Dissolution-Precipitation in Porous Media: Experiments and Modelling

The evolution of porosity due to dissolution/precipitation processes of minerals and the associated change of transport parameters are of major interest for natural geological environments and engineered underground structures. We designed a reproducible and fast to conduct 2D experiment, which is flexible enough to investigate several process couplings implemented in the numerical code OpenGeosys-GEM (OGS-GEM). We investigated advective-diffusive transport of solutes, effect of liquid phase density on advective transport, and kinetically controlled dissolution/precipitation reactions causing porosity changes. In addition, the system allowed to investigate the influence of microscopic (pore scale) processes on macroscopic (continuum scale) transport. A Plexiglas tank of dimension 10 × 10 cm was filled with a 1 cm thick reactive layer consisting of a bimodal grain size distribution of celestite (SrSO4) crystals, sandwiched between two layers of sand. A barium chloride solution was injected into the tank causing an asymmetric flow field to develop. As the barium chloride reached the celestite region, dissolution of celestite was initiated and barite precipitated. Due to the higher molar volume of barite, its precipitation caused a porosity decrease and thus also a decrease in the permeability of the porous medium. The change of flow in space and time was observed via injection of conservative tracers and analysis of effluents. In addition, an extensive post-mortem analysis of the reacted medium was conducted. We could successfully model the flow (with and without fluid density effects) and the transport of conservative tracers with a (continuum scale) reactive transport model. The prediction of the reactive experiments initially failed. Only the inclusion of information from post-mortem analysis gave a satisfactory match for the case where the flow field changed due to dissolution/precipitation reactions. We concentrated on the refinement of post-mortem analysis and the investigation of the dissolution/precipitation mechanisms at the pore scale. Our analytical techniques combined scanning electron microscopy (SEM) and synchrotron X-ray micro-diffraction/micro-fluorescence performed at the XAS beamline (Swiss Light Source). The newly formed phases include an epitaxial growth of barite micro-crystals on large celestite crystals (epitaxial growth) and a nano-crystalline barite phase (resulting from the dissolution of small celestite crystals) with residues of celestite crystals in the pore interstices. Classical nucleation theory, using well-established and estimated parameters describing barite precipitation, was applied to explain the mineralogical changes occurring in our system. Our pore scale investigation showed limits of the continuum scale reactive transport model. Although kinetic effects were implemented by fixing two distinct rates for the dissolution of large and small celestite crystals, instantaneous precipitation of barite was assumed as soon as oversaturation occurred. Precipitation kinetics, passivation of large celestite crystals and metastability of supersaturated solutions, i.e. the conditions under which nucleation cannot occur despite high supersaturation, were neglected. These results will be used to develop a pore scale model that describes precipitation and dissolution of crystals at the pore scale for various transport and chemical conditions. Pore scale modelling can be used to parameterize constitutive equations to introduce pore-scale corrections into macroscopic (continuum) reactive transport models. Microscopic understanding of the system is fundamental for modelling from the pore to the continuum scale.

The influence of saliva on the dissolution of calcium fluoride after application of different fluoride gels in vitro

OBJECTIVE: To determine the formation and dissolution of calcium fluoride on the enamel surface after application of two fluoride gel-saliva mixtures. METHOD AND MATERIALS: From each of 80 bovine incisors, two enamel specimens were prepared and subjected to two different treatment procedures. In group 1, 80 specimens were treated with a mixture of an amine fluoride gel (1.25% F-; pH 5.2; 5 minutes) and human saliva. In group 2, 80 enamel blocks were subjected to a mixture of sodium fluoride gel (1.25% F; pH 5.5; 5 minutes) and human saliva. Subsequent to fluoride treatment, 40 specimens from each group were stored in human saliva and sterile water, respectively. Ten specimens were removed after each of 1 hour, 24 hours, 2 days, and 5 days and analyzed according to potassium hydroxide-soluble fluoride. RESULTS: Application of amine fluoride gel resulted in a higher amount of potassium hydroxide-soluble fluoride than did sodium fluoride gel 1 hour after application. Saliva exerted an inhibitory effect according to the dissolution rate of calcium fluoride. However, after 5 days, more than 90% of the precipitated calcium fluoride was dissolved in the amine fluoride group, and almost all potassium hydroxide-soluble fluoride was lost in the sodium fluoride group. Calcium fluoride apparently dissolves rapidly, even at almost neutral pH. CONCLUSION: Considering the limitations of an in vitro study, it is concluded that highly concentrated fluoride gels should be applied at an adequate frequency to reestablish a calcium fluoride-like layer.

