Analytical solutions of linear finite- and small-strain one-dimensional consolidation

Analytical solutions are presented for linear finite-strain one-dimensional consolidation of initially unconsolidated soil layers with surcharge loading for both one- and two-way drainage. These solutions complement earlier solutions for initially unconsolidated soil layers without surcharge and initially normally consolidated soil layers with surcharge. Small-strain solutions for the consolidation of initially unconsolidated soil layers with surcharge loading are also presented, and the relationship between the earlier solutions for initially unconsolidated soil without surcharge and the corresponding small-strain solutions, which was not addressed in the earlier work, is clarified. The new solutions for initially unconsolidated soil with surcharge loading can be applied to the analysis of low stress consolidation tests and to the partial validation of numerical solutions of non-linear finite-strain consolidation. They also clarify a formerly perplexing aspect of finite-strain solution charts first noted in numerical solutions. Copyright (C) 2004 John Wiley Sons, Ltd.

On the chaotic rotation of a liquid-filled gyrostat via the Melnikov-Holmes-Marsden integral

The non-linear motions of a gyrostat with an axisymmetrical, fluid-filled cavity are investigated. The cavity is considered to be completely filled with an ideal incompressible liquid performing uniform rotational motion. Helmholtz theorem, Euler's angular momentum theorem and Poisson equations are used to develop the disturbed Hamiltonian equations of the motions of the liquid-filled gyrostat subjected to small perturbing moments. The equations are established in terms of a set of canonical variables comprised of Euler angles and the conjugate angular momenta in order to facilitate the application of the Melnikov-Holmes-Marsden (MHM) method to investigate homoclinic/heteroclinic transversal intersections. In such a way, a criterion for the onset of chaotic oscillations is formulated for liquid-filled gyrostats with ellipsoidal and torus-shaped cavities and the results are confirmed via numerical simulations. (c) 2006 Elsevier Ltd. All rights reserved.

Emergent temporal behaviour and collaborative work

Although collaboration manifestly takes place in time, the role of time in shaping the behaviour of collaborations, and collaborative systems, is not well understood. Time is more than clock-time or the subjective experience of time; its effects on systems include differential rates of change of system elements, temporally non-linear behaviour and phenomena such as entrainment and synchronization. As a system driver, it generates emergent effects shaping systems and their behaviour. In the paper we present a systems view of time, and consider the implications of such a view through the case of collaborative development of a new university timetabling system. Teasing out the key temporal phenomena using the notion of temporal trajectories helps us understand the emergent temporal behaviour and suggests a means for improving outcomes.

Bifurcations in the non-dissipative dynamics of ultra cold atoms

We present the first experimental observation of several bifurcations in a controllable non-linear Hamiltonian system. Dynamics of cold atoms are used to test predictions of non-linear, non-dissipative Hamiltonian dynamics.

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.

Computationally efficient solution techniques for adsorption problems involving steep gradients in bidisperse particles

A piecewise uniform fitted mesh method turns out to be sufficient for the solution of a surprisingly wide variety of singularly perturbed problems involving steep gradients. The technique is applied to a model of adsorption in bidisperse solids for which two fitted mesh techniques, a fitted-mesh finite difference method (FMFDM) and fitted mesh collocation method (FMCM) are presented. A combination (FMCMD) of FMCM and the DASSL integration package is found to be most effective in solving the problems. Numerical solutions (FMFDM and FMCMD) were found to match the analytical solution when the adsorption isotherm is linear, even under conditions involving steep gradients for which global collocation fails. In particular, FMCMD is highly efficient for macropore diffusion control or micropore diffusion control. These techniques are simple and there is no limit on the range of the parameters. The techniques can be applied to a variety of adsorption and desorption problems in bidisperse solids with non-linear isotherm and for arbitrary particle geometry.

Relative dispersions of intra-albumin transit times across rat and elasmobranch perfused livers, and implications for intra- and inter-species scaling of hepatic clearance using microsomal data

It is recognized that vascular dispersion in the liver is a determinant of high first-pass extraction of solutes by that organ. Such dispersion is also required for translation of in-vitro microsomal activity into in-vivo predictions of hepatic extraction for any solute. We therefore investigated the relative dispersion of albumin transit times (CV2) in the livers of adult and weanling rats and in elasmobranch livers. The mean and normalized variance of the hepatic transit time distribution of albumin was estimated using parametric non-linear regression (with a correction for catheter influence) after an impulse (bolus) input of labelled albumin into a single-pass liver perfusion. The mean +/- s.e. of CV2 for albumin determined in each of the liver groups were 0.85 +/- 0.20 (n = 12), 1.48 +/- 0.33 (n = 7) and 0.90 +/- 0.18 (n = 4) for the livers of adult and weanling rats and elasmobranch livers, respectively. These CV2 are comparable with that reported previously for the dog and suggest that the CV2 Of the liver is of a similar order of magnitude irrespective of the age and morphological development of the species. It might, therefore, be justified, in the absence of other information, to predict the hepatic clearances and availabilities of highly extracted solutes by scaling within and between species livers using hepatic elimination models such as the dispersion model with a CV2 of approximately unity.

