A finite difference scheme for inviscid flows with non-eequilibrium chemistry and internal energy

A finite difference scheme is presented for the inviscid terms of the equations of compressible fluid dynamics with general non-equilibrium chemistry and internal energy.

Results of an aqueous source term model for a radiological risk assessment of the Drigg LLW Site, U.K.

A radionuclide source term model has been developed which simulates the biogeochemical evolution of the Drigg low level waste (LLW) disposal site. The DRINK (DRIgg Near field Kinetic) model provides data regarding radionuclide concentrations in groundwater over a period of 100,000 years, which are used as input to assessment calculations for a groundwater pathway. The DRINK model also provides input to human intrusion and gaseous assessment calculations through simulation of the solid radionuclide inventory. These calculations are being used to support the Drigg post closure safety case. The DRINK model considers the coupled interaction of the effects of fluid flow, microbiology, corrosion, chemical reaction, sorption and radioactive decay. It represents the first direct use of a mechanistic reaction-transport model in risk assessment calculations.

State estimation using model order reduction for unstable systems

The problem of state estimation occurs in many applications of fluid flow. For example, to produce a reliable weather forecast it is essential to find the best possible estimate of the true state of the atmosphere. To find this best estimate a nonlinear least squares problem has to be solved subject to dynamical system constraints. Usually this is solved iteratively by an approximate Gauss–Newton method where the underlying discrete linear system is in general unstable. In this paper we propose a new method for deriving low order approximations to the problem based on a recently developed model reduction method for unstable systems. To illustrate the theoretical results, numerical experiments are performed using a two-dimensional Eady model – a simple model of baroclinic instability, which is the dominant mechanism for the growth of storms at mid-latitudes. It is a suitable test model to show the benefit that may be obtained by using model reduction techniques to approximate unstable systems within the state estimation problem.

On the distortion of turbulence by a progressive surface wave

