CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[snappyHexMesh] SnappyHexMesh in Parallel problem

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By swifty

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 22, 2012, 20:00
Unhappy SnappyHexMesh in Parallel problem
  #1
New Member
 
Join Date: Apr 2012
Location: UK
Posts: 3
Rep Power: 14
swifty is on a distinguished road
Hello,

I am trying to run snappyHexMesh in parallel using openFOAM 2.1.0 to mesh an aerofoil. I had it working without any problems in serial first, I only have problems now I am trying in parallel. I start by running blockMesh to generate a base mesh, then I decompose the mesh using scotch for 2 processors. Then change the decompose dictionary to ptscotch and run snappyHexMesh. This fails with the following error,

[0] processorPolyPatch::calcGeometry : Writing my 6925 faces to OBJ file "/home/aswift/OpenFOAM/aswift-2.1.0/run/coord_seligFmt/aerofoil/processor0/procBoundary0to1_faces.obj"
[1] processorPolyPatch::calcGeometry : Writing my 6925 faces to OBJ file "/home/aswift/OpenFOAM/aswift-2.1.0/run/coord_seligFmt/aerofoil/processor1/procBoundary1to0_faces.obj"
[1] processorPolyPatch::calcGeometry : Dumping cell centre lines between corresponding face centres to OBJ file"/home/aswift/OpenFOAM/aswift-2.1.0/run/coord_seligFmt/aerofoil/processor1/procBoundary1to0_faceCentresConnections.obj"
[0] processorPolyPatch::calcGeometry : Dumping cell centre lines between corresponding face centres to OBJ file"/home/aswift/OpenFOAM/aswift-2.1.0/run/coord_seligFmt/aerofoil/processor0/procBoundary0to1_faceCentresConnections.obj"
[1]
[1]
[1] --> FOAM FATAL ERROR:
[1] face 597 area does not match neighbour by 0.103929% -- possible face ordering problem.
patchrocBoundary1to0 my area:1.33687e-05 neighbour area:1.33549e-05 matching tolerance:6.7218e-10
Mesh face:2108231 vertices:4((0.535348 -0.114329 0.03125) (0.539191 -0.114329 0.03125) (0.539191 -0.110848 0.03125) (0.535348 -0.110851 0.03125))
If you are certain your matching is correct you can increase the 'matchTolerance' setting in the patch dictionary in the boundary file.
Rerun with processor debug flag set for more information.
[1]
[1] From function processorPolyPatch::calcGeometry()
[1] in file meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C at line 239.
[1]
FOAM parallel run exiting
[1]

To try and solve this error I have tried adding a renumberMesh command before snappyHexMesh in my run file, but that doesn't help. I have tried changing the matchTolerance value in the processor0 and 1 folders, but the value appears to be overwritten.

I have attached the snappyHexMesh log and the files I use for the run.

I would be grateful for any help.

Regards

Swifty
Attached Files
File Type: txt log.snappyHexMesh.txt (24.5 KB, 32 views)
File Type: gz aerofoil.tar.gz (11.2 KB, 18 views)

Last edited by swifty; April 22, 2012 at 20:10. Reason: extra information
swifty is offline   Reply With Quote

Old   April 23, 2012, 07:12
Default
  #2
kid
Senior Member
 
cfdkid
Join Date: Mar 2009
Posts: 133
Rep Power: 17
kid is on a distinguished road
If snappyHexMesh runs in serial what is need to run it in parallel. use decomposePar and proceed with solver. By running snappyHexMesh in parallel you might be making multiple instances of same mesh on each processor.
If anything you want to run in parallel it should be solver and not snappyHexMesh.
kid is offline   Reply With Quote

Old   April 23, 2012, 20:21
Default
  #3
New Member
 
Join Date: Apr 2012
Location: UK
Posts: 3
Rep Power: 14
swifty is on a distinguished road
I have gone back to the basics and I am doing the grid generation in serial and then run the solver in parallel. I found a problem with the decomposition of the mesh, to make it work I had to delete ccx ccy ccz cellLevel pointLevel from 0. I also added the keyword

structured yes;

to the decomposeParDict. I can now run simpleFoam in parallel.
wht likes this.
swifty is offline   Reply With Quote

Old   June 25, 2012, 03:26
Default
  #4
Senior Member
 
Join Date: Dec 2011
Posts: 111
Rep Power: 19
haakon will become famous soon enough
I have the exactly same problem with snappyhexMesh as swifty. I try to make a mesh with sHM in parallel, and it fails with the same error message. In serial everything seems to be OK.

If anyone can help me/us with this I would highly appreciate that.
haakon is offline   Reply With Quote

Old   September 26, 2012, 09:37
Default
  #5
Senior Member
 
Eloïse
Join Date: Jul 2012
Location: Trondheim, Norway
Posts: 113
Rep Power: 14
Eloise is on a distinguished road
If you still have this issue, go to have a look at the following thread:
http://www.cfd-online.com/Forums/ope...on-method.html
Eloise is offline   Reply With Quote

