ALBERT

Description: A coherent-vortex analysis is made of a computational solution for the free decay of homogeneous, Charney-isotropic geostrophic turbulence at large Reynolds number. The method of analysis is a vortex detection and measurement algorithm that we call a vortex census. The census demonstrates how, through non-conservative interactions among closely approaching vortices, the vortex population evolves towards fewer, larger, sparser, and more weakly deformed vortices. After emergence from random initial conditions and a further period of population adjustment, there is a period of approximately self-similar temporal evolution in the vortex statistics. This behaviour is consistent with a mean-vortex scaling theory based on the conservation of energy, vortex extremum, and vortex aspect ratio. This period terminates as the population approaches a late-time non-turbulent end-state vortex configuration. The end state develops out of merger and alignment interactions among like-sign vortices, and even during the scaling regime, local clusters of nearly aligned vortices are common.

Description: When two drops of radius R touch, surface tension drives an initially singular motion which joins them into a bigger drop with smaller surface area. This motion is always viscously dominated at early times. We focus on the early-time behaviour of the radius rm of the small bridge between the two drops. The flow is driven by a highly curved meniscus of length 2πrm and width Δ ≪ rm around the bridge, from which we conclude that the leading-order problem is asymptotically equivalent to its two-dimensional counterpart. For the case of inviscid surroundings, an exact two-dimensional solution (Hopper 1990) shows that Δ ∝ rm3 and rm ∼ (tγ/πη) ln [tγ/(ηR)]; and thus the same is true in three dimensions. We also study the case of coalescence with an external viscous fluid analytically and, for the case of equal viscosities, in detail numerically. A significantly different structure is found in which the outer-fluid forms a toroidal bubble of radius Δ ∝ rm3/2 at the meniscus and rm ∼ (tγ/4πη) ln [tγ/(ηR)]. This basic difference is due to the presence of the outer-fluid viscosity, however small. With lengths scaled by R a full description of the asymptotic flow for rm(t) ≪ 1 involves matching of lengthscales of order rm2 rm3/2, rm, 1 and probably rm7/4.

Description: In this paper we investigate the linear stability of detonations in which the underlying steady one-dimensional solutions are of the pathological type. Such detonations travel at a minimum speed, which is greater than the Chapman-Jouguet (CJ) speed, have an internal frozen sonic point at which the thermicity vanishes, and the unsupported wave is supersonic (i.e. weak) after the sonic point. Pathological detonations are possible when there are endothermic or dissipative effects present in the system. We consider a system with two consecutive irreversible reactions A→B→C, with an Arrhenius form of the reaction rates and the second reaction endothermic. We determine analytical asymptotic solutions valid near the sonic pathological point for both the one-dimensional steady equations and the equations for linearized perturbations. These are used as initial conditions for integrating the equations. We show that, apart from the existence of stable modes, the linear stability of the pathological detonation is qualitatively the same as for CJ detonations for both one- and two-dimensional disturbances. We also consider the stability of overdriven detonations for the system. We show that the frequency of oscillation for one-dimensional disturbances, and the cell size based on the wavenumber with the highest group velocity for two-dimensional disturbances, are both very sensitive to the detonation speed for overdriven detonations near the pathological speed. This dependence on the degree of overdrive is quite different from that obtained when the unsupported detonation is of the CJ type.

Description: This paper considers the role of long finite-amplitude Rossby waves in determining the evolution of flow along a rapidly rotating channel with an uneven floor. The Rossby waves travel on a potential vorticity interface in a channel with a cross-channel step change in depth, where step position varies slowly along the channel. A nonlinear wave equation is derived describing the evolution of the potential vorticity interface. To leading order this is the hydraulic equation derived by Haynes, Johnson & Hurst (1993). Dispersion appears at the next order. Various solution regimes are identified. As well as slowly varying hydraulic solutions, two further types of steady solutions appear: approach-controlled flows and twin supercritical leaps. Both these solutions are characterized by leaps between supercritical branches of the hydraulic function. It is shown how the position and size of these 'supercritical leaps' can be determined within the context of hydraulic theory. To fully resolve the internal structure of these leaps dispersive effects must be included and leaps are shown to correspond to kink soliton solutions of the steady unforced problem. It is also shown that increasing dispersion (decreasing topographic length scale) causes the loss of the subcritical solution branch in some subcritical flows. The only candidate for a steady solution in these regimes is then an approach-controlled flow. Integrations of initial value problems show that in general flows evolve towards the dispersive form of the solution predicted by hydraulic theory, at least near the topographic perturbation. However, in those subcritical flows where sufficiently large dispersion causes the subcritical branch to disappear, unsteady integrations evolve to approach-controlled flows even when the dispersion is sufficiently small that the subcritical branch still exists.

