|
[Sponsors] |
November 1, 2016, 23:53 |
[simpleFoam] Floating point exception
|
#1 |
New Member
Siyu Chen
Join Date: Sep 2016
Location: Houghton, USA
Posts: 7
Rep Power: 9 |
Hi, I am Siyu, a beginner.
I want to simulate the flow through aggregates in 2D using simpleFoam solver. The aggregates are scanned from X-ray CT. Then I import just two pictures into 3D Slicer, and export as STL file. I used blockMesh, snappyHexMesh, and checkMesh. But when I type "simpleFoam", this error coming out. According to some of the information through Google and on this forum, I thought it may be (1) Boundary condition; (2) Meshing. Actually, I have revised in these two aspects for days. Problem still can NOT be resolved. Please help me think about it. Thanks. All of the files before meshing are available at (I uploaded on Google Drive): https://drive.google.com/open?id=0B4...GcwQTdHdmptWk0 Some details as following: ----------------'Geometry'----------------- Set up face as the inlet, bottom face as the outlet, no value for left and right face, empty for front and back faces. https://drive.google.com/open?id=0B4...01PRFBRcHlWWmc ----------------boundary condition----------------- Code:
// p dimensions [0 2 -2 0 0 0 0]; internalField uniform 0; boundaryField { Agg { //type zeroGradient; type fixedValue; value uniform 0; } inlet { type zeroGradient; //type fixedValue; //value uniform 1e-3; } outlet { //type fixedFluxPressure; type fixedValue; value uniform 0; } left { type zeroGradient; } right { type zeroGradient; } frontAndBack { type empty; } } Code:
// U dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { Agg { type noSlip; //type fixedValue; //value uniform (0 0 0); } inlet { //type zeroGradient; type fixedValue; value uniform (0 -0.001 0); } outlet { type inletOutlet; //value uniform 0; inletValue uniform (0 0 0); } left { type noSlip; //type fixedValue; //value uniform (0 0 0); } right { type noSlip; //type fixedValue; //value uniform (0 0 0); } frontAndBack { type empty; } } Code:
Create time Create polyMesh for time = 0 Time = 0 Mesh stats points: 1825571 faces: 3949696 internal faces: 3216923 cells: 1085217 faces per cell: 6.60386 boundary patches: 6 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 728927 prisms: 58108 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 0 polyhedra: 298182 Breakdown of polyhedra by number of faces: faces number of cells 5 42486 6 90848 7 714 8 2476 9 90836 10 1421 11 1044 12 49539 13 100 14 394 15 18045 16 2 17 5 18 271 21 1 Checking topology... Boundary definition OK. ***Total number of faces on empty patches is not divisible by the number of cells in the mesh. Hence this mesh is not 1D or 2D. 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 100 202 ok (non-closed singly connected) outlet 100 202 ok (non-closed singly connected) left 100 202 ok (non-closed singly connected) right 100 202 ok (non-closed singly connected) frontAndBack 116838 141393 ok (non-closed singly connected) Agg 615535 661836 multiply connected (shared edge) <<Writing 153 conflicting points to set nonManifoldPoints Checking geometry... Overall domain bounding box (-100 -100 0) (0 0 1) Mesh has 2 geometric (non-empty/wedge) directions (1 1 0) Mesh has 2 solution (non-empty) directions (1 1 0) ***Number of edges not aligned with or perpendicular to non-empty directions: 1054845 <<Writing 1500798 points on non-aligned edges to set nonAlignedEdges Boundary openness (7.90745e-18 -1.62165e-17 4.64804e-15) OK. Max cell openness = 1.34114e-15 OK. Max aspect ratio = 25.2667 OK. Minimum face area = 1.05879e-22. Maximum face area = 1.00726. Face area magnitudes OK. ***Zero or negative cell volume detected. Minimum negative volume: -2.52891e-07, Number of negative volume cells: 31 <<Writing 31 zero volume cells to set zeroVolumeCells Mesh non-orthogonality Max: 179.57 average: 19.4536 *Number of severely non-orthogonal (> 70 degrees) faces: 604. ***Number of non-orthogonality errors: 20. <<Writing 624 non-orthogonal faces to set nonOrthoFaces ***Error in face pyramids: 1184 faces are incorrectly oriented. <<Writing 1171 faces with incorrect orientation to set wrongOrientedFaces ***Max skewness = 167.391, 144 highly skew faces detected which may impair the quality of the results <<Writing 144 skew faces to set skewFaces Coupled point location match (average 0) OK. Failed 5 mesh checks. End fvSolution Code:
solvers { p { solver PCG; preconditioner DIC; tolerance 1e-07; relTol 0.1; } U { solver PBiCG; preconditioner DILU; tolerance 1e-07; relTol 0; } } SIMPLE { nNonOrthogonalCorrectors 0; residualControl { p 1e-5; U 1e-5; } pRefCell 0; pRefValue 0; } relaxationFactors { p 0.15; U 0.94; } fvSchemes Code:
ddtSchemes { default steadyState; } gradSchemes { default Gauss linear; grad(p) Gauss linear; } divSchemes { default none; div(phi,U) bounded Gauss limitedLinearV 1; div((nuEff*dev2(T(grad(U))))) Gauss linear; // div(phid,p) Gauss limitedLinear 1; // div(phiU,p) Gauss linear; // div(U) Gauss linear; } laplacianSchemes { default none; laplacian(nuEff,U) Gauss linear corrected; laplacian((1|A(U)),p) Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } fluxRequired { default no; p ; } Code:
... Time = 28 DILUPBiCG: Solving for Ux, Initial residual = 0.576676, Final residual = 9.89667e-08, No Iterations 11 DILUPBiCG: Solving for Uy, Initial residual = 0.589817, Final residual = 2.19584e-08, No Iterations 12 DICPCG: Solving for p, Initial residual = 0.909012, Final residual = 258.481, No Iterations 1001 time step continuity errors : sum local = 3.31355e+20, global = 4.27055e+16, cumulative = 4.27046e+16 ExecutionTime = 300.04 s ClockTime = 306 s Time = 29 DILUPBiCG: Solving for Ux, Initial residual = 0.930444, Final residual = 3.52896e-08, No Iterations 15 DILUPBiCG: Solving for Uy, Initial residual = 0.940062, Final residual = 1.46116e-08, No Iterations 13 DICPCG: Solving for p, Initial residual = 0.987867, Final residual = 7531.96, No Iterations 1001 time step continuity errors : sum local = 4.12972e+30, global = 3.83275e+27, cumulative = 3.83275e+27 ExecutionTime = 363.53 s ClockTime = 370 s Time = 30 DILUPBiCG: Solving for Ux, Initial residual = 0.938886, Final residual = 6.16514e-08, No Iterations 32 DILUPBiCG: Solving for Uy, Initial residual = 0.905041, Final residual = 4.36084e-08, No Iterations 28 #0 Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" #2 ? in "/lib64/libc.so.6" #3 Foam::DICPreconditioner::calcReciprocalD(Foam::Field<double>&, Foam::lduMatrix const&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" #4 Foam::DICPreconditioner::DICPreconditioner(Foam::lduMatrix::solver const&, Foam::dictionary const&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" #5 Foam::lduMatrix::preconditioner::addsymMatrixConstructorToTable<Foam::DICPreconditioner>::New(Foam::lduMatrix::solver const&, Foam::dictionary const&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" #6 Foam::lduMatrix::preconditioner::New(Foam::lduMatrix::solver const&, Foam::dictionary const&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" #7 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" #8 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so" #9 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/bin/simpleFoam" #10 Foam::fvMatrix<double>::solve() in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/bin/simpleFoam" #11 ? in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/bin/simpleFoam" #12 __libc_start_main in "/lib64/libc.so.6" #13 ? in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/bin/simpleFoam" Floating point exception Thanks for your help! |
|
Tags |
simplefoam; mesh |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
icoFoam floating point exception (8) | leizhao512 | OpenFOAM Running, Solving & CFD | 7 | November 1, 2018 11:43 |
Floating point exception error after icoFoam | Deepacfd | OpenFOAM Bugs | 1 | February 10, 2017 05:35 |
A floating point exception has occurred: floating point exception [Overflow]. | starlight | STAR-CCM+ | 4 | May 4, 2016 09:08 |
A floating point exception - SEM Model | yansheng | STAR-CCM+ | 1 | April 4, 2016 04:57 |
Floating point exception from twoPhaseEulerFoam | openfoammaofnepo | OpenFOAM Running, Solving & CFD | 1 | March 19, 2016 13:56 |