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

[Other] preserveFaceZones constraint on a multi region mesh

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

Like Tree1Likes
  • 1 Post By ViktorKL

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 29, 2018, 05:23
Default preserveFaceZones constraint on a multi region mesh
  #1
New Member
 
Viktor Klüber
Join Date: Jan 2018
Posts: 10
Rep Power: 8
ViktorKL is on a distinguished road
Hello,

I am modelling a multi region mesh consisting of a solid and a fluid part in OpenFoam.
It is supposed to depict a duct structure.
The interfaces between solid and fluid part (INTERFf & INTERFs) are of "regionCouple" boundary type and transfered to zones by (eg. for solid region):

Code:
setSet -batch solid.setSet -region solid
setsToZones -noFlipMap -region solid
with solid.setSet:
Code:
faceSet solidFluid new patchToFace INTERFs
The same goes for the fluid region.

For using a customized conjugate solver for conform meshes on this model, owner and neighbour cells of the defined zone must be processed on the same processor.
So, I wanted to use the preserveFaceZones constraint in decomposeParDict like explained in:
https://github.com/OpenFOAM/OpenFOAM...composeParDict

I used:
Code:
 
preserveFaceZones (fluidSolid solidFluid);
in both the decomposeParDict of the fluid region and the decomposeParDict of the solid one.

When now decomposing the mesh automatically (eg. using simple, hierachical or scotch) the preseveFaceZone does not have any influence on the decomposition.
The single decomposed meshes are still inconsistently distributed.
As an example the zones for one processor are depicted on the following picture, where red shows the solid and blue the fluid zone.
For my purpose, these zones should be overlaying.



I executed the decomposition first with fe4, then of4, fe32 and later of5 without any improvement.
Unfortunately, I also did not find any instructions or similiar topic on a forum explaining the preserveFaceZones constraint in detail. Therfore I opened this thread.

Is there any possibility using an automatic method for my purpose?
Maybe I did some severe mistakes in defining the zones or executing the decomposition.
Right now I am doing the decomposition manually with the setFieldsDict.
Thus, an automized decomposition would save me a lot of time in the future.


Thank you for your help in advance.


Kind regards,
Viktor
Attached Images
File Type: jpeg preserveFaceZones.jpeg (107.2 KB, 233 views)
ViktorKL is offline   Reply With Quote

Old   October 30, 2018, 08:15
Default
  #2
Member
 
Saleh Abuhanieh
Join Date: Nov 2017
Posts: 83
Rep Power: 9
Saleh Abuhanieh is on a distinguished road
Hi Viktor,


Did you solve it?
Saleh Abuhanieh is offline   Reply With Quote

Old   October 30, 2018, 08:24
Default
  #3
New Member
 
Viktor Klüber
Join Date: Jan 2018
Posts: 10
Rep Power: 8
ViktorKL is on a distinguished road
Hey,
Unfortunately no! I am still reliant on the manual decomposition of my multi-region meshes via the setFieldDict scheme.
Best regards,
Viktor
Saleh Abuhanieh likes this.
ViktorKL is offline   Reply With Quote

Old   October 30, 2018, 09:55
Default
  #4
Member
 
Saleh Abuhanieh
Join Date: Nov 2017
Posts: 83
Rep Power: 9
Saleh Abuhanieh is on a distinguished road
Hi Vektor,


Could you please explain how you are decomposing manually?


I am facing the same problem, I opened a thread while ago:


Segmentation fault initiated from a region-to-region BC
Saleh Abuhanieh is offline   Reply With Quote

Old   May 21, 2019, 12:28
Default
  #5
New Member
 
Gaoqiang Yang
Join Date: Apr 2018
Posts: 14
Rep Power: 8
shanyeyun is on a distinguished road



Hi,

Have you solved the problem?
May I ask that how do you view mesh in each processor? After I decompose the mesh, when I tried to view in paraFoam it failed.


Thanks a lot
shanyeyun is offline   Reply With Quote

Old   May 21, 2019, 12:45
Default
  #6
New Member
 
Viktor Klüber
Join Date: Jan 2018
Posts: 10
Rep Power: 8
ViktorKL is on a distinguished road
Quote:
Originally Posted by shanyeyun View Post



Hi,

Have you solved the problem?
May I ask that how do you view mesh in each processor? After I decompose the mesh, when I tried to view in paraFoam it failed.


Thanks a lot
Unfortunately I haven't found a solution yet.
You can view single segements of the decomposed mesh by executing parafoam in the respective processor# subfolder.
ViktorKL is offline   Reply With Quote

Reply

Tags
constraints, decomposepardict, preservefacezones

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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem with chtMultiregionFoam radiation boundary condition baran_foam OpenFOAM Running, Solving & CFD 10 December 17, 2019 18:36
How I can introduce my power heat (W) in chtMultiRegionFoam? aminem OpenFOAM Pre-Processing 32 August 29, 2019 03:23
[Other] conformed FSI mesh for unstructured fluid region ashish.svm OpenFOAM Meshing & Mesh Conversion 10 August 2, 2019 09:40
[snappyHexMesh] Creating multiple multiple cell zones with snappyHexMesh - a newbie in deep water! divergence OpenFOAM Meshing & Mesh Conversion 0 January 23, 2019 05:17
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 07:20


All times are GMT -4. The time now is 05:24.