CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Floating point exception error (http://www.cfd-online.com/Forums/openfoam-solving/71116-floating-point-exception-error.html)

Alan December 15, 2009 21:05

Floating point exception error
 
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::printStack(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 December 26, 2009 01:38

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

kerstin December 27, 2009 08:10

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

louisgag December 27, 2009 11:32

Maybe you could try to modify the pressure initial conditions to

Code:

    inlet         
    {

        type            zeroGradient;

    }

    outlet         
    {
        type            fixedValue;
        value          uniform 0;
    }


louisgag December 27, 2009 11:49

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.

matejfor December 28, 2009 03:54

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

Alan December 28, 2009 19:01

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

mihir1310 January 13, 2010 20:30

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::printStack(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::pureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam::perfectGas> > > > >::calculate() in "/home/mihirsamel/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#4 Foam::ePsiThermo<Foam::pureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam::perfectGas> > > > >::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;


sheikh nasir February 8, 2012 04:31

Floating error
 
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

jdiorio April 6, 2012 13:33

Quote:

Originally Posted by louisgag (Post 240991)
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.

louisgag April 6, 2012 14:02

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


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