Effects of pH and acid concentration on erosive dissolution of enamel, dentine, and compressed hydroxyapatite

The aims of this study were to determine the effects of pH and acid concentration on the dissolution of enamel, dentine, and compressed hydroxyapatite (HA) in citric acid solutions (15.6 and 52.1 mmol l(-1) ; pH 2.45, 3.2, and 3.9), using a pH-stat system. After an initial adjustment period, the dissolution rates of enamel and HA were constant, while that of dentine decreased with time. The dissolution rate increased as the pH decreased, and this was most marked for enamel. To compare substrates, the rate of mineral dissolution was normalized to the area occupied by mineral at the specimen surface. For a given acid concentration, the normalized dissolution rate of HA was always less than that for either dentine or enamel. The dissolution rate for dentine mineral was similar to that for enamel at pH 2.45 and greater at pH 3.2 and pH 3.9. The concentration of acid significantly affected the enamel dissolution rate at pH 2.45 and pH 3.2, but not at pH 3.9, and did not significantly affect the dissolution rates of dentine or HA at any pH. The variation in response of the dissolution rate to acid concentration/buffer capacity with respect to pH and tissue type might complicate attempts to predict erosive potential from solution composition.

Biogeochemical processes in a clay formation in situ experiment: Part F - Reactive transport modelling

Reactive transport modelling was used to simulate solute transport, thermodynamic reactions, ion exchange and biodegradation in the Porewater Chemistry (PC) experiment at the Mont Terri Rock Laboratory. Simulations show that the most important chemical processes controlling the fluid composition within the borehole and the surrounding formation during the experiment are ion exchange, biodegradation and dissolution/precipitation reactions involving pyrite and carbonate minerals. In contrast, thermodynamic mineral dissolution/precipitation reactions involving alumo-silicate minerals have little impact on the fluid composition on the time-scale of the experiment. With the accurate description of the initial chemical condition in the formation in combination with kinetic formulations describing the different stages of bacterial activities, it has been possible to reproduce the evolution of important system parameters, such as the pH, redox potential, total organic C. dissolved inorganic C and SO(4) concentration. Leaching of glycerol from the pH-electrode may be the primary source of organic material that initiated bacterial growth, which caused the chemical perturbation in the borehole. Results from these simulations are consistent with data from the over-coring and demonstrate that the Opalinus Clay has a high buffering capacity in terms of chemical perturbations caused by bacterial activity. This buffering capacity can be attributed to the carbonate system as well as to the reactivity of clay surfaces.

Kinetic assessment of persistent halogenated xenobiotics in cell culture models: comparison of mono- and poly-halogenated compounds

