E. 一様流の中に置かれた円柱のまわりの流れ

(しばらく工事中)

「ソフトウェアセミナー：FreeFem++による有限要素プログラミング -- 中級編-- (2016/2/11, 12)」 で公開された NS-cylinder.edp

-- リンク切れになっている。

Internet Archive の WayBack Machine で読むことはできるが、 それによると、資料そのものは、もともとパリ第6大学 Laboratoire Jacques-Louis Lions のサーバーに置いてあった。 現在はそのサーバーも無くなったようだが、代替サーバーがあるようで

https://www.ljll.fr/~suzukia/FreeFempp-tutorial-JSIAM2016/

に置いてある。

Atsushi Suzuki 氏による 解説PDF “Finite element programming by FreeFem++ -- intermediate course” のpp. 44-52 に書かれているものである。

(なお、2022年度の応用数値解析特論の授業資料 「応用数値解析特論 第 12 回 ～Navier-Stokes 方程式に対する有限要素法 (3)～ 非定常 Navier-Stokes 方程式」 は、このプログラムの解説をしている。)

curl -O https://www.ljll.fr/~suzukia/FreeFempp-tutorial-JSIAM2016/EDP/NS-cylinder.edp
FreeFem++ NS-cylinder.edp

時間発展する系なので、動画を見るのが分かりやすいであろう。

https://m-katsurada.sakura.ne.jp/ana2026/NS-cylinder.mp4

(追跡する時間は、オリジナルのプログラムより長くしてあります。)