Solute transport by surface runoff from low-angle slopes: theory and application

The removal of chemicals in solution by overland how from agricultural land has the potential to be a significant source of chemical loss where chemicals are applied to the soil surface, as in zero tillage and surface-mulched farming systems. Currently, we lack detailed understanding of the transfer mechanism between the soil solution and overland flow, particularly under field conditions. A model of solute transfer from soil solution to overland flow was developed. The model is based on the hypothesis that a solute is initially distributed uniformly throughout the soil pore space in a thin layer at the soil surface. A fundamental assumption of the model is that at the time runoff commences, any solute at the soil surface that could be transported into the soil with the infiltrating water will already have been convected away from the area of potential exchange. Solute remaining at the soil surface is therefore not subject to further infiltration and may be approximated as a layer of tracer on a plane impermeable surface. The model fitted experimental data very well in all but one trial. The model in its present form focuses on the exchange of solute between the soil solution and surface water after the commencement of runoff. Future model development requires the relationship between the mass transfer parameters of the model and the time to runoff: to be defined. This would enable the model to be used for extrapolation beyond the specific experimental results of this study. The close agreement between experimental results and model simulations shows that the simple transfer equation proposed in this study has promise for estimating solute loss to surface runoff. Copyright (C) 2000 John Wiley & Sons, Ltd.

Computation of bifurcation boundaries for power systems: A new Delta-plane method

This paper is devoted to the problems of finding the load flow feasibility, saddle node, and Hopf bifurcation boundaries in the space of power system parameters. The first part contains a review of the existing relevant approaches including not-so-well-known contributions from Russia. The second part presents a new robust method for finding the power system load flow feasibility boundary on the plane defined by any three vectors of dependent variables (nodal voltages), called the Delta plane. The method exploits some quadratic and linear properties of the load now equations and state matrices written in rectangular coordinates. An advantage of the method is that it does not require an iterative solution of nonlinear equations (except the eigenvalue problem). In addition to benefits for visualization, the method is a useful tool for topological studies of power system multiple solution structures and stability domains. Although the power system application is developed, the method can be equally efficient for any quadratic algebraic problem.

Finite element analysis of heat transfer and mineralization in layered hydrothermal systems with upward throughflow

We use the finite element method to solve the coupled problem between convective pore-fluid flow, heat transfer and mineralization in layered hydrothermal systems with upward throughflow. In particular, we present the improved rock alteration index (IRAI) concept for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in the systems. To validate the numerical method used in the computation, analytical solutions to a benchmark problem have been derived. After the numerical method is validated, it is used to investigate the pattern of pore-fluid Aom, the distribution of temperature and the mineralization pattern of gold minerals in a layered hydrothermal system with upward throughflow. The related numerical results have demonstrated that the present concept of IRAI is useful and applicable for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in hydrothermal systems. (C) 2000 Elsevier Science S.A. All rights reserved.

Identifying regional dust transport pathways: Application of kinematic trajectory modelling to a trans-Tasman case

PI kinematic trajectory model is used to investigate potential pathways of dust transport from Australia to New Zealand. Historically, these have been assumed to follow rather direct west-east trajectories spanning 2 to 3 days, often resulting in red snow events in the Southern Alps of New Zealand. However, results from the present study which examined the route taken by air parcels originating in southern Australia during dust storms on 24 and 25 May 1994, indicate that trans-Tasman dust transport trajectories are more diverse than previously thought, and display considerable variation during single events. These mon divergent pathways tie in more closely with aeolian dust sedimentation patterns identified by ocean coring in the Tasman Sea, and may account for the deposition of Australian dust on sub-Antarctic islands located well south of the Australian continent. Copyright 2000 John Wiley Sons, Ltd.

Decoherence and coherent population transfer between two coupled systems

We show that an arbitrary system described by two dipole moments exhibits coherent superpositions of internal states that can be completely decoupled fi om the dissipative interactions (responsible for decoherence) and an external driving laser field. These superpositions, known as dark or trapping states, can he completely stable or can coherently interact with the remaining states. We examine the master equation describing the dissipative evolution of the system and identify conditions for population trapping and also classify processes that can transfer the population to these undriven and nondecaying states. It is shown that coherent transfers are possible only if the two systems are nonidentical, that is the transitions have different frequencies and/or decay rates. in particular, we find that the trapping conditions can involve both coherent and dissipative interactions, and depending on the energy level structure of the system, the population can be trapped in a linear superposition of two or more bare states, a dressed state corresponding to an eigenstate of the system plus external fields or, in some cases. in one of the excited states of the system. A comprehensive analysis is presented of the different processes that are responsible for population trapping, and we illustrate these ideas with three examples of two coupled systems: single V- and Lambda-type three-level atoms and two nonidentical tao-level atoms, which are known to exhibit dark states. We show that the effect of population trapping does not necessarily require decoupling of the antisymmetric superposition from the dissipative interactions. We also find that the vacuum-induced coherent coupling between the systems could be easily observed in Lambda-type atoms. Our analysis of the population trapping in two nonidentical atoms shows that the atoms can be driven into a maximally entangled state which is completely decoupled from the dissipative interaction.