We evaluated the suitability of single and multiple cell type cultures as model systems to characterise cellular kinetics of highly lipophilic compounds with potential ecotoxicological impact. Confluent mono-layers of human skin fibroblasts, rat astrocytoma C6 cells, non-differentiated and differentiated mouse 3T3 cells were kept in culture medium supplemented with 10% foetal calf serum. For competitive uptake experiments up to four different cell types, grown on glass sectors, were exposed for 3h to (14)C-labelled model compounds, dissolved either in organic solvents or incorporated into unilamellar lecithin liposomes. Bromo-, or chloro-benzenes, decabromodiphenylether (DBP), and dichlorodiphenyl ethylene (DDE) were tested in rather high concentration of 20 microM. Cellular toxicity was low. Compound levels were related to protein, DNA, and triglyceride contents. Cellular uptake was fast and dependent on physico-chemical properties of the compounds (lipophilicity, molecular size), formulation, and cell type. Mono-halogenated benzenes showed low and similar uptake levels (=low accumulation compounds). DBP and DDE showed much higher cellular accumulations (=high accumulation compounds) except for DBP in 3T3 cells. Uptake from liposomal formulations was mostly higher than if compounds were dissolved in organic solvents. The extent of uptake correlated with the cellular content of triglycerides, except for DBP. Uptake competition between different cell types was studied in a sectorial multi-cell culture model. For low accumulation compounds negligible differences were found among C6 cells and fibroblasts. Uptake of DDE was slightly and that of DBP highly increased in fibroblasts. Well-defined cell culture systems, especially the sectorial model, are appropriate to screen for bioaccumulation and cytotoxicity of (unknown) chemical entities in vitro.

Efficacy of telavancin against penicillin-resistant pneumococci and Staphylococcus aureus in a rabbit meningitis model and determination of kinetic parameters

The penetration of telavancin was 2% into inflamed meninges and ca. 1 per thousand into noninflamed meninges after two intravenous injections (30 mg/kg of body weight). In experimental meningitis, telavancin was significantly superior to vancomycin combined with ceftriaxone against a penicillin-resistant pneumococcal strain. Against a methicillin-sensitive staphylococcal strain, telavancin was slightly but not significantly superior to vancomycin.

Substrate binding tunes conformational flexibility and kinetic stability of an amino acid antiporter

We used single molecule dynamic force spectroscopy to unfold individual serine/threonine antiporters SteT from Bacillus subtilis. The unfolding force patterns revealed interactions and energy barriers that stabilized structural segments of SteT. Substrate binding did not establish strong localized interactions but appeared to be facilitated by the formation of weak interactions with several structural segments. Upon substrate binding, all energy barriers of the antiporter changed thereby describing the transition from brittle mechanical properties of SteT in the unbound state to structurally flexible conformations in the substrate-bound state. The lifetime of the unbound state was much shorter than that of the substrate-bound state. This leads to the conclusion that the unbound state of SteT shows a reduced conformational flexibility to facilitate specific substrate binding and a reduced kinetic stability to enable rapid switching to the bound state. In contrast, the bound state of SteT showed an increased conformational flexibility and kinetic stability such as required to enable transport of substrate across the cell membrane. This result supports the working model of antiporters in which alternate substrate access from one to the other membrane surface occurs in the substrate-bound state.

Modeling the marine aragonite cycle: changes under rising carbon dioxide and its role in shallow water CaCO3 dissolution

The marine aragonite cycle has been included in the global biogeochemical model PISCES to study the role of aragonite in shallow water CaCO3 dissolution. Aragonite production is parameterized as a function of mesozooplankton biomass and aragonite saturation state of ambient waters. Observation-based estimates of marine carbonate production and dissolution are well reproduced by the model and about 60% of the combined CaCO3 water column dissolution from aragonite and calcite is simulated above 2000 m. In contrast, a calcite-only version yields a much smaller fraction. This suggests that the aragonite cycle should be included in models for a realistic representation of CaCO3 dissolution and alkalinity. For the SRES A2 CO2 scenario, production rates of aragonite are projected to notably decrease after 2050. By the end of this century, global aragonite production is reduced by 29% and total CaCO3 production by 19% relative to pre-industrial. Geographically, the effect from increasing atmospheric CO2, and the subsequent reduction in saturation state, is largest in the subpolar and polar areas where the modeled aragonite production is projected to decrease by 65% until 2100.

Effects of buffering properties and undissociated acid concentration on dissolution of dental enamel in relation to pH and acid type

