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

icoFoam - Floating point exception (core dumped)

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

Like Tree2Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   July 18, 2013, 05:44
Default
  #21
Member
 
Amin
Join Date: May 2013
Posts: 76
Rep Power: 5
Mirage12 is on a distinguished road
i do not really want to a final steady-state solution, that is why i am trying a transient solver .

but nimasam suggested, that i can use simpleFoam for for laminar flow, that is why i try to add the temperature in the simpleFoam solver.
Mirage12 is offline   Reply With Quote

Old   August 4, 2013, 11:49
Default
  #22
New Member
 
Gabriel Boucher
Join Date: Jul 2013
Posts: 23
Rep Power: 5
Alfalfa is on a distinguished road
Hello everyone, I also have a problem with my courant number.

In order to practice with OpenFOAM (since I'm a beginner), I tried to create a simple nozzle with blockMesh. To run the case, I use the elbow simulation with icoFoam. I changed the meshing to use the one I created and run blockMesh. but when I run the simulation, I get courant number of the magintue e289 on the very first step... I tried a coarse meshing, but nothing changed, I tried really small time step, still nothing and now I don't know what to do .

Here is my meshing:

Code:
convertToMeters 0.1;

vertices
(
    (-25 0 0.5)        //0
    (0 -10 0.5)        //1
    (0 10 0.5)        //2
    (-25 5 0.5)        //3
    (-25 0 -0.5)    //4
    (0 -10 -0.5)    //5
    (0 10 -0.5)        //6
    (-25 5 -0.5)    //7


    
);

blocks
(
    hex (0 1 2 3 4 5 6 7) (25 25 1) simpleGrading (1 1 1)

);

edges
(
);

boundary
(
    inlet
    {
        type patch;
        faces
        (
            (0 3 7 4)
        );
    }
    outlet
    {
        type patch;
        faces
        (
        (1 2 6 5)
        );
    }
    upperWall
    {
        type wall;
        faces
        (
            (3 2 6 7)
        );
    }
    lowerWall
    {
        type wall;
        faces
        (
            (0 1 5 4)
        );
    }
    frontAndBack
    {
        type empty;
        faces
        (
            (0 1 2 3)
            (4 5 6 7)
        );
    }
);

mergePatchPairs
(
);
Here's the controlDict :

Code:
application     icoFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         0.1;

deltaT          0.001;

writeControl    timeStep;

writeInterval   20;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;
I run OpenFOAM 2.2.0 on Ubuntu 12.04 LTS.

Any help would be appreciated. I try to learn OpenFoam by modifying the tutorial and this is my first (pathetic) attempt.

Thanks
Alfalfa is offline   Reply With Quote

Old   August 4, 2013, 12:04
Default
  #23
New Member
 
Gabriel Boucher
Join Date: Jul 2013
Posts: 23
Rep Power: 5
Alfalfa is on a distinguished road
Oh nevermind, I found my error using checkMesh . My mesh orientation was not good.
Alfalfa is offline   Reply With Quote

Old   April 16, 2016, 05:58
Default Floating point exception(core dumped)
  #24
New Member
 
sagar
Join Date: Apr 2016
Posts: 1
Rep Power: 0
shaggadude is on a distinguished road
Hello Foamers,
I am new to openfoam.I have created my own solver adding temperature equation (only advection-diffusion with velocity) to the icoFoam solver. to be specific ,code is
solve
( fvm::div(phi,T)
- fvm::laplacian(DT,T)
);
I am getting this error

Courant Number mean: 0.0248669 max: 2.20364
smoothSolver: Solving for Ux, Initial residual = 0.202777, Final residual = 1.21499e+270, No Iterations 1000
smoothSolver: Solving for Uy, Initial residual = 0.351417, Final residual = 8.3534e+269, No Iterations 1000
DILUPBiCG: Solving for T, Initial residual = 0.00762719, Final residual = 8.13923e-06, No Iterations 23
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<d ouble>&, Foam::Field<double> const&) const at ??:?
#6 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#7 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#8 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#9 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#10 ? at ??:?
#11 ? at ??:?
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13 ? at ??:?
Floating point exception (core dumped)

After some debugging I fund out that lowering value of nu to 0.01 I am getting the output but as I lower the values to 0.001 it starts giving this output.I want to simulate it with nu=1e^-6 .What is wrong with this?I am attaching the files of the fvsolution,geometry etc.
Attached Files
File Type: zip trialexample.zip (73.9 KB, 1 views)
shaggadude is offline   Reply With Quote

Old   April 19, 2016, 10:17
Default same probleme with turbulence model
  #25
New Member
 
Mehdi
Join Date: Apr 2016
Posts: 5
Rep Power: 2
MehdiCFD is on a distinguished road
helo everyone i have the same probleme, i would like to simulate an airflow around an airfoil i used rhosimplefoam like a solver and when i'm compiling its apear this message:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 3.0.1-d8a290b55d28
Exec : rhoSimpleFoam
Date : Apr 19 2016
Time : 15:10:30
Host : "mehdi-HP-250-G4-Notebook-PC"
PID : 7408
Case : /home/mehdi/OpenFOAM/mehdi-3.0.1/run/airfoil/airfoil_rhoSimpleFoam
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


