evapPhaseChangeFOAM
 

i have developed a solver based on interFOAM, to solve energy equation and besides consider mass transfer between two phases, this solver works correctly for one dimensional case (stephan phase change problem)
but for two case studies, temperature at interface behaves strangely
now any suggestion, cooperation or idea will be helpful

P.S
developed files and case studies are available in attachment
+
some descriptions :) can be found here:

http://www.4shared.com/document/-eBG...OF_method.html

Hi.nimasam,im intersted in this problem,too.

have you solved this problem already?
Would you tell me the way to resolve this problem?

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

Hello Nimasam,
i have the same problem: phase change and evaporation due to heat transfer...
OF not have a default solver for this.
have you any suggestion for me to solve the problem in your solver?

Thank you

my solver solves evaporation with VOF interface tracking method, if you look at here:
http://www.cfd-online.com/Forums/blo...-openfoam.html
you will find some little description about how does it works and present situation of solver.

Thanks
 

Hi nimasam,
thanks for the reply.
do you share the correct solver code?
i have some problems with c++ programmation.

Hi
send me your email

many tanks
 

thanks nimasam for your understanding...
my email is: seranto2001@yahoo.it

Hello nimasam,
please re-send me the solver code,
the mail not arrived....
Thanks
Antonello

Hi Nima Sam

I need to simulate boiling flow in OF but I could not find any solver for that. From this forum I found that you have developed your own code. Is it for boiling flow? May I use it?

Thanks

Hi, nimasam. I have always confused the implementation in VOF model. Can you send me your final code. My email is nygbook@gmail.com. Thank you!

Hi Nima Sam.

I'm also working on this topic and I would also be interested in having your source code, if possible, in order to test it on experimental data I have.

My email adress is julien.commanay at free.fr

Regards.

Dear Nima Sam,
Could you please send me the solver? my email id is mrityunjay111@gmail.com.
Thanks in advance
MK

Dear Nima Sam,
Could you please send me the solver to zcarral@yahoo.es
Is it applied to bubble growing due to evaporation?
Thanks in advance

Dear Nima Sam,

I'm studying the instability of a boiling front and I'm really interested in the source code you developed. Could you please send me the OF solver ?
My email is psapin@imft.fr

Thanks,

Paul

Hi nimasam,
Will you please send me your code, my email: ahmmedshakil@hotmail.com

I would like to take a look at the code too ( albert.passalacqua@gmail.com)

PS. In the description you say you need level set. In reality, the level set alone would not be enough, due to conservation errors. There are combined level set - vof methods, and but frankly they become easily too complicated, they can lead to inconsistent representations, and and they are expensive (in some case the conservation problem is not resolved, and typically you need adaptive refinement anyways). I would suggest you take a look at the work of Popinet (Gerris code). They developed a geometric VOF with height functions/paraboloid interpolation to calculate the curvature. Combined with adaptive grid refinement it seems more promising than coupling level set and vof :-)

Best,

Hi nimasam,

neoxwolf@yandex.com

I would be pleased if you could send your code. I will compare with my experimental results.

thank you

Hello, Nima Sam.

Would you be so kind to send your code to smirnovskysaha@gmail.com?
Thanks in advance!

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!

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

solver
 

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!

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

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!

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

source code
 

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

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

code
 

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

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.

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.

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.

Thank, I look that

Hi Nima Sam

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

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

Help
 

Hi everybody
I have a primary problem:mad:
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?:)

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

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é.

how add source term and T.eqn for evap
 

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...:):confused::)

Hello! Have you seen my post?

Hi Alexsander
Thank you for Reply
I see now...:D:D
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....:confused: