Multiple floating objects

CKH · September 12, 2011, 10:27

Dear Foamers,

My eventual goal is to simulate multiple ships floating side-by-side in waves.

As a first step, I have modified the floatingObject example under the interDyMFoam sub-directory to model 2 boxes. I have tried to refine my mesh and adjust the time steps, but my two attempts seemed to stop at the same point in time regardless of the mesh refinement.

The error message that I am getting is as follows:

Code:

MULES: Solving for alpha1
Liquid phase volume fraction = 0.503994516429  Min(alpha1) = -2.545080202e-15  Max(alpha1) = 1.07101310392e+292
[0] #0  Foam::error::printStack(Foam::Ostream&) in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #1  Foam::sigFpe::sigHandler(int) in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #2   in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3  Foam::PhiScheme<double, Foam::interfaceCompressionLimiter>::limiter(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libinterfaceProperties.so"
[0] #4  Foam::limitedSurfaceInterpolationScheme<double>::weights(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #5  Foam::surfaceInterpolationScheme<double>::interpolate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #6  Foam::fv::gaussConvectionScheme<double>::interpolate(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMI am trying to simulate esh> const&) const in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #7  Foam::fv::gaussConvectionScheme<double>::flux(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #8  
[0]  in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[0] #9  
[0]  in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[0] #10  
[0]  in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[0] #11  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #12  
[0]  in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[ceeckh:01715] *** Process received signal ***
[ceeckh:01715] Signal: Floating point exception (8)
[ceeckh:01715] Signal code:  (-6)
[ceeckh:01715] Failing at address: 0x3e8000006b3
[ceeckh:01715] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x33d80) [0x7ff0ebd8bd80]
[ceeckh:01715] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7ff0ebd8bd05]
[ceeckh:01715] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x33d80) [0x7ff0ebd8bd80]
[ceeckh:01715] [ 3] /home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libinterfaceProperties.so(_ZNK4Foam9PhiSchemeIdNS_27interfaceCompressionLimiterEE7limiterERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEE+0x562) [0x7ff0ebb47ca2]
[ceeckh:01715] [ 4] /home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZNK4Foam33limitedSurfaceInterpolationSchemeIdE7weightsERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEE+0x25) [0x7ff0ee778195]
[ceeckh:01715] [ 5] /home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZNK4Foam26surfaceInterpolationSchemeIdE11interpolateERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEE+0x35) [0x7ff0ee76d965]
[ceeckh:01715] [ 6] /home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZNK4Foam2fv21gaussConvectionSchemeIdE11interpolateERKNS_14GeometricFieldIdNS_13fvsPatchFieldENS_11surfaceMeshEEERKNS3_IdNS_12fvPatchFieldENS_7volMeshEEE+0x3c) [0x7ff0ee76580c]
[ceeckh:01715] [ 7] /home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZNK4Foam2fv21gaussConvectionSchemeIdE4fluxERKNS_14GeometricFieldIdNS_13fvsPatchFieldENS_11surfaceMeshEEERKNS3_IdNS_12fvPatchFieldENS_7volMeshEEE+0x15) [0x7ff0ee76eca5]
[ceeckh:01715] [ 8] interDyMFoam() [0x47a445]
[ceeckh:01715] [ 9] interDyMFoam() [0x47a56d]
[ceeckh:01715] [10] interDyMFoam() [0x42cf85]
[ceeckh:01715] [11] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xff) [0x7ff0ebd76eff]
[ceeckh:01715] [12] interDyMFoam() [0x4243b9]
[ceeckh:01715] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 1715 on node ceeckh exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

I have also attached my blockMesh file and the dictionary files under the system directory. Really appreciate if anyone could help me with this problem.

Thank you.

kiehian

jhoepken · September 13, 2011, 04:26

Hi,

From what I can see in the very short log file excerpt, I am guessing that your time step is way to large. Your alpha is exploding. Do you use the CFL based time step control?

P.S.: A longer excerpt of the logfile would have been helpful.

CKH · September 13, 2011, 05:15

Hi Jens,

Thanks for helping me out with this.

As this is based on the example case of dam break with a floating box (floatingObject), I am not too sure what is a good value for U in the Courant number.

However, the time step I have used is much smaller than if I assume the same U value as the example and recalculate the required time step for the mesh I am using. (I used order of 1E-5 vs a CFL-based order of 1E-3).

Below is a longer excerpt of the log file, which includes the last few stable time steps at the set Courant number of 0.25:

Code:

Interface Courant Number mean: 0.000325787599653 max: 0.179187841883
Courant Number mean: 0.00423727237639 max: 0.25531975441
deltaT = 4.87923599541e-06
Time = 0.0180437

Centre of mass: (-0.350072782873 0.0651312159514 0.199485635643)
Linear velocity: (-0.143284476665 78.636476626 2.34574296499)
Angular velocity: (423.010487298 1.35901090149 -1.96111987984)
Centre of mass: (0.350087441965 0.0684517136853 -0.202485695229)
Linear velocity: (0.138516646667 84.4349192533 -3.20615213875)
Angular velocity: (-439.76316529 0.902495680635 -0.661229840877)
GAMG:  Solving for cellDisplacementx, Initial residual = 0.00445716240619, Final residual = 6.38886137841e-06, No Iterations 3
GAMG:  Solving for cellDisplacementy, Initial residual = 0.00450350518389, Final residual = 6.4389446049e-06, No Iterations 3
GAMG:  Solving for cellDisplacementz, Initial residual = 0.00298902386253, Final residual = 6.17529067946e-06, No Iterations 3
Execution time for mesh.update() = 2.56 s
time step continuity errors : sum local = 1.45043905361e-12, global = -2.69333372684e-14, cumulative = 2.26640809904e-06
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 2.62647669918e-06, No Iterations 9
time step continuity errors : sum local = 3.81386055149e-18, global = -8.34194182978e-19, cumulative = 2.26640809904e-06
MULES: Solving for alpha1
Liquid phase volume fraction = 0.504072648962  Min(alpha1) = -2.94594100549e-15  Max(alpha1) = 0.99999999999
MULES: Solving for alpha1
Liquid phase volume fraction = 0.50404131509  Min(alpha1) = -2.90139095429e-15  Max(alpha1) = 0.99999999999
MULES: Solving for alpha1
Liquid phase volume fraction = 0.504009991105  Min(alpha1) = -2.85744731133e-15  Max(alpha1) = 0.99999999999
GAMG:  Solving for p_rgh, Initial residual = 0.00381675881406, Final residual = 3.0326763017e-05, No Iterations 3
time step continuity errors : sum local = 6.12634251692e-09, global = 3.40911217402e-10, cumulative = 2.26674901026e-06
GAMGPCG:  Solving for p_rgh, Initial residual = 0.00323918316439, Final residual = 4.77938073206e-09, No Iterations 5
time step continuity errors : sum local = 9.36394656984e-13, global = 2.95758071703e-15, cumulative = 2.26674901321e-06
smoothSolver:  Solving for epsilon, Initial residual = 0.00105894903597, Final residual = 2.53065650016e-07, No Iterations 3
smoothSolver:  Solving for k, Initial residual = 0.00139861690848, Final residual = 2.97634249492e-07, No Iterations 3
ExecutionTime = 16249.63 s  ClockTime = 16414 s

