Truncation errors in finite difference models for solute transport equation with first-order reaction

The truncation errors associated with finite difference solutions of the advection-dispersion equation with first-order reaction are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation and a temporally and spatially weighted parametric approach is used for differentiating among the various finite difference schemes. The numerical truncation errors are defined using Peclet and Courant numbers and a new Sink/Source dimensionless number. It is shown that all of the finite difference schemes suffer from truncation errors. Tn particular it is shown that the Crank-Nicolson approximation scheme does not have second order accuracy for this case. The effects of these truncation errors on the solution of an advection-dispersion equation with a first order reaction term are demonstrated by comparison with an analytical solution. The results show that these errors are not negligible and that correcting the finite difference scheme for them results in a more accurate solution. (C) 1999 Elsevier Science B.V. All rights reserved.

Therapeutic efficacy of zeta-cypermethrin pour-on for the treatment of biting and sucking lice in cattle under field conditions

Objective To assess the efficacy of zeta-cypermethrin pour-on to control cattle lice. Design Five field trials in south-eastern Australia. Procedure Zeta-cypermethrin pour-on, deltamethrin pour-on and pour-on vehicle were applied to groups of 10 cattle. Lice were counted before treatment and 14, 28, 42 and 56 days after treatment. Results Zeta-cypermethrin pour-on given at 2.5 mg/kg was equivalent to, or marginally more effective than a deltamethrin pour-on at 0.75 mg/kg. It eliminated B bovis and H eurysternus and gave good control of L vituli and S capillatus. Zeta-cypermethrin at 1 mg/kg gave good control of B bovis and H eurysternus but was not satisfactory against L vituli and S capillatus. Conclusion Zeta-cypermethrin pour-on, given at 2.5 mg/kg, is an effective treatment for cattle lice control. Zeta-cypermethrin, and other synthetic pyrethroid pour-ons, are the treatment of choice to control B bovis.

Influence of field-based nutrient enrichment on the photobiology of the giant clam Tridacna maxima

Nutrients were added to 12 microatolls in One Tree Island lagoon every low tide for 13 mo to an initial concentration of 10 mu M (ammonium, N) and 2 mu M (phosphate, P). These concentrations remained above background for 2 to 3 h after addition. The addition of ammonium (N and NI-P but not P alone) significantly increased P, (gross photosynthesis) P,, (net photosynthesis) and R (respiration) per unit wet-tissue weight and cc (photosynthetic efficiency) in Tridacna maxima after 3 mo nutrient enrichment. These responses to small and transient changes in ammonium concentrations suggest that symbiotic clams are not nutrient-replete, and that even subtle changes in nutrients can have a measurable effect on photosynthesis. The same clams did not show significant differences in photosynthetic parameters 6 mo after the beginning of nutrient enrichment, suggesting that their previous responses had either been seasonal or that symbiotic clams such as T. maxima are able to adjust their photophysiology following external changes in nutrient concentrations.

Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method. Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately. Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.

Effect of disorder on quantum phase transitions in anisotropic XY spin chains in a transverse field

We present some exact results for the effect of disorder on the critical properties of an anisotropic XY spin chain in a transverse held. The continuum limit of the corresponding fermion model is taken and in various cases results in a Dirac equation with a random mass. Exact analytic techniques can then be used to evaluate the density of states and the localization length. In the presence of disorder the ferromagnetic-paramagnetic or Ising transition of the model is in the same universality class as the random transverse field Ising model solved by Fisher using a real-space renormalization-group decimation technique (RSRGDT). If there is only randomness in the anisotropy of the magnetic exchange then the anisotropy transition (from a ferromagnet in the x direction to a ferromagnet in the y direction) is also in this universality class. However, if there is randomness in the isotropic part of the exchange or in the transverse held then in a nonzero transverse field the anisotropy transition is destroyed by the disorder. We show that in the Griffiths' phase near the Ising transition that the ground-state energy has an essential singularity. The results obtained for the dynamical critical exponent, typical correlation length, and for the temperature dependence of the specific heat near the Ising transition agree with the results of the RSRODT and numerical work. [S0163-1829(99)07125-8].

Finite element modelling of reactive mass transport problems in fluid-saturated porous media

We use the finite element method to solve reactive mass transport problems in fluid-saturated porous media. In particular, we discuss the mathematical expression of the chemical reaction terms involved in the mass transport equations for an isothermal, non-equilibrium chemical reaction. It has turned out that the Arrhenius law in chemistry is a good mathematical expression for such non-equilibrium chemical reactions especially from the computational point of view. Using the finite element method and the Arrhenius law, we investigate the distributions of PH (i.e. the concentration of H+) and the relevant reactive species in a groundwater system. Although the main focus of this study is on the contaminant transport problems in groundwater systems, the related numerical techniques and principles are equally applicable to the orebody formation problems in the geosciences. Copyright (C) 1999 John Wiley & Sons, Ltd.

