Subgrid-Scale Stress Modelling Based on the Square of the Velocity Gradient Tensor

A new subgrid scale model is proposed for Large Eddy Simulations in complex geometries. This model which is based on the square of the velocity gradient tensor accounts for the effects of both the strain and the rotation rate of the smallest resolved turbulent fluctuations. Moreover it recovers the proper y ³ near-wall scaling for the eddy viscosity without requiring dynamic procedure. It is also shown from a periodic turbulent pipe flow computation that the model can handle transition.     

KOLMOGOROV BEHAVIOR OF NEAR-WALL TURBULENCE AND ITS APPLICATION IN TURBULENCE MODELING

The near-wall behavior of turbulence is re-examined in a way different from that proposed by Hanjalic and Launder¹ and followers^2,3,4,5. It is shown that at a certain distance from the wall, all energetic large eddies will reduce to Kolmogorov eddies (the smallest eddies in turbulence). All the important wall parameters, such as friction velocity, viscous length scale, and mean strain rate at the wall, are characterised by Kolmogorov microscales. According t o this Kolmogorov behavior of near-wall turbulence, the turbulence quantities, such as turbulent kinetic energy, dissipation rate, etc. at the location where the large eddies become “Kolmogorov” eddies, can be estimated by using both direct numerical simulation (DNS) data and asymptotic analysis of near-wall turbulence. This information will provide useful boundary conditions for the turbulent transport equations. As a n example, the concept is incorporated in the standard κ - εmodel which is then applied t o channel and boundary layer flows. Using appropriate boundary conditions (based on Kolmogorov behaviour of near-wall turbulence), there is no need for any wall-modification to the κ - ε equations (including model constants). Results compare very well with the DNS and experimental data.     

An experimental contribution to near-wall measurements by means of a special laser Doppler anemometry technique

 A new experimental technique for the investigation of near-wall turbulence using laser Doppler anemometry is presented, which allows an accurate measurement of the flow field very close to the wall, with good resolution and a high data rate. Such a technique is tested in a fully developed turbulent flow (with Reynolds numbers between 4,300 and 67,000) by carrying out a careful statistical analysis of the streamwise and wall-normal velocity components within the near-wall region, at distances from the wall ranging from approximately y ⁺ = 1 to y ⁺ = 100. The velocity profiles, Reynolds stresses and higher-order moments of the two-dimensional boundary layer are presented. The results, which are in agreement with the most recent data in the literature, testify the validity of the proposed experimental solution. Moreover, the accuracy of the results allows the friction velocity to be calculated as the intercept at the wall of the best linear fit of the total stress profile; in this way, an unambiguous examination of the normalized statistics is possible. Received: 17 April 2001 / Accepted: 15 August 2001     

Experimental Investigation of the Log-Law for an Adverse Pressure Gradient Turbulent Boundary Layer Flow at Re θ = 10000

We present an experimental investigation of a turbulent boundary layer flow at a significant adverse pressure gradient at Reynolds number Re _θ ?=?10000 using large field PIV. The testcase is designed to start from a zero pressure gradient flow at Re _θ ?=?8000 with a distinct log-law region following a slowly rising adverse pressure gradient. This allows to reveal a breakdown of the log-law under the effect of the adverse pressure gradient. The region described by the log-law is progressively reduced in terms of y ^?+? and then joins into a modified log-law which gives a good fit to the data up to at least y/δ ₉₉?≈?0.2. The scaling in the overlap region is demonstrated using the mean velocity slope diagnostic function, enabled due to the high quality of the PIV data. Locally, the velocity profile is measured down to the wall using long-range microscopic PIV with particle tracking velocimetry to determine the wall shear stress directly in the adverse pressure gradient region.     

DNS of a spatially developing turbulent boundary layer with passive scalar transport

