数値シミュレーション技法基礎セミナー第3回熱伝導と行列解法.pdfVIP

数値シミュレーション技法基礎セミナー 2007/09/07 数値シミュレーション技法基礎セミナー 第３回　熱伝導と行列解法 担当：中久喜伴益 (理学研究科?地球惑星システム学専攻) 2007 年９月７日 数値シミュレーション技法基礎セミナー 2007/09/07 目次 １．偏微分方程式の数値解法とシミュレーション 1.1 偏微分方程式の型 1.2 物質の輸送を表す微分方程式 1.3 偏微分方程式の数値解法と離散化の方法 1.4 定常解を解く問題：境界値問題 ２．離散化方程式の求め方～熱伝導 (拡散) 方程式を例に 2.1 拡散方程式 2.2 有限差分法 2.3 有限体積法 ３．熱伝導 (拡散) 方程式の差分方程式 3.1 一次元熱伝導方程式 3.2 二次元熱伝導方程式 3.3 時間積分の方法 ４．プレートの冷却のシミュレーション 4.1 プレートの温度構造 4.2 変数変換と数値積分による解法 4.3 差分法による解法 4.4 差分法における境界条件の与え方 4.5 プレート冷却シミュレーションのアルゴリズム ５．波動の数値解析 5.1 高階の時間微分の扱い 5.2 波動の数値シミュレーションのアルゴリズム ６．連立一次方程式の解法 6.1 直接法と反復法 6.2 ガウスの消去法 6.3 対称行列の場合 6.4 特殊な係数行列を持つ連立一次方程式の解法 6.5 三重対角行列のガウスの消去法 6.6 帯対称行列に対する修正コレスキー法 6.7 反復法による連立一次方程式の解法 6.8 緩和法によるポワソン方程式の解法 ７．参考文献など 7.1 参考文献 7.2 参考資料 数値シミュレーション技法基礎セミナー 2007/09/07 １．偏微分方程式の数値解法とシミュレーション 1.1 偏微分方程式の型 楕円型：ポワソン方程式など ?2φ = ρ 境界値問題。離散化して得られる連立一次方程式を解く。 放物型：熱伝導方程式 (拡散方程式)、移流 (対流)?拡散方程式（非圧縮粘性流れ） Dφ = κ?2φ Dt 初期値問題。現在のタイムステップの値を使って、繰り返し直後の値を求める。陽解法 では連立一次方程式を解くことなしに直接値を決められる。他に比べて解くのが楽。 双曲型：移流方程式（超音速流れ）、波動方程式 ?φ +u ? ?φ = 0 ?t 初期値問題。数値的に不安定になりやすく解法に工夫が必要。 1.2 物質の輸送を表す微分方程式 生成?消滅 dφ = aφ dt 拡散 ?φ = κ?2φ ?t 移流?拡散 Dφ 2 ?φ = κ? φ , +u ? ?φ = κ?2φ Dt ?t 1.3 偏微分方程式の数値解法と離散化の方法 空間の偏微分を離散化する方法には主に以下３つの方法がある。 ?有限差分法 (FDM) 物体を点の集まりと考え、テイラー展開に基づいて、微分を差分に置き換え、差分方程 式を作る。 ?有限体積法 (FVM) 物体を小さな箱（有限体積）に分け、体積での保存則を満たすように積分で表した保存 則の式を離散化して差分方程式に置き換える。 ?有限要素法 (FEM) 物体を小さな要素に分け、要素内部で求める関数を２次式などの次数の低い多項式で近 似して表す。多項式の係数を微分方程式に適合するように最小二乗法などで決める。 他にも、粒子法やスペクトル法がある。時間方向の離散化にはどの方法で差分法を用い る。 数値シミュレーション技法基礎セミナー 2007/09/07 1.4 定常解を解く問題：境界値問題 定常状態の解を求める場合には時間微分をゼロと置いた式、 ? 2T 0 = κ ?x2 を離散化する。この場合は、解が境界条件だけで決まる。離散化