[ssyadav]

Hello Nimasam

Kindly send me the OpenFoam based code to simulate boiling with OpenFoam. My email id. is ssyadav25@gmail.com

Thanks in advance.

Shyam Sunder

[elisabet]

Nima and Alexander,

I'm a little bit confused about the state of the art of your solvers. I'ld really appreciate you could send me your codes and, if possible, comments. I'll try to understand the differences and limitations.

I'm currently interested in steam condensation, so any suggestion is welcome!
Thank you in advance,

Elisabet (elisabet.masdelesvalls@gits.ws)

[sahas]

Hello, Elisabet!

Currently I have no deal with simulation of evaporation/condensation using VOF method. Some information about evaporation model implementation can be found in this thread. Still remaining problem is taking correctly into account the latent heat. As mentioned in the literature, the possible solution is to introduce weight function in such a way that the integral from the volume latent heat source multiplied by this weight function gives the value of correct latent heat power. Unfortunately, I do not have much information about this.

[serena_bears]

Hi Nimasan,

could you please send me your final code? Now I'm working on droplet condensation problem. Thank you very much.

My email is LONG0032@e.ntu.edu.sg

[tju_shq]

Hi Nimasan,

could you please send me your final code? Now I'm working on droplet condensation problem. Thank you very much.

My email is shiquan@tju.edu.cn

[asla9796]

Hello, I am new to CFD and is trying to simulate the same problem of droplet evaporation. This link posted in the thread which contains the files for solver appears to be empty now. Please do look into this and do share the code again if possible.

My mail ID is mm14b007@smail.iitm.ac.in

Thank you for your time

[gaza]

Hi
for now my code does work
I will share my code when it will work
see to OpenFOAM 1606 there is interEvaporatingCondensingFoam

[Ank09]

Hi Nimasam!
I'm working on the droplet evaporation topic and I'm very interested in your code. Could you please send me a copy of your final code?
Thanks a lot! My best wishes!
My email: andyking0928@foxmail.com

[GauthamKrishnan]

Hi Nimasam!
I'm working on the droplet evaporation and I'm very interested in your code. Could you please send me a copy of your final code+ case file
Thanks a lot! My best wishes!
My email: [email]gkriz66@gmail.com

[nimasam]

Dear foamer

The latest version of my solver (phaseChangeHeatFoam), several test cases and published papers are available in github.

https://github.com/NimaSam/phaseChangeHeatFoam/


please inform me about possible bugs.


Best Regards

[ssyadav]

Dear Nima

Thank you very much for posting your solver. I tried compiling it in OpenFOAM 2.2.2, I am getting some errors which are as follows:

1) smoothInterfaceProperties.C:27:49: fatal error: alphaContactAngleFvPatchScalarField.H: No such file or directory
#include "alphaContactAngleFvPatchScalarField.H"

I could resolve this error by changing the path in the "smoothInterfaceProperties/Make/options" file

2) /usr/bin/ld: cannot find -ltwoPhaseInterfaceProperties

This library is included in "phaseChangeHeatFoam/Make/options" but I could not locate it my OpenFOAM installation.

Please help.

Also how I can extend your solver to include Conjugate heat transfer in a solid, particularly the boundary condition at the solid-fluid interface, i.e. the incompressible version of "compressible::turbulentTemperatureCoupledBaffleMi xed"

Thanks.

Shyam

[nimasam]

Dear Shyam

1- for your second question, i guess the name of library has been changed in newer version of OpenFOAM, i dont have OpenFOAM-2.2.0 to check it my self. but i guess it became libtwoPhaseProperties or libinterfaceProperties

2- you can see this paper

[Elham]

Originally Posted by nimasam

Dear Shyam

1- for your second question, i guess the name of library has been changed in newer version of OpenFOAM, i dont have OpenFOAM-2.2.0 to check it my self. but i guess it became libtwoPhaseProperties or libinterfaceProperties

2- you can see this paper

