理論・主要公式
速度場 $\mathbf{F}=(P,Q)$ を各点で評価し、流れ粒子を $\dot x=P,\ \dot y=Q$(RK4)で移流させる。
ソース:$\mathbf{F}=(Qx,\,Qy)$、渦:$\mathbf{F}=(-Qy,\,Qx)$、サドル:$\mathbf{F}=(Qx,\,-Qy)$
発散 $\nabla\!\cdot\!\mathbf{F}=\partial_x P+\partial_y Q$、回転 $\nabla\!\times\!\mathbf{F}=\partial_x Q-\partial_y P$ をプローブで中心差分により実測。
既知解の検証:渦($Q\!=\!1$)→ 回転 $=2$・発散 $=0$、ソース($Q\!=\!1$)→ 発散 $=2$・回転 $=0$。
💬 解説ダイアログ
🙋
「ソース」って流体がどこから湧いてくるんですか?実際にはありえなくないですか?
🎓
2次元の理想化された数学モデルだ。ただ実用的な近似として、例えば地面に穴が開いて水が噴き出す「プランジングジェット」や、細管からの噴流などがソースに近い流れを示す。CAEでは「ポテンシャル流れ」として、粘性のない速い流れの概算に使われる。
🙋
「サドル」プリセットの流線がX字に見えますね。これは何ですか?
🎓
$\mathbf{F}=(Qx,-Qy)$ の双曲場だ。x方向に押し出し、y方向に引き込む。よどみ点(流れがぶつかって速度0になる点)まわりの局所流れがまさにこの形で、力学系の不動点安定性解析の基本例にもなっている。発散も回転も0なのをプローブで確かめてみてほしい。
🙋
「渦」と「ソース」で発散・回転の数値がどう変わるんですか?
🎓
渦 $\mathbf{F}=(-Qy,Qx)$ は回転$=2Q$・発散$=0$、ソース $\mathbf{F}=(Qx,Qy)$ は発散$=2Q$・回転$=0$になる。$Q=1$ なら渦の回転は2、ソースの発散は2だ。プローブの読み値が理論値と一致するのを確認できる。
本ツールでは、2次元非圧縮性流体の流れを、基本解の重ね合わせによりモデル化する。各基本流は流れ関数 \(\psi\) または速度ポテンシャル \(\phi\) で記述され、速度場 \(\boldsymbol{v} = (u, v)\) は \(\boldsymbol{v} = \nabla \phi\) または \(u = \partial \psi / \partial y\), \(v = -\partial \psi / \partial x\) により算出される。例えば、強さ \(\Gamma\) の渦は \(\psi = \frac{\Gamma}{2\pi} \ln r\) で与えられ、その周囲に循環流を生じる。また、強さ \(Q\) のソースは \(\phi = \frac{Q}{2\pi} \ln r\) により放射状の流れを生成する。これらの基本解を線形結合した全流れ関数 \(\Psi = \sum \psi_i\) により、複雑な流線パターンが形成される。流線は \(\Psi = \text{const}\) の等高線として定義され、ユーザーがキャンバス上に指定した出発点から数値積分により追跡される。このモデルにより、渦対やよどみ点などの物理現象を直感的に理解できる。
産業での実際の使用例
航空機業界では、ボー風力発電事業者が翼周りの空気流れ解析の初期検討に本ツールを活用。翼断面(2次元)のソースと渦の重ね合わせで揚力発生メカニズムを即座に可視化し、設計案の絞り込みに貢献。自動車業界では、トヨタがラジエーターグリル周りの冷却風の流線を確認し、整流板の配置最適化に利用。また、風力発電では、風力タービンメーカーがブレード断面周りの流れを可視化し、失速特性の予測に役立てている。
研究・教育での活用
大学の流体力学講義では、東大や東工大がポテンシャル流理論の理解促進に導入。学生がソース・シンク・渦のパラメータを変更しながら流線パターンを観察し、重ね合わせの原理を直感的に習得。研究では、宇宙機関が翼後縁の渦放出パターンの予備解析に使用し、実験計画の効率化に貢献。
CAE解析との連携や実務での位置付け
本ツールは、本格的なCAE(例:CAEソフト Fluent、OpenFOAM)による3次元解析の前段階として位置づけられる。設計初期段階で流れの大まかな傾向を把握し、計算負荷の高いCAEの解析条件(境界条件やメッシュ生成の方針)を迅速に決定。実務では、CAEの結果検証やトラブルシューティングの際に、基本流の重ね合わせで現象を分解し、原因特定を支援する補助ツールとして機能する。
「流線は流体の粒子の軌跡(流跡線)と常に一致する」と思いがちですが、実際には定常流(時間変化しない流れ)でのみ一致します。非定常流では、流線は瞬間瞬間の速度場の接線を結んだものであり、粒子の実際の移動経路とは異なります。特に渦や一様流を時間変化させて重ね合わせる場合、この違いに注意が必要です。
「ソース(湧き出し)とシンク(吸い込み)を重ねれば必ず双極子(ダブレット)になる」と思いがちですが、実際には両者の強さが等しく、かつ無限に近づけた極限でのみ双極子とみなせます。有限の距離で配置した場合は単なるソース・シンク対となり、流線の形状も異なります。ツール上で距離を変えながら挙動を確認する際は、この近似関係を意識してください。
「キャンバス上で流線の出発点を増やせば、すべての流れの様子を完全に把握できる」と思いがちですが、実際には流線はあくまで瞬間の速度場の可視化であり、出発点の位置によっては重要な循環や淀み点を捉えきれない場合があります。特に渦の中心付近や流れの分岐点では、出発点のわずかな違いで流線の描画結果が大きく変わることがあるため、複数の視点から配置を試すことが重要です。
具体的な計算例
渦プリセット $\mathbf{F}=(-Qy,Qx)$ を選び、強さ Q=1 にしてカーソルを任意の点へ動かすと、回転 ∇×F は理論値 2Q=2、発散 ∇·F は 0 と表示されます。次にソースプリセット $\mathbf{F}=(Qx,Qy)$・Q=1 に切り替えると、発散が 2Q=2、回転が 0 になります。これは中心差分による数値計算が解析的な既知解と一致することを示します。サドルプリセットでは発散・回転ともに 0 で、流線がX字(双曲線)を描きます。