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

evapPhaseChangeFOAM

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

Like Tree14Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 31, 2016, 10:07
Default boiling in openFOAM with evapPhaseChangeFOAM
  #81
New Member
 
Shyam Sunder
Join Date: Sep 2015
Posts: 27
Rep Power: 9
ssyadav is on a distinguished road
Hello Nimasam

Kindly send me the OpenFoam based code to simulate boiling with OpenFoam. My email id. is ssyadav25@gmail.com

Thanks in advance.

Shyam Sunder
ssyadav is offline   Reply With Quote

Old   May 15, 2016, 06:44
Default VOF solver for condensation
  #82
Member
 
Elisabet Mas de les Valls
Join Date: Mar 2009
Location: Barcelona, Spain
Posts: 64
Rep Power: 15
elisabet is on a distinguished road
Nima and Alexander,

I'm a little bit confused about the state of the art of your solvers. I'ld really appreciate you could send me your codes and, if possible, comments. I'll try to understand the differences and limitations.

I'm currently interested in steam condensation, so any suggestion is welcome!
Thank you in advance,

Elisabet (elisabet.masdelesvalls@gits.ws)
elisabet is offline   Reply With Quote

Old   May 16, 2016, 16:29
Default
  #83
Member
 
Alexander
Join Date: Mar 2009
Posts: 49
Rep Power: 15
sahas is on a distinguished road
Hello, Elisabet!

Currently I have no deal with simulation of evaporation/condensation using VOF method. Some information about evaporation model implementation can be found in this thread. Still remaining problem is taking correctly into account the latent heat. As mentioned in the literature, the possible solution is to introduce weight function in such a way that the integral from the volume latent heat source multiplied by this weight function gives the value of correct latent heat power. Unfortunately, I do not have much information about this.
sahas is offline   Reply With Quote

Old   May 30, 2016, 22:51
Default
  #84
New Member
 
LONG JIAO
Join Date: May 2016
Posts: 1
Rep Power: 0
serena_bears is on a distinguished road
Hi Nimasan,

could you please send me your final code? Now I'm working on droplet condensation problem. Thank you very much.

My email is LONG0032@e.ntu.edu.sg
serena_bears is offline   Reply With Quote

Old   July 15, 2016, 10:21
Default
  #85
New Member
 
Shen shiquan
Join Date: Jul 2016
Location: The State Key Laboratory of Engines (Tianjin University)
Posts: 12
Rep Power: 8
tju_shq is on a distinguished road
Hi Nimasan,

could you please send me your final code? Now I'm working on droplet condensation problem. Thank you very much.

My email is shiquan@tju.edu.cn
tju_shq is offline   Reply With Quote

Old   September 9, 2016, 16:48
Exclamation Link to solver files not working
  #86
New Member
 
Aslamah Rahman
Join Date: Sep 2016
Posts: 1
Rep Power: 0
asla9796 is on a distinguished road
Hello, I am new to CFD and is trying to simulate the same problem of droplet evaporation. This link posted in the thread which contains the files for solver appears to be empty now. Please do look into this and do share the code again if possible.

My mail ID is mm14b007@smail.iitm.ac.in

Thank you for your time

Last edited by asla9796; September 9, 2016 at 16:50. Reason: forgot to share id
asla9796 is offline   Reply With Quote

Old   September 11, 2016, 11:50
Default
  #87
Senior Member
 
Przemek
Join Date: Jun 2011
Posts: 245
Rep Power: 14
gaza is on a distinguished road
Hi
for now my code does work
I will share my code when it will work
see to OpenFOAM 1606 there is interEvaporatingCondensingFoam
__________________
best regards
pblasiak
gaza is offline   Reply With Quote

Old   September 20, 2016, 12:20
Default
  #88
New Member
 
AndyKing
Join Date: Dec 2013
Location: China
Posts: 1
Rep Power: 0
Ank09 is on a distinguished road
Hi Nimasam!
I'm working on the droplet evaporation topic and I'm very interested in your code. Could you please send me a copy of your final code?
Thanks a lot! My best wishes!
My email: andyking0928@foxmail.com
Ank09 is offline   Reply With Quote

Old   April 1, 2017, 13:16
Default
  #89
New Member
 
Gautham Krishnan
Join Date: Feb 2016
Posts: 4
Rep Power: 8
GauthamKrishnan is on a distinguished road
Hi Nimasam!
I'm working on the droplet evaporation and I'm very interested in your code. Could you please send me a copy of your final code+ case file
Thanks a lot! My best wishes!
My email: [email]gkriz66@gmail.com
GauthamKrishnan is offline   Reply With Quote

