View Issue Details

IDProjectCategoryView StatusLast Update
0001229OpenFOAMBugpublic2014-03-20 09:58
Reporteruser889Assigned Touser4 
PrioritynormalSeverityfeatureReproducibilityalways
Status resolvedResolutionunable to reproduce 
PlatformLinuxOSCent OS 6.3OS Version(please specify)
Summary0001229: Baffle creation does not work in snappyHexMeshDict (entry faceType baffle;)
DescriptionHello,

this Feature http://www.openfoam.org/version2.2.0/snappyHexMesh.php Baffle and Boundary Creation (the entry faceType baffle; in snappyHexMeshDict) seems not to work with the Version Build : 2.3.x-f5bc84928a6d. With Build: 2.2.x-f9a78f7f7596 it works fine.
You can test it maybe with the propeller tutorial.
If I use version Build: 2.2.x-f9a78f7f7596 and make a checkMesh after meshing I get 2 Regions as I aspected. If I use Build : 2.3.x-f5bc84928a6d I get only 1 Region in checkMesh.log and I can’t also use the Workflow for generating AMI Faces as described in the propeller Tutorial.
Here are some abstracts of my different snappyHexMesh log files... don’t compare the numbers they are different ... as you see 2.3.x seems to recognize the keyword faceType baffle; but I don’t get my two region mesh.

Build : 2.2.x-f9a78f7f7596
Morphing phase
--------------
Converting zoned faces into baffles ...
For zone faceZ_Sliding_box_leicht_groesserAMI found patches Sliding_box_leicht_groesserAMI_AMI and Sliding_box_leicht_groesserAMI_AMI_slave
Edge intersection testing:
    Number of edges : 40541007
    Number of edges to retest : 7085820
    Number of intersected edges : 4162658
Created 1021381 baffles in = 15.67 s


Surface : Sliding_box_leicht_groesserAMI
    faces to become baffle : 2042762
    points to duplicate : 0

Snapping to features in 10 iterations ...
Constructing mesh displacer ...

Build : 2.3.x-f5bc84928a6d


Morphing phase
--------------

Converting zoned faces into baffles ...
For zone faceZ_Sliding_box_leicht_groesserAMI found patches Sliding_box_leicht_groesserAMI_AMI_fein and Sliding_box_leicht_groesserAMI_AMI_fein_slave
Edge intersection testing:
    Number of edges : 19252367
    Number of edges to retest : 3583678
    Number of intersected edges : 2304343
Created 517848 baffles in = 11.21 s


Surface : Sliding_box_leicht_groesserAMI
    faces to become baffle : 1035704
    points to duplicate : 0

Snapping to features in 10 iterations ...
Constructing mesh displacer ...

without keywordd faceType baffle

Build : 2.2.x-f9a78f7f7596

Morphing phase
--------------

Converting zoned faces into baffles ...
For zone faceZ_Sliding_box_leicht_groesserAMI found patches Sliding_box_leicht_groesserAMI_AMI and Sliding_box_leicht_groesserAMI_AMI_slave
Edge intersection testing:
    Number of edges : 20093094
    Number of edges to retest : 3612834
    Number of intersected edges : 2347178
Created 520934 baffles in = 10.29 s



Snapping to features in 10 iterations...
Constructing mesh displacer...
Steps To ReproducePropeller Tutorial
mesh wirh 2.2.x-f9a78f7f7596 --> checkMesh --> 2 Regions reported

Propelle Tutorial

mesh with 2.3.x-f5bc84928a6d --> checkMesh --> 1 Region reported
TagsNo tags attached.

Activities

user889

2014-03-19 09:13

  ~0002959

sorry another remark ... maybe it is not a bug but I investigate the following:

meshed and calculated always with Build : 2.2.x-f9a78f7f7596

1. mesh snappy with entry faceType baffle; mesh has around 7 mio cells
2. checkMesh --> two Regions

3. not sure if it is right make one try with mergeOrSplitBaffles -split and one without --> must I do split Baffles?????


