ABSTRACTArtificial compressibility method is extended to solve the low Mach variable density Navier–Stokes equations for the simulation of axisymmetric laminar radiative diffusion flames. A combustion model based on the generalized state relationships for species mass fractions and temperature is adopted and the radiation heat loss is estimated as a fraction of heat of combustion. Conventional finite difference method along with a total variation diminishing scheme is used for the spatial discretization on nonuniform half-staggered grid layouts. To reach the steady-state condition, a three-step low-storage explicit Runge–Kutta method is used. The accuracy of the proposed method is reported for confined and unconfined diffusion flames.