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

Floating point exception error

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 15, 2009, 21:05
Default Floating point exception error
  #1
New Member
 
Alan Harrland
Join Date: Mar 2009
Posts: 21
Rep Power: 8
Alan is on a distinguished road
Hi,

I am trying to run a scenario of a three dimensional wedge in supersonic flow, using rhoCentralFoam. The solution will iterate for a while, but then it will crash and I get the following error:

#0 Foam::error:rintStack(Foam::Ostream&) in "/home/cdoolan/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/cdoolan/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) in "/home/cdoolan/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/cdoolan/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/rhoCentralFoam"
#5 main in "/home/cdoolan/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/rhoCentralFoam"
#6 __libc_start_main in "/lib/libc.so.6"
#7 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/i386/elf/start.S:122

Floating point exception

I have tried turning off the floating point exception catching, which allows the solver to iterate until the end time, but I cannot view in post processing as the solution space contains NaN values. I am building my mesh in Gmsh, and importing into OpenFOAM (with which I have no problems).

My initial condition files are below:

p:
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.6                                   |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 1;

boundaryField
{
    inlet           
    {
        type            fixedValue;
        value           uniform 1;
    }

    outlet          
    {
        type            zeroGradient;
    }

    obstacle
    {
        type            zeroGradient;
    }

    defaultFaces
    {
    type        empty;
    }
}

// ************************************************************************* //
T:
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.6                                   |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 1;

boundaryField
{
    inlet           
    {
        type            fixedValue;
        value           uniform 1;
    }

    outlet          
    {
        type            zeroGradient;
    }

    obstacle
    {
        type            zeroGradient;
    }

    defaultFaces    
    {
        type            empty;
    }
}

// ************************************************************************* //
U:
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.6                                   |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    inlet           
    {
        type            fixedValue;
        value           uniform (5 0 0);
    }

    outlet          
    {
        type            zeroGradient;
    }

    obstacle
    {
        type            slip;
    }

    defaultFaces    
    {
        type            empty;
    }
}

// ************************************************************************* //
and my boundary file:
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.6                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

4
(
    inlet
    {
        type            patch;
        nFaces          156;
        startFace       14250;
    }
    outlet
    {
        type            patch;
        nFaces          894;
        startFace       14406;
    }
    obstacle
    {
        type            patch;
        nFaces          114;
        startFace       15300;
    }
    defaultFaces
    {
        type            empty;
        nFaces          0;
        startFace       15414;
    }
)

// ************************************************************************* //
I can email the working directory (its 600kb) if that would help, but it is too large to upload here (97kb limit).

Thanks, any help would be much appreciated, as I am pulling my hair out here, trying to figure this out.

Thanks,
Alan
Alan is offline   Reply With Quote

Old   December 26, 2009, 01:38
Default
  #2
New Member
 
Alan Harrland
Join Date: Mar 2009
Posts: 21
Rep Power: 8
Alan is on a distinguished road
Is there anyone here willing to have a look at this for me? I can't for the life of me figure it out. It would be greatly appreciated as I am pulling my hair out.

Thanks in advance,
Alan
Alan is offline   Reply With Quote

Old   December 27, 2009, 08:10
Wink
  #3
New Member
 
Kerstin Heinen
Join Date: Mar 2009
Location: Ludwigshafen, Germany
Posts: 27
Rep Power: 8
kerstin is on a distinguished road
Hi Alan,
I would say even as advanced cfd user it is nearly impossible to give the answer to your problem...

Just some advice how I would proceed if I were you:

-Start writing out results very often and look at intermediate results. If crashing is due to bad mesh quality you may see it in the intermediate results , that some single cells get large values which are then later "exploding"

- check your direction vector for velocity in the boundary condition, does it really match to the geometry?

- how do you initialise the solution? If this is a turbulent flow, don't start from scratch...You can run a "potentialFoam" before or reduce your velocity magnitude, to have less severe starting conditions for the solver...

