[Shif]

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

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

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

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

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:

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:

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:

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:

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:

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

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