A direct numerical simulation (DNS) of a spatially developing turbulent boundary layer over a flat plate under zero pressure gradient (ZPG) has been carried out. The evolution of several passive scalars with both isoscalar and isoflux wall boundary condition are computed during the simulation. The Navier–Stokes equations as well as the scalar transport equation are solved using a fully spectral method. The highest Reynolds number based on the free-stream velocity U_∞ and momentum thickness θ is Re_θ=830, and the molecular Prandtl numbers are 0.2, 0.71 and 2. To the authors’ knowledge, this Reynolds number is to date the highest with such a variety of scalars. A large number of turbulence statistics for both flow and scalar fields are obtained and compared when possible to existing experimental and numerical simulations at comparable Reynolds number. The main focus of the present paper is on the statistical behaviour of the scalars in the outer region of the boundary layer, distinctly different from the channel-flow simulations. Agreements as well as discrepancies are discussed while the influence of the molecular Prandtl number and wall boundary conditions is also highlighted. A Pr scaling for various quantities is proposed in outer scalings. In addition, spanwise two-point correlation and instantaneous fields are employed to investigate the near-wall streak spacing and the coherence between the velocity and the scalar fields. Probability density functions (PDF) and joint probability density functions (JPDF) are shown to identify the intermittency both near the wall and in the outer region of the boundary layer. The present simulation data will be available online for the research community.     

Calculation of complex near-wall turbulent flows with a low-Reynolds-number k-ϵ model

An improved low-Reynolds-number k-? model has been formulated and tested against a range of DNS (direct numerical simulation) and experimental data for channel and complex shear layer flows. The model utilizes a new form of damping function adopted to account for both wall proximity effects and viscosity influences and a more flexible damping argument based on the gradient of the turbulent kinetic energy on the wall. Additionally, the extra production of the inhomogeneous part of the viscous dissipation near a wall has been added to the dissipation equation with significantly improved results. The proposed model was successfully applied to the calculation of a range of wall shear layers in zero, adverse and favourable pressure gradients as well as backward-facing-step separated flows.     

BOUNDARY TREATMENT AND AN EFFICIENT PRESSURE ALGORITHM FOR INTERNAL TURBULENT FLOWS USING THE PDF METHOD

The generalized Langevin model, which is used to model the motion of stochastic particles in the velocity–composition joint probability density function (PDF) method for reacting turbulent flows, has been extended to incorporate solid wall effects. Anisotropy of Reynolds stresses in the near-wall region has been addressed. Numerical experiments have been performed to demonstrate that the forces in the near-wall region of a turbulent flow cause the stochastic particles approachi ng a solid wall to reverse their direction of motion normal to the wall and thereby, leave the near-wall layer. This new boundary treatment has subsequently been implemented in a full-scale problem to prove its validity. The test problem considered here is that of an isothermal, non-reacting turbulent flow in a two-dimensional channel with plug inflow and a fixed back-pressure. An efficient pressure correction method, developed in the spirit of the PISO algorithm, has been implemented. The pressure correction strategy is easy to implement and is completely consistent with the time- marching scheme used for the solution of the Lagrangian momentum equations. The results show remarkable agreement with both k–ϵ and algebraic Reynolds stress model calculations for the primary velocity. The secondary flow velocity and the turbulent moments are in better agreement with the algebraic Reynolds stress model predictions than the k– ϵ predictions. © 1997 by John Wiley & Sons, Ltd.     

Features of the turbulent wall pressure field on a long towed cylinder

Turbulent wall pressure fluctuation correlation functions were measured in water on a towed cylindrical model of length 129.8 m and diameter 3.8 cm for steady speeds ranging from 6.2 to 15.5 m/s. The drag on the model was measured with a strut-mounted load cell to provide estimates of the momentum thickness and friction velocity that are used for scaling the correlation functions. Very high momentum thickness Reynolds numbers Re_θ were achieved, and varied from 4.8 × 10⁵ to 1.1 × 10⁶. The ratio of boundary layer thickness to cylinder radius was approximately 24, which is an order of magnitude greater than previous laboratory investigations. The ratio of momentum thickness to viscous length scale is significantly greater than for flat plate cases at comparable Re_θ. A similarity scaling is shown to be more effective than outer or inner boundary layer scalings for collapsing the correlation functions. Comparisons with the early streamwise and transverse correlation measurements of Willmarth and Yang are favorable, and show consistent trends of a more rapid loss of correlated energy for cylindrical turbulent boundary layers than for flat plate cases. Convection velocities are also presented and shown to collapse well with separation scaled on outer variables. A simple model that relates the peak of the correlation function to the average coherence levels is shown to be valid for spatial separations less than the approximate momentum thickness.     

A method for estimating wall friction in turbulent wall-bounded flows