A model of specification-based testing of interactive systems

In this paper we present a model of specification-based testing of interactive systems. This model provides the basis for a framework to guide such testing. Interactive systems are traditionally decomposed into a functionality component and a user interface component; this distinction is termed dialogue separation and is the underlying basis for conceptual and architectural models of such systems. Correctness involves both proper behaviour of the user interface and proper computation by the underlying functionality. Specification-based testing is one method used to increase confidence in correctness, but it has had limited application to interactive system development to date.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Coupled experimental and thermodynamic modeling studies for metallurgical smelting and coal combustion slag systems

An extensive research program focused on the characterization of various metallurgical complex smelting and coal combustion slags is being undertaken. The research combines both experimental and thermodynamic modeling studies. The approach is illustrated by work on the PbO-ZnO-Al2O3-FeO-Fe2O3-CaO-SiO2 system. Experimental measurements of the liquidus and solidus have been undertaken under oxidizing and reducing conditions using equilibration, quenching, and electron probe X-ray microanalysis. The experimental program has been planned so as to obtain data for thermodynamic model development as well as for pseudo-ternary Liquidus diagrams that can be used directly by process operators. Thermodynamic modeling has been carried out using the computer system FACT, which contains thermodynamic databases with over 5000 compounds and evaluated solution models. The FACT package is used for the calculation of multiphase equilibria in multicomponent systems of industrial interest. A modified quasi-chemical solution model is used for the liquid slag phase. New optimizations have been carried out, which significantly improve the accuracy of the thermodynamic models for lead/zinc smelting and coal combustion processes. Examples of experimentally determined and calculated liquidus diagrams are presented. These examples provide information of direct relevance to various metallurgical smelting and coal combustion processes.

Robust Mixture Modelling Using the t Distribution

Normal mixture models are being increasingly used to model the distributions of a wide variety of random phenomena and to cluster sets of continuous multivariate data. However, for a set of data containing a group or groups of observations with longer than normal tails or atypical observations, the use of normal components may unduly affect the fit of the mixture model. In this paper, we consider a more robust approach by modelling the data by a mixture of t distributions. The use of the ECM algorithm to fit this t mixture model is described and examples of its use are given in the context of clustering multivariate data in the presence of atypical observations in the form of background noise.

Numerical modelling of single-layer folding: clarification of an issue regarding the possible effect of computer codes and the influence of initial irregularities

The influence of initial perturbation geometry and material propel-ties on final fold geometry has been investigated using finite-difference (FLAC) and finite-element (MARC) numerical models. Previous studies using these two different codes reported very different folding behaviour although the material properties, boundary conditions and initial perturbation geometries were similar. The current results establish that the discrepancy was not due to the different computer codes but due to the different strain rates employed in the two previous studies (i.e. 10(-6) s(-1) in the FLAC models and 10(-14) s(-1) in the MARC models). As a result, different parts of the elasto-viscous rheological field were bring investigated. For the same material properties, strain rate and boundary conditions, the present results using the two different codes are consistent. A transition in Folding behaviour, from a situation where the geometry of initial perturbation determines final fold shape to a situation where material properties control the final geometry, is produced using both models. This transition takes place with increasing strain rate, decreasing elastic moduli or increasing viscosity (reflecting in each case the increasing influence of the elastic component in the Maxwell elastoviscous rheology). The transition described here is mechanically feasible but is associated with very high stresses in the competent layer (on the order of GPa), which is improbable under natural conditions. (C) 2000 Elsevier Science Ltd. All rights reserved.

Mechanisms of substitution reactions on cyclometallated platinum(IV) complexes: "Quasi-labile" systems

The substitution reactions of SMe2 by phosphines (PMePh2, PEtPh2, PPh3, P(4-MeC6H4)(3), P(3-MeC6H4)(3), PCy3) on Pt-IV complexes having a cyclometalated imine ligand, two methyl groups in a cis-geometrical arrangement, a halogen, and a dimethyl sulfide as ligands, [Pt(CN)(CH3)(2)(X)(SMe2)], have been studied as a function of temperature, solvent, and electronic and steric characteristics of the phosphines and the X and CN ligands. In all cases, a limiting dissociative mechanism has been found, where the dissociation of the SMe2 ligand corresponds to the rate-determining step. The pentacoordinated species formed behaves as a true pentacoordinated Pt-IV compound in a steady-state concentration, given the solvent independence of the rate constant. The X-ray crystal structures of two of the dimethyl sulfide complexes and a derivative of the pentacoordinate intermediate have been determined. Differences in the individual rate constants for the entrance of the phosphine ligand can only be estimated as reactivity ratios. In all cases an effect of the phosphine size is detected, indicating that an associative step takes place from the pentacoordinated intermediate. The nature of the (CN) imine and X ligands produces differences in the dimethyl sulfide dissociation reactions rates, which can be quantified by the corresponding DeltaS double dagger values (72, 64, 48, 31, and 78 J K-1 mol(-1) for CN/X being C6H4CHNCH2C6H5/Br, C6H4CHNCH2-(2,4,6-(CH3)(3))C6H2/Br, C6H4CHNCH2C6H5/Cl, C6Cl4CHNCH2C6H5/Cl, and C6W4CH2NCHC6H5/ Pr, respectively). As a whole, the donor character of the coordinated C-aromatic and X atoms have the greatest influence on the dissociativeness of the rate-determining step.

