[論文レビュー] A Two Stage CVT / Eikonal Convection Mesh Deformation Approach for Large Nodal Deformations

Stephan Schmidt|arXiv (Cornell University)|Nov 27, 2014

Computer Graphics and Visualization Techniques参考文献 8被引用数 7

ひとこと要約

本稿では、表面パラメータライゼーションに重心ボロノイタイルレーション（CVT）を、ボリュームメッシュ変形にイクォン型レイトレーシングを組み合わせることで、大規模なノード形状最適化における2段階メッシュ変形手法を提案する。この手法は、CVTステップを形状微分カーネル内に埋め込み、レイの方向をイクォン方程式で計算することで、リメッシュなしでメッシュ品質を維持し、2次元および3次元ナビエ–ストークス流れにおける最適化を成功裏に実行した。それぞれエネルギー目的関数は8.0%および32.8%低下した。

ABSTRACT

A two step mesh deformation approach for large nodal deformations, typically arising from non-parametric shape optimization, fluid-structure interaction or computer graphics, is considered. Two major difficulties, collapsed cells and an undesirable parameterization, are overcome by considering a special form of ray tracing paired with a centroid Voronoi reparameterization. The ray direction is computed by solving an Eikonal equation. With respect to the Hadamard form of the shape derivative, both steps are within the kernel of the objective and have no negative impact on the minimizer. The paper concludes with applications in 2D and 3D fluid dynamics and automatic code generation and manages to solve these problems without any remeshing. The methodology is available as a FEniCS shape optimization add-on at http://www.mathematik.uni-wuerzburg.de/~schmidt/femorph.

研究の動機と目的

• 形状最適化中の大規模ノード変形におけるメッシュ劣化と劣悪なパラメータライゼーションを解決すること。
• 流体構造連成問題および非パラメトリック形状最適化におけるメッシュ品質を維持し、リメッシュを回避すること。
• 最適性を損なわせることなく、表面メッシュ品質制御を形状微分フレームワークに統合すること。
• 複雑で大規模な変形に対しても、安定したボリュームメッシュ変形を実現するため、イクォン型レイトレーシングを用いること。
• 2次元および3次元非圧縮性ナビエ–ストークス形状最適化において、リメッシュなしで本手法の有効性を実証すること。

提案手法

• 2段階の手法：まず、重心ボロノイタイルレーション（CVT）に基づく接線補正ステップにより、メッシュの悪さを最小化することで、表面ノードの分布を改善する。
• CVTステップはハダマール形状微分のカーネル内に埋め込まれており、解の最適性に影響を与えない。
• 次に、境界から体積内へのレイの方向を計算するためにイクォン方程式を解き、安定したレイトレーシングに基づくメッシュ変形を実現する。
• 風向きをイクォン解から得た方向に従い、対流拡散方程式を解くことでボリュームメッシュ変形を達成し、境界移動の滑らかな伝播を保証する。
• 自動有限要素離散化にはFEniCSを用い、ニュートンソルバを適用し、安定化された降下法のためのヘッセ行列の低次元近似にはラプラース＝ベルトラミ作用素を用いる。
• 各形状更新ステップで射影を用いて幾何的および体積制約を強制する。

実験結果

リサーチクエスチョン

• RQ1 リメッシュなしの手法は、形状最適化における大規模ノード変形でもメッシュ品質を維持できるか？
• RQ2 形状微分の最適性を損なわせることなく、表面パラメータライゼーションをどのように改善できるか？
• RQ3 イクォン型レイトレーシングは、複雑で大規模な変形に対しても、安定したボリュームメッシュ変形を提供できるか？
• RQ4 CVTに基づく接線補正は、形状微分カーネルを保持したまま、メッシュの悪さをどの程度低減できるか？
• RQ5 本2段階手法は、リメッシュなしで2次元および3次元ナビエ–ストークス流れにおける流体エネルギー散逸の最適化を成功裏に実行できるか？

主な発見

• 3次元ではエネルギー目的関数を266.1282から244.9740に低下させ、相対低下率8.0%を達成したが、リメッシュは一切行わず。
• 2次元ではエネルギー目的関数を195.6071から131.3547に低下させ、相対低下率32.8%を達成し、収束性の高さが示された。
• CVTに基づく接線補正ステップは、メッシュの悪さを効果的に低減させ、表面パラメータライゼーションを改善したが、細胞接続性は変更しなかった。
• イクォン型レイトレーシングにより、大規模なノード変位に対しても安定的かつ滑らかなボリュームメッシュ変形が実現した。
• 最適化プロセス全体がハダマール形状微分のカーネル内に留まり、最適性が保持され、メッシュ感度の複雑さの問題を回避した。
• 本手法は2次元および3次元非圧縮性ナビエ–ストークス問題に成功裏に適用され、全過程で良好に条件付けられたボリューム要素を維持した。