A cell by cell anisotropic adaptive mesh ALE scheme for the numerical solution of the Euler equations

In this paper a cell by cell anisotropic adaptive mesh technique is added to an existing staggered mesh Lagrange plus remap finite element ALE code for the solution of the Euler equations. The quadrilateral finite elements may be subdivided isotropically or anisotropically and a hierarchical data structure is employed. An efficient computational method is proposed, which only solves on the finest level of resolution that exists for each part of the domain with disjoint or hanging nodes being used at resolution transitions. The Lagrangian, equipotential mesh relaxation and advection (solution remapping) steps are generalised so that they may be applied on the dynamic mesh. It is shown that for a radial Sod problem and a two-dimensional Riemann problem the anisotropic adaptive mesh method runs over eight times faster.

Consistent Dirichlet Boundary Conditions for Numerical Solution of Moving Boundary Problems

We consider the imposition of Dirichlet boundary conditions in the finite element modelling of moving boundary problems in one and two dimensions for which the total mass is prescribed. A modification of the standard linear finite element test space allows the boundary conditions to be imposed strongly whilst simultaneously conserving a discrete mass. The validity of the technique is assessed for a specific moving mesh finite element method, although the approach is more general. Numerical comparisons are carried out for mass-conserving solutions of the porous medium equation with Dirichlet boundary conditions and for a moving boundary problem with a source term and time-varying mass.

GDASH: a grid-enabled program for structure solution from powder diffraction data

The simulated annealing approach to structure solution from powder diffraction data, as implemented in the DASH program, is easily amenable to parallelization at the individual run level. Very large scale increases in speed of execution can therefore be achieved by distributing individual DASH runs over a network of computers. The GDASH program achieves this by packaging DASH in a form that enables it to run under the Univa UD Grid MP system, which harnesses networks of existing computing resources to perform calculations.

MDASH: a multi-core-enabled program for structure solution from powder diffraction data

The simulated annealing approach to structure solution from powder diffraction data, as implemented in the DASH program, is easily amenable to parallelization at the individual run level. Modest increases in speed of execution can therefore be achieved by executing individual DASH runs on the individual cores of CPUs.

Voronoi, Delaunay, and Block-Structured Mesh Refinement for Solution of the Shallow-Water Equations on the Sphere

Alternative meshes of the sphere and adaptive mesh refinement could be immensely beneficial for weather and climate forecasts, but it is not clear how mesh refinement should be achieved. A finite-volume model that solves the shallow-water equations on any mesh of the surface of the sphere is presented. The accuracy and cost effectiveness of four quasi-uniform meshes of the sphere are compared: a cubed sphere, reduced latitude–longitude, hexagonal–icosahedral, and triangular–icosahedral. On some standard shallow-water tests, the hexagonal–icosahedral mesh performs best and the reduced latitude–longitude mesh performs well only when the flow is aligned with the mesh. The inclusion of a refined mesh over a disc-shaped region is achieved using either gradual Delaunay, gradual Voronoi, or abrupt 2:1 block-structured refinement. These refined regions can actually degrade global accuracy, presumably because of changes in wave dispersion where the mesh is highly nonuniform. However, using gradual refinement to resolve a mountain in an otherwise coarse mesh can improve accuracy for the same cost. The model prognostic variables are height and momentum collocated at cell centers, and (to remove grid-scale oscillations of the A grid) the mass flux between cells is advanced from the old momentum using the momentum equation. Quadratic and upwind biased cubic differencing methods are used as explicit corrections to a fast implicit solution that uses linear differencing.

The influence of mineral solubility and soil solution concentration on the toxicity of copper to Eisenia fetida Savigny

