Main CfdTm Forum Saturne Aster ATAAC TWikiEditRaw editAttachPrint version

TWiki> Main Web>HybridModel (2010-08-13, FlavienBillard)

TWiki> Main Web>HybridModel (2010-08-13, FlavienBillard)

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.

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.

(4) |

(5) |

(6) |

(7) |

(8) |

(9) |

(10) |

(11) |

(12) |

(13) |

(14) |

(15) |

Latex rendering error!! dvi file was not created.

Topic revision: r11 - 2010-08-13 - 21:45:28 - FlavienBillard

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.