[Filippo_azz]

Good evening, everyone. I have been working with multiphase simulations using interFoam. Currently, I am attempting to incorporate temperature field analysis into my simulations. I have successfully integrated the temperature (T) equation into interFoam (specifically working with OpenFOAM version 8) quite effortlessly. However, re-running the simulation to include these changes is quite time-consuming.
In my specific case, the fluid properties are not affected by the temperature field. Therefore, I am exploring ways to calculate the temperature field using the velocity (U) and volume fraction (Alpha) fields already evaluated in the previus simulations. To achieve this, I've adapted the scalarTransportFoam solver to read the alpha field. This adaptation is necessary to accurately compute the property values. I've modified the solver's equation to the following:

c*(fvm::ddt(rho, T) + fvm::div(rhoPhi, T)) - fvm::laplacian(k, T)

However, I am encountering an issue. Although the simulation initiates, it yields non-physical and odd results, suggesting that I might be missing a crucial element in my approach.
Could anyone provide insights or guidance on this matter? Any assistance would be greatly appreciated. Thank you very much.