Absorption spectrum of a strongly driven atom in a detuned squeezed vacuum

We present numerical and analytical results for the Mollow probe absorption spectrum of a coherently driven two-level system in a narrow bandwidth squeezed vacuum field. The spectra are calculated for the case where the Rabi frequency of the driving field is much larger than the natural linewidth and the squeezed vacuum carrier frequency is detuned from the driving laser frequency. The driving laser is on resonance. We show that in a detuned squeezed vacuum the standard Mellow features are each split into triplets. The central components of each triplet are weakly dependent on the squeezing phase but the sidebands strongly depend on the phase and can have dispersive or absorptive/emissive profiles. We also derive approximate analytical expressions for the spectral features and find that the multi-peak structure of the spectrum can be interpreted either via the eigenfrequencies of a generalized Floquet Hamiltonian or in terms of three-photon transitions between dressed stales involving a probe field photon and a correlated photon pair from the squeezed vacuum field.

The crack tip strain field of AISI 4340 - Part I - Measurement technique

This is the first paper in a study on the influence of the environment on the crack tip strain field for AISI 4340. A stressing stage for the environmental scanning electron microscope (ESEM) was constructed which was capable of applying loads up to 60 kN to fracture-mechanics samples. The measurement of the crack tip strain field required preparation (by electron lithography or chemical etching) of a system of reference points spaced at similar to 5 mu m intervals on the sample surface, loading the sample inside an electron microscope, image processing procedures to measure the displacement at each reference point and calculation of the strain field. Two algorithms to calculate strain were evaluated. Possible sources of errors were calculation errors due to the algorithm, errors inherent in the image processing procedure and errors due to the limited precision of the displacement measurements. Estimation of the contribution of each source of error was performed. The technique allows measurement of the crack tip strain field over an area of 50 x 40 mu m with a strain precision better than +/- 0.02 at distances larger than 5 mu m from the crack tip. (C) 1999 Kluwer Academic Publishers.

The crack tip strain field of AISI 4340 - Part II - Experimental measurements

Crack tip strain maps have been measured for AISI 4340 high strength steel. No significant creep was observed. The measured values of CTOD were greater than expected from the HRR model. Crack tip branching was observed in every experiment. The direction of crack branching was in the same direction as a major ridge'' of epsilon(yy) strain, which in turn was in the same direction as predicted by the HRR model. Furthermore, the measured magnitudes of the epsilon(y)y strain in this same direction were in general greater than the values predicted by the HRR model. This indicates more plasticity in the crack tip region than expected from the HRR model. This greater plasticity could be related to the larger than expected CTOD values. The following discrepancies between the measured strain fields for AISI 4340 and the HRR predictions are noteworthy: (1) The crack branching. (2) Values of CTOD significantly higher than predicted by HRR. (3) The major ridge'' of epsilon(yy) strain an angle of about 60 degrees with the direction of overall propagation of the fatigue precrack, in which the measured magnitudes of the epsilon(yy) strain were greater than the values predicted by the HRR model. (4) Asymmetric shape of the plastic zone as measured by the epsilon(yy) strain. (5) Values of shear strain gamma(xy) significantly higher than predicted by the HRR model. (C) 1999 Kluwer Academic Publishers.

The crack tip strain field of AISI 4340 - Part III - Hydrogen influence

This paper studied the influence of hydrogen and water vapour environments on the plastic behaviour in the vicinity of the crack tip for AISI 4340. Hydrogen and water vapour (at a pressure of 15 Torr) significantly increased the crack tip opening displacement. The crack tip strain distribution in 15 Torr hydrogen was significantly different to that measured in vacuum. In the presence of sufficient hydrogen, the plastic zone was larger, was elongated in the direction of crack propagation and moreover there was significant creep. These observations support the hydrogen enhanced localised plasticity model for hydrogen embrittlement in this steel. The strain distribution in the presence of water vapour also suggests that SCC in AISI 4340 occurs via the hydrogen enhanced localised plasticity mechanism. (C) 1999 Kluwer Academic Publishers.

Formal derivation of finite state machines for class testing

