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.     

Inflow boundary condition for DNS of turbulent boundary layers on supersonic blunt cones

For direct numerical simulation (DNS) of turbulent boundary layers, gen- eration of an appropriate inflow condition needs to be considered. This paper proposes a method, with which the inflow condition for spatial-mode DNS of turbulent boundary layers on supersonic blunt cones with different Mach numbers, Reynolds numbers and wall temperature conditions can be generated. This is based only on a given instant flow field obtained by a temporal-mode DNS of a turbulent boundary layer on a flat plate. Effectiveness of the method is shown in three typical examples by comparing the results with those obtained by other methods.     

DNS,LES and RANS of turbulent heat transfer in boundary layer with suddenly changing wall thermal conditions

The objectives of this study are to investigate a thermal field in a turbulent boundary layer with suddenly changing wall thermal conditions by means of direct numerical simulation (DNS), and to evaluate predictions of a turbulence model in such a thermal field, in which DNS of spatially developing boundary layers with heat transfer can be conducted using the generation of turbulent inflow data as a method. In this study, two types of wall thermal condition are investigated using DNS and predicted by large eddy simulation (LES) and Reynolds-averaged Navier–Stokes equation simulation (RANS). In the first case, the velocity boundary layer only develops in the entrance of simulation, and the flat plate is heated from the halfway point, i.e., the adiabatic wall condition is adopted in the entrance, and the entrance region of thermal field in turbulence is simulated. Then, the thermal boundary layer develops along a constant temperature wall followed by adiabatic wall. In the second case, velocity and thermal boundary layers simultaneously develop, and the wall thermal condition is changed from a constant temperature to an adiabatic wall in the downstream region. DNS results clearly show the statistics and structure of turbulent heat transfer in a constant temperature wall followed by an adiabatic wall. In the first case, the entrance region of thermal field in turbulence can be also observed. Thus, both the development and the entrance regions in thermal fields can be explored, and the effects upstream of the thermal field on the adiabatic region are investigated. On the other hand, evaluations of predictions by LES and RANS are conducted using DNS results. The predictions of both LES and RANS almost agree with the DNS results in both cases, but the predicted temperature variances near the wall by RANS give different results as compared with DNS. This is because the dissipation rate of temperature variance is difficult to predict by the present RANS, which is found by the evaluation using DNS results.     

Lagrangian model on the turbulent motion of small solid particle in turbulent boundary layer flows

.Intr0ductionSurfaceerosionofmaterialbysolid-particleimpactisanimportantprobleminmultiphaseflowindustriaIdevicesandthecharacteristicsoftheparticIe'smotioninaturbulentboundarylayerflowisthebaseofthestudyofthematerialsurfaceerosion.Manycalculationmodelshave…     

A theoretical mode for the instability of turbulent boundary layer over compliant surface

A theoretical model for the instability of turbulent boundary layer over compliant surfaces is described. The investigation of instability is carried out from a time-asymptotic space-time perspective that classifies instabilities as either convective or absolute. Results are compared against experimental observations of surface waves on elastic and viscoelastic compliant layers.     

Bubble migration in a turbulent boundary layer

The wall void peaking distribution observed in an upward turbulent bubbly boundary layer along a flat plate is generated by bubbles that move towards the plate, come into contact with the wall and then slide along it. This transverse ‘migration’ has been studied using flow visualization, high speed video and particle tracking techniques to measure the trajectories of mono-disperse air bubbles at very low void fractions. Investigations have been performed at four Reynolds numbers in the range 280 < Re_θ < 3000, covering both the laminar and turbulent regimes, with mono-disperse bubbles of mean equivalent diameter between 2 mm and 6 mm. Lagrangian statistics calculated from hundreds of trajectories show that the migration only occurs in the turbulent regime and for bubble diameters below some critical value: 3.5 mm < d_eqcrit < 4 mm. Above this size (We > 3), the interface deformation is such that bubbles do not remain at the wall, even when they are released at the surface. Also, bubble migration is shown to be non-systematic, to have a non-deterministic character in the sense that trajectories differ significantly, to increase with Reynolds number and to take place on a short time scale. A series of experiments with isolated bubbles demonstrates that these results are not influenced by bubble–bubble interactions and confirm that two-way coupling in the flow is limited. Flow visualizations show that the migration originates with the capture of bubbles inside the large turbulent structures of the boundary layer (‘bulges’). The bubbles begin to move towards the wall as they cross these structures, and the point at which they reach the wall is strongly correlated with the position of the deep ‘valleys’ which separate the turbulent ‘bulges’. The analysis of the mean Lagrangian trajectories of migrating bubbles confirms these observations. Firstly, the average time of migration calculated from these trajectories coincides with the mean transit time of the bubbles across the structures. Secondly, once the trajectories have been scaled by this transit time and the boundary layer thickness δ, they all have the same form in the region y/δ < 0.4, independent of the Reynolds number.     

