CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Problem with FloatingObject (https://www.cfd-online.com/Forums/openfoam-solving/96415-problem-floatingobject.html)

Leech January 22, 2012 14:41
Problem with FloatingObject
 
Hi,

i used the floatingObject tutorial as basic case and changed some stuff. (Other box geometrie, no water pillar crashing, but waves by groovy BC, etc.)
blockMesh, topoSet, subsetMesh and setFields run without any errors and they create the geometrie i wanted. But when i run interDyMFoam this is what happens:

Quote:
Parallel processing using SYSTEMOPENMPI with 3 processors
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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, global = 0, cumulative = 0
GAMGPCG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0, global = 0, cumulative = 0
Courant Number mean: 0 max: 0

Starting time loop

Interface Courant Number mean: 0 max: 0
Courant Number mean: 0 max: 0
deltaT = 1.19990400768e-05
Time = 1.199904007679385782754455380683111798134632408618 927001953125e-05

Centre of mass: (0.5 0.5 0.5)
Linear velocity: (0 0 -0.0061535137797)
Angular velocity: (-1.23699579461 1.23699579461 -3.03140649339e-16)
GAMG: Solving for cellDisplacementx, Initial residual = 0, Final residual = 0, No Iterations 0
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.88 s
time step continuity errors : sum local = 0, global = 0, cumulative = 0
GAMGPCG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0, global = 0, cumulative = 0
MULES: Solving for alpha1
Liquid phase volume fraction = 0.499999999998 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.499999999998 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.499999999998 Min(alpha1) = 0 Max(alpha1) = 1
GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 0.0023053244698, No Iterations 1
time step continuity errors : sum local = 3.17665412683e-08, global = -9.92537010983e-10, cumulative = -9.92537010983e-10
GAMGPCG: Solving for p_rgh, Initial residual = 0.000819301605459, Final residual = 2.66295853907e-09, No Iterations 6
time step continuity errors : sum local = 1.04743489104e-13, global = -1.25880765395e-14, cumulative = -9.9254959906e-10
smoothSolver: Solving for epsilon, Initial residual = 6.93749749838e-05, Final residual = 3.93397656793e-09, No Iterations 1
smoothSolver: Solving for k, Initial residual = 1, Final residual = 5.646373754e-09, No Iterations 2
ExecutionTime = 23.27 s ClockTime = 23 s

Interface Courant Number mean: 0 max: 0
Courant Number mean: 1.46703267158e-09 max: 1.59921765751e-06
deltaT = 1.43971203686e-05
Time = 2.639616044539331232037827901226734184092492796480 655670166015625e-05

Centre of mass: (0.5 0.5 0.499999805108)
Linear velocity: (-5.64965799146e-07 -2.75751900641e-07 -0.0127947704547)
Angular velocity: (-2.55623059431 2.55623845484 3.58503208285e-05)
GAMG: Solving for cellDisplacementx, Initial residual = 1, Final residual = 4.94059419397e-06, No Iterations 10
GAMG: Solving for cellDisplacementy, Initial residual = 1, Final residual = 4.94059419395e-06, No Iterations 10
GAMG: Solving for cellDisplacementz, Initial residual = 1, Final residual = 4.82235690693e-06, No Iterations 10
Execution time for mesh.update() = 8.34 s
time step continuity errors : sum local = 1.2567710507e-13, global = -1.51038793094e-14, cumulative = -9.92564702939e-10
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 3.1516056274e-06, No Iterations 8
time step continuity errors : sum local = 3.9608467269e-19, global = -7.56805235012e-20, cumulative = -9.92564703015e-10
MULES: Solving for alpha1
Liquid phase volume fraction = 0.499992190839 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.499995996431 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.499999802097 Min(alpha1) = 0 Max(alpha1) = 1
GAMG: Solving for p_rgh, Initial residual = 0.0107100284668, Final residual = 6.25337728779e-05, No Iterations 8
time step continuity errors : sum local = 4.02883959303e-09, global = -2.09767644717e-09, cumulative = -3.09024115018e-09
GAMGPCG: Solving for p_rgh, Initial residual = 3.51711714171e-05, Final residual = 1.91859209088e-09, No Iterations 4
time step continuity errors : sum local = 5.54714114581e-12, global = 7.93573595888e-13, cumulative = -3.08944757659e-09
smoothSolver: Solving for epsilon, Initial residual = 0.00016822924594, Final residual = 1.59679720409e-07, No Iterations 1
smoothSolver: Solving for k, Initial residual = 0.945174262496, Final residual = 5.39143318514e-09, No Iterations 3
ExecutionTime = 54.84 s ClockTime = 55 s

Interface Courant Number mean: 0 max: 0
Courant Number mean: 0.00097852148895 max: 0.0132106618024
deltaT = 1.7275549307e-05
Time = 4.367170975238391560308351402319715361954877153038 97857666015625e-05

Centre of mass: (0.499999999979 0.49999999999 0.499999599455)
Linear velocity: (-6.56540194232 -6.57855838451 76.3356295205)
Angular velocity: (16214.9562388 -16213.5539856 -1.60771114785)
GAMG: Solving for cellDisplacementx, Initial residual = 0.320308797821, Final residual = 8.14063828051e-06, No Iterations 8
GAMG: Solving for cellDisplacementy, Initial residual = 0.320332721129, Final residual = 8.14124007992e-06, No Iterations 8
GAMG: Solving for cellDisplacementz, Initial residual = 0.320358747561, Final residual = 8.1156648521e-06, No Iterations 8
Execution time for mesh.update() = 7.28 s
time step continuity errors : sum local = 6.65618595413e-12, global = 9.52233462908e-13, cumulative = -3.08849534313e-09
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 3.60728334458e-06, No Iterations 9
time step continuity errors : sum local = 2.40107292951e-17, global = 8.12773516175e-18, cumulative = -3.088495335e-09
MULES: Solving for alpha1
Liquid phase volume fraction = 0.499991767373 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.499995800386 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.499999833481 Min(alpha1) = 0 Max(alpha1) = 1
GAMG: Solving for p_rgh, Initial residual = 0.000364083509199, Final residual = 2.65613705406e-06, No Iterations 7
time step continuity errors : sum local = 1.12988147847e-08, global = 5.44795910026e-09, cumulative = 2.35946376526e-09
GAMGPCG: Solving for p_rgh, Initial residual = 0.000665537320881, Final residual = 2.69318919171e-09, No Iterations 6
time step continuity errors : sum local = 1.54092547837e-13, global = -1.34854857419e-14, cumulative = 2.35945027977e-09
smoothSolver: Solving for epsilon, Initial residual = 0.000190031852915, Final residual = 1.93766873312e-07, No Iterations 1
smoothSolver: Solving for k, Initial residual = 0.363429700421, Final residual = 9.19435886496e-07, No Iterations 2
ExecutionTime = 88.26 s ClockTime = 88 s

Interface Courant Number mean: 1.10514975191e-06 max: 0.0108393115476
Courant Number mean: 0.00101659386078 max: 0.0137133881986
deltaT = 2.07292261075e-05
Time = 6.440093585985957164684323483783146002679131925106 048583984375e-05

Centre of mass: (0.499700600756 0.499700000784 0.503480997046)
Linear velocity: (-15.1047541567 -15.1316708777 166.219646585)
Angular velocity: (35428.0375917 -35356.1342217 -36.814993058)
GAMG: Solving for cellDisplacementx, Initial residual = 0.999973212739, Final residual = 6.56364965277e-06, No Iterations 10
GAMG: Solving for cellDisplacementy, Initial residual = 0.999973329662, Final residual = 6.56385657994e-06, No Iterations 10
GAMG: Solving for cellDisplacementz, Initial residual = 0.999919316933, Final residual = 6.46238168198e-06, No Iterations 10
Execution time for mesh.update() = 8.37 s
time step continuity errors : sum local = 1.8229498367e-13, global = -1.59536359068e-14, cumulative = 2.35943432614e-09
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 1.00000544823, No Iterations 101
time step continuity errors : sum local = 1.82295976855e-13, global = -1.59532959623e-14, cumulative = 2.35941837284e-09
MULES: Solving for alpha1
Liquid phase volume fraction = 0.558964333137 Min(alpha1) = -2.85833788947e+297 Max(alpha1) = 8.00055428652e+265
[0] #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[0] #1 Foam::sigFpe::sigHandler(int)[1] #0 Foam::error::printStack(Foam::Ostream&)[2] #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[0] #2 Uninterpreted:
[0] #3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[1] #1 in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[2] #1 Foam::sigFpe::sigHandler(int)Foam::sigFpe::sigHand ler(int) in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[0] #4 in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[1] #2 Uninterpreted:
[1] #3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[2] #2 Uninterpreted:
[2] #3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&)void Foam::multiply<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[1] #4 in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[2] #4 in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/interDyMFoam"
[0] #5 void Foam::multiply<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&)void Foam::multiply<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&)Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::operator*<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/interDyMFoam"
[1] #5 in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/interDyMFoam"
[2] #5 in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/interDyMFoam"
[0] #6 Foam::fv::gaussConvectionScheme<double>::flux(Foam ::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) constFoam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::operator*<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&)Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::operator*<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/interDyMFoam"
[2] #6 Foam::fv::gaussConvectionScheme<double>::flux(Foam ::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/interDyMFoam"
[1] #6 in "/opt/openfoam200/platforms/linuxGccDPOpt/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 "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
[2] #7
in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
[1] #7

[0] in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/interDyMFoam"
[0] #8 [2] in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/interDyMFoam"
[2] #8
[1] in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/interDyMFoam"
[1] #8
[0] in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/interDyMFoam"
[0] #9 __libc_start_main
[2] in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/interDyMFoam"
[2] #9 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
[0] #10 [1] in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/interDyMFoam"
[1] #9 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
[2] #10 in "/lib/i386-linux-gnu/libc.so.6"
[1] #10


[0] in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/interDyMFoam"
[mstrlnx:18577] *** Process received signal ***
[mstrlnx:18577] Signal: Floating point exception (8)
[mstrlnx:18577] Signal code: (-6)
[mstrlnx:18577] Failing at address: 0x4891
[mstrlnx:18577] [ 0] [0xb77f140c]
[mstrlnx:18577] [ 1] [0xb77f1424]
[mstrlnx:18577] [ 2] /lib/i386-linux-gnu/libc.so.6(gsignal+0x4f) [0xb57bec8f]
[mstrlnx:18577] [ 3] /opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6sigFpe10sigHandlerEi+0x61) [0xb5e38551]
[mstrlnx:18577] [ 4] [0xb77f1400]
[mstrlnx:18577] [ 5] /opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZN4Foam8multiplyERNS_5FieldIdEERKN S_5UListIdEES6_+0x20) [0xb5d79f10]
[mstrlnx:18577] [ 6] interDyMFoam(_ZN4Foam8multiplyINS_13fvsPatchFieldE NS_11surfaceMeshEEEvRNS_14GeometricFieldIdT_T0_EER KS6_S9_+0x46) [0x80a8476]
[mstrlnx:18577] [ 7] interDyMFoam(_ZN4FoammlINS_13fvsPatchFieldENS_11su rfaceMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEE RKS7_RKS8_+0x327) [0x80cd4e7]
[mstrlnx:18577] [ 8] /opt/openfoam200/platforms/linuxGccDPOpt/lib/libfiniteVolume.so(_ZNK4Foam2fv21gaussConvectionSc hemeIdE4fluxERKNS_14GeometricFieldIdNS_13fvsPatchF ieldENS_11surfaceMeshEEERKNS3_IdNS_12fvPatchFieldE NS_7volMeshEEE+0x56) [0xb6a92476]
[mstrlnx:18577] [ 9] interDyMFoam() [0x8090498]
[mstrlnx:18577] [10] interDyMFoam() [0x806d365]
[mstrlnx:18577] [11] /lib/i386-linux-gnu/libc.so.6(__libc_start_main+0xf3) [0xb57aa113]
[mstrlnx:18577] [12] interDyMFoam() [0x8065961]
[mstrlnx:18577] *** End of error message ***
[1] in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/interDyMFoam"
[mstrlnx:18578] *** Process received signal ***
[mstrlnx:18578] Signal: Floating point exception (8)
[mstrlnx:18578] Signal code: (-6)
[mstrlnx:18578] Failing at address: 0x4892
[mstrlnx:18578] [ 0] [0xb775440c]
[mstrlnx:18578] [ 1] [0xb7754424]
[mstrlnx:18578] [ 2] /lib/i386-linux-gnu/libc.so.6(gsignal+0x4f) [0xb5721c8f]
[mstrlnx:18578] [ 3] /opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6sigFpe10sigHandlerEi+0x61) [0xb5d9b551]
[mstrlnx:18578] [ 4] [0xb7754400]
[mstrlnx:18578] [ 5] /opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZN4Foam8multiplyERNS_5FieldIdEERKN S_5UListIdEES6_+0x20) [0xb5cdcf10]
[mstrlnx:18578] [ 6] interDyMFoam(_ZN4Foam8multiplyINS_13fvsPatchFieldE NS_11surfaceMeshEEEvRNS_14GeometricFieldIdT_T0_EER KS6_S9_+0x46) [0x80a8476]
[mstrlnx:18578] [ 7] interDyMFoam(_ZN4FoammlINS_13fvsPatchFieldENS_11su rfaceMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEE RKS7_RKS8_+0x327) [0x80cd4e7]
[mstrlnx:18578] [ 8] /opt/openfoam200/platforms/linuxGccDPOpt/lib/libfiniteVolume.so(_ZNK4Foam2fv21gaussConvectionSc hemeIdE4fluxERKNS_14GeometricFieldIdNS_13fvsPatchF ieldENS_11surfaceMeshEEERKNS3_IdNS_12fvPatchFieldE NS_7volMeshEEE+0x56) [0xb69f5476]
[mstrlnx:18578] [ 9] interDyMFoam() [0x8090498]
[mstrlnx:18578] [10] interDyMFoam() [0x806d365]
[mstrlnx:18578] [11] /lib/i386-linux-gnu/libc.so.6(__libc_start_main+0xf3) [0xb570d113]
[mstrlnx:18578] [12] interDyMFoam() [0x8065961]
[mstrlnx:18578] *** End of error message ***
[2] in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/interDyMFoam"
[mstrlnx:18579] *** Process received signal ***
[mstrlnx:18579] Signal: Floating point exception (8)
[mstrlnx:18579] Signal code: (-6)
[mstrlnx:18579] Failing at address: 0x4893
[mstrlnx:18579] [ 0] [0xb77d840c]
[mstrlnx:18579] [ 1] [0xb77d8424]
[mstrlnx:18579] [ 2] /lib/i386-linux-gnu/libc.so.6(gsignal+0x4f) [0xb57a5c8f]
[mstrlnx:18579] [ 3] /opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6sigFpe10sigHandlerEi+0x61) [0xb5e1f551]
[mstrlnx:18579] [ 4] [0xb77d8400]
[mstrlnx:18579] [ 5] /opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZN4Foam8multiplyERNS_5FieldIdEERKN S_5UListIdEES6_+0x20) [0xb5d60f10]
[mstrlnx:18579] [ 6] interDyMFoam(_ZN4Foam8multiplyINS_13fvsPatchFieldE NS_11surfaceMeshEEEvRNS_14GeometricFieldIdT_T0_EER KS6_S9_+0x46) [0x80a8476]
[mstrlnx:18579] [ 7] interDyMFoam(_ZN4FoammlINS_13fvsPatchFieldENS_11su rfaceMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEE RKS7_RKS8_+0x327) [0x80cd4e7]
[mstrlnx:18579] [ 8] /opt/openfoam200/platforms/linuxGccDPOpt/lib/libfiniteVolume.so(_ZNK4Foam2fv21gaussConvectionSc hemeIdE4fluxERKNS_14GeometricFieldIdNS_13fvsPatchF ieldENS_11surfaceMeshEEERKNS3_IdNS_12fvPatchFieldE NS_7volMeshEEE+0x56) [0xb6a79476]
[mstrlnx:18579] [ 9] interDyMFoam() [0x8090498]
[mstrlnx:18579] [10] interDyMFoam() [0x806d365]
[mstrlnx:18579] [11] /lib/i386-linux-gnu/libc.so.6(__libc_start_main+0xf3) [0xb5791113]
[mstrlnx:18579] [12] interDyMFoam() [0x8065961]
[mstrlnx:18579] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 18577 on node mstrlnx exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------
sorry its a bit long :) (By the way the groovy BC is not active, there are no waves, just a water surface with the floatingThing)
It crashes and writes floating point exception. So i tried the refine the mesh (up to 2 million cells) and to reduce the time step (down to 0.0001 s). Both didnt helped.
What i can also see in the log is that the velocity of the floatingObject goes very high! So i guess its just falling through the water.. (not so much of a floatingObject ;) )
Any idea why the object isnt bound to the surface? What further information you need? Shall i post some of the case files here?

