Finite element modelling of three-dimensional convection problems in fluid-saturated porous media heated from below

We use the finite element method to model three-dimensional convective pore-fluid flow in fluid-saturated porous media when they are heated from below. In particular, we employ the particle-tracking technique to mimic the trajectories of particles in three-dimensional fluid flow problems. The related numerical results demonstrated that: (1) The progressive asymptotic approach procedure, which was previously developed for the finite element modelling of two-dimensional convective pore-fluid flow problems, is equally applicable to the finite element modelling of three-dimensional convective pore-fluid flow in fluid-saturated porous media heated from below. (2) The perturbation of gravity at different planes has a significant effect on the pattern of three-dimensional convective pore-fluid flow and therefore, may influence the pattern of orebody formation and mineralization in three-dimensional hydrothermal systems. Copyright (C) 2001 John Wiley & Sons, Ltd.

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.

Ischaemic preconditioning in the rat brain: a longitudinal magnetic resonance imaging (MRI) study

Ischaemic preconditioning in rats was studied using MRI. Ischaemic preconditioning was induced, using an intraluminal filament method, by 30 min middle cerebral artery occlusion (MCAO), and imaged 24 h later. The secondary insult of 100 min MCAO was induced 3 days following preconditioning and imaged 24 and 72 h later. Twenty four hours following ischaemic preconditioning most rats showed small sub-cortical hyperintense regions not seen in sham-preconditioned rats. Twenty-four hours and 72 h following the secondary insult preconditioned animals showed significantly smaller lesions (24 h = 112 +/- 31 mm(3), mean +/- standard error; 72 h = 80 +/- 35 mm(3)) which were confined to the striatum, than controls (24 h = 234 +/- 32 mm(3), p = 0.026; 72 h = 275 +/- 37 mm(3), p = 0.003). In addition during Lesion maturation from 24 to 72 h post-secondary MCAO, preconditioned rats displayed an average reduction in lesion size as measured by MRI whereas sham-preconditioned rats displayed increases in lesion size; this is the first report of such differential lesion volume evolution in cerebral ischaemic preconditioning. Copyright (C) 2001 John Wiley & Sons, Ltd.

Processing of urinary pheromones in Antechinus stuartii (Marsupialia: Dasyuridae): Functional magnetic resonance imaging of the brain

Little is known of the neural mechanisms of marsupial olfaction. However, functional magnetic resonance imaging (fMRI) has made it possible to visualize dynamic brain function in mammals without invasion. In this study, central processing of urinary pheromones was investigated in the brown antechinus, Antechinus stuartii, using fMRI. Images were obtained from 18 subjects (11 males, 7 females) in response to conspecific urinary olfactory stimuli. Significant indiscriminate activation occurred in the accessory olfactory bulb, entorhinal, frontal, and parietal cortices in response to both male and female urine. The paraventricular nucleus of hypothalamus, ventrolateral thalamic nucleus, and medial preoptic area were only activated in response to male urine. Results of this MRI study indicate that projections of accessory olfactory system are activated by chemo-sensory cues. Furthermore, it appears that, based on these experiments, urinary pheromones may act on the hypothalamo-pituitary-adrenocortical axis via the paraventricular nucleus of the hypothalamus and may play an important role in the unique life history pattern of A. stuartii. Finally, this study has demonstrated that fMRI may be a powerful tool for investigations of olfactory processes in mammals.

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.

Magnetic-field-induced superconductivity in layered organic molecular crystals with localized magnetic moments

The synthetic organic compound λ(BETS)2FeCl4 undergoes successive transitions from an antiferromagnetic insulator to a metal and then to a superconductor as a magnetic field is increased. We use a Hubbard-Kondo model to clarify the role of the Fe3+ magnetic ions in these phase transition. In the high-field regime, the magnetic field acting on the electron spins is compensated by the exchange field He due to the magnetic ions. This suggests that the field-induced superconducting state is the same as the zero-field superconducting state which occurs under pressure or when the Fe3+ ions are replaced by non-magnetic Ga3+ ions. We show how Hc can be extracted from the observed splitting of the Shybnikov-de Haas frequencies. Furthermore, we use this method of extracting He to predict the field range for field-induced superconductivity in other materials. We also show that at high fields the spin fluctuations of the localized spins are not important.

Effect of material anisotropy on the onset of convective flow in three-dimensional fluid-saturated faults