We describe a simple method for estimating turbulent boundary layer wall friction using the fit of measured velocity data to a boundary layer model profile that extends the logarithmic profile all the way to the wall. Two models for the boundary layer profile are examined, the power-series interpolation scheme of Spalding and the Musker profile which is based on the eddy viscosity concept. The performance of the method is quantified using recent experimental data in zero pressure gradient flat-plate turbulent boundary layers, and favorable pressure gradient turbulent boundary layers in a pipe, for which independent measurements of wall shear are also available. Between the two model profiles tested, the Musker profile performs much better than the Spalding profile. Results show that the new procedure can provide highly accurate estimates of wall shear with a mean error of about 0.5% in friction velocity, or 1% in shear stress, an accuracy that is comparable to that from independent direct measurements of wall shear stress. An important advantage of the method is its ability to provide accurate estimates of wall shear not only based on many data points in a velocity profile but also very sparse data points in the velocity profile, including only a single data point such as that originating from a near-wall probe.     

Numerical prediction of turbulent boundary layer noise from a sharp-edged flat plate

An efficient hybrid uncorrelated wall plane waves–boundary element method (UWPW-BEM) technique is proposed to predict the flow-induced noise from a structure in low Mach number turbulent flow. Reynolds-averaged Navier-Stokes equations are used to estimate the turbulent boundary layer parameters such as convective velocity, boundary layer thickness, and wall shear stress over the surface of the structure. The spectrum of the wall pressure fluctuations is evaluated from the turbulent boundary layer parameters and by using semi-empirical models from literature. The wall pressure field underneath the turbulent boundary layer is synthesized by realizations of uncorrelated wall plane waves (UWPW). An acoustic BEM solver is then employed to compute the acoustic pressure scattered by the structure from the synthesized wall pressure field. Finally, the acoustic response of the structure in turbulent flow is obtained as an ensemble average of the acoustic pressures due to all realizations of uncorrelated plane waves. To demonstrate the hybrid UWPW-BEM approach, the self-noise generated by a flat plate in turbulent flow with Reynolds number based on chord Re_c = 4.9 × 10⁵ is predicted. The results are compared with those obtained from a large eddy simulation (LES)-BEM technique as well as with experimental data from literature.     

A New Low-Reynolds-Number One-Equation Model of Turbulence

In this study, we propose a new Low-Reynolds-Number (LRN)one-equation model, which is derived from an LRN two-equation(k-ε) model. The derivation of the transport equation, in principle, is based on the assumption that the turbulent structure parameter remains constant. However, the relation for the turbulent structure parameter a ₁(=|− |/k) is modified to account for near-wall turbulence. As a result, the present one-equation model contains a term which takes the near-wall limiting behavior explicitly into account. Thus, the present model provides the correct wall-limiting behavior of turbulence in the vicinity of the wall and can be applied to the analysis of heat transfer. The validity of the present model is tested in channel flows, boundary layer flows with and without pressure gradient, plane wall jet, and flow with separation and reattachment. The calculated results showed good agreement with the direct numerical simulation (DNS) and experimental data. This revised version was published online in July 2006 with corrections to the Cover Date.     

Near-wall turbulent characteristics along very long thin circular cylinders

Understanding the salient physics within the turbulent boundary layer of towed thin cylinders is paramount to the Navy sonar array communities. However, the required long array length to achieve wide acoustic aperture creates unique and consistent flow characteristics that suggest simplified tangential forcing expressions suitable for design purposes. One well-known fact is that the majority of the array surface experiences very thick turbulent boundary layers (TBL) and large Reynolds numbers. The resultant statistics are most commonly dependent on the inner and outer length scales. Herein, we resolve the near-wall TBL structure under those flow conditions by large-eddy simulation. The turbulent mean-flow statistics showed near-wall consistency using only inner scaling. But both inner and outer variables were found necessary to properly scale the turbulent fluctuations. An expression for the tangential wall-friction coefficient (C_t) indicates two distinct flow regimes as characterized by the near-wall turbulent flow structure. The respective parameters appear independent of the outer length scale. Thickening (or thinning) the cylinder near their common threshold (defined by a radius-based Reynolds number) transitioned the turbulent character between the two regimes.     

An experimental study of a two-dimensional plane turbulent wall jet

 Laser-Doppler measurements were conducted in a plane turbulent wall jet at a Reynolds number based on inlet velocity, Re ₀, of 9600. The initial development as well as the fully developed flow was studied. Special attention was given to the near-wall region, including the use of small measuring volumes and the application of specific near-wall data corrections, so that wall shear stresses were determined directly from the mean velocity gradient at the wall using only data below y ⁺=4. It was possible to resolve the inner peak in the streamwise turbulence intensity as well as the inner (negative) peak in the shear stress. Limiting values of (u′)⁺ and uv ⁺ were determined. Turbulence data from the outer region of the flow were compared to earlier hot wire measurements and large differences in the normal turbulence intensity and the shear stress were found. These differences can be attributed to high turbulence intensity effects on the hot-wires. Received: 17 October 1996 / Accepted: 8 December 1997     