Interface Courant Number mean: 0.000309976335649 max: 0.170755924658
Courant Number mean: 0.00408901570098 max: 0.252552720149
deltaT = 4.80567967387e-06
Time = 0.0180485

Centre of mass: (-0.350073471848 0.0655091405472 0.199496919991)
Linear velocity: (-0.143450542713 78.6455253187 2.35034888503)
Angular velocity: (423.024493198 1.35852694248 -1.97080261532)
Centre of mass: (0.350088108055 0.0688575278015 -0.202501124668)
Linear velocity: (0.13869169085 84.4542457509 -3.21511392728)
Angular velocity: (-439.832644812 0.906936640071 -0.650025314925)
GAMG:  Solving for cellDisplacementx, Initial residual = 0.0043653668749, Final residual = 6.29435478485e-06, No Iterations 3
GAMG:  Solving for cellDisplacementy, Initial residual = 0.00441179954511, Final residual = 6.41083269994e-06, No Iterations 3
GAMG:  Solving for cellDisplacementz, Initial residual = 0.00293315398511, Final residual = 6.22896204948e-06, No Iterations 3
Execution time for mesh.update() = 2.56 s
time step continuity errors : sum local = 9.22278154598e-13, global = 2.91299410271e-15, cumulative = 2.26674901613e-06
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 4.95868725448e-06, No Iterations 8
time step continuity errors : sum local = 4.57477513469e-18, global = -6.5832569033e-19, cumulative = 2.26674901613e-06
MULES: Solving for alpha1
Liquid phase volume fraction = 0.504076048219  Min(alpha1) = -2.81480134182e-15  Max(alpha1) = 0.99999999999
MULES: Solving for alpha1
Liquid phase volume fraction = 0.504043012626  Min(alpha1) = -1.09755604631e-12  Max(alpha1) = 0.99999999999
MULES: Solving for alpha1
Liquid phase volume fraction = 0.504009986925  Min(alpha1) = -2.73122412697e-15  Max(alpha1) = 0.99999999999
GAMG:  Solving for p_rgh, Initial residual = 0.00384916493515, Final residual = 3.37872200267e-05, No Iterations 3
time step continuity errors : sum local = 6.49569788618e-09, global = 2.73567741513e-10, cumulative = 2.26702258387e-06
GAMGPCG:  Solving for p_rgh, Initial residual = 0.00335223305572, Final residual = 5.54940824135e-09, No Iterations 5
time step continuity errors : sum local = 1.05775382703e-12, global = 1.24722265671e-15, cumulative = 2.26702258511e-06
smoothSolver:  Solving for epsilon, Initial residual = 0.001023876895, Final residual = 2.66025571663e-07, No Iterations 3
smoothSolver:  Solving for k, Initial residual = 0.00139380629685, Final residual = 3.12120916037e-07, No Iterations 3
ExecutionTime = 16263.06 s  ClockTime = 16428 s

Interface Courant Number mean: 0.00032901121076 max: 0.252398496212
Courant Number mean: 0.00422058905311 max: 0.31475605604
deltaT = 3.8060983017e-06
Time = 0.0180523

Centre of mass: (-0.350074018154 0.065808489468 0.199505874105)
Linear velocity: (-0.143592936418 78.6536849538 2.35453712191)
Angular velocity: (423.038897196 1.35814474038 -1.97859681296)
Centre of mass: (0.35008863626 0.0691790053441 -0.202513378632)
Linear velocity: (0.138846405826 84.4713104606 -3.22307091531)
Angular velocity: (-439.895920779 0.910585338043 -0.640608809966)
GAMG:  Solving for cellDisplacementx, Initial residual = 0.00344389380666, Final residual = 5.50298289881e-06, No Iterations 3
GAMG:  Solving for cellDisplacementy, Initial residual = 0.00348180695973, Final residual = 5.56835601745e-06, No Iterations 3
GAMG:  Solving for cellDisplacementz, Initial residual = 0.00232298891502, Final residual = 5.92286399503e-06, No Iterations 3
Execution time for mesh.update() = 2.57 s
time step continuity errors : sum local = 8.37741028636e-13, global = 9.87800259381e-16, cumulative = 2.2670225861e-06
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 6.76592070776e-06, No Iterations 8
time step continuity errors : sum local = 5.6685560109e-18, global = 5.56846492987e-19, cumulative = 2.2670225861e-06
MULES: Solving for alpha1
Liquid phase volume fraction = 0.504054939957  Min(alpha1) = -2.69876398288e-15  Max(alpha1) = 0.99999999999
MULES: Solving for alpha1
Liquid phase volume fraction = 0.504032457732  Min(alpha1) = -2.66665258023e-15  Max(alpha1) = 0.99999999999
MULES: Solving for alpha1
Liquid phase volume fraction = 0.504009984315  Min(alpha1) = -2.63488693047e-15  Max(alpha1) = 0.99999999999
GAMG:  Solving for p_rgh, Initial residual = 0.00993855133339, Final residual = 3.87047122965e-05, No Iterations 3
time step continuity errors : sum local = 4.78330720513e-09, global = 3.09696673749e-10, cumulative = 2.26733228278e-06
GAMGPCG:  Solving for p_rgh, Initial residual = 0.00377928621112, Final residual = 5.10720921935e-09, No Iterations 5
time step continuity errors : sum local = 6.28297520456e-13, global = -2.27077377704e-14, cumulative = 2.26733226007e-06
smoothSolver:  Solving for epsilon, Initial residual = 0.000799159002912, Final residual = 9.04288901952e-08, No Iterations 3
smoothSolver:  Solving for k, Initial residual = 0.00106847871696, Final residual = 1.08432620282e-07, No Iterations 3
ExecutionTime = 16276.51 s  ClockTime = 16441 s

