渦励振（VIV）解析
理論と物理
渦励振現象の物理的背景
海洋ライザーや煙突が風で揺れるのは単なる風荷重じゃないって聞いたんですが、どういうメカニズムなんですか？
円柱状の構造物に流れが当たると、後流にカルマン渦列が交互に放出される。この渦放出に伴って構造物に周期的な揚力変動が作用し、構造物が振動する。これが渦励振（Vortex-Induced Vibration, VIV）だよ。
渦放出周波数 $f_s$ はStrouhal数 $St$ で特徴づけられる。
ここで $D$ は円柱直径、$U$ は一様流速だ。Re数が $300 < Re < 3 \times 10^5$ の範囲では $St \approx 0.2$ がよく使われる。
渦の放出周波数が構造物の固有振動数に近づくとまずいんですよね？
その通り。$f_s$ が構造物の固有振動数 $f_n$ に近づくとロックイン（lock-in）現象が発生する。渦放出周波数が構造物の振動数に同期し、振幅が急激に増大するんだ。
この領域では換算速度 $U^ = U/(f_n D)$ がおよそ $4 < U^ < 8$ の範囲に対応する。
支配方程式
流体と構造が相互作用するので、両方の方程式を連立する必要があるんですか？
流体側はNavier-Stokes方程式を解く。非圧縮性を仮定すると、
構造側は運動方程式で表される。単純化した1自由度モデルだと、
ここで $F_L(t)$ は流体から受ける揚力だ。質量比 $m^* = m / (\rho_f D^2 L)$ と減衰比 $\zeta$ が応答を支配する。
この二つの方程式が界面条件でつながるんですね。
流体-構造界面で速度の適合条件と力の平衡条件を課す。
ALE（Arbitrary Lagrangian-Eulerian）法を使って移動メッシュを扱うのが標準的なアプローチだ。
カルマン渦の「規則正しさ」——乱流の中の秩序
流体が円柱を過ぎるとき、Re数40〜200程度の層流域で現れるカルマン渦列は自然界の「秩序ある乱れ」の代表例だ。渦の放出周波数はStrouhal数 St = fD/U ≈ 0.2 で驚くほど安定している。これは円柱直径が1mmでも1mでも、流速を正規化すれば同じ無次元数になるという「相似則」の美しさだ。煙突・海洋ライザー管・橋梁のハンガーケーブル、さらには人体の血管まで、円柱状構造物があればどこでもカルマン渦が現れる。工学的には「厄介もの」として扱われる渦励振だが、渦の規則性を利用した「渦流量計」（流量を渦の周波数で計測）として、化学・石油プラントで広く使われている。
各項の物理的意味
- 構造-熱連成項：温度変化による熱膨張が構造変形を誘発し、変形が温度場に影響する。$\sigma = D(\varepsilon - \alpha \Delta T)$。【日常の例】夏に線路のレールが伸びて隙間が狭くなる——温度上昇→熱膨張→応力発生の典型例。電子基板がはんだ付け後に反るのも、異なる材料の熱膨張率差による。エンジンのシリンダーブロックは高温部と低温部の温度差で熱応力が発生し、最悪の場合亀裂に至る。
- 流体-構造連成（FSI）項：流体圧力・せん断力が構造を変形させ、構造変形が流体領域を変化させる双方向の相互作用。【日常の例】強風で吊り橋のケーブルが振動する（渦励振）——風の力が構造を揺らし、揺れた構造が風の流れを変え、さらに振動が増幅する。心臓の血流と血管壁の弾性変形、航空機の翼のフラッタ（空力弾性不安定性）も典型的なFSI問題。片方向のみの連成で済む場合もあるが、変形が大きい場合は双方向連成が必須。
- 電磁-熱連成項：ジュール発熱 $Q = J^2/\sigma$ が温度上昇を引き起こし、温度変化が電気抵抗を変化させるフィードバックループ。【日常の例】電気ストーブのニクロム線は電流が流れると発熱（ジュール熱）して赤くなる——温度が上がると抵抗が変わり、電流分布も変化する。IHクッキングヒーターの渦電流発熱、送電線の温度上昇による弛み増加もこの連成の例。
- データ転写項：異なる物理場間のメッシュ不一致を補間で解決。【日常の例】天気予報で「気温のデータ」と「風のデータ」を合わせて体感温度を計算するとき、それぞれの観測地点が異なれば補間が必要——CAEの連成解析でも、構造メッシュとCFDメッシュは一般に一致しないため、界面でのデータ転写（補間）精度が結果の信頼性に直結する。
仮定条件と適用限界
- 弱連成仮定（片方向連成）：一方の物理場が他方に影響するが逆は無視可能な場合に有効
- 強連成が必要なケース：FSIでの大変形、電磁-熱連成での温度依存性が強い場合
- 時間スケールの分離：各物理場の特性時間が大きく異なる場合、サブサイクリングで効率化可能
- 界面条件の整合性：連成界面でのエネルギー・運動量保存が数値的に満たされることを確認
- 適用外ケース：3つ以上の物理場が同時に強く連成する場合、モノリシック手法が必要になることがある
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 熱膨張係数 $\alpha$ | 1/K | 鋼: 約12×10⁻⁶、アルミ: 約23×10⁻⁶ |
| 連成界面力 | N/m²（圧力）またはN（集中力） | 流体側と構造側で力の釣り合いを確認 |
| データ転写誤差 | 無次元（%） | 補間精度はメッシュ密度比に依存。5%以下が目安 |
数値解法と実装
ALE法によるメッシュ移動
構造物が動くと流体メッシュも変形させないといけないですよね。どうやって処理するんですか？
ALE法ではメッシュ速度 $\mathbf{w}$ を導入して、Navier-Stokes方程式の移流項を修正する。
メッシュの移動はラプラス方程式やバネアナロジーで決定する。大変形時にはメッシュリメッシングが必要になるよ。
リメッシングのタイミングってどう決めるんですか？
要素品質指標（アスペクト比やスキューネス）をモニタリングして、閾値を下回ったら自動リメッシングを実行する。Ansys FluentではDynamic Mesh機能で設定できる。STAR-CCM+のモーフィングメッシュも同様の仕組みだ。
連成アルゴリズム
流体と構造をどういう順番で解くかって重要なんですか？
非常に重要だ。大きく分けて2種類ある。
弱連成（Loose/Weak coupling）: 各タイムステップで流体→構造を1回ずつ解く。計算は速いが、質量比が小さい（$m^* < 5$程度の）系では不安定になりやすい。
強連成（Strong coupling）: 各タイムステップ内でサブイテレーションを行い、界面条件を収束させる。安定だが計算コストが高い。
強連成では、Aitken緩和やIQN-ILS（Interface Quasi-Newton Inverse Least Squares）法で収束を加速する。
海洋ライザーのように質量比が小さい場合は強連成が必須ということですね。
その通り。水中構造物は $m^* \approx 1$ になることもあるから、弱連成ではadded mass effectによる数値不安定が起きる。
乱流モデルの選択
VIV解析ではRANSとLESのどちらを使うべきですか？
渦放出を正確に捕えるにはLESやDESが望ましい。RANSの $k$-$\omega$ SST モデルでは2D渦放出のタイミングは捉えられるが、3D的な渦構造やスパン方向の相関が表現できない。
| 乱流モデル | 渦放出精度 | 計算コスト | 推奨Re範囲 |
|---|---|---|---|
| URANS (k-omega SST) | 中 | 低 | Re < 10^4 |
| DES/DDES | 高 | 中 | 10^4 < Re < 10^6 |
| LES (Wall-Resolved) | 非常に高 | 非常に高 | Re < 10^5 |
| LES (Wall-Modeled) | 高 | 高 | Re < 10^6 |
「ロック・イン」現象——構造が渦を支配する瞬間
VIV解析で最も重要な概念がロック・イン（Lock-in）だ。通常、渦の放出周波数は流速に比例して変化する。しかし流速が構造の固有振動数に近づくと、渦の周波数が構造に「引きずられて」固有振動数に張り付いてしまう。この状態がロック・インで、流速が変わっても渦は同じ周波数で放出し続け、振動が持続・増大する。ロック・インの速度範囲（ロック・イン幅）は構造の質量比と減衰比によって変わる。軽くて減衰が小さい構造（鋼製スリムな煙突など）は広いロック・イン域を持ち、より危険だ。VIV解析では「ロック・イン域で何次モードがどの流速で共振するか」のマッピングが設計の核心になる。
モノリシック法
全物理場を1つの連立方程式系として同時に解く。強い連成に対して安定だが、実装が複雑でメモリ消費が大きい。
パーティション法（分離反復法）
各物理場を独立に解き、界面でデータ交換。実装が容易で既存ソルバーを活用可能。弱い連成に適する。
界面データ転写
最近傍法（最も簡単だが精度低い）、射影法（保存的）、RBF補間（メッシュ非一致に強い）。保存性と精度のバランスが重要。
サブイタレーション
各連成ステップ内で十分な反復を行い、界面条件の整合性を確保。残差基準は各物理場の典型値に基づいてスケーリング。
Aitken緩和
連成反復の緩和係数を自動調整。過緩和による発散を防止し、収束を加速する適応的手法。
安定性条件
added mass効果（流体-構造連成で構造密度≈流体密度の場合）に注意。不安定な場合はロビン型界面条件やIQN-ILS法を適用。
Aitken緩和のたとえ
Aitken緩和は「シーソーのバランス取り」に似ている。一方が強く押しすぎると反対側が跳ね上がり、その反動でまた強く押しすぎる——この振動を抑えるために、押す力を自動的に調整するのがAitken緩和。連成反復が振動して収束しないとき、前回の修正量を見て次の修正量を自動調整する適応的手法。
実践ガイド
解析モデルの構築手順
実際にVIV解析を始めるとき、どういう手順で進めればいいですか？
まずジオメトリの準備だ。円柱直径 $D$ に対して、流入方向に $20D$ 以上、後流方向に $40D$ 以上、横方向に $20D$ 以上の計算領域を確保する。これがJournal of Fluids and Structuresなどで推奨されている最低限のサイズだ。
メッシュはどのくらい細かくすればいいんですか？
円柱表面の最初のセル高さは壁面 $y^+$ に依存する。URANSなら $y^+ \approx 1$、LESでも $y^+ < 1$ が必要だ。円柱周りに最低200分割以上を確保し、O型メッシュで境界層を解像する。
| メッシュパラメータ | URANS推奨 | LES推奨 |
|---|---|---|
| 円柱周方向分割数 | 200以上 | 360以上 |
| 壁面第1層高さ | $y^+ \approx 1$ | $y^+ < 1$ |
| 境界層内の層数 | 20以上 | 30以上 |
| 後流域解像度 | $D/20$ | $D/40$ |
| スパン方向分割（3D） | - | $\pi D / 20$ 以下 |
時間刻みとCourant数
非定常解析だと時間刻みの設定も重要ですよね。
渦放出1周期あたり最低200タイムステップ以上が必要だ。具体的には、
最大Courant数はURANSで $CFL < 5$、LESでは $CFL < 1$ を目標にする。
定常状態に達するまでどのくらい回す必要がありますか？
初期過渡を除外するために最低20周期、統計量を取得するのにさらに20周期以上流す。lift coefficientの時刻歴をモニタリングして、振幅が安定したことを確認するんだ。
メッシュ収束確認
メッシュの妥当性はどうやって検証するんですか？
3水準以上のメッシュ密度でStrouhal数、平均抗力係数 $\bar{C}_D$、揚力係数の振幅 $C_{L,rms}$ を比較する。Richardson外挿法でグリッド収束指数（GCI）を算出するのが定量的だ。
ここで $F_s = 1.25$（安全係数）、$r$ はメッシュ細分化率、$p$ は収束の次数だよ。