Structure of turbulent flow over arrays of cubical roughness

The structure of turbulent flow over large roughness consisting of regular arrays of cubical obstacles is investigated numerically under constant pressure gradient conditions. Results are analysed in terms of first- and second-order statistics, by visualization of instantaneous flow fields and by conditional averaging. The accuracy of the simulations is established by detailed comparisons of first- and second-order statistics with wind-tunnel measurements. Coherent structures in the log region are investigated. Structure angles are computed from two-point correlations, and quadrant analysis is performed to determine the relative importance of Q2 and Q4 events (ejections and sweeps) as a function of height above the roughness. Flow visualization shows the existence of low-momentum regions (LMRs) as well as vortical structures throughout the log layer. Filtering techniques are used to reveal instantaneous examples of the association of the vortices with the LMRs, and linear stochastic estimation and conditional averaging are employed to deduce their statistical properties. The conditional averaging results reveal the presence of LMRs and regions of Q2 and Q4 events that appear to be associated with hairpin-like vortices, but a quantitative correspondence between the sizes of the vortices and those of the LMRs is difficult to establish; a simple estimate of the ratio of the vortex width to the LMR width gives a value that is several times larger than the corresponding ratio over smooth walls. The shape and inclination of the vortices and their spatial organization are compared to recent findings over smooth walls. Characteristic length scales are shown to scale linearly with height in the log region. Whilst there are striking qualitative similarities with smooth walls, there are also important differences in detail regarding: (i) structure angles and sizes and their dependence on distance from the rough surface; (ii) the flow structure close to the roughness; (iii) the roles of inflows into and outflows from cavities within the roughness; (iv) larger vortices on the rough wall compared to the smooth wall; (v) the effect of the different generation mechanism at the wall in setting the scales of structures.

Calculating the ocean's mean dynamic topography from a mean sea surface and a geoid

In principle the global mean geostrophic surface circulation of the ocean can be diagnosed by subtracting a geoid from a mean sea surface (MSS). However, because the resulting mean dynamic topography (MDT) is approximately two orders of magnitude smaller than either of the constituent surfaces, and because the geoid is most naturally expressed as a spectral model while the MSS is a gridded product, in practice complications arise. Two algorithms for combining MSS and satellite-derived geoid data to determine the ocean’s mean dynamic topography (MDT) are considered in this paper: a pointwise approach, whereby the gridded geoid height field is subtracted from the gridded MSS; and a spectral approach, whereby the spherical harmonic coefficients of the geoid are subtracted from an equivalent set of coefficients representing the MSS, from which the gridded MDT is then obtained. The essential difference is that with the latter approach the MSS is truncated, a form of filtering, just as with the geoid. This ensures that errors of omission resulting from the truncation of the geoid, which are small in comparison to the geoid but large in comparison to the MDT, are matched, and therefore negated, by similar errors of omission in the MSS. The MDTs produced by both methods require additional filtering. However, the spectral MDT requires less filtering to remove noise, and therefore it retains more oceanographic information than its pointwise equivalent. The spectral method also results in a more realistic MDT at coastlines. 1. Introduction An important challenge in oceanography is the accurate determination of the ocean’s time-mean dynamic topography (MDT). If this can be achieved with sufficient accuracy for combination with the timedependent component of the dynamic topography, obtainable from altimetric data, then the resulting sum (i.e., the absolute dynamic topography) will give an accurate picture of surface geostrophic currents and ocean transports.

The fall speeds of sub-100 mu m ice crystals

Estimates for the sedimentation rate of realistic ice crystals at sizes smaller than 100 µm are presented. These calculations, which exploit new results for the capacitance of ice crystals, are compared with laboratory studies and found to be in good agreement. The results highlight a weakness in contemporary ice particle fall speed parametrizations for very small crystals, which can lead to sedimentation rates being overestimated by a factor of two. The theoretical approach applied here may also be useful for calculating the sedimentation rate and mobility of non-spherical aerosol particles.

Language processing with dynamic fields

