Finite element modelling of rock alteration and metamorphic process in hydrothermal systems

We use the finite element method to simulate the rock alteration and metamorphic process in hydrothermal systems. In particular, we consider the fluid-rock interaction problems in pore-fluid saturated porous rocks. Since the fluid rock interaction takes place at the contact interface between the pore-fluid and solid minerals, it is governed by the chemical reaction which usually takes place very slowly at this contact interface, from the geochemical point of view. Due to the relative slowness of the rate of the chemical reaction to the velocity of the pore-fluid flow in the hydrothermal system to be considered, there exists a retardation zone, in which the conventional static theory in geochemistry does not hold true. Since this issue is often overlooked by some purely numerical modellers, it is emphasized in this paper. The related results from a typical rock alteration and metamorphic problem in a hydrothermal system have shown not only the detailed rock alteration and metamorphic process, but also the size of the retardation zone in the hydrothermal system. Copyright (C) 2001 John Wiley & Sons, Ltd.

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.

The VITATOPS (vitamins to prevent stroke) Trial: Rationale and design of an international, large, simple, randomised trial of homocysteine-lowering multivitamin therapy in patients with recent transient ischaemic attack or stroke

Background: Epidemiological studies suggest that raised plasma concentrations of total homocysteine (tHcy) may be a common, causal and treatable risk factor for atherothromboembolic ischaemic stroke. Although tHcy can be lowered effectively with small doses of folic acid, vitamin B-12 and vitamin B-6, it is not known whether lowering tHcy, by means of multivitamin therapy, can prevent stroke and other major atherothromboembolic vascular events. Purpose: To determine whether vitamin supplements (folic acid 2 mg, B-6 25 Mg, B-12 500 mug) reduce the risk of stroke, and other serious vascular events, in patients with recent stroke or transient ischaemic attacks of the brain or eye (TIA). Methods: An international, multi-centre, randomised, double-blind, placebo-controlled clinical trial. Results: As of November 2001, more than 1,400 patients have been randomised from 10 countries in four continents. Conclusion: VITATOPS aims to recruit and follow up 8,000 patients between 2000 and 2004, and provide a reliable estimate of the safety and effectiveness of dietary supplementation with folic acid, vitamin B-12, and vitamin B-6 in reducing recurrent serious vascular events among a wide range of patients with TIA and stroke. Copyright (C) 2002 S. Karger AG, Basel.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. 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.

Improving maximum walking distance in early peripheral arterial disease: Randomised controlled trial

The purpose of this study was to determine the impact of increased physical activity and cessation of smoking on the natural history of early peripheral arterial disease, We conducted a randomised controlled trial in Perth, Western Australia, involving 882 men with early peripheral arterial disease identified via population-based screening using the Edinburgh Claudication Questionnaire and the ankle:brachial index. Members of the control group (n = 441) received usual care from their general practitioner while members of the intervention group (n = 441) were allocated to a stop smoking and keep walking regime - a combined community-based intervention of cessation of smoking (where applicable) and increased physical activity. Postal follow-up occurred at two and 12 months post-entry into the trial. The main outcome of interest was maximum walking distance. There were no statistically significant differences in the characteristics of the intervention and usual care groups at recruitment. Follow-up information at two and 12 months was available for 85% and 84% of participants, respectively. At 12 months, more men allocated to the intervention group had improved their maximum walking distance (23% vs 15%; chi(2) = 9.74, df = 2, p = 0.008). In addition, more men in the intervention group reported walking more than three times per week for recreation (34% vs 25%, p = 0.01). Although not statistically significant, more men in the intervention group who were smokers when enrolled in the trial had stopped smoking (12% vs 8%, p = 0.43). It is concluded that referral of older patients with intermittent claudication to established physiotherapy programs in the community can increase levels of physical activity and reduce disability related to peripheral arterial disease. A combination of simple and safe interventions that are readily available in the community through physiotherapists and general practitioners has the potential to improve early peripheral arterial disease.

