CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   interDyMFoam in parallel (http://www.cfd-online.com/Forums/openfoam-solving/118096-interdymfoam-parallel.html)

idefix May 21, 2013 02:34

interDyMFoam in parallel
 
Hello,

I get the following error massage by running interDyMFoam in parallel:
Quote:

user@pc53:~/OpenFOAM/user-2.1.1/run/surface_stripping/parameterstudie/HC3/2.1.1/turbulenzmodell-testen/dynamicMesh/01test> [0] #0 Foam::error::printStack(Foam::Ostream&) [1] #0 Foam::error::printStack(Foam::Ostream&)[2] #0 Foam::error::printStack(Foam::Ostream&)[3] #0 Foam::error::printStack(Foam::Ostream&) in "/home/user/OpenFOAM/OpenFOAM -2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #1 Foam::sigFpe::sigHandler(int) in "/home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[2] #1 Foam::sigFpe::sigHandler(int) in "/home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #1 Foam::sigFpe::sigHandler(int) in "/home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #1 Foam::sigFpe::sigHandler(int) in "/home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[2] #2 in "/home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #2 in "/home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #2 in "/home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #2 in "/lib64/libc.so.6"
[2] #3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double , Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Fo am::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/lib64/libc.so.6"
[1] #3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double , Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Fo am::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/lib64/libc.so.6"
[3] #3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double , Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Fo am::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/lib64/libc.so.6"
[0] #3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double , Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Fo am::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVo lume.so"
[2] #4 in "/home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[1] #4 in "/home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[3] #4 in "/home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #4 at MULES.C:0
[3] #5 at MULES.C:0
[2] #5 Foam::MULES::explicitSolve(Foam::GeometricField<do uble, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> cons t&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double)Foam::MULES::explicitSolve(Foam::GeometricF ield<double, Foam::fvPatchField, Foam::vo lMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) a t MULES.C:0
[0] #5 Foam::MULES::explicitSolve(Foam::GeometricField<do uble, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> cons t&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) at MULES.C:0
[1] #5 Foam::MULES::explicitSolve(Foam::GeometricField<do uble, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> cons t&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolum e.so"
[2] #6 in "/home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #6 in "/home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/
lib/libfiniteVolume.so"
[3] #6 in "/home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[1] #6
[2] in "/home/user/OpenFOAM/user-2.1.1/platforms/linux64GccDPOpt/bin/VOFinterDyMFoam"
[2] #7 __libc_start_main

[0] in "/home/user/OpenFOAM/user-2.1.1/platforms/linux64GccDPOpt/bin/VOFinterDyMFoam"
[0] #7 __libc_start_main in "/lib64/libc.so.6"
[2] #8 [1] in "/home/user/OpenFOAM/user-2.1.1/platforms/linux64GccDPOpt/bin/VOFinterDyMFoam"
[1] #7 __libc_start_main[3] in "/home/user/OpenFOAM/user-2.1.1/platforms/linux64GccDPOpt/bin/VOFinterDyMFoam"
[3] #7 __libc_start_main
in "/lib64/libc.so.6"
[0] #8 in "/lib64/libc.so.6"
[1] #8 in "/lib64/libc.so.6"
[3] #8 [2] at /usr/src/packages/BUILD/glibc-2.11.2/csu/../sysdeps/x86_64/elf/start.S:116
[pc48:30685] *** Process received signal ***
[pc48:30685] Signal: Floating point exception (8)
[pc48:30685] Signal code: (-6)
[pc48:30685] Failing at address: 0x216000077dd


[pc48:30685] [ 0] /lib64/libc.so.6(+0x32a60) [0x7fb0fade6a60]
[pc48:30685] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x7fb0fade69e5]
[pc48:30685] [ 2] /lib64/libc.so.6(+0x32a60) [0x7fb0fade6a60]
[pc48:30685] [ 3] /home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES7limiterINS_17geo metricOneFieldENS_9zeroFieldES3_EEvRN S_5FieldIdEERKT_RKNS_14GeometricFieldIdNS_12fvPatc hFieldENS_7volMeshEEERKNSA_IdNS_13fvsPatchFieldENS _11surfaceMeshEEESK_RKT0_RKT1_ddi+0x1078) [0x7fb0fd9c7998]
[pc48:30685] [ 4] /home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(+0x806b6e) [0x7fb0fd9b9b6e]
[pc48:30685] [ 5] /home/user/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES13explicitSolveER NS_14GeometricFieldIdNS_12fvPatchFiel dENS_7volMeshEEERKNS1_IdNS_13fvsPatchFieldENS_11su rfaceMeshEEERS8_dd+0x80) [0x7fb0fd9b9e10]
[pc48:30685] [ 6] VOFinterDyMFoam() [0x434288]
[pc48:30685] [ 7] /lib64/libc.so.6(__libc_start_main+0xfd) [0x7fb0fadd2b7d]
[pc48:30685] [ 8] VOFinterDyMFoam() [0x425049]
[pc48:30685] *** End of error message ***