The fine structure of (v, 3) directed triple systems: v equivalent to 2 (mod 3)

The fine structure of a directed triple system of index lambda is the vector (c(1), c(2),...,C-lambda), where c(i) is the number of directed triples appearing precisely i times in the system. We determine necessary and sufficient conditions for a vector to be the fine structure of a directed triple system of index 3 for upsilon = 2 (mod 3).

Experimental study of phase equilibria in the systems Fe-Zn-O and Fe-Zn-Si-O at metallic iron saturation

The reported experimental work on the systems Fe-Zn-O and Fe-Zn-Si-O in equilibrium with metallic iron is part of a wider research program that combines experimental and thermodynamic computer modeling techniques to characterize zinc/lead industrial slags and sinters in the system PbO-ZnO-SiO2-CaO-FeO-Fe2O3. Extensive experimental,investigations using high-temperature equilibration and quenching techniques followed by electron probe X-ray microanalysis (EPMA) were carried out. Special experimental; procedures were developed to enable accurate measurements in these ZnO-containing systems to be performed in equilibrium with metallic iron; The systems Fe-Zn-O and FeZn-Si-O were experimentally investigated in equilibrium with metallic iron in the temperature ranges 900 degreesC to 1200 degreesC (1173 to 1473 K) and from 1000 degreesC to 1350 degreesC (1273 to 1623 K), respectively. The liquidus surface in the system Fe-Zn-Si-O in equilibrium with metallic iron was characterized in the composition ranges 0 to 33 wt pet ZnO and 0 to 40 wt pet SiO2. The wustite (Fe,Zn)O, zincite (Zn,Fe)O, willemite (Zn,Fe)(2)SiO4, arid fayalite: (Fe,Zn)(2)SiO4 solid solutions in equilibrium with metallic iron were measured.

Hybrid modelling of coupled pore fluid-solid deformation problems

A hybrid formulation for coupled pore fluid-solid deformation problems is proposed. The scheme is a hybrid in the sense that we use a vertex centered finite volume formulation for the analysis of the pore fluid and a particle method for the solid in our model. The pore fluid formally occupies the same space as the solid particles. The size of the particles is not necessarily equal to the physical size of materials. A finite volume mesh for the pore fluid flow is generated by Delaunay triangulation. Each triangle possesses an initial porosity. Changes of the porosity are specified by the translations of the mass centers of particles. Net pore pressure gradients are applied to the particle centers and are considered in the particle momentum balance. The potential of our model is illustrated by means of a simulation of coupled fracture and fluid flow developed in porous rock under biaxial compression condition.

A model for predicting the future incidence of coronary heart disease within percentiles of coronary heart disease risk

Background We present a method (The CHD Prevention Model) for modelling the incidence of fatal and nonfatal coronary heart disease (CHD) within various CHD risk percentiles of an adult population. The model provides a relatively simple tool for lifetime risk prediction for subgroups within a population. It allows an estimation of the absolute primary CHD risk in different populations and will help identify subgroups of the adult population where primary CHD prevention is most appropriate and cost-effective. Methods The CHD risk distribution within the Australian population was modelled, based on the prevalence of CHD risk, individual estimates of integrated CHD risk, and current CHD mortality rates. Predicted incidence of first fatal and nonfatal myocardial infarction within CHD risk strata of the Australian population was determined. Results Approximately 25% of CHD deaths were predicted to occur amongst those in the top 10 percentiles of integrated CHD risk, regardless of age group or gender. It was found that while all causes survival did not differ markedly between percentiles of CHD risk before the ages of around 50-60, event-free survival began visibly to differ about 5 years earlier. Conclusions The CHD Prevention Model provides a means of predicting future CHD incidence amongst various strata of integrated CHD risk within an adult population. It has significant application both in individual risk counselling and in the identification of subgroups of the population where drug therapy to reduce CHD risk is most cost-effective. J Cardiovasc Risk 8:31-37 (C) 2001 Lippincott Williams & Wilkins.