Finite element modelling of temperature gradient driven rock alteration and mineralization in porous rock masses

We present finite element simulations of temperature gradient driven rock alteration and mineralization in fluid saturated porous rock masses. In particular, we explore the significance of production/annihilation terms in the mass balance equations and the dependence of the spatial patterns of rock alteration upon the ratio of the roll over time of large scale convection cells to the relaxation time of the chemical reactions. Special concepts such as the gradient reaction criterion or rock alteration index (RAI) are discussed in light of the present, more general theory. In order to validate the finite element simulation, we derive an analytical solution for the rock alteration index of a benchmark problem on a two-dimensional rectangular domain. Since the geometry and boundary conditions of the benchmark problem can be easily and exactly modelled, the analytical solution is also useful for validating other numerical methods, such as the finite difference method and the boundary element method, when they are used to dear with this kind of problem. Finally, the potential of the theory is illustrated by means of finite element studies related to coupled flow problems in materially homogeneous and inhomogeneous porous rock masses. (C) 1998 Elsevier Science S.A. All rights reserved.

Molecular development of the olfactory nerve pathway

There are, at least, two major questions concerning the molecular development of the olfactory nerve pathway. First, what are the molecular cues responsible for guiding axons from the nasal cavity to the olfactory bulb? Second, what is the molecular basis of axon targeting to specific glomeruli once axons reach the olfactory bulb? Studies in the primary olfactory pathway have focused on the role of the extracellular matrix and ensheathing cells in establishing an initial substrate for growth of pioneer axons between the periphery and brain. The primary axons also express a multitude of cell adhesion molecules that regulate fasciculation of axons and hence may play a role in fascicle formation in the olfactory nerve. Although the olfactory neuroepithelium principally consists of a morphologically homogeneous class of primary olfactory neurons, there are numerous subpopulations of olfactory neurons expressing chemically distinct phenotypes. In particular, numerous subpopulations have been characterized by expression of unique carbohydrate residues and olfactory receptor proteins. Some of these molecules have recently been implicated in axon guidance and targeting to specific glomeruli.

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.

Effects of beta-escin and saponin on the transverse-tubular system and sarcoplasmic reticulum membranes of rat and toad skeletal muscle

Mechanically skinned skeletal muscle fibres from rat and toad were exposed to the permeabilizing agents beta-escin and saponin. The effects of these agents on the sealed transverse tubular system (t-system) and sarcoplasmic reticulum (SR) were examined by looking at changes in the magnitude of the force responses to t-system depolarization, the time course of the fluorescence of fura-2 trapped in the sealed t-system, and changes in the magnitude of caffeine-induced contractures following SR loading with Ca2+ under defined conditions. In the presence of 2 mu g ml(-1) beta-escin and saponin, the response to t-system depolarization was not completely abolished, decreasing to a plateau, and a large proportion of fura-2 remained in the sealed t-system. At 10 mu g ml(-1), both agents abolished the ability of both rat and toad preparations to respond to t-system depolarization after 3 min of exposure, but a significant amount of fura-2 remained in sealed t-tubules even after exposure to 100 mu g ml(-1) beta-escin and saponin for 10 min. beta-Escin took longer than saponin to reduce the t-system depolarizations and fura-2 content of the sealed t-system to a similar level. The ability of the SR to load Ca2+ was reduced to a lower level after treatment with beta-escin than saponin. This direct effect on the SR occurred at much lower concentrations for rat (2 mu g ml(-1) beta-escin and 10 mu g ml(-1) saponin) than toad (10 mu g ml(-1) beta-escin and 150 mu g ml(-1) saponin). The reverse order in sensitivities to beta-escin and saponin of t-system and SR membranes indicates that the mechanisms of action of beta-escin and saponin are different in the two types of membrane. In conclusion, this study shows that: (1) beta-escin has a milder action on the surface membrane than saponin; (2) beta-escin is a more potent modifier of SR function; (3) simple permeabilization of membranes is not sufficient to explain the effects of beta-escin and saponin on muscle membranes; and (4) the t-system network within muscle fibres is not a homogeneous compartment.

Thylakoid membrane architecture

