二次流れ
二次流れの理論基礎
概要
二次流れって翼列の中で主流と違う方向に流れるやつですよね？
そう。翼間流路内で主流方向と垂直な速度成分を持つ流れを二次流れと呼ぶ。端壁（ハブ・シュラウド）の境界層が翼列の圧力差で横に曲げられることで発生する。
主な渦構造
どんな渦構造があるんですか？
代表的な渦構造を挙げよう。
| 渦の名称 | 発生メカニズム | 影響 |
|---|---|---|
| 通路渦 (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の古典理論の対応が今も研究テーマになっている。
二次流れの数値計算手法
渦の同定手法
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検証」という段階的精緻化アプローチが多く採用されている。
二次流れの実務適用
端壁コンタリング
二次流れ損失を減らす方法はありますか？
非軸対称端壁コンタリング（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予測はどうやって実験検証するんですか？
代表的な実験手法を挙げよう。
| 手法 | 測定対象 | 特徴 |
|---|---|---|
| 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より若干精度が高いという第三者評価がある。結論として、使用するツールが決まった後はそのツールの乱流モデル実装に合わせた実験ベースの精度検証データを蓄積することが長期的な設計信頼性の基盤になる。
二次流れの先端研究
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が定量評価に有用だ。
関連トピック
なった
詳しく
報告