Finite element modelling of reactive mass transport problems in fluid-saturated porous media

We use the finite element method to solve reactive mass transport problems in fluid-saturated porous media. In particular, we discuss the mathematical expression of the chemical reaction terms involved in the mass transport equations for an isothermal, non-equilibrium chemical reaction. It has turned out that the Arrhenius law in chemistry is a good mathematical expression for such non-equilibrium chemical reactions especially from the computational point of view. Using the finite element method and the Arrhenius law, we investigate the distributions of PH (i.e. the concentration of H+) and the relevant reactive species in a groundwater system. Although the main focus of this study is on the contaminant transport problems in groundwater systems, the related numerical techniques and principles are equally applicable to the orebody formation problems in the geosciences. Copyright (C) 1999 John Wiley & Sons, Ltd.

Modelling solute transport in structured soils: Performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO3-, Cl-, PO43-) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed tone-way ANOVA, p < 0.001) indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale. (C) 2001 Elsevier Science Ltd. All rights reserved.

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.

Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. (c) 2004 Elsevier Ltd. All rights reserved.

Supersaturated solutions evaluated with an in vitro stratum corneum tape stripping technique

The flux of a compound across a membrane from any formulation, whether it contains penetration enhancers or not, is limited by its saturated solubility in the vehicle. Under such conditions the concentration of the permeant in the outer layers of the stratum corneum is also saturated. Consequently, when the permeation of a drug from a supersaturated solution leads to enhanced penetration, the concentration of the drug in the outer layers of the membrane is also supersaturated. Therefore, the stratum corneum may possess antinucleant properties which inhibit or retard the crystallisation process. In this study, the enhanced in vitro permeation of supersaturated solutions of piroxicam across human skin in diffusion cells was demonstrated. The amount of permeant in the stratum corneum was determined using a tape stripping technique. Supersaturated solutions up to four degrees of saturation were investigated which produced a linear relationship between the degree of saturation and the amount of piroxicam in the stratum corneum (R-2 = 0.970). Furthermore, the amount of piroxicam in the viable layers of the skin also increased with increasing degree of saturation. An analysis of the results suggested that enhanced penetration across human skin from supersaturated solutions of piroxicam may occur as a result of the antinucleating ability of the intercellular lipids of the stratum corneum. (C) 1997 Elsevier Science B.V.

Pharmacotherapy of the ion transport defect in cystic fibrosis

1. More than 1300 different mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) cause cystic fibrosis (CF), a disease characterized by deficient epithelial Cl- secretion and enhanced Na+ absorption. The clinical course of the disease is determined by the progressive lung disease. Thus, novel approaches in pharmacotherapy are based primarily on correction of the ion transport defect in the airways. 2. The current therapeutic strategies try to counteract the deficiency in Cl- secretion and the enhanced Na+ absorption. A number of compounds have been identified, such as genistein and xanthine derivatives, which directly activate mutant CFTR. Other compounds may activate alternative Ca2+-activated Cl- channels or basolateral K+ channels, which supply the driving force for Cl- secretion. Apart from that, Na+ channel blockers, such as phenamil and benzamil, are being explored, which counteract the hyperabsorption of NaCl in CF airways. 3. Clinical trials are under way using purinergic compounds such as the P2Y(2) receptor agonist INS365. Activation of P2Y(2) receptors has been found to both activate Cl- secretion and inhibit Na+ absorption. 4. The ultimate goal is to recover Cl- channel activity of mutant CFTR by either enhancing synthesis and expression of the protein or by activating silent CFTR Cl- channels. Strategies combining these drugs with compounds facilitating Cl- secretion and inhibiting Na+ absorption in vivo may have the best chance to counteract the ion transport defect in cystic fibrosis.

Solution techniques for transport problems involving steep concentration gradients: application to noncatalytic fluid solid reactions

