後流（ウェイク）
後流（ウェイク）の理論基礎
概要
先生、後流って物体の後ろにできる流れのことですよね？
その通り。物体が流体中を運動する（あるいは流れの中に物体がある）とき、物体の下流側に形成される速度欠損領域を後流（wake）と呼ぶ。航空機の後流乱気流、風力発電のウェイク干渉、自動車の空気抵抗、スポーツ選手のスリップストリームなど、応用範囲は極めて広い。
後流の構造
物体から十分離れた遠方後流を考えよう。後流の速度プロファイルは、一様流速度 $U_\infty$ から速度欠損 $u_{def}(x, y)$ を引いた形になる。
遠方後流では $u_{def} \ll U_\infty$ と仮定できるので、線形化された境界層方程式が適用できる。
自己相似解
後流にも自己相似解があるんですか？
遠方後流の自己相似性は Townsend の理論で確立されている。
2D後流（円柱などの2D物体）:
軸対称後流（球などの3D物体）:
噴流は $u_c \propto x^{-1}$ でしたよね。後流のほうが減衰が遅いんですね。
そうだ。後流は一様流の中の「穴」なので、噴流のような強い自己誘導がなく、拡散が遅い。
運動量積分と抗力
後流から物体の抗力が分かるって聞いたんですが。
これは非常に重要な関係だ。後流の運動量欠損を積分すると物体の抗力が得られる。
2D物体の場合（単位スパン長さあたり）:
軸対称物体の場合:
この関係はJones の公式とも呼ばれ、風洞実験で物体の抗力を非接触で測定する方法の基礎だ。後流のピトー管トラバースから速度分布を求め、運動量積分で抗力を算出する。
力を直接測らなくても抗力が分かるんですね。
そうだ。運動量保存則の帰結であり、CFDでもコントロールボリュームの運動量収支から抗力を計算する手法として使われる。壁面の圧力・摩擦力の積分と一致するはずだ。
後流の安定性
後流の安定性解析も重要ですか？
後流プロファイルの安定性解析から、カルマン渦列の特性を説明できる。後流速度プロファイルに対する時間的安定性解析で、
- 反対称モード（sinuous mode）: 渦列の蛇行運動。カルマン渦に対応
- 対称モード（varicose mode）: 後流幅の脈動。通常は反対称モードより不安定性が弱い
Monkewitz (1988) は、後流が「絶対不安定」になる条件（速度欠損が十分大きいとき）を示した。絶対不安定な後流は自励振動を起こし、上流からの擾乱なしでもカルマン渦列が自発的に形成される。
トラック隊列走行と後流の省エネ——スリップストリームの計算
大型トラックが高速道路で隊列を組んで走ると、後続車は前車の後流（ウェイク）の中に入り、空気抵抗が20〜30%減少します。これが「スリップストリーム」効果で、燃費改善の観点から自動運転隊列走行として実用化が進んでいます。後流の速度欠損がどこまで回復するか（後流回復長さ）はRe数と物体の形状に依存し、トラック同士の最適車間距離の設計に直結します。CFDで「前車後流の速度分布を後続車の入口境界条件として使う」連成解析が採用されており、2台、3台と連結が増えるほど後続車の省エネ効果が増すことが計算で確認されています。
後流（ウェイク）の数値計算手法
数値手法の選択
後流のCFDにはどんな手法を使うんですか？
後流の解析は、物体自体の解析と後流域の解析の2つの側面がある。
| 目的 | 手法 | 備考 |
|---|---|---|
| 物体近傍の剥離と近傍後流 | RANS / DES / LES | 壁面解像が必要 |
| 遠方後流の拡散と回復 | RANS / LES | 広い計算領域が必要 |
| 後流の安定性解析 | DNS + Floquet / BiGlobal | 基底状態の精密計算が前提 |
| 風力発電のウェイク干渉 | Actuator Line/Disk + LES | 風車をモデル化して後流に集中 |
後流域のメッシュ設計
後流域のメッシュで注意することは何ですか？
後流は下流に行くほど拡がるので、メッシュもそれに追従させる必要がある。
- 物体直後: 最も細かいメッシュ。再循環領域の長さ程度まで物体表面メッシュと同等の解像度
- 中間後流（$5D\text{--}20D$）: 渦構造が崩壊していく過程。流れ方向にメッシュを緩やかに粗くする（成長率 $< 1.1$）
- 遠方後流（$> 20D$）: 自己相似領域。後流幅に対して少なくとも10セル以上を配置
- 横方向: 後流幅の3倍以上の領域を確保
数値拡散で後流が早く消えてしまう問題はどう対処しますか？
後流の速度欠損は下流で非常に小さくなるため、数値拡散の影響を受けやすい。対策は、
1. 高次精度スキーム: 最低2次精度。LESなら中心差分系
2. メッシュの等方性: 流れ方向にセルを引き延ばしすぎない。アスペクト比 $< 5$
3. 十分な解像度: 速度欠損が $1\%$ 以下の領域でも、欠損プロファイルを分解できるメッシュ
4. AMR (Adaptive Mesh Refinement): 渦度や速度勾配に基づいて動的にメッシュを細かくする
運動量積分法による抗力計算
CFDで後流から抗力を計算する方法を教えてください。
物体から十分下流の断面（例えば $10D$ 下流）で速度分布を取得し、運動量積分を実行する。
2つ目の圧力項は遠方では小さいが、物体近傍の断面では無視できない。この方法で得た抗力と、壁面の圧力・摩擦力の直接積分で得た抗力が一致することを確認するのが、CFD検証の良いプラクティスだ。
OpenFOAM での後流解析
OpenFOAM で後流の統計量を取得するにはどうすればいいですか？
fieldAverage function object を使って時間平均場を計算する。
```
functions
{
fieldAverage1
{
type fieldAverage;
libs ("libfieldFunctionObjects.so");
writeControl writeTime;
fields
(
U { mean on; prime2Mean on; base time; }
p { mean on; prime2Mean on; base time; }
);
}
}
```
これで UMean（時間平均速度）と UPrime2Mean（Reynolds応力テンソル）が出力される。後流の速度欠損プロファイルは UMean から $U_\infty$ を引いて得られる。
後流の「自己相似性」——どこまで行っても同じ形をした流れ
後流（ウェイク）の面白い性質の一つが「自己相似性」です。物体から十分離れた遠後流域では、速度欠損プロファイルがガウス分布に近い形に収束し、スケールを変えると前の断面と同じ形になります。この自己相似解はGaussian far-wake solutionとして理論的に求まり、測定位置を変えても相似パラメータでスケーリングすれば一本の曲線に乗ります。数値計算で後流を正確に追跡するには「ドメインが十分長い（物体後方に直径の50〜100倍）」ことが必要です。後流を数本の断面でしか評価しないと、自己相似解への収束が確認できないまま解析が終わります。
後流（ウェイク）の実務適用
実務での後流解析
実務で後流解析が必要になるのはどんな場面ですか？
主な応用場面を挙げよう。
| 応用 | 関心のある量 | Re範囲 | 典型的な手法 |
|---|---|---|---|
| 航空機の後流乱気流 | 翼端渦の強度と持続距離 | $10^7$ | RANS + 渦追跡 |
| 風力発電のウェイク干渉 | 速度欠損率、乱流付加 | $10^7$ | Actuator Line + LES |
| 自動車の空力 | 後流域の圧力回復、$C_D$ | $10^6$ | RANS (SST) / DES |
| 建築物周りの風環境 | 歩行者レベルの風速 | $10^6\text{--}10^7$ | RANS (k-ε) / LES |
| 海洋構造物（VIV） | 渦放出周波数、横力 | $10^4\text{--}10^6$ | URANS / DES |
風力発電ウェイクの解析
風力発電のウェイク干渉って大きな問題なんですか？
非常に大きな問題だ。上流の風車のウェイクにより、下流の風車が受ける風速が $20\text{--}40\%$ 低下する。出力はほぼ $U^3$ に比例するので、速度の低下は発電量に3乗で効く。
ウェイク解析の手法は3段階ある。
1. 工学モデル: Jensen/Frandsen モデル。ウェイクの拡がりを線形仮定。ファームレイアウト最適化に使用
2. Actuator Disk/Line + RANS: 風車をソース項でモデル化。個々の風車のウェイク予測
3. Actuator Line + LES: 翼素をラインソースで表現し、後流の乱流構造まで解像。最も正確だが計算コスト大
Jensen モデルはどんな式ですか？
最もシンプルなウェイクモデルだ。
ここで $C_T$ はスラスト係数、$k_{wake}$ はウェイク拡がり率（陸上で $\approx 0.075$、洋上で $\approx 0.04$）、$D$ はロータ直径だ。
後流の回復距離
後流はどのくらい下流で回復するんですか？
物体の種類と流れの条件による。目安はこうだ。
| 物体 | 後流が $95\%$ 回復する距離 | 備考 |
|---|---|---|
| 円柱 | $50\text{--}100D$ | 亜臨界域 |
| 球 | $30\text{--}50D$ | 3Dなので回復が速い |
| 風力タービン | $10\text{--}20D$ | 大気乱流による混合で加速 |
| 航空機（翼端渦） | 数km〜十数km | 非常に持続性が高い |
航空機の後流はそんなに長いんですか。
翼端渦は渦核の径が小さく渦度が集中しているため、粘性散逸に非常に長い時間がかかる。大型機の後流は後続機にとって深刻な危険となるため、管制間隔の設定に直結する。Crow 不安定性による翼端渦の蛇行と崩壊の予測が活発に研究されている。
レイノルズの実験（1883年）——乱流発見の瞬間
オズボーン・レイノルズは、管内の水にインクを流す実験で「層流から乱流への遷移」を発見しました。流速を上げていくと、インクの線がある瞬間にグチャグチャに乱れる。この劇的な瞬間を、レイノルズは数学的に $Re = \rho uD/\mu$ という無次元数で表現した。100年以上経った今も、CFDエンジニアが最初に確認するのはこのレイノルズ数です。
後流（ウェイク）のソフトウェア比較
ツール別の特徴
後流解析に適したCFDツールはどれですか？
応用分野に応じたツール選択を整理しよう。
| 応用 | 推奨ツール | 理由 |
|---|---|---|
| 自動車空力 | Fluent, STAR-CCM+, PowerFLOW | DES/LES対応、自動メッシュ |
| 風力発電ウェイク | SOWFA (OpenFOAM), EllipSys3D | Actuator Line実装済み |
| 航空機後流 | OVERFLOW, FUN3D, elsA | NASA/航空宇宙向け高精度コード |
| 汎用後流解析 | OpenFOAM, Fluent | 広い適用範囲 |
| 建築風工学 | WindSim, Fluent, OpenFOAM | 大気境界層のモデリング |
STAR-CCM+ での自動車後流解析
自動車の後流解析はSTAR-CCM+が多いと聞きます。
STAR-CCM+ は自動車CFDでの実績が豊富だ。特徴は、
- 自動メッシュ: ポリヘドラルメッシュ + プリズム層の自動生成。複雑な車体形状に対応
- IDDES: Improved Delayed DES。後流の大規模渦構造を直接解像
- Surface Wrapper: 汚いCADデータのクリーンアップ。実務で重要な機能
- Design Manager: 形状パラメータの自動最適化
自動車の $C_D$ 予測では、後流域の圧力回復が鍵だ。DESを使うと、RANSでは再現できない後流の非定常構造（rear pillar vortex, C-pillar vortex等）を捉えることができ、$C_D$ の予測精度が向上する。
OpenFOAM (SOWFA) での風力発電ウェイク
SOWFA って何ですか？
Simulator fOr Wind Farm Applications の略で、NREL（米国国立再生可能エネルギー研究所）が開発した OpenFOAM ベースのツールだ。大気境界層のLESにActuator Line モデルを組み込んで、ウィンドファーム全体のウェイク干渉を解析できる。
主な特徴は、
- 大気安定度（安定/中立/不安定）の模擬
- コリオリ力と地表面粗度の効果
- Actuator Line/Disk モデルのロータ表現
- 先行風車のウェイクが後続風車に与える影響の定量評価
Fluent での後流解析設定
Fluent で物体後流を解析する場合の推奨設定を教えてください。
自動車のAhmed body（Re $\sim 10^6$）を例にすると、
- メッシュ: Poly-hexcore。後流域に Body of Influence で追加のリファインメント
- 乱流モデル: SST $k$-$\omega$ DDES（後流の非定常構造を解像）
- 空間離散化: Bounded Central Differencing（LES領域）、Second Order Upwind（RANS領域）
- 時間刻み: $\Delta t \approx 5 \times 10^{-4} L/U_\infty$
- 統計収集: $100 L/U_\infty$ 以上の統計平均時間
Ahmed body って後端のスラント角度で後流構造がガラッと変わるんですよね。
その通り。スラント角 $\phi = 25°$ で後流構造がC-pillar渦支配からフルスカッシュバック型に急変し、$C_D$ が不連続的に変化する。これはCFDベンチマークとして非常に有名だ。
後流CFDツール比較——SimScaleのクラウドCFDは本格後流解析に使えるか
後流（Wake）解析ツールの選択肢は近年急速に広がった。SimScaleはブラウザベースのクラウドCFDで、OpenFOAMバックエンドを使いながら専門的なHPC環境なしにLES解析が実行できる。自動車後流の簡易解析やド・ラフト走行の概念検討レベルでは実用的だが、高精度な後流回復距離の予測には解像度が不足するケースがある。一方NUMECA FINE/Openはオープンソースとの連携が容易で、風力後流のSALES（Semi-Analytical LES）モデルを内蔵する。商用ツールのリファレンスとしてはOpenFOAM + windFarmerの組み合わせが研究機関での標準的な選択だ。ツール選定は計算目的の精度要件とライセンスコストのバランスで決まる。
後流（ウェイク）の先端研究
後流の絶対/対流不安定性
後流の安定性理論について詳しく教えてください。
Huerre & Monkewitz (1985, 1990) のフレームワークが標準だ。
- 対流不安定性: 擾乱が下流に流されながら成長。上流に情報が伝わらない。外部擾乱のアンプリファイア
- 絶対不安定性: 擾乱が発生点に留まりつつ成長。自励振動を引き起こす。カルマン渦列の発生源
後流プロファイルのパラメータ $\Lambda = (U_{min} - U_\infty)/(U_{min} + U_\infty)$ が臨界値 $\Lambda_{cr} \approx -0.9$ より小さいとき（速度欠損が十分大きいとき）、後流は絶対不安定になる。
円柱の直後は絶対不安定で、十分下流では対流不安定になるわけですね。
その通り。絶対不安定な領域が「グローバル振動子（global oscillator）」として機能し、カルマン渦列の周波数を決定する。これが Strouhal 数の物理的な起源だ。
後流の渦追跡手法
航空機の翼端渦のような長距離後流はどう計算するんですか？
フルNavier-Stokes計算では計算領域が膨大になるため、ハイブリッド手法が使われる。
- 渦粒子法（Vortex Particle Method）: 渦度場を離散渦粒子で表現。格子なしで渦の移流と拡散を追跡
- Temporal approach: 航空機の直後の断面を初期条件として、時間発展で渦の減衰を追跡（空間座標を時間に変換）
- RANS/LES → 渦粒子法の引き継ぎ: 航空機近傍はCFD、遠方場は渦粒子法で効率的に計算
Proper Orthogonal Decomposition (POD)
後流の解析でPODが使われると聞きました。
PODはスナップショット（瞬時流れ場の集合）からエネルギー的に最適なモード分解を行う手法だ。後流解析では、
第1モードが渦放出の対称/反対称パターンに対応し、寄与率（エネルギー比）が $50\%$ 以上を占めることが多い。低次元モデル（ROM）の構築や流れ制御の設計に活用される。
風力発電のウェイク最新研究
風力発電のウェイク研究の最前線はどうなっていますか？
以下のトピックが活発だ。
- ウェイクステアリング: 上流の風車のヨー角を故意にずらし、ウェイクを横に逸らして下流風車への影響を軽減。Howland et al. (2019, Nature Energy)
- Dynamic Induction Control: ロータの推力を時間変動させて後流の混合を促進
- LiDAR計測との統合: 前方LiDARでウェイクをリアルタイム検出し、制御にフィードバック
- 機械学習によるウェイクモデル: GAN やオートエンコーダで高速ウェイク予測
ウェイクステアリングは実際のウィンドファームで使われているんですか？
Siemens Gamesa や Vestas が実機テストを行い、ファーム全体の発電量が $1\text{--}3\%$ 向上する結果が報告されている。一見小さいが、大規模洋上風力で年間数十億円の差になりうる。
風力発電機後流の重合効果——オフショア風力ファームのCFD最適化
洋上風力発電所（オフショアウィンドファーム）では、上流タービンの後流（Wake）が下流タービンに当たることで発電量が15〜40%低下する「Wake Effect」が大きな問題だ。後流の回復距離は大気乱流強度と風車間距離（D: ロータ径）に依存し、通常8D〜10D程度を要する。RANS-CFDによるファームレベル解析は数十台のタービンを同時に扱えるが、個々の後流構造は近似的にしか捉えられない。最前線ではLES解析によるタービン間の複雑な渦干渉の解明と、機械学習（GP回帰）を組み合わせたファーム最適化（ヨー制御、ティルティング）が実用化段階に入っており、配置最適化だけで年間発電量を5〜10%増加させる実績が報告されている。
後流（ウェイク）のトラブル対応
よくあるトラブル
後流解析のCFDでよくある問題を教えてください。
1. 後流が早く消える（速度欠損の過小評価）
原因と対策:
| 原因 | 対策 |
|---|---|
| 数値拡散（1次風上スキーム） | 2次精度以上。LES部分は中心差分系 |
| 後流域のメッシュが粗い | 後流幅に対して最低10セル。Body of Influence を活用 |
| メッシュのアスペクト比が大きい | 後流域のアスペクト比を5以下に。等方的なメッシュが理想 |
| RANSの渦粘性が過大 | DES/DDESに切り替え。後流のLES領域で渦粘性を下げる |
2. 後流が非対称（物理的に対称であるべき場合）
対称物体なのに後流が偏るんですが。
確認ポイント:
- メッシュの対称性: 非構造メッシュの場合、対称面でメッシュが非対称になることがある。対称面に合わせてメッシュを生成する
- 非定常性: 実は物理的に非対称かもしれない。カルマン渦列は瞬時には非対称だ。時間平均を取ると対称になるはず
- ヒステリシス: 高Re数では複数の安定状態が共存する場合がある。Ahmed bodyの $\phi = 25°$ 近傍など
3. 壁面からの剥離位置がずれる
後流の構造は剥離位置に支配されるため、壁面の剥離位置を正しく予測することが最も重要だ。
対策:
- $y^+$ の確認: 壁関数使用時は $30 < y^+ < 300$、壁面解像時は $y^+ < 1$
- 遷移モデル: 層流境界層からの遷移位置が重要な場合（球の抗力危機など）、$\gamma$-$Re_\theta$ 遷移モデルを使用
- Grid Induced Separation (GIS): DDESを使ってRANS領域の境界層をLESに切り替える際のメッシュ起因の非物理的剥離を防止
4. DES/DDESでの灰色領域問題
DES って万能じゃないんですか？
DES の最大の問題は「灰色領域」だ。RANS領域とLES領域の遷移領域で、どちらのモデルも十分に機能しない。
症状: 剥離直後の渦が解像されず、本来あるべき乱流構造が発達しない。Modelled Stress Depletion (MSD) とも呼ばれる。
対策:
- DDES（Delayed DES）を使用して境界層内のLES化を防ぐ
- IDDES（Improved DDES）でさらに改善
- シールド関数のキャリブレーションを確認
- メッシュの粗密遷移を滑らかにする（急激なメッシュ変化はNG）
5. 運動量積分と壁面力積分の不一致
後流の運動量積分から求めた抗力と、壁面の圧力・摩擦力の積分から求めた抗力が一致しない場合、
- 運動量積分の断面位置が不適切: 物体から離れすぎると後流が計算領域外に出る。近すぎると圧力項が大きく不確実
- 出口境界条件の影響: 出口が近すぎると運動量が保存されない
- メッシュ収束していない: 両方の値がメッシュを細かくしても一致に向かわない場合は、より根本的な問題がある
両者の一致は解の信頼性の良いチェックですね。
その通りだ。壁面積分と後流の運動量積分の一致は「conservation check」として非常に有用だ。差が $5\%$ 以内であれば許容範囲、$1\%$ 以内なら非常に良い結果だ。
後流CFDの収束不良——非定常渦放出をRANSで定常解析する矛盾
円柱後流のCFD解析で「残差が収束しない」という問題は、実はCFDの問題ではなく物理の問題であることが多い。Re>46では円柱後流にカルマン渦列が発生し、流れ場は本質的に非定常だ。これを定常RANS（Steady-State RANS）で解こうとすると、ソルバは安定解を見つけられず残差が振動し続ける。対策は2つ：①非定常RANS（URANS）に切り替え時間平均量として評価する、②SAS（Scale-Adaptive Simulation）やDESを用いて大規模渦を直接計算する。実務では「残差が振動 = 物理的に非定常 = 定常解析の仮定が成立しない」のサインと認識し、解析手法の根本から見直すことが肝要だ。
なった
詳しく
報告