Showing papers on "Turbulence published in 1989"

XFOIL: An Analysis and Design System for Low Reynolds Number Airfoils

Abstract: Calculation procedures for viscous/inviscid analysis and mixed-inverse design of subcritical airfoils are presented. An inviscid linear-vorticity panel method with a Karman-Tsien compressiblity correction is developed for direct and mixed-inverse modes. Source distributions superimposed on the airfoil and wake permit modeling of viscous layer influence on the potential flow. A two-equation lagged dissipation integral method is used to represent the viscous layers. Both laminar and turbulent layers are treated, with an e 9-type amplification formulation determinining the transition point. The boundary layer and transition equations are solved simultaneously with the inviscid flowfield by a global Newton method. The procedure is especially suitable for rapid analysis of low Reynolds number airfoil flows with transitional separation bubbles. Surface pressure distributions and entire polars are calculated and compared with experimental data. Design procedure examples are also presented.

Parameterization of Orography-Induced Turbulence in a Mesobeta--Scale Model

Abstract: The possibility of extending existing techniques for turbulence parameterization in the planetary boundary layer to attitude, orography-induced turbulence events is examined. Starting from a well-tested scheme, we show that it is possible to generalize the specification method of the length scales, with no deterioration of the scheme performance in the boundary layer. The new scheme is implemented in a two-dimensional version of a limited-area, numerical model used for the simulation of mesobeta-scale atmospheric flows. Three well-known cases of orographically induced turbulence are studied. The comparison with observations and former studies shows a satisfactory behavior of the new scheme.

Bifurcation theory of poloidal rotation in tokamaks: A model for L-H transition.

Abstract: It is shown that the poloidal momentum balance equation in tokamaks has bifurcated solutions. The poloidal flow velocity ${U}_{p}$ can suddenly become more positive when the ion collisionality decreases. The corresponding radial electric field ${E}_{r}$ becomes more negative and hence suppresses the turbulent fluctuations. Thus, plasma confinement is improved. The theory is employed to explain the L-H transition observed in tokamaks.

Coherent structure of the convective boundary layer derived from large-eddy simulations

Abstract: Turbulence in the convective boundary layer (CBL) uniformly heated from below and topped by a layer of uniformly stratified fluid is investigated for zero mean horizontal flow using large-eddy simulations (LES). The Rayleigh number is effectively infinite, the Froude number of the stable layer is 0.09 and the surface roughness height relative to the height of the convective layer is varied between 10−6 and 10−2. The LES uses a finite-difference method to integrate the three-dimensional grid-volume-averaged Navier–Stokes equations for a Boussinesq fluid. Subgrid-scale (SGS) fluxes are determined from algebraically approximated second-order closure (SOC) transport equations for which all essential coefficients are determined from the inertial-range theory. The surface boundary condition uses the Monin–Obukhov relationships. A radiation boundary condition at the top of the computational domain prevents spurious reflections of gravity waves. The simulation uses 160 × 160 × 48 grid cells. In the asymptotic state, the results in terms of vertical mean profiles of turbulence statistics generally agree very well with results available from laboratory and atmospheric field experiments. We found less agreement with respect to horizontal velocity fluctuations, pressure fluctuations and dissipation rates, which previous investigations tend to overestimate. Horizontal spectra exhibit an inertial subrange. The entrainment heat flux at the top of the CBL is carried by cold updraughts and warm downdraughts in the form of wisps at scales comparable with the height of the boundary layer. Plots of instantaneous flow fields show a spoke pattern in the lower quarter of the CBL which feeds large-scale updraughts penetrating into the stable layer aloft. The spoke pattern has also been found in a few previous investigations. Small-scale plumes near the surface and remote from strong updraughts do not merge together but decay while rising through large-scale downdraughts. The structure of updraughts and downdraughts is identified by three-dimensional correlation functions and conditionally averaged fields. The mean circulation extends vertically over the whole boundary layer. We find that updraughts are composed of quasi-steady large-scale plumes together with transient rising thermals which grow in size by lateral entrainment. The skewness of the vertical velocity fluctuations is generally positive but becomes negative in the lowest mesh cells when the dissipation rate exceeds the production rate due to buoyancy near the surface, as is the case for very rough surfaces. The LES results are used to determine the root-mean-square value of the surface friction velocity and the mean temperature difference between the surface and the mixed layer as a function of the roughness height. The results corroborate a simple model of the heat transfer in the surface layer.

