Showing papers in "Journal of Fluid Mechanics in 1974"
TL;DR: In this article, Spark shadow pictures and measurements of density fluctuations suggest that turbulent mixing and entrainment is a process of entanglement on the scale of the large structures; some statistical properties of the latter are used to obtain an estimate of entrainedment rates, and large changes of the density ratio across the mixing layer were found to have a relatively small effect on the spreading angle.
Abstract: Plane turbulent mixing between two streams of different gases (especially nitrogen and helium) was studied in a novel apparatus Spark shadow pictures showed that, for all ratios of densities in the two streams, the mixing layer is dominated by large coherent structures High-speed movies showed that these convect at nearly constant speed, and increase their size and spacing discontinuously by amalgamation with neighbouring ones The pictures and measurements of density fluctuations suggest that turbulent mixing and entrainment is a process of entanglement on the scale of the large structures; some statistical properties of the latter are used to obtain an estimate of entrainment rates Large changes of the density ratio across the mixing layer were found to have a relatively small effect on the spreading angle; it is concluded that the strong effects, which are observed when one stream is supersonic, are due to compressibility effects, not density effects, as has been generally supposed
3,339 citations
TL;DR: In this article, it was shown that Kolmogorov's third hypothesis is logically inconsistent, save under assumptions that are extreme and unlikely, and a widely used justification of lognormality due to Yaglom and based on probabilistic argument involving a self-similar cascade, will also be discussed.
Abstract: Kolmogorov’s “third hypothesis” asserts that in intermittent turbulence the average \( \bar \varepsilon \) of the dissipation e, taken over any domain D, is ruled by the lognormal probability distribution. This hypothesis will be shown to be logically inconsistent, save under assumptions that are extreme and unlikely. A widely used justification of lognormality due to Yaglom and based on probabilistic argument involving a self-similar cascade, will also be discussed. In this model, lognormality indeed applies strictly when D is “an eddy,” typically a three-dimensional box embedded in a self-similar hierarchy, and may perhaps remain a reasonable approximation when D consists of a few such eddies. On the other hand, the experimental situation is better described by considering averages taken over essentially one-dimensional domains D.
1,567 citations
TL;DR: A mixing layer is formed by bringing two streams of water, moving at different velocities, together in a lucite-walled channel as mentioned in this paper, where dye is injected between the two streams just before they are brought together, marking the vorticitycarrying fluid.
Abstract: A mixing layer is formed by bringing two streams of water, moving at different velocities, together in a lucite-walled channel. The Reynolds number, based on the velocity difference and the thickness of the shear layer, varies from about 45, where the shear layer originates, to about 850 at a distance of 50 cm. Dye is injected between the two streams just before they are brought together, marking the vorticity-carrying fluid. Unstable waves grow, and fluid is observed to roll up into discrete two-dimensional vortical structures. These turbulent vortices interact by rolling around each other, and a single vortical structure, with approximately twice the spacing of the former vortices, is formed. This pairing process is observed to occur repeatedly, controlling the growth of the mixing layer. A simple model of the mixing layer contains, as the important elements controlling growth, the degree of non-uniformity in the vortex train and the ‘lumpiness’ of the vorticity field.
1,335 citations
TL;DR: In this article, the Segre-Silberberg effect of inertia-induced lateral migration of a neutrally buoyant rigid sphere in a Newtonian fluid is studied theoretically for simple shear flow and for two-dimensional Poiseuille flow.
Abstract: The familiar Segre-Silberberg effect of inertia-induced lateral migration of a neutrally buoyant rigid sphere in a Newtonian fluid is studied theoretically for simple shear flow and for two-dimensional Poiseuille flow. It is shown that the spheres reach a stable lateral equilibrium position independent of the initial position of release. For simple shear flow, this position is midway between the walls, whereas for Poiseuille flow, it is 0·6 of the channel half-width from the centre-line. Particle trajectories are calculated in both cases and compared with available experimental data. Implications for the measurement of the rheological properties of a dilute suspension of spheres are discussed.
790 citations
TL;DR: In this article, a fluid-fluid interface that joins a solid surface forms a common line and if the common line moves along the solid, a mutual displacement process is involved and is studied here.
Abstract: A fluid-fluid interface that joins a solid surface forms a common line. If the common line moves along the solid, a mutual displacement process is involved and is studied here. Some simple experiments motivate the formulation of the basic assumption of the analysis. The basic assumption is a formalization of the idea that the fluid-fluid interface rolls on or unrolls off the solid. This forms an axiom for the mostly kinematical analysis that follows. The predictions are tested through a series of qualitative experiments. The role of the no-slip boundary condition at the solid surface is discussed.
785 citations
TL;DR: In this article, the change in energy dissipation due to a small hump on a body in a uniform steady flow is calculated, and the result is used in conjunction with the variational methods of optimal control to obtain the optimality conditions for four minimum-drag problems of fluid mechanics.
Abstract: In this paper, the change in energy dissipation due to a small hump on a body in a uniform steady flow is calculated. The result is used in conjunction with the variational methods of optimal control to obtain the optimality conditions for four minimum-drag problems of fluid mechanics. These conditions imply that the unit-area profile of smallest drag has a front end shaped like a wedge of angle 90°.
699 citations
TL;DR: In this article, a review of the geophysical information and the fluid dynamics of convection in a Boussinesq fluid of infinite Prandtl number is presented and analyzed in terms of simple physical models.
Abstract: Plate tectonics provides a remarkably accurate kinematic description of the motion of the earth's crust but a fully dynamical theory requires an understanding of convection in the mantle. Thus the properties of plates and of the mantle must be related to a systematic study of convection. This paper reviews both the geophysical information and the fluid dynamics of convection in a Boussinesq fluid of infinite Prandtl number. Numerical experiments have been carried out on several simple two-dimensional models, in which convection is driven by imposed horizontal temperature gradients or else by heating either internally or from below. The results are presented and analysed in terms of simple physical models. Although the computations are highly idealized and omit variation of viscosity and other major features of mantle convection, they can be related to geophysical measurements. In particular, the external gravity field depends on changes in surface elevation; this suggests an observational means of investigating convection in the upper mantle.
662 citations
TL;DR: In this paper, the Strouhal number as a function of Reynolds number measured by Moller (1938) has been confirmed using water flow and the lower critical Reynolds number, first reported by Cometta (1957), was found to be Re = 6 × 103.
Abstract: Vortex shedding from spheres has been studied in the Reynolds number range 400 < Re < 5 × 106. At low Reynolds numbers, i.e. up to Re = 3 × 103, the values of the Strouhal number as a function of Reynolds number measured by Moller (1938) have been confirmed using water flow. The lower critical Reynolds number, first reported by Cometta (1957), was found to be Re = 6 × 103. Here a discontinuity in the relationship between the Strouhal and Reynolds numbers is obvious. From Re = 6 × 103 to Re = 3 × 105 strong periodic fluctuations in the wake flow were observed. Beyond the upper critical Reynolds number (Re = 3.7 × 105) periodic vortex shedding could not be detected by the present measurement techniques.The hot-wire measurements indicate that the signals recorded simultaneously at different positions on the 75° circle (normal to the flow) show a phase shift. Thus it appears that the vortex separation point rotates around the sphere. An attempt is made to interpret this experimental evidence.
477 citations
TL;DR: In this paper, a simple model for the experimentally observed instability of the vortex ring to azimuthal bending waves of wavelength comparable with the core size is presented, and short-wave instabilities are discussed for both the ring and the vortex pair.
Abstract: A simple model for the experimentally observed instability of the vortex ring to azimuthal bending waves of wavelength comparable with the core size is presented. Short-wave instabilities are discussed for both the vortex ring and the vortex pair. Instability for both the ring and the pair is predicted to occur whenever the self-induced rotation of waves on the filament passes through zero. Although this does not occur for the first radial bending mode of a vortex filament, it is shown to be possible for bending modes with a more complex radial structure with at least one node at some radius within the core. The previous work of Widnall & Sullivan (1973) is discussed and their experimental results are compared with the predictions of the analysis presented here.
391 citations
TL;DR: In this paper, a hot-film anemometer measurement was carried out in a fully developed turbulent channel flow with a viscous sublayer and the results showed that the streamwise velocity fluctuations decreased at a higher rate than the mean velocity; in the region y+ [lsim ] 0·1, these fluctuations vanished at the same rate as the average velocity.
Abstract: Hot-film anemometer measurements have been carried out in a fully developed turbulent channel flow. An oil channel with a thick viscous sublayer was used, which permitted measurements very close to the wall. In the viscous sublayer between y+ ≃ 0·1 and y+ = 5, the streamwise velocity fluctuations decreased at a higher rate than the mean velocity; in the region y+ [lsim ] 0·1, these fluctuations vanished at the same rate as the mean velocity.The streamwise velocity fluctuations u observed in the viscous sublayer and the fluctuations (∂u/∂y)0 of the gradient at the wall were almost identical in form, but the fluctuations of the gradient at the wall were found to lag behind the velocity fluctuations with a lag time proportional to the distance from the wall. Probability density distributions of the streamwise velocity fluctuations were measured. Furthermore, measurements of the skewness and flatness factors made by Kreplin (1973) in the same flow channel are discussed. Measurements of the normal velocity fluctuations v at the wall and of the instantaneous Reynolds stress −ρuv were also made. Periods of quiescence in the − ρuv signal were observed in the viscous sublayer as well as very active periods where ratios of peak to mean values as high as 30:1 occurred.
380 citations
TL;DR: In this article, the stability of two-dimensional convection rolls with respect to three-dimensional disturbances is analyzed, and it is found that convection roll are unstable for Prandtl numbers less than about 5, where the instability is caused by momentum advection terms in the equations of motion.
Abstract: Steady solutions in the form of two-dimensional rolls are obtained for convection in a horizontal layer of fluid heated from below as a function of the Rayleigh and Prandtl numbers. Rigid boundaries of infinite heat conductivity are assumed. The stability of the two-dimensional rolls with respect to three-dimensional disturbances is analysed. It is found that convection rolls are unstable for Prandtl numbers less than about 5 with respect to an oscillatory instability investigated earlier by Busse (1972) for the case of free boundaries. Since the instability is caused by the momentum advection terms in the equations of motion the Rayleigh number for the onset of instability increases strongly with Prandtl number. Good agreement with various experimental observations is found.
TL;DR: In this paper, Stokes' infinitesimal-wave expansion for steady progressive free-surface waves has been extended to high order using a computer to perform the coefficient arithmetic, which is valid for any finite value of the wavelength and solutions of high accuracy can be obtained for most values of the wave height and water depth.
Abstract: Stokes' infinitesimal-wave expansion for steady progressive free-surface waves has been extended to high order using a computer to perform the coefficient arithmetic. Stokes' expansion has been found to be incapable of yielding the highest wave for any value of the water depth since convergence is limited by a square-root branch-point some distance short of the maximum. By reformulating the problem using a different independent parameter, the highest waves are obtained correctly. Series summation and analytic continuation are facilitated by the use of Pade approximants. The method is valid in principle for any finite value of the wavelength and solutions of high accuracy can be obtained for most values of the wave height and water depth.
TL;DR: In this article, the authors examined consistency and uniqueness questions raised by both the 1941 and 1962 Kolmogorov inertial-range theories and showed that even scalar nonlinear cascade processes need not yield asymptotic log-normality.
Abstract: Consistency and uniqueness questions raised by both the 1941 and 1962 Kolmogorov inertial-range theories are examined. The 1941 theory, although unlikely from the viewpoint of vortex-stretching physics, is not ruled out just because the dissipation fluctuates; but self-consistency requires that dissipation fluctuations be confined to dissipation-range scales by a spacewise mixing mechanism. The basic idea of the 1962 theory is a self-similar cascade mechanism which produces systematically increasing intermittency with a decrease of scale size. This concept in itself requires neither the third Kolmogorov hypothesis (log-normality of locally averaged dissipation) nor the first hypothesis (universality of small-scale statistics as functions of scale-size ratios and locally averaged dissipation). It does not even imply that the inertial range exhibits power laws. A central role for dissipation seems arbitrary since conservation alone yields no simple relation between the local dissipation rate and the corresponding proper inertial-range quantity: the local rate of energy transfer. A model rate equation for the evolution of probability densities is used to illustrate that even scalar nonlinear cascade processes need not yield asymptotic log-normality. The approximate experimental support for Kolmogorov's hypothesis takes on added significance in view of the wide variety of a priori admissible alternatives.If the Kolmogorov law is asymptotically valid, it is argued that the value of μ depends on the details of the nonlinear interaction embodied in the Navier–Stokes equation and cannot be deduced from overall symmetries, invariances and dimensionality. A dynamical equation is exhibited which has the same essential invariances, symmetries, dimensionality and equilibrium statistical ensembles as the Navier–Stokes equation but which has radically different inertial-range behaviour.
TL;DR: In this paper, a 3: 1 symmetric expansion in a duct with an aspect ratio of 9·2: 1 downstream of the expansion is reported. But the velocity profiles were in good agreement with those obtained by solving the two-dimensional momentum equation, although there were substantial threedimensional effects in the vicinity of the separation regions.
Abstract: Flow visualization and laser-anemometry measurements are reported in the flow downstream of a plane 3: 1 symmetric expansion in a duct with an aspect ratio of 9·2: 1 downstream of the expansion. The flow was found to be markedly dependent on Reynolds number, and strongly three-dimensional even well away from the channel corners except at the lowest measurable velocities. The measurements at a Reynolds number of 56 indicated that the separation regions behind each step were of equal length. Symmetric velocity profiles existed from the expansion to a fully developed, parabolic profile far downstream, although there were substantial three-dimensional effects in the vicinity of the separation regions. The velocity profiles were in good agreement with those obtained by solving the two-dimensional momentum equation. At a Reynolds number of 114, the two separation regions were of different lengths, leading to asymmetric velocity profiles; three dimensional effects were much more pronounced. At a Reynolds number of 252, a third separation zone was found on one wall, downstream of the smaller of the two separation zones adjacent to the steps. As at the lower Reynolds numbers, the flow was very stable. At higher Reynolds numbers the flow became less stable and periodicity became increasingly important in the main stream; this was accompanied by a highly disturbed fluid motion in the separation zones, as the flow tended towards turbulence.
TL;DR: In this paper, the stability of small travelling-wave disturbances in the flow over a flat plate is discussed and an iterative method is used to generate an asymptotic series solution in inverse powers of the Reynolds number Rx = Ux/v to the power one half.
Abstract: The stability of small travelling-wave disturbances in the flow over a flat plate is discussed. An iterative method is used to generate an asymptotic series solution in inverse powers of the Reynolds number Rx = Ux/v to the power one half. The neutral-stability boundaries given by the first two terms of this series are obtained and compared with experimental data. It is shown that the parallel flow approximation leads to a valid solution at very large Reynolds numbers.
TL;DR: The inviscid stability of swirling flows with mean velocity profiles similar to that obtained by Batchelor (1964) for a trailing vortex from an aircraft is studied with respect to infinitesimal non-axisymmetric disturbances as discussed by the authors.
Abstract: The inviscid stability of swirling flows with mean velocity profiles similar to that obtained by Batchelor (1964) for a trailing vortex from an aircraft is studied with respect to infinitesimal non-axisymmetric disturbances. The flow is characterized by a swirl parameter q involving the ratio of the magnitude of the maximum swirl velocity to that of the maximum axial velocity. It is found that, as the swirl is continuously increased from zero, the disturbances die out quickly for a small value of q if n = 1 (n is the azimuthal wavenumber of the Fourier disturbance of type exp{i(αx + nϕ − αct)}); but for negative values of n, the amplification rate increases and then decreases, falling to negative values at q slightly greater than 1·5 for n = −1. The maximum amplification rate increases for increasingly negative n up to n = −6 (the highest mode investigated), and corresponds to q ≃ 0·85. The applicability of these results to attempts at destabilizing vortices is briefly discussed.
TL;DR: In this paper, the problem of natural convection in a cavity of small aspect ratio with differentially heated end walls is considered, and it is shown by use of matched asymptotic expansions that the flow consists of two distinct regimes : a parallel flow in the core region and a second, non-parallel flow near the ends of the cavity.
Abstract: The problem of natural convection in a cavity of small aspect ratio with differentially heated end walls is considered. It is shown by use of matched asymptotic expansions that the flow consists of two distinct regimes : a parallel flow in the core region and a second, non-parallel flow near the ends of the cavity. A solution valid at all orders in the aspect ratio A is found for the core region, while the first several terms of the appropriate asymptotic expansion are obtained for the end regions. Parametric limits of validity for the parallel flow structure are discussed. Asymptotic expressions for the Nusselt number and the single free parameter of the parallel flow solution, valid in the limit as A → 0, are derived.
TL;DR: In this article, the effect of surface roughness on the flow past spheres has been investigated over the Reynolds number range 5 × 104 < Re < 6 × 106, where the Strouhal number for each of the various roughness conditions was equal to its value for a smooth sphere.
Abstract: The effect of surface roughness on the flow past spheres has been investigated over the Reynolds number range 5 × 104 < Re < 6 × 106. The drag coefficient has been determined as a function of the Reynolds number for five surface roughnesses. With increasing roughness parameter the critical Reynolds number decreases. At the same time the transcritical drag coefficient rises, having a maximum value of 0·4.The vortex shedding frequency has been measured under subcritical flow conditions. It was found that the Strouhal number for each of the various roughness conditions was equal to its value for a smooth sphere. Beyond the critical Reynolds number no prevailing shedding frequency could be detected by the measurement techniques employed.The drag coefficient of a sphere under the blockage conditions 0·5 < ds/dt < 0·92 has been determined over the Reynolds number range 3 × 104 < Re < 2 × 106. Increasing blockage causes an increase in both the drag coefficient and the critical Reynolds number. The characteristic quantities were referred to the flow conditions in the smallest cross-section between sphere and tube. In addition the effect of the turbulence level on the flow past a sphere under various blockage conditions was studied.
TL;DR: In this paper, a finite element program suitable for solving incompressible, viscous free surface problems in steady axisymmetric or plane flows is presented. But the authors do not consider the non-Newtonian flow, non-zero Reynolds numbers, and transient flow.
Abstract: : The authors discuss the creation of a finite element program suitable for solving incompressible, viscous free surface problems in steady axisymmetric or plane flows. For convenience in extending program capability to non-Newtonian flow, non-zero Reynolds numbers, and transient flow, a Galerkin formulation of the governing equations is chosen, rather than an extremum principle. The resulting program is used to solve the Newtonian die-swell problem for creeping jets free of surface tension constraints. The authors conclude that a Newtonian jet expands about 13%, in substantial agreement with experiments made with both small finite Reynolds numbers and small ratios of surface tension to viscous forces. The solutions to the related stick-slip problem and the tube inlet problem, both of which also contain stress singularities, are also given. (Modified author abstract)
TL;DR: In this paper, it was shown that the super-critically stable, finite amplitude, long, monochromatic wave obtained by Lin (1969, 1970, 1971) is stable to side-band disturbances under modal interaction if the bandwidth is less in magnitude than to the ratio of the amplitude to the film thickness.
Abstract: The nonlinear stability of a viscous film flowing steadily down an inclined plane is investigated by the method of multiple scales. It is shown that the super-critically stable, finite amplitude, long, monochromatic wave obtained by Lin (1969, 1970, 1971) is stable to side-band disturbances under modal interaction if the bandwidth is less in magnitude than to the ratio of the amplitude to the film thickness. Near the upper branch of the linear neutral-stability curve where the amplification rate ci is O(e2), the nonlinear evolution of initially infinitesimal waves of a finite bandwidth is shown to obey the Landau-Stuart equation, Near the lower branch of the neutral curve, the nonlinear evolution is stronger. An equation is derived for describing this strong nonlinear development of relatively long waves. In practice, disturbance of this type clusters in the form of a hump which cannot be constructed only by the first few harmonics.
TL;DR: In this article, it was shown that the surface wind drift in the ocean substantially reduces the maximum wave height ξx and wave orbital velocity that can be attained before breaking, where q is the magnitude of the surface drift at the point where the wave profile crosses the mean water level and c is the wave speed.
Abstract: It is shown that the surface wind drift in the ocean substantially reduces the maximum wave height ξx and wave orbital velocity that can be attained before breaking. If q is the magnitude of the surface drift at the point where the wave profile crosses the mean water level and c is the wave speed, then
\[
\zeta_{\max} = \frac{c^2}{2g}\bigg(1-\frac{q}{c}\bigg)^2.
\]
Incipient breaking in a steady wave train is characterized by the occurrence of stagnation points at wave crests, but not necessarily by discontinuities in slope. After breaking, there is in the mean flow a stagnation point relative to the wave profile near the crest of the broken wave, on one side of which the water tumbles forward and behind which it recedes more smoothly to the rear. Some simple flow visualization studies indicate the general extent of the wake behind the breaking region.
TL;DR: In this paper, the Korteweg-de Vries (KdV) equation is tested experimentally as a model for moderate amplitude waves propagating in one direction in relatively shallow water of uniform depth.
Abstract: The Korteweg-de Vries (KdV) equation is tested experimentally as a model for moderate amplitude waves propagating in one direction in relatively shallow water of uniform depth. For a wide range of initial data, comparisons are made between the asymptotic wave forms observed and those predicted by the theory in terms of the number of solitons that evolve, the amplitude of the leading soliton, the asymptotic shape of the wave and other qualitative features. The KdV equation is found to predict accurately the number of evolving solitons and their shapes for initial data whose asymptotic characteristics develop in the test section of the wave tank. The accuracy of the leading-soliton amplitudes computed by the KdV equation could not be conclusively tested owing to the viscous decay of the measured wave amplitudes; however, a procedure is presented for estimating the decay in amplitude of the leading wave. Computations suggest that the KdV equation predicts the amplitude of the leading soliton to within the expected error due to viscosity (12%) when the non-decayed amplitudes are less than about a quarter of the water depth. Indeed, agreement to within about 20% is observed over the entire range of experiments examined, including those with initial data for which the non-decayed amplitudes of the leading soliton exceed half the fluid depth.
TL;DR: In this paper, a calculation procedure for three-dimensional parabolic flows is applied to predict the velocity and temperature fields in helically coiled pipes, where the curvature produces a secondary flow and causes departures from the symmetric velocity profile of Poiseuille flow.
Abstract: A calculation procedure for three-dimensional parabolic flows is applied to predict the velocity and temperature fields in helically coiled pipes. The curvature produces a secondary flow and causes departures from the symmetric velocity profile of Poiseuille flow. Predictions are presented of flow and heat transfer in the developing and fully developed regions. Comparisons of the developing and fully developed velocity profiles with experimental data exhibit good agreement. The development of the wall temperature for the case of axially uniform heat flux with an isothermal periphery has been compared with experimental data and the agreement is good. Predictions for fully developed temperature profiles and heat-transfer coefficients also exhibit good agreement with experimental data. Effects of the Dean number on the friction factor and heat transfer are presented.
TL;DR: In this paper, the strength and spacing of vortices in the wake of a circular cylinder have been obtained for conditions under which the body undergoes lateral vibrations, and an inverse relation between the initial circulation and the length of the vortex formation region was obtained for cylinder oscillations of up to 50% of a diameter, at vibration frequencies both above and below the Strouhal shedding frequency.
Abstract: The strength (initial circulation) and spacing of vortices in the wake of a circular cylinder have been obtained for conditions under which the body undergoes lateral vibrations. The vibrations of the cylinder were at all times synchronized with those in the wake, thereby suppressing the natural Strouhal frequency in favour of a common synchronized or ‘locked-in’ frequency for the body-wake system. All experiments were performed at a Reynolds number of 144 or 190. An inverse relation between the initial circulation K and the length lF of the vortex formation region was obtained for cylinder oscillations of up to 50% of a diameter, at vibration frequencies both above and below the Strouhal shedding frequency. The initial circulation K of the vortices was increased by as much as 65%, at lF = 1·6 diameters, from the stationary-cylinder value of K corresponding to lF = 3·2d. An increase in the rate A of vorticity generation of 80% from the stationary-cylinder wake value was obtained with the cylinder vibrating at 30% of a diameter and 110% of the Strouhal frequency. Both flow-visualization and hot-wire results show that the lateral spacing of the vortex street decreases as the vibration amplitude of the cylinder is increased, but that the longitudinal vortex spacing is independent of changes in amplitude. The longitudinal spacing, however, varies inversely with the vibration frequency. The street approaches a single line of vortices of alternating sign as the amplitude of vibration approaches values near a full cylinder diameter, and secondary vortex formation at these large amplitudes is associated with the vanishing lateral spacing of the street. Observation of the wake has elucidated the mechanism of vortex formation; the entrainment processes in the formation region have been observed at small intervals over a cycle of the cylinder's motion.
TL;DR: In this article, a two-dimensional stability analysis of the flow in a straight alluvial channel has been carried out, using the vorticity transport equation, and an attempt has been made to account for the influence of gravity on bed-load transport, and this turned out to change the stability quite significantly.
Abstract: A two-dimensional stability analysis of the flow in a straight alluvial channel has been carried out, using the vorticity transport equation. In the analysis an attempt has been made to account for the influence of gravity on bed-load transport, and this turned out to change the stability quite significantly.In the case of instability, the further growth of the dunes has been investigated using a second-order approximation, This nonlinear theory explains the experimental fact that the dunes very soon become asymmetric.
TL;DR: In this article, Batchelor's (1959) constant-strain dissipation spectrum and a joint probability distribution of the scalar field and its space derivatives were analyzed for convection of a sparse distribution of sheets of passive scalar in a random straining field whose correlation scale is large compared with the sheet size.
Abstract: The stretching of line elements, surface elements and wave vectors by a random, isotropic, solenoidal velocity field in D dimensions is studied. The rates of growth of line elements and (D – 1)-dimensional surface elements are found to be equal if the statistics are invariant to velocity reversal. The analysis is applied to convection of a sparse distribution of sheets of passive scalar in a random straining field whose correlation scale is large compared with the sheet size. This is Batchelor's (1959) κ−1 spectral regime. Some exact analytical solutions are found when the velocity field varies rapidly in time. These include the dissipation spectrum and a joint probability distribution that describes the simultaneous effect of Stretching and molecular diffusivity κ on the amplitude profile of a sheet. The latter leads to probability distributions of the scalar field and its space derivatives. For a growing κ−1 range at zero κ, these derivatives have essentially lognormal statistics. In the steady-state κ−1 regime at κ > 0, intermittencies measured by moment ratios are much smaller than for lognormal statistics, and they increase less rapidly with the order of the derivative than in the κ = 0 case. The κ > 0 distributions have singularities a t zero amplitude, due to a background of highly diffused sheets. The results do not depend strongly on D. But as D → ∞, temporal fluctuations in the stretching rates become negligible and Batchelor's (1959) constant-strain dissipation spectrum is recovered.
TL;DR: In this article, the averaged equations of slow flow in random arrays of fixed spheres are developed as a hierarchy of integro-differential equations, and an iteration procedure is described for obtaining the mean drag in the case of small volume concentration c. The results for the mean flow confirm Childress' terms in clogc and c, but indicate that for practical values of c numerical evaluation of integrals is needed, rather than expansion in powers of c and log c.
Abstract: The averaged equations of slow flow in random arrays of fixed spheres are developed as a hierarchy of integro-differential equations, and an iteration procedure is described for obtaining the mean drag in the case of small volume concentration c. The leading approximation is that given by Brinkman's model of flow past a single fixed sphere, in which the effects of all other spheres are treated as a Darcy resistance. The higher approximations take account of the modification to the mean flow, particularly in the near field, due to the localized nature of the actual resistance. Thus the second approximation finds the change due to a second sphere, and averages over all its possible positions. The result for the mean drag confirms Childress’ terms in clogc and c (apart from an arithmetical correction to the latter), but indicates that for practical values of c numerical evaluation of integrals is needed, rather than expansion in powers of c and log c. The last section of the paper develops the corresponding results for flow through random arrays of fixed parallel circular cylinders.
TL;DR: In this paper, the decay of two-dimensional, homogeneous, isotropic, incompressible turbulence is investigated both by means of numerical simulation (in spectral as well as in grid-point form), and theoretically by use of the direct-interaction approximation and the test-field model.
Abstract: The decay of two-dimensional, homogeneous, isotropic, incompressible turbulence is investigated both by means of numerical simulation (in spectral as well as in grid-point form), and theoretically by use of the direct-interaction approximation and the test-field model. The calculations cover the range of Reynolds numbers 50 ≤ RL ≤ 100. Comparison of spectral methods with finite-difference methods shows that one of the former with a given resolution is equivalent in accuracy to one of the latter with twice the resolution. The numerical simulations at the larger Reynolds numbers suggest that earlier reported simulations cannot be used in testing inertial-range theories. However, the large-scale features of the flow field appear to be remarkably independent of Reynolds number.The direct-interaction approximation is in satisfactory agreement with simulations in the energy-containing range, but grossly underestimates enstrophy transfer at high wavenumbers. The latter failing is traced to an inability to distinguish between convection and intrinsic distortion of small parcels of fluid. The test-field model on the other hand appears to be in excellent agreement with simulations at all wavenumbers, and for all Reynolds numbers investigated.
TL;DR: In this paper, the authors consider the motion of the mass of fluid ejected through a sharp-edged orifice by the motion a piston and show that the vorticity formed by viscous forces within the separated flow at the sharp edge rolls up to form a concentrated vortex which consists of a core of very fine scale turbulence surrounded by a co-moving bubble of much larger scale turbulence.
Abstract: We consider the motion of the mass of fluid ejected through a sharp-edged orifice by the motion of a piston. The vorticity formed by viscous forces within the separated flow at the sharp edge rolls up to form a concentrated vortex which, after a development period, consists of a core of very fine scale turbulence surrounded by a co-moving bubble of much larger scale turbulence. This bubble entrains outer fluid, mixes with it, and deposits the majority into a wake together with some small fraction of the total vorticity of the ring. Enough fluid is retained to account for the slow growth of the whole fluid mass. A theory which takes account of both the growth process and the loss of vorticity is proposed. By comparison with experimental measurements we have determined that the entrainment coefficient for turbulent vortex rings has a value equal to 0.011 ± 0.001, while their effective drag coefficient is 0.09 ± 0.01. These results seem to be independent of Reynolds number to within experimental accuracy.
TL;DR: In this paper, internal waves of the fundamental mode propagating into a shoaling region have been studied experimentally in a continuously stratified fluid and the waves divide into three classes depending upon the ratio of the bottom slope γ to the wave-characteristic slope c.
Abstract: Internal waves of the fundamental mode propagating into a shoaling region have been studied experimentally in a continuously stratified fluid. The waves divide into three classes depending upon the ratio of the bottom slope γ to the wave-characteristic slope c. For γ/c 1, the waves are inhomogeneous and have complex spatial dependence.