[ce73stargazer]

Hi there

I was running a two-phase fluid problem, however, it crashed shortly after it started and I have no idea why it happened, any suggestions here?

Below is the output messages

Courant Number mean: 0.00693316 max: 0.249238
Interface Courant Number mean: 0.00020949 max: 0.0426987
deltaT = 0.00083725
Time = 0.209107

PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 0.00120734, Final residual = 0.019386, No Iterations 1000
Phase-1 volume fraction = 1250.59 Min(alpha.water) = -199775 Max(alpha.water) = 224852
MULES: Correcting alpha.water
Phase-1 volume fraction = 1250.59 Min(alpha.water) = -63308.6 Max(alpha.water) = 79178.7
Relaxing time: 0.01 s
DILUPBiCG: Solving for Ux, Initial residual = 0.0597176, Final residual = 2.96317e-11, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.00743806, Final residual = 1.08024e-10, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 0.00285687, Final residual = 9.85841e-08, No Iterations 24
GAMG: Solving for p_rgh, Initial residual = 7.48022e-06, Final residual = 9.49453e-08, No Iterations 4
time step continuity errors : sum local = 2.1612e-07, global = 8.21126e-09, cumulative = 1.3899e-06
GAMG: Solving for p_rgh, Initial residual = 8.79021e-06, Final residual = 9.01058e-08, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 1.74997e-07, Final residual = 8.3287e-08, No Iterations 1
time step continuity errors : sum local = 1.89999e-07, global = 2.11415e-08, cumulative = 1.41104e-06
GAMG: Solving for p_rgh, Initial residual = 6.32173e-07, Final residual = 8.35427e-08, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 1.03803e-07, Final residual = 9.82178e-09, No Iterations 41
time step continuity errors : sum local = 2.23999e-08, global = 4.81646e-10, cumulative = 1.41152e-06
ExecutionTime = 164.97 s ClockTime = 165 s

Courant Number mean: 0.00694366 max: 0.248626
Interface Courant Number mean: 0.000200139 max: 0.0413379
deltaT = 0.000840938
Time = 0.209948

PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 0.00118213, Final residual = 0.019572, No Iterations 1000
Phase-1 volume fraction = 1305.51 Min(alpha.water) = -62596.1 Max(alpha.water) = 78556.6
MULES: Correcting alpha.water
Phase-1 volume fraction = 1305.51 Min(alpha.water) = -126116 Max(alpha.water) = 143771
Relaxing time: 0.01 s
DILUPBiCG: Solving for Ux, Initial residual = 0.0600228, Final residual = 1.23929e-10, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.00764295, Final residual = 8.54626e-11, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 0.00292367, Final residual = 3.5649e+92, No Iterations 1000
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#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 ??:?
#4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::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 ??:?
#5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8 main at ??:?
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10 ? at ??:?
Floating point exception (core dumped)

And here is the output for checkMesh

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 4.0
Exec : checkMesh
Date : Sep 02 2016
Time : 14:06:46
Host : "stargazer"
PID : 7811
Case : /home/zhifang/OpenFoamTests/of40/waveFlume
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
points: 61200
internal points: 0
faces: 120449
internal faces: 59251
cells: 29950
faces per cell: 6
boundary patches: 6
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 29950
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0

Checking topology...
Boundary definition OK.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology
front 29950 30600 ok (non-closed singly connected)
back 29950 30600 ok (non-closed singly connected)
inlet 50 102 ok (non-closed singly connected)
beach 50 102 ok (non-closed singly connected)
atmosphere 599 1200 ok (non-closed singly connected)
bottom 599 1200 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (0 -2 -0.015) (59.9 0.2 0.015)
Mesh has 2 geometric (non-empty/wedge) directions (1 1 0)
Mesh has 2 solution (non-empty) directions (1 1 0)
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (1.27986e-19 1.22547e-16 1.56297e-13) OK.
Max cell openness = 2.40253e-16 OK.
Max aspect ratio = 99.942 OK.
Minimum face area = 3.00174e-05. Maximum face area = 0.0100058. Face area magnitudes OK.
Min volume = 3.00174e-06. Max volume = 0.000300174. Total volume = 2.79128. Cell volumes OK.
Mesh non-orthogonality Max: 6.2735 average: 0.886639
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.0816914 OK.
Coupled point location match (average 0) OK.

Mesh OK.

End

I did the test with OpenFOAM4.0, however, i tried this test with other OF version (OF240 and FoamExtend 3.1), it crashed in similar ways. What I see is really crazy here is the alpha.water value, which I think should be a value between 0 to 1, I have no idea why it gets to such crazy number. In contrast, at the start of the simulation, it looks fine to me (shown below)

Starting time loop

Courant Number mean: 0 max: 0
Interface Courant Number mean: 0 max: 0
deltaT = 0.0011976
Time = 0.0011976

PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 0, Final residual = 0, No Iterations 0
Phase-1 volume fraction = 0.871242 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.871242 Min(alpha.water) = 0 Max(alpha.water) = 1
Relaxing time: 0.02 s
DILUPBiCG: Solving for Ux, Initial residual = 0.000847825, Final residual = 8.08268e-15, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1.15755e-06, Final residual = 1.76858e-17, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 9.39922e-08, No Iterations 19
GAMG: Solving for p_rgh, Initial residual = 4.84878e-11, Final residual = 4.84878e-11, No Iterations 0
time step continuity errors : sum local = 3.05765e-15, global = -3.37343e-16, cumulative = -3.37343e-16
GAMG: Solving for p_rgh, Initial residual = 1.55096e-08, Final residual = 1.55096e-08, No Iterations 0
GAMG: Solving for p_rgh, Initial residual = 1.55096e-08, Final residual = 1.55096e-08, No Iterations 0
time step continuity errors : sum local = 9.7804e-13, global = 3.90545e-14, cumulative = 3.87172e-14
GAMG: Solving for p_rgh, Initial residual = 1.55101e-08, Final residual = 1.55101e-08, No Iterations 0
GAMG: Solving for p_rgh, Initial residual = 1.55101e-08, Final residual = 1.10642e-09, No Iterations 1
time step continuity errors : sum local = 6.97712e-14, global = 2.67358e-14, cumulative = 6.5453e-14
ExecutionTime = 0.36 s ClockTime = 0 s

Courant Number mean: 7.99602e-08 max: 6.41771e-05
Interface Courant Number mean: 0 max: 0
deltaT = 0.00143023
Time = 0.00262784

PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 3.94481e-09, Final residual = 3.94481e-09, No Iterations 0
Phase-1 volume fraction = 0.871242 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.871242 Min(alpha.water) = 0 Max(alpha.water) = 1
Relaxing time: 0.01 s
DILUPBiCG: Solving for Ux, Initial residual = 0.0821968, Final residual = 3.37494e-13, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.0995908, Final residual = 8.26784e-12, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 0.00466376, Final residual = 7.75137e-08, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 7.49742e-08, Final residual = 7.49742e-08, No Iterations 0
time step continuity errors : sum local = 7.0052e-12, global = -3.54341e-12, cumulative = -3.47796e-12
GAMG: Solving for p_rgh, Initial residual = 1.98265e-07, Final residual = 2.43121e-08, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 2.42912e-08, Final residual = 2.42912e-08, No Iterations 0
time step continuity errors : sum local = 2.27165e-12, global = -1.0442e-12, cumulative = -4.52216e-12
GAMG: Solving for p_rgh, Initial residual = 2.42823e-08, Final residual = 2.42823e-08, No Iterations 0
GAMG: Solving for p_rgh, Initial residual = 2.42823e-08, Final residual = 9.56305e-09, No Iterations 1
time step continuity errors : sum local = 8.94311e-13, global = -4.50834e-13, cumulative = -4.97299e-12
ExecutionTime = 0.48 s ClockTime = 0 s

Courant Number mean: 1.95193e-07 max: 3.24846e-05
Interface Courant Number mean: 0 max: 0
deltaT = 0.00170148
Time = 0.00432932

PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 1.05581e-08, Final residual = 4.42952e-07, No Iterations 1000
Phase-1 volume fraction = 0.871243 Min(alpha.water) = 0 Max(alpha.water) = 1.02223
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.871243 Min(alpha.water) = 0 Max(alpha.water) = 1.02223
Relaxing time: 0 s
DILUPBiCG: Solving for Ux, Initial residual = 0.0632514, Final residual = 6.06682e-14, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.0673685, Final residual = 1.29071e-12, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 0.00626382, Final residual = 6.36146e-08, No Iterations 6
GAMG: Solving for p_rgh, Initial residual = 5.24492e-08, Final residual = 5.24492e-08, No Iterations 0
time step continuity errors : sum local = 8.47279e-12, global = -7.67813e-12, cumulative = -1.26511e-11
GAMG: Solving for p_rgh, Initial residual = 2.60459e-07, Final residual = 2.13011e-08, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 2.1271e-08, Final residual = 2.1271e-08, No Iterations 0
time step continuity errors : sum local = 3.44117e-12, global = -2.89041e-12, cumulative = -1.55415e-11
GAMG: Solving for p_rgh, Initial residual = 2.12575e-08, Final residual = 2.12575e-08, No Iterations 0
GAMG: Solving for p_rgh, Initial residual = 2.12575e-08, Final residual = 9.22819e-09, No Iterations 1
time step continuity errors : sum local = 1.49291e-12, global = -1.3672e-12, cumulative = -1.69087e-11
ExecutionTime = 1.4 s ClockTime = 1 s

Anyone can help? thanks in advance.