課題2 (もしかすると二つ出して、選択問題にするかも)

2次元渦無し非圧縮流の定常流で、流体の占める領域 $\Omega$ と、その境界 $\Gamma=\rd\Omega$ での流速の法線成分 $v_n:=\bm{v} \cdot\bm{n}$ が分かっている場合に、速度ポテンシャル $\phi$ 、 流れ関数 $\psi$ を計算して、等ポテンシャル線、流線、速度場を 可視化せよ。 領域$\Omega$ と境界値(流速の法線成分) $v_n$ は、 自分で興味のあるもの、自分の都合の良いものを選んで良い (後の注意を読んでおくこと)。

$\phi$ は、ポテンシャル問題

$$\Laplacian\phi=0 \Omega$$ (1)

$$\frac{\rd\phi}{\rd\bm{n}}=v_n \Gamma$$ (2)

の解である。

ポテンシャル問題 (1), (2) を解いて、等ポテンシャル線と速度場 $\bm{v}$ を求めるサンプル・プログラム potential2d-v0.edp を公開してある。

ターミナルで次のようにして入手する

curl -O http://nalab.mind.meiji.ac.jp/~mk/complex2/potential2d-v0.edp

大筋は、 $\Omega$ と $v_n$ を自分が決めたものにするようにプログラムを書き換えれば良い。 (弱形式は変更する必要がない。)

流線の書き方には色々なやり方がある (一つくらいノーヒントの問を入れておくことにする)。 選んだ問題によっては、分かりやすい図が描けるように調整が必要な場合もある。