|
[Sponsors] |
[snappyHexMesh] SnappyHexMesh in Parallel problem |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 22, 2012, 20:00 |
SnappyHexMesh in Parallel problem
|
#1 |
New Member
Join Date: Apr 2012
Location: UK
Posts: 3
Rep Power: 14 |
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 Last edited by swifty; April 22, 2012 at 20:10. Reason: extra information |
|
April 23, 2012, 07:12 |
|
#2 |
Senior Member
cfdkid
Join Date: Mar 2009
Posts: 133
Rep Power: 17 |
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. |
|
April 23, 2012, 20:21 |
|
#3 |
New Member
Join Date: Apr 2012
Location: UK
Posts: 3
Rep Power: 14 |
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. |
|
June 25, 2012, 03:26 |
|
#4 |
Senior Member
Join Date: Dec 2011
Posts: 111
Rep Power: 19 |
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. |
|
September 26, 2012, 09:37 |
|
#5 |
Senior Member
Eloïse
Join Date: Jul 2012
Location: Trondheim, Norway
Posts: 113
Rep Power: 14 |
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 |
|
February 24, 2015, 09:31 |
|
#6 |
New Member
Join Date: Feb 2015
Posts: 18
Rep Power: 11 |
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.... |
|
September 17, 2015, 07:37 |
|
#7 | |
Member
Join Date: Dec 2014
Posts: 50
Rep Power: 11 |
Quote:
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! |
||
September 17, 2015, 12:47 |
|
#8 |
Senior Member
Eloïse
Join Date: Jul 2012
Location: Trondheim, Norway
Posts: 113
Rep Power: 14 |
Hi Harak,Try using simple decomposition method instead of scotch or pscotch for snappyHexMesh.
Regards, Eloïse |
|
September 17, 2015, 13:02 |
|
#9 | |
Member
Join Date: Dec 2014
Posts: 50
Rep Power: 11 |
Quote:
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! |
||
September 17, 2015, 13:26 |
|
#10 |
Senior Member
Eloïse
Join Date: Jul 2012
Location: Trondheim, Norway
Posts: 113
Rep Power: 14 |
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 |
|
November 6, 2015, 05:40 |
snappyHexMesh parallel run problem
|
#11 |
New Member
jaydeep
Join Date: Sep 2015
Location: Pune, Maharashtra, India
Posts: 7
Rep Power: 11 |
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 |
|
|
|
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 |