Interface Courant Number mean: 0.000229833000784 max: 0.302370957161
Courant Number mean: 0.00309451410083 max: 0.302370957161
deltaT = 3.13814065273e-06
Time = 0.0180554

Centre of mass: (-0.350074468952 0.066055327944 0.199513269145)
Linear velocity: (-0.143713194496 78.6591757349 2.35756092727)
Angular velocity: (423.057900844 1.359646774 -1.98905854207)
Centre of mass: (0.350089072193 0.0694441117524 -0.202523504097)
Linear velocity: (0.1389756274 84.4833121253 -3.2287992894)
Angular velocity: (-439.94364224 0.913700977288 -0.632466847919)
GAMG:  Solving for cellDisplacementx, Initial residual = 0.00282989255253, Final residual = 4.59766250124e-06, No Iterations 3
GAMG:  Solving for cellDisplacementy, Initial residual = 0.00286157907736, Final residual = 4.60982988026e-06, No Iterations 3
GAMG:  Solving for cellDisplacementz, Initial residual = 0.00191149708391, Final residual = 4.53120679811e-06, No Iterations 3
Execution time for mesh.update() = 2.57 s
time step continuity errors : sum local = 5.18033387764e-13, global = -1.87226050676e-14, cumulative = 2.26733224135e-06
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 8.27971214728e-06, No Iterations 7
time step continuity errors : sum local = 4.28989587445e-18, global = 2.57836271406e-19, cumulative = 2.26733224135e-06
MULES: Solving for alpha1
Liquid phase volume fraction = 0.504056385492  Min(alpha1) = -0.00133561642123  Max(alpha1) = 0.999999999989
MULES: Solving for alpha1
Liquid phase volume fraction = 0.5040331605  Min(alpha1) = -0.00018638629599  Max(alpha1) = 0.999999999989
MULES: Solving for alpha1
Liquid phase volume fraction = 0.504009952958  Min(alpha1) = -2.55791225148e-15  Max(alpha1) = 0.999999999989
GAMG:  Solving for p_rgh, Initial residual = 0.0103214760728, Final residual = 3.67270530086e-05, No Iterations 3
time step continuity errors : sum local = 3.12773541036e-09, global = 1.49081191687e-10, cumulative = 2.26748132254e-06
GAMGPCG:  Solving for p_rgh, Initial residual = 0.00379524481323, Final residual = 7.96472187558e-09, No Iterations 5
time step continuity errors : sum local = 7.1170746306e-13, global = -3.17949346268e-14, cumulative = 2.26748129074e-06
smoothSolver:  Solving for epsilon, Initial residual = 0.000692502328127, Final residual = 3.30601940636e-07, No Iterations 3
smoothSolver:  Solving for k, Initial residual = 0.00096642910504, Final residual = 7.24130195228e-07, No Iterations 3
ExecutionTime = 16289.61 s  ClockTime = 16455 s

Interface Courant Number mean: 0.000237780909322 max: 0.630649717903
Courant Number mean: 0.00286284794595 max: 1.60259184612
deltaT = 4.89419863457e-07
Time = 0.0180559

Centre of mass: (-0.350074539319 0.0660938260995 0.1995144235)
Linear velocity: (-0.143783194091 78.6609425431 2.35873797312)
Angular velocity: (423.069269931 1.3619925901 -1.99633914247)
Centre of mass: (0.35008914024 0.0694854617637 -0.202525085422)
Linear velocity: (0.139047057032 84.4882663814 -3.23126008249)
Angular velocity: (-439.965971177 0.915393153408 -0.627434148581)
GAMG:  Solving for cellDisplacementx, Initial residual = 0.000462900041698, Final residual = 5.84277791375e-06, No Iterations 2
GAMG:  Solving for cellDisplacementy, Initial residual = 0.000465965617342, Final residual = 5.77636772713e-06, No Iterations 2
GAMG:  Solving for cellDisplacementz, Initial residual = 0.000328786048385, Final residual = 6.52068720835e-06, No Iterations 2
Execution time for mesh.update() = 2.34 s
time step continuity errors : sum local = 1.10996862176e-13, global = -4.95869179953e-15, cumulative = 2.26748128578e-06
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 9.11419305921e-06, No Iterations 7
time step continuity errors : sum local = 1.01170753357e-18, global = 5.83478574024e-21, cumulative = 2.26748128578e-06
MULES: Solving for alpha1
Liquid phase volume fraction = 0.504004938521  Min(alpha1) = -2.55383589479e-15  Max(alpha1) = 0.999999999989
MULES: Solving for alpha1
Liquid phase volume fraction = 0.504007410302  Min(alpha1) = -2.54976550153e-15  Max(alpha1) = 0.999999999989
MULES: Solving for alpha1
Liquid phase volume fraction = 0.504009919507  Min(alpha1) = -2.5457010644e-15  Max(alpha1) = 0.999999999989
GAMG:  Solving for p_rgh, Initial residual = 0.180956363452, Final residual = 0.0012560699121, No Iterations 2
time step continuity errors : sum local = 3.52412976029e-09, global = 9.27564793752e-11, cumulative = 2.26757404226e-06
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0165923879666, Final residual = 1.61534301916e-09, No Iterations 6
time step continuity errors : sum local = 7.32155643598e-15, global = -2.92801907078e-16, cumulative = 2.26757404197e-06
smoothSolver:  Solving for epsilon, Initial residual = 0.000187859833342, Final residual = 1.05268296273e-07, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.000309447981548, Final residual = 4.17949801455e-07, No Iterations 2
ExecutionTime = 16302.76 s  ClockTime = 16468 s