--------------------------------------------------------------------------
mpirun noticed that process rank 2 with PID 30685 on node pc48 exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

The last time steps in my log-file are:
Quote:

Interface Courant Number mean: 2.18646334687e-05 max: 0.254532086074
Courant Number mean: 0.000239346723633 max: 1.65895015383
deltaT = 4.20235726211e-11
Time = 1.87186135528

Selected 173 cells for refinement out of 75979.
Refined from 75979 to 77190 cells.
Selected 0 split points out of a possible 8295.
Execution time for mesh.update() = 0.309999999998 s
time step continuity errors : sum local = 1.00039221839e-06, global = -3.483097512e-09, cumulative = -0.018811995198
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 4.08790127103e-06, No Iterations 7
time step continuity errors : sum local = 4.08950462107e-12, global = 9.21422024607e-13, cumulative = -0.0188119951971
MULES: Solving for alpha1
Phase-1 volume fraction = 0.416386179111 Min(alpha1) = -1.5267257622e-26 Max(alpha1) = 1.00000115258
MULES: Solving for alpha1
Phase-1 volume fraction = 0.416386178858 Min(alpha1) = -8.91381120585e-27 Max(alpha1) = 1.00000115258
MULES: Solving for alpha1
Phase-1 volume fraction = 0.416386178604 Min(alpha1) = -1.08596070544e-26 Max(alpha1) = 1.00000115258
GAMG: Solving for p_rgh, Initial residual = 0.244905115864, Final residual = 0.00121327413155, No Iterations 5
time step continuity errors : sum local = 1.32713444941e-07, global = 2.84927591875e-08, cumulative = -0.0188119667043
GAMGPCG: Solving for p_rgh, Initial residual = 3.76384153428e-08, Final residual = 8.68458829731e-10, No Iterations 1
time step continuity errors : sum local = 7.32584014989e-06, global = 1.0541322618e-06, cumulative = -0.018810912572
smoothSolver: Solving for epsilon, Initial residual = 0.996385458551, Final residual = 2.11894825636e-06, No Iterations 1000
bounding epsilon, min: -4.9480321014e+16 max: 5.0654485474e+21 average: 1.72580986552e+18
smoothSolver: Solving for k, Initial residual = 0.985758981367, Final residual = 0.00712526678361, No Iterations 1000
bounding k, min: -5.04700160787e-05 max: 2.20748593107e+15 average: 2.01192671257e+12
ExecutionTime = 153713.62 s ClockTime = 154490 s

Interface Courant Number mean: 1.41025652279e-05 max: 0.0848984222164
Courant Number mean: 0.000467736127991 max: 9.6991387301
deltaT = 2.16635588945e-12
--> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 1010
Increased the timePrecision from 12 to 13 to distinguish between timeNames at time 1.87186135528
Time = 1.871861355284

Selected 2 cells for refinement out of 77190.
Refined from 77190 to 77204 cells.
Selected 0 split points out of a possible 8296.
Execution time for mesh.update() = 0.25 s
time step continuity errors : sum local = 1.35646528805e-06, global = 5.46496372006e-08, cumulative = -0.0188108579224
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 4.98161826634e-06, No Iterations 6
time step continuity errors : sum local = 6.75739225988e-12, global = -6.96771494063e-13, cumulative = -0.0188108579231
MULES: Solving for alpha1
Phase-1 volume fraction = 0.416386178542 Min(alpha1) = -1.05025651157e-27 Max(alpha1) = 1.00000115261
MULES: Solving for alpha1
Phase-1 volume fraction = 0.416386178367 Min(alpha1) = -2.21912630123e-27 Max(alpha1) = 1.00000115264
MULES: Solving for alpha1
Phase-1 volume fraction = 0.416386178192 Min(alpha1) = -1.87440792842e-27 Max(alpha1) = 1.00000115266
--> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 1010
Increased the timePrecision from 13 to 14 to distinguish between timeNames at time 1.87186135528
GAMG: Solving for p_rgh, Initial residual = 0.0522052818922, Final residual = 0.000491766590874, No Iterations 3
time step continuity errors : sum local = 0.00817760277586, global = -0.00106629138053, cumulative = -0.0198771493036
GAMGPCG: Solving for p_rgh, Initial residual = 0.0157662082368, Final residual = 9.9904851778e-10, No Iterations 7
time step continuity errors : sum local = 1.06749920518e-09, global = 3.4918147483e-11, cumulative = -0.0198771492687
smoothSolver: Solving for epsilon, Initial residual = 0.022241455999, Final residual = 1.00845657085e-10, No Iterations 1
bounding epsilon, min: -4.26554509282e+20 max: 2.44885187134e+27 average: 6.35003618774e+22
smoothSolver: Solving for k, Initial residual = 1.13116151847e-06, Final residual = 2.87138632584e-07, No Iterations 1
bounding k, min: -325452149996 max: 9.33124200194e+19 average: 2.19327730332e+15
ExecutionTime = 153714.73 s ClockTime = 154491 s

