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

[snappyHexMesh] Cyclic BC with snappyHexMesh crashes in multiple processors

Register Blogs Community New Posts Updated Threads Search

Like Tree4Likes
  • 3 Post By wyldckat
  • 1 Post By wyldckat

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 25, 2013, 10:05
Default Cyclic BC with snappyHexMesh crashes in multiple processors
  #1
New Member
 
Join Date: Apr 2013
Posts: 11
Rep Power: 13
jgil9 is on a distinguished road
Hi all,
I am new to Open Foam and Hope you can give me a little hand, my research requires me to model a steady state turbulent case using cyclic boundaries in parallel. Seemed to me that the case motorBike was a good start, and though it handles the cyclic BC in a single processor snappy crashes in multiple processors. I have tried preservePatchs but no luck. I have also read that cyclicAMI could be a solution... I tried and though I get past snappy it crashes in potentialFoam. I have been stuck here for some time and would really appreciate some help.
Thanks In Advanced
jgil9 is offline   Reply With Quote

Old   August 7, 2013, 08:11
Default Help with CyclicAMI
  #2
New Member
 
Join Date: Apr 2013
Posts: 11
Rep Power: 13
jgil9 is on a distinguished road
Hi Bruno,
I have been trying to Apply Cyclic or CyclicAMI boundaries to the motorBike case with no success. I a farily new to OpenFoam and for my research, I need to apply Cyclic boundaries in a decomposed domain for a k-omeaga model, therefore motorbike fit the bill except for the cyclic boundaries. I appreciate the guidance.
Sincerely
JC
jgil9 is offline   Reply With Quote

Old   August 16, 2013, 08:02
Default
  #3
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,982
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Hi JC,

I've moved these posts of yours to a single thread, to make it easier to answer and to keep the discussion in a single thread

Have a look into this thread: http://www.cfd-online.com/Forums/ope...tml#post392721 - namely post #9. There you'll see an example case on how you can add cyclic patches, after meshing with snappyHexMesh.
In other words, you first should mesh without cyclic patches and only after the mesh is complete is when you can add the cyclic patches to it.

Best regards,
Bruno
Cherif, acgnipper and arsenis like this.
__________________
wyldckat is offline   Reply With Quote

Old   August 31, 2013, 05:34
Default Thanks!
  #4
New Member
 
Join Date: Apr 2013
Posts: 11
Rep Power: 13
jgil9 is on a distinguished road
Thanks Bruno, A couple of weeks or so after I posted the initial question I stumbled on an alternative solution. I was able to run snappy with cyclic boundaries but on a single processor and only decomposing before running the solvers. I will take a look at your suggestion though!
Do you then suggest to blockMesh, topoSet, decompose, run snappy in parallel, reconstructParMesh and the createPatch?

Thanks for the help, sorry for delayed response.
JC

Last edited by jgil9; August 31, 2013 at 08:37. Reason: Adding procedure question
jgil9 is offline   Reply With Quote

Old   August 31, 2013, 09:32
Default
  #5
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,982
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Hi JC,

Quote:
Originally Posted by jgil9 View Post
Do you then suggest to blockMesh, topoSet, decompose, run snappy in parallel, reconstructParMesh and the createPatch?
If you run this command:
Code:
createPatch -help
you'll see that you can also run createPatch in parallel, which means that you can run it before reconstructParMesh

Best regards,
Bruno
__________________
wyldckat is offline   Reply With Quote

Old   September 1, 2013, 02:25
Default Case
  #6
New Member
 
Join Date: Apr 2013
Posts: 11
Rep Power: 13
jgil9 is on a distinguished road
Bruno,
Thanks for the tip , I have set up the case and ran it, but I have been getting an error in createPatch, Face 3616 specified in set is not an external face of the mesh. I believe I have followed the example case you recommended. I don't mean to impose and I and I a very great-full for the time you have already devoted to this thread but I attaching the case file and I would greatly appreciate your input.
Thanks
JC
Attached Files
File Type: gz motorBike.tar.gz (11.8 KB, 12 views)
jgil9 is offline   Reply With Quote

