Construction of three-phase data to model multiphase flow in porous media: Comparing an optimization approach to the finite element approach

Multiphase flow modelling is a major issue in the assessment of groundwater pollution. Three-phase flows are commonly governed by mathematical models that associate a pressure equation with two saturation equations. These equations involve a number of secondary variables that reflect the fluid behaviour in a porous medium. To improve the computational efficiency of multiphase flow simulators, several simplified reformulations of three-phase flow equations have been proposed. However, they require the construction of new secondary variables adapted to the reformulated flow equations. In this article, two different approaches are compared to quantify these variables. A numerical example is given for a typical fine sand.     

Numerical simulation of percussive drilling

This paper presents a novel dynamical model to analyze the long‐term response of a percussive drilling system. This departs from existing approaches that usually consider a single activation and bit/rock interaction cycle for the analysis of the process performance. The proposed model integrates the axial dynamics of an elastic piston and an elastic drill bit, a motion‐dependent pressure law to drive the piston, and a generalized bit/rock interaction law representative of the dynamic indentation taking place at the bit/rock interface. It applies to down‐the‐hole percussive drilling as well as top‐hole, with minor modifications. The model does not account for the angular motion or the hole cleaning, however. The model is first formulated mathematically; then, a finite‐dimensional approximation is proposed for computations. Numerical analyses of the model response, for a low‐size down‐the‐hole percussive system, follow. The period‐1 stationary response for the reference configuration is studied in detail, and parametric analyses assessing the influence on the rate of penetration of the bit/rock interaction parameters, the feed force, and the percussive activation parameters are conducted. These analyses reveal that the multiscale nature of the process is well captured by the model and recover expected trends for the influence of the parameters. They also suggest that a significant increase of the penetration rate can be achieved by increasing the percussive frequency. Copyright © 2015 John Wiley & Sons, Ltd.     

Sites occupancy of Fe3+ in Garfield Nontronite: A spectroscopic study

New data on the structure of Garfield nontronite have been produced by the use of different spectroscopic techniques: Mössbauer spectroscopy, optical spectroscopy, X-ray absorption edges and EXAFS and NMR. The tetrahedral iron content is found to be no higher than 1 percent. All iron atoms belong to the octahedral sheet excluding the possibility of the presence of non crystallized phases. Some ambiguities remain about the coherence of the octahedral sheet because of the presence of two doublets in the Mössbauer spectrum and at least two lineshapes in NMR spectra of OH which correspond to different environments.     

Synthesis of monticellite–forsterite and merwinite–forsterite symplectites in the CaO–MgO–SiO<Subscript>2</Subscript> model system: influence of temperature and water content on microstructure evolution

Homogeneous single crystals of synthetic monticellite with the composition \({\text{Ca}}_{0.88}{\text{Mg}}_{1.12}{\text{SiO}}_4\) (Mtc I) were annealed in a piston-cylinder apparatus at temperatures between 1000 and \(1200\,^{\circ }\hbox {C}\), pressures of 1.0–1.4 GPa, for run durations from 10 min to 24 h and applying bulk water contents ranging from 0.0 to 0.5 wt% of the total charge. At these conditions, Mtc I breaks down to a fine-grained, symplectic intergrowth. Thereby, two types of symplectites are produced: a first symplectite type (Sy I) is represented by an aggregate of rod-shaped forsterite immersed in a matrix of monticellite with end-member composition (Mtc II), and a second symplectite type (Sy II) takes the form of a lamellar merwinite–forsterite intergrowth. Both symplectites may form simultaneously, where the formation of Sy I is favoured by the presence of water. Sy I is metastable with respect to Sy II and is successively replaced by the latter. For both symplectite types, the characteristic spacing of the symplectite phases is independent of run duration and is only weeakly influenced by the water content, but it is strongly temperature dependent. It varies from about 400 nm at \(1000\,^{\circ }\hbox {C}\) to 1200 nm at \(1100\,^{\circ }\hbox {C}\) in Sy I, and from 300 nm at \(1000\,^{\circ }\hbox {C}\) to 700 nm at \(1200\,^{\circ }\hbox {C}\) in Sy II. A thermodynamic analysis reveals that the temperature dependence of the characteristic spacing of the symplectite phases is due to a relatively high activation energy for chemical segregation by diffusion within the reaction front as compared to the activation energy for interface reactions at the reaction front. The temperature dependence of the characteristic lamellar spacing and the temperature-time dependence of overall reaction progress have potential for applications in geo-thermometry and geo-speedometry.     

Hot Spots of Nitrification in the Elbe Estuary and Their Impact on Nitrate Regeneration

Estuaries act as an organic matter and nutrient filter in the transition between the land, rivers and the ocean. In the past, high nutrient and organic carbon load and low oxygen concentration made the Elbe River estuary (NW Europe) a sink for dissolved inorganic nitrogen. A recent reduction in loads and subsequent recovery of the estuary changed its biogeochemical function, so that nitrate is no longer removed on its transition towards the coastal North Sea. Nowadays in the estuary, nitrification appears to be a significant nitrate source. To quantify nitrification and determine actively nitrifying regions in the estuary, we measured the concentrations of ammonium, nitrite and nitrate, the dual stable isotopes of nitrate and net nitrification rates in the estuary on five cruises from August 2012 to August 2013. The nitrate concentration increased markedly downstream of the port of Hamburg in summer and spring, accompanied by a decrease of nitrate isotope values that was clearest in summer exactly at the location where nitrate concentration started to increase. Ammonium and nitrite peaked in the Hamburg port region (up to 18 and 8 μmol L^?1, respectively), and nitrification rates in this region were up to 7 μmol L^?1 day^?1. Our data show that coupled re-mineralization and nitrification are significant internal nitrate sources that almost double the estuary’s summer nitrate concentration. Furthermore, we find that the port of Hamburg is a hot spot of nitrification, whereas the maximum turbidity zone (MTZ) only plays a subordinate role in turnover of nitrate.     