Previous work on generating state machines for the purpose of class testing has not been formally based. There has also been work on deriving state machines from formal specifications for testing non-object-oriented software. We build on this work by presenting a method for deriving a state machine for testing purposes from a formal specification of the class under test. We also show how the resulting state machine can be used as the basis for a test suite developed and executed using an existing framework for class testing. To derive the state machine, we identify the states and possible interactions of the operations of the class under test. The Test Template Framework is used to formally derive the states from the Object-Z specification of the class under test. The transitions of the finite state machine are calculated from the derived states and the class's operations. The formally derived finite state machine is transformed to a ClassBench testgraph, which is used as input to the ClassBench framework to test a C++ implementation of the class. The method is illustrated using a simple bounded queue example.

Paramagnetic limiting of the upper critical field of the layered organic superconductor kappa-(BEDT-TTF)(2)Cu(SCN)(2)

We report detailed measurements of the interlayer magnetoresistance of the layered organic superconductor kappa-(BEDT-TTF)(2)Cu(SCN)(2) for temperatures down to 0.5 K and fields up to 30 T. The upper critical field is determined from the resistive transition for a wide range of temperatures and field directions. For magnetic fields parallel to the layers, the upper critical field increases approximately linearly with decreasing temperature. The upper critical field at low temperatures is compared to the Pauli paramagnetic limit, at which singlet superconductivity should be destroyed by the Zeeman splitting of the electron spins. The measured value is comparable to a value for the paramagnetic limit calculated from thermodynamic quantities but exceeds the limit calculated from BCS theory. The angular dependence of the upper critical field shows a cusplike feature for fields close to the layers, consistent with decoupled layers.

Response of a two-level atom to a narrow-bandwidth squeezed-vacuum excitation

Using the coupled-system approach we calculate the optical spectra of the fluorescence and transmitted fields of a two-level atom driven by a squeezed vacuum of bandwidths smaller than the natural atomic linewidth. We find that in this regime of squeezing bandwidths the spectra exhibit unique features, such as a hole burning and a three-peak structure, which do not appear for a broadband excitation. We show that the features are unique to the quantum nature of the driving squeezed vacuum field and donor appear when the atom is driven by a classically squeezed field. We find that a quantum squeezed-vacuum field produces squeezing in the emitted fluorescence field which appears only in the squeezing spectrum while there is no squeezing in the total field. We also discuss a nonresonant excitation and find that depending on the squeezing bandwidth there is a peak or a hole in the spectrum at a frequency corresponding to a three-wave-mixing process. The hole appears only for a broadband excitation and results from the strong correlations between squeezed-vacuum photons.

Transport properties of strongly correlated metals: A dynamical mean-field approach

The temperature dependence of the transport properties of the metallic phase of a frustrated Hubbard model on the hypercubic lattice at half-filling is calculated. Dynamical mean-held theory, which maps the Hubbard model onto a single impurity,Anderson model that is solved self-consistently, and becomes exact in the limit of large dimensionality, is used. As the temperature increases there is a smooth crossover from coherent Fermi liquid excitations at low temperatures to incoherent excitations at high temperatures. This crossover leads to a nonmonotonic temperature dependence for the resistance, thermopower, and Hall coefficient, unlike in conventional metals. The resistance smoothly increases from a quadratic temperature dependence at low temperatures to large values which can exceed the Mott-Ioffe-Regel value ha/e(2) (where a is a lattice constant) associated with mean free paths less than a lattice constant. Further signatures of the thermal destruction of quasiparticle excitations are a peak in the thermopower and the absence of a Drude peak in the optical conductivity. The results presented here are relevant to a wide range of strongly correlated metals, including transition metal oxides, strontium ruthenates, and organic metals.

Finite element modelling of dissipative structures for nonequilibrium chemical reactions in fluid-saturated porous media

We use the finite element method to model and predict the dissipative structures of chemical species for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. In particular, we explore the conditions under which dissipative structures of the species may exist in the Brusselator type of nonequilibrium chemical reaction. Since this is the first time the finite element method and related strategies have been used to study the chemical instability problems in a fluid-saturated porous medium, it is essential to validate the method and strategies before they are put into application. For this purpose, we have rigorously derived the analytical solutions for dissipative structures of chemical species in a benchmark problem, which geometrically is a square. Comparison of the numerical solutions with the analytical ones demonstrates that the proposed numerical method and strategy are robust enough to solve chemical instability problems in a fluid-saturated porous medium. Finally, the related numerical results from two application examples indicate that both the regime and the magnitude of pore-fluid flow have significant effects on the nature of the dissipative structures that developed for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. The motivation for this study is that self-organization under conditions of pore-fluid flow in a porous medium is a potential mechanism of the orebody formation and mineralization in the upper crust of the Earth. (C) 2000 Elsevier Science S.A. All rights reserved.