[vince_44]

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

[colinB]

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

[vince_44]

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

[colinB]

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

[asoltoon]

Hello Eric,

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

Regards, Ali

[Ralph M]

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

[asoltoon]

Hi Ralph,

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

Regards, Ali

[parkh32]

Originally Posted by bouclette

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)

[Ralph M]

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

[parkh32]

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

[asoltoon]

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

[Ralph M]

@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

[parkh32]

Hi Ralph,

thanks for your reply, I will check BCs again.

hs//

[asoltoon]

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

[vince_44]

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

[parkh32]

Originally Posted by Ralph M

@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//

[sony anggara]

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]

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]

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

[Ralph M]

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