[論文レビュー] A Chebyshev-based rectangular-polar integral solver for scattering by general geometries described by non-overlapping patches

Oscar P. Bruno, Emmanuel Garza|arXiv (Cornell University)|Jul 4, 2018

Electromagnetic Scattering and Analysis参考文献 19被引用数 9

ひとこと要約

本稿では、一般の非重複する論理的四角形パッチで記述される幾何形状に対する音響散乱を扱う、チェビシェフコロケーションおよびフェージェール求積を用いた高次精度の長方形極座標積分ソルバを提案する。表面依存の長方形極座標変換と事前計算済みの特異積分、およびチェビシェフ展開を組み合わせることで、1波長あたり数ポイントの最小限のメッシュでスペクトル精度を達成し、開口部および閉じた表面、複雑なCAD生成幾何形状を含む、堅牢で高速かつ高精度な解法を実現する。

ABSTRACT

This paper introduces a high-order-accurate strategy for integration of singular kernels and edge-singular integral densities that appear in the context of boundary integral equation formulations of the problem of acoustic scattering. In particular, the proposed method is designed for use in conjunction with geometry descriptions given by a set of arbitrary non-overlapping logically-quadrilateral patches---which makes the algorithm particularly well suited for treatment of CAD-generated geometries. Fejér's first quadrature rule is incorporated in the algorithm, to provide a spectrally accurate method for evaluation of contributions from far integration regions, while highly-accurate precomputations of singular and near-singular integrals over certain "surface patches" together with two-dimensional Chebyshev transforms and suitable surface-varying "rectangular-polar" changes of variables, are used to obtain the contributions for singular and near-singular interactions. The overall integration method is then used in conjunction with the linear-algebra solver GMRES to produce solutions for sound-soft open- and closed-surface scattering obstacles, including an application to an aircraft described by means of a CAD representation. The approach is robust, fast, and highly accurate: use of a few points per wavelength suffices for the algorithm to produce far-field accuracies of a fraction of a percent, and slight increases in the discretization densities give rise to significant accuracy improvements.

研究の動機と目的

• 境界積分方程式における音響散乱の設定で生じる特異核および近似特異核に対する高次精度積分法の開発。
• 非重複する論理的四角形パッチで記述される複雑な幾何形状に対する、効率的かつ高精度な散乱問題の解法を可能にすること。
• 最小限の離散化ポイント数（1波長あたり）でスペクトル収束と高精度を達成し、計算コストを低減すること。
• GMRESに基づくソルバに統合することで、開口部および閉じた表面における散乱問題の堅牢でスケーラブルな解法を実現すること。
• 実世界の工学的幾何形状、特に航空機やグライダーのCADモデルへの適用可能性を示すこと。

提案手法

• 本手法はハイブリッド・ニストルム・コロケーション戦略を採用：遠方領域の滑らかな被積分関数に対しては、フェージェールの第一求積則を用い、スペクトル精度を達成する。
• 特異および近似特異的相互作用に対しては、2次元チェビシェフ変換を用いて、核統合密度の高精度事前計算を実施する。
• 核特異点および幾何的端縁に近づくように、表面依存の長方形極座標変換を適用し、解像度と精度を向上させる。
• この座標変換により特異点で微分が消えるようにし、特別な求積則を必要とせずに、点の集約による高次収束を可能にする。
• 密度のチェビシェフ多項式展開により、微分の効率的計算および幾何的量の算出が可能となり、電磁気的境界積分方程式でも有用である。
• フェージェール求積とチェビシェフ変換の間で共通の格子点を活用することで、離散変換の計算が簡素化され、オーバーヘッドが低減される。

実験結果

リサーチクエスチョン

• RQ1 一般の非重複パッチベース幾何形状における境界積分方程式の特異核に対して、高次精度積分法を開発可能か？
• RQ2 エッジおよび幾何的特異点を扱いながら、1波長あたり最小限のポイント数でスペクトル精度を達成するにはどうすればよいか？
• RQ3 幾何的再パラメータライゼーションやメッシュ生成を必要とせずに、複雑なCAD生成幾何形状にどの程度適用可能か？
• RQ4 遠方場の精度および解の誤差の観点から、前方演算子および完全な散乱ソルバの収束特性はいかなるものか？
• RQ5 本手法を電磁気的散乱および音響散乱以外の境界積分方程式設定へ拡張可能か？

主な発見

• 1波長あたり数ポイントで遠方場の精度が1パーセント未満に達し、高い効率性と堅牢性を示した。
• k=1の2×2×2の立方体に対して、エッジ固有の座標変換を適用した場合、遠方場誤差が桁違いに低減され、p=2で顕著な改善が得られた。
• k=10π（λ=0.2）の半径1の円板に対して、本手法はポアソンスポット現象を正確に再現し、物理的妥当性および微細構造の解像度を確認した。
• 開口部および閉じた表面の両方において、滑らかで特異な被積分関数に対して、理論的期待と整合的な高次収束率を示した。
• 87パッチからなるCAD生成グライダーに対して、本手法は全場および強度分布を正確に計算し、実世界の工学的幾何形状への適用可能性を示した。
• チェビシェフに基づく事前計算と長方形極座標変換の併用により、過剰なメッシュ細分化を伴わず、高精度を達成できることが実証された。