[fmerk]

Hello everyone,

I'm currently working on a spillway simulation where I want to set my inlet boundary as a constant flow height (h=3m). I want the flow to be pressure-driven by a hydrostatic pressure distribution at the named flow height. I chose the interFoam solver.

I divided the inlet patch into a inletWater (alpha=1) and inletAir (alpha=0). Now, I'd like to set up a hydrostatic pressure for the inletWater patch:

inletWater
{

type groovyBC;
valueExpresion "9810.*(3-pos().z)+100000";
fractionExpression "1";
rho rho;
psi none;
gamma 1;
p0 uniform 100000;
value uniform 100000;

}

Checking this setup (and many alternatives) does not provide a hydrostatic pressure distribution, but a single value for the whole patch (=100000).


I also prescribed a water phase with a pressure distribution for the internal field using funkySetFields which works successfully.

expressions
(
pressureWater
{
field p_rgh; //field to initialise
expression "9810.*(2-pos().z)+100000";
condition "(pos().z<2) && (pos().z>=0)"; //ranging from 2 meter to 0 meter
keepPatches 1;
}

);

How can I setup my inlet to have a constant hydrostatic pressure for a given flow height?

I'm probably missing something fundamental. Since I'm googling since weeks now, I'm getting frustrated.



Kind regards,

fmerk