CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Wiki > LES filters

LES filters

From CFD-Wiki

Revision as of 09:15, 3 January 2012 by Peter (Talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

In Large eddy simulation (LES) only the large scale motions of the flow are solved for by filtering out the small and universal eddies. In practical applications of some SGS models, implicit filtering is done by the grid itself and programmers need not worry about the filtering operation. The values of velocity on the grid are the filtered values of velocity. However, for some SGS models, such as the Dynamic subgrid-scale model an explicit filtering step is required to compute the SGS stress tensor. Additionally, in the theoretical analysis of LES, filtering a function is defined as convoluting the function with a filtering kernel, just as is typically done in electrical engineering.

Some of the commonly used filters are defined below. In all cases,  \Delta is the filter width,  G(x) is the filtering kernel in physical space and \widehat{G(k)} is the filtering kernel in Fourier-wavenumber space.

Box filter

The Box filter is the same as the "grid filter" whereby the filter cuts off the values of the function beyond a half filter width away.


G(x) = \frac{1}{\Delta} H\left( \frac{1}{2}\Delta -|x| \right)

where H is the Heaviside function,


\widehat{G(k)} = \frac{\sin \left (\frac{1}{2} k \Delta \right)}{ \frac{1}{2} k \Delta}

Gaussian filter

The Gaussian filter is a normalized Gaussian function. The Fourier transform of a Gaussian function is also a Gaussian, hence the G(x) and \widehat{G(k)} have very similar forms,


G(x) = \left( \frac{6}{\pi \Delta^2} \right)^{1/2} e^{ \frac{-6 x^2}{\Delta^2}}



\widehat{G(k)} = e^{ \frac{-k^2 \Delta^2}{24} }
My wiki