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

Free Surface Ship Flow

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

Like Tree18Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   January 18, 2013, 10:18
Default
  #201
Member
 
vincent
Join Date: Apr 2011
Posts: 41
Rep Power: 7
vince_44 is on a distinguished road
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
Attached Files
File Type: doc checkMesh.doc (17.5 KB, 14 views)
File Type: doc SnappyHexMesh.doc (41.0 KB, 13 views)
vince_44 is offline   Reply With Quote

Old   January 18, 2013, 10:31
Default
  #202
Senior Member
 
Join Date: Aug 2010
Location: Groningen, The Netherlands
Posts: 216
Rep Power: 11
colinB is on a distinguished road
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
Ralph M likes this.

Last edited by colinB; January 18, 2013 at 10:47.
colinB is offline   Reply With Quote

Old   January 18, 2013, 10:53
Default
  #203
Member
 
vincent
Join Date: Apr 2011
Posts: 41
Rep Power: 7
vince_44 is on a distinguished road
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
vince_44 is offline   Reply With Quote

Old   January 21, 2013, 04:42
Default
  #204
Senior Member
 
Join Date: Aug 2010
Location: Groningen, The Netherlands
Posts: 216
Rep Power: 11
colinB is on a distinguished road
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
colinB is offline   Reply With Quote

Old   January 27, 2013, 11:49
Default
  #205
Member
 
Ali Khalifesoltani
Join Date: Mar 2011
Location: Esfahan, Iran
Posts: 52
Rep Power: 7
asoltoon is on a distinguished road
Hello Eric,

I think the your link has some problems. It doesn't open. Could you possibly post the right link?

Regards, Ali
asoltoon is offline   Reply With Quote

Old   January 27, 2013, 15:40
Default
  #206
Senior Member
 
Ralph Moolenaar
Join Date: Aug 2010
Location: 's-Hertogenbosch, the Netherlands
Posts: 120
Rep Power: 7
Ralph M is on a distinguished road
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
Ralph M is offline   Reply With Quote

Old   January 28, 2013, 02:19
Default
  #207
Member
 
Ali Khalifesoltani
Join Date: Mar 2011
Location: Esfahan, Iran
Posts: 52
Rep Power: 7
asoltoon is on a distinguished road
Hi Ralph,

I will take a look at the tutorial.
thanks for your reply.

Regards, Ali
asoltoon is offline   Reply With Quote

Old   February 3, 2013, 04:03
Default
  #208
New Member
 
hs
Join Date: Mar 2012
Posts: 22
Rep Power: 6
parkh32 is on a distinguished road
Quote:
Originally Posted by bouclette View Post
I take it from the log file.
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)
parkh32 is offline   Reply With Quote

Old   February 3, 2013, 07:18
Default
  #209
Senior Member
 
Ralph Moolenaar
Join Date: Aug 2010
Location: 's-Hertogenbosch, the Netherlands
Posts: 120
Rep Power: 7
Ralph M is on a distinguished road
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
Ralph M is offline   Reply With Quote

Old   February 3, 2013, 16:26
Default
  #210
New Member
 
hs
Join Date: Mar 2012
Posts: 22
Rep Power: 6
parkh32 is on a distinguished road
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
parkh32 is offline   Reply With Quote

Old   February 4, 2013, 11:31
Default
  #211
Member
 
Ali Khalifesoltani
Join Date: Mar 2011
Location: Esfahan, Iran
Posts: 52
Rep Power: 7
asoltoon is on a distinguished road
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
asoltoon is offline   Reply With Quote

Old   February 4, 2013, 15:48
Default
  #212
Senior Member
 
Ralph Moolenaar
Join Date: Aug 2010
Location: 's-Hertogenbosch, the Netherlands
Posts: 120
Rep Power: 7
Ralph M is on a distinguished road
@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
Ralph M is offline   Reply With Quote

Old   February 4, 2013, 16:10
Default
  #213
New Member
 
hs
Join Date: Mar 2012
Posts: 22
Rep Power: 6
parkh32 is on a distinguished road
Hi Ralph,

thanks for your reply, I will check BCs again.

hs//
parkh32 is offline   Reply With Quote

Old   February 5, 2013, 10:46
Default
  #214
Member
 
Ali Khalifesoltani
Join Date: Mar 2011
Location: Esfahan, Iran
Posts: 52
Rep Power: 7
asoltoon is on a distinguished road
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
asoltoon is offline   Reply With Quote

Old   February 6, 2013, 07:37
Default
  #215
Member
 
vincent
Join Date: Apr 2011
Posts: 41
Rep Power: 7
vince_44 is on a distinguished road
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
vince_44 is offline   Reply With Quote

Old   February 8, 2013, 19:30
Default
  #216
New Member
 
hs
Join Date: Mar 2012
Posts: 22
Rep Power: 6
parkh32 is on a distinguished road
Quote:
Originally Posted by Ralph M View Post
@Hs: did you also checked your BC's? When checkmesh is okay the error must be somewhere else....

Cheers,

Ralph
Hi Ralph,

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//
parkh32 is offline   Reply With Quote

Old   February 19, 2013, 06:39
Default Ship On Free Surface
  #217
New Member
 
Join Date: Feb 2013
Posts: 6
Rep Power: 5
sony anggara is on a distinguished road
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
sony anggara is offline   Reply With Quote

Old   February 19, 2013, 06:52
Default Ship On Free Surface
  #218
New Member
 
Join Date: Feb 2013
Posts: 6
Rep Power: 5
sony anggara is on a distinguished road
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
sony anggara is offline   Reply With Quote

Old   February 19, 2013, 06:53
Default Ship On Free Surface
  #219
New Member
 
Join Date: Feb 2013
Posts: 6
Rep Power: 5
sony anggara is on a distinguished road
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
sony anggara is offline   Reply With Quote

Old   February 19, 2013, 06:54
Default
  #220
Senior Member
 
Ralph Moolenaar
Join Date: Aug 2010
Location: 's-Hertogenbosch, the Netherlands
Posts: 120
Rep Power: 7
Ralph M is on a distinguished road
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
Ralph M 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
Free-Surface Ship Flow - Boundary Conditions James Date CFX 1 February 19, 2013 06:42
ship free-surface analysis Andrea Mercuri CD-adapco 0 September 28, 2004 11:01
Free Surface Flow for Ship sam FLUENT 6 October 24, 2003 05:29
viscous free surface flow past a ship hull lololo Main CFD Forum 0 June 12, 2002 23:02
meshing for surface ship flow boris FLUENT 0 April 24, 2002 20:27


All times are GMT -4. The time now is 09:40.