ポテンシャル流れシミュレーター 戻る
流体力学・CFD

ポテンシャル流れシミュレーター

円柱まわりの2次元ポテンシャル流れを流線・圧力係数カラーマップで可視化。マグヌス効果(回転する円柱による揚力)も体験。ダランベールのパラドックスを確認。

流れタイプ
パラメータ
一様流速 U∞ (m/s) 1.0
円柱半径 R 1.0
表示設定
解析結果
0.00
揚力係数 CL
0.00
抗力係数 CD
−3.00
Cp (上端)
−3.00
Cp (下端)

理論メモ

流れ関数(円柱 + 循環):
$\Psi = U_\infty r\!\left(1-\frac{R^2}{r^2}\right)\!\sin\theta - \frac{\Gamma}{2\pi}\ln r$

揚力(Kutta-Joukowski):
$L = \rho\, U_\infty\, \Gamma$

圧力係数:
$C_p = 1 - \left(\frac{V}{U_\infty}\right)^2$
Cp低(高速)
Cp高(低速)
キャンバスをクリックすると局所速度とCpを表示

上図: 流線 + Cpカラーマップ  |  下図: 円柱表面のCp分布(θ = 0° → 360°)

ポテンシャル流れシミュレーターとは

🧑‍🎓
このシミュレーターで「循環Γ」を増やすと、円柱の周りの流れがぐるっと回り始めますね。これが何の役に立つんですか?
🎓
ざっくり言うと、これが「揚力」を生み出す仕組みなんだ。例えば、回転する野球のボールがカーブする「マグヌス効果」は、この循環によって生じる圧力差が原因なんだよ。上のスライダーでΓを大きくすると、円柱の上の流速が上がり(青い部分)、圧力が下がるのがわかるね。これで上向きの力が働くんだ。
🧑‍🎓
え、そうなんですか!でも、流れが完全に対称だと、円柱を流れる流体は左右で打ち消し合って、抗力がゼロになるって聞きました。これが「ダランベールのパラドックス」ですよね?
🎓
その通り!理想的なポテンシャル流れでは、圧力分布が前後対称になるから抗力がゼロになるんだ。実際の流体には粘性があるからそうはならないんだけどね。シミュレーターで「循環Γ」を0に戻して、円柱の表面の色(圧力係数$C_p$)を見てごらん。前と後ろでほとんど同じ色(圧力)になってるだろ?これがパラドックスの可視化だよ。
🧑‍🎓
なるほど!じゃあ「流量Q」のスライダーは何に使うんですか?マイナスにすると吸い込まれるみたいですが。
🎓
これは「点源」や「点吸込み」を表現しているんだ。実務では、例えばエンジンの吸気口や換気扇の吸い込み口付近の流れを、このような単純化したモデルでまず理解するのに使うんだ。Qをプラス(源)にすると円柱から四方に流れが出ていく様子が、Cpカラーマップを見ると円柱表面の圧力がどう変わるかがわかる。いろいろ触ってみて、流線のパターンがどう変わるか観察してみよう。

物理モデルと主要な数式

このシミュレーションの基礎は、一様流中の円柱の周りの「流れ関数」です。これに循環(自由渦)や源・吸込みの効果を重ね合わせています。

$$\Psi = U_\infty \left(r - \frac{R^2}{r}\right)\sin\theta - \frac{\Gamma}{2\pi}\ln r + \frac{Q}{2\pi}\theta$$

ここで、$\Psi$は流れ関数[m²/s]、$U_\infty$は一様流速[m/s]、$r, \theta$は円柱中心からの極座標、$R$は円柱半径[m]、$\Gamma$は循環強度[m²/s]、$Q$は流量[m²/s](正で源、負で吸込み)です。流線は$\Psi$=一定の線として描画されます。

流れから生じる最も重要な力である「揚力」は、クッタ・ジュコーフスキーの定理で与えられます。また、表面の圧力分布は圧力係数$C_p$で評価します。

$$L = \rho\, U_\infty\, \Gamma \quad , \quad C_p = 1 - \left(\frac{V}{U_\infty}\right)^2$$

$L$は単位長さあたりの揚力[N/m]、$\rho$は流体密度[kg/m³]です。$C_p$は表面速度$V$を用いて定義され、$C_p$が低い(青)ほど圧力が低く、高い(赤)ほど圧力が高いことを示します。揚力はこの圧力差の積分結果です。

実世界での応用

航空工学(翼理論):薄翼理論の基礎となる考え方です。循環を付加することで翼に働く揚力を説明し、クッタ・ジュコーフスキーの定理は翼設計の根幹をなします。シミュレーターでΓを変えて揚力が変化する様子は、その直接的な体験です。

スポーツボールの流体解析:サッカーのバナナシュートや野球のカーブ、卓球のトップスピンは全て「マグヌス効果」、すなわち回転による循環の発生が原因です。ポテンシャル流れモデルは、この効果を理想化して理解する第一歩となります。

