
[Sponsors] 
March 26, 2015, 17:51 

#61 
New Member
Mostafa Mobli
Join Date: Feb 2015
Posts: 8
Rep Power: 9 

April 23, 2015, 14:22 

#62 
New Member
Yaoyina
Join Date: Sep 2014
Posts: 7
Rep Power: 10 
Hello, Nima Sam.
Would you be so kind to send your code to yaoyn14@mails.tsinghua.edu.cn? Thanks in advance! 

May 26, 2015, 16:09 

#63 
New Member
Chad
Join Date: Sep 2011
Posts: 16
Rep Power: 13 
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. 

May 27, 2015, 17:11 

#64 
Senior Member
Join Date: May 2011
Posts: 231
Rep Power: 14 

July 9, 2015, 00:50 

#65 
New Member
Join Date: Dec 2012
Posts: 19
Rep Power: 12 
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 

August 7, 2015, 10:36 

#66 
Senior Member
Przemek
Join Date: Jun 2011
Posts: 245
Rep Power: 14 
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.cfdonline.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 

August 7, 2015, 10:46 

#67  
Senior Member
Join Date: May 2011
Posts: 231
Rep Power: 14 
I do not think, there is correct one!
what is wrong with your results? Quote:


August 11, 2015, 04:43 

#68 
Senior Member
Przemek
Join Date: Jun 2011
Posts: 245
Rep Power: 14 
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 025 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.cfdonline.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 

August 11, 2015, 05:50 

#69 
Senior Member
Join Date: May 2011
Posts: 231
Rep Power: 14 
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, 

August 11, 2015, 07:54 

#70 
Senior Member
Przemek
Join Date: Jun 2011
Posts: 245
Rep Power: 14 
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.cfdonline.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 

August 11, 2015, 08:04 

#71  
Senior Member
Przemek
Join Date: Jun 2011
Posts: 245
Rep Power: 14 
Quote:
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 

August 11, 2015, 08:07 

#72  
Senior Member
Join Date: May 2011
Posts: 231
Rep Power: 14 
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:


August 11, 2015, 08:12 

#73 
Senior Member
Przemek
Join Date: Jun 2011
Posts: 245
Rep Power: 14 
I created my solver based on interPhaseChangeFoam. In VOF energy equation is solved for both phases.
__________________
best regards pblasiak 

August 11, 2015, 08:18 

#74 
Senior Member
Join Date: May 2011
Posts: 231
Rep Power: 14 
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... 

August 11, 2015, 08:30 

#75  
Senior Member
Przemek
Join Date: Jun 2011
Posts: 245
Rep Power: 14 
Quote:
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 

August 11, 2015, 08:31 

#76 
Senior Member
Join Date: May 2011
Posts: 231
Rep Power: 14 
no I am not...I do not have it..could you please send it to me ?


August 11, 2015, 08:44 

#77  
Senior Member
Przemek
Join Date: Jun 2011
Posts: 245
Rep Power: 14 
Quote:
http://fluid.itcmp.pwr.wroc.pl/~pblasiak/download.html it is for OF1.6.x
__________________
best regards pblasiak 

September 21, 2015, 02:24 
Hi Nimasam

#78 
Member
Akr
Join Date: Apr 2015
Location: India
Posts: 53
Rep Power: 9 

December 23, 2015, 07:54 

#79 
Member
thomas
Join Date: Jul 2014
Posts: 50
Rep Power: 10 
Hi Nimasan,
could you please send me your code? My email is thomas.schaub.h@gmail.com Thank you very much in advance! 

December 29, 2015, 05:14 

#80 
New Member
Li Cong
Join Date: Dec 2015
Posts: 1
Rep Power: 0 
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!


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

