
[Sponsors] 
Multiphase flow and Phase change due to heat transferevaporation 

LinkBack  Thread Tools  Display Modes 
June 2, 2009, 13:58 

#41 
Member
Rachel Vogl
Join Date: Jun 2009
Posts: 48
Rep Power: 9 
Dear All,
I would like to know if this solver is available for testing/preview/further development? I am interested in a similar problem involving phase change and chemical reactions. 

June 9, 2009, 04:55 

#42 
Senior Member
Sandy Lee
Join Date: Mar 2009
Posts: 213
Rep Power: 10 
Dear all,
I have three questions about interPhaseChangeFoam solver: 1) In OF, it seems that the model fomulation in the file Kunz.C is different from the literature, for example, in OF, it is: Kunz::mDotAlphal(): ( mcCoeff_*sqr(limitedAlpha1) *max(p  pSat(), p0_)/max(p  pSat(), 0.01*pSat()), mvCoeff_*min(p  pSat(), p0_) ) Kunz::mDotP(): ( mcCoeff_*sqr(limitedAlpha1)*(1.0  limitedAlpha1) *pos(p  pSat())/max(p  pSat(), 0.01*pSat()), (mvCoeff_)*limitedAlpha1*neg(p  pSat()) ) where, limitedAlpha1 = min(max(alpha1_,scalar(0)),scalar(1)); mcCoeff_ = Cc_*rho2()/tInf_; mvCoeff_ = Cv_*rho2()/(0.5*rho1()*sqr(UInf_)*tInf_); However, in literature, it is: Mprod=(Cv/0.5*UInf2*tInf2)rho2/rho1*gamma*min[0,ppv] Mdest=(Cc/tInf)rho2*gamma2[1gamma] Are they same? What is the differences? 2) In the UEqn.H of interPhaseChangeFoam,  UEqn.H  surfaceScalarField muf = twoPhaseProperties>muf() + fvc::interpolate(rho*turbulence>nuSgs()); fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U)  fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)  fvm::laplacian(muf, U)  (fvc::grad(U) & fvc::grad(muf)) // fvc::div(muf*(fvc::interpolate(dev2(fvc::grad(U))) & mesh.Sf())) ); UEqn.relax(); if (momentumPredictor) { solve ( UEqn == fvc::reconstruct ( ( fvc::interpolate(interface.sigmaK())*fvc::snGrad(g amma)  ghf*fvc::snGrad(rho)  fvc::snGrad(pd) ) * mesh.magSf() ) ); }  What is the meanings ' UEqn.relax() ' ? In PISO, it still need some relaxation factors? wrong ... 3) If I want to simulate the full wet flows of a hydrofoil firstly, how should I specify the parameters? I just need to set Cc_ and Cv_ 0, right? But I still could not get the convergent results, why? Could you help me to explain them? Thank you very much. Regards, Sandy sandy.lee37@gmail.com 

June 21, 2009, 08:27 

#43 
Senior Member
Sandy Lee
Join Date: Mar 2009
Posts: 213
Rep Power: 10 
Hi man,
Now, I don't want to fight with the Liquid Hydrogen or flash evaporation again. I just want to make the foil caviting, but why my gamma equation can not be solved by interPhaseChangeFoam. I choose the oneEqEddy in the Les turbulent model and Kunz's interPhaseChange model during my simulating. I could not find what is the matter with them. Please give a finger. Thank you very much. Sandy 

June 30, 2009, 11:59 

#44 
Senior Member
isabel
Join Date: Apr 2009
Location: Spain
Posts: 171
Rep Power: 9 
In the interPhaseChange solver. In gammaEqnSubCycle.H, we have this line:
MULES::implicitSolve(oneField(),gamma,phi,phiGamma ,Sp,Su,1,0) Does anybody know what means " oneField() " ? And Sp, Su ? 

June 30, 2009, 12:38 

#45 
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,953
Rep Power: 41 
The name means "Questions about this class should only be asked once (and not in multiple unrelated threads) otherwise people will stop answering"
Another explanation can be found here: http://foam.sourceforge.net/doc/Doxy...1oneField.html The second explanation is the right one, but I like the first one better 

July 1, 2009, 08:28 

