2Dベクトル場可視化 戻る
解析ツール

2次元ベクトル場・流線可視化ツール

一様流・ソース・シンク・渦・サドル・双極子といった基本場を切り替え、流れ粒子をリアルタイムに移流(RK4)。矢印・流線・カラーマップを重ね、カーソルを動かすとプローブ点の |F|・発散・回転をその場で計測します。

流れのプリセット

場の強さ Q1.0
矢印グリッド密度15

一時停止中はスライダーを動かすと結果が即座に更新されます。

計算結果
プローブ |F|
発散 ∇·F
回転 ∇×F
粒子平均速度
ベクトル場と流線
理論・主要公式

速度場 $\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だ。プローブの読み値が理論値と一致するのを確認できる。

よくある質問

Q. ベクトル場の「発散」と「回転」とは何ですか?
A. 発散(div)はある点での流体の「湧き出し度」で、∇·v=∂u/∂x+∂v/∂yで計算します。非圧縮流では∇·v=0(質量保存)。回転(rot/curl)は「渦の強さ」で、ω=∂v/∂x-∂u/∂yで計算します。ポテンシャル流れでは∇×v=0(非回転)が条件です。
Q. 電磁場とベクトル場の関係は?
A. 電場Eや磁場Bもベクトル場であり、マクスウェル方程式はベクトル解析(∇·E、∇×E等)で記述されます。正電荷の周りの電場はソースと、磁石の周りは双極子とそっくりな形をしています。流体力学と電磁気学は数学的に深く対応しています。
Q. 有限要素法(FEM)ではどう使いますか?
A. FEMの結果は各節点での速度・応力・温度ベクトルとして出力されます。このようなベクトル場可視化(矢印表示・流線表示)はParaViewやEnSightなどの後処理ソフトで行います。流線を見ることで流れの構造(再循環域・渦など)を把握できます。
Q. 流線・流跡線・煙跡線の違いは?
A. 流線(streamline)はある瞬間の速度場の接線曲線、流跡線(pathline)は1つの流体粒子の軌跡、煙跡線(streakline)は同じ点を通過した全粒子の軌跡です。定常流(時間変化なし)では3つは一致しますが、非定常流では異なります。

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の結果検証やトラブルシューティングの際に、基本流の重ね合わせで現象を分解し、原因特定を支援する補助ツールとして機能する。

よくある誤解と注意点

「流線は流体の粒子の軌跡(流跡線)と常に一致する」と思いがちですが、実際には定常流(時間変化しない流れ)でのみ一致します。非定常流では、流線は瞬間瞬間の速度場の接線を結んだものであり、粒子の実際の移動経路とは異なります。特に渦や一様流を時間変化させて重ね合わせる場合、この違いに注意が必要です。

「ソース(湧き出し)とシンク(吸い込み)を重ねれば必ず双極子(ダブレット)になる」と思いがちですが、実際には両者の強さが等しく、かつ無限に近づけた極限でのみ双極子とみなせます。有限の距離で配置した場合は単なるソース・シンク対となり、流線の形状も異なります。ツール上で距離を変えながら挙動を確認する際は、この近似関係を意識してください。

「キャンバス上で流線の出発点を増やせば、すべての流れの様子を完全に把握できる」と思いがちですが、実際には流線はあくまで瞬間の速度場の可視化であり、出発点の位置によっては重要な循環や淀み点を捉えきれない場合があります。特に渦の中心付近や流れの分岐点では、出発点のわずかな違いで流線の描画結果が大きく変わることがあるため、複数の視点から配置を試すことが重要です。

使い方ガイド

  1. プリセット(一様流/ソース/シンク/渦/サドル/双極子)を選び、場の種類を切り替えます。流れ粒子がRK4で移流し、矢印・流線がリアルタイム更新されます
  2. 「場の強さ Q」スライダー(0.2〜3)で場の強度を、「矢印グリッド密度」スライダー(6〜26本)で矢印の本数を調整します
  3. 「流れ粒子アニメーション」「速度ベクトル矢印」「流線」のチェックボックスで表示を切り替え、流れのパターンを比較します
  4. カーソルをキャンバス上で動かすとプローブ点が移動し、|F|・発散 ∇·F・回転 ∇×F がライブ計測されます。「速度大きさ(カラーマップ)」で高速域と低速域を色分けできます

具体的な計算例

渦プリセット $\mathbf{F}=(-Qy,Qx)$ を選び、強さ Q=1 にしてカーソルを任意の点へ動かすと、回転 ∇×F は理論値 2Q=2、発散 ∇·F は 0 と表示されます。次にソースプリセット $\mathbf{F}=(Qx,Qy)$・Q=1 に切り替えると、発散が 2Q=2、回転が 0 になります。これは中心差分による数値計算が解析的な既知解と一致することを示します。サドルプリセットでは発散・回転ともに 0 で、流線がX字(双曲線)を描きます。

実務での注意点

  1. 「双極子」プリセットは等強度のソース(+x側)とシンク(−x側)の対で、両者を無限に近づけた極限がダブレット(円柱まわり流れの構成要素)に一致します
  2. 「矢印グリッド密度」を上げる(最大26本)ほど描画負荷が増えるため、確認後は15本前後での運用が快適です
  3. 本ツールは各点で定義された定常ベクトル場を扱い、粒子は場に沿って移流します。場そのものが時間変化する非定常流(渦の合体など)は表現しません