[xucong1228]

Hello Foamers
I am new to Openfoam and I am doing a simulation about drainage,you can see model I use from the attachment

length=5.5m
width=3m
depth=0.1m

The left side is inlet (Q=20 L/s)
The right side is outlet and there is a gully on the bottom.(we can assume the gully is the second outlet)
The top floor is atmosphere.

I used the LTSInterFoam for simulation because I heard this Foam is faster than interFoam.However, I found some very strange result after simulation:

1.As you can see in the attachment, the water cannot get into the gully.
2.The flow rate of inlet is not equal to the flow rate of two outlets.

I think it should be the problem of boundary condition, but I don't know which one it is.
So I listed my files U/P_pgh/alpha, I hope someone could help me.

P_pgh

dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type           buoyantPressure;
        value          uniform 0;
        
    }

    outlet
    {
        type           zeroGradient;
    }

    outlet2
    {
        type           zeroGradient;
    }

    sides
    {
        type           buoyantPressure;
        value          uniform 0;
    }

    atmosphere
    {
        type           totalPressure;
        p0              uniform 0;
        U               U;
        phi             phi;
        rho             rho;
        psi             none;
        gamma           1;
        value           uniform 0;
    }

    lowerwall
    {
        type            buoyantPressure;
        value          uniform 0;
    }
}

U

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type            flowRateInletVelocity;
        volumetricFlowRate  constant 20e-3;           
        rhoInlet            1000;
    }

    outlet
    {
        type            zeroGradient;
    }

    outlet2
    {
        type            zeroGradient;                  
    }

    sides
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }

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

    lowerwall
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
}

alpha

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            inletOutlet;
        inletValue     uniform 1;
        value           uniform 1;
    }

    outlet
    {
        type            zeroGradient;
    }

    sides
    {
        type            zeroGradient;
    }

    atmosphere
    {
        type            inletOutlet;
        inletValue     uniform 0;
        value           uniform 0;
    }

    lowerwall
    {
        type            zeroGradient;
    }

    outlet2
    {
       type            zeroGradient;
    }
}

Lastly,can anyone tell me how to calculate the flow rate of water in this multiphase flow.I used the 'calculator' and 'integrate variable' from paraview.However i don't think it is precise.
Any help would be great, thank you so much!