Interface Courant Number mean: 6.32874578052e-05 max: 1.6630627622
Courant Number mean: 0.000736934887108 max: 1.6630627622
deltaT = 7.3567436812e-08
Time = 0.018056

Centre of mass: (-0.350074549897 0.0660996129943 0.199514597035)
Linear velocity: (-0.143783847451 78.6608037959 2.35882283859)
Angular velocity: (423.08140202 1.36346322286 -1.99905764201)
Centre of mass: (0.35008915047 0.0694916773827 -0.202525323155)
Linear velocity: (0.139063042666 84.4876274801 -3.23105709637)
Angular velocity: (-439.967177522 0.916019853021 -0.625638103411)
GAMG:  Solving for cellDisplacementx, Initial residual = 8.3298614334e-05, Final residual = 7.76464457012e-06, No Iterations 2
GAMG:  Solving for cellDisplacementy, Initial residual = 8.36912304455e-05, Final residual = 7.54171790949e-06, No Iterations 2
GAMG:  Solving for cellDisplacementz, Initial residual = 6.72628065284e-05, Final residual = 7.15960384245e-06, No Iterations 2
Execution time for mesh.update() = 2.34 s
time step continuity errors : sum local = 1.10054377168e-15, global = -4.40126768938e-17, cumulative = 2.26757404193e-06
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 6.55950473684e-06, No Iterations 7
time step continuity errors : sum local = 7.48728609849e-21, global = 1.02060308413e-21, cumulative = 2.26757404193e-06
--> FOAM Warning :
    From function Time::operator++()
    in file db/Time/Time.C at line 937
    Increased the timePrecision from 6 to 7 to distinguish between timeNames at time 0.0180559336059
MULES: Solving for alpha1
Liquid phase volume fraction = 0.503994516429  Min(alpha1) = -2.545080202e-15  Max(alpha1) = 1.07101310392e+292
[0] #0  Foam::error::printStack(Foam::Ostream&) in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #1  Foam::sigFpe::sigHandler(int) in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #2   in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3  Foam::PhiScheme<double, Foam::interfaceCompressionLimiter>::limiter(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libinterfaceProperties.so"
[0] #4  Foam::limitedSurfaceInterpolationScheme<double>::weights(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #5  Foam::surfaceInterpolationScheme<double>::interpolate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #6  Foam::fv::gaussConvectionScheme<double>::interpolate(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #7  Foam::fv::gaussConvectionScheme<double>::flux(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #8 
[0]  in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[0] #9 
[0]  in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[0] #10 
[0]  in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[0] #11  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #12 
[0]  in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[ceeckh:01715] *** Process received signal ***
[ceeckh:01715] Signal: Floating point exception (8)
[ceeckh:01715] Signal code:  (-6)
[ceeckh:01715] Failing at address: 0x3e8000006b3
[ceeckh:01715] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x33d80) [0x7ff0ebd8bd80]
[ceeckh:01715] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7ff0ebd8bd05]
[ceeckh:01715] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x33d80) [0x7ff0ebd8bd80]
[ceeckh:01715] [ 3] /home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libinterfaceProperties.so(_ZNK4Foam9PhiSchemeIdNS_27interfaceCompressionLimiterEE7limiterERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEE+0x562) [0x7ff0ebb47ca2]
[ceeckh:01715] [ 4] /home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZNK4Foam33limitedSurfaceInterpolationSchemeIdE7weightsERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEE+0x25) [0x7ff0ee778195]
[ceeckh:01715] [ 5] /home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZNK4Foam26surfaceInterpolationSchemeIdE11interpolateERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEE+0x35) [0x7ff0ee76d965]
[ceeckh:01715] [ 6] /home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZNK4Foam2fv21gaussConvectionSchemeIdE11interpolateERKNS_14GeometricFieldIdNS_13fvsPatchFieldENS_11surfaceMeshEEERKNS3_IdNS_12fvPatchFieldENS_7volMeshEEE+0x3c) [0x7ff0ee76580c]
[ceeckh:01715] [ 7] /home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZNK4Foam2fv21gaussConvectionSchemeIdE4fluxERKNS_14GeometricFieldIdNS_13fvsPatchFieldENS_11surfaceMeshEEERKNS3_IdNS_12fvPatchFieldENS_7volMeshEEE+0x15) [0x7ff0ee76eca5]
[ceeckh:01715] [ 8] interDyMFoam() [0x47a445]
[ceeckh:01715] [ 9] interDyMFoam() [0x47a56d]
[ceeckh:01715] [10] interDyMFoam() [0x42cf85]
[ceeckh:01715] [11] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xff) [0x7ff0ebd76eff]
[ceeckh:01715] [12] interDyMFoam() [0x4243b9]
[ceeckh:01715] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 1715 on node ceeckh exited on signal 8 (Floating point exception).

Thanks once again.

elvis · September 13, 2011, 05:58

Hi,

are you trying s.th. like http://powerlab.fsb.hr/ped/kturbo/Op...F/twoBoats.mpg ? I believe there is a extend http://openfoam-extend.git.sourcefor...enFOAM-1.6-ext =>tutorials/basic/sixDOFSolver/sixDOFmotion with bodyA & bodyB
but I am not sure

CKH · September 13, 2011, 06:04

Hi Elvis,

Thanks for pointing out the two boats case. I have looked at the extend version and actually, my blockMeshDict file is based on the case. However, instead of subjecting the two boxes to a steady flow, what I am plan to do is to subject them to waves. So as a first step, I tried to modify the floatingObject case to see if it works for two objects and if it does, my next step is to implement some form of wave input.

jhoepken · September 13, 2011, 08:11

About your crashing simulation, you might want to check the following:

Gravity pointing into the correct direction?
Field (alpha) initialised correctly?
Are there bad cells in your mesh?

