# evapPhaseChangeFOAM

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

September 25, 2017, 12:02
#101
Senior Member

Elham
Join Date: Oct 2009
Posts: 184
Rep Power: 16
Quote:
 Originally Posted by nimasam Dear Elham Please re-read my points, you will figure out how it is derived. equation 18: 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: then right hand side of equation 18 becomes: so alpha1*div(U) is appearing here now consider boiling and condensation separately, for boiling the source term is Sb=S*alpha1 so for condensation Sb=S*(1.0-alpha1), then:
Sorry but I am still confused.

Base on

Quote:
 ddt(alpha1)+div(U,alpha1)+Sp*alpha1+Su=0
For boiling it should be :

Quote:
 Sp=- div(U) and Su= -S

 September 25, 2017, 12:50 #102 Senior Member   Nima Samkhaniani Join Date: Sep 2009 Location: Tehran, Iran Posts: 1,266 Blog Entries: 1 Rep Power: 24 i said Sb=S*alpha1, if you consider Su=S, then you missed alpha1 in equation. about negative signs, it is better to refer the code to see the exact formulation of MULES, i may missed the signs in post __________________ My Personal Website (http://nimasamkhaniani.ir/) Telegram channel (https://t.me/cfd_foam)

September 26, 2017, 10:11
#103
Senior Member

Elham
Join Date: Oct 2009
Posts: 184
Rep Power: 16
Quote:
 Originally Posted by nimasam i said Sb=S*alpha1, if you consider Su=S, then you missed alpha1 in equation. about negative signs, it is better to refer the code to see the exact formulation of MULES, i may missed the signs in post
Dear Nima,

Thanks for all kind answers. I got the idea.You have separated the general source term into two parts; condensation and boiling so that S=Sb+Sc. Then for the boiling phase you just have Sb as source term and for condensation just Sc.

Cheers,

Elham

November 4, 2017, 14:35
#104
New Member

Cláudio Corrêa
Join Date: Jun 2017
Location: Brazil
Posts: 14
Rep Power: 8
Quote:
 Originally Posted by nimasam 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
I'm interested in this problem too. I have trying to solve the problem of water vapor condensation in vertical pipe. My problem is three-dimensional. Could you help me solve this problem? I'm used OpenFoam 4.1

 November 5, 2017, 06:21 #105 Senior Member   Nima Samkhaniani Join Date: Sep 2009 Location: Tehran, Iran Posts: 1,266 Blog Entries: 1 Rep Power: 24 Dear Cláudio Corrêaread this post, i upload a solver based on OpenFOAM_2.2.0 which considers both evaporation and condensation and it is suitable for your work __________________ My Personal Website (http://nimasamkhaniani.ir/) Telegram channel (https://t.me/cfd_foam)

 December 6, 2017, 00:42 constructor classes #106 Senior Member   Elham Join Date: Oct 2009 Posts: 184 Rep Power: 16 Dear Nima, I am wondering why you defined two different word classes, word& type and word& alpha1Name in phaseChangeHeatFoam? In constructor of phaseChangeTwoPhaseMixture.C you have: Foam:haseChangeTwoPhaseMixture:haseChangeTwoPh aseMixture ( const word& type, const volVectorField& U, const surfaceScalarField& phi, const word& alpha1Name ) I would appreciate if you let me know. Cheers, Elham

 February 14, 2018, 03:28 anonymous term in energy equation #107 Senior Member   Elham Join Date: Oct 2009 Posts: 184 Rep Power: 16 Dear Niam, Regarding to phaseChangeHeatFoam the energy equation (25) in your paper is : ddT + div (UT) - lapplacian (k/rhoC, T) = -1/rhoC *mdotTriplePrime*Hlg there is a term in TEqn.H in line 19 Code: `- fvm::Sp(fvc::div(phi),T)` which is not corresponding to any term in energy equation. The whole left side of TEqn is as following: Code: ``` fvScalarMatrix TEqn ( fvm::ddt(T) + fvm::div(phi, T) - fvm::Sp(fvc::div(phi),T) - Foam::fvm::laplacian( k/rhoC_ , T,"laplacian(alphaEff,T)") /* + ( fvc::div(fvc::absolute(phi, U), p) + fvc::ddt(rho, KE) + fvc::div(rhoPhi, KE) )/rhoC */ );``` Would you please explain a bout it? Regards, Elham

 February 14, 2018, 05:27 #108 New Member   Saicharan Join Date: Jan 2018 Location: Bangalore, India Posts: 29 Rep Power: 8 Hi Nima Sam. Could you kindly send me your solver to the email ID saicharanb56@gmail.com ? I would really appreciate it. I am trying to simulate sessile droplet evaporation and I believe your solver will be ideal for me. I am, however, new to OpenFOAM and do not possess enough skill to code my own solver. Thanks in advance.

 February 20, 2018, 01:10 #109 Senior Member   Nima Samkhaniani Join Date: Sep 2009 Location: Tehran, Iran Posts: 1,266 Blog Entries: 1 Rep Power: 24 Dear Saicharan please look this post __________________ My Personal Website (http://nimasamkhaniani.ir/) Telegram channel (https://t.me/cfd_foam)

