An algorithm for isentropic flow

An efficient algorithm is presented for the solution of the equations of isentropic gas dynamics with a general convex gas law. The scheme is based on solving linearized Riemann problems approximately, and in more than one dimension incorporates operator splitting. In particular, only two function evaluations in each computational cell are required. The scheme is applied to a standard test problem in gas dynamics for a polytropic gas

An efficient flux difference splitting algorithm for unsteady duct flows of a real gas

A numerical scheme is presented for the solution of the Euler equations of compressible flow of a real gas in a single spatial coordinate. This includes flow in a duct of variable cross-section, as well as flow with slab, cylindrical or spherical symmetry, as well as the case of an ideal gas, and can be useful when testing codes for the two-dimensional equations governing compressible flow of a real gas. The resulting scheme requires an average of the flow variables across the interface between cells, and this average is chosen to be the arithmetic mean for computational efficiency, which is in contrast to the usual “square root” averages found in this type of scheme. The scheme is applied with success to five problems with either slab or cylindrical symmetry and for a number of equations of state. The results compare favourably with the results from other schemes.

An efficient shock capturing algorithm for compressible flows in a duct of variable cross-section

A numerical scheme is presented for the solution of the Euler equations of compressible flow of a gas in a single spatial co-ordinate. This includes flow in a duct of variable cross-section as well as flow with slab, cylindrical or spherical symmetry and can prove useful when testing codes for the two-dimensional equations governing compressible flow of a gas. The resulting scheme requires an average of the flow variables across the interface between cells and for computational efficiency this average is chosen to be the arithmetic mean, which is in contrast to the usual ‘square root’ averages found in this type of scheme. The scheme is applied with success to five problems with either slab or cylindrical symmetry and a comparison is made in the cylindrical case with results from a two-dimensional problem with no sources.

An algorithm for the shallow water equations with body fitted meshes

An efficient algorithm based on flux difference splitting is presented for the solution of the two-dimensional shallow water equations in a generalised coordinate system. The scheme is based on solving linearised Riemann problems approximately and in more than one dimension incorporates operator splitting. The scheme has good jump capturing properties and the advantage of using body-fitted meshes. Numerical results are shown for flow past a circular obstruction.

An isentropic flow algorithm with body fitted meshes

An efficient algorithm based on flux difference splitting is presented for the solution of the three-dimensional equations of isentropic flow in a generalised coordinate system, and with a general convex gas law. The scheme is based on solving linearised Riemann problems approximately and in more than one dimension incorporates operator splitting. The algorithm requires only one function evaluation of the gas law in each computational cell. The scheme has good shock capturing properties and the advantage of using body-fitted meshes. Numerical results are shown for Mach 3 flow of air past a circular cylinder. Furthermore, the algorithm also applies to shallow water flows by employing the familiar gas dynamics analogy.

A numerical scheme for two-dimensional, open channel flows with non-rectangular geometries

An algorithm based on flux difference splitting is presented for the solution of two-dimensional, open channel flows. A transformation maps a non-rectangular, physical domain into a rectangular one. The governing equations are then the shallow water equations, including terms of slope and friction, in a generalized coordinate system. A regular mesh on a rectangular computational domain can then be employed. The resulting scheme has good jump capturing properties and the advantage of using boundary/body-fitted meshes. The scheme is applied to a problem of flow in a river whose geometry induces a region of supercritical flow.

The surface geometry of carbonmonoxide and hydrogen co-adsorbed on Ni 111

A combination of photoelectron spectroscopy, temperature programmed desorption and low energy electron diffraction structure determinations have been applied to study the p(2 x 2) structures of pure hydrogen and co-adsorbed hydrogen and CO on Ni {111}. In agreement with earlier work atomic hydrogen is found to adsorb on fcc and hcp sites in the pure layer with H-Ni bond lengths of 1.74Angstrom. The substrate interlayer distances, d(12) = 2.05Angstrom and d(23) = 2.06Angstrom, are expanded with respect to clean Ni {111} with buckling of 0.04Angstrom in the first layer. In the co-adsorbed phase Co occupies hcp sites and only the hydrogen atoms on fcc sites remain on the surface. d(12) is even further expanded to 2.08Angstrom with buckling in the first and second layer of 0.06 and 0.02Angstrom, respectively. The C-O, C-Ni, and H-Ni bond lengths are within the range of values also found for the pure adsorbates.

The interplay between geometry, electronic structure, and reactivity of Cu-Ni bimetallic (111) surfaces

The mutual influence of surface geometry (e.g. lattice parameters, morphology) and electronic structure is discussed for Cu-Ni bimetallic (111) surfaces. It is found that on flat surfaces the electronic d-states of the adlayer experience very little influence from the substrate electronic structure which is due to their large separation in binding energies and the close match of Cu and Ni lattice constants. Using carbon monoxide and benzene as probe molecules, it is found that in most cases the reactivity of Cu or Ni adlayers is very similar to the corresponding (111) single crystal surfaces. Exceptions are the adsorption of CO on submonolayers of Cu on Ni(111) and the dissociation of benzene on Ni/Cu(111) which is very different from Ni(111). These differences are related to geometric factors influencing the adsorption on these surfaces.

