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

Floating point exception while running simpleFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 24, 2022, 15:40
Default Floating point exception while running simpleFoam
  #1
New Member
 
Shif
Join Date: Jan 2022
Posts: 3
Rep Power: 4
Shif is on a distinguished road
Hello.

I am trying to do a simulation of a wind turbine blade (1/3 of the whole rotor) and have referred to the other thread on this forum for help every time I have encountered any problem.
Currently, I am encountering a "Floating Point Exception" when I run simpleFoam.

I have gone through some of the threads where others have faced the same issue and the first diagnosis recommended was to check their mesh. My checkMesh seems to have no problems, If it does have any issues that i have failed to identify please let me know
Code:
Create time

Create mesh for time = 0

Time = 0

Mesh stats
    points:           388192
    internal points:  328012
    faces:            4346366
    internal faces:   4226014
    cells:            2143095
    faces per cell:   4
    boundary patches: 6
    point zones:      0
    face zones:       1
    cell zones:       1

Overall number of cells of each type:
    hexahedra:     0
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    2143095
    polyhedra:     0

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
    cyclic1             7493     3860     ok (non-closed singly connected)
    cyclic2             7493     3860     ok (non-closed singly connected)
    inlet               1331     719      ok (non-closed singly connected)
    inlet_top           6034     3125     ok (non-closed singly connected)
    outlet              1319     713      ok (non-closed singly connected)
    blade               96682    48343    ok (closed singly connected)

Checking faceZone topology for multiply connected surfaces...
    FaceZone            Faces    Points   Surface topology
    interior-fluid      4226014  388190   multiply connected (shared edge)
  <<Writing 388175 conflicting points to set nonManifoldPoints

Checking basic cellZone addressing...
    CellZone            Cells        Points       Volume       BoundingBox
    cellMRFZone         2143095      388192       5.15191e+06  (-135 -116.913 -180) (0 116.913 90)

Checking geometry...
    Overall domain bounding box (-135 -116.913 -180) (0 116.913 90)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (-3.38768e-17 -3.01073e-16 -5.20473e-16) OK.
    Max cell openness = 2.94847e-16 OK.
    Max aspect ratio = 7.90672 OK.
    Minimum face area = 0.000464667. Maximum face area = 91.3371.  Face area magnitudes OK.
    Min volume = 5.23826e-06. Max volume = 281.275.  Total volume = 5.15191e+06.  Cell volumes OK.
    Mesh non-orthogonality Max: 63.8832 average: 16.2646
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 0.767696 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End
My controlDict file is
Code:
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     simpleFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         1000;

deltaT          1;

writeControl    timeStep;

writeInterval   50;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

functions
{
    forces
    {
        type forces;
        libs ("libforces.so");
        writeControl timeStep;
        writeInterval 50;
        patches (blade);
//        pname p;
//        Uname U;
        rho rhoInf;
        log true;
        rhoInf 1;
        CofR (0 0 0);
    }
}
fvSchemes
Code:
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
	default         steadyState;
}

gradSchemes
{
	default	Gauss linear;
}


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

laplacianSchemes
{
	default	Gauss linear corrected;
}

snGradSchemes
{
	default	corrected;
}

wallDist
{
	method	meshWave;
}


// ************************************************************************* //
fvSolution
Code:
solvers
{
    p
    {
        solver          GAMG;
        tolerance       1e-7;
        relTol          0;
        smoother        GaussSeidel;
		cacheAgglomeration no;
		maxIter	500;
    }

    "(U|k|omega)"
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        tolerance       1e-7;
        relTol          0;
        nSweeps         2;
    }
}
SIMPLE
{
    nNonOrthogonalCorrectors 0;
    consistent      yes;
	pRefCell	0;
	pRefValue	0;
	residualControl
	{
	p	1e-6;
	"(U|k|omega)"	1e-6;
	}
}

relaxationFactors
{
    fields
    {
        p               0.3;
    }
    equations
    {
        "(U|k|omega)"   0.7;
    }
}
my boundary conditions are as follows:
k:

