[Thamali]

Hi,
I am doing a multiphase simulation as a single phase.
This is about gasification/combustion of a single particle. And I am resolving inside the particle too.
My particle is not moving. Therefore, it only involves in enthalpy by reaction and changes to the gas phase species and mass by reactions.
I assume the particle and gas inside the particle is in thermal equilibrium.
I use a binary function(gamma or whatever which is '0' and '1') to identify the two different phases in governing equations.
My problem is when I use OpenFoam thermo packages and if I solve enthalpy of gas phase equation (by considering solid phase terms as sources) the equation does not converge.
Therefore, I tried to solve temperature equation for thermo.T() and it does not captured by thermo package since I do not solve enthalp and it considers the 'enthalpy' change is zero and thermo.correct() gives same constant temperature over and over during simulations.
I thought I can do a trick like solve the 'T_otherEqn' and calculate

thermo.he()=cp*T_other;

and then thermo.T() will change based on the he().

thermo.T() should be =T_other.

But, thermo.T() calcluated by correct() function is different from T_other.

Could somebody explain this to me?
Or what are my other options plsss?
If I don't use available thermo packages, coding willl become complicated and lot of work to do.

Thanks a lot in advance..!