Simulation and Modelling of Turbulent Trailing-Edge Flow

Computations of turbulent trailing-edge flow have been carried out at a Reynolds number of 1000 (based on the free-stream quantities and the trailing-edge thickness) using an unsteady 3D Reynolds-Averaged Navier–Stokes (URANS) code, in which two-equation (k–ε) turbulence models with various low-Re near wall treatments were implemented. Results from a direct numerical simulation (DNS) of the same flow are available for comparison and assessment of the turbulence models used in the URANS code. Two-dimensional URANS calculations are carried out with turbulence mean properties from the DNS used at the inlet; the inflow boundary-layer thickness is 6.42 times the trailing-edge thickness, close to typical turbine blade flow applications. Many of the key flow features observed in DNS are also predicted by the modelling; the flow oscillates in a similar way to that found in bluff-body flow with a von Kármán vortex street produced downstream. The recirculation bubble predicted by unsteady RANS has a similar shape to DNS, but with a length only half that of the DNS. It is found that the unsteadiness plays an important role in the near wake, comparable to the modelled turbulence, but that far downstream the modelled turbulence dominates. A spectral analysis applied to the force coefficient in the wall normal direction shows that a Strouhal number based on the trailing-edge thickness is 0.23, approximately twice that observed in DNS. To assess the modelling approximations, an a priori analysis has been applied using DNS data for the key individual terms in the turbulence model equations. A possible refinement to account for pressure transport is discussed.

[1]  Paul A. Durbin,et al.  A Procedure for Using DNS Databases , 1998 .

[2]  T. G. Thomas,et al.  Structure and energetics of a turbulent trailing edge flow , 2003 .

[3]  Klaus Bremhorst,et al.  A Modified Form of the k-ε Model for Predicting Wall Turbulence , 1981 .

[4]  P. Durbin SEPARATED FLOW COMPUTATIONS WITH THE K-E-V2 MODEL , 1995 .

[5]  P. Spalart Direct simulation of a turbulent boundary layer up to Rθ = 1410 , 1988, Journal of Fluid Mechanics.

[6]  T. Shih,et al.  New time scale based k-epsilon model for near-wall turbulence , 1993 .

[7]  M. Alam,et al.  Two-Equation Turbulence Modelling of a Transitional Separation Bubble , 2000 .

[8]  K. Hannemann,et al.  Numerical simulation of the absolutely and convectively unstable wake , 1989 .

[9]  A. Roshko On the drag and shedding frequency of two-dimensional bluff bodies , 1954 .

[10]  P. R. Voke,et al.  Turbulent Simulation of a Flat Plate Boundary Layer and Near Wake , 1997 .

[11]  Brian Launder,et al.  Contribution towards a Reynolds-stress closure for low-Reynolds-number turbulence , 1976, Journal of Fluid Mechanics.

[12]  W. Dawes,et al.  Predicting turbulent flow in a staggered tube bundle , 1999 .

[13]  Z. Yang,et al.  A Galilean and tensorial invariant k-epsilon model for near wall turbulence , 1993 .

[14]  William N. Dawes,et al.  The practical application of solution-adaption to the numerical simulation of complex turbomachinery problems , 1992 .

[15]  J. C. Rotta,et al.  Turbulent boundary layers in incompressible flow , 1962 .

[16]  Neil D. Sandham,et al.  Direct Numerical Simulation of Turbulent Flow over a Rectangular Trailing Edge , 2001 .

[17]  V. C. Patel,et al.  Turbulence models for near-wall and low Reynolds number flows - A review , 1985 .