CKH · September 14, 2011, 07:14

Hi Jens,

Thanks once again for your pointers. I've looked into my input files and found that my gravity was orientated along the wrong axis wrt the boats/alpha. However, after fixing that problem and further refinement of my mesh, the simulation still crashes somewhere at the start (around 0.1s). checkMesh says ok and visualising in paraFoam seems fine too.

Attached is a jpg of the alpha1 field at 0s. I'm now trying to re-run with a finer time-step and also with turbulence turned off to see if I can get the simulation to run properly.

However, if you have further insights, please let me know!

Thanks!

kiehian

jhoepken · September 14, 2011, 17:30

Again, a new log file would help a lot

CKH · September 15, 2011, 00:20

Apologies, I was thinking of running the simulation with smaller time step up and running first. However, even with smaller time steps, both the kEpsilon and laminar simulation still do not work...

Below is an extract of the log file from the turbulent run, after the Courant number started exploding (the full log is too long, so I've put it in the attachments). Also attached are the files from my system / 0 and constant directories.

Code:

Interface Courant Number mean: 0.000216627976107 max: 0.373212770233
Courant Number mean: 0.0054433824093 max: 0.373212770233
deltaT = 1.72109219044e-10
Time = 0.01444356487

Centre of mass: (-0.349814340732 -0.125696851454 0.196298253316)
Linear velocity: (12649.2884767 -1046153.3357 -294487.425168)
Angular velocity: (-12316886.264 -54713.605122 -180578.983616)
Centre of mass: (0.34981707595 -0.124358692651 -0.196568327016)
Linear velocity: (-13810.1784143 -823605.271104 218608.055854)
Angular velocity: (9409918.23348 -78418.1421569 197937.084767)
GAMG:  Solving for cellDisplacementx, Initial residual = 0.00914785781768, Final residual = 3.37125821226e-06, No Iterations 4
GAMG:  Solving for cellDisplacementy, Initial residual = 0.00552863013895, Final residual = 8.90932047055e-06, No Iterations 3
GAMG:  Solving for cellDisplacementz, Initial residual = 0.005588990035, Final residual = 3.24807476154e-06, No Iterations 4
Execution time for mesh.update() = 3.47 s
time step continuity errors : sum local = 4.13760397787e-10, global = -3.54569889692e-12, cumulative = -4.14669791128e-06
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 5.59830547694e-06, No Iterations 5
time step continuity errors : sum local = 2.31635707417e-15, global = -2.50480722521e-17, cumulative = -4.1466979113e-06
MULES: Solving for alpha1
Liquid phase volume fraction = 0.500498189336  Min(alpha1) = -5.38274365806e-18  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.500534512433  Min(alpha1) = -4.03241934233e-20  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.000412929757808, Final residual = 2.02051992133e-06, No Iterations 5
time step continuity errors : sum local = 9.4694440486e-08, global = -3.84738504627e-08, cumulative = -4.18517176177e-06
GAMGPCG:  Solving for p_rgh, Initial residual = 0.000120911404203, Final residual = 5.06119181433e-09, No Iterations 4
time step continuity errors : sum local = 3.3542482784e-10, global = -1.14461031922e-12, cumulative = -4.18517290638e-06
smoothSolver:  Solving for epsilon, Initial residual = 0.075115075581, Final residual = 2.99738298272e-07, No Iterations 7
smoothSolver:  Solving for k, Initial residual = 0.0990449776516, Final residual = 4.83857085574e-07, No Iterations 9
ExecutionTime = 7134.3 s  ClockTime = 7356 s

Interface Courant Number mean: 0.000265597013331 max: 0.46347609714
Courant Number mean: 0.00660038424558 max: 0.46347609714
deltaT = 9.2836080667e-11
Time = 0.01444356497

Centre of mass: (-0.349813644548 -0.125761293767 0.196278901707)
Linear velocity: (3932.20152874 -387209.391939 -116728.243164)
Angular velocity: (-4320019.60806 -17361.820749 -58813.3686743)
Centre of mass: (0.349816233099 -0.124407234488 -0.196555340918)
Linear velocity: (-4956.92370943 -300428.028058 76014.0091893)
Angular velocity: (3204649.11753 -27042.7989216 70933.421513)
GAMG:  Solving for cellDisplacementx, Initial residual = 0.00246801333602, Final residual = 6.28707160732e-06, No Iterations 3
GAMG:  Solving for cellDisplacementy, Initial residual = 0.00145454425153, Final residual = 4.92632806556e-06, No Iterations 3
GAMG:  Solving for cellDisplacementz, Initial residual = 0.00148012964054, Final residual = 6.56529965109e-06, No Iterations 3
Execution time for mesh.update() = 3.31 s
time step continuity errors : sum local = 1.80928869248e-10, global = -6.17405253258e-13, cumulative = -4.18517352378e-06
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 7.89014358714e-06, No Iterations 5
time step continuity errors : sum local = 1.42755490437e-15, global = -1.05503818049e-16, cumulative = -4.18517352389e-06
MULES: Solving for alpha1
Liquid phase volume fraction = 0.500526484771  Min(alpha1) = -1.14193700815e-20  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.500537597551  Min(alpha1) = -2.80569412702e-22  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.000725206106519, Final residual = 3.81567763131e-06, No Iterations 6
time step continuity errors : sum local = 7.49250057179e-08, global = 4.02202015722e-08, cumulative = -4.14495332231e-06
GAMGPCG:  Solving for p_rgh, Initial residual = 0.000207894126332, Final residual = 5.59511538545e-09, No Iterations 4
time step continuity errors : sum local = 1.92976789873e-10, global = 3.40385082818e-12, cumulative = -4.14494991846e-06
smoothSolver:  Solving for epsilon, Initial residual = 0.011047608263, Final residual = 3.1070352165e-07, No Iterations 5
bounding epsilon, min: -510464.980869 max: 1.25574249239e+19 average: 1.05141718141e+15
smoothSolver:  Solving for k, Initial residual = 0.016215050954, Final residual = 9.9945203402e-07, No Iterations 93
bounding k, min: -164901.521304 max: 2.3334191091e+12 average: 1462769853.14
ExecutionTime = 7150.81 s  ClockTime = 7373 s

Interface Courant Number mean: 6.95458888e-05 max: 0.149152602651
Courant Number mean: 0.00173898180367 max: 0.149152602651
deltaT = 1.08396666746e-10
Time = 0.01444356507

Centre of mass: (-0.349813604946 -0.125769994463 0.196276191003)
Linear velocity: (2049.47828255 -680546.86927 -172941.572078)
Angular velocity: (-7478075.25696 -6598.23653197 -44771.4877576)
Centre of mass: (0.349816142596 -0.124415687175 -0.196554024332)
Linear velocity: (-9073.05806102 -591770.188542 147068.65847)
Angular velocity: (6482924.26627 -49343.8084364 133078.931852)
GAMG:  Solving for cellDisplacementx, Initial residual = 0.000295146175863, Final residual = 5.97019978886e-06, No Iterations 2
GAMG:  Solving for cellDisplacementy, Initial residual = 0.000218583751145, Final residual = 4.22958199331e-06, No Iterations 2
GAMG:  Solving for cellDisplacementz, Initial residual = 0.000198654046882, Final residual = 6.55616581993e-06, No Iterations 2
Execution time for mesh.update() = 3.02 s
time step continuity errors : sum local = 2.2532231667e-10, global = 3.97438238407e-12, cumulative = -4.14494594408e-06
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 4.51009197064e-06, No Iterations 6
time step continuity errors : sum local = 1.01622430033e-15, global = 1.66777362241e-17, cumulative = -4.14494594406e-06
MULES: Solving for alpha1
Liquid phase volume fraction = 0.500538350555  Min(alpha1) = -4.36546661836e-21  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.500540610469  Min(alpha1) = -4.35752954849e-21  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.000292468857355, Final residual = 2.23564579873e-06, No Iterations 4
time step continuity errors : sum local = 1.17614932829e-07, global = -4.08291430287e-08, cumulative = -4.18577508709e-06
GAMGPCG:  Solving for p_rgh, Initial residual = 0.00138218498679, Final residual = 1.84084382503e-09, No Iterations 5
time step continuity errors : sum local = 2.81116504523e-10, global = -3.61968945516e-12, cumulative = -4.18577870678e-06
smoothSolver:  Solving for epsilon, Initial residual = 0.00345703998198, Final residual = 5.09397214168e-07, No Iterations 12
bounding epsilon, min: -1793882202.49 max: 5.93459191412e+19 average: 2.74891994338e+15
smoothSolver:  Solving for k, Initial residual = 0.956101528916, Final residual = 0.455462713576, No Iterations 1000
bounding k, min: -49197.00595 max: 5.17537636805e+15 average: 126756438419
ExecutionTime = 7191.28 s  ClockTime = 7413 s

Interface Courant Number mean: 4.6025422324e-05 max: 5.74951952632
Courant Number mean: 0.00094935070414 max: 74.7795030671
deltaT = 3.62387627292e-13
Time = 0.01444356508

Centre of mass: (-0.349813603593 -0.125770458619 0.196276074721)
Linear velocity: (-199870.392266 -1281286.18559 -321764.799104)
Angular velocity: (-14508169.6235 1064848.82357 2146686.30104)
Centre of mass: (0.349816136322 -0.124416087817 -0.196553922142)
Linear velocity: (-17332.0760682 -1106798.10145 282366.944374)
Angular velocity: (12470116.3121 -96147.2441149 253658.466637)
GAMG:  Solving for cellDisplacementx, Initial residual = 4.08774120991e-05, Final residual = 6.11355013872e-06, No Iterations 1
GAMG:  Solving for cellDisplacementy, Initial residual = 2.72076026644e-05, Final residual = 3.89219989053e-06, No Iterations 1
GAMG:  Solving for cellDisplacementz, Initial residual = 3.78022826952e-05, Final residual = 6.7544189797e-06, No Iterations 1
Execution time for mesh.update() = 2.76 s
time step continuity errors : sum local = 9.39818042371e-13, global = -1.21012085843e-14, cumulative = -4.18577871888e-06
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 9.30087569536e-06, No Iterations 5
time step continuity errors : sum local = 8.74112085685e-18, global = -9.61087795481e-19, cumulative = -4.18577871888e-06
MULES: Solving for alpha1
Liquid phase volume fraction = 0.500540239206  Min(alpha1) = -2.22896482145e-21  Max(alpha1) = 1
--> FOAM Warning : 
    From function Time::operator++()
    in file db/Time/Time.C at line 937
    Increased the timePrecision from 10 to 11 to distinguish between timeNames at time 0.0144435650752
MULES: Solving for alpha1
Liquid phase volume fraction = 0.500540589846  Min(alpha1) = -1.24863124537e-23  Max(alpha1) = 1
--> FOAM Warning : 
    From function Time::operator++()
    in file db/Time/Time.C at line 937
    Increased the timePrecision from 11 to 12 to distinguish between timeNames at time 0.0144435650754
GAMG:  Solving for p_rgh, Initial residual = 0.543802456097, Final residual = 0.00350015483377, No Iterations 4
time step continuity errors : sum local = 1.33565722745e-08, global = -1.47244315827e-09, cumulative = -4.18725116204e-06
GAMGPCG:  Solving for p_rgh, Initial residual = 1.38260922845e-06, Final residual = 6.14744253446e-09, No Iterations 3
time step continuity errors : sum local = 0.23336264053, global = -0.00361142303455, cumulative = -0.00361561028571
smoothSolver:  Solving for epsilon, Initial residual = 0.999994172508, Final residual = 7.77978833782e-07, No Iterations 100
smoothSolver:  Solving for k, Initial residual = 0.68037087876, Final residual = 0.0738085489472, No Iterations 1000
ExecutionTime = 7232.47 s  ClockTime = 7454 s

Interface Courant Number mean: 3.63949734232 max: 994593.962793
Courant Number mean: 77.1411855051 max: 25125249.7756
deltaT = 3.60581119121e-21
Time = 0.0144435650752

Centre of mass: (-0.349813603593 -0.125770458619 0.196276074721)
Linear velocity: (-26538645334.4 64962531.0391 -30616967.994)
Angular velocity: (64724043.6967 148919030194 248829311320)
Centre of mass: (0.349816136322 -0.124416087817 -0.196553922142)
Linear velocity: (-17658.4200507 -1110505.36576 282661.546484)
Angular velocity: (12503044.3746 -98176.7068542 258394.570317)
GAMG:  Solving for cellDisplacementx, Initial residual = 1.68520346229e-05, Final residual = 4.40230661425e-06, No Iterations 1
GAMG:  Solving for cellDisplacementy, Initial residual = 1.01648094939e-05, Final residual = 2.82058347291e-06, No Iterations 1
GAMG:  Solving for cellDisplacementz, Initial residual = 1.94435624256e-05, Final residual = 4.96067629504e-06, No Iterations 1
Execution time for mesh.update() = 2.76 s
time step continuity errors : sum local = 2.32199323754e-09, global = -3.59342002862e-11, cumulative = -0.00361561032164
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 6.34199512265e-06, No Iterations 5
time step continuity errors : sum local = 1.47260697833e-14, global = -1.79753525488e-15, cumulative = -0.00361561032165
--> FOAM Warning : 
    From function Time::operator++()
    in file db/Time/Time.C at line 937
    Increased the timePrecision from 12 to 13 to distinguish between timeNames at time 0.0144435650752
MULES: Solving for alpha1
[0] #0  Foam::error::printStack(Foam::Ostream&)[1] #0  Foam::error::printStack(Foam::Ostream&) in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #1  Foam::sigFpe::sigHandler(int) in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #1  Foam::sigFpe::sigHandler(int) in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #2   in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #2   in "/lib/x86_64-linux-gnu/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, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/lib/x86_64-linux-gnu/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, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[1] #4  void Foam::MULES::explicitSolve<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::zeroField const&, Foam::zeroField const&, double, double) in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #4  void Foam::MULES::explicitSolve<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::zeroField const&, Foam::zeroField const&, double, double) in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[1] #5  Foam::MULES::explicitSolve(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #5  Foam::MULES::explicitSolve(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[1] #6   in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #6  

[1]  in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[1] #7  __libc_start_main[0]  in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[0] #7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #8   in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #8  

[1]  in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[ceeckh:06398] *** Process received signal ***
[ceeckh:06398] Signal: Floating point exception (8)
[ceeckh:06398] Signal code:  (-6)
[ceeckh:06398] Failing at address: 0x3e8000018fe
[ceeckh:06398] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x33d80) [0x7fca9b968d80]
[ceeckh:06398] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7fca9b968d05]
[ceeckh:06398] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x33d80) [0x7fca9b968d80]
[ceeckh:06398] [ 3] /home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES7limiterINS_17geometricOneFieldENS_9zeroFieldES3_EEvRNS_5FieldIdEERKT_RKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNSA_IdNS_13fvsPatchFieldENS_11surfaceMeshEEESK_RKT0_RKT1_ddi+0x11e0) [0x7fca9e34e810]
[ceeckh:06398] [ 4] /home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES13explicitSolveINS_17geometricOneFieldENS_9zeroFieldES3_EEvRKT_RNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNS7_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERSE_RKT0_RKT1_dd+0x2a0) [0x7fca9e3501a0]
[ceeckh:06398] [ 5] /home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES13explicitSolveERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERS8_dd+0x21) [0x7fca9e342491]
[ceeckh:06398] [ 6] interDyMFoam() [0x42d0fe]
[ceeckh:06398] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xff) [0x7fca9b953eff]
[ceeckh:06398] [ 8] interDyMFoam() [0x4243b9]
[ceeckh:06398] *** End of error message ***
[0]  in "/home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/bin/interDyMFoam"
[ceeckh:06397] *** Process received signal ***
[ceeckh:06397] Signal: Floating point exception (8)
[ceeckh:06397] Signal code:  (-6)
[ceeckh:06397] Failing at address: 0x3e8000018fd
[ceeckh:06397] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x33d80) [0x7f36bc39dd80]
[ceeckh:06397] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7f36bc39dd05]
[ceeckh:06397] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x33d80) [0x7f36bc39dd80]
[ceeckh:06397] [ 3] /home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES7limiterINS_17geometricOneFieldENS_9zeroFieldES3_EEvRNS_5FieldIdEERKT_RKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNSA_IdNS_13fvsPatchFieldENS_11surfaceMeshEEESK_RKT0_RKT1_ddi+0x11e0) [0x7f36bed83810]
[ceeckh:06397] [ 4] /home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES13explicitSolveINS_17geometricOneFieldENS_9zeroFieldES3_EEvRKT_RNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNS7_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERSE_RKT0_RKT1_dd+0x2a0) [0x7f36bed851a0]
[ceeckh:06397] [ 5] /home/kh/OpenFOAM/OpenFOAM-2.0.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES13explicitSolveERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERS8_dd+0x21) [0x7f36bed77491]
[ceeckh:06397] [ 6] interDyMFoam() [0x42d0fe]
[ceeckh:06397] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xff) [0x7f36bc388eff]
[ceeckh:06397] [ 8] interDyMFoam() [0x4243b9]
[ceeckh:06397] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 1 with PID 6398 on node ceeckh exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------
2 total processes killed (some possibly by mpirun during cleanup)

