A unifying moving mesh method for curves, surfaces, and domains based on mesh equidistribution and alignment
Min Zhang, Weizhang Huang
TL;DR
This work introduces a unifying moving mesh method that simultaneously handles curves, surfaces, and domains as $m$-dimensional objects in $\mathbb{R}^d$ by exploiting mesh equidistribution and alignment within a metric tensor framework. It derives a gradient-flow based MMPDE for mesh movement, provides closed-form analytical velocity expressions via scalar-by-matrix differentiation, and proves nonsingularity of the evolving mesh. The approach does not require an analytical parameterization of the geometry, instead leveraging normals and tangents (or curvature) computed from the mesh to keep points on $S$ while controlling concentration through a curvature- or identity-based metric. Numerical experiments on a range of 2D and 3D curves and surfaces demonstrate robust, singularity-free mesh movement and curvature-driven concentration patterns, illustrating practical applicability to geometric PDEs and simulations on general manifolds.
Abstract
A unifying moving mesh method is developed for general $m$-dimensional geometric objects in $d$-dimensions ($d \ge 1$ and $1\le m \le d$) including curves, surfaces, and domains. The method is based on mesh equidistribution and alignment and does not require the availability of an analytical parametric representation of the underlying geometric object. Mathematical characterizations of shape and size of $m$-simplexes and properties of corresponding edge matrices and affine mappings are derived. The equidistribution and alignment conditions are presented in a unifying form for $m$-simplicial meshes. The equation for mesh movement is defined based on the moving mesh PDE approach, and suitable projection of the nodal mesh velocities is employed to ensure the mesh points stay on the underlying geometric object. The analytical expression for the mesh velocities is obtained in a compact matrix form. The nonsingularity of moving meshes is proved. Numerical results for curves ($m=1$) and surfaces ($m=2$) in two and three dimensions are presented to demonstrate the ability of the developed method to move mesh points without causing singularity and control their concentration.