To quantify the relationships between buffering properties and acid erosion and hence improve models of erosive potential of acidic drinks, a pH-stat was used to measure the rate of enamel dissolution in solutions of citric, malic and lactic acids, with pH 2.4-3.6 and with acid concentrations adjusted to give buffer capacities (β) of 2-40 (mmol·l(-1))·pH(-1) for each pH. The corresponding undissociated acid concentrations, [HA], and titratable acidity to pH 5.5 (TA5.5) were calculated. In relation to β, the dissolution rate and the strength of response to β varied with acid type (lactic > malic ≥ citric) and decreased as pH increased. The patterns of variation of the dissolution rate with TA5.5 were qualitatively similar to those for β, except that increasing pH above 2.8 had less effect on dissolution in citric and malic acids and none on dissolution in lactic acid. Variations of the dissolution rate with [HA] showed no systematic dependence on acid type but some dependence on pH. The results suggest that [HA], rather than buffering per se, is a major rate-controlling factor, probably owing to the importance of undissociated acid as a readily diffusible source of H(+) ions in maintaining near-surface dissolution within the softened layer of enamel. TA5.5 was more closely correlated with [HA] than was β, and seems to be the preferred practical measure of buffering. The relationship between [HA] and TA5.5 differs between mono- and polybasic acids, so a separate analysis of products according to predominant acid type could improve multivariate models of erosive potential.

The effect of monoalkyl phosphates and fluoride on dissolution of hydroxyapatite, and interactions with saliva

The aims were to investigate the effect of monoalkyl phosphates (MAPs) and fluoride on dissolution rate of native and saliva-coated hydroxyapatite (HA). Fluoride at 300 mg/l (as NaF) inhibited dissolution of native HA by 12%, while potassium and sodium dodecyl phosphates (PDP, SDP), at 0.1% or higher, inhibited dissolution by 26-34%. MAPs, but not fluoride, also showed persistence of action. MAPs at 0.5% and fluoride at 300 mg/l were then tested separately against HA pre-treated with human saliva for 2 or 18 h. Agents were applied with brushing to half the specimens, and without brushing to the other half. In control (water-treated) specimens, pre-treatment of HA with human saliva reduced dissolution rate on average by 41% (2 h) and 63% (18 h). Brushing did not have a statistically significant effect on dissolution rate of saliva-coated specimens. In brushed specimens, fluoride significantly increased the inhibition due to 2- or 18-hour saliva pre-treatment. It is hypothesised that brushing partially removes the salivary film and allows KOH-soluble calcium fluoride formation at the surfaces of HA particles. Inhibition was reduced by PDP in 2-hour/non-brushed specimens and in 18-hour/brushed specimens. PDP did not affect dissolution rates in the remaining groups and SDP did not affect dissolution rate in any group. Possible reasons for these variable results are discussed. The experiments show that pre-treatment with saliva can significantly modify results of tests on potential anti-erosive agents and it is recommended that saliva pre-treatment should be a routine part of testing such agents.

Kinetic modeling of sodium in the lunar exosphere

Our knowledge about the lunar environment is based on a large volume of ground-based, remote, and in situ observations. These observations have been conducted at different times and sampled different pieces of such a complex system as the surface-bound exosphere of the Moon. Numerical modeling is the tool that can link results of these separate observations into a single picture. Being validated against previous measurements, models can be used for predictions and interpretation of future observations results. In this paper we present a kinetic model of the sodium exosphere of the Moon as well as results of its validation against a set of ground-based and remote observations. The unique characteristic of the model is that it takes the orbital motion of the Moon and the Earth into consideration and simulates both the exosphere as well as the sodium tail self-consistently. The extended computational domain covers the part of the Earth’s orbit at new Moon, which allows us to study the effect of Earth’s gravity on the lunar sodium tail. The model is fitted to a set of ground-based and remote observations by tuning sodium source rate as well as values of sticking, and accommodation coefficients. The best agreement of the model results with the observations is reached when all sodium atoms returning from the exosphere stick to the surface and the net sodium escape rate is about 5.3 × 1022 s−1.