|
[Sponsors] |
January 18, 2013, 10:18 |
|
#201 |
Member
vincent
Join Date: Apr 2011
Posts: 45
Rep Power: 15 |
Hi Colin
Thanks for your quite reply, the error message is following: [1] #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [1] #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [1] #2 in "/lib/x86_64-linux-gnu/libc.so.6" [1] #3 Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [1] #4 Foam::GaussSeidelSmoother::smooth(Foam::Field<doub le>&, Foam::Field<double> const&, unsigned char, int) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [1] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [1] #6 [1] in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/LTSInterFoam" [1] #7 [1] in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/LTSInterFoam" [1] #8 [1] in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/LTSInterFoam" [1] #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" [1] #10 [1] in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/LTSInterFoam" [fluxiz-xeon2670:07523] *** Process received signal *** [fluxiz-xeon2670:07523] Signal: Floating point exception (8) [fluxiz-xeon2670:07523] Signal code: (-6) [fluxiz-xeon2670:07523] Failing at address: 0x3e800001d63 [fluxiz-xeon2670:07523] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7f8e36fea4a0] [fluxiz-xeon2670:07523] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7f8e36fea425] [fluxiz-xeon2670:07523] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7f8e36fea4a0] [fluxiz-xeon2670:07523] [ 3] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam19GaussSeidelSmoother6smoot hERKNS_4wordERNS_5FieldIdEERKNS_9lduMatrixERKS5_RK NS_10FieldFieldIS4_dEERKNS_8UPtrListIKNS_17lduInte rfaceFieldEEEhi+0x1f4) [0x7f8e3801fa34] [fluxiz-xeon2670:07523] [ 4] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam19GaussSeidelSmoother6smoo thERNS_5FieldIdEERKS2_hi+0x2d) [0x7f8e3801ffdd] [fluxiz-xeon2670:07523] [ 5] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam12smoothSolver5solveERNS_5 FieldIdEERKS2_h+0x40a) [0x7f8e38019efa] [fluxiz-xeon2670:07523] [ 6] LTSInterFoam() [0x45a698] [fluxiz-xeon2670:07523] [ 7] LTSInterFoam() [0x45b080] [fluxiz-xeon2670:07523] [ 8] LTSInterFoam() [0x42f901] [fluxiz-xeon2670:07523] [ 9] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7f8e36fd576d] [fluxiz-xeon2670:07523] [10] LTSInterFoam() [0x4334cd] [fluxiz-xeon2670:07523] *** End of error message *** -------------------------------------------------------------------------- mpirun noticed that process rank 1 with PID 7523 on node fluxiz-xeon2670 exited on signal 8 (Floating point exception). -------------------------------------------------------------------------- I atteched the checkMesh and snappyHexMesh files. Regards Vince |
|
January 18, 2013, 10:31 |
|
#202 |
Senior Member
Join Date: Aug 2010
Location: Groningen, The Netherlands
Posts: 216
Rep Power: 19 |
What are your computer specifications?
Anything under 40 GB of RAM might cause the Error! You have 37M cells and you need approximately 1 GB of RAM per 0.8M cells To further check this: what was the number of cells in the previous calculations? what are your computer specifications? Do you get a smaller ratio of cells per GB RAM with the earlier calculations than with the latest calculation, then this is for sure the error source. I hope this is the solution regards Colin Last edited by colinB; January 18, 2013 at 10:47. |
|
January 18, 2013, 10:53 |
|
#203 |
Member
vincent
Join Date: Apr 2011
Posts: 45
Rep Power: 15 |
I use a computer with a bi-xeon E5, 95Go of RAM. So I don't think there is a problem with RAM.
In this case, when I change the scheme of div(rho*phi, U) to linuearUpwind, the calculation run without problem but the result is very bad. So what's happen with the limitedLinear scheme? regards Vince |
|
January 21, 2013, 04:42 |
|
#204 |
Senior Member
Join Date: Aug 2010
Location: Groningen, The Netherlands
Posts: 216
Rep Power: 19 |
Dear Vince,
sorry, when it comes to judge on how the schemes impair the calculation, my knowledge is rather limited and I cannot help you. However it would be also interesting to me what the solution looks like to me. regards Colin |
|
January 27, 2013, 11:49 |
|
#205 |
Member
Ali Khalifesoltani
Join Date: Mar 2011
Location: Esfahan, Iran
Posts: 56
Rep Power: 15 |
Hello Eric,
I think the your link has some problems. It doesn't open. Could you possibly post the right link? Regards, Ali |
|
January 27, 2013, 15:40 |
|
#206 |
Senior Member
Ralph Moolenaar
Join Date: Aug 2010
Location: 's-Hertogenbosch, the Netherlands
Posts: 120
Rep Power: 16 |
Dear Ali,
I'm not sure if Eric will read your post, but the latest version of OF21 has a Wigley hull tutorial included in the LTSInterFOAM solver I think. Is that what you're looking for? Ralph
__________________
CFD for marine applications? Go to http://www.marinecfd.com/ and join the OF Ship Hydromechanics Group: http://www.cfd-online.com/Forums/gro...mechanics.html |
|
January 28, 2013, 02:19 |
|
#207 |
Member
Ali Khalifesoltani
Join Date: Mar 2011
Location: Esfahan, Iran
Posts: 56
Rep Power: 15 |
Hi Ralph,
I will take a look at the tutorial. thanks for your reply. Regards, Ali |
|
February 3, 2013, 04:03 |
|
#208 |
New Member
hs
Join Date: Mar 2012
Posts: 22
Rep Power: 14 |
Hi Ben,
I'v just started interDyMFoam. blockMesh and sHM were completed but interDyMFoam got errors as below. Could you give me some guide for me, please? // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Selecting dynamicFvMesh dynamicMotionSolverFvMesh Selecting motion solver: displacementLaplacian Selecting motion diffusion: inverseDistance Reading field p_rgh Reading field alpha1 Reading field U Reading/calculating face flux field phi Reading transportProperties Selecting incompressible transport model Newtonian Selecting incompressible transport model Newtonian Selecting turbulence model type RASModel Selecting RAS turbulence model kEpsilon kEpsilonCoeffs { Cmu 0.09; C1 1.44; C2 1.92; sigmaEps 1.3; } Reading g Calculating field g.h PIMPLE: Operating solver in PISO mode time step continuity errors : sum local = 0.250418602578, global = -0.250418602578, cumulative = -0.250418602578 GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 6.73986681523e-07, No Iterations 5 time step continuity errors : sum local = 1.68778803235e-07, global = 1.97915736808e-08, cumulative = -0.250418582786 Courant Number mean: 4.08539047445 max: 206.187674774 Starting time loop Interface Courant Number mean: 0 max: 0 Courant Number mean: 0.00990637845405 max: 0.499970113419 deltaT = 0.00242483026188 Time = 0.00242483 sixDoFRigidBodyMotion constraints converged in 1 iterations Constraint force: (-2.15121188445e-05 -4.59773475343e-05 0) Constraint moment: (0 0 0) Centre of mass: (0 0 0) Linear velocity: (1.16435796365e-09 2.48855499246e-09 -0.00848442071477) Angular velocity: (7.64408796861e-08 1.78424900329e-06 8.65034058421e-07) GAMG: Solving for cellDisplacementx, Initial residual = 1, Final residual = 4.29473890619e-06, No Iterations 9 GAMG: Solving for cellDisplacementy, Initial residual = 0, Final residual = 0, No Iterations 0 GAMG: Solving for cellDisplacementz, Initial residual = 0, Final residual = 0, No Iterations 0 Execution time for mesh.update() = 3.17 s time step continuity errors : sum local = 4.76660924674e-10, global = 5.58948732343e-11, cumulative = -0.25041858273 GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 5.58784483893e-06, No Iterations 5 time step continuity errors : sum local = 0.288013547787, global = 0.288012658057, cumulative = 0.037594075327 MULES: Solving for alpha1 Phase-1 volume fraction = -3.08729405329e+294 Min(alpha1) = -1.51690089806e+296 Max(alpha1) = 2.81618849836e+296 #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::PhiScheme<double, Foam::interfaceCompressionLimiter>::limiter(Foam:: GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libinterfaceProperties.so" #4 Foam::limitedSurfaceInterpolationScheme<double>::w eights(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #5 Foam::surfaceInterpolationScheme<double>::interpol ate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #6 Foam::fv::gaussConvectionScheme<double>::interpola te(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #7 Foam::fv::gaussConvectionScheme<double>::flux(Foam ::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #8 at interDyMFoam.C:0 #9 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/interDyMFoam" #10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #11 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/interDyMFoam" Floating point exception (core dumped) |
|
February 3, 2013, 07:18 |
|
#209 |
Senior Member
Ralph Moolenaar
Join Date: Aug 2010
Location: 's-Hertogenbosch, the Netherlands
Posts: 120
Rep Power: 16 |
Hey Hs,
You can see in the output that alpha1 is not behaving well, the min/max value are way beyond the range of [0-1]. Did you used setFields to 'fill' your domain and what are you trying to simulate? What is your output of "checkMesh -latestTime"? Cheers, Ralph
__________________
CFD for marine applications? Go to http://www.marinecfd.com/ and join the OF Ship Hydromechanics Group: http://www.cfd-online.com/Forums/gro...mechanics.html |
|
February 3, 2013, 16:26 |
|
#210 |
New Member
hs
Join Date: Mar 2012
Posts: 22
Rep Power: 14 |
Hi Ralph,
I am trying to simulate the 2DoF wigleyHull in interDyMFoam. I'v modified files in "floatingObject" and run the case. The output of "checkMesh" is as below. Thanks, hs// // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create polyMesh for time = 0 Time = 0 Mesh stats points: 209138 faces: 574173 internal faces: 545493 cells: 182759 boundary patches: 5 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 173693 prisms: 688 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 0 polyhedra: 8378 Checking topology... Boundary definition OK. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces ... Patch Faces Points Surface topology inlet 1422 1565 ok (non-closed singly connected) outlet 1110 1247 ok (non-closed singly connected) atmosphere 19416 19813 ok (non-closed singly connected) sides 4160 4503 ok (non-closed singly connected) hull_wall 2572 2801 ok (non-closed singly connected) Checking geometry... Overall domain bounding box (-1.5 -0.75 -0.5000000357) (2.5 0.75 0.0399) Mesh (non-empty, non-wedge) directions (1 1 1) Mesh (non-empty) directions (1 1 1) Boundary openness (-3.52092954661e-16 2.36648403091e-16 -4.55752380682e-14) OK. Max cell openness = 3.21304649248e-16 OK. Max aspect ratio = 7.81781560787 OK. Minumum face area = 1.68314246315e-05. Maximum face area = 0.00540716576367. Face area magnitudes OK. Min volume = 9.41898154846e-08. Max volume = 0.000270208268515. Total volume = 3.23398498218. Cell volumes OK. Mesh non-orthogonality Max: 52.1678045455 average: 8.40504498863 Non-orthogonality check OK. Face pyramids OK. Max skewness = 1.34743015352 OK. Coupled point location match (average 0) OK. Mesh OK. End |
|
February 4, 2013, 11:31 |
|
#211 |
Member
Ali Khalifesoltani
Join Date: Mar 2011
Location: Esfahan, Iran
Posts: 56
Rep Power: 15 |
Hello guys,
I am simulating a planing boat with LTSInterFoam. I have simulated the steady problem with 0 DOF, but I want to simulate the history of craft's motion with 2 or more DOF. My question is how to implement the rigid body motion to the craft? Is the solver a good choice for transient problem or should I couple the solver with another rigid body motion solver? any response is appreciated, Best regards, Ali |
|
February 4, 2013, 15:48 |
|
#212 |
Senior Member
Ralph Moolenaar
Join Date: Aug 2010
Location: 's-Hertogenbosch, the Netherlands
Posts: 120
Rep Power: 16 |
@Hs: did you also checked your BC's? When checkmesh is okay the error must be somewhere else....
@Ali: Check google for navyFoam of shipFoam...both for OF16. Certainly use interFoam as a basis since LTS isn't transient! Cheers, Ralph
__________________
CFD for marine applications? Go to http://www.marinecfd.com/ and join the OF Ship Hydromechanics Group: http://www.cfd-online.com/Forums/gro...mechanics.html |
|
February 4, 2013, 16:10 |
|
#213 |
New Member
hs
Join Date: Mar 2012
Posts: 22
Rep Power: 14 |
Hi Ralph,
thanks for your reply, I will check BCs again. hs// |
|
February 5, 2013, 10:46 |
|
#214 |
Member
Ali Khalifesoltani
Join Date: Mar 2011
Location: Esfahan, Iran
Posts: 56
Rep Power: 15 |
Hi Ralph,
I searched the web and found a little documents about navyFoam that was not so useful, and some materials about shipFoam containing a set of folders(0 , constant, system) but it appears that the solver itself is not in the package. could you possibly give me the related materials about the shipFoam and navyFoam if you have or guide me how to find them? Thanks in advance, Regards, Ali |
|
February 6, 2013, 07:37 |
|
#215 |
Member
vincent
Join Date: Apr 2011
Posts: 45
Rep Power: 15 |
Hi all
Just a little question, maybe stupid but I would like to know if somebody have made some test with differents quality mesh on STL CAD? There is a dependence between the number and quality cells on STL CAD and final result in OpenFOAM? Thanks a lot sheers Vince |
|
February 8, 2013, 19:30 |
|
#216 | |
New Member
hs
Join Date: Mar 2012
Posts: 22
Rep Power: 14 |
Quote:
I checked the BCs and ran again. interDyMFoam was crashed again after 20minutes running. I attached the log file, Could you give some guide for me? Thanks, hs// |
||
February 19, 2013, 06:39 |
Ship On Free Surface
|
#217 |
New Member
Join Date: Feb 2013
Posts: 6
Rep Power: 13 |
I am Sony, from Indonesia
now i am doing my Master thesis project by using CFD ansys 13 version. i intend to generate free surface simulation to predict ship resistance. the boundary condition in tutorial in term "free surface over bump" does not show good result. the simulation is always not finish. any body can help for explaining boundary condition used in free surface ship simulation in CFX? and may be also showing physical model that can be adopt? Thank you Regard Sony |
|
February 19, 2013, 06:52 |
Ship On Free Surface
|
#218 |
New Member
Join Date: Feb 2013
Posts: 6
Rep Power: 13 |
I am Sony, from Indonesia
now i am doing my master thesis project by using CFX ansys 13 version. i intend to generate free surface simulation to predict ship resistance. the boundary condition in tutorial in term "free surface over bump" does not show good result. the simulation is always not finish. any body can help for explaining boundary condition used in free surface ship simulation in CFX? and may be also showing physical model that can be adopt? Thank you Regard Sony |
|
February 19, 2013, 06:53 |
Ship On Free Surface
|
#219 |
New Member
Join Date: Feb 2013
Posts: 6
Rep Power: 13 |
I am Sony, from Indonesia
now i am doing my thesis project by using CFX ansys 13 version. i intend to generate free surface simulation to predict ship resistance. the boundary condition in tutorial in term "free surface over bump" does not show good result. the simulation is always not finish. any body can help for explaining boundary condition used in free surface ship simulation in CFX? and may be also showing physical model that can be adopt? Thank you Regard Sony |
|
February 19, 2013, 06:54 |
|
#220 |
Senior Member
Ralph Moolenaar
Join Date: Aug 2010
Location: 's-Hertogenbosch, the Netherlands
Posts: 120
Rep Power: 16 |
Hi Sony,
Right now you're posting in the OpenFOAM part of the forum...maybe you have more luck finding the answer that you need in the CFX forum? Could you also remove one (of your two) post above? Thanks! Ralph
__________________
CFD for marine applications? Go to http://www.marinecfd.com/ and join the OF Ship Hydromechanics Group: http://www.cfd-online.com/Forums/gro...mechanics.html |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Free-Surface Ship Flow - Boundary Conditions | James Date | CFX | 1 | February 19, 2013 06:42 |
ship free-surface analysis | Andrea Mercuri | Siemens | 0 | September 28, 2004 12:01 |
Free Surface Flow for Ship | sam | FLUENT | 6 | October 24, 2003 06:29 |
viscous free surface flow past a ship hull | lololo | Main CFD Forum | 0 | June 13, 2002 00:02 |
meshing for surface ship flow | boris | FLUENT | 0 | April 24, 2002 21:27 |