bluuebunny/crossref_metadata_2025_split huggingface dataset:
An investigation of the stability of a magnetohydrodynamic shear layer
Authors: Malcolm D. G. 
Published: March 2006

This paper presents an experimental investigation of the conditions of instability and post-critical flow in an electrically driven shear layer which has been previously studied both theoretically and experimentally. The shear layer forms as a primarily azimuthal flow in an electrically conducting fluid between the edges of two parallel and coaxial circular electrodes, mounted in insulating planes, when a current is passed between them in the presence of a strong, axially applied magnetic field. The aximuthal velocity at any radial position in the layer is maximum near the electrode planes and decreases to zero at the centreplane between the electrodes. Radial and axial secondary velocity components apparently exist, though smaller in magnitude than the azimuthal velocity throughout most of the layer.The experiments employ a miniature hot-film sensor as an instantaneous velocity indicator in mercury. The following experimental results are reported.
\begin{eqnarray*}
& R_c = 2420,\quad 150 \leqslant M \leqslant 350,\\
& R_c = 345\,M^{\frac{1}{3}},\quad 350\leqslant M \leqslant 650,
\end{eqnarray*}
where Rc is the critical Reynolds number of the shear layer, based upon the maximum velocity in the layer and the layer thickness, and M is the Hartmann number based upon the distance from the electrodes to the centreplane. For values of M below about 300, oscillograms of the hot-film signal show that the layer goes unstable in a chaotic fashion. For M ≥ 300, however, the instabilities grow to produce a wave-like post-critical flow of steady amplitude. The waves are apparently parallel to the applied magnetic field and drift in the azimuthal direction at a constant proportion (10–15%) of the critical velocity. The maximum intensity of the peak to peak fluctuations from the measured mean azimuthal velocity appears to be about 0·05. These fluctuations decrease to zero at the centreplane. By an approximate measurement of wavelength it is determined that the integral number of waves in the post-critical flow increases with M.The physical mechanism governing the onset of instability has not been established. The manner in which hydrodynamic and electromagnetic forces interact with the destabilizing disturbances cannot be fully understood without further theoretical analysis.
The stress system in a suspension of force-free particles
Authors: Batchelor G. K. 
Published: March 2006

The purpose of the paper is to consider in general terms the properties of the bulk stress in a suspension of non-spherical particles, on which a couple (but no force) may be imposed by external means, immersed in a Newtonian fluid. The stress is sought in terms of the instantaneous particle orientations, and the problem of determining these orientations from the history of the motion is not considered. The bulk stress and bulk velocity gradient in the suspension are defined as averages over an ensemble of realizations, these averages being equal to integrals over a suitably chosen volume of ambient fluid and particles together when the suspension is statistically homogeneous. Without restriction on the type of particle or the concentration or the Reynolds number of the motion, the contribution to the bulk stress due to the presence of the particles is expressed in terms of integrals involving the stress and velocity over the surfaces of particles together with volume integrals not involving the stress. The antisymmetric part of this bulk stress is equal to half the total couple imposed on the particles per unit volume of the suspension. When the Reynolds number of the relative motion near one particle is small, a suspension of couple-free particles of constant shape is quasi-Newtonian; i.e. the dependence of the bulk stress on bulk velocity gradient is linear. Two significant features of a suspension of non-spherical particles are (1) that this linear relation is not of the Newtonian form and (2) that the effect of exerting a couple on the particles is not confined to the generation of an antisymmetrical part of the bulk stress tensor. The role of surface tension at the particle boundaries is described.In the case of a dilute suspension the contributions to the bulk stress from the various particles are independent, and the contributions arising from the bulk rate of strain and from the imposed couple are independent for each particle. Each particle acts effectively as a force doublet (i.e. equal and opposite adjoining ‘Stokeslets’) whose tensor strength determines the disturbance flow far from the particle and whose symmetrical and antisymmetrical parts are designated as a stresslet and a couplet. The couplet strength is determined wholly by the externally imposed couple on the particle; but the stresslet strength depends both on the bulk rate of strain and, for a non-spherical particle, on the rate of rotation of the particle relative to the fluid resulting from the imposed couple. The general properties of the stress system in a dilute suspension are illustrated by the specific and complete results which may be obtained for rigid ellipsoidal particles by use of the work by Jeffery (1922).
The Cameron—Martin—Wiener method in turbulence and in Burgers’ model: general formulae, and application to late decay
Authors: Kahng W-H. , Siegel A. 
Published: March 2006

