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

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

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

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

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

理論・主要公式
流れ関数(円柱 + 循環):
$\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$

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

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

よくある質問

理想流体(非粘性・非圧縮)のポテンシャル流れでは、一様流中の物体に抗力が働かないという理論的結果を指します。本シミュレーターで循環Γ=0、源Q=0に設定すると、円柱前後の圧力係数分布が対称になり、抗力がゼロとなる様子を視覚的に確認できます。
循環強度Γ(ガンマ)の値を正または負に設定してください。Γを大きくすると、円柱周りの流線の非対称性が強まり、揚力が発生します。同時に圧力係数カラーマップで、円柱の上下で圧力差が生じる様子を確認できます。
Qを正の値に設定すると円柱から流体が湧き出す「源」、負の値に設定すると流体が吸い込まれる「吸込み」を表現します。これにより、吹出しや吸込みがある場合の流線の変化や、物体表面の圧力分布への影響を可視化できます。
バグではありません。このシミュレーションはポテンシャル流れを仮定しており、円柱内部も計算領域に含まれています。実際の流体では円柱内部に流れは侵入しませんが、本ツールでは理論上の流れ関数を全域で計算しているため、内部にも流線が描画されます。

実世界での応用

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

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

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

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

よくある誤解と注意点

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

第二に、パラメータ設定のオーダー(桁)に注意です。例えば、一様流速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の値を変えると流線パターンが劇的に変わりますが、これは「特異点」が持つ数学的な強さの表れで、現実の広がりを持つ吸込み口の「一次近似」として解釈するのがコツです。

使い方ガイド

  1. 流速U(m/s)をスライダーまたは数値入力で0.5~10m/sの範囲で設定します
  2. 円柱半径R(mm)を1~50mmで調整し、レイノルズ数Re=ρUD/μの値を確認します
  3. 循環Γ(m²/s)を-10~+10の範囲で変更してマグヌス効果による揚力変化を観察します
  4. 格子解像度G(分割数)を増加させて計算精度を向上させ、圧力係数Cpの分布をリアルタイム更新します
  5. 流線とCpカラーマップから上端・下端の圧力係数、揚力係数CL、抗力係数CDを読み取ります

具体的な計算例

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

実務での注意点