Most cited
This page lists all time most cited articles for this title. Please use the publication date filters on the left if you would like to restrict this list to recently published content, for example to articles published in the last three years. The number of times each article was cited is displayed to the right of its title and can be clicked to access a list of all titles this article has been cited by.
- Cited by 853
Maximal deformation of an impacting drop
- CHRISTOPHE CLANET, CÉDRIC BÉGUIN, DENIS RICHARD, DAVID QUÉRÉ
-
- Published online by Cambridge University Press:
- 11 October 2004, pp. 199-208
-
- Article
- Export citation
-
We first study the impact of a liquid drop of low viscosity on a super-hydrophobic surface. Denoting the drop size and speed as $D_{0}$ and $U_{0}$, we find that the maximal spreading $D_{\hbox{\scriptsize\it max}}$ scales as $D_{0}\hbox{\it We}^{1/4}$ where We is the Weber number associated with the shock ($\hbox{\it We}\,{\equiv}\,\rho U_{0}^2 D_{0}/\sigma$, where $\rho$ and $\sigma$ are the liquid density and surface tension). This law is also observed to hold on partially wettable surfaces, provided that liquids of low viscosity (such as water) are used. The law is interpreted as resulting from the effective acceleration experienced by the drop during its impact. Viscous drops are also analysed, allowing us to propose a criterion for predicting if the spreading is limited by capillarity, or by viscosity.
- Cited by 849
A universal time scale for vortex ring formation
- Part of:
- MORTEZA GHARIB, EDMOND RAMBOD, KARIM SHARIFF
-
- Published online by Cambridge University Press:
- 10 April 1998, pp. 121-140
-
- Article
- Export citation
-
The formation of vortex rings generated through impulsively started jets is studied experimentally. Utilizing a piston/cylinder arrangement in a water tank, the velocity and vorticity fields of vortex rings are obtained using digital particle image velocimetry (DPIV) for a wide range of piston stroke to diameter (L/D) ratios. The results indicate that the flow field generated by large L/D consists of a leading vortex ring followed by a trailing jet. The vorticity field of the leading vortex ring formed is disconnected from that of the trailing jet. On the other hand, flow fields generated by small stroke ratios show only a single vortex ring. The transition between these two distinct states is observed to occur at a stroke ratio of approximately 4, which, in this paper, is referred to as the ‘formation number’. In all cases, the maximum circulation that a vortex ring can attain during its formation is reached at this non-dimensional time or formation number. The universality of this number was tested by generating vortex rings with different jet exit diameters and boundaries, as well as with various non-impulsive piston velocities. It is shown that the ‘formation number’ lies in the range of 3.6–4.5 for a broad range of flow conditions. An explanation is provided for the existence of the formation number based on the Kelvin–Benjamin variational principle for steady axis-touching vortex rings. It is shown that based on the measured impulse, circulation and energy of the observed vortex rings, the Kelvin–Benjamin principle correctly predicts the range of observed formation numbers.
- Cited by 848
Scaling in thermal convection: a unifying theory
- SIEGFRIED GROSSMANN, DETLEF LOHSE
-
- Published online by Cambridge University Press:
- 25 March 2000, pp. 27-56
-
- Article
- Export citation
-
A systematic theory for the scaling of the Nusselt number Nu and of the Reynolds number Re in strong Rayleigh–Bénard convection is suggested and shown to be compatible with recent experiments. It assumes a coherent large-scale convection roll (‘wind of turbulence’) and is based on the dynamical equations both in the bulk and in the boundary layers. Several regimes are identified in the Rayleigh number Ra versus Prandtl number Pr phase space, defined by whether the boundary layer or the bulk dominates the global kinetic and thermal dissipation, respectively, and by whether the thermal or the kinetic boundary layer is thicker. The crossover between the regimes is calculated. In the regime which has most frequently been studied in experiment (Ra [lsim ] 1011) the leading terms are Nu ∼ Ra1/4Pr1/8, Re ∼ Ra1/2Pr−3/4 for Pr [lsim ] 1 and Nu ∼ Ra1/4Pr−1/12, Re ∼ Ra1/2Pr−5/6 for Pr [gsim ] 1. In most measurements these laws are modified by additive corrections from the neighbouring regimes so that the impression of a slightly larger (effective) Nu vs. Ra scaling exponent can arise. The most important of the neighbouring regimes towards large Ra are a regime with scaling Nu ∼ Ra1/2Pr1/2, Re ∼ Ra1/2Pr−1/2 for medium Pr (‘Kraichnan regime’), a regime with scaling Nu ∼ Ra1/5Pr1/5, Re ∼ Ra2/5Pr−3/5 for small Pr, a regime with Nu ∼ Ra1/3, Re ∼ Ra4/9Pr−2/3 for larger Pr, and a regime with scaling Nu ∼ Ra3/7Pr−1/7, Re ∼ Ra4/7Pr−6/7 for even larger Pr. In particular, a linear combination of the ¼ and the 1/3 power laws for Nu with Ra, Nu = 0.27Ra1/4 + 0.038Ra1/3 (the prefactors follow from experiment), mimics a 2/7 power-law exponent in a regime as large as ten decades. For very large Ra the laminar shear boundary layer is speculated to break down through the non-normal-nonlinear transition to turbulence and another regime emerges. The theory presented is best summarized in the phase diagram figure 2 and in table 2.
- Cited by 848
The determination of the bulk stress in a suspension of spherical particles to order c2
- G. K. Batchelor, J. T. Green
-
- Published online by Cambridge University Press:
- 29 March 2006, pp. 401-427
-
- Article
- Export citation
-
An exact formula is obtained for the term of order c2 in the expression for the bulk stress in a suspension of force-free spherical particles in Newtonian ambient fluid, where c is the volume fraction of the spheres and c [Lt ] 1. The particles may be of different sizes, and composed of either solid or fluid of arbitrary viscosity. The method of derivation circumvents the familiar obstacle, of non-absolutely convergent integrals representing the effect of all pair interactions in which one specified particle takes part, by the judicious use of a certain quantity which is affected by the presence of distant particles in a similar way and whose mean value is known exactly. The bulk stress is in general of non-Newtonian form and depends on the statistical properties of the suspension which in turn are dependent on the type of bulk flow.
The formula contains two functions which are parameters of the flow field due to two spherical particles immersed in fluid in which the velocity gradient is uniform at infinity. One of them, p(r, t), represents the probability density for the vector r separating the centres of the two particles. The variation of p(r, t) for a moving material point in r-space due to hydrodynamic action is found in terms of a function q(r), and this gives p(r, t) explicitly over the whole of the region of r-space occupied by trajectories of one particle centre relative to another which come from infinity. In a region of closed trajectories, steady-state hydrodynamic action alone does not determine the relation between the values of p (r, t) for different material points. The function q(r) is singular when the spheres touch, and the contribution of nearly-touching spheres to the bulk stress is evidently important. Approximate numerical values of all the relevant functions are presented for the case of rigid spherical particles of uniform size.
In the case of steady pure straining motion of the suspension, all trajectories in r-space come from infinity, the suspension has isotropic structure and the stress behaviour can be represented (to order c2) in terms of an effective viscosity ${\mathop\mu\limits^{*}}$. It is estimated from the available numerical data that for a suspension of identical rigid spherical particles \[ {\mathop\mu\limits^{*}}/\mu = 1 + 2.5c + 7.6c^2, \] the error bounds on the coefficient of c2 being about ∓ 0.8. In the important case of steady simple shearing motion, there is a region of closed trajectories of one sphere centre relative to another, of infinite volume. The stress system is here not of Newtonian form, and numerical results are not obtainable until the probability, density p(r, t) can be made determinate in the region of closed trajectories by the introduction of some additional physical process, such as three-sphere encounters or Brownian motion, or by the assumption of some particular initial state.
In the analogous problem for an incompressible solid suspension it may be appropriate to assume that for many methods of manufacture p(r, t) is uniform over the accessible part of r-space, in which event the solid suspension has ‘Newtonian’ elastic behaviour and the ratio of the effective shear modulus to that of the matrix is estimated to be 1 + 2·5c + 5·2c2 for a suspension of identical rigid spheres.
- Cited by 837
The motion of bubbles in a vertical temperature gradient
- N. O. Young, J. S. Goldstein, M. J. Block
-
- Published online by Cambridge University Press:
- 28 March 2006, pp. 350-356
-
- Article
- Export citation
-
It has been observed experimentally that small bubbles in pure liquids can be held stationary or driven downwards by means of a sufficiently strong negative temperature gradient in the vertical direction. This effect is demonstrated to be due to the stresses resulting from the thermal variation of surface tension at the bubble surface. The flow field within and around the bubble is derived, and an expression for the magnitude of the temperature gradient required to hold the bubble stationary is obtained. This expression is verified experimentally.
- Cited by 835
The three-dimensional nature of boundary-layer instability
- P. S. Klebanoff, K. D. Tidstrom, L. M. Sargent
-
- Published online by Cambridge University Press:
- 28 March 2006, pp. 1-34
-
- Article
- Export citation
-
An experimental investigation is described in which principal emphasis is given to revealing the nature of the motions in the non-linear range of boundary-layer instability and the onset of turbulence. It has as its central purpose the evaluation of existing theoretical considerations and the provision of a sound physical model which can be taken as a basis for a theoretical approach. The experimental method consisted of introducing, in a two-dimensional boundary layer on a flat plate at ‘incompressible’ speeds, three-dimensional disturbances under controlled conditions using the vibrating-ribbon technique, and studying their growth and evolution using hot-wire methods. It has been definitely established that longitudinal vortices are associated with the non-linear three-dimensional wave motions. Sufficient data were obtained for an evaluation of existing theoretical approaches. Those which have been considered are the generation of higher harmonics, the interaction of the mean flow and the Reynold stress, the concave streamline curvature associated with the wave motion, the vortex loop and the non-linear effect of a three-dimensional perturbation. It appears that except for the latter they do not adequately describe the observed phenomena. It is not that they are incorrect or may not play a role in some aspect of the local behaviour, but from the over-all point of view the results suggest that it is the non-linear effect of a three-dimensional perturbation which dominates the behaviour. A principal conclusion to be drawn is that a new perspective, one that takes three-dimensionality into account, is required in connexion with boundary-layer instability. It is demonstrated that the actual breakdown of the wave motion into turbulence is a consequence of a new instability which arises in the aforementioned three-dimensional wave motion. This instability involves the generation of ‘hairpin’ eddies and is remarkably similar in behaviour to ‘inflexional’ instability. It is also shown that the results obtained from the study of controlled disturbances are equally applicable to ‘natural’ transition.
- Cited by 833
Scaling of hard thermal turbulence in Rayleigh-Bénard convection
- Bernard Castaing, Gemunu Gunaratne, François Heslot, Leo Kadanoff, Albert Libchaber, Stefan Thomae, Xiao-Zhong Wu, Stéphane Zaleski, Gianluigi Zanetti
-
- Published online by Cambridge University Press:
- 26 April 2006, pp. 1-30
-
- Article
- Export citation
-
An experimental study of Rayleigh-Bénard convection in helium gas at roughly 5 K is performed in a cell with aspect ratio 1. Data are analysed in a ‘hard turbulence’ region (4 × 107 < Ra < 6 × 1012) in which the Prandtl number remains between 0.65 and 1.5. The main observation is a simple scaling behaviour over this entire range of Ra. However the results are not the same as in previous theories. For example, a classical result gives the dimensionless heat flux, Nu, proportional to $Ra^{\frac{1}{3}}$ while experiment gives an index much closer to $\frac{2}{7}$. A new scaling theory is described. This new approach suggests scaling indices very close to the observed ones. The new approach is based upon the assumption that the boundary layer remains in existence even though its Rayleigh number is considerably greater than unity and is, in fact, diverging. A stability analysis of the boundary layer is performed which indicates that the boundary layer may be stabilized by the interaction of buoyancy driven effects and a fluctuating wind.
- Cited by 830
Kinetic theory representation of hydrodynamics: a way beyond the Navier–Stokes equation
- XIAOWEN SHAN, XUE-FENG YUAN, HUDONG CHEN
-
- Published online by Cambridge University Press:
- 27 February 2006, pp. 413-441
-
- Article
- Export citation
-
We present in detail a theoretical framework for representing hydrodynamic systems through a systematic discretization of the Boltzmann kinetic equation. The work is an extension of a previously proposed formulation. Conventional lattice Boltzmann models can be shown to be directly derivable from this systematic approach. Furthermore, we provide here a clear and rigorous procedure for obtaining higher-order approximations to the continuum Boltzmann equation. The resulting macroscopic moment equations at each level of the systematic discretization give rise to the Navier–Stokes hydrodynamics and those beyond. In addition, theoretical indications to the order of accuracy requirements are given for each discrete approximation, for thermohydrodynamic systems, and for fluid systems involving long-range interactions. All these are important for complex and micro-scale flows and are missing in the conventional Navier–Stokes order descriptions. The resulting discrete Boltzmann models are based on a kinetic representation of the fluid dynamics, hence the drawbacks in conventional higher-order hydrodynamic formulations can be avoided.
- Cited by 826
Coherent structures and turbulence
- A. K. M. Fazle Hussain
-
- Published online by Cambridge University Press:
- 21 April 2006, pp. 303-356
-
- Article
- Export citation
-
This is a personal statement on the present state of understanding of coherent structures, in particular their spatial details and dynamical significance. The characteristic measures of coherent structures are discussed, emphasizing coherent vorticity as the crucial property. We present here a general scheme for educing structures in any transitional or fully turbulent flow. From smoothed vorticity maps in convenient flow planes, this scheme recognizes patterns of the same mode and parameter size, and then phase-aligns and ensemble-averages them to obtain coherent structure measures. The departure of individual realizations from the ensemble average denotes incoherent turbulence. This robust scheme has been used to educe structures from velocity data using a rake of hot wires as well as direct numerical simulations and can educe structures using newer measurement techniques such as digital image processing. Our recent studies of coherent structures in several free shear flows are briefly reviewed. Detailed data in circular and elliptic jets, mixing layers, and a plane wake reveal that incoherent turbulence is produced at the ‘saddles’ and then advected to the ‘centres’ of the structures. The mechanism of production of turbulence in shear layers is the stretching of longitudinal vortices or ‘ribs’ which connect the predominantly spanwise ‘rolls’; the ribs induce spanwise contortions of rolls and cause mixing and dissipation, mostly at points where they connect with rolls. We also briefly discuss the role of coherent structures in aerodynamic noise generation and argue that the structure breakdown process, rather than vortex pairing, is the dominant mechanism of noise generation. The ‘cut-and-connect’ interaction of coherent structures is proposed as a specific mechanism of aerodynamic noise generation, and a simple analytical model of it shows that it can provide acceptable predictions of jet noise. The coherent-structures approach to turbulence, apart from explaining flow physics, has also enabled turbulence management via control of structure evolution and interactions. We also discuss some new ideas under investigation: in particular, helicity as a characteristic property of coherent structures.
- Cited by 824
Note on a paper of John W. Miles
- Louis N. Howard
-
- Published online by Cambridge University Press:
- 28 March 2006, pp. 509-512
-
- Article
- Export citation
-
The theorem X established by Miles in the preceding paper is here given a simpler and more general proof. Some further theoretical results concerning the stability of heterogeneous shear flows are also presented, in particular a demonstration that the complex wave velocity of any unstable mode must lie in a certain semicircle.
- Cited by 818
Radiation stress and mass transport in gravity waves, with application to ‘surf beats’
- M. S. Longuet-Higgins, R. W. Stewart
-
- Published online by Cambridge University Press:
- 28 March 2006, pp. 481-504
-
- Article
- Export citation
-
This paper studies the second-order currents and changes in mean surface level which are caused by gravity waves of non-uniform amplitude. The effects are interpreted in terms of the radiation stresses in the waves.
The first example is of wave groups propagated in water of uniform mean depth. The problem is solved first by a perturbation analysis. In two special cases the second-order currents are found to be proportional simply to the square of the local wave amplitude: (a) when the lengths of the groups are large compared to the mean depth, and (b) when the groups are all of equal length. Then the surface is found to be depressed under a high group of waves and the mass-transport is relatively negative there. In case (a) the result can be simply related to the radiation stresses, which tend to expel fluid from beneath the higher waves.
The second example considered is the propagation of waves of steady amplitude in water of gradually varying depth. Assuming no loss of energy by friction or reflexion, the wave amplitude must vary horizontally, to maintain the flux of energy constant; it is shown that this produces a negative tilt in the mean surface level as the depth diminishes. However, if the wave height is limited by breaking, the tilt is positive. The results are in agreement with some observations by Fairchild.
Lastly, the propagation of groups of waves from deep to shallow water is studied. As the mean depth decreases, so the response of the fluid to the radiation stresses tends to increase. The long waves thus generated may be reflected as free waves, and so account for the 'surf beats’ observed by Munk and Tucker.
Generalle speaking, the changes in mean sea level produced by ocean waves are comparable with those due to horizontal wind stress. It may be necessary to allow for the wave stresses in calculating wind stress coefficients.
- Cited by 815
A simple dynamical model of intermittent fully developed turbulence
- Uriel Frisch, Pierre-Louis Sulem, Mark Nelkin
-
- Published online by Cambridge University Press:
- 12 April 2006, pp. 719-736
-
- Article
- Export citation
-
We present a phenomenological model of intermittency called the β-model and related to the Novikov-Stewart (1964) model. The key assumption is that in scales ∼ l02−n only a fraction βn of the total space has an appreciable excitation. The model, the idea of which owes much to Kraichnan (1972, 1974), is dynamical in the sense that we work entirely with inertial-range quantities such as velocity amplitudes, eddy turn-over times and energy transfer. This gives more physical insight than the traditional approach based on probabilistic models of the dissipation.
The β-model leads in an elementary way to the concept of the self-similarity dimension D, a special case of Mandelbrot's (1974, 1976) ‘fractal dimension’. For three-dimensional turbulence, the correction B to the $\frac{5}{3}$ exponent of the energy spectrum is equal to $\frac{1}{3}(3 - D)$ and is related to the exponent μ of the dissipation correlation function by $B = \frac{1}{3}\mu $ (0.17 for the currently accepted value). This is a borderline case of the Mandelbrot inequality $B \leqslant \frac{1}{3}\mu $. It is shown in the appendix that this inequality may be derived from the Navier-Stokes equation under the strong, but plausible, assumption that the inertial-range scaling laws for second- and fourth-order moments have the same viscous cut-off.
The predictions of the β-model for the spectrum and for higher-order statistics are in agreement with recent conjectures based on analogies with critical phenomena (Nelkin 1975) but generally diasgree with the 1962 Kolmogorov lognormal model. However, the sixth-order structure function 〈δv6(l)〉 and the dissipation correlation function 〈ε(r)ε(r + 1)〉 are related by \[ \langle \delta v^6(l)\rangle /l^2\sim\langle \epsilon({\bf r})\epsilon({\bf r}+{\bf l})\rangle \] in both models. We conjecture that this relation is model independent.
Finally, some possible directions for further numerical and experimental work on intermittency are indicated.
- Cited by 804
Wave formation in laminar flow down an inclined plane
- T. Brooke Benjamin
-
- Published online by Cambridge University Press:
- 28 March 2006, pp. 554-573
-
- Article
- Export citation
-
This paper deals theoretically with a problem of hydrodynamic stability characterized by small values of the Reynolds number R. The primary flow whose stability is examined consists of a uniform laminar stream of viscous liquid running down an inclined plane under the action of gravity, being bounded on one side by a free surface influenced by surface tension. The problem thus has a direct bearing on the properties of thin liquid films such as have important uses in chemical engineering.
Numerous experiments in the past have shown that in flow down a wall the stream is noticeably agitated by waves except when R is quite small; on a vertical water film, for instance, waves may be observed until R is reduced to some value rather less than 10. The present treatment is accordingly based on methods of approximation suited to fairly low values of R, and thereby avoids the severe mathematical difficulties usual in stability problems at high R. The formulation of the problem resembles that given by Yih (1954); but the method of solution differs from his, and the respective results are in conflict. In particular, there is dis-agreement over the matter of the stability of a strictly vertical stream at very small R. In contrast with the previous conclusions, it is shown here that the flow is always unstable: that is, a class of undamped waves exists for all finite values of R. However, the rates of amplification of unstable waves are shown to become very small when R is made fairly small, and their wavelengths to become very large; this provides a satisfactory explanation for the apparent absence of waves in some experimental observations, and also for the wide scatter among existing estimates of the ‘quasi-critical’ value of R below which waves are undetectable. In view of the controversial nature of these results, emphasis is given to various points of agreement between the present work and the established theory of roll waves; the latter theory gives a clear picture of the physical mechanism of wave formation on gravitational flows, and in its light the results obtained here appear entirely reasonable.
The conditions governing neutral stability are worked out to the third order in a parameter which is shown to be small; but a less accurate approximation is then justified as an adequate basis for an easily workable theory providing a ready check with experiment, This theory is used to predict the value of R at which observable waves should first develop on a vertical water film, and also the length and velocity of the waves. These three predictions are compared with the experimental results found by Binnie (1957), and are substantially confirmed.
- Cited by 801
Oblique and parallel modes of vortex shedding in the wake of a circular cylinder at low Reynolds numbers
- C. H. K. Williamson
-
- Published online by Cambridge University Press:
- 26 April 2006, pp. 579-627
-
- Article
- Export citation
-
Two fundamental characteristics of the low-Reynolds-number cylinder wake, which have involved considerable debate, are first the existence of discontinuities in the Strouhal-Reynolds number relationship, and secondly the phenomenon of oblique vortex shedding. The present paper shows that both of these characteristics of the wake are directly related to each other, and that both are influenced by the boundary conditions at the ends of the cylinder, even for spans of hundreds of diameters in length. It is found that a Strouhal discontinuity exists, which is not due to any of the previously proposed mechanisms, but instead is caused by a transition from one oblique shedding mode to another oblique mode. This transition is explained by a change from one mode where the central flow over the span matches the end boundary conditions to one where the central flow is unable to match the end conditions. In the latter case, quasi-periodic spectra of the velocity fluctuations appear; these are due to the presence of spanwise cells of different frequency. During periods when vortices in neighbouring cells move out of phase with each other, ‘vortex dislocations’ are observed, and are associated with rather complex vortex linking between the cells. However, by manipulating the end boundary conditions, parallel shedding can be induced, which then results in a completely continuous Strouhal curve. It is also universal in the sense that the oblique-shedding Strouhal data (Sθ) can be collapsed onto the parallel-shedding Strouhal curve (S0) by the transformation, S0 = Sθ/cosθ, where θ is the angle of oblique shedding. Close agreement between measurements in two distinctly different facilities confirms the continuous and universal nature of this Strouhal curve. It is believed that the case of parallel shedding represents truly two-dimensional shedding, and a comparison of Strouhal frequency data is made with several two-dimensional numerical simulations, yielding a large disparity which is not clearly understood. The oblique and parallel modes of vortex shedding are both intrinsic to the flow over a cylinder, and are simply solutions to different problems, because the boundary conditions are different in each case.
- Cited by 798
Some specific features of atmospheric tubulence
- A. M. Oboukhov
-
- Published online by Cambridge University Press:
- 28 March 2006, pp. 77-81
-
- Article
- Export citation
-
The spectrum of atmospheric turbulence is very broad by comparison with spectra in wind tunnels. We introduce the notion of small-scale and large-scale turbulence. Small-scale turbulence consists of a set of disturbances, the scales of which do not exceed the distance to the wall and for which the hypothesis of three-dimensional isotropy is valid in a certain rough approximation. Large-scale turbulence is essentially anisotropic; the horizontal scale in the atmosphere is much larger than the vertical one, the latter being confined to a certain characteristic height H. The horizontal scale varies widely according to the external conditions and characteristics of the medium.
- Cited by 792
Collapse of an initially spherical vapour cavity in the neighbourhood of a solid boundary
- Milton S. Plesset, Richard B. Chapman
-
- Published online by Cambridge University Press:
- 29 March 2006, pp. 283-290
-
- Article
- Export citation
-
Vapour bubble collapse problems lacking spherical symmetry are solved here using a numerical method designed especially for these problems. Viscosity and compressibility in the liquid are neglected. Two specific cases of initially spherical bubbles collapsing near a plane solid wall were simulated: a bubble initially in contact with the wall, and a bubble initially half its radius from the wall at the closest point. It is shown that the bubble develops a jet directed towards the wall rather early in the collapse history. Free surface shapes and velocities are presented at various stages in the collapse. Velocities are scaled like (Δp/ρ)½ where ρ is the density of the liquid and Δp is the constant difference between the ambient liquid pressure and the pressure in the cavity. For \[ \Delta p/\rho = 10^6 {\rm cm}^2/\sec^2 \approx 1\, \hbox{atm/density of water} \] the jet had a speed of about 130m/sec in the first case and 170m/sec in the second when it struck the opposite side of the bubble. Such jet velocities are of a magnitude which can explain cavitation damage. The jet develops so early in the bubble collapse history that compressibility effects in the liquid and the vapour are not important.
- Cited by 789
Calculations of the development of an undular bore
- D. H. Peregrine
-
- Published online by Cambridge University Press:
- 28 March 2006, pp. 321-330
-
- Article
- Export citation
-
If a long wave of elevation travels in shallow water it steepens and forms a bore. The bore is undular if the change in surface elevation of the wave is less than 0·28 of the original depth of water. This paper describes the growth of an undular bore from a long wave which forms a gentle transition between a uniform flow and still water. A physical account of its development is followed by the results of numerical calculations. These use finite-difference approximations to the partial differential equations of motion. The equations of motion are of the same order of approximation as is necessary to derive the solitary wave. The results are in general agreement with the available experimental measurements.
- Cited by 787
Turbulent entrainment in stratified flows
- T. H. Ellison, J. S. Turner
-
- Published online by Cambridge University Press:
- 28 March 2006, pp. 423-448
-
- Article
- Export citation
-
When a fluid which is lighter than its surroundings is emitted by a source under a sloping roof (or a heavier fluid from a source on a sloping floor), it may flow as a relatively thin turbulent layer. The motion of this layer is governed by the rate at which it entrains the ambient fluid. A theory is presented in which it is assumed that the entrainment is proportional to the velocity of the layer multiplied by an empirical function, E(Ri), of the overall Richardson number for the layer defined by Ri = g(ρa - ρ) h/ρaV2. This theory predicts that in most practical cases the layer will rapidly attain an equilibrium state in which Ri does not vary with distance downstream, and the gravitational force on the layer is just balanced by the drag due to entrainment together with friction on the floor or roof.
Two series of laboratory experiments are described from which E(Ri) can be determined. In the first, the spread of a surface jet of fluid lighter than that over which it is flowing is measured; in the second, a study is made of the flow of a heavy liquid down the sloping floor of a channel. These experiments show that E falls off rapidly as Ri increases and is probably negligible when Ri is more than about 0·8.
The theoretical and experimental results allow predictions to be made of flow velocities once the rate of supply of density difference is known. An estimate is also given of the uniform velocity which the ambient fluid must possess in order to cause the motion of the layer to be reversed.
- Cited by 785
Turbulence: the filtering approach
- M. Germano
-
- Published online by Cambridge University Press:
- 26 April 2006, pp. 325-336
-
- Article
- Export citation
-
Explicit or implicit filtered representations of chaotic fields like spectral cut-offs or numerical discretizations are commonly used in the study of turbulence and particularly in the so-called large-eddy simulations. Peculiar to these representations is that they are produced by different filtering operators at different levels of resolution, and they can be hierarchically organized in terms of a characteristic parameter like a grid length or a spectral truncation mode. Unfortunately, in the case of a general implicit or explicit filtering operator the Reynolds rules of the mean are no longer valid, and the classical analysis of the turbulence in terms of mean values and fluctuations is not so simple.
In this paper a new operatorial approach to the study of turbulence based on the general algebraic properties of the filtered representations of a turbulence field at different levels is presented. The main results of this analysis are the averaging invariance of the filtered Navier—Stokes equations in terms of the generalized central moments, and an algebraic identity that relates the turbulent stresses at different levels. The statistical approach uses the idea of a decomposition in mean values and fluctuations, and the original turbulent field is seen as the sum of different contributions. On the other hand this operatorial approach is based on the comparison of different representations of the turbulent field at different levels, and, in the opinion of the author, it is particularly fitted to study the similarity between the turbulence at different filtering levels. The best field of application of this approach is the numerical large-eddy simulation of turbulent flows where the large scale of the turbulent field is captured and the residual small scale is modelled. It is natural to define and to extract from the resolved field the resolved turbulence and to use the information that it contains to adapt the subgrid model to the real turbulent field. Following these ideas the application of this approach to the large-eddy simulation of the turbulent flow has been produced (Germano et al. 1991). It consists in a dynamic subgrid-scale eddy viscosity model that samples the resolved scale and uses this information to adjust locally the Smagorinsky constant to the local turbulence.
- Cited by 785
Solute transport in heterogeneous porous formations
- Gedeon Dagan
-
- Published online by Cambridge University Press:
- 20 April 2006, pp. 151-177
-
- Article
- Export citation
-
Solute transport in porous formations is governed by the large-scale heterogeneity of hydraulic conductivity. The two typical lengthscales are the local one (of the order of metres) and the regional one (of the order of kilometres). The formation is modelled as a random fixed structure, to reflect the uncertainty of the space distribution of conductivity, which has a lognormal probability distribution function. A first-order perturbation approximation, valid for small log-conductivity variance, is used in order to derive closed-form expressions of the Eulerian velocity covariances for uniform average flow. The concentration expectation value is determined by using a similar approximation, and it satisfies a diffusion equation with time-dependent apparent dispersion coefficients. The longitudinal coefficients tend to constant values in both two- and three-dimensional flows only after the solute body has travelled a few tens of conductivity integral scales. This may be an exceedingly large distance in many applications for which the transient stage prevails. Comparison of theoretical results with recent field experimental data is quite satisfactory.
The variance of the space-averaged concentration over a volume V may be quite large unless the lengthscale of the initial solute body or of V is large compared with the conductivity integral scale. This condition is bound to be obeyed for transport at the local scale, in which case the concentration may be assumed to satisfy the ergodic hypothesis. This is not generally the case at the regional scale, and the solute concentration is subjected to large uncertainty. The usefulness of the prediction of the concentration expectation value is then quite limited and the dispersion coefficients become meaningless.
In the second part of the study, the influence of knowledge of the conductivity and head at a set of points upon transport is examined. The statistical moments of the velocity and concentration fields are computed for a subensemble of formations and for conditional probability distribution functions of conductivity and head, with measured values kept fixed at the set of measurement points. For conditional statistics the velocity is not stationary, and its mean and variance vary throughout the space, even if its unconditional mean and variance are constant. The main aim of the analysis is to examine the reduction of concentration coefficient of variation, i.e. of its uncertainty, by conditioning. It is shown that measurements of transmissivity on a grid of points can be effective in reducing concentration variance, provided that the distance between the points is smaller than two conductivity integral scales. Head conditioning has a lesser effect upon variance reduction.