CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

Multiple floating objects

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By CKH

Reply
 
LinkBack Thread Tools Display Modes
Old   September 12, 2011, 09:27
Default Multiple floating objects
  #1
CKH
New Member
 
Kie Hian
Join Date: Aug 2011
Location: Singapore
Posts: 19
Rep Power: 6
CKH is on a distinguished road
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
Attached Files
File Type: gz blockMeshDict.tar.gz (2.5 KB, 24 views)
File Type: gz system.tar.gz (1.7 KB, 26 views)
CKH is offline   Reply With Quote

Old   September 13, 2011, 03:26
Default
  #2
Senior Member
 
Jens Höpken
Join Date: Apr 2009
Location: Duisburg, Germany
Posts: 156
Rep Power: 8
jhoepken is on a distinguished road
Send a message via Skype™ to jhoepken
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.
jhoepken is offline   Reply With Quote

Old   September 13, 2011, 04:15
Default
  #3
CKH
New Member
 
Kie Hian
Join Date: Aug 2011
Location: Singapore
Posts: 19
Rep Power: 6
CKH is on a distinguished road
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.
CKH is offline   Reply With Quote

Old   September 13, 2011, 04:58
Default
  #4
Senior Member
 
Elvis
Join Date: Mar 2009
Location: Sindelfingen, Germany
Posts: 579
Blog Entries: 5
Rep Power: 13
elvis is on a distinguished road
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
elvis is offline   Reply With Quote

Old   September 13, 2011, 05:04
Default
  #5
CKH
New Member
 
Kie Hian
Join Date: Aug 2011
Location: Singapore
Posts: 19
Rep Power: 6
CKH is on a distinguished road
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.
CKH is offline   Reply With Quote

Old   September 13, 2011, 07:11
Default
  #6
Senior Member
 
Jens Höpken
Join Date: Apr 2009
Location: Duisburg, Germany
Posts: 156
Rep Power: 8
jhoepken is on a distinguished road
Send a message via Skype™ to jhoepken
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?
jhoepken is offline   Reply With Quote

Old   September 14, 2011, 06:14
Default
  #7
CKH
New Member
 
Kie Hian
Join Date: Aug 2011
Location: Singapore
Posts: 19
Rep Power: 6
CKH is on a distinguished road
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
Attached Images
File Type: jpg floatObj2.jpg (12.1 KB, 128 views)
CKH is offline   Reply With Quote

Old   September 14, 2011, 16:30
Default
  #8
Senior Member
 
Jens Höpken
Join Date: Apr 2009
Location: Duisburg, Germany
Posts: 156
Rep Power: 8
jhoepken is on a distinguished road
Send a message via Skype™ to jhoepken
Again, a new log file would help a lot
jhoepken is offline   Reply With Quote

Old   September 14, 2011, 23:20
Default
  #9
CKH
New Member
 
Kie Hian
Join Date: Aug 2011
Location: Singapore
Posts: 19
Rep Power: 6
CKH is on a distinguished road
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!
Attached Files
File Type: gz 0.org.tar.gz (1.3 KB, 33 views)
File Type: gz system.tar.gz (1.8 KB, 25 views)
File Type: gz constant_blockMeshDict.tar.gz (3.1 KB, 31 views)
File Type: gz log.mpirun_part1.tar.gz (59.4 KB, 16 views)
File Type: gz log.mpirun_part2.tar.gz (72.5 KB, 11 views)
chun likes this.
CKH is offline   Reply With Quote

Old   September 16, 2011, 02:41
Default
  #10
Senior Member
 
Jens Höpken
Join Date: Apr 2009
Location: Duisburg, Germany
Posts: 156
Rep Power: 8
jhoepken is on a distinguished road
Send a message via Skype™ to jhoepken
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.
jhoepken is offline   Reply With Quote

Old   September 21, 2011, 23:13
Default
  #11
CKH
New Member
 
Kie Hian
Join Date: Aug 2011
Location: Singapore
Posts: 19
Rep Power: 6
CKH is on a distinguished road
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).
CKH is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 21:16.