Shif |
March 24, 2022 14:40 |
Floating point exception while running simpleFoam
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
|