Description: The transient evolution of the bubble-size probability density functions resulting from the breakup of an air bubble injected into a fully developed turbulent water flow has been measured experimentally using phase Doppler particle sizing (PDPA) and image processing techniques. These measurements were used to determine the breakup frequency of the bubbles as a function of their size and of the critical diameter Dc defined as Dc = 1.26 (σ/ρ)3/5∈-2/5, where ∈ is the rate of dissipation per unit mass and per unit time of the underlying turbulence. A phenomenological model is proposed showing the existence of two distinct bubble size regimes. For bubbles of sizes comparable to Dc, the breakup frequency is shown to increase as (σ/ρ)-2/5 ∈3/5 √D/Dc-1, while for large bubbles whose sizes are greater than 1.63Dc, it decreases with the bubble size as ∈1/3 D-2/3. The model is shown to be in good agreement with measurements performed over a wide range of bubble sizes and turbulence intensities.

Description: Direct numerical simulations of turbulence resulting from Kelvin-Helmholtz instability in stratified shear flow are used to examine the geometry of the dissipation range in a variety of flow regimes. As the buoyancy and shear Reynolds numbers that quantify the degree of isotropy in the dissipation range increase, alignment statistics evolve from those characteristic of parallel shear flow to those found previously in studies of stationary, isotropic, homogeneous turbulence (e.g. Ashurst et al. 1987; She et al. 1991; Tsinober et al. 1992). The analysis yields a limiting value for the mean compression rate of scalar gradients that is expected to be characteristic of all turbulent flows at sufficiently high Reynolds number. My main focus is the value of the constant q that appears in both the Batchelor (1959) and Kraichnan (1968) theoretical forms for the passive scalar spectrum. Taking account of the effects of time-dependent strain, I propose a revised estimate of q, denoted qe, which appears to agree with spectral shapes derived from simulations and observations better than do previous theoretical estimates. The revised estimate is qe = 7.3±0.4, and is expected to be valid whenever the buoyancy Reynolds number exceeds O(102). The Kraichnan (1968) spectral form, in which effects of intermittency are accounted for, provides a better fit to the DNS results than does the Batchelor (1959) form.

Description: A methodology termed the 'filtered mass density function' (FMDF) is developed and implemented for large-eddy simulation (LES) of variable-density chemically reacting turbulent flows at low Mach numbers. This methodology is based on the extension of the 'filtered density function' (FDF) scheme recently proposed by Colucci et al (1998) for LES of constant-density reacting flows. The FMDF represents the joint probability density function of the subgrid-scale (SGS) scalar quantities and is obtained by solution of its modelled transport equation. In this equation, the effect of chemical reactions appears in a closed form and the influences of SGS mixing and convection are modelled. The stochastic differential equations (SDEs) which yield statistically equivalent results to those of the FMDF transport equation are derived and are solved via a Lagrangian Monte Carlo scheme. The consistency, convergence, and accuracy of the FMDF and the Monte Carlo solution of its equivalent SDEs are assessed. In non-reacting flows, it is shown that the filtered results via the FMDF agree well with those obtained by the 'conventional' LES in which the finite difference solution of the transport equations of these filtered quantities is obtained. The advantage of the FMDF is demonstrated in LES of reacting shear flows with non-premixed reactants. The FMDF results are appraised by comparisons with data generated by direct numerical simulation (DNS) and with experimental measurements. In the absence of a closure for the SGS scalar correlations, the results based on the conventional LES are significantly different from those obtained by DNS. The FMDF results show a closer agreement with DNS. These results also agree favourably with laboratory data of exothermic reacting turbulent shear flows, and portray several of the features observed experimentally.

Description: The objective of this study is to develop a method of controlling vortex shedding behind a bluff body using control theory. A suboptimal feedback control procedure for local sensing and local actuation is developed and applied to the flow behind a circular cylinder. The location of sensors for feedback is limited to the cylinder surface and the control input from actuators is the blowing and suction on the cylinder surface. Three different cost functionals to be minimized (J1 and J2) or maximized (J3) are investigated: J1 is proportional to the pressure drag of the cylinder, J2 is the square of the difference between the target pressure (inviscid flow pressure) and real flow pressure on the cylinder surface, and J3 is the square of the pressure gradient on the cylinder surface, respectively. Given the cost functionals, the flow variable to be measured by the sensors and the control input from the actuators are determined from the suboptimal feedback control procedure. Several cases for each cost functional have been numerically simulated at Re = 100 and 160 to investigate the performance of the control algorithm. For all actuations, vortex shedding becomes weak or disappears, and the mean drag and drag/lift fluctuations significantly decrease. For a given magnitude of the blowing/suction, reducing J2 provides the largest drag reduction among the three cost functionals.