二次流れ
理論と物理
概要
二次流れって翼列の中で主流と違う方向に流れるやつですよね？
そう。翼間流路内で主流方向と垂直な速度成分を持つ流れを二次流れと呼ぶ。端壁（ハブ・シュラウド）の境界層が翼列の圧力差で横に曲げられることで発生する。
主な渦構造
どんな渦構造があるんですか？
代表的な渦構造を挙げよう。
| 渦の名称 | 発生メカニズム | 影響 |
|---|---|---|
| 通路渦 (Passage Vortex) | 端壁BLが翼間圧力差で巻き上がる | 主要な二次流れ損失源 |
| 馬蹄渦 (Horseshoe Vortex) | 翼前縁で端壁BLが分岐 | SS側が通路渦に合流 |
| コーナー渦 (Corner Vortex) | 翼面-端壁の交差部で発生 | 剥離を誘発 |
| チップ漏れ渦 (Tip Leakage Vortex) | 翼端隙間からの漏れ流れ | 効率低下の主因（回転翼） |
| スクレイパー渦 | シュラウド壁面の相対運動 | 遷音速段で顕著 |
通路渦と馬蹄渦の関係を教えてください。
馬蹄渦は翼前縁で2つに分かれる。圧力面側（PS leg）は隣の翼に向かい、吸い込み面側（SS leg）はそのまま通路渦に巻き込まれて増強される。この合流した渦が通路渦の主体だ。
二次流れ損失の定量化
二次流れ損失はどのくらいですか？
タービン翼列の全損失の30～50%が二次流れに起因すると言われている。CFDではエントロピー生成率で可視化するのが有効だ。
この量をCFD-Postで計算して体積積分すれば、翼面損失・端壁損失・チップ漏れ損失を分離して評価できる。
ターボ機械の二次流れ理論——Hawthorne（1955年）と馬蹄渦の体系化
ターボ機械翼列の「二次流れ（Secondary Flow）」を理論的に整理したのは英国のW.R. Hawthorne（1955年）だ。Hawthorneは翼の前縁前方で入射境界層渦度が馬蹄形に分裂・伸長し「馬蹄渦（Horseshoe Vortex）」を形成する過程を渦度輸送方程式で記述した。この理論は翼列の端壁損失のメカニズムを初めて定量的に説明し、ターボ機械設計における端壁処理の重要性を示した先駆けとなった。Hawthorne自身はケンブリッジ大学の工学教授として複数世代の航空宇宙エンジニアを育て、後の現代ターボCFDの基礎を担う研究者を多数輩出した。彼の二次流れ理論は現代のCFDで数値的に検証されており、馬蹄渦の形状・強さのCFD予測とHawthorneの古典理論の対応が今も研究テーマになっている。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 速度 $u$ | m/s | 入口条件で体積流量から換算する際、断面積の単位に注意 |
| 圧力 $p$ | Pa | ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 |
| 密度 $\rho$ | kg/m³ | 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C |
| 粘性係数 $\mu$ | Pa·s | 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 |
| レイノルズ数 $Re$ | 無次元 | $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 |
| CFL数 | 無次元 | $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結 |
数値解法と実装
渦の同定手法
CFDの結果から渦構造をどうやって抽出するんですか？
複数の渦同定法がある。
| 手法 | 定義 | 特徴 |
|---|---|---|
| Q基準 | 渦度テンソルの大きさ > ひずみ速度テンソルの大きさ | 最も広く使用、CFD-Postに標準搭載 |
| λ2基準 | 圧力ヘッシアンの第2固有値が負 | せん断効果を除去、より正確 |
| ヘリシティ | $H = \mathbf{v} \cdot \boldsymbol{\omega}$ | 渦の回転方向を判別可能 |
| 壁面限界流線 | 壁面せん断応力の方向 | 剥離線・付着線の同定 |
Q基準の等値面を表示するのが一番手軽ですか？
そう。CFD-PostでQ=正の値の等値面を全圧損失係数や渦度で着色すると、通路渦やチップ漏れ渦の3D構造が一目で分かる。
メッシュ要件
二次流れを正確に予測するにはどのくらいのメッシュが必要ですか？
端壁付近のメッシュ密度がカギだ。
- 端壁y+: < 1（Low-Re SSTモデル使用時）
- 端壁プリズム層: 15～20層
- 翼面-端壁の交差部: メッシュ細分化（コーナー渦を捕える）
- 通路中央のスパン方向: 40セル以上
翼面だけでなく端壁のメッシュ品質が二次流れの予測精度を左右する。
TurboGridで端壁のメッシュを細かくするにはどうしますか？
TurboGridのBoundary Layer Refinementで端壁（Hub/Shroud）に対して専用のプリズム層を設定する。翼面と端壁の両方に独立にy+制御が可能だ。
乱流モデルの影響
二次流れの予測に乱流モデルの違いは大きいですか？
SST k-omegaとk-epsilonで通路渦の位置やサイズに有意な差が出る。SSTのほうが端壁近傍の逆圧力勾配を正確に捉えるから、通路渦の強さと位置がより実験に近い。LESならさらに非定常の渦構造まで解像できるが、計算コストが2桁以上増える。
ターボ二次流れのCFD数値手法——コーナー渦の予測とSSTモデルの精度
ターボ機械の翼エンドウォール付近に形成される「コーナー渦（Corner Vortex）」の予測精度は使用する乱流モデルに強く依存する。標準k-εモデルはせん断力の等方性を仮定するため、強い曲率と圧力勾配が同時に働くエンドウォール周辺で渦の強さを著しく過小評価する（実験比40〜60%の過小評価）。SSTモデルはk-εとk-ωを切り替えるブレンディング関数でエンドウォール付近の精度を改善するが、それでも二次流れの渦の位置を±5〜10%の範囲でずらすことがある。最高精度には差分レイノルズ応力モデル（DRSM）かLESが必要だが、設計サイクルでの計算コストが問題だ。実務では「SSTで二次流れの傾向を把握し、最終設計のみLES検証」という段階的精緻化アプローチが多く採用されている。
風上差分（Upwind）
1次風上: 数値拡散が大きいが安定。2次風上: 精度向上するが振動のリスク。高レイノルズ数流れでは必須。
中心差分（Central Differencing）
2次精度だが、Pe数 > 2で数値振動が発生。低レイノルズ数の拡散支配流れに適する。
TVDスキーム（MUSCL、QUICK等）
リミッタ関数により数値振動を抑制しつつ高精度を維持。衝撃波や急勾配の捕捉に有効。
有限体積法 vs 有限要素法
FVM: 保存則を自然に満足。CFDの主流。FEM: 複雑形状・マルチフィジックスに有利。SPH等のメッシュフリー法も発展中。
CFL条件（クーラン数）
陽解法: CFL ≤ 1が安定条件。陰解法: CFL > 1でも安定だが、精度と反復回数に影響。LES: CFL ≈ 1を推奨。物理的意味: 1タイムステップで情報が1セル以上進まないこと。
残差モニタリング
連続の式・運動量・エネルギーの各残差が3〜4桁低下で収束と判断。質量保存の残差は特に重要。
緩和係数
圧力: 0.2〜0.3、速度: 0.5〜0.7が一般的な初期値。発散する場合は緩和係数を下げる。収束後は上げて加速。
非定常計算の内部反復
各タイムステップ内で定常解に収束するまで反復。内部反復数: 5〜20回が目安。残差がタイムステップ間で変動する場合は時間刻みを見直す。
SIMPLE法のたとえ
SIMPLE法は「交互に調整する」手法。まず速度を仮に求め（予測ステップ）、その速度で質量保存が満たされるよう圧力を補正し（補正ステップ）、補正された圧力で速度を修正する——このキャッチボールを繰り返して正解に近づく。2人で棚を水平にする作業に似ている：片方が高さを合わせ、もう片方がバランスを取り、これを交互に繰り返す。
風上差分のたとえ
風上差分は「川の流れに立って上流の情報を重視する」手法。川の中にいる人が下流を見ても水の出所は分からない——上流の情報が下流を決めるという物理を反映した離散化手法。精度は1次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
端壁コンタリング
二次流れ損失を減らす方法はありますか？
非軸対称端壁コンタリング（Non-Axisymmetric Endwall Contouring）が最も効果的な手法の一つだ。端壁の形状を翼間で凸凹させることで、端壁付近の圧力分布を変え、二次流れを抑制する。
どのくらい損失を減らせるんですか？
翼列の二次流れ損失を10～30%低減、段効率で0.5～1.5ポイント改善が報告されている。
CFDによる端壁最適化
CFDで端壁コンタリングの最適化はできますか？
できる。端壁形状をフーリエ級数やスプライン曲面でパラメータ化し、CFDベースの最適化を行う。
1. 端壁を翼間の周方向×軸方向のグリッドで離散化（5×5～10×10制御点）
2. 各点の径方向変位を設計変数とする
3. 目的関数: 全圧損失係数の最小化 or 段効率の最大化
4. 遺伝的アルゴリズムやベイズ最適化で探索
設計変数が多くて大変そうですね。
25～100変数になるから、代理モデル（Kriging）を使って評価回数を数百回に抑えるのが実用的だ。FINE/DesignやoptiSLangにはこうした代理モデルベース最適化が実装されている。
チップ漏れ損失の低減
チップ漏れ損失を減らす方法は？
| 手法 | 原理 | 効果 |
|---|---|---|
| チップクリアランス縮小 | 漏れ流量の直接低減 | 最も効果的だが製造・運用の制約 |
| 翼端スクイーラーチップ | リム構造で漏れ流れを絞る | 0.5～1%効率改善 |
| ウィングレット | 翼端に突起を付けて漏れ渦を抑制 | 0.5～1.5%効率改善 |
| ケーシングトリートメント | シュラウドに溝を刻んで流れを制御 | 主に安定性向上 |
軸流タービンの翼端クリアランス——クリアランスギャップ流れのCFD設計
軸流タービンのロータブレード先端とケーシング間の「翼端クリアランス（Tip Clearance）」は設計上重要なパラメータだ。クリアランスが大きいと「翼端漏れ流れ（Tip Leakage Flow）」が圧力面から吸い込み面へ通過し、大きな漏れ渦（Tip Leakage Vortex）を形成して段効率を低下させる。一般にクリアランスを設計値の1%から2%に増やすと段効率が0.5〜1ポイント低下する経験則がある。CFD解析では製造公差による翼端形状のバリエーション（±0.1mm）を含めた感度解析が設計信頼性評価に重要で、SquealerチップやA型先端冷却穴の有無が漏れ流れ量に与える影響を事前評価する。三菱重工エネルギーが公開した研究では、Squealerチップの最適形状をCFD最適化で決定し、翼端損失を20%低減した成果が発表されている。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
実験的検証手法
二次流れのCFD予測はどうやって実験検証するんですか？
代表的な実験手法を挙げよう。
| 手法 | 測定対象 | 特徴 |
|---|---|---|
| 5孔ピトーブ管トラバース | 3D速度場+全圧・静圧 | 翼列出口の定番手法 |
| PIV (粒子画像速度計) | 2D/3D瞬時速度場 | 渦構造の可視化に優れる |
| 油膜法 | 壁面限界流線 | 剥離線の同定 |
| 感圧塗料 (PSP) | 翼面静圧分布 | 非定常圧力分布も可能 |
5孔プローブの結果とCFDは合いますか？
全圧損失の周方向分布は概ね一致するが、通路渦コアの位置がスパン方向に1～2%ずれることがある。これはRANSの乱流粘性が渦の拡散を過大評価するためだ。LESで改善される。
CFDベンチマーク
二次流れ予測の標準ベンチマーク問題はありますか？
Durham大学のLinear Cascade（低速タービン翼列）やAAA委員会のテストケースがよく使われる。NUMECA や Ansys が公開しているバリデーションケースも参考になる。
ソルバー間比較
ソルバーによって二次流れの予測に差が出ますか？
同じメッシュ・同じ乱流モデルならソルバー間の差は小さい。むしろメッシュ品質と乱流モデルの選択のほうが支配的だ。ただしMixing Planeの実装が異なると、多段での二次流れの累積に差が出ることがある。
ターボ二次流れCFDの専用ツール——CFX vs Fine/Turboの乱流モデル実装差
ターボ機械の二次流れ予測でANSYS CFXとNUMECA FINE/Turboの乱流モデル実装を比較すると、CFXはSST（Shear Stress Transport）モデルのターボ機械実績が豊富で翼列解析のデファクト標準として機能している。SSTモデルの圧力勾配補正項の実装がわずかに異なるため、エンドウォール二次流れの予測では両ツールで効率予測が0.5〜1ポイント異なるケースが報告されている。FINE/TurboはCurvature-corrected SSTを標準で提供しており、旋回が強い流路（ラジアルコンプレッサー）での二次流れ予測でCFXより若干精度が高いという第三者評価がある。結論として、使用するツールが決まった後はそのツールの乱流モデル実装に合わせた実験ベースの精度検証データを蓄積することが長期的な設計信頼性の基盤になる。
選定で最も重要な3つの問い
- 「何を解くか」：二次流れに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
RANSの限界
RANSだと二次流れのどこが不正確ですか？
通路渦やチップ漏れ渦の非定常変動は RANS では時間平均されてしまう。実際にはこれらの渦は蛇行や崩壊を繰り返しており、瞬時の損失分布は時間平均と大きく異なる。
LESの適用
LESで二次流れを解析した例はありますか？
学術的にはWall-Resolved LESでタービン翼列の二次流れを解析した研究がある。チップ漏れ渦の蛇行周波数やコーナー剥離の非定常振動が初めて解明された。ただし計算コストは RANS の $10^4$ 倍以上だ。
実務的なアプローチ: SASとSDES
実務で使えるスケール解像手法はありますか？
SAS の結果は RANS とどう違いますか？
端壁付近のエントロピー生成分布がより集中した構造を示す。RANS では数値拡散で「ぼやけた」渦がSASではシャープなコアとして表現される。損失の絶対値は RANS と大差ないが、損失の発生メカニズムの理解が深まる。
トラブルシューティング
メッシュ起因の偽渦
二次流れの解析でメッシュが原因の偽の渦構造が出ることはありますか？
ある。特に以下のケースで注意が必要だ。
1. 端壁O-gridの結合部: TurboGridのO-gridと主流H-gridの接合面でメッシュの非直交性が高いと、偽の渦が発生する
2. 粗すぎるスパン方向メッシュ: 通路渦のサイズに対してセル数が不足すると、渦が分裂や消滅する
3. Mixing Plane面直後: 周方向情報が平均化された後に擬似的な渦構造が再形成される
確認方法: 2水準以上のメッシュで渦構造が定性的に変わらないことを確認する。
損失の評価面
二次流れ損失を評価する断面はどこに設定すべきですか？
翼列出口から下流に翼弦の0.5～1.0倍の距離に設定するのが標準だ。近すぎるとウェイクが完全に発達しておらず、遠すぎると混合損失と二次流れ損失が分離しにくくなる。
損失分離の手法
翼面損失と二次流れ損失を分離できますか？
古典的にはスパン中央の翼面損失を2D損失と見なし、全損失との差を二次流れ損失とする方法がある。CFDではエントロピー生成率を体積積分する方法がより正確だ。
| 損失成分 | 評価方法 |
|---|---|
| 翼面境界層損失 | 翼面近傍のエントロピー生成積分 |
| 端壁損失 | 端壁近傍のエントロピー生成積分 |
| チップ漏れ損失 | チップ隙間領域のエントロピー生成積分 |
| 混合損失 | 翼列出口以降のエントロピー生成積分 |
この分離をやると、どこに改善余地があるか分かりやすいですね。
そう。損失の分離は最適化の方向を決める重要な指針だ。端壁損失が支配的なら端壁コンタリング、チップ漏れが支配的ならスクイーラーチップ、という具合に設計改善の優先順位が決まる。
二次流れCFDの可視化で渦が見えない——Q基準とλ₂基準の選択と閾値
ターボ機械CFDで「渦構造可視化をしたが二次流れの渦が見えない」問題は、閾値設定の誤りが原因であることが多い。渦の可視化にはQ基準（Q>0の領域）またはλ₂基準（λ₂<0の領域）が標準的だが、閾値の選択が難しい。Q基準のQを大きく設定すると強い主渦だけが見え、弱い二次流れ渦が消える。逆に小さすぎると渦でない剪断層も着色されてノイジーになる。推奨手順：①全域のQ値の最大値Qmaxを求め、②0.01×Qmax〜0.1×Qmaxの範囲で閾値を変えて動的に確認する。ParaViewのisovolume機能とアニメーションで閾値をスイープすることで渦の強度分布を把握できる。また渦軸を線表示するVortex Core Filterが定量評価に有用だ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——二次流れの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告