|
[Sponsors] |
Floating point exception while running simpleFoam |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 | |
New Member
Shif
Join Date: Jan 2022
Posts: 3
Rep Power: 4 ![]() |
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 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); } } 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; } // ************************************************************************* // 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; } } 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; } } // ************************************************************************* // 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; } } // ************************************************************************* // 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; } } // ************************************************************************* // 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; } } // ************************************************************************* // 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; } } // ************************************************************************* // Quote:
Thank you |
||
![]() |
![]() |
![]() |
Tags |
floating point exception, mrfsimplefoam, wind blade |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Floating Point exception (core dumped) chtMultiRegionFoam | ayu12 | OpenFOAM Running, Solving & CFD | 1 | October 12, 2021 11:31 |
Floating point exception (core dumped), running a new solver | Mahyar Javidi | OpenFOAM Running, Solving & CFD | 6 | April 7, 2018 12:43 |
Floating point exception (yes...) | sandri_92 | OpenFOAM Running, Solving & CFD | 0 | October 23, 2016 14:42 |
Floating point exception (core dumped) for GAMG solver | yuhou1989 | OpenFOAM Running, Solving & CFD | 2 | March 24, 2015 19:28 |
simpleFoam Floating point exception error -help | sudhasran | OpenFOAM Running, Solving & CFD | 3 | March 12, 2012 16:23 |