[論文レビュー] Analysis of the first variation and a numerical gradient flow for integral Menger curvature

Tobias Hermes|arXiv (Cornell University)|Aug 27, 2014

Geometric Analysis and Curvature Flows参考文献 59被引用数 7

ひとこと要約

本稿は閉曲線に対する積分Menger曲率の最初の変種を分析し、三角多項式近似と適応的時間刻み法を用いた数値勾配流れを開発する。円がスケーリングされたエネルギーの定常点であることを証明し、広範な数値実験を通じて、特に大きな $ p $ に対してエネルギーを効率的に最小化すること、自己貫通を回避し、対称的かつ理想に近い配置に収束することを示している。

ABSTRACT

In this thesis, we consider the knot energy "integral Menger curvature" which is the triple integral over the inverse of the classic circumradius of three distinct points on the given knot to the power $p\in [2,\infty)$. We prove the existence of the first variation for a subset of a certain fractional Sobolev space if $p>3$ and for a subset of a certain Hölder space otherwise. We also discuss how fractional Sobolev and Hölder spaces can be generalised for $2π$-periodic, closed curves. Since this energy is not invariant under scaling, we additionally consider a rescaled version of the energy, where we take the energy to the power one over $p$ and multiply by the length of the curve to a certain power. We prove that a circle is at least a stationary point of the rescaled energy. Furthermore, we show that in general a functional with a scaling behaviour like $E(rγ)=r^αE(γ)$, $α\in\mathbb{R}$ cannot have stationary points unless $α=0$. Consequently "integral Menger curvature" for $p eq 3$ can be used as a subsidiary condition for a Lagrange-multiplier rule. We consider a numerical gradient flow for the rescaled energy. For this purpose we use trigonometric polynomials to approximate the knots and the trapezoidal rule for numerical integration, which is very efficient in this case. Moreover, we derive a suitable representation of the first variation. We present an algorithm for the adaptive choice of time step size and for the redistribution of the Fourier coefficients. After discussing the full discretization we present a wide collection of example flows.

研究の動機と目的

• 分数次ソボレフ空間および Hölder 空間において、$ p > 3 $ および $ p \leq 3 $ のそれぞれに対して、積分Menger曲率の最初の変種の存在を確立すること。
• エネルギーのスケーリング挙動を分析し、唯一 $ \alpha = 0 $ のときのみ関数的エネルギーが定常点を持つことを示し、ラグランジュ乗数枠組みにおける制約としての妥当性を裏付けること。
• 三角多項式と台形則を用いて、スケーリングされたエネルギーのための数値勾配流れを開発すること。
• 収束性と安定性を向上させるために、適応的時間刻み法およびフーリエ係数再配分アルゴリズムを設計すること。
• さまざまな Knot 型（無結、三葉結び、8の字結びなど）が異なる $ p $-値のもとでどのように進化するかを数値的に示すこと。

提案手法

• 曲線上の3点の外接円半径の幾何的解析を用いて最初の変種を導出し、数値計算に適した表現を得る。
• 閉曲線のパラメータ化に三角多項式を用い、スペクトル精度とフーリエに基づく効率的計算を実現する。
• 数値積分に台形則を適用し、周期性と滑らかな関数に対する高い収束速度を活用する。
• 安定性と効率性を維持するため、エネルギーおよび曲率の変化に基づいて $ \tau $ を調整する適応的時間刻み戦略を実装する。
• 進化する曲線における解像度の維持とクラスタリングの防止を目的として、フーリエ係数再配分アルゴリズムを導入する。
• 空間的スペクトル近似、時間的刻み法、適応的制御を組み合わせた完全離散化手法を構築し、広範な数値例によって検証する。

実験結果

リサーチクエスチョン

• RQ1 分数次ソボレフ空間または Hölder 空間において、積分Menger曲率の最初の変種がどのような条件下で存在するか。
• RQ2 スケーリングされた積分Menger曲率エネルギーに定常点が存在するか。円はそのような点の一つであるか。
• RQ3 スケーリング挙動 $ E(r\gamma) = r^\alpha E(\gamma) $ が定常点の存在に与える影響は。
• RQ4 非局所的かつ滑らかでないエネルギーに対して、安定的かつ効率的な勾配流れを実現するための数値手法は何か。
• RQ5 さまざまな Knot 型が勾配流れのもとでどのように進化するか。また、$ p $ を変化させた場合、どのような配置に近づくか。

主な発見

• $ p > 3 $ の場合、最初の変種は分数次ソボレフ空間の部分集合に存在する。$ p \leq 3 $ の場合、それは Hölder 空間の部分集合に存在する。
• 円はスケーリングされたエネルギーの定常点であることが、解析的導出および数値的検証によって確認された。
• $ p = 3 $ のスケーリングエネルギーは、再配分なしで 1,000,000 ステップ後に $ 5_1 $ 結びに対して対称的配置に収束し、エネルギーは約 42.1562 に安定する。
• $ p = 50 $ の場合、エネルギーは顕著に減少し、自己貫通を回避する。$ 5_2 $ 結びは 50,000 ステップ後にエネルギー $ \approx 43.87673 $ に達する。
• 再配分アルゴリズムにより収束が加速される。$ p = 50 $ の $ 5_1 $ 結びにおいて、再配分ありでは 50,000 ステップで改善された配置に到達するが、再配分なしでは 1,000,000 ステップ以上を要する。
• 流れは初期の対称性を保存するが、丸め誤差のため破れ、特に大きな $ p $ に対してはより低いエネルギーの配置に到達する。