Scope for further analytical solutions for constant flux infiltration into a semi-infinite soil profile or redistribution in a finite soil profile

[1] We attempt to generate new solutions for the moisture content form of the one-dimensional Richards' [1931] equation using the Lisle [1992] equivalence mapping. This mapping is used as no more general set of transformations exists for mapping the one-dimensional Richards' equation into itself. Starting from a given solution, the mapping has the potential to generate an infinite number of new solutions for a series of nonlinear diffusivity and hydraulic conductivity functions. We first seek new analytical solutions satisfying Richards' equation subject to a constant flux surface boundary condition for a semi-infinite dry soil, starting with the Burgers model. The first iteration produces an existing solution, while subsequent iterations are shown to endlessly reproduce this same solution. Next, we briefly consider the problem of redistribution in a finite-length soil. In this case, Lisle's equivalence mapping is generalized to account for arbitrary initial conditions. As was the case for infiltration, however, it is found that new analytical solutions are not generated using the equivalence mapping, although existing solutions are recovered.

Participation in recreational physical activity: Why do socioeconomic groups differ?

This qualitative study explored how influences on recreational physical activity (RPA) were patterned by socioeconomic position. Face-to-face interviews were conducted with 10 males and 10 females in three socioeconomic groups (N = 60). Influences salient across all groups included previous opportunities, physical health, social assistance, safety, environmental aesthetics and urban design, physical and health benefits, and barriers of self-consciousness, low skill, and weather/time of year. Influences more salient to the high socioeconomic group included social benefits, achieving a balanced lifestyle, and the barrier of an unpredictable lifestyle. Influences more salient to the high and mid socioeconomic groups included efficacy, perceived need, activity demands, affiliation, emotional benefits, and the barrier of competing demands. Influences more salient to the low socioeconomic group included poor health and barriers of inconvenient access and low personal functioning. Data suggest that efforts to increase RPA in the population should include both general and socioeconomically targeted strategies.

Finite-difference time-domain-based studies of MRI pulsed field gradient-induced eddy currents inside the human body

In modern magnetic resonance imaging (MRI), patients are exposed to strong, rapidly switching magnetic gradient fields that, in extreme cases, may be able to elicit nerve stimulation. This paper presents theoretical investigations into the spatial distribution of induced current inside human tissues caused by pulsed z-gradient fields. A variety of gradient waveforms have been studied. The simulations are based on a new, high-definition, finite-difference time-domain method and a realistic inhomogeneous 10-mm resolution human body model with appropriate tissue parameters. it was found that the eddy current densities are affected not only by the pulse sequences but by many parameters such as the position of the body inside the gradient set, the local biological material properties and the geometry of the body. The discussion contains a comparison of these results with previous results found in the literature. This study and the new methods presented herein will help to further investigate the biological effects caused by the switched gradient fields in a MRI scan. (C) 2002 Wiley Periodicals, Inc.

Texture alignment in simple shear

We illustrate the flow behaviour of fluids with isotropic and anisotropic microstructure (internal length, layering with bending stiffness) by means of numerical simulations of silo discharge and flow alignment in simple shear. The Cosserat theory is used to provide an internal length in the constitutive model through bending stiffness to describe isotropic microstructure and this theory is coupled to a director theory to add specific orientation of grains to describe anisotropic microstructure. The numerical solution is based on an implicit form of the Material Point Method developed by Moresi et al. [1].

Mental rotation processes in mathematically gifted boys: An fMRI investigation

