Constant Mean Curvature Surfaces from Discrete Harmonic Maps

一句话总结

把”曲面有恒定平均曲率 (CMC) 当且仅当其高斯映射调和”这条经典光滑理论直接离散化，用一次 S2 值调和映射求解加一次泊松求解就能构造出离散 CMC 曲面，并以曲面法向（而非顶点位置）作为设计变量。

研究背景

• 领域现状：恒定平均曲率曲面（CMC，平均曲率处处为常数 \(H\) 的曲面）在物理、微分几何和建筑中都很有价值，因为它们兼具美学与结构优势（均匀法向载荷下各方向主应力相等，适合做充气结构或网壳）。已有的离散方法主要分两派：一派在三角网格上做体积约束下的面积最小化，另一派用离散可积系统做显式构造。
• 核心痛点：这两派各自离散了光滑理论的不同侧面。面积最小化类方法直接以顶点位置为变量，容易出现网格退化、需要重网格化，并且会破坏与输入域之间的参数化对应；可积系统类方法（如把问题转成 S3 中极小曲面）与 CMC 曲面几何的联系间接、不直观，难以做交互建模。此外，非最小 (\(H \neq 0\)) 的一般 CMC 曲面的构造、离散化和设计仍是难点。
• 本文 idea：回到 Ruh–Vilms 的经典结论——浸入曲面是 CMC 当且仅当它的高斯映射 \(N: M \to S^2\) 调和。于是把设计变量从”位置”换成”法向”：先求一个 S2 值调和映射，再对它做一次”共轭”（解泊松方程）得到常高斯曲率曲面，最后取单位法向偏移即得一对 CMC 曲面。整个流程只需要离散 Dirichlet 能量最小化 + 一次泊松求解，几乎逐字对应光滑理论。

方法

整体框架：光滑理论给出一条三步链——从 S2 上的调和映射 \(N\) 出发，先构造它的共轭映射 \(g\)（满足 \(\mathrm{d}g = N \times \star\,\mathrm{d}N\)，是一张常高斯曲率曲面），再沿法向做单位偏移 \(f_\pm = g \pm N\)，得到一对共形参数化的 CMC 曲面。本文把这条链逐环节离散到三角网格（或其对偶）上。

flowchart LR
  A["边界法向 N0"] --> B["Dirichlet 能量最小化<br/>得离散调和映射 N: V→S2"]
  B --> C["解对偶泊松方程<br/>得共轭曲面 g（常高斯曲率）"]
  C --> D["单位法向偏移 f± = A(g) ± N"]
  D --> E["一对离散 CMC 曲面"]

关键设计：

• 以高斯映射为主变量。变量是单位法向 \(N\) 而非顶点位置。这带来两个好处：一是高斯映射是尺度无关的，”先定形状再定尺寸”符合设计直觉；二是在高斯映射上施加 Dirichlet 边界条件，等价于沿边界指定 CMC 曲面的切平面，从而天然支持自由边界问题与切平面控制。构造出的曲面自动是共形参数化的，方便把设计域上的纹理、UV 直接传到结果曲面。

• 离散调和映射与共轭。离散 Dirichlet 能量取 \(E(N) = \tfrac{1}{2}\sum_{ij} w_{ij}\,\lVert N_i - N_j\rVert^2\)，其中 \(w_{ij}\) 是常规余切权。关键观察是：这里的余切权只依赖固定的输入度量、在优化中不变，比”用曲面自身诱导、随迭代变化的权重去最小化面积”要简单得多。调和条件（每个内部顶点处 \(N_i \times \sum_j w_{ij} N_j = 0\)）在几何上等价于对偶多边形逐顶点”闭合”，闭合的对偶边向量正是共轭曲面的微分 \(\varphi\)，对它做（对偶）泊松求解即得 \(g\)。由于 \(g\) 定义在对偶顶点（面）上而 \(N\) 在原始顶点上，需要用面积加权平均把二者搬到同一套组合结构上再做偏移。

• Möbius 不变的鲁棒 Dirichlet 能量。在对偶网格上，标准的对偶余切拉普拉斯在非本征 Delaunay 网格上会出现负权甚至无穷权（两三角形四点近似共圆时），数值很脆。本文从离散 Willmore 能量导出一个新的、Möbius 不变的对偶 Dirichlet 能量 \(\psi^{\text{möb}}_{ij}\)，它只依赖固定度量下的角度量，且即便在非 Delaunay 网格上也恒为非负，是最鲁棒的对偶能量选择。实验中它在对偶网格上的平均曲率逼近误差比对偶余切离散低约 20%。