Thanks a lot!
Have a nice evening
Leech

Leech January 23, 2012 10:17
Hi,

i found this here: http://www.cfd-online.com/Forums/ope...-than-1-a.html

There they found that there is a impulse in prgh crashing the case. I tested my floating obejct case (the original one) and it crashed. Then i tested the floating object case on my new setup machine (Foam 2.1.0) and there it runs. Then i used the 2.1.0 files and again changed the mesh. But still it keeps crashing after the velocity goes extremly high and then alpha1 goes high..

Any suggestions?

Thanks!


Second question: In the linked Thread they solved the problem by fixing the body. I want to fix the body to move in x and y direction. How can i do this?

Leech January 25, 2012 10:48
I tried now the following:


Quote:
floatingObject
{
type sixDoFRigidBodyDisplacement;
centreOfMass (4 4 3.5);
momentOfInertia (0.1 0.1 0.1);
mass 10;
rhoInf 1; // needed only for solvers solving for kinematic pressure
report on;
value uniform (0 0 0);
velocity (0 0 0);
acceleration (0 0 0);

restraints
{
verticalSpring1
{
sixDoFRigidBodyMotionRestraint linearSpring;

linearSpringCoeffs
{
anchor (4 4 0);
refAttachmentPt (4 4 3.5);
stiffness 4000;
damping 2;
restLength 4;
}
}
}

constraints
{
maxIterations 500;

fixedAxis1
{
sixDoFRigidBodyMotionConstraint fixedAxis;
tolerance 1e-06;
relaxationFactor 0.7;
fixedAxisCoeffs
{
axis ( 0 0 1 );
}
}

fixedAxis2
{
sixDoFRigidBodyMotionConstraint fixedAxis;
tolerance 1e-06;
relaxationFactor 0.7;
fixedAxisCoeffs
{
axis ( 0 1 0 );
}
}

fixedAxis3
{
sixDoFRigidBodyMotionConstraint fixedAxis;
tolerance 1e-06;
relaxationFactor 0.7;
fixedAxisCoeffs
{
axis ( 1 0 0 );
}
}
}
}

