网格生成 07:计算共形几何主线

共形几何到底提供什么

共形几何不直接替你“画格子”。它先把网格的方向、奇异点和布局转译成黎曼面上的内在对象,使存在性、自由度和全局障碍能够被定理描述。

flowchart LR
    A[quad / hex 目标] --> B[平坦锥度量]
    B --> C[有限 holonomy]
    C --> D[交叉场与测度叶状结构]
    D --> E[亚纯二次/四次微分]
    E --> F[除子与 Abel-Jacobi]
    F --> G[自然坐标]
    G --> H[轨线 / T-layout / 网格]

第一层:平坦锥度量

把每个理想四边形看成单位正方形,除奇异点外曲率为零。内部价数vv 的总角为vπ/2v\pi/2,曲率集中量为

Kv=2π−vπ2.K_v=2\pi-v\frac{\pi}{2}.

因此四边网格首先对应一个带锥奇点的平坦度量。离散 Ricci/Yamabe 流可在共形类内调边长,使顶点曲率逼近给定目标;但得到平坦锥度量后,还需控制 holonomy、边界对齐和轨线有限性。

第二层:holonomy

沿闭路平行移动一个方向,回来后可能旋转。四边网格允许的旋转只能是

Hol(γ)∈{eikπ/2:k=0,1,2,3}.\mathrm{Hol}(\gamma)\in \{e^{ik\pi/2}:k=0,1,2,3\}.

Gauss–Bonnet 只约束总曲率;holonomy 是沿每条闭路的全局约束。前者成立不推出后者成立。

第三层:二次与四次微分

局部写四次微分

Φ=ϕ(z) dz4.\Phi=\phi(z)\,dz^4.

换坐标w=w(z)w=w(z) 时

ϕ~(w)=ϕ(z(w))(dzdw)4,\widetilde\phi(w)=\phi(z(w))\left(\frac{dz}{dw}\right)^4,

所以Φ\Phi 是坐标无关对象。选取局部四次根并积分得到自然坐标

ζ=∫Φ1/4.\zeta=\int \Phi^{1/4}.

在ζ\zeta 平面中,水平线和竖直线天然正交;拉回曲面后成为候选四边网格线。

零点、极点与价数

若z=0z=0 附近

Φ=zm dz4,\Phi=z^m\,dz^4,

则

ζ∝z1+m/4.\zeta\propto z^{1+m/4}.

绕原点一圈时,自然坐标中的总角为

2π(1+m4)=(m+4)π2.2\pi\left(1+\frac m4\right)=\frac{(m+4)\pi}{2}.

因此对应四边网格价数

v=m+4.v=m+4.

m=−1m=-1 对应三价点,m=1m=1 对应五价点。对闭亏格gg 曲面,四次微分除子的次数满足

∑imi=deg⁡(K4)=4(2g−2)=8g−8,\sum_i m_i=\deg(K^4)=4(2g-2)=8g-8,

与组合恒等式∑(vi−4)=8g−8\sum(v_i-4)=8g-8 完全一致。

Abel–Jacobi 为什么比总次数更强

给定候选除子

D=∑imipi,D=\sum_i m_i p_i,

次数正确只说明局部曲率账目平衡。要成为某个亚纯四次微分的除子,还需

D∼4K,D\sim 4K,

即D−4KD-4K 是主除子。用 Abel–Jacobi 映射表述,就是其像在 Jacobian variety 的周期格中为零。直觉上:局部奇异点“电荷”总和正确后,还必须让它们在所有2g2g 个独立闭路方向上的全局周期同时闭合。

Riemann–Roch 告诉我们什么

对除子DD,亚纯函数空间

L(D)={f≠0:(f)+D≥0}∪{0}L(D)=\{f\neq0:(f)+D\ge0\}\cup\{0\}

的维数满足

ℓ(D)−ℓ(K−D)=deg⁡D−g+1.\ell(D)-\ell(K-D)=\deg D-g+1.

它回答的不是“某个优化器能否找到解”,而是满足极点约束的全局函数空间有多少自由度。研究前沿正在把这种空间描述用于“生成一族网格”,再在其中优化几何和仿真目标。

测度叶状结构与 Strebel 微分

二次微分

q=ψ(z) dz2q=\psi(z)\,dz^2

的水平轨线满足q(γ˙2)>0q(\dot\gamma^2)>0。当几乎所有非临界水平轨线闭合时,可把曲面分解成有限个圆柱区域;Strebel 微分正适合控制这种分解。图值调和映射把曲面映到度量图,等值集形成闭合叶片,从而给出圆柱分解与结构化布局。

一条理论到算法的完整链

  1. 检查曲面拓扑,构造同调基;
  2. 选择奇异点/曲线系统与目标指标;
  3. 计算调和 1-形式、周期矩阵或图值调和映射;
  4. 用 Abel–Jacobi 检查/优化除子;
  5. 用 Ricci 流或代数表示构造平坦度量/亚纯微分;
  6. 在切开曲面上积分自然坐标;
  7. 检查 deck transformation 与整数周期;
  8. 提取临界轨线、T-mesh 或格点网格;
  9. 优化几何并验证 Jacobian、边界和求解误差。

理论保证的边界

  • 共形不保证面积均匀;
  • 存在亚纯微分不保证所有轨线形成有限网格布局;
  • 除子合法不保证离散算法数值稳定;
  • 曲面 quad layout 可行不自动保证体内 all-hex 有效;
  • 论文对特定拓扑类的定理不能直接扩为任意 CAD 装配体。

核心论文

自测

  1. 为什么∑mi=8g−8\sum m_i=8g-8 仍不够?
  2. 自然坐标为什么产生正交网格线?
  3. Ricci 流获得平坦锥度量后还缺哪些条件?

答案:它只检查次数,未检查除子的线性等价/周期;自然坐标把微分化为dζ4d\zeta^4,横纵坐标线正交;还缺有限 holonomy、边界对齐、整数周期、有限轨线与局部单射。

文献深读

离散共形底座、Surface Foliation 与 QMG I–III 的理论依赖,以及 2026 年二次微分/特殊函数拓展,见计算共形几何与大连理工研究谱系。