[Loekatoni]

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