[論文レビュー] Strong form mesh-free $hp$-adaptive solution of linear elasticity problem

Mitja Jančič, Gregor Kosec|arXiv (Cornell University)|Oct 13, 2022

Advanced Numerical Analysis Techniques被引用数 4

ひとこと要約

本稿では、径数関数生成有限差分法（RBF-FD）を用いた、線形弾性問題を解くためのメッシュフリーhp適応型数値法を提示する。本手法は、イノベーティブな陰・陽（IMEX）誤差指標を用い、反復的に解法を行い、局所的な高次再評価により誤差を推定し、修正されたテキサス3段階戦略を用いてメッシュの細分化を判断し、ノードの再配置と近似次数の調整を実施する。2次元および3次元の接触問題において、非適応的手法に比べ優れた収束性と精度を達成した。

ABSTRACT

We present an algorithm for $hp$-adaptive collocation-based mesh-free numerical analysis of partial differential equations. Our solution procedure follows a well-established iterative solve-estimate-mark-refine paradigm. The solve phase relies on the Radial Basis Function-generated Finite Differences (RBF-FD) using point clouds generated by advancing front node positioning algorithm that supports variable node density. In the estimate phase, we introduce an Implicit-Explicit (IMEX) error indicator, which assumes that the error relates to the difference between the implicitly obtained solution (from the solve phase) and a local explicit re-evaluation of the PDE at hand using a higher order approximation. Based on the IMEX error indicator, the modified Texas Three Step marking strategy is used to mark the computational nodes for $h$-, $p$- or $hp$-(de-)refinement. Finally, in the refine phase, nodes are repositioned and the order of the method is locally redefined using the variable order of the augmenting monomials according to the instructions from the mark phase. The performance of the introduced $hp$-adaptive method is first investigated on a two-dimensional Peak problem and further applied to two- and three-dimensional contact problems. We show that the proposed IMEX error indicator adequately captures the global behaviour of the error in all cases considered and that the proposed $hp$-adaptive solution procedure significantly outperforms the non-adaptive approach. The proposed $hp$-adaptive method stands for another important step towards a fully autonomous numerical method capable of solving complex problems in realistic geometries without the user intervention.

研究の動機と目的

• 実際の幾何形状における複雑な線形弾性問題を解くための、完全に自動化されたメッシュフリーhp適応型数値法の開発。
• 散乱ノード点群と適応的細分化戦略を活用することで、メッシュ生成および細分化におけるユーザーの干渉を排除すること。
• h細分化（ノード再配置）とp細分化（変動近似次数）を統合することにより、解の精度と収束効率を向上させること。
• ピーク問題および解析解のない3次元の摩擦疲労接触問題を含むベンチマーク問題において、提案手法の妥当性を検証すること。
• 線形弾性問題におけるグローバル誤差の挙動を、高勾配領域および特異点付近で的確に捉えることができる、新しいIMEXベース誤差指標の導入と評価。

提案手法

• 本手法は、可変ノード密度を持つ点群を生成する前進フロントノード配置アルゴリズムを用いたRBF-FDに基づき、solve-estimate-mark-refineの反復的パラダイムを採用する。
• 解法フェーズでは、多尺度基底関数と可変ステンシルサイズを用いたRBF-FDにより、強形式での偏微分作用素の近似を実施する。
• 誤差推定フェーズでは、陰的に解かれた解と、局所的に陽的に高次再評価されたPDEとの差分を計算するIMEX誤差指標を導入する。
• マークフェーズでは、IMEX誤差指標に基づき、ノードをh-細分化、p-細分化、またはhp-細分化のいずれに分類するかを決定する、修正されたテキサス3段階戦略を用いる。
• 細分化フェーズでは、前進フロントアルゴリズムを用いてノードを再配置し、マークの結果に従い、局所的に多項式を追加することで近似次数を調整する。
• 本手法は、メッシュトポロジーを必要とせず、統一されたフレームワーク内でh細分化（ノード追加／削除）とp細分化（次数上昇）を両立させることで、hp適応性を実現する。

実験結果

リサーチクエスチョン

• RQ1 IMEXベース誤差指標は、メッシュフリーRBF-FD解における線形弾性問題のグローバル誤差挙動を、効果的に捉えることができるか？
• RQ2 提案されたhp適応フレームワークは、非適応的RBF-FD手法に比べ、顕著に収束速度と精度を向上させるか？
• RQ3 特異点や局所的応力集中を示す接触問題のような問題に対して、本手法はどの程度の性能を示すか？
• RQ4 ユーザー定義のメッシュ生成や手動での細分化なしに、複雑な3次元幾何形状においても高い精度を達成できるか？
• RQ5 可変ノード密度と可変近似次数の組み合わせは、急峻な勾配および境界層を的確に解像するのに有効であるか？

主な発見

• IMEX誤差指標は、2次元ピーク問題および3次元摩擦疲労接触問題を含むすべてのテストケースにおいて、グローバル誤差挙動を的確に捉えた。
• 提案されたhp適応型手法は、非適応的RBF-FDに比べ、特に高応力勾配領域および特異点付近で顕著に高速な収束と高い精度を達成した。
• 摩擦疲労接触問題では、0.2067 mm（ドメイン幅の200分の1）の接触領域を高解像度で解像し、微小スケール特徴のモデル化能力を示した。
• 初期段階の反復で特異点付近で近似次数を自動的に増加させ、その後にh細分化により局所的場の記述を向上させた。これは、適応的細分化に関する既存の文献と整合的であった。
• 可変ノード密度と適応的次数選択の導入により、計算リソースのより効率的な配分が実現され、全体的な計算コストが低減した。
• 角ノードを削除して数値的安定性を確保した場合でも、本手法は頑健で安定した挙動を示し、実世界の問題への実用的妥当性を裏付けた。