- make a "checkMesh"...with errors in the output there my experience is, that you will not have a stable run in most cases...

Good luck for this "hunt". It wouldn't be CFD, if it would never happen that one has exactly the same problems, you have now.
Kerstin
kerstin is offline   Reply With Quote

Old   December 27, 2009, 11:32
Default
  #4
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 208
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Maybe you could try to modify the pressure initial conditions to

Code:
    inlet           
    {

        type            zeroGradient;

    }

    outlet          
    {
        type            fixedValue;
        value           uniform 0;
    }
louisgag is offline   Reply With Quote

Old   December 27, 2009, 11:49
Default
  #5
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 208
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
I also have a FPE error that keeps coming back and I can't figure out what initiates it...

Code:
Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider
Evaluation of GGI weighting factors:
  Largest slave weighting factor correction : 0.000117224699752 average: 0.000117210464361
  Largest master weighting factor correction: 4.4408920985e-16 average: 1.09561482693e-16

smoothSolver:  Solving for Ux, Initial residual = 4.18366053939e-05, Final residual = 3.99534809799e-11, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.000183494576854, Final residual = 1.58095734727e-10, No Iterations 5
PCG:  Solving for p, Initial residual = 0.00115761064277, Final residual = 9.69845796082e-07, No Iterations 285
PCG:  Solving for p, Initial residual = 9.99140357817e-05, Final residual = 9.84892398469e-07, No Iterations 16
time step continuity errors : sum local = 1.61074737955e-11, global = -1.38664307317e-13, cumulative = -2.44455566515e-08
PCG:  Solving for p, Initial residual = 2.7137207064e-05, Final residual = 9.65820558314e-07, No Iterations 41
PCG:  Solving for p, Initial residual = 1.15863553234e-05, Final residual = 9.31153715832e-07, No Iterations 3
time step continuity errors : sum local = 1.5228622047e-11, global = -5.55997394556e-14, cumulative = -2.44456122512e-08
PBiCG:  Solving for omega, Initial residual = 1.98154832644e-05, Final residual = 3.46778670816e-08, No Iterations 1
PBiCG:  Solving for k, Initial residual = 8.33751071488e-05, Final residual = 2.53152769991e-07, No Iterations 1
ExecutionTime = 16216.54 s  ClockTime = 16307 s

Courant Number mean: 0.0394310396556 max: 0.499997648141 velocity magnitude: 124.35829345
deltaT = 9.57280153069e-06
Time = 0.316468

Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider
Evaluation of GGI weighting factors:
  Largest slave weighting factor correction : 2.36652718848e-05 average: 2.36504740638e-05
  Largest master weighting factor correction: 4.4408920985e-16 average: 1.09561482693e-16

smoothSolver:  Solving for Ux, Initial residual = 4.18646624848e-05, Final residual = 4.00745405967e-11, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.000183554079661, Final residual = 1.58722402612e-10, No Iterations 5
PCG:  Solving for p, Initial residual = 0.00113267387872, Final residual = 9.99882920123e-07, No Iterations 149
PCG:  Solving for p, Initial residual = 9.82117938326e-05, Final residual = 9.95912006942e-07, No Iterations 15
time step continuity errors : sum local = 1.62857195533e-11, global = 1.21882881021e-12, cumulative = -2.44443934224e-08
PCG:  Solving for p, Initial residual = 2.63425356134e-05, Final residual = 9.5916336339e-07, No Iterations 39
PCG:  Solving for p, Initial residual = 1.13546165489e-05, Final residual = 9.1377403705e-07, No Iterations 3
time step continuity errors : sum local = 1.4942541397e-11, global = 1.21024107242e-12, cumulative = -2.44431831814e-08
PBiCG:  Solving for omega, Initial residual = 1.97786652847e-05, Final residual = 3.4608657236e-08, No Iterations 1
PBiCG:  Solving for k, Initial residual = 8.33789411068e-05, Final residual = 2.53398294466e-07, No Iterations 1
bounding k, min: -0.122378016374 max: 237.949650915 average: 11.0237179132
ExecutionTime = 16216.96 s  ClockTime = 16308 s

