Filtering in the context of large eddy simulation (LES) is a mathematical operation intended to remove a range of small scales from the solution to the Navier-Stokes equations. Because the principal difficulty in simulating turbulent flows comes from the wide range of length and time scales, this operation makes turbulent flow simulation cheaper by reducing the range of scales that must be resolved. The LES filter operation is low-pass, meaning it filters out the scales associated with high frequencies.
Homogeneous filters[]
The same DNS velocity field filtered using a
box filter and
Δ = L/32
The same DNS velocity field filtered using a
box filter and
Δ = L/16
Definition in physical space[]
The low-pass filtering operation used in LES can be applied to a spatial and temporal field, for example
. The LES filter operation may be spatial, temporal, or both. The filtered field, denoted with a bar, is defined as:[1][2]
![\overline {\phi ({\boldsymbol {x}},t)}=\displaystyle {\int _{{-\infty }}^{{\infty }}}\int _{{-\infty }}^{{\infty }}\phi ({\boldsymbol {r}},t^{{\prime }})G({\boldsymbol {x}}-{\boldsymbol {r}},t-t^{{\prime }})dt^{{\prime }}d{\boldsymbol {r}},](https://wikimedia.org/api/rest_v1/media/math/render/svg/5277270a6ba30dfc857a8bdbf6c1a6f1d21e1048)
where
is a convolution kernel unique to the filter type used. This can be written as a convolution operation:
![\overline{\phi} = G \star \phi .](https://wikimedia.org/api/rest_v1/media/math/render/svg/9a0b7f599c3da3b4507ad0fcbd74b14e8f94afc7)
The filter kernel
uses cutoff length and time scales, denoted
and
respectively. Scales smaller than these are eliminated from
Using this definition, any field
may be split up into a filtered and sub-filtered (denoted with a prime) portion, as
![\phi = \bar{\phi} + \phi^{\prime} .](https://wikimedia.org/api/rest_v1/media/math/render/svg/786d4c15c126a71ddef29b0ce657c806d3d1c96b)
This can also be written as a convolution operation,
![\phi ^{{\prime }}=\left(1-G\right)\star \phi .](https://wikimedia.org/api/rest_v1/media/math/render/svg/cb88be48ba9eed4a2418be7fa72e49a3d90a5dc6)
Definition in spectral space[]
The filtering operation removes scales associated with high frequencies, and the operation can accordingly be interpreted in Fourier space. For a scalar field
the Fourier transform of
is
a function of
the spatial wave number, and
the temporal frequency.
can be filtered by the corresponding Fourier transform of the filter kernel, denoted
![\overline {{\hat {\phi }}}({\boldsymbol {k}},\omega )={\hat {\phi }}({\boldsymbol {k}},\omega ){\hat {G}}({\boldsymbol {k}},\omega )](https://wikimedia.org/api/rest_v1/media/math/render/svg/29f2adf21a406327b83ee28ffe32ef3d07cc9d2e)
or,
![\overline {{\hat {\phi }}}={\hat {G}}{\hat {\phi }}.](https://wikimedia.org/api/rest_v1/media/math/render/svg/2f498eced7a60bb1b5125d277e0856dc0e84db90)
The filter width
has an associated cutoff wave number
and the temporal filter width
also has an associated cutoff frequency
The unfiltered portion of
is:
![{\hat {\phi ^{{\prime }}}}=(1-{\hat {G}}){\hat {\phi }}.](https://wikimedia.org/api/rest_v1/media/math/render/svg/e962307a54b415a9b999e88d19f4c1f876ee4e8c)
The spectral interpretation of the filtering operation is essential to the filtering operation in large eddy simulation, as the spectra of turbulent flows is central to LES subgrid-scale models, which reconstruct the effect of the sub-filter scales (the highest frequencies). One of the challenges in subgrid modeling is to effectively mimic the cascade of kinetic energy from low to high frequencies. This makes the spectral properties of the implemented LES filter very important to subgrid modeling efforts.
Homogeneous filter properties[]
Homogeneous LES filters must satisfy the following set of properties when applied to the Navier-Stokes equations.[1]
- 1. Conservation of constants
- The value of a filtered constant must be equal to the constant,
![\overline {a}=a,](https://wikimedia.org/api/rest_v1/media/math/render/svg/cb3f0d35352a5bc493761f9e5fbe4a4f1c3aa97d)
- which implies,
![\int _{{-\infty }}^{{\infty }}\int _{{-\infty }}^{{\infty }}G({\boldsymbol {\xi }},t^{{\prime }})d^{3}{\boldsymbol {\xi }}dt^{{\prime }}=1.](https://wikimedia.org/api/rest_v1/media/math/render/svg/1ceeb36ac3cc3740d1be791144f56616263490e4)
- 2. Linearity
![\overline {\phi +\psi }=\overline {\phi }+\overline {\psi }.](https://wikimedia.org/api/rest_v1/media/math/render/svg/fee280ea53c87db62c5dfeb3b7b75e0968e5c190)
- 3. Commutation with derivatives
![\overline {{\frac {\partial \phi }{\partial s}}}={\frac {\partial \overline {\phi }}{\partial s}},\qquad s={\boldsymbol {x}},t.](https://wikimedia.org/api/rest_v1/media/math/render/svg/17d71d469735716652b6e1b911246caa06d43847)
- If notation is introduced for operator commutation
for two arbitrary operators
and
, where
![[f,g]\phi =f\circ g(\phi )-g\circ f(\phi )=f(g(\phi ))-g(f(\phi )),](https://wikimedia.org/api/rest_v1/media/math/render/svg/c16f02f51de3e6ea23348440769390f4ee0e1ea5)
- then this third property can be expressed as
![\left[G\star ,{\frac {\partial }{\partial s}}\right]=0.](https://wikimedia.org/api/rest_v1/media/math/render/svg/3ffc45070eaa4ca1d242e99c812ab5cb6d39c772)
Filters satisfying these properties are generally not Reynolds operators, meaning, first:
![{\displaystyle {\begin{array}{rcl}{\overline {\overline {\phi }}}&\neq &{\overline {\phi }},\\G\star G\star \phi =G^{2}\star \phi &\neq &G\star \phi ,\end{array}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e845ffe7c3e16133454537552ec0742194b2a380)
and second,
![\overline {\phi ^{{\prime }}}=G\star (1-G)\star \phi \neq 0.](https://wikimedia.org/api/rest_v1/media/math/render/svg/5445610a8264bea36323629699a76ee70f6f989d)
Inhomogeneous filters[]
Implementations of filtering operations for all but the simplest flows are inhomogeneous filter operations. This means that the flow either has non-periodic boundaries, causing problems with certain types of filters, or has a non-constant filter width
, or both. This prevents the filter from commuting with derivatives, and the commutation operation leads to several additional error terms:
![{\begin{array}{rcl}\left[{\frac {\partial }{\partial {\boldsymbol {x}}}},G\star \right]\phi &=&{\frac {\partial }{\partial {\boldsymbol {x}}}}\left(G\star \phi \right)-G\star {\frac {\partial \phi }{\partial {\boldsymbol {x}}}}\\&=&{\frac {\partial }{\partial {\boldsymbol {x}}}}\int _{{\Omega }}G({\boldsymbol {x}}-{\boldsymbol {r}},\Delta ({\boldsymbol {x}},t))\phi ({\boldsymbol {r}},t)d{\boldsymbol {r}}-G\star {\frac {\partial \phi }{\partial {\boldsymbol {x}}}}\\&=&\left({\frac {\partial G}{\partial \Delta }}\star \phi \right){\frac {\partial \Delta }{\partial x}}+\int _{{d\Omega }}G(x-r,\Delta (x,t))\phi (r,t){\boldsymbol {n}}dS\end{array}},](https://wikimedia.org/api/rest_v1/media/math/render/svg/8e57745aa7084e0d0022df94ab931e87a98445e8)
where
is the vector normal to the surface of the boundary
and
[1]
The two terms both appear due to inhomogeneities. The first is due to the spatial variation in the filter size
while the second is due to the domain boundary. Similarly, the commutation of the filter
with the temporal derivative leads to an error term resulting from temporal variation in the filter size,
![\left[{\frac {\partial }{\partial t}},G\star \right]=\left({\frac {\partial G}{\partial \Delta }}\star \phi \right){\frac {\partial \Delta }{\partial t}}.](https://wikimedia.org/api/rest_v1/media/math/render/svg/924c69bb805ddb80ca5ad23979fb9d5c77c964ac)
Several filter operations which eliminate or minimize these error terms have been proposed.[citation needed]
Classic large eddy simulation filters[]
![[icon]](//upload.wikimedia.org/wikipedia/commons/thumb/1/1c/Wiki_letter_w_cropped.svg/20px-Wiki_letter_w_cropped.svg.png) | This section needs expansion with: Proper alignment of plots. You can help by . (January 2020) |
Turbulent energy spectrum and effect of filtering operations
[3]
There are three filters ordinarily used for spatial filtering in large eddy simulation. The definition of
and
and a discussion of important properties, is given.[2]
Box filter[]
Box filter in physical and spectral space
The filter kernel in physical space is given by:
![G({\boldsymbol {x}}-{\boldsymbol {r}})={\begin{cases}{\frac {1}{\Delta }},&{\text{if}}\left|{\boldsymbol {x}}-{\boldsymbol {r}}\right|\leq {\frac {\Delta }{2}},\\0,&{\text{otherwise}}.\end{cases}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/465ce596d57ef747917fe6008ac33c74c67560c8)
The filter kernel in spectral space is given by:
![{\hat {G}}({\boldsymbol {k}})={\frac {\sin {({\frac {1}{2}}k\Delta )}}{{\frac {1}{2}}k\Delta }}.](https://wikimedia.org/api/rest_v1/media/math/render/svg/cb165eadcbe9ebbfeab268e9faa64142fcaa3e2f)
Gaussian filter in physical and spectral space
Gaussian filter[]
The filter kernel in physical space is given by:
![G({\boldsymbol {x}}-{\boldsymbol {r}})=\left({\frac {6}{\pi \Delta ^{{2}}}}\right)^{{{\frac {1}{2}}}}\exp {\left(-{\frac {6({\boldsymbol {x-r}})^{2}}{\Delta ^{2}}}\right)}.](https://wikimedia.org/api/rest_v1/media/math/render/svg/5ea85de7269b1a65df8af44d656df9b4fc5952d6)
The filter kernel in spectral space is given by:
![{\hat {G}}({\boldsymbol {k}})=\exp {\left(-{\frac {{\boldsymbol {k}}^{2}\Delta ^{2}}{24}}\right)}.](https://wikimedia.org/api/rest_v1/media/math/render/svg/d53c50f37e45028da8a6a72a7f3f9dc26c3558f6)
Sharp spectral filter in physical and spectral space
Sharp spectral filter[]
The filter kernel in physical space is given by:
![G({\boldsymbol {x}}-{\boldsymbol {r}})={\frac {\sin {(\pi ({\boldsymbol {x-r}})/\Delta )}}{\pi ({\boldsymbol {x-r}})}}.](https://wikimedia.org/api/rest_v1/media/math/render/svg/afe726113d59e09d4318716ee40528cb03ece794)
The filter kernel in spectral space is given by:
![{\hat {G}}({\boldsymbol {k}})=H\left(k_{c}-\left|k\right|\right),\qquad k_{c}={\frac {\pi }{\Delta }}.](https://wikimedia.org/api/rest_v1/media/math/render/svg/45acf3b6955557603ca79cb84512ee14d4c8e6e0)
See also[]
References[]