[prinux]

Hello Foamers,

I'm relatively new to OpenFOAM and have been tasked to extend interPhaseChangeFoam for solving multiphase simulation with heat and mass transfer. I have added an EEqn.H to the solver files and it compiles without issues.

My question is : How do I initialize temperatures for the different phases in my simulation. I am able to add temperature as a boundary condition to the walls of the domain.

Previously, I had tried to do the same using setFieldsDict, but all it does is initialize the temperature field for the phases, but doesn't get solved for the next time step and the temperature field reverts to the internal field value defined in 0/T.orig file. Do I have to add a T.liquid or a T.vapour file in the 0 folder of my case files, like in multiPhaseEulerFoam, or have I missed something while extending the solver?

I'm using OpenFOAM v8 on Ubuntu 20.04.

Thanks,
Prithvi