ALBERT

Beschreibung: This note provides some explanation of the fact that, contrary to the requirements of local isotropy, the skewness S of the streamwise temperature derivative ∂θ/∂x1 has been observed to be a non-zero constant of magnitude of about unity in high-Reynolds-number and high-Péclet-number turbulent shear flows. Measurements in slightly heated homogeneous shear flows and in unsheared grid turbulence suggest that S is non-zero only when the mean shear dU1/dx2 and the mean temperature gradient dT/dx2 are both non-zero. The sign of S is given by –sgn (dU1/dx2).sgn (dT/dx2). For fixed dU1/dx2, S is of the form tanh (αdT/dx2), α being a constant, while for fixed dT/dx2, it is of the form S/S* = 1 − β1 exp (− β2τ), where S* is a characteristic value of S, β1 and β2 are positive constants, and τ can be interpreted as a ‘total strain’. The derivative skewness data in other (inhomogeneous) shear flows are also compatible with the latter relation. Predictions from a simplified transport equation for [formula omitted], derived in the light of the present experimental observations, are in reasonable agreement with the measured values of S. A possible physical mechanism maintaining S is discussed. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: When pure solvent is separated from a solution of non-zero concentration Cb by a semi-permeable membrane, permeable to solvent (water) but not to solute, water flows osmotically across the membrane towards the solution. Its velocity J is given by J = PΔC, where P is a constant and ΔC is the concentration difference across the membrane. Because the osmotic flow advects solute away from the membrane, ΔC is usually less than Cb, by a factor γ which depends on the thickness of and flow in a concentration boundary layer. In this paper the layer is analysed on the assumption that the stirring motions in the bulk solution, which counter the osmotic advection, can be represented as two-dimensional stagnation-point flow. The steady-state results are compared with those of the standard physiological model in which the layer has a given thickness δ and the osmotic advection is countered only by diffusion. It turns out that the standard theory, although mechanistically inadequate, accurately predicts the value of γ over a wide range of values of the governing parameter β = PCbδ/D (where D is the solute diffusivity) if δ is given by where ν is the kinematic viscosity of the fluid and α is the stirring parameter. The final approach to the steady state is also analysed, and it is shown to be achieved in a time scale (D/ν)1/3/αk′ where k′ is a dimensionless number whose dependence on β is computed. Moreover, if β exceeds a certain critical value (≈ 10), the approach to the steady state is not monotonic but takes the form of a damped oscillation (in practice, however, β is unlikely to rise significantly above 1). The theory is extended to the case where the solute concentration is non-zero on both sides of the membrane and in that case it is shown that J is bounded as β → ∞. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: A rational asymptotic theory describing the perturbed flow in a turbulent boundary layer encountering a small two-dimensional hump is presented. The theory is valid in the limit of very high Reynolds number in the case of an aerodynamically smooth surface, or in the limit of small drag coefficient in the case of a rough surface. The method of matched asymptotic expansions is used to obtain a multiple-structured flow, along the general lines of earlier laminar studies. The leading-order velocity perturbations are shown to be precisely the inviscid, irrotational, potential flow solutions over most of the domain. The Reynolds stresses are found to vary across a thin layer adjacent to the surface, and display a singular behaviour near the surface which needs to be resolved by an even thinner wall layer. The Reynolds stress perturbations are calculated by means of a second-order closure model, which is shown to be the minimum level of sophistication capable of describing these variations. The perturbation force on the hump is also calculated, and its order of magnitude is shown to depend on the level of turbulence closure; a cruder turbulence model gives rise to spuriously large forces. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: Steady potential flow around a two-dimensional bubble with surface tension, either free or attached to a wall, is considered. The results also apply to a liquid drop. The flow and the bubble shape are determined as functions of the contact angle β and the dimensionless pressure ratio γ = (pb − ps)/½ρU2. Here pb is the pressure in the bubble, ps = p∞ + ½ρU2 is the stagnation pressure, p∞ is the pressure at infinity, ρ is the fluid density and U is the velocity at infinity. The surface tension σ determines the dimensions of the bubble, which are proportional to 2σ/ρU2. As γ tends to ∞, the bubble surface tends to a circle or circular arc, and as γ decreases the bubble elongates in the direction normal to the flow. When γ reaches a certain value γ0(β), opposite sides of the bubble touch each other. The problem is formulated as an integrodifferential equation for the bubble surface. This equation is discretized and solved numerically by Newton's method. Bubble profiles, the bubble area, the surface energy and the kinetic energy are presented for various values of β and γ. In addition a perturbation solution is given for γ large when the bubble is nearly a circular arc, and a slender-body approximation is presented for β ∼ ½π and γ ∼ γ0(β), when the bubble is slender. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: The evolution of the shape of a slender inviscid drop in an axisymmetric straining motion is studied at low Reynolds numbers. It is found that the shape equation can be solved by polynommals with time-dependent coefficients. A global stability result can be used to show simply that only one possible equilibrium is stable. It is further shown that if the slender drop starts with a long-wavelength waist then it cannot go to this stable equilibrium and must either extend indefinitely or burst. In the class of trinomial shapes, it is shown that the drop either bursts or goes to the stable equilibrium, depending on whether or not the initial shape has a waist. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: Conditions are found for the appearance of non-uniform progressive waves of permanent form from a long-wave modulation of a finite-amplitude Stokes wave on deep water. The waveheight at which the modulated waves can occur is a very slowly decreasing function of the modulation wavelength for values up to 150 times the original wavelength. Some qualitative remarks are made about the problem of determining the stability of the new waves. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: Based on the parabolic approximation, a refraction—diffraction model for linear water waves is developed. With the assumption that the water depth (refraction index) is slowly varying, the model equation describes the forward-scattered wavefield. Two examples are considered in particular: (i) wave diffraction by a long thin barrier on a uniform slope, and (ii) wave convergence over a semicircular step shoal. For the former problem, a similarity solution in terms of Fresnel integrals is obtained for the wavefield in the neighbourhood of the shadow boundary. For the latter problem, the resulting Schrödinger equation is solved numerically. The wavefield near the caustics as well as in the shadow region is obtained and compared with experimental data. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: The coherent structure dynamics in the near field of a circular jet has been experimentally explored by inducing ‘stable’ vortex pairing through controlled excitation (see Zaman & Hussain 1980) and applying phase-averaging techniques. Hot-wire measurements were made in a 7·62 cm air jet with laminar exit boundary layer at the Reynolds number ReD = 3·2 × 104, excited at the Strouhal number StD = 0·85. At a particular phase during the pairing process, spatial distributions of the phase-average longitudinal and lateral velocity perturbations (〈u)〉, 〈v〉), vorticity, streamlines, the coherent and background Reynolds stresses and turbulence intensities have been educed. These data have been obtained for four different locations occupied by the vortices at the same phase (preceding, during, and following the pairing event), in the region 0 〈 x/D 〈 5. Spatial distributions of these measures at four successive phases during the pairing process are also educed in an attempt to further understand the vortex-pairing dynamics. The flow physics is discussed on the basis of measurements over the physical extent of the vortical structures, phase-locked to specific phases of the pairing event and thus do not involve use of the Taylor hypothesis. The computed pseudostream functions at particular phases are compared with the corresponding streamlines drawn by the method of isoclines. Transition of the vortices is examined on the basis of vorticity diffusion, the superimposed random fluctuation field intensities and Reynolds stress and phase-locked circumferential correlation measurements. The peak vorticity drops rapidly owing to transition and interaction of the vortices during pairing but, farther downstream, the decay can be attributed to destruction of the coherent vorticity by the background turbulence Reynolds stress, especially at the locations of the latter's ‘saddle points’. Controlled excitation enhances the initial circumferential coherence of the vortical structures, but is ineffective in delaying turbulent breakdown near the end of the potential core; the breakdown appears to occur through evolution of the circumferential lobe structures. The coherent structure Reynolds stress is found to be much larger than the background turbulence Reynolds stress for 0 〈 x/D ≲ 3, but these two are comparable near the end of the jet potential core. The zone average of the coherent structure Reynolds stress over the cross-section of the merging vortex pair is much larger than that over a single vortical structure either before or after the completion of pairing. During the pairing process, such average correlations are found to be the largest at an early phase of the process while entrainment, turbulent breakdown as well as rapid diffusion of vorticity occur at a later phase. The regions of alternate positive and negative coherent Reynolds stresses associated with the structures and their interactions help explain ‘negative production’. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: Fully developed periodic flow (with non-zero mean) of a Newtonian fluid in a rigid curved tube has been investigated both numerically and experimentally. Results are reported for the mean friction factor, the amplitude ratio and phase angle between flow rate and pressure drop, the axial velocity profile, and the wall shear stress distribution. The numerical results (obtained by a finite difference method) are restricted to rather slow flows (mean Dean number [formula omitted]), while the experimental results (extracted from instantaneous flow rate-pressure drop measurements) extend up to [formula omitted]. A ‘resonant’ interaction between the axial and secondary flows at intermediate frequencies appears to be a characteristic feature of periodic flow in a curved tube. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: This paper reports on an extensive experimental study of the flows due to under-expanded axisymmetric jets impinging on flat plates. The range of plate locations extends to a point where the jet is just subsonic but the main emphasis is on the behaviour in the first shock cell. Plate inclinations from 90° to 30° were investigated by means of comprehensive surface pressure measurements and shadowgraph pictures. Wherever possible, the main features of the results have been reconstructed using inviscid analyses of the wave interactions. The flows are shown to be extremely complex due to the local structure of the free jet and, particularly, due to interactions between shock waves in the free jet and those created by the plate. In the near field, these interactions tend to be the controlling factors but at larger distances from the nozzle, mixing effects become increasingly important. The maximum pressure on the plate when it is inclined can be very much larger than when the plate is perpendicular, owing to the possibility of high pressure recoveries through multiple shock systems. Correlations are presented for some of the main features on perpendicular plates and it is shown that the integrated pressure loads for both normal and inclined plates can be predicted well by a simple momentum balance. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: When a pipe, connected to reservoirs at both ends, is subjected to an oscillating compression and expansion, fluid is alternately squeezed into and out of the reservoirs, causing a considerable amount of dissipation. This and two other related geometries (involving a symmetrical channel, and a pair of circular disks) are analysed for their flow patterns and energy dissipation, as a function of frequency. It is found that, under certain circumstances, the impedance due to transverse flow can greatly exceed the acoustical impedance (due to longitudinal flow). © 1980, Cambridge University Press. All rights reserved.