• 数值求解与边界处理。调和映射用热流（Dirichlet 能量的梯度下降）求解，采用 Sobolev \(H^1\) 梯度（预条件后收敛更快），并可在临界点附近切换牛顿法加速；每步用 S2 上的指数映射把切向更新拉回球面，并用回溯线搜索保证能量下降。对边界顶点，其周围对偶边向量一般不闭合，通过补充新的边/面让其闭合以得到规整的边界曲线。曲面的形状算子用最低阶 Regge 有限元逼近，用于评估平均曲率与主曲率方向。

设计层面，作者还给出三类变分问题：切平面边界条件（含”曲面切于/正交于某平面”的对称约束）、CMC 曲面逼近给定输入（把输入的高斯映射投影到 S2 调和映射空间）、以及周期 CMC 曲面（在非单连通域上通过上同调周期向量、并用牛顿法优化共形模数来闭合周期，构造 Delaunay 结点面、Wente 环面等）。

实验结果

主实验是网格细化下的收敛性与逼近精度验证，对比三种 Dirichlet 能量离散：原始余切、对偶余切、以及本文的 Möbius 能量 \(E^{\text{möb}}\)。核心结论是平均曲率误差随细化线性收敛、共形畸变二次收敛，且 \(E^{\text{möb}}\) 在对偶网格上比对偶余切更准更稳。

量度	收敛/表现	备注
平均曲率 \(L^2\)、\(L^\infty\) 误差	随网格细化线性收敛	离散平均曲率聚集在 \(1/2\) 附近，标准差约 0.004
拟共形畸变 \(Q\) 的 \(L^2\) 误差	二次收敛	粗网格上已满足 \(Q \approx 1\)
对偶网格平均曲率误差	\(E^{\text{möb}}\) 比对偶余切低约 20%	且避免了余切权求逆带来的数值不稳

其余实验用文字概括：与 Surface Evolver（Brakke）及 Pan 等人方法做一致性对比时，因求解的边值问题不同（切平面 vs 顶点位置），通过”用本文方法生成曲面→提取边界曲线与体积→喂给对方算法”的方式对齐，结果几何一致（相对包围盒对角线的 Hausdorff 距离约 \(10^{-3}\)）；不同的是本文保持了与输入域的共形对应，而对方方法或因重网格化丢失对应、或因切向滑动累积参数化畸变。此外在高度各向异性、不规则三角化上，本文无需重网格或预处理即可重建出与规则网格一致的 CMC 曲面。

亮点与局限

• 亮点：
  • 概念极简且忠于光滑理论——只用”Dirichlet 能量最小化 + 泊松求解”就复刻了 CMC 与调和高斯映射的经典对应，可读性和可实现性都很好。
  • 以法向/高斯映射为设计变量，天然支持切平面边界条件、共形参数化、以及面积-体积方法难以触及的对称与周期 CMC 曲面（结点面、Wente 环面、空间填充的蜂窝状结构等）。
  • 提出的 Möbius 不变对偶 Dirichlet 能量在非 Delaunay 网格上恒非负，显著提升数值鲁棒性与精度。
• 局限：
  • 构造只能得到 CMC-\(\tfrac{1}{2}\) 曲面，其他常数值需靠整体缩放得到（形状与高斯映射不变），本质上是一个归一化选择。
  • 非单连通域上共轭曲面的可积性受 1-形式上同调阻碍，闭合曲面需要额外解周期闭合条件（优化共形模数），对多亏格情形更复杂。
  • 输出曲面完全继承输入三角化的质量：不规则/各向异性网格的瑕疵会被共形地保留到结果中，方法本身不改善网格质量。
  • 与既有方法无法直接同题比较（边值问题不同），只能做间接一致性验证。

延伸思考

• 把设计变量从位置换到高斯映射，是”用规范/对偶量而非原始量参数化几何”的又一例证，和用中曲率、共形因子、连接形式等做主变量的思路一脉相承；这类换元往往能把非线性约束（如体积约束、单位法向约束）转成更规整的调和/泊松问题。
• 文中留了几个开放方向：其他 Dirichlet 能量离散（均值权、Friedel 等人的重标度权）会诱导出怎样的离散 CMC 概念尚未探索；对更一般多亏格闭曲面的周期闭合、以及把该表示接入交互式建筑找形与网壳优化流程，都值得跟进。
• 与建筑几何的结合点很实在：CMC 曲面在均匀法向载荷下主应力各向同性、适合做曲率对齐网壳，本文的共形参数化和切平面控制正好对接玻璃-钢结构的面板划分与找形需求。