[pixarzhang]

I ran the "3Dwaves" case. the only change I made is set Wind to (1 0 0 ). the boundary field Uy of wind is almost 0，but the internal field of wind near boundary is very differ！
the figure is attached.


I think it's not the solver's but the boundary condition's problem, when i adopted groovyBC in interFoam case, I also came into the discontinuity problem.

Does anyone know how to solve the problem?


the boundary field,the upside Uy is 0.

the internal Field,the upside Uy is not 0,especially some part have value up to 0.4 or low to -1

My object concern on the wind field upon the wave. I hope to get an natural wind field , not the corrected field with Utarget.

the size of my computational zone is x*y*z=10m*7m*2m.at first, i set a relaxationZone from the inlet boundary.
but after several seconds.the Uy diffused as:

DILUPBiCG:  Solving for Uy, Initial residual = 0.00810394, Final residual = 1.05103e+19, No Iterations 1001

then deltaT becomes very little,the computation collapsed at last.do you know where goes wrong, is the case can't use this solver or the BC is not proper?

the relaxationZone is set as

relaxationNames     ( inlet );

wind                ( 1 0 0 );

initializationName  inlet;

inletCoeffs
{
    waveType            stokesFirst;
    Tsoft               0;
    depth               2.75;
    period              1;
    direction           ( 1 0 0 );
    phi                 0;
    height              0.2;
    waveNumber          (3.1415926 0 0);
    omega               6.28319;

    relaxationZone
    {
        relaxationScheme    Spatialwater;
        relaxationShape     Rectangular;
        beachType           Empty;
        relaxType           INLET;
        startX              ( 0 0.7 0 );
        endX                ( 5 0.7 2 );
        orientation         ( 1 0 0 );
    }
}

the initial U field boundary is set as

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }
    
    outlet
    {
        type            zeroGradient;
    }

    cyclic2
    {
        type            cyclic;
    }

    bottom
    {
        type            slip;
    }
    
    atmosphere
    {
        type            slip;
    }
    
    defaultFaces
    {
        type            zeroGradient;
    }
}

alpha.water field boundary

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }
    outlet
    {
        type            zeroGradient;
    }
    cyclic2
    {
        type            cyclic;
    }

    bottom
    {
        type            fixedValue;
        value           uniform 1;
    }

    atmosphere
    {
        type            fixedValue;
        value           uniform 0;
    }

    defaultFaces
    {
        type            empty;
    }
}

p_rhg field boundary is set as

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }
    outlet
    {
        type            zeroGradient;
    }
    bottom
    {
        type            slip;
    }
    atmosphere
    {
        type            zeroGradient;
    }
    cyclic2_half0
    {
        type            cyclic;
    }
    cyclic2_half1
    {
        type            cyclic;
    }
}

below is the U field at 10.6s, the diffusion occured over the wave crest out of the relaxationZone (the relaxationZone is at the right part)