[Barlo]

Hello everyone,
I'm trying to simulate the injection of methane into a constant volume combustion chamber with the solver XiFoam, using non-homogeneous mixture as a thermodynamic model. I found a problem: (that is established since the early moments of simulation) a strange behaviour of temperature (see picture below). Initially I thought it was due to problems of mesh or the stability of the discretization schemes, but, after some simulations on a hexahedral homogeneous grid, and using bounded schemes the problem is not solved yet. So, perhpas the problem lies in how the solver can extracts temperature from enthalpy (density is also affected obviously), but I'm not able to understand how cp and R are treated as functions of the mixture fraction (if it is actually taken into account in XiFoam). Does anyone have an idea about it? Is there anyone who has already ran into this problem? I enclose all files as set.
Thanking you in advance,
best regards,
Lorenzo Bartolucci

[Barlo]

Any suggestion?