Dear Nima,

I cannot compile your code due to the following error:

In file included from phaseChangeHeatFoam.C:57:0:
createFields.H: In function ‘int main(int, char**)’:
createFields.H:106:72: error: no matching function for call to ‘Foam::smoothInterfaceProperties::smoothInterfaceProperties(Foam::volScalarField&, Foam::volVectorField&, Foam::phaseChangeTwoPhaseMixture&)’

Do you have any idea to fix it?

Cheers,

Elham

[nimasam]

1-which version of openfoam, do u use?
2-you should first compile smoothInterfaceProperties using following command:

wmake libso

then you can compile the solver

[Elham]

Originally Posted by nimasam

1-which version of openfoam, do u use?
2-you should first compile smoothInterfaceProperties using following command:

then you can compile the solver

The library is compiled but still have the error. I use OpenFOAM 2.3.0.

[nimasam]

Dear Elham it's been written for OpenFOAM_220, for upper versions of openfoam, you may need to change some of libraries name, because the basic solver which is interfoam is changed a lot through openfoam 2.2 to openfoam 3.0

[Elham]

Dear Nima,

I have read your paper "The evaluation of the diffuse interface method for phase change simulations using OpenFOAM" and the phaseChangeHeatFoam code. To derive Sp you have used the right hand side of equation 18 but I am still wondering how you have derived su , as part of source term?
I would appreciate if you give me some clues.

Cheers,

Elham

[nimasam]

Several points to clarify coding:
1)as i remember the formulation of MULES is like that:

ddt(alpha1)+div(U,alpha1)+Sp*alpha1+Su=0

2- convection term in equation 18, then:

U & grad(alpha1)=div(U,alpha1)-alpha.div(U)

3-the mass transfer source term is multiplying with void fraction, means:

mDotBoiling=mdot*alpha1;
mDotCondensation=mdot*(1-alpha1);

[Elham]

Originally Posted by nimasam

Several points to clarify coding:
1)as i remember the formulation of MULES is like that:

2- convection term in equation 18, then:

3-the mass transfer source term is multiplying with void fraction, means:

Dear Nima,

I am confused about deriving the source terms. Equation 18 is as following:

ddt(alpha1)+U& grad(alpha1) +grad(alpha1(1-alpha1)Uc) = -mDot(1/rho1-alpha1(1/rho1-1/rho2)))

and in eq 16:

div (U) = mDot(1/rho2-1/rho1)

substituting eq 16 into eq 18 and rearranging:

ddt(alpha1) + grad(U,alpha1) + grad (alpha1(1-alpha1)Uc) = -mDot/rho1

So we should have:

sp = 0 and su = mDot/rho1

But in alphaEqn.C:

sp=mDot(1/rho1-alpha1(1/rho-1/rho2))
su=divU*alpha1 + mDot(1/rho1-alpha1(1/rho-1/rho2))

Thnaks for any help.

Elham

[nimasam]

Dear Elham
Please re-read my points, you will figure out how it is derived.
equation 18:

dt(alpha1)+U& grad(alpha1) +grad(alpha1(1-alpha1)Uc) = -mDot(1/rho1-alpha1(1/rho1-1/rho2)))

lets consider S=-mDot(1/rho1-alpha1(1/rho1-1/rho2)))
then source terms for boiling and condensations would be:
Sb=S*alpha1
Sc=S*(1-alpha1)
Also i mentioned in previous post that:

U & grad(alpha1)=div(U,alpha1)-alpha1*div(U)

then right hand side of equation 18 becomes:

dt(alpha1)+div(U,alpha1)-alpha.div(U) +grad(alpha1(1-alpha1)Uc)

so alpha1*div(U) is appearing here
now consider boiling and condensation separately,
for boiling the source term is Sb=S*alpha1
so

su=alpha1*div(U)
sp=S

for condensation Sb=S*(1.0-alpha1), then:

su=alpha1*div(U)+S
sp=-S