Near-wall flow in a three-dimensional boundary layer on the endwall of a 30° bend

 Turbulence measurements are reported on the three-dimensional turbulent boundary layer along the centerline of the flat endwall in a 30° bend. Profiles of mean velocities and Reynolds stresses were obtained down to y ⁺≈2 for the mean flow and y ⁺≈8 for the turbulent stresses. Mean velocity data collapsed well on a simple law-of-the-wall based on the magnitude of the resultant velocity. The turbulence intensity and turbulent shear stress magnitude both increased with increased three-dimensionality. The ratio of these two quantities, the a ₁ structure parameter, decreased in the central regions of the boundary layer and showed profile similarity for y ⁺<50. The shear stress vector angle lagged behind the velocity gradient vector angle in the outer region of the boundary layer, however there was an indication that the shear stress vector tends to lead the velocity gradient vector close to the wall. Received: 16 July 1996/Accepted: 14 July 1997     

Large eddy simulation of compressible channel flow

Astract The present study is a contribution to the analysis of wall-bounded compressible flows, including a special focus on wall modeling for compressible turbulent boundary layer in a plane channel. large eddy simulation (LES) of fully developed isothermal channel flows at Re = 3,000 and Re = 4,880 with a sufficient mesh refinement at the wall are carried out in the Mach number range 0.3 ≤ M ≤ 3 for two different source term formulations: first the classical extension of the incompressible configuration by Coleman et al. (J. Fluid Mech. 305:159–183, 1995), second a formulation presently derived to model both streamwise pressure drop and streamwise internal energy loss in a spatially developed compressible channel flow. It is shown that the second formulation is consistent with the spatial problem and yields a much stronger cooling effect at the wall than the classical formulation. Based on the present LES data bank, compressibility and low Reynolds number effects are analysed in terms of coherent structure and statistics. A study of the universality of the structure of the turbulence in non-hypersonic compressible boundary layers (M≤5) is performed in reference to Bradshaw (Annu. Rev. Fluid. Mech. 9:33–54, 1977). An improvement of the van Driest transformation is proposed; it accounts for both density and viscosity changes in the wall layer. Consistently, a new integral wall scaling (y ^c+) which accounts for strong temperature gradients at the wall is developed for the present non-adiabatic compressible flow. The modification of the strong Reynolds analogy proposed by Huang et al. (J. Fluid Mech. 305:185–218, 1995) to model the correlation between velocity and temperature for non-adiabatic wall layers is assessed on the basis of a Crocco–Busemann relation specific to channel flow. The key role of the mixing turbulent Prandtl number Pr _m is pointed out. Results show very good agreement for both source formulations although each of them involve a very different amount of energy transfer at the wall. The present work was performed within the framework of the French–German research initiative “large eddy simulation of complex flows’ (UR 507). The computing resources were provided by IDRIS-France. The authors gratefully acknowledge the financial support from the Centre National de la Recherche Scientifique (CNRS), the Centre d’été Mathématique de Recherche Avancée en Calcul Scientifique (CEMRACS) and the Direction Générale de l’Armement (DGA/D4S).     

Turbulent separated reattached flow in a two-dimensional curved-wall diffuser

A turbulent separation-reattachment flow in a two-dimensional asymmetrical curved-wall diffuser is studied by a two-dimensional laser doppler velocimeter. The turbulent boundary layer separates on the lower curved wall under strong pressure gradient and then reattaches on a parallel channel. At the inlet of the diffuser, Reynolds number based on the diffuser height is 1.2×10⁵ and the velocity is 25.2m/s. The results of experiments are presented and analyzed in new defined streamline-aligned coordinates. The experiment shows that after Transitory Detachment Reynolds shear stress is negative in the near-wall backflow region. Their characteristics are approximately the same as in simple turbulent shear layers near the maximum Reynolds shear stress. A scale is formed using the maximum Reynolds shear stresses. It is found that a Reynolds shear stress similarity exists from separation to reattachment and the Schofield-Perry velocity law exists in the forward shear flow. Both profiles are used in the experimental work that leads to the design of a new eddy-viscosity model. The length scale is taken from that developed by Schofield and Perry. The composite velocity scale is formed by the maximum Reynolds shear stress and the Schofield-Perry velocity scale as well as the edge velocity of the boundary layer. The results of these experiments are presented in this paper.     

