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.

H-reflex modulation during passive lengthening and shortening of the human triceps surae

1. The present study investigated the effects of lengthening and shortening actions on IT-reflex amplitude. H-reflexes were evoked in the soleus (SOL) and medial gastroenemius (MG) of human subject, during passive isometric, lengthening and shortening actions performed at angular velocities of 0, +/-2, +/-5 and +/- 15 deg s(-1). 2. H-reflex amplitude, in froth SOL and MG were significantly depressed during passive lengthening actions and facilitated during passive shortening actions, when compared with the isometric R-reflex amplitude. 3. Four experiments were performed in which the latencies front the onset of movement to delivery of the stimulus were altered. Passive H-reflex modulation during lengthening actions was found tee begin at latencies of less than 60 ms suggesting that this inhibition was due to peripheral and/or spinal mechanisms. 4. It is postulated that, the H-reflex modulation seen in the present study is related to the tunic discharge of muscle spindle afferents and the consequent effects of transmission within the la pathway. Inhibition of the H-reflex at less than 60 ms after the onset of muscle lengthening may he attributed to several mechanisms, which cannot be distinguished using the current protocol. These may include the inability to evoke volleys in la fibres that are refractory following muscle spindle discharge during; rapid muscle lengthening, a reduced probability of transmitter release front the presynaptic terminal (homosynaptic post.-activation depression) and presynaptic inhibition of la afferents from plantar flexor agonists. Short latency facilitation of the H-reflex may be attributed to temporal summation of excitatory postsynaptic potentials arising from muscle spindle afferents during rapid muscle lengthening. At longer latencies, presynaptic inhibition of Ia afferents cannot be excluded as a potential inhibitory mechanism.

Finite element modelling of heat transfer through permeable cracks in hydrothermal systems with upward throughflow

We use the finite element method to model the heat transfer phenomenon through permeable cracks in hydrothermal systems with upward throughflow. Since the finite element method is an approximate numerical method, the method must be validated before it is used to soh,e any new, kind of problem. However, the analytical solution, which can be used to validate the finite element method and other numerical methods, is rather limited in the literature, especially, for the problem considered here. Keeping this in mind, we have derived analytical solutions for the temperature distribution along the vertical axis of a crack in a fluid-saturated porous layer. After the finite element method is validated by comparing the numerical solution with the analytical solution for the same benchmark problem, it is used to investigate the pore-fluid flow and heat transfer in layered hydrothermal systems with vertical permeable cracks. The related analytical and numerical results have demonstrated that vertical cracks are effective and efficient members to transfer heat energy from the bottom section to the top section in hydrothermal systems with upward throughflow.

Computer simulations of coupled problems in geological and geochemical systems

The finite element method is used to simulate coupled problems, which describe the related physical and chemical processes of ore body formation and mineralization, in geological and geochemical systems. The main purpose of this paper is to illustrate some simulation results for different types of modelling problems in pore-fluid saturated rock masses. The aims of the simulation results presented in this paper are: (1) getting a better understanding of the processes and mechanisms of ore body formation and mineralization in the upper crust of the Earth; (2) demonstrating the usefulness and applicability of the finite element method in dealing with a wide range of coupled problems in geological and geochemical systems; (3) qualitatively establishing a set of showcase problems, against which any numerical method and computer package can be reasonably validated. (C) 2002 Published by Elsevier Science B.V.

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.

Analysis of steady-state heat transfer through mid-crustal vertical cracks with upward throughflow in hydrothermal systems

We conduct a theoretical analysis of steady-state heat transfer problems through mid-crustal vertical cracks with upward throughflow in hydrothermal systems. In particular, we derive analytical solutions for both the far field and near field of the system. In order to investigate the contribution of the forced advection to the total temperature of the system, two concepts, namely the critical Peclet number and the critical permeability of the system, have been presented and discussed in this paper. The analytical solution for the far field of the system indicates that if the pore-fluid pressure gradient in the crust is lithostatic, the critical permeability of the system can be used to determine whether or not the contribution of the forced advection to the total temperature of the system is negligible. Otherwise, the critical Peclet number should be used. For a crust of moderate thickness, the critical permeability is of the order of magnitude of 10(-20) m(2), under which heat conduction is the overwhelming mechanism to transfer heat energy, even though the pore-fluid pressure gradient in the crust is lithostatic. Furthermore, the lower bound analytical solution for the near field of the system demonstrates that the permeable vertical cracks in the middle crust can efficiently transfer heat energy from the lower crust to the upper crust of the Earth. Copyright (C) 2002 John Wiley Sons, Ltd.

Properties of the vertebrate skeletal muscle tubular system as a sealed compartment

