[论文解读] Numerical Methods for Biomembranes: conforming subdivision methods versus non-conforming PL methods

Jingmin Chen, Thomas Yu+4 more|arXiv (Cornell University)|Jan 28, 2019

Advanced Numerical Analysis Techniques参考文献 64被引用 4

一句话总结

本文比較了在生物膜建模中用於求解幾何變分問題（特別是Willmore能量最小化）的符合子結構曲面（SS）方法與非符合分段線性（PL）方法。論文證明基於離散Willmore能量的PL方法無法收斂至平滑的Willmore最小化解，而基於SS的方法則因高階精度與共形性而成功收斂，進而提出一種基於共形幾何的PL方法正則化方法。

ABSTRACT

The Canham-Helfrich-Evans models of biomembranes consist of a family of geometric constrained variational problems. In this article, we compare two classes of numerical methods for these variational problems based on piecewise linear (PL) and subdivision surfaces (SS). Since SS methods are based on spline approximation and can be viewed as higher order versions of PL methods, one may expect that the only difference between the two methods is in the accuracy order. In this paper, we prove that a numerical method based on minimizing any one of the `PL Willmore energies' proposed in the literature would fail to converge to a solution of the continuous problem, whereas a method based on minimization of the bona fide Willmore energy, well-defined for SS but not PL surfaces, succeeds. Motivated by this analysis, we propose also a regularization method for the PL method based on techniques from conformal geometry. We address a number of implementation issues crucial for the efficiency of our solver. A software package called Wmincon accompanies this article, provides parallel implementations of all the relevant geometric functionals. When combined with a standard constrained optimization solver, the geometric variational problems can then be solved numerically. To this end, we realize that some of the available optimization algorithms/solvers are capable of preserving symmetry, while others manage to break symmetry; we explore the consequences of this observation.

研究动机与目标

• 分析分段線性（PL）與子結構曲面（SS）方法在求解生物膜幾何變分問題時的收斂行為。
• 識別為何基於PL的Willmore能量最小化即使理論上一致，仍無法收斂至平滑解。
• 提出一種基於共形幾何的正則化技術，以改善PL方法的收斂性質。
• 開發並釋出一個並行軟體套件Wmincon，用於在PL與SS曲面上計算幾何泛函。
• 探討應用於Willmore型問題的優化求解器中對稱性保持與破壞的現象。

提出的方法

• 將Willmore能量形式化為連續幾何泛函，並為PL與SS曲面定義離散類比。
• 證明任何一致的PL Willmore能量最小化均無法收斂至平滑Willmore最小化解，原因在於缺乏共形性與不充分的正則性。
• 提出一種基於共形參量化與離散共形幾何的PL方法正則化，以穩定收斂行為。
• 在名為Wmincon的並行C++軟體套件中實現幾何泛函（Willmore、面積、體積、平均曲率）。
• 使用標準約束優化求解器搭配Wmincon，數值求解Canham-Helfrich-Evans問題。
• 分析優化求解器的對稱性行為，區分對稱性保持與對稱性破壞的算法。

实验结果

研究问题

• RQ1 為何基於PL的離散Willmore能量最小化即使理論一致，仍無法收斂至平滑Willmore最小化解？
• RQ2 能否基於共形幾何的正則化恢復PL方法在Willmore能量最小化中的收斂性？
• RQ3 不同優化求解器如何影響幾何變分問題中解的對稱性？
• RQ4 共形性與高階逼近在確保收斂至平滑解中扮演何種角色？
• RQ5 所提出的正則化技術能否推廣至其他非符合有限元方法在幾何PDE中的應用？

主要发现

• 任何定義於PL曲面上的一致離散Willmore能量均無法收斂至平滑Willmore最小化解，如命題3.9與3.12所證明。
• 子結構曲面方法因其高階逼近與共形性，成功收斂至平滑Willmore最小化解，而PL方法則無法做到。
• 所提出的基於共形幾何的PL方法正則化技術恢復了收斂行為，並提升了數值穩定性。
• 軟體套件Wmincon可有效實現PL與SS曲面上幾何泛函的並行計算，促進變分問題的數值求解。
• 部分優化求解器能保持對稱性，收斂至對稱最小化解；而其他求解器則破壞對稱性，收斂至非對稱解，影響數值結果的可靠性。
• 理論分析表明，即使對於虧格0的情況，PL曲面的空間過於龐大且不規則，無法支持一致收斂至平滑Willmore最小化解。