[mlubej]

Hello nimasam.
I have tried to develop my solver for boiling for dimensioning a heat exchanger with phase change for academic purposes. I would be really glad if i could look into your final code, I have problems with temperature distribution in steam phase...

my email is mlubej@siol.net, I am a phd student of chemical engineering.
Best regards, Martin

[KrisC]

Dear Nima Sam,

I'm working on the optimization of an evaporator for liquified petroleum gas. For this, I need to model the evaporation of LPG as a consequence of heat transfer.

Can you please sent me your solver?
My email is: kcoumans (at) hotmail.com

Thank you very much!

[Jmmymiller]

Hi nimasam:

I am working on the heat pipe simulation. Is you code can deal with condensation also?
Anyway, I would like to see the code! Please send it to me.
jimmymiller38200@gmail.com

Thanks

[Chandler]

Hi, nimasam! I am also interested in this problem. but i have little progress.
would you please send me your final code and case? and i would like to read the your paper related to this solver.
my email: zengqingyun89@163.com
Thank you!

[mecman]

Dear Nima Sam,
Could you please send me the solver?m.Bahreini1990@gmail.com
Regards,

[ageorg]

Dear nimasam,

I am also working in boiling of refrigerants in microchannels. I would also be very interested in your code. My email is anastasios.georgoulas@gmail.com. If you send me the code please send me also your full details (name, surname, affiliation etc) for citation and acknowledgement purposes in any potential future publication of my research group.

Thank you very much in advance,

ageorg

[wond]

Originally Posted by houkensjtu

Hi nimasam!

I had tried your code above in the attached file.
Actually because I am using of 2.1.x, so what I did is to read your code and apply all the changes by myself to interFoam. After some effort I compiled the solver successfully.
Unfortunately, I found when I run the StephanProblem, there were two issues:

1. temperature went extremely high near the interface, which caused extremely high velocity which is surely not correct.
2. pressure inside the domain went extremely high and cause some unrealistic velocity

would u please send me a copy of your final code? I am now stuck and don't know where was the problem.

my email:

houkensjtu@gmail.com

Thank u!

Hi, Houken

haver you find out the reson of the high temperature and pressure? I have carried out the same model by commercial software and the results seemed to be unreasonable which is same with yours.

thanks!

Wond

[Annabel]

Originally Posted by nimasam

i solved the problem, if u like to received final code send me ur email

hey nimasam,
i am newbie in openfoam. i would like to get the code of evaporation problem, because i also have a urgent one
yaqileee@gmail.com

[sahas]

Hello everyone!

If you are interesting in VOF-solver with evaporation I can give some useful (I hope ) information because we realized it for the problem of drop evaporation.

First of all we chose interPhaseChangeFoam as the starting point and add temperature equation there. To avoid problems with divergence during numerical solving one may write temperature equation in non-conservative form in the following way:

fvScalarMatrix TEqn
(
rhoC*fvm::ddt(T) + rhoC*fvm::div(phi, T) - fvm::laplacian(K , T) + L*Y
);

Here:
volScalarField K = alpha1*K1 + (scalar(1)-alpha1)*K2;
volScalarField rhoC = alpha1*rho1*C1 + (scalar(1) - alpha1)*rho2*C2;

K1, K2, C1 and C2 - heat transfer and heat capacity coefficients for both phases (by the way we consider only one-component gas phase - air without steam, it is lack of our model).
L - latent evaporation heat;
Y - evaporated mass in the volume (will be defined below).

Next, it is need to modify equation for alpha1 (alphaEqn.H). In the solver interPhaseChangeFoam there are two source terms named Su and Sp. Su answer for the source term which calculated explicitly, Sp - implicitly, so the next equation is solving:
(... left part for alpha1 equation ...) = Sp*alpha1 + Su.
In our model it is necessary to calculate evaporation implicitly, so Su should be equal 0 and Sp should be calculated as Sp = -Y/rho1/alpha1.