Beschreibung: The mechanism of the creation of secondary vortices behind an impulsively started circular cylinder is analysed in this paper by a higher order of accuracy numerical method. This is a combination of second-order and fourth-order compact finite difference schemes to resolve complete unsteady Navier–Stokes equations. The fourth-order compact scheme is used to calculate the Poisson equation of the stream function and the second-order alternating direction implicit scheme to resolve the vorticity transport equation. In particular, the growth of primary and secondary vortices with time is analysed for Reynolds numbers equal to 300, 550 and 1000. A single secondary vortex first appears at a Reynolds number equal to 300 on the surface of the cylinder. At R = 550, this creation is found numerically at dimensionless time t about 2·85, and this single secondary vortex is transformed into a pair of secondary vortices at t about 5. For R = 1000, two single vortices can be observed at t about 2·5, one near the separation point and another more important, easily identified in flow structure. These secondary vortices are transformed into a pair of secondary vortices at t about 4·5. A numerical analysis of the influence of the grid systems and the time step is also given. All numerical results presented here are compared with experimental visualizations. The comparison is found satisfactory. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: Measurements of velocity fluctuations in the wake of a thin non-lifting aerofoil are presented: in one set of experiments the flow was symmetrical, while in the other the upper surface of the aerofoil was roughened to increase the surface shear stress. Measurements were confined to the near wake, where the disturbed region lies within the inner layers of the original boundary layers; thus the boundary-layer thickness is not a relevant length scale. The practical relevance of the experiment is to the prediction of flow over aerofoils, where only the initial region of the wake significantly affects the aerofoil pressure distribution. Temperature-conditioned sampling techniques were used, one boundary layer at a time being heated so that fluid from each boundary layer could be traced within the wake. In contrast to the behaviour of merging shear layers in ducts and jets, the wake interaction involves significant fine-scale mixing; the results reveal a three-layer structure, with a fine-scale inner wake of mixed fluid separating two layers in which structural changes are confined to the region of time sharing, or internal intermittency, between mixed and unmixed fluid. The implications of the results for calculation methods are discussed. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: Decaying turbulence in neutral and stratified fluids has been studied experimentally for relatively high mesh Reynolds numbers and long time-histories. The neutral case indicates an initial period decay law, q2∝ t−1, through non-dimensional time which is considerably longer than previous measurements at the same mesh Reynolds number (Re = 48260). The stratified experiment resulted in a decay rate virtually identical to that of the neutral case through Wgt/M = 275. However the decay rate sharply decreased after this time when the field of turbulence was replaced by internal gravity waves. A critical Richardson number marks the transition from the turbulence to an internal gravity wave domain. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: A generalized eddy-viscosity function νT, is introduced in order to express the Reynolds stress in an incompressible dusty gas as a linear combination of the Kronecker and rate-of-strain tensors. On the basis of Saffman's dusty-gas model a transport equation for the eddy viscosity is derived from the general turbulence energy equations, thereby introducing two additional functions, the specific turbulence kinetic energy E1, and a scale variable s. In order to determine the three variables modified Prandtl–Wieghardt relation among them is accepted and a transport equation for s is postulated in the same manner as in the clean-gas turbulence transport model (firstly proposed by Harlow & Nakayama 1967) but with the inclusion of an additional term accounting for the dust particles stabilizing action. We are considering values of loading (mass ratio of particles) of order of unity, with particle/gas density ratios of order of 103 and volume concentrations of the order of 10−3, so that particle–particle interactions are neglected. Supposing that the particles nearly follow the gas motion, following well at large scales and poorly at small, an application of the theory to problem of numerical calculations of the dusty-gas parameters such as mean velocity profile of turbulent pipe flow is given. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: Measurements of the thermal and velocity structure of the near-surface mixing layer of a freshwater lake in moderate wind conditions from fixed or mobile arrays of sensors reveal large-scale coherent structures consisting of narrow fronts across which both the temperature and the horizontal component of the current increase. These fronts are generally transverse to the wind direction and are inclined to the vertical, and appear to be similar to fronts, reported as temperature ‘ramps’, in the near-surface atmospheric boundary layer. The time derivatives of the temperature are skewed in a sense consistent with observations in laboratory and atmospheric boundary layers, and of a magnitude consistent with measurements in the latter. Evidence is presented to show that bubbles generated by breaking waves are carried down in the large-scale pattern of flow associated with the fronts in the mixing layer. The presence of a Langmuir circulation associated with wind rows has not been established in these experiments. The relevance of the observations to the ocean mixing layer is discussed. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: The convective instability of a layer of fluid heated from below is studied on the assumption that the flux of heat through the boundaries is unaffected by the motion in the layer. It is shown that when the heat flux is above the critical value for the onset of convection, motion takes place on a horizontal scale much greater than the layer depth. Following Childress & Spiegel (1980) the disparity of scales is exploited in an expansion scheme that results in a nonlinear evolution equation for the leading-order temperature perturbation. This equation which does not depend on the vertical co-ordinate, is solved analytically where possible and numerically where necessary; most attention is concentrated on solutions representing two-dimensional rolls. It is found that for any given heat flux a continuum of steady solutions is possible for all wave numbers smaller than a given cut off. Stability analysis reveals, however, that each mode is unstable to one of longer wavelength than itself, so that any long box will eventually contain a single roll, even though the most rapidly growing mode on linear theory has much shorter wavelength. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: A model is presented of the two-dimensional boundary-layer and interior flow in a rectangular box resulting from the application of a quadratic temperature variation on its lower surface. The other walls are insulating. It is shown that a similarity form exists for the narrow, thermocline-like layer near the lower surface, and that this can satisfy all known consistency conditions with the interior, together with either laminar or turbulent side-wall regions. The interior temperature and Nusselt number are shown to be insensitive to Prandtl number, and to be primarily functions of the horizontal Rayleigh number RaL. Specifically, the interior temperature, relative to the coldest applied value, is 60% of the total applied temperature range. The Nusselt number is predicted to vary as [formula omitted] for a box with unit aspect ratio. The dynamics of the side-wall region, and the details of the imposed temperature variation appear to be unimportant in determining the overall buoyancy exchange. The solutions are compared with numerical and observational results, and generally good agreement is found. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: Energy stability theory has been formulated for two-dimensional buoyancy–thermocapillary convection in a layer with a free surface. The theory yields a critical Rayleigh number RE for which R 〈 RE is a sufficient condition for stability of the layer. RE emerges from the variational formulation as an eigenvalue of a nonlinear system of Euler–Lagrange equations. For the case of small capillary number (large mean surface tension) explicit values are obtained for RE. The analogous linear-theory results for this case are obtained in terms of a critical Rayleigh number RL. These are compared. It is found that the existence of the deformable interface can lead to a stabilization relative to the case of a planar interface. This result is explained in physical terms. The energy theory is then generalized to include general flow problems having three-dimensional disturbances, non-Newtonian bulk fluids and general interfacial mechanics such as surface viscosity and elasticity. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: Linear stability analysis is applied to the problem of a density-stratified fluid contained in an inclined slot being subjected to a lateral temperature gradient. Stability equations are solved using the Galerkin technique with 12 terms in the truncated expansion series. Within the range of θ considered, |θ| 〈 75°, critical instability was found to be of the stationary type. Results of critical thermal Rayleigh numbers and wavenumbers at all inclination angles are in good agreement with the experimental results obtained earlier (Paliwal & Chen 1980). Contrary to intuition, these results show that the system is more stable when the lower wall is heated. This is shown to be the result of the increased vertical solute gradient in the steady state prior to the onset of instabilities when the heating is from below. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: Results from idealized ocean models indicate that equatorially trapped baroclinic waves incident on an eastern boundary may be partially transmitted north and south along the coast as boundary-trapped internal Kelvin waves. The offshore scale of the coastal internal Kelvin waves is the internal Rossby radius of deformation δR, which decreases as the Coriolis parameter f increases. The effect of the presence of a continental slope of width Ls, along a north–south oriented coastline, on the poleward propagation of coastal trapped internal Kelvin waves is studied in a two-layer β-plane model. The waves propagate from regions near the equator where δR 〉 Ls to mid-latitudes where δR 〈 Ls. It is assumed that f varies slowly on the alongshore scale of the waves L, that L ≫ Ls, and that either the topographic slope is weak or that the upper-layer depth is small compared to the lower-layer depth. All of the coastal trapped waves present in the model are non-dispersive. For most values of f, the cross-shelf eigen-functions consist of the internal Kelvin wave and an infinite set of continental shelf waves whose vertical structure depends on δR/Ls. For δR/Ls ≫ 1, the shelf waves are bottom trapped while for δR/Ls ≪ 1 they are barotropic. The wave speeds Cn of the shelf waves vary linearly with f, whereas the wave speed c0 of the internal Kelvin wave is independent of f. As f increases through critical values fCn, where Cn approaches C0, the phase speeds and the eigenfunctions vary so that the eigenfunctions represent a different type of wave on either side of fCn. In the slowly-varying approximation, the alongshore energy flux in each eigenfunction is a constant. It follows that an internal Kelvin wave which has a wavelength short enough that the slowly-varying approximation remains valid and which propagates poleward from the equatorial region where f 〈 fC1 will transform into a shelf wave, at values of f near fC1, and will continue propagation poleward in that form. As a result, coastal trapped baroclinic disturbances may be able to propagate efficiently from the equatorial region to mid-latitudes where they may take the form of barotropic shelf waves. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: It is shown that all parallel inviscid shear flows of constant density are unstable to a wide class of initial infinitesimal three-dimensional disturbances in the sense that, according to linear theory, the kinetic energy of the disturbance will grow at least as fast as linearly in time. This can occur even when the disturbance velocities are bounded, because the streamwise length of the disturbed region grows linearly with time. This finding may have implications for the observed tendency of turbulent shear flows to develop a longitudinal streaky structure. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: We consider the onset of Rayleigh–Bénard convection from random fluctuations arising within a fluid. In the specific case in which the fluctuations are thermodynamically determined, we reduce the problem to a random initial value problem for the Fourier modes. For the case of weak nonlinear convection, it is possible to truncate the number of modes and this truncated set is solved both by a Monte Carlo technique and by moment methods for various Rayleigh numbers. We find three stages in the evolution of ordered convection from random fluctuations which correspond to time intervals in which the fluctuations and the nonlinearity have different degrees of importance. It is shown that no simple moment truncation method will succeed and that the time for onset of convection is a mean over a distribution of times for which members of an ensemble exhibit appreciable convective transport. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: The thermodynamics of irreversible processes is normally limited to processes that can be adequately described by linear constitutive relations, like those of Fourier and Newton in a simple gas. In this paper we use thermodynamic arguments to derive the (nonlinear) Burnett equations for a monatomic gas, thus avoiding the complicated kinetic theory by which the equations were discovered and which somewhat obscures the origin of the various terms in the equations. Expressions are given for the entropy, its flux and its production rate correct to second-order in Knudsen number. The theory involves five phenomenological parameters, and as there are eleven coefficients in the second-order terms of Burnett's equations, we are able to deduce several necessary constraints between these coefficients. Compact forms for the equations are found that clarify their physical significance. The general method we have developed is applicable to media other than simple gases. In a final section we use our theory of Burnett's equations to draw some general conclusions concerning the second law of thermodynamics. It is shown that the Clausius-Duhem inequality holds only for the linear theory of constitutive relations; and that axiomatic generalizations of the inequality to nonlinear processes – common in continuum mechanics–fail because the vital distinction between reversible and irreversible processes is not made. © 1980, Cambridge University Press. All rights reserved.

