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

桂田 [4] の該当箇所
https://m-katsurada.sakura.ne.jp/labo/text/welcome-to-freefem/node12.html

日本応用数理学会の 「ソフトウェアセミナー：FreeFem++による有限要素プログラミング - 中級編- (2016/2/11, 12)」 で鈴木厚先生が解説された資料を拝借している。

このセミナーの資料は (セミナー当時とは違って) 以下から入手できる。

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

応用数値解析の過去の講義で、このプログラムを解説している。そのときのスライド資料

https://m-katsurada.sakura.ne.jp/ana2022/ANA12_1219_handout.pdf#page=6

(方程式等はこの資料を見て下さい。)

ターミナルでこうすれば実行できる

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/NSstream.mp4 流線
• https://m-katsurada.sakura.ne.jp/ana2026/NSvp.mp4 速度場と圧力場