网格生成 07:计算共形几何主线
共形几何到底提供什么
共形几何不直接替你“画格子”。它先把网格的方向、奇异点和布局转译成黎曼面上的内在对象,使存在性、自由度和全局障碍能够被定理描述。
flowchart LR
A[quad / hex 目标] --> B[平坦锥度量]
B --> C[有限 holonomy]
C --> D[交叉场与测度叶状结构]
D --> E[亚纯二次/四次微分]
E --> F[除子与 Abel-Jacobi]
F --> G[自然坐标]
G --> H[轨线 / T-layout / 网格]第一层:平坦锥度量
把每个理想四边形看成单位正方形,除奇异点外曲率为零。内部价数 的总角为,曲率集中量为
因此四边网格首先对应一个带锥奇点的平坦度量。离散 Ricci/Yamabe 流可在共形类内调边长,使顶点曲率逼近给定目标;但得到平坦锥度量后,还需控制 holonomy、边界对齐和轨线有限性。
第二层:holonomy
沿闭路平行移动一个方向,回来后可能旋转。四边网格允许的旋转只能是
Gauss–Bonnet 只约束总曲率;holonomy 是沿每条闭路的全局约束。前者成立不推出后者成立。
第三层:二次与四次微分
局部写四次微分
换坐标 时
所以 是坐标无关对象。选取局部四次根并积分得到自然坐标
在 平面中,水平线和竖直线天然正交;拉回曲面后成为候选四边网格线。
零点、极点与价数
若 附近
则
绕原点一圈时,自然坐标中的总角为
因此对应四边网格价数
对应三价点, 对应五价点。对闭亏格 曲面,四次微分除子的次数满足
与组合恒等式 完全一致。
Abel–Jacobi 为什么比总次数更强
给定候选除子
次数正确只说明局部曲率账目平衡。要成为某个亚纯四次微分的除子,还需
即 是主除子。用 Abel–Jacobi 映射表述,就是其像在 Jacobian variety 的周期格中为零。直觉上:局部奇异点“电荷”总和正确后,还必须让它们在所有 个独立闭路方向上的全局周期同时闭合。
Riemann–Roch 告诉我们什么
对除子,亚纯函数空间
的维数满足
它回答的不是“某个优化器能否找到解”,而是满足极点约束的全局函数空间有多少自由度。研究前沿正在把这种空间描述用于“生成一族网格”,再在其中优化几何和仿真目标。
测度叶状结构与 Strebel 微分
二次微分
的水平轨线满足。当几乎所有非临界水平轨线闭合时,可把曲面分解成有限个圆柱区域;Strebel 微分正适合控制这种分解。图值调和映射把曲面映到度量图,等值集形成闭合叶片,从而给出圆柱分解与结构化布局。
一条理论到算法的完整链
- 检查曲面拓扑,构造同调基;
- 选择奇异点/曲线系统与目标指标;
- 计算调和 1-形式、周期矩阵或图值调和映射;
- 用 Abel–Jacobi 检查/优化除子;
- 用 Ricci 流或代数表示构造平坦度量/亚纯微分;
- 在切开曲面上积分自然坐标;
- 检查 deck transformation 与整数周期;
- 提取临界轨线、T-mesh 或格点网格;
- 优化几何并验证 Jacobian、边界和求解误差。
理论保证的边界
- 共形不保证面积均匀;
- 存在亚纯微分不保证所有轨线形成有限网格布局;
- 除子合法不保证离散算法数值稳定;
- 曲面 quad layout 可行不自动保证体内 all-hex 有效;
- 论文对特定拓扑类的定理不能直接扩为任意 CAD 装配体。
核心论文
- Surface foliation theory, 2017
- Metric-based quad meshing, 2019
- Meromorphic quartic differentials and Abel–Jacobi, 2020
- Optimizing singularities by Abel–Jacobi, 2021
自测
- 为什么 仍不够?
- 自然坐标为什么产生正交网格线?
- Ricci 流获得平坦锥度量后还缺哪些条件?
答案:它只检查次数,未检查除子的线性等价/周期;自然坐标把微分化为,横纵坐标线正交;还缺有限 holonomy、边界对齐、整数周期、有限轨线与局部单射。
文献深读
离散共形底座、Surface Foliation 与 QMG I–III 的理论依赖,以及 2026 年二次微分/特殊函数拓展,见计算共形几何与大连理工研究谱系。