[Aritra Roy Choudhury]

Dear community,
I am running a reacting simulation for premixed H2-air gas. I am using k-epsilon turbulence model, PaSR model with Cmix = 0.001, and ode solver for chemistry solving.
When running in parallel (in reactingFoam in OpenFOAM-v2312) I am observing discontinuous mass fraction values at the different processor domains. I also tried using OpenFOAM-v2112 with no use. This behaviour is not there when running in serial. This behaviour is only visible for some minor species and not for U, P and T fields. I have attached the contours for reference. Please have a look at them and if anyone has faced the same issue please let me know how did you fixed it.
Thank you

[GeekCFD]

Hi Aritra,
Have you performed the Grid Independence Study for your mesh? The mesh is likely the issue.
The other issue can be a 'relatively large' deltaT value in controlDict.

Can you please share the case setup for your case?
It is otherwise very difficult to comment on anything in particular.

-----------------
Thanks,
GeekCFD