乱流モデル

乱流の程度を測りたい場合に、必要なものは全ての乱流運動のデータではなく、大まかな流れである場合が多い。実験や理論から、大きな渦は流れ場の形状の影響を強く受けるが、小さな渦は影響が少ない。そのため、工学的な用途から乱流を計算する場合、平均的な流れに関わる大きな渦の挙動を求めることが重要である。この解くべき流れを、ある種の平均操作を施すことによって導き、流れとそれ以外の乱流渦との間の関係を与える乱流モデルが必要になる。

物理学の未解決問題

乱流（特にその内部構造）の振る舞いを記述する理論上のモデルを構築することは可能か？

乱流を含む流れの解析手法は以下のように分類される^[1]。リストの後ろの方ほど厳密な取り扱いだが、計算負荷も増大する。

摩擦係数をレイノルズ数の関数として与えたり、熱伝達率（ヌセルト数）をレイノルズ数とプラントル数の関数として与えるなどの方法。流体力学の初歩段階で教えられ、大変有効な手法だが、単純流れにしか適用できない。コンピュータを用いた数値解析の必要はない。

運動方程式をある方向に積分し、残された方向への常微分方程式に変換する方法。

RANSに代表される、時間平均、空間平均あるいはアンサンブル平均して得られる方程式に基づく方法。

乱流流れを平均成分と変動成分の和として表すが、それらは切り離すことができず、乱れの効果は乱流応力という形で流れの方程式に現れる。乱流モデルではこれを近似することが必要となる。乱流モデルを用いる解析手法で対象とする大きな渦の変動を抽出するには何らかの特別な操作が必要となるが、これには大きく分けて二つの方法があり、一つはレイノルズ（Reynolds）平均と呼ばれる時間的および空間的平均操作に基づくRANSと呼ばれる方法で、もう一つは空間フィルタ操作に基づくLESと呼ばれる方法である。

定常流れ(もしくは非常にゆったりとした非定常流れ)だけを解くには、乱流運動の全てがモデル化の対象となり、(非圧縮性流体では)連続の式とナビエ-ストークス方程式に対してレイノルズ平均が施される。このレイノルズ平均が施された方程式をRANS（Reynolds-Averaged Navier-Stokes equations, RANS）と呼ぶ。

速度成分の2点相関に対する方程式、あるいはそのフーリエ変換を用いる方法。一様乱流以外にはほとんど用いられない。

乱流の比較的大きな構造を直接計算の対象とし、それより細かい乱れに対してモデル化を行う計算をラージエディシミュレーション（英語版）（Large-Eddy Simulation, LES）と呼ぶ。LESの方程式には、ナビエ-ストークス方程式に空間的に大小のスケールで分離するためにフィルターをかける。このことから、LESにおけるフィルターは格子平均と呼ばれることもある。

上記に挙げた乱流モデルを用いずに、乱流の全ての運動に対してナビエ-ストークス方程式を解く方法を乱流の直接数値シミュレーション（Direct Numerical Simulation, DNS）と呼ぶ。必要な計算格子の数はレイノルズ数とともに急速に増大するため、幾何形状が単純な低レイノルズ数流れにのみ適用可能である。

• 数値流体力学編集委員会　編『数値流体力学シリーズ３　乱流解析』東京大学出版会、1995年。 
• 梶島岳夫『乱流の数値シミュレーション』養賢堂、1999年。 　# 改訂版2014年、ISBN 978-4-84250526-8。
• 村上周三東京大学出版会、2000年。 
• 社団法人　土木学会　応用力学委員会　編『いまさら聞けない計算力学の常識』丸善、2008年。 
• Marcel Lesieur、Olivier Méteis、Pierre Comte：「乱流のシミュレーション：LESによる数値計算と可視化」、森北出版、ISBN 978-4-627-67331-1 (2010年9月10日)。
  • 原著：Large-Eddy Simulations of Turbulence, Cambrige Univ. Press (2005).
• 金田行雄（編）：「乱流の計算科学：乱流解明のツールとしての大規模数値シミュレーション」、共立出版、ISBN 978-4-320-12270-3 (2012年7月30日)。