CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   problem with new OF170 using floatingObject tutorial example (http://www.cfd-online.com/Forums/openfoam-solving/77595-problem-new-of170-using-floatingobject-tutorial-example.html)

 anmartin June 28, 2010 14:14

problem with new OF170 using floatingObject tutorial example

Hello,

Now I have installed the new OF170, when I try to run the floatingObject example from multiphase/interDyMFoam tutorial I obtain the following error:

Interface Courant Number mean: 0.00868469717822 max: 3.07048076829
Courant Number mean: 0.0214761703961 max: 3.07048076829
deltaT = 0.000368438312644
Time = 1.2731

Centre of mass: (0.373595449653 0.351882909424 0.482406400186)
Linear velocity: (-0.357961120598 -0.419740647472 -0.322074163028)
Angular velocity: (1.07106504184 -0.903173170713 -0.0110838047329)
GAMG: Solving for cellDisplacementx, Initial residual = 0.0056044966347, Final residual = 5.8395571122e-06, No Iterations 5
GAMG: Solving for cellDisplacementy, Initial residual = 0.00549399203972, Final residual = 5.71456061284e-06, No Iterations 5
GAMG: Solving for cellDisplacementz, Initial residual = 0.00997063228984, Final residual = 3.85256246304e-06, No Iterations 6
Execution time for mesh.update() = 0.7 s
time step continuity errors : sum local = 3.84344042594e-10, global = 2.49126712936e-11, cumulative = -0.000105905828067
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 2.81917055336e-06, No Iterations 7
time step continuity errors : sum local = 1.08353248665e-15, global = 2.55464254723e-16, cumulative = -0.000105905828067
MULES: Solving for alpha1
Liquid phase volume fraction = 0.533377598497 Min(alpha1) = -3.16898215629e+294 Max(alpha1) = 2.40723995432
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 in "/lib/libc.so.6"
#5 Foam::tmp<Foam::GeometricField<Foam::outerProduct< Foam::Vector<double>, double>::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<double>(Foam::GeometricField<doubl e, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/opt/openfoam170/lib/linux64GccDPOpt/libinterfaceProperties.so"
#6 Foam::tmp<Foam::GeometricField<Foam::outerProduct< Foam::Vector<double>, double>::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<double>(Foam::GeometricField<doubl e, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam170/lib/linux64GccDPOpt/libinterfaceProperties.so"
#7 Foam::LimitedScheme<double, Foam::vanLeerLimiter<Foam::NVDTVD>, Foam::limitFuncs::magSqr>::limiter(Foam::Geometric Field<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam170/lib/linux64GccDPOpt/libfiniteVolume.so"
#8 Foam::limitedSurfaceInterpolationScheme<double>::w eights(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam170/lib/linux64GccDPOpt/libinterfaceProperties.so"
#9 Foam::surfaceInterpolationScheme<double>::interpol ate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam170/lib/linux64GccDPOpt/libinterfaceProperties.so"
#10 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/openfoam170/lib/linux64GccDPOpt/libfiniteVolume.so"
#11 Foam::fv::gaussConvectionScheme<double>::flux(Foam ::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam170/lib/linux64GccDPOpt/libfiniteVolume.so"
#12
in "/opt/openfoam170/applications/bin/linux64GccDPOpt/interDyMFoam"
#13
in "/opt/openfoam170/applications/bin/linux64GccDPOpt/interDyMFoam"
#14
in "/opt/openfoam170/applications/bin/linux64GccDPOpt/interDyMFoam"
#15 __libc_start_main in "/lib/libc.so.6"
#16
in "/opt/openfoam170/applications/bin/linux64GccDPOpt/interDyMFoam"
Floating point exception

I will be very grateful if someone could help me

best regards

 anmartin July 1, 2010 01:09

loatingObject tutorial case & MULES:alpha1 greater than 1

Hello,
I still continues trying to run the floating tutorial case.

Anybody know what it is possible that alpha start to increase and produce simulation crash?

Ex:
Interface Courant Number mean: 0.00333864 max: 1.10327
Courant Number mean: 0.00852379 max: 1.10327
deltaT = 0.000195507
Time = 1.3364602

Centre of mass: (0.381972 0.360547 0.467263)
Linear velocity: (-0.165585 -0.204445 -0.203064)
Angular velocity: (-0.805702 0.711156 0.00404295)
GAMG: Solving for cellDisplacementx, Initial residual = 0.00244061, Final residual = 6.89207e-06, No Iterations 4
GAMG: Solving for cellDisplacementy, Initial residual = 0.00239952, Final residual = 6.77586e-06, No Iterations 4
GAMG: Solving for cellDisplacementz, Initial residual = 0.00378275, Final residual = 9.78091e-06, No Iterations 4
Execution time for mesh.update() = 1 s
time step continuity errors : sum local = 4.2295e-10, global = 4.07878e-11, cumulative = -9.58168e-05
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 9.20557e-06, No Iterations 8
time step continuity errors : sum local = 3.8935e-15, global = 3.69278e-18, cumulative = -9.58168e-05
MULES: Solving for alpha1
Liquid phase volume fraction = 0.536127 Min(alpha1) = -
4.94881e-15 Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.535844 Min(alpha1) = -0.159738 Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.535648 Min(alpha1) = -0.312417 Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.535528 Min(alpha1) = -0.421065 Max(alpha1) = 1

GAMG: Solving for p_rgh, Initial residual = 0.000486513, Final residual = 3.42267e-06, No Iterations 2
time step continuity errors : sum local = 9.36507e-08, global = 3.03639e-09, cumulative = -9.58138e-05
GAMGPCG: Solving for p_rgh, Initial residual = 0.000118452, Final residual = 2.18228e-09, No Iterations 4
time step continuity errors : sum local = 5.93232e-11, global = -4.55011e-12, cumulative = -9.58138e-05
DILUPBiCG: Solving for epsilon, Initial residual = 0.00973365, Final residual = 4.66459e-10, No Iterations 3
DILUPBiCG: Solving for k, Initial residual = 0.0139425, Final residual = 1.52868e-09, No Iterations 3
ExecutionTime = 200.59 s ClockTime = 201 s

Interface Courant Number mean: 0.00586782 max: 2.06191
Courant Number mean: 0.00939326 max: 2.36691
deltaT = 1.89614e-05
Time = 1.3364792

Centre of mass: (0.381969 0.360543 0.467259)
Linear velocity: (-0.165586 -0.204465 -0.203104)
Angular velocity: (-0.805398 0.710926 0.00409667)
GAMG: Solving for cellDisplacementx, Initial residual = 0.00408855, Final residual = 5.87695e-06, No Iterations 5
GAMG: Solving for cellDisplacementy, Initial residual = 0.00400199, Final residual = 5.74196e-06, No Iterations 5
GAMG: Solving for cellDisplacementz, Initial residual = 0.00562898, Final residual = 8.24578e-06, No Iterations 5
Execution time for mesh.update() = 1.07 s
time step continuity errors : sum local = 5.7535e-12, global = -4.41295e-13, cumulative = -9.58138e-05
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 6.28997e-06, No Iterations 9
time step continuity errors : sum local = 3.61893e-17, global = -3.48269e-18, cumulative = -9.58138e-05
MULES: Solving for alpha1
Liquid phase volume fraction = 0.533907 Min(alpha1) = -6.94436e+293 Max(alpha1) = 1.61898

#0 Foam::error::printStack(Foam::Ostream&) in "/home/aml/OpenFOAM/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/aml/OpenFOAM/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 in "/lib/libc.so.6"
#5 Foam::tmp<Foam::GeometricField<Foam::outerProduct< Foam::Vector<double>, double>::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<double>(Foam::GeometricField<doubl e, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/home/aml/OpenFOAM/openfoam170/lib/linux64GccDPOpt/libinterfaceProperties.so"
#6 Foam::tmp<Foam::GeometricField<Foam::outerProduct< Foam::Vector<double>, double>::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<double>(Foam::GeometricField<doubl e, Foam::fvPatchField, Foam::volMesh> const&) in "/home/aml/OpenFOAM/openfoam170/lib/linux64GccDPOpt/libinterfaceProperties.so"
#7 Foam::LimitedScheme<double, Foam::vanLeerLimiter<Foam::NVDTVD>, Foam::limitFuncs::magSqr>::limiter(Foam::Geometric Field<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/aml/OpenFOAM/openfoam170/lib/linux64GccDPOpt/libfiniteVolume.so"
#8 Foam::limitedSurfaceInterpolationScheme<double>::w eights(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/aml/OpenFOAM/openfoam170/lib/linux64GccDPOpt/libinterfaceProperties.so"
#9 Foam::surfaceInterpolationScheme<double>::interpol ate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/aml/OpenFOAM/openfoam170/lib/linux64GccDPOpt/libinterfaceProperties.so"
#10 Foam::fv::gaussConvectionScheme<double>::interpola te(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/aml/OpenFOAM/openfoam170/lib/linux64GccDPOpt/libfiniteVolume.so"
#11 Foam::fv::gaussConvectionScheme<double>::flux(Foam ::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/aml/OpenFOAM/openfoam170/lib/linux64GccDPOpt/libfiniteVolume.so"
#12
in "/home/aml/OpenFOAM/openfoam170/applications/bin/linux64GccDPOpt/interDyMFoam"
#13
in "/home/aml/OpenFOAM/openfoam170/applications/bin/linux64GccDPOpt/interDyMFoam"
#14
in "/home/aml/OpenFOAM/openfoam170/applications/bin/linux64GccDPOpt/interDyMFoam"
#15 __libc_start_main in "/lib/libc.so.6"
#16
in "/home/aml/OpenFOAM/openfoam170/applications/bin/linux64GccDPOpt/interDyMFoam"

best

 All times are GMT -4. The time now is 02:12.