三个概念先分开

重网格 remeshing：已有离散几何上改变采样与连接，目标可能只是几何逼近和元素规则性；
各向异性网格 anisotropic mesh：单元在一个度量场中近似等边，在欧氏空间中可以很细长；
自适应 mesh adaptation：根据数值解或误差估计反复修改网格，是“求解—估计—改网格”的闭环。

三者常共享局部 split/collapse/swap/smooth 操作，但目标函数和正确性标准不同。

统一语言：黎曼度量场

令 $M(x)$ 为对称正定矩阵。边 $e$ 的度量长度可写为

\ell_M(e)=\int_0^1 \sqrt{e^\mathsf T M(x_0+te)e}\,dt.

在 $M$ 中近似单位、等边的网格，在欧氏空间会沿误差变化慢的方向拉长。若 $H$ 是数值解 Hessian，许多方法把其特征向量作为网格主方向，再对特征值做截断、归一化与复杂度控制。关键难点不是写出 $M$ ，而是：误差模型是否适合目标 PDE， $M$ 是否平滑，边界几何是否允许目标各向异性，以及生成器能否真正实现该度量。

1999–2008：从均匀采样到曲面变分重网格

级别	工作	核心思想	限制/后续问题
A	Centroidal Voronoi Tessellations, 1999	生成元等于 Voronoi 区域质心，统一采样、量化与优化	Lloyd 只保证能量下降到局部驻点
A	Interactive Geometry Remeshing, 2002	在参数域/几何信号上交互控制采样与近似	参数化畸变会进入采样误差
A	Isotropic Surface Remeshing, 2003	用边长阈值、flip、平滑和投影形成经典局部操作循环	投影与特征保护是主要工程边界
A	Anisotropic Polygonal Remeshing, 2003	用曲率/误差张量驱动方向和密度	曲率估计噪声会直接污染度量
B	Approximated Centroidal Voronoi Diagrams for Uniform Polygonal Mesh Coarsening, 2004	以近似 CVT 进行曲面简化和均匀化	目标以几何采样为主，不直接等同于 PDE 最优
B	Harmonic Functions for Quadrilateral Remeshing, 2005	用调和标量场组织四边形布局	临界点与拓扑结构决定 patch 质量
A	Recent Advances in Remeshing of Surfaces, 2008	总结参数化、Voronoi、局部操作与特征保持	作为历史地图读，不能覆盖后来的鲁棒/学习方法

2004–2016：连续网格与各向异性理论成熟

Persson–Strang 的 DistMesh 展示了另一条极简路线：把节点当作受力系统，并通过隐式几何函数约束边界。这类方法特别适合教学，因为“尺寸函数—边长力—投影”关系直接，但工业 CAD 的拓扑与特征恢复并不因此消失。

级别	工作	贡献	应核对的假设
A	Persson–Strang, A Simple Mesh Generator in MATLAB, 2004	用弹簧松弛和符号距离函数构造二维/三维示例	距离函数质量、收敛准则与拓扑复杂度
A	Loseille–Alauzet, Continuous Mesh Framework I, 2011	把离散网格极限化为连续度量对象，便于推导最优度量	连续最优如何被有限单元生成器实现
B	Optimal 3D Highly Anisotropic Mesh Adaptation, 2009	推进三维强各向异性度量构造与适配	极端特征值比下的稳健局部操作
A	Three-Dimensional Adaptive Domain Remeshing, 2014	将 level set、移动边界与表面/体重网格联成闭环	几何重构误差和物理量转移误差
A	A Decade of Progress on Anisotropic Mesh Adaptation, 2016	系统总结度量、误差模型、局部操作与工业瓶颈	综述中的理论模型与具体软件实现差异

自适应闭环真正多了什么

flowchart LR
    A[当前网格] --> B[PDE 求解]
    B --> C[误差/伴随/Hessian 估计]
    C --> D[构造并归一化度量 M]
    D --> E[split collapse swap smooth]
    E --> F[边界投影与有效性检查]
    F --> G[场量保守转移]
    G --> A

单看最终网格会漏掉至少三个误差源：估计器误差、网格实现度量的误差、解在两张网格之间转移的误差。移动界面还会多出几何演化误差。

软件和近年可复现路线

级别	入口	适合研究什么	注意事项
C	MMG publications	三角形/四面体的各向同性与各向异性局部重网格	逐项记录 Hausdorff、hmin/hmax、gradation 与 level-set 参数
C	MMG2D metric adaptation tutorial	从标量/张量度量到实际网格	输出“看起来细长”不是度量符合度的证明
B	Tetrahedral Remeshing in Large-Scale Simulation, 2022	大规模局部修改、并行化和数值模拟耦合	分区接口与可重复性
A	Anisotropic hp-Mesh Adaptation Methods, 2022	把 $h$ 、 $p$ 和方向选择放进统一误差框架	需要有限元误差估计基础
B	Goal-Oriented Mesh Adaptation in Firedrake, 2025	在现代 PDE 软件中连接伴随误差与自适应	软件栈、自动微分和具体方程假设

论文比较不能只报最小角

各向异性网格故意在欧氏空间制造小角或高长宽比。评价应至少包括：

Q_M(K),\qquad \|M-M_h\|,\qquad \eta(u_h),\qquad N,\qquad t_{mesh},\qquad t_{solve}.

其中 $Q_M$ 是度量空间质量， $\eta$ 是误差估计， $N$ 是复杂度预算。公平比较要求相同自由度或相同总成本，而不是一张网格固定单元数、另一张任意加密。

与计算共形几何的交叉点

各向异性度量给出“局部应该多长、多宽、朝哪里”；共形/方向场方法给出“全局坐标线怎样闭合、在哪里允许奇异”。真正困难的是兼容：一个由 PDE Hessian 推出的方向场，未必满足可积性、整数周期或 Abel–Jacobi 条件。

可研究的问题包括：

将误差度量投影到可积分的交叉场/微分空间时，损失多少最优性？
怎样把度量复杂度约束写入 quad/hex 的整数变量？
奇异点位置应由几何曲率、PDE 误差还是两者共同决定？
自适应修改后怎样保持无缝参数化和特征对齐？

建议复现实验

选一个解析解已知、Hessian 强各向异性的二维 Poisson 问题：

用均匀三角网格作基线；
从精确 Hessian 构造理想度量，再与恢复 Hessian 比较；
用 MMG 适配到相同顶点预算；
统计度量边长分布、插值误差、有限元误差和求解时间；
改变度量截断比，观察极端各向异性何时让生成器或求解器失稳。

返回文献地图 · 对应教学正文