<< 4 次精度中心差分による音波計算(1) :
x 方向の音波の伝播(音波減衰あり) >>

[deepconv/sample]

概要

x4 次精度中心差分演算モジュールを用いてx 平面内を伝播する音波の計算を行う. クーラン数を 0.5 程度の値で固定し, 音波減衰項の係数を適当に変えて行う. なお, 微分演算モジュールが異なることを除けば計算設定は音波計算(6)と同じである.

計算設定

プログラムソース arare2
重力 なし
温度 温度一様(298.15 K)
音速 347 m/s
粘性 なし
定圧比熱 1004 J K^-1 kg^-1
積分時間 40 秒
計算領域 水平(x) : 10000 m
鉛直(z) : 10000 m
クーラン数 0.52
格子点間隔 水平(x) : 200 m
鉛直(z) : 200 m
時間間隔 0.3 秒
境界条件 水平(x) : 周期境界条件
鉛直(z) : 固定境界条件
数値解法 水平(x) : 陽解法
鉛直(z) : 陰解法



[HOME]

 

音波減衰項
1.0d-6 * dx^2 / dt
arare_dim_ss-R50_TL3_TS03_sw6_a1.0d-6_center4.nc
arare_dim_fs-R50_TL3_TS03_sw6_a1.0d-6_center4.nc
arare_dim_sf-R50_TL3_TS03_sw6_a1.0d-6_center4.nc



Exner_sw6_a1.0d-6_center4
エクスナー関数の時間変化


VelX_sw6_a1.0d-6_center4
水平速度の時間変化


音波減衰項
2.5d-7 * dx^2 / dt
arare_dim_ss-R50_TL3_TS03_sw6_a5.0d-7_center4.nc
arare_dim_fs-R50_TL3_TS03_sw6_a5.0d-7_center4.nc
arare_dim_sf-R50_TL3_TS03_sw6_a5.0d-7_center4.nc



Exner_sw6_a2.5d-7_center4
エクスナー関数の時間変化


VelX_sw6_a2.5d-7_center4
水平速度の時間変化


音波減衰項
5.0d-7 * dx^2 / dt
arare_dim_ss-R50_TL3_TS03_sw6_a2.5d-7_center4.nc
arare_dim_fs-R50_TL3_TS03_sw6_a2.5d-7_center4.nc
arare_dim_sf-R50_TL3_TS03_sw6_a2.5d-7_center4.nc



Exner_sw6_a5.0d-7_center4
エクスナー関数の時間変化


VelX_sw6_a5.0d-7_center4
水平速度の時間変化

   
Last update: Thu Apr 14 16:02:08 JST 2005
Copyright © dcmodel 2005 (ODAKA Masatsugu)