Boundary condition for multidomain geometry

I want to solve my PDE problem on multidomain geometry. I imported mesh definition to Matlab and I tried to specify multidomain geometry using geometryFromMesh(model,nodes,elements,ElementIDToRegionID) I get the geometry with appropriate number of cells, but several faces that I don't understand where from they comes. In particular, my geometry has three cells/materials and six material interfaces (material1/material2, material1/material3, material2/material3, material1/exterior, material2/exterior and material3/exterior). However, Matlab somehow defines 39 faces.

I would like to know:

• How Matlab define those faces?
• How I can get info which nodes are in any particular face?
• How to force a Matlab to define faces reasonably for my problem (e.g. on different material interface), so that I can specify boundary conditions?

Answer 1

[Ravi Kumar]

Hi Urban,

Would you be able to share your code and data? Geometric faces are identified by analyzing how drastic the normal changes between adjacent element facets.

You can find nodes associated with a face using findNodes(model.Mesh,'region','Face',fID) function.

Regards,

Ravi

Answer 2

[Urban Simoncic]

Thank you for the answer. Function findNodes solves my current problem, but does not guarantee that slightly different problem can be solved in this way. I have defined my surfaces as a list of triangles and each such surface will have different boundary condition. Matlab creates far more faces than is the number of my surface (39 vs. 6) and I can check which faces are contained in each of my surfaces with the function findNodes. However, it may happen in some other case that one particular face in partially contained in one of my surfaces and partially in another. For this reason I would generally like to have an option to redefine faces, because otherwise I may not be able to define boundary conditions as I want.

[Ravi Kumar]

Currenlty, finding the required faces might be the only option. An example data would help me figure out if there are any workaround.

Answer 3

[Urban Simoncic]

Some my example data is attached. Anyway, if Matlab is smart enough that each face is always on single interface between two materials (not split on multiple interfaces), than finding the required faces will work for me.

[Ravi Kumar]

You might want to refine the mesh that you used in geometryFromMesh call. When I plot the geometry, there seems to a bunch of spurious faces at the bottom as shown in the figure below. These can lead to problem during analysis.