
[Sponsors] 
January 14, 2013, 14:02 

#21 
New Member
Martin Lubej
Join Date: Jul 2012
Posts: 5
Rep Power: 13 
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 

January 23, 2013, 10:46 
solver

#22 
New Member
Kris Coumans
Join Date: Oct 2012
Location: Netherlands
Posts: 1
Rep Power: 0 
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! 

January 29, 2013, 14:10 

#23 
New Member
KuanLin Lee
Join Date: Jul 2012
Posts: 2
Rep Power: 0 
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 

April 12, 2013, 00:16 

#24 
New Member
007
Join Date: Oct 2012
Posts: 9
Rep Power: 13 
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! 

July 4, 2013, 12:19 

#25 
Member
Mohammad Bahreini
Join Date: Dec 2012
Posts: 36
Rep Power: 13 
Dear Nima Sam,
Could you please send me the solver?m.Bahreini1990@gmail.com Regards, 

July 11, 2013, 05:01 
source code

#26 
Member
Anastasios
Join Date: Mar 2009
Posts: 34
Rep Power: 17 
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 

July 17, 2013, 00:59 

#27  
New Member
Join Date: Jun 2013
Posts: 18
Rep Power: 12 
Quote:
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 

July 25, 2013, 04:45 
code

#28  
New Member
Join Date: Jul 2013
Posts: 1
Rep Power: 0 
Quote:
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 

July 25, 2013, 06:45 

#29 
Member
Alexander
Join Date: Mar 2009
Posts: 49
Rep Power: 17 
Hello everyone!
If you are interesting in VOFsolver 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 nonconservative 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 onecomponent 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 VOFinterface 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,1E8) ); 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. 

July 29, 2013, 06:32 

#30 
New Member
Y.LANCE
Join Date: Feb 2013
Posts: 7
Rep Power: 13 
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 07:35. 

July 29, 2013, 08:43 

#31 
Member
Alexander
Join Date: Mar 2009
Posts: 49
Rep Power: 17 
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 nonconservative form. But during equation transformation some additional terms arise: rho*(Cp1Cp2)*T*Y/rho1Cp*(rho1rho2)*T*Y/rho1. Fortunately these terms are negligible comparing with L*Y so they can be omitted. 

July 29, 2013, 08:51 

#32 
New Member
Y.LANCE
Join Date: Feb 2013
Posts: 7
Rep Power: 13 
Thank, I look that


August 8, 2013, 02:29 

#33 
New Member
Josefine Wilms
Join Date: Sep 2009
Posts: 10
Rep Power: 16 

August 8, 2013, 02:46 

#34 
Senior Member
Mohammad Shakil Ahmmed
Join Date: Oct 2012
Location: AUS
Posts: 137
Rep Power: 14 
Hi sahas,
Can you upload any test case of your solver? 

August 8, 2013, 10:04 
Help

#35 
New Member
Join Date: Dec 2011
Posts: 12
Rep Power: 14 
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? 

August 8, 2013, 10:22 

#36 
Member
Alexander
Join Date: Mar 2009
Posts: 49
Rep Power: 17 

August 10, 2013, 04:34 

#37 
New Member
Jose
Join Date: Oct 2012
Posts: 6
Rep Power: 13 
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é. 

September 6, 2013, 02:31 
how add source term and T.eqn for evap

#38 
New Member
sasan
Join Date: Sep 2013
Posts: 28
Rep Power: 12 
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... 

September 6, 2013, 04:41 

#39  
Member
Alexander
Join Date: Mar 2009
Posts: 49
Rep Power: 17 
Quote:
Hello! Have you seen my post? 

September 6, 2013, 04:53 

#40  
New Member
sasan
Join Date: Sep 2013
Posts: 28
Rep Power: 12 
Quote:
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.... 

Tags 
boiling, evaporation, interfoam, phase change 
Thread Tools  Search this Thread 
Display Modes  