In order to investigate the effect of material anisotropy on convective instability of three-dimensional fluid-saturated faults, an exact analytical solution for the critical Rayleigh number of three-dimensional convective flow has been obtained. Using this critical Rayleigh number, effects of different permeability ratios and thermal conductivity ratios on convective instability of a vertically oriented three-dimensional fault have been examined in detail. It has been recognized that (1) if the fault material is isotropic in the horizontal direction, the horizontal to vertical permeability ratio has a significant effect on the critical Rayleigh number of the three-dimensional fault system, but the horizontal to vertical thermal conductivity ratio has little influence on the convective instability of the system, and (2) if the fault material is isotropic in the fault plane, the thermal conductivity ratio of the fault normal to plane has a considerable effect on the critical Rayleigh number of the three-dimensional fault system, but the effect of the permeability ratio of the fault normal to plane on the critical Rayleigh number of three-dimensional convective flow is negligible.

Effects of hot intrusions on pore-fluid flow and heat transfer in fluid-saturated rocks

We use the finite element method to solve coupled problems between pore-fluid flow and heat transfer in fluid-saturated porous rocks. In particular, we investigate the effects of both the hot pluton intrusion and topographically driven horizontal flow on the distributions of the pore-flow velocity and temperature in large-scale hydrothermal systems. Since general mineralization patterns are strongly dependent on distributions of both the pore-fluid velocity and temperature fields, the modern mineralization theory has been used to predict the general mineralization patterns in several realistic hydrothermal systems. The related numerical results have demonstrated that: (1) The existence of a hot intrusion can cause an increase in the maximum value of the pore-fluid velocity in the hydrothermal system. (2) The permeability of an intruded pluton is one of the sensitive parameters to control the pore-fluid flow, heat transfer and ore body formation in hydrothermal systems. (3) The maximum value of the pore-fluid velocity increases when the bottom temperature of the hydrothermal system is increased. (4) The topographically driven flow has significant effects on the pore-fluid flow, temperature distribution and precipitation pattern of minerals in hydrothermal systems. (5) The size of the computational domain may have some effects on the pore-fluid flow and heat transfer, indicating that the size of a hydrothermal system may affect the pore-fluid flow and heat transfer within the system. (C) 2003 Elsevier Science B.V. All rights reserved.

Convective instability of 3-D fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the convective instability of 3-D fluid-saturated geological fault zones when they are heated uniformly from below. In particular, we have derived exact analytical solutions for the critical Rayleigh numbers of different convective flow structures. Using these critical Rayleigh numbers, three interesting convective flow structures have been identified in a geological fault zone system. It has been recognized that the critical Rayleigh numbers of the system have a minimum value only for the fault zone of infinite length, in which the corresponding convective flow structure is a 2-D slender-circle flow. However, if the length of the fault zone is finite, the convective flow in the system must be 3-D. Even if the length of the fault zone is infinite, since the minimum critical Rayleigh number for the 2-D slender-circle flow structure is so close to that for the 3-D convective flow structure, the system may have almost the same chance to pick up the 3-D convective flow structures. Also, because the convection modes are so close for the 3-D convective flow structures, the convective flow may evolve into the 3-D finger-like structures, especially for the case of the fault thickness to height ratio approaching zero. This understanding demonstrates the beautiful aspects of the present analytical solution for the convective instability of 3-D geological fault zones, because the present analytical solution is valid for any value of the ratio of the fault height to thickness. Using the present analytical solution, the conditions, under which different convective flow structures may take place, can be easily determined.

Thermal Characteristics of Polybia scutellaris Nests (Hymenoptera: Vespidae) Using Computational Fluid Dynamics: A Possible Adaptation to Tropical Climates

Polybia scutellaris constructs huge nests characterized by numerous spinal projections on the surface. We investigated the thermal characteristics of P scutellaris nests in order to determine whether their nest temperature is homeothermically maintained and whether the spines play a role in the thermoregulation of the nests. In order to examine these hypotheses, we measured the nest temperature in a active nest and in an abandoned nest. The temperature in the active nest was almost stable at 27 degrees C, whereas that of the abandoned nest varied with changes in the ambient temperature, suggesting that nest temperature was maintained by the thermogenesis of colony individuals. In order to predict the thermal properties of the spines, a numerical simulation was employed. To construct a 3D-model of a P scutellaris nest, the nest architecture was simplified into an outer envelope and the surface spines, for both of which the initial temperature was set at 27 degrees C. The physical properties of the simulated nest were regarded to be those of wood since the nest of this species is constructed from plant materials. When the model was exposed to cool air (12 degrees C), the temperature was lower in the models with more spines. On the other hand, when the nest was heated (42 degrees C), the temperature increase was smaller in models with more spines. It is suggested that the spines act as a heat radiator, not as an insulator, against the changes in ambient temperature.

Simultaneous analysis of parabens in cosmetic products by stir bar sorptive extraction and liquid chromatography

