直接数値シミュレーション(ちょくせつすうちシミュレーション、: Direct Numerical Simulation、DNS)とは、流れ(層流乱流)を数値的に解析するCFD手法の1つであり、以下の基礎方程式をそのまま(モデル化なしで)解き、流れに含まれる全ての大きさの渦をシミュレートすることである。

ここでp圧力、ρは密度、νは 粘性係数 、λは 第2粘性係数 である。各項はそれぞれ、

である。

特徴

一般にDNSは、全ての大きさの渦を解像できるように十分細かい計算格子を用意する必要とするため、必要となる計算格子は膨大なものとなり、実用上は解析困難な場合が多い。特に、高レイノルズ数の流れにおいては最小渦のスケールが微細となるためDNSの適用は困難である。そのため、実際には乱流モデルを用いて解析されることが多い。

格子スケールとレイノルズ数の関係

DNSを行うには、全ての渦を計算できるように細かい計算格子を必要とする。実際に必要な格子点数は、次元解析から調べることができる。乱流の最小渦の大きさは、 コルモゴロフ (Kolmogorov)のマイクロスケールηと呼ばれ、DNSではこの大きさまで細かく解像する必要がある。ここで、格子点の間隔をl と表すと次式となる。

ここで次元解析より、

と与えられる。 粘性散逸率 ε [m2/s3] は次式から与える。

これらより格子点数N を求めると、

と表される。

ここでRe は、流れを特徴付ける無次元数のレイノルズ数であり、次式で表される。

ここで、U は代表速度 [m/s]、L は代表長さ [m] であり、流れ場のスケールとして平均速度と物体の大きさで表される。レイノルズ数Re は流れの乱れ具合を表すもので、この値が大きい程流れは強く乱れたものである。上式より、レイノルズ数が大きいほど、計算に必要な格子点数は増える。

たとえば、人が速度 1 m/s程度で歩くことを考えれば、必要な格子点数はN 〜109 程度となり、その数値解析は容易ではない。


この記事では、Creative Commons Attribution-Share-Alike License 3.0の下に公開されているWikipediaの記事直接数値シミュレーションの資料を使用しています。