scispace - formally typeset
Search or ask a question

Showing papers on "Turbulence published in 1996"


Journal ArticleDOI
TL;DR: A scheme for the characterization of the quality of direct turbulence measurements using the Navier-Stokes equations, which are the basis of direct surface-based measurements of turbulent fluxes via the eddy correlation method.

1,456 citations


Journal ArticleDOI
TL;DR: In this paper, the Smagorinsky eddy-viscosity model is combined with a spatially averaged dynamic model for complex-geometry inhomogeneous flows, and a new dynamic model formulation is introduced that combines advantages of the statistical and local approaches.
Abstract: The dynamic model for large-eddy simulation of turbulence samples information from the resolved velocity field in order to optimize subgrid-scale model coefficients. When the method is used in conjunction with the Smagorinsky eddy-viscosity model, and the sampling process is formulated in a spatially local fashion, the resulting coefficient field is highly variable and contains a significant fraction of negative values. Negative eddy viscosity leads to computational instability and as a result the model is always augmented with a stabilization mechanism. In most applications the model is stabilized by averaging the relevant equations over directions of statistical homogeneity. While this approach is effective, and is consistent with the statistical basis underlying the eddy-viscosity model, it is not applicable to complex-geometry inhomogeneous flows. Existing local formulations, intended for inhomogeneous flows, are most commonly stabilized by artificially constraining the coefficient to be positive. In this paper we introduce a new dynamic model formulation, that combines advantages of the statistical and local approaches. We propose to accumulate the required averages over flow pathlines rather than over directions of statistical homogeneity. This procedure allows the application of the dynamic model with averaging to in-homogeneous flows in complex geometries. We analyse direct numerical simulation data to document the effects of such averaging on the Smagorinsky coefficient. The characteristic Lagrangian time scale over which the averaging is performed is chosen based on measurements of the relevant Lagrangian autocorrelation functions, and on the requirement that the model be purely dissipative, guaranteeing numerical stability when coupled with the Smagorinsky model. The formulation is tested in forced and decaying isotropic turbulence and in fully developed and transitional channel flow. In homogeneous flows, the results are similar to those of the volume-averaged dynamic model, while in channel flow, the predictions are slightly superior to those of the spatially (planar) averaged dynamic model. The relationship between the model and vortical structures in isotropic turbulence, as well as ejection events in channel flow, is investigated. Computational overhead is kept small (about 10% above the CPU requirements of the spatially averaged dynamic model) by using an approximate scheme to advance the Lagrangian tracking through first-order Euler time integration and linear interpolation in space.

1,149 citations


Journal ArticleDOI
TL;DR: In this article, a large-eddy simulation (LES) formalism, along with various subgrid-scale models developed since Smagorinsky's model, is presented, with an emphasis on the generation of coherent vortices.
Abstract: The paper presents large-eddy simulation (LES) formalism, along with the various subgrid-scale models developed since Smagorinsky’s model. We show how Kraichnan’s spectral eddy viscosity may be implemented in physical space, yielding the structure-function model. Recent developments of this model that allow the eddy viscosity to be inhibited in transitional regions are discussed. We present a dynamic procedure, where a double filtering allows one to dynamically determine the subgrid-scale model constants. The importance of backscatter effects is discussed. Alternatives to the eddy-viscosity assumption, such as scale- similarity models, are considered. Pseudo-direct simulations in which numerical diffusion replaces subgrid transfers are mentioned. Various applications of LES to incompressible and compressible turbulent flows are given, with an emphasis on the generation of coherent vortices

1,101 citations


Book ChapterDOI
TL;DR: In this paper, the authors argue that the active turbulence and coherent motions near the top of a vegetation canopy are patterned on a plane mixing layer, because of instabilities associated with the characteristic strong inflection in the mean velocity profile.
Abstract: This paper argues that the active turbulence and coherent motions near the top of a vegetation canopy are patterned on a plane mixing layer, because of instabilities associated with the characteristic strong inflection in the mean velocity profile. Mixing-layer turbulence, formed around the inflectional mean velocity profile which develops between two coflowing streams of different velocities, differs in several ways from turbulence in a surface layer. Through these differences, the mixing-layer analogy provides an explanation for many of the observed distinctive features of canopy turbulence. These include: (a) ratios between components of the Reynolds stress tensor; (b) the ratio K H /K M of the eddy diffusivities for heat and momentum; (c) the relative roles of ejections and sweeps; (d) the behaviour of the turbulent energy balance, particularly the major role of turbulent transport; and (e) the behaviour of the turbulent length scales of the active coherent motions (the dominant eddies responsible for vertical transfer near the top of the canopy). It is predicted that these length scales are controlled by the shear length scale L s = U(h)/U′(h) (where h is canopy height, U(z) is mean velocity as a function of height z, and U′ = dU/dz). In particular, the streamwise spacing of the dominant canopy eddies Λ x = mL s , with m = 8.1. These predictions are tested against many sets of field and wind-tunnel data. We propose a picture of canopy turbulence in which eddies associated with inflectional instabilities are modulated by larger-scale, inactive turbulence, which is quasi-horizontal on the scale of the canopy.

