[論文レビュー] Towards Unstructured Mesh Generation Using the Inverse Poisson Problem

Guy Bunin|ArXiv.org|Feb 17, 2008

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

ひとこと要約

本稿では、平面領域における非構造的四角形メッシュの生成のための新しい手法を提示する。これは逆ポアソン（IP）問題を解くことで実現される。メッシュ生成は、不規則な頂点がポアソン方程式における点源に対応するグローバル最適化問題として定式化され、解がユーザーが定義した解像度および方向性を持つ適合メッシュ場を生成する。主な貢献は、滑らかな遷移と良好な形状のセルを有する高品質で境界に沿ったメッシュを生成する包括的なメッシュ生成スキームの確立である。

ABSTRACT

A novel approach to unstructured quadrilateral mesh generation for planar domains is presented. Away from irregular vertices, the resulting meshes have the properties of nearly conformal grids. The technique is based on a theoretical relation between the present problem, and the inverse Poisson (IP) problem with point sources. An IP algorithm is described, which constructs a point-source distribution, whose sources correspond to the irregular vertices of the mesh. Both the background theory and the IP algorithm address the global nature of the mesh generation problem. The IP algorithm is incorporated in a complete mesh generation scheme, which also includes an algorithm for creating the final mesh. Example results are presented and discussed.

研究の動機と目的

• 高品質なセル形状と滑らかな解像度遷移を保証する理論的根拠に基づいたグローバルな非構造的四角形メッシュ生成アプローチの開発。
• メッシュ生成のグローバルな性質を、メッシュ特性がポテンシャル場から生じる逆ポアソン（IP）問題に結びつけることによる定式化。
• IP解による不規則頂点の配置を組み込んだメッシュ生成パイプラインの構築。
• 境界または重要な点におけるユーザー定義の入力によって、セル品質と方向性を制御する境界に沿ったメッシュの生成。

提案手法

• メッシュ生成を逆ポアソン問題の解法として定式化する：目標となるポテンシャル場 φ が与えられたとき、∇²φ = ρ を満たす源分布 ρ を再構成する。
• メッシュ内の不規則な頂点は、ポアソン方程式における点源（δ関数）としてモデル化され、その位置と次数（k値）はIP解によって決定される。
• ポテンシャル φ は局所的メッシュ解像度の対数として解釈され、その勾配はクロスフィールドおよび測地線曲線を介してメッシュの方向性を定義する。
• 有理関数補間法を用いてIP問題を解き、境界を含むユーザー指定点における必要な解像度および方向性を近似する。
• 最終的なメッシュ生成アルゴリズムにより、連続的フィールドから有効な四角形メッシュを生成し、接続性とセル形状の品質を保証する。
• 特異点（不規則頂点）におけるスターゲオデシックを、k値から導かれる角度シフトを用いて計算し、特異点周囲のメッシュエッジの適切な整合性を確保する。

実験結果

リサーチクエスチョン

• RQ1 非構造的四角形メッシュ生成のグローバルな性質を、偏微分方程式フレームワークを用いてどのように形式化できるか。
• RQ2 逆ポアソン問題が、メッシュ内での不規則頂点の最適配置および種別を効果的に特定するために有効に利用できるか。
• RQ3 逆ポアソンフレームワーク内での点源分布を用いて、ユーザー定義の解像度および方向性を持つ適合メッシュ場をどの程度正確に近似できるか。
• RQ4 得られた連続的フィールドを、境界に沿った有効で高品質な四角形メッシュに信頼性高く変換する方法は何か。

主な発見

• 本手法は、不規則な頂点が必要な箇所にのみ配置され、それ以外は非常に規則的な領域を持つ境界に沿った四角形メッシュを効果的に生成した。
• 図6の例では、粗いメッシュ（1898セル）の最小β値が0.576、平均β値が0.945に達し、良好な形状のセルであることが示された。
• 解像度を2倍にした（7630セル）場合、最小β値は0.627に上昇し、平均は0.975に達し、解像度が細かくなるにつれてセル品質が向上した。
• 図8では、粗いメッシュの最小β値が0.506、平均が0.933であったが、細かいメッシュではそれぞれ0.562および0.958に向上し、一貫した改善が確認された。
• 本手法は異なるドメインにおいてもメッシュ品質を維持しており、β値は常に0.5以上、平均は0.95～0.98の範囲に保たれ、近似的に正方形のセル形状を示した。
• 解像度スケーリングに対して本手法は頑健であり、φに定数を加えて解像度を2倍にした場合でも源分布は変化せず、メッシュ構造が保持された。