A mathematical model for estimating the extent of solute- and water-flux heterogeneity in multiple sample percolation experiments

Multiple sampling is widely used in vadose zone percolation experiments to investigate the extent in which soil structure heterogeneities influence the spatial and temporal distributions of water and solutes. In this note, a simple, robust, mathematical model, based on the beta-statistical distribution, is proposed as a method of quantifying the magnitude of heterogeneity in such experiments. The model relies on fitting two parameters, alpha and zeta to the cumulative elution curves generated in multiple-sample percolation experiments. The model does not require knowledge of the soil structure. A homogeneous or uniform distribution of a solute and/or soil-water is indicated by alpha = zeta = 1, Using these parameters, a heterogeneity index (HI) is defined as root 3 times the ratio of the standard deviation and mean. Uniform or homogeneous flow of water or solutes is indicated by HI = 1 and heterogeneity is indicated by HI > 1. A large value for this index may indicate preferential flow. The heterogeneity index relies only on knowledge of the elution curves generated from multiple sample percolation experiments and is, therefore, easily calculated. The index may also be used to describe and compare the differences in solute and soil-water percolation from different experiments. The use of this index is discussed for several different leaching experiments. (C) 1999 Elsevier Science B.V. All rights reserved.

Modeling the fluid-flow-induced stress and collapse in a dendritic network

An analytical approach to the stress development in the coherent dendritic network during solidification is proposed. Under the assumption that stresses are developed in the network as a result of the friction resisting shrinkage-induced interdendritic fluid flow, the model predicts the stresses in the solid. The calculations reflect the expected effects of postponed dendrite coherency, slower solidification conditions, and variations of eutectic volume fraction and shrinkage. Comparing the calculated stresses to the measured shear strength of equiaxed mushy zones shows that it is possible for the stresses to exceed the strength, thereby resulting in reorientation or collapse of the dendritic network.

The Heisenberg antiferromagnet on an anisotropic triangular lattice: linear spin-wave theory

We consider the effect of quantum spin fluctuations on the ground-state properties of the Heisenberg antiferromagnet on an anisotropic triangular lattice using linear spin-wave (LSW) theory. This model should describe the magnetic properties of the insulating phase of the kappa-(BEDT-TTF)(2)X family of superconducting molecular crystals. The ground-state energy, the staggered magnetization, magnon excitation spectra, and spin-wave velocities are computed as functions of the ratio of the antiferromagnetic exchange between the second and first neighbours, J(2)/J(1). We find that near J(2)/J(1) = 0.5, i.e., in the region where the classical spin configuration changes from a Neel-ordered phase to a spiral phase, the staggered magnetization vanishes, suggesting the possibility of a quantum disordered state. in this region, the quantum correction to the magnetization is large but finite. This is in contrast to the case for the frustrated Heisenberg model on a square lattice, for which the quantum correction diverges logarithmically at the transition from the Neel to the collinear phase. For large J(2)/J(1), the model becomes a set of chains with frustrated interchain coupling. For J(2) > 4J(1), the quantum correction to the magnetization, within LSW theory, becomes comparable to the classical magnetization, suggesting the possibility of a quantum disordered state. We show that, in this regime, the quantum fluctuations are much larger than for a set of weakly coupled chains with non-frustrated interchain coupling.

Hypoxia stimulates cerebral blood flow in the estuarine crocodile (Crocodylus porosus)

The effect of N-2 respiration on cerebral blood flow (CBF) velocity on the dorsal surface of cerebellum was examined in the estuarine crocodile, Crocodylus porosus, using epi-illumination microscopy. Twelve minutes of N-2 respiration resulted in a 126% increase in CBF velocity. N-2 respiration had no effect on blood pressure, indicating an underlying cerebral vasodilation. In addition, heart rate increased significantly. Systemic injections of aminophylline and the NO synthase (NOS) inhibitor nitro-L-arginine (L-NA) did not affect the hypoxia induced increase in CBF. We conclude that C. porosus responds to hypoxia with adenosine and nitric oxide (NO) independent cerebral vasodilation, and that this is likely to be a mechanism protecting the brain from energy deficiency during prolonged dives. (C) 1999 Elsevier Science Ireland Ltd. All rights reserved.

Targeting local tissues by transdermal application: Understanding drug physicochemical properties that best exploit protein binding and blood flow effects

The targeting of topically applied drug molecules into tissues below a site of application requires an understanding of the complex interrelationships between the drug, its formulation, the barrier properties of the skin, and the physiological processes occurring below the skin that are responsible for drug clearance from the site, tissue, and/or systemic distribution and eventual elimination. There is still a certain amount of controversy over the ability of topically applied drugs to penetrate into deeper tissues by diffusion or whether this occurs by redistribution in the systemic circulation. The major focus of our work in this area has been in determining how changes in drug structure and physicochemical properties, such as protein binding and lipophilicity, affect drug clearance into the local dermal microcirculation and lymphatics, as well as subsequent distribution into deeper tissues below an application site. The present study outlines our recent thinking on the drug molecule optimal physical attributes, in terms of plasma and tissue partitioning behaviour, that offer the greatest potential for deep tissue targeting. Drug Dev. Res. 46:309-315, 1999. (C) 1999 Wiley-Liss, Inc.

A novel two-chain proteinase inhibitor generated by circularization of a multidomain precursor protein

