[samtheman23]

All,

I built a '2D' numerical wave tank (1 element wide) with a fixed, truncated cylinder in the middle (static mesh).

I know I am breaking protocol by not posting everything to make it reproducible but I was wondering if anyone had encountered a similar phenomenon.

Here is a video sped up 10x: https://youtu.be/o-_OqZrtPJk
When I initialize my domain the pressure (and velocity is 0) everywhere. But on the first time step you see these 'pressure waves' above the MWL emanate from the cylinder outward.

Here is my waveProperties:

seaLevel            0;

relaxationNames     ( inlet outlet );

initializationName  outlet;

inletCoeffs
{
    waveType            stokesFirst;
    Tsoft               20;
    depth               95;
    period              8;
    direction           ( 1 0 0 );
    phi                 0;
    height              4;
    waveNumber          (0.0628806 0 0);
    omega               0.785398;

    relaxationZone
    {
        relaxationScheme    Spatial;
        relaxationShape     Rectangular;
        beachType           Empty;
        relaxType           INLET;
        startX              ( -650 0 0 );
        endX                ( -620 0 0.1 );
        orientation         ( 1 0 0 );
    }
}

outletCoeffs
{
    waveType            potentialCurrent;
    U                   ( 0 0 0 );
    Tsoft               20;

    relaxationZone
    {
        relaxationScheme    Spatial;
        relaxationShape     Rectangular;
        beachType           Empty;
        relaxType           OUTLET;
        startX              ( 640 0 0 );
        endX                ( 650 0 0.1 );
        orientation         ( 1 0 0 );
    }
}

Here is my 0/p_rgh:

boundaryField
{
    cylinder
    {
        type            zeroGradient;
    }
    front
    {
        type            empty;
    }
    bottom
    {
        type            zeroGradient;
    }
    inlet
    {
        type            zeroGradient;
    }
    back
    {
        type            empty;
    }
    outlet
    {
        type            zeroGradient;
    }
    atmosphere
    {
        type            totalPressure;
        rho             none;
        psi             none;
        gamma           1;
        p0              uniform 0;
        value           uniform 0;
    }
}

Here is my 0/U:

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

internalField   uniform (0 0 0);

boundaryField
{
    cylinder
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    front
    {
        type            empty;
    }
    bottom
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    inlet
    {
        type            waveVelocity;
        refValue        uniform (0 0 0);
        refGradient     uniform (0 0 0);
        valueFraction   uniform 1;
        value           uniform (0 0 0);
    }
    back
    {
        type            empty;
    }
    outlet
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    atmosphere
    {
        type            pressureInletOutletVelocity;
        value           uniform (0 0 0);
    }
}

My only guess is that there is some flux condition on the inlet/outlet. I also know that a common fix is an AR of 1 at the waterline. Or maybe going to StokesSecond? Any other ideas?