流線 — CAE用語解説

良い着眼点だ。これらは時間の捉え方が根本的に異なる。流線は、ある瞬間に、その場の速度ベクトルに接するように引かれた仮想的な線だ。一方、流跡線はある1つの流体粒子が時間とともに通る軌跡、流脈線はある瞬間に特定の点を通過する全ての粒子が描く線だ。定常流では三つは一致するが、非定常流では異なる。例えば、飛行機の翼周りの定常解析では流線を見るが、ヘリコプターのローターが回転する非定常解析では流跡線を見る必要がある。

速度場を

その通りだ。支配方程式である非圧縮性ナビエ-ストークス方程式は、

核心的な質問だ。主要なアルゴリズムは2つある。1つ目はオイラー法の応用で、最も単純だ。あるシード点

2つ目で、標準的に使われるのが4次/5次ルンゲ-クッタ法（RK45）だ。これは、1ステップで複数回速度を評価し、高次の精度を達成する。例えば、Ansys FluentやOpenFOAMの`streamLine`関数は、デフォルトでこのRK法を採用している。ステップサイズ

これが重要なポイントで、補間が必要になる。FVM（有限体積法）ベースのソフトでは、セル中心の値から節点への補間（例えば線形補間）を行い、節点値から任意の点での値を求める（形状関数を使った内挿）。Ansys CFXでは要素内で二次補間を行う。補間の精度が、流線の滑らかさと物理的正確さに直結する。粗いメッシュで細かい渦を追おうとすると、補間誤差で流線がガタガタになったり、存在しない構造が見えたりする。

主な終了条件は4つだ。1) 計算領域の境界に到達。2) 速度が極めて小さな領域（

目的に応じて戦略が変わる。例えば、流入条件の均一性を見たいなら、流入境界面上に等間隔で数十点配置する。剥離渦を可視化したいなら、剥離が起こる翼や車体のエッジの直後に線状に配置する。Ansys Fluentの「Line/Rake」シード機能がこれだ。循環域（デッドゾーン）を探すには、領域全体にランダムに数百点配置する「Point Cloud」が有効だ。重要なのは「見たい物理現象の上流にシードを置く」という原則だ。

速度の大きさは基本的だが、それだけでは不十分だ。実務では、圧力（特に全圧損失を見る）、渦度（渦の強さ）、温度（熱の輸送経路）、乱流エネルギー（混合の度合い）で色分けすることが多い。例えば、ポンプ内部流れでは、流線を圧力で色分けすると、キャビテーションが発生する低圧域を流線が通過しているか一目でわかる。色の範囲は、物理的に意味のある値に設定すべきだ。例えば温度なら、流入温度と壁面温度の間にする。

必ずしも計算誤差ではない。可視化ソフトがメッシュを表示せず、流線だけを表示している場合によくある。しかし、壁面を貫通する流線は明らかな問題だ。原因は主に二つ。1) 壁面境界条件（通常は無滑り条件）が正しく設定されていない。2) メッシュが粗すぎて、壁面近傍の速度勾配を解像できておらず、補間によって壁面でゼロでない速度が計算されてしまう。後者を確認するには、壁面のせん断応力やy+の分布をチェックする必要がある。

流線そのものから直接数値を出すのは難しい。流線は定性的評価の王者だ。定量的評価には、流線の元になった速度場を使う。例えば、熱交換器のフィン間流れの均一性を評価するなら、下流の特定断面における速度分布の標準偏差や最大速度/平均速度比を計算する。流線は、その不均一なパターンがどこで、なぜ発生しているのかを視覚的に理解するための手がかりとして使う。可視化（流線）と定量データ（グラフ、コンター）の組み合わせが鉄則だ。

大きな違いがある。まずシード設定だ。Fluentは従来から「Point」、「Line」、「Surface」、「Volume」など多様だが、設定画面がやや分かりにくい。Star-CCM+は「パーツ」と呼ばれる幾何学的要素（点、線、面、体）に直接シードを設定する直感的なインターフェースで、特に「Field Function」を使って複雑なシード分布を関数で定義できるのが強みだ。

両者とも高次のルンゲ-クッタ法を実装しているが、並列計算への対応に差がある。Fluentは流線計算自体は並列化されているが、非常に多数（数万本）の流線を生成する場合、Post処理段階で重くなる傾向がある。Star-CCM+はシミュレーションとポスト処理が統合された環境で、流線の生成も効率的に並列分散処理される。1000万セルモデルで数万本の流線を生成する場合、Star-CCM+の方が高速なことが多い。

ParaViewは非常に強力だ。FluentやOpenFOAMのデータを読み込み、「Stream Tracer」フィルタを使う。アルゴリズムはRK2, RK4, RK45から選べ、シードも多様に設定できる。商用ソフトに引けを取らないが、インターフェースの学習コストが高い点と、超大型データ（数十GB以上）の扱いでメモリ管理に気を使う点が難点だ。VisItも同様だが、ParaViewほどメジャーではない。研究機関や予算制約のあるプロジェクトではParaViewが第一選択肢になる。

COMSOLの特徴は、流線を「粒子追跡」モジュールの一部として扱う点だ。これにより、流線に沿った任意の変数（温度、化学種濃度、電位など）の積分値を容易に計算できる。例えば、マイクロ流体チップ内の特定経路での反応物質の平均滞留時間を、流線に沿って計算できる。これは純粋なCFDソフトにはない強みで、マルチフィジックス解析に特化している。ただし、大規模乱流解析のような専用CFDソフトの性能には及ばない。

最も多い原因は二つ。第一に、出口境界条件だ。流出（Outflow）や圧力出口（Pressure Outlet）では、流線はその面で計算を終了する。これが「切れる」正しい挙動だ。問題は、流線が途中で消える場合で、これはその位置のセルデータが欠損（NaNや異常値）している可能性が高い。第二に、メッシュ品質だ。極端に細長いセルや歪んだセルでは、補間計算が失敗し、速度ベクトルが求まらなくなる。対策としては、まず該当箇所のメッシュを確認し、セル品質（Skewness > 0.95 など）を改善する。

メッシュを細かくするのは根本解決になるが、計算コストがかかる。まず試すべきは流線計算ソルバーの設定変更だ。ステップサイズ

これはよくある。原因は主にデータの補間方法とシード点の厳密な位置だ。Fluentは独自の内部補間法を使い、ParaViewはVTKライブラリを通じてデータを読み、そのライブラリの補間法を使う。シード点の座標値も、丸め誤差レベルで異なる可能性がある。信用すべきは大まかな流れのパターン（剥離の有無、大きな渦の位置）だ。細部の違いに一喜一憂する必要はない。ただし、圧力や速度の極値が流線上で大きく異なる場合は、データのエクスポート/インポート過程で欠損や変換誤差が生じている可能性を疑うべきだ。

それで正しい。非定常流では、流線はあくまで「瞬間の速度場」に基づく。時間が経てば速度場自体が変わるので、流線の形も連続的には変化しない。むしろ、流線が滑らかに変形するアニメーションを見せられたら、それは流跡線を表示しているか、または間違った可視化をしている可能性が高い。非定常流の時間発展を見たいなら、粒子追跡（Particle Tracking） 機能を使って流跡線をアニメーション化するか、複数の時刻の流線を並べて比較するのが正しい方法だ。Ansys EnsightやTechPlotは、このような時系列可視化に強い。