[WaterHammer1985]

I am trying to adapt the weirOverflow tutorial case (interFoam) to a similar domain with a mesh converted from fluent. Schematic of domain is below. Case files are located here https://ufile.io/w2xaw (couldn't upload below).



When I setFields, the alpha.water looks good but the pressure field in the domain does not look good. The water depth is 0.54m (5300Pa) but the pressure upstream of the weir is 5.3e+09Pa (see image). As time progresses in the simulation, the area of high pressure shows odd results (see image).



I've run checkMesh and it looks ok to me. I'm using the same boundary conditions as the weirOverflow tutorial (note that inlet_air is being modeled as a wall).



If anyone has advice on how I should be troubleshooting this it would be greatly appreciated. Thanks!

















0/U

FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#include        "include/initialConditions"

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

internalField   uniform (0 0 0);

boundaryField
{
    inlet_flow
    {
        type            variableHeightFlowRateInletVelocity;
        flowRate        $inletFlowRate;
        alpha           alpha.water;
        value           uniform (0 0 0);
    }

    outlet
    {
        type            zeroGradient;
    }

    (sides|bottom|inlet_air|vent)
    {
        type            noSlip;
    }

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

    defaultFaces
    {
        type            empty;
    }
}

0/p_rgh

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#include        "include/initialConditions"

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

internalField   uniform $pressure;

boundaryField
{
    inlet_flow
    {
        type            zeroGradient;
    }

    outlet
    {
        type            zeroGradient;
    }

    (sides|inlet_air|bottom|vent)
    {
        type            fixedFluxPressure;
    }
    top
    {
        type            totalPressure;
        p0              uniform 0;
    }

    defaultFaces
    {
        type            empty;
    }
}

0/k

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#include        "include/initialConditions"

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

internalField   uniform $turbulentKE;

boundaryField
{
    inlet_flow
    {
        type            fixedValue;
        value           $internalField;
    }

    outlet
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }

    (sides|inlet_air|bottom|vent)
    {
        type            kqRWallFunction;
        value           $internalField;
    }

    top
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }

    defaultFaces
    {
        type            empty;
    }
}

0/epsilon

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#include        "include/initialConditions"

dimensions      [0 2 -3 0 0 0 0];

internalField   uniform $turbulentEpsilon;

boundaryField
{
    inlet_flow
    {
        type            fixedValue;
        value           $internalField;
    }

    outlet
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }

    (sides|inlet_air|bottom|vent)
    {
        type            epsilonWallFunction;
        value           $internalField;
    }
   
    top
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }

    defaultFaces
    {
        type            empty;
    }
}

0/include/initialconditions

inletFlowRate        0.007506;
pressure             0;
turbulentKE          4.14e-03;
turbulentEpsilon     4.39e-05;

constant/setFieldsDict

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      setFieldsDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

defaultFieldValues
(
    volScalarFieldValue alpha.water 0
);

regions
(
    boxToCell
    {
        box (-100 -50 -100) (100 1.75 0.54);

        fieldValues
        (
            volScalarFieldValue alpha.water 1
            volVectorFieldValue U (0 0.139 0)
        );
    }
);