Old   February 24, 2015, 09:31
Default
  #6
New Member
 
Join Date: Feb 2015
Posts: 18
Rep Power: 11
jojosaxo is on a distinguished road
Hello everyone,

I'm new to Openfoam and I'm trying to run snappyHexMesh in parallel.
As some people above I got some strange message about facing area no-matching between processor.

FOAM FATAL ERROR:
[1] face 597 area does not match neighbour by 0.103929% -- possible face ordering problem.
patchrocBoundary1to0 my area:1.33687e-05 neighbour area:1.33549e-05 matching tolerance:6.7218e-10

I'm actually doing a mesh around a Naca aifroil.

This issue actually doesn't occur with a 3D case I'm working on.

Does anyone know where this problem is coming from and how to solve it???

Thanks for your help

See Attached my file....
Attached Files
File Type: txt decomposeParDict.txt (1.1 KB, 32 views)
File Type: txt snappyHexMeshDict.txt (10.5 KB, 9 views)
jojosaxo is offline   Reply With Quote

Old   September 17, 2015, 07:37
Default
  #7
Member
 
Join Date: Dec 2014
Posts: 50
Rep Power: 11
Harak is on a distinguished road
Quote:
Originally Posted by swifty View Post
I have gone back to the basics and I am doing the grid generation in serial and then run the solver in parallel. I found a problem with the decomposition of the mesh, to make it work I had to delete ccx ccy ccz cellLevel pointLevel from 0. I also added the keyword

structured yes;

to the decomposeParDict. I can now run simpleFoam in parallel.
Hey swifty,

I now it's been a long time but are you able to explain this a little further? I would really appreciate it!
What did you delete and where did you delete it?

Thanks!
Harak is offline   Reply With Quote

Old   September 17, 2015, 12:47
Default
  #8
Senior Member
 
Eloïse
Join Date: Jul 2012
Location: Trondheim, Norway
Posts: 113
Rep Power: 14
Eloise is on a distinguished road
Hi Harak,Try using simple decomposition method instead of scotch or pscotch for snappyHexMesh.
Regards,
Eloïse
Eloise is offline   Reply With Quote

Old   September 17, 2015, 13:02
Default
  #9
Member
 
Join Date: Dec 2014
Posts: 50
Rep Power: 11
Harak is on a distinguished road
Quote:
Originally Posted by Eloise View Post
Hi Harak,Try using simple decomposition method instead of scotch or pscotch for snappyHexMesh.
Regards,
Eloïse
Hey Eloise,

thanks for your quick reply.

How would you decompose this geometry using simple?



I would like to split it in 32 part, because I've 32 cores available. I'm thinking of 4 subdomains to both left and right sides and 2 subdomains to the top...so 4*4*2=32. Is that possible or is there a better way you can provide?

Thanks a lot!
Harak is offline   Reply With Quote

Old   September 17, 2015, 13:26
Default
  #10
Senior Member
 
Eloïse
Join Date: Jul 2012
Location: Trondheim, Norway
Posts: 113
Rep Power: 14
Eloise is on a distinguished road
Tricky indeed I guess it's best to avoid having empty partitions, while trying to have equivalent number of cells per partitions... I don't know how many cells you have, but it might not be best to use as many partitions pas possible. Try first (4,4,1), just to see if the decomposition method is working. If it does, you can then vary the decomposition per direction and see what works faster for you.
Regards,
Eloïse
Eloise is offline   Reply With Quote

Old   November 6, 2015, 05:40
Default snappyHexMesh parallel run problem
  #11
New Member
 
jaydeep
Join Date: Sep 2015
Location: Pune, Maharashtra, India
Posts: 7
Rep Power: 11
jaydeepKhajure is on a distinguished road
Hi all,

I need little help with snappyHexMesh. It's a 3D geometry, placed in the middle of the domain. When I run snappyHexMesh on single processor, it goes well. But, when I run it in parallel, object I placed in refinement box or in domain goes missing. It does not include the patch. I have checked my snappyHexMeshDict, decomposeParDict and blockMeshDict with motorbike case, everything seems fine.

Can anybody look into this ?

Regards,

Jaydeep
jaydeepKhajure is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[snappyHexMesh] Running snappyHexMesh in parallel creates new time directories hconel OpenFOAM Meshing & Mesh Conversion 5 September 27, 2022 16:20
snappyHexMesh in parallel leads to segmentation fault ChrisHa OpenFOAM Pre-Processing 1 January 14, 2019 11:04
Problem with foam-extend 4.0 ggi parallel run Metikurke OpenFOAM Running, Solving & CFD 1 December 6, 2018 16:51
SnappyHexMesh OF-1.6-ext crashes on a parallel run norman1981 OpenFOAM Bugs 5 December 7, 2011 13:48
[snappyHexMesh] Parallel mesh generation using snappyHexMesh aki_yafuji OpenFOAM Meshing & Mesh Conversion 0 December 25, 2010 04:49


All times are GMT -4. The time now is 13:55.