[Lee97]

Hello, I am trying to generate a mesh for the 3D finite square wing with sharp trailing edge.

I generated mesh pointwise and exported it for OpenFOAM.

But, when I run "chekMesh -allGeometry -allTopology"

It causes 9 errors.

And, When I run simpleFoam, the iteration conducted 10 iterations, and the calculation stopped.

But, When I run in ANSYS Fluent, it runs calculate and seems to converge.

In my opinion, in ParaView, main problems(skewness, non-orthogonal) are caused near the wing surface and trailing edge.

Anyone has solutions?


It's my checkMesh

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 4.1
Exec : /opt/openfoam4/platforms/linux64GccDPInt32Opt/bin/checkMesh -allGeometry -allTopology
Date : Sep 20 2022
Time : 21:54:11
Host : "anode0"
PID : 135522
Case : /home/users/DongCheon/3D_Finite_wing/NACA0015_Wing_SA_simple
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Enabling all (cell, face, edge, point) topology checks.

Enabling all geometry checks.

Time = 0

Mesh stats
points: 1888492
faces: 7620588
internal faces: 7551414
cells: 2936697
faces per cell: 5.16635
boundary patches: 4
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 1556638
prisms: 167058
wedges: 0
pyramids: 144880
tet wedges: 0
tetrahedra: 1068121
polyhedra: 0

Checking topology...
Boundary definition OK.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Topological cell zip-up check OK.
Face-face connectivity OK.
<<Writing 451 cells with two non-boundary faces to set twoInternalFacesCells
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology Bounding box
Airfoil 43466 42108 ok (closed singly connected) (0 -0.0750026 0) (1 0.0750026 6.6)
FarField 16120 8216 ok (non-closed singly connected) (-15 -25 -23) (45 25 27)
Inlet 4802 2480 ok (non-closed singly connected) (-15 -25 -23) (-15 25 27)
Outlet 4786 2472 ok (non-closed singly connected) (45 -25 -23) (45 25 27)

Checking geometry...
Overall domain bounding box (-15 -25 -23) (45 25 27)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (5.37758e-16 4.14489e-16 -3.41606e-16) OK.
***High aspect ratio cells found, Max aspect ratio: 1.86952e+145, number of cells 255312
<<Writing 255312 cells with high aspect ratio to set highAspectRatioCells
Minimum face area = 4.26742e-11. Maximum face area = 2.8502. Face area magnitudes OK.
***Zero or negative cell volume detected. Minimum negative volume: -1.09265e-23, Number of negative volume cells: 1
<<Writing 1 zero volume cells to set zeroVolumeCells
Mesh non-orthogonality Max: 90 average: 20.8956
*Number of severely non-orthogonal (> 70 degrees) faces: 216149.
***Number of non-orthogonality errors: 1.
<<Writing 216150 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
***Max skewness = 13.7631, 3631 highly skew faces detected which may impair the quality of the results
<<Writing 3631 skew faces to set skewFaces
Coupled point location match (average 0) OK.
***Error in face tets: 156 faces with low quality or negative volume decomposition tets.
<<Writing 146 faces with low quality or negative volume decomposition tets to set lowQualityTetFaces
*Edges too small, min/max edge length = 5.99044e-06 3.09297, number too small: 631977
<<Writing 631977 points on short edges to set shortEdges
*There are 1 faces with concave angles between consecutive edges. Max concave angle = 53.5414 degrees.
<<Writing 1 faces with concave angles to set concaveFaces
Face flatness (1 = flat, 0 = butterfly) : min = 0.143127 average = 0.999959
*There are 5 faces with ratio between projected and actual area < 0.8
Minimum ratio (minimum flatness, maximum warpage) = 0.143127
<<Writing 5 warped faces to set warpedFaces
Cell determinant (wellposedness) : minimum: 0 average: 0.796436
***Cells with small determinant (< 0.001) found, number of cells: 819951
<<Writing 819951 under-determined cells to set underdeterminedCells
***Concave cells (using face planes) found, number of cells: 254
<<Writing 254 concave cells to set concaveCells
Face interpolation weight : minimum: 0 average: 0.456055
***Faces with small interpolation weight (< 0.05) found, number of faces: 547
<<Writing 547 faces with low interpolation weights to set lowWeightFaces
Face volume ratio : minimum: -1.05158 average: 0.845126
***Faces with small volume ratio (< 0.01) found, number of faces: 85
<<Writing 85 faces with low volume ratio cells to set lowVolRatioFaces

Failed 9 mesh checks.

End