Modeling of individual coherent structures in wall region of a turbulent boundary layer

Models for individual coherent structures in the wall region of a turbulent boundary layer are proposed. Method of numerical simulations is used to follow the evolution of the structures. It is found that the proposed model does bear many features of coherent structures found in experiments. Project supported by the National Natural Science Foundation of China (Grant No. 19732005) and National Climbing Project.     

Effect of pressure gradient on coherent structures in a turbulent boundary layer

By using the idea of resonant triad of the theory of hydrodynamic stability, the effect of pressure gradient on coherent structures in a turbulent boundary layer is investigated. The favorable pressure gradient suppresses the generation of the coherent structure, while the adverse pressure gradient has the opposite effect. The scale, form, as well as the propagation speed of the coherent structures are different from those with zero pressure gradient. The theoretical results are, in general, in agreement with those found from experiments. From the calculated probability density curve of the circulation differences of the nearly streamwise vortex pairs, it is found that the adverse pressure gradient makes the vortex pair more symmetric. Project supported by the National Natural Science Foundation of China.     

The heat transfer problem in the wall region of a turbulent boundary layer

The problem of heat transfer in the wall region of a turbulent boundary layer has been investigated. The resonant triad in the theory of hydrodynamic stability was used to obtain the velocity field induced by the coherent structures in the wall region of the turbulent boundary layer, while the small scale turbulence was represented by a simple model. By such a new approach of modeling, the 3-D temperature field is calculated, the mean temperature profile in the wall region and the Nusselt number characterizing the heat flux, which was found to be in good agreement with the experimental observations are obtained. The instantaneous temperature field had streaky structures, thus offering a mechanism for their generation found in numerical simulations. Project supported by the National Natural Science Foundation of China (Grant No. 19132011).     

Self-similar turbulent boundary layer in pressure gradient

Self-similar flows in a turbulent boundary layer when the free-stream velocity is specified as a power function of longitudinal coordinate are investigated. The self-similar formulation not only simplifies solving of the problem by reducing the equations of motion to ordinary differential equations but also provides a mean for formulating closure conditions. It is shown that for the class of flows under consideration that depend on three governing parameters the dimensionless mixing length is a function of the normalised distance from the wall and the exponent in the law specifying the free-stream velocity distribution in the outer region and a universal function of local Reynolds number in the wall region, the latter corollary being true even when the skin friction vanishes. In calculations this function is set to be independent of pressure gradient, which gives the results very close to experimental data. There exist four different self-similar flow regimes. Each regime is related to its similarity parameter, one of which is the well-known Clauser equilibrium parameter and the other three are established for the first time. In case of adverse pressure gradient when the exponent lies within certain limits, which depend on Reynolds number, the problem has two solutions with different values of the boundary layer thickness and skin friction, which points out the possibility of hysteresis in near-separating flow. Separation occurs not at the minimum value of the exponent that corresponds to the strongest adverse pressure gradient but at a higher one whose dependence on Reynolds number is calculated in the paper. The results of the theory are in good agreement with experimental data. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Transitional intermittency in a rough wall boundary layer

The transitional intermittency factor distributions were evaluated from the wall friction records acquired in the rough wall boundary layer under free stream turbulence. Important characteristics of by-pass transition were derived from the factor distributions. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Near-separating,self-similar turbulent boundary layer

The self-similar flow in a turbulent boundary layer, which is in a state close to separation as a result of the effect of adverse pressure gradient, is investigated. Such a boundary layer has a triple-deck asymptotic structure. Between outer and near-wall regions above the logarithmic sublayer, i. e. the constant-stress layer, an intermediate region — the gradient sublayer — is formed, where the shear stress varies linearly due to adverse longitudinal pressure gradient. In the external part of the gradient sublayer, the velocity profile obeys the square-root law. The velocity profile obtained from the solution for the outer region satisfies a slip condition on the wall. The slip value decreases as the similarity parameter increases and vanishes at the value of Ω = 0.0911, which corresponds to separation, here δ^* is the displacement thickness, and U and U′ are the free-stream velocity and its derivative with respect to the longitudinal coordinate. In this case, the exponent m in the law specifying the free-stream self-similar velocity distribution increases, with separation occurring not at the minimal value of m = −1/3, which corresponds to the strongest adverse pressure gradient, but at the value m = 0.228. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The problem of the turbulent boundary layer

