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

buoyantSimpleFoam and watertank

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

Like Tree25Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 18, 2022, 08:15
Default p_rgh and non-uniform density
  #101
New Member
 
Fatih Tüzün
Join Date: Dec 2022
Posts: 1
Rep Power: 0
FatihT is on a distinguished road
Hello,
I'm facing a similar problem with non-uniform density over a water tank with supercritical water. The density is in a range of 85 to 1012 kg/m³.

I did the initialization with codeStream. (see below)

If you calculate the density from the pressure difference at the very bottom,
(p-p_rgh)/(gh) density corresponds to the local density.

Is there a way to use domainIntegrate to integrate the density up to the current height?
Or something like: p_hyd += rho*g*(h-h[-1])

buoyantFoam/pEqn.H:

// For steady closed-volume compressible cases adjust the pressure level
// to obey overall mass continuity
if (adjustMass && !thermo.incompressible())
{
p += (initialMass - fvc::domainIntegrate(thermo.rho()))
/fvc::domainIntegrate(psi);
p_rgh = p - rho*gh - pRef;
p_rgh.correctBoundaryConditions();
}
Attached Images
File Type: jpg prgh.jpg (52.4 KB, 12 views)
FatihT 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



All times are GMT -4. The time now is 08:25.