Female reproductive tissues of the ornamental tobacco amass high levels of serine proteinase inhibitors (PIs) for protection against pests and pathogens. These PIs are produced from a precursor protein composed of six repeats each with a protease reactive site. Here we show that proteolytic processing of the precursor generates five single-chain PIs and a remarkable two-chain inhibitor formed by disulfide-bond Linkage of Nand C-terminal peptide fragments. Surprisingly, PI precursors adopt this circular structure regardless of the number of inhibitor domains, suggesting this bracelet-like conformation is characteristic of the widespread potato inhibitor II (Pot II) protein family.

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.

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.

Two-wave-mixing induced self-focusing in an active photorefractive oscillator

We derive a nonlinear wave equation for a signal beam which is coupled to a pump beam by two-wave-mixing in a photorefractive crystal. This equation describes self-focusing of the signal beam. We compare two-wave-mixing induced spatial self-focusing of single-pass experiments in a diffusion-type photorefractive crystal and of a photorefractive oscillator using the same crystal. We observe that the nonlinear refractive index change in the oscillator is decreased while increasing resonator losses.

Submarine groundwater discharge and associated chemical input to a coastal sea

This paper presents a theoretical model of flow and chemical transport processes in subterranean estuaries (unconfined brackish groundwater aquifers at the ocean-land interface). The model shows that groundwater circulation and oscillating flow, caused by wave setup and tide, may constitute up to 96% of submarine groundwater discharge (SGWD) compared with 4% due to the net groundwater discharge. While these local flow processes do not change the total amount of land-derived chemical input to the ocean over a long period (e.g., yearly), they induce fluctuations of the chemical transfer rate as the aquifer undergoes saltwater intrusion. This may result in a substantial increase in chemical fluxes to the ocean over a short period (e.g., monthly and by a factor of 20 above the averaged level), imposing a possible threat to the marine environment. These results are essentially consistent with the experimental findings of Moore [1996] and have important implications for coastal resources management.

Increasing the efficiency of solute leaching: impacts of flow interruption with drainage of the "preferential flow paths"

Most soils contain preferential flow paths that can impact on solute mobility. Solutes can move rapidly down the preferential flow paths with high pore-water velocities, but can be held in the less permeable region of the soil matrix with low pore-water velocities, thereby reducing the efficiency of leaching. In this study, we conducted leaching experiments with interruption of the flow and drainage of the main flow paths to assess the efficiency of this type of leaching. We compared our experimental results to a simple analytical model, which predicts the influence of the variations in concentration gradients within a single spherical aggregate (SSA) surrounded by preferential flow paths on leaching. We used large (length: 300 mm, diameter: 216 mm) undisturbed field soil cores from two contrasting soil types. To carry out intermittent leaching experiments, the field soil cores were first saturated with tracer solution (CaBr2), and background solution (CaCl2) was applied to mimic a leaching event. The cores were then drained at 25- to 30-cm suction to empty the main flow paths to mimic a dry period during which solutes could redistribute within the undrained region. We also conducted continuous leaching experiments to assess the impact of the dry periods on the efficiency of leaching. The flow interruptions with drainage enhanced leaching by 10-20% for our soils, which was consistent with the model's prediction, given an optimised equivalent aggregate radius for each soil. This parameter quantifies the time scales that characterise diffusion within the undrained region of the soil, and allows us to calculate the duration of the leaching events and interruption periods that would lead to more efficient leaching. Application of these methodologies will aid development of strategies for improving management of chemicals in soils, needed in managing salts in soils, in improving fertiliser efficiency, and in reclaiming contaminated soils. (C) 2000 Elsevier Science B.V. All rights reserved.

Two-dimensional nonlinear dynamics of evanescent-wave guided atoms in a hollow fiber

We describe the classical and quantum two-dimensional nonlinear dynamics of large blue-detuned evanescent-wave guiding cold atoms in hollow fiber. We show that chaotic dynamics exists for classic dynamics, when the intensity of the beam is periodically modulated. The two-dimensional distributions of atoms in (x,y) plane are simulated. We show that the atoms will accumulate on several annular regions when the system enters a regime of global chaos. Our simulation shows that, when the atomic flux is very small, a similar distribution will be obtained if we detect the atomic distribution once each the modulation period and integrate the signals. For quantum dynamics, quantum collapses, and revivals appear. For periodically modulated optical potential, the variance of atomic position will be suppressed compared to the no modulation case. The atomic angular momentum will influence the evolution of wave function in two-dimensional quantum system of hollow fiber.

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.

A two-dimensional analytical solution of groundwater responses to tidal loading in an estuary and ocean

Previous studies on tidal dynamics of coastal aquifers have focussed on the inland propagation of oceanic tides in the cross-shore direction, a configuration that is essentially one-dimensional. Aquifers at natural coasts can also be influenced by tidal waves in nearby estuaries, resulting in a more complex behaviour of head fluctuations in the aquifers. We present an analytical solution to the two-dimensional depth-averaged groundwater flow equation for a semi-infinite aquifer subject to oscillating head conditions at the boundaries. The solution describes the tidal dynamics of a coastal aquifer that is adjacent to a cross-shore estuary. Both the effects of oceanic and estuarine tides on the aquifer are included in the solution. The analytical prediction of the head fluctuations is verified by comparison with numerical solutions computed using a standard finite-difference method. An essential feature of the present analytical solution is the interaction between the cross- and along-shore tidal waves in the aquifer area near the estuary's entry. As the distance from the estuary or coastline increases, the wave interaction is weakened and the aquifer response is reduced, respectively, to the one-dimensional solution for oceanic tides or the solution of Sun (Sun H. A two-dimensional analytical solution of groundwater response to tidal loading in an estuary, Water Resour Res 1997;33:1429-35) for two-dimensional non-interacting tidal waves. (C) 2000 Elsevier Science Ltd. All rights reserved.