1,094 citations


Book
01 Jan 1996
TL;DR: In this paper, the authors present solutions for nutrient transfer to osmotrophs in the full range of flow regimes for which solutions have been published, and extend some of those solutions to new parameter domains and flow environments.
Abstract: We present solutions for nutrient transfer to osmotrophs in the full range of flow regimes for which solutions have been published, and we extend some of those solutions to new parameter domains and flow environments. These regimes include stagnant water; steady, unitorm flow arising from swimming or sinking; steady shear flows; and fluctuating shear from dissipation of turbulence, as well as the combined effects of turbulence-induced shear and swimming or sinking. Solutions for nutrient fluxes cannot be carried over from one flow regime to another. In all cases, however, mass transfer increases with cell size and with flow velocity. Cell shape becomes particularly important at high flow velocities. For steady, uniform flow arising from sinking or swimming, we find asymptotic analytic and numerical solutions from the engineering literature superior to those in more common use within oceanography. These engineering solutions suggest flow effects an order of magnitude smaller than commonly supposed. A cell radius near 20 Ilm is needed before swimming or sinking can be expected to increase the flux of nutrients, such as nitrate or phosphate, substantially (by ~50%) over the stagnant-water case. We find sound asymptotic solutions for the case oflinear shear and supplement them with numerical solutions of our own to cover the range of cell sizes and shear rates of interest for phytoplankton. We extend them further to cover viscous shears from dissipating turbulence for cells smaller than the Kolmogorov scale (order of 1-6mm in the ocean). Our analysis suggests turbulence effects an order of magnitude greater than previously postulated, with a cell size of 60~lm needed to experience substantial gain. Cell rotation, whether induced by the propulsion mechanism in swimming or passively by shear across the cell perimeter, will reduce the rate of nutrient transfer relative to a non-rotating cell unless the axis of rotation parallels the direction of flow. Although in calm water dinoflagellates by swimming are able to increase nutrient uptake, in strong turbulence they may not be able to maintain a rotational axis parallel to the direction of swimming or the direction of shear, resulting in a relative reduction in flux. Conversely, large chains of diatoms and filamentous cyanobacteria that span the radius of the smallest vortices are best able to take advantage of turbulence. Despite these deductions from a diversity of analytic and numerical solutions, unequivocal data to test the contribution of advection to nutrient acquisition by phytoplankton are scarce owing, in large part, to the inability to visualize, record and thus mimic fluid motions in the vicinities of cells in natural flows.

697 citations


Journal ArticleDOI
TL;DR: In this article, two new methods for distinguishing two-dimensional (2D) turbulence from slab turbulence are applied to Helios magnetometer data, and they indicate that solar wind magnetic turbulence possesses a dominant (∼85 % by energy) 2D component.
Abstract: Two new methods for distinguishing two-dimensional (2D) turbulence from slab turbulence are applied to Helios magnetometer data. Two-component models with varying slab and 2D ingredients are considered. Both methods indicate that solar wind magnetic turbulence possesses a dominant (∼85 % by energy) 2D component. The presence of such a large 2D component provides a natural solution to the long-standing problem of “too small” cosmic ray mean free paths derived from quasilinear scattering theory when using the slab model.

618 citations


Journal ArticleDOI
TL;DR: In this article, a review of numerical models for turbulent fluid-particle flows is presented, which is structured according to the turbulence models used for the continuous phase: turbulence energy-dissipation models, large eddy simulations, direct numerical simulations, and discrete vortex models.
Abstract: Numerical models for turbulent fluid-particle flows are reviewed. The two approaches typically used for modelling the dispersed (particle) phase are the trajectory and two-fluid formulations, while volume- averaged models are most common for the continuous (fluid) phase. The review is structured according to the turbulence models used for the continuous phase: turbulence energy-dissipation models, large- eddy simulations, direct numerical simulations, and discrete vortex models. The applications of these models to simulate particle dispersion due to fluid turbulence and the adjustments to the models to account for the modulation of the carrier phase turbulence by the particles are addressed.