Old   June 7, 2017, 05:41
Default phaseChangeHeatFoam
  #90
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,265
Blog Entries: 1
Rep Power: 23
nimasam is on a distinguished road
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
__________________
Telegram channel (https://telegram.me/openfoam4Iranian)
My Weblog in Persian(http://openfoam.blogfa.com/)
My Personal Website (http://nimasamkhaniani.ir/)
nimasam is offline   Reply With Quote

Old   June 7, 2017, 13:41
Default phaseChangeHeatFoam compilation error
  #91
New Member
 
Shyam Sunder
Join Date: Sep 2015
Posts: 27
Rep Power: 9
ssyadav is on a distinguished road
Dear Nima

Thank you very much for posting your solver. I tried compiling it in OpenFOAM 2.2.2, I am getting some errors which are as follows:

1) smoothInterfaceProperties.C:27:49: fatal error: alphaContactAngleFvPatchScalarField.H: No such file or directory
#include "alphaContactAngleFvPatchScalarField.H"

I could resolve this error by changing the path in the "smoothInterfaceProperties/Make/options" file

2) /usr/bin/ld: cannot find -ltwoPhaseInterfaceProperties

This library is included in "phaseChangeHeatFoam/Make/options" but I could not locate it my OpenFOAM installation.

Please help.

Also how I can extend your solver to include Conjugate heat transfer in a solid, particularly the boundary condition at the solid-fluid interface, i.e. the incompressible version of "compressible::turbulentTemperatureCoupledBaffleMi xed"

Thanks.

Shyam
ssyadav is offline   Reply With Quote

Old   June 8, 2017, 02:46
Default
  #92
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,265
Blog Entries: 1
Rep Power: 23
nimasam is on a distinguished road
Dear Shyam

1- for your second question, i guess the name of library has been changed in newer version of OpenFOAM, i dont have OpenFOAM-2.2.0 to check it my self. but i guess it became libtwoPhaseProperties or libinterfaceProperties

2- you can see this paper
__________________
Telegram channel (https://telegram.me/openfoam4Iranian)
My Weblog in Persian(http://openfoam.blogfa.com/)
My Personal Website (http://nimasamkhaniani.ir/)
nimasam is offline   Reply With Quote

Old   August 31, 2017, 13:13
Default
  #93
Senior Member
 
Elham
Join Date: Oct 2009
Posts: 184
Rep Power: 15
Elham is on a distinguished road
Quote:
Originally Posted by nimasam View Post
Dear Shyam

1- for your second question, i guess the name of library has been changed in newer version of OpenFOAM, i dont have OpenFOAM-2.2.0 to check it my self. but i guess it became libtwoPhaseProperties or libinterfaceProperties

2- you can see this paper
Dear Nima,

I cannot compile your code due to the following error:

Code:
In file included from phaseChangeHeatFoam.C:57:0:
createFields.H: In function ‘int main(int, char**)’:
createFields.H:106:72: error: no matching function for call to ‘Foam::smoothInterfaceProperties::smoothInterfaceProperties(Foam::volScalarField&, Foam::volVectorField&, Foam::phaseChangeTwoPhaseMixture&)’
Do you have any idea to fix it?

Cheers,

Elham
Elham is offline   Reply With Quote

Old   August 31, 2017, 14:44
Default
  #94
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,265
Blog Entries: 1
Rep Power: 23
nimasam is on a distinguished road
1-which version of openfoam, do u use?
2-you should first compile smoothInterfaceProperties using following command:
Quote:
wmake libso
then you can compile the solver
__________________
Telegram channel (https://telegram.me/openfoam4Iranian)
My Weblog in Persian(http://openfoam.blogfa.com/)
My Personal Website (http://nimasamkhaniani.ir/)
nimasam is offline   Reply With Quote

Old   September 1, 2017, 09:53
Default
  #95
Senior Member
 
Elham
Join Date: Oct 2009
Posts: 184
Rep Power: 15
Elham is on a distinguished road
Quote:
Originally Posted by nimasam View Post
1-which version of openfoam, do u use?
2-you should first compile smoothInterfaceProperties using following command:

then you can compile the solver
The library is compiled but still have the error. I use OpenFOAM 2.3.0.
Elham is offline   Reply With Quote

Old   September 2, 2017, 02:19
Default
  #96
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,265
Blog Entries: 1
Rep Power: 23
nimasam is on a distinguished road
Dear Elham it's been written for OpenFOAM_220, for upper versions of openfoam, you may need to change some of libraries name, because the basic solver which is interfoam is changed a lot through openfoam 2.2 to openfoam 3.0
__________________
Telegram channel (https://telegram.me/openfoam4Iranian)
My Weblog in Persian(http://openfoam.blogfa.com/)
My Personal Website (http://nimasamkhaniani.ir/)
nimasam is offline   Reply With Quote

Old   September 22, 2017, 13:06
Default
  #97
Senior Member
 
Elham
Join Date: Oct 2009
Posts: 184
Rep Power: 15
Elham is on a distinguished road
Dear Nima,

I have read your paper "The evaluation of the diffuse interface method for phase change simulations using OpenFOAM" and the phaseChangeHeatFoam code. To derive Sp you have used the right hand side of equation 18 but I am still wondering how you have derived su , as part of source term?
I would appreciate if you give me some clues.

Cheers,

Elham
Elham is offline   Reply With Quote

Old   September 23, 2017, 05:01
Default
  #98
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,265
Blog Entries: 1
Rep Power: 23
nimasam is on a distinguished road
Several points to clarify coding:
1)as i remember the formulation of MULES is like that:
Quote:
ddt(alpha1)+div(U,alpha1)+Sp*alpha1+Su=0
2- convection term in equation 18, then:
Quote:
U & grad(alpha1)=div(U,alpha1)-alpha.div(U)
3-the mass transfer source term is multiplying with void fraction, means:
Quote:
mDotBoiling=mdot*alpha1;
mDotCondensation=mdot*(1-alpha1);
__________________
Telegram channel (https://telegram.me/openfoam4Iranian)
My Weblog in Persian(http://openfoam.blogfa.com/)
My Personal Website (http://nimasamkhaniani.ir/)
nimasam is offline   Reply With Quote

Old   September 25, 2017, 02:10
Default
  #99
Senior Member
 
Elham
Join Date: Oct 2009
Posts: 184
Rep Power: 15
Elham is on a distinguished road
Quote:
Originally Posted by nimasam View Post
Several points to clarify coding:
1)as i remember the formulation of MULES is like that:

2- convection term in equation 18, then:

3-the mass transfer source term is multiplying with void fraction, means:
Dear Nima,

I am confused about deriving the source terms. Equation 18 is as following:

ddt(alpha1)+U& grad(alpha1) +grad(alpha1(1-alpha1)Uc) = -mDot(1/rho1-alpha1(1/rho1-1/rho2)))

and in eq 16:

div (U) = mDot(1/rho2-1/rho1)

substituting eq 16 into eq 18 and rearranging:

ddt(alpha1) + grad(U,alpha1) + grad (alpha1(1-alpha1)Uc) = -mDot/rho1

So we should have:

sp = 0 and su = mDot/rho1

But in alphaEqn.C:

sp=mDot(1/rho1-alpha1(1/rho-1/rho2))
su=divU*alpha1 + mDot(1/rho1-alpha1(1/rho-1/rho2))

Thnaks for any help.

Elham
Elham is offline   Reply With Quote

Old   September 25, 2017, 12:21
Default
  #100
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,265
Blog Entries: 1
Rep Power: 23
nimasam is on a distinguished road
Dear Elham
Please re-read my points, you will figure out how it is derived.
equation 18:
Quote:
dt(alpha1)+U& grad(alpha1) +grad(alpha1(1-alpha1)Uc) = -mDot(1/rho1-alpha1(1/rho1-1/rho2)))
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:
Quote:
U & grad(alpha1)=div(U,alpha1)-alpha1*div(U)
then right hand side of equation 18 becomes:
Quote:
dt(alpha1)+div(U,alpha1)-alpha.div(U) +grad(alpha1(1-alpha1)Uc)
so alpha1*div(U) is appearing here
now consider boiling and condensation separately,
for boiling the source term is Sb=S*alpha1
so
Quote:
su=alpha1*div(U)
sp=S
for condensation Sb=S*(1.0-alpha1), then:
Quote:
su=alpha1*div(U)+S
sp=-S
__________________
Telegram channel (https://telegram.me/openfoam4Iranian)
My Weblog in Persian(http://openfoam.blogfa.com/)
My Personal Website (http://nimasamkhaniani.ir/)
nimasam is offline   Reply With Quote

Reply

Tags
boiling, evaporation, interfoam, phase change

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 08:57.