Confocal scanning laser microscopic observations were made on live chloroplasts in intact cells and on mechanically isolated, intact chloroplasts. Chlorophyll fluorescence was imaged to observe thylakoid membrane architecture. C-3 plant species studied included Spinacia oleracea L., Spathiphyllum sp. Schott, cv. 'Mauna Loa', and Pisum sativum L. C-4 plants were also investigated: Saccharum officinarum L., Sorghum bicolor L. Moench, Zea mays L. and Panicum miliaceum L. Some Spinacia chloroplasts were treated with 3-(3,4-dichlorophenyl)-1,1-dimethylurea (DCMU) to enhance or sodium dithionite (SD) to reduce the photosystem II fluorescence signal. Confocal microscopy images of C-3 chloroplasts differed from electron microscopy pictures because they showed discrete spots of bright fluorescence with black regions between them. There was no evidence of fluorescence from stroma thylakoids. The thylakoid membrane system at times appeared to be string-like, with brightly fluorescing grana lined up like beads. C-4 bundle sheath chloroplasts were imaged from three different types of C-4 plants. Saccharum and Sorghum bundle sheath chloroplasts showed homogeneous fluorescence and were much dimmer than mesophyll chloroplasts. Zea had rudimentary grana, and dim, homogeneous intergrana fluorescence was visualised. Panicum contained thylakoids similar in appearance and string-like arrangement to mesophyll chloroplasts. Isolated Pisum chloroplasts, treated with a drop of 5 mM MgCl2 showed a thylakoid membrane system which appeared to be unravelling. Spongy mesophyll chloroplasts of Spinacia treated with 5 mM sodium dithionite showed a granal thylakoid system with distinct regions of no fluorescence. A time-series experiment provided evidence of dynamic membrane rearrangements over a period of half an hour.

Asymmetric MRI magnet design using a hybrid numerical method

This paper describes a hybrid numerical method for the design of asymmetric magnetic resonance imaging magnet systems. The problem is formulated as a field synthesis and the desired current density on the surface of a cylinder is first calculated by solving a Fredholm equation of the first kind. Nonlinear optimization methods are then invoked to fit practical magnet coils to the desired current density. The field calculations are performed using a semi-analytical method. A new type of asymmetric magnet is proposed in this work. The asymmetric MRI magnet allows the diameter spherical imaging volume to be positioned close to one end of the magnet. The main advantages of making the magnet asymmetric include the potential to reduce the perception of claustrophobia for the patient, better access to the patient by attending physicians, and the potential for reduced peripheral nerve stimulation due to the gradient coil configuration. The results highlight that the method can be used to obtain an asymmetric MRI magnet structure and a very homogeneous magnetic field over the central imaging volume in clinical systems of approximately 1.2 m in length. Unshielded designs are the focus of this work. This method is flexible and may be applied to magnets of other geometries. (C) 1999 Academic Press.

A hybrid, inverse approach to the design of magnetic resonance imaging magnets

This paper describes a hybrid numerical method of an inverse approach to the design of compact magnetic resonance imaging magnets. The problem is formulated as a field synthesis and the desired current density on the surface of a cylinder is first calculated by solving a Fredholm equation of the first, kind. Nonlinear optimization methods are then invoked to fit practical magnet coils to the desired current density. The field calculations are performed using a semi-analytical method. The emphasis of this work is on the optimal design of short MRI magnets. Details of the hybrid numerical model are presented, and the model is used to investigate compact, symmetric MRI magnets as well as asymmetric magnets. The results highlight that the method can be used to obtain a compact MRI magnet structure and a very homogeneous magnetic field over the central imaging volume in clinical systems of approximately 1 m in length, significantly shorter than current designs. Viable asymmetric magnet designs, in which the edge of the homogeneous region is very close to one end of the magnet system are also presented. Unshielded designs are the focus of this work. This method is flexible and may be applied to magnets of other geometries. (C) 2000 American Association of Physicists in Medicine. [S0094-2405(00)00303-5].

Grazers and diggers: Exploitation competition and coexistence among foragers with different feeding strategies on a single resource

