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 26, 2015, 17:51
Default
  #61
New Member
 
Mostafa Mobli
Join Date: Feb 2015
Posts: 8
Rep Power: 9
moblimostafa is on a distinguished road
Dear Nima

Could you send me your your code?
mostafamobli@gmail.com
moblimostafa is offline   Reply With Quote

Old   April 23, 2015, 14:22
Default
  #62
New Member
 
Yaoyina
Join Date: Sep 2014
Posts: 7
Rep Power: 10
yaoyina is on a distinguished road
Hello, Nima Sam.

Would you be so kind to send your code to yao-yn14@mails.tsinghua.edu.cn?
Thanks in advance!
yaoyina is offline   Reply With Quote

Old   May 26, 2015, 16:09
Default
  #63
New Member
 
Chad
Join Date: Sep 2011
Posts: 16
Rep Power: 13
chaolian is on a distinguished road
Hi, Nima

Thanks for your contribution. I am interested in your solver.
Could you please send your solver to me: chaolian@mtu.edu

Thanks in advance.

chaolian is offline   Reply With Quote

Old   May 27, 2015, 17:11
Default
  #64
Senior Member
 
Join Date: May 2011
Posts: 231
Rep Power: 14
Kanarya is on a distinguished road
Hi Nima,

I am also interested in your code
kah.recep@gmail.com

thanks in advance!!
Kanarya is offline   Reply With Quote

Old   July 9, 2015, 00:50
Default
  #65
New Member
 
Join Date: Dec 2012
Posts: 19
Rep Power: 12
smraniaki is on a distinguished road
Hello Nimasam

Thank you very much for your contribution to the world of OpenFoam. I am also interested to see your code. Would you please email it to rifel2002@yahoo.com
Smran
smraniaki is offline   Reply With Quote

Old   August 7, 2015, 10:36
Default
  #66
Senior Member
 
Przemek
Join Date: Jun 2011
Posts: 245
Rep Power: 14
gaza is on a distinguished road
Hi Nima Sam,
I compiled your solver but the results are wrong.
I used the fluids properties the same as Bitan Shu. You can see here

http://www.cfd-online.com/Forums/ope...pvofhardt.html

I am also trying to do solver for evaporation/condensation but my results are also wrong.
Can anybody upload correct solver for evaporation/condensation purposes?
__________________
best regards
pblasiak
gaza is offline   Reply With Quote

Old   August 7, 2015, 10:46
Default
  #67
Senior Member
 
Join Date: May 2011
Posts: 231
Rep Power: 14
Kanarya is on a distinguished road
I do not think, there is correct one!
what is wrong with your results?

Quote:
Originally Posted by gaza View Post
Hi Nima Sam,
I compiled your solver but the results are wrong.
I used the fluids properties the same as Bitan Shu. You can see here

http://www.cfd-online.com/Forums/ope...pvofhardt.html

I am also trying to do solver for evaporation/condensation but my results are also wrong.
Can anybody upload correct solver for evaporation/condensation purposes?
gaza likes this.
Kanarya is offline   Reply With Quote

Old   August 11, 2015, 04:43
Default
  #68
Senior Member
 
Przemek
Join Date: Jun 2011
Posts: 245
Rep Power: 14
gaza is on a distinguished road
Hi Kanarya,
In attachement

http://fluid.itcmp.pwr.wroc.pl/~pblasiak/download.html

there is stefan problem case with settings as in Bitan Shu conference paper:

http://pl.scribd.com/doc/28869783/Pa...0080214#scribd

Additionally there is file interfacePosStefProbl_analSol.dat where are analitycal results for this case for time 0-25 s. For 25 s interface position should be 3.8 mm but numerical value is much smaller (see figure in the attachement).

I created my own solver. Results are close to analytical but still incorrect.
I uploaded this with damBreak case here

http://www.cfd-online.com/Forums/ope...-modeling.html

In Nima Sam I also added additional term as in my solver in TEqn.H to avoid unboundedness in temperature field.

