[hcen]

Hi everyone,

I am a newbie to OpenFoam, and am now trying to simulate a compressed air-water two phase pipe flow.

I've a mesh created from ICEM with only hex cells. The solver have been crashing and I have no clue what is leading to that. The model is pre-filled with water and has one compressed air inlet and two outlets. Below is the message from terminal. Could anyone take a look and give me a few advises? Thx!

Courant Number mean: 0.00344514823 max: 0.102423844
Max Ur Courant Number = 0.809298941
deltaT = 3.73662022e-06
Time = 0.125824067

PIMPLE: iteration 1
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.0454209699 Min(alpha.air) = 2.97973812e-316 Max(alpha.air) = 1
Constructing momentum equations
smoothSolver: Solving for e.air, Initial residual = 0.000370432209, Final residual = 4.56780784e-09, No Iterations 2
smoothSolver: Solving for e.water, Initial residual = 3.72096286e-07, Final residual = 5.83454128e-11, No Iterations 1
min T.air -51.5335211
min T.water 208.307104
GAMG: Solving for p_rgh, Initial residual = 0.000283639836, Final residual = 1.04823569e-10, No Iterations 2
PIMPLE: iteration 2
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.0454209586 Min(alpha.air) = -3.97078944e-17 Max(alpha.air) = 1
Constructing momentum equations
[9] #0 Foam::error:rintStack(Foam::Ostream&)[11] #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[9] #1 Foam::sigFpe::sigHandler(int) at ??:?
[11] #1 Foam::sigFpe::sigHandler(int) at ??:?
[9] #2 ? at ??:?
[11] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[9] #3 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[11] #3 ? in "/lib/x86_64-linux-gnu/libm.so.6"
[9] #4 pow in "/lib/x86_64-linux-gnu/libm.so.6"
[11] #4 pow in "/lib/x86_64-linux-gnu/libm.so.6"
[9] #5 Foam:ow(Foam::Field<double>&, Foam::UList<double> const&, double const&) in "/lib/x86_64-linux-gnu/libm.so.6"
[11] #5 Foam:ow(Foam::Field<double>&, Foam::UList<double> const&, double const&) at ??:?
[9] #6 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:ow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensioned<double> const&) at ??:?
[9] #7 Foam::dragModels::SchillerNaumann::CdRe() const at ??:?
[11] #6 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:ow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensioned<double> const&) at ??:?
[9] #8 Foam::dragModel::Ki() const at ??:?
[11] #7 Foam::dragModels::SchillerNaumann::CdRe() const at ??:?
[9] #9 Foam::dragModel::K() const at ??:?
[11] #8 Foam::dragModel::Ki() const at ??:?
[9] #10 Foam::BlendedInterfacialModel<Foam::dragModel>::K( ) const at ??:?
[9] #11 Foam::twoPhaseSystem::Kd() const at ??:?
[11] #9 Foam::dragModel::K() const at ??:?
[11] #10 Foam::BlendedInterfacialModel<Foam::dragModel>::K( ) const at ??:?
[11] #11 Foam::twoPhaseSystem::Kd() const at ??:?
[9] #12 ? at ??:?
[11] #12 at ??:?
[9] #13 __libc_start_main? in "/lib/x86_64-linux-gnu/libc.so.6"
[9] #14 ? at ??:?
[tande-Precision-Tower-7910:08931] *** Process received signal ***
[tande-Precision-Tower-7910:08931] Signal: Floating point exception (8)
[tande-Precision-Tower-7910:08931] Signal code: (-6)
[tande-Precision-Tower-7910:08931] Failing at address: 0x3e8000022e3
[tande-Precision-Tower-7910:08931] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x35250) [0x7f1e66adb250]
[tande-Precision-Tower-7910:08931] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x37) [0x7f1e66adb1c7]
[tande-Precision-Tower-7910:08931] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x35250) [0x7f1e66adb250]
[tande-Precision-Tower-7910:08931] [ 3] /lib/x86_64-linux-gnu/libm.so.6(+0x12f0c) [0x7f1e67099f0c]
[tande-Precision-Tower-7910:08931] [ 4] /lib/x86_64-linux-gnu/libm.so.6(pow+0x1c) [0x7f1e670ab44c]
[tande-Precision-Tower-7910:08931] [ 5] /home/tande/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam3powERNS_5FieldIdEERKNS_5UL istIdEERKd+0x38) [0x7f1e67f03348]
[tande-Precision-Tower-7910:08931] [ 6] /home/tande/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so(_ZN4Foam3powINS _12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14Geometr icFieldIdT_T0_EEEERKS7_RKNS_11dimensionedIdEE+0x28 e) [0x7f1e6befb2ae]
[tande-Precision-Tower-7910:08931] [ 7] /home/tande/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam10dragModels15SchillerNaumann4CdReEv+0x1d0) [0x7f1e6b0634a0]
[tande-Precision-Tower-7910:08931] [ 8] /home/tande/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam9dragModel2KiEv+0x132) [0x7f1e6b02a1f2]
[tande-Precision-Tower-7910:08931] [ 9] /home/tande/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam9dragModel1KEv+0x2e) [0x7f1e6b029ffe]
[tande-Precision-Tower-7910:08931] [10] /home/tande/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam23Blende dInterfacialModelINS_9dragModelEE1KEv+0x4a3) [0x7f1e6b450453]
[tande-Precision-Tower-7910:08931] [11] /home/tande/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam14twoPha seSystem2KdEv+0x32) [0x7f1e6b436042]
[tande-Precision-Tower-7910:08931] [12] twoPhaseEulerFoam() [0x43510b]
[tande-Precision-Tower-7910:08931] [13] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0) [0x7f1e66ac6ac0]
[tande-Precision-Tower-7910:08931] [14] twoPhaseEulerFoam() [0x4496a9]
[tande-Precision-Tower-7910:08931] *** End of error message ***
at ??:?
[11] #13 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[11] #14 ? at ??:?
[tande-Precision-Tower-7910:08933] *** Process received signal ***
[tande-Precision-Tower-7910:08933] Signal: Floating point exception (8)
[tande-Precision-Tower-7910:08933] Signal code: (-6)
[tande-Precision-Tower-7910:08933] Failing at address: 0x3e8000022e5
[tande-Precision-Tower-7910:08933] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x35250) [0x7fcddfb32250]
[tande-Precision-Tower-7910:08933] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x37) [0x7fcddfb321c7]
[tande-Precision-Tower-7910:08933] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x35250) [0x7fcddfb32250]
[tande-Precision-Tower-7910:08933] [ 3] /lib/x86_64-linux-gnu/libm.so.6(+0x12f0c) [0x7fcde00f0f0c]
[tande-Precision-Tower-7910:08933] [ 4] /lib/x86_64-linux-gnu/libm.so.6(pow+0x1c) [0x7fcde010244c]
[tande-Precision-Tower-7910:08933] [ 5] /home/tande/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam3powERNS_5FieldIdEERKNS_5UL istIdEERKd+0x38) [0x7fcde0f5a348]
[tande-Precision-Tower-7910:08933] [ 6] /home/tande/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so(_ZN4Foam3powINS _12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14Geometr icFieldIdT_T0_EEEERKS7_RKNS_11dimensionedIdEE+0x22 e) [0x7fcde4f5224e]
[tande-Precision-Tower-7910:08933] [ 7] /home/tande/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam10dragModels15SchillerNaumann4CdReEv+0x1d0) [0x7fcde40ba4a0]
[tande-Precision-Tower-7910:08933] [ 8] /home/tande/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam9dragModel2KiEv+0x132) [0x7fcde40811f2]
[tande-Precision-Tower-7910:08933] [ 9] /home/tande/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam9dragModel1KEv+0x2e) [0x7fcde4080ffe]
[tande-Precision-Tower-7910:08933] [10] /home/tande/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam23Blende dInterfacialModelINS_9dragModelEE1KEv+0x4a3) [0x7fcde44a7453]
[tande-Precision-Tower-7910:08933] [11] /home/tande/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam14twoPha seSystem2KdEv+0x32) [0x7fcde448d042]
[tande-Precision-Tower-7910:08933] [12] twoPhaseEulerFoam() [0x43510b]
[tande-Precision-Tower-7910:08933] [13] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0) [0x7fcddfb1dac0]
[tande-Precision-Tower-7910:08933] [14] twoPhaseEulerFoam() [0x4496a9]
[tande-Precision-Tower-7910:08933] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 9 with PID 8931 on node tande-Precision-Tower-7910 exited on signal 8 (Floating point exception).

Everything seemed fine to me. The residuals are small, reached with only 1 or 2 iterations. I just didn't understand what was causing the problem. Any suggestion would be appreciated!!!

Best,
HC