Courant Number mean: 0.0394313290225 max: 0.499998490876 velocity magnitude: 124.35851925
deltaT = 9.57280153069e-06
Time = 0.316477

Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider
Evaluation of GGI weighting factors:
  Largest slave weighting factor correction : 0.000160794100478 average: 0.000160779324961
  Largest master weighting factor correction: 4.4408920985e-16 average: 1.22708860616e-16

smoothSolver:  Solving for Ux, Initial residual = 4.18648802413e-05, Final residual = 4.01555017906e-11, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.000183664619068, Final residual = 1.5925205909e-10, No Iterations 5
PCG:  Solving for p, Initial residual = 0.0012822507928, Final residual = 9.74487717197e-07, No Iterations 314
PCG:  Solving for p, Initial residual = 0.000183503994712, Final residual = 9.54017668329e-07, No Iterations 27
time step continuity errors : sum local = 1.56040920447e-11, global = 2.71683168776e-14, cumulative = -2.4443156013e-08
PCG:  Solving for p, Initial residual = 4.54649995932e-05, Final residual = 9.7803002942e-07, No Iterations 92
PCG:  Solving for p, Initial residual = 1.87289344062e-05, Final residual = 8.32829759138e-07, No Iterations 5
time step continuity errors : sum local = 1.36217381462e-11, global = -7.85213463584e-13, cumulative = -2.44439412265e-08
PBiCG:  Solving for omega, Initial residual = 1.97895544736e-05, Final residual = 3.43292848071e-08, No Iterations 1
PBiCG:  Solving for k, Initial residual = 8.33872916715e-05, Final residual = 2.5269138733e-07, No Iterations 1
ExecutionTime = 16217.56 s  ClockTime = 16308 s

Courant Number mean: 0.0394316251525 max: 0.499999326467 velocity magnitude: 124.358741937
deltaT = 9.57280153069e-06
Time = 0.316487

Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider
Evaluation of GGI weighting factors:
  Largest slave weighting factor correction : 0.000285528627237 average: 0.000285513904142
  Largest master weighting factor correction: 4.4408920985e-16 average: 9.49532850008e-17

smoothSolver:  Solving for Ux, Initial residual = 4.18150647223e-05, Final residual = 4.02129097578e-11, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.000183919043323, Final residual = 1.59663449673e-10, No Iterations 5
PCG:  Solving for p, Initial residual = 0.00118091660555, Final residual = 9.97058782232e-07, No Iterations 280
PCG:  Solving for p, Initial residual = 0.000103443967517, Final residual = 8.8633487677e-07, No Iterations 12
time step continuity errors : sum local = 1.4493480812e-11, global = 2.1001888378e-13, cumulative = -2.44437312076e-08
PCG:  Solving for p, Initial residual = 2.80136398208e-05, Final residual = 9.62893451566e-07, No Iterations 42
PCG:  Solving for p, Initial residual = 1.11499919963e-05, Final residual = 9.28004174336e-07, No Iterations 3
time step continuity errors : sum local = 1.51747798538e-11, global = 4.16825424229e-13, cumulative = -2.44433143822e-08
PBiCG:  Solving for omega, Initial residual = 1.98065657161e-05, Final residual = 3.40451890596e-08, No Iterations 1
PBiCG:  Solving for k, Initial residual = 8.33920881206e-05, Final residual = 2.52098556426e-07, No Iterations 1
ExecutionTime = 16218.13 s  ClockTime = 16309 s