Thank you once again!

jhoepken · September 16, 2011, 03:41

I've never used the floatingBodies in this release, so please correct me if I am wrong: Your log file shows the linear and angular velocities of both bodies, right? The linear velocities seem to be unrealistically high. Is it possible, that the motion of the bodies cause the high velocities? Maybe you can try a simulation on the same case but fixinig all DOF of each body. If this does not break, the high velocities and the crash of the simulation may be related to the motion of the bodies.
I hope this helps.

CKH · September 22, 2011, 00:13

Just an update on the case. I recreated a mesh that's closer to what I'm looking at (more ship-hull-like) rather than two square blocks.

Jens: I saw the high velocities too and was thinking that it may be due to my inputs under 6dofrigidbodymotions. I did what you've suggested, i.e. to fix the two bodies and the case seems to run ok (although it took the past few days of parallel running on my 4 core desktop!)

I am not entirely sure if this is correct, but another possible reason that the original (two floating bodies) case was crashing could be due to the fact that I've used a uniform mesh (grading of 1) throughout, when I should use a finer grading closer to the hull.

I'm gonna to carry on working towards simulating two floating vessels in waves and will try to post more updates of my progress (and of course problems).

Ferdinand · March 7, 2016, 13:26

Hi CKH;

You could solve the problem ?. What was the mistake? Could you give more information about this.

