🌐 JP | 🇬🇧 EN | Last sync: 2025-11-16
基礎数理道場 > 偏微分方程式と境界値問題 > 第5章
🎯 学習目標
- 有限差分法（FDM）の基礎と様々なスキームを理解する
- 有限要素法（FEM）の理論と実装を習得する
- 時間積分スキーム（陽解法・陰解法）の特性を学ぶ
- 安定性解析と収束性の理論的基礎を理解する
- メッシュ生成と要素の選択を学ぶ
- 2次元・3次元問題への拡張を理解する
- スパース行列の効率的な扱い方を習得する
- プロセスシミュレーション（熱処理、反応拡散）への応用を実装する
📖 数値解法の基礎
数値解法の分類
有限差分法（Finite Difference Method, FDM） :
- 微分を差分近似で置き換える
- 構造格子で実装が容易
- 複雑形状への適用が困難
有限要素法（Finite Element Method, FEM） :
- 変分原理に基づく弱形式を利用
- 非構造格子で複雑形状に対応
- 要素内での補間が高精度
有限体積法（Finite Volume Method, FVM） :
- 保存則を積分形式で扱う
- 流体力学で広く使用
- 質量・エネルギー保存が厳密
安定性と収束性
安定性（Stability） : 数値誤差が時間発展で発散しない条件
CFL条件（Courant-Friedrichs-Lewy） : 波動方程式の安定性条件
\[ C = c \frac{\Delta t}{\Delta x} \leq C_{\text{max}} \]
収束性（Convergence） : メッシュ幅 \(\Delta x \to 0\) で真の解に近づく性質
一貫性（Consistency） : 差分式が微分方程式に収束する性質
Laxの等価定理 : 一貫性 + 安定性 ⇒ 収束性
まとめ
- 有限差分法 は実装が容易だが、複雑形状への適用が困難。FTCS, BTCS, Crank-Nicolsonの特性を理解することが重要
- 有限要素法 は変分原理に基づき、非構造格子で複雑形状に対応。線形三角形要素が基本
- 安定性と収束性 は数値解法の信頼性を保証する基本概念。CFL条件やLaxの等価定理が重要
- 適応的メッシュ細分化 により、効率的に高精度解が得られる
- 時間依存問題 は半離散化後、陽解法または陰解法で時間積分する
- 非線形問題 はNewton-Raphson法などの反復法で解く
- スパース行列 の効率的な扱いが大規模問題の鍵
- プロセスシミュレーション（焼入れ、熱応力）など、材料科学への実用的応用が豊富
← 第4章: 変分法と最適化（準備中） シリーズTOPへ
免責事項
- 本コンテンツは教育・研究・情報提供のみを目的としており、専門的な助言(法律・会計・技術的保証など)を提供するものではありません。
- 本コンテンツおよび付随するCode examplesは「現状有姿(AS IS)」で提供され、明示または黙示を問わず、商品性、特定目的適合性、権利非侵害、正確性・完全性、動作・安全性等いかなる保証もしません。
- 外部リンク、第三者が提供するデータ・ツール・ライブラリ等の内容・可用性・安全性について、作成者および東北大学は一切の責任を負いません。
- 本コンテンツの利用・実行・解釈により直接的・間接的・付随的・特別・結果的・懲罰的損害が生じた場合でも、適用法で許容される最大限の範囲で、作成者および東北大学は責任を負いません。
- 本コンテンツの内容は、予告なく変更・更新・提供停止されることがあります。
- 本コンテンツの著作権・ライセンスは明記された条件(例: CC BY 4.0)に従います。当該ライセンスは通常、無保証条項を含みます。