Some efficient solution techniques for solving models of noncatalytic gas-solid and fluid-solid reactions are presented. These models include those with non-constant diffusivities for which the formulation reduces to that of a convection-diffusion problem. A singular perturbation problem results for such models in the presence of a large Thiele modulus, for which the classical numerical methods can present difficulties. For the convection-diffusion like case, the time-dependent partial differential equations are transformed by a semi-discrete Petrov-Galerkin finite element method into a system of ordinary differential equations of the initial-value type that can be readily solved. In the presence of a constant diffusivity, in slab geometry the convection-like terms are absent, and the combination of a fitted mesh finite difference method with a predictor-corrector method is used to solve the problem. Both the methods are found to converge, and general reaction rate forms can be treated. These methods are simple and highly efficient for arbitrary particle geometry and parameters, including a large Thiele modulus. (C) 2001 Elsevier Science Ltd. All rights reserved.

Influence of seaward boundary condition on contaminant transport in unconfined coastal aquifers

Contaminant transport in coastal aquifers is complicated partly due to the conditions at the seaward boundary including seawater intrusion and tidal variations of sea level. Their inclusion in modelling this system will be computationally expensive. Therefore, it will be instructive to investigate the consequence of simplifying the seaward boundary condition by neglecting the seawater density and tidal variations in numerical predictions of contaminant transport in this zone. This paper presents a comparison of numerical predictions for a simplified seaward boundary condition with experimental results for a corresponding realistic one including a saltwater interface and tidal variations. Different densities for contaminants are considered. The comparison suggests that the neglect of the seawater intrusion and tidal variations does not affect noticeably the overall migration rate of the plume before it reaches the saltwater interface. However, numerical prediction shows that a more dense contaminant travels further seaward and part of the solute mass exits under the sea if the seawater density is not included. This is not consistent with the experimental result, which shows that the contaminant travels upwards to the shoreline along the saltwater interface. Neglect of seawater density, therefore, will result in an underestimation of the exit rate of solute mass around the coastline and fictitious migration paths under the seabed. For a less dense contaminant, neglect of seawater density has little effect on numerical prediction of migration paths. (C) 2001 Elsevier Science B.V. All rights reserved.

Cloning and gastrointestinal expression of rat hephaestin: relationship to other iron transport proteins

The membrane-bound ceruloplasmin homolog hephaestin plays a critical role in intestinal iron absorption. The aims of this study were to clone the rat hephaestin gene and to examine its expression in the gastrointestinal tract in relation to other genes encoding iron transport proteins. The rat hephaestin gene was isolated from intestinal mRNA and was found to encode a protein 96% identical to mouse hephaestin. Analysis by ribonuclease protection assay and Western blotting showed that hephaestin was expressed at high levels throughout the small intestine and colon. Immunofluorescence localized the hephaestin protein to the mature villus enterocytes with little or no expression in the crypts. Variations in iron status had a small but nonsignificant effect on hephaestin expression in the duodenum. The high sequence conservation between rat and mouse hephaestin is consistent with this protein playing a central role in intestinal iron absorption, although its precise function remains to be determined.

Influence of nano-structure on electrolytic properties in CeO2 based system

Doped ceria (CeO2) compounds are fluorite type oxides that show oxygen ionic conductivity higher than yttria stabilized zirconia, in oxidizing atmosphere. In order to improve the conductivity, the effective index was suggested to maximize the oxygen ionic conductivity in doped CeO2 based oxides. In addition, the true microstructure of doped CeO2 was observed at atomic scale for conclusion of conduction mechanism. Doped CeO2 had small domains (10-50 nm) with ordered structure in a grain. It is found that the electrolytic properties strongly depended on the nano-structural feature at atomic scale in doped CeO2 electrolyte.

The effect of soft tissue deloading tape on thoracic spine pressure pain thresholds in asymptomatic subjects

The application of tape to deload soft tissue is used in the management of thoracic spine pain. A reported clinical feature of this treatment is reduced tenderness of the spine during postero-anterior mobilizations. A randomized, single blind, placebo controlled, repeated measures design study was employed to investigate the effects of deloading tape on pressure pain threshold measurements at the level of the T7 spinous process in an asymptomatic group of 24 subjects. Pressure pain thresholds were assessed prior to and following the application of deloading tape, placebo sham tape and no-tape control conditions. All subjects received all three conditions in a randomized order on three separate days. Differences between the pre- and post-measurements were used as indicators of change in a subject's pressure pain threshold. No significant change in pressure pain threshold measurements was found between conditions. In summary, this study demonstrated that deloading tape applied to the level of the T7 spinous process did not significantly change pressure pain threshold measurements in asymptomatic subjects, raising the possibility that any pain relieving effect may well be conditional upon pain being present. (C) 2002 Elsevier Science Ltd. All rights reserved.

