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.

Testing the influence of small crystals on ice size spectra using Doppler lidar observations

Analysis of the vertical velocity of ice crystals observed with a 1.5micron Doppler lidar from a continuous sample of stratiform ice clouds over 17 months show that the distribution of Doppler velocity varies strongly with temperature, with mean velocities of 0.2m/s at -40C, increasing to 0.6m/s at -10C due to particle growth and broadening of the size spectrum. We examine the likely influence of crystals smaller than 60microns by forward modelling their effect on the area-weighted fall speed, and comparing the results to the lidar observations. The comparison strongly suggests that the concentration of small crystals in most clouds is much lower than measured in-situ by some cloud droplet probes. We argue that the discrepancy is likely due to shattering of large crystals on the probe inlet, and that numerous small particles should not be included in numerical weather and climate model parameterizations.

Industrial radioactive barite scale: suppression of radium uptake by introduction of competing ions

Incorporation of radioactive isotopes during the formation of barite mineral scale is a widespread phenomenon occurring within the oil, mining and process industries. In a series of experiments radioactive barite/celestite solid solutions (SSBarite-Celcstite) have been synthesized under controlled conditions by the counter diffusion of Ra-226, Ba2+, Sr24+ and SO42- ions through a porous medium (silica gel), to investigate inhibiting effects in Ra uptake associated with the introduction of a competing ion (Sr2+). From characterization studies, the particle size and the morphology of the crystals appear to be related to the initial [Sr]/[Ba] molar ratio of the starting solution. Typically, systems richer in Sr produce smaller sized crystals and clusters characterized by a lower degree of order. The activity introduced to the system is mainly incorporated in the crystals generated from the barite/celestite solid solution as suggested by the activity profiles of the hydrogel columns analysed by gamma-spectrometry. There is a relationship between the initial [Sr]/[Ba] molar ratio of the starting solution and the activity exhibited by the synthesized crystals. An effective inhibition of the Ra-226 uptake during formation of the crystals (SSBarite-Celestite) was obtained through the introduction of a competing ion (Sr2+): the higher the initial [Sr]/[Ba] molar ratio of the starting solution, the lower the intensity of the activity peak in the crystals. (C) 2003 Published by Elsevier Ltd.

Radar scattering by aggregate snowflakes

The radar scattering properties of realistic aggregate snowflakes have been calculated using the Rayleigh-Gans theory. We find that the effect of the snowflake geometry on the scattering may be described in terms of a single universal function, which depends only on the overall shape of the aggregate and not the geometry or size of the pristine ice crystals which compose the flake. This function is well approximated by a simple analytic expression at small sizes; for larger snowflakes we fit a curve to Our numerical data. We then demonstrate how this allows a characteristic snowflake radius to be derived from dual wavelength radar measurements without knowledge of the pristine crystal size or habit, while at the same time showing that this detail is crucial to using such data to estimate ice water content. We also show that the 'effective radius'. characterizing the ratio of particle volume to projected area, cannot be inferred from dual wavelength radar data for aggregates. Finally, we consider the errors involved in approximating snowflakes by 'air-ice spheres', and show that for small enough aggregates the predicted dual wavelength ratio typically agrees to within a few percent, provided some care is taken in choosing the radius of the sphere and the dielectric constant of the air-ice mixture; at larger sizes the radar becomes more sensitive to particle shape, and the errors associated with the sphere model are found to increase accordingly.

Dynamic scaling of bred vectors in spatially extended chaotic systems

We unfold a profound relationship between the dynamics of finite-size perturbations in spatially extended chaotic systems and the universality class of Kardar-Parisi-Zhang (KPZ). We show how this relationship can be exploited to obtain a complete theoretical description of the bred vectors dynamics. The existence of characteristic length/time scales, the spatial extent of spatial correlations and how to time it, and the role of the breeding amplitude are all analyzed in the light of our theory. Implications to weather forecasting based on ensembles of initial conditions are also discussed.

A rotating disc voltammetry study of the 1,8-dihydroxyanthraquinone mediated reduction of colloidal indigo

