[jiadongw]

Hello Foamers,
I’m working on the waveExampleSolitary tutorial by OpenFoam v1806.
But the tutorial crash when I change the mesh. The new mesh is conducted by ICEM.
The computational domain is attached.

The error message is:

Courant Number mean: 0.00142094 max: 0.666531
Interface Courant Number mean: 0.000296555 max: 0.0831518
deltaT = 0.00152778
Time = 3.37456

PIMPLE: iteration 1
Updating Boussinesq wave model for patch INLET
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.510157  Min(alpha.water) = -0.000917064  Max(alpha.water) = 1.13209
Updating Boussinesq wave model for patch INLET
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.510156  Min(alpha.water) = -0.00167893  Max(alpha.water) = 1.13282
Updating Boussinesq wave model for patch INLET
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.510154  Min(alpha.water) = -0.00231185  Max(alpha.water) = 1.13354
Updating Boussinesq wave model for patch INLET
Updating shallowWaterAbsorption wave model for patch OUTLET
DICPCG:  Solving for p_rgh, Initial residual = 0.00680415, Final residual = 0.00046163, No Iterations 14
time step continuity errors : sum local = 1.32835e-06, global = -2.59204e-08, cumulative = -3.63653e-06
[0] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[0] #1  Foam::sigFpe::sigHandler(int) at ??:?
[0] #2  ? in /lib/x86_64-linux-gnu/libc.so.6
[0] #3  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
[0] #4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
[0] #5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[0] #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[0] #7  Foam::fvMatrix<double>::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?
[0] #8  Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?
[0] #9  ? in ~/OpenFOAM/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/bin/interFoam
[0] #10  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
[0] #11  ? in ~/OpenFOAM/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/bin/interFoam
[ryou-virtual-machine:07043] *** Process received signal ***
[ryou-virtual-machine:07043] Signal: Floating point exception (8)
[ryou-virtual-machine:07043] Signal code:  (-6)
[ryou-virtual-machine:07043] Failing at address: 0x3e800001b83
[ryou-virtual-machine:07043] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7ff1dceddf20]
[ryou-virtual-machine:07043] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xc7)[0x7ff1dcedde97]
[ryou-virtual-machine:07043] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7ff1dceddf20]
[ryou-virtual-machine:07043] [ 3] /home/ryou/OpenFOAM/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5scaleERNS_5FieldIdEES3_RKNS_9lduMatrixERKNS_10FieldFieldIS1_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEERKS2_h+0xaf)[0x7ff1de318d7f]
[ryou-virtual-machine:07043] [ 4] /home/ryou/OpenFOAM/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0xab9)[0x7ff1de31cbf9]
[ryou-virtual-machine:07043] [ 5] /home/ryou/OpenFOAM/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x6f5)[0x7ff1de31ee65]
[ryou-virtual-machine:07043] [ 6] /home/ryou/OpenFOAM/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x187)[0x7ff1e108cf07]
[ryou-virtual-machine:07043] [ 7] /home/ryou/OpenFOAM/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE24solveSegregatedOrCoupledERKNS_10dictionaryE+0x3a8)[0x7ff1e0c28418]
[ryou-virtual-machine:07043] [ 8] /home/ryou/OpenFOAM/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZNK4Foam6fvMesh5solveERNS_8fvMatrixIdEERKNS_10dictionaryE+0x23)[0x7ff1e0bcb6a3]
[ryou-virtual-machine:07043] [ 9] interFoam(+0x4c36a)[0x55a72125736a]
[ryou-virtual-machine:07043] [10] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7ff1dcec0b97]
[ryou-virtual-machine:07043] [11] interFoam(+0x560ea)[0x55a7212610ea]
[ryou-virtual-machine:07043] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 0 on node ryou-virtual-machine exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

I guess this problem is due to wrong definition of setFieldsDict, because the water.alpha of initial inlet seems to be weird, as attached figure shown.
The setFieldsDict is defined as using the new mesh:

defaultFieldValues
(
    volScalarFieldValue alpha.water 0
);

regions
(
    boxToCell
    {
        box (0 0 -0.1513274595) (15 1 0.1513274595);
        fieldValues
        (
            volScalarFieldValue alpha.water 1
        );
    }
);

Any suggestions are appreciated.
Thank you.

[shixiangyu]

i think u can change CFL <0.5 ， could tell me how to plot the first pic. did u uesd paraview?