ALBERT

Description: The following radiocarbon date list contains dates of samples measured since our previous list (R, 1979, v 21, p 131-137). As before, age calculations are based on the Libby half-life 5570 ± 30 yr and reported in years before 1950. The modern standard is 0.95 of the activity of NBS oxalic acid. Sample pretreatment, combustion, and counting technique are essentially the same as described in R, 1971, v 13, p 135-140, supplemented by new techniques for groundwater processing (R, 1979, v 21, p 131-137) and for soil sample treatment (R, 1977, v 19, p 465-475).

Description: This date list includes most of the archaeologic and geologic samples dated in this laboratory since publication of our last date list (R, 1978, v 20, p 210–233), as well as some samples dated previously which lacked adequate sample information.

Description: Samples in this paper are a continuation of the Florida geologic samples, reported earlier (R, 1978, v 20, p 436–440). Each specimen analyzed in this list was a carbonate shell. All dates are from single shells as opposed to means of multiple shell samples.

Description: The Menlo Park radiocarbon laboratory of the US Geological Survey began operation in January 1976 in a new building with a counting room 9.8m below the ground surface. Shielding above the counters consists of 30cm lead and 8.2m earth. The ionizing cosmic ray flux at the counters is 0.28 cpm/cm2. Background levels obtained thus far are summarized in Table 1.