592 citations


Journal ArticleDOI
TL;DR: In this article, a formalism for analyzing errors in nonlinear problems is developed in the context of finite difference approximations for the Navier?Stokes equations when the flow is fully turbulent.

590 citations


Journal ArticleDOI
TL;DR: In this article, a cubic relation between the strain and vorticity tensor and the stress tensor was proposed, which does much better than a conventional eddy-viscosity scheme in capturing effects of streamline curvature over a range of flows.

578 citations


Journal ArticleDOI
TL;DR: In this paper, a numerical model has been built to simulate the small-scale atmospheric flows within the urban canopy, based on the lower atmosphere classical assumptions and the standard κ-ϵ” two-equation turbulence model.

557 citations


Journal ArticleDOI
TL;DR: In this article, a theory was developed to correct mean-velocity profiles for the effects of wandering and to provide complete quantitative estimates of its amplitude and contributions to Reynolds stress fields.
Abstract: Experiments have been performed on the tip vortex trailing from a rectangular NACA 0012 half-wing. Preliminary studies showed the vortex to be insensitive to the introduction of a probe and subject only to small wandering motions. Meaningful velocity measurements could therefore be made using hot-wire probes.Detailed analysis of the effects of wandering was performed to properly reveal the flow structure in the core region and to give confidence in measurements made outside the core. A theory has been developed to correct mean-velocity profiles for the effects of wandering and to provide complete quantitative estimates of its amplitude and contributions to Reynolds stress fields. Spectral decomposition was found to be the most effective method of separating these contributions from velocity fluctuations due to turbulence.Outside the core the flow structure is dominated by the remainder of the wing wake which winds into an ever-increasing spiral. There is no large region of axisymmetric turbulence surrounding the core and little sign of turbulence generated by the rotational motion of the vortex. Turbulence stress levels vary along the wake spiral in response to the varying rates of strain imposed by the vortex. Despite this complexity, the shape of the wake spiral and its turbulent structure reach an approximately self-similar form.On moving from the spiral wake to the core the overall level of velocity fluctuations greatly increases, but none of this increase is directly produced by turbulence. Velocity spectra measured at the vortex centre scale in a manner that implies that the core is laminar and that velocity fluctuations here are a consequence of inactive motion produced as the core is buffeted by turbulence in the surrounding spiral wake. Mean-velocity profiles through the core show evidence of a two-layered structure that dies away with distance downstream.

Journal ArticleDOI
TL;DR: In this article, a general classification of two-phase flows and a number of possible ways to formulate two-fluid models are discussed, and a general procedure to develop such a model is presented.

Journal ArticleDOI
TL;DR: In this article, the statistical uncertainty associated with the sampling of random processes such as those which occur in turbulence research are given, and formulas based on normal distribution assumptions and on any general distribution shape are given for means, variances, Reynolds stresses, correlation coefficients, homogeneous and mixed turbulent triple products and fourth order turbulence moments.
Abstract: Methods for calculating the statistical uncertainty associated with the sampling of random processes such as those which occur in turbulence research are given. In particular, formulas based on normal distribution assumptions and on any general distribution shape are given for means, variances, Reynolds stresses, correlation coefficients, homogeneous and mixed turbulent triple products and fourth order turbulence moments. In addition, two resampling algorithms, the “bootstrap” and “jackknife”, are presented and compared using actual turbulence data. The availability of these methods will allow turbulence data to be presented with statistical uncertainty error bars on all turbulence quantities.

