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.