The main question is about Y calculation.
Y should be "evaporated mass in the volume" whereas when we consider evaporation we have only "evaporated mass from the square" (named S).
It can be shown that:
Y = S*|grad(alpha1)|

Quantity S is evaluating according to the evaporation model. We use model of B. Abramzon and W.A. Sirignano and define S as:

S = D*Sh*rho_gas*ln(1+Bm)/R,

where D - diffusion coefficient, Sh - modified Sherwood number, rho_gas=rho2 - gas density, Bm - mass Spalding number (it is a function of temperature), R - radius of VOF-interface curvature (for the model details you can ask me or look at articles ). The main problem is the calculation of R, let me know if you need more information about it.

So, in alphaEqn.H you need the next code:

volScalarField vDotvAlphal (-Y/rho1/max(alpha1,1E-8) );
volScalarField Sp
(
IOobject ("Sp", runTime.timeName(), mesh), vDotvAlphal
);

volScalarField Su should be equal 0.

Unfortunately, our code is written in "ugly" manner and I do not want to paste it here (I am ashamed for it). But if you really need it simply write here and I post it in a separate message.

[Yoann]

Hello sahas.

I work on VOF solver with evaporation and I'm interested by your code.
If is possible I would like to see your code.

Thanks.

[sahas]

Hello Yoann!

I attach tar.gz archive with the solver. I should note that there is a directory phaseChangeTwoPhaseMixtures which is superfluous for present solver but I feel too lazy to clean the dependencies

The new/modified files are: interPhaseChangeFoam.C, TEqn.H, alphaEqn.H and createFields.H.

Also I beg my pardon if the problems with solver running arise since the man responsible for the code is absent now. Nevertheless do not hesitate to ask me about the code

Apropos, I should make some notes about temperature equation.
Right equation is in conservative form:
fvScalarMatrix TEqn
(
fvm::ddt(rhoC,T) + fvm::div(rhoC*phi, T) - fvm::laplacian(K , T) + L*Y
);

But as I mentioned before there are difficulties with its numerical solving. So it is better to use non-conservative form. But during equation transformation some additional terms arise: -rho*(Cp1-Cp2)*T*Y/rho1-Cp*(rho1-rho2)*T*Y/rho1. Fortunately these terms are negligible comparing with L*Y so they can be omitted.

[Yoann]

Thank, I look that

[Fine]

Hi Nima Sam

Could you please send me your code? My email: fine.wilms@gmail.com

[ahmmedshakil]

Hi sahas,
Can you upload any test case of your solver?

[adampaya]

Hi everybody
I have a primary problem
I installed the PGI visual Fortran successfully but
i don't know how could i run a program you know
first it should be compiled then Run it
can anybody explain the sequence of steps needed to Run a program in PGI Visual Fortran I mean?

[sahas]

Originally Posted by ahmmedshakil

Hi sahas,
Can you upload any test case of your solver?

I can do it only in a few weeks because the man who is responsible for the code is absent now.

[llidito]

Hi Nima Sam,

I would like to take a look at the code. Could you please send it to me?

michaelcrespos@hotmail.com

Thanks and nice Saturday,
José.

[Sasy]

Hi every body
i want simulate evaporation in vertical tube.i use interFoam solver,and i know i shoud adding T.eq and mass source...but i dont have idea...
any body know?gudie me plz...

[sahas]

Originally Posted by Sasy

Hi every body
i want simulate evaporation in vertical tube.i use interFoam solver,and i know i shoud adding T.eq and mass source...but i dont have idea...
any body know?gudie me plz...

Hello! Have you seen my post?

[Sasy]

Originally Posted by sahas

Hello! Have you seen my post?

Hi Alexsander
Thank you for Reply
I see now...
is this solver that you are attach,simulate the evaporation at the interface between fluid and gas by increasing Temp?is this solver have source term for phase change?
i dont work with interPHasechangeFoam solver....but i think this solver simulate Cavitation with reduce pressure....