CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Pre-Processing

Multi Region setup

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

LinkBack Thread Tools Display Modes
Old   April 20, 2014, 13:13
Default Multi Region setup
Tony Ladd
Join Date: Aug 2013
Posts: 34
Rep Power: 4
tladd is on a distinguished road
I was wondering if someone could help me with the setup of a multiregion problem. I am trying to solve a Laplace equation in two domains, with continuity of pressure across the interface and a jump in the normal gradient (Darcy flow with two different permeabilities). I was using the conjugateCavity tutorial from foam-extend-3.0/tutorials/coupled/conjugateHeatSimpleFoam/conjugateCavity/ as a starting point and it worked without any problems.

For my own case file I have the mesh set up (passing checkMesh in both regions) and I can view it in ParaFoam - all seems OK. However if I try to view the pressure field in the slave region ("exterior") paraFoam crashes:



request for polyMesh interior from objectRegistry stefan failed
available objects of type polyMesh are


From function objectRegistry::lookupObject<Type>(const word&) const
in file lnInclude/objectRegistryTemplates.C at line 139.

FOAM aborting


The reference to interior comes from the "shadowRegion interior" statement in the boundary file for the exterior region. It does not have a matching statement anywhere because the interior region is the master mesh. In conjugateCavity the statement is "shadowRegion region0". Again there is no matching statement anywhere in the case tree. I thought perhaps this was a special name and I tried that but it did not work either. I read a couple of tutorials: the report by Magnusson and the slides by Craven from OpenFOAM 6 but they both used different case setups and parameter lists than the tutorial in foam-extend-3.0. So I am a bit stuck and would appreciate any help or hints.

I am attaching the case directory. You get the current mesh setup by typing Allrun - I don't have a solver yet.

On another note - I was looking at the multiRegionHeater tutorial as well. In that case it seems you set the whole mesh up in polymesh (rather than a master polymesh and subregion polymesh dir(s)) and then add the subregion definitions with setSet. Are these equally good alternatives or are there reasons for choosing one set up over another.

Sorry for the long post

Attached Files
File Type: gz stefan.tar.gz (2.8 KB, 6 views)
tladd is offline   Reply With Quote

Old   April 21, 2014, 09:22
Default Follow up
Tony Ladd
Join Date: Aug 2013
Posts: 34
Rep Power: 4
tladd is on a distinguished road
I got paraFoam to show the initial pressure field.

Apparently the mistake was in constant/exterior/polymesh/boundary. I had the parameter "attached=on" in the regionCouple patch (as inMagnusson's report) but in the conjugateCavity tutorial it is set to off

Not sure what the difference is. Any enlightenment?
tladd is offline   Reply With Quote


Thread Tools
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
Unstabil Simulation with chtMultiRegionFoam mbay101 OpenFOAM Running, Solving & CFD 13 December 28, 2013 14:12
Multi region meshing & recovering the original patch names fluidpath OpenFOAM Native Meshers: snappyHexMesh and Others 4 May 19, 2013 19:13
Using starToFoam clo OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 33 September 26, 2012 04:04
StarToFoam error Kart OpenFOAM Meshing & Mesh Conversion 1 February 4, 2010 05:38
Import gmsh msh to Foam adorean Open Source Meshers: Gmsh, Netgen, CGNS, ... 24 April 27, 2005 08:19

All times are GMT -4. The time now is 04:59.