The OECD 14 d earthworm acute toxicity test was used to determine the toxicity of copper added as copper nitrate (Cu(NO3)(2)), copper sulphate (CuSO4) and malachite (Cu-2(OH)(2)(CO3)) to Eisenia fetida Savigny. Cu(NO3)(2), and CuSO4 were applied in both an aqueous (aq) and solid (s) form, Cu-2(OH)(2)(CO3) was added as a solid. Soil solution was extracted by centrifugation, and analysed for copper. Two extractants [0.01 M CaCl2 and 0.005 M diethylenetriminpentaacetic acid (DTPA)] were used as a proxy of the bioavailable copper fraction in the soil. For bulk soil copper content the calculated copper toxicity decreased in the order nitrate > sulphide > carbonate, the same order as decreasing solubility of the metal compounds. For Cu(NO3)(2) and CuSO4, the LC50s obtained were not significantly different when the compound was added in solution or solid form. There was a significant correlation between the soil solution copper concentration and the percentage earthworm mortality for all 3 copper compounds (P less than or equal to 0.05) indicating that the soil pore water copper concentration is important for determining copper availability and toxicity to E. fetida. In soil avoidance tests the earthworms avoided the soils treated with Cu(NO3)(2) (aq and s) and CuSO4 (aq and s), at all concentrations used (110-8750 mug Cu g(-1), and 600-8750 mug Cu g(-1) respectively). In soils treated with Cu-2(OH2)CO3, avoidance behaviour was exhibited at all concentrations greater than or equal to3500 mug Cu g(-1). There was no significant correlation between the copper extracted by either CaCl2 or DTPA and percentage mortality. These two extractants are therefore not useful indicators of copper availability and toxicity to E. fetida.

Naturally occurring radioactive material (NORM) from a former phosphoric acid processing plant

In recent years there has been an increasing awareness of the radiological impact of non-nuclear industries that extract and/or process ores and minerals containing naturally occurring radioactive material (NORM). These industrial activities may result in significant radioactive contamination of (by-) products, wastes and plant installations. In this study, scale samples were collected from a decommissioned phosphoric acid processing plant. To determine the nature and concentration of NORM retained in pipe-work and associated process plant, four main areas of the site were investigated: (1) the 'Green Acid Plant', where crude acid was concentrated; (2) the green acid storage tanks; (3) the Purified White Acid (PWA) plant, where inorganic impurities were removed; and (4) the solid waste, disposed of on-site as landfill. The scale samples predominantly comprise the following: fluorides (e.g. ralstonite); calcium sulphate (e.g. gypsum); and an assemblage of mixed fluorides and phosphates (e.g. iron fluoride hydrate, calcium phosphate), respectively. The radioactive inventory is dominated by U-238 and its decay chain products, and significant fractionation along the series occurs. Compared to the feedstock ore, elevated concentrations (<= 8.8 Bq/g) of U-238 Were found to be retained in installations where the process stream was rich in fluorides and phosphates. In addition, enriched levels (<= 11 Bq/g) of Ra-226 were found in association with precipitates of calcium sulphate. Water extraction tests indicate that many of the scales and waste contain significantly soluble materials and readily release radioactivity into solution. (c) 2005 Elsevier Ltd. All rights reserved.

Displacement of Zn through acidic light-textured soils

