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

problem with new OF170 using floatingObject tutorial example

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 28, 2010, 14:14
Default problem with new OF170 using floatingObject tutorial example
  #1
Member
 
angel
Join Date: May 2009
Location: Spain
Posts: 46
Rep Power: 8
anmartin is on a distinguished road
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:rintStack(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"
#3 Foam::fv::gaussGrad<double>::grad(Foam::GeometricF ield<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam170/applications/bin/linux64GccDPOpt/interDyMFoam"
#4 Foam::fv::gaussGrad<double>::grad(Foam::GeometricF ield<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam170/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::tmp<Foam::GeometricField<Foam:uterProduct< 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:uterProduct< 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 is offline   Reply With Quote

Old   July 1, 2010, 01:09
Default loatingObject tutorial case & MULES:alpha1 greater than 1
  #2
Member
 
angel
Join Date: May 2009
Location: Spain
Posts: 46
Rep Power: 8
anmartin is on a distinguished road
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:rintStack(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"
#3 Foam::fv::gaussGrad<double>::grad(Foam::GeometricF ield<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/aml/OpenFOAM/openfoam170/applications/bin/linux64GccDPOpt/interDyMFoam"
#4 Foam::fv::gaussGrad<double>::grad(Foam::GeometricF ield<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/aml/OpenFOAM/openfoam170/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::tmp<Foam::GeometricField<Foam:uterProduct< 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:uterProduct< 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
anmartin 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 on Fluent Tutorial: Horizontal Film Boilig Feng FLUENT 2 April 13, 2013 05:34
newbie problem with cavity tutorial miki OpenFOAM Running, Solving & CFD 8 September 2, 2012 15:22
Help! Compiled UDF problem 4 Wave tank tutorial Shane FLUENT 1 September 3, 2010 02:32
Missing movingWall value in some OF-1.6.x tutorial cases. aunola OpenFOAM Bugs 1 September 13, 2009 17:55
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 19:13


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