June 2, 2016, 19:55
|
Mesh created using Matlab
|
#1
|
Member
Zhifang Hu
Join Date: Sep 2015
Location: Brisbane
Posts: 54
Rep Power: 10
|
Hi there
I wrote some matlab codes to help me creating mesh for OpenFOAM, however, it does not work very well at this point. I am wondering if anyone could give me some suggestions based the errors showed on checkMesh.log (see below)
Quote:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.1 |
| \\ / A nd | Web: http://www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 3.1-1dd681f6e943
Exec : checkMesh
Date : Jun 01 2016
Time : 11:35:51
Host : ubuntu14.04
PID : 8716
CtrlDict : /home/zhifang/foam/foam-extend-3.1/etc/controlDict
Case : /home/zhifang/OpenFoamTests/FoamExtendTests/OLAFoam/baseWaveFlume
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create polyMesh for time = 0
From function void polyMesh::initMesh()
in file meshes/polyMesh/polyMeshInitMesh.C at line 81
Truncating neighbour list at 7 for backward compatibility
Time = 0
Mesh stats
all points: 24
live points: 24
all faces: 29
live faces: 29
internal faces: 7
cells: 6
boundary patches: 6
point zones: 0
face zones: 0
cell zones: 0
Overall number of cells of each type:
hexahedra: 6
prisms: 0
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 Area [m^2] Surface topology
front 6 12 20 ok (non-closed singly connected)
back 6 12 20 ok (non-closed singly connected)
inlet 3 8 0.3 ok (non-closed singly connected)
beach 3 8 0.3 ok (non-closed singly connected)
atmosphere 2 6 0.06 ok (non-closed singly connected)
bottom 2 6 0.06 ok (non-closed singly connected)
Checking geometry...
This is a 2-D mesh
Overall domain bounding box (0 -0.015 -5) (2 0.015 5)
Mesh (non-empty, non-wedge) directions (1 0 1)
Mesh (non-empty) directions (1 0 1)
Mesh (non-empty, non-wedge) dimensions 2
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (0 0 0) Threshold = 1e-06 OK.
***High aspect ratio cells found, Max aspect ratio: 1.72667e+10, number of cells 6 Threshold = 1000
Writing 6 cells with high aspect ratio to set highAspectRatioCells
Minumum face area = 0.03. Maximum face area = 5. Face area magnitudes OK.
***Zero or negative cell volume detected. Minimum negative volume: -0.15, Number of negative volume cells: 6
Writing 6 zero volume cells to set zeroVolumeCells
Mesh non-orthogonality Max: 180 average: 180 Threshold = 70
***Number of non-orthogonality errors: 7.
Writing 7 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 36 faces are incorrectly oriented.
Writing 29 faces with incorrect orientation to set wrongOrientedFaces
Max skewness = 3.33429e-14 OK.
Failed 4 mesh checks.
End
|
Last edited by ce73stargazer; June 3, 2016 at 00:44.
|
|
|