Summary Existing measurements of low-speed turbulent surface friction on a flat plate, in the absence of pressure gradient and roughness, are shown to be consistent with a simple analysis based on functional similarity in the velocity profile. In particular, the fully developed turbulent boundary layer is found to be unique within the accuracy of the experimental data, with uniqueness defined as the existence of a definite correspondence between local friction coefficient and momentum thickness Reynolds number. The relationships known as the law of the wall and the velocity defect law are found to describe the turbulent velocity profiles accurately for a considerable range of Reynolds numbers, and an effort is made to clarify the physical significance of these formulae. Finally, the proper definition of a length Reynolds number is discussed in terms of the asymptotic local properties of the ideal boundary layer, and numerical values for ideal mean and local friction coefficients are tabulated against Reynolds numbers based on momentum thickness and on distance from the leading edge.

---

Zusammenfassung Es wird gezeigt, dass vorhandene Messungen der turbulenten Wandschubspannung an der glatten ebenen Platte in inkompressibler Strömung ohne Druckgradient durch eine einfache Berechnung in Übereinstimmung gebracht werden können. Die Rechnung beruht auf einer funktionellen Ähnlichkeit der Geschwindigkeitsverteilung. Es wird im besonderen gefunden, dass die vollentwickelte turbulente Grenzschicht innerhalb der Messgenauigkeit einem eindeutigen Zusammenhang zwischen dem örtlichen Reibungskoeffizienten und der Reynoldsschen Zahl, bezogen auf die Impulsdicke, folgt. Die Beziehungen, die als Wandgesetz und Mittengesetz bekannt sind, beschreiben die Geschwindigkeitsverteilung genau innerhalb eines erheblichen Bereiches Reynoldsscher Zahlen, und es wird versucht, den physikalischen Inhalt dieser Gesetzmässigkeiten zu vertiefen. Abschliessend wird eine zweckmässige Definition der auf Plattenlänge bezogenen Reynoldsschen Zahl diskutiert, die auf dem asymptotischen örtlichen Zustand der idealen Grenzschicht beruht. Rechenwerte der idealen, mittleren und örtlichen Reibungskoeffizienten, bezogen auf beide obigen Definitionen der Reynoldsschen Zahl, werden tabelliert.

     

Direct simulation of the turbulent boundary layer on a plate

A numerical method for the integration of three-dimensional Navier–Stokes equations for compressible fluid as applied to direct numerical simulation is proposed. By way of example, the boundary layer on a plate is simulated. The computations were carried out for Re_θ = 1500. The computational grid consisted of a half billion nodes. The flow region includes the laminar, transitional, and turbulent zones. The numerically obtained distributions of average velocity, friction, and pulsations are compared with experimental data and available numerical solutions.     

Three-dimensional turbulent boundary layer on wind turbine blades

The momentum integral technique for predicting the boundary–layer growth in three–dimensional flow has been extended to include the entrainment equation as the closure model. Special attention has been devoted to those terms in the differential equations that change the boundary–layer structure from that of cvasi two–dimensional steady flow, at outboard locations, to a tree–dimensional flow pattern. It is concluded that the stall is delayed due to the boundary–layer reattachment at inboard sections in conjunction with the onset of a spanwise vortex like structure. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical simulation of the turbulent boundary layer equations via a Runge-Kutta algorithm

This paper addresses a hybrid computational procedure for the step-by-step calculation of momentum transfer in turbulent boundary layer flows along flat plates. The proposed procedure relies on a modified method of lines wherein transversal discretizations are carried out by a “control volume” being infinitesimal in the streamwise direction and finite in the transversal direction of the fluid flow. Using mixing length theory and coarse intervals in the transversal direction, the resulting system of ordinary differential equations of first order may be readily integrated on a personal computer utilizing a fourth-order Runge-Kutta algorithm. In general, a maximum number of sixteen lines is necessary at the trailing edge of the flat plate for a typical calculation. As a consequence, computing time and storage for each run were very small when compared to other finite-difference methods. Furthermore, to validate the hybrid procedure involving the method of lines and control volumes (MOLCV), comparisons with experimental data have been done in terms of both velocity distributions and local skin friction coefficients. For all cases tested, the proposed methodology predicts the growth of the boundary layer of gases correctly.     