We apply the Cameron—Martin—Wiener (formerly ‘Wiener—Hermite’) expansion of a random velocity field to the analytical study of turbulence. The kernels of this expansion contain all statistical information about the ensemble. Complete expressions are derived for constructing statistical quantities in terms of the kernels, and for the equations of motion of the kernels. We rigorously prove the Gaussian trend of the velocity field of the Navier—Stokes equation in the very late stage when the non-linear term is neglected. Then-dependence (nis the order of derivative) of the flatness factor, minus three for derivatives of the velocity field, shows a rapid increase withnin this stage.The late decay problem of the Burgers model of turbulence is studied analytically with a view to obtaining suggestive guidelines for fitting the non-linear aspects of the model turbulence. We can divide the energy spectrum density into two parts, the larger of which is a kind of steady solution, which we call the ‘equilibrium state’, which remains self-similar in time in terms of an appropriate variable. The deviation from this ‘equilibrium solution’ satisfies the Kármán—Howarth equation. As initial velocity field, we take two particular cases: (a) a pure Gaussian, and (b) a non-Gaussian velocity field. With these two cases a detailed spectral analysis has been obtained. The energy spectrum deviation from ‘equilibrium’ declines exponentially to zero for all wave-numbers. The Gaussian case shows that the flatness factor minus three increases rapidly withn, while the non-Gaussian case does not show any marked dependence onn.
Gravity wave refraction by islands
Authors: Lautenbacher Conrad C. 
Published: March 2006

The refractive influence on tsunami run-up of the offshore bottom topography of islands is analyzed. Shallow water theory is used to treat problems in which the geometry resembles that of individual Hawaiian islands and in which the incident wave is plane and monochromatic.Mathematically, the differential equation for long-wave propagation is converted into an integral equation to which numerical methods are applied. Results of practical importance include the run-up on island coastal areas. The results are used in conjunction with earlier one-dimensional analyses to estimate the total tsunami run-up.
The turbulent boundary layer on a rotating cylinder in an axial stream
Authors: Cham T-S. , Head M. R. 
Published: March 2006

Calculations are presented of the development of the turbulent boundary layer on a rotating cylinder with its axis along the stream, and the results are compared with measurements made by Parr. With the choice of a suitable co-ordinate system, the boundary layer downstream of the nose of the cylinder approximates closely to a condition of two-dimensionality, and a simple integral method of solution can be applied. The only evidence of three-dimensionality lies in the destabilizing effect of rotation on the turbulence structure of the layer and an analysis of this instability has been made which relates changes in mixing length and entrainment to an instability parameter in the form of a Richardson number.In spite of the changes in shear stress distribution and entrainment brought about by rotation, mean velocity profiles and skin friction values are found to be related to H and Rθ in the same way as for the stationary flat plate.
A theory of turbulence in stratified fluids
Authors: Long Robert R. 
Published: March 2006

An effort is made to understand turbulence in fluid systems like the oceans and atmosphere in which the Richardson number is generally large. Toward this end, a theory is developed for turbulent flow over a flat plate which is moved and cooled in such a way as to produce constant vertical fluxes of momentum and heat. The theory indicates that in a co-ordinate system fixed in the plate the mean velocity increases linearly with heightzabove a turbulent boundary layer and the mean density decreases asz3, so that the Richardson number is large far from the plate. Near the plate, the results reduce to those of Monin & Obukhov.Thecurvatureof the density profile is essential in the formulation of the theory. When the curvature is negative, a volume of fluid, thoroughly mixed by turbulence, will tend to flatten out at a new level well above the original centre of mass, thereby transporting heat downward. When the curvature is positive a mixed volume of fluid will tend to fall a similar distance, again transporting heat downward. A well-mixed volume of fluid will also tend to rise when the density profile is linear, but this rise is negligible on the basis of the Boussinesq approximation. The interchange of fluid of different, mean horizontal speeds in the formation of the turbulent patch transfers momentum. As the mixing in the patch destroys the mean velocity shear locally, kinetic energy is transferred from mean motion to disturbed motion. The turbulence can arise in spite of the high Richardson number because the precise variations of mean density and mean velocity mentioned above permit wave energy to propagate from the turbulent boundary layer to the whole region above the plate. At the levels of reflexion, where the amplitudes become large, wave-breaking and turbulence will tend to develop.The relationship between the curvature of the density profile and the transfer of heat suggests that the density gradient near the level of a point of inflexion of the density curve (in general cases of stratified, shearing flow) will increase locally as time goes on. There will also be a tendency to increase the shear through the action of local wave stresses. If this results in a progressive reduction in Richardson number, an ultimate outbreak of Kelvin–Helmholtz instability will occur. The resulting sporadic turbulence will transfer heat (and momentum) through the level of the inflexion point. This mechanism for the appearance of regions of low Richardson number is offered as a possible explanation for the formation of the surfaces of strong density and velocity differences observed in the oceans and atmosphere, and for the turbulence that appears on these surfaces.
On stagnation-point conditions in non-equilibrium inviscid blunt-body flows
Authors: Vinokur Marcel 
Published: March 2006