Colloidal indigo is reduced to an aqueous solution of leuco-indigo in a mediated two-electron process converting the water-insoluble dye into the water-soluble leuco form. The colloidal dye does not interact directly with the electrode surface, and to employ an electrochemical process for this reduction, the redox mediator 1,8-dihydroxyanthraquinone (1,8-DHAQ) is used to transfer electrons from the electrode to the dye. The mediated reduction process is investigated at a (500-kHz ultrasound-assisted) rotating disc electrode, and the quantitative analysis of voltammetric data is attempted employing the Digisim numerical simulation software package. At the most effective temperature, 353 K, the diffusion coefficient for 1,8-DHAQ is (0.84 +/- 0.08)x10(-9) m(2) s(-1), and it is shown that an apparently kinetically controlled reaction between the reduced form of the mediator and the colloidal indigo occurs within the diffusion layer at the electrode surface. The apparent bimolecular rate constant k (app)=3 mol m(-3) s(-1) for the rate law d[leuco-indigo]/dt = k(app) x [mediator] x [indigo] is determined and attributed to a mediator diffusion controlled dissolution of the colloid particles. The average particle size and the number of molecules per particles are estimated from the apparent bimolecular rate constant and confirmed by scanning electron microscopy.

The mechanism of bacterial indigo reduction

The reduction of water-insoluble indigo by the recently isolated moderate thermophile, Clostridium isatidis, has been studied with the aim of developing a sustainable technology for industrial indigo reduction. The ability to reduce indigo was not shared with C. aurantibutyricum, C. celatum and C. papyrosolvens, but C. papyrosolvens could reduce indigo carmine (5,5-indigosulfonic acid), a soluble indigo derivative. The supernatant from cultures of C. isatidis, but not from cultures of the other bacteria tested, decreased indigo particle size to one-tenth diameter. Addition of madder powder, anthraquinone-2,6-disulfonic acid, and humic acid all stimulated indigo reduction by C. isatidis. Redox potentials of cultures of C. isatidis were about 100 mV more negative than those of C. aurantibutyricum, C. celatum and C. papyrosolvens, and reached –600 mV versus the SCE in the presence of indigo, but potentials were not consistently affected by the addition of the quinone compounds, which probably act by modifying the surface of the bacteria or indigo particles. It is concluded that C. isatidis can reduce indigo because (1) it produces an extracellular factor that decreases indigo particle size, and (2) it generates a sufficiently reducing potential.

Particulate matter mass concentration (PM10) under different ventilation methods in classrooms

Recently, studies have shown that the classroom environment is very important for students' health and performance. Thus, the evaluation of indoor air quality (IAQ) in a classroom is necessary to ensure students' well-being. In this paper the emphasis is on airborne concentration of particulate matter (PM) in adult education rooms. The mass concentration of PM10 particulates was measured in two classrooms under different ventilation methods in the University of Reading, UK, during the winter period of 2008. In another study the measurement of the concentration of particles was accompanied with measurements of CO2 concentration in these classrooms but this study is the subject of another publication. The ambient PM10, temperature, relative humidity, wind speed and direction, and rainfall events were monitored as well. In general, this study showed that outdoor particle concentrations and outdoor meteorological parameters were identified as significant factors influencing indoor particle concentration levels. Ventilation methods showed significant effects on air change rate and on indoor/outdoor (I/O) concentration ratios. Higher levels of indoor particulates were seen during occupancy periods. I/O ratios were significantly higher when classrooms were occupied than when they were unoccupied, indicating the effect of both people presence and outdoor particle concentration levels. The concentrations of PM10 indoors and outdoors did not meet the requirements of WHO standards for PM10 annual average.

Efficient on the fly calculation of time correlation functions in computer simulations

Time correlation functions yield profound information about the dynamics of a physical system and hence are frequently calculated in computer simulations. For systems whose dynamics span a wide range of time, currently used methods require significant computer time and memory. In this paper, we discuss the multiple-tau correlator method for the efficient calculation of accurate time correlation functions on the fly during computer simulations. The multiple-tau correlator is efficacious in terms of computational requirements and can be tuned to the desired level of accuracy. Further, we derive estimates for the error arising from the use of the multiple-tau correlator and extend it for use in the calculation of mean-square particle displacements and dynamic structure factors. The method described here, in hardware implementation, is routinely used in light scattering experiments but has not yet found widespread use in computer simulations.

