松元亮治（千叶大理）-AstrophysicsLaboratory.PDFVIP

差分法の基礎 松元亮治（千葉大理） 本章の目標 線形及び非線形の１変数波動方程式を例 にして差分法の基礎を学ぶ 差分法の数値的安定性を学ぶ 解の単調性を維持する解法として一次精 度風上差分法を導入する 本章の参考書：流体力学の数値計算法 藤井孝蔵 著、東京大学出版会（１９９４） 差分近似 格子間隔： Δx、Δｙ 格子点における物理量： 中心差分 上式から下式を引くと 中心差分 各種差分式 式(1.1)+(1.2) 線形スカラー移流方程式 厳密解 ＦＴＣＳスキーム 時間について前進差分、空間について中心差分 クーラン数 FTCS スキームは陽解法 １次元波動伝播シミュレーション のアルゴリズム １．メッシュ生成：座標値ｘ[j]をセットする ２．初期条件：初期値u[j]をセットする ３．時刻ｔが決められた時刻に達するまで以 下を繰り返す (a)Δｔ後の値を差分式にもとづいて計算 (b)境界の値を境界条件から決める (c) 時刻をΔｔだけ進める FTCSスキームによる計算例 クーラン数=0.25 FTCS スキームは数値的に不安定 フォンノイマンの安定性解析 ＦＴＣＳスキームの差分方程式に次式を代入 uj n のプロフィール プロフィールの変化 ＦＴＣＳスキームの安定性 複素数に拡張 増幅率 常に不安定！ Lax-Friedrich スキーム Stable when クーラン条件 クーラン条件の意味 1タイムステップに波が1メッシュ以上 伝わってはいけないという条件 Lax-Friedrich スキームによる 計算結果 クーラン数=0.25 1次精度風上差分法 風上差分法の安定性 Stable when １次精度風上差分による計算例 クーラン数 =0.25 クーラン数 =0.8 高次精度スキームの例: Lax-Wendroff Method テイラー展開 を代入 中心差分を用いると 2step Lax-Wendroff method Lax-Wendroff 法による計算例 クーラン数 =0.25 クーラン数 = 0.80 Godunovの定理 In 1D scalar advection equation,any scheme higher than 2nd order accuracy with form Cannot maintain the monotonicity of the solution 1st Order Upwind Scheme preserves the monotonicity 人工粘性 数値振動を抑制するために、 以下の拡散項を加える 速度勾配に比例した拡散係数の例 保存形 流体力学、磁気流体力学の方程式は 以下の保存形に書くことができる ここで 保存形式の物理的意味 まで空間積分