CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Open Source Meshers: Gmsh, Netgen, CGNS, ...

False name of internal faceZone/faceSet by Gmsh

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   November 29, 2012, 09:01
Default False name of internal faceZone/faceSet by Gmsh
  #1
Senior Member
 
Hisham's Avatar
 
Hisham El Safti
Join Date: Apr 2011
Location: Braunschweig, Germany
Posts: 247
Blog Entries: 10
Rep Power: 8
Hisham is on a distinguished road
Dear Foamers,

I try to make an internal faceZone/faceSet in a Gmsh geometry. I defined a physical surface of the internal surfaces in question. I run
Code:
gmsh -3 file.geo; gmshToFoam file.msh
. By viewing the output in paraFoam, the faceZone/faceSet exists but has a different name than the physical surface. It has the name of a cellZone that is not even adjacent to the surfaces.

Question time:
1. Did anyone face the same problem or has anyone produced internal faceZones without this problem?
2. Is there a workaround?

Best regards,
Hisham
Hisham is offline   Reply With Quote

Old   November 29, 2012, 10:25
Default
  #2
Senior Member
 
Hisham's Avatar
 
Hisham El Safti
Join Date: Apr 2011
Location: Braunschweig, Germany
Posts: 247
Blog Entries: 10
Rep Power: 8
Hisham is on a distinguished road
It seems that there is a bug in gmshToFoam. Changing:
Code:
01058         forAll(zoneFaces, zoneI)
01059         {
01060             if (zoneFaces[zoneI].size())
01061             {
01062                 label physReg = zoneToPhys[zoneI];
01063 
01064                 Map<word>::const_iterator iter = physicalNames.find(physReg);
01065 
01066                 word zoneName = "faceZone_" + name(zoneI);
01067                 if (iter != physicalNames.end())
01068                 {
01069                     zoneName = iter();
01070                 }
01071 
01072                 Info<< "Writing zone " << zoneI << " to faceZone "
01073                     << zoneName << " and faceSet"
01074                     << endl;
to:

Code:
01058         forAll(zoneFaces, zoneI)
01059         {
01060             if (zoneFaces[zoneI].size())
01061             {
01062                 label physReg = patchToPhys[zoneI];
01063 
01064                 Map<word>::const_iterator iter = physicalNames.find(physReg);
01065 
01066                 word zoneName = "faceZone_" + name(zoneI);
01067                 if (iter != physicalNames.end())
01068                 {
01069                     zoneName = iter();
01070                 }
01071 
01072                 Info<< "Writing zone " << zoneI << " to faceZone "
01073                     << zoneName << " and faceSet"
01074                     << endl;
solves the problem for me!

Regards,
Hisham
Hisham is offline   Reply With Quote

Reply

Tags
faceset, facezone, gmsh, gmshtofoam, internal face

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Internal faces from gmsh how to create patches in OpenFoam podallaire Open Source Meshers: Gmsh, Netgen, CGNS, ... 27 April 25, 2012 21:24
2D Mesh Generation Tutorial for GMSH aeroslacker Open Source Meshers: Gmsh, Netgen, CGNS, ... 12 January 19, 2012 04:52
Internal boundary problem stefanke OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 2 May 24, 2011 02:21
Import problem ARC Open Source Meshers: Gmsh, Netgen, CGNS, ... 0 February 27, 2010 11:56
CFX4.3 -build analysis form Chie Min CFX 5 July 12, 2001 23:19


All times are GMT -4. The time now is 16:48.