Journal ArticleDOI
TL;DR: In this article, the authors studied the Taylor-microscale Reynolds number in a wind tunnel with triangular winged grid bars and determined a linear dependence between n1 and C1*: C1* = 4.5 − 2.4n1.
Abstract: Using an active grid devised by Makita (1991), shearless decaying turbulence is studied for the Taylor-microscale Reynolds number, Rλ, varying from 50 to 473 in a small (40 × 40 cm2 cross-section) wind tunnel. The turbulence generator consists of grid bars with triangular wings that rotate and flap in a random way. The value of Rλ is determined by the mean speed of the air (varied from 3 to 14 m s–1) as it passes the rotating grid, and to a lesser extent by the randomness and rotation rate of the grid bars. Our main findings are as follows. A weak, not particularly well-defined scaling range (i.e. a power-law dependence of both the longitudinal (u) and transverse (v) spectra, F11(k1) and F22(k1) respectively, on wavenumber k1) first appears at Rλ ∼ 50, with a slope, n1, (for the u spectrum) of approximately 1.3. As Rλ was increased, n1 increased rapidly until Rλ ∼ 200 where n ∼ 1.5. From there on the increase in n1 was slow, and even by Rλ = 473 it was still significantly below the Kolmogorov value of 1.67. Over the entire range, 50 [les ] Rλ [les ] 473, the data were well described by the empirical fit: . Using a modified form of the Kolmogorov similarity law: F11(k1) = C1*e2/3k1–5/3(k 1η)5/3–n1 where e is the turbulence energy dissipation rate and η is the Kolmogorov microscale, we determined a linear dependence between n1 and C1*: C1* = 4.5 – 2.4n1. Thus for n1 = 5/3 (which extrapolation of our results suggests will occur in this flow for Rλ ∼ 104), C1* = 0.5, the accepted high-Reynolds-number value of the Kolmogorov constant. Analysis of the p.d.f. of velocity differences Δu(r) and Δv(r) where r is an inertial subrange interval, conditional dissipation, and other statistics showed that there was a qualitative difference between the turbulence for Rλ 200 (strong turbulence). For the latter, the p.d.f.s of Δu(r) and Δv(r) had super Gaussian tails and the dissipation (both of the u and v components) conditioned on Δu(r) and Δv(r) was a strong function of the velocity difference. For Rλ 200 are consistent with the predictions of the Kolmogorov refined similarity hypothesis (and make a distinction between the dynamical and kinematical contributions to the conditional statistics). They have much in common with similar statistics done in shear flows at much higher Rλ, with which they are compared.

Journal ArticleDOI
TL;DR: In this paper, a large-eddy simulation was used to study mixing of turbulent, coannular jets discharging into a sudden expansion, which resembles that of a coaxial jet-combustor, and the goal of the calculation was to gain some insight into the phenomena leading to lean blowout in such combustion devices.
Abstract: Large-eddy simulation (LES) was used to study mixing of turbulent, coannular jets discharging into a sudden expansion. This geometry resembles that of a coaxial jet-combustor, and the goal of the calculation was to gain some insight into the phenomena leading to lean blow-out (LBO) in such combustion devices. This is a first step in a series of calculations, where the focus is on the fluid dynamical aspects of the mixing process in the combustion chamber. The effects of swirl, chemical reactions and heat release were not taken into account. Mixing of fuel and oxidizer was studied by tracking a passive scalar introduced in the central jet. The dynamic subgrid-scale (DM) model was used to model both the subgrid-scale stresses and the subgrid-scale scalar flux. The Reynolds number was 38000, based on the bulk velocity and diameter of the combustion chamber. Mean velocities and Reynolds stresses are in good agreement with experimental data. Animated results clearly show that intermittent pockets of fuel-rich fluid (from the central jet) are able to cross the annular jet, virtually undiluted, into the recirculation zone. Most of the fuel-rich fluid is, however, entrained into the recirculation zone near the instantaneous reattachment point. Fuel trapped in the recirculation zone is, for the most part, entrained back into the step shear layer close to the base of the burner.

Journal Article
TL;DR: In this article, high-frequency (5 Hz) in situ measurements of flow speed were collected in Spartina alternijlora, Juncus roemerianus, and Di-tichlis spicata canopies using hot-film anemometry sensor arrays.
Abstract: The transport of particulate and dissolved matter on the surface of coastal marshes is controlled by the hydrodynamic characteristics of over-marsh flows. High-frequency (5 Hz) in situ measurements of flow speed were collected in Spartina alternijlora, Juncus roemerianus, and Di:-tichlis spicata canopies using hot-film anemometry sensor arrays. These data indicate that mean flow speed, turbulence intensity, and the shape of the vertical speed profile are influenced by variations in plant morphology and stem density. Mean flow speed and turbulence intensity are inversely related to stem density and to distance from the creek edge. Flow energies decrease by about one order of magnitude when flows encounter the vegetated marsh surface and continue to decrease as vegetation density increa,ses. Turbulent flow energy also decays exponentially with increasing distance from the creek edge. Reductions in flow speed coupled with energy decay provide a hydrologic mechanism for sediment deposition patterns commonly observed in marsh systems. Suspended matter transport is also affected by plant-flow interactions. Vertical flow structure is strongly influenced by canopy morphology (plant type and plant shape). Plant-flow interactions result in vertical speed profiles whose shapes deviate from the logarithmic profile typical in free-stream conditions and in the development of transitional flow regimes (i.e. neither laminar nor fully turbulent).