Old   September 1, 2013, 07:17
Default
  #7
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,982
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Hi JC,

The problem with creating the sets before using snappyHexMesh is that the locations of the face indexes will be completely lost, as shown in the attached image.
I used:
Code:
reconstructParMesh -constant
after snappyHexMesh ran in parallel.
Then ran paraFoam and chose to see only the sets.

Given how you are defining the base mesh, you have two possibilities:
  1. Either you change "blockMeshDict" to create the 6 walls as patches. Then after sHM, use createPatch to change these patches to walls.
  2. Or you can use topoSet and createPatch before sHM for creating those patches. Then after sHM, use createPatch again to change the patches to cyclics.
Best regards,
Bruno
Attached Images
File Type: jpg Screenshot from 2013-09-01 11:10:10.jpg (47.2 KB, 83 views)
Cherif likes this.
__________________
wyldckat is offline   Reply With Quote

Old   September 1, 2013, 09:28
Default
  #8
New Member
 
Join Date: Apr 2013
Posts: 11
Rep Power: 13
jgil9 is on a distinguished road
Hi Bruno,
Thanks! Still have not been able to crack this problem. Tried both suggestions but with no luck, I am attaching the case with suggestion no. 1. As per your recommendation I changed blockMeshDict and then tired to create new cyclic patches after running sHM, but have been running into a point matching problem (though I turned off pointsync, and/or relaxed the matchTolerance) I also played around with the separation vector.
Thanks JC
Attached Files
File Type: gz motorBike.tar.gz (12.0 KB, 10 views)
jgil9 is offline   Reply With Quote

Old   September 1, 2013, 14:26
Default
  #9
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,982
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Hi JC,

OK, not much time for seeing your case, so here are two examples:
Best regards,
Bruno
__________________
wyldckat is offline   Reply With Quote

Old   September 1, 2013, 18:50
Default
  #10
New Member
 
Join Date: Apr 2013
Posts: 11
Rep Power: 13
jgil9 is on a distinguished road
Hi Bruno,
I can't thank you enough for your help... really! And I know you must be an extremely busy person (I respect that and try not to waste your time... I had read and tried to use cyclicAMI before, and though that using your suggestion of createPatch might do the trick). Using ciclicAMI the run gets passed sHM and patchSummary which is great but crashed in potentialFoam with a continuity error

So only when you have time... here is the case file for cyclicAMI

A very great-full forum user!
JC
Attached Files
File Type: gz motorBikeCyclicAMI.tar.gz (12.1 KB, 20 views)

Last edited by jgil9; September 2, 2013 at 07:37.
jgil9 is offline   Reply With Quote

Old   September 7, 2013, 15:53
Default
  #11
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,982
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Hi JC,

Unfortunately something strange is going on with this mesh.
Attached is what should work, but doesn't.

I suggest that you step back a few steps and start by trying to mesh without the motorbike, because I think the bike is introducing several issues into the mesh, leading createPatch to not be able to create the cyclic patches.

Best regards,
Bruno
Attached Files
File Type: gz motorBikeCyclic_doesnt_work.tar.gz (10.2 KB, 10 views)
__________________
wyldckat 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] SnappyHexMesh crashes if nRelaxIter>0 kishore96 OpenFOAM Meshing & Mesh Conversion 11 August 31, 2017 08:42
[snappyHexMesh] Tutorial crashes: snappyHexMesh floating point exception. jasv OpenFOAM Meshing & Mesh Conversion 4 May 10, 2016 03:55
[snappyHexMesh] Using snappyHexMesh for multiple enclosed regions richard_vega OpenFOAM Meshing & Mesh Conversion 0 November 13, 2014 15:28
multiple .stl in snappyHexMesh jki OpenFOAM Pre-Processing 1 August 2, 2014 14:07
OpenFOAM static build on Cray XT5 asaijo OpenFOAM Installation 9 April 6, 2011 13:21


All times are GMT -4. The time now is 10:23.