bluuebunny/crossref_metadata_2025 huggingface dataset:
Steady compressible swirl flows with closed streamlines at high Reynolds numbers
Authors: Guderley Karl Gottfried
Published: April 2006

In steady axisymmetric flows in a closed swirl chamber one can distinguish between the swirl flow proper, with components normal to the meridian plane, and a secondary flow whose components lie in the meridian plane. One can trace the motion of a particle within the meridian plane. The closed path so obtained will be called a streamline, to be parametrized by a stream function ϕ. One can distinguish between the flow in a boundary layer, where the velocity gradient is large, and a core flow, where the velocity and temperature gradients are relatively small. The present article is concerned with only the core flow. In high Reynolds number flows in which the streamlines are not closed there are three quantities which are constant along the streamlines: the total enthalpy (the right-hand side of Bernoulli's equation), the entropy and the moment of momentum of the particles with respect to the axis of symmetry. These are determined by conditions in the entrance cross-section. In flows with closed streamlines these quantities are ultimately determined by the cumulative effects of viscosity and heat conductivity. Conditions expressing cumulative effects enter the analysis as integrabiliby conditions necessary for the existence of a second approximation in a development of the flow field with respect to the reciprocal of the Reynolds number. They are re-expressed as the balance equations for energy, entropy and angular momentum which are to be satisfied on all surfaces ϕ = constant. One thus obtains an algorithm which leads from expressions for the total enthalpy, the entropy and the angular momentum as functions of ϕ to the residuals in the balance equations, also computed as functions of ϕ. The functions with which this algorithm starts must be chosen such that the residuals in the balance equations become zero. The secondary flow can be arbitrarily slow only if the Prandtl number is ½. At the centre of the secondary motion the balance equations are linearly dependent. This fact introduces an additional free parameter which allows one to compute secondary flows with different speeds. The linearized algorithm has the character of a Fredholm integral equation. This suggests an iterative solution similar to a Neumann series. The particles experience periodic changes of state which can be discussed as thermodynamic cycles. Such an analysis shows that the heat inputs occur on average at lower temperatures than the heat outputs. Besides the work that maintains the swirl motion, which is provided by shear-force components normal to the meridian plane, one therefore needs additional work, provided by the shear-force components within the meridian plane, to maintain a secondary motion. Responsible for this state of affairs is the fact that particles which do not quite maintain adiabatic temperatures move within a field with a large pressure gradient caused by the swirl component. This makes the flow sensitive to disturbances in the energy balance.
A numerical model of the air flow above water waves. Part 2
Authors: Gent P. R.
Published: April 2006

Further results from the nonlinear numerical model of the air flow in a deep turbulent boundary layer above water waves described in Gent & Taylor (1976) are presented. The results are calculated with the surface roughness z0 both constant and varying with position along the wave. With the form used when z0 varies, the fractional rate |ζ| of energy transfer per radian advance in phase due to the working of the pressure forces is larger than for z0 constant both when the transfer is from wind to waves and when it is from waves to wind. The latter case occurs when the waves are travelling faster than, or against, the wind. The energy transfer rates are compared with other theoretical predictions and with recent field observations.
Unsteady flows in a semi-infinite contracting or expanding pipe
Authors: Uchida Shigeo , Aoki Hiroshi
Published: April 2006

Physiological pumps produce flows by alternate contraction and expansion of the vessel. When muscles start to squeeze its wall the valve at the upstream end is closed and that at the downstream end is opened, and the fluid is pumped out in the downstream direction. These systems can be modelled by a semi-infinite pipe with one end closed by a compliant membrane which prevents only axial motion of the fluid, leaving radial motion completely unrestricted. In the present paper an exact similar solution of the Navier–Stokes equation for unsteady flow is a semi-infinite contracting or expanding circular pipe is calculated and reveals the following characteristics of this type of flow. In a contracting pipe the effects of viscosity are limited to a thin boundary layer attached to the wall, which becomes thinner for higher Reynolds numbers. In an expanding pipe the flow adjacent to the wall is highly retarded and eventually reverses at Reynolds numbers above a critical value. The pressure gradient along the axis of pipe is favourable for a contracting wall, while it is adverse for an expanding wall in most cases. These solutions are valid down to the state of a completely collapsed pipe, since the nonlinearity is retained in full. The results of the present theory may be applied to the unsteady flow produced by a certain class of forced contractions and expansions of a valved vein or a thin bronchial tube.
A note on the westward drift of the earth's magnetic field
Authors: Wood W. W.
Published: April 2006

A model of the earth's liquid core is assumed in which the underlying magnetic field and velocity are zonal and axially symmetric. Alfvén waves that vary as ei(kϕ−σt) are considered, where ϕ is the angle of longitude. Buoyancy and Coriolis forces Ω × U are included.

For a wide class of basic states and regions of flow, it is shown that roughly as many of the waves with a given k [ges ] 2 propagate eastwards as propagate westwards. All these waves are neutrally stable. The class of basic states is restricted by certain inequalities involving their velocity, magnetic field and entropy gradient.

it is observed that the known equivalence (Malkus 1967) between Alfvén waves with frequencies σ [Lt ] Ω and inertial waves with frequencies σp which are O(Ω) still holds when buoyancy forces are present. The equivalence requires σp2 to be real. If σp is pure imaginary, as is possible (though perhaps uncommon) in an unstably stratified medium, then the corresponding Alfvén wave is not neutrally stable and travels westwards.
Formation of the hexagonal pattern on the surface of a ferromagnetic fluid in an applied magnetic field
Authors: Gailītis A.
Published: April 2006

When a ferromagnetic fluid with a horizontal free surface is subjected to a uniform vertical applied magnetic field B0, it is known (Cowley & Rosensweig 1967) that the surface may be unstable when the field strength exceeds a certain critical value Bc. In this paper we consider, by means of an energy minimization principle, the possible forms that the surface may then take. Under the assumption that |μ − 1| [Lt ] 1 (where μ is the magnetic permeability of the fluid), it is shown that when B0 is near to Bc there are three equilibrium configurations for the surface: (i) flat surface, (ii) stationary hexagonal pattern, (iii) stationary square pattern. Configuration (i) is stable for B0 < Bc, (ii) is stable for B0 > Bc and B0−Bc sufficiently small, and (iii) is stable for some higher values of B0. In each configuration the fluid is static, and the surface is in equilibrium under the joint action of gravity, surface tension, and magnetic forces. The amplitude of the surface perturbation in cases (ii) and (iii) is calculated, and hysteresis effects associated with increase and decrease of B0 are discussed.
On steady flow in a partially filled rotating cylinder
Authors: Gans Roger F.
Published: April 2006

The flow in a partially filled cylinder rotating at right angles to the earth's gravity is found under the assumptions of rapid rotation and small viscosity. Effects of viscosity, nonlinear interaction and finite container length are included.
A method-of-characteristics calculation of coronary blood flow
Authors: Rumberger John A. , Nerem Robert M.
Published: April 2006

The one-dimensional unsteady flow equations for flow in an elastic tube have been solved by employing the method of characteristics and used to predict the development of the flow and pressure wave forms in the left coronary artery system of the horse. Input data to this model include in vivo measurements of wave speed and aortic-root pressure, both of which were carried out in horse experiments. In addition, estimates of vessel taper and fluid losses due to branching were established from plastic casts of the horse coronary arteries and from in vivo flowmeter data. The calculated results have confirmed earlier experiments in revealing a relatively large systolic component of flow in the major epicardial vessels. However, the calculations indicate that, once within the myocardium, the systolic flow component quickly diminishes and the diastolic flow component becomes increasingly important. The pressure pulse does not peak as observed in the aorta, but rather rounds out with systolic pressures decreasing slowly and diastolic pressures decreasing more rapidly with distance from the coronary ostium. The presence of relatively large amplitude, low frequency waves (of the order of 5–10 Hz), which were observed mainly during diastole in horse experiments, has also been confirmed by the computer calculations. Similar calculations carried out for conditions simulating the dog and human coronary systems indicate that such oscillations become higher in frequency and lower in amplitude with decreasing animal size.
On the growth, limiting thickness and dominant eddy scale of turbulent shearing layers in the atmosphere
Authors: Lessen Martin , Barcilon Albert , Butler Thomas E.
Published: April 2006

The growth of an unbounded, density-stratified, turbulent shearing layer in the presence of a gravity field is studied using the postulate of marginal instability. It is found, for a similarity mean velocity and density distribution, that after a rapid initial growth rate the growth slows asymptotically to zero as the Richardson number approaches a value of 1/4. Furthermore, the theory predicts a constant dominant turbulent eddy scale in all but the initial stages of growth of the turbulent shearing layer.

Both the general growth characteristics and the constant dominant turbulent eddy scale predicted by the theory are confirmed by experimental data.
Nonlinear stability of parallel flows with subcritical Reynolds numbers. Part 1. An asymptotic theory valid for small amplitude disturbances
Authors: Itoh Nobutake
Published: April 2006

A strict distinction is made between the two fundamental assumptions in the Stuart-Watson theory of nonlinear stability, one of which is that the amplitude of disturbance is sufficiently small, while the other is that the damping or amplification rate for an infinitesimal disturbance is small. This distinction leads to classification of the nonlinear stability theory into two asymptotic theories: the theory based on the first assumption can be applied to subcritical flows with Reynolds numbers away from the neutral curve, even to flows with no neutral curve, such as plane Couette flow or pipe Poiseuille flow, while the theory based on the second assumption is available only for Reynolds numbers and wavenumbers in the neighbourhood of the neutral curve. In the theory based on the first assumption the concept of trajectories in phase space, together with the method of eigenfunction expansion, is introduced in order to display nonlinear behaviour of the disturbance amplitude and to provide the most rational definition of the Landau constant available for classification of the behaviour patterns.
Nonlinear stability of parallel flows with subcritical Reynolds numbers. Part 2. Stability of pipe Poiseuille flow to finite axisymmetric disturbances
Authors: Itoh Nobutake
Published: April 2006

A theoretical study is presented of the spatial stability of flow in a circular pipe to small but finite axisymmetric disturbances. The disturbance is represented by a Fourier series with respect to time, and the truncated system of equations for the components up to the second-harmonic wave is derived under a rational assumption concerning the magnitudes of the Fourier components. The solution provides a relation between the damping rate and the amplitude of disturbance. Numerical calculations are carried out for Reynolds numbers R between 500 and 4000 and βR [les ] 5000, β being the non-dimensional frequency. The results indicate that the flow is stable to finite disturbances as well as to infinitesimal disturbances for all values of R and βR concerned.
Grid turbulence near a moving wall
Authors: Thomas N. H. , Hancock P. E.
Published: April 2006

Decaying grid turbulence was passed over a wall moving at the stream speed. For the high Reynolds number of the experiment, the field due to the wall constraint on the normal component of the velocity fluctuations is found to extend further into the flow than the influence of the viscous boundary condition on the tangential-component fluctuations. Measurements of the variances, length scales and spectra of the three velocity components of the turbulence are compared with the results of a previous experiment and with the theoretical predictions for an idealization of the flow. A simple model for some departures from the theory is proposed.
Stability of time-periodic flows in a circular pipe
Authors: Yang W. H. , Yih Chia-Shun
Published: April 2006

The stability of time-periodic flows in a circular pipe is investigated. The disturbance is assumed to be axially symmetric and to have a small amplitude, so that the governing differential equation is linear. Calculations are carried out for the first ten modes for a range of values of the frequency of the primary motion, of the wavenumber of the disturbance, and of the Reynolds number of the primary flow. In the ranges of the parameters for which the calculations have been carried out, the flows are found to be stable and, as for Stokes flows (von Kerczek & Davis 1974), it is conjectured that the flows under study here are stable for all frequencies and all Reynolds numbers.
The accelerated fully rough turbulent boundary layer
Authors: Coleman Hugh W. , Moffat Robert J. , Kays William M.
Published: April 2006

The behaviour of a fully rough turbulent boundary layer subjected to favourable pressure gradients both with and without blowing was investigated experimentally using a porous test surface composed of densely packed spheres of uniform size. Measurements of profiles of mean velocity and the components of the Reynolds-stress tensor are reported for both unblown and blown layers. Skin-friction coefficients were determined from measurements of the Reynolds shear stress and mean velocity.

An appropriate acceleration parameterKrfor fully rough layers is defined which is dependent on a characteristic roughness dimension but independent of molecular viscosity. For a constant blowing fractionFgreater than or equal to zero, the fully rough turbulent boundary layer reaches an equilibrium state whenKris held constant. Profiles of the mean velocity and the components of the Reynolds-stress tensor are then similar in the flow direction and the skin-friction coefficient, momentum thickness, boundary-layer shape factor and the Clauser shape factor and pressure-gradient parameter all become constant.

Acceleration of a fully rough layer decreases the normalized turbulent kinetic energy and makes the turbulence field much less isotropic in the inner region (forFequal to zero) compared with zero-pressure-gradient fully rough layers. The values of the Reynolds-shear-stress correlation coefficients, however, are unaffected by acceleration or blowing and are identical with values previously reported for smooth-wall and zero-pressure-gradient rough-wall flows. Increasing values of the roughness Reynolds number with acceleration indicate that the fully rough layer does not tend towards the transitionally rough or smooth-wall state when accelerated.
The slow motion of a rigid particle in a second-order fluid
Authors: Brunn P.
Published: April 2006

The purpose of the present paper is to reach some general conclusions on the motion of rigid particles in a homogeneous shear flow of a viscoelastic fluid. Under the basic assumption of nearly Newtonian slow flow, the creeping-motion equations for a second-order fluid with characteristic time constants κ0(2) and κ0(11) can be employed. It is shown that the κ0(2) contributions to the hydrodynamic force F and couple G depend upon the hydrodynamic force, couple and stresslet which act upon the particle in a Newtonian fluid (termed F(1), G(1) and S(1), respectively). Since this relation involves time derivatives of F(1) and G(1), a little reflexion is needed to realize that the modification of the classical Stokes law for steady translation in a quiescent fluid can have no κ0(2) term. Since no results of such generality are possible for the κ0(11) contributions we focus attention on transversely isotropic particles. Employing the concept of material tensors, the symmetry of such particles dictates the form these tensors adopt. This alone is sufficient to show that sedimentation in a quiescent fluid is accompanied by a change in orientation until a stable terminal orientation is attained. Depending upon the type of particle only one of the two orientations, axis of symmetry parallel or perpendicular to the external force, is stable. Another result concerns two-dimensional shear flow, for which we show that the symmetry axis has to drift through various Jeffery orbits until an equilibrium orientation is reached. While the orbits C = 0 and C = ∞ are equilibrium orbits for every transversely isotropic particle there may be a third such preferred orbit, which we denote by C*. In order for these orbits to be stable certain restrictions have to hold, showing that the orbits C = 0 and C* cannot both be stable. For the special case of a rigid tridumbbell of axis ratio s the orbit C* does not exist. If s > 1 the drift for this particle is into the orbit C = 0 while for s < 1 it is into the orbit C = ∞. This agrees qualitatively quite well with experimental results obtained for rods and disks. No quantitative comparison is possible; the particle shape influences the result quantitatively owing to its effect on the combination of the fluid parameters κ0(2) and κ0(11).
Convective instability of stably stratified water in the ocean
Authors: Leibovich S.
Published: April 2006

A recent theoretical description of interactions between surface waves and currents in the ocean is extended to allow density stratification. The interaction leads to a convective instability even when the density stratification is statically stable. An unspecified random surface wave field is permitted provided that it is statistically stationary.

The instability can be traced to torques produced by variations of a ‘vortex force’. Non-diffusive instabilities produced by this mechanism in water of infinite depth are explored in detail for arbitrary distributions of the destabilizing force. Stability is determined by an eigenvalue problem formally identical to that determining normal modes of infinitesimal internal waves in fluid with a density profile that is not monotone and thereby has a statically unstable region. Some tentative remarks are offered about the problem when dissipation is allowed.

Application of the present theory to Langmuir circulations is discussed. Also, according to the present theory, internal wave propagation should be modified by the vortex force arising from the interaction between the surface waves and the current.
An approximate theory for the streaming motion past axisymmetric bodies at Reynolds numbers of from 1 to approximately 100
Authors: Kolansky Michael S. , Weinbaum Sheldon , Pfeffer Robert
Published: April 2006

In Weinbaum et al. (1976) a simple new pressure hypothesis is derived which enables one to take account of the displacement interaction, the geometrical change in streamline radius of curvature and centrifugal effects in the thick viscous layers surrounding two-dimensional bluff bodies in the intermediate Reynolds number range O(1) < Re < O(102) using conventional Prandtl boundary-layer equations. The new pressure hypothesis states that the streamwise pressure gradient as a function of distance from the forward stagnation point on the displacement body is equal to the wall pressure gradient as a function of distance along the original body. This hypothesis is shown to be equivalent to stretching the streamwise body co-ordinate in conventional first-order boundary-layer theory. The present investigation shows that the same pressure hypothesis applies for the intermediate Reynolds number flow past axisymmetric bluff bodies except that the viscous term in the conventional axisymmetric boundary-layer equation must also be modified for transverse curvature effects O(δ) in the divergence of the stress tensor. The approximate solutions presented for the location of separation and the detailed surface pressure and vorticity distribution for the flow past spheres, spheroids and paraboloids of revolution at various Reynolds numbers in the range O(1) < Re < O(102) are in good agreement with available numerical Navier–Stokes solutions.
Internal gravity waves in a time-varying stratification
Authors: Rotunno Richard
Published: April 2006

The influence of slow time variations of the Brun-Väisälä frequencyNupon the energy of internal gravity waves is investigated. It is found that, when time variations inNare produced by a mean deformation field (reversible mean state), the wave energy can vary in either direct or inverse proportion, depending on the wavenumber orientation. WhenNchanges owing to a certain type of irreversible process, the wave energy varies with only inverse proportionality.

The nocturnal planetary boundary layer (NPBL) provides an example whereN=N(z, t). The full initial/boundary-value problem for anN(z, t) similar to the climatological mean for theNPBLis solved.
The excitation of damped waves diffracted over a submerged circular sill
Authors: Pite H. D.
Published: April 2006

A mathematical model of rectilinear surface waves incident on a submerged circular sill is developed. The class of waves considered consists of those for which the ratio of wavelength to water depth is large but which do not necessarily belong in the longwave category. A friction damping term is introduced into the equations of motion and the solutions obtained for the regions over the sill and in the ocean are matched by assuming a continuous surface and energy flux at the sill edge. The results show large reductions in the Q-factor of the resonance peaks brought about by friction damping. It is also found that, except at low frequency, a large number of overlapping resonance peaks which are out of phase with one another occupy a relatively narrow frequency band such that these resonance peaks effectively cancel one another. Experiments were performed to determine the friction constant used in the equations of motion and, using this friction constant, the theoretical results of wave resonance are verified.
The refraction of head seas by a long ship. Part 2. Waves of long wavelength
Authors: Ursell F.
Published: April 2006

It is known that head seas cannot travel without deformation along a cylinder of full constant cross-section, and recent calculations have indicated that the wave amplitude near the cylinder ultimately decreases as the waves travel along the cylinder, i.e. that the waves are refracted away from the axis of the cylinder. It was assumed in these calculations that the cross-section was a half-immersed circle of radius a of the same order as the wavelength 2π/K, but the method can probably be adapted to arbitrary full constant cross-sections. (There is however another calculation which indicates that for a thin ship the wave amplitude ultimately increases.) In the present paper these calculations are extended. The circular section is again studied but it is now supposed that the wavenumber Ka may be small. Uniformly valid expressions for the wave potential are obtained which show that for small Ka the refraction becomes significant only when Kx (the dimensionless distance along the cylinder) is so large that the product (Kx)½ v0(Ka) is also large; here the function v0(Ka) ∼ 2Ka arises in the solution of a certain eigenvalue problem. (The uniformly valid expressions also suggest an interpretation of the thin-ship calculation which resolves the apparent inconsistency.) The same method is applied to the waves generated by a pulsating source on an infinite cylinder, and similar results are obtained.
The motion of two-dimensional vortex pairs in a ground effect
Authors: Barker Steven J. , Crow Steven C.
Published: April 2006

A new technique for generating a pair of line vortices in the laboratory has been developed. The mean flow of these vortices is highly two-dimensional, although most of the flow field is turbulent. This two-dimensionality permits the study of vortex motions in the absence of the Crow mutual induction instability and other three-dimensional effects. The vortices are generated in a water tank of dimensions 15 × 122 × 244 cm. They propagate vertically and their axes span the 15 cm width of the tank. One wall of the tank is transparent, and the flow is visualized using fluorescein dye. High speed photography is used to study both the transition to turbulence during the vortex formation process and the interaction of the turbulent vortices with a simulated ground plane.

Transition occurs first in an annular region surrounding the core of each vortex, starting with a shear-layer instability on the rolled-up vortex sheet. The turbulent region then grows both radially inwards and radially outwards until the entire recirculation cell is turbulent. A ‘relaminarization’ of the vortex core appears to take place somewhat later.

The interaction of the vortex pair with the ground plane does not follow the predictions of potential-flow theory for line vortices. Although the total circulation is apparently conserved, the vortices remain at a larger distance from the ground than is expected and eventually ‘rebound’ or move away from the ground. Differences between a free-surface boundary condition and a smooth or rough ground plane are discussed. The ground-plane interaction is qualitatively very similar to that of aircraft trailing vortices observed in recent flight tests.
The scattering of short surface waves by a cylinder
Authors: Alker G.
Published: April 2006

The scattering of short surface waves by a partially immersed cylinder is considered. The cylinder is taken to be circular and to pass through two fixed points on the free surface with its centre on or above the free surface. Of particular interest is the behaviour of the solution as the cross-section of the immersed part of the cylinder approaches a semicircle. The method of matched asymptotic expansions is used.
Simple rotational flows
Authors: Shercliff J. A.
Published: April 2006

The paper develops and discusses some new additions to the available stock of analytical solutions of the nonlinear equations of fluid motion. The motions are steady, two-dimensional and devoid of viscous or other rotational forces (although such forces must have been significant during any starting process). The fluid density is constant.

The solutions are in two groups, referred respectively to Cartesian and polar co-ordinates. In both the stream function is of separable form, i.e. expressible as a product of two functions, each dependent on one co-ordinate. A remarkable variety of motions is revealed. Those that are most significant physically are described as bends (rapid transitions from one rectilinear flow to another) or as loops (closed, non-circular, vortex-type flows). The effects of boundary layers at walls or instability are not explored.

The paper closes with a mention of some preliminary experiments on loop flows in which all streamlines are ellipses and some discussion of the applicability of bend flows. Generalizations to axisymmetric flows and compressible flows are also mentioned briefly.
Experiments on the structure of turbulence in fully developed pipe flow. Part 2. A statistical procedure for identifying ‘bursts’ in the wall layers and some characteristics of flow during bursting periods
Authors: Heidrick T. R. , Banerjee S. , Azad R. S.
Published: April 2006

This paper is the second of a pair describing two-point velocity measurements in fully developed pipe flow. A method of processing hot-film anemometer signals to identify intervals of high energy production (‘bursts’) in wall turbulence is presented. The method uses filtered cross-stream spatial derivatives of the axial velocity fluctuations. It is demonstrated to be more sensitive to ‘bursts’ than several other methods of indentification. The bursts identified in this manner are shown to have similar characteristics to those observed in visual studies.

The technique has been applied to the wall region of turbulent pipe flow. Mean burst rates have been obtained at various distances from the wall for three Reynolds numbers. It is shown that the mean burst rate cannot be reliably obtained from a previously used technique based on the autocorrelation of the axial velocity fluctuations.

On the basis of our experiments, the mean burst rate and the turbulent shear stress have been found to vary similarly with distance from the wall. In the region near the wall where the shear stress is constant the mean burst rate is independent of the kinematic viscosity.

Some characteristics of the velocity fluctuations during burst intervals have been studied. All the bursts began with a relative minimum in the axial velocity fluctuations followed by a peak in the cross-stream spatial derivative. A second peak always occurred midway through the burst. The sequence of events is somewhat similar to that in the last stage of laminar-to-turbulent transition.
On the theory of solitary Rossby waves
Authors: Redekopp L. G.
Published: April 2006

The evolution of long, finite amplitude Rossby waves in a horizontally sheared zonal current is studied. The wave evolution is described by the Korteweg–de Vries equation or the modified Korteweg-de Vries equation depending on the atmospheric stratification. In either case, the cross-stream modal structure of these waves is given by the long-wave limit of the neutral eigensolutions of the barotropic stability equation. Both non-singular and singular eigensolutions are considered and the appropriate analysis is developed to yield a uniformly valid description of the motion in the critical-layer region where the wave speed matches the flow velocity. The analysis demonstrates that coherent, propagating, eddy structures can exist in stable shear flows and that these eddies have peculiar interaction properties quite distinct from the traditional views of turbulent motion.
Stochastic closure for nonlinear Rossby waves
Authors: Holloway Greg , Hendershott Myrl C.
Published: April 2006

An extension of the turbulence ‘test-field model’ (Kraichnan 1971a) is given for two-dimensional flow with Rossby-wave propagation. Such a unified treatment of waves and turbulence is necessary for flows in which the relative strength of nonlinear terms depends upon the length scale considered. We treat the geophysically interesting case in which long, fast Rossby waves propagate substantially without interaction while short Rossby waves are thoroughly dominated by advection. We recover the observations of Rhines (1975) that the tendency of two-dimensional flow to organize energy into larger scales of motion is inhibited by Rossby waves and that an initially isotropic flow develops anisotropy preferring zonal motion. The anisotropy evolves to an equilibrium functional dependence on the isotropic part of the flow spectrum. Theoretical results are found to be in quantitative agreement with numerical flow simulations.
Growth and equilibrium of short gravity waves in a wind-wave tank
Authors: Plant W. J. , Wright J. W.
Published: April 2006

Temporal and spatial development of short gravity waves in a linear wind-wave tank has been measured for wind speeds up to 15 m/s using microwave Doppler spectrometry. Surface waves of wavelength 4·1 cm, 9·8 cm, 16·5 cm and 36 cm were observed as a function of fetch, wind speed and wind duration. The waves grew exponentially from inception until they were about 10 dB smaller than their maximum height, and the temporal growth and spectral transport (spatial growth) rates were about equal when the wave amplitude was sufficiently small. The amplitude of a short gravity wave of fixed wavelength was found to decrease substantially at winds, fetches or durations greater than those at which the short gravity wave was approximately the dominant wave; such phenomena are sometimes referred to as overshoot. The dominant short gravity wave was observed to reach a maximum amplitude which depended only on wavelength, showing that wave breaking induced by an augmented wind drift cannot be the primary limitation to the wave height. Waves travelling against the wind were observed for wavelengths of 9·8 cm, 16·5 cm and 36 cm and were shown to be generated by the air flow at low wind speeds.

Measured initial growth rates for 16·5 cm and 36 cm waves were greater than expected, suggesting the existence of a growth mechanism in addition to direct transfer from the wind via linear instability of the boundary-layer flow. Initial temporal growth rates and spectral transport rates were compared to yield an experimental determination of the magnitude of the sum of nonlinear interactions and dissipation in short gravity waves. If the steady-state energy input in the neighbourhood of the dominant wave occurs at the measured initial temporal growth rates, then most of the energy input is locally dissipated; relatively little is advected away. Calculated gravitycapillary nonlinear energy transfer rates match those determined from initial growth rates for 9·8 cm waves and the gravity–capillary wave interaction continues to be significant for waves as long as 16·5 cm. For longer waves the gravity–capillary interaction is too small to bring the short gravity wave to a steady state when it is the dominant wave of the wind-wave system.
Three-dimensional natural convection in a box: a numerical study
Authors: Mallinson G. D. , Davis G. De Vahl
Published: April 2006

The solution of the steady-state Navier–Stokes equations in three dimensions has been obtained by a numerical method for the problem of natural convection in a rectangular cavity as a result of differential side heating. In the past, this problem has generally been treated as though it were two-dimensional. The solutions explore the three-dimensional motion generated by the presence of no-slip adiabatic end walls. For Ra = 104, the three-dimensional motion is shown to be the result of the inertial interaction of the rotating flow with the stationary walls together with a contribution arising from buoyancy forces generated by longitudinal temperature gradients. The inertial effect is inversely dependent on the Prandtl number, whereas the thermal effect is nearly constant. For higher values of Ra, multiple longitudinal flows develop which are a delicate function of Ra, Pr and the cavity aspect ratios.
Laminar boundary-layer flow past a two-dimensional slot
Authors: Sobey Ian J.
Published: April 2006

Observations by Cornhill ' Roach (1976) of sudanophilic lesions in the vicinity of intercostal arteries in rabbit aortas have shown that lesions develop on the downstream side of the associated ostia. There is considerable conjecture as to the role which varying levels of wall shear stress play in the development of such lesions; Cornhill ' Roach implicate high wall shear stress levels. We develop a consistent model of steady boundary-layer flow past a side slot assuming that there is Stokes flow in the side slot and that the main body of the boundary layer remains undisturbed. Our results show that increased levels of wall shear stress occur both upstream and downstream of the slot. If the withdrawal of fluid through the side slot is sufficiently great there may be a stagnation point on the downstream side of the slot. The wall shear stress level near the slot depends on both normal and transverse motions at the mouth of the slot. Indeed, very near the slot, on a length scale comparable with the slot width, the wall shear stress level depends only on the transverse motions at the mouth of the slot.
Nonlinear deep-water waves: theory and experiment. Part 2. Evolution of a continuous wave train
Authors: Lake Bruce M. , Yuen Henry C. , Rungaldier Harald , Ferguson Warren E.
Published: April 2006

Results of an experimental investigation of the evolution of a nonlinear wave train on deep water are reported. The initial stage of evolution is found to be characterized by exponential growth of a modulational instability, as was first discovered by Benjamin ' Feir. At later stages of evolution it is found that the instability does not lead to wave-train disintegration or loss of coherence. Instead, the modulation periodically increases and decreases, and the wave train exhibits the Fermi–Pasta–Ulam recurrence phenomenon. Results of an earlier study of nonlinear wave packets by Yuen ' Lake, in which solutions of the nonlinear Schrödinger equation were shown to provide quantitatively correct descriptions of the properties of nonlinear wave packets, are applied to describe the experimentally observed wave-train phenomena. A comparison between the laboratory data and numerical solutions of the nonlinear Schrödinger equation for the long-time evolution of nonlinear wave trains is given.
A note on some nonlinear water-wave experiments and the comparison of data with theory
Authors: Lake Bruce M. , Yuen Henry C.
Published: April 2006

The problem of the instability of a uniform, nonlinear, deep-water wave train to infinitesimal long-wave perturbations, first studied by Benjamin ' Feir (1967) and Benjamin (1967), is re-examined. It is found that the apparent discrepancy between the experimental and theoretical growth rates of the instability is associated with the experimental generation of waves which do not have the Stokes wave profiles assumed in the theory. Experimental and theoretical results relating the initial wave steepness and the most unstable long-wave perturbation are used to obtain a correction factor, which is found to account for the mismatch in wave forms and which resolves the discrepancy in growth rates. The results illustrate that, when theory is compared with experiments in which the values of certain higher-order (nonlinear) quantities must be deduced from measurements of first-order quantities, great care must be taken to ascertain that the experimental conditions and the theoretical assumptions are indeed compatible to the required order.
Stability of time-dependent double-diffusive convection in an inclined slot
Authors: Chen C. F. , Sandford R. D.
Published: April 2006

The two-dimensional motion of a stably stratified fluid containing two solutes with different molecular diffusivities in an inclined slot has recently been examined by Chen (1975, hereafter referred to as I). The two solutes have continuous opposing gradients with the slower-diffusing one more dense at the bottom. It is found that, in the steady state, there exists a slow upward flow along the slope driven by the slight buoyancy difference near the wall, not unlike the solution found by Wunsch (1970) and Phillips (1970) for a single solute. For the time-dependent flow resulting from switching on the diffusivities at t = 0, there may be a flow reversal near the wall depending on the relative magnitude of λ and τ (where λ is the ratio of the density gradient and τ−1 is the ratio of the diffusivity of the faster-diffusing solute T to that of the slower-diffusing one S). By examining the distributions of S and T across the slot, it becomes apparent that in cases with flow reversal double-diffusive instability is likely to occur.

In this paper, we examine the stability of time-dependent double-diffusive convection in an inclined slot both analytically and experimentally. The time-dependent perturbation equations are numerically integrated starting with an initial distribution of small random disturbances in the vorticity. The growth or decay of the kinetic energy of the perturbations serves to indicate whether the flow is unstable or stable. The results show that the flow becomes more unstable (a) with increasing λ at a given angle of inclination with respect to the vertical and (b) with increasing angle of inclination at a given value of λ. Experiments were carried out in a 2[sdot ]54 cm wide slot using sugar and salt solutions at angles of inclination of 30°, 45° and 60°. Results obtained confirm the trends predicted by the analysis. Good agreement was obtained between the predicted and the experimental values of the critical wavelength for the case λ = 0[sdot ]7.
The effect of Brownian motion on the bulk stress in a suspension of spherical particles
Authors: Batchelor G. K.
Published: April 2006

The effect of Brownian motion of particles in a statistically homogeneous suspension is to tend to make uniform the joint probability density functions for the relative positions of particles, in opposition to the tendency of a deforming motion of the suspension to make some particle configurations more common. This smoothing process of Brownian motion can be represented by the action of coupled or interactive steady ‘thermodynamic’ forces on the particles, which have two effects relevant to the bulk stress in the suspension. Firstly, the system of thermodynamic forces on particles makes a direct contribution to the bulk stress; and, secondly, thermodynamic forces change the statistical properties of the relative positions of particles and so affect the bulk stress indirectly. These two effects are analysed for a suspension of rigid spherical particles. In the case of a dilute suspension both the direct and indirect contributions to the bulk stress due to Brownian motion are of order ø2, where ø([Lt ] 1) is the volume fraction of the particles, and an explicit expression for this leading approximation is constructed in terms of hydrodynamic interactions between pairs of particles. The differential equation representing the effects of the bulk deforming motion and the Brownian motion on the probability density of the separation vector of particle pairs in a dilute suspension is also investigated, and is solved numerically for the case of relatively strong Brownian motion. The suspension has approximately isotropic structure in this case, regardless of the nature of the bulk flow, and the effective viscosity representing the stress system to order ϕ2 is found to be
\[
\mu^{*} = \mu(1+2.5\phi + 6.2\phi^2).
\]
The value of the coefficient of ø2 for steady pure straining motion in the case of weak Brownian motion is known to be 7[sdot ]6, which indicates a small degree of ‘strain thickening’ in the ø2-term.
Waves on water jets
Authors: Hoyt J. W. , Taylor J. J.
Published: April 2006

By the use of high-speed photography, instabilities occurring in high Reynolds number water jets discharging into air have been made visible. These instabilities include the axisymmetric mode accompanying the transition from laminar to turbulent flow at the nozzle exit, spray formation as a culmination of the axisymmetric disturbances, and, further downstream, helical disturbances which result in the entire jet assuming a helical form. The final disruption of the jet is due to amplification of the helical waves. It is further shown that the amplification of the helical disturbances is due in part to aerodynamic form drag, since jets discharging into surrounding air moving at the same speed as the jet remain relatively stable, compared with the case when the jet is discharged into stagnant air.
The diffusion of scalar and vector fields by homogeneous stationary turbulence
Authors: Knobloch Edgar
Published: April 2006

An exact Eulerian formulation of the problem of diffusion of passive scalar and vector fields by a turbulent velocity field is obtained. It is shown that, in the short autocorrelation time limit, the diffusion equation is exact for any turbulence. For non-zero autocorrelation times the form of the first few correction terms to the diffusion equation is found. As a result of these corrections the diffusion of scalar, divergence-free and curl-free vector fields will be different. The calculations use the Kubo–Van Kampen–Terwiel technique and are carried out for zero ordinary diffusivity and for homogeneous, stationary, isotropic, incompressible, helical turbulence.
Long surface waves incident on a submerged horizontal plate
Authors: Siew P. F. , Hurley D. G.
Published: April 2006

A train of surface gravity waves of wavelength λ in a channel of depth H is incident on a horizontal plate of length l that is submerged to a depth c. Under the assumption that both λ and l are large compared with H, the method of matched asymptotic expansions is used to show that, to first order, the reflexion coefficient R and the transmission coefficient T are given by
\[
R = \chi \left\{\frac{\sigma l}{(gH)^{\frac{1}{2}}}\sin\frac{\sigma l}{(gc)^{\frac{1}{2}}}-2\bigg(\frac{c}{H}\bigg)^{\frac{1}{2}}\bigg(1-\cos\frac{\sigma l}{(gc)^{\frac{1}{2}}}\bigg)\right\}
\]
and
\[
T =\chi\left\{2i\left[\sin\frac{\sigma l}{(gc)^{\frac{1}{2}}}+\frac{\sigma l}{b}\bigg(\frac{c}{g}\bigg)^{\frac{1}{2}}\right]\right\}
\]
where
\begin{eqnarray*}
\chi &=& 1\left/ \left\{2\bigg(\frac{c}{H}\bigg)^{\frac{1}{2}}\bigg(1-\cos\frac{\sigma l}{(gc)^{\frac{1}{2}}}\bigg)+\frac{\sigma l}{b}\bigg(\frac{H}{g}\bigg)^{\frac{1}{2}}\bigg(1+\frac{c}{H}\bigg)\sin\frac{\sigma l}{(gc)^{\frac{1}{2}}}\right.\right.\\
&&\left. +2i\bigg(\sin\frac{\sigma l}{(gc)^{\frac{1}{2}}}+\frac{\sigma l}{b}\bigg(\frac{c}{g}\bigg)^{\frac{1}{2}}\cos\frac{\sigma l}{(gc)^{\frac{1}{2}}}\bigg)\right\},
\end{eqnarray*}
σ is the angular frequency and g the acceleration due to gravity.
On Hamilton's principle for surface waves
Authors: Miles John W.
Published: April 2006

The boundary-value problem for irrotational surface waves is derived from a variational integral I with the Lagrangian density [Lscr ] = Ξ ηt - [Hscr ] where Ξ (X, t) is the value of the velocity potential at the free surface, y = η(x, t), and [Hscr ] is the energy density in x space. [Hscr ] then is expressed as a functional of Ξ and η, qua canonical variables, by solving a reduced boundary-value problem for the potential, after which the requirement that I be stationary with respect to independent variations of Ξ and η yields a pair of evolution equations for Ξ and η. The Fourier expansions Ξ = pn(t)ϕn*(x) and η = qn(t)ϕn(x), where {ϕn} is an orthogonal set of basis functions, reduce I to Hamilton's action integral, in which the complex amplitudes pn and qn appear as canonically conjugate co-ordinates, and yield canonical equations for pn and qn that are the spectral transforms of the evolution equations for Ξ and η. The evolution equations are reduced (asymptotically) to partial differential equations for pn and qn by expanding [Hscr ] in powers of α = a/d and β = (d/l)2, where a and l are scales of amplitude and wavelength. Explicit third approximations are developed for β = O(α).
A note regarding ‘On Hamilton's principle for surface waves’
Authors: Milder D. Michael
Published: April 2006

The canonical form of the equations for the free-surface elevation and potential of an irrotational fluid is more than a coincidence. The elevation is a ‘generalized coordinate’ field sufficient to define the system Lagrangian without explicit reference to the motion of the fluid interior. The Lagrangian and the associated field equations are complete and self-contained in the two-dimensional surface co-ordinates, but non-local (integro-differential) in form; the canonical equations derived by Miles are just the Hamiltonian counterparts.