Journal ArticleDOI
TL;DR: In this paper, hot-wire measurements were conducted in the very near wake (x/d⩽10) of a circular cylinder at a Reynolds number based on cylinder diameter.
Abstract: Hot-wire measurements were conducted in the very near wake (x/d⩽10) of a circular cylinder at a Reynolds number based on cylinder diameter, Re d of 3900. Measurements of the streamwise velocity component with the use of single sensor hot-wire probes were found to be inaccurate for such flowfields where high flow angles are present. An X-array probe provided detailed streamwise and lateral velocity component statistics. Frequency spectra of these two velocity components are also presented. Measurements with a 4-sensor hot-wire probe confirmed that the very near wake region is dominantly two-dimensional, thus validating the accuracy of the present X-array data.

Journal ArticleDOI
TL;DR: In this article, a filter-structure-function (FSF) model is proposed for the simulation of a quasi-incompressible boundary layer developing spatially over an adiabatic flat plate, with a low level of upstream forcing.
Abstract: It is well known that subgrid models such as Smagorinsky's cannot be used for the spatially growing simulation of the transition to turbulence of flat-plate boundary layers, unless large-amplitude perturbations are introduced at the upstream boundary: they are over-dissipative, and the flow simulated remains laminar. This is also the case for the structure-function model (SF) of Metais & Lesieur (1992). In the present paper we present a sequel to this model, the filtered-structure-function (FSF) model. It consists of removing the large-scale fluctuations of the field before computing its second-order structure function. Analytical arguments confirm the superiority of the FSF model over the SF model for large-eddy simulations of weakly unstable transitional flows. The FSF model is therefore used for the simulation of a quasi-incompressible (M∞ = 0.5) boundary layer developing spatially over an adiabatic flat plate, with a low level of upstream forcing. With the minimal resolution 650 × 32 × 20 grid points covering a range of streamwise Reynolds numbers Rex1 e [3.4 × 105, 1.1 × 106], transition is obtained for 80 hours of time-processing on a CRAY 2 (whereas DNS of the whole transition takes about ten times longer). Statistics of the LES are found to be in acceptable agreement with experiments and empirical laws, in the laminar, transitional and turbulent parts of the domain. The dynamics of low-pressure and high-vorticity distributions is examined during transition, with particular emphasis on the neighbourhood of the critical layer (defined here as the height of the fluid travelling at a speed equal to the phase speed of the incoming Tollmien–Schlichting waves). Evidence is given that a subharmonic-type secondary instability grows, followed by a purely spanwise (i.e. time-independent) mode which yields peak-and-valley splitting and transition to turbulence. In the turbulent region, flow visualizations and local instantaneous profiles are provided. They confirm the presence of low- and high-speed streaks at the wall, weak hairpins stretched by the flow and bursting events. It is found that most of the vorticity is produced in the spanwise direction, at the wall, below the high-speed streaks. Isosurfaces of eddy viscosity confirm that the FSF model does not perturb transition much, and acts mostly in the vicinity of the hairpins.

Journal ArticleDOI
TL;DR: In this article, it was shown that wave packets belonging to the inertial range exhibit properties intermediate between those of weak and strong turbulent cascades, and that successive collisions are correlated since wave packets are distorted by diverging field lines.
Abstract: Kraichnan (1965) proposed that MHD turbulence occurs as a result of collisions between oppositely directed Alfv\'en wave packets. Recent work has generated some controversy over the nature of non linear couplings between colliding Alfv\'en waves. We find that the resolution to much of the confusion lies in the existence of a new type of turbulence, intermediate turbulence, in which the cascade of energy in the inertial range exhibits properties intermediate between those of weak and strong turbulent cascades. Some properties of intermediate MHD turbulence are: (i) in common with weak turbulent cascades, wave packets belonging to the inertial range are long lived; (ii) however, components of the strain tensor are so large that, similar to the situation in strong turbulence, perturbation theory is not applicable; (iii) the breakdown of perturbation theory results from the divergence of neighboring field lines due to wave packets whose perturbations in velocity and magnetic fields are localized, but whose perturbations in displacement are not; (iv) 3--wave interactions dominate individual collisions between wave packets, but interactions of all orders $n\geq 3$ make comparable contributions to the intermediate turbulent energy cascade; (v) successive collisions are correlated since wave packets are distorted as they follow diverging field lines; (vi) in common with the weak MHD cascade, there is no parallel cascade of energy, and the cascade to small perpendicular scales strengthens as it reaches higher wave numbers; (vii) For an appropriate weak excitation, there is a natural progression from a weak, through an intermediate, to a strong cascade.

