高阶网格不是“在线性边上多放几个点”

对参考单元$\hat K$，高阶几何映射写成

$$x(\xi)=\sum_{i=1}^{n_p}N_i^{(p)}(\xi)x_i.$$

角点有效并不保证整个曲边单元有效。真正条件是

$$\det J(\xi)>0,\qquad \forall \xi\in\hat K.$$

当边界节点被投影到 CAD 曲面，内部高阶节点也会移动，$\det J$ 是一个高次多项式，可能在采样点之间变负。因此这条文献线围绕三个问题展开：怎样弯曲、怎样检查、怎样在保持几何精度时修复。

2002–2012：曲边生成从经验形变走向优化

2013–2018：untangling、变分框架与工业几何

2018–2026：工具链、局部操作与证书

有效性检查为何不能只采样

在 Bernstein–Bézier 基下，$\det J$ 可表示为控制系数组合。若所有相关 Bernstein 系数为正，可以得到全域正性的充分证书；若有负系数，不能立刻断言无效，可递归细分参数域以收紧界。

这形成三级检查：

1. 快速 quadrature/sample 检查，用于筛除明显问题；
2. Bernstein 包围与自适应细分，用于严格证书；
3. 若证书失败，进入带 barrier 的优化或局部拓扑修改。

必须区分“在若干积分点$\det J>0$”和“对整个参考单元$\det J>0$”。前者只是有限采样证据。

边界层是独立难题

CFD 边界层要求法向极细、切向较粗的棱柱/六面体层。若第一层高度为$h_0$、增长率为$r$、层数为$n$，总厚度为

$$H=h_0\frac{r^n-1}{r-1}.$$

但几何可行性还受局部曲率半径、邻近壁面距离和凹角控制。沿法向推进到达 medial axis 前，偏移曲面可能相交；在尖角附近，法向甚至不唯一。因而边界层算法必须包含碰撞检测、局部停止/合并、特征分类和体网格过渡。

补充近年入口：

• Robust Boundary Layer Generation for Complex Marine Geometries, 2023
• Anisotropic Boundary-Layer Surface Meshing, 2026

与求解器的真实联系

高阶几何并不必然提高解精度。至少要共同考虑：

$$\text{几何阶次 }p_g, \quad \text{解空间阶次 }p_u, \quad \text{积分阶次 }p_q, \quad \kappa(J), \quad \min \det J.$$

几何映射畸变会放大梯度、影响积分误差和矩阵条件数。论文比较若只报最小 Jacobian，却不报 PDE 误差、积分规则和自由度，就无法证明方法对高阶仿真更好。

推荐复现实验

1. 从同一个线性三角/四面体网格生成$p=2,3,4$ 曲边网格；
2. 在圆/球和高曲率 CAD 上投影边界节点；
3. 对比弹性形变、Winslow/目标矩阵和 barrier 优化；
4. 同时做采样 Jacobian 与 Bernstein 证书；
5. 在相同自由度下运行 Poisson/弹性问题，记录误差、条件数和总时间；
6. 制造薄间隙与凹角，记录边界层首先在哪个约束上失败。

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