Courant Number mean: 0.0394319159774 max: 0.500000147171 velocity magnitude: 124.358960763
deltaT = 9.57170436146e-06
Time = 0.316496

Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider
Evaluation of GGI weighting factors:
  Largest slave weighting factor correction : 0.000397873642995 average: 0.000397859002476
  Largest master weighting factor correction: 4.4408920985e-16 average: 1.18326401309e-16

smoothSolver:  Solving for Ux, Initial residual = 4.17820269512e-05, Final residual = 4.02754195966e-11, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.000183916215376, Final residual = 1.60180223308e-10, No Iterations 5
PCG:  Solving for p, Initial residual = 0.00116209995737, Final residual = 9.91222118143e-07, No Iterations 96
PCG:  Solving for p, Initial residual = 9.7432771035e-05, Final residual = 8.65531722367e-07, No Iterations 12
time step continuity errors : sum local = 1.41496821471e-11, global = -1.1742061533e-12, cumulative = -2.44444885883e-08
PCG:  Solving for p, Initial residual = 2.64247799783e-05, Final residual = 9.26829060343e-07, No Iterations 41
PCG:  Solving for p, Initial residual = 1.09427368381e-05, Final residual = 9.06631309017e-07, No Iterations 3
time step continuity errors : sum local = 1.48215256106e-11, global = -1.20279142129e-12, cumulative = -2.44456913798e-08
PBiCG:  Solving for omega, Initial residual = 1.98363312799e-05, Final residual = 3.4374674519e-08, No Iterations 1
PBiCG:  Solving for k, Initial residual = 8.33873645469e-05, Final residual = 2.51784635071e-07, No Iterations 1
ExecutionTime = 16218.51 s  ClockTime = 16309 s

Courant Number mean: 0.0394276758498 max: 0.499943657681 velocity magnitude: 124.359179205
deltaT = 9.57170436146e-06
Time = 0.316506

Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider
Evaluation of GGI weighting factors:
  Largest slave weighting factor correction : 0.000497859420717 average: 0.000497844892901
  Largest master weighting factor correction: 4.4408920985e-16 average: 8.47275466161e-17

smoothSolver:  Solving for Ux, Initial residual = 4.17905810822e-05, Final residual = 4.03898968375e-11, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.000184025898617, Final residual = 1.60739458371e-10, No Iterations 5
PCG:  Solving for p, Initial residual = 0.00112936194147, Final residual = 9.74360513696e-07, No Iterations 274
PCG:  Solving for p, Initial residual = 9.76900241161e-05, Final residual = 9.7094946801e-07, No Iterations 11
time step continuity errors : sum local = 1.58701620076e-11, global = 1.80901108399e-13, cumulative = -2.44455104787e-08
PCG:  Solving for p, Initial residual = 2.57632315604e-05, Final residual = 9.37704584393e-07, No Iterations 41
PCG:  Solving for p, Initial residual = 1.09362657999e-05, Final residual = 8.88228716267e-07, No Iterations 3
time step continuity errors : sum local = 1.45180616517e-11, global = 3.24254714623e-13, cumulative = -2.44451862239e-08
PBiCG:  Solving for omega, Initial residual = 1.9899688921e-05, Final residual = 3.67212803247e-08, No Iterations 1
PBiCG:  Solving for k, Initial residual = 8.33898572745e-05, Final residual = 2.51738401993e-07, No Iterations 1
bounding k, min: -0.802415643715 max: 237.717244127 average: 11.0287760469
ExecutionTime = 16219.02 s  ClockTime = 16310 s

Courant Number mean: 0.0394279428752 max: 0.499944485086 velocity magnitude: 124.359399357
deltaT = 9.57170436146e-06
Time = 0.316516

Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider
Evaluation of GGI weighting factors:
  Largest slave weighting factor correction : 0.000585499543768 average: 0.000585485158807
  Largest master weighting factor correction: 4.4408920985e-16 average: 9.64141047701e-17

