A numerical study of pore-fluid, thermal and mass flow iu fluid-saturated porous rock basins

We present a numerical methodology for the study of convective pore-fluid, thermal and mass flow in fluid-saturated porous rock basins. lit particular, we investigate the occurrence and distribution pattern of temperature gradient driven convective pore-fluid flow and hydrocarbon transport in the Australian North West Shelf basin. The related numerical results have demonstrated that: (1) The finite element method combined with the progressive asymptotic approach procedure is a useful tool for dealing with temperature gradient driven pore-fluid flow and mass transport in fluid-saturated hydrothermal basins; (2) Convective pore-fluid flow generally becomes focused in more permeable layers, especially when the layers are thick enough to accommodate the appropriate convective cells; (3) Large dislocation of strata has a significant influence off the distribution patterns of convective pore;fluid flow, thermal flow and hydrocarbon transport in the North West Shelf basin; (4) As a direct consequence of the formation of convective pore-fluid cells, the hydrocarbon concentration is highly localized in the range bounded by two major faults in the basin.

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].

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.

Visual function: How spiders find the right rock to crawl under

Animals that go on hunting expeditions face the problem of finding the way home at the end of the day. A group of hunting spiders has now been added to the list of animals that use the celestial pattern of polarized light as a compass for navigation. (C) 1999 Elsevier Science Ltd. All rights reserved.

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.

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.

Strain-dependent fluid flow defined through rock mass classification schemes

Strain-dependent hydraulic conductivities are uniquely defined by an environmental factor, representing applied normal and shear strains, combined with intrinsic material parameters representing mass and component deformation moduli, initial conductivities, and mass structure. The components representing mass moduli and structure are defined in terms of RQD (rock quality designation) and RMR (rock mass rating) to represent the response of a whole spectrum of rock masses, varying from highly fractured (crushed) rock to intact rock. These two empirical parameters determine the hydraulic response of a fractured medium to the induced-deformations The constitutive relations are verified against available published data and applied to study one-dimensional, strain-dependent fluid flow. Analytical results indicate that both normal and shear strains exert a significant influence on the processes of fluid flow and that the magnitude of this influence is regulated by the values of RQD and RMR.

Toxicity of biorational insecticides to Helicoverpa spp. (Lepidoptera : Noctuidae) and predators in cotton field

Field trials on upland cotton (Gossypium hirstum L.) during its reproductive phase were used to assess the toxicity of several biorational pesticides and chemicals to Helicoverpa armigera (Hubner) and H. puntigera Wallengren, as well as major predators at Dalby, Queensland, Australia. Moderate rate-dependent control was obtained in plots treated with neem (Azadirachta indica A. Juss) seed extract-azadirachtin (Aza) at rates of 30, 60 and 90 g/ha. Plots treated with Talstar EC (bifenthrin) applications achieved the best results, followed by treatment with alternation of chemicals (methomyl, bifenthrin, thiodicarb and endosulfan) and biorational insecticides (neem oil, azadirachtin and Bacillus thuringiensis kurstaki var. Berliner). Predators, including lady beetles, lacewings, spiders and predatory bugs, were insensitive to Aza, tooseendanin (Tsdn) and BT applications. In contrast, chemicals were very destructive of predators. All treatments provided some protection from infestation of H. armigera and H. puntigera. The effect of Aza on Helicoverpa spp. was reflected in a relatively higher yield of seed cotton harvested from Aza-treated plots compared with the control, but chemical control achieved significantly higher yields than any other treatment.

Prediction of rock alteration patterns: a potential tool in mineral exploration

Hydrothermal alteration of a quartz-K-feldspar rock is simulated numerically by coupling fluid flow and chemical reactions. Introduction of CO2 gas generates an acidic fluid and produces secondary quartz, muscovite and/or pyrophyllite at constant temperature and pressure of 300 degrees C and 200 MPa. The precipitation and/or dissolution of the secondary minerals is controlled by either mass-action relations or rate laws. In our simulations the mass of the primary elements are conserved and the mass-balance equations are solved sequentially using an implicit scheme in a finite-element code. The pore-fluid velocity is assumed to be constant. The change of rock volume due to the dissolution or precipitation of the minerals, which is directly related to their molar volume, is taken into account. Feedback into the rock porosity and the reaction rates is included in the model. The model produces zones of pyrophyllite quartz and muscovite due to the dissolution of K-feldspar. Our model simulates, in a simplified way, the acid-induced alteration assemblages observed in various guises in many significant mineral deposits. The particular aluminosilicate minerals produced in these experiments are associated with the gold deposits of the Witwatersrand Basin.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Heat flow for the Yang-Mills-Higgs field and the Hermitian Yang-Mills-Higgs metric

For a parameter lambda > 0, we study a type of vortex equations, which generalize the well-known Hermitian-Einstein equation, for a connection A and a section phi of a holomorphic vector bundle E over a Kahler manifold X. We establish a global existence of smooth solutions to heat flow for a self-dual Yang-Mills-Higgs field on E. Assuming the lambda -stability of (E, phi), we prove the existence of the Hermitian Yang-Mills-Higgs metric on the holomorphic bundle E by studying the limiting behaviour of the gauge flow.

A numerical study of flexural buckling of foliated rock slopes

The occurrence of foliated rock masses is common in mining environment. Methods employing continuum approximation in describing the deformation of such rock masses possess a clear advantage over methods where each rock layer and each inter-layer interface (joint) is explicitly modelled. In devising such a continuum model it is imperative that moment (couple) stresses and internal rotations associated with the bending of the rock layers be properly incorporated in the model formulation. Such an approach will lead to a Cosserat-type theory. In the present model, the behaviour of the intact rock layer is assumed to be linearly elastic and the joints are assumed to be elastic-perfectly plastic. Condition of slip at the interfaces are determined by a Mohr-Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformations. The model is incorporated into the finite element program AFENA and validated against an analytical solution of elementary buckling problems of a layered medium under gravity loading. A design chart suitable for assessing the stability of slopes in foliated rock masses against flexural buckling failure has been developed. The design chart is easy to use and provides a quick estimate of critical loading factors for slopes in foliated rock masses. It is shown that the model based on Euler's buckling theory as proposed by Cavers (Rock Mechanics and Rock Engineering 1981; 14:87-104) substantially overestimates the critical heights for a vertical slope and underestimates the same for sub-vertical slopes. Copyright (C) 2001 John Wiley & Sons, Ltd.

A new formation model for M32: a threshed early-type spiral galaxy?

The origin of M32, the closest compact elliptical galaxy (cE), is a long-standing puzzle of galaxy fort-nation in the Local Group. Our N-body/smoothed particle hydrodynamics simulations suggest a new scenario in which the strong tidal field of M31 can transform a spiral galaxy into a compact elliptical galaxy. As a low-luminosity spiral galaxy plunges into the central region of M31, most of the outer stellar and gaseous components of its disk are dramatically stripped as a result of M31's tidal field. The central bulge component, on the other hand, is just weakly influenced by the tidal field, owing to its compact configuration, and retains its morphology. M31's strong tidal field also induces rapid gas transfer to the central region, triggers a nuclear starburst, and consequently forms the central high-density and more metal-rich stellar populations with relatively young ages. Thus, in this scenario, M32 was previously the bulge of a spiral galaxy tidally interacting with M31 several gigayears ago. Furthermore, we suggest that cE's like M32 are rare, the result of both the rather narrow parameter space for tidal interactions that morphologically transform spiral galaxies into cE's and the very short timescale (less than a few times 10(9) yr) for cE's to be swallowed by their giant host galaxies (via dynamical friction) after their formation.