[kjetil]

Originally Posted by ngj

I am not aware of any other solver. However, if you are can share your setup, I will try to give it a look and see, if I can spot where this un-physical pressure originates from.

Thanks. Well, I paste here what I believe is most crucial. If needed the entire case tree can be found online here http://folk.ntnu.no/kjetilbi/openfoam/nr5-auto2/ .

Content of 0/U:

internalField uniform (0 0 0);

boundaryField

{
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
pipewall
{
type fixedValue;
value uniform (0 0 0);
}
}

Content of 0/P:

internalField uniform 0;

boundaryField
{
inlet
{
type fixedValue;
value uniform 0;
}
outlet
{
type zeroGradient;
}
pipewall
{
type zeroGradient;

}
}

And a snippet from transportProperties:

// liq
phase2
{
transportModel Newtonian;
nu nu [0 2 -1 0 0 0 0] 1e-6;
rho rho [1 -3 0 0 0 0 0] 1000;
rho0 rho0 [1 -3 0 0 0 0 0] 1000;
psi psi [ 0 -2 2 0 0 ] 1e-05;
}

// gas
phase1
{
transportModel Newtonian;
nu nu [0 2 -1 0 0 0 0] 1.2e-5;
rho rho [1 -3 0 0 0 0 0] 1.2;
rho0 rho0 [1 -3 0 0 0 0 0] 1.2;
psi psi [ 0 -2 2 0 0 ] 1e-05;
}

pMin pMin [ 1 -1 -2 0 0 0 0 ] 0;

sigma sigma [ 1 0 -2 0 0 0 0 ] 0.07;

Basically, I only "know" the outlet pressure - which should be 1 atm. The inlet is not pressurized from outside, but instead the (negative) hydrostatic pressure is supposed "pull the water" from the outside reservoir.

I would very much appreciate any comments on how to improve this approach.

Thanks.