A rapid-distortion model is developed to investigate the interaction of weak turbulence with a monochromatic irrotational surface water wave. The model is applicable when the orbital velocity of the wave is larger than the turbulence intensity, and when the slope of the wave is sufficiently high that the straining of the turbulence by the wave dominates over the straining of the turbulence by itself. The turbulence suffers two distortions. Firstly, vorticity in the turbulence is modulated by the wave orbital motions, which leads to the streamwise Reynolds stress attaining maxima at the wave crests and minima at the wave troughs; the Reynolds stress normal to the free surface develops minima at the wave crests and maxima at the troughs. Secondly, over several wave cycles the Stokes drift associated with the wave tilts vertical vorticity into the horizontal direction, subsequently stretching it into elongated streamwise vortices, which come to dominate the flow. These results are shown to be strikingly different from turbulence distorted by a mean shear flow, when `streaky structures' of high and low streamwise velocity fluctuations develop. It is shown that, in the case of distortion by a mean shear flow, the tendency for the mean shear to produce streamwise vortices by distortion of the turbulent vorticity is largely cancelled by a distortion of the mean vorticity by the turbulent fluctuations. This latter process is absent in distortion by Stokes drift, since there is then no mean vorticity. The components of the Reynolds stress and the integral length scales computed from turbulence distorted by Stokes drift show the same behaviour as in the simulations of Langmuir turbulence reported by McWilliams, Sullivan & Moeng (1997). Hence we suggest that turbulent vorticity in the upper ocean, such as produced by breaking waves, may help to provide the initial seeds for Langmuir circulations, thereby complementing the shear-flow instability mechanism developed by Craik & Leibovich (1976). The tilting of the vertical vorticity into the horizontal by the Stokes drift tends also to produce a shear stress that does work against the mean straining associated with the wave orbital motions. The turbulent kinetic energy then increases at the expense of energy in the wave. Hence the wave decays. An expression for the wave attenuation rate is obtained by scaling the equation for the wave energy, and is found to be broadly consistent with available laboratory data.

Use of airborne remote sensing to detect riverside Brassica rapa to aid in risk assessment of transgenic crops

High resolution descriptions of plant distribution have utility for many ecological applications but are especially useful for predictive modelling of gene flow from transgenic crops. Difficulty lies in the extrapolation errors that occur when limited ground survey data are scaled up to the landscape or national level. This problem is epitomized by the wide confidence limits generated in a previous attempt to describe the national abundance of riverside Brassica rapa (a wild relative of cultivated rapeseed) across the United Kingdom. Here, we assess the value of airborne remote sensing to locate B. rapa over large areas and so reduce the need for extrapolation. We describe results from flights over the river Nene in England acquired using Airborne Thematic Mapper (ATM) and Compact Airborne Spectrographic Imager (CASI) imagery, together with ground truth data. It proved possible to detect 97% of flowering B. rapa on the basis of spectral profiles. This included all stands of plants that occupied >2m square (>5 plants), which were detected using single-pixel classification. It also included very small populations (<5 flowering plants, 1-2m square) that generated mixed pixels, which were detected using spectral unmixing. The high detection accuracy for flowering B. rapa was coupled with a rather large false positive rate (43%). The latter could be reduced by using the image detections to target fieldwork to confirm species identity, or by acquiring additional remote sensing data such as laser altimetry or multitemporal imagery.

Skill Rating by Bayesian Inference

Systems Engineering often involves computer modelling the behaviour of proposed systems and their components. Where a component is human, fallibility must be modelled by a stochastic agent. The identification of a model of decision-making over quantifiable options is investigated using the game-domain of Chess. Bayesian methods are used to infer the distribution of players’ skill levels from the moves they play rather than from their competitive results. The approach is used on large sets of games by players across a broad FIDE Elo range, and is in principle applicable to any scenario where high-value decisions are being made under pressure.

Low density lipoprotein undergoes oxidation within lysosomes in cells

The oxidised low density lipoprotein (LDL) hypothesis of atherosclerosis proposes that LDL undergoes oxidation in the interstitial fluid of the arterial wall. We have shown that aggregated (vortexed) nonoxidised LDL was taken up by J774 mouse macrophages and human monocyte-derived macrophages and oxidised intracellularly, as assessed by the microscopic detection of ceroid, an advanced lipid oxidation product. Confocal microscopy showed that the ceroid was located in the lysosomes. To confirm these findings, J774 macrophages were incubated with acetylated LDL, which is internalised rapidly to lysosomes, and then incubated (chase incubation) in the absence of any LDL. The intracellular levels of oxysterols, measured by HPLC, increased during the chase incubation period, showing that LDL must have been oxidised inside the cells. Furthermore, we found that this oxidative modification was inhibited by lipid-soluble antioxidants, an iron chelator taken up by fluid-phase pinocytosis and the lysosomotropic drug chloroquine, which increases the pH of lysosomes. The results indicate that LDL oxidation can occur intracellularly, most probably within lysosomes.

Low density lipoprotein undergoes oxidation within lysosomes in cells

The oxidized low density lipoprotein (LDL) hypothesis of atherosclerosis proposes that LDL undergoes oxidation in the interstitial fluid of the arterial wall. We have shown that aggregated (vortexed) nonoxidized LDL was taken up by J774 mouse macrophages and human monocyte-derived macrophages and oxidized intracellularly, as assessed by the microscopic detection of ceroid, an advanced lipid oxidation product. Confocal microscopy showed that the ceroid was located in the lysosomes. To confirm these findings, J774 macrophages were incubated with acetylated LDL, which is internalized rapidly to lysosomes, and then incubated (chase incubation) in the absence of any LDL. The intracellular levels of oxysterols, measured by HPLC, increased during the chase incubation period, showing that LDL must have been oxidized inside the cells. Furthermore, we found that this oxidative modification was inhibited by lipid-soluble antioxidants, an iron chelator taken up by fluid-phase pinocytosis and the lysosomotropic drug chloroquine, which increases the pH of lysosomes. The results indicate that LDL oxidation can occur intracellularly, most probably within lysosomes.

Prooxidant and antioxidant properties of human serum ultrafiltrates toward LDL: important role of uric acid

Oxidized LDL is present within atherosclerotic lesions, demonstrating a failure of antioxidant protection. A normal human serum ultrafiltrate of M-r below 500 was prepared as a model for the low M-r components of interstitial fluid, and its effects on LDL oxidation were investigated. The ultrafiltrate (0.3%, v/v) was a potent antioxidant for native LDL, but was a strong prooxidant for mildly oxidized LDL when copper, but not a water-soluble azo initiator, was used to oxidize LDL. Adding a lipid hydroperoxide to native LDL induced the antioxidant to prooxidant switch of the ultrafiltrate. Uric acid was identified, using uricase and add-back experiments, as both the major antioxidant and prooxidant within the ultrafiltrate for LDL. The ultrafiltrate or uric acid rapidly reduced Cu2+ to Cu+. The reduction of Cu2+ to Cu+ may help to explain both the antioxidant and prooxidant effects observed. The decreased concentration of Cu2+ would inhibit tocopherol-mediated peroxidation in native LDL, and the generation of Cu+ would promote the rapid breakdown of lipid hydroperoxides in mildly oxidized LDL into lipid radicals. The net effect of the low M-r serum components would therefore depend on the preexisting levels of lipid hydroperoxides in LDL.jlr These findings may help to explain why LDL oxidation occurs in atherosclerotic lesions in the presence of compounds that are usually considered to be antioxidants.

The surface structure of BaO on Pt(111): (2x2)-reconstructed BaO(111)

The surface structure of BaO(111) has been determined using STM and computer modelling. The BaO(111) surface was prepared in thin film form on Pt(111) and presents a surface with twice the lattice parameter expected for that of the bulk termination, i.e. a (2 x 2) reconstruction. Computer modelling indicates that the bulk termination is unstable, but that the (2 x 2) reconstructed BaO(111) surface has a low surface energy and is hence a stable surface reconstruction. The (2 x 2) reconstruction consists of small, three-sided pyramids with (100) oriented sides and either oxygen or barium ions at the apices. Less regular surface reconstructions containing the same pyramids are almost equally stable, indicating that we may also expect less regular regions to appear with a fairly random distribution of these surface species. The simulations further suggest that a regular (4 x 4) reconstruction built up of bigger pyramids is even more energetically favourable, and some evidence is found for such a structure in the STM. (c) 2006 Elsevier B.V. All rights reserved.

An HMV mapping study of a pulsating jet system

A new system for the generation of hydrodynamic modulated voltammetry (HMV) is presented. This system consists of an oscillating jet produced through the mechanical vibration of a large diaphragm. The structure of the cell is such that a relatively small vibration is transferred to a large fluid flow at the jet outlet. Positioning of an electrode (Pt, 0.5 mm or 25 mu m diameter) over the exit of this jet enables the detection of the modulated flow of liquid. While this flow creates modest mass transfer rates (time averaged similar to 0.015 cm s(-1)) it can also be used to create a HMV system where a 'lock-in' approach is adopted to investigate the redox chemistry in question. This is demonstrated for the Fe(CN)(6)(3-/4-) redox system. Here 'lock-in' to the modulated hydrodynamic signal is achieved through the deployment of bespoke software. The apparatus and procedure is shown to produce a simple and efficient way to obtain the desired signal. In addition the spatial variation of the HMV signal, phase correction and time averaged current with respect to the jet orifice is presented. (C) 2008 Elsevier B.V. All rights reserved.

Approximate Riemann solutions of the two-dimensional shallow-water equations

A finite-difference scheme based on flux difference splitting is presented for the solution of the two-dimensional shallow-water equations of ideal fluid flow. A linearised problem, analogous to that of Riemann for gasdynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearised problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second-order scheme which avoids non-physical, spurious oscillations. An extension to the two-dimensional equations with source terms, is included. The scheme is applied to a dam-break problem with cylindrical symmetry.