WALEモデル
理論と物理
概要
先生、WALEモデルって産業LESで一番人気って聞いたんですけど、どういうモデルですか？
WALE（Wall-Adapting Local Eddy-viscosity）モデルはNicoud-Ducros (1999) が提案したLES SGSモデルだ。速度勾配テンソルの2乗テンソル $g_{ij}^2 = \bar{g}_{ik}\bar{g}_{kj}$ のトレースレス対称部分 $S_{ij}^d$ を使うことで、壁面近傍で自動的に $\nu_{\text{sgs}} \to 0$ となる。Van Driest減衰関数も壁面距離も不要だ。
支配方程式
具体的な式を教えてください。
SGS渦粘性は次式で計算される。
ここで $S_{ij}^d$ は速度勾配テンソルの2乗の対称トレースレス部分だ。
モデル定数は $C_w = 0.325$（一般的にはSmagorinskyの $C_s = 0.1$ に相当）。
WALEの壁面挙動
なぜ壁面で自動的にゼロになるんですか？
壁面近傍では速度場が $u \sim y$（線形プロファイル）となる。このとき $\bar{S}_{ij} \sim O(1)$ だが $S_{ij}^d \sim O(y)$ となるため、WALEの分子は $O(y^3)$、分母は $O(1)$ となり、$\nu_{\text{sgs}} \sim y^3$ で減衰する。これは壁面近傍の $\nu_t \sim y^3$ という理論的要件と一致する。
比較すると、Smagorinskyモデルでは $\nu_{\text{sgs}} \sim |\bar{S}| \sim O(1)$（壁面で減衰しない）ため、Van Driest減衰 $f = 1 - \exp(-y^+/A^+)$ が必要だった。
| モデル | 壁面での $\nu_{\text{sgs}}$ | 追加処理 |
|---|---|---|
| Smagorinsky | $O(1)$（減衰しない） | Van Driest減衰必須 |
| Dynamic Smagorinsky | $O(y^3)$（自動） | 不要 |
| WALE | $O(y^3)$（自動） | 不要 |
WALEは動的Smagorinskyと同じ壁面挙動を、テストフィルタの計算なしで実現してるんですね。だから産業用途で人気なのか。
WALEモデルが「壁面にもvan Driest不要」を実現した仕組み
従来のSmagorinskyモデルでは壁面付近での過剰散逸を防ぐためにvan Driest減衰関数が必要でした。1999年にNicoud & Ducrosが提案したWALE（Wall-Adapting Local Eddy-viscosity）モデルは、速度勾配テンソルの不変量を巧みに組み合わせることで「壁面近傍で渦粘性が自然に $y^3$ オーダーでゼロに近づく」性質を実現しています。van Driest関数が不要というのは見た目以上に大きな利点で、複雑形状でも壁面距離の計算が不要になるのです。
各項の物理的意味
- 時間項 $\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$。時間刻みの安定性に直結 |
数値解法と実装
実装の詳細
WALEモデルの実装で注意すべき点はありますか？
$S_{ij}^d$ の計算には速度勾配テンソル $\bar{g}_{ij} = \partial \bar{u}_i / \partial x_j$ の9成分と、その2乗 $\bar{g}_{ij}^2 = \bar{g}_{ik}\bar{g}_{kj}$ が必要だ。計算コストはSmagorinskyモデルと大差ない。
| 計算ステップ | 必要な演算 |
|---|---|
| 1. $\bar{g}_{ij}$ の計算 | 速度勾配テンソル（9成分） |
| 2. $\bar{g}_{ij}^2$ の計算 | テンソル積（9成分） |
| 3. $S_{ij}^d$ の計算 | 対称化＋トレース除去（6成分） |
| 4. $\nu_{\text{sgs}}$ の計算 | スカラー演算 |
ソルバー別の設定
各ソルバーでの設定方法を教えてください。
WALEはすべての主要CFDソルバーで標準サポートされている。設定が簡単で、特別なパラメータ調整が不要な点が産業用途で高く評価されている。
グリッド幅 $\Delta$ の定義
WALEモデルでのグリッド幅はどう定義しますか？
一般的には体積等価幅 $\Delta = V^{1/3}$ が使われる。構造格子では $\Delta = (\Delta_x \Delta_y \Delta_z)^{1/3}$ だ。WALEモデルは壁面距離に依存しないため、$\Delta$ の定義に対するロバスト性が高い。
WALEはシンプルな実装でありながら壁面挙動が正しく、定数のチューニングも不要。産業LESのデフォルトSGSモデルとして最適な選択ですね。
WALEの $S_{ij}^d$ テンソル——なぜ「速度勾配の二乗」を使うのか
WALEモデルの肝は $\mathbf{g}^2 = \mathbf{g} \cdot \mathbf{g}$（速度勾配テンソルの2乗）の対称偏差成分 $S_{ij}^d$ です。なぜわざわざ二乗するのか？ 速度勾配テンソルを一乗のまま使うと壁面近傍で渦粘性が $y$ に比例して消えますが、固体壁面の物理では $y^3$ オーダーで消えるのが正しい。速度勾配の二乗を使うことで、この正しいスケーリングを自動的に得られるという仕組みです。シンプルな修正で大きな改善を得た好例として、乱流モデリングのテキストでよく紹介されます。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
適用範囲
WALEモデルはどういう場面で使うべきですか？
| 適用場面 | 理由 |
|---|---|
| 自動車空力LES | 壁面近傍の挙動が正しく、設定が簡単 |
| 建築物周りの風環境解析 | メッシュ非均一性にロバスト |
| 混合・撹拌のLES | 壁面距離不要で複雑形状に対応 |
| DES/DDESのLES領域のSGSモデル | 壁面依存性がないため相性が良い |
メッシュ要件
WALEモデルに必要なメッシュ解像度は？
標準的なWall-Resolved LESの要件と同じだ。
| パラメータ | 推奨値（壁単位） |
|---|---|
| $y^+$ (第1層) | < 1 |
| $\Delta x^+$ (流れ方向) | 20〜50 |
| $\Delta z^+$ (スパン方向) | 10〜20 |
| 境界層内の層数 | 15〜25 |
$C_w = 0.325$ は変える必要がありますか？
通常は変更不要だ。ただし非常に粗いメッシュ（Wall-Modeled LES的な使い方）では $C_w$ を若干大きくすることで数値安定性が改善する場合がある。逆にDNS級のメッシュでは $C_w$ を小さくして過剰散逸を抑えることもある。ただしこれは例外的なケースだ。
他のSGSモデルとの比較
WALEと他のSGSモデルの結果の差はどのくらいですか？
十分に解像されたLES（Grid-Resolved LES）では、SGSモデルの選択は結果にほとんど影響しない。SGSモデルが重要になるのは解像度が限界的な場合だ。
| メッシュ解像度 | SGSモデルの影響 |
|---|---|
| 十分に細かい（$k_{\max}\eta > 1$） | ほぼ無視できる |
| 適度（$k_{\max}\eta \sim 0.5$） | 10%程度の差 |
| 粗い（$k_{\max}\eta < 0.3$） | 大きな差（モデル依存性高い） |
つまり十分なメッシュ解像度があればSGSモデルの選択は二の次で、メッシュの品質が最も重要だと。
WALEモデルが「建物周辺風環境」で支持される理由
都市風環境のLESでWALEモデルが好まれる背景には実務的な理由があります。ビル間の複雑な形状では壁面距離 $y^+$ の計算が煩雑になり、van Driest減衰関数の適用範囲が曖昧になりがち。WALEモデルはその必要がないため「とにかく流して結果を見たい」という実務ニーズにフィットします。環境アセスメント案件では風速比や乱流強度の評価精度より計算スピードが優先されることも多く、適度な精度と低コストのバランスがWALEの強みです。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
ソルバー別の実装特性
各ソルバーのWALE実装に違いはありますか？
モデル式自体は同じだが、速度勾配の計算方法やグリッド幅の定義に微差がある。
| 特性 | Fluent | STAR-CCM+ | OpenFOAM |
|---|---|---|---|
| 勾配計算 | Green-Gauss or LSQ | Hybrid LSQ | Gauss linear |
| $\Delta$ の定義 | $V^{1/3}$ | $V^{1/3}$ | cubeRootVol |
| 数値散逸の低減 | BCD推奨 | Blended推奨 | filteredLinear推奨 |
| DES/DDESとの組合せ | 可能 | 可能 | 可能 |
WALEをDDESのSGSモデルとして使えるんですか？
できる。DDESのLES領域のSGSモデルとしてWALEを指定する構成は一般的だ。Fluentでは DDES + WALE SGS の組合せが選択可能。壁面距離に依存しないWALEはDDESとの相性が特に良い。
LES SGSモデルの選択フロー
LESのSGSモデルを選ぶ判断フローを教えてください。
1. 産業LESで定数チューニングを避けたい → WALE
LESのSGSモデルを選ぶ判断フローを教えてください。
1. 産業LESで定数チューニングを避けたい → WALE
2. 遷移流を含むLESで最高精度が必要 → Dynamic Smagorinsky
3. 純粋せん断流（Couette流等）でゼロ散逸が必要 → $\sigma$ モデル
4. 既存のベンチマーク結果と比較したい → Smagorinsky ($C_s = 0.1$)
5. 簡単な予備検討 → Smagorinsky
迷ったらWALE、精度追求なら動的Smagorinsky、というのが現実的な判断ですね。
WALEモデル定数 $C_w = 0.5$ はどこから来たか
WALEモデルのモデル定数 $C_w$ は多くの実装で0.5が採用されていますが、「なぜ0.5なのか？」と問われると意外と答えが難しい。Nicoud & Ducros（1999）の原論文では等方性乱流のDNSデータとの比較から導出されていますが、Smagorinskyの $C_s$ と同じく万能ではありません。Fluent、OpenFOAM、StarCCM+いずれも0.5をデフォルトにしており「業界標準」になっていますが、剪断流れや回転流れでは調整が必要なケースもある——という認識を持っておくのが大切です。
選定で最も重要な3つの問い
- 「何を解くか」：WALEモデルに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
WALE-ABLモデル
WALEモデルの発展版はありますか？
大気境界層（ABL）への適用に特化したWALE-ABLモデルが提案されている。大気LESでは地表面粗さの効果が重要で、標準WALEに粗さパラメータを追加する修正が行われている。
WALEとAMR（適応格子細分化）
WALEとAMR（Adaptive Mesh Refinement）の組合せは？
WALEは壁面距離に依存しないため、AMRでメッシュが局所的に細分化されても安定に動作する。Smagorinskyモデル+Van Driest減衰はAMRで壁面距離が変化すると問題が生じるが、WALEではそのリスクがない。
非ニュートン流体への拡張
非ニュートン流体のLESにもWALEは使えますか？
Rudman-Blackburn (2006) らが非ニュートン流体（べき乗則流体、ビンガム流体等）へのWALEの拡張を提案している。SGS粘性に実効粘度を組み込む修正が行われる。ポリマー溶液や血液流れのLESで活用されている。
$\sigma$ モデルやVremanモデルとの理論的比較
WALE、$\sigma$ モデル、Vremanモデルはどれも壁面で自動的にゼロになるんですよね。違いは何ですか？
数学的な性質が異なる。
| 性質 | WALE | $\sigma$ | Vreman |
|---|---|---|---|
| 純粋せん断で $\nu_{\text{sgs}} = 0$ | いいえ | はい | はい |
| 純粋回転で $\nu_{\text{sgs}} = 0$ | はい | はい | はい |
| 2Dの軸対称拡大で $\nu_{\text{sgs}} = 0$ | いいえ | はい | はい |
| 壁面での減衰 $O(y^3)$ | はい | はい | はい |
$\sigma$ モデルとVremanモデルは理論的にはWALEより優れた性質を持つが、実際のLES計算での結果の差は小さいことが多い。WALEの実績の豊富さと全ソルバーでの標準対応が、産業用途での優位性だ。
WALEモデルと燃焼LESの相性
大型ガスタービン燃焼器のLESでWALEモデルが採用されるケースが増えています。燃焼器内部は複雑形状の壁面が多く、壁面距離を求めること自体が重い前処理になる。動的モデルは負値のクリッピングロジックを燃焼ソルバーと統合するのが難しいという現場の声もあります。WALEはそのどちらの問題も回避できるシンプルさが武器。燃焼反応モデルとの結合では「流体モデルはシンプルに、化学モデルに計算資源を集中させる」という設計思想が合理的という判断です。
トラブルシューティング
よくある問題と対策
WALEモデルで問題が起きることはありますか？
1. SGS粘性が過小
症状: WALEの $\nu_{\text{sgs}}$ が非常に小さく、事実上のno-model LESになっている
原因: メッシュが細かすぎてSGS成分がほぼゼロ、または分母の $|\bar{S}|$ が支配的で $\nu_{\text{sgs}}$ が抑制
対策: これは問題ではなく正常な動作。メッシュが十分に細かければSGSモデルの寄与は小さくて当然だ。エネルギースペクトルで解像度を確認し、十分ならそのまま使えばよい。
2. 速度勾配の数値誤差
速度勾配の計算精度が重要って聞いたんですが。
原因: 非構造格子の歪んだセルで速度勾配の計算精度が低下し、$S_{ij}^d$ に誤差が混入
対策:
- Least Squares法で勾配を計算（Green-Gaussより精度が高い）
- セルのスキューネスを0.7以下に抑える
- メッシュの急激なサイズ変化を避ける
3. LES結果がRANSと変わらない
症状: 時間平均結果がSST k-omegaのRANS結果とほぼ同じ
原因: メッシュが粗すぎてLESが渦を解像できていない。SGS粘性がRANSの渦粘性と同程度になっている
対策:
- SGS粘性比 $\nu_{\text{sgs}}/\nu$ のコンター図を確認。$\nu_{\text{sgs}}/\nu > 10$ ならLES領域のメッシュが粗すぎる
- Pope criterion: 解像された乱流エネルギーが全体の80%以上か確認
- メッシュを細かくするかDDESに切替え
WALEモデル自体はほとんどトラブルがなく、問題の大半はメッシュ品質や解像度に起因するんですね。SGSモデルよりメッシュが重要という原則が再確認できました。
WALEとWMなのに結果が変わらない——そのときのチェックポイント
「WALEモデルに切り替えたのにSmagorinskyと結果がほぼ同じ」という相談がフォーラムでよく見られます。多くの場合、原因はSGSモデルの違いを上回る「格子依存性」です。格子が粗すぎると、どのSGSモデルを使っても大スケール渦の解像が不足して似たような結果になります。まずエネルギースペクトルのロールオフ位置を確認し、格子のカットオフ周波数が慣性域にあるかをチェックする——これがWALEのパフォーマンスを正しく評価する前提条件です。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——WALEモデルの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告