課題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$ を自分が決めたものにするようにプログラムを書き換えれば良い。 (弱形式は変更する必要がない。)

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