Beschreibung: An electrically conducting Boussinesq fluid is confined between two rigid horizontal planes. The fluid is permeated by a strong uniform horizontal magnetic field and the entire system rotates rapidly about a vertical axis. By heating the fluid from below and cooling it from above the system becomes unstable to small perturbations when the adverse temperature gradient becomes sufficiently large. Attention is restricted to small values of the Ekman number E and the ratio q of the thermal and magnetic diffusivities (see (1.2) and (1.3) below). In this parameter range marginal convection is steady and its character depends on the relative sizes of the Coriolis and Lorentz forces as measured by the parameter λ (see (1.1) below). When λ ≥ 2/3½, motion consists of a single roll, whose axis is perpendicular to the applied magnetic field. On the other hand, when λ 〈 2/3½, two distinct rolls are possible: the axis of each roll lies oblique but with equal angle to the applied magnetic field. Only the latter case is discussed here. Once the Rayleigh number R exceeds its critical value Rc only one of the two sets of single rolls remains stable, while its squared amplitude increases linearly with R – Rc. For certain values of the parameters λ and τ (see (1.6) below) a second critical value may be isolated at which the system becomes unstable to unidirectional geostrophic flow perturbations aligned with the applied magnetic field. The instability sets in as either a steady or oscillatory shear flow dependent on the values taken by λ and τ. In both cases, after the secondary instability sets in, the roll amplitude remains largely insensitive to further increase in the Rayleigh number with the consequence that the geostrophic flow is stabilized. The amplitude of the shear, on the other hand, increases with R, adjusting its magnitude to ensure stability of the convection rolls. © 1980, Cambridge University Press. All rights reserved.