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

[snappyHexMesh] Problem with boundaries with sHM in parallel running

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 24, 2019, 07:56
Default Problem with boundaries with sHM in parallel running
  #1
New Member
 
Loek
Join Date: Oct 2018
Posts: 14
Rep Power: 2
Loekatoni is on a distinguished road
Dear Foamers,


I got a question regarding using snappyHexMesh in parallel. So when I snap my mesh I create a boundary which I call 'bubble'. If I don't run in parallel, than this can be used. But now I'm running the sHM in parallel and a couple of new boundaries are introduces (the places where the mesh is cut to become split over processors). Now I got a couple of new boundary patches and these overwrite my 'bubble' boundary patch. So now if I want to run my solver in parallel the solver cannot find the bubble patchField for p or U.



The code inside the p boundary are given for the non-parallel snapped mesh and the parallel snapped mesh:


__________________________________________________ _____

The non-parallel case:


The commands I run are:
surfaceCheck constant/triSurface/bubble.stl
cp -r 0.orig 0
surfaceOrient constant/triSurface/bubble.stl "(1e10 1e10 1e10)" constant/triSurface/bubble.stl
surfaceFeatureExtract
surfaceFeatureConvert constant/triSurface/bubble.eMesh edges.vtk
blockMesh
snappyHexMesh -overwrite
extrudeMesh
-> which conclude in the boundaries for p shown below:



dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
stillWall
{
type zeroGradient;
}
movingWall
{
type zeroGradient;
}
cylinderCenter
{
type zeroGradient;
}
bubbleWall
{
type zeroGradient;
}
front
{
type empty;
}
back
{
type empty;
}
bubble
{
type zeroGradient;
}
}
_____________________________________________


The parallel case


The commands I run:
surfaceCheck constant/triSurface/bubble.stl
cp -r 0.orig 0
surfaceOrient constant/triSurface/bubble.stl "(1e10 1e10 1e10)" constant/triSurface/bubble.stl
surfaceFeatureExtract
surfaceFeatureConvert constant/triSurface/bubble.eMesh edges.vtk
blockMesh
decomposePar
mpirun -np 4 snappyHexMesh -overwrite -parallel

mpirun -np 4 extrudeMesh -parallel




-> which conclude in the boundaries for p shown below:



dimensions [0 2 -2 0 0 0 0];


internalField uniform 0;

boundaryField
{
stillWall
{
type zeroGradient;
}
movingWall
{
type zeroGradient;
}
cylinderCenter
{
type zeroGradient;
}
bubbleWall
{
type zeroGradient;
}
front
{
type empty;
}
back
{
type empty;
}
procBoundary0to1
{
type processor;
value uniform 0;
}
procBoundary0to2
{
type processor;
value uniform 0;
}
}
_________________________________________


You see in the second case the boundaryPatch for bubble is gone.

Does anyone know how to counter this problem?


Sincerely,
Loek
Loekatoni is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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
Problem with foam-extend 4.0 ggi parallel run Metikurke OpenFOAM Running, Solving & CFD 1 December 6, 2018 15:51
Error running openfoam in parallel fede32 OpenFOAM Programming & Development 5 October 4, 2018 16:38
Problem for parallel running chtMultiregionFoam Giancarlo_IngChimico OpenFOAM Programming & Development 4 December 21, 2016 03:26
Problem in Running OpenFoam in Parallel himanshu28 OpenFOAM Running, Solving & CFD 1 July 11, 2013 09:19
Problem running a parallel fluent job on local machine via mpd highhopes FLUENT 0 March 3, 2011 05:07


All times are GMT -4. The time now is 17:18.