ITPASS 実習レポート７月２７日出題分

数値計算実習課題

氏名　寺岡拓也

問 1-A 風上差分

スクリプト URL

http://itpass.scitec.kobe-u.ac.jp/~tteraoka/report02/scripts/quiz1a.txt

使い方

1. chmod a+x quiz1a.rb
2. ./quiz1a.rb を入力する
3. 移流速度を聞かれるので、入力する

結果アニメーション

● c が0より大きい　

● c が0より小さい

参考•資料

納多さんのサンプルスクリプトを利用しました。

共同作業者

特になし。

工夫したこと, 感想

プログラム実行の時好きな移流速度を入力できるようにした。

問 1-B 安定なパラメータ条件の確認

パラメータ設定

　　・ 格子点数 m = 320
　　・ 系の大きさ L = 2π
　　・ 移流速度 c = 1
　　・ 全ステップ数 nt = 300 大の時　3000　小の時　
　　・ 出力間隔 nt_out = 5　大の時　30　小の時　
　　・ 時間間隔　(ⅰ) Δx / c より少しに大きい場合　Δt = 0.0197
　 　　　 　　　 (ⅱ) 上に同じ　小さい場合　Δt = 0.0195

結果アニメーション

(ⅰ) 少し大きい場合

(ⅱ) 少し小さい場合

参考•資料

前問と同じく、納多さんのスクリプトを利用しました。

共同作業者

特になし。

工夫したこと, 感想

・アニメを見やすくするパラメータの設定がイマイチわかりませんでした。