What do you think about Nima Sam solver and my?
Do you know what is wrong?
Did you validate Nima Sam solver on damBreak case?
__________________
best regards
pblasiak
gaza is offline   Reply With Quote

Old   August 11, 2015, 05:50
Default
  #69
Senior Member
 
Join Date: May 2011
Posts: 231
Rep Power: 14
Kanarya is on a distinguished road
Hi Gaza,

Which solver is yours?
Do you consider that your water or vapour is reached in saturation ?
If not, you need source and sink term for both phases...
For instance, if water is cool enough, it needs to gain energy in order to evaporate for this reason you need source term...

Best regards,
Kanarya is offline   Reply With Quote

Old   August 11, 2015, 07:54
Default
  #70
Senior Member
 
Przemek
Join Date: Jun 2011
Posts: 245
Rep Power: 14
gaza is on a distinguished road
Hi Kanarya,
My solver is interEvapCondPhaseChangeFoam (first link). The second link is a test case for solver
evaPhaseChangeFoamU (solver created by Nima Sam). As you can see in the figure in attachement
numeirical results is not equal analytical, which should be 3.8mm for 25 s for fluid properties used (see Bitan Shu paper). Why do you think that Nima Sam solver works well?

I did not consider that water is in saturation condition because I do not know how to provide this or
temperature at interface equal to saturation temperature. I only used initial conditions for water close to saturation temperature.

Here you can see how to compile my solver

http://www.cfd-online.com/Forums/ope...-modeling.html

If you run stefan problem with my solver you can see that results are close to analytical but still not well. It is clear in damBreak case where you can see the difference in comparison to interFoam.
Do you know how can I improve my solver?
__________________
best regards
pblasiak
gaza is offline   Reply With Quote

Old   August 11, 2015, 08:04
Default
  #71
Senior Member
 
Przemek
Join Date: Jun 2011
Posts: 245
Rep Power: 14
gaza is on a distinguished road
Quote:
Originally Posted by Kanarya View Post
Hi Gaza,

Which solver is yours?
Do you consider that your water or vapour is reached in saturation ?
If not, you need source and sink term for both phases...
For instance, if water is cool enough, it needs to gain energy in order to evaporate for this reason you need source term...

Best regards,

And one more question: what do you mean "you need source and sink term for both phases"?
My TEqn looks like

fvScalarMatrix TEqn
(
fvm::ddt(rhoCp, T)
+ fvm::div(rhoPhiCpf, T)
- fvm::Sp(fvc::ddt(rhoCp) + fvc::div(rhoPhiCpf), T) // additional term to avoid uboundedness ot T
- fvm::laplacian(kappaf, T)
==
fvm::Sp(hEvap*(mCond + mVap)/T, T)
);

hEvap - enthalpy of vaporisation
mCond - mass due to condensation
mVap - mass due to evaporation

mCond and mVap are calculated using phase change model. In this case Ganapathy.
__________________
best regards
pblasiak
gaza is offline   Reply With Quote

Old   August 11, 2015, 08:07
Default
  #72
Senior Member
 
Join Date: May 2011
Posts: 231
Rep Power: 14
Kanarya is on a distinguished road
As I can see here you did not specified for which phase you are solving T equation right?so there is not volume fraction of alpha or sth?or I am wrong?

Quote:
Originally Posted by gaza View Post
And one more question: what do you mean "you need source and sink term for both phases"?
My TEqn looks like

fvScalarMatrix TEqn
(
fvm::ddt(rhoCp, T)
+ fvm::div(rhoPhiCpf, T)
- fvm::Sp(fvc::ddt(rhoCp) + fvc::div(rhoPhiCpf), T) // additional term to avoid uboundedness ot T
- fvm::laplacian(kappaf, T)
==
fvm::Sp(hEvap*(mCond + mVap)/T, T)
);

hEvap - enthalpy of vaporisation
mCond - mass due to condensation
mVap - mass due to evaporation

