CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Problem with FloatingObject

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 22, 2012, 15:41
Default Problem with FloatingObject
  #1
Member
 
Pierre
Join Date: Sep 2010
Posts: 57
Rep Power: 6
Leech is on a distinguished road
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:rintStack(Foam::Ostream&) in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[0] #1 Foam::sigFpe::sigHandler(int)[1] #0 Foam::error:rintStack(Foam::Ostream&)[2] #0 Foam::error:rintStack(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:perator*<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:perator*<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:perator*<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 is offline   Reply With Quote

Old   January 23, 2012, 11:17
Default
  #2
Member
 
Pierre
Join Date: Sep 2010
Posts: 57
Rep Power: 6
Leech is on a distinguished road
Hi,

i found this here: OpenFoam.1.7.x floatingObject tutorial case & MULES:alpha1 greater than 1

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?

Last edited by Leech; January 23, 2012 at 12:16.
Leech is offline   Reply With Quote

Old   January 25, 2012, 11:48
Default
  #3
Member
 
Pierre
Join Date: Sep 2010
Posts: 57
Rep Power: 6
Leech is on a distinguished road
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?
Leech is offline   Reply With Quote

Old   January 26, 2012, 12:42
Default
  #4
New Member
 
Mark Beal
Join Date: Feb 2011
Posts: 24
Rep Power: 5
msbealo is on a distinguished road
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
__________________
Dynamic Fluid Design
www.dynamic-fluid-design.com
msbealo is offline   Reply With Quote

Old   January 26, 2012, 15:46
Default hello openfoamers
  #5
Member
 
Jignesh
Join Date: Aug 2011
Location: India
Posts: 68
Rep Power: 5
jignesh_thaker2007 is on a distinguished road
Send a message via Yahoo to jignesh_thaker2007
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
jignesh_thaker2007 is offline   Reply With Quote

Old   January 26, 2012, 17:05
Default
  #6
Member
 
Pierre
Join Date: Sep 2010
Posts: 57
Rep Power: 6
Leech is on a distinguished road
Quote:
Originally Posted by msbealo View Post
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...
Leech is offline   Reply With Quote

Old   January 30, 2012, 10:15
Default
  #7
New Member
 
Mark Beal
Join Date: Feb 2011
Posts: 24
Rep Power: 5
msbealo is on a distinguished road
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
__________________
Dynamic Fluid Design
www.dynamic-fluid-design.com
msbealo is offline   Reply With Quote

Old   January 30, 2012, 10:18
Default
  #8
New Member
 
Mark Beal
Join Date: Feb 2011
Posts: 24
Rep Power: 5
msbealo is on a distinguished road
Quote:
Originally Posted by jignesh_thaker2007 View Post
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
__________________
Dynamic Fluid Design
www.dynamic-fluid-design.com
msbealo is offline   Reply With Quote

Old   March 26, 2012, 18:01
Default disable time step adjustment fixes the floatingObject tutorial
  #9
nsf
Member
 
Nicolas Edh
Join Date: Mar 2010
Location: Uppsala, Sweden
Posts: 74
Rep Power: 6
nsf is on a distinguished road
Quote:
Originally Posted by Leech View Post
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
nsf is offline   Reply With Quote

Old   March 29, 2012, 05:44
Default
  #10
Member
 
Pierre
Join Date: Sep 2010
Posts: 57
Rep Power: 6
Leech is on a distinguished road
Quote:
Originally Posted by nsf View Post
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.
Leech is offline   Reply With Quote

Old   March 29, 2012, 16:24
Default
  #11
nsf
Member
 
Nicolas Edh
Join Date: Mar 2010
Location: Uppsala, Sweden
Posts: 74
Rep Power: 6
nsf is on a distinguished road
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!
nsf 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
Problem with interFoam; Wave/wiggle alpha1 behavior JonW OpenFOAM 3 February 23, 2013 21:41
UDF compiling problem Wouter Fluent UDF and Scheme Programming 6 June 6, 2012 05:43
natural convection problem for a CHT problem Se-Hee CFX 2 June 10, 2007 07:29
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 20:13
Is this problem well posed? Thomas P. Abraham Main CFD Forum 5 September 8, 1999 15:52


All times are GMT -4. The time now is 08:17.