[czhongrong]

Hi, my simulation is diverging. I've tried increasing nNonOrthogonalCorrectors up to 10, and also introducing residualControl in fvSolution, but nothing helps. My checkMesh log is as follows.

The boundary conditions are similar to my other case, which runs perfectly. I have only changed the velocity at the inlet and with a slightly different geometry.

May I know if the "number of severely non-orthogonal faces" is the cause? Otherwise, how can I solve the diverging problem?

Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
points: 983340
faces: 2838800
internal faces: 2729200
cells: 928000
boundary patches: 7
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 928000
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
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
inlet 1600 1701 ok (non-closed singly connected)
outlet 1600 1701 ok (non-closed singly connected)
upperWall 34400 34911 ok (non-closed singly connected)
ground_wall 33600 34101 ok (non-closed singly connected)
cylinder 15200 15390 ok (non-closed singly connected)
lateral_front 11600 12140 ok (non-closed singly connected)
lateral_back 11600 12140 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-0.625 -0.5 0) (1.875 0.5 1.5)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (2.649697203e-16 -2.956581033e-14 -6.494254509e-15) OK.
Max cell openness = 2.890069133e-14 OK.
Max aspect ratio = 728.8145848 OK.
Minumum face area = 2.451588289e-08. Maximum face area = 0.004215342807. Face area magnitudes OK.
Min volume = 4.596728041e-10. Max volume = 7.903767764e-05. Total volume = 3.745717595. Cell volumes OK.
Mesh non-orthogonality Max: 88.46293676 average: 19.90242632
*Number of severely non-orthogonal faces: 110720.
Non-orthogonality check OK.
<<Writing 110720 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
Max skewness = 3.314142436 OK.
Coupled point location match (average 0) OK.

Mesh OK.

End