4. calculate both (with split Baffles and without) with pimpleDyMFoam on 16 cores per timestep ca. 43seconds
min AMI weight around 1 but max AMI weight arround 2 !? (for both (with and without splitBaffles see below)

Other procedure (also same 2.2.x version and workstation)
1. mesh same with snappy without entry faceType baffle; --> exactly the same size as the above grid (7 mio cells)
2. checkMesh --> one Region
3. procedure as explain in propeller tutorial so
3.1. create AMI Faces with topoSet
3.2. removeRedudantZones
3.3. createBaffles
3.4. split Baffels
4. checkMesh --> two Region
5. calculate with pimpleDyMFoam on 16 cores per timestep ca. 36 seconds

without keyword faceType baffle; in snappyHexMeshDict
 
Courant Number mean: 8.64905971732e-05 max: 2.59737784981
Time = 5e-05

solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 5e-05 transformation: ((-0.00779796735578 0.000817395792011 5.55111512313e-17) (0.999988144161 (0 0 -0.00486944938111)))
AMI: Creating addressing and weights between 255196 source faces and 255196 target faces
AMI: Patch source weights min/max/average = 0.9871475696, 1.07806140172, 1.00033048769
AMI: Patch target weights min/max/average = 0.989115709357, 1.07463380688, 1.00032654688

Mesh stats
    points: 8869498
    faces: 23829621
    internal faces: 22296551
    cells: 7501199
    faces per cell: 6.14917321884
    boundary patches: 103
    point zones: 0
    face zones: 3
    cell zones: 1

Overall number of cells of each type:
    hexahedra: 6871916
    prisms: 134330
    wedges: 0
    pyramids: 0
    tet wedges: 1835
    tetrahedra: 16
    polyhedra: 493102

Checking geometry...
    Overall domain bounding box (-1.00609661216 -1.69330000877 -0.417299985886) (1.02999997139 1.20097994804 1.18283840849)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (-3.70160124945e-15 3.21047477424e-15 -5.01686200307e-15) OK.
    Max cell openness = 4.25289656852e-16 OK.
    Max aspect ratio = 11.7877514808 OK.
    Minimum face area = 1.37630026846e-09. Maximum face area = 0.00038455321401. Face area magnitudes OK.
    Min volume = 8.45801299421e-12. Max volume = 5.85978575478e-06. Total volume = 5.53804682654. Cell volumes OK.
    Mesh non-orthogonality Max: 64.9593120247 average: 7.84795082726
    Non-orthogonality check OK.
    Face pyramids OK.
 ***Max skewness = 18.6291767531, 145 highly skew faces detected which may impair the quality of the results
  <<Writing 145 skew faces to set skewFaces
    Coupled point location match (average 0) OK.

Failed 1 mesh checks.

End


with keyword faceType baffle; in snappyHexMeshDict

Courant Number mean: 8.64867181354e-05 max: 2.61435366216
Time = 5e-05

solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 5e-05 transformation: ((-0.00779796735578 0.000817395792011 5.55111512313e-17) (0.999988144161 (0 0 -0.00486944938111)))
AMI: Creating addressing and weights between 255205 source faces and 255205 target faces
AMI: Patch source weights min/max/average = 0.986455891728, 2.18858144453, 1.0002513297
AMI: Patch target weights min/max/average = 0.987297347862, 2.06349439625, 1.00024762667

Mesh stats
    points: 8835046
    faces: 24033593
    internal faces: 22053960
    cells: 7501199
    faces per cell: 6.1440248419
    boundary patches: 106
    point zones: 0
    face zones: 1
    cell zones: 1

Overall number of cells of each type:
    hexahedra: 6840674
    prisms: 161282
    wedges: 0
    pyramids: 0
    tet wedges: 2250
    tetrahedra: 16
    polyhedra: 496977

Checking geometry...
    Overall domain bounding box (-1.00609661216 -1.69330000877 -0.417299985886) (1.02999997139 1.20097994804 1.18283840849)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (-9.75893179571e-17 8.4309853923e-15 -2.406533543e-16) OK.
    Max cell openness = 4.37988882248e-16 OK.
    Max aspect ratio = 15.0406628612 OK.
    Minimum face area = 1.37641030639e-09. Maximum face area = 0.00038455321401. Face area magnitudes OK.
    Min volume = 8.45734875078e-12. Max volume = 5.85978575478e-06. Total volume = 5.53804742597. Cell volumes OK.
    Mesh non-orthogonality Max: 64.9586414022 average: 7.85492136799
    Non-orthogonality check OK.
    Face pyramids OK.
 ***Max skewness = 18.6394512432, 147 highly skew faces detected which may impair the quality of the results
  <<Writing 147 skew faces to set skewFaces
    Coupled point location match (average 0) OK.

Failed 1 mesh checks.

End

Finalising parallel run

user4

2014-03-19 12:17

  ~0002961

I am running 23x with

   faceType boundary;

and get 2 regions. Same when running with

   faceType baffle;


If you use 'faceType internal' you'll have to duplicate the faces ('createBaffles') and points, if you use 'baffle' you'll only need to duplicate the points (mergeOrSplitBaffles -split), if 'boundary' this is already done for you.

user889

2014-03-19 15:33

  ~0002963

hmmm ... i can not reproduce it ... I have the one Case were I always have the error ... in other Cases now not but I can't see the difference in the Moment ...

propeller Case works fine whether 2.2.x or 2.3.0 or 2.3.x ...

BUT can you update tutorial Case propeller ... propellerTip.obj is missing and also in snappyHexMeshDict there stands propellerTip.obj.gz ... I copied it from 2.2.0

Then the Allruns are running perfectly!

Thank you very much for your time

user4

2014-03-20 09:58

  ~0002965

The propellerTip.obj.gz is stored in the tutorials/resources directory and gets copied by the ./Allrun.pre script.

I'll close this bug - if you encounter problems open a new one.

Issue History

Date Modified Username Field Change
2014-03-19 07:30 user889 New Issue
2014-03-19 09:13 user889 Note Added: 0002959
2014-03-19 12:17 user4 Note Added: 0002961
2014-03-19 15:33 user889 Note Added: 0002963
2014-03-20 09:58 user4 Note Added: 0002965
2014-03-20 09:58 user4 Status new => resolved
2014-03-20 09:58 user4 Fixed in Version => 2.3.x
2014-03-20 09:58 user4 Resolution open => unable to reproduce
2014-03-20 09:58 user4 Assigned To => user4