But it still crashes. The body still begins to move along the z-axis?
I got 4 constraints? How can it move?

msbealo January 26, 2012 11:42
Leech,

I suggest running the interFoam solver first until it reached a steady state, and then start using the interDymFoam solver. In the past I've had experiences of large forces in the z direction at the beginning of the simulation which lead to large velocities and displacement of the object. The object might not be able to recover from these initial kicks.

Regards,

Mark Beal

jignesh_thaker2007 January 26, 2012 14:46
hello openfoamers
 
I tried to solve the problem of 3-D pipe flow and i want to use the K-Epsilon turbulence model. it is two phase flow problem. water and air is fed at inlet with different velocity and inside pipe it mix with each other and coming out to the another hand. so i want to plot different regimes like slug, plug, stratified, wavy etc. so plz tel me what is the inlet and outlet boundary condition for K, epsilon, alpha, pressure and velocity. and also tel me is there any roll of nut and nutilda? and tel me which solver i used rasinterFoam, interFoam or interDyMFoam?

plz help me i am so confused.

i am waiting for reply

jignesh thaker

Leech January 26, 2012 16:05
Quote:
Originally Posted by msbealo (Post 341337)
Leech,

I suggest running the interFoam solver first until it reached a steady state, and then start using the interDymFoam solver. In the past I've had experiences of large forces in the z direction at the beginning of the simulation which lead to large velocities and displacement of the object. The object might not be able to recover from these initial kicks.

