CFD Online Discussion Forums

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

oldirtybrindley March 7, 2014 19:31
floating point exception
 
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 March 7, 2014 20:06
[
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;
}

// ************************************************************************* //

Lieven March 8, 2014 16:56
Hi olddirtybrindly,

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

Cheers,

L


All times are GMT -4. The time now is 18:31.