[Taataa]

Hi all,

I am modeling a two phase (air and water) flow in a container with compressed air at the bottom and a column of water on top. The top boundary is open to atmosphere and the rest are walls. I set the initial condition for the pressure field using setFields as follows (two times atmospheric pressure):

boxToCell
    {
        box (0 0 -1) (0.1 0.5 1);
        fieldValues
        (
    		volScalarFieldValue p_rgh 303975
        );
    }

I am using multiphaseEulerFoam solver and set the top boundary conditions for U and p_rgh as follows:

p_rgh:
    atmosphere
    {
        type            totalPressure;
        p0              $internalField;
        U               U.air;
        phi             phi.air;
        value           $internalField;
    }

U:
    atmosphere
    {
        type            pressureInletOutletVelocity;
        value           uniform (0 0 0);
    }

The problem is that no matter what value I set for the initial pressure the results are the same. A bubble will form and rises. The water column should move upward due to effect of compressed air at the bottom. Any suggestions is appreciated.