Journal ArticleDOI
TL;DR: The available data describing the incompressible zero-pressure gradient boundary layer are reviewed in association with a range of studies which are either new or, to date, not very generally available as discussed by the authors.

Book
01 Jan 1996
TL;DR: In this article, the authors investigated the air bubble entrainment in free-surface turbulent shear flows and developed an analysis of the air entraining processes in free surface flows.
Abstract: In high velocity water flows, large quantities of air bubbles are entrained at the free-surfaces. Practical applications are found in Chemical, Civil, Environmental, Mechanical, Mining and Nuclear Engineering. Air-water flows are observed in small-scale as well as large-scale flow situations. E.g., thin circular jets used as mixing devices in chemical plants (Qw ~ 0.001 L/s, diameter ~ 1 mm), and spillway flows (Qw > 10,000 m3/s, flow thickness over 10 m). In each case, however, the interactions between the entrained air bubbles and the turbulence field are significant. This monograph investigates the "air bubble entrainment in free-surface turbulent shear flows". It develops an analysis of the air entrainment processes in free-surface flows. The air-water flows are investigated as homogeneous mixtures with variable density. The variations of fluid density result from the non-uniform air bubble distributions and the turbulent diffusion process. Several types of air-water free-surface flows are studied : plunging jet flows (Part II), open channel flows (Part III), and turbulent water jets discharging into air (Part IV). Each configuration can be characterised as a high-velocity free-surface flow with turbulent shear layer and large air bubble content. Experimental observations confirm the conceptual idea that the air-water mixture behaves as a homogeneous compressible fluid. The monograph presents numerous and recent experimental investigations with mean velocities up to 57 m/s and mean air contents up to 70%. The analysis of experimental studies provides new information on the air-water flow field : air bubble distributions, air-water velocity profiles, air bubble sizes and bubble-turbulence interactions. The results show a strong similarity between all the flow patterns. In each case the distributions of air concentration (i.e. void fraction) can be approximated by a simple advective diffusion theory. New analysis is developed for each flow configuration and compared successfully with model and prototype data. The velocity distributions in air-water flows have the same shape as for monophase flows. However the presence of air bubbles modifies some turbulence characteristics while the turbulence controls the mechanism of bubble breakup. The book presents new useful information for design engineers and research-and-development scientists who need a better understanding of the fluid mechanics of air-water flows. Both qualitative and quantitative information are provided. In some cases the limits of our knowledge are pointed out. The book consists of five parts. Part I introduces the topic and its relevance, develops a dimensional analysis and discusses the air-water gas transfer process. In each subsequent part, the distributions of air content and air-water velocity are described. The results are grouped as : plunging jet flows (Part II), open channel flows (Part III) and high-velocity water jets discharging into the atmosphere (Part IV). In Part V, an analogy between the various types of air-water flows is developed. In the appendices, tables of physical and chemical properties of fluids are provided in appendix A. The report presents results expressed in SI Units. A table of unit conversions is given in appendix B. Estimates of bubble rise velocity are discussed in appendix C. Appendix D develops sound celerity calculations in two-phase flows. Appendices E, G, H and I present complete calculations of the air bubble diffusion process. Boundary layer characteristics and jet trajectory calculations are detailed in appendices F and J respectively. Appendix K defines bubble size distribution characteristic parameters. Observations by LEONARDO DA VINCI are recounted in appendix L. 'Errare Humanum Est'. Appendix M presents a correction form. Readers who find an error or mistake are welcome to record the error on the page and to send a copy to the author. At the beginning of the book, the reader will find the table of contents, a list of symbols, a glossary and an album of colourful photographs of 'white waters'.

Journal ArticleDOI
TL;DR: In this paper, a new method to model the effect of the solid boundaries on the rest of the flowfield in large-eddy simulations is proposed, where the filtered Navier-Stokes equations are solved up to the first computational point from there to the wall, a simplified set of equations is solved, and an estimate of the instantaneous wall shear stress required to impose boundary conditions is obtained.
Abstract: A new method to model the effect of the solid boundaries on the rest of the flowfield in large-eddy simulations is proposed The filtered Navier-Stokes equations are solved up to the first computational point From there to the wall, a simplified set of equations is solved, and an estimate of the instantaneous wall shear stress required to impose boundary conditions is obtained Computations performed for the plane channel, square duct, and the rotating channel flow cases gave improved results compared with existing models The additional computing time required by the model is on the order of 10-15% of the overall computing time The mean flow quantities and low-order statistics, which are of primary interest in engineering calculations, are in very good agreement with the reference data available in the literature