In order to gain understanding of the movement of pollutant metals in soil. the chemical mechanisms involved in the transport of zinc were studied. The displacement of zinc through mixtures of sand and cation exchange resin was measured to validate the methods used for soil. With cation exchange capacities of 2.5 and 5.0 cmol(c) kg(-1). 5.6 and 8.4 pore volumes of 10 mM CaCl2, respectively, were required to displace a pulse of ZnCl2. A simple Burns-type model (Wineglass) using an adsorption coefficient (K-d) determined by fitting a straight line relationship to an adsorption isotherm gave a good fit to the data (K-d=0.73 and 1.29 ml g(-1), respectively). Surface and subsurface samples of an acidic sandy loam (organic matter 4.7 and 1.0%. cation exchange capacity (CEC) 11.8 and 6.1 cmol(c) kg(-1) respectively) were leached with 10 mM calcium chloride. nitrate and perchlorate. With chloride. the zinc pulse was displaced after 25 and 5 pore volumes, respectively. The Kd values were 6.1 and 2.0 ml g(-1). but are based on linear relationships fitted to isotherms which are both curved and show hysteresis. Thus. a simple model has limited value although it does give a general indication of rate of displacement. Leaching with chloride and perchlorate gave similar displacement and Kd values, but slower movement occurred with nitrate in both soil samples (35 and 7 pore volumes, respectively) which reflected higher Kd values when the isotherms were measured using this anion (7.7 and 2.8 ml g(-1) respectively). Although pH values were a little hi-her with nitrate in the leachates, the differences were insufficient to suggest that this increased the CEC enough to cause the delay. No increases in pH occurred with nitrate in the isotherm experiments. Geochem was used to calculate the proportions of Zn complexed with the three anions and with fulvic acid determined from measurements of dissolved organic matter. In all cases, more than 91% of the Zn was present as Zn2+ and there were only minor differences between the anions. Thus, there is an unexplained factor associated with the greater adsorption of Zn in the presence of nitrate. Because as little as five pore volumes of solution displaced Zn through the subsurface soil, contamination of ground waters may be a hazard where Zn is entering a light-textured soil, particularly where soil salinity is increased. Reductions in organic matter content due to cultivation will increase the hazard. (C) 2004 Elsevier B.V. All rights reserved.

Uptake Pathways of Polycyclic Aromatic Hydrocarbons in White Clover

An understanding of the primary pathways of plant uptake of organic pollutants is important to enable the risks from crops grown on contaminated soils to be assessed. A series of experiments were undertaken to quantify the importance of the pathways of contamination and the Subsequent transport within the plant using white clover plants grown in solution culture. Root uptake was primarily an absorption process, but a component of the contamination was a result of the transpiration flux to the shoot for higher Solubility compounds. The root contamination can be easily predicted using a simple relationship with K-OW, although if a composition model was used based on lipid content, a significant under prediction of the contamination was observed. Shoot uptake was driven by the transpiration stream flux which was related to the solubility of the individual PAH rather than the K-OW. However, the experiment was over a short duration, 6 days, and models based on K-OW may be better for crops grown in the field where the vegetation will approach equilibrium and transpiration cannot easily be measured, A significant fraction of the shoot contamination resulted from aerial deposition derived from volatilized PAH. This pathway was more significant for compounds approaching log K-OA > 9 and log K-AW < -3. The shoot uptake pathways need further investigation to enable them to be modeled separately, There was no evidence of significant systemic transport of the PAR so transfer outside the transpiration stream is likely to be limited.

Zinc deficiency in selected cultivars of wheat and barley as tested in solution culture

The relative zinc (Zn) efficiencies of 33 wheat and 3 barley cultivars were determined by growing them in chelate-buffered culture solutions. Zn efficiency, determined by growth in a Zn-deficient solution relative to that in a medium containing an adequate concentration of Zn, was found to vary between 10% and 63% among the cultivars tested. Out of the 36 cultivars tested, 12 proved to be Zn efficient, 10 were Zn inefficient, and the remaining 14 varieties were classed as intermediate. The most Zn-efficient cultivars included Bakhtawar, Gatcher S61, Wilgoyne, and Madrigal, and the most Zn inefficient included Durati, Songlen, Excalibur, and Chakwal-86. Zn-efficient cultivars accumulated greater amounts of Zn in their shoots than inefficient cultivars, but the correlation between shoot Zn and shoot dry matter production was poor. All the cultivars accumulated higher concentrations of iron (Fe), copper (Cu), manganese (Mn), and phosphorus (P) at deficient levels of Zn, compared with adequate Zn concentrations. The Zn-inefficient cultivars accumulated higher concentrations of these other elements compared to efficient cultivars.