In non-equilibrium inviscid blunt-body flows, the state of the gas at the stagnation point is known to be in thermodynamic equilibrium for all finite relaxation times. The dependence of this state on the non-equilibrium processes and body geometry is investigated for the most general conditions. The stagnation-point state is always found to be in a narrow range bounded on one side by the state obtained in an equilibrium flow. The other bound, called the frozen limit, is far removed from the state obtained in an identically frozen flow (infinite relaxation times). For certain state variables, the frozen-limit value lies outside the range determined by frozen and equilibrium flow. Significant errors are found in several published predictions of the stagnation-point state, resulting from the non-analytic approach to equilibrium in nearly frozen flow.The two bounds on the pressure are expressible in terms of the normal shock density ratios for equilibrium and frozen flow. The actual pressure for an arbitrary flow situation is found to depend only on the shock nose radius and the relation between density and time in the relaxation zone behind a normal shock wave. If the density law is given by a single relaxation model, a closed form expression for the pressure is obtained. The analysis is carried out for both plane and axisymmetric flows, and is also valid for non-equilibrium free-stream conditions.
Thermal boundary-layer theory near the stagnation point in three-dimensional fluctuating flow
Authors: Ghoshal S. , Ghoshal A. 
Published: March 2006

The equations of motion and energy governing a three-dimensional fluctuating flow of an incompressible fluid in the vicinity of a stagnation point on a regular surface have been integrated analytically. The velocity of the oncoming flow relative to the body oscillates in magnitude but not in direction.It has also been shown that the analysis of Lighthill for the two-dimensional fluctuating flow may be extended to the three-dimensional flow (both chordwise and spanwise), namely for each point on the body there is a critical frequency ω0 such that for frequencies ω > ω0 the oscillations are to a close approximation ordinary ‘shear waves’, unaffected by the mean flow; the phase advance in the skin friction is then 45°. For frequencies ω < ω0 the oscillations may be closely approximated by the sum of two parts: one quasi-steady part and the other proportional to the acceleration of the oncoming stream. The phase advance in the skin friction is then tan−1 (ω/ω0).
An optical study of turbulence
Authors: Roe G. E. 
Published: March 2006

This paper describes theoretical and experimental work carried out at the Cavendish Laboratory of the University of Cambridge. The main object of the work was to develop a new technique for measuring the structure of fluid turbulence.A parallel beam of light is passed through the turbulent region, containing refractive index fluctuations, and analyzed on exit by gratings of periodic transmissivity. Two forms of analysis yield (a) the spatial power spectrum of the refractive index fluctuations in the turbulence, and (b) the velocity distribution within the beam aperture. The method does not disturb the fluid physically, does not depend on the existence of a mean flow velocity, and works well in liquids.One of the limitations of this single-beam method is that it produces information averaged along the path length of the beam in the turbulence, and to overcome this a cross-beam technique, using two beams intersecting at right-angles, has been developed in theory. This method gives the spatial power spectrum of the refractive index fluctuations, as does the single beam method, but the results are characteristic only of the volume of intersection of the beams.The paper first discusses the theory of the single-beam and crossed-beam techniques, and then experimental results obtained with the single-beam method.The turbulent region investigated was a rectangular tank of water, heated from below and cooled from above, producing convective turbulence of high Rayleigh number (4·1 × 108), a system difficult to analyze by conventional methods of measurement, such as the hot-wire anemometer.Spectral density functions (power spectra) of refractive index, and hence in this case temperature fluctuations, have been measured, as have velocity distributions. Statistical analysis of the results also gives useful information about the Eulerian time scale of the turbulent field.
A study of free jet impingement. Part 2. Free jet turbulent structure and impingement heat transfer
Authors: Donaldson Coleman Dup. , Snedeker Richard S. , Margolis David P. 
Published: March 2006

An experimental study of jet impingement is completed with the presentation of the measured turbulent characteristics of the circular subsonic jet and the heat transfer rates measured when this jet impinges normal to a flat plate. The data suggest that for impingement very close to the stagnation point, the heat transfer can be computed by applying a turbulent correction factor to the laminar value calculated for a flow having the same pressure distribution as that present in the impingement region. The correction factor is found to be a function of the axial distance and not of Reynolds number. Farther away, the measurements agree well with the heat transfer estimated using the method of Rosenbaum & Donaldson (1967). At large distances from the stagnation point, the heat transfer falls off in inverse proportion with the distance.The documentation of the turbulent jet flow field includes measurements of the radial and axial velocity fluctuations and their spectra, as well as the radial distribution of turbulent shear$\overline{w^{\prime}u^{\prime}}$. In addition, measurements of the turbulence near the stagnation point and the total pressure fluctuation at the stagnation point are presented.
Finite-amplitude stability of pipe flow
Authors: Davey A. , Nguyen H. P. F. 
Published: March 2006