smoothSolver:  Solving for Ux, Initial residual = 4.17858218617e-05, Final residual = 4.05074168974e-11, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.000184094269779, Final residual = 1.61390224942e-10, No Iterations 5
PCG:  Solving for p, Initial residual = 0.00112989517074, Final residual = 9.91499532231e-07, No Iterations 152
PCG:  Solving for p, Initial residual = 9.69060456898e-05, Final residual = 9.08043230719e-07, No Iterations 12
time step continuity errors : sum local = 1.48405476369e-11, global = -9.5042254985e-13, cumulative = -2.44461366465e-08
PCG:  Solving for p, Initial residual = 2.59766736557e-05, Final residual = 9.42002427652e-07, No Iterations 41
PCG:  Solving for p, Initial residual = 1.1185836139e-05, Final residual = 9.04423006686e-07, No Iterations 3
time step continuity errors : sum local = 1.47813793588e-11, global = -1.04539352081e-12, cumulative = -2.444718204e-08
PBiCG:  Solving for omega, Initial residual = 1.98655390897e-05, Final residual = 3.54224756962e-08, No Iterations 1
PBiCG:  Solving for k, Initial residual = 8.3438685774e-05, Final residual = 2.51623930413e-07, No Iterations 1
ExecutionTime = 16219.45 s  ClockTime = 16310 s

Courant Number mean: 0.0394282050904 max: 0.499945308896 velocity magnitude: 124.35961817
deltaT = 9.57170436146e-06
Time = 0.316525

Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider
Evaluation of GGI weighting factors:
  Largest slave weighting factor correction : 0.000660805915266 average: 0.000660791703336
  Largest master weighting factor correction: 3.33066907388e-16 average: 6.42760698467e-17

smoothSolver:  Solving for Ux, Initial residual = 4.17822084748e-05, Final residual = 4.06156755246e-11, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.000184161482139, Final residual = 1.619696451e-10, No Iterations 5
PCG:  Solving for p, Initial residual = 0.00111502624192, Final residual = 9.94145230431e-07, No Iterations 280
PCG:  Solving for p, Initial residual = 9.81284451381e-05, Final residual = 9.83206605956e-07, No Iterations 11
time step continuity errors : sum local = 1.6066266089e-11, global = -1.6661252468e-13, cumulative = -2.44473486525e-08
PCG:  Solving for p, Initial residual = 2.56735821773e-05, Final residual = 9.78307634024e-07, No Iterations 43
PCG:  Solving for p, Initial residual = 1.10964929717e-05, Final residual = 8.87408739619e-07, No Iterations 3
time step continuity errors : sum local = 1.45008481624e-11, global = 1.02944097058e-14, cumulative = -2.44473383581e-08
PBiCG:  Solving for omega, Initial residual = 1.98451187283e-05, Final residual = 2.70582444337e-08, No Iterations 1
PBiCG:  Solving for k, Initial residual = 8.34347962654e-05, Final residual = 2.51644588879e-07, No Iterations 1
ExecutionTime = 16219.97 s  ClockTime = 16311 s

Courant Number mean: 0.0394284702228 max: 0.499946124724 velocity magnitude: 124.359835338
deltaT = 9.57170436146e-06
Time = 0.316535

Initializing the GGI interpolator between master/shadow patches: insideSlider/outsideSlider
Evaluation of GGI weighting factors:
  Largest slave weighting factor correction : 0.000723788760952 average: 0.000723774752211
  Largest master weighting factor correction: 3.33066907388e-16 average: 9.20316454623e-17