Regards,

Mark Beal
Thanks Mark. I'll try that out tomorrow and will tell if it works.
Anyway i go some strange problems: Just to be shure i copied the floatingObject tutorial fresh and had it running on 4 cores. And after 1,8s (case time) it crashes. But when i run it single core it goes through without crashing. I am getting really confused about this floating object thing. Sometimes it works, sometimes not, sometimes on multicore, sometimes not...

msbealo January 30, 2012 09:15
Leech,

Yeah.. that's an annoying problem and I'm not sure of the solution. To me, there appears to be a problem with the interDymFoam solver when used on multiple cores. If you find a solution (or the cause) then I'd be very interested to know it. I'm not working on any dynamic mesh problems right now so I can't offer a solution.

Regards,

Mark

msbealo January 30, 2012 09:18
Quote:
Originally Posted by jignesh_thaker2007 (Post 341362)
I tried to solve the problem of 3-D pipe flow..... so plz tel me what is the inlet and outlet boundary condition for K, epsilon, alpha, pressure and velocity..... and tel me which solver i used rasinterFoam, interFoam or interDyMFoam?

jignesh thaker

Jignesh,

To start I'd work with some pre-existing tutorials and try to adapt them to your case. You might need to combine a number of different cases to get what you need. I suggest you start a new thread and ask again.

