[論文レビュー] A Multiscale Method for Porous Microstructures

Donald L. Brown, Daniel Peterseim|arXiv (Cornell University)|Nov 7, 2014

Advanced Mathematical Modeling in Engineering参考文献 16被引用数 11

ひとこと要約

本稿では、穴あき領域におけるPoincaré定数の厳密な解析を通じて、最適収束率を確立することで、切断補正関数を用いた局所直交分解（LOD）を用いた多スケール有限要素法を、多孔質微細構造におけるノイマン問題の解法に適用する。フィラメント構造では微細構造依存の境界を示し、孤立粒子では一様な境界を示す。

ABSTRACT

In this paper we develop a multiscale method to solve problems in complicated porous microstructures with Neumann boundary conditions. By using a coarse-grid quasi-interpolation operator to define a fine detail space and local orthogonal decomposition, we construct multiscale corrections to coarse-grid basis functions with microstructure. By truncating the corrector functions we are able to make a computationally efficient scheme. Error results and analysis are presented. A key component of this analysis is the investigation of the Poincaré constants in perforated domains as they may contain micro-structural information. Using a constructive method originally developed for weighted Poincaré inequalities, we are able to obtain estimates on Poincaré constants with respect to scale and separation length of the pores. Finally, two numerical examples are presented to verify our estimates.

研究の動機と目的

• 細部の特徴が顕著に現れる複雑で高周期的な多孔質微細構造におけるノイマン問題を解く課題に対処すること。
• グローバルに細かいスケールで解くことなく、微細構造の影響を捉えられる計算効率の高い多スケールスキームを開発すること。
• 穴あき領域におけるPoincaré定数が、孔の大きさや分離長などの微細構造パラメータにどのように依存するかを分析すること。
• 穴あき領域におけるグローバルに支持されたおよび局所化された多スケール基底関数の両方の誤差推定を確立すること。
• 重み付きPoincaré不等式の技術を用いて、多孔質媒体におけるPoincaré定数を推定する構成的手段を提供すること。

提案手法

• 粗い穴あきグリッドを構築し、局所L^2射影に基づく準補間作用素を用いて細かいスケール空間を定義する。
• 局所直交分解（LOD）を適用して関数空間を粗いスケール成分と細かいスケール成分に分割し、基底関数の多スケール補正を可能にする。
• 計算効率を確保しながらも精度を維持するために、補正関数を局所的なパッチに切り詰める。
• ガラーキン直交性と減衰推定を用いて、補正関数の切り詰めによって生じる誤差を評価する。
• 粗いメッシュサイズHと切り詰めパラメータmに依存する誤差境界を導出し、収束率O(H)およびmに関して指数関数的減衰を示す。
• 重み付きPoincaré不等式の手法を用いて、穴あき領域におけるPoincaré定数を推定する構成的手順を開発する。

実験結果

リサーチクエスチョン

• RQ1 多スケール法は、振動係数ではなく幾何的特徴（例：孔）を有する微細構造に対処するためにどのように適合可能か？
• RQ2 穴あき領域におけるPoincaré定数は、孔の大きさや分離長などの微細構造パラメータにどのように依存するか？
• RQ3 局所直交分解フレームワークは、厳密な誤差制御が可能な穴あき領域におけるノイマン問題へどのように拡張可能か？
• RQ4 微細構造は補正関数の減衰率および局所化スキームの収束にどのように影響を与えるか？
• RQ5 Poincaré定数が微細構造に対して一様に有界となる条件は何か、また幾何的特徴に敏感になる条件は何か？

主な発見

• 粗いグリッド近似における収束率は最適であり、O(H)であり、切り詰めパラメータmに関して指数関数的減衰を示す。
• 穴あき領域におけるPoincaré定数は微細構造の幾何に依存する：非常に細いフィラメント構造では劣化するが、孤立粒子では一様に有界のまま保たれる。
• Poincaré定数を推定する構成的手法により、特に網状フィラメント微細構造において、分離長および孔の大きさに比例する明示的境界が得られる。
• 局所化された多スケールスキームは、グローバル版と同等の収束率を維持し、誤差は切り詰めパラメータmおよびメッシュサイズHによって制御される。
• 数値例は、H、m、および微細構造パラメータに関して理論的収束率を確認している。
• 本手法で用いられる準補間作用素は、ロバストで効果的であり、さまざまな微細構造配置において性能が検証されている。