はじめての音響数値シミュレーションプログラミングガイドについて

[yomogi]

はじめまして。現在、「はじめての音響数値シミュレーションプログラミングガイド」でＦＤＴＤ法を勉強させていただいております。
ＦＤＴＤ法のプログラムの内、プログラム５－１３（ｐ１１３）と５－２３（ｐ１１９）に関しての質問です。

５－２３は境界条件として、インピーダンス境界を適応していると思うのですが、５－１３は特に境界条件についての記載がありませんでした。
５－１３はどのような境界条件になっているのでしょうか。
プログラムを実行してみたところ減衰することなく波が伝わっているようでした。

よろしくお願いいたします。

[Hoshi Kazuma]

はじめまして，星と申します。

プログラム５－１３（ｐ１１３）は，特に境界条件を設定してませんね。

この場合，粒子速度 (プログラムではUx, Uy) が更新されるかポイントです。

プログラム５－１３では，

Ux[0,  y]（四角の下辺）
Ux[41,y]（四角の上辺）
Uy[x,  0]（四角の左辺）
Uy[x,31]（四角の右辺）

の全てが0のまま更新されません。よって四辺は，全て空気の粒子が動かない辺，すなわち剛壁とみなされます。

プログラム中の
<pre>
for n in range(5000)
</pre>
繰り返し数5000を500から1000まで，徐々に変化させながら音波の進行を見ると，固定端反射していることが確認できます。
是非お試しください。