February 25, 2018, 12:22
#110
Senior Member

Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,266
Blog Entries: 1
Rep Power: 24
Quote:
 Originally Posted by Elham Dear Niam, Regarding to phaseChangeHeatFoam the energy equation (25) in your paper is : ddT + div (UT) - lapplacian (k/rhoC, T) = -1/rhoC *mdotTriplePrime*Hlg there is a term in TEqn.H in line 19 Code: `- fvm::Sp(fvc::div(phi),T)` which is not corresponding to any term in energy equation. The whole left side of TEqn is as following: Code: ``` fvScalarMatrix TEqn ( fvm::ddt(T) + fvm::div(phi, T) - fvm::Sp(fvc::div(phi),T) - Foam::fvm::laplacian( k/rhoC_ , T,"laplacian(alphaEff,T)") /* + ( fvc::div(fvc::absolute(phi, U), p) + fvc::ddt(rho, KE) + fvc::div(rhoPhi, KE) )/rhoC */ );``` Would you please explain a bout it? Regards, Elham
fvc::div(phi)=0, so theoretically it is zero, but numerically it is not, and it help numerical stability
__________________
My Personal Website (http://nimasamkhaniani.ir/)
Telegram channel (https://t.me/cfd_foam)

 February 25, 2018, 20:32 #111 Senior Member   Elham Join Date: Oct 2009 Posts: 184 Rep Power: 16 Thanks. Elham

July 20, 2018, 10:12
#112
Senior Member

Przemek
Join Date: Jun 2011
Posts: 249
Rep Power: 15
Quote:
 Originally Posted by Elham Dear Niam, Regarding to phaseChangeHeatFoam the energy equation (25) in your paper is : ddT + div (UT) - lapplacian (k/rhoC, T) = -1/rhoC *mdotTriplePrime*Hlg there is a term in TEqn.H in line 19 Code: `- fvm::Sp(fvc::div(phi),T)` which is not corresponding to any term in energy equation. The whole left side of TEqn is as following: Code: ``` fvScalarMatrix TEqn ( fvm::ddt(T) + fvm::div(phi, T) - fvm::Sp(fvc::div(phi),T) - Foam::fvm::laplacian( k/rhoC_ , T,"laplacian(alphaEff,T)") /* + ( fvc::div(fvc::absolute(phi, U), p) + fvc::ddt(rho, KE) + fvc::div(rhoPhi, KE) )/rhoC */ );``` Would you please explain a bout it? Regards, Elham

Dear Elham,
It is also written in the documentation
https://openfoam.org/release/2-2-0/n...s-boundedness/

According to this link the same can be achieved with
setting divergence scheme as bounded in fvScheme file, e.g.
Code:
` div(phi,T)      bounded Gauss upwind;`
__________________
best regards
pblasiak

 May 18, 2022, 08:14 #113 Member   hari charan Join Date: Sep 2021 Location: India,hyderabad Posts: 96 Rep Power: 4 Hi nimasam, Could you tell me how you solved your problem if not possible to send .If you can send you can mail me to this mail : saicharan662000@gmail.com Thanks and regards

February 24, 2023, 15:17
Ancient code
#114
New Member

Michael Jensen
Join Date: May 2022
Posts: 27
Rep Power: 3
Quote:
 Originally Posted by 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

To anybody like me finding this post all these years later, it uses openfoam 2.2.0, for which installation instructions exist only to up only to ubuntu 13.10, which is, in this case, unsupported and obsolete. If you follow those instructions, they will not succeed. Because the files they depend on won't be there. I have contacted the author about this and been directed to.. I can't recall where, (..I could track it down..) but it was one of the documentation pages you'll find if you search for evaporation/concensation/openfoam. So this is merely a confirmation of what you surely already expect. I can vouch that the author is very friendly if you are reasonable with whatever request you might make :-)

Best regards, Mike :-)

 Tags boiling, evaporation, interfoam, phase change