Deforming a mesh

[dokken]

If you are only perturbing boundary nodes, you will need to apply smoothing to the interior nodes to make sure that you Get a valid deformation, see

Changing the shape in the surface normal directions

The simplest mesh deformation can be described as follows: import dolfinx from mpi4py import MPI from typing import List import numpy as np def deform_mesh(V, bcs: List[dolfinx.fem.DirichletBCMetaClass]): mesh = V.mesh u = dolfinx.fem.Function(V) dolfinx.fem.petsc.set_bc(u.vector, bcs) deformation_array = u.x.array.reshape((-1, mesh.geometry.dim)) mesh.geometry.x[:, :mesh.geometry.dim] += deformation_array mesh = dolfinx.mesh.create_unit_square(MPI.COMM_WORLD, 10, 10) …