reactingFoam, correct solvers / discretization schemes
 

Hi everyone,

I am using reactingFoam to solve a flow inside a circular tube: injection of a flow of air at ambient temperature through a centre tube, with a surrounding co-axial flow of hot gases. The simulation is in laminar regime and i use chemical species O2, CO2, H2O, and N2 with no chemical reactions.

The species and enthalpy seem to converge, however the axial velocity and pressure residuals oscillate greatly, and I am trying to solve this but I am not being able to.

Can anyone advise me on this? It seems like the default solvers / discretization schemes that reactingFoam uses are not the correct ones for flows highly dependent on the velocity field.

Any help is much appreciated!

Thanks in advance :)