Best regards

xoitx · March 21, 2016, 14:05

Hey
Did you manage to solve it?
Thanks

cfd lover · February 20, 2018, 04:17

I finally solved it. In my case the problem was about k and epsilon files. The internal field and inlet value was not accurate. Although i have solved this error and the case is running now, it has not converged yet.

FloatingSolar · February 20, 2019, 10:08

Hey,

CFD lover, would you be willing to upload your files? I am also trying to simulate multiple floating objects, but I can't get it to work.

Kind regards,
Daan

February 20, 2018, 04:17	Floating point exception error	#14
cfd lover New Member Ali Mohammadi Join Date: Oct 2017 Posts: 15 Rep Power: 9	I finally solved it. In my case the problem was about k and epsilon files. The internal field and inlet value was not accurate. Although i have solved this error and the case is running now, it has not converged yet.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
how to set periodic boundary conditions	Ganesh	FLUENT	15	November 18, 2020 07:09
OpenFOAM static build on Cray XT5	asaijo	OpenFOAM Installation	9	April 6, 2011 13:21
Inserting Boundary Layers in multiple objects in Gambit for use in FLUENT	3quarters	FLUENT	0	July 16, 2010 09:55
Multiple domain/Multiple material problems... possible or not?	Amiga500	OpenFOAM	1	April 14, 2010 16:32
How to simulate floating objects on the sea?	RX	Siemens	2	December 16, 2006 17:31

