Subdomains are currently supported only in 2D

1

Hi I am trying to generate a 3D mesh with subdomains using mshr and CSG, but when I try to implement subdomains I get the error:

*** Error:   Unable to setting subdomain.
*** Reason:  Subdomains are currently supported only in 2D.
*** Where:   This error was encountered inside CSGGeometry.cpp.
*** Process: 0
*** 
*** DOLFIN version: 2019.1.0
*** Git changeset:  74d7efe1e84d65e9433fd96c50f1d278fa3e3f3f

Is there someway to overcome this issue? or some beta version of FEniCS that set subdomains in 3D?

2

Hi, this is an issue with mshr which is not actively maintained nor well documented. For complex meshing, we suggest using external mesher software like Gmsh for instance

3

Hi, according to the python tutorial there are many applications relying on the subdomain class. Are there any elegant alternatives to replace the subdomain class in the 3D cases of (i) partial integration region/boundary, i.e. dx/ds (subdomain_markers), (ii) region-constant parameter, (iii) subdomain refinement?

4

As bleyerj said. It is only mshr that is limited, and subdomains can be used with 3D geometries in fenics as They are exported from external software