Adsorption and dissociation of benzene on bimetallic surfaces - the influence of surface geometry and electronic structure

This topical review discusses the influence of the surface geometry (e.g. lattice parameters and termination) and electronic structure of well-defined bimetallic surfaces on the adsorption and dissociation of benzene. The available data can be divided into two categories with combinations of non-transition metals and transition metals on the one side and combinations of two transition metals on the other. The main effect of non-transition metals in surface alloys is site blocking which can suppress chemisorption and dissociation of the molecules completely. When two transition metals are combined, the effects are less dramatic. They mainly affect the strength of the chemisorption bond and the degree of dissociation due to electronic and template effects.

The surface geometry of carbon monoxide and oxygen co-adsorbed on Ni{111}

Low energy electron diffraction (LEED) structure determinations have been performed for the p(2 x 2) structures of pure oxygen and oxygen co-adsorbed with CO on Ni{111}. Optimisation of the non-geometric parameters led to very good agreement between experimental and theoretical IV-curves and hence to a high accuracy in the structural parameters. In agreement with earlier work atomic oxygen is found to adsorb on fee sites in both structures. In the co-adsorbed phase CO occupies atop sites. The positions of the substrate atoms are almost identical, within 0.02 Angstrom, in both structures, implying that the interaction with oxygen dominates the arrangement of Ni atoms at the surface.

The TAMORA algorithm: satellite rainfall estimates over West Africa using multi-spectral SEVIRI data

A multi-spectral rainfall estimation algorithm has been developed for the Sahel region of West Africa with the purpose of producing accumulated rainfall estimates for drought monitoring and food security. Radar data were used to calibrate multi-channel SEVIRI data from MSG, and a probability of rainfall at several different rain-rates was established for each combination of SEVIRI radiances. Radar calibrations from both Europe (the SatPrecip algorithm) and Niger (TAMORA algorithm) were used. 10 day estimates were accumulated from SatPrecip and TAMORA and compared with kriged gauge data and TAMSAT satellite rainfall estimates over West Africa. SatPrecip was found to produce large overestimates for the region, probably because of its non-local calibration. TAMORA was negatively biased for areas of West Africa with relatively high rainfall, but its skill was comparable to TAMSAT for the low-rainfall region climatologically similar to its calibration area around Niamey. These results confirm the high importance of local calibration for satellite-derived rainfall estimates. As TAMORA shows no improvement in skill over TAMSAT for dekadal estimates, the extra cloud-microphysical information provided by multi-spectral data may not be useful in determining rainfall accumulations at a ten day timescale. Work is ongoing to determine whether it shows improved accuracy at shorter timescales.

A cluster-based implementation of a fault-tolerant parallel reduction algorithm using swarm-array computing

Recent research in multi-agent systems incorporate fault tolerance concepts. However, the research does not explore the extension and implementation of such ideas for large scale parallel computing systems. The work reported in this paper investigates a swarm array computing approach, namely ‘Intelligent Agents’. In the approach considered a task to be executed on a parallel computing system is decomposed to sub-tasks and mapped onto agents that traverse an abstracted hardware layer. The agents intercommunicate across processors to share information during the event of a predicted core/processor failure and for successfully completing the task. The agents hence contribute towards fault tolerance and towards building reliable systems. The feasibility of the approach is validated by simulations on an FPGA using a multi-agent simulator and implementation of a parallel reduction algorithm on a computer cluster using the Message Passing Interface.

An analysis of the impact of R&D on productivity using Bayesian model averaging with a reversible jump algorithm.

A Bayesian Model Averaging approach to the estimation of lag structures is introduced, and applied to assess the impact of R&D on agricultural productivity in the US from 1889 to 1990. Lag and structural break coefficients are estimated using a reversible jump algorithm that traverses the model space. In addition to producing estimates and standard deviations for the coe¢ cients, the probability that a given lag (or break) enters the model is estimated. The approach is extended to select models populated with Gamma distributed lags of di¤erent frequencies. Results are consistent with the hypothesis that R&D positively drives productivity. Gamma lags are found to retain their usefulness in imposing a plausible structure on lag coe¢ cients, and their role is enhanced through the use of model averaging.

The effects of variation in snow properties on snow mass estimation using the Chang algorithm

Estimating snow mass at continental scales is difficult but important for understanding landatmosphere interactions, biogeochemical cycles and Northern latitudes’ hydrology. Remote sensing provides the only consistent global observations, but the uncertainty in measurements is poorly understood. Existing techniques for the remote sensing of snow mass are based on the Chang algorithm, which relates the absorption of Earth-emitted microwave radiation by a snow layer to the snow mass within the layer. The absorption also depends on other factors such as the snow grain size and density, which are assumed and fixed within the algorithm. We examine the assumptions, compare them to field measurements made at the NASA Cold Land Processes Experiment (CLPX) Colorado field site in 2002–3, and evaluate the consequences of deviation and variability for snow mass retrieval. The accuracy of the emission model used to devise the algorithm also has an impact on its accuracy, so we test this with the CLPX measurements of snow properties against SSM/I and AMSR-E satellite measurements.