A mathematical model is presented that describes a system where two consumer species compete exploitatively for a single renewable resource. The resource is distributed in a patchy but homogeneous environment; that is, all patches are intrinsically identical. The two consumer species are referred to as diggers and grazers, where diggers deplete the resource within a patch to lower densities than grazers. We show that the two distinct feeding strategies can produce a heterogeneous resource distribution that enables their coexistence. Coexistence requires that grazers must either move faster than diggers between patches or convert the resources to population growth much more efficiently than diggers. The model shows that the functional form of resource renewal within a patch is also important for coexistence. These results contrast with theory that considers exploitation competition for a single resource when the resource is assumed to be well mixed throughout the system.

Habitat structure and egg distributions in the processionary caterpillar Ochrogaster lunifer: lessons for conservation and pest management

1. The spatial and temporal distribution of eggs laid by herbivorous insects is a crucial component of herbivore population stability, as it influences overall mortality within the population. Thus an ecologist studying populations of an endangered butterfly can do little to increase its numbers through habitat management without knowledge of its egg-laying patterns across individual host-plants under different habitat management regimes. At the other end of the spectrum, a knowledge of egg-laying behaviour can do much to control pest outbreaks by disrupting egg distributions that lead to rapid population growth. 2. The distribution of egg batches of the processionary caterpillar Ochrogaster lunifer on acacia trees was monitored in 21 habitats during 2 years in coastal Australia. The presence of egg batches on acacias was affected by host-tree 'quality' (tree size and foliar chemistry that led to increased caterpillar survival) and host-tree 'apparency' (the amount of vegetation surrounding host-trees). 3. In open homogeneous habitats, more egg batches were laid on high-quality trees, increasing potential population growth. In diverse mixed-species habitats, more egg batches were laid on low-quality highly apparent trees, reducing population growth and so reducing the potential for unstable population dynamics. The aggregation of batches on small apparent trees in diverse habitats led to outbreaks on these trees year after year, even when population levels were low, while site-wide outbreaks were rare. 4. These results predict that diverse habitats with mixed plant species should increase insect aggregation and increase population stability. In contrast, in open disturbed habitats or in regular plantations, where egg batches are more evenly distributed across high-quality hosts, populations should be more unstable, with site-wide outbreaks and extinctions being more common. 5. Mixed planting should be used on habitat regeneration sites to increase the population stability of immigrating or reintroduced insect species. Mixed planting also increases the diversity of resources, leading to higher herbivore species richness. With regard to the conservation of single species, different practices of habitat management will need to be employed depending on whether a project is concerned with methods of rapidly increasing the abundance of an endangered insect or concerned with the maintenance of a stable, established insect population that is perhaps endemic to an area. Suggestions for habitat management in these different cases are discussed. 6. Finally, intercropping can be highly effective in reducing pest outbreaks, although the economic gains of reduced pest attack may be outweighed by reduced crop yields in mixed-crop systems.

Water wave-driven seepage in marine sediments

The action of water waves moving over a porous seabed drives a seepage flux into and out of the marine sediments. The volume of fluid exchange per wave cycle may affect the rate of contaminant transport in the sediments. In this paper, the dynamic response of the seabed to ocean waves is treated analytically on the basis of pore-elastic theory applied to a porous seabed. The seabed is modelled as a semi-infinite, isotropic, homogeneous material. Most previous investigations on the wave-seabed interaction problem have assumed quasi-static conditions within the seabed, although dynamic behaviour often occurs in natural environments. Furthermore, wave pressures used in the previous approaches were obtained from conventional ocean wave theories: which are based on the assumption of an impermeable rigid seabed. By introducing a complex wave number, we derive a new wave dispersion equation, which includes the seabed characteristics (such as soil permeability, shear modulus, etc.). Based on the new closed-form analytical solution, the relative differences of the wave-induced seabed response under dynamic and quasi-static conditions are examined. The effects of wave and soil parameters on the seepage flux per wave cycle are also discussed in detail. (C) 2000 Elsevier Science Ltd. All rights reserved.

A method for the design of MRI radiofrequency coils based on triangular and pulse basis functions

This paper presents a numerical technique for the design of an RF coil for asymmetric magnetic resonance imaging (MRI) systems. The formulation is based on an inverse approach where the cylindrical surface currents are expressed in terms of a combination of sub-domain basis functions: triangular and pulse functions. With the homogeneous transverse magnetic field specified in a spherical region, a functional method is applied to obtain the unknown current coefficients. The current distribution is then transformed to a conductor pattern by use of a stream function technique. Preliminary MR images acquired using a prototype RF coil are presented and validate the design method. (C) 2002 Elsevier Science B.V. All rights reserved.