#46 
Senior Member
Sandy Lee
Join Date: Mar 2009
Posts: 213
Rep Power: 10 
You are great! Gschaider. Sure, we should study from busy ants from onefield to another and others in everyday, right? And should have multihead but one .
In fact, I do feel confused to use this solver in these days . In Kunz.C and phaseChangeTwoPhaseMixture.C of this solver, the variables mDotAlphal() and vDotAlphal() are const, so if I want to change the values Cc and Cv of the file transportProperties from 0 (namely, the fullwet flows) to 30000 and 900000 (namely, Kunz's model), the code will report error infomations. I would like to know how to avoid this problems? Thanks. 

July 2, 2009, 22:54 

#47 
Senior Member
Sandy Lee
Join Date: Mar 2009
Posts: 213
Rep Power: 10 
Hi isabel, Sp and Su are the direct variables of vDotAlphal(), however, they are just indirectly connect with mDotAlphal().I think it is maybe a key to this solver.


July 7, 2009, 05:06 

#48 
Senior Member
Sandy Lee
Join Date: Mar 2009
Posts: 213
Rep Power: 10 
Sp = (vDotvAlphal  vDotcAlphal) * gamma , and it will be solved implicitly.
Su = vDotcAlphal, and it is an explicit term in this equation. But, why Su also includes the term divU*gamma in gammaEqu.H ? If the MULES::implicitSolver was chose, this term should be deleted, right? Who knew it? Please help me out. Thanks. 

July 16, 2009, 05:54 

#49 
Senior Member
isabel
Join Date: Apr 2009
Location: Spain
Posts: 171
Rep Power: 9 
I have studied the tutorial "Solve cavitating flow arroun a 2d hydrofoil using a user modified version of interPhaseChangeFoam", cut I don't understand what these mean:
vDotvAlphal vDotcalphal vDotvP vDotcP Does anybody know what these mean? 

July 16, 2009, 09:22 

#50  
Senior Member
Sandy Lee
Join Date: Mar 2009
Posts: 213
Rep Power: 10 
Quote:
vDotcalphal = =(Cc*rho2/tInf)*gamma^2*{1.0/rho1  gamma*(1.0/rho1  1.0/rho2)} vDotvP = {Cv*rho2/(0.5*rho1*sqr(UInf)*tInf)}*gamma*(1.0/rho1  1.0/rho2) vDotcP = ={(Cc*rho2/tInf)*gamma^2*(1.0  gamma)/(ppSat)} *(1.0/rho1  1.0/rho2) Am I right? Last edited by sandy; July 20, 2009 at 00:49. 

July 17, 2009, 02:56 

#51 
Senior Member
isabel
Join Date: Apr 2009
Location: Spain
Posts: 171
Rep Power: 9 
Thank you very much, sandy.
By the way, what you asked before about what is the meaning ' UEqn.relax() ' In PISO. This line is to adjust the coefficients so as to incorporate the underrelaxation coefficient, such that a solution to UEqn == source will produce a partlyrelaxed version of U. It is explained in this link: Trying to figure out the details of simpleFoam 

July 20, 2009, 00:57 

#52 
Senior Member
Sandy Lee
Join Date: Mar 2009
Posts: 213
Rep Power: 10 
I think, in UEqn.H, the red and blod terms which is
fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U)  fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)  fvm::laplacian(muf, U)  (fvc::grad(U) & fvc::grad(muf)) // fvc::div(muf*(fvc::interpolate(dev2(fvc::grad(U))) & mesh.Sf())) ); should be deleted, because they destroyed the original PDE. Am I right? 

August 5, 2009, 09:12 

#53  
New Member
Eric
Join Date: Jul 2009
Location: Belgium
Posts: 3
Rep Power: 9 
Dear Hamed,
have you got any progress in your simulation? I'm simulating something very similar, liquid propane release into atmospheric area. I started with dieselFoam with evaporation model on with no success, then I found your thread. If there is any experience you can share with me, it will be highly appreciated. Regards, Eric Quote:


August 6, 2009, 05:05 
compressibleLesMixtureInterPhaseChangeFoam

#54 
Member
Hamed Aghajani
Join Date: Mar 2009
Location: London, UK
Posts: 77
Rep Power: 9 
Dear Eric,
I skipped tracing vaporized liquid at that time, and haven't returned back yet. To do that, I think, there should be an extra equation for vaporized liquid. Be cause of miscible behavior of vaporized liquid and surrounding gas, say, Air, the new solver should have capability to act as "reactingFoam" do and calculate the mixture thermoPhysical property in that way. Beside this it is better to add compressibilty and turbulence modelling as well. I have no specific Idea how to do All this. Seeking for every comment and help Hamed 

August 6, 2009, 05:10 

#55 
Member
Rachel Vogl
Join Date: Jun 2009
Posts: 48
Rep Power: 9 
Hello,
Have you looked at rhoReactingFoam ? It has compressibility effect with multiphase reacting flows. I am trying to use, however yet to figure out the details. The absence of tutorial is making it tough to create a case to test this solver. 

August 6, 2009, 05:10 

#56 
Member
Rachel Vogl
Join Date: Jun 2009
Posts: 48
Rep Power: 9 
Forgot to mention it is a new solver in OF1.6


August 6, 2009, 05:13 

#57 
New Member
Eric
Join Date: Jul 2009
Location: Belgium
Posts: 3
Rep Power: 9 
Dear Hamed,
thanks for your reply and suggestion. I'll continue my simulation and see how it goes for the time being. I'll update you guys later if I made any progress. have a nice day! Eric 

August 6, 2009, 05:14 

#58 
New Member
Eric
Join Date: Jul 2009
Location: Belgium
Posts: 3
Rep Power: 9 
Hi Rachel,
Haven't update to 1.6 yet, I'm using 1.5.x. So have no idea how the new solver works! Probablly will take a look at it tomorrow. Cheers, Eric 

August 6, 2009, 06:05 
Fantastic, rhoReactingFoam!

#59 
Member
Hamed Aghajani
Join Date: Mar 2009
Location: London, UK
Posts: 77
Rep Power: 9 
Rachel,
Thank you very much for your clue, I try to shift to 1.6 and before that, OpenSUSE 11.1 . Please Keep us update of any further progress, B4N Hamed 

August 6, 2009, 06:15 

#60 
Senior Member
Sandy Lee
Join Date: Mar 2009
Posts: 213
Rep Power: 10 
Hi Rachel, are there some cases or tutorials about interPhaseChangeFoam solver in OF_1.6.


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Modelling Supersonic TwoPhase Flow with Phase Change  wes  OpenFOAM Running, Solving & CFD  8  April 26, 2016 07:21 
About phase change heat and mass transfer  Michael  FLUENT  2  February 13, 2011 02:49 
Two phase flow with phase change  Ahmad AlZoubi  CFX  1  November 26, 2008 04:59 
Twophase flow in Tjunction, multiphase of DPM?  Tony  FLUENT  2  July 8, 2008 01:26 
how to deal with phasechange heat exchanger?  cherry  FLUENT  1  April 16, 2002 21:59 