SIMPLE: convergence criteria
field p tolerance 0.01
field U tolerance 0.0001
field e tolerance 0.001
field "(k|epsilon|omega)" tolerance 0.001

Reading thermophysical properties

Selecting thermodynamics package
{
type hePsiThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
sigmak 1;
sigmaEps 1.3;
}

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/home/mehdi/OpenFOAM/OpenFOAM-3.0.1/platforms/linuxGccDPInt32Opt/lib/libOpenFOAM.so"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam:perator/(Foam::tmp<Foam::Field<double> > const&, Foam::tmp<Foam::Field<double> > const&) at ??:?
#5 Foam::fluidThermo::nu(int) const at ??:?
#6 Foam::TurbulenceModel<Foam::geometricOneField, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>, Foam::compressibleTurbulenceModel, Foam::fluidThermo>::nu(int) const at ??:?
#7 Foam::nutkWallFunctionFvPatchScalarField::calcNut( ) const at ??:?
#8 Foam::nutWallFunctionFvPatchScalarField::updateCoe ffs() at ??:?
#9 Foam::fvPatchField<double>::evaluate(Foam::UPstrea m::commsTypes) at ??:?
#10 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::evaluate() at ??:?
#11 Foam::RASModels::kEpsilon<Foam::EddyDiffusivity<Fo am::ThermalDiffusivity<Foam::CompressibleTurbulenc eModel<Foam::fluidThermo> > > >::correctNut() at ??:?
#12 Foam::RASModels::kEpsilon<Foam::EddyDiffusivity<Fo am::ThermalDiffusivity<Foam::CompressibleTurbulenc eModel<Foam::fluidThermo> > > >::kEpsilon(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&, Foam::word const&) at ??:?
#13 Foam::RASModel<Foam::EddyDiffusivity<Foam::Thermal Diffusivity<Foam::CompressibleTurbulenceModel<Foam ::fluidThermo> > > >::adddictionaryConstructorToTable<Foam::RASModels ::kEpsilon<Foam::EddyDiffusivity<Foam::ThermalDiff usivity<Foam::CompressibleTurbulenceModel<Foam::fl uidThermo> > > > >::New(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) at ??:?
#14 Foam::RASModel<Foam::EddyDiffusivity<Foam::Thermal Diffusivity<Foam::CompressibleTurbulenceModel<Foam ::fluidThermo> > > >::New(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) at ??:?
#15 ? at ??:?
#16 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#17 ? at ??:?
Exception en point flottant (core dumped)
mehdi@mehdi-HP-250-G4-Notebook-PC:~/OpenFOAM/mehdi-3.0.1/run/airfoil/airfoil_rhoSimpleFoam$
if someone can help me please it's for my master project, thank you
MehdiCFD is offline   Reply With Quote

Old   April 26, 2016, 12:03
Default
  #26
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18
immortality is on a distinguished road
Quote:
Originally Posted by shaggadude View Post
Hello Foamers,
I am new to openfoam.I have created my own solver adding temperature equation (only advection-diffusion with velocity) to the icoFoam solver. to be specific ,code is
solve
( fvm::div(phi,T)
- fvm::laplacian(DT,T)
);
I am getting this error

Courant Number mean: 0.0248669 max: 2.20364
smoothSolver: Solving for Ux, Initial residual = 0.202777, Final residual = 1.21499e+270, No Iterations 1000
smoothSolver: Solving for Uy, Initial residual = 0.351417, Final residual = 8.3534e+269, No Iterations 1000
DILUPBiCG: Solving for T, Initial residual = 0.00762719, Final residual = 8.13923e-06, No Iterations 23
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<d ouble>&, Foam::Field<double> const&) const at ??:?
#6 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#7 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#8 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#9 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#10 ? at ??:?
#11 ? at ??:?
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13 ? at ??:?
Floating point exception (core dumped)

After some debugging I fund out that lowering value of nu to 0.01 I am getting the output but as I lower the values to 0.001 it starts giving this output.I want to simulate it with nu=1e^-6 .What is wrong with this?I am attaching the files of the fvsolution,geometry etc.
Hi
your courant number seems too high.may it cause the run to be unstable. reduce it to be lower than 1 by lowering time step. may it be solved.
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Reply

Tags
(core dumped), floating point exception, icofoam

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
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 06:36
simpleFoam Floating point exception error -help sudhasran OpenFOAM Running, Solving & CFD 3 March 12, 2012 17:23
Pipe flow in settlingFoam floating point exception jochemvandenbosch OpenFOAM Running, Solving & CFD 4 February 16, 2012 04:24
block-structured mesh for t-junction Robert@cfd ANSYS Meshing & Geometry 20 November 11, 2011 05:59
Finished simulation doesn't start: floating point exception [Divide by zero] MaxCFD STAR-CCM+ 3 June 26, 2011 10:31


All times are GMT -4. The time now is 00:48.