A directional ghost-cell immersed boundary method for low Mach number reacting flows with interphase heat and mass transfer

•General Robin (Neumann) BC are accurately imposed and satisfied.•Neither polynomial assumptions or complex interpolations are required.•Low Mach number variable-density flows are considered.•Reaction-induced interfacial phenomena can be accurately resolved.•Second-order accuracy is preserved for all the implementations.