Importance of climate,forest fires and human population size in the Holocene boreal forest composition change in northern Europe

The relative importance of climate, forest fires and human population size on long‐term boreal forest composition were statistically investigated at regional and local scales in Fennoscandia. We employ pollen data from lakes, reflecting regional vegetation, and small forest hollows, reflecting local vegetation, from Russia, Finland and Sweden to reconstruct the long‐term forest composition. As potential drivers of the Holocene forest dynamics we consider climate, generated from a climate model and oxygen isotope data, past forest fires generated from sedimentary charcoal data and human population size derived from radiocarbon dated archaeological findings. We apply the statistical method of variation partitioning to assess the relative importance of these environmental variables on long‐term boreal forest composition. The results show that climate is the main driver of the changes in Holocene boreal forest composition at the regional scale. However, at the local scale the role of climate is relatively small. In general, the importance of forest fires is low both at regional and local scales. The fact that both climate and forest fires explain relatively small proportions of variation in long‐term boreal vegetation in small forest hollow records demonstrates the complexity of factors affecting stand‐scale forest dynamics. The relative importance of human population size was low in both the prehistorical and the historical time periods. However, this is the first time that this type of data has been used to statistically assess the importance of human population size on boreal vegetation and the spatial representativeness of the data may cause bias to the analysis.     

Glassy orthopyroxene granodiorites of the Pannonian Basin: tracers of ultra-high-temperature deep-crustal anatexis triggered by Tertiary basaltic volcanism

Glassy orthopyroxene granodiorite-tonalite (named pincinite after type locality) was described from basaltic lapilli tuffs of the Pliocene maar near Pinciná village in the Slovakian part of the Pannonian Basin. Two pincinite types exhibit a qualitatively similar mineral composition (quartz, An_20–55 plagioclase, intergranular silicic glass with orthopyroxene and ilmenite, ±K-feldspar), but strongly different redox potential and formation PT conditions. Peraluminous pincinite is reduced (6–7% of total iron as Fe³⁺ in corundum-normative intergranular dacitic glass) and contains ilmenite with 8–10 mol% Fe₂O₃ and orthopyroxene dominated by ferrosilite. High-density (up to 0.85 g/cm³) primary CO₂ inclusions with minor H₂, CH₄, H₂S, CO and N₂ (<2 mol% total) are present in Qtz and Plg. Equilibrium PT conditions inferred from the intergranular Opx–Ilm–Glass assemblage and fluid density correspond to 1,170±50°C, 5.6±0.4 kbar, respectively. Metaluminous pincinite is more oxidised (25–27% of total iron as Fe³⁺ in diopside-normative intergranular glass of rhyolite–trachyte–dacite composition) and contains Fe₂O₃-rich ilmenite (17–29 mol%) associated with enstatite. Fluid inclusions are composed of CO₂–H₂O mixtures with up to 38 mol% H₂O. Raman spectroscopy revealed H₂S along with dominant CO₂ in the carbonic phase. Equilibrium PT parameters for the intergranular Opx–Ilm–Glass assemblage correspond to 740±15°C, 2.8±0.1 kbar, respectively. Reducing gas species (<2 mol% total) in the CO₂-inclusions of the peraluminous pincinite resulted from hydrogen diffusion due to fH₂ gradient imposed during decrease of redox potential from the log fO₂ values near QFM during Qtz + Plg growth, to QFM-2 incidental to the superimposed Opx + Ilm assemblage in the intergranular melt. The decrease in oxygen fugacity was recorded also in the metaluminous pincinite, where log fO₂ values changed from ~QFM + 2.6 to QFM + 0.4, but hydrogen diffusion did not occur. Absence of OH-bearing minerals, major and trace element abundances (e.g. REE 300–320, Nb 55–57, Th 4–31, Zr 240–300 ppm, FeO_tot/MgO up to 11), and Sr–O isotope ratios in the pincinites are diagnostic of high-temperature anorogenic magmas originated by dehydration melting of biotite in quartz-feldspathoid crust (⁸⁷Sr/⁸⁶Sr>0.705–0.706, ¹⁸O>9 V-SMOW) around alkali basalt reservoir in depths between 17 and 20 km, and around late stage derivatives of the basalt fractionation, intruding the crust up to depths of 10–11 km. Low water activity in the pincinite parental melt was caused by CO₂-flux from the Tertiary basaltic reservoirs and intrusions. The anatexis leads to generation of a melt-depleted granulitic crust beneath the Pannonian Basin, and the pincinites are interpreted as equivalents of igneous charnockites and enderbites quenched at temperatures above solidus and unaffected by sub-solidus re-equilibration and metamorphic overprint.