smoothSolver:  Solving for Ux, Initial residual = 4.17873817625e-05, Final residual = 4.07204477989e-11, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.000184210804886, Final residual = 1.62641723107e-10, No Iterations 5
PCG:  Solving for p, Initial residual = 0.00110876280886, Final residual = 9.61326185025e-07, No Iterations 204
PCG:  Solving for p, Initial residual = 9.64025215336e-05, Final residual = 9.69402485222e-07, No Iterations 11
time step continuity errors : sum local = 1.58389919616e-11, global = -1.25756269618e-12, cumulative = -2.44485959208e-08
PCG:  Solving for p, Initial residual = 2.5342476828e-05, Final residual = 9.85823081162e-07, No Iterations 41
PCG:  Solving for p, Initial residual = 1.10396708963e-05, Final residual = 9.27978145016e-07, No Iterations 3
time step continuity errors : sum local = 1.51621484794e-11, global = -1.19091666363e-12, cumulative = -2.44497868375e-08
[louis-dell:03845] *** Process received signal ***
[louis-dell:03845] Signal: Floating point exception (8)
[louis-dell:03845] Signal code:  (-6)
[louis-dell:03845] Failing at address: 0x3e800000f05
[louis-dell:03846] *** Process received signal ***
[louis-dell:03846] Signal: Floating point exception (8)
[louis-dell:03846] Signal code:  (-6)
[louis-dell:03846] Failing at address: 0x3e800000f06
[louis-dell:03845] [ 0] /lib/libc.so.6 [0x7f46d1640530]
[louis-dell:03845] [ 1] /lib/libc.so.6(gsignal+0x35) [0x7f46d16404b5]
[louis-dell:03845] [ 2] /lib/libc.so.6 [0x7f46d1640530]
[louis-dell:03845] [ 3] /home/louis/OpenFOAM/OpenFOAM-1.5-dev/lib/linux64GccDPOpt/libOpenFOAM.so(_ZNK4Foam5PBiCG5solveERNS_5FieldIdEERKS2_h+0xbd8) [0x7f46d255f418]
[louis-dell:03845] [ 4] /home/louis/OpenFOAM/OpenFOAM-1.5-dev/lib/linux64GccDPOpt/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE5solveERNS_7IstreamE+0x13f) [0x7f46d327b71f]
[louis-dell:03845] [ 5] /home/louis/OpenFOAM/OpenFOAM-1.5-dev/lib/linux64GccDPOpt/libincompressibleRASModels.so(_ZN4Foam5solveIdEENS_9lduMatrix17solverPerformanceERKNS_3tmpINS_8fvMatrixIT_EEEE+0x4e) [0x7f46d3e7fc0e]
[louis-dell:03845] [ 6] /home/louis/OpenFOAM/OpenFOAM-1.5-dev/lib/linux64GccDPOpt/libincompressibleRASModels.so(_ZN4Foam14incompressible9RASModels9kOmegaSST7correctEv+0x123f) [0x7f46d3ebe79f]
[louis-dell:03845] [ 7] turbDyMFoam [0x41a488]
[louis-dell:03845] [ 8] /lib/libc.so.6(__libc_start_main+0xfd) [0x7f46d162babd]
[louis-dell:03845] [ 9] turbDyMFoam [0x415169]
[louis-dell:03845] *** End of error message ***
[louis-dell:03846] [ 0] /lib/libc.so.6 [0x7ff0e7f99530]
[louis-dell:03846] [ 1] /lib/libc.so.6(gsignal+0x35) [0x7ff0e7f994b5]
[louis-dell:03846] [ 2] /lib/libc.so.6 [0x7ff0e7f99530]
[louis-dell:03846] [ 3] /home/louis/OpenFOAM/OpenFOAM-1.5-dev/lib/linux64GccDPOpt/libOpenFOAM.so(_ZNK4Foam5PBiCG5solveERNS_5FieldIdEERKS2_h+0xbd8) [0x7ff0e8eb8418]
[louis-dell:03846] [ 4] /home/louis/OpenFOAM/OpenFOAM-1.5-dev/lib/linux64GccDPOpt/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE5solveERNS_7IstreamE+0x13f) [0x7ff0e9bd471f]
[louis-dell:03846] [ 5] /home/louis/OpenFOAM/OpenFOAM-1.5-dev/lib/linux64GccDPOpt/libincompressibleRASModels.so(_ZN4Foam5solveIdEENS_9lduMatrix17solverPerformanceERKNS_3tmpINS_8fvMatrixIT_EEEE+0x4e) [0x7ff0ea7d8c0e]
[louis-dell:03846] [ 6] /home/louis/OpenFOAM/OpenFOAM-1.5-dev/lib/linux64GccDPOpt/libincompressibleRASModels.so(_ZN4Foam14incompressible9RASModels9kOmegaSST7correctEv+0x123f) [0x7ff0ea81779f]
[louis-dell:03846] [ 7] turbDyMFoam [0x41a488]
[louis-dell:03846] [ 8] /lib/libc.so.6(__libc_start_main+0xfd) [0x7ff0e7f84abd]
[louis-dell:03846] [ 9] turbDyMFoam [0x415169]
[louis-dell:03846] *** End of error message ***
mpirun noticed that job rank 0 with PID 3845 on node louis-dell exited on signal 8 (Floating point exception). 
1 additional process aborted (not shown)
From this output all I can tell is that it occurs all of a sudden right before calculating k and omega.. Usually around the same timestep. However, if I start from the last saved state, the FPE will usually take longer to appear.. Anyone has any hints for me?

