数値モデリング演習(2007年度後期)
|
本演習の目的
地球流体的題材をもとに, 実際に方程式を数値的に解くことで数値モデリングとは何かを体得する. 気候モデリング特論とリンクしており, 講義の 第5回目までが終了してから演習を行う.
担当
教員: 渡部(第1回),中村(第2, 3回)
TA: 古関, 山本
第1回 11/19(月) 13:00-15:00
内容
以下の1次元移流(減衰)方程式
を, 初期条件
のもとで数値的に解くことを考える. 差分の設定は以下のようにする.
N | J | Δt | Δx | tの範囲 | xの範囲 |
---|---|---|---|---|---|
200 | 50 | 1.0 | 2.0 | 0,200 | 0,100 |
適当なスキームで数値解を求め, 真の解と比較して用いた差分スキームの特性を 理解する.
α | γ | c | - |
---|---|---|---|
減衰係数 | 初期値のx幅 | 移流速度 | コメント |
0. | 20 | 1.0 | デフォルト |
0. | 5 | 1.0 | 数値分散がよく見える |
0. | 20 | 2.005 | CFLが満たされない |
0.01 | 20 | 1.0 | 移流減衰解 |
第3回 未定 (月曜 13:00-15:00)