Dynamical structure of the turbulent boundary layer on rough surface

The experiments on structure of turbulent boundary layer on the plane rough wall without pressure gradient are presented. Sand roughness of the wall is considered. Measurements are carried out using Time-Resolved PIV technique in planes parallel and perpendicular to the wall. The results on rough wall are compared with the base case of boundary layer on smooth wall. Hairpin vortices have been detected. Topology and typical size of those structures substantially differ in the cases in question. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The boundary layer on a wall nearly parallel to a magnetic field

This paper calculates the structure of the boundary layer which occurs at the wall of a nonconducting duct in fully-developed laminar flow of a conducting fluid when the wall is inclined at a small angle to the uniform transverse magnetic field and when the Hartmann numberM is high. When is large compared withM ^–1/2 the layer splits into a Hartmann layer on the inclined wall and a weak shear layer emanating from the duct corner along field lines. These two layers cause equal reductions in the total flow rate.

---

Zusammenfassung Diese Arbeit behandelt die Struktur der Grenzschicht, die an der Wand eines nicht leitfähigen Kanals entsteht die eine vollentwickelte Laminarströmung einer leitfähigen Flüssigkeit führt, wobei die Wand einen kleinen Winkel mit dem gleichförmigen transversalen Magnetfeld bildet und die Hartmannsche ZahlM groß ist. Für große Werte von im Vergleich zuM ^–1/2 teilt sich die Grenzschicht in einen Hartmannsche Schicht an der geneigten Wand und eine schwache Scherschicht, die von der Ecke des Kanals entlang der Feldlinien ausgeht. Die beiden Schichten verursachen die gleichen Verminderungen in der Durchflußmenge.

     

Formulation of wall boundary conditions in turbulent flow computations on unstructured meshes

Features of the formulation and numerical implementation of wall boundary conditions in turbulent flow computations on unstructured meshes are discussed. A method is proposed for implementing weak wall boundary conditions for a finite-volume discretization of the Reynolds-averaged Navier-Stokes equations on unstructured meshes. The capabilities of the approach are demonstrated in several gasdynamic simulations in comparison with the method of near-wall functions. The influence of the near-wall resolution on the accuracy of the computations is analyzed, and the grid dependence of the solution is compared in the case of the near-wall function method and weak boundary conditions.     

A Lagrangian simulation of particle deposition in a turbulent boundary layer in the presence of thermophoresis

Knowledge of particle deposition in turbulent flows is often required in engineering situations. Examples include fouling of turbine blades, plate-out in nuclear reactors and soot deposition. Thus it is important for numerical simulations to be able to predict particle deposition. Particle deposition is often principally determined by the forces acting on the particles in the boundary layer. The particle tracking facility in the CFD code uses the eddy lifetime model to simulate turbulent particle dispersion, no specific boundary layer being modelled. The particle tracking code has been modified to include a boundary layer. The non-dimensional yplus, y⁺, distance of the particle from the wall is determined and then values for the fluid velocity, fluctuating fluid velocity and eddy lifetime appropriate for a turbulent boundary layer used. Predictions including the boundary layer have been compared against experimental data for particle deposition in turbulent pipe flow. The results giving much better agreement. Many engineering problems also involve heat transfer and hence temperature gradients. Thermophoresis is a phenomena by which small particles experience a force in the opposite direction to the temperature gradient. Thus particles will tend to deposit on cold walls and be repulsed by hot walls. The effect of thermophoresis on the deposition of particles can be significant. The modifications of the particle tracking facility have been extended to include the effect of thermophoresis. A preliminary test case involving the deposition of particles in a heated pipe has been simulated. Comparison with experimental data from an extensive experimental programme undertaken at ISPRA, known as STORM (Simplified Tests on Resuspension Mechanisms), has been made.     

