[論文レビュー] Pre-Shape Calculus: Foundations and Application to Mesh Quality Optimization

Daniel Luft, Volker Schulz|arXiv (Cornell University)|Dec 16, 2020

Computational Geometry and Mesh Generation参考文献 48被引用数 4

ひとこと要約

本稿では、形状最適化にメッシュ品質最適化を統合するpre-shape微分が両方の法線方向および接線方向に作用することにより、古典的形状最適化を拡張する統一的理論枠組み「Pre-Shape Calculus」を提案する。この手法により、形状変形を防ぎつつ、最小限の計算コストで表面および体積メッシュ品質を向上させる接線勾配降下法を用いて、形状とメッシュ品質を同時に最適化できる。

ABSTRACT

Deformations of the computational mesh arising from optimization routines usually lead to decrease of mesh quality or even destruction of the mesh. We propose a theoretical framework using pre-shapes to generalize classical shape optimization and calculus. We define pre-shape derivatives and derive according structure and calculus theorems. In particular, tangential directions are featured in pre-shape derivatives, in contrast to classical shape derivatives featuring only normal directions. Techniques from classical shape optimization and -calculus are shown to carry over to this framework. An optimization problem class for mesh quality is introduced, which is solvable by use of pre-shape derivatives. This class allows for simultaneous optimization of classical shape objectives and mesh quality without deteriorating the classical shape optimization solution. The new techniques are implemented and numerically tested for 2D and 3D.

研究の動機と目的

• PDE制約付き問題における解の信頼性を損なう可能性がある形状最適化中のメッシュ品質の劣化を解消すること。
• パラメータ化（メッシュ）をpre-shapeとして取り扱うことで、古典的形状最適化とメッシュ品質向上を統合する理論的枠組みを構築すること。
• 最適形状を維持しながら形状最適化を正則化するpre-shapeパラメータ化追跡問題のクラスを定式化すること。
• 元の形状最適化プロセスに干渉せずに、pre-shape勾配降下法を用いてメッシュ品質最適化を数値的に実装できることを可能にすること。
• 2次元および3次元において、表面および体積メッシュ品質向上の両方に対して、フレームワークの実現可能性と有効性を示すこと。

提案手法

• 参照多様体の微分同相写像としてpre-shapeを定義し、形状とメッシュパラメータ化を同時に取り扱えるようにする。
• 法線成分および接線成分を含むpre-shape微分を導入し、古典的形状微分を一般化する。
• pre-shape微分の構造定理および微分法則（物質微分の公式や正規化条件を含む）を導出する。
• 表面および体積上の所望の要素体積分布を目的関数として持つpre-shapeパラメータ化追跡問題を構築する。
• 形状変形を避けるために、pre-shape微分の接線成分のみを用いた勾配降下法を実装する。
• 局所的な接線正規直交基底を計算するためのグラム・シュミット直交化法を用い、メッシュ品質追跡のためのヤコビアン行列式を計算する。

実験結果

リサーチクエスチョン

• RQ1 統一的微分法の枠組み内で、形状最適化中のメッシュ品質劣化を形式的にどのように扱えるか？
• RQ2 法線および接線変動を含むpre-shape微分の数学的構造は何か？
• RQ3 元の形状最適化問題における最適形状を変えることなく、pre-shape微分を用いてメッシュ品質を最適化できるか？
• RQ4 pre-shape微分を数値的に実装するにはどうすればよいか？ これにより、最小限の計算コストでメッシュ品質向上を達成できるか？
• RQ5 pre-shape微分の接線成分のみで、形状を保持したまま所望のメッシュ品質に収束を導けるか、その範囲はどの程度か？

主な発見

• pre-shape微分法の枠組みは、接線変動を組み込むことで、形状変形を伴わずメッシュ品質最適化を可能にする古典的形状微分法の一般化に成功した。
• 数値実験の結果、pre-shape微分の接線成分のみを用いることで、目的関数および勾配ノルムの収束が達成され、元の球形の形状が保持された。
• 3次元球体テストケースでは、48回の反復および1256.78秒の計算時間で収束し、所望の非一様な表面要素体積分布が達成された。
• 最適形状が不変のままであることが確認され、メッシュ品質最適化の過程ですら球形幾何学が保持された。
• 実装は追加計算コストが最小限であり、形状最適化プロセスに干渉することなく、メッシュ品質最適化が達成された。
• フレームワークは任意の次元に適用可能であり、表面および体積メッシュに適用可能であり、2次元および3次元における数値的妥当性が確認された。