Numerical modelling of contaminant transport in coastal aquifers

This paper employs a two-dimensional variable density flow and transport model to investigate the transport of a dense contaminant plume in an unconfined coastal aquifer. Experimental results are also presented to show the contaminant plume in a freshwater-seawater flow system. Both the numerical and experimental results suggest that the neglect of the seawater interface does not noticeably affect the horizontal migration rate of the plume before it reaches the interface. However, the contaminant will travel further seaward and part of the solute mass will exit under the sea if the higher seawater density is not included. If the seawater density is included, the contaminant will travel upwards towards the beach along the freshwater-saltwater interface as shown experimentally. Neglect of seawater density, therefore, will result in an underestimate of solute mass rate exiting around the coastline. (C) 2002 IMACS. Published by Elsevier Science B.V. All rights reserved.

Synaptic vesicle transport and synaptic membrane transporter sites in excitatory amino acid nerve terminals in Alzheimer disease

The apparent L-[H-3]glutamate uptake rate (v') was measured in synaptic vesicles isolated from cerebral cortex synaptosomes prepared from autopsied Alzheimer and non-Alzheimer dementia cases, and age-matched controls. The initial synaptosome preparations exhibited similar densities of D-[H-3]aspartate membrane binding sites (B-MAX values) in the three groups. In control brain the temporal cortex D-[H-3]aspartate B-MAX was 132% of that in motor cortex, parallel with the L- [H-3]glutamate v' values (temporal = 139% of motor; NS). Unlike D- [H-3]aspartate B-MAX values, L- [H-3]glutamate v' values were markedly and selectively lower in Alzheimer brain preparations than in controls, particularly in temporal cortex. The difference could not be attributed to differential effects of autopsy interval or age at death. Non-Alzheimer dementia cases resembled controls. The selective loss of vesicular glutamate transport is consistent with a dysfunction in the recycling of transmitter glutamate.

Abnormal extracellular matrix protein transport associated with increased apoptosis of vascular smooth muscle cells in Marfan syndrome and bicuspid aortic valve thoracic aortic aneurysm

Background - Marfan syndrome (MS) is a genetic disorder caused by a mutation in the fibrillin gene FBN1. Bicuspid aortic valve (BAV) is a congenital heart malformation of unknown cause. Both conditions are associated with ascending aortic aneurysm and premature death. This study examined the relationship among the secretion of extracellular matrix proteins fibrillin, fibronectin, tenascin, and vascular smooth muscle cell (VSMC) apoptosis. The role of matrix metalloproteinase (MMP)- 2 in VSMC apoptosis was studied in MS aneurysm. Methods and Results - Aneurysm tissue was obtained from patients undergoing surgery ( MS: 4 M, 1 F, age 27 - 45 years; BAV: 3 M, 2 F, age 28 - 65 years). Normal aorta from subjects with nonaneurysm disease was also collected ( 4 M, 1 F, age 23 - 93 years). MS and BAV aneurysm histology showed areas of cystic medial necrosis (CMN) without inflammatory infiltrate. Immunohistochemical study of cultured MS and BAV VSMC showed intracellular accumulation and reduction of extracellular distribution of fibrillin, fibronectin, and tenascin. Western blot showed no increase in expression of fibrillin, fibronectin, or tenascin in MS or BAV VSMC and increased expression of MMP-2 in MS VSMCs. There was 4-fold increase in loss of cultured VSMC incubated in serum-free medium for 24 hours in both MS ( 27 +/- 8%) and BAV ( 32 +/- 14%) compared with control ( 7 +/- 5%). Conclusions - In MS and BAV there is alteration in both the amount and quality of secreted proteins and an increased degree of VSMC apoptosis. Up-regulation of MMP-2 might play a role in VSMC apoptosis in MS VSMC. The findings suggest the presence of a fundamental cellular abnormality in BAV thoracic aorta, possibly of genetic origin.