[nico.vanesch]

Dear all,

I am trying to reproduce/validate the result of the paper of Khater et al (2019):

https://www.nature.com/articles/s41598-019-40069-9

In this research, an incompressible non-isotherm multiphase solver was used (using ANSYS FLUENT) with a variable density, viscosity and surface tension. I've already compiled the interTempFoam succesfully and it runs smoothly.
So I'm looking for something like this:




I already looked into integrating the thermoPhysical library into the solver, but these are incompatible since these are only available for compressible solvers.
If anyone can lead me the way to the answer, I would be very grateful.

Best,

//UPDATE
never mind, found it myself.
You can use:
http://www.tfd.chalmers.se/~hani/kur...nFoam%20v2.pdf
to make the viscosity temperature dependent.

Then, you can use the temperatureDependent function to make a temperature dependent surface tension:
sigma
{
type temperatureDependent;
sigma polynomial
((a 1)
(b 1)
);
with a and b as in the equations above.

Then, for the density I recalculated the polynomial above to fit to the bousineqq approximation and used the code from the beginning of this post. The only thing that I know still have to figure out, is to make the Prandtl number to be dependent of the temperature dependent viscosity. Anyone have suggestions on this?

[tonnykz]

Hi Nico,
I have also added temperature equation to OF-v1812 using information provided in this sub Adding the Energy Equation to interFoam (OF 2.4.0).
But results for temperature equation diverging.
Which version of OpenFOAM you have used?
Is there anything trick to make it work on simple damTempBreak tutorial?
Thank you,
Best regards

Tonnykz

[AUbuntu]

Hello Guys,


I see that this thread is very old. I am interested to know if any solver has been developed in OpenFOAM to be used for solving energy equation for interFoam during this 10 years period similar to interTempFoam.
As I followed this thread and used different schemes such as harmonic and vanLeer scheme to suppress the diffusion at the water-air interface. vanLeer gives me the best results however there is still some unrealistic increase in temperature at the water-air interface. Does anyone have any idea how to completely solve this issue?


Regards

[tonnykz]

Originally Posted by AUbuntu

Hello Guys,


I see that this thread is very old. I am interested to know if any solver has been developed in OpenFOAM to be used for solving energy equation for interFoam during this 10 years period similar to interTempFoam.
As I followed this thread and used different schemes such as harmonic and vanLeer scheme to suppress the diffusion at the water-air interface. vanLeer gives me the best results however there is still some unrealistic increase in temperature at the water-air interface. Does anyone have any idea how to completely solve this issue?


Regards

Hi AUbuntu,
I have switched to the formulation of energy equation as in this article https://github.com/MahdiNabil/CFD-PC
I don't see same divergence with using it.
Best,
Tonnykz

[AUbuntu]

Originally Posted by tonnykz

Hi AUbuntu,
I have switched to the formulation of energy equation as in this article https://github.com/MahdiNabil/CFD-PC
I don't see same divergence with using it.
Best,
Tonnykz

Hello Tonnykz,

Many thanks for your reply.
So, based on my understanding in the paper they used the energy equation in the form of enthalpy (h)
Is this the only change that you have considered?
I saw they also mentioned using different surface-tension force model. Could you please also kindly let me know if you have also considered that?

regards

[tonnykz]

Originally Posted by AUbuntu

Hello Tonnykz,

Many thanks for your reply.
So, based on my understanding in the paper they used the energy equation in the form of enthalpy (h) instead of internal energy (u).
Is this the only change that you have considered?
I saw they also mentioned using different surface-tension force model. Could you please also kindly let me know if you have also considered that?

regards

Hi AUbuntu,
Yes, they are using enthalpy and that what I have implemented.
But they have thermal phase change implemented apart from special surface-tension models. Latter i did not consider. But also adopted thermal phase change models they provide. You can use noPhaseChange model if you want to model transient energy equation transport due to convection and diffusion.
Best,
Tonnykz