Transport of particles in an atmospheric turbulent boundary layer

A program incorporating the parallel code of large eddy simulation (LES) and particle transportation model is developed to simulate the motion of particles in an atmospheric turbulent boundary layer (ATBL). A model of particles of 100-micrometer order coupling with large scale ATBL is proposed. Two typical cases are studied, one focuses on the evolution of particle profile in the ATBL and the landing displacement of particles, whereas the other on the motion of particle stream.The English text was polished by Yunming Chen.     

A kind of fast changing coherent structure in a turbulent boundary layer

Coherent structures of a turbulent boundary layer were investigated by hydrogen bubble method. A kind of fast changing structure was observed. That is a spot in which all the hydrogen bubbles vanish much faster than in other regions. This investigation verified that dark-spot is formed by a strong sweep from outer layer. Inside a dark-spot the local instantaneous flow speed might be four times of its neighboring high-speed streaks. Comparing with the low/high speed streaks, both dark-spot and the vortical structures around it are changing very fast. Around dark-spot intensive shear layers are formed and indications of the generation of small-scale structures could be observed. The project supported by the National Natural Science Foundation of China (No. 19672004) and the National Climbing Project     

Higher-order and length-scale statistics from DNS of a decelerated planar wall-bounded turbulent flow

A DNS database is presented to document third- and fourth-order moments and their budgets for fully developed plane-channel flow and for strained plane-channel flow. The effect of straining has a similar effect on statistics, such as the skewness and flatness of velocities, as in an adverse-pressure-gradient (APG) boundary layer. In addition to higher-order statistics, some modeling implications are also described, including issues related to decomposition of the velocity–pressure gradient correlations, the assumption of dissipation isotropy and a fourth-order turbulence time scale. An analysis of two-point correlations along the inhomogeneous direction is made to include the effects of APG straining. This reveals the advantage of a one-point length-scale based on wall-normal velocity fluctuations as an alternative to the traditionally used length-scale based on turbulence kinetic energy. The present study should prove to be useful for turbulence modelers in need of data to develop and assess higher-than-second-order Reynolds-averaged Navier–Stokes closure models.     

Visualization of flow structures in a turbulent boundary layer using a new technique

In this paper, an experimental investigation on the flow structures in a turbulent boundary layer employing a special laser light sheet-Hydrogen bubble flow visualization technique is described. It is observed that the high/low speed streaks are directly related to the hairpin or horseshoe-like vortices. This observation can give a better understanding of the physical mechanism in the turbulent boundary layer. Fluid Mechanic Institute, BUAA     

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.     

Multiple distortion of a supersonic turbulent boundary layer

Two experiments were performed to study the response of a supersonic turbulent boundary layer to successive distortions. In the first experiment (Case 1), the flow passed over a forward-facing ramp formed by 20° compression corner followed by a 20° expansion corner located about 4_o downstream, where _o is the incoming boundary layer thickness. In the second experiment (Case 2), the forward-facing ramp was constructed of curved compression and expansion surfaces with the same turning angles and total step height as in Case 1. The radii of curvature for the compression and expansion surfaces were equal to 12_o. In both experiments, the flow relaxation was observed over a distance equal to 12_o. In this relaxation region, the mean and turbulent flow behavior of the boundary layer was measured. The mean velocity profile was found to be altered by the distortion. Recovery of the profile began near the wall and occurred rapidly, but in the outer part of the boundary layer, recovery proceeded slowly. Turbulence measurements revealed a dramatic reduction in the turbulence shear stress and a progressively decaying streamwise Reynolds stress profile.     

