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

floating point exception

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 7, 2014, 20:31
Default floating point exception
  #1
New Member
 
J. Brindley
Join Date: Mar 2014
Posts: 11
Rep Power: 4
oldirtybrindley is on a distinguished road
Not sure what's going on. Can't get simpleFoam to run

Code:
Time = 1

smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.081410493, No Iterations 4
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.061596729, No Iterations 3
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.091684187, No Iterations 3
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.0059667553, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.40650934, Final residual = 0.0012374181, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.036894522, Final residual = 0.00027957842, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.0062564286, Final residual = 4.8770264e-05, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.0013606213, Final residual = 1.2695641e-05, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.00035198339, Final residual = 3.117196e-06, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.00010876533, Final residual = 1.0124916e-06, No Iterations 3
GAMG:  Solving for p, Initial residual = 3.7107364e-05, Final residual = 3.2700493e-07, No Iterations 3
GAMG:  Solving for p, Initial residual = 1.4011258e-05, Final residual = 1.2603484e-07, No Iterations 3
GAMG:  Solving for p, Initial residual = 5.5796704e-06, Final residual = 3.7850784e-08, No Iterations 6
GAMG:  Solving for p, Initial residual = 2.3340101e-06, Final residual = 6.442118e+119, No Iterations 1000
time step continuity errors : sum local = 1.2632859e+123, global = 4.4873194e+101, cumulative = 4.4873194e+101
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4  Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8  Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) at ??:?
#9  Foam::incompressible::RASModels::kOmegaSST::correct() at ??:?
#10  
 at ??:?
#11  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12  
 at ??:?
Floating point exception (core dumped)
oldirtybrindley is offline   Reply With Quote

Old   March 7, 2014, 21:06
Default
  #2
New Member
 
J. Brindley
Join Date: Mar 2014
Posts: 11
Rep Power: 4
oldirtybrindley is on a distinguished road
[
Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.3.0-f5222ca19ce6
Exec   : checkMesh
Date   : Mar 07 2014
Time   : 20:04:09
Host   : "AVA-411296"
PID    : 26818
Case   : /home/kathyb/OpenFOAM/kathyb-2.3.0/run/SplitTipTesting/AOA/simplesolvefixed(donotuse)
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

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

Create polyMesh for time = 0

Time = 0

Mesh stats
    points:           30927
    faces:            72559
    internal faces:   60772
    cells:            21783
    faces per cell:   6.1208741
    boundary patches: 5
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     15749
    prisms:        1429
    wedges:        0
    pyramids:      0
    tet wedges:    128
    tetrahedra:    2
    polyhedra:     4475
    Breakdown of polyhedra by number of faces:
        faces   number of cells
            4   1186
            5   814
            6   740
            7   218
            8   90
            9   678
           10   86
           11   47
           12   401
           13   3
           14   6
           15   150
           18   56

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
                   Patch    Faces   Points                  Surface topology
                   inlet        0        0                        ok (empty)
                  outlet        0        0                        ok (empty)
                   walls       40       67  ok (non-closed singly connected)
               frontBack        0        0                        ok (empty)
                    Foil    11747    13554  ok (non-closed singly connected)

Checking geometry...
    Overall domain bounding box (-0.020220271 -0.010812593 -0.00023150215) (0.020260984 0.010818352 0.15183358)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (-5.8212676e-17 4.5742309e-17 1.775036e-17) OK.
    Max cell openness = 2.9738872e-16 OK.
    Max aspect ratio = 21.38308 OK.
    Minimum face area = 6.6009305e-09. Maximum face area = 1.3716404e-05.  Face area magnitudes OK.
    Min volume = 1.1250083e-12. Max volume = 2.2923998e-08.  Total volume = 2.8560307e-05.  Cell volumes OK.
    Mesh non-orthogonality Max: 62.160689 average: 13.151348
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 2.6576963 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

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

solvers
{
    p
    {
        solver           GAMG;
        tolerance        0.02;
        relTol           0;
        smoother         GaussSeidel;
        nPreSweeps       0;
        nPostSweeps      2;
        cacheAgglomeration on;
        agglomerator     faceAreaPair;
        nCellsInCoarsestLevel 10;
        mergeLevels      1;
    }


    U
    {
        solver           PBiCG;
	preconditioner	 DILU;
        tolerance        1e-6;
        relTol		0.1;
    }

    k
    {
        solver           PBiCG;
        preconditioner	 DILU;
        tolerance        1e-6;
        relTol		0.1;
    }

    omega
    {
                solver           PBiCG;
preconditioner	 DILU;
        tolerance        1e-6;
        relTol		0.1;
    }

	cellDisplacement
	{
	solver	GAMG;
	tolerance	1e-5;
	relTol		0;
	smoother	GaussSeridel;
	cacheAgglomeration	true;
	nCellsInCoursestLevel	10;
	agglomerator	faceAreaPair;
	mergeLevels	1;

	}


}

SIMPLE
{
    nNonOrthogonalCorrectors 0;

    pReffCell 0;
    pRefValue 0;
    residualControl
    {
        p 1e-4;
	U 1e-4;
    }
}

potentialFlow
{
    nNonOrthogonalCorrectors 10;
}

relaxationFactors
{
    fields
    {
        p               0.3;
    }
    equations
    {
        U               0.7;
        k               0.7;
        omega           0.7;
    }
}

cache
{
    grad(U);
}

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

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         Gauss linear;
    grad(p)         Gauss linear;
    grad(U)	    Gauss linear;
    grad(k)    	    cellLimited Gauss linear 1;
    grad(omega)    	    cellLimited Gauss linear 1;
}

divSchemes
{
    default         none;
    div(phi,U)      bounded Gauss upwind;;
    div(phi,k)      bounded Gauss upwind;
    div(phi,omega)  bounded Gauss upwind;
    div((nuEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear limited corrected 0.5;
    laplacian(DomegaEff,omega)		Gauss linear limited 0.33;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         limited 0.33;
}

fluxRequired
{
    default         no;
    p;
}

// ************************************************************************* //
oldirtybrindley is offline   Reply With Quote

Old   March 8, 2014, 17:56
Default
  #3
Senior Member
 
Lieven
Join Date: Dec 2011
Location: Leuven, Belgium
Posts: 297
Rep Power: 15
Lieven will become famous soon enough
Hi olddirtybrindly,

Do you get the same crash when you run the case with pisofoam?

Cheers,

L
Lieven 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
Floating point exception with pimpleDyMFoam ebah6 OpenFOAM Running, Solving & CFD 8 May 13, 2013 09:01
Inlet Velocity Profile BC - Floating Point exception during solution initialization Janshi STAR-CCM+ 4 March 14, 2012 11:21
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


All times are GMT -4. The time now is 02:36.