non orthogonal faces - numerical setup
 

Dear all,

i'm working on a multiphase flow around an half hull using the solvers LTSInterFoam and interFoam.
The mesh that i'm using is provided by a commercial mesher and then imported into OpenFOAM format. It passed all checks, but there is a warning regarding the non orthogonal faces. Here is the complete result of the checkMesh:

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

Create polyMesh for time = 0

Time = 0

Mesh stats
points: 1760011
faces: 5047515
internal faces: 4945735
cells: 1644541
boundary patches: 9
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 1558359
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 86182

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
TRANSOM 1283 1372 ok (non-closed singly connected)
TOP 2080 2187 ok (non-closed singly connected)
INLET 4394 4542 ok (non-closed singly connected)
SIDE 14000 14258 ok (non-closed singly connected)
BOTTOM 160 189 ok (non-closed singly connected)
OUTLET 4394 4542 ok (non-closed singly connected)
SYMM 18304 18985 ok (non-closed singly connected)
HULL 55525 56820 ok (non-closed singly connected)
DECK 1640 1787 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-3.3 -2.21696e-06 -5.775) (13.2 4.95 1.65)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (1.33454e-15 -7.00367e-16 -2.58466e-15) OK.
Max cell openness = 3.48471e-16 OK.
Max aspect ratio = 74.0582 OK.
Minumum face area = 3.83174e-07. Maximum face area = 0.683414. Face area magnitudes OK.
Min volume = 4.08448e-10. Max volume = 0.561563. Total volume = 606.336. Cell volumes OK.
Mesh non-orthogonality Max: 86.568 average: 11.4283
*Number of severely non-orthogonal faces: 1294.
Non-orthogonality check OK.
<<Writing 1294 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
Max skewness = 3.26948 OK.
Coupled point location match (average 0) OK.

Mesh OK.



I've attached also some images about the location of the non orthogonal faces. Hope this could help!
My question is: considering this mesh, which numerical setup (fvScheme, fvSolution) could be optimal in order to reach convergence in terms of numerical scheme and solvers?
It seems that growing up the velocity at inlet boundary lets the divergence happen.
I should be able to make simulation at high Froude number (up to 1.5), so with quite high value of velocity.

Thanks in advance to everyone for your support!

Yours sincerely,

Andrea