为什么直边、各向同性网格不够

三类场景会迅速暴露基础网格的局限：

• 高阶 FEM/DG：解空间是高次多项式，直边几何误差可能成为主误差；
• 边界层/激波：不同方向尺度相差几个数量级；
• 高曲率 CAD：粗直边近似会错过真实边界与法向。

高阶曲边单元

$p$ 阶等参单元同时用高阶形函数表达几何和解：

$$x(\xi)=\sum_{i=1}^{N_p}X_iN_i^{(p)}(\xi), \qquad u_h(\xi)=\sum_{i=1}^{N_p}U_iN_i^{(p)}(\xi).$$

把边界中间节点投影到 CAD 后，内部 Jacobian 可能变负。高阶网格生成因此至少包含：

1. 生成有效直边拓扑网格；
2. 将边界高阶节点贴到几何；
3. 通过弹性、调和、优化或 PDE 变形传播到内部；
4. 以 Jacobian 下界为约束修复无效单元；
5. 用高阶积分验证几何和求解误差。

各向异性度量的单位单元思想

在每点给 SPD 矩阵$M(x)$。边$e$ 的度量长度

$$\ell_M(e)=\int_0^1\sqrt{e^TM(x(t))e}\,dt.$$

理想网格希望所有边在度量空间中长度约为 1，并让单元在变换$y=M^{1/2}x$ 后接近正则。二维度量单元面积目标满足

$$|K|\sqrt{\det M_K}\approx \text{常数}.$$

这同时编码物理尺寸、拉伸方向和长宽比。

从 Hessian 到度量并不唯一

若只近似标量场$u$，可从恢复 Hessian$H(u)$ 构造$M$；但实际问题还会遇到：

• 多变量误差：速度、压力、温度的度量如何相交；
• 非光滑解：Hessian 恢复噪声很大；
• 几何度量：曲率要求与 PDE 误差要求如何合并；
• 复杂度约束：给定节点数预算怎样归一化$M$；
• 可实现性：度量变化太快时任何网格都难以同时满足。

SPD 度量的“交集”不是逐元素取最大值，因为特征向量可能不同；应使用椭球包含或广义特征分解定义。

边界层网格

粘性流近壁区域常用棱柱/六面体层。法向第$k$ 层厚度可按几何级数

$$h_k=h_0r^k, \qquad H_n=\sum_{k=0}^{n-1}h_0r^k =h_0\frac{r^n-1}{r-1}.$$

$h_0$ 由目标$y^+$、摩擦速度和黏性尺度确定，不是任意“第一层很薄”。工程困难在凹角碰撞、窄缝层重叠、尖锐特征、层终止与和外部 tet/hex 网格过渡。

曲面度量与共形结构怎样配合

共形参数化保持局部角度，但各向异性度量需要指定方向与尺度。两者可在统一语言下理解：

$$M=\lambda R \begin{bmatrix} e^{2s}&0\\0&e^{-2s} \end{bmatrix}R^T.$$

• $\lambda$ 控制面积尺度；
• $R$ 控制主方向；
• $s$ 控制各向异性；
• $s=0$ 时退化为共形缩放。

因此研究机会不在“共形或各向异性二选一”，而在怎样让 Beltrami/准共形结构、方向场与误差度量全局兼容。

高阶有效性是连续约束

对高阶映射，要求

$$\min_{\xi\in\hat K}\det J(\xi)>0.$$

这是半无限约束。常见近似包括：

• 高阶 quadrature/采样点检测；
• 把$\det J$ 写入 Bernstein 基并用凸包性质给下界；
• barrier 能量$-\log\det J$；
• 逐步弯曲并回溯；
• 基于弹性或超弹性变形的 untangling。

采样没有发现负值，只是“未发现失败”；系数界或解析证书才更接近严格保证。

研究实验应该怎样做

至少比较四个轴：

2026 的相关信号

大连理工 2026 年公开学术报告讨论了黎曼度量驱动的 P2 曲边网格，强调路径保持、各向异性捕捉、受约束 Delaunay 和曲边交换。这反映出前沿已从“生成合法一次单元”转向“几何阶次、度量方向和有效性同时控制”。报告信息适合做线索，结论仍应回到正式论文和可复现实验核验。

自测

1. 为什么把中边节点投影到 CAD 可能翻转单元？
2. 度量空间等边在欧氏空间中可能是什么形状？
3. 采样 Jacobian 全正为何不是严格证明？

答案：边界弯曲会改变单元内部高阶映射；可能是按特征方向拉长的细长单元；采样点之间仍可能出现负值。

文献深读

度量场、误差驱动和动态重网格见重网格、各向异性与自适应文献谱系；曲边有效性、变分优化和边界层见高阶曲边网格与边界层文献谱系。