In this paper we present some results concerning the stability of flow in a circular pipe to small but finite axisymmetric disturbances. The flow is unstable if the amplitude of a disturbance exceeds a critical value, the equilibrium amplitude, which we have calculated for a wide range of wave-numbers and Reynolds numbers. For large values of the Reynolds number, R, and for a real value of the wave-number, α, we indicate that the energy density of a critical disturbance is of order c2i, where −ααci is the damping rate of the associated infinitesimal disturbance. The energy, per unit length of the pipe, of a critical disturbance which is concentrated near the axis of the pipe is of order R−2, and the wave-number α is of order R1/3 For a critical disturbance which is concentrated near the wall of the pipe the energy is of order $R^{-\frac{3}{2}}$ and α is of order R½. This suggests that non-linear instability is most likely to be caused by a ‘centre’ mode rather than by a ‘wall’ mode. The wall mode solution is also essentially the solution for the problem of plane Couette flow when αR is large. We compare it with the true solution.In an appendix Dr A. E. Gill indicates how some of the results of this paper may be inferred from a simple scale analysis.
Convection in a rotating annulus uniformly heated from below
Authors: Davies-Jones Robert P. , Gilman Peter A. 
Published: March 2006

We present a linear stability analysis, to second order in initial amplitude, of Bénard convection of a Boussinesq fluid in a thin rotating annulus for modest Taylor numbers T ([les ] 104). The work is motivated in part by the desire to study further a mechanism for maintaining, through horizontal Reynolds stresses induced in the convection, the sun's ‘equatorial acceleration’, which has been demonstrated for a rotating convecting spherical shell by Busse & Durney. The annulus is assumed to have stress free, perfectly conducting top and bottom (which allows separation of the equations) and non-conducting non-slip sides. A laboratory experiment which fulfills these conditions (except perhaps the free bottom) is being developed with H. Snyder.We study primarily annuli with gap-width to depth ratios a of order unity. The close, non-slip side-walls produce a number of effects not present in the infinite plane case, including overstability at high Prandtl numbers P, and multiple minima in Rayleigh number R on the stability boundary. The latter may give rise to vacillation. The eigenfunctions for stationary convection for a = 2, T [lsim ] 2000 clearly show momentum of the same sense as the rotation is transported from the inner to the outer half of the annulus, corresponding to transport toward equatorial latitudes on the sphere. The complete second-order solutions for the induced circulations indeed give faster rotation in the outer half, except for large P (> 102), in which case thermal stresses dominate. At all P, this differential rotation is qualitatively a thermal wind. Overstable convective cells, and stationary cells at higher T, induce more complicated differential rotations.
Calculation of boundary-layer development using the turbulent energy equation: compressible flow on adiabatic walls
Authors: Bradshaw P. , Ferriss D. H. 
Published: March 2006

The basic method described by Bradshaw, Ferriss & Atwell (1967) is extended to compressible flow in two-dimensional boundary layers in arbitrary pressure gradient (excluding shock waves and expansion fans) by invoking Morkovin's hypothesis (Favre 1964) that the turbulence structure is unaffected by compressibility. Using the same empirical functions as in incompressible flow, skin friction in zero pressure gradient is predicted to within 3% of Spalding & Chi's (1964) correlation for free-stream Mach numbers less than 5. Comparisons with experiments in pressure gradient are restricted by the lack of data, but, since Morkovin's hypothesis does not depend on pressure gradient, methods which use it (of which the present method seems to be the first) can be checked fairly adequately by comparisons with data in zero pressure gradient.No ‘compressibility transformations’ are needed, although the Crocco relation is used, provisionally, for the temperature: since the calculations take only about 20% longer than in incompressible flow, Morkovin's hypothesis does as much as any transformation could do. It is pointed out that, in supersonic flow, surface curvature which is large enough to induce a significant longitudinal pressure gradient is also large enough to have a very significant effect on the turbulence structure.
Wave forces on a circular dock
Authors: Garrett C. J. R. 
Published: March 2006

The scattering of surface gravity waves by a circular dock is considered in order to determine the horizontal and vertical forces and torque on the dock. An incident plane wave is expanded in Bessel functions, and for each mode the problem is formulated in terms of the potential on the cylindrical surface containing the dock and extending to the bottom. The solution is shown to have phase independent of depth and so may be obtained from an infinite set of real equations, which are solved numerically by Galerkin's method. The convergence of the solution is discussed, and some numerical results are presented.This problem has been investigated previously by Miles & Gilbert (1968) by a different method, but their work contained errors.