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

Multiphase flow and Phase change due to heat transferevaporation

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

Like Tree1Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   June 2, 2009, 13:58
Default
  #41
Member
 
Rachel Vogl
Join Date: Jun 2009
Posts: 48
Rep Power: 8
Rachel is on a distinguished road
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.
Rachel is offline   Reply With Quote

Old   June 9, 2009, 04:55
Default
  #42
Senior Member
 
Sandy Lee
Join Date: Mar 2009
Posts: 207
Rep Power: 9
sandy is on a distinguished road
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,p-pv]
Mdest=(Cc/tInf)rho2*gamma2[1-gamma]

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


sandy is offline   Reply With Quote

Old   June 21, 2009, 08:27
Default
  #43
Senior Member
 
Sandy Lee
Join Date: Mar 2009
Posts: 207
Rep Power: 9
sandy is on a distinguished road
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
sandy is offline   Reply With Quote

Old   June 30, 2009, 11:59
Default
  #44
Senior Member
 
isabel
Join Date: Apr 2009
Location: Spain
Posts: 171
Rep Power: 8
isabel is on a distinguished road
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 ?
isabel is offline   Reply With Quote

Old   June 30, 2009, 12:38
Default
  #45
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,912
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by isabel View Post
Does anybody know what means " oneField() " ?
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
gschaider is offline   Reply With Quote

Old   July 1, 2009, 08:28
Default
  #46
Senior Member
 
Sandy Lee
Join Date: Mar 2009
Posts: 207
Rep Power: 9
sandy is on a distinguished road
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 full-wet 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.
sandy is offline   Reply With Quote

Old   July 2, 2009, 22:54
Default
  #47
Senior Member
 
Sandy Lee
Join Date: Mar 2009
Posts: 207
Rep Power: 9
sandy is on a distinguished road
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.
sandy is offline   Reply With Quote

Old   July 7, 2009, 05:06
Default
  #48
Senior Member
 
Sandy Lee
Join Date: Mar 2009
Posts: 207
Rep Power: 9
sandy is on a distinguished road
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.
sandy is offline   Reply With Quote

Old   July 16, 2009, 05:54
Default
  #49
Senior Member
 
isabel
Join Date: Apr 2009
Location: Spain
Posts: 171
Rep Power: 8
isabel is on a distinguished road
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?
isabel is offline   Reply With Quote

Old   July 16, 2009, 09:22
Default
  #50
Senior Member
 
Sandy Lee
Join Date: Mar 2009
Posts: 207
Rep Power: 9
sandy is on a distinguished road
Quote:
Originally Posted by isabel View Post
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?
vDotvAlphal = (Cv*rho2/(0.5*rho1*sqr(UInf)*tInf))*min(p - pSat, 0)*{1.0/rho1 - gamma*(1.0/rho1 - 1.0/rho2)}

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)/(p-pSat)} *(1.0/rho1 - 1.0/rho2)

Am I right?

Last edited by sandy; July 20, 2009 at 00:49.
sandy is offline   Reply With Quote

Old   July 17, 2009, 02:56
Default
  #51
Senior Member
 
isabel
Join Date: Apr 2009
Location: Spain
Posts: 171
Rep Power: 8
isabel is on a distinguished road
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 partly-relaxed version of U.
It is explained in this link:

Trying to figure out the details of simpleFoam


isabel is offline   Reply With Quote

Old   July 20, 2009, 00:57
Default
  #52
Senior Member
 
Sandy Lee
Join Date: Mar 2009
Posts: 207
Rep Power: 9
sandy is on a distinguished road
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?
sandy is offline   Reply With Quote

Old   August 5, 2009, 09:12
Default
  #53
New Member
 
Eric
Join Date: Jul 2009
Location: Belgium
Posts: 3
Rep Power: 7
eric_zh is on a distinguished road
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:
Originally Posted by haghajani View Post
Dear All,

How can I trace the vaporized liquid (changed phase from 1->0) in to air.
I am trying to study the effect of flash evaporation in high pressure liquid Hydrogen release. I used interPhaseChangeFoam for a simple case and the results seems reasonable. Now I want also see the vapor of hydrogen inside surrounding air.

Best regards,
Hamed Aghajani
hamed.aghajani@gmail.com
h.aghajani@kingston.ac.uk
eric_zh is offline   Reply With Quote

Old   August 6, 2009, 05:05
Lightbulb compressibleLesMixtureInterPhaseChangeFoam
  #54
Member
 
Hamed Aghajani
Join Date: Mar 2009
Location: London, UK
Posts: 77
Rep Power: 8
haghajani is on a distinguished road
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
haghajani is offline   Reply With Quote

Old   August 6, 2009, 05:10
Default
  #55
Member
 
Rachel Vogl
Join Date: Jun 2009
Posts: 48
Rep Power: 8
Rachel is on a distinguished road
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.
Rachel is offline   Reply With Quote

Old   August 6, 2009, 05:10
Default
  #56
Member
 
Rachel Vogl
Join Date: Jun 2009
Posts: 48
Rep Power: 8
Rachel is on a distinguished road
Forgot to mention it is a new solver in OF-1.6
Rachel is offline   Reply With Quote

Old   August 6, 2009, 05:13
Default
  #57
New Member
 
Eric
Join Date: Jul 2009
Location: Belgium
Posts: 3
Rep Power: 7
eric_zh is on a distinguished road
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
eric_zh is offline   Reply With Quote

Old   August 6, 2009, 05:14
Default
  #58
New Member
 
Eric
Join Date: Jul 2009
Location: Belgium
Posts: 3
Rep Power: 7
eric_zh is on a distinguished road
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
eric_zh is offline   Reply With Quote

Old   August 6, 2009, 06:05
Thumbs up Fantastic, rhoReactingFoam!
  #59
Member
 
Hamed Aghajani
Join Date: Mar 2009
Location: London, UK
Posts: 77
Rep Power: 8
haghajani is on a distinguished road
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
haghajani is offline   Reply With Quote

Old   August 6, 2009, 06:15
Default
  #60
Senior Member
 
Sandy Lee
Join Date: Mar 2009
Posts: 207
Rep Power: 9
sandy is on a distinguished road
Hi Rachel, are there some cases or tutorials about interPhaseChangeFoam solver in OF_1.6.
sandy is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 00:14.