渦放出 — CAE用語解説
渦放出
渦放出とは何か
渦放出って流れの中に物体を置くと渦が交互に出てくる現象ですよね？
理論と物理
渦放出の基本概念
「渦放出」って、教科書で見たんですけど、具体的にどんな現象なんですか？単に渦ができることとは違うんですか？
良い質問だ。渦放出（Vortex Shedding）は、流体が非流線形物体（円柱など）の背後で、周期的に交互に渦を生成・放出する現象だ。単なる渦形成と決定的に違うのはその「周期性」と「規則性」にある。例えば直径
交互に、というのはどういうパターンなんですか？なぜ規則正しく起こるんですか？
物体の上下（または左右）で圧力分布が不安定になり、交互に渦が剥離する。これが「カルマン渦列」を形成する。規則性が生まれるメカニズムは、剥離点の移動と渦の成長・放出のフィードバックにある。レイノルズ数が約300から3×10^5の範囲で、この周期的な剥離が顕著に観測される。煙突や海底パイプラインの振動問題の根源だ。
支配方程式はナビエ-ストークス方程式ですよね？渦放出を予測するために、特別な項や条件は必要ないんですか？
その通り、非圧縮性ナビエ-ストークス方程式が基礎だ。特別な項は不要だが、非定常計算が絶対条件だ。時間発展項
数値解法と実装
CFDでの離散化とソルバー設定
CFDで渦放出をシミュレーションする時、一番気をつけるべき離散化のポイントは何ですか？
主に3点だ。第一に「時間離散」。二次精度以上の陰的または陽的スキームが必須だ。例えば、Ansys Fluentではbounded second order implicitがよく使われる。第二に「空間離散」。移流項には二次風上差分（QUICKやSecond Order Upwind）が必要で、一次風上では数値拡散で渦が消えてしまう。第三に「時間刻み幅
メッシュはどのくらい細かくすればいいんでしょうか？境界層の解像も必要ですか？
渦放出の主なメカニズムは物体背後の中～遠方の剥離領域にある。だから、壁面解像RANS（Reynolds-Averaged Navier-Stokes）を使う場合、
乱流モデルはどう選べばいいですか？LES（Large Eddy Simulation）じゃないとダメですか？
目的による。ストローハル数などの大まかな周波数を知りたいだけなら、非定常RANS（URANS）のSST k-ωモデルでも可能だ。しかし、渦の詳細な3次元構造や遷移現象を捉えるにはLES、または分離渦シミュレーション（DES）が望ましい。ただし、LESは計算コストが桁違いに高い。実務では、ISO 13715などの風荷重評価ではURANS結果をベースに安全率を掛けることが多い。
実践ガイド
解析ワークフローと検証
実際に円柱の渦放出解析を始めるとしたら、最初に何から手を付けるべきですか？
まずは「2次元、低レイノルズ数」のベンチマークケースからだ。具体的には、直径1m、流速1 m/s、動粘性係数1×10^-6 m²/sとしてレイノルズ数Re=10^6とするのではなく、Re=100や1000といった層流域から始める。この方が渦のパターンが明確で、計算も軽い。メッシュ依存性調査を必ず行い、リフト係数（C_L）の振幅と周波数がメッシュを細かくしても収束することを確認する。
結果が合っているかどうか、どうやって検証すればいいですか？実験データが無い場合が多いですよね。
実験データがなくても、文献値と比較できる。例えば、円柱の2次元流れでRe=100なら、ストローハル数Stは約0.164、時間平均抗力係数C_Dは約1.35±0.05、リフト係数C_Lの振幅は約0.33という確立された数値実験結果がある。自分のシミュレーションでこれらの値（特にSt）が±5%以内に収まっているか確認する。また、後流の渦パターンが「カルマン渦列」の可視化画像と定性的に一致するかも重要だ。
解析時間はどのくらい見積もればいいですか？定常状態になるまで計算する必要がありますか？
少なくとも「初期過渡現象が消え、周期的定常状態に達した後」で、さらに10周期分以上のデータを取得する必要がある。例えば、ストローハル数から周期T=1/fを求め、初期過渡に5T、データ取得に10T分の物理時間を計算する計画を立てる。時間刻み幅が0.01Tなら、合計で1500ステップ以上だ。モニタポイントでのリフト係数の時系列を確認し、振幅と周期が一定になるまで計算を続ける。
ソフトウェア比較
各ソルバーの特徴と適用例
Ansys FluentとOpenFOAMで渦放出解析をやる場合、大きな違いは何ですか？
根本的な違いは「ソルバーアルゴリズムの選択肢と初期設定」だ。Fluentは圧力基盤のソルバーがデフォルトで、SIMPLE/PISO系列のアルゴリズムが強力に最適化されている。非定常計算の設定もGUIから直感的。一方、OpenFOAMはコードベースなので、例えば`pimpleFoam`（PIMPLEアルゴリズム）や`icoFoam`（非圧縮性層流）などソルバーを選択し、`fvSchemes`ファイルで数値スキームを全て自分で定義する必要がある。柔軟性は高いが、設定ミスのリスクも高い。
Abaqus/CFDやCOMSOL Multiphysicsのようなマルチフィジックスソフトではどうですか？
それらは「流体-構造連成（FSI）」解析に強い。単なる渦放出の可視化なら専用CFDソフトで十分だが、渦放出による「流れ誘起振動（VIV）」を解析する場合はこれらのツールが本領を発揮する。COMSOLでは「Laminar Flow」物理場と「Solid Mechanics」物理場を「Two-Way Coupling」で結合し、メッシュ変形を考慮できる。Abaqus/CFDも同様に、構造解析で培われた非線形接触や大変形解析と流体を連成させられる。ただし、純粋なCFDの機能面や計算速度ではFluentやSTAR-CCM+に一歩譲る部分もある。
STAR-CCM+の特徴は何ですか？
STAR-CCM+は「オールインワン」のワークフローと、非常に強力な「メッシング」機能が売りだ。特にポリヘドラルメッシュは、複雑な形状でも比較的少ないセル数で高品質なメッシュを生成できる。渦放出のような剥離流れでは、後流領域に自動で適応メッシュ細分化（AMR）を設定するのが簡単で、ユーザの手間を大幅に減らせる。また、シミュレーションの運用手法（Design Manager）が整っており、パラメータスタディや最適化と組み合わせやすい。自動車のアウターパネルや航空機の脚周りの流れなど、産業界で広く使われている。
トラブルシューティング
よくあるエラーと対策
非定常計算を始めたのに、いつまで経っても周期的な渦が出てきません。リフト係数の時系列がノイズのようで、きれいな正弦波になりません。なぜですか？
最も多い原因は3つだ。1. **時間刻み幅が大きすぎる**：1周期を10ステップ以下で解像しようとしている。周期Tを予測して、
計算が発散してしまいます。特に、PISOソルバーを使っている時に圧力の修正で発散することが多いです。
これは古典的な問題だ。非定常剥離流れは圧力と速度の強い連成を持つ。対策は、1. **アンダーリラクセーション係数を下げる**：圧力と運動量の緩和係数をデフォルト（0.3, 0.7など）からより小さい値（0.1, 0.5）に下げて安定性を優先する。2. **非定常計算の最初は時間刻み幅を極端に小さくする**：最初の数ステップは
文献値とストローハル数が10%以上ずれています。メッシュと時間刻みは十分細かくしたつもりです。
メッシュと時間刻み以外の要因を疑うべきだ。第一に**計算領域のサイズ**。特に流出境界が物体に近すぎないか？円柱の場合、流出境界は少なくとも下流20D以上は離す。側方境界も15D以上離し、ブロッケージ効果を避ける。第二に**乱流モデルの影響**。URANSのSST k-ωモデルは、渦粘性が高すぎて渦の減衰が早まり、周波数がずれることがある。第三に**2次元近似の限界**。現実の渦放出は3次元的で、特に高レイノルズ数ではスパン方向に変動がある。2次元解析は本質的に近似なので、ずれが生じることは覚悟する必要がある。
計算コストが高すぎて、実用的なサイズのモデルで非定常計算ができません。何か対策は？
産業設計では完全な非定常計算を避ける方法もある。1. **スペクトル解析の利用**：定常RANS計算で得た流れ場を基に、線形安定性解析や固有値解析を行い、支配的な周波数を推定する（いわゆる「渦放出の線形理論」）。2. **経験式・規格の利用**：建築分野ではAIJ（日本建築学会）推薦の風荷重規定や、ISO 13715などに基づき、ストローハル数と設計風速から強制力の周波数を算定する。3. **ハイブリッドRANS/LESの利用**：全域をLESするのではなく、剥離予想領域だけをLES解像度とするDES（Detached Eddy Simulation）やSAS（Scale Adaptive Simulation）モデルを使う。これらはRANSとLESの中間的コストで、後流の剥離渦を捉えられる。
関連トピック
なった
詳しく
報告