› ポテンシャル流れシミュレーター 戻る
流体解析

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

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

流れパターン
パラメータ
一様流速 U∞1.0 m/s
循環 Γ0.0 m²/s
代表長さ a1.00 m
表示
ライブ数値
1.00
一様流速 U∞ [m/s]
0.0
循環 Γ [m²/s]
2
淀み点の数
2.00
瞬間最大流速 |V| [m/s]
0.00
揚力係数 CL
0.0
揚力 L=ρU∞Γ [N/m]
−3.00
Cp 上端
−3.00
Cp 下端
ポテンシャル流れ・粒子アニメーション
高速・低圧 (Cp<0)
低速・高圧 (Cp→1)

トレーサー粒子が速度場を移流。色=圧力係数Cp、黄=淀み点、緑=揚力ベクトル。プリセットで一様流・湧き出し・円柱・自由渦・ランキン卵形を切替。

理論・主要公式
複素ポテンシャルの重ね合わせ(一様流+二重子+渦):
$w(z)=U_\infty z+\dfrac{\mu}{2\pi z}-\dfrac{i\Gamma}{2\pi}\ln z,\quad u-iv=\dfrac{dw}{dz}$

円柱半径(二重子強さ μ):$a=\sqrt{\dfrac{\mu}{2\pi U_\infty}}$

流れ関数(円柱+循環):$\Psi = U_\infty r\!\left(1-\frac{a^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$

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

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

よくある質問

理想流体(非粘性・非圧縮)のポテンシャル流れでは、一様流中の物体に抗力が働かないという理論的結果を指します。本シミュレーターの「円柱(二重子)」で循環Γ=0にすると、円柱の前後で圧力係数の分布が対称になり(前後の淀み点で Cp=1、上下端で Cp=−3)、圧力による抗力がゼロになる様子を確認できます。
循環強度Γ(ガンマ)の値を正または負に設定してください。Γを大きくすると、円柱周りの流線の非対称性が強まり、揚力が発生します。同時に圧力係数カラーマップで、円柱の上下で圧力差が生じる様子を確認できます。
Q は「湧き出し」と「ランキン卵形」のプリセットで使います。Q を正にすると一点から流体が湧き出す「源」、負にすると吸い込まれる「吸込み」になります。一様流と湧き出しを重ねるとランキン半物体、湧き出しと吸い込みを組み合わせるとランキン卵形(閉じた物体)のまわりの流れになり、物体の大きさや淀み点の位置が Q と U∞ で決まる様子を確認できます。
円柱の内部は速度ゼロとして扱い、流線・トレーサー粒子・圧力色は円柱の外側だけに描きます。複素ポテンシャルの式そのものは円柱の内側でも値を持ちますが、物理的な意味はないため表示していません。

実世界での応用

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

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

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

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

よくある誤解と注意点

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

第二に、循環の大きさと淀み点の位置に注意です。円柱の淀み点は、循環が小さいうちは表面上を下(または上)へ移動し、Γ が 4πU∞a を超えると表面から離れて流れの中に1つだけ現れ、円柱のまわりを周回する閉じた流線ができます。U∞=1 m/s、a=1 m なら境目は Γ≈12.6 m²/s で、Γ=15 m²/s では淀み点は中心から約 1.8 m の位置に移ります。これはポテンシャル流れとして正しい解ですが、実際の回転円柱でここまでの循環を得るには非常に速い回転が必要です。まずは Γ を U∞a と同じくらい(例えば U∞=5 m/s、a=1 m なら 5 m²/s 前後)から試すと、表面上の淀み点の移動と揚力の関係が見やすくなります。

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

使い方ガイド

  1. 一様流速 U∞ を 0.2〜5 m/s の範囲で設定します
  2. 代表長さ a(円柱では半径、ランキン卵形では湧き出しと吸い込みの間隔の半分)を 0.4〜2.0 m で調整します
  3. 「円柱(二重子)」で循環Γを −18〜+18 m²/s の範囲で変え、マグヌス効果による揚力の変化を観察します
  4. 流れパターン(一様流/湧き出し/円柱(二重子)/自由渦/ランキン卵形)を切り替えて、流線と Cp カラーマップの変化を比較します
  5. 上端・下端の圧力係数、揚力係数 CL、揚力 L、淀み点の数を読み取ります(抗力はダランベールのパラドックスにより常に0なので表示しません)

具体的な計算例

半径R=1.0m(直径D=2.0m)の円柱に流速U=5m/s(大気中ρ=1.225kg/m³)の空気流が作用する場合:レイノルズ数Re=ρUD/μ≈6.8×10⁵となります。無循環(Γ=0)では理想流体のダランベールのパラドックスに従いCD=0ですが、Γ=+5m²/sを加えると、クッタ・ジュコーフスキーの定理 L'=ρUΓ=1.225·5·5≈30.6N/m の揚力が発生し、CL=Γ/(U·R)=1.0となります。同時にCp(上端)≈-3.65、Cp(下端)≈-2.38へと変化し、流線の非対称化を確認できます。

実務での注意点