mCond and mVap are calculated using phase change model. In this case Ganapathy.
Kanarya is offline   Reply With Quote

Old   August 11, 2015, 08:12
Default
  #73
Senior Member
 
Przemek
Join Date: Jun 2011
Posts: 245
Rep Power: 14
gaza is on a distinguished road
Quote:
Originally Posted by Kanarya View Post
As I can see here you did not specified for which phase you are solving T equation right?so there is not volume fraction of alpha or sth?or I am wrong?
I created my solver based on interPhaseChangeFoam. In VOF energy equation is solved for both phases.
__________________
best regards
pblasiak
gaza is offline   Reply With Quote

Old   August 11, 2015, 08:18
Default
  #74
Senior Member
 
Join Date: May 2011
Posts: 231
Rep Power: 14
Kanarya is on a distinguished road
yes but if you look to the solver alpha.H it is only solving for one phase not like two fluid model, if you check the twoPhaseEulerFoam you will see...
what is you final gaol...what do you want to simulate with this solver ...probably two fluid model suits you better...
Quote:
Originally Posted by gaza View Post
I created my solver based on interPhaseChangeFoam. In VOF energy equation is solved for both phases.
Kanarya is offline   Reply With Quote

Old   August 11, 2015, 08:30
Default
  #75
Senior Member
 
Przemek
Join Date: Jun 2011
Posts: 245
Rep Power: 14
gaza is on a distinguished road
Quote:
Originally Posted by Kanarya View Post
yes but if you look to the solver alpha.H it is only solving for one phase not like two fluid model, if you check the twoPhaseEulerFoam you will see...
what is you final gaol...what do you want to simulate with this solver ...probably two fluid model suits you better...
yes but I want to use one fluid approach because it is less computationally expensive.
I think that alpha equation in my solver is solved in wrong way but I do not know how to fix it.

Are you using Nima Sam solver? What do you think about results of it?
__________________
best regards
pblasiak
gaza is offline   Reply With Quote

Old   August 11, 2015, 08:31
Default
  #76
Senior Member
 
Join Date: May 2011
Posts: 231
Rep Power: 14
Kanarya is on a distinguished road
no I am not...I do not have it..could you please send it to me ?
Quote:
Originally Posted by gaza View Post
yes but I want to use one fluid approach because it is less computationally expensive.
I think that alpha equation in my solver is solved in wrong way but I do not know how to fix it.

Are you using Nima Sam solver? What do you think about results of it?
Kanarya is offline   Reply With Quote

Old   August 11, 2015, 08:44
Default
  #77
Senior Member
 
Przemek
Join Date: Jun 2011
Posts: 245
Rep Power: 14
gaza is on a distinguished road
Quote:
Originally Posted by Kanarya View Post
no I am not...I do not have it..could you please send it to me ?
yes of course, here it is (third link):

http://fluid.itcmp.pwr.wroc.pl/~pblasiak/download.html

it is for OF1.6.x
__________________
best regards
pblasiak
gaza is offline   Reply With Quote

Old   September 21, 2015, 02:24
Default Hi Nimasam
  #78
Member
 
Akr
Join Date: Apr 2015
Location: India
Posts: 53
Rep Power: 9
NightWing is on a distinguished road
Can i get the final code???


email: arunkraj03@gmail.com
NightWing is offline   Reply With Quote

Old   December 23, 2015, 07:54
Default
  #79
Member
 
thomas
Join Date: Jul 2014
Posts: 50
Rep Power: 10
thomas. is on a distinguished road
Hi Nimasan,

could you please send me your code?

My email is thomas.schaub.h@gmail.com

Thank you very much in advance!
thomas. is offline   Reply With Quote

Old   December 29, 2015, 05:14
Smile
  #80
LCC
New Member
 
Li Cong
Join Date: Dec 2015
Posts: 1
Rep Power: 0
LCC is on a distinguished road
Hello, nimasam. I have some problems about two phase change progress. Can you send me your final code? My email is 18600156862@163.com. Thank you very much!
LCC 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 04:01.