[論文レビュー] Scalar boundary value problems on junctions of thin rods and plates. I. Asymptotic analysis and error estimates

Renata Bunoiu, Giuseppe Cardone+1 more|arXiv (Cornell University)|Mar 24, 2014

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

ひとこと要約

本稿では、薄板と薄い棒の接合部における、異なる材料特性を有するスカラー境界値問題について、漸近解析と誤差推定を展開する。問題は小さな厚さパラメータ $ h $ と伝導度対比 $ \gamma_j(h) = \gamma_j h^{-\alpha} $ を用いてモデル化される。ノイマンおよびディリクレ境界条件の下で、解の収束速度を非一様重み付きソボレフノルムで導出し、$ \alpha = 1 $ の場合に $ O(|\ln h|^{-1}) $、$ \alpha = 0 $ の場合に $ O(h^{3/2}) $ の収束率を示す。境界層と次元削減問題を含む明示的な漸近展開が得られる。

ABSTRACT

We derive asymptotic formulas for the solutions of the mixed boundary value problem for the Poisson equation on the union of a thin cylindrical plate and several thin cylindrical rods. One of the ends of each rod is set into a hole in the plate and the other one is supplied with the Dirichlet condition. The Neumann conditions are imposed on the whole remaining part of the boundary. Elements of the junction are assumed to have contrasting properties so that the small parameter, i.e. the relative thickness, appears in the differential equation, too, while the asymptotic structures crucially depend on the contrastness ratio. Asymptotic error estimates are derived in anisotropic weighted Sobolev norms.

研究の動機と目的

• 薄板と薄い棒の接合部における、異なる材料特性を有するポアソン方程式の解の漸近的挙動を分析すること。
• 解の漸近的近似に関する、非一様重み付きソボレフノルムにおける厳密な誤差推定を導出すること。
• 対比パラメータ $ \alpha $（0 または 1）が漸近的構造および収束速度に与える影響を特定すること。
• 特に $ \alpha = 1 $ の場合に、接合部および棒の端部に境界層が存在するかを同定し、その形を特定すること。
• 漸近展開と伝達条件を用いて、3次元から1次元および2次元モデルへの次元削減を正当化すること。

提案手法

• 薄板 $ \Omega_\bullet(h) $ と薄い棒 $ \Omega_j(h) $ からなる接合部におけるポアソン方程式の混合境界値問題を定式化し、主に境界でノイマン条件を、棒の端部でディリクレ条件を課す。
• 相対的厚さを表す小さなパラメータ $ h $ と、異なる材料の伝導度をモデル化する対比パラメータ $ \gamma_j(h) = \gamma_j h^{-\alpha} $ を導入する。
• 一致する漸近展開を用いて解 $ u(h,x) $ の漸近展開を構築し、$ \alpha = 1 $ の場合に境界層補正項 $ A_j(\ln h) \widetilde{\mathbf{W}}^j $ を含める。
• 極限問題として、板領域 $ \omega_0 $ における2次元ディリクレ問題と、各棒 $ (0,l_j) $ における1次元問題を導出し、それぞれの解 $ U_0^\# $ と $ U_j^\# $ を得る。
• 重み付きソボレフノルムとコンドラチェフ理論を用いて特異性を制御し、$ L^2 $ および $ H^1 $-型ノルムにおける誤差推定を導出する。
• エネルギー推定と非一様不等式を用いて、3次元解が次元削減された2次元および1次元モデルに収束することを正当化する。

実験結果

リサーチクエスチョン

• RQ1 ポアソン方程式の解の漸近的挙動は、対比パラメータ $ \alpha $ にどのように依存するか？
• RQ2 特に $ \alpha = 1 $ の場合に、接合部および棒の端部付近での解の構造はどのようなものか？境界層はどのように形成されるか？
• RQ3 漸近的近似の収束速度 $ h $ における評価は何か？また、$ \alpha = 0 $ と $ \alpha = 1 $ の場合でどのように異なるか？
• RQ4 非一様重み付きソボレフノルムにおける誤差推定は、幾何的特異性と材料対比が解に与える影響をどのように反映するか？
• RQ5 3次元問題は効果的に2次元と1次元の結合モデルに還元可能か？この還元の精度はいかほどか？

主な発見

• $ \alpha = 1 $ の場合、解は $ L^2 $ および $ H^1 $ で、$ O(|\ln h|^{-1}) $ の収束速度で次元削減された2次元-1次元モデルに強く収束し、誤差推定は $ \| \nabla_x (u_0 - U_0^0 - \sum_j \chi_{0j} A_j \widetilde{\mathbf{W}}^j) \|_{L^2(\Omega_\bullet)} + h^{-1} \sum_j \| \nabla_x (u_j - U_j^0) \|_{L^2(\Omega_j)} \leq c h \mathcal{N} $ となる。
• $ \alpha = 0 $ の場合、収束速度は $ O(h^{3/2}) $ であり、誤差推定は $ \| \nabla_x (u_0 - U_0^\#) \|_{L^2(\Omega_\bullet)} + \sum_j \| \nabla_x (u_j - U_j^\#) \|_{L^2(\Omega_j)} \leq c h^{3/2} \mathcal{N} $ となる。
• 棒の端部におけるディリクレ条件は、グリーン関数に対数特異性を引き起こし、板および棒における $ H^1 $-ノルムで $ O(|\ln h|^{-1}) $ の収束をもたらす。
• 境界層補正項 $ A_j(\ln h) \widetilde{\mathbf{W}}^j $ は $ \alpha = 1 $ の場合に不可欠であり、接合部および棒の端部付近の特異的挙動を捉えている。
• 極限問題は、$ \omega_0 $ における2次元ディリクレ問題と、各棒における1次元混合問題であり、解 $ U_0^\# \in H^2(\omega_0) $ および $ U_j^\# \in H^2(0,l_j) $ を得る。
• 板 $ \omega_0 \times (0,1) $ における $ H^1 $-ノルムでの収束は $ O(|\ln h|^{-1/2}) $ であり、勾配における対数特異性の影響を反映している。