Description: Recent experimental work has suggested that the relative fractionation of 14C to 13C may differ from the accepted value of b = 2. In order to explore the implications of this possibility, the standard formulae for correcting radiocarbon dates for fractionation effects are rederived, but without making any of the usual assumptions or approximations. A generalized dating equation is derived (where ASN and AON are normalized sample and standard activities, β is a factor which reflects changes in atmospheric 13C and 14C content, {RST(o)/RST}b accounts for post-depositional changes in sample 13C ratio, and tcal is calendar age in years before ad 1950. The errors in calculated ages which might arise from different b values are estimated and shown to be small relative to other dating uncertainties. The effect of b ≠ 2 may be important in the calibration of radiocarbon dates using tree-ring samples of known age. A theoretical analysis suggests that b ≠ 2 effects may result in a correlation between age anomaly (ie, the difference between radiocarbon age and calendar age) and sample 13C data. However, an analysis of published data reveals no meaningful correlation. This result, while not eliminating the possibility that b ≠ 2, highlights its unimportance even in this high-precision application of radiocarbon dating.

Description: The carbon content, pH and 14C concentration of humic acids were determined for three soil series of Arctic and Subarctic ecosystems. The measured 14C ages were interpreted in the light of an equilibrium model of humus formation and of mineralization processes in recent soils, and the coefficient of renovation, Kr, was calculated for humic acids. The comparison of Kr for series formed under different climatic conditions suggested that global warming could accelerate decomposition of soil organic matter and possibly increase productivity of ecosystems of the Arctic region.

Description: Δ14C records are reported for post-bomb corals from three sites in the tropical Atlantic Ocean. In corals from 18°S in the Brazil Current, Δ14C values increased from ca. −58% in the early 1950s to +138% by 1974, then decreased to 110‰ by 1982. Shorter records from 8ºS off Brazil and from the Cape Verde Islands (17°N) showed initially higher Δ14C values before 1965 than those at 18ºS, but showed lower rates of increase of Δ14C during the early 1960s. There is general agreement between the coral results and Δ14C of dissolved inorganic carbon (DIC) measured in seawater previously for locations in the tropical Atlantic Ocean. Δ14C values at our tropical ocean sites increased at a slower rate than those observed previously in the temperate North Atlantic (Florida and Bermuda), owing to the latter's proximity to the bomb 14C input source in the northern, hemisphere. Model results show that from 1960–1980 the Cape Verde coral and selected DIG Δ14C values from the North Equatorial Current agree with that calculated for the North Atlantic based on an isopycnal mixing model with a constant water mass renewal rate between surface and subsurface waters. This is in contrast to Δ14C values in Bermuda corals that showed higher post-bomb values than those predicted using a constant water mass renewal rate, hence indicating that ventilation in the western north Atlantic Ocean had decreased by a factor of 3 during the 1960s and 1970s (Druffel 1989).

Description: Although soil contains about three times the amount of carbon present in the preindustrial atmosphere, determining how perturbations (e.g., changing land use, CO2 fertilization, changing climate and anthropogenic nitrogen deposition) alter soil carbon storage and influence atmospheric CO2 levels has proved elusive. Not knowing the soil carbon turnover times causes part of this uncertainty. I outline a strategy for using radiocarbon measurements to estimate soil organic matter turnover times and inventories in native soil. The resulting estimates of carbon exchange produce reasonable agreement with measurements of CO2 fluxes from soil. Furthermore, derivatives of the model are used to explore soil carbon dynamics of cultivated and recovering soil. Because the models can reproduce observed soil 14C measurements in native, cultivated, and recovering ecosystems (i.e., the underlying assumptions appear reasonable), the native model was modified to estimate the potential rate of additional carbon storage because of CO2 fertilization. This process may account for 45–65% of the “missing CO2 sink.”

Description: Recently, in an intercomparison of the Hohenheim German oak chronology (Becker 1993) and the Göttingen chronology (Leuschner and Delorme 1988), an error was discovered in the former (Leuschner, in preparation). Due to an error in adding sections at 5241 bc, 41 yr are missing in the published Hohenheim chronology. After correction of the error, the two chronologies synchronize over their entire common length, back to 7200 bc.

Description: Extensive measurements of radiocarbon have been used in New Zealand since the mid-1960s to follow carbon (C) movement and turnover in soils. We present here unpublished radiocarbon (14C) measurements on a range of eight New Zealand soils with details of the sites, ecosystems, climates, soil descriptions and associated analytical data. An overview is also given of published 14C measurements on soils, and the use of these measurements to model soil C turnover.

Description: We report here the remainder of the Hamburg University dates on thin soil layers (HAM 1652–3129).

Description: In order to verify sediment trap samples as indicators of upper ocean 14C concentrations, particulate inorganic radiocarbon (PICΔ14C) collected by time-series sediment traps in the Sea of Okhotsk and the Bering Sea was measured by accelerator mass spectrometry (AMS). All of the PICΔ14C measurements were 〈 0‰, in contrast to GEOSECS 14C data in the upper ocean from the northwestern North Pacific. This difference is attributed to the upwelling of deepwater that contains low Δ14C of dissolved inorganic carbon (DICΔ14C) and to the decrease over time of surface DICΔ14C owing to the decrease of atmospheric Δ14C values. In addition, PICΔ14C values showed significant seasonal variability: PICΔ14C collected in the fall was the greatest (-22‰ on average), whereas PICΔ14C collected in winter showed an average minimum of −48‰. It is likely that this difference was caused by changes in mixed layer thickness. Although some uncertainties remain, further study on PICΔ14C will enable us to estimate seasonal variability in DICΔ14C and air-sea CO2 exchange rate.

Description: For several undisturbed sites in Germany, 14C data are reported for soil organic matter (SOM) (4 sites), soil CO2 (10 sites) and dissolved organic carbon (DOC) (1 site). With the assumption of a fast degradable component (lifetime ca. 1 yr) and a slow degradable component (lifetime ca. 100 yr), a range between 0.6 and 1.6 mm yr-1 has been determined for the downward migration rates of soil organic carbon at the sampling sites from the soil 14C data. The soil CO2 measurements show that in deciduous forests the fast degradable component is ca. 60% and the slow degradable component is ca. 40% of the SOM. In coniferous forests this ratio is reversed. The 14C results for DOC could not be explained with the assumption of a first-order decay process. The removal of soil organic carbon by DOC is of minor importance for the estimation of carbon budgets for the investigated site.

Description: Procedures and equipment have been described in previous date lists. Except as otherwise indicated, wood, charcoal, and peat samples are pretreated with dilute NaOH-Na4P2O7 and dilute H3PO4 before conversion to the counting gas methane; marls and lake cores are treated with acid only. Very calcareous materials are treated with HCl instead of H3PO4.

Description: Most of the 14C measurements reported here were made between October 1979 and October 1980. Equipment, measurement, and treatment of samples are as reported previously (R, 1968, v 10, p 36-37; 1976, v 18, p 290; 1980, v 22, p 1045).

Description: Routine radiocarbon analyses were last reported for the Laboratory of Isotope Geochemistry at the University of Arizona in 1971 (Haynes, Grey, and Long, 1971), and a special date list on packrat middens appeared in 1978 (Mead, Thompson, and Long, 1978). This list presents results obtained from our gas proportional counting facility before its major renovation and before the addition of a liquid scintillation counting system. The characteristics of these new systems will be described in the next date list.

Description: The results presented in this date list are from the same sampling program previously reported (R, 1978, v 20, p 436–440). Each sample reported herein was geologic in origin and, in most instances, represents a single shell date.

Description: The extension of the bristlecone pine chronology, which is, thus far, the only complete tree-ring-dated series that extends beyond 6000 bc and which is used for the calibration of radiocarbon dates, has been going on over the past 22 years, that is, ever since Schulman (1956; 1958) established the extremely old ages of the bristlecone pine. The few people in this world who search for very old wood samples to extend a dendrochronology may appreciate a method of establishing quickly the approximate ages of the samples.

Description: The following list consists entirely of dates for archaeologic samples from the British Isles, mainly measured over the period from May 1979 to June 19801. The dates were obtained by liquid scintillation counting of benzene using the laboratory procedures outlined in previous lists (see, eg, BM-VIII, R, 1976, v 18, p 16).

Description: This list contains most of the measurements made during 1979, since our last list (R, 1979, v 21, p 180–185). The laboratory procedures were outlined previously (R, 1968, v 10, p 29–34; R, 1971, v 13, p 29–31). The collagen extraction follows Longin (1970).

Description: The following radiocarbon dates are a partial list of samples measured for a variety of projects and materials since May 1980. Chemical and counting procedures remain the same as indicated in R, v 20, p 274–282.

Description: This list contains the results of measurements made during 1978 and 1979. Ages were computed on the radiocarbon half-life of 5568 ± 30 years. Statistical errors quoted herein are 1σ counting errors which include sample, background, and standard. The age limit reported is calculated on the basis of 3σ activity. The σ13C values in table 1 were measured by Krueger Enterprises, Inc, Geochron Laboratories Division1 and calculated based on Craig's Peedee Belemnite (PDB) limestone standard (Craig, 1957). The total alkalinity as bicarbonate values reported in table 1 was determined using techniques described by Brown, Skougstad, and Fishman (1970). Unless otherwise stated, all samples were collected and submitted by personnel of the US Geological Survey.

Description: The dates presented below are from some important archaeologic and Quaternary sites. All the dates are in years bp, based on τ½ = 5568 years. When converting archaeologic dates into the ad/bc scale, 1950 should be used as the base year. The dates are not corrected for 13C fractionation. All the dates older than 10,000 years have been given with 2σ errors.

Description: Measurements of cosmogenic nuclides made in situ in the Earth's surface are being used to help resolve a wide range of geologic and chronologic questions. Cosmogenic nuclides (3He, 10Be, 14C, 21Ne, 26Al 36C1 are presently used) can reveal rock exposure history information leading to estimates of timing of surface forming events, rates and styles of erosion, and timing and durations of episodes of burial. Depending on the problems being tackled, a significant source of error (±10–25%) for any cosmogenic nuclide method is the present uncertainty in the spatial and temporal variability of the rates of production of these in-situ nuclides.

Description: We present a compilation of 14C soil dates measured at the University of Hamburg through 1984 (HAM-1597).

Description: At the University of Miami Tritium Laboratory and the University of Washington Quaternary Isotope Laboratory, more than 1000 large-volume Pacific Ocean radiocarbon samples were measured for the WOCE program. Here we present a comprehensive data set, and a brief discussion of our findings.

Description: Bomb 14C from nuclear tests in the atmosphere has proved to be a particularly useful tool in the study of the carbon cycle. We provide here a ca. 30-yr time series of 14C concentrations in the atmosphere between 28°N and 71°N and in the ocean surface between 45°S and 45°N. More recently (since 1990), a north-south profile also has been obtained for 14C in the surface waters of the Atlantic Ocean. The measurements were performed using the conventional technique of beta counting or large samples (4 to 5 liter CO2) in CO2 proportional counters. These data show that the 14C concentration in the atmosphere is leveling off with a time constant of 0.055 yr-1, and is now approaching that of the ocean surface at lower latitudes.Additional tracer studies have been concerned especially with the penetration of bomb 14C into the deep ocean. The Norwegian and Greenland seas are of interest as a sink for atmospheric CO2 and also a source of water for the deep Atlantic Ocean. During the last five years, several 14C depth profiles have been measured from the Fram Strait (79°N) to south of Iceland (62°N), using the AMS technique available at the University of Arizona AMS Facility. We considered it important to repeat and compare a few of the profiles with those produced by the GEOSECS expedition in 1972 and the TTO expedition in 1981. The profiles show that water descending to the deep Atlantic Ocean is originating mainly from intermediate and surface depths in the Nordic Seas. However, the ventilation rate of the Norwegian Sea deepwater is too slow to be an important component in the transfer of water over the Greenland-Scotland Ridge.

Description: The Geological Survey of Canada (GSC) has developed, over the past decade, a user-oriented database, Date Locator File, of Canadian samples dated by the 14C technique. This database presently contains 〉3500 soil and soil-related dates. The primary category in this suite of dates is peat, as a large portion of the Canadian landscape is covered with this type of organic soil. The data is available gratis to all researchers in a large variety of formats from simple lists to complex tables for inclusion in publications. The site localities can also be plotted on base maps suitable for publication. The database is actively augmented on an ongoing basis, but to continue to be relevant, it depends largely on the altruism of the scientific community.

Description: The following list contains samples of geologic interest that were processed from February 1974 through May 1980 at the Illinois State Geological Survey (ISGS) Radiocarbon Dating Laboratory. The archaeologic samples processed during the same period will be published in our next date list. The benzene liquid scintillation technique was used following laboratory procedures previously reported by Coleman (1973; 1974).

Description: The following radiocarbon dates are a partial list of samples measured for a variety of projects and materials since August 1980. Chemical and counting procedures remain the same as indicated in R, v 20, p 274-282.

Description: Results reported here are for samples of geologic or geographic context completed during the period 1975 to 1978.

Description: Measurements have continued with the same proportional counter system, pretreatment procedure, methane preparation and measurement, and calculation, as described previously (R, 1970, v 12, p 298–318). Uncertainties quoted are single standard deviations originating from standard, sample, background counting rates, and half-life. No 13C/12C ratios were measured. Sample descriptions have been prepared in cooperation with submitters.

Description: Acid hydrolysis is used to fractionate the soil organic carbon pool into relatively slow- and fast-cycling compartments on soils from Arizona, the Great Plains states and Michigan collected for carbon isotope tracer studies related to soil carbon sequestration, for studies of shifts in C3/C4 vegetation, and for “pre-bomb” soil-carbon inventories. Prior to hydrolysis, soil samples are first treated with cold 0.5–1N HCl to remove soil carbonates if necessary. Samples are then dispersed in a concentrated NaCl solution (ρ≍1.2 g cm-3) and floated plant fragments are skimmed off the surface. After rinsing and drying, all remaining recognizable plant fragments are picked from the soil under 20x magnification. Plant-free soils, and hot, 6N HCl acid-hydrolysis residue and hydrolyzate fractions are analyzed for carbon content, δ13C and 14C age, and the carbon distribution is verified within 1–2% by stable-carbon isotope mass balance. On average, the recalcitrant residue fraction is 1800 yr older and 2.6% more 13C-depleted than total soil organic carbon. A test of hydrolysis with fresh plant fragments produced as much as 71–76% in the acid-hydrolysis residue pool. Thus, if plant fragments are not largely removed prior to hydrolysis, the residue fraction may date much younger than it actually is.

Description: In this paper we consider the development of shear dispersion following the introduction of a diffusing tracer substance into a tube or duct containing flowing fluid, with emphasis on the characterization of the temporal variation of concentration at a fixed axial position. Asymptotic results are derived by assuming that the distance downstream of the point of tracer introduction, appropriately non-dimensionalized, is large. First, we consider the central moments of the temporal concentration variation, including their dependence on transverse position and on the initial transverse distribution of tracer. The moments for finite Péclet number are expressed in terms of their infinite-Péclet-number counterparts, and the latter are given explicitly for Poiseuille flow. Then, assuming the Péclet number is infinite, we derive an approximate solution for the Green's function expressing tracer concentration following its introduction at an arbitrary point within the tube. The solution is expressed in terms of three numerically evaluated functions of a dimensionless time variable, with parametric dependence on the distance downstream of the point of tracer release. The method is illustrated by calculation of the approximate solution for dispersion in Poiseuille flow. Unlike previous approximations, the present solution is uniformly asymptotic and represents the tails of the concentration distribution as well as the approximately Gaussian central part; in these three regions, simpler analytic forms of the approximation are given. Comparison with previous computational solutions suggests the present approximation remains reasonably accurate even at quite short distances from the point where tracer is released.

Description: The influence of different boundary conditions applied in the contact line region on the outer meniscus shape is analysed by means of a finite-element numerical simulation of the steady movement of a liquid-gas meniscus in a capillary tube. The free-surface steady shape is obtained by solving the unsteady creeping-flow approximation of the Navier-Stokes equations starting from some initial shape. Comparisons of the outer solutions obtained using two different inner models, together with that published by Lowndes (1980), indicate the relative insensitivity of the outer solution to the type of model utilized in the contact line region.

Description: The stability of steady axisymmetric convection in cylinders heated from below and insulated laterally is investigated numerically using a mixed finite-difference/Chebyshev collocation method to solve the base flow and the linear stability equations. Linear stability boundaries are given for radius to height ratios Γ from 0.9 to 1.56 and for Prandtl numbers Pr = 0.02 and Pr = 1. Depending on Γ and Pr, the azimuthal wavenumber of the critical mode may be m = 1,2,3, or 4. The dependence of the critical Rayleigh number on the aspect ratio and the instability mechanisms are explained by analysing the energy transfer to the critical modes for selected cases. In addition to these results the onset of buoyant convection in liquid bridges with stress-free conditions on the cylindrical surface is considered. For insulating thermal boundary conditions, the onset of convection is never axisymmetric and the critical azimuthal wavenumber increases monotonically with Γ. The critical Rayleigh number is less then 1708 for most aspect ratios.

Description: Since Bödewadt's (1940) seminal work on the boundary layer flow produced by a fluid in solid-body rotation over a stationary disk of infinite radius there has been much interest in determining the stability of such flows. To date, it appears that there is no theoretical study of the stability of Bödewadt's self-similar solution to perturbations that are not self-similar. Experimental studies have been compromised due to the difficulty in establishing these steady flows in the laboratory. Savaş. (1983, 1987) has studied the end wall boundary layers of flow in a circular cylinder following impulsive spin-down. During the first few radians of rotation, the endwall boundary layers have a structure very similar to Bödewadt layers. For certain conditions, Savaş has observed a series of axisymmetric waves travelling radially inwards in the endwall boundary layers. The conjecture is that these waves represent a mode of instability of the Bödewadt layer. Within a few radians of rotation however, the centrifugal instability of the sidewall layer dominates the spin-down process and the endwall waves are difficult to examine further. Here, the impulsive spin-down problem is examined numerically for Savaş' (1983, 1987) conditions and good agreement with his experiments is achieved. New experimental results are also presented, which include quantitative space-time information regarding the axisymmetric waves. These agree well with both the numerics and the earlier experimental work. Further, a related problem is considered numerically. This flow is also initially in solid-body rotation, but only the endwalls are impulsively stopped, keeping the sidewall rotating. This results in a flow virtually identical to the usual spin-down flow for the first few radians of rotation, except in the immediate vicinity of the sidewall. The sidewall layer is no longer centrifugally unstable and the circular waves on the endwalls are observed without the influence of the sidewall instability.

Description: It is well known that subgrid models such as Smagorinsky's cannot be used for the spatially growing simulation of the transition to turbulence of flat-plate boundary layers, unless large-amplitude perturbations are introduced at the upstream boundary: they are over-dissipative, and the flow simulated remains laminar. This is also the case for the structure-function model (SF) of Métais & Lesieur (1992). In the present paper we present a sequel to this model, the filtered-structure-function (FSF) model. It consists of removing the large-scale fluctuations of the field before computing its second-order structure function. Analytical arguments confirm the superiority of the FSF model over the SF model for large-eddy simulations of weakly unstable transitional flows. The FSF model is therefore used for the simulation of a quasi-incompressible (M∞ = 0.5) boundary layer developing spatially over an adiabatic flat plate, with a low level of upstream forcing. With the minimal resolution 650 × 32 × 20 grid points covering a range of streamwise Reynolds numbers Rex1 ε [3.4 × 105, 1.1 × 106], transition is obtained for 80 hours of time-processing on a CRAY 2 (whereas DNS of the whole transition takes about ten times longer). Statistics of the LES are found to be in acceptable agreement with experiments and empirical laws, in the laminar, transitional and turbulent parts of the domain. The dynamics of low-pressure and high-vorticity distributions is examined during transition, with particular emphasis on the neighbourhood of the critical layer (defined here as the height of the fluid travelling at a speed equal to the phase speed of the incoming Tollmien–Schlichting waves). Evidence is given that a subharmonic-type secondary instability grows, followed by a purely spanwise (i.e. time-independent) mode which yields peak-and-valley splitting and transition to turbulence. In the turbulent region, flow visualizations and local instantaneous profiles are provided. They confirm the presence of low- and high-speed streaks at the wall, weak hairpins stretched by the flow and bursting events. It is found that most of the vorticity is produced in the spanwise direction, at the wall, below the high-speed streaks. Isosurfaces of eddy viscosity confirm that the FSF model does not perturb transition much, and acts mostly in the vicinity of the hairpins.

Description: The interaction between a zero-pressure-gradient turbulent boundary layer and a pair of strong, common-flow-down, streamwise vortices with a sizeable velocity deficit is studied by large-eddy simulation. The subgrid-scale stresses are modelled by a localized dynamic eddy-viscosity model. The results agree well with experimental data. The vortices drastically distort the boundary layer, and produce large spanwise variations of the skin friction. The Reynolds stresses are highly three-dimensional. High levels of kinetic energy are found both in the upwash region and in the vortex core. The two secondary shear stresses are significant in the vortex region, with magnitudes comparable to the primary one. Turbulent transport from the immediate upwash region is partly responsible for the high levels of turbulent kinetic energy in the vortex core; its effect on the primary stress (úú) is less significant. The mean velocity gradients play an important role in the generation of (úú) in all regions, while they are negligible in the generation of turbulent kinetic energy in the vortex core. The pressure - strain correlations are generally of opposite sign to the production terms except in the vortex core, where they have the same sign as the production term in the budget of (úú). The results highlight the limitations of the eddy-viscosity assumption (in a Reynolds-averaged context) for flows of this type, as well as the excessive diffusion predicted by typical turbulence models.

Description: The paper addresses the mathematical modelling of the formation of a pointed drop in a four-roller mill, observed by Taylor (1934) in the Cavendish Laboratory. Since the experiments were carried out with drops of small diameter compared to the mill size, the method of matched asymptotic expansions is applicable. A two-dimensional Stokes flow generated by the rotating rollers in the mill but with no drop effect (outer problem) is computed numerically by a boundary-element method. The local expansion of that flow at the centre of the mill, where the drop is to be positioned, is used as a far field for the flow around the drop in unbounded fluid (inner problem). Employing a plane-flow model and using complex-variable techniques, the explicit solutions previously obtained by the author are adapted to the inner problem. It is proved that, with an increasing rotation rate of the rollers, the drop does develop two apparent cusps on the interface, and its shapes have striking similarities with Taylor's experiments. Response diagrams showing the drop distortion versus the elongational strain demonstrate that these are one-to-one function of each other if the drop diameter is greater than a critical value determined by the size of the mill but cease to be one-to-one otherwise. This behaviour is identified with a sudden transition from a rounded drop to a cusped one at a critical strain.

Description: The effects of molecular diffusivities of heat and mass on the counter-gradient scalar and momentum transfer in strongly stable stratification are experimentally investigated in unsheared and sheared stratified water mixing-layer flows downstream of turbulence-generating grids. Experiments are carried out in two kinds of stably stratified water flows. In the case of thermal stratification, the difference between the turbulent fluxes of an active scalar (heat with the Prandtl number of Pr ≈ 6) and a passive scalar (mass with the Schmidt number of Sc ≈ 600) is investigated. In the case of salt stratification, the effects of the molecular diffusion of the active scalar (salt) with a very high Schmidt number of Sc ≈ 600 on the counter-gradient scalar transfer is studied. Comparisons of the effects of molecular diffusivities are also made between thermally stratified water and air (Pr ≈ 0.7) flows. Further, the effects of mean shear on the counter-gradient scalar and momentum transfer are investigated for both stratified cases. Instantaneous temperature, concentration and streamwise and vertical velocities are simultaneously measured using a combined technique with a resistance thermometer, a laser-induced fluorescence method, and a laser-Doppler velocimeter with high spatial resolution. Turbulent scalar fluxes, joint probability density functions, and cospectra are estimated.The results of the first case show that both active heat and passive mass develop counter-gradient fluxes but that the counter-gradient flux of passive mass is about 10% larger than that of active heat, mostly due to molecular diffusion effects at small scales. The counter-gradient scalar transfer mechanism in stable stratification can be explained by considering the relative balance between the available potential energy and the turbulent kinetic energy as in Schumann (1987). In thermally and salt-stratified water mixing-layer flows with the active scalars of high Prandtl and Schmidt numbers, the buoyancy-induced motions with finger-like structures first contribute to the counter-gradient scalar fluxes at small scales, and then the large-scale motions, which bring fluid back to its original levels, generate the counter-gradient fluxes at large scales. The contribution of the small-scale motions to the counter-gradient fluxes in stratified water flows is quite different from that in stratified air flows. The higher Prandtl or Schmidt number of the active scalar generates both the stronger buoyancy effects and the longer time-oscillation period of the counter-gradient scalar fluxes. The time-oscillation occurs at large scales but the counter-gradient fluxes at small scales persist without oscillating. The mean shear acts to reduce the counter-gradient scalar and momentum transfer at large scales, and therefore the counter-gradient fluxes in sheared stratified flows can be seen only in very strong stratification. The behaviour of the counter-gradient momentum flux in strong stratification is quite similar to that of the counter-gradient scalar flux.

Description: The improvement of heat transfer conditions in liquid-metal magnetohydrodynamic (MHD) flows is of prime importance for self-cooled fusion blanket design concepts. For many years the research was based on stationary inertialess assumptions since it was expected that time-dependent inertial flows would be suppressed by strong electromagnetic damping, especially in the extreme range of fusion relevant parameters. In the present analysis the stationary inertialess assumptions are abandoned. Nevertheless, the classical ideas usually used to obtain inertialess asymptotic solutions are drawn on. The basic inertial equations are reduced to a coupled two-dimensional problem by analytical integration along magnetic field lines. The magnetic field is responsible for a quasi-two-dimensional flow; the non-uniform distribution of the wall conductivity creates a wake-type profile, the MHD effect reducing to a particular forcing and friction. The solution for the two-dimensional variables, the field aligned component of vorticity, the stream function, and the electric potential are obtained by numerical methods. In a flat channel with non-uniform electrical wall conductivity, time-dependent solutions similar to the Kármán vortex street behind bluff bodies are possible. The onset of the vortex motion, i.e. the critical Reynolds number depends strongly on the strength of the magnetic field expressed by the Hartmann number. Stability analyses in viscous hydrodynamic wakes often use the approximation of a unidirectional flow which does not take into account the spatial evolution of the wake. The present problem exhibits a wake-type basic flow, which does not change along the flow path. It represents, therefore, an excellent example to which the simple linear analysis on the basis of Orr-Sommerfeld-type equations applies exactly. Once unstable, the flow first exhibits a regular time periodic vortex pattern which is rearranged further downstream. One can observe an elongation, pairing, or sometimes more complex merging of vortices. All these effects lead to larger flow structures with lower frequencies. The possibility for a creation and maintenance of time-dependent vortex-type flow pattern in MHD flows is demonstrated.

Description: Experiments have been performed, in capillary tubes, on the displacement of a viscous fluid (glycerine) by a less viscous one (a glycerine-water mixture) with which it is miscible in all proportions. A diagnostic measure of the amount of viscous fluid left behind on the tube wall has been found, for both vertical and horizontal tubes, as a function of the Péclet (Pe) and Atwood (At) numbers, as well as a parameter that is a measure of the relative importance of viscous and gravitational effects. The asymptotic value of this diagnostic quantity, for large Pe and an At of unity, has been found to agree with that found in immiscible displacements, while the agreement with the numerical results of Part 2 (Chen & Meiburg 1966), over the whole range of At, is very good. At values of the average Pe greater than 1000 a sharp interface existed so that it was possible to make direct comparisons between the present results and a prior experiment with immiscible fluids, in particular an effective surface tension at the diffusing interface could be evaluated. The effect of gravity on the amount of viscous fluid left on the tube wall has been investigated also, and compared with the results of Part 2. A subsidiary experiment has been performed to measure both the average value of the diffusion coefficient between pure glycerine and several glycerine-water mixtures, in order to be able to calculate a representative Péclet number for each experiment, and the local value as a function of the local concentration of glycerine, in the dilute limit.

Description: Thermoacoustic refrigeration occurs in periodic flow in a duct with heat transfer within the fluid and to the tube. This study considers the periodic limit cycle with large pressure oscillations that is obtained in a tube when prescribed, phase-shifted, periodic velocities at the tube ends, at frequencies lower than acoustic eigenmodes, sweep a length comparable to the tube length. The temperature differences between the two ends are of arbitrary magnitude, heat transfer in the transverse direction within the fluid is assumed to be very effective and the thermal mass of the wall is large. The geometry is two-dimensional, axisymmetric, and conduction is accounted for, not only in the fluid, but also with and within the tube wall. A perturbation solution valid in a local near-isothermal limit determines the equilibrium longitudinal temperature profile that is reached at the periodic regime, the pressure field including longitudinal gradients, and the longitudinal enthalpy flux. Results are presented for tubes open at both ends and also with one end closed. In the latter case, a singularity occurs in the temperature at the closed end, with behaviour identical to Rott's result for acoustic flow with small pressure amplitude. Other new results obtained for tubes open at both ends show that when velocities at both ends are in opposite phase, internal singularities in the temperature profiles may occur.

Description: Following the study in Part 1 of cross-flow and other non-symmetric effects on vortex/wave interactions in boundary layers, the present Part 2 applies the ideas of Part 1 and related works to an incident axisymmetric flow supplemented by a small swirl or azimuthal velocity. This is with a view to possibly increasing understanding of vortex breakdown. The wave components involved are predominantly inviscid Rayleigh-like ones. The presence of the swirl leads to extra features and complications associated mainly with extra logarithmic contributions but for the dominant interactions essentially the same equations as in Part 1 are found. These dominant nonlinear interactions must be based on azimuthal wavenumbers of ± 1 in the case of the Squire jet with swirl. In contrast to Part 1, which consisted mainly of an analysis of the quasi-bounded solutions, a representative set of numerical solutions of the full integro-differential amplitude equations is presented, for realistic axial and swirl velocity profiles. The work points also to the influence of further increases in the incident swirl.

Description: An asymptotic analysis for the long-time unsteady laminar far wake of a bluff body due to a step change in its travelling velocity from U1 to U2 is presented. For U1 ≥ 0 and U2 〉 0, the laminar wake consists of a new wake of volume flux Q2 corresponding to the current velocity U2, an old wake of volume flux Q1 corresponding to the original velocity U1, and a transition zone that connects these two wakes. The transition zone acts as a sink (or a source) of volume flux (Q2 - Q1) and is moving away from the body at speed U2. Streamwise diffusion is negligible in the new and old wakes but a matched asymptotic expansion that retains the streamwise diffusion is required to determine the vorticity transport in the transition zone. A source of volume flux Q2 located near the body needs to be superposed on the unsteady wake to form the global flow field around the body. The asymptotic predictions for the unsteady wake velocity, unsteady wake vorticity, and the global flow field around the body agree well with finite difference solutions for flow over a sphere at finite Reynolds numbers. The long-time unsteady flow structures due to a sudden stop (U2 = 0) and an impulsive reverse (U1 U2 〈 0) of the body are analysed in detail based on the asymptotic solutions for the unsteady wakes and the finite difference solutions. The elucidation of the long-time behaviour of such unsteady flows provides a framework for understanding the long-time particle dynamics at finite Reynolds number.

Description: A theory and simulation code are developed to study non-steady sources as means to control sonic booms of supersonic aircraft. A key result is that the source of sonic boom pressure is not confined to the length of the aircraft but occupies an extensive segment of the flight path. An aircraft in non-steady flight functions as a synthetic aperture antenna, generating complex acoustic waves with no simple relation to instantaneous volume or lift distributions. The theory applies linear acoustics to slender non-steady sources but requires no farfield approximation. The solution for pressure contains a term not seen in Whitham's theory for sonic booms of distant supersonic aircraft. The term describes a pressure field that decays algebraically behind the Mach cone and, in the case of steady flight, integrates to a ground load equal to the weight of the aircraft. The algebraic term is separate from those that describe the sonic boom. Two non-steady source phenomena are evaluated: periodic velocity changes (surge), and periodic longitudinal lift redistribution (slosh). Surge can attenuate a sonic boom and covert it into prolonged weak reverberation, but accelerations needed to produce the phenomenon seem too large for practical use. Slosh may be practical and can alter sonic booms but does not, on average, result in boom attenuation. The conclusion is that active sonic boom abatement is possible in theory but maybe not practical.

Description: The evolution of long waves generated by a pressure disturbance acting on an initially unperturbed free surface in a channel of finite depth is analysed. Both off-critical and transcritical conditions are considered in the context of the fully nonlinear inviscid problem. The solution is achieved by using an accurate boundary integral approach and a time-stepping procedure for the free-surface dynamics. The discussion emphasizes the comparison between the present results and those provided by both the Boussinesq and the related Korteweg-de Vries model. For small amplitudes of the forcing, the predictions of the asymptotic theories are essentially confirmed. However, for finite intensities of the disturbance, several new features significantly affect the physical results. In particular, the interaction among different wave components, neglected in the Korteweg-de Vries approximation, is crucial in determining the evolution of the wave system. A substantial difference is indeed observed between the solutions of the Korteweg-de Vries equation and those of both the fully nonlinear and the Boussinesq model. For increasing dispersion and fixed nonlinearity the agreement between the Boussinesq and fully nonlinear description is lost, indicating a regime where dispersion becomes dominant. Consistently with the long-wave modelling, the transcritical regime is characterized by an unsteady flow and a periodic emission of forward-running waves. However, also in this case, quantitative differences are observed between the three models. For larger amplitudes, wave steepening is almost invariably observed as a precursor of the localized breaking commonly detected in the experiments. The process occurs at the crests of either the trailing or the upstream-emitted wave system for Froude numbers slightly sub- and super-critical respectively.

Description: Experiments and simulations of a travelling wave state of incompressible Newtonian flow in a curved duct of square cross-section are presented. The travelling wave mode develops from the well-documented steady four-cell flow state and is characterized by oscillations of the two Dean vortices near the centre of the outer wall. The oscillations were induced by a carefully positioned pin at 5° from the inlet of the curved section along the symmetry line of the cross-section. It was shown that the travelling wave state is characteristic for curved duct flow and that the pin made it possible to observe the oscillations within the 270° long curved duct. Travelling waves were observed at flow rates above Dn = 170 (Dn = Re/(R/a)l/2, where Re is the Reynolds number, R is the radius of curvature of the duct and a is the duct dimension. The curvature ratio, R/a, is 15.1). If no other disturbances are imposed, the oscillations are the result of the selective amplification of random disturbances in the flow, leading to a broad frequency spectrum. The travelling wave was found to lock in to an imposed periodic disturbance at a selected frequency. The flow structure of the locked state was investigated in detail, using flow visualization and a one-component laser Doppler anemometer to measure streamwise or spanwise velocities. Direct numerical simulations using the package CFDS-FLOW3D are in very good agreement with the experiments and confirm the existence of a fully developed, streamwise-periodic travelling wave state. The inflow region between the two Dean vortices, which transports low-speed fluid away from the outer wall, creates strongly inflectional spanwise profiles of the streamwise velocity. Similarities with twisting vortices in a curved channel and sinuous oscillations of Görtler vortices show that the travelling waves observed here result from a secondary shear instability of these spanwise inflectional profiles.

Description: The anisotropic behaviour of density-gradient fluctuations in stably stratified grid turbulence and the consequences for simplified (isotropic) estimates of scalar dissipation rates χ were experimentally studied in a thermally stratified wind tunnel at moderate Reynolds numbers (Reλ ≃ 20). Strong stable stratifications were attained, with Brunt-Väisälä frequency N as high as 4 rad s-1. The correlation method was used to estimate the mean-square cross-stream and streamwise density gradients. Cross-stream gradients were measured using two cold wires. The mean-square vertical gradients were found to become larger than the streamwise gradients by as much as a factor of 2.2 for the largest dimensionless buoyancy times (Nt = 7). This corresponds to a 40% error in the scalar dissipation estimates based on ∂θ/∂x alone, and assuming the validity of the isotropic relations. Gradient spectral relations show that this buoyancy-induced anisotropy persists at all length scales. Better closure of the scalar variance balance was attained than in previously reported measurements by other researchers. This is attributed to our use of cold-wire temperature sensors having larger length-to-diameter ratio than used in the previous measurements.

Description: A two-dimensional Boussinesq equation, utt-uxx+ 3(u2)xx - uxxxx - uyy = 0, is introduced to describe the propagation of gravity waves on the surface of water, in particular the head-on collision of oblique waves. This equation combines the two-way propagation of the classical Boussinesq equation with the (weak) dependence on a second spatial variable, as occurs in the two-dimensional Korteweg-de Vries (2D KdV) (or KPII) equation. Exact and general solitary-wave, two-soliton and resonant solutions are obtained from the Hirota bilinear form of the equation. The existence of a distributed-soliton solution is investigated, but it is shown that this is not a possibility. However the connection with the classical 2D KdV equation (which does possess such a solution) is explored via a suitable parametric representation of the dispersion relation. A three-soliton solution is also constructed, but this exists only if an auxiliary constraint among the six parameters is satisfied; thus the two-dimensional Boussinesq equation is not one of the class of completely integrable equations, confirming the analysis of Hietarinta (1987). This constraint is automatically satisfied for the classical Boussinesq equation (which is completely integrable). Graphical reproductions of some of the solutions of the two-dimensional Boussinesq equations are also presented.

Description: Individual travelling cavitation bubbles were examined as they interacted with the flow over a two-dimensional hydrofoil. Each bubble was produced from a single nucleus created upstream of the hydrofoil, and the flow near the hydrofoil was visualized using particle imaging velocimetry (PIV). Travelling bubbles were observed to generate a local region of turbulence as they passed close to an unstable laminar boundary layer. By producing a locally turbulent region, the bubbles could temporarily sweep away a portion of attached cavitation at the foil midchord. Also, the bubbles were observed to strongly interact with a turbulent boundary layer, producing local regions of patch cavitation.

Description: The problem of gas motion in a tube closed at one end and driven at the other by an oscillating piston is studied theoretically. When the piston vibrates with a finite amplitude at the first acoustic resonance frequency, periodic shock waves are generated, travelling back and forth in the tube. A perturbation method, based on a small Mach number, M and a global mass conservation condition, is employed to formulate a solution of the problem in the form of two standing waves separated by a jump (shock front). By expanding the equations of motion in a series of a small parameter ε = M1/2, all hydrodynamic properties are predicted with an accuracy to second-order terms, i.e. to ε2. It is found that the first-order solution coincides with the previous theories of Betchov (1958) and Chester (1964), while additional terms predict a non-homogeneous time-averaged pressure along the tube. This prediction compares favourably with experimental results from the literature. The importance of the phenomenon is discussed in relation to different transport processes in resonance tubes.

Description: Two-dimensional mixing driven by an instability mechanism which is concentrated near one of the boundaries is considered, with particular application to Langmuir circulations driven by a wave spectrum. The question of how to define the equivalent of the Rayleigh number is attacked using the energy balance equations and simple truncated models of the instability. Given a particular horizontal wavelength for the disturbance, the strength of the forcing on the cells, and thus the growth rate, is determined by a tradeoff between maximizing the depth-averaged forcing and maximizing the depth of penetration. As a result of this tradeoff, long-wavelength cells grow more slowly, but penetrate more deeply and have a larger equivalent Rayleigh number. At finite amplitude, these long-wavelength cells come to dominate the flow field. The depth of penetration of, and density transport accomplished by, Langmuir cells is considered as a function of the mean stratification and diffusion. An application to oceanic mixed layers is considered assuming the Mellor-Yamada 2 1/2-level turbulence closure model to define the background level of turbulent mixing. For many realistic cases, Langmuir cells are predicted to dominate the vertical transport of momentum and density.

Description: A three-dimensional computer simulation of a concentrated emulsion in shear flow has been developed for low-Reynolds-number finite-capillary-number conditions. Numerical results have been obtained using an efficient boundary integral formulation with periodic boundary conditions and up to twelve drops in each periodically replicated unit cell. Calculations have been performed over a range of capillary numbers where drop deformation is significant up to the value where drop breakup is imminent. Results have been obtained for dispersed-phase volume fractions up to 30% and dispersed- to continuous-phase viscosity ratios in the range of 0 to 5. The results reveal a complex rheology with pronounced shear thinning and large normal stresses that is associated with an anisotropic microstructure that results from the alignment of deformed drops in the flow direction. The viscosity of an emulsion is only a moderately increasing function of the dispersed-phase volume fraction, in contrast to suspensions of rigid particles or undeformed drops. Unlike rigid particles, deformable drops do not form large clusters.

Description: A numerical study of the superharmonic instabilities of short-crested waves on water of finite depth is performed in order to measure their time scales. It is shown that these superharmonic instabilities can be significant - unlike the deep-water case - in parts of the parameter regime. New resonances associated with the standing wave limit are studied closely. These instabilities 'contaminate' most of the parameter space, excluding that near two-dimensional progressive waves; however, they are significant only near the standing wave limit. The main result is that very narrow bands of both short-crested waves 'close' to two-dimensional standing waves, and of well developed short-crested waves, perturbed by superharmonic instabilities, are unstable for depths shallower than approximately a non-dimensional depth d = 1; the study is performed down to depth d = 0.5 beyond which the computations do not converge sufficiently. As a corollary, the present study predicts that these very narrow sub-domains of shortcrested wave fields will not be observable, although most of the short-crested wave fields will be.

Description: In a series of laboratory experiments the growth of double-diffusive salt fingers from an initial configuration of two homogeneous reservoirs with salt in the lower and sugar in the upper layer was investigated. For most of the experiments the stability ratio was between 2.5 and 3, where the latter value is at the upper limit (the ratio of salt to sugar diffusivities) for which fingers can exist. In these experiments long slender fingers are generated at the interface. Essentially all theories or physical bases for models of salt fingers presuppose such a configuration of long fingers. Our measurements show that the length of fingers at high stability ratio increases with time like i'/2, with a coefficient that is consistent with the diffusive spread of the faster diffusing component (salt). When the initial stability ratio is closer to unity, fingers penetrate into the reservoirs very rapidly carrying with them large anomalies of salt and sugar which give rise to convective overturning of the reservoirs. The convection sweeps away the ends of the fingers, and when it is intense enough (as it is when the sugar anomaly is large) it can reduce the finger height to a value less than the width. After this initial phase the finger length grows linearly with time as has been found in previous studies. These results show that salt fingers can evolve in quite different ways depending on the initial stability ratio and must cast doubt on the use of simple similarity arguments to parameterize the heat and salt fluxes produced by fingers.

Description: We show that Kolmogorov's (1941b) inertial-range law for the third-order structure function can be derived from a dynamical equation including pressure terms and mean flow gradient terms. A new inertial-range law, relating the two-point pressure-velocity correlation to the single-point pressure-strain tensor, is also derived. This law shows that the two-point pressure-velocity correlation, just like the third-order structure function, grows linearly with the separation distance in the inertial range. The physical meaning of both this law and Kolmogorov's law is illustrated by a Fourier analysis. An inertial-range law is also derived for the third-order velocity-enstrophy structure function of two-dimensional turbulence. It is suggested that the second-order vorticity structure function of two-dimensional turbulence is constant and scales with εω 2/3 in the enstrophy inertial range, εω being the enstrophy dissipation. Owing to the constancy of this law, it does not imply a Fourier-space inertial-range law, and therefore it is not equivalent to the k-1 law for the enstrophy spectrum, suggested by Kraichnan (1967) and Batchelor (1969).

Description: Simulations of decaying two-dimensional turbulence suggest that the one-point vorticity density has the self-similar form ρω ∼ t f(ωt) implied by Batchelor's (1969) similarity hypothesis, except in the tails. Specifically, similarity holds for |ω| 〈 ωm, while ρω falls off rapidly above. The upper bound of the similarity range, ωm, is also nearly conserved in high-Reynolds-number hyperviscosity simulations and appears to be related to the average amplitude of the most intense vortices (McWilliams 1990), which was an important ingredient in the vortex scaling theory of Carnevale et al. (1991). The universal function f also appears to be hyperbolic, i.e. f(x) ∼ c/2|x|1+qc, for |x| 〉 x*, where qc = 0.4 and x* = 70, which along with the truncated similarity form implies a phase transition in the vorticity moments matrix presented between the self-similar 'background sea' and the coherent vortices. Here cq and c are universal. Low-order moments are therefore consistent with Batchelor's similarity hypothesis whereas high-order moments are similar to those predicted by Carnevale et al. (1991). A self-similar but less well-founded expression for the energy spectrum is also proposed. It is also argued that ωs = x*/t represents 'mean sea-level', i.e. the (average) threshold separating the vortices and the sea, and that there is a spectrum of vortices with amplitudes in the range (ωs, ωm). The total area occupied by vortices is also found to remain constant in time, with losses due to mergers of large-amplitude vortices being balanced by gains due to production of weak vortices. By contrast, the area occupied by vortices above a fixed threshold decays in time as observed by McWilliams (1990).

Description: Motivated by the study of blood flow in the major coronary arteries, which are situated on the outer surface of the pumping heart, we analyse flow of an incompressible Newtonian fluid in a tube whose curvature varies both along the tube and with time. Attention is restricted to the case in which the tube radius is fixed and its centreline moves in a plane. Nevertheless, the governing equations are very complicated, because the natural coordinate system involves acceleration, rotation and deformation of the frame of reference, and their derivation forms a major part of the paper. Then they are applied to two particular, relatively simple examples: a tube of uniform but time-dependent curvature; and a sinuous tube, representing a small-amplitude oscillation about a straight pipe. In the former case the curvature is taken to be small and to vary by a small amount, and the solution is developed as a triple power series in mean curvature ratio δ0, curvature variation ε and Dean number D. In the latter case the Reynolds number is taken to be large and a linearized solution for the perturbation to the flow in the boundary layer at the tube wall is obtained, following Smith (1976a). In each case the solution is taken far enough that the first non-trivial effects of the variable curvature can be determined. Results are presented in terms of the oscillatory wall shear stress distribution and, in the uniform curvature case, the contribution of steady streaming to the mean wall shear stress is calculated. Estimation of the parameters for the human heart indicates that the present results are not directly applicable, but point the way for future work.

Description: The aim of the paper is to present a new transport process which is likely to have great importance for understanding the internal constitution of the stars.In order to set the problem in context, we first give a short presentation of the physical properties of the Sun and stars, described usually under the names Standard Solar Model or Standard Stellar Models (SSM). Next we show that an important shortcoming of SSM is that they do not explain the age dependence of the lithium deficiency of stars of known age: stars of galactic clusters and the Sun. It was suggested a long time ago that the presence of a macroscopic diffusion process in the radiative zone should be assumed, below the surface convective zone of solar-like stars. It is then possible for the lithium present in the convective zone to be carried to the thermonuclear burning level below the convective zone. The first assumption was that differential rotation generates turbulence and therefore that a turbulent diffusion process takes place. However, this model predicts a lithium abundance which is strongly rotation dependent, contrary to the observations. Furthermore, as the diffusion coefficient is large all over the radiative zone, it prevents the possibility of gravitational separation by diffusion and consequently leads to the impossibility of explaining the difference in helium abundance between the surface and the centre of the Sun. The consequence is obviously that we need to take into account another physical process.Stars having a mass M 〈 1.3M[odot ] have a convective zone which begins close to the stellar surface and extends down to a depth which is an appreciable fraction of the stellar radius. In the convective zone, strong stochastic motions carry, at least partially, heat transfer. These motions do not vanish at the lower boundary and generate internal waves into the radiative zone. These random internal waves are at the origin of a diffusion process which can be considered as responsible for the diffusive transport of lithium down to the lithium burning level. This is certainly not the only physical process responsible for lithium deficiency in main sequence stars, but its properties open the way to a completely consistent analysis of lithium deficiency.The model of generation of gravity waves is based on a model of heat transport in the convective zone by diving plumes. The horizontal component of the turbulent motion at the boundary of the convective zone is assumed to generate the horizontal motion of internal waves. The result is a large horizontal component of the diffusion coefficient, which produces in a short time an horizontally uniform chemical composition. It is known that gravity waves, in the absence of any dissipative process, cannot generate vertical mixing. Therefore, the vertical component of the diffusion coefficient is entirely dependent on radiative damping. It decreases quickly in the radiative zone, but is large enough to be responsible for lithium burning.Owing to the radial dependence of velocity amplitude, the diffusion coefficient increases when approaching the stellar centre. However, very close to the centre, nonlinear dissipative and radiative damping of internal waves become large and the diffusion coefficient vanishes at the very centre.

Description: Strongly nonlinear vortex-Tollmien-Schlichting-wave interaction equations are derived for the case where the undisturbed motion represents the developing flow in a circular pipe. The effect upon the equations of moving the wave input position further downstream is investigated and the development of the flow is found to be accelerated by increasing the size of the wave disturbance. Numerical solutions of the three-dimensional interaction equations are presented and indicate that the form of interaction considered here appears to promote the three-dimensionality as the flow develops downstream. It is shown that one of the interactions considered here can develop within an initially two-dimensional Blasius boundary layer.

Description: It is well known that in any conservative system that admits resonant triad interactions, a uniform (test) wavetrain that participates in a single triad is unstable if it has the highest frequency in the triad, and neutrally stable otherwise. We show that this result changes significantly in the presence of coupled triads : with coupling, the test wave can be unstable to a high-frequency perturbation. The coupling sends energy from the (weak) high-frequency source into particular low-frequency waves that grow even though they had zero amplitudes initially. This mechanism thereby selects these low-frequency waves from the spectrum of low-frequency waves available for triad interactions. Moreover, the instability persists in the presence of weak damping, provided the wave amplitudes exceed two thresholds. First, the initial amplitude of the test wavetrain must be large enough for the instability to dominate the damping. Secondly, the (small) initial amplitudes of the high-frequency perturbations must exceed a threshold in order for the low-frequency waves to grow to a prescribed amplitude.

Description: The spectral mechanisms of the differential diffusion of pairs of passive scalars with different molecular diffusivities are studied in stationary isotropic turbulence, using direct numerical simulation data at Taylor-scale Reynolds number up to 160 on 1283 and 2563 grids. Of greatest interest are the roles of nonlinear triadic interactions between different scale ranges of the velocity and scalar fields in the evolution of spectral coherency between the scalars, and the effects of mean scalar gradients. Analysis of single-scalar spectral transfer (extending the results of a previous study) indicates a robust local forward cascade behaviour at high wavenumbers, which is strengthened by both high diffusivity and mean gradients. This cascade is driven primarily by moderately non-local interactions in which two small-scale scalar modes are coupled via a lower-wavenumber velocity mode near the peak of the energy dissipation spectrum. This forward cascade is coherent, tending to increase the coherency between different scalars at high wavenumbers but to decrease it at lower wavenumbers. However, at early times coherency evolution at high wavenumbers is dominated by de-correlating effects due to a different type of non-local triad consisting of two scalar modes with a moderate scale separation and a relatively high-wavenumber velocity mode. Consequently, although the small-scale motions play little role in spectral transfer, they are responsible for the rapid de-correlation observed at early times. At later times both types of competing triadic interactions become important over a wider wavenumber range, with increased relative strength of the coherent cascade, so that the coherency becomes slow-changing. When uniform mean scalar gradients are present, a stationary state develops in the coherency spectrum as a result of a balance between a coherent mean gradient contribution (felt within about 1 eddy-turnover time) and the net contribution from scale interactions. The latter is made less de-correlating because of a strengthened coherent forward cascade, which is in turn caused by uniform mean gradients acting as a primarily low-wavenumber source of scalar fluctuations with the same spectral content as the velocity field.

Description: The effect of wall cooling on hypersonic boundary-layer separation near a compression ramp is considered. Two cases are identified corresponding to the value of the average Mach number M̄ across the upstream boundary layer approaching the compression ramp. The flow is referred to as supercritical for M̄ 〉 1 and subcritical for M̄ 〈 1. The interaction is described by triple-deck theory, and numerical results are given for both cases for various ramp angles and levels of wall cooling. The effect of wall cooling on the absolute instability described recently by Cassel, Ruban & Walker (1995) for an uncooled wall is of particular interest; a stabilizing effect is observed for supercritical boundary layers, but a strong destabilizing influence occurs in the subcritical case. Wall cooling also influences the location and size of the separated region. For supercritical flow, progressive wall cooling reduces the size of the recirculating-flow region, the separation point moves downstream, and upstream influence is diminished. In contrast for the subcritical case downstream influence is reduced with increased cooling. In either situation, a sufficient level of wall cooling eliminates separation altogether for the ramp angles considered. The present numerical results closely confirm the strong wall cooling theory of Kerimbekov, Ruban & Walker (1994).

Description: A number of new experiments have been performed on the rise of air bubbles in clean mixtures of distilled water and pure, reagent grade, glycerine covering a range of the relevant parameter, the Morton number, Mo = gv4ρ3/σ3, of 1013. Here g is the acceleration due to gravity, v the kinematic viscosity, ρ the density and σ the surface tension of the mixture. In these careful measurements several scaling regimes have been found that have not been discussed before in the extensive literature on the subject. The transitions between these regimes have been delineated and attempts made to discuss the dynamical processes that might be important in each of them.

Description: A potential flow model has been formulated for scallop swimming. Under the small-disturbance approximation, the problem of the unsteady flow past the wing-like configuration of a scallop is separated into two linear sub-problems: the steady lifting problem and the unsteady symmetric thickness problem. The latter is associated with the expansion and contraction of the boundary surface of the scallop due to the shell opening and closing. A quasi-two-dimensional analytical solution of the thickness problem was obtained to give the time-dependent fluid forces acting on the outer surfaces of the shells. In addition to the added-mass effect, which has been widely accepted in the hydrodynamics of aquatic locomotion, there are two other mechanisms in the fluid reaction: flow-induced pseudo-elasticity and pseudo-viscosity. The pseudoelasticity provides a force proportional to the gape angle displacement, and will assist shell opening but resist shell closing. The pseudo-viscosity force is proportional to the angular velocity of the gape, and benefits both shell opening and closing. Their roles are discussed through comparison with those of shell inertia, hinge ligament elasticity and hinge damping. At 10°C the hinge damping in the scallop was found to be almost compensated by the flow pseudo-viscosity. The unsteady fluid reaction may have a significant effect on the operation of the dynamic swimming system of scallops.

Description: The linear stability of elliptic pipe flow is considered for finite aspect ratios thereby bridging the gap between the small-aspect-ratio analysis of Davey & Salwen (1994) and the large-aspect-ratio asymptotics of Hocking (1977). The flow is found to become linearly unstable above an aspect ratio of about 10.4 to the spanwise-modulated analogue of the Orr-Sommerfeld mode to which plane Poiseuille flow first loses stability. This disturbance is found to possess a series of intense vortices along its critical layer at lateral stations far removed from the central minor axis. The critical Reynolds number appears to fall from infinity as the aspect ratio increases above 10.4, ultimately approaching Hocking's (1977) asymptotic result at much larger aspect ratios.

Description: In this paper, we have identified a new mechanism which can promote rapid growth of three-dimensional disturbances. The mechanism involves the interaction between a planar mode and an oblique mode, or a pair of oblique modes, which are phase-locked in the sense that they have the same phase speed. This allows a powerful nonlinear interaction to take place within the common critical layer(s). The disturbance is not required to form a subharmonic resonant triad, and hence the mechanism operates under much less restrictive conditions than does subharmonic resonance (although it is somewhat less powerful). We show that the quadratic interaction between the planar mode and the oblique modes drives an exceptionally large forced mode with the difference frequency, which in turn interacts with the planar mode to contribute the dominant nonlinear effect. This interaction can cause the oblique modes to grow super-exponentially provided that their magnitude is sufficiently small. As a result of the super-exponential growth, the oblique mode may soon become strong enough to produce a feedback effect on the planar mode, so that the interactions eventually become fully coupled. This subsequent stage takes slightly different forms depending on whether a single or a pair of oblique modes is present. Both cases are investigated. Particular attention is paid to symmetric plane shear layers, e.g. a planar wake or jet, for which subharmonic resonance of sinuous modes is inactive.

Description: We have investigated the lift force on a small isolated particle which is attached to a flat smooth surface and embedded within the viscous sublayer of the turbulent boundary layer over this surface. We have developed a novel experimental technique with which it is possible to measure both the mean and fluctuating lift force by gluing the particle on top of a silicium cantilever. The deflection of this cantilever is measured with a focused laser beam. The sensitivity of the focus detection system allows us to measure a lift force with an average value around 10-8N and with a standard deviation of approximately 5% of the mean. This means that our device is at least a factor of 100 more sensitive than previous devices and at the same time able to measure the lift forces on smaller particles. Data for the mean lift force (FL+) as a function of the particle radius (a+), where both parameters have been non-dimensionalized with the kinematic viscosity v and the friction velocity u*, are obtained in the range 0.3 〈 a+ 〈 2. The data support the relationship: FL+ = (56.9 + 1.1) (a+)1.87±0.04. Also results on the fluctuating lift force have been obtained. We find that the ratio of the r.m.s. to the mean lift force is approximately 2.8.

Description: Gravity-driven granular flow of slightly frictional particles down an inclined, bumpy chute is studied. A modified kinetic model which includes the frictional energy loss effects is used, and the boundary conditions for a bumpy wall with small friction are derived by ensuring the balance of momentum and energy. At the free surface, the condition of vanishing of the solid volume fraction is used. The mean velocity, the fluctuation kinetic energy and the solid volume fraction profiles are evaluated. It is shown that steady granular gravity flow down a bumpy frictional chute could be achieved at arbitrary inclination angles. The computational results also show that the slip velocity may vary considerably depending on the granular layer height, the surface boundary roughness, the friction coefficient and the inclination angles. The model predictions are compared with the existing experimental and simulation data, and good agreement is observed. In particular, the model can well predicate the features of the variation of solid volume fraction and fluctuation energy profiles for different particle-wall friction coefficients and wall roughnesses.

Description: Experimental investigations in the three-dimensional boundary layer of a swept flat plate with the pressure gradient induced from outside are aimed at enhancing knowledge of the transition process in the presence of pure crossflow instability. The development of disturbances is characterized by the occurrence of both stationary and travelling instability modes, by early nonlinear development and by complex dependence upon the environmental conditions. Experiments under natural conditions of transition showed a good correspondence of the identified modes with those predicted by local linear stability theory. The disturbance growth, however, is generally overpredicted. Controlled excitation of crossflow vortices allowing measurements closer to the linear range of amplification confirmed this result. Nonlinear effects such as interaction between stationary disturbances and base flow and between travelling and stationary modes have already been observed when the naturally excited instabilities become of measurable size. The most striking feature of the disturbance development is the complex dependence on initial conditions. Experiments under systematically varied environments showed that surface roughness represents the key parameter responsible for the initiation of stationary crossflow vortices. In contrast to two-dimensional boundary layers, free-stream turbulence influences the transition process indirectly. Only for turbulence levels Tu 〉 0.2% and smooth surfaces do the travelling instability waves dominate. The location of the final breakdown of laminar flow is clearly determined by the saturation amplitude of crossflow vortices. The receptivity to sound, two-dimensional surface roughness and non-uniformities of the test-section mean flow was found to be very weak.