風力発電・タービン設計の初期検討:ブレード周りの複雑な流れを、まずは非粘性のポテンシャル流れとしてモデル化し、おおよその圧力分布や流線を把握するために用いられます。特に、循環の概念は渦流れを理解する上で重要です。

建築・環境風工学の基礎検討:建物や構造物周りの風の流れを大まかに把握する際に、ポテンシャル流れによる簡易モデルが使われることがあります。源・吸込みのモデルは、換気口付近の流れのイメージづくりに役立ちます。

よくある誤解と注意点

このシミュレーターを使い始めるときに、特に初学者がハマりがちな落とし穴がいくつかあります。まず第一に、「ポテンシャル流れは『完全な』流体のモデルだ」という点を忘れないでください。粘性を無視しているので、現実の流れとは根本的に異なります。例えば、循環Γを加えて揚力が発生する様子を見て「これで翼の揚力が全部説明できる!」と早合点する人がいますが、実際の翼では粘性によって生じる「循環」と「後縁条件」が複雑に絡み合っています。このツールはその核となる物理概念を抽出した「理想化された第一歩」だと捉えましょう。

第二に、パラメータ設定のオーダー(桁)に注意です。例えば、一様流速U∞を1 m/s、循環Γを10 m²/sなどと、無茶な組み合わせにすると、流線がぐちゃぐちゃになって物理的に意味のない結果になります。目安としては、Γの大きさはU∞ × R(半径)のオーダー程度から始めるのが安全です。例えばU∞=5 m/s、R=1 mなら、Γは5 m²/s前後から試してみると、きれいな非対称流れが観察できます。

第三に、「点源・吸込み(Q)」は数学的な抽象化ツールであることを理解しましょう。実務で「流量1 m²/sの吸込み」と言われてもピンときませんよね?これは、例えば直径10cmの円形ダクトから一定量の空気を吸い込む場合、その効果を単一の「点」で代表させたものと考えるのです。Qの値を変えると流線パターンが劇的に変わりますが、これは「特異点」が持つ数学的な強さの表れで、現実の広がりを持つ吸込み口の「一次近似」として解釈するのがコツです。

関連する工学分野

この一見シンプルな円柱まわりの流れの計算は、実は様々な先端工学分野の基礎として顔を出します。まず風力発電のブレード設計です。ブレードの各断面は「翼型」ですが、それを一次元化して円柱と類似の循環流れでモデル化する「揚力線理論」という手法があります。シミュレーターでΓを変えて揚力が変わるのを見る体験は、ブレードにかかるトルクをどう制御するかを考える第一歩です。

次に、船舶・潜水艦のプロペラ(スクリュー)設計です。プロペラが起こす流れは、一様流に循環と「らせん状の渦」を加えた複雑なものですが、その基本要素はここで学ぶ重ね合わせの原理と同じです。特に、プロペラ後方にできる「渦列」の予測には、ポテンシャル流れの考え方が応用されます。

さらに意外なところでは、マイクロ流体デバイス(Lab-on-a-chip)の流路設計にも通じます。極小スケールでは、粘性力が支配的になることが多いですが、複数の流入孔(点源)を配置して所望の流れパターンを作り出す設計思想は、このシミュレーターでQを複数配置したらどうなるか?と考えることそのものです。また、建築分野の風環境解析では、建物を単純な円柱や直方体に置き換えたポテンシャル流れの計算が、初期の風況予測に使われることがあります。

発展的な学習のために

このシミュレーターに慣れてきたら、次のステップとして「複素ポテンシャル」の世界に踏み込むことを強くお勧めします。実は、先ほどの流れ関数Ψと、もう一つの「速度ポテンシャルΦ」を組み合わせた複素数W(z)=Φ+iΨ(ここでz=x+iy)を考えると、数学的に驚くほど見通しが良くなります。一様流、循環、点源は全て、$$W(z) = U_\infty z + \frac{\Gamma}{2\pi i}\ln z + \frac{Q}{2\pi}\ln z$$のような複素関数の和で表現できます。これを使えば、円柱だけでなく、「ジュコーフスキー変換」という手法で翼型(ジャークー翼など)の周りの流れを計算できるようになります。これが「薄翼理論」への架け橋です。

学習の順序としては、1) このシミュレーターで感覚をつかむ → 2) 複素ポテンシャルと等角写像の基礎を学ぶ → 3) ジュコーフスキー翼を導出し、その圧力分布を計算してみる → 4) 現実の粘性の影響を学ぶため、「境界層理論」や「渦法」といった粘性を考慮した数値解法(CFD)の入門書に進む、という流れがスムーズです。このツールは、流体力学の美しい理想世界「ポテンシャル理論」の入口です。ここで得た「重ね合わせの原理」と「循環の物理的意味」の直感は、その後より複雑なCFDの結果を解釈する際の、確かな羅針盤となってくれるでしょう。