Finite element modelling of reactive fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins

We present the finite element simulations of reactive mineral carrying fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins. In particular we explore the mixing of reactive sulfide and sulfate fluids and the relevant patterns of mineralization for Load, zinc and iron minerals in the regime of temperature-gradient-driven convective flow. Since the mineralization and ore body formation may last quite a long period of time in a hydrothermal basin, it is commonly assumed that, in the geochemistry, the solutions of minerals are in an equilibrium state or near an equilibrium state. Therefore, the mineralization rate of a particular kind of mineral can be expressed as the product of the pore-fluid velocity and the equilibrium concentration of this particular kind of mineral Using the present mineralization rate of a mineral, the potential of the modern mineralization theory is illustrated by means of finite element studies related to reactive mineral-carrying fluids mixing problems in materially homogeneous and inhomogeneous porous rock basins.

On the theory of mixing-frequency electron spin-echo envelope modulation spectroscopy

It has recently been stated that the parametrization of the time variables in the one-dimensional (I-D) mixing-frequency electron spin-echo envelope modulation (MIF-ESEEM) experiment is incorrect and hence the wrong frequencies for correlated nuclear transitions are predicted. This paper is a direct response to such a claim, its purpose being to show that the parametrization in land 2-D MIF-ESEEM experiments possesses the same form as that used in other 4-pulse incrementation schemes and predicts the same correlation frequencies. We show that the parametrization represents a shearing transformation of the 2-D time-domain and relate the resulting frequency domain spectrum to the HYSCORE spectrum in terms of a skew-projection. It is emphasized that the parametrization of the time-domain variables may be chosen arbitrarily and affects neither the computation of the correct nuclear frequencies nor the resulting resolution. The usefulness or otherwise of the MIF parameters \gamma\ > 1 is addressed, together with the validity of the original claims of the authors with respect to resolution enhancement in cases of purely homogeneous and inhomogeneous broadening. Numerical simulations are provided to illustrate the main points.

Interpretation of negative values of the interaction parameter in the adsorption equation through the effects of surface layer heterogeneity

The problem of the negative values of the interaction parameter in the equation of Frumkin has been analyzed with respect to the adsorption of nonionic molecules on energetically homogeneous surface. For this purpose, the adsorption states of a homologue series of ethoxylated nonionic surfactants on air/water interface have been determined using four different models and literature data (surface tension isotherms). The results obtained with the Frumkin adsorption isotherm imply repulsion between the adsorbed species (corresponding to negative values of the interaction parameter), while the classical lattice theory for energetically homogeneous surface (e.g., water/air) admits attraction alone. It appears that this serious contradiction can be overcome by assuming heterogeneity in the adsorption layer, that is, effects of partial condensation (formation of aggregates) on the surface. Such a phenomenon is suggested in the Fainerman-Lucassen-Reynders-Miller (FLM) 'Aggregation model'. Despite the limitations of the latter model (e.g., monodispersity of the aggregates), we have been able to estimate the sign and the order of magnitude of Frumkin's interaction parameter and the range of the aggregation numbers of the surface species. (C) 2004 Elsevier B.V All rights reserved.

Unitary R-matrices for topological quantum computing

The main problem with current approaches to quantum computing is the difficulty of establishing and maintaining entanglement. A Topological Quantum Computer (TQC) aims to overcome this by using different physical processes that are topological in nature and which are less susceptible to disturbance by the environment. In a (2+1)-dimensional system, pseudoparticles called anyons have statistics that fall somewhere between bosons and fermions. The exchange of two anyons, an effect called braiding from knot theory, can occur in two different ways. The quantum states corresponding to the two elementary braids constitute a two-state system allowing the definition of a computational basis. Quantum gates can be built up from patterns of braids and for quantum computing it is essential that the operator describing the braiding-the R-matrix-be described by a unitary operator. The physics of anyonic systems is governed by quantum groups, in particular the quasi-triangular Hopf algebras obtained from finite groups by the application of the Drinfeld quantum double construction. Their representation theory has been described in detail by Gould and Tsohantjis, and in this review article we relate the work of Gould to TQC schemes, particularly that of Kauffman.