Some new results of the observation on the coherent structures of turbulent boundary layer in flow with adverse pressure gradient

By a suitable manipulation of hydrogen bubble generation, some new results were obtained: (1) The long-streaks are generated along the interfaces between low and high-speed streaks. The long-streaks are generally stretching and are moving faster than its neighboring high-speed streaks. The hydrogen bubbles in long-streaks have longer life. (2) The stream-wise vortices are also generated along the interfaces. The project is supported by the National Natural Science Foundation of China.     

A mechanism for excitation of coherent structures in wall region of a turbulent boundary layer

Introduction Themechanismforthegenerationofcoherentstructuresinthewallregionofaturbulent boundarylayerhasalwaysbeeninconcernandinvestigated.AccordingtoTsujimotoand Miyake[1],thecharacteristicsofturbulenceinthewallregionweremainlydeterminedbythe generationandevolutionofcoherentstructures,notbythesmall_scaleturbulence.However, excitationsfromregionofy >60werefoundtobenecessary,otherwisethewallregionwould degeneratetolaminarflow.Therefore,theinvestigationofthemechanismthathowcoherent structuresi…     

The structures of turbulent boundary layer in the vicinity of separation

Large coherent structures of turbulent boundary layer in the vicinity of separation were observed in a water channel by the hydrogen bubble method. Motion pictures of the de views were taken. The features of the instantaneous velocity profiles, the large transverse and streamwise vortices were discussed.     

Instability waves in the wall region of a turbulent boundary layer on a flat plate

Coherent structures and the bursting phenomena in the wall region of a turbulent boundary layer play a very important role in determining the characteristics of the boundary layer. Yet the nature and the origin of the coherent structures are unclear until now. In this paper, nonlinear stability calculations for the wall region of a turbulent boundary layer have been made. It was found that there do exist instability waves which may be responsible for the coherent structures. The project is supported by the National Natural Science Foundation of China.     

Numerical prediction of turbulent flow structure generated by a synthetic cross‐jet into a turbulent boundary layer

A detailed numerical study using large‐eddy simulation (LES) and unsteady Reynolds‐averaged Navier–Stokes (URANS) was undertaken to investigate physical processes that are engendered in the injection of a circular synthetic (zero‐net mass flux) jet in a zero pressure gradient turbulent boundary layer. A complementary study was carried out and was verified by comparisons with the available experimental data that were obtained at corresponding conditions with the aim of achieving an improved understanding of fluid dynamics of the studied processes. The computations were conducted by OpenFOAM C++, and the physical realism of the incoming turbulent boundary layer was secured by employing random field generation algorithm. The cavity was computed with a sinusoidal transpiration boundary condition on its floor. The results from URANS computation and LES were compared and described qualitatively and quantitatively. There is a particular interest for acquiring the turbulent structures from the present numerical data. The numerical methods can capture vortical structures including a hairpin (primary) vortex and secondary structures. However, the present computations confirmed that URANS and LES are capable of predicting current flow field with a more detailed structure presented by LES data as expected. Copyright © 2011 John Wiley & Sons, Ltd.     

Computation of shock wave/turbulent boundary layer interactions using a two-equation model with compressibility corrections

In this paper computational results for two different types of shock wave / turbulent boundary layer interaction flows are presented. It is shown that upstream effects of the shock induced separation cannot be reproduced by Wilcox's (1991) k--model, whereas downstream of the interaction, predictions of pressure distribution and skin friction are acceptable. The inclusion of the compressible part of the dissipation rate and the pressure dilatation in the model has noticeable, but not dramatic effects on wall pressure and skin friction in the selected flow cases.     

Bursting frequency in turbulent boundary layers

Employing laser Doppler anemometry and VITA techniques, the bursting frequency in turbulent boundary layers has been measured over the Reynolds-number range 320 to 1470. The result indicates that the mean and non-dimensional bursting frequency scaled with the variables appropriate for the wall region was constant and independent of Reynoids number. When the same data are plotted using the outer variables of boundary layer to normalize the bursting frequency, the non-dimensional frequency increases as the Reynolds number increases. This is in agreement with the results of Blackwelder et al. (1983) who used hot wire anemometry and VITA technique. The project is supported by the National Natural Science Foundation of China