TWiki> Main Web>TWikiUsers>RichardHoward>RichardHowardLESwithQ (2009-12-28, DominiqueLaurence)
The theory behind an LES model adapted for both high and low Reynolds number flows is explained here

## A general LES model

Before attempting to develop any kind of LES model one MUST do a DNS.

The reason for this is that only DNS can demonstrate without any doubt that the numerical method used is capable of resolving local turbulent oscillations. We can define a DNS as a simulation in which the turbulence kinetic energy of the flow is conserved without the introduction of any kind of turbulence modelling ( i.e. the turbulence kinetic energy statistics balance). Typically a CFD code will have been written to conserve momentum explicitly and so a DNS is the only real test that the system is truly capable of capturing the local spatial and time oscillations associated with "natural" turbulence.

It is not important to study different Reynolds numbers as a single DNS is enough to demonstrate the viability of the numerical scheme.

## The objective of an LES

Viewed in a different way a DNS shows us how much oscillation the numerical scheme is capable of dealing with. This means that if a given cell encounters an oscillation that is stronger it must be reduced or "smoothed" to a level that the numercial scheme can cope with. This "smoothing" process is what the LES model must do.

One of the key problems with LES modelling (or indeed any type of turbulence modelling) is that the behaviour of the "smoothing" process changes significantly with the amount of work it has to do.

For example a tiny grid cell very close to a wall doesn't require the same sort of modelling as a large grid cell that covers an entire boundary layer.

## Some features necessary in an LES model

There are several features that are important for the turbulent subgrid viscosity eddy viscosity of an LES model:

• on a no slip surface

• in the centre of a vortex

It should be noted that there is no physical reason for the turbulent subgrid viscosity to always be positive. It is common practice to prevent negative subgrid viscosity from occuring to ensure numerical stability but this has no physical justification.

In order to gain a better understanding of how LES models work it is useful to look at the classic Smagorinsky model.

## How much smoothing do we have to do in an LES?

The answer is: it depends on the flow !!

Next to walls however there is a relatively good non-dimensional property which can be used to quite accurately describe the flow behaviour.

This quantity is where For example near the wall where and the popular empiric profile of Reichardt can also be used.

The main problem with is that it is inherently related to the distance from the wall whereas LES (and indeed any kind of turbulence modelling) should be only related to the local flow behaviour.

A local version can be given by

so

## a local version of The is a useful parameter to analyse LES turbulent flow field.

N.B. For models The old version was developed to make use of this parameter. The main difficulty with the usage of this parameter within a model seems to be its Reynolds number dependence.

Another quantity that has a far more subtle Reynolds number dependence is related to the Smagorinsky turbulent eddy viscosity. where the reference unity values of the friction velocity and and length scale are used to non dimensionalise the parameter.

## Boundary layer turbulence as a paradigm for inhomogeneous turbulence

DNS data can be used to plot and compare it with different powers of From this we can get an idea of how the dependence of should change with .

For example ,  ,  ,  , This sort of way at looking at suggests that we require at least 5 different models for or in another sense the sensitivity of the turbulent viscosity is different at different distances from the wall.

In order to obtain the y dependence we can make use of and the approximate assumption that . In addition we can convert the wall associated dependence to a local flow dependence ,  ,  ,  , Since the dimensions of are the same as : and and also has no dependence, we can use with to provide the necessary dimensions for  ,  ,  ,  , For larger values of a turbulent viscosity associated with high Reynolds numbers needs to be developed.

## High Reynolds number LES

The most common way of dealing with the wall shear stress in high Reynolds number wall flows is to use a power law or log law.

We must remember that we are trying to model part of the momentum equation for the near wall grid and the model is represented as In effect at high Reynolds numbers because the the viscosity contribution becomes negligeable.

The friction velocity can either be obtained from a power law e.g. which can be written as (Werner Wengle, , )

or a log law which can be written as ( , )

In both cases the value of the friction velocity can be obtained as a function of the local grid velocity, wall distance and the molecular viscosity: and this is the value that is included in the momentum equation.

For the power law shown above the relation is so This suggests that at high Reynolds numbers the turbulent viscosity should be modelled as and ## Modelling high Reynolds number oscillations

The wall laws and power laws discussed above demonstrate that the turbulent viscosity needs to have some molecular viscosity dependence at high Reynolds numbers. ,  , Latex rendering error!! dvi file was not created.

Current Tags:
create new tag
, view all tags
Topic revision: r19 - 2009-12-28 - 16:26:19 - DominiqueLaurence
Main Web
22 Mar 2019

Manchester CfdTm
Code_Saturne

ATAAC
KNOO

DESider
FLOMANIA