3章流体および磁気流体力学方程式の風上差分.pdfVIP

61 第3章 流体および磁気流体力学方程式の風 上差分 花輪知幸（名大理） 3.1 はじめに 私たちは初日に移流方程式や Burgers 方程式の数値解法として、風上差分法を習った。風上 差分法は衝撃波の取り扱いに優れていることや、非物理的な数値振動を起こさない (=TVD 条 件を遵守する) ことを学んだ。昨日は、流体力学方程式や磁気流体力学方程式は波動方程式の 集合体—システム方程式—であることを学んだ。この 2 日間の学習成果を総合すると、風上差 分法は流体力学方程式や磁気流体力学方程式の数値解法としても有効であろうと容易に想像で きる。実際、衝撃波を伴う (磁気) 流れを解析する方法として多くの数値計算コードに採用され ている。 「流束は風上で評価せよ」という風上差分法の原理 (概念) は単純であるが、これを実際に流 体力学方程式や磁気流体力学方程式に応用することは簡単ではない。最初に習った線型波動方 程式は 1 成分で波の位相速度が一定であったのに対し、流体力学方程式は連立であり波の速度 も時間や場所によって異なるためである。これらの違いをまとめたのが表 3.1 である。 名称 変数 自由度 線型/非線型 方程式 移流方程式 u 1 線型 3.1 Burgers 方程式 u 1 非線型 3.2, 3.3 Maxwell 方程式 Ey, Ez, Bx, By 4 線型 3.4, 3.5 流体力学方程式 ρ, v, P 3 非線型 3.6 磁気流体力学方程式 ρ, vx, vy, vz, By, Bz 7 非線型 3.7 表 3.1: さまざまな波動方程式の比較。変数の数は独立変数を (x, t) とした 1 次元の場合で数 えた。 移流方程式 ?u ?u + c = 0 (3.1) ?t ?x Burgers 方程式 [波の速度を明示した形式 (上) と流束を明示した形式 (下) ?u ?u + u = 0 (3.2) ?t ?x 62 第 3 章 流体および磁気流体力学方程式の風上差分
 ?u ? u2 + = 0 (3.3) ?t ?x 2 (真空中の)Maxwell 方程式 [ベクトル形式と成分ごとに分解した形式]