# CavitatingFoam psil psiv

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

 March 10, 2008, 06:22 Hi! Can anyone explain me the #1 New Member   Marco Zecchi Join Date: Mar 2009 Location: Italy Posts: 8 Rep Power: 14 Hi! Can anyone explain me the meaning of "psiv" and "psil" in cavitatingFoam/nozzle2D/constant/thermodynamicProperties/? Thank you!

 March 10, 2008, 12:09 Marco, As in most OpenFOAM ap #2 Member   David P. Schmidt Join Date: Mar 2009 Posts: 71 Rep Power: 14 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

 March 10, 2008, 13:30 Thank you for your quickly rep #3 New Member   Marco Zecchi Join Date: Mar 2009 Location: Italy Posts: 8 Rep Power: 14 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

 March 12, 2008, 11:13 Marco, The equations for ma #4 Member   David P. Schmidt Join Date: Mar 2009 Posts: 71 Rep Power: 14 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

 July 31, 2009, 10:07 #5 New Member   varun Join Date: Jul 2009 Location: Braunschweig, Germany Posts: 28 Rep Power: 14 hi, I was wondering how to calculate 'psil' and 'psiv'? does anyone know the value of psiv for a vapour like air? pls help!

 October 7, 2009, 06:10 #6 New Member   Jml Join Date: Mar 2009 Posts: 23 Rep Power: 14 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

December 9, 2009, 15:39
#7
New Member

Sukanta Rakshit
Join Date: Jun 2009
Posts: 16
Rep Power: 14
Quote:
 Originally Posted by jml 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)

June 3, 2014, 08:06
#8
Member

Join Date: Jul 2011
Posts: 54
Rep Power: 12
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 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.

February 17, 2020, 17:23
Calculating psiL and psiV
#9
Member

Join Date: Mar 2019
Posts: 79
Rep Power: 4
Quote:
 Originally Posted by A_Pete 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