Thanks,

-Louis

PS: I will provide the full logfile if anyone is interested.

Last edited by louisgag; December 27, 2009 at 22:14.
louisgag is offline   Reply With Quote

Old   December 28, 2009, 03:54
Default
  #6
Member
 
matej forman
Join Date: Mar 2009
Location: Brno, Czech Republic
Posts: 92
Rep Power: 8
matejfor is on a distinguished road
Alan,

I'd support the idea with mesh checking!
Although I've never done any supersonic flow computations, I'm not very sure with your boundary cond. First if the obstacle is wall, it should have not a patch, but wall keyword in the boundary file, but I'm not sure if it makes difference. What is more important I think you need some more advanced BC tuned for supersonic pressure then zeroGradient. Looking at the tutorials for your solver I see there are some pressureTransmittive BCs.
For example I do not understand your pressure at inlet = 1;
Why 1?

My guess would be identify the closest tutorial to your case and have a look at the tutorial settings for BC.

good luck

matej
matejfor is offline   Reply With Quote

Old   December 28, 2009, 19:01
Default
  #7
New Member
 
Alan Harrland
Join Date: Mar 2009
Posts: 21
Rep Power: 8
Alan is on a distinguished road
Thanks for the replies guys. I think the problem may lie in my mesh. There are some highly skew elements along the wedge's surface, and when the shockwave reaches these elements, the simulation crashes.

I don't think there should be a problem with the BC's, as it is based on other supersonic tutorials in the rhoCentralFoam solver, and I have done other simulations without problems. I will try and play around with them though to see if it helps.

Thanks again for the responses.

Alan
Alan is offline   Reply With Quote

Old   January 13, 2010, 20:30
Default
  #8
Member
 
Mihir
Join Date: Mar 2009
Posts: 40
Rep Power: 8
mihir1310 is on a distinguished road
Im trying to solve a free-jet problem by rhoCentralFoam [OF 1.5]. My domain represents a close chamber in which the jet enters at Supersonic conditions.
My code crashes due to a FPE , I am getting -ve pressure & density at the outlet .

Quote:
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/mihirsamel/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/mihirsamel/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib/libc.so.6"
#3 Foam::ePsiThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::calculate() in "/home/mihirsamel/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#4 Foam::ePsiThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::correct() in "/home/mihirsamel/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#5 main in "/home/mihirsamel/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhoCentralFoam"
#6 __libc_start_main in "/lib/libc.so.6"
#7 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
Floating point exception
this is obviously due to solving the h equation. I do not understand what is causing this error . ill post my BC;s & the thermoPhysicalProperties directory .

Code:
dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 42000;