Journal ArticleDOI
TL;DR: In this article, the evolution of the energy density in solar wind fluctuations at MHD scales is derived on the basis of transport theories appropriate to a radially expanding solar wind, which represent a departure from the well-known WKB description, including the effects of mixing, driving by stream-stream interactions (compression and shear) and interstellar pick-up ions as well as nonisotropic MHD turbulence.
Abstract: On the basis of transport theories appropriate to a radially expanding solar wind, new results for the evolution of the energy density in solar wind fluctuations at MHD scales are derived. The models, which represent a departure from the well-known WKB description, include the effects of “mixing”, driving by stream-stream interactions (compression and shear) and interstellar pick-up ions as well as non-isotropic MHD turbulence. Magnetometer data from Voyager 1 and 2 and Pioneer 11 are compared to the turbulence-based models and close agreement is found between theory and data for a reasonable choice of parameters.

Book ChapterDOI
TL;DR: In this paper, the theoretical basis for random flight models for the trajectories of tracer particles in turbulence is reviewed, and their application to calculate dispersion in the principal types of atmospheric turbulence (stratified, vertically-inhomogeneous, Gaussian or non-Gaussian turbulence in the surface layer and above) is surveyed.
Abstract: We review the theoretical basis for, and the advantages of, random flight models for the trajectories of tracer particles in turbulence. We then survey their application to calculate dispersion in the principal types of atmospheric turbulence (stratified, vertically-inhomogeneous, Gaussian or non-Gaussian turbulence in the surface layer and above), and show that they are especially suitable for some problems (e.g., quantifying ground emissions).

Journal ArticleDOI
TL;DR: In this article, a high-speed video system was used to study the interaction between sediment particles and turbulence in the wall region of an open channel flow with both smooth and transitionally rough beds.
Abstract: A high-speed video system was used to study the interaction between sediment particles and turbulence in the wall region of an open channel flow with both smooth and transitionally rough beds. In smooth flows, particles immersed within the viscous sublayer were seen to accumulate along low-speed wall streaks; apparently due to the presence of quasi-streamwise vortices in the wall region. Larger particles did not tend to group along streaks, however their velocity was observed to respond to the streaky structure of the flow velocity in the wall region. In transitionally rough flows particle sorting was not observed. Coherent flow structures in the form of shear layers typically observed in the near-wall region interacted with sediment particles lying on the channel bottom, resulting in the particles being entrained into suspension. Although there has been some speculation that this process would not be effective in entraining particles totally immersed in the viscous sublayer, the results obtained demonstrate the opposite. The entrainment mechanism appears to be the same independent of the roughness condition of the bottom wall, smooth or transitionally rough. In the latter case, however, hiding effects tend to preclude the entrainment of particles with sizes finer than that of the roughness elements. The analysis of particle velocity during entrainment shows that the streamwise component tends to be much smaller than the local mean flow velocity, while the vertical component tends to be much larger than the local standard deviation of the vertical flow velocity fluctuations, which would indicate that such particles are responding to rather extreme flow ejection events.

Journal ArticleDOI
TL;DR: In this paper, a large eddy simulation (LES) of the incompressible Navier-Stokes equations was performed at channel flow Reynolds numbers, Reτ, of 180 and 644 (based on friction velocity and channel half width); subgridscale stresses were parametrized using the Lagrangian dynamic eddy viscosity model.
Abstract: Particle transport in fully‐developed turbulent channel flow has been investigated using large eddy simulation (LES) of the incompressible Navier–Stokes equations. Calculations were performed at channel flow Reynolds numbers, Reτ, of 180 and 644 (based on friction velocity and channel half width); subgrid‐scale stresses were parametrized using the Lagrangian dynamic eddy viscosity model. Particle motion was governed by both drag and gravitational forces and the volume fraction of the dispersed phase was small enough such that particle collisions were negligible and properties of the carrier flow were not modified. Material properties of the particles used in the simulations were identical to those in the DNS calculations of Rouson and Eaton [Proceedings of the 7th Workshop on Two‐Phase Flow Predictions (1994)] and experimental measurements of Kulick et al. [J. Fluid Mech. 277, 109 (1994)]. Statistical properties of the dispersed phase in the channel flow at Reτ=180 are in good agreement with the DNS; reas...