Regards,

Mark

nsf March 26, 2012 17:01
disable time step adjustment fixes the floatingObject tutorial
 
Quote:
Originally Posted by Leech (Post 341376)
Thanks Mark. I'll try that out tomorrow and will tell if it works.
Anyway i go some strange problems: Just to be shure i copied the floatingObject tutorial fresh and had it running on 4 cores. And after 1,8s (case time) it crashes. But when i run it single core it goes through without crashing. I am getting really confused about this floating object thing. Sometimes it works, sometimes not, sometimes on multicore, sometimes not...
Hello,

I've had similar problem twice now. The first time was with a separate case with interFoam. I experienced unexplainable crashes.

The second time was when I tried the floatingObject tutorial. The simulation crashes after 0.73s. Both parallel and serial.

In both cases I tried different fvSchemes but it didn't help.

What solved the issue for me was to disable time step adjustment. In the floatingObject tutorial I chose the standard time step of 0.01s and disabled time adjustment. After that the simulations completed without issue.

I should perhaps say I'm running 2.1.x but there hasn't been any relevant change to the tutorial files since the first release of 2.0.x.

Best

Nicolas

Leech March 29, 2012 04:44
Quote:
Originally Posted by nsf (Post 351596)
Hello,

I've had similar problem twice now. The first time was with a separate case with interFoam. I experienced unexplainable crashes.