boundaryField
{
    wall
    {
        type            slip;
    }
    axis
    {
        type            empty;
    }
    outlet
    {
        type            nonReflective;
        refValue        uniform 6600;
        gamma           1.4;
        phi             phi;
        psi             psi;
        valueFraction   uniform 0.2;
        refGradient     uniform 0.2;
        value           uniform 6600;
    }
    inlet
    {
        type            totalPressure;
        U               U;
        phi             phi;
        rho             none;
        psi             none;
        gamma           1.4;
        p0              uniform 101325;
        value           uniform 101325;
    }
    frontAndBackPlanes
    {
        type            empty;
    }
    front
    {
        type            wedge;
    }
    back
    {
        type            wedge;
    }
}
U
Code:
dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (260.39 0 0);

boundaryField
{
    wall
    {
        type            slip;
    }
    axis
    {
        type            empty;
    }
    outlet
    {
        type            zeroGradient;
    }
    inlet
    {
        type            fixedValue;
        value           uniform (315 0 0);
    }
    frontAndBackPlanes
    {
        type            empty;
    }
    front
    {
        type            wedge;
    }
    back
    {
        type            wedge;
    }
}
Code:
thermoType      hThermo<pureMixture<sutherlandTransport<specieThermo<hConstThermo<perfectGas>>>>>;

mixture         perfectGas 1 28.9 1004.3 200.0 1.458e-6 110.4;

rhoMin          rhoMin     [1 -3 0 0 0]     0.001;

Last edited by mihir1310; January 13, 2010 at 21:02.
mihir1310 is offline   Reply With Quote

Old   February 8, 2012, 04:31
Default Floating error
  #9
Member
 
Join Date: Jan 2012
Posts: 58
Rep Power: 5
sheikh nasir is on a distinguished road
Hello
i am getting floating error:invalid number. I am working on train travelling in tunnel. Can any body help me. My email is sheikhnasir39@gmail.com

thanks
sheikh nasir is offline   Reply With Quote

Old   April 6, 2012, 13:33
Default
  #10
New Member
 
Join Date: Jan 2010
Posts: 23
Rep Power: 7
jdiorio is on a distinguished road
Quote:
Originally Posted by louisgag View Post
I also have a FPE error that keeps coming back and I can't figure out what initiates it...

From this output all I can tell is that it occurs all of a sudden right before calculating k and omega.. Usually around the same timestep. However, if I start from the last saved state, the FPE will usually take longer to appear.. Anyone has any hints for me?

Thanks,

-Louis

PS: I will provide the full logfile if anyone is interested.
Hi Louis,

Do you have time-step adjustment on?

I noticed a similar behavior when using a moving mesh with GGI doing cell layer/addition. Residuals were fine and then suddenly an FPE. I could restart from the latest time step and it would make it past that point, but then at some later time it would happen.

I eventually set an maximum time-step that was pretty low, and the problem went away - no more need to baby sit and keep restarting. Don't know if that helps.
jdiorio is offline   Reply With Quote

Old   April 6, 2012, 14:02
Default
  #11
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 208
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Hi James,

it has been quite a while since I found a way out of this problem, so I don't recall exactly how I did it. Thanks for the tip, I will keep it in mind if I stumble against such a problem again.


-Louis
louisgag 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
Installation of Netgen in SuSE Linux 92 edvardsenpriv Open Source Meshers: Gmsh, Netgen, CGNS, ... 23 January 16, 2009 07:12
Problem with compile the setParabolicInlet ivanyao OpenFOAM Running, Solving & CFD 6 September 5, 2008 20:50
How to get the max value of the whole field waynezw0618 OpenFOAM Running, Solving & CFD 4 June 17, 2008 05:07
Gmsh and samplesurface touf Open Source Meshers: Gmsh, Netgen, CGNS, ... 2 December 10, 2007 03:27
Floating point error Nabil FLUENT 2 January 19, 2007 06:30


All times are GMT -4. The time now is 23:25.