A penalization method for the simulation of weakly compressible reacting gas-particle flows with general boundary conditions

•Weakly compressible approximation combined with volume penalization method.•Dirichlet, Neumann and Robin type boundary conditions at particle surface.•Model validation & application to finite or infinite rate volume & surface reaction.•Efficient combination external and intraparticle field calculation.•Effects of compressibility systematically evaluated.