The programming language python in earth system simulations

-scale vary from a planetary scale and million years for convection problems to 100km and 10 years for fault systems simulations. Various techniques are in use to deal with the time dependency (e.g. Crank-Nicholson), with the non-linearity (e.g. Newton-Raphson) and weakly coupled equations (e.g. non-linear Gauss-Seidel). Besides these high-level solution algorithms discretization methods (e.g. finite element method (FEM), boundary element method (BEM)) are used to deal with spatial derivatives. Typically, large-scale, three dimensional meshes are required to resolve geometrical complexity (e.g. in the case of fault systems) or features in the solution (e.g. in mantel convection simulations). The modelling environment escript allows the rapid implementation of new physics as required for the development of simulation codes in earth sciences. Its main object is to provide a programming language, where the user can define new models and rapidly develop high-level solution algorithms. The current implementation is linked with the finite element package finley as a PDE solver. However, the design is open and other discretization technologies such as finite differences and boundary element methods could be included. escript is implemented as an extension of the interactive programming environment python (see www.python.org). Key concepts introduced are Data objects, which are holding values on nodes or elements of the finite element mesh, and linearPDE objects, which are defining linear partial differential equations to be solved by the underlying discretization technology. In this paper we will show the basic concepts of escript and will show how escript is used to implement a simulation code for interacting fault systems. We will show some results of large-scale, parallel simulations on an SGI Altix system. Acknowledgements: Project work is supported by Australian Commonwealth Government through the Australian Computational Earth Systems Simulator Major National Research Facility, Queensland State Government Smart State Research Facility Fund, The University of Queensland and SGI.

Towards identifying the new structures formed on the γ-radiolysis of Ultem

Ultem irradiated up to 10.0 MGy has been analysed using C-13, H-1 and D-2 proton-carbon and proton-proton correlation NMR spectroscopy to shed light on the formation of new structures. Chemical shifts and correlation data were used to determine the structure or partial structures of several new components. The spectra indicated the presence of new groups and structures involving the isopropylidene group, the imide ring, and hydrogen-abstraction reactions. Possible pathways for formation of the new structures are proposed and the G-values for their formation have been estimated. (C) 2003 Elsevier Science Ltd. All rights reserved.

A new approach to the synthesis of nanocrystal conjugated polymer composites

A novel one pot process has been developed for the preparation of PbS nanocrystals in the conjugated polymer poly 2-methoxy,5-(2 ethyl-hexyloxy-p-phenylenevinylene) (MEH-PPV). Current techniques for making such composite materials rely upon synthesizing the nanocrystals and conducting polymer separately, and subsequently mixing them. This multi-step technique has two serious drawbacks: templating surfactant must be removed before mixing, and co-solvent incompatibility causes aggregation. In our method, we eliminate the need for an initial surfactant by using the conducting polymer to terminate and template nanocrystal growth. Additionally, the final product is soluble in a single solvent. We present materials analysis which shows PbS nanocrystals can be grown directly in a conducting polymer, the resulting composite is highly ordered and nanocrystal size can be controlled.