[show]

Hi, there,

An error was appeared when I run the waveFlume, the error was appeared during the calculation and the reminder of the error is as follow,

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigSegv::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::sampledSurfaceElevation::sampleIntegrateAndW rite(Foam::sampledSurfaceElevation::fieldGroup<dou ble>&) at ??:?
#4 Foam::OutputFilterFunctionObject<Foam::sampledSurf aceElevation>::execute(bool) at ??:?
#5 Foam::functionObjectList::execute(bool) at ??:?
#6 Foam::Time::run() const at ??:?
#7 main at ??:?
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9 ? at ??:?

However, when I delete the following function in controlDict/system, the error will disappeared.
functions
{
#includeIfPresent "../waveGaugesNProbes/surfaceElevationAnyName_controlDict";
}

But, I want to get the surface elevation, so I can't delete it. Anyone can help me to solve this problem? Thanks a lot!

[ngj]

Good morning Shengchao,

Are you getting this error, when you are running the standard wave flume test or have you made modifications?

My best guess is that your wave gauges are outside of the computational domain.

Kind regards,

Niels

[show]

Hi, Niels,

Thank you for your quick reply, I modify the original waveFlume, In waveProperties/Constant the solitary wave is used

relaxationNames (inlet outlet);

initializationName outlet;

inletCoeffs
{
    waveType    solitaryFirst;
    height    0.08;
    depth       0.2;
    direction   (1 0 0);
    x0          (1 0 0);
    
    // 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      (26 0.0 -1);
        endX        (30 0.0  1);
        orientation      (1.0 0.0 0.0);
    }
};

In blockMeshDict, a rectangle bridge model is put in the wave flume, the code is as follows,

convertToMeters 0.001;

vertices        
(
    (  0 -200 0 )
    ( 15000 -200 0 )
    ( 15324 -200 0 )
    ( 30000 -200 0 )
    (  0 0 0 )
    ( 15000 0 0 )
    ( 15324 0 0 )
    ( 30000 0 0 )
    (  0 40 0 )
    ( 15000 40 0 )
    ( 15324 40 0 )
    ( 30000 40 0 )
    (  0 200 0 )
    ( 15000 200 0 )
    ( 15324 200 0 )
    ( 30000 200 0 )
    
                 
    (  0 -200 0.1 )
    ( 15000 -200 0.1 )
    ( 15324 -200 0.1 )
    ( 30000 -200 0.1 )
    (  0 0 0.1 )
    ( 15000 0 0.1 )
    ( 15324 0 0.1 )
    ( 30000 0 0.1 )
    (  0 40 0.1 )
    ( 15000 40 0.1 )
    ( 15324 40 0.1 )
    ( 30000 40 0.1 )
    (  0 200 0.1 )
    ( 15000 200 0.1 )
    ( 15324 200 0.1 )
    ( 30000 200 0.1 )                
);

blocks          
(
    hex (0 1 5 4 16 17 21 20) ( 300 50 1 ) simpleGrading (0.5 1 1)
    hex (4 5 9 8 20 21 25 24) ( 300 10 1 ) simpleGrading (0.5 1 1)
    hex (8 9 13 12 24 25 29 28) ( 300 40 1 ) simpleGrading (0.5 1 1)
    hex (1 2 6 5 17 18 22 21) ( 20 50 1 ) simpleGrading (1 1 1)
    hex (9 10 14 13 25 26 30 29) ( 20 40 1 ) simpleGrading (1 1 1)
    hex (2 3 7 6 18 19 23 22) ( 300 50 1 ) simpleGrading (2 1 1)
    hex (6 7 11 10 22 23 27 26) ( 300 10 1 ) simpleGrading (2 1 1)
    hex (10 11 15 14 26 27 31 30) ( 300 40 1 ) simpleGrading (2 1 1)
);

edges           
(
);

patches         
(
    patch inlet
    (
        (0 16 20 4)
        (4 20 24 8)
        (8 24 28 12)
    )
    wall bottom 
    (
        (0 1 17 16)
        (1 2 18 17)
        (2 3 19 18)
    )
    patch outlet
    (
        (3 19 23 7)
        (7 23 27 11)
        (11 27 31 15)
    )

    patch atmosphere 
    (
        (12 13 29 28)
        (13 14 30 29)
        (14 15 31 30)
    )

    wall bridgemodel
    (
        (5 21 25 9)
        (6 22 26 10)
        (5 6 22 21)
        (9 10 26 25)
    )

    empty frontBack
    (
        (0 1 5 4)
        (4 5 9 8)
        (8 9 13 12)
        (1 2 6 5)
        (9 10 14 13)
        (2 3 7 6)
        (6 7 11 10)
        (10 11 15 14)
        (16 17 21 20)
        (20 21 25 24)
        (24 25 29 28)
        (17 18 22 21)
        (25 26 30 29)
        (18 19 23 22)
        (22 23 27 26)
        (26 27 31 30)
    )
);

mergePatchPairs
(
);

and the probeDefinitions/constant is

surfaceElevationAnyName
{
    type waveGauge;

    pointDistribution lineDistribution;
    N 3;
    linestart (12.0 -0.15 0.05);
    lineend   (14.0 -0.15 0.05);
    add       ( 0 0.3 0);
    axis      y;
    stretch 1.;
}

the wave gauge is in the computational domain.

Thanks a lot!

[ngj]

Hallo,

Yes, my initial guess was correct: your wave gauges are outside of the computational domain.

Good luck,

Niels

[show]

Hi, Neils,

The problem was solved, Thank you for your information!

cheers,

shengchao