September 13, 2011, 04:26		#2
jhoepken Senior Member Jens Höpken Join Date: Apr 2009 Location: Duisburg, Germany Posts: 159 Rep Power: 17	Hi, From what I can see in the very short log file excerpt, I am guessing that your time step is way to large. Your alpha is exploding. Do you use the CFL based time step control? P.S.: A longer excerpt of the logfile would have been helpful.

September 13, 2011, 05:58		#4
elvis Senior Member Elvis Join Date: Mar 2009 Location: Sindelfingen, Germany Posts: 620 Blog Entries: 6 Rep Power: 24	Hi, are you trying s.th. like http://powerlab.fsb.hr/ped/kturbo/Op...F/twoBoats.mpg ? I believe there is a extend http://openfoam-extend.git.sourcefor...enFOAM-1.6-ext =>tutorials/basic/sixDOFSolver/sixDOFmotion with bodyA & bodyB but I am not sure

September 13, 2011, 06:04		#5
CKH New Member Kie Hian Join Date: Aug 2011 Location: Singapore Posts: 19 Rep Power: 15	Hi Elvis, Thanks for pointing out the two boats case. I have looked at the extend version and actually, my blockMeshDict file is based on the case. However, instead of subjecting the two boxes to a steady flow, what I am plan to do is to subject them to waves. So as a first step, I tried to modify the floatingObject case to see if it works for two objects and if it does, my next step is to implement some form of wave input.

September 13, 2011, 08:11		#6
jhoepken Senior Member Jens Höpken Join Date: Apr 2009 Location: Duisburg, Germany Posts: 159 Rep Power: 17	About your crashing simulation, you might want to check the following: Gravity pointing into the correct direction? Field (alpha) initialised correctly? Are there bad cells in your mesh?

September 14, 2011, 17:30		#8
jhoepken Senior Member Jens Höpken Join Date: Apr 2009 Location: Duisburg, Germany Posts: 159 Rep Power: 17	Again, a new log file would help a lot

September 16, 2011, 03:41		#10
jhoepken Senior Member Jens Höpken Join Date: Apr 2009 Location: Duisburg, Germany Posts: 159 Rep Power: 17	I've never used the floatingBodies in this release, so please correct me if I am wrong: Your log file shows the linear and angular velocities of both bodies, right? The linear velocities seem to be unrealistically high. Is it possible, that the motion of the bodies cause the high velocities? Maybe you can try a simulation on the same case but fixinig all DOF of each body. If this does not break, the high velocities and the crash of the simulation may be related to the motion of the bodies. I hope this helps.

September 22, 2011, 00:13		#11
CKH New Member Kie Hian Join Date: Aug 2011 Location: Singapore Posts: 19 Rep Power: 15	Just an update on the case. I recreated a mesh that's closer to what I'm looking at (more ship-hull-like) rather than two square blocks. Jens: I saw the high velocities too and was thinking that it may be due to my inputs under 6dofrigidbodymotions. I did what you've suggested, i.e. to fix the two bodies and the case seems to run ok (although it took the past few days of parallel running on my 4 core desktop!) I am not entirely sure if this is correct, but another possible reason that the original (two floating bodies) case was crashing could be due to the fact that I've used a uniform mesh (grading of 1) throughout, when I should use a finer grading closer to the hull. I'm gonna to carry on working towards simulating two floating vessels in waves and will try to post more updates of my progress (and of course problems). Ferdinand likes this.

March 7, 2016, 13:26		#12
Ferdinand New Member Luis Fernando Join Date: Nov 2013 Location: Perú Posts: 19 Rep Power: 13	Hi CKH; You could solve the problem ?. What was the mistake? Could you give more information about this. Best regards

March 21, 2016, 14:05		#13
xoitx New Member Alpha Beta Join Date: Mar 2016 Posts: 28 Rep Power: 10	Hey Did you manage to solve it? Thanks

February 20, 2019, 10:08		#15
FloatingSolar New Member Daan Join Date: Feb 2019 Posts: 4 Rep Power: 7	Hey, CFD lover, would you be willing to upload your files? I am also trying to simulate multiple floating objects, but I can't get it to work. Kind regards, Daan