An asymptotic approach for shock-wave/transpired turbulent boundary layer interactions

Summary In this work perturbation techniques are used to study the problem of the interaction between a shock wave and a transpired turbulent boundary layer at transonic speeds. In the case considered here, the Mach number is assumed to be high enough for the sonic line to penetrate deep into the boundary layer so that it ends close to the wall. The flow region is divided into a region of strong interaction and regions of weak interaction. For the regions of weak interaction, upstream and downstream of the shock, a classical two-deck structure is assumed to hold for the boundary layer. Solutions chosen for these regions must account for the effects of blowing or suction. The strong interaction region on the other hand is shown to consist of three decks. A detailed analysis of the whole flow field is carried out and solutions valid in the double limit as Reynolds number tends to infinity and Mach number tends to one are proposed. Solutions of adjacent layers are shown to match so providing a smooth solution for the entire flow region. The analysis yields solutions for the pressure and skin-friction profiles.

---

Zusammenfassung In dieser Veröffentlichung wird mit Methoden der Störungsrechnung das Problem der Wechselwirkung zwischen einer Schockwelle und einer turbulenten Grenzschicht mit Blas- und Saugeffekten beim Übergang zu Überschallgeschwindigkeiten untersucht. In dem hier betrachteten Fall ist die Machzahl so hoch, daß die Schockwellenfront tief in die Grenzschicht eindringt und nahe bei der Wand endet. Das Strömungsgebiet wird in ein Gebiet starker Wechselwirkung und in ein Gebiet schwacher Wechselwirkung eingeteilt. Für die Gebiete schwacher Wechselwirkung (vor und hinter der Schockwellenfront) wird angenommen, daß die Grenzschicht eine klassische zwei-Deck Struktur besitzt. Lösungen, die für diese Gebiete gewählt werden, müssen die Blas- oder Saugeffekte erklären. Es wird gezeigt, daß das Gebiet starker Wechselwirkung dagegen aus drei Decks besteht. Eine ausführliche Berechnung des gesamten Strömungsfeldes wird durchgeführt, und es werden Lösungen angegeben, die gültig sind für den doppelten Grenzfall, daß die Reynoldszahl gegen unendlich und die Machzahl gegen eins strebt. Es zeigt sich, daß die Lösungen für benachbarte Schichten aneinander passen, und daß sich damit eine stetige Lösung für das gesamte Strömungsgebiet ergibt. Die Lösungen für die Druck- und Wandreibungsprofile werden bestimmt.

     

Laws of the wall for velocity and temperature in supersonic turbulent boundary layers

Scaling laws for velocity and temperature profiles in the near-wall region of sub- and supersonic turbulent boundary layers have been developed, which allow us to represent velocity and temperature profiles in compressible gas stream in terms of those in an incompressible boundary layer. They are obtained as asymptotic expansions of the solutions to the Reynolds equations in a small parameter — the Mach number based on the friction velocity and gas enthalpy on the wall. The leading term of the expansion for velocity corresponds to known Van Driest's formula. However, the obtained solution contains additional terms of order unity, which explains the contradiction between Van Driest's formula and experimental data. The law of the wall for temperature, which has been formulated for the first time, has an analogous structure. Besides the von Kármán constant and the turbulent Prandtl number in the logarithmic region, known for incompressible flow, the obtained relations contain three new universal constants, which do not depend on gas molecular properties and the specific heat ratio. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Scaling laws for turbulent boundary layer flows in favorable and adverse pressure gradient

An analytical theory is proposed to describe incompressible plane and axisymmetric turbulent boundary layer flows in favorable and adverse pressure gradients for near-equilibrium conditions. Scaling laws for the mean velocity, the Reynolds stress components, and the skin friction have been established. A universal friction law makes it possible to represent the skin friction distributions corresponding to different Reynolds numbers and pressure gradients in terms of a function of one variable. The theory is based on general physical assumptions and does not involve any special hypotheses on the nature of the turbulent motion. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)