[論文レビュー] Interpolants of Lattice Functions for the Analysis of Atomistic/Continuum Multiscale Methods

Christoph Ortner, Alexander V. Shapeev|arXiv (Cornell University)|Apr 17, 2012

Advanced Mathematical Modeling in Engineering参考文献 10被引用数 19

ひとこと要約

本稿では、ブラーヴェ格子上の格子関数に対して一般化された（準）補間関数のクラスを導入し、離散的な原子スケールデータの正確な連続的表現を可能にする。ノード補間関数における2次精度および準補間関数における3次精度といった、原子スケール／連続体マルチスケール法の解析に不可欠な主要な近似特性を確立し、ソボレフ空間およびレベーグ空間における厳密な誤差評価を提示する。

ABSTRACT

We introduce a general class of (quasi-)interpolants of functions defined on a Bravais lattice, and establish several technical results for these interpolants that are crucial ingredients in the analysis of atomistic models and atomistic/continuum multi-scale methods.

研究の動機と目的

• マルチスケールモデリングにおける離散的格子関数の連続的補間関数の一般枠組みを構築すること。
• 数値解析のための連続関数による離散的原子スケールデータの近似において高次精度を保証すること。
• ノード補間関数および準補間関数のソボレフ空間およびレベーグ空間における厳密な近似誤差推定を確立すること。
• 安定で可逆的かつ局所的にサポートされる補間関数を通じて、原子スケール／連続体結合法の解析を支援すること。
• マルチスケール法の数値解析における既存結果を統一的かつ拡張的に扱える関数空間設定および技術的ツールを提供すること。

提案手法

• コンパクトに台を持つ $ W^{1,\infty} $ の基底関数 $ \bar{\zeta} $ を用い、$ \bar{v}(x) = \sum_{\xi \in \mathbb{Z}^d} v(\xi) \bar{\zeta}(x - \xi) $ としてノード補間関数を定義。この $ \bar{\zeta} $ は補間条件および多項式再現条件を満たす。
• 畳み込みにより $ \tilde{v}(x) = (\bar{\zeta} \ast \bar{v})(x) $ として準補間関数を導入し、より高い正則性を保ち、3次精度を達成可能にする。
• 基底関数 $ \bar{\zeta} $ に構造的仮定（Z1–Z4）を課す：正則性、コンパクトな台、アフィン関数の多項式再現、整数点におけるノード補間。
• 結合密度補題および畳み込み恒等式を用いて、滑らかな関数に対する補間関数による近似誤差推定を導出する。
• 逆推定および局所的ノルム同値性を確立し、$ \|\nabla^j(\tilde{J}v - v)\|_{L^p} $ を $ \|\nabla^k v\|_{L^p} $ の項で評価する。ここで $ \tilde{J} $ は準補間関数空間への射影作用素。
• 任意の立方多項式が、ある離散関数 $ w \in \mathscr{U} $ を用いて $ \tilde{w} $ として正確に表現可能であることを証明し、準補間関数空間の近似の一貫性を保証する。

実験結果

リサーチクエスチョン

• RQ1 離散的格子関数を、物理的および数学的性質を保持したまま、一貫的かつ安定的に連続関数に写像する方法は何か？
• RQ2 基底関数 $ \bar{\zeta} $ にどのような条件を課すことで、ノード補間関数では2次精度、準補間関数では3次精度が達成されるか？
• RQ3 補間関数の近似誤差は、ソボレフノルムにおいて一様に有界に抑えられるか？また、その依存関係は、元の関数の正則性にどのように依存するか？
• RQ4 補間関数フレームワークは、原子スケール／連続体マルチスケール法における収束性および一貫性の解析にどのように応用できるか？
• RQ5 補間演算が可逆となる条件は何か？また、その可逆性はマルチスケール結合スキームの安定性にどのように影響するか？

主な発見

• ノード補間関数 $ \bar{v} $ は $ W^{1,\infty}_{\rm loc}(\mathbb{R}^d; \mathbb{R}^m) $ に属し、局所的リプシッツ連続性と安定性を保証する。
• 準補間関数 $ \tilde{v} $ は $ W^{3,\infty}_{\rm loc}(\mathbb{R}^d; \mathbb{R}^m) $ に属し、より高い正則性と滑らかな関数に対する3次精度を保証する。
• 任意の $ v \in W^{k,p}_{\rm loc}(\mathbb{R}^d) $ に対して $ k \leq 4 $ のとき、誤差 $ \|\nabla^j(\tilde{J}v - v)\|_{L^p} $ は $ C_{\tilde{J}} \|\nabla^k v\|_{L^p} $ で有界である。これは最適収束率を示す。
• 任意の立方多項式 $ p $ は、ある $ w \in \mathscr{U} $ を用いて $ \tilde{w} $ として正確に表現可能であり、準補間関数空間の一貫性が保証される。
• 補間関数 $ \tilde{J} $ は準補間関数空間への射影作用素であり、その近似誤差は関数の $ k $ 階ソボレフ半ノルムに依存して制御される。
• 基底関数 $ \bar{\zeta} $ をQ1基底関数（区分的線形）として選ぶと、$ \bar{v} $ はテンソル積線形Bスプライン、$ \tilde{v} $ は立方Bスプラインとして得られ、このフレームワークの一般性および標準的スプライン関数との関連性が裏付けられる。