For the simulation of low Mach number reactive flows with significant density changes a modified lattice-BGK model in combination with the conventional convective-diffusion solvers for equations of temperature and species is proposed. Together with boundary-fitting conditions and local grid refinement the scheme enables the accurate consideration of low Mach number combustion in complex geometry as the flows around porous burners or droplets combustion.