CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

evapPhaseChangeFOAM

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree7Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   January 14, 2013, 15:02
Default
  #21
New Member
 
Martin Lubej
Join Date: Jul 2012
Posts: 5
Rep Power: 4
mlubej is on a distinguished road
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
mlubej is offline   Reply With Quote

Old   January 23, 2013, 11:46
Default solver
  #22
New Member
 
Kris Coumans
Join Date: Oct 2012
Location: Netherlands
Posts: 1
Rep Power: 0
KrisC is on a distinguished road
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!
KrisC is offline   Reply With Quote

Old   January 29, 2013, 15:10
Default
  #23
New Member
 
Kuan-Lin Lee
Join Date: Jul 2012
Posts: 2
Rep Power: 0
Jmmymiller is on a distinguished road
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
Jmmymiller is offline   Reply With Quote

Old   April 12, 2013, 01:16
Default
  #24
New Member
 
Zeng Qingyun
Join Date: Oct 2012
Posts: 8
Rep Power: 4
Chandler is on a distinguished road
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!
Chandler is offline   Reply With Quote

Old   July 4, 2013, 13:19
Default
  #25
Member
 
Mohammad Bahreini
Join Date: Dec 2012
Posts: 35
Rep Power: 4
mecman is on a distinguished road
Dear Nima Sam,
Could you please send me the solver?m.Bahreini1990@gmail.com
Regards,
mecman is offline   Reply With Quote

Old   July 11, 2013, 06:01
Default source code
  #26
New Member
 
Anastasios
Join Date: Mar 2009
Posts: 25
Rep Power: 7
ageorg is on a distinguished road
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
ageorg is offline   Reply With Quote

Old   July 17, 2013, 01:59
Default
  #27
New Member
 
Join Date: Jun 2013
Posts: 16
Rep Power: 3
wond is on a distinguished road
Quote:
Originally Posted by houkensjtu View Post
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
wond is offline   Reply With Quote

Old   July 25, 2013, 05:45
Default code
  #28
New Member
 
Join Date: Jul 2013
Posts: 1
Rep Power: 0
Annabel is on a distinguished road
Quote:
Originally Posted by nimasam View Post
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
Annabel is offline   Reply With Quote

Old   July 25, 2013, 07:45
Default
  #29
Member
 
Alexander
Join Date: Mar 2009
Posts: 37
Rep Power: 7
sahas is on a distinguished road
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.
nimasam likes this.
sahas is offline   Reply With Quote

Old   July 29, 2013, 07:32
Default
  #30
New Member
 
Y.LANCE
Join Date: Feb 2013
Posts: 7
Rep Power: 3
Yoann is on a distinguished road
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.

Last edited by Yoann; July 29, 2013 at 08:35.
Yoann is offline   Reply With Quote

Old   July 29, 2013, 09:43
Default
  #31
Member
 
Alexander
Join Date: Mar 2009
Posts: 37
Rep Power: 7
sahas is on a distinguished road
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.
Attached Files
File Type: gz evapInterFoam.gz (25.8 KB, 74 views)
sahas is offline   Reply With Quote

Old   July 29, 2013, 09:51
Default
  #32
New Member
 
Y.LANCE
Join Date: Feb 2013
Posts: 7
Rep Power: 3
Yoann is on a distinguished road
Thank, I look that
Yoann is offline   Reply With Quote

Old   August 8, 2013, 03:29
Default
  #33
New Member
 
Josefine Wilms
Join Date: Sep 2009
Posts: 8
Rep Power: 7
Fine is on a distinguished road
Hi Nima Sam

Could you please send me your code? My email: fine.wilms@gmail.com
Fine is offline   Reply With Quote

Old   August 8, 2013, 03:46
Default
  #34
Senior Member
 
Mohammad Shakil Ahmmed
Join Date: Oct 2012
Location: AU
Posts: 109
Rep Power: 5
ahmmedshakil is on a distinguished road
Hi sahas,
Can you upload any test case of your solver?
ahmmedshakil is offline   Reply With Quote

Old   August 8, 2013, 11:04
Default Help
  #35
New Member
 
Join Date: Dec 2011
Posts: 11
Rep Power: 5
adampaya is on a distinguished road
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?
adampaya is offline   Reply With Quote

Old   August 8, 2013, 11:22
Default
  #36
Member
 
Alexander
Join Date: Mar 2009
Posts: 37
Rep Power: 7
sahas is on a distinguished road
Quote:
Originally Posted by ahmmedshakil View Post
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.
sahas is offline   Reply With Quote

Old   August 10, 2013, 05:34
Default
  #37
New Member
 
Jose
Join Date: Oct 2012
Posts: 6
Rep Power: 4
llidito is on a distinguished road
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ť.
llidito is offline   Reply With Quote

Old   September 6, 2013, 03:31
Default how add source term and T.eqn for evap
  #38
New Member
 
sasan
Join Date: Sep 2013
Posts: 28
Rep Power: 3
Sasy is on a distinguished road
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...
Sasy is offline   Reply With Quote

Old   September 6, 2013, 05:41
Default
  #39
Member
 
Alexander
Join Date: Mar 2009
Posts: 37
Rep Power: 7
sahas is on a distinguished road
Quote:
Originally Posted by Sasy View Post
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?
sahas is offline   Reply With Quote

Old   September 6, 2013, 05:53
Default
  #40
New Member
 
sasan
Join Date: Sep 2013
Posts: 28
Rep Power: 3
Sasy is on a distinguished road
Quote:
Originally Posted by sahas View Post
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....
Sasy is offline   Reply With Quote

Reply

Tags
boiling, evaporation, interfoam, phase change

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 20:46.