Journal ArticleDOI
TL;DR: In this article, an evaluation of four one-equation eddy viscosity-transport turbulence closure models as applied to three-dimensional shock wave/boundary-layer interactions is presented.
Abstract: An evaluation of four one-equation eddy viscosity-transport turbulence closure models as applied to three-dimensional shock wave/boundary-layer interactions is presented herein. Comparisons of two versions of the Baldwin-Barth model, an approach of Edwards and McRae, and a modified form of the Spalart-Allmaras model are presented for two test cases, one involving Mach 8 flow over a flat plate/sharp fin apparatus and the other involving Mach 3 flow over a cylinder-offset-cone geometry. Strengths and weaknesses of the one-equation approaches are highlighted through direct comparison with experimental data, and the effect of grid refinement is examined.

Journal ArticleDOI
TL;DR: In this paper, the effect of compressibility on mixing layers was investigated using direct numerical simulation databases and it was found that the dilatational contribution to dissipation is negligible even when eddy shocklets are observed in the flow.
Abstract: Direct numerical simulation databases have been used to study the effect of compressibility on mixing layers. The simulations cover convective Mach numbers from 0.2 to 1.2 and all contain a fully resolved turbulent energy cascade to small spatial scales. Statistical information is extracted from the databases to determine reasons for the reduced growth rate that is observed as the convective Mach number is increased. It is found that the dilatational contribution to dissipation is negligible even when eddy shocklets are observed in the flow. Also pressure-dilatation is not found to be significant. Using an accurate relation between the momentum thickness growth rate and the production of turbulence kinetic energy together with integrated equations for the Reynolds stress tensor it is shown that reduced pressure fluctuations are responsible for the changes in growth rate via the pressure–strain term. A deterministic model for the required pressure fluctuations is given based on the structure of variable-density vortices and the assumption that the limiting eddies are sonic. Simple anisotropy considerations are used to close the averaged equations. Good agreement with turbulence statistics obtained from the simulations is found.

Journal ArticleDOI
TL;DR: In this paper, the authors measured velocity and concentration profiles as well as flow resistance inside and outside the sheet-flow layer of a movable bed in steady currents using four kinds of sediment.
Abstract: Sheet-flow layer of movable bed in steady currents has been investigated experimentally, using four kinds of sediment. Velocity and concentration profiles as well as flow resistance were measured inside and outside the sheet-flow layer. The range of fall-velocity–to–friction-velocity ratio ( w /U f) achieved in the study is 0.2 0.8–1). In the suspension-mode sheet flow, the flow resistance is found to depend not only on the Shields parameter but also on a parameter involving w . Measured velocity profiles are found to follow the logarithmic law near the bed, outside the sheet-flow layer, whereas inside the sheet-flow layer, they satisfy a power law. Measured concentration profiles, on the other hand, indicate a linear variation with the distance from the bed inside the sheet flow-layer near the bed, whereas away from the bed, the familiar Rouse distribution is sat...

Journal ArticleDOI
TL;DR: In this article, a selfconsistent model for stellar turbulent convection is presented, which is similar in spirit to the CM model (Canuto \& Mazzitelli 1991) since it accounts for the full spectrum of the turbulent eddies rather than only one eddy.
Abstract: We present a selfconsistent model for stellar turbulent convection which is similar in spirit to the CM model (Canuto \& Mazzitelli 1991) since it accounts for the full spectrum of the turbulent eddies rather than only one eddy, as done in the mixing length theory (MLT). The model differs from the CM model in the treatment of the rate of energy input $n_s(k)$ from the source that generates the turbulence. In the present model, $n_s(k)$ is controlled by both the {\it source} and {\it the turbulence} it ultimately generates, thus ensuring a selfconsistent modeling of the turbulence. This improves the CM model in which $n_s(k)$ was taken to be equal to the growth rate of the {\it linear} unstable convective modes. However, since the formulation of a selfconsistent treatment is far from simple, we were forced to use a representation of the nonlinear interactions less complete than the one in the CM model. The ensuing equations were solved numerically for a wide range of convective efficiencies. The results are the convective flux, the mean square turbulent velocity, the root mean squared turbulent pressure and the turbulent viscosity. We implemented the model in the ATON stellar structure code and computed the evolution of a solar model. The results are generally similar to those of the CM model and thus quite different from the MLT. The present model requires a smaller overshoot into the upper radiative zone than does the CM model, in accord with recent empirical estimates. Application to POP II stars and comparison with the very metal-poor globular cluster M68 yields an age in the range $11\div 12$ Gyr. This is somewhat younger than the CM age, which in turn is younger than the corresponding MLT age, a result of possible cosmological interest.