Mechanisms of burst release from pH-responsive polymeric microparticles.

Microencapsulation of drugs into preformed polymers is commonly achieved through solvent evaporation techniques or spray drying. We compared these encapsulation methods in terms of controlled drug release properties of the prepared microparticles and investigated the underlying mechanisms responsible for the “burst release” effect. Using two different pH-responsive polymers with a dissolution threshold of pH 6 (Eudragit L100 and AQOAT AS-MG), hydrocortisone, a model hydrophobic drug, was incorporated into microparticles below and above its solubility within the polymer matrix. Although, spray drying is an attractive approach due to rapid particle production and relatively low solvent waste, the oil-in-oil microencapsulation method is superior in terms of controlled drug release properties from the microparticles. Slow solvent evaporation during the oil-in-oil emulsification process allows adequate time for drug and polymer redistribution in the microparticles and reduces uncontrolled drug burst release. Electron microscopy showed that this slower manufacturing procedure generated non-porous particles whereas thermal analysis and X-ray diffractometry showed that drug loading above the solubility limit of the drug in the polymer generated excess crystalline drug on the surface of the particles. Raman spectral mapping illustrated that drug was homogeneously distributed as a solid solution in the particles when loaded below saturation in the polymer with consequently minimal burst release.

Efficient non-linear data assimilation in geophysical fluid dynamics

New ways of combining observations with numerical models are discussed in which the size of the state space can be very large, and the model can be highly nonlinear. Also the observations of the system can be related to the model variables in highly nonlinear ways, making this data-assimilation (or inverse) problem highly nonlinear. First we discuss the connection between data assimilation and inverse problems, including regularization. We explore the choice of proposal density in a Particle Filter and show how the ’curse of dimensionality’ might be beaten. In the standard Particle Filter ensembles of model runs are propagated forward in time until observations are encountered, rendering it a pure Monte-Carlo method. In large-dimensional systems this is very inefficient and very large numbers of model runs are needed to solve the data-assimilation problem realistically. In our approach we steer all model runs towards the observations resulting in a much more efficient method. By further ’ensuring almost equal weight’ we avoid performing model runs that are useless in the end. Results are shown for the 40 and 1000 dimensional Lorenz 1995 model.

Radar scattering from ice aggregates using the horizontally aligned oblate spheroid approximation

The assumed relationship between ice particle mass and size is profoundly important in radar retrievals of ice clouds, but, for millimeter-wave radars, shape and preferred orientation are important as well. In this paper the authors first examine the consequences of the fact that the widely used ‘‘Brown and Francis’’ mass–size relationship has often been applied to maximumparticle dimension observed by aircraftDmax rather than to the mean of the particle dimensions in two orthogonal directions Dmean, which was originally used by Brown and Francis. Analysis of particle images reveals that Dmax ’ 1.25Dmean, and therefore, for clouds for which this mass–size relationship holds, the consequences are overestimates of ice water content by around 53% and of Rayleigh-scattering radar reflectivity factor by 3.7 dB. Simultaneous radar and aircraft measurements demonstrate that much better agreement in reflectivity factor is provided by using this mass–size relationship with Dmean. The authors then examine the importance of particle shape and fall orientation for millimeter-wave radars. Simultaneous radar measurements and aircraft calculations of differential reflectivity and dual-wavelength ratio are presented to demonstrate that ice particles may usually be treated as horizontally aligned oblate spheroids with an axial ratio of 0.6, consistent with them being aggregates. An accurate formula is presented for the backscatter cross section apparent to a vertically pointing millimeter-wave radar on the basis of a modified version of Rayleigh–Gans theory. It is then shown that the consequence of treating ice particles as Mie-scattering spheres is to substantially underestimate millimeter-wave reflectivity factor when millimeter-sized particles are present, which can lead to retrieved ice water content being overestimated by a factor of 4.h

