decomposePar -region (OF 2.1) without splitting regions in separate parts
 

I am running a simulation with a multiregion solver (a modification of chtMultiRegionFoam) and want to run some of my calculations in parallel.

I know I can use decomposePar -region regionName to decompose the different regions separately, but this gives quiet an illogical mesh division because 1 region is almost surrounded by another region but gets split over all processors regardless. (See the image of the full mesh with the 3 regions coloured and a snapshot of what part is set to processor0).

My question is: can anybody tell me how I can make sure that the decomposition is done based on position in the mesh, not based on position in the region? For this test I just used the simple decomposition method.

Hi Michiel,
I'm just facing the same problem in these days...
Have you found the solution?

Thanks in advance,

Andrea