The functional brain organisation of mathematically gifted adolescents may be different from those of average mathematical ability. In this study we used fMRI to examine the neural circuitry that mediates the performance of mathematically gifted boys and average ability controls while engaged in mental rotation. Eight math gifted male adolescents and five average ability male adolescents were presented 18 control and 18 mental rotation trials in two separate blocks. Participants selected one of four test stimuli to match the target stimulus by pressing one of four fibreoptic buttons. The control task required a simple 'best match' for the target stimulus. EPI scans were acquired on a 3-T MR scanner and a fixed effects statistical analysis (SPM99) was used to identify areas of significant activation in the rotation tasks, for the two groups. The results indicate that during mental rotation both groups activate the parietal lobes bilaterally, though to different levels. Moreover, the math gifted are uniformly bilateral in their pattern of activation, and engage some anterior regions not found in those of average ability. These regions include bilateral prefrontal cortex and the right anterior cingulate, which may serve to heighten concentration, and to optimise the pre-planning of purposeful actions.

A spectral clustering algorithm for manufacturing cell formation

A graph clustering algorithm constructs groups of closely related parts and machines separately. After they are matched for the least intercell moves, a refining process runs on the initial cell formation to decrease the number of intercell moves. A simple modification of this main approach can deal with some practical constraints, such as the popular constraint of bounding the maximum number of machines in a cell. Our approach makes a big improvement in the computational time. More importantly, improvement is seen in the number of intercell moves when the computational results were compared with best known solutions from the literature. (C) 2009 Elsevier Ltd. All rights reserved.

Estimating the potential refolding yield of recombinant proteins expressed as inclusion bodies

Recombinant protein production in bacteria is efficient except that insoluble inclusion bodies form when some gene sequences are expressed. Such proteins must undergo renaturation, which is an inefficient process due to protein aggregation on dilution from concentrated denaturant. In this study, the protein-protein interactions of eight distinct inclusion-body proteins are quantified, in different solution conditions, by measurement of protein second virial coefficients (SVCs). Protein solubility is shown to decrease as the SVC is reduced (i.e., as protein interactions become more attractive). Plots of SVC versus denaturant concentration demonstrate two clear groupings of proteins: a more aggregative group and a group having higher SVC and better solubility. A correlation of the measured SVC with protein molecular weight and hydropathicity, that is able to predict which group each of the eight proteins falls into, is presented. The inclusion of additives known to inhibit aggregation during renaturation improves solubility and increases the SVC of both protein groups. Furthermore, an estimate of maximum refolding yield (or solubility) using high-performance liquid chromatography was obtained for each protein tested, under different environmental conditions, enabling a relationship between yield and SVC to be demonstrated. Combined, the results enable an approximate estimation of the maximum refolding yield that is attainable for each of the eight proteins examined, under a selected chemical environment. Although the correlations must be tested with a far larger set of protein sequences, this work represents a significant move beyond empirical approaches for optimizing renaturation conditions. The approach moves toward the ideal of predicting maximum refolding yield using simple bioinformatic metrics that can be estimated from the gene sequence. Such a capability could potentially screen, in silico, those sequences suitable for expression in bacteria from those that must be expressed in more complex hosts. (C) 2004 Wiley Periodicals, Inc.

Loss of Ergodicity in the Transition from Annealed to Quenched Disorder in a Finite Kinetic Ising Model

We consider a kinetic Ising model which represents a generic agent-based model for various types of socio-economic systems. We study the case of a finite (and not necessarily large) number of agents N as well as the asymptotic case when the number of agents tends to infinity. The main ingredient are individual decision thresholds which are either fixed over time (corresponding to quenched disorder in the Ising model, leading to nonlinear deterministic dynamics which are generically non-ergodic) or which may change randomly over time (corresponding to annealed disorder, leading to ergodic dynamics). We address the question how increasing the strength of annealed disorder relative to quenched disorder drives the system from non-ergodic behavior to ergodicity. Mathematically rigorous analysis provides an explicit and detailed picture for arbitrary realizations of the quenched initial thresholds, revealing an intriguing ""jumpy"" transition from non-ergodicity with many absorbing sets to ergodicity. For large N we find a critical strength of annealed randomness, above which the system becomes asymptotically ergodic. Our theoretical results suggests how to drive a system from an undesired socio-economic equilibrium (e. g. high level of corruption) to a desirable one (low level of corruption).