Validation of ACE-FTS satellite data in the upper troposphere/lower stratosphere (UTLS) using non-coincident measurements

CO, O3, and H2O data in the upper troposphere/lower stratosphere (UTLS) measured by the Atmospheric Chemistry Experiment Fourier Transform Spectrometer(ACE-FTS) on Canada’s SCISAT-1 satellite are validated using aircraft and ozonesonde measurements. In the UTLS, validation of chemical trace gas measurements is a challenging task due to small-scale variability in the tracer fields, strong gradients of the tracers across the tropopause, and scarcity of measurements suitable for validation purposes. Validation based on coincidences therefore suffers from geophysical noise. Two alternative methods for the validation of satellite data are introduced, which avoid the usual need for coincident measurements: tracer-tracer correlations, and vertical tracer profiles relative to tropopause height. Both are increasingly being used for model validation as they strongly suppress geophysical variability and thereby provide an “instantaneous climatology”. This allows comparison of measurements between non-coincident data sets which yields information about the precision and a statistically meaningful error-assessment of the ACE-FTS satellite data in the UTLS. By defining a trade-off factor, we show that the measurement errors can be reduced by including more measurements obtained over a wider longitude range into the comparison, despite the increased geophysical variability. Applying the methods then yields the following upper bounds to the relative differences in the mean found between the ACE-FTS and SPURT aircraft measurements in the upper troposphere (UT) and lower stratosphere (LS), respectively: for CO ±9% and ±12%, for H2O ±30% and ±18%, and for O3 ±25% and ±19%. The relative differences for O3 can be narrowed down by using a larger dataset obtained from ozonesondes, yielding a high bias in the ACEFTS measurements of 18% in the UT and relative differences of ±8% for measurements in the LS. When taking into account the smearing effect of the vertically limited spacing between measurements of the ACE-FTS instrument, the relative differences decrease by 5–15% around the tropopause, suggesting a vertical resolution of the ACE-FTS in the UTLS of around 1 km. The ACE-FTS hence offers unprecedented precision and vertical resolution for a satellite instrument, which will allow a new global perspective on UTLS tracer distributions.

Intergenerational wealth mobility in rural Bangladesh

Unique residential history data with retrospective information on parental assets are used to study household wealth mobility in 141 villages in rural Bangladesh. Regression estimates of father–son correlations and analyses of intergenerational transition matrices show substantial persistence in wealth even when we correct for measurement errors in parental wealth. We do not find wealth mobility to be higher between periods of a person's life than between generations. We find that the process of household division plays an important role: sons who splinter off from the father's household experience greater (albeit downward) mobility in wealth. Despite significant occupational mobility across generations, its contribution to wealth mobility, net of human capital attainment of individuals, appears insignificant. Low wealth mobility in our data is primarily explained by intergenerational persistence in educational attainment.

On the statistical modeling of persistence in total ozone anomalies

Geophysical time series sometimes exhibit serial correlations that are stronger than can be captured by the commonly used first‐order autoregressive model. In this study we demonstrate that a power law statistical model serves as a useful upper bound for the persistence of total ozone anomalies on monthly to interannual timescales. Such a model is usually characterized by the Hurst exponent. We show that the estimation of the Hurst exponent in time series of total ozone is sensitive to various choices made in the statistical analysis, especially whether and how the deterministic (including periodic) signals are filtered from the time series, and the frequency range over which the estimation is made. In particular, care must be taken to ensure that the estimate of the Hurst exponent accurately represents the low‐frequency limit of the spectrum, which is the part that is relevant to long‐term correlations and the uncertainty of estimated trends. Otherwise, spurious results can be obtained. Based on this analysis, and using an updated equivalent effective stratospheric chlorine (EESC) function, we predict that an increase in total ozone attributable to EESC should be detectable at the 95% confidence level by 2015 at the latest in southern midlatitudes, and by 2020–2025 at the latest over 30°–45°N, with the time to detection increasing rapidly with latitude north of this range.