<< 4 次精度中心差分による音波計算(2) :
xz 方向の音波の伝播(音波減衰あり) >>
[deepconv/sample]

概要

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

計算設定

プログラムソース arare2
重力 なし
温度 温度一様(298.15 K)
音速 347 m/s
粘性 なし
定圧比熱 1004 J K^-1 kg^-1
積分時間 200 秒
計算領域 水平(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
0 -- 50 秒までのエクスナー関数の時間変化
arare_dim_ss-R50_TL3_TS03_sw8_a1.0d-6_center4.nc
arare_dim_fs-R50_TL3_TS03_sw8_a1.0d-6_center4.nc
arare_dim_sf-R50_TL3_TS03_sw8_a1.0d-6_center4.nc

Exner_sw8_a1.0d-6_center4-t0
t = 0 秒
Exner_sw8_a1.0d-6_center4-t10
t = 10 秒
Exner_sw8_a1.0d-6_center4-t20
t = 20 秒

Exner_sw8_a1.0d-6_center4-t30
t = 30 秒
Exner_sw8_a1.0d-6_center4-t40
t = 40 秒
Exner_sw8_a1.0d-6_center4-t50
t = 50 秒
z = 5000 m での x-t 図
描画の技術的問題から t の範囲は 0--50 となっている.

Exner_sw8_a1.0d-6_center4-xt
VelX_sw8_a1.0d-6_center4-xt
VelZ_sw8_a1.0d-6_center4-xt
x = 5000 m での z-t 図
描画の技術的問題から t の範囲は 0--50 となっている.

Exner_sw8_a1.0d-6_center4-zt
VelX_sw8_a1.0d-6_center4-zt
VelZ_sw8_a1.0d-6_center4-zt
音波減衰項
5.0d-7 * dx^2 / dt
0 -- 50 秒までのエクスナー関数の時間変化
arare_dim_ss-R50_TL3_TS03_sw8_a5.0d-7_center4.nc
arare_dim_fs-R50_TL3_TS03_sw8_a5.0d-7_center4.nc
arare_dim_sf-R50_TL3_TS03_sw8_a5.0d-7_center4.nc

Exner_sw8_a5.0d-7_center4-t0
t = 0 秒
Exner_sw8_a5.0d-7_center4-t10
t = 10 秒
Exner_sw8_a5.0d-7_center4-t20
t = 20 秒

Exner_sw8_a5.0d-7_center4-t30
t = 30 秒
Exner_sw8_a5.0d-7_center4-t40
t = 40 秒
Exner_sw8_a5.0d-7_center4-t50
t = 50 秒
z = 5000 m での x-t 図
描画の技術的問題から t の範囲は 0--50 となっている.

Exner_sw8_a5.0d-7_center4-xt
VelX_sw8_a5.0d-7_center4-xt
VelZ_sw8_a5.0d-7_center4-xt
x = 5000 m での z-t 図
描画の技術的問題から t の範囲は 0--50 となっている.

Exner_sw8_a5.0d-7_center4-zt
VelX_sw8_a5.0d-7_center4-zt
VelZ_sw8_a5.0d-7_center4-zt
音波減衰項
2.5d-7 * dx^2 / dt
0 -- 50 秒までのエクスナー関数の時間変化
arare_dim_ss-R50_TL3_TS03_sw8_a2.5d-7_center4.nc
arare_dim_fs-R50_TL3_TS03_sw8_a2.5d-7_center4.nc
arare_dim_sf-R50_TL3_TS03_sw8_a2.5d-7_center4.nc

Exner_sw8_a2.5d-7_center4-t0
t = 0 秒
Exner_sw8_a2.5d-7_center4-t10
t = 10 秒
Exner_sw8_a2.5d-7_center4-t20
t = 20 秒

Exner_sw8_a2.5d-7_center4-t30
t = 30 秒
Exner_sw8_a2.5d-7_center4-t40
t = 40 秒
Exner_sw8_a2.5d-7_center4-t50
t = 50 秒
z = 5000 m での x-t 図
描画の技術的問題から t の範囲は 0--50 となっている.

Exner_sw8_a2.5d-7_center4-xt
VelX_sw8_a2.5d-7_center4-xt
VelZ_sw8_a2.5d-7_center4-xt
x = 5000 m での z-t 図
描画の技術的問題から t の範囲は 0--50 となっている.

Exner_sw8_a2.5d-7_center4-zt
VelX_sw8_a2.5d-7_center4-zt
VelZ_sw8_a2.5d-7_center4-zt
     
Last update: Thu Apr 14 16:03:50 JST 2005
Copyright © dcmodel 2005 (ODAKA Masatsugu)