[gentodin]

Hi all, I am trying to implement my temperature equation to OpenFOAM.

I have an additional mass update that changes the cells at every timestep from gas to solid phase hence a change in density.

After implementing what I think is correct, the temperature reading in the gas phase goes higher and higher until the simulation cannot solve anymore due to floating-point error, (is this what u call diverge?).

Nevertheless, I would like to check with you guys if you happen to know what are things to look out for when implementing such a complex system?

[gentodin]

Hi Anybody knows what are the issues?

For example,
I am updating
forAll(alpha1,celli)
{

if(some condition)
{
alpha1[celli] = 1.0;
T[celli] = 1900;

}
}


What are the things I need to lookout for when I update the Volume fraction alpha1?