Code:
FoamFile
{
    version     2.0;
    format      binary;
    arch        "LSB;label=32;scalar=64";
    class       volScalarField;
    location    "0";
    object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0.06;

boundaryField
{
    blade
    {
        type            kqRWallFunction;
        value           $internalField;
    }
    inlet
    {
        type            fixedValue;
        value           $internalField;
    }
    inlet_top
    {
        type            fixedValue;
        value           $internalField;
    }
    outlet
    {
        type            zeroGradient;
    }
    cyclic1
    {
        type            cyclicAMI;
    }
	  cyclic2
    {
        type            cyclicAMI;
    }
}


// ************************************************************************* //
nut:

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

internalField   uniform 0;

boundaryField
{
    blade
    {
        type            nutkWallFunction;
        value           $internalField;
    }
    inlet
    {
        type            calculated;
        value           $internalField;
    }
    inlet_top
    {
        type            calculated;
        value           $internalField;
    }
    outlet
    {
        type            calculated;
        value           $internalField;
    }
    cyclic1
    {
        type            cyclicAMI;
    }
	  cyclic2
    {
        type            cyclicAMI;
    }
}


// ************************************************************************* //
omega:

Code:
FoamFile
{
    version     2.0;
    format      binary;
    arch        "LSB;label=32;scalar=64";
    class       volScalarField;
    location    "0";
    object      omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 4.5939;

boundaryField
{
    blade
    {
        type            omegaWallFunction;
        value           $internalField;
    }
    inlet
    {
        type            fixedValue;
        value           $internalField;
    }
    inlet_top
    {
        type            fixedValue;
        value           $internalField;
    }
    outlet
    {
        type            zeroGradient;
    }
    cyclic1
    {
        type            cyclicAMI;
    }
  	cyclic2
    {
        type            cyclicAMI;
    }
}


// ************************************************************************* //
p:

Code:
FoamFile
{
    version     2.0;
    format      binary;
    arch        "LSB;label=32;scalar=64";
    class       volScalarField;
    location    "0";
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{
    blade
    {
        type            zeroGradient;
    }
    inlet
    {
        type            zeroGradient;
    }
    inlet_top
    {
        type            zeroGradient;
    }
    outlet
    {
        type            fixedValue;
        value           $internalField;
    }
    cyclic1
    {
        type            cyclicAMI;
    }
	cyclic2
    {
        type            cyclicAMI;
    }
}


// ************************************************************************* //
U:

Code:
FoamFile
{
    version     2.0;
    format      binary;
    arch        "LSB;label=32;scalar=64";
    class       volVectorField;
    location    "0";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform (0 0 -10);

boundaryField
{
    blade
    {
        type            noSlip;
    }
    inlet
    {
        type            fixedValue;
        value           $internalField;
    }
    inlet_top
    {
        type            fixedValue;
        value           $internalField;
    }
    outlet
    {
        type            zeroGradient;
    }
    cyclic1
    {
        type            cyclicAMI;
    }
    cyclic2
    {
        type            cyclicAMI;
    }
}

// ************************************************************************* //
Also one thing i noticed was when i run simpleFoam -dry-run i get this warning, I tried looking up for a solution but i couldnt find much resources
Quote:
--> FOAM Warning :
From void Foam::faceZone::checkAddressing() const
in file meshes/polyMesh/zones/faceZone/faceZone.C at line 187
Illegal face index 21 outside range 0..20
could this warning be the root of the problem? I could not find any solution to fix this.

Thank you
Shif is offline   Reply With Quote

Reply

Tags
floating point exception, mrfsimplefoam, wind blade

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Floating Point exception (core dumped) chtMultiRegionFoam ayu12 OpenFOAM Running, Solving & CFD 1 October 12, 2021 12:31
Floating point exception (core dumped), running a new solver Mahyar Javidi OpenFOAM Running, Solving & CFD 6 April 7, 2018 13:43
Floating point exception (yes...) sandri_92 OpenFOAM Running, Solving & CFD 0 October 23, 2016 15:42
Floating point exception (core dumped) for GAMG solver yuhou1989 OpenFOAM Running, Solving & CFD 2 March 24, 2015 20:28
simpleFoam Floating point exception error -help sudhasran OpenFOAM Running, Solving & CFD 3 March 12, 2012 17:23


All times are GMT -4. The time now is 09:27.