Adaptive meshing becomes truly powerful when refinement is expressed not as a set of discrete rules, but as a continuous geometric field defined over the domain. Metric based refinement provides exactly this: a Riemannian metric tensor field that encodes the desired element size, stretching, and orientation at every point. It is the mathematical foundation that unifies curvature driven refinement, physics driven anisotropy, and error based adaptation into a single framework.

At each point \( x \) in the domain, a symmetric positive definite matrix \( M(x) \) defines a local inner product:

\[ |v|_M = \sqrt{v^T M(x) v} \]

This replaces Euclidean distance with a metric weighted distance, effectively reshaping space.

If \( M(x) \) has eigenvalues \( \lambda_1, \lambda_2, \lambda_3 \) and eigenvectors \( e_1, e_2, e_3 \), then:

• \( e_i \) gives the preferred refinement direction
• \( \lambda_i^{-1/2} \) gives the desired element size along that direction

In regions where curvature or solution gradients are high, the metric increases the local “density” of space, forcing the mesher to generate smaller elements. In smooth or low gradient regions, the metric relaxes, allowing larger elements.

The mesher’s goal becomes simple:

generate elements that are unit sized in the metric space, even if they are anisotropic in Euclidean space.

This is the core idea that turns refinement into a continuous, mathematically controlled process.

• Principal curvature
• Curvature gradients
• Feature lines and sharp edges

These ensure that the mesh conforms to the underlying shape with appropriate resolution.

• Hessians of velocity, pressure, temperature, or stress
• Boundary layer gradients
• Vorticity or strain rate tensors

These capture directional physics and produce anisotropic elements aligned with flow or stress fields.

• A posteriori error estimators
• Jump indicators
• Residual based metrics

These refine the mesh where the numerical solution is under resolved.

In practice, engineering workflows often combine geometric and physical metrics:

\[ M(x) = \alpha M_{\text{geometry}}(x) + \beta M_{\text{physics}}(x) \]

This produces meshes that respect both shape and simulation accuracy.

Once the metric field is defined, the mesher operates in a transformed space where:

• isotropic elements in metric space → become anisotropic elements in Euclidean space
• smooth metric variation → produces smooth refinement transitions
• metric eigenvectors → become element orientation
• metric eigenvalues → become element stretching and size

This is why metric based refinement naturally produces:

• elongated elements in boundary layers
• compressed elements near curvature hotspots
• directional refinement along stress paths
• smooth gradation between refinement zones

The metric field acts as a continuous control law for the mesher.

A typical workflow looks like this:

1. Compute geometric or physical signals
   (curvature, Hessians, error indicators)
2. Construct the metric field
   including normalization for target error
3. Generate or adapt the mesh
   ensuring elements are unit sized in metric space
4. Validate and iterate
   using solver feedback or updated metrics

This loop converges to a mesh that is both efficient and numerically accurate.

• Captures fillets, blends, and small radii
• Preserves sharp edges and functional geometry
• Improves stress prediction and fatigue analysis

• Produces boundary layer anisotropy
• Aligns elements with flow direction
• Ensures solver stability and accuracy

• Concentrates refinement around stress concentrations
• Supports iterative solve refine loops
• Achieves error equidistributed meshes

• Provides smooth, curvature aware surface discretization
• Avoids over refinement in flat regions

Within TheMeshProject, the metric field is the bridge between:

• geometric understanding (curvature, features, segmentation)
• numerical requirements (gradients, stresses, error)
• mesh generation (anisotropy, size control, smooth transitions)

It transforms refinement from a set of heuristics into a continuous, mathematically grounded field that governs element quality and solver performance.

This unified view is what allows TheMeshProject to produce meshes that are both visually coherent and computationally optimized.