[論文レビュー] High-Order Quadrature on Multi-Component Domains Implicitly Defined by Multivariate Polynomials

Robert I. Saye|arXiv (Cornell University)|May 19, 2021

Advanced Numerical Analysis Techniques参考文献 81被引用数 28

ひとこと要約

本稿では、多次元多項式によって暗黙的に定義された複雑なマルチコンポonent領域における統合のための高次数の数値積分法を提示する。幾何学的形状を多価関数として再定式化し、1次元のガウス＝ルジャンドル則または双曲正接・双曲正 secant（tanh-sinh）則を適用することで、h-細分において最大22次精度を達成し、q-細分においてはほぼ指数関数的収束を実現する。特異点やトポロジーの複雑さに対しても自動的に処理可能である。

ABSTRACT

A high-order quadrature algorithm is presented for computing integrals over curved surfaces and volumes whose geometry is implicitly defined by the level sets of (one or more) multivariate polynomials. The algorithm recasts the implicitly defined geometry as the graph of an implicitly defined, multi-valued height function, and applies a dimension reduction approach needing only one-dimensional quadrature. In particular, we explore the use of Gauss-Legendre and tanh-sinh methods and demonstrate that the quadrature algorithm inherits their high-order convergence rates. Under the action of $h$-refinement with $q$ fixed, the quadrature schemes yield an order of accuracy of $2q$, where $q$ is the one-dimensional node count; numerical experiments demonstrate up to 22nd order. Under the action of $q$-refinement with the geometry fixed, the convergence is approximately exponential, i.e., doubling $q$ approximately doubles the number of accurate digits of the computed integral. Complex geometry is automatically handled by the algorithm, including, e.g., multi-component domains, tunnels, and junctions arising from multiple polynomial level sets, as well as self-intersections, cusps, and other kinds of singularities. A variety of numerical experiments demonstrates the quadrature algorithm on two- and three-dimensional problems, including: randomly generated geometry involving multiple high-curvature pieces; challenging examples involving high degree singularities such as cusps; adaptation to simplex constraint cells in addition to hyperrectangular constraint cells; and boolean operations to compute integrals on overlapping domains.

研究の動機と目的

• 多次元多項式のゼロレベル集合によって定義される曲面および体積領域における積分のための高次数積分スキームの開発。
• メッシュ生成やヒューリスティックな分割なしに、マルチコンポonent領域、トンネル、接合部、尖点、自己交差などの複雑な幾何形状を処理すること。
• 高次数収束を示す1次元積分則を活用することで、高精度かつ高効率な積分を実現すること。
• 非滑らか関数が界面で不連続であっても、暗黙的定義領域の個々の連結成分を自動的に検出・統合し、精度を保持すること。
• 2次元および3次元の両方で動作する次元に依存しないフレームワークを提供すること。幾何的前処理は最小限に抑える。

提案手法

• 暗黙的定義領域を多価関数として再定式化し、座標方向に沿った1次元積分に統合問題を還元する。
• ガウス＝ルジャンドル則およびtanh-sinh則を含む高次数1次元積分則を、高さ関数に適用して領域上の積分を計算する。
• 各積分セルに対して、多項式系の実数解を区間(0,1)内で特定するため、de Casteljauの分割、ニュートン法にバイセクションの安全措置を組み合わせたハイブリッド根の探索戦略を用いる。
• 積分重みを定数関数を正確に統合できるように正規化し、tanh-sinh則の近似指数的収束特性を保持する。
• 本手法は次元に依存せず、体積積分と表面積分の両方に均等に適用可能であり、それぞれd次元でO(q^d)およびO(q^{d-1})のノード数を有する。
• トポロジーの複雑さ（複数のコンポonent、交差、特異点など）が、多項式レベル集合から自然に生じ、追加のユーザー介入なしに処理可能である。

実験結果

リサーチクエスチョン

• RQ1 1次元積分則のみを用いて、複雑なトポロジーを持つ暗黙的定義領域に対し高次数積分を達成できるか？
• RQ2 このような領域において、h-細分（メッシュサイズ）およびq-細分（積分則の次数）に対する収束速度はどのように変化するか？
• RQ3 明示的なメッシュ生成や分割なしに、複数の連結成分、接合部、特異点（尖点など）を自動で検出し、統合できるか？
• RQ4 ガウス＝ルジャンドル則およびtanh-sinh則を用いることで、曲面界面や幾何的特異点が存在する状況でも高次数精度がどの程度維持されるか？
• RQ5 ハイブリッド根の探索戦略（de Casteljauの分割、ニュートン法、固有値ソルバーの組み合わせ）は、高次多項式系の界面交差をどのようにして高精度かつ安定して計算可能にするか？

主な発見

• 固定qにおけるh-細分において、数値実験により最大22次収束を示し、近似的に2q次の精度を達成する。
• 固定幾何形状におけるq-細分において、収束はほぼ指数的であり、qを2倍にすると計算された積分値の有効数字が概ね2倍に増加する。
• 特異点や自己交差を含む高カーブチャーラの特徴を持つ領域に対しても、特別な処理やメッシュ細分化なしに積分を正しく計算可能である。
• 複数の多項式レベル集合から生じる複数の連結成分、トンネル、接合部を、ユーザー定義の分解なしに自動で検出し、統合可能である。
• de Casteljauの分割、ニュートン法、固有値ソルバーを組み合わせたハイブリッド根の探索戦略は、95％以上の根が10回未満のニュートン反復でフルマシン精度に達し、高い精度と安定性を実現した。
• 正規化されたtanh-sinh積分則を用いることで、定数関数の正確な統合が保証され、滑らかな被積分関数に対しても近似指数的収束特性が維持される。