Confocal imaging of impermeant fluorescent dyes trapped in the tubular (t-) system of skeletal muscle fibres of rat and cane toad was used to examine changes in the morphology of the t-system upon mechanical skinning, the time course of dye loss from the sealed t-systern in mechanically skinned fibres and the influence of rapid application and removal of glycerol on the morphology of the sealed t-system. In contrast to intact fibres, which have a t-systern open to the outside, the sealed t-systern of toad mechanically skinned fibres consistently displayed local swellings (vesicles). The occurrence of vesicles in the sealed t-system of rat-skinned fibres was infrequent. Application and removal of 200-400 mM glycerol to the sealed t-system did not produce any obvious changes in its morphology. The dyes fluo-3, fura-2 and Oregon green 488 were lost from the sealed t-system of toad fibres at different rates suggesting that the mechanism of organic anion transport across the tubular wall was not by indiscriminate bulk transport. The rate of fluo-3 and fura-2 loss from the sealed t-system of rat fibres was greater in rat than in toad fibres and could be explained by differences in surface area: volume ratio of the t-system in the two fibre types. Based on the results presented here and on other results from this laboratory, an explanation is given for the formation of numerous vesicles in toad-skinned fibres and lack of vesicle formation in rat-skinned fibres. This explanation can also help with better understanding the mechanism responsible for vacuole formation in intact fibres. (C) 2002 Elsevier Science Ltd. All rights reserved.

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.

An equivalent algorithm for simulating thermal effects of magma intrusion problems in porous rocks

An equivalent algorithm is proposed to simulate thermal effects of the magma intrusion in geological systems, which are composed of porous rocks. Based on the physical and mathematical equivalence, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with a physically equivalent heat source. From the analysis of an ideal solidification model, the physically equivalent heat source has been determined in this paper. The major advantage in using the proposed equivalent algorithm is that the fixed finite element mesh with a variable integration time step can be employed to simulate the thermal effect of the intruded magma solidification using the conventional finite element method. The related numerical results have demonstrated the correctness and usefulness of the proposed equivalent algorithm for simulating the thermal effect of the intruded magma solidification in geological systems. (C) 2003 Elsevier 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.

Theoretical and numerical analyses of convective instability in porous media with temperature-dependent viscosity

Exact analytical solutions of the critical Rayleigh numbers have been obtained for a hydrothermal system consisting of a horizontal porous layer with temperature-dependent viscosity. The boundary conditions considered are constant temperature and zero vertical Darcy velocity at both the top and bottom of the layer. Not only can the derived analytical solutions be readily used to examine the effect of the temperature-dependent viscosity on the temperature-gradient driven convective flow, but also they can be used to validate the numerical methods such as the finite-element method and finite-difference method for dealing with the same kind of problem. The related analytical and numerical results demonstrated that the temperature-dependent viscosity destabilizes the temperature-gradient driven convective flow and therefore, may affect the ore body formation and mineralization in the upper crust of the Earth. Copyright (C) 2003 John Wiley Sons, Ltd.

Active and passive smoking and the risk of subarachnoid hemorrhage - An international population-based case-control study

Background and Purpose - This study was undertaken to better clarify the risks associated with cigarette smoking and subarachnoid hemorrhage (SAH). Methods - The study included 432 incident cases of SAH frequency matched to 473 community SAH-free controls to determine dose-dependent associations of active and passive smoking ( at home) and smoking cessation with SAH. Results - Compared with never smokers not exposed to passive smoking, the adjusted odds ratio for SAH among current smokers was 5.0 (95% confidence interval [CI], 3.1 to 8.1); for past smokers, 1.2 ( 95% CI, 0.8 to 2.0); and for passive smokers, 0.9 ( 95% CI, 0.6 to 1.5). Current and lifetime exposures showed a clear dose-dependent effect, and risks appeared more prominent in women and for aneurysmal SAH. Approximately 1 in 3 cases of SAH could be attributed to current smoking, but risks decline quickly after smoking cessation, even among heavy smokers. Conclusions - A strong positive association was found between cigarette smoking and SAH, especially for aneurysmal SAH and women, which is virtually eliminated within a few years of smoking cessation. Large opportunities exist for preventing SAH through smoking avoidance and cessation programs.

Ethical preferences for influencing superiors: A 41-society study

With a 41-society sample of 9990 managers and professionals, we used hierarchical linear modeling to investigate the impact of both macro-level and micro-level predictors on subordinate influence ethics. While we found that both macro-level and micro-level predictors contributed to the model definition, we also found global agreement for a subordinate influence ethics hierarchy. Thus our findings provide evidence that developing a global model of subordinate ethics is possible, and should be based upon multiple criteria and multilevel variables. Journal of International Business Studies (2009) 40, 1022-1045. doi:10.1057/jibs.2008.109

Surgical management of lung cancer in Western Australia in 1996 and its outcomes

Background: All cases of lung cancer diagnosed in Western Australia in 1996 in which surgery was the primary treatment, were reviewed. Reported herein are the characteristics of the patients, the treatment outcomes and a comparison of the management undertaken with that recommended by international guidelines. Methods: All patients with a new diagnosis of lung cancer in Western Australia in the calendar year of 1996 were identified using two different population-based registration systems: the Western Australian (WA) Cancer Registry and the WA Hospital Morbidity Data System. A structured questionnaire on the diagnosis and management was completed for each case. Date of death was determined through the WA Cancer Registry. Results: Six hundred and sixty-eight patients with lung cancer were identified; 132 (20%) were treated with surgery. Lobectomy was the most frequently performed procedure (71%), followed by pneumonectomy (19%). Major complications affected 23% of patients. Postoperative mortality was 6% (3% lobectomy, 12% pneumonectomy). At 5 years the absolute survival was as follows for stage I, II, IIIA, IIIB, respectively: 51%, 45%, 12%, 5%. Conclusions: Investigations and choice of surgery in WA in 1996 reflect current international guidelines. The survival of patients with resectable lung cancer remains unsatisfactory.