Interface Courant Number mean: 0.000202964601196 max: 2.35952666146
Courant Number mean: 0.0136565407562 max: 249.050219909
deltaT = 4.34923504633e-15
Time = 1.8718613552835

Selected 0 cells for refinement out of 77204.
Selected 172 split points out of a possible 8296.
Unrefined from 77204 to 76000 cells.
Execution time for mesh.update() = 0.349999999977 s
time step continuity errors : sum local = 3.46465509528e-07, global = -2.39190192038e-09, cumulative = -0.0198771516606
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 3.59797605787e-06, No Iterations 7
time step continuity errors : sum local = 1.24657460846e-12, global = 8.14608391133e-14, cumulative = -0.0198771516605
--> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 1010
Increased the timePrecision from 14 to 15 to distinguish between timeNames at time 1.87186135528
MULES: Solving for alpha1
Phase-1 volume fraction = 0.416386173663 Min(alpha1) = -5.79771957911e-31 Max(alpha1) = 1.00000115266
--> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 1010
Increased the timePrecision from 15 to 16 to distinguish between timeNames at time 1.87186135528
MULES: Solving for alpha1
Phase-1 volume fraction = 0.416386152159 Min(alpha1) = -5.79771957904e-31 Max(alpha1) = 1.00000115266
MULES: Solving for alpha1
Phase-1 volume fraction = 0.416386130654 Min(alpha1) = -7.12738923027e-31 Max(alpha1) = 1.00000115266
GAMG: Solving for p_rgh, Initial residual = 0.94516627313, Final residual = 0.00895500854434, No Iterations 2
time step continuity errors : sum local = 9.62003610005e-07, global = -3.95529727442e-09, cumulative = -0.0198771556158
GAMGPCG: Solving for p_rgh, Initial residual = 0.00644022734333, Final residual = 4.24042829899e-09, No Iterations 7
time step continuity errors : sum local = 5.93317777085e-08, global = 4.51542234236e-10, cumulative = -0.0198771551643
smoothSolver: Solving for epsilon, Initial residual = 0.999996273117, Final residual = 2.73070503443e-10, No Iterations 1
bounding epsilon, min: -6.04982593122e+25 max: 4.90595707702e+33 average: 3.47291805011e+29
smoothSolver: Solving for k, Initial residual = 3.99587395154e-06, Final residual = 3.12890370849e-06, No Iterations 1000
ExecutionTime = 153716.96 s ClockTime = 154494 s

Interface Courant Number mean: 0.00942374410484 max: 1674.29348036
Courant Number mean: 0.170367047615 max: 4175.18098897
deltaT = 2.59765393424e-19
Time = 1.871861355283515

Selected 172 cells for refinement out of 76000.
Refined from 76000 to 77204 cells.
Selected 0 split points out of a possible 8296.
Execution time for mesh.update() = 0.309999999998 s
time step continuity errors : sum local = 1.39878280483e-07, global = 3.38689799327e-09, cumulative = -0.0198771517774
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 9.72804698248e-06, No Iterations 6
time step continuity errors : sum local = 1.36074248473e-12, global = 2.11250369829e-14, cumulative = -0.0198771517774
--> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 1010
Increased the timePrecision from 16 to 17 to distinguish between timeNames at time 1.87186135528
My grid has two cells in the third dimension, the other two dimensions have a lot of cells.

The grid is a rectangular channel. Bottom and top are walls.

My boundary conditions are:
timeVaryingMappedFixedValue for gas and liquid
inletOutlet for outlet
fixedValue for wall
symmetryPlane for frontAndBack

Has anybone an idea what is going wrong?

Thanks a lot
Idefix


All times are GMT -4. The time now is 19:41.