查看“︁大涡模拟”︁的源代码

[[File:LES_Turbulent_Velocity_Field.png|link=//upload.wikimedia.org/wikipedia/commons/thumb/8/85/LES_Turbulent_Velocity_Field.png/300px-LES_Turbulent_Velocity_Field.png|right|thumb|300x300px|[[湍流]]气体速度场的大涡模拟。]]
'''大涡模拟'''('''Large eddy simulation,''' '''LES''') 是用于[[计算流体力学|计算流体动力学]]的[[湍流]]数学模型。它最初由约瑟夫·斯玛格林斯基于 1963 年提出，用于模拟大气气流， <ref name="Smagorinsky_1963">{{Cite journal |last=Smagorinsky |first=Joseph |date=March 1963 |title=General Circulation Experiments with the Primitive Equations |url=https://archive.org/details/sim_monthly-weather-review_1963-03_91_3/page/99 |journal=Monthly Weather Review |volume=91 |issue=3 |page=99–164 |bibcode=1963MWRv...91...99S |doi=10.1175/1520-0493(1963)091<0099:GCEWTP>2.3.CO;2 |doi-access=free}}</ref>并由迪尔多夫(1970)推广。 <ref name="Deardorff_1970">{{Cite journal |last=Deardorff |first=James |year=1970 |title=A numerical study of three-dimensional turbulent channel flow at large Reynolds numbers |url=https://archive.org/details/sim_journal-of-fluid-mechanics_1970-04-13_41_2/page/n214 |journal=[[Journal of Fluid Mechanics]] |volume=41 |issue=2 |page=453–480 |bibcode=1970JFM....41..453D |doi=10.1017/S0022112070000691 |s2cid=121884175}}</ref> LES 目前广泛应用于多个研究领域，包括[[燃烧]]、 <ref name="Pitsch_2006">{{Cite journal |last=Pitsch |first=Heinz |year=2006 |title=Large-Eddy Simulation of Turbulent Combustion |url=http://caltechconf.library.caltech.edu/145/1/Abstract1.pdf |journal=Annual Review of Fluid Mechanics |volume=38 |issue=1 |page=453–482 |bibcode=2006AnRFM..38..453P |doi=10.1146/annurev.fluid.38.050304.092133 |s2cid=5487815 |access-date=2022-10-25 |archive-date=2023-03-07 |archive-url=https://web.archive.org/web/20230307163055/http://caltechconf.library.caltech.edu/145/1/Abstract1.pdf |dead-url=no }}</ref>声学、 <ref name="Wagner_2007">{{Cite book|last=Wagner|first=Claus|last2=Hüttl|first2=Thomas|last3=Sagaut|first3=Pierre|title=Large-Eddy Simulation for Acoustics|publisher=Cambridge University Press|year=2007|isbn=978-0-521-87144-0}}</ref>和大气边界层模拟。 <ref name="Sullivan_1994">{{Cite journal |last=Sullivan |first=Peter P. |last2=McWilliams |first2=James C. |last3=Moeng |first3=Chin-Hoh |year=1994 |title=A subgrid-scale model for large-eddy simulation of planetary boundary-layer flows |journal=Boundary-Layer Meteorology |volume=71 |page=247–276 |bibcode=1994BoLMe..71..247S |citeseerx=10.1.1.463.6006 |doi=10.1007/BF00713741 |issn=0006-8314 |s2cid=53051046 |number=3}}</ref>

通过数值求解[[纳维-斯托克斯方程|Navier-Stokes 方程]]来模拟湍流需要解决非常广泛的时间和长度尺度，所有这些都会影响流场。这样的分辨率可以通过直接数值模拟(DNS) 来实现，但 DNS 的计算成本很高，难以模拟具有复杂几何形状或流动配置的实际工程系统，例如湍流喷射、泵、车辆和起落架。

LES 背后的主要思想是通过 Navier-Stokes 方程的[[濾波器|低通滤波]]忽略最小长度尺度来降低计算成本。这种可以被视为时间和空间平均的低通滤波，有效地从数值解中去除了小尺度信息。然而，这些信息并不是无关紧要的，它对流场的影响必须被建模，由此衍生小尺度可以发挥重要作用的问题的活跃研究领域，例如近壁流， <ref name="piomelli-02">{{Cite journal |last=Piomelli, Ugo |last2=Elias Balaras |year=2002 |title=Wall-layer models for large-eddy simulations. |journal=Annual Review of Fluid Mechanics |volume=34 |issue=34 |page=349–374 |bibcode=2002AnRFM..34..349P |doi=10.1146/annurev.fluid.34.082901.144919}}</ref> <ref name="spalart-09">{{Cite journal |last=Spalart, P. R. |year=2009 |title=Detached-eddy simulation. |journal=Annual Review of Fluid Mechanics |volume=41 |issue=1 |page=181–202 |bibcode=2009AnRFM..41..181S |doi=10.1146/annurev.fluid.010908.165130}}</ref>反应流、 <ref name="Pitsch_2006" />和多相流。 <ref name="fox-12">{{Cite journal |last=Fox, R. O. |year=2012 |title=Large-eddy-simulation tools for multiphase flows. |journal=Annual Review of Fluid Mechanics |volume=44 |issue=1 |page=47–76 |bibcode=2012AnRFM..44...47F |doi=10.1146/annurev-fluid-120710-101118}}</ref>


== 过滤器定义和属性 ==
[[File:DNS_Velocity_Field.png|link=//upload.wikimedia.org/wikipedia/commons/thumb/7/77/DNS_Velocity_Field.png/300px-DNS_Velocity_Field.png|right|thumb|300x300px|由均匀衰减湍流的直接数值模拟(DNS) 产生的速度场。域大小为<math>L^3</math> .]]
[[File:DNS_Filtered_Velocity_Field_Small.png|link=//upload.wikimedia.org/wikipedia/commons/thumb/d/d4/DNS_Filtered_Velocity_Field_Small.png/300px-DNS_Filtered_Velocity_Field_Small.png|right|thumb|300x300px|使用盒式过滤器过滤的相同 DNS 速度字段和<math>\Delta=L/32</math> .]]
[[File:DNS_Filtered_Velocity_Field_Large.png|link=//upload.wikimedia.org/wikipedia/commons/thumb/8/89/DNS_Filtered_Velocity_Field_Large.png/300px-DNS_Filtered_Velocity_Field_Large.png|right|thumb|300x300px|使用盒式过滤器过滤的相同 DNS 速度字段和<math>\Delta=L/16</math> .]]
LES 滤波器可应用于空间和时间场<math>\phi(\boldsymbol{x},t)</math>并执行空间滤波操作、时间滤波操作或两者一并进行。过滤后的流场加了上划线，定义为： <ref name="Pope_2000"/> <ref name="Sagaut_2006">{{Cite book|last=Sagaut, Pierre|title=Large Eddy Simulation for Incompressible Flows|publisher=Springer|year=2006|edition=Third|isbn=978-3-540-26344-9}}</ref>

: <math>
\overline{\phi(\boldsymbol{x},t)} = \displaystyle{
\int_{-\infty}^{\infty}} \int_{-\infty}^{\infty} \phi(\boldsymbol{r},\tau) G(\boldsymbol{x}-\boldsymbol{r},t - \tau) d\tau d \boldsymbol{r}
</math>

其中<math>G</math>是滤波器卷积核。这也可以写成：

: <math>
\overline{\phi} = G \star \phi .
</math>

过滤器内核<math>G</math>有一个相关的截止长度尺度<math>\Delta</math>和截止时间尺度<math>\tau_{c}</math> .小于这些的尺度将被消掉。使用上述过滤器定义，任何流场<math>\phi</math>可以分为过滤和子过滤（用素数表示）部分，如

: <math>
\phi = \bar{\phi} + \phi^{\prime} .
</math>

需要注意的是，大涡模拟滤波操作不满足雷诺算子的性质。

== 过滤的控制方程 ==
LES的控制方程是通过过滤控制流场的[[偏微分方程]]得到的<math>\rho \boldsymbol{u}(\boldsymbol{x},t)</math> .不可压缩和可压缩 LES 控制方程之间存在差异，这导致了新滤波操作的定义。

=== 不可压缩流动 ===
对于不可压缩流动，对[[連續性方程式|连续性方程]]和 Navier-Stokes 方程进行滤波，得到滤波后的不可压缩连续性方程，

: <math>
\frac{ \partial \bar{u_i} }{ \partial x_i } = 0
</math>

和过滤后的 Navier-Stokes 方程，

: <math>
\frac{ \partial \bar{u_i} }{ \partial t }
+ \frac{ \partial }{ \partial x_j } \left( \overline{ u_i u_j } \right) 
= - \frac{1}{\rho} \frac{ \partial \overline{p} }{ \partial x_i } 
+ \nu \frac{\partial}{\partial x_j} \left( \frac{ \partial \bar{u_i} }{ \partial x_j } + \frac{ \partial \bar{u_j} }{ \partial x_i } \right)
= - \frac{1}{\rho} \frac{ \partial \overline{p} }{ \partial x_i } 
+ 2 \nu \frac{\partial}{\partial x_j} \bar{S}_{ij},
</math>

其中<math>\bar{p}</math>是过滤后的压力场和<math>\bar{S}_{ij}</math>是使用过滤速度评估的应变率张量。[[非線性系統|非线性]]滤波平流项<math>\overline{u_i u_j}</math>是LES建模困难的主要原因。它需要未过滤的速度场，这是未知的，因此必须对其进行建模。下面的分析说明了非线性带来的困难，即它导致大小尺度之间的相互作用，防止尺度分离。

过滤后的平流项可以按照莱昂纳德(1975)<ref name="Leonard_1974">{{Cite book|last=Leonard|first=A.|title=Energy cascade in large-eddy simulations of turbulent fluid flows|work=Advances in Geophysics A|year=1975|isbn=9780120188185|series=Advances in Geophysics|volume=18|pages=237–248|bibcode=1975AdGeo..18..237L|doi=10.1016/S0065-2687(08)60464-1}}</ref>拆分为：

: <math>
\overline{u_i u_j} = \tau_{ij} + \overline{u}_i \overline{u}_j
</math>

其中<math>\tau_{ij}</math>是残余应力张量，因此过滤后的 Navier-Stokes 方程变为

: <math>
\frac{ \partial \bar{u_i} }{ \partial t }
+ \frac{ \partial }{ \partial x_j } \left( \overline{u}_i \overline{u}_j \right) 
= - \frac{1}{\rho} \frac{ \partial \overline{p} }{ \partial x_i } 
+ 2 \nu \frac{\partial}{\partial x_j} \bar{S}_{ij}
- \frac{ \partial \tau_{ij} }{ \partial x_j }
</math>


与残余应力张量<math>\tau_{ij}</math>对所有未封闭的项进行分组。 Leonard 将这个应力张量分解为<math>\tau_{ij} = L_{ij} + C_{ij} + R_{ij}</math>并为每个项提供物理解释。 <math>L_{ij} = \overline{ \bar{u}_{i} \bar{u}_{j} } - \bar{u}_{i} \bar{u}_{j}</math> 是伦纳德张量，代表大尺度之间的相互作用；<math>R_{ij} = \overline{u^{\prime}_{i} u^{\prime}_{j}}</math> 为类雷诺应力项，表示子滤波器尺度 (SFS) 之间的相互作用；<math>C_{ij} = \overline{\bar{u}_{i} u^{\prime}_{j}} + \overline{\bar{u}_{j} u^{\prime}_{i}}
</math> 是克拉克张量， <ref name="Clark">{{Cite journal |last=Clark |first=R. |last2=Ferziger |first2=J. |last3=Reynolds |first3=W. |author-link3=William Craig Reynolds |year=1979 |title=Evaluation of subgrid-scale models using an accurately simulated turbulent flow |url=https://archive.org/details/sim_journal-of-fluid-mechanics_1979-03-09_91_1/page/n2 |journal=[[Journal of Fluid Mechanics]] |volume=91 |page=1–16 |bibcode=1979JFM....91....1C |doi=10.1017/S002211207900001X |s2cid=120228458}}</ref>表示大尺度和小尺度之间的跨尺度相互作用。 <ref name="Leonard_1974" />建模未闭合项<math>\tau_{ij}</math>是亚网格尺度（SGS）模型需要解析的项。亚网格应力张量<math>\tau_{ij}</math>的存在，使这变得具有挑战性必须考虑所有尺度之间的相互作用，包括过滤尺度与未过滤尺度。

被动标量<math>\phi</math>的滤波控制方程 ，例如混合分数或温度，可以写成

: <math>
\frac{ \partial \overline{\phi} }{ \partial t }
+ \frac{\partial}{\partial x_j} \left( \overline{u}_j \overline{\phi} \right)
= \frac{\partial \overline{J_{\phi}} }{\partial x_j} 
+ \frac{ \partial q_j }{ \partial x_j }
</math>

其中<math>J_{\phi}</math>是<math>\phi</math>的扩散通量，<math>q_j</math>是标量<math>\phi</math>的子滤波器通量。过滤后的扩散通量<math>\overline{J_{\phi}}</math>是未封闭的，除非假设它具有特定的形式，例如梯度扩散模型<math>J_{\phi} = D_{\phi} \frac{ \partial \phi }{ \partial x_i }</math> . <math>q_j</math>的定义形式类似<math>\tau_{ij}</math> ：

: <math>
q_j = \bar{\phi} \overline{u}_j - \overline{\phi u_j}
</math>

并且可以类似地分解为各种尺度之间相互作用的贡献。这种子过滤器通量也需要子过滤器模型。

==== 推导 ====
使用[[爱因斯坦求和约定]]，笛卡尔坐标中不可压缩流体的 Navier-Stokes 方程为

: <math> \frac{\partial u_i}{\partial x_i} = 0 </math>
: <math> \frac{\partial u_i}{\partial t} + \frac{\partial u_iu_j}{\partial x_j}
= - \frac{1}{\rho} \frac{\partial p}{\partial x_i}
+ \nu \frac{\partial^2 u_i}{\partial x_j \partial x_j}.
</math>

过滤动量方程导致

: <math> \overline{\frac{\partial u_i}{\partial t}} + \overline{\frac{\partial u_iu_j}{\partial x_j}}
= - \overline{\frac{1}{\rho} \frac{\partial p}{\partial x_i}}
+ \overline{\nu \frac{\partial^2 u_i}{\partial x_j \partial x_j}}.
</math>

如果我们假设过滤和微分对等，那么

: <math> \frac{\partial \bar{u_i}}{\partial t} + \overline{\frac{\partial u_iu_j}{\partial x_j}}
= - \frac{1}{\rho} \frac{\partial \bar{p}}{\partial x_i}
+ \nu \frac{\partial^2 \bar{u_i}}{\partial x_j \partial x_j}.
</math>

这个方程模拟了过滤变量的时间变化<math>\bar{u_i}</math> .由于未过滤的变量<math>u_i</math>未知，无法直接计算<math>\overline{\frac{\partial u_iu_j}{\partial x_j}}</math> 。然而，<math> \frac{\partial \bar{u_i}\bar{u_j}}{\partial x_j}</math>是已知的。进行了替换：

: <math> \frac{\partial \bar{u_i}}{\partial t} + \frac{\partial \bar{u_i}\bar{u_j}}{\partial x_j}
= - \frac{1}{\rho} \frac{\partial \bar{p}}{\partial x_i}
+ \nu \frac{\partial^2 \bar{u_i}}{\partial x_j \partial x_j}
- \left(\overline{ \frac{\partial u_iu_j}{\partial x_j}} - \frac{\partial \bar{u_i}\bar{u_j}}{\partial x_j}\right).
</math>

令<math>\tau_{ij} = \overline{u_i u_j} - \bar{u}_{i} \bar{u}_{j}</math> 。得到的方程组是 LES 方程：

: <math> \frac{\partial \bar{u_i}}{\partial t} + \bar{u_j} \frac{\partial \bar{u_i}}{\partial x_j}
= - \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\tau_{ij}}{\partial x_j}.
</math>

=== 可压缩控制方程 ===
对于可压缩流动的控制方程，每个方程都从质量守恒开始过滤。这给出了：

: <math>
\frac{\partial \overline{\rho}}{\partial t} + \frac{ \partial \overline{u_i \rho} }{\partial x_i} = 0
</math>

这导致了一个额外的子过滤器项。然而，希望避免必须对质量守恒方程的子过滤器尺度进行建模。出于这个原因，Favre <ref name="Favre_1983">{{Cite journal |last=Favre, Alexandre |year=1983 |title=Turbulence: space-time statistical properties and behavior in supersonic flows |journal=Physics of Fluids A |volume=23 |issue=10 |page=2851–2863 |bibcode=1983PhFl...26.2851F |doi=10.1063/1.864049}}</ref>提出了一种密度加权滤波操作，称为 Favre 滤波，定义为任意量<math>\phi</math>为：

: <math>
\tilde{\phi} = \frac{ \overline{\rho \phi} }{ \overline{\rho} }
</math>

在不可压缩的限制下，它变成了正常的过滤操作。这使得质量守恒方程化为：

: <math>
\frac{\partial \overline{\rho}}{\partial t} + \frac{ \partial \overline{\rho} \tilde{u_i} }{ \partial x_i } = 0.
</math>

然后可以将此概念扩展到编写可压缩流的 Favre 过滤动量方程。根据弗雷曼： <ref name="Vreman_1995">{{Cite journal |last=Vreman, Bert |last2=Geurts, Bernard |last3=Kuerten, Hans |year=1995 |title=Subgrid-modelling in LES of compressible flow |url=https://research.utwente.nl/en/publications/subgridmodelling-in-les-of-compressible-flow(7c54958d-ebdd-4422-bf9d-c17052984a68).html |journal=[[Applied Scientific Research]] |volume=45 |issue=3 |page=191–203 |doi=10.1007/BF00849116}}</ref>

: <math>
\frac{ \partial \overline{\rho} \tilde{u_i} }{ \partial t }
+ \frac{ \partial \overline{\rho} \tilde{u_i} \tilde{u_j} }{ \partial x_j }
+ \frac{ \partial \overline{p} }{ \partial x_i }
- \frac{ \partial \tilde{\sigma}_{ij} }{ \partial x_j }
= - \frac{ \partial \overline{\rho} \tau_{ij}^{r} }{ \partial x_j }
+ \frac{ \partial }{ \partial x_j } \left( \overline{\sigma}_{ij} - \tilde{\sigma}_{ij} \right)
</math>

其中<math>\sigma_{ij}</math>是剪应力张量，对于牛顿流体，由下式给出：

: <math>
\sigma_{ij} = 2 \mu(T) S_{ij} - \frac{2}{3} \mu(T) \delta_{ij} S_{kk}
</math>

<math>\frac{ \partial }{\partial x_j} \left( \overline{\sigma}_{ij} - \tilde{\sigma}_{ij} \right)</math>这一项表示评估粘度的子过滤器粘性贡献<math>\mu(T)</math>使用 Favre 过滤温度<math>\tilde{T}</math> . Favre 滤波动量场的亚网格应力张量：<math>
\tau_{ij}^{r} = \widetilde{ u_i \cdot u_j } - \tilde{u_i} \tilde{u_j}
</math>

以此类推，莱昂纳德分解也可以写成滤波三重积的残余应力张量<math>\overline{\rho \phi \psi}</math> 。三重乘积可以使用 Favre 过滤运算符重写为<math>\overline{\rho} \widetilde{\phi \psi}</math> ，这是一个未封闭的项（它需要场<math>\phi</math>和<math>\psi</math>的信息, 当只有<math>\tilde{\phi}</math>和<math>\tilde{\psi}</math>两项是已知的时候）。它可以以类似于上面提到的<math>\overline{u_i u_j}</math>的方式分解，这导致子滤波器应力张量<math>\overline{\rho} \left( \widetilde{\phi \psi} - \tilde{\phi} \tilde{\psi} \right)</math> 。这个子过滤器项可以分为三种类型的相互作用的贡献：伦达德张量<math>L_{ij}</math>, 代表解析尺度之间的相互作用；克拉克张量<math>C_{ij}</math> ，表示已解析和未解析的尺度之间的相互作用；和雷诺张量<math>R_{ij}</math> ，它表示未解析的尺度之间的相互作用。 <ref name="Sagaut_2009">{{Cite book|last=Garnier, E.|last2=Adams, N.|last3=Sagaut, P.|title=Large eddy simulation for compressible flows|year=2009|publisher=Springer|isbn=978-90-481-2818-1|doi=10.1007/978-90-481-2819-8|url=http://cds.cern.ch/record/1339029}}</ref>

:

=== 滤波动能方程 ===
除了过滤的质量和动量方程之外，过滤动能方程可以提供额外的见解。可以过滤动能场以产生总过滤动能：

: <math>
\overline{E} = \frac{1}{2} \overline{ u_i u_i }
</math>

过滤后的总动能可以分解为两项：过滤后的速度场的动能<math>E_f</math> ,

: <math>
E_f = \frac{1}{2} \overline{u_i} \, \overline{u_i}
</math>

和剩余动能<math>k_r</math> ,

: <math>
k_r = \frac{1}{2} \overline{ u_i u_i } - \frac{1}{2} \overline{u_i} \, \overline{u_i} = \frac{1}{2} \tau_{ii}^{r}
</math>

这样<math>\overline{E} = E_f + k_r</math> .

守恒方程为<math>E_f</math>可以通过将滤波后的动量传输方程乘以<math>\overline{u_i}</math>产生：

: <math>
\frac{\partial E_f}{\partial t} 
+ \overline{u_j} \frac{\partial E_f}{\partial x_j} 
+ \frac{1}{\rho} \frac{\partial \overline{u_i} \bar{p} }{ \partial x_i }
+ \frac{\partial \overline{u_i} \tau_{ij}^{r}}{\partial x_j} 
- 2 \nu \frac{ \partial \overline{u_i} \bar{S_{ij}} }{ \partial x_j }
= 
- \epsilon_{f} 
- \Pi
</math>

其中<math>\epsilon_{f} = 2 \nu \bar{S_{ij}} \bar{S_{ij}}</math>是粘性应力对过滤速度场的动能的耗散，并且<math>\Pi = -\tau_{ij}^{r} \bar{S_{ij}}</math>表示动能的子过滤尺度 (SFS) 耗散。

左边的项代表运输，右边的项是耗散动能的汇项。 <ref name="Pope_2000"/>

<math>\Pi</math> SFS 耗散项特别令人感兴趣，因为它代表能量从大分辨尺度到小未分辨尺度的转移。一般， <math>\Pi</math>将能量从大尺度转移到小尺度。然而，瞬间<math>\Pi</math>可以是积极的''或''消极的，这意味着它也可以作为源项<math>E_f</math>, 滤波后的速度场的动能。能量从未解析到已解析尺度的传递称为'''反向散射'''（同样，能量从已解析到未解析尺度的传递称为'''前向散射'''）。 <ref name="Piomelli_1991">{{Cite journal |last=Piomelli, U. |last2=Cabot, W. |last3=Moin, P. |author-link3=Parviz Moin |last4=Lee, S. |year=1991 |title=Subgrid-scale backscatter in turbulent and transitional flows |url=https://semanticscholar.org/paper/b954cf23648b5a86c716490fd8637bea54325841 |journal=Physics of Fluids A |volume=3 |issue=7 |page=1766–1771 |bibcode=1991PhFl....3.1766P |doi=10.1063/1.857956 |s2cid=54904570}}</ref>

== LES的数值方法 ==
大涡模拟涉及使用[[计算流体力学|计算流体动力学]]求解离散滤波控制方程。 LES 从域大小中解析尺度<math>L</math>过滤器尺寸<math>\Delta</math> ，因此必须解决很大一部分高波数湍流波动。这需要高阶数值方案，或者如果使用低阶数值方案，则需要精细的网格分辨率。 波普 <ref name="Pope_2000">{{Cite book|title=Turbulent Flows|year=2000|publisher=Cambridge University Press|last=Pope, S. B.}}</ref>的第 13 章解决了网格分辨率有多精细的问题<math>\Delta x</math>需要解析过滤的速度场<math>\overline{u}(\boldsymbol{x})</math> . Ghosal <ref name="Ghosal_1996">{{Cite journal |last=Ghosal, S. |date=April 1996 |title=An analysis of numerical errors in large-eddy simulations of turbulence |journal=[[Journal of Computational Physics]] |volume=125 |issue=1 |page=187–206 |bibcode=1996JCoPh.125..187G |doi=10.1006/jcph.1996.0088}}</ref>发现，对于低阶离散化方案，例如在有限体积方法中使用的那些，截断误差可以与子滤波器尺度贡献相同，除非滤波器宽度<math>\Delta</math>比网格间距<math>\Delta x</math>大得多。虽然偶数阶方案具有截断误差，但它们是非耗散的， <ref name="Leveque_1992">{{Cite book|title=Numerical Methods for Conservation Laws|last=Randall J. Leveque|year=1992|publisher=Birkhäuser Basel|edition=2nd|isbn=978-3-7643-2723-1}}</ref>并且由于子滤波器尺度模型是耗散的，偶数阶方案不会像耗散方案那样强烈地影响子滤波器尺度模型的贡献。

=== 过滤器实现 ===
大涡模拟中的滤波操作可以是隐式的，也可以是显式的。隐式过滤认识到子过滤器比例模型将以与许多数值方案相同的方式消散。通过这种方式，可以假设网格或数值离散化方案是 LES 低通滤波器。虽然这充分利用了网格分辨率，并消除了计算子滤波器比例模型项的计算成本，但很难确定与一些数值问题相关的 LES 滤波器的形状。此外，截断误差也可能成为问题。 <ref name="Grinstein_2007">{{Cite book|title=Implicit large eddy simulation|url=https://archive.org/details/isbn_9780521869829|last=Grinstein, Fernando|last2=Margolin, Len|last3=Rider, William|year=2007|publisher=Cambridge University Press|isbn=978-0-521-86982-9}}</ref>

在显式滤波中， LES 滤波器应用于离散的 Navier-Stokes 方程，提供明确定义的滤波器形状并减少截断误差。然而，显式过滤比隐式过滤需要更精细的网格，并且计算量与<math>(\Delta x)^4</math>成正比。Sagaut (2006) 的第 8 章更详细地介绍了 LES 的数值问题。 <ref name="Sagaut_2006" />

== 大涡模拟的边界条件 ==
入口边界条件对LES的精度影响很大，对于LES入口条件的处理是一个复杂的问题。理论上，一个良好的 LES 边界条件应包含以下特征： <ref name="Li_P">{{Cite journal |last=Li |first=P. |last2=Eckels |first2=S. |last3=Mann |first3=G. |last4=Zhang |first4=N. |year=2018 |title=A Method of Measuring Turbulent Flow Structures With Particle Image Velocimetry and Incorporating Into Boundary Conditions of Large Eddy Simulations |journal=Journal of Fluids Engineering |publisher=ASME International |volume=140 |issue=7 |at=071401-071401-11 |doi=10.1115/1.4039256}}</ref>

(1) 提供准确的流动特性信息，即速度和湍流；

(2) 满足 Navier-Stokes 方程和其他物理；

(3)易于实施和适应不同情况。

目前，为 LES 生成入口条件的方法大致分为 Tabor 等人分类的两类： <ref>{{Cite journal |last=Tabor |first=G. R. |last2=Baba-Ahmadi |first2=M. H. |year=2010 |title=Inlet conditions for large eddy simulation: a review |journal=Computers & Fluids |volume=39 |issue=4 |page=553&ndash;567 |doi=10.1016/j.compfluid.2009.10.007}}</ref>

产生湍流入口的第一种方法是根据具体情况合成它们，例如傅里叶技术、原理正交分解（POD）和涡流方法。合成技术试图在入口处构建具有合适的类湍流特性的湍流场，并使其易于指定湍流参数，例如湍流动能和湍流耗散率。此外，使用随机数生成的入口条件在计算上并不昂贵。然而，该方法存在一个严重的缺陷。合成的湍流不满足由 Navier-Stokes 方程控制的流体流动的物理结构。 <ref name="Li_P"/>

第二种方法涉及一个单独的前体计算，以生成一个湍流数据库，该数据库可以引入到入口处的主要计算中。数据库（有时称为“库”）可以通过多种方式生成，例如循环域、预先准备好的库和内部映射。然而，前驱体模拟产生湍流流入的方法需要很大的计算能力。

研究人员检查了各种合成和前体计算的应用，发现入口湍流越真实，LES 预测结果就越准确。 <ref name="Li_P"/>

== 模型未解析的尺度 ==
为了讨论未解析尺度的建模，首先必须对未解析尺度进行分类。它们分为两组：'''已解决的子过滤器尺度'''(SFS) 和'''子网格尺度'''(SGS)。

解析的子滤波器尺度表示波数大于截止波数的尺度<math>k_c</math> ，但其影响被过滤器抑制。仅当使用波空间中的非局部滤波器（例如盒式或高斯滤波器）时，才存在已解析的子滤波器尺度。这些解析的子过滤器尺度必须使用过滤器重建来建模。

子网格比例是小于截止滤波器宽度<math>\Delta</math>的任何比例。SGS 模型的形式取决于过滤器的实现。如LES 数值方法部分所述，如果考虑隐式 LES，则不实施 SGS 模型，并且假设离散化的数值效应模拟未解决的湍流运动的物理特性。

=== 子网格比例模型 ===
如果没有普遍有效的湍流描述，则在构建和应用 SGS 模型时必须利用经验信息，并辅以伽利略不变性<ref name="Pope_2000"/>等基本物理约束。<ref name="Meneveau_2010">{{Cite journal |last=Meneveau, C. |year=2010 |title=Turbulence: Subgrid-Scale Modeling |journal=Scholarpedia |volume=5 |issue=1 |page=9489 |bibcode=2010SchpJ...5.9489M |doi=10.4249/scholarpedia.9489 |doi-access=free}}</ref>存在两类 SGS 模型；第一类是'''功能模型'''，第二类是'''结构模型'''。一些模型可能被归类为两者。

==== 功能（涡粘）模型 ====
功能模型比结构模型更简单，只关注以物理上正确的速率耗散能量。这些基于人工涡流粘度方法，其中湍流的影响集中在湍流粘度中。该方法将亚网格尺度上的动能耗散视为类似于分子扩散。在这种情况下，偏斜部分<math>\tau_{ij}</math>被建模为：

: <math>
\tau_{ij}^r - \frac{1}{3} \tau_{kk} \delta_{ij} = -2 \nu_\mathrm{t} \bar{S}_{ij}
</math>

其中<math>\nu_\mathrm{t}</math>是湍流涡流粘度和<math>\bar{S}_{ij} = \frac{1}{2} \left( \frac{\partial \bar{u}_i }{\partial x_j} + \frac{\partial \bar{u}_j}{ \partial x_i} \right)</math>是应变率张量。

根据量纲分析，涡流粘度的单位必须为<math>\left[ \nu_\mathrm{t} \right] = \frac{\mathrm{m^2}}{\mathrm{s}}</math> 。大多数涡流粘度 SGS 模型将涡流粘度建模为特征长度尺度和特征速度尺度的乘积。

===== Smagorinsky-Lilly 模型 =====
第一个成功开发的 SGS 模型是 Smagorinsky-Lilly SGS 模型，它由Smagorinsky <ref name="Smagorinsky_1963"/>开发并用于 Deardorff 的第一个 LES 模拟。 <ref name="Deardorff_1970"/>它将涡流粘度建模为：

: <math>\nu_\mathrm{t} 
= C \Delta^2\sqrt{2\bar{S}_{ij}\bar{S}_{ij}} 
= C \Delta^2 \left| \bar{S} \right|
</math>

其中<math>\Delta</math>是网格大小，<math>C</math>是一个常数。

该方法假设小尺度的能量产生和耗散处于平衡状态——即， <math>\epsilon = \Pi</math> .

===== 动态模型（Germano 等） =====

:

Germano 等 <ref name="Germano_1991">{{Cite journal |last=Germano, M. |last2=Piomelli, U. |last3=Moin, P. |author-link3=Parviz Moin |last4=Cabot, W. |year=1991 |title=A dynamic subgrid‐scale eddy viscosity model |url=https://semanticscholar.org/paper/35659a343e691dc3440a1e6c2da0240dad61761a |journal=[[Physics of Fluids|Physics of Fluids A]] |volume=3 |issue=7 |page=1760–1765 |bibcode=1991PhFl....3.1760G |doi=10.1063/1.857955 |s2cid=55719851}}
</ref>使用 Smagorinsky 模型确定了许多研究，每个研究都发现了不同的 Smagorinsky 常数值<math>C</math>针对不同的流量配置。为了试图为 SGS 模型制定一种更通用的方法，Germano 等提出了一个动态 Smagorinsky 模型，它使用了两个过滤器：一个网格 LES 过滤器，表示为<math>\overline{f}</math> ，以及一个测试 LES 滤波器<math>\hat{f}</math>，用于任何湍流的场<math>f</math> 。测试过滤器的尺度大于网格过滤器，并在 LES 表示的已经平滑的场上增加了对湍流场的额外平滑。将测试滤波器应用于 LES 方程（通过将“网格”滤波器应用于 Navier-Stokes 方程获得）会产生一组新的方程，它们的形式相同但 SGS 应力项<math>\tau_{ij} = \overline{u_{i} u_{j}} - \bar{u}_{i} \bar{u}_{j}</math>被<math>T_{ij} = \widehat{\overline{u_{i} u_{j}}} - \hat{\bar{u}}_{i} \hat{\bar{u}}_{j}</math>所替代。 Germano等注意到即使两者都由于存在未解析的尺度而无法精确计算，这两个张量之间仍然存在一个精确的关系。这种关系被称为 Germano 恒等式：

<math>
L_{ij} = T_{ij} - \hat{\tau}_{ij}.
</math>

这里<math> L_{ij} = \widehat{\bar{u}_{i} \bar{u}_{j}} - \widehat{\bar{u}_{i}} \widehat{\bar{u}_{j}}</math>可以直接计算得出，因为它只涉及过滤的速度和测试过滤的操作。同一性的意义在于，如果假设湍流是自相似的，那么网格和测试级别的 SGS 应力具有相同的形式<math>\tau_{ij} - (\tau_{kk}/3)\delta_{ij} = - 2 C \Delta^{2} |\bar{S}_{ij}| \bar{S}_{ij}</math>和<math>T_{ij} - (T_{kk}/3)\delta_{ij} = - 2 C \hat{\Delta}^{2} |\hat{\bar{S}}_{ij}| \hat{\bar{S}}_{ij}</math> ，则 Germano 恒等式提供了一个方程，从中可以得到 Smagorinsky 系数<math>C</math> （不再是“常数”）可能得以推出。

为了做到这一点，需要在原始推导中引入x两个额外的步骤。首先，假设即使<math>C</math>原则上是可变的，变化足够慢，可以从过滤操作中移出：<math> \widehat{C (.)} = C \widehat{(.)} </math> 。其次，既然<math>C</math>是一个标量，Germano 恒等式与一个二阶张量（选择应变张量的速率）联系起来，以将其转换为一个标量方程，从而推出<math>C</math>。 Lilly <ref name="Lilly_1992">{{Cite journal |last=Lilly |first=D. K. |year=1992 |title=A proposed modification of the Germano subgrid-scale closure method |journal=Physics of Fluids A |volume=4 |issue=3 |page=633–636 |bibcode=1992PhFlA...4..633L |doi=10.1063/1.858280}}</ref>发现了一种从张量恒等式中推出 <math>C</math>的更好的方法。他指出，Germano 恒等式需要在空间中的每个点上满足单个量的九个方程（其中只有五个是独立的） <math>C</math>值 。因此对<math>C</math>的推算条件多余了。他提议<math>C</math>通过最小化残差使用最小二乘拟合来确定。由此得：<math>
C = \frac{ L_{ij} m_{ij} }{ m_{kl} m_{kl} }.
</math>

其中<math>
m_{ij} = \alpha_{ij} - \widehat{\beta}_{ij} 
</math>。为简洁起见<math> \alpha_{ij} = - 2 \hat{\Delta}^{2} | \hat{\bar{S}} | \hat{\bar{S}}_{ij} </math>, <math> \beta_{ij} = - 2 \Delta^2 | \bar{S} | \bar{S}_{ij} </math>在 LES 模拟中实施该模型的最初尝试被证明是不成功的。首先，计算出的系数根本不像假设的那样“缓慢变化”，而且变化与任何其他湍流场一样多。其次，计算<math>C</math>可以是积极的，也可以是消极的。后一个事实本身不应被视为缺点，因为使用过滤的 DNS 字段的先验测试表明，本地子网格耗散率<math> - \tau_{ij} \bar{S}_{ij}</math>即使在流体域上的积分始终为正，表示大尺度上的能量净耗散，在湍流场中的积分几乎与正数一样可能为负数。与涡粘度的严格正相反，正值的轻微优势导致观察到的净耗散。这种从小到大的能量“反向散射”确实对应于 Smagorinsky 模型中的负 C 值。然而，发现 Germano-Lilly 公式不能产生稳定的计算。通过在同质方向上平均分子和分母（在流动中存在这样的方向），采用了一种特别的方法

:

: <math>
C = \frac{ 
\left\langle L_{ij} m_{ij} \right\rangle
}{ 
\left\langle m_{kl} m_{kl} \right\rangle
}.
</math>

当平均涉及足够大的统计样本时，计算<math>C</math>是正的（或至少极少负值）稳定的计算是可能的。简单地将负值设置为零（一个称为“剪裁”的过程），无论有或没有平均，也会提高计算稳定性。 Meneveau 提出<ref name="Meneveauetal_1996">{{Cite journal |last=Meneveau, C. |last2=Lund, T. S. |last3=Cabot, W. H. |year=1996 |title=A Lagrangian dynamic subgrid-scale model of turbulence |journal=J. Fluid Mech. |volume=319 |issue=1 |page=353–385 |bibcode=1996JFM...319..353M |doi=10.1017/S0022112096007379 |hdl=2060/19950014634 |s2cid=122183534 |hdl-access=free}}</ref>对具有指数衰减“记忆”的拉格朗日流体轨迹进行平均。这可以应用于缺乏均匀方向的问题，并且如果进行平均的有效时间足够长，但又不足以消除感兴趣的空间不均匀性，那么这可以是稳定的。

Lilly 对 Germano 方法的修改，然后是统计平均或综合去除负粘度区域似乎是临时的，即使它可以“工作”。 Ghosal 等人提出了一种称为“动态定位模型”（DLM）的最小二乘最小化过程的替代公式。 <ref name="Ghosal_1995">{{Cite journal |last=Ghosal |first=S. |last2=Lund |first2=T.S. |last3=Moin |first3=P. |last4=Akselvoll |first4=K. |year=1995 |title=A dynamic localization model for large-eddy simulation of turbulent flows |url=https://archive.org/details/sim_journal-of-fluid-mechanics_1995-03-10_286/page/n234 |journal=Journal of Fluid Mechanics |volume=286 |issue= |page=229–255 |doi=10.1017/S0022112095000711 |s2cid=124586994}}</ref>在这种方法中，首先定义一个量

: <math>
E_{ij} = L_{ij} - T_{ij} + \hat{\tau}_{ij}
</math>

张量<math>\tau_{ij}</math>和<math>T_{ij}</math>替换为适当的 SGS 模型。然后，该张量表示子网格模型未能在每个空间位置符 Germano 恒等式的数量。在Lilly的方法中，<math>C</math>可以脱离帽符算子：

: <math>
 \widehat{C (.)} = C \widehat{(.)}
</math>

使得<math>E_{ij}</math>成为一个<math>C</math>的代数函数，然后通过要求确定<math>E_{ij} E_{ij}</math>被认为是<math>C</math>的函数具有最小的可能值。然而，由于<math>C</math>如此获得的结果与湍流中的任何其他波动量一样可变，原始假设<math>C</math>不能后天证明。在 DLM 方法中，通过不调用从测试过滤操作中删除<math>C</math>的步骤来避免这种不一致。相反，人们通过数量定义了整个流的全局误差：

: <math>
 E [ C ] = \int E_{ij} E_{ij} dV
</math>

其中积分范围在整个流体体积上。这个全局误差<math>E[C(x,y,z,t)]</math>是空间变化函数<math>C(x,y,z,t)</math> 的函数（这里的瞬时<math>t</math>是固定的，只是作为一个参数出现），<math>C</math>的结果需要让误差最小化。这个变分问题的解决方法是<math>C</math>必须满足第二类 Fredholm 积分方程

: <math> 
C (\boldsymbol{x}) = f ( \boldsymbol{x} ) + 
\int K(\boldsymbol{x}, \boldsymbol{y}) C ( \boldsymbol{y} ) d\boldsymbol{y}
</math>

其中函数<math>K(\boldsymbol{x}, \boldsymbol{y})</math>和<math>f ( \boldsymbol{x} )</math>根据解析的项<math>L_{ij},\alpha_{ij},\beta_{ij}</math>定义，因此在每个时间步长和整个流体域的积分范围都是已知的。积分方程通过迭代程序进行数值求解，如果与预处理方案一起使用，则发现收敛通常很快。尽管这种变分方法消除了 Lilly 方法中固有的不一致性，但<math>C(x,y,z,t)</math>从积分方程获得的仍然显示出与负粘度相关的不稳定性。这可以通过坚持来解决<math>E[C]</math>受约束最小化<math>C(x,y,z,t) \geq 0</math> .这导致了一个<math>C</math>的非线性方程：

: <math> 
C (\boldsymbol{x}) = \left[ f ( \boldsymbol{x} ) + 
\int K(\boldsymbol{x}, \boldsymbol{y}) C ( \boldsymbol{y} ) d\boldsymbol{y} \right]_{+}
</math>

这里后缀 + 表示“正数部分”，即， <math> x_{+} = (x + |x|)/2 </math> 。尽管这表面上看起来像“剪裁”，但它不是一个临时方案，而是约束变分问题的真正解决方案。这种 DLM(+) 模型被发现是稳定的，并且对于强制和衰减的各向同性湍流、通道流动和各种其他更复杂的几何形状产生了出色的结果。如果流动恰好具有均匀的方向（让我们说方向 x 和 z），那么可以引入假定 <math> C = C(y,t) </math> 。然后，变分方法立即产生 Lilly 的结果，对均匀方向进行平均，无需对先前结果进行临时修改。

DLM(+) 模型的一个缺点是它没有描述反向散射，而反向散射在分析 DNS 数据时被认为是真实的“事物”。开发了两种方法来解决这个问题。由于 Carati 等人的一种方法。 <ref name="Carati_1995">{{Cite journal |last=Carati |first=D. |last2=Ghosal |first2=S. |last3=Moin |first3=P. |year=1995 |title=On the representation of backscatter in Dynamic Localization models |journal=Physics of Fluids |volume=7 |issue=3 |page=606–616 |doi=10.1063/1.868585}}</ref>类似于朗道的脉动流体力学理论，增加了一个由涨落耗散定理确定的振幅的脉动力。在第二种方法中，有人注意到任何“反向散射”能量出现在解析尺度中，只是以亚网格尺度中的能量为代价。 DLM 可以通过一种简单的方式进行修改，以考虑到这一物理事实，从而在本质上稳定的同时允许反向散射。 在DLM 的这个 k 方程版本中，DLM(k) 替换在 Smagorinsky 涡流粘度模型中的<math> \Delta | \bar{S} | </math>一项，以<math> \sqrt{k} </math>作为适当的速度尺度。确定<math>C</math>的流程保持与“无约束”版本相同，除了张量<math> \alpha_{ij} = - 2 \hat{\Delta} \sqrt{K} \hat{\bar{S}}_{ij} </math>, <math> \beta_{ij} = - 2 \hat{\Delta} \sqrt{k} \bar{S}_{ij} </math>（其中子测试尺度动能 K 与子网格尺度动能 k 的关系为<math> K = k + L_{ii}/2 </math> ）（接着追踪 Germano 恒等式）。为了确定 k，我们现在使用传输方程

: <math> 
\frac{\partial k}{\partial t} + u_{j} \frac{\partial k}{\partial x_{j}} = 
- \tau_{ij} \bar{S}_{ij} - \frac{C_{*}}{\Delta} k^{3/2} + 
\frac{\partial }{\partial x_j} \left( D \Delta \sqrt{k} \frac{\partial k }{\partial x_j} \right) 
+ \nu \frac{\partial^{2} k }{\partial x_j \partial x_j}
</math>

其中<math>\nu</math>是运动粘度，<math>C_{*},D</math>是分别代表动能耗散和扩散的正系数。这些可以按照 DLM(+) 中的约束最小化的动态过程来确定。这种方法虽然比 DLM(+) 实施起来更昂贵，但被发现是稳定的，并且与测试的各种流的实验数据有很好的一致性。此外，DLM(k) 在数学上不可能导致计算不稳定，因为大尺度和 SGS 能量的总和不会因构造而增加。这两种结合反向散射的方法都运作良好。与 DLM(+) 相比，它们产生的模型耗散稍小，性能有所提高。 DLM(k) 模型还产生子网格动能，这可能是一个感兴趣的物理量。这些改进是在模型实现的成本有所增加的情况下实现的。

动态模型起源于 1990 年[[史丹佛大學|斯坦福大学]]湍流研究中心(CTR) 的[https://ctr.stanford.edu/ctr-summer-program 暑期项目] {{Wayback|url=https://ctr.stanford.edu/ctr-summer-program |date=20221029130930 }}。一系列“CTR-Tea”研讨会庆祝了湍流建模这一重要里程碑的[https://ctr.stanford.edu/event-type/ctr-tea 30 周年] {{Wayback|url=https://ctr.stanford.edu/event-type/ctr-tea |date=20221030050059 }}。

== 另见 ==

* [[流体力学]]
* 伽利略不变性——某些类型滤波器的重要属性
* [[雷诺平均纳维－斯托克斯方程|雷诺平均 Navier-Stokes 方程]]
* [[湍流]]

== 延伸阅读 ==

* Heus, T.;&nbsp;van Heerwaarden, C. C.;Jonker, H. J. J.;&nbsp;Pier Siebesma, A.;&nbsp;Axelsen, S. DOI: [https://www.geosci-model-dev.net/3/415/2010/gmd-3-415-2010.html 10.5194/gmd-3-415-2010] {{Wayback|url=https://www.geosci-model-dev.net/3/415/2010/gmd-3-415-2010.html |date=20161012172556 }}. ISSN: 1991-9603.

== 参考资料 ==
<references group="" responsive="1"></references>
[[Category:计算流体力学]]
[[Category:亂流]]
[[Category:流体力学]]
[[Category:流体动力学]]
[[Category:偏微分方程]]