This tutorial is part of viskex.

Plot boundaries mesh tags of a 3D dolfinx mesh

import dolfinx.mesh
import mpi4py.MPI

import viskex

import common_dolfinx as common  # isort: skip

Generate meshes of the unit cube by dividing each edge of the cube in 6 segments, using either a tetrahedral or hexahedral mesh.

cube_tetra = dolfinx.mesh.create_unit_cube(mpi4py.MPI.COMM_WORLD, 6, 6, 6, dolfinx.mesh.CellType.tetrahedron)
cube_hexa = dolfinx.mesh.create_unit_cube(mpi4py.MPI.COMM_WORLD, 6, 6, 6, dolfinx.mesh.CellType.hexahedron)

Mark boundaries according to the $(x, y, z)$ position of the vertices of the mesh:

• if $x > 0$ and $x \leq \frac{1}{3}$, assign boundary mesh tag equal to $1$;
• if $x \geq \frac{2}{3}$ and $x < 1$, assign boundary mesh tag equal to $2$;
• if $x = 0$, assign boundary mesh tag equal to $3$;
• if $x = 1$, assign boundary mesh tag equal to $4$;
• otherwise, leave the boundary mesh tag unassigned.

cube_tetra_subdomains = common.mark_subdomains(cube_tetra)
cube_hexa_subdomains = common.mark_subdomains(cube_hexa)
cube_tetra_boundaries = common.mark_boundaries(cube_tetra, cube_tetra_subdomains)
cube_hexa_boundaries = common.mark_boundaries(cube_hexa, cube_hexa_subdomains)

We plot the resulting boundaries.

viskex.dolfinx.plot_mesh_tags(cube_tetra, cube_tetra_boundaries, "boundaries")

error: XDG_RUNTIME_DIR is invalid or not set in the environment.
MESA: error: ZINK: failed to choose pdev
glx: failed to create drisw screen

viskex.dolfinx.plot_mesh_tags(cube_hexa, cube_hexa_boundaries, "boundaries")

error: XDG_RUNTIME_DIR is invalid or not set in the environment.
MESA: error: ZINK: failed to choose pdev
glx: failed to create drisw screen

We can also plot all boundaries with tag equal to $4$ (first and second plot), $2$ (third and fourth plot) or both (fifth and sixth plot), which will be displayed with colors. The rest of the mesh facets are still included in the plot, but are colored in gray.

viskex.dolfinx.plot_mesh_tags(
    cube_tetra, cube_tetra_boundaries, "boundaries_4", viskex.utils.values_in([4]))

error: XDG_RUNTIME_DIR is invalid or not set in the environment.
MESA: error: ZINK: failed to choose pdev
glx: failed to create drisw screen

viskex.dolfinx.plot_mesh_tags(
    cube_hexa, cube_hexa_boundaries, "boundaries_4", viskex.utils.values_in([4]))

error: XDG_RUNTIME_DIR is invalid or not set in the environment.
MESA: error: ZINK: failed to choose pdev
glx: failed to create drisw screen

viskex.dolfinx.plot_mesh_tags(
    cube_tetra, cube_tetra_boundaries, "boundaries_2", viskex.utils.values_in([2]))

error: XDG_RUNTIME_DIR is invalid or not set in the environment.
MESA: error: ZINK: failed to choose pdev
glx: failed to create drisw screen

viskex.dolfinx.plot_mesh_tags(
    cube_hexa, cube_hexa_boundaries, "boundaries_2", viskex.utils.values_in([2]))

error: XDG_RUNTIME_DIR is invalid or not set in the environment.
MESA: error: ZINK: failed to choose pdev
glx: failed to create drisw screen

viskex.dolfinx.plot_mesh_tags(
    cube_tetra, cube_tetra_boundaries, "boundaries_24", viskex.utils.values_in([2, 4]))

error: XDG_RUNTIME_DIR is invalid or not set in the environment.
MESA: error: ZINK: failed to choose pdev
glx: failed to create drisw screen

viskex.dolfinx.plot_mesh_tags(
    cube_hexa, cube_hexa_boundaries, "boundaries_24", viskex.utils.values_in([2, 4]))

error: XDG_RUNTIME_DIR is invalid or not set in the environment.

MESA: error: ZINK: failed to choose pdev
glx: failed to create drisw screen

What's next?

• next step: plot a scalar field
• repeat this step with firedrake