underworldcode / underworld3

https://underworldcode.github.io/underworld3/
Other
21 stars 10 forks source link

Identify inner/outer surface tags of spherical shell in gmsh #182

Closed gthyagi closed 7 months ago

gthyagi commented 7 months ago

In this line https://github.com/underworldcode/underworld3/blob/c953b6fde0ea8526cba87b4d2a720032c9ce5bcb/src/underworld3/meshing.py#L570

The default surface tags are not (1, 2), where 1 designates the outer surface and 2 designates the inner surface. However, the current generation process has resulted with 2 assigned to the inner surface and 3 to the outer surface.

How to check inner/outer surface?

gmsh.model.get_bounding_box(2, 2)

Output: (-1.2200001, -1.2200001, -1.2200001, 1.2200001, 1.2200001, 1.2200001)

We can use above value and model inner/outer radius to identity inner/outer surfaces.

gthyagi commented 7 months ago

refer to #184