A sensitive and precise stir bar sorptive extraction (SBSE) combined with LC (SBSE/LC) analysis is described for simultaneous determination of methyl, ethyl, propyl, and butyl parabens in commercial cosmetic products in agreement with the European Union Cosmetics Directive 76/768/EEC. Important factors in the optimization of SB SE efficiency are discussed, such as time and temperature of extraction, pH, and ionic strength of the sample, matrix effects, and liquid desorption conditions by different modes (magnetic stirring, ultrasonic). The LOQs of the SBSE/LC method ranged from 30 to 200 ng/mg, with linear response over a dynamic range, from the LOQ to 2.5 mu g/mg, with a coefficient of determination higher than 0.993. The interday precision of the SBSE/LC method presented a coefficient of variation lower than 5%. The effectiveness of the proposed method was proven for analysis of commercial cosmetic products such as body creams, antiperspirant creams, and sunscreens.

Cost effectiveness of adding magnetic resonance imaging to the usual management of suspected scaphoid fractures

Objective: To determine the cost effectiveness of a magnetic resonance imaging scan (MRI) within 5 days of injury compared with the usual management of occult scaphoid fracture. Methods: All patients with suspected scaphoid fractures in five hospitals were invited to participate in a randomised controlled trial of usual treatment with or without an MRI scan. Healthcare costs were compared, and a cost effectiveness analysis of the use of MRI in this scenario was performed. Results: Twenty eight of the 37 patients identified were randomised: 17 in the control group, 11 in the MRI group. The groups were similar at baseline and follow up in terms of number of scaphoid fractures, other injuries, pain, and function. Of the patients without fracture, the MRI group had significantly fewer days immobilised: a median of 3.0 (interquartile range 3.0-3.0) v 10.0 (7-12) in the control group (p = 0.006). The MRI group used fewer healthcare units (median 3.0, interquartile range 2.0-4.25) than the control group (5.0, 3.0-6.5) (p = 0.03 for the difference). However, the median cost of health care in the MRI group ($594.35 AUD, $551.35-667.23) was slightly higher than in the control group ($428.15, $124.40-702.65) (p = 0.19 for the difference). The mean incremental cost effectiveness ratio derived from this simulation was that MRI costs $44.37 per day saved from unnecessary immobilisation (95% confidence interval $4.29 to $101.02). An illustrative willingness to pay was calculated using a combination of the trials measure of the subjects' individual productivity losses and the average daily earnings. Conclusions: Use of MRI in the management of occult scaphoid fracture reduces the number of days of unnecessary immobilisation and use of healthcare units. Healthcare costs increased non-significantly in relation to the use of MRI in this setting. However, when productivity losses are considered, MRI may be considered cost effective, depending on the individual case.

Thermal analysis, nuclear magnetic resonance spectroscopy, and impedance spectroscopy of N,N-dimethyl-pyrrolidinium iodide: An ionic solid exhibiting rotator phases

N,N-dimethyl-pyrrolidinium iodide has been investigated using differential scanning calorimetry, nuclear magnetic resonance (NMR) spectroscopy, second moment calculations, and impedance spectroscopy. This pyrrolidinium salt exhibits two solid-solid phase transitions, one at 373 K having an entropy change, Delta S, of 38 J mol(-1) K-1 and one at 478 K having Delta S of 5.7 J mol(-1) K-1. The second moment calculations relate the lower temperature transition to a homogenization of the sample in terms of the mobility of the cations, while the high temperature phase transition is within the temperature region of isotropic tumbling of the cations. At higher temperatures a further decrease in the H-1 NMR linewidth is observed which is suggested to be due to diffusion of the cations. (C) 2005 American Institute of Physics.

Magnetic field effects on the conductivity of organic bipolar and unipolar devices at room temperature

Magnetic field effects on the conductivity of different types of organic devices: undoped and dye doped aluminium (III) 8-hydroxyquinoline (Alq(3))-based organic light emitting diodes (OLEDs), electron-only Alq(3)-based diodes, and a hole-only N,N`-diphenyl-N,N`-bis(1-naphthyl)1,1`-biphenyl-4,4`-diamine (alpha-NPD)-based diode were studied at room temperature. Only negative magnetoresistance (MR) was observed for the Alq(3)-based devices. The addition of a rubrene dye in Alq(3)-based OLEDs quenches the MR by a factor of 5. The alpha-NPD hole-only device showed only positive MR. Our results are discussed with respect to the actual models for MR in organic semiconductors. Our results are in good agreement with the bipolaron model. (C) 2009 Elsevier B.V. All rights reserved.