Turbulent oscillatory boundary layers at high Reynolds numbers

Abstract: This study deals with turbulent oscillatory boundary-layer flows over both smooth and rough beds. The free-stream flow is a purely oscillating flow with sinusoidal velocity variation. Mean and turbulence properties were measured mainly in two directions, namely in the streamwise direction and in the direction perpendicular to the bed. Some measurements were made also in the transverse direction. The measurements were carried out up to Re = 6 × 106 over a mirror-shine smooth bed and over rough beds with various values of the parameter a/ks covering the range from approximately 400 to 3700, a being the amplitude of the oscillatory free-stream flow and ks the Nikuradse's equivalent sand roughness. For smooth-bed boundary-layer flows, the effect of Re is discussed in greater detail. It is demonstrated that the boundary-layer properties change markedly with Re. For rough-bed boundary-layer flows, the effect of the parameter a/ks is examined, at large values (O(103)) in combination with large Re.

Lagrangian statistics from direct numerical simulations of isotropic turbulence

Abstract: A comprehensive study is reported of the Lagrangian statistics of velocity, acceleration, dissipation and related quantities, in isotropic turbulence. High-resolution direct numerical simulations are performed on 643 and 1283 grids, resulting in Taylor-scale Reynolds numbers Rλ in the range 38-93. The low-wavenumber modes of the velocity field are forced so that the turbulence is statistically stationary. Using an accurate numerical scheme, of order 4000 fluid particles are tracked through the computed flow field, and hence time series of Lagrangian velocity and velocity gradients are obtained.The results reported include: velocity and acceleration autocorrelations and spectra; probability density functions (p.d.f.'s) and moments of Lagrangian velocity increments; and p.d.f.'s, correlation functions and spectra of dissipation and other velocity-gradient invariants. It is found that the acceleration variance (normalized by the Kolmogorov scales) increases as R½λ - a much stronger dependence than predicted by the refined Kolmogorov hypotheses. At small time lags, the Lagrangian velocity increments are distinctly non-Gaussian with, for example, flatness factors in excess of 10. The enstrophy (vorticity squared) is found to be more intermittent than dissipation, having a standard-deviation-to-mean ratio of about 1.5 (compared to 1.0 for dissipation). The acceleration vector rotates on a timescale about twice the Kolmogorov scale, while the timescales of acceleration magnitude, dissipation and enstrophy appear to scale with the Lagrangian velocity timescale.

A Global Climatology of Albedo, Roughness Length and Stomatal Resistance for Atmospheric General Circulation Models as Represented by the Simple Biosphere Model (SiB)

Abstract: Components of the Simple Biosphere Model (SiB) of Sellers et al. (1986) were used to generate global monthly fields of surface albedo (0.4-4.0 microns), roughness length and minimum surface (stomatal) resistance. SiB consists of three submodels which describe the roles of radiative transfer, turbulent transfer and surface resistance in determining the energy balance of the vegetated land surface. These three submodels were detached from SiB and used on the SiB parameter set (total and green leaf area index, leaf angle orientation, canopy dimensions, etc.) to calculate global monthly fields of albedo, roughness length and minimum stomatal resistance at 1 x 1 deg resolution. Time series of various parameters are also displayed for each vegetation type for specified grid points. The SiB results compare reasonably well with appropriate measurements obtained from the literature and have the additional merit of being mutually consistent; the three submodels use many common parameters, which ensures that, for each grid area, the calculated surface properties are closely interrelated as is the case in nature. The derived fields provide a check on the operation of the submodels and the correctness of the parameter set. They can also be used as prescribed fields for GCMs that do not have biophysically based land surface parameterizations.

Turbulent Boundary-Layer Separation

Abstract: This article summarizes our present understanding of the physical behavior of two-dimensional turbulent separated flows, which occur due to adverse pressure gradients around streamlined and bluff bodies. The physical behavior of turbulence is flow dependent, so detailed experimental infor­ mation is needed for understanding such flows and modeling their physics for calculation methods. An earlier review (Simpson 1 985) discussed in much detail prior experimental and computational work, and this was followed by an updated review of calculation methods only (Simpson 1 987). Here additional recent references are added to those cited in the two other works. By separation, we mean the entire process of departure or breakaway, or the breakdown of boundary-layer flow. An abrupt thickening of the rotational-flow region next to a wall and significant values of the normal­ to-wall velocity component must accompany breakaway, or otherwise this region would not have any significant interaction with the free-stream flow. This unwanted interaction causes a reduction in the performance of the flow device of interest (e.g. a loss of lift on an airfoil or a loss of pressure rise in a diffuser). It is too narrow a view to use vanishing surface shearing stress or flow reversal as the criterion for separation. Only in steady two-dimensional flow do these conditions usually accompany separation. In unsteady two­ dimensional flow the surface shear stress can change sign with flow reversal without the occurrence of breakaway_ Conversely, the breakdown of the boundary-layer concept can occur before any flow reversal is encountered. In three-dimensional flow the rotational layer can depart without the

Reynolds-number effects on the structure of a turbulent channel flow

Abstract: A high resolution, two component laser-Doppler anemometer has been used for turbulence measurements at a high data rate in a channel flow of water. Measurements of the velocity components in the stream direction and in a direction normal to the wall are reported over the Reynolds number range of 3000–40000. The combination of high spatial resolution and high data rates enabled accurate reconstruction of time dependent velocity traces. Long-time statistical averages of these signals clearly show that profiles of the dimensionless turbulence quantities such as turbulence intensities and Reynolds stress are strongly Reynolds-number dependent over a large part of the channel flow. For instance, in the Reynolds-number range of this investigation, it is shown that the fluctuating turbulence quantities do not scale with wall variables even as close as 15 viscous lengths from the wall. The velocity traces and associated power spectra exposed two phenomena which may explain the Reynolds number dependencies.

Characteristic-eddy decomposition of turbulence in a channel

Abstract: The proper orthogonal decomposition technique (Lumley's decomposition) is applied to the turbulent flow in a channel to extract coherent structures by decomposing the velocity field into characteristic eddies with random coefficients. In the homogeneous spatial directions, a generaliztion of the shot-noise expansion is used to determine the characteristic eddies. In this expansion, the Fourier coefficients of the characteristic eddy cannot be obtained from the second-order statistics. Three different techniques are used to determine the phases of these coefficients. They are based on: (1) the bispectrum, (2) a spatial compactness requirement, and (3) a functional continuity argument. Results from these three techniques are found to be similar in most respects. The implications of these techniques and the shot-noise expansion are discussed. The dominant eddy is found to contribute as much as 76 percent to the turbulent kinetic energy. In both 2D and 3D, the characteristic eddies consist of an ejection region straddled by streamwise vortices that leave the wall in the very short streamwise distance of about 100 wall units.

Advances in turbulence

Abstract: Contents: The State of Turbulence Research.- Contributions of Numerical Simulation Data Bases to the Physics, Modeling, and Measurement of Turbulence.- The Self-Preservation of Turbulent Flows and Its Relation to Initial Conditions and Coherent Structures.- Engineering Turbulence Models.- Chaos and the Onset of Turbulence.- Advances in Turbulence Measurement Techniques.- The Role of Smoke Visualization and Hot-Wire Anemometry in the Study of Transition.- Index.

Observation of organized structure in turbulent flow within and above a forest canopy

Abstract: Ramp patterns of temperature and humidity occur coherently at several levels within and above a deciduous forest as shown by data gathered with up to seven triaxial sonic anemometer/thermometers and three Lyman-alpha hygrometers at an experimental site in Ontario, Canada. The ramps appear most clearly in the middle and upper portion of the forest. Time/height cross-sections of scalar contours and velocity vectors, developed from both single events and ensemble averages of several events, portray details of the flow structures associated with the scalar ramps. Near the top of the forest they are composed of a weak ejecting motion transporting warm and/or moist air out of the forest followed by strong sweeps of cool and/or dry air penetrating into the canopy. The sweep is separated from the ejecting air by a sharp scalar microfront. At approximately twice the height of the forest, ejections and sweeps are of about equal strength.

Bar-generating cross-shore flow mechanisms on a beach

Abstract: Random waves normally incident on a dissipative beach induce a variety of cross-shore flows, such as asymmetric oscillatory flow, wave grouping-induced long-wave flow, breaking-induced turbulent flow, and momentum decay-induced undertow. These flows are identified, analyzed and hindcasted in a set of laboratory experiments with the aim of revealing the role of each of the flow mechanisms in the two-dimensional case of bar generation on a beach.

Transport of Passive Scalars in a Turbulent Channel Flow

Abstract: A direct numerical simulation of a turbulent channel flow with three passive scalars at different molecular Prandtl numbers is performed. Computed statistics including the turbulent Prandtl numbers are compared with existing experimental data. The computed fields are also examined to investigate the spatial structure of the scalar fields. The scalar fields are highly correlated with the streamwise velocity; the correlation coefficient between the temperature and the streamwise velocity is as high as 0.95 in the wall region. The joint probability distributions between the temperature and velocity fluctuations are also examined; they suggest that it might be possible to model the scalar fluxes in the wall region in a manner similar to the Reynolds stresses.

Drag and Reconfiguration of Broad Leaves in High Winds

Abstract: Drag was measured and changes of configuration noted as a variety of leaves, leaflets, and clusters were subjected to turbulent winds of 10 and 20 m s-1. Leaves with acute bases and short petioles had the highest surface-specific drag, fluttered erratically and, most commonly, tore. Leaves with lobed bases and long petioles had lower drag, fluttered little and reconfigured into increasingly acute cones. Pinnately compound leaves had the lowest drag and formed cylinders with leaflets layered alternately. For all but individual white oak leaves, drag coefficients (based on original surface area) decreased with increasing wind speed. Single leaves of white poplar were unstable at all speeds but resisted damage even at 30 m s"1; clusters formed stable cones. These results are contrasted with the behaviour of flags in wind and are related to wind-throw in trees.

Aerosol particle deposition in numerically simulated channel flow

Abstract: The trajectories of rigid spherical particles in a turbulent channel flow are computed using a pseudospectral computer program to simulate the three‐dimensional, time‐dependent flow field. It is assumed that the channel is vertical so that gravity cannot directly cause the deposition of particles on the walls. The particles are assumed to be sufficiently small and widely separated so that their influence on the fluid velocity field can be ignored. It is found that when the particles are assigned random initial locations with initial velocities that are equal to the local fluid velocity, the particles tend to accumulate in the viscous sublayer. At the edge of the viscous sublayer, the particles that deposit on the wall typically possess normal components of velocity that are comparable in magnitude to the intensity of the normal component of the velocity in the core of the channel (i.e., of the order of magnitude of the friction velocity). A shear‐induced lift force having the form derived by Saffman for laminar flow is found to have virtually no effect on particle trajectories, except within the viscous sublayer where it plays a significant role both in the inertial deposition of particles and in the accumulation of trapped particles. The Reynolds number of the particles that deposit does not remain small compared with unity.

Evaluation of turbulent transport and dissipation closures in second-order modeling

Abstract: We show that the turbulence statistics from our (96)3 large-eddy-simulation (LES) studies of a convective boundary layer are in excellent agreement with those from the Deardorff–Willis laboratory convection tank. Using these LES data, we evaluate contemporary parameterizations for turbulent transport and dissipation in second-order closure models of the convective boundary layer. The gradient-diffusion parameterization for turbulent transport fares poorly, due in large part to the direct influence of buoyancy. This leads to poor predictions of the vertical profiles of some turbulence statistics. We also find that the characteristic length scales for the mechanical and thermal dissipation rates typically used in second-order closure models are a factor of 2–3 too small; this leads to underpredictions of turbulence kinetic energy levels. Finally, we find that the flux and variance budgets for conservative scalars are substantially different in top-down and bottom-up diffusion. In order to reproduce...

New approximate boundary conditions for large eddy simulations of wall-bounded flows

Abstract: Two new approximate boundary conditions have been applied to the large eddy simulation of channel flow with and without transpiration. These new boundary conditions give more accurate results than those previously in use, and allow significant reduction of the required CPU time over simulations in which no‐slip conditions are applied. Mean velocity profiles and turbulence intensities compare well both with experimental data and with the results of resolved simulations. The influence of the approximate boundary conditions remains confined near the point of application and does not affect the turbulence statistics in the core of the flow.

Spectral approach to non-isotropic turbulence subjected to rotation

Abstract: The non-isotropic effects of solid-body rotation on homogeneous turbulence are investigated in this paper. A spectral formalism using eigenmodes introduces the spectral Coriolis effects more easily and leads to simpler expressions for the integral quadratic terms which come mostly from classical two-point closures. The analysis is then applied to a specific eddy damped quasi-normal Markovian model, which includes the inertial waves regime in the evaluation of triple correlations. This procedure allows for a departure from isotropy by external rotation effects. When started with rigorously isotropic initial data, the various trends observed on the Reynolds stresses and the integral lengthscales remain in accordance with the results from direct simulations; moreover they reflect a very specific spectral angular distribution. Such an angular dependence allows a drain of spectral energy from the parallel to the normal wave vectors (with respect to the rotation axis), and thus appears consistent with a trend toward two-dimensionality.

Near-wall structure of a turbulent boundary layer with riblets

Abstract: A detailed wind tunnel study has been carried out on the near-wall turbulence structure over smooth and riblet wall surfaces under zero pressure gradient. Time-average quantities as ‘well as conditionally sampled profiles were obtained using hotwire/film anemometry, along with a simultaneous flow visualization using the smoke-wire technique and a sheet of laser light. The experimental results indicated a significant change of the structure in the turbulent boundary layer near the riblet surface. The change was confined within a small volume of the flow close to the wall surface. A conceptual model for the sequence of the bursts was then proposed based on an extensive study of the flow visualization, and was supported by the results of conditionally sampled velocity fields. A possible mechanism of turbulent drag reduction by riblets is discussed.

On the structure of pressure fluctuations in simulated turbulent channel flow

Abstract: Pressure fluctuations in a turbulent channel flow are investigated by analyzing a database obtained from a direct numerical simulation. Detailed statistics associated with the pressure fluctuations are presented. Characteristics associated with the rapid (linear) and slow (nonlinear) pressure are discussed. It is found that the slow pressure fluctuations are larger than the rapid pressure fluctuations throughout the channel except very near the wall, where they are about the same magnitude. This is contrary to the common belief that the nonlinear source terms are negligible compared to the linear source terms. Probability density distributions, power spectra, and two-point correlations are examined to reveal the characteristics of the pressure fluctuations. The global dependence of the pressure fluctuations and pressure-strain correlations are also examined by evaluating the integral associated with Green's function representations of them. In the wall region where the pressure-strain terms are large, most contributions to the pressure-strain terms are from the wall region (i.e., local), whereas away from the wall where the pressure-strain terms are small, contributions are global. Structures of instantaneous pressure and pressure gradients at the wall and the corresponding vorticity field are examined.

Second-moment closure and its use in modelling turbulent industrial flows

Abstract: Second-moment turbulence models focus directly on the transport equations for the Reynolds stresses rather than supposing the stress and strain fields to be directly linked via an eddy viscosity. This elaboration enables the effects of complex strains and force fields on the turbulence structure to be better captured. The paper summarizes the principal modelling strategies adopted for the unknown processes in these equations and presents the forms that have been found most useful in engineering calculations. Methods adopted for overcoming significant problems of numerical instability and lack of convergence compared with eddy-viscosity-based schemes are also presented. Applications involving momentum and heat transfer in complex flows are drawn from the advanced technology sectors of the power generation and aircraft industries.

Structure of the turbulent flow field under breaking waves in the surf zone

Abstract: The structure of turbulence and its role in the breaking wave dynamics within the surf zone have been investigated through laboratory experiments using several flow visualization techniques and a fibre-optic LDV system. The results indicate that there exists a characteristic structure of large-scale eddies referred to here as ‘horizontal eddies’ and ‘obliquely descending eddies’, which has a significant role in the generation of Reynolds stress and thus affects the deformation of the mean flow field. The experiments also reveal that these eddies caused by the wave breaking bring a large amount of vorticity (with non-zero average) into otherwise almost irrotational velocity fields, resulting in the generation of vorticity-related mean flow fields as well as turbulence (vorticity-containing velocity fluctuation). This means that the breaking waves in the surf zone can be regarded as pseudowaves which consist of irrotational velocity components as ‘wave motion’ and appreciable amounts of rotational mean velocity components as ‘eddying motion’ (with non-zero mean vorticity) together with turbulence. It is found that the generation of the mean rotational velocity component due to wave breaking causes considerable increase in mass and momentum transport, as compared with ordinary non-breaking waves, and thus a decrease in wave height.

Direct numerical simulation of stratified homogeneous turbulent shear flows

Abstract: The exact time-dependent three-dimensional Navier-Stokes and temperature equations are integrated numerically to simulate stably stratified homogeneous turbulent shear flows at moderate Reynolds numbers whose horizontal mean velocity and mean temperature have uniform vertical gradients. The method uses shear-periodic boundary conditions and a combination of finite-difference and pseudospectral approximations. The gradient Richardson number Ri is varied between 0 and 1. The simulations start from isotropic Gaussian fields for velocity and temperature both having the same variances. The simulations represent approximately the conditions of the experiment by Komori et al. (1983) who studied stably stratified flows in a water channel (molecular Prandtl number Pr = 5). In these flows internal gravity waves build up, superposed by hot cells leading to a persistent counter-gradient heat-flux (CGHF) in the vertical direction, i.e. heat is transported from lower-temperature to higher-temperature regions. Further, simulations with Pr = 0.7 for air have been carried out in order to investigate the influence of the molecular Prandtl number. In these cases, no persistent CGHF occurred. This confirms our general conclusion that the counter-gradient heat flux develops for strongly stable flows (Ri [approximate] 0.5–1.0) at sufficiently large Prandtl numbers (Pr = 5). The flux is carried by hot ascending, as well as cold descending turbulent cells which form at places where the highest positive and negative temperature fluctuations initially existed. Buoyancy forces suppress vertical motions so that the cells degenerate to two-dimensional fossil turbulence. The counter-gradient heat flux acts to enforce a quasi-static equilibrium between potential and kinetic energy. Previously derived turbulence closure models for the pressure-strain and pressure-temperature gradients in the equations for the Reynolds stress and turbulent heat flux are tested for moderate-Reynolds-number flows with strongly stable stratification (Ri = 1). These models overestimate the turbulent interactions and underestimate the buoyancy contributions. The dissipative timescale ratio for stably stratified turbulence is a strong function of the Richardson number and is inversely proportional to the molecular Prandtl number of the fluid.

Reynolds number and end‐wall effects on a lid‐driven cavity flow

Abstract: A series of experiments has been conducted in a lid‐driven cavity of square cross section (depth=width=150 mm) for Reynolds numbers (Re, based on lid speed and cavity width) between 3200 and 10 000, and spanwise aspect ratios (SAR) between 0.25:1 and 1:1. Flow visualization using polystyrene beads and two‐dimensional laser‐Doppler anemometer (LDA) measurements have shed new light on the momentum transfer processes within the cavity. This paper focuses on the variation, with Re and SAR, of the mean and the rms velocities profiles, as well as the ∼(U’V’) profile, along the horizontal and vertical centerlines in the symmetry plane. In addition, the contribution of the large‐scale ‘‘organized structures,’’ and the high‐frequency ‘‘turbulent’’ velocity fluctuations to the total rms is examined. At low Re, the organized structures account for most of the energy contained in the flow irrespective of SAR. As the Re increases, however, so does the energy content of the higher frequency fluctuations. This trend is ...