We construct a mapping from complex recursive linguistic data structures to spherical wave functions using Smolensky's filler/role bindings and tensor product representations. Syntactic language processing is then described by the transient evolution of these spherical patterns whose amplitudes are governed by nonlinear order parameter equations. Implications of the model in terms of brain wave dynamics are indicated.

The acoustic signature of fluid flow in complex porous media

Effective medium approximations for the frequency-dependent and complex-valued effective stiffness tensors of cracked/ porous rocks with multiple solid constituents are developed on the basis of the T-matrix approach (based on integral equation methods for quasi-static composites), the elastic - viscoelastic correspondence principle, and a unified treatment of the local and global flow mechanisms, which is consistent with the principle of fluid mass conservation. The main advantage of using the T-matrix approach, rather than the first-order approach of Eshelby or the second-order approach of Hudson, is that it produces physically plausible results even when the volume concentrations of inclusions or cavities are no longer small. The new formulae, which operates with an arbitrary homogeneous (anisotropic) reference medium and contains terms of all order in the volume concentrations of solid particles and communicating cavities, take explicitly account of inclusion shape and spatial distribution independently. We show analytically that an expansion of the T-matrix formulae to first order in the volume concentration of cavities (in agreement with the dilute estimate of Eshelby) has the correct dependence on the properties of the saturating fluid, in the sense that it is consistent with the Brown-Korringa relation, when the frequency is sufficiently low. We present numerical results for the (anisotropic) effective viscoelastic properties of a cracked permeable medium with finite storage porosity, indicating that the complete T-matrix formulae (including the higher-order terms) are generally consistent with the Brown-Korringa relation, at least if we assume the spatial distribution of cavities to be the same for all cavity pairs. We have found an efficient way to treat statistical correlations in the shapes and orientations of the communicating cavities, and also obtained a reasonable match between theoretical predictions (based on a dual porosity model for quartz-clay mixtures, involving relatively flat clay-related pores and more rounded quartz-related pores) and laboratory results for the ultrasonic velocity and attenuation spectra of a suite of typical reservoir rocks. (C) 2003 Elsevier B.V. All rights reserved.

The role of gypsum in the reactions of phosphate with soils

The aim of this study was to examine the mechanisms by which gypsum increases the sorption of fertilizer-P in soils of and and semi-arid regions. Either gypsum or soil (Usher from the UK; pH 7.8, 7% organic matter, 21% CaCO3: Yasouj from Iran; pH 8.2, 1.4% OM, 18% CaCO3: Ghanimeh from Saudi Arabia; pH 7.8, 1% OM, 26% CaCO3, 13% gypsum) was shaken for 24 It with KH2PO4 solutions in 10 mM CaCl2. With gypsum, grinding increased sorption by a factor of about 3, and increase in pH from 5.6 to 7.5 greatly increased sorption. Scanning electron micrographs (SEM) and EDX quantitative analysis showed that small crystals of gypsum disappeared and roughly spherical particles of dicalcium phosphate (DCPD) were formed. Analysis of equilibrium Solutions showed, using GEOCHEM, that octa-calcium phosphate (OCP) coated the DCPD. For the soils, sorption was in the order Ghanimeh > Yasouj > Usher. Removal of gypsum from Ghanimeh reduced sorption, with precipitated gypsum having a greater effect than gypsum mixed physically with the soil. Addition to Usher had no effect. SEM and EDX could not be used in the soil matrix, but solubility analysis again showed that solutions were close to equilibrium with OCP. Usher was unresponsive to added gypsum, presumably because of its small sorption capacity and high organic matter content. In Ghanimeh and Yasouj soils, gypsum increased sorption by being a source of readily available Ca2+ (c) 2005 Elsevier B.V. All rights reserved.

Campanian crustacean burrow system from Israel with brood and nursery chambers representing communal organization

