An algorithm for solving variable density incompressible flows is presented. This algorithm is implemented on a finite volume based fully implicit solver employing a non-staggered hexahedral structured grid arrangement. Momentum interpolation is used for pressure-velocity coupling and the resulting Poisson equation for pressure is modified to a Helmholtz equation for better convergence. This algorithm is used for solving problems involving the variation of density due to temperatu re, flu id mixing and turbulent combustion flows. Copyright © 2007 by ASME.