[yangzhou]

Hi Foamers,

I am a newcomer of waves2Foam. I have confronted some problems in generating waves(OpenFoam 4.0, Ubuntu 17.04). I am getting approximately 5%-6% reduction on the peak and valley of the wave, the error is little bit more than I expected.
1.JPG
2.JPG
Wave properties input

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.5                                   |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      environmentalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

seaLevel    0.00;

// A list of the relaxation zones in the simulation. The parameters are given
// in <name>Coeffs below.
relaxationNames (inlet outlet);

initializationName outlet;

inletCoeffs
{
    // Wave type to be used at boundary "inlet" and in relaxation zone "inlet"
    waveType    stokesFirst;  
    
    // Ramp time of 2 s
    Tsoft       2;

    // Water depth at the boundary and in the relaxation zone
    depth       2;

    // Wave period
    period      2.0;

    // Phase shift in the wave
    phi         0.000000;

    // Wave number vector, k. 
    direction  (1.0 0.0 0.0);

    // Wave height
    height      0.04;
    
    // Specifications on the relaxation zone shape and relaxation scheme
    relaxationZone
    {
        relaxationScheme Spatial;
        relaxationShape  Rectangular;
        beachType        Empty;
    
        relaxType   INLET;
        startX      (0 0.0 -1);
        endX        (5 0.0  1);
        orientation      (1.0 0.0 0.0);
    }
};

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

    relaxationZone
    {
        relaxationScheme Spatial;
        relaxationShape  Rectangular;
        beachType        Empty;    

        relaxType   OUTLET;
        startX      (13 0.0 -1);
        endX        (18 0.0  1);
        orientation      (1.0 0.0 0.0);
    }
};



// ************************************************************************* //


Blockmesh Dictionary
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.5                                   |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

The blockmesh dic

convertToMeters 1;

vertices        
(
    (  0 -2 0 )
    ( 18 -2 0 )
    (  0  1 0 )
    ( 18  1 0 )
                 
    (  0 -2 0.1 )
    ( 18 -2 0.1 )
    (  0  1 0.1 )
    ( 18  1 0.1 )                    
);

blocks          
(
    hex (0 1 3 2 4 5 7 6) ( 150 450 1 ) simpleGrading (1 1 1)
);

edges           
(
);

patches         
(
    patch inlet
    (
        (0 4 6 2)
    )
    wall bottom 
    (
        (0 1 5 4)
    )
    patch outlet
    (
        (1 5 7 3)
    )

    patch atmosphere 
    (
        (2 3 7 6)
    )
    
    empty frontBack
    (
        (0 1 3 2)
        (4 5 7 6)
    )
);

mergePatchPairs
(
);

// ************************************************************************* //

Can anybody guide me how to narrow the gap? Many thanks

[MehdiCFD]

me too i obtained the same range of error with Airy (stokesFirst) about 6%, but you have to refine your mesh bellow sea surface and verify tha you don't make a brutal grading and let a large distance between the end of the refined mesh and the max amplitude size.
me i'm using AR=9 delta z=0.0033m delta x,y=0.03m near the surface (H=0.1 m i use 0.35 m refined mesh bellow surface).

i hope that will be usefull for you and good luck.
Best regards,
Mehdi

[yangzhou]

Dear Mehdi,
I have tried different mesh method by divide the wavelength（λ/20 for example）and the amplitude of the wave (H/6). By the results is not satisfying. It seems different mesh method do not affect the results of wave generation.
Capture.JPG
Any tips would be helpful to me.
Best regards,
Yang

Originally Posted by MehdiCFD

me too i obtained the same range of error with Airy (stokesFirst) about 6%, but you have to refine your mesh bellow sea surface and verify tha you don't make a brutal grading and let a large distance between the end of the refined mesh and the max amplitude size.
me i'm using AR=9 delta z=0.0033m delta x,y=0.03m near the surface (H=0.1 m i use 0.35 m refined mesh bellow surface).

i hope that will be usefull for you and good luck.
Best regards,
Mehdi

[MehdiCFD]

Dear Yangzhou,

so it's not a mesh problem, may be you have to use turbulence modelling if your Reynolds number is high then 10^5 then you should use it to modelise the production of turbulent energy may be it will be decreasing your relative error, i suggeste you to use K-w SST turbulence model
if it that doen't resolve your problem you have to look into your fvSchemes div Schemes, the Gauss linearUpwind Schemes of div(phi.U) is more efficient and less dessipative then limitedLinear schemes, you have to read different paper that use waves2Foam solver.

Best regards,
Mehdi