レイノルズ平均ナビエ-ストークス方程式 （英 : Reynolds-averaged Navier-Stokes equation, 略称: RANS ）とは、 時間平均化 [1] オズボーン・レイノルズ が提唱したレイノルズ分解（英語版 ） が方程式の前提にあり、レイノルズ分解によって流れの瞬間物理量は時間平均値と変動量に分けられる [2] 乱流 を記述するために用いられる。乱流特性に関する知識に基づく近似を用いることで、ナビエ-ストークス方程式 の近似時間平均解を与えることができる。
定常過程 における非圧縮性ニュートン流体 のRANS方程式は直交座標系 においてアインシュタインの縮約記法 を用いて次のように表される。
ρ
u
¯
j
∂
u
¯
i
∂
x
j
=
ρ
f
¯
i
+
∂
∂
x
j
[
−
p
¯
δ
i
j
+
μ
(
∂
u
¯
i
∂
x
j
+
∂
u
¯
j
∂
x
i
)
−
ρ
u
i
′
u
j
′
¯
]
{\displaystyle \rho {\bar {u}}_{j}{\frac {\partial {\bar {u}}_{i}}{\partial x_{j}}}=\rho {\bar {f}}_{i}+{\frac {\partial }{\partial x_{j}}}\left[-{\bar {p}}\delta _{ij}+\mu \left({\frac {\partial {\bar {u}}_{i}}{\partial x_{j}}}+{\frac {\partial {\bar {u}}_{j}}{\partial x_{i}}}\right)-\rho {\overline {u_{i}^{\prime }u_{j}^{\prime }}}\right]}
この方程式の左辺は、平均流量の不安定さと平均流による対流に起因する流体要素の平均運動量の変化を表す。平均物体力、平均圧力場に起因する等方性応力、粘性応力、変動する速度場に起因するレイノルズ応力（英語版 ） と呼ばれる見かけの応力
ρ
u
i
′
u
j
′
¯
{\displaystyle \rho {\overline {u_{i}^{\prime }u_{j}^{\prime }}}}
非線形のレイノルズ応力項は、RANS方程式を解くために追加のモデリングを必要とし、 それは多くの異なる乱流モデル を創出する。時間平均演算子
.
¯
{\displaystyle {\overline {.}}}
レイノルズ演算子（英語版 ） である。
RANS方程式の導出 ある瞬間のナビエ-ストークス方程式からRANS方程式を導出するのに必要な基本的ツールは、レイノルズ分解である。レイノルズ分解とは流れの成分（例えば流速
u
{\displaystyle u}
u
¯
{\displaystyle {\bar {u}}}
u
′
{\displaystyle u^{\prime }}
(
u
′
¯
=
0
)
{\displaystyle ({\bar {u^{\prime }}}=0)}
u
(
x
,
t
)
=
u
¯
(
x
)
+
u
′
(
x
,
t
)
{\displaystyle u({\boldsymbol {x}},t)={\bar {u}}({\boldsymbol {x}})+u^{\prime }({\boldsymbol {x}},t)}
x
=
(
x
,
y
,
z
)
{\displaystyle {\boldsymbol {x}}=(x,y,z)}
いくつかの文献 [3]
¯
{\displaystyle {\bar {\!}}}
u
¯
{\displaystyle {\bar {u}}}
U
{\displaystyle U}
u
′
{\displaystyle u^{\prime }}
u
{\displaystyle u}
u
,
u
¯
,
u
′
{\displaystyle u,{\bar {u}},u^{\prime }}
非圧縮粘性流体のナビエ-ストークス方程式をテンソル表記 により表すと次のようになる。
∂
u
i
∂
x
i
=
0
{\displaystyle {\frac {\partial u_{i}}{\partial x_{i}}}=0}
∂
u
i
∂
t
+
u
j
∂
u
i
∂
x
j
=
f
i
−
1
ρ
∂
p
∂
x
i
+
ν
∂
2
u
i
∂
x
j
∂
x
j
{\displaystyle {\frac {\partial u_{i}}{\partial t}}+u_{j}{\frac {\partial u_{i}}{\partial x_{j}}}=f_{i}-{\frac {1}{\rho }}{\frac {\partial p}{\partial x_{i}}}+\nu {\frac {\partial ^{2}u_{i}}{\partial x_{j}\partial x_{j}}}}
f
i
{\displaystyle f_{i}}
次にそれぞれの瞬間物理量を平均値と変動量に分けることで、結果以下の式となる。 [4]
∂
u
i
¯
∂
x
i
=
0
{\displaystyle {\frac {\partial {\bar {u_{i}}}}{\partial x_{i}}}=0}
∂
u
i
¯
∂
t
+
u
′
j
∂
u
′
i
∂
x
j
¯
=
f
i
¯
−
1
ρ
∂
p
¯
∂
x
i
+
ν
∂
2
u
i
¯
∂
x
j
∂
x
j
{\displaystyle {\frac {\partial {\bar {u_{i}}}}{\partial t}}+{\overline {{{u^{\prime }}_{j}}{\frac {\partial {u^{\prime }}_{i}}{\partial x_{j}}}}}={\bar {f_{i}}}-{\frac {1}{\rho }}{\frac {\partial {\bar {p}}}{\partial x_{i}}}+\nu {\frac {\partial ^{2}{\bar {u_{i}}}}{\partial x_{j}\partial x_{j}}}}
∂
u
i
∂
x
i
=
∂
u
i
¯
∂
x
i
+
∂
u
i
′
∂
x
i
=
0
{\displaystyle {\frac {\partial u_{i}}{\partial x_{i}}}={\frac {\partial {\bar {u_{i}}}}{\partial x_{i}}}+{\frac {\partial u_{i}^{\prime }}{\partial x_{i}}}=0}
∂
u
i
¯
∂
t
+
u
j
¯
∂
u
i
¯
∂
x
j
=
f
i
¯
−
1
ρ
∂
p
¯
∂
x
i
+
ν
∂
2
u
i
¯
∂
x
j
∂
x
j
−
∂
u
i
′
u
j
′
¯
∂
x
j
{\displaystyle {\frac {\partial {\bar {u_{i}}}}{\partial t}}+{\bar {u_{j}}}{\frac {\partial {\bar {u_{i}}}}{\partial x_{j}}}={\bar {f_{i}}}-{\frac {1}{\rho }}{\frac {\partial {\bar {p}}}{\partial x_{i}}}+\nu {\frac {\partial ^{2}{\bar {u_{i}}}}{\partial x_{j}\partial x_{j}}}-{\frac {\partial {\overline {u_{i}^{\prime }u_{j}^{\prime }}}}{\partial x_{j}}}}
さらに変形すると以下のようになる。
ρ
∂
u
i
¯
∂
t
+
ρ
u
j
¯
∂
u
i
¯
∂
x
j
=
ρ
f
i
¯
+
∂
∂
x
j
[
−
p
¯
δ
i
j
+
2
μ
S
i
j
¯
−
ρ
u
i
′
u
j
′
¯
]
{\displaystyle \rho {\frac {\partial {\bar {u_{i}}}}{\partial t}}+\rho {\bar {u_{j}}}{\frac {\partial {\bar {u_{i}}}}{\partial x_{j}}}=\rho {\bar {f_{i}}}+{\frac {\partial }{\partial x_{j}}}\left[-{\bar {p}}\delta _{ij}+2\mu {\bar {S_{ij}}}-\rho {\overline {u_{i}^{\prime }u_{j}^{\prime }}}\right]}
ここで
S
i
j
{\displaystyle S_{ij}}
S
i
j
=
1
2
(
∂
u
i
∂
x
j
+
∂
u
j
∂
x
i
)
{\displaystyle S_{ij}={\frac {1}{2}}{\biggl (}{\frac {\partial u_{i}}{\partial x_{j}}}+{\frac {\partial u_{j}}{\partial x_{i}}}{\biggr )}}
最後に、時間での積分により、時間依存性が削除されるため時間微分項を消去する。
ρ
u
j
¯
∂
u
i
¯
∂
x
j
=
ρ
f
i
¯
+
∂
∂
x
j
[
−
p
¯
δ
i
j
+
2
μ
S
i
j
¯
−
ρ
u
i
′
u
j
′
¯
]
{\displaystyle \rho {\bar {u_{j}}}{\frac {\partial {\bar {u_{i}}}}{\partial x_{j}}}=\rho {\bar {f_{i}}}+{\frac {\partial }{\partial x_{j}}}\left[-{\bar {p}}\delta _{ij}+2\mu {\bar {S_{ij}}}-\rho {\overline {u_{i}^{\prime }u_{j}^{\prime }}}\right]}
脚注
^ 変数
x
{\displaystyle x}
X
¯
{\displaystyle {\bar {X}}}
X
¯
=
lim
T
→
∞
1
T
∫
t
0
t
0
+
T
x
d
t
.
{\displaystyle {\bar {X}}=\lim _{T\to \infty }{\frac {1}{T}}\int _{t_{0}}^{t_{0}+T}x\,dt.}
時間平均値が一義的に定義されるためには、平均値
X
¯
{\displaystyle {\bar {X}}}
t
0
{\displaystyle t_{0}}
^ レイノルズ分解 CAE用語集｜ ソフトウェアクレイドル 、2018年10月15日閲覧。 ^ Tennekes, H.; Lumley, J. L. (1992). A first course in turbulence (14. print. ed.). Cambridge, Mass. [u.a.]: MIT Press. ISBN 978-0-262-20019-6 . ^ それぞれの瞬間物理量を平均値と変動量に分けることで、以下の式となる。
∂
(
u
i
¯
+
u
i
′
)
∂
x
i
=
0
{\displaystyle {\frac {\partial \left({\bar {u_{i}}}+u_{i}^{\prime }\right)}{\partial x_{i}}}=0}
∂
(
u
i
¯
+
u
i
′
)
∂
t
+
(
u
j
¯
+
u
j
′
)
∂
(
u
i
¯
+
u
i
′
)
∂
x
j
=
(
f
i
¯
+
f
i
′
)
−
1
ρ
∂
(
p
¯
+
p
′
)
∂
x
i
+
ν
∂
2
(
u
i
¯
+
u
i
′
)
∂
x
j
∂
x
j
.
{\displaystyle {\frac {\partial \left({\bar {u_{i}}}+u_{i}^{\prime }\right)}{\partial t}}+\left({\bar {u_{j}}}+u_{j}^{\prime }\right){\frac {\partial \left({\bar {u_{i}}}+u_{i}^{\prime }\right)}{\partial x_{j}}}=\left({\bar {f_{i}}}+f_{i}^{\prime }\right)-{\frac {1}{\rho }}{\frac {\partial \left({\bar {p}}+p^{\prime }\right)}{\partial x_{i}}}+\nu {\frac {\partial ^{2}\left({\bar {u_{i}}}+u_{i}^{\prime }\right)}{\partial x_{j}\partial x_{j}}}.}
さらに方程式を時間平均すると次式が得られる
∂
(
u
i
¯
+
u
i
′
)
∂
x
i
¯
=
0
{\displaystyle {\overline {\frac {\partial \left({\bar {u_{i}}}+u_{i}^{\prime }\right)}{\partial x_{i}}}}=0}
∂
(
u
i
¯
+
u
i
′
)
∂
t
¯
+
(
u
j
¯
+
u
j
′
)
∂
(
u
i
¯
+
u
i
′
)
∂
x
j
¯
=
(
f
i
¯
+
f
i
′
)
¯
−
1
ρ
∂
(
p
¯
+
p
′
)
∂
x
i
¯
+
ν
∂
2
(
u
i
¯
+
u
i
′
)
∂
x
j
∂
x
j
¯
.
{\displaystyle {\overline {\frac {\partial \left({\bar {u_{i}}}+u_{i}^{\prime }\right)}{\partial t}}}+{\overline {\left({\bar {u_{j}}}+u_{j}^{\prime }\right){\frac {\partial \left({\bar {u_{i}}}+u_{i}^{\prime }\right)}{\partial x_{j}}}}}={\overline {\left({\bar {f_{i}}}+f_{i}^{\prime }\right)}}-{\frac {1}{\rho }}{\overline {\frac {\partial \left({\bar {p}}+p^{\prime }\right)}{\partial x_{i}}}}+\nu {\overline {\frac {\partial ^{2}\left({\bar {u_{i}}}+u_{i}^{\prime }\right)}{\partial x_{j}\partial x_{j}}}}.}
また、非線形項を以下のように簡略化することができる
u
i
u
i
¯
=
(
u
i
¯
+
u
i
′
)
(
u
i
¯
+
u
i
′
)
¯
=
u
i
¯
u
i
¯
+
u
i
¯
u
i
′
+
u
i
′
u
i
¯
+
u
i
′
u
i
′
¯
=
u
i
¯
u
i
¯
+
u
i
′
u
i
′
¯
{\displaystyle {\overline {u_{i}u_{i}}}={\overline {\left({\bar {u_{i}}}+u_{i}^{\prime }\right)\left({\bar {u_{i}}}+u_{i}^{\prime }\right)}}={\overline {{\bar {u_{i}}}{\bar {u_{i}}}+{\bar {u_{i}}}u_{i}^{\prime }+u_{i}^{\prime }{\bar {u_{i}}}+u_{i}^{\prime }u_{i}^{\prime }}}={\bar {u_{i}}}{\bar {u_{i}}}+{\overline {u_{i}^{\prime }u_{i}^{\prime }}}}
![{\displaystyle \rho {\bar {u}}_{j}{\frac {\partial {\bar {u}}_{i}}{\partial x_{j}}}=\rho {\bar {f}}_{i}+{\frac {\partial }{\partial x_{j}}}\left[-{\bar {p}}\delta _{ij}+\mu \left({\frac {\partial {\bar {u}}_{i}}{\partial x_{j}}}+{\frac {\partial {\bar {u}}_{j}}{\partial x_{i}}}\right)-\rho {\overline {u_{i}^{\prime }u_{j}^{\prime }}}\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/55f5ce9608974c0a3fb0bae737071f65bdd1bb6e)
![{\displaystyle \rho {\frac {\partial {\bar {u_{i}}}}{\partial t}}+\rho {\bar {u_{j}}}{\frac {\partial {\bar {u_{i}}}}{\partial x_{j}}}=\rho {\bar {f_{i}}}+{\frac {\partial }{\partial x_{j}}}\left[-{\bar {p}}\delta _{ij}+2\mu {\bar {S_{ij}}}-\rho {\overline {u_{i}^{\prime }u_{j}^{\prime }}}\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/55d9b40bc4afefd0522f20e117de4068287863d5)
![{\displaystyle \rho {\bar {u_{j}}}{\frac {\partial {\bar {u_{i}}}}{\partial x_{j}}}=\rho {\bar {f_{i}}}+{\frac {\partial }{\partial x_{j}}}\left[-{\bar {p}}\delta _{ij}+2\mu {\bar {S_{ij}}}-\rho {\overline {u_{i}^{\prime }u_{j}^{\prime }}}\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1abc38e1557474fc40b63b261909fe3ded4d45b0)