The second time was when I tried the floatingObject tutorial. The simulation crashes after 0.73s. Both parallel and serial.

In both cases I tried different fvSchemes but it didn't help.

What solved the issue for me was to disable time step adjustment. In the floatingObject tutorial I chose the standard time step of 0.01s and disabled time adjustment. After that the simulations completed without issue.

I should perhaps say I'm running 2.1.x but there hasn't been any relevant change to the tutorial files since the first release of 2.0.x.

Best

Nicolas

Hi!

Thanks for that tip! It seems to help at least to get it running a bit longer.
If you havnt seen it yet, i did a bug report on the crashes concerning interDyMFoam. Link:
http://www.openfoam.org/mantisbt/view.php?id=417

Ill continue testing.

nsf March 29, 2012 15:24
Hi,

You could also try to select a time step that insures that the interface courant number is below 0.3. This is suggested inte this report

[URL="http://publications.lib.chalmers.se/cpl/record/index.xsql?pubid=100901"]http://publications.lib.chalmers.se/cpl/record/index.xsql?pubid=100901[/URL

If you google the author you can find the report as a pdf.

Regarding the floatingObject tutorial I could run it in parallel on ubuntu 10.04, with adjustable run time. This was on a three year old amd phenom II.

The floatingObject crashed on ubuntu 11.10 on a intel i5 sandy bridge. I've just recompiled on the sandy bridge machine and am going to try to run in in serial once more.
(the other case was to big to run in serial and was crashing on SLES10, on 2.0.x)

If I still experience crashes I'll post back in your bug report.

Good luck with your simulation!


All times are GMT -4. The time now is 22:51.