A multiscale approach to hybrid RANS/LES wall modeling within a high-order discontinuous Galerkin scheme using function enrichment

We present a novel approach to hybrid Reynolds-averaged Navier-Stokes (RANS)/ large eddy simulation (LES) wall modeling based on function enrichment, which overcomes the common problem of the RANS-LES transition and enables coarse meshes near the boundary. While the concept of function enrichment as an efficient discretization technique for turbulent boundary layers has been proposed in an earlier article by Krank & Wall (A new approach to wall modeling in LES of incompressible flow via function enrichment. J Comput Phys. 2016;316:94-116), the contribution of this work is a rigorous derivation of a new multiscale turbulence modeling approach and a corresponding discontinuous Galerkin discretization scheme. In the near-wall area, the Navier-Stokes equations are explicitly solved for an LES and a RANS component in one single equation. This is done by providing the Galerkin method with an independent set of shape functions for each of these two methods; the standard high-order polynomial basis resolves turbulent eddies, where the mesh is sufficiently fine and the enrichment automatically computes the ensemble-averaged flow if the LES mesh is too coarse. As a result of the derivation, the RANS model is applied solely to the RANS degrees of freedom, which effectively prevents the typical issue of a log-layer mismatch in attached boundary layers. As the full Navier-Stokes equations are solved in the boundary layer, spatial refinement gradually yields wall-resolved LES with exact boundary conditions. Numerical tests show the outstanding characteristics of the wall model regarding grid independence, superiority compared to equilibrium wall models in separated flows, and achieve a speed-up by two orders of magnitude compared to wall-resolved LES.     

Laminarisation and re-transition of a turbulent boundary layer subjected to favourable pressure gradient

 Experimental results are reported for the response of an initially turbulent boundary layer (Re _θ≈1700) to a favourable pressure gradient with a peak value of K≡(−υ/ρU ³ _E ) dp/dx equal to 4.4×10^-6. In the near-wall region of the boundary layer (y/δ<0.1) the turbulence intensity u′ scales roughly with the free-stream velocity U _E until close to the location where K is a maximum whereas in the outer region u′ remains essentially frozen. Once the pressure gradient is relaxed, the turbulence level increases throughout the boundary layer until K falls to zero when the near wall u′ levels show a significant decrease. The intermittency γ is the clearest indicator of a fundamental change in the turbulence structure: once K exceeds 3×10^-6, the value of γ in the immediate vicinity of the wall γ _s falls rapidly from unity, reaches zero at the location where K again falls below 3×10^-6 and then rises back to unity. Although γ is practically zero throughout the boundary layer in the vicinity of γ _s =0, the turbulence level remains high. The explanation for what appears to be a contradiction is that the turbulent frequencies are too low to induce turbulent mixing. The mean velocity profile changes shape abruptly where K exceeds 3×10^-6. Values for the skin friction coefficient, based upon hot-film measurements, peak at the same location as K and fall to a minimum close to the location where K drops back to zero. Received: 28 January 1998/Accepted: 8 April 1998     

Transitionally rough zero pressure gradient turbulent boundary layers

Near-wall measurements are performed to study the effects of surface roughness and viscous shear stresses on the transitionally rough regime (5 < k ⁺ < 70) of a zero pressure gradient turbulent boundary layer. The x-dependence is known from the eleven consecutive measurements in the streamwise direction, which allows for the computation of the streamwise gradients in the boundary layer equations. Thus, the skin friction is computed from the integrated boundary layer equation with errors of 3 and 5% for smooth and rough, respectively. It is found that roughness destroys the viscous layer near the wall, thus, reducing the contribution of the viscous stress in the wall region. As a result, the contribution in the wall shear stress due to form drag increases, while the viscous stress decreases. This yields Reynolds number invariance in the skin friction as k ⁺ increases into the fully rough regime. Furthermore, the roughness at the wall reduces the high peak of the streamwise component of the Reynolds stress in the near-wall region. However, for the Reynolds wall-normal and shear stress components, its contribution is not significantly altered for sand grain roughness.