|
[Sponsors] |
April 1, 2011, 07:04 |
fvScheme and NonOrthogonal Correctors
|
#1 |
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 16 |
Hi all,
I have to run a simulation with a mesh of prisms and hexas. This is my checkMesh: Create polyMesh for time = 0 Time = 0 Mesh stats points: 101178 internal points: 0 faces: 259884 internal faces: 118566 cells: 69762 boundary patches: 6 point zones: 0 face zones: 0 cell zones: 1 Overall number of cells of each type: hexahedra: 29640 prisms: 40122 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 0 polyhedra: 0 Checking topology... Boundary definition 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 obstacle 1560 3120 ok (non-closed singly connected) bottom 49 100 ok (non-closed singly connected) right 68 138 ok (non-closed singly connected) top 49 100 ok (non-closed singly connected) left 68 138 ok (non-closed singly connected) defaultFaces 139524 101178 ok (non-closed singly connected) Checking geometry... Overall domain bounding box (0 0 0) (0.00447 0.0062 5e-05) Mesh (non-empty, non-wedge) directions (1 1 0) Mesh (non-empty) directions (1 1 0) All edges aligned with or perpendicular to non-empty directions. Boundary openness (1.63471e-19 -2.90979e-18 -4.90414e-20) OK. Max cell openness = 2.19738e-16 OK. Max aspect ratio = 13.3858 OK. Minumum face area = 8.1356e-11. Maximum face area = 7.04863e-09. Face area magnitudes OK. Min volume = 4.0678e-15. Max volume = 3.16076e-13. Total volume = 1.21253e-09. Cell volumes OK. Mesh non-orthogonality Max: 31.2939 average: 4.74797 Non-orthogonality check OK. Face pyramids OK. Max skewness = 0.495827 OK. Mesh OK. What are the best setting for fvScheme and non-orthogonal corrector? Is the option "limited 0.5" a good choise? Something like this: Code:
gradSchemes { default cellLimited Gauss linear 1; } laplacianSchemes { default Gauss linear limited 0.5; } snGradSchemes { default limited 0.5; } And for the divergence what i have to choose? (Iam using interFoam) divSchemes { div(rho*phi,U) ?? div(phi,alpha) ?? div(phirb,alpha) ?? } and finally for the non orthogonal corrector in fvSolution, what is best choise? Thanks for any help andrea |
|
|
|