CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Flouting point exception twoPhaseEulerFoam - I set delta t (http://www.cfd-online.com/Forums/openfoam/88384-flouting-point-exception-twophaseeulerfoam-i-set-delta-t.html)

cristina87 May 16, 2011 09:28

Flouting point exception twoPhaseEulerFoam - I set delta t
 
Hi,

I receive the Flouting point exception when and I don't know why because I calculated delta t according the formula from tutorial and I set it in controlDict. I have the velocity on OX direction - from 0 to 12 so 13 is the d dimension and I have 20 cells on OX direction => delta x = 0.65. The velocity is (1 0 0) so I set my delta t to 0.65. The Co number I set to be max 0.2. I attached also my case.

Starting time loop

Courant Number mean: 0 max: -0
Max Ur Courant Number = -0
deltaT = 0.65
Time = 0.65

DILUPBiCG: Solving for alpha, Initial residual = 1, Final residual = 6.93889e-18, No Iterations 1
Dispersed phase volume fraction = 7.38636e+295 Min(alpha) = 0 Max(alpha) = 1.1375e+298
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam170/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam170/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::fv::gaussGrad<double>::grad(Foam::GeometricF ield<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam170/lib/linuxGccDPOpt/libfiniteVolume.so"
#4 Foam::fv::gaussGrad<double>::grad(Foam::GeometricF ield<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam170/lib/linuxGccDPOpt/libfiniteVolume.so"
#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/applications/bin/linuxGccDPOpt/twoPhaseEulerFoam"
#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/applications/bin/linuxGccDPOpt/twoPhaseEulerFoam"
#7 Foam::LimitedScheme<double, Foam::Limited01Limiter<Foam::limitedLinearLimiter< Foam::NVDTVD> >, Foam::limitFuncs::magSqr>::limiter(Foam::Geometric Field<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam170/lib/linuxGccDPOpt/libfiniteVolume.so"
#8 Foam::limitedSurfaceInterpolationScheme<double>::w eights(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam170/lib/linuxGccDPOpt/libfiniteVolume.so"
#9 Foam::surfaceInterpolationScheme<double>::interpol ate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam170/lib/linuxGccDPOpt/libfiniteVolume.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/linuxGccDPOpt/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/linuxGccDPOpt/libfiniteVolume.so"
#12
in "/opt/openfoam170/applications/bin/linuxGccDPOpt/twoPhaseEulerFoam"
#13
in "/opt/openfoam170/applications/bin/linuxGccDPOpt/twoPhaseEulerFoam"
#14 __libc_start_main in "/lib/libc.so.6"
#15
in "/opt/openfoam170/applications/bin/linuxGccDPOpt/twoPhaseEulerFoam"
Floating point exception

Cristina

freemankofi May 26, 2011 15:35

Cristina,
I think you might probably have an issue with your case setup itself and not the dt based on the error.

Check your BCs, transportproperties, etc to make sure that everything is alright.

BTW, NO test case file was attached to your posting.


All times are GMT -4. The time now is 14:03.