# The Hybrid Approach

In many hybrid approaches the main challenge resides in coupling the two different velocity fields arising from RANS (statistically averaged) and LES (filtered) . This is often done by applying a matching criteria, i.e. assuming the same turbulent viscosity at the interface, the same kinetic energy or dissipation etc. This poses a problem since these values represent totally different properties of different velocity fields. Instead of a single one velocity field to couple RANS and LES, the model presented here allows an overlap of both fields, with the RANS model driving the near wall RANS velocity field, without damping in any way the dissipation of resolved fluctuations. Many sub-grid models assume that the flow contains an inertial sub range and hence the sub-grid motions can be assumed to be isotropic. This is true only if the grid is small enough for the anisotropy introduced by the mean shear to be neglected. At high Reynolds numbers, the refinement of the grid becomes too costly, therefore restricting the LES method to low Reynolds numbers flows. As the solid boundary is approached, the mean shear becomes high enough to introduce anisotropy across a range of diminishing scales. It is then necessary for the model to represent at the same time subgrid-scale contributions to the mean shear stress and isotropic dissipation effects.

## Modelling

The instantaneous velocity can be decomposed as

where is the averaged velocity and is the fluctuating component.

Schumann (4) proposed to split the residual stress tensor into a "locally isotropic" part and an "inhomogeneous" part. The isotropic part is proportional to the fluctuating strain and does not affect the mean flow equations but determines the rate of energy dissipation. The inhomogeneous part is proportional to the mean strain and controls the shear stress and mean velocity profile:

where denotes ensemble averaging of the filtered equations. The viscosities and are based on fluctuating and mean strains. The isotropic part of the residual stress tensor has a zero time mean value. By refining the grid the residual stresses must tend to zero, therefore the inhomogeneous part must have a grid dependence parameter in the turbulent viscosity . Schumann (4) used a mixing length model for with the length scale computed as , where is a constant that is difficult to prescribe for all types of flows. (4) and (5) tried to derive a theoretical value for the constant but they were forced to introduced corrective constants to agree with a range of experiments. (6) used the same principle of splitting the residual stress but in their mixing length model, they use the spanwise size of the cell as the length scale. They argue that for the near wall region in a channel flow, the important structures are streaks that are finely spaced on the spanwise direction. Therefore a coarse resolution in the spanwise direction would lead to larger eddies and a thicker viscous sublayer.

Sullivan et al. (7) developed a similar approach for planetary boundary layer flows but chose to match the Monin-Obukhov similarity theory (8). Bagget (9) used a similar approach to compare two hybrid models, one "Schumann-like" and one "DES-like" but found excessive streamwise fluctuations leading to streaks that were much too large.

In the context of hybrid LES-RANS, a blending function, , can be used to introduce a smooth transition between the resolved and the ensemble averaged turbulence parts. In the present study the total residual stress is written as:

 (1)

In this way the averaged stress would be:

which is just the RANS stress,and the total shear stress would be . It is therefore necessary that the blending function tends to one in the region where is resolved correctly and to zero in the region near the wall where the shear stress is under resolved due to the coarse grid. The total rate of transfer of energy from the filtered motions to the residual scales is given by (assuming that (10))

which shows how the RANS viscosity contributes to dissipation in association with the mean velocity only, i.e. the resolved turbulent stresses are free to develop independently from the RANS viscosity.

## The turbulent viscosities models.

For the isotropic viscosity , (4) used a model based on the sub-grid energy. (6) used the standard (12) model based on the fluctuating strain. Here the later approach is used:

with . In the frame of unstructured codes, the filter width is taken as twice the cell volume ().

In this study, the elliptic relaxation model $\varphi-f$ of (11) is used to calculate the RANS viscosity. This model solves for the ratio used in the turbulent viscosity as:

 (2)

where . For the channel flow calculations presented here, the choice of the RANS models hardly makes any difference but the elliptic relaxation method has been shown to perform well on separating and impinging flows and the aim of the present case is only to show the robustness of the coupling even with a sophisticated RANS model.

The blending function has been parametrised by the ratio of the turbulent length scale to the filter width:

 (3)

Here are empirical constants. These values were chosen to match the shear stress profile based on channel flow results at with DNS data. When using the model (equation (2)), the wall distance is not a desirable parameter and the blending function can be formulated using . The blending function has been devised to connect the two length scales smoothly so its value is close to zero near the wall and unity far from it. Similar functions have been used in other hybrid approaches (See (13), (14) or (15)). Although the function in equation (3) is totally empirical, it has been tested for a range of Reynolds numbers and grids and gave satisfactory results (not presented here). The function allows a higher contribution from the LES part as the grid is refined. Different coefficients have been used in the optimisation of the blending function, but the results are not greatly affected but nevertheless are always better than standard LES on the same mesh. In equation ((1)) the averaged velocity has been calculated as a running average with an averaging window of about 10 times the eddy turnover time. Although it is possible also to use plane averaging in the case of the channel flow, this was not done in order to keep the formulation applicable for 3D flows where no plane averaging is possible.

## References

 (4)

 (5)

 (6)

 (7)

 (8)

 (9)

 (10)

 (11)

 (12)

 (13)

 (14)

 (15)

Latex rendering error!! dvi file was not created.

Current Tags:
create new tag
, view all tags
Topic revision: r11 - 2010-08-13 - 21:45:28 - FlavienBillard
Main Web
17 Dec 2017

Manchester CfdTm
Code_Saturne

ATAAC
KNOO

DESider
FLOMANIA