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

CavitatingFoam psil psiv

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

Like Tree4Likes
  • 4 Post By A_Pete

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 10, 2008, 06:22
Default Hi! Can anyone explain me the
  #1
New Member
 
Marco Zecchi
Join Date: Mar 2009
Location: Italy
Posts: 8
Rep Power: 14
kyubi is on a distinguished road
Hi!
Can anyone explain me the meaning of "psiv" and "psil" in cavitatingFoam/nozzle2D/constant/thermodynamicProperties/?

Thank you!
kyubi is offline   Reply With Quote

Old   March 10, 2008, 12:09
Default Marco, As in most OpenFOAM ap
  #2
Member
 
David P. Schmidt
Join Date: Mar 2009
Posts: 71
Rep Power: 14
schmidt_d is on a distinguished road
Marco,
As in most OpenFOAM applications, psi in cavitatingFoam is the derivative of density with respect to pressure. Here, the v and l appended to psi refer to the vapor and liquid, respectively.
-DPS
schmidt_d is offline   Reply With Quote

Old   March 10, 2008, 13:30
Default Thank you for your quickly rep
  #3
New Member
 
Marco Zecchi
Join Date: Mar 2009
Location: Italy
Posts: 8
Rep Power: 14
kyubi is on a distinguished road
Thank you for your quickly reply David!

I'm also trying to figure out which are the governing equations used in cavitatingFoam, but I'm not able to understand it only from the source code...
Can you help me again?

Marco
kyubi is offline   Reply With Quote

Old   March 12, 2008, 11:13
Default Marco, The equations for ma
  #4
Member
 
David P. Schmidt
Join Date: Mar 2009
Posts: 71
Rep Power: 14
schmidt_d is on a distinguished road
Marco,

The equations for mass & momentum are the standard single phase, variable density versions, since cavitatingFoam is using a pseudo-fluid approach. The user can choose various compressibility models (see the compressibilityModels sub-directory in the source). I am most familiar with Wallis' model. The derivation is given in his text book, which is now out of print. The compressibility model sets the value of psi as a function of mixture fraction, which feeds into the pressure solution. Unfortunately, I find that you can run into serious stability problems when the density ratio between liquid and vapor is extreme.

I once tried an approach where you analytically integrate Walllis' relationship. It worked well for very high speed flow and small length scales, but required specialized non-reflecting BC's.

@Article{Schmidt1999b,
author = {D.P. Schmidt and C.J. Rutland and M.L. Corradini},
title = "{A Fully Compressible Model of Cavitating Flow}",
journal = {Atomization and Sprays},
year = {1999},
volume = {9},
}

The openFOAM implementation is more general and can run at low Mach #.

-David
schmidt_d is offline   Reply With Quote

Old   July 31, 2009, 10:07
Default
  #5
New Member
 
varun
Join Date: Jul 2009
Location: Braunschweig, Germany
Posts: 28
Rep Power: 14
vw.cfd is on a distinguished road
hi,

I was wondering how to calculate 'psil' and 'psiv'?
does anyone know the value of psiv for a vapour like air?

pls help!
vw.cfd is offline   Reply With Quote

Old   October 7, 2009, 06:10
Default
  #6
jml
New Member
 
Jml
Join Date: Mar 2009
Posts: 23
Rep Power: 14
jml is on a distinguished road
Hello Schmidt, I have seen that you talk about the Wallis model. I'm setting the thermodynamicproperties file, but I have three parameters of density. żDo you know what's the difference between rholSat, rhovSat and rhoMin?

Thanks
jml is offline   Reply With Quote

Old   December 9, 2009, 15:39
Default
  #7
New Member
 
Sukanta Rakshit
Join Date: Jun 2009
Posts: 16
Rep Power: 14
srakshit is on a distinguished road
Quote:
Originally Posted by jml View Post
Hello Schmidt, I have seen that you talk about the Wallis model. I'm setting the thermodynamicproperties file, but I have three parameters of density. żDo you know what's the difference between rholSat, rhovSat and rhoMin?

Thanks
rholSat --> Liquid Density at Saturation Point (i.e. at equillibrium)

rhovSat --> Vapour Density at Saturation Point (i.e. at equillibrium)

rhoMin --> very low value of density used to keep the density positive (value of 0.1 is fine)
srakshit is offline   Reply With Quote

Old   June 3, 2014, 08:06
Default
  #8
Member
 
Join Date: Jul 2011
Posts: 54
Rep Power: 12
A_Pete is on a distinguished road
I know it's an old thread, but since I was looking for this myself I felt like I should at least answer to this post, after having found the answer.

Quote:
Originally Posted by vw.cfd View Post
hi,

I was wondering how to calculate 'psil' and 'psiv'?
does anyone know the value of psiv for a vapour like air?

pls help!

'psil' can be calculated as psil = 1/a_l^2 , where a_l is the sonic speed of the liquid phase. Same goes for 'psiv' and the appropriate sonic speed of the vapour phase.
A_Pete is offline   Reply With Quote

Old   February 17, 2020, 17:23
Default Calculating psiL and psiV
  #9
Member
 
Join Date: Mar 2019
Posts: 79
Rep Power: 4
mm66 is on a distinguished road
Quote:
Originally Posted by A_Pete View Post
I know it's an old thread, but since I was looking for this myself I felt like I should at least answer to this post, after having found the answer.




'psil' can be calculated as psil = 1/a_l^2 , where a_l is the sonic speed of the liquid phase. Same goes for 'psiv' and the appropriate sonic speed of the vapour phase.
Hi Foamers,

I was going through the tutorial for cavitatingFoam. It mentions [0 -2 2 0 0] as the unit for both psiL and psiV. So it is:
Code:
rho * isentropic compressibility (Bs)
where Bs*rho is the inverse of (speed of sound)^2:

Code:
Bs*rho = 1/c^2
If isothermal compressibility (Bt) is instead available, then Bs can be calculated using:

Code:
Cp/Cv = Bt/Bs
I appreciate if you please correct me if I am wrong here.

Regards,
MJ
mm66 is offline   Reply With Quote

Reply

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
CavitatingFoam Equations morteza OpenFOAM Running, Solving & CFD 13 October 5, 2020 13:15
CavitatingFoam for Diesel injector nozzles javier_salvador OpenFOAM Running, Solving & CFD 0 July 16, 2008 10:37
CavitatingFoam with a moving mesh mass conservation error idosil OpenFOAM Running, Solving & CFD 3 November 27, 2007 18:27
Problem with cavitatingFoam clo OpenFOAM Running, Solving & CFD 1 June 7, 2007 12:53


All times are GMT -4. The time now is 05:26.