Programming internet based DAI applications in Qu-Prolog

This paper presents the unique collection of additional features of Qu-Prolog, a variant of the Al programming language Prolog, and illustrates how they can be used for implementing DAI applications. By this we mean applications comprising communicating information servers, expert systems, or agents, with sophisticated reasoning capabilities and internal concurrency. Such an application exploits the key features of Qu-Prolog: support for the programming of sound non-clausal inference systems, multi-threading, and high level inter-thread message communication between Qu-Prolog query threads anywhere on the internet. The inter-thread communication uses email style symbolic names for threads, allowing easy construction of distributed applications using public names for threads. How threads react to received messages is specified by a disjunction of reaction rules which the thread periodically executes. A communications API allows smooth integration of components written in C, which to Qu-Prolog, look like remote query threads.

Willingness to pay for conservation of the Asian elephant in Sri Lanka: A contingent valuation study

Results from a CVM survey of willingness to pay for the conservation of the Asian elephant of a sample of urban residents in three selected housing schemes in Colombo, the capital of Sri Lanka, are reported. Face– to–face surveys were conducted using an interview schedule. A non-linear logit regression model was constructed to analyse the respondents’ responses for the payment principle questions and to identify the factors that influence their responses. We investigate whether urban residents’ WTP for the conservation of elephants is sufficient to compensate farmers for the damage caused by elephants, and consequently to raise farmers’ tolerance of the presence of elephants on the farming fields. We find that beneficiaries (the urban residents) could compensate losers (the farmers in the HEC affected areas) and be better off than in the absence of elephants in Sri Lanka. This suggests that there is a strong economic case for the conservation of the wild elephant population in Sri Lanka. However, we have insufficient data to determine Sri Lanka’s optimal elephant population in the Kaldor-Hicks sense.

The net benefit of saving the Asian Elephant: A policy and contingent valuation Study

Reports results from a contingent valuation survey of willingness to pay for the conservation of the Asian elephant of a sample of urban residents living in three selected housing schemes in Colombo, the capital of Sri Lanka. Face–to–face surveys were conducted using an interview schedule. A non-linear logit regression model is used to analyse the respondents’ responses for the payment principle questions and to identify the factors that influence their responses. We investigate whether urban residents’ willingness to pay for the conservation of elephants is sufficient to compensate farmers for the damage caused by elephants. We find that the beneficiaries (the urban residents) could compensate losers (the farmers in the areas affected by human–elephant conflict) and be better off than in the absence of elephants in Sri Lanka. Therefore, there is a strong economic case for the conservation of the wild elephant population in Sri Lanka. However, we have insufficient data to determine the optimal level of this elephant population in the Kaldor-Hicks sense. Nevertheless, the current population of elephant in Sri Lanka is Kaldor-Hicks preferable to having none.

Quantitative analysis of the time courses of enzyme-catalyzed reactions

The catalytic properties of enzymes are usually evaluated by measuring and analyzing reaction rates. However, analyzing the complete time course can be advantageous because it contains additional information about the properties of the enzyme. Moreover, for systems that are not at steady state, the analysis of time courses is the preferred method. One of the major barriers to the wide application of time courses is that it may be computationally more difficult to extract information from these experiments. Here the basic approach to analyzing time courses is described, together with some examples of the essential computer code to implement these analyses. A general method that can be applied to both steady state and non-steady-state systems is recommended. (C) 2001 academic Press.

Flexural waves transmitted by rectangular piezoceramic transducers

Rectangular piezoceramic transducers are widely used in ultrasonic evaluation and health monitoring techniques and structural vibration control applications. In this paper the flexural waves excited by rectangular transducers adhesively attached to isotropic plates are investigated. In view of the difficulties in developing accurate analytical models describing the transfer characteristics of the transducer due to the complex electromechanical transduction processes and transducer-structure interactions involved, a combined theoretical-experimental approach is developed. A multiple integral transform method is used to describe the propagation behaviour of the waves in the plates, while a heterodyne Doppler laser vibrometer is employed as a non-contact receiver device. This combined theoretical-experimental approach enables the efficient characterization of the electromechanical transfer properties of the piezoelectric transducer which is essential for the development of optimized non-destructive evaluation systems. The results show that the assumption of a uniform contact pressure distribution between the transducer and the plate can accurately predict the frequency spectrum and time domain response signals of the propagating waves along the main axes of the rectangular transmitter element.

Parametric signal regeneration

We present an ultra-high bandwidth all-optical digital signal regeneration device concept utilising non-degenerate parametric interaction in a one-dimensional waveguide. Performance is analysed in terms of re-amplification, re-timing, and re-shaping (including centre frequency correction) of time domain multiplexed signals. Bandwidths of 10-100 THz are achievable. (C) 2001 Published by Elsevier Science B.V.