Main CfdTm Forum Saturne Aster ATAAC TWikiEditRaw editAttachPrint version

TWiki> Main Web>TWikiUsers>RichardHoward>RichardHowardLESwithQ (2009-12-28, DominiqueLaurence)

TWiki> Main Web>TWikiUsers>RichardHoward>RichardHowardLESwithQ (2009-12-28, DominiqueLaurence)

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.

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.

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.

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

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.

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.

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

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.

Topic revision: r19 - 2009-12-28 - 16:26:19 - DominiqueLaurence

Code_Saturne

KNOO

FLOMANIA

Photo Wall

Upcoming Events

Add Event

Copyright © by the contributing authors. Unless noted otherwise, all material on this web site is the property of the contributing authors.