Phosphorite-filled crustacean burrows associated with a Campanian-age omission surface in the north-western Negev are described. The phosphatic burrow casts weather out displaying scratches (bioglyphs) and two types of local swellings (chambers), which are flattened normal to the course of the burrow. The more abundant chamber type is a flattened spheroid (diameter 45-50 mm) or a flattened, highly prolate ellipsoid of larger dimensions, with bioglyphs. The other type is a flattened spheroid (diameter 45 mm), gently rounded on the upper side and flat on the base. Rings of elevations on the cast (representing moats) form interconnected circlets, each capped by about eight rounded hemispherical tubercles (4 x 4 mm) (pits on original), the whole forming a discrete network. The first type of chamber may have hosted the young (nursery chamber) and/or stored food. The second type of cast replicates a chamber with a pitted floor, which may have formed a brood chamber for 60-70 spherical eggs, each about 3 mm in diameter. Brood chambers in crustacean burrow systems were previously suspected, but only at burrow terminations. The interpreted K-type breeding strategy, brood care and associated functions require a high degree of social organization, none of which has been observed in extant crustaceans, but all occur within social insects.

Solvent Effects on the Formation of Nanoparticles and Multilayered Coatings Based on Hydrogen-Bonded Interpolymer Complexes of Poly(acrylic acid) with Homo- and Copolymers of N-Vinyl Pyrrolidone

The formation of hydrogen-bonded interpolymer complexes between poly(acrylic acid) and poly(N-vinyl pyrrolidone) as well as amphiphilic copolymers of N-vinyl pyrrolidone with vinyl propyl ether has been studied in aqueous and organic solutions. It was demonstrated that introduction of vinyl propyl ether units into the macromolecules of the nonionic polymer enhances their ability to form complexes in aqueous solutions due to more significant contribution of hydrophobic effects. The complexation was found to be a multistage process that involves the formation of primary polycomplex particles, which further aggregate to form spherical nanoparticles. Depending on the environmental factors (pH, solvent nature), these nanoparticles may either form stable colloidal solutions or undergo further aggregation, resulting in precipitation of interpolymer complexes. In organic solvents, the intensity of complex formation increases in the following order: methanol < ethanol < isopropanol < dioxane. The multilayered coatings were developed using layer-by-layer deposition of interpolymer complexes on glass surfaces. It was demonstrated that the solvent nature affects the efficiency of coating deposition.

Droplet shape of an anisotropic liquid

We investigate how a droplet of a complex liquid is modified by its internal nanoscale structure. As the liquid passes from an isotropic disordered state to an anisotropic layered morphology, the droplet shape switches from a smooth spherical cap to a terraced hyperbolic profile, which can be modeled as a stack of thin concentric circular disks with a repulsion between adjacent disk edges. Our ability to resolve the detailed shape of these defect-free droplets offers a unique opportunity to explore the underlying physics.

Converting the nanodomains of a diblock-copolymer thin film from spheres to cylinders with an external electric field

We investigate the ability of an applied electric field to convert the morphology of a diblock-copolymer thin film from a monolayer of spherical domains embedded in the matrix to cylindrical domains that penetrate through the matrix. As expected, the applied field increases the relative stability of cylindrical domains, while simultaneously reducing the energy barrier that impedes the transition to cylinders. The effectiveness of the field is enhanced by a large dielectric contrast between the two block-copolymer components, particularly when the low-dielectric contrast component forms the matrix. Furthermore, the energy barrier is minimized by selecting sphere-forming diblock copolymers that are as compositionally symmetric as possible. Our calculations, which are the most quantitatively reliable to date, are performed using a numerically precise spectral algorithm based on self-consistent-field theory supplemented with an exact treatment for linear dielectric materials.

Consistent approximate models of the global atmosphere: shallow, deep, hydrostatic, quasi-hydrostatic and non-hydrostatic

We study global atmosphere models that are at least as accurate as the hydrostatic primitive equations (HPEs), reviewing known results and reporting some new ones. The HPEs make spherical geopotential and shallow atmosphere approximations in addition to the hydrostatic approximation. As is well known, a consistent application of the shallow atmosphere approximation requires omission of those Coriolis terms that vary as the cosine of latitude and of certain other terms in the components of the momentum equation. An approximate model is here regarded as consistent if it formally preserves conservation principles for axial angular momentum, energy and potential vorticity, and (following R. Müller) if its momentum component equations have Lagrange's form. Within these criteria, four consistent approximate global models, including the HPEs themselves, are identified in a height-coordinate framework. The four models, each of which includes the spherical geopotential approximation, correspond to whether the shallow atmosphere and hydrostatic (or quasi-hydrostatic) approximations are individually made or not made. Restrictions on representing the spatial variation of apparent gravity occur. Solution methods and the situation in a pressure-coordinate framework are discussed. © Crown copyright 2005.

The structure of crystallisable copolymers of L-lactide, epsilon-caprolactone and glycolide

We apply a new X-ray scattering approach to the study of melt-spun filaments of tri-block and random terpolymers prepared from lactide, caprolactone and glycolide. Both terpolymers contain random sequences, in both cases the overall fraction of lactide units is similar to 0.7 and C-13 and H-1 NMR shows the lactide sequence length to be similar to 9-10. A novel representation of the X-ray fibre pattern as series of spherical harmonic functions considerably facilitates the comparison of the scattering from the minority crystalline phase with hot drawn fibres prepared from the poly(L-lactide) homopolymer. Although the fibres exhibit rather disordered structures we show that the crystal structure is equivalent to that displayed by poly(L-lactide) for both the block and random terpolymers. There are variations in the development of a two-phase structure which reflect the differences in the chain architectures. There is evidence that the random terpolymer includes non-lactide units in to the crystal interfaces to achieve a well defined two-phase structure. (c) 2005 Published by Elsevier Ltd.

Rotational spectra induced by vibrations

Mizushima and Venkateswarlu showed in 1953 that certain molecules have the property that excited vibrational states may possess rotational spectra even when the rotational spectrum of the ground vibrational state is forbidden by symmetry. We call such a spectrum a vibrationally induced rotational spectrum, and have made a systematic examination of the point groups which permit such behaviour. We also give formulae for the approximate line frequencies and intensities in these spectra, and discuss some of the problems involved in observing them. The spectra can only arise from degenerate vibrational states, and are of three possible types: i) symmetric top perpendicular spectra, shown by molecules belonging to the point groups Dnh, Dn and Cnh, where n is odd; (ii) symmetric top parallel spectra, shown by molecules belonging to Dnd and S2n, where n is even; and (iii) spherical top spectra, shown by molecules belonging to T or Td. Excited vibrational states of polar molecules of point groups Cnv or Cn, where n is odd, may also possess vibrationally induced perpendicular components of type (i), in addition to their ordinary parallel spectra. In addition to the above limitations on the point groups there are, in general, limitations on the symmetry species of the degenerate vibrational states.

Wave functions for the methane molecule

If the potential field due to the nuclei in the methane molecule is expanded in terms of a set of spherical harmonics about the carbon nucleus, only the terms involving s, f, and higher harmonic functions differ from zero in the equilibrium configuration. Wave functions have been calculated for the equilibrium configuration, first including only the spherically symmetric s term in the potential, and secondly including both the s and the f terms. In the first calculation the complete Hartree-Fock S.C.F. wave functions were determined; in the second calculation a variation method was used to determine the best form of the wave function involving f harmonics. The resulting wave functions and electron density functions are presented and discussed

Vibrational intensities in methane

The absorption intensities of the two infra-red active vibrations in methane have been obtained from a perturbation calculation on the equilibrium wave functions derived in the preceding paper. The perturbation field is the change in the potential field due to the nuclei which results from moving the nuclei in the vibrational coordinate concerned, and a simplified form of second order perturbation theory, developed by Pople and Schofield, is used for the calculation. The main approximation involved is the neglect of f and higher harmonics in the spherical harmonic expansion of the nuclear field. The resulting dipole moment derivatives are approximately three times larger than the experimental values, but they show qualitative features and sign relationships which are significant.