A robust mesh optimisation method for multiphase porous media flows

Flows of multiple fluid phases are common in many subsurface reservoirs. Numerical simulation of these flows can be challenging and computationally expensive. Dynamic adaptive mesh optimisation and related approaches, such as adaptive grid refinement can increase solution accuracy at reduced computational cost. However, in models or parts of the model domain, where the local Courant number is large, the solution may propagate beyond the region in which the mesh is refined, resulting in reduced solution accuracy, which can never be recovered. A methodology is presented here to modify the mesh within the non-linear solver. The method allows efficient application of dynamic mesh adaptivity techniques even with high Courant numbers. These high Courant numbers may not be desired but a consequence of the heterogeneity of the domain. Therefore, the method presented can be considered as a more robust and accurate version of the standard dynamic mesh adaptivity techniques.     

Fine-scale volume heterogeneity measurements in sand

As part of the effort to characterize the acoustic environment during the high frequency sediment acoustics experiment (SAX99), fine-scale variability of sediment density was measured by an in situ technique and by core analysis. The in situ measurement was accomplished by a newly developed instrument that measures sediment conductivity. The conductivity measurements were conducted on a three-dimensional (3-D) grid, hence providing a set of data suited for assessing sediment spatial variability. A 3-D sediment porosity matrix is obtained from the conductivity data through an empirical relationship (Archie's Law). From the porosity matrix, sediment bulk density is estimated from known average grain density. A number of cores were taken at the SAX99 site, and density variations were measured using laboratory techniques. The power spectra were estimated from both techniques and were found to be appropriately fit by a power-law. The exponents of the horizontal one-dimensional (1-D) power-law spectra have a depth-dependence and range from 1.72 to 2.41. The vertical 1-D spectra have the same form, but with an exponent of 2.2. It was found that most of the density variability is within the top 5 mm of the sediment, which suggests that sediment volume variability will not have major impact on acoustic scattering when the sound frequency is below 100 kHz. At higher frequencies, however, sediment volume variability is likely to play an important role in sound scattering.     

Composition,age and origin of two proterozoic diorite—tonalite complexes in the Arabian shield

Large amounts of diorite—tonalite magma were intruded into the island-arc successions of the southern Arabian shield between ca. 900 and 700 Ma ago. Major oxide, trace element, rare earth (REE) and isotopic data are presented for two plutons exemplifying older and younger members of this plutonic phase. The Thurrat pluton, which was emplaced into virtually unmetamorphosed volcanics of sequence B, has yielded a 10-point Rb-Sr isochron indicating emplacement 744 ± 22 Ma ago and an initial ${}^{87}\text{Sr}{}^{86}\text{Sr}$ ratio of 0.70281. It consists of gabbros, diorites, quartz diorites, tonalites and low-Al₂O₃ trondhjemites which are compositionally primitive, with depleted lithophile element contents and flat REE patterns with a negative Eu anomaly in the trondhjemites. The magma was most probably mantle-derived, and analogies with other trondhjemitic plutons suggest that it was probably emplaced in an island-arc setting. The Bidah pluton, which was emplaced into a compositionally very immature succession of metamorphosed volcanics and volcaniclastics of sequence C, has yielded a nine-point near-isochron (MSWD = 2.86) indicating an age of 901 ± 37Ma and an initial ${}^{87}\text{Sr}{}^{86}\text{Sr}$ ratio of 0.70246. This date is accepted as the age of emplacement. The rocks are compositionally primitive gabbros, diorites, quartz diorites, tonalites and granodiorites with depleted lithophile element contents, and flat REE patterns with a negative Eu anomaly in the more siliceous components. The country-rock volcanics were formed in an immature island-arc environment, and the composition of the Bidah pluton is suggestive of a mainly mantle-derived magma emplaced into that arc. The data for these two plutons therefore provide further evidence that most of the rock material added to the Arabian shield between 900 and 700 Ma ago was derived from the mantle.     

Geochemical dispersion of heavy metals via organic complexing: A laboratory study of copper,lead, zinc,and nickel behaviour at a simulated sediment-water boundary

Data are presented in this study from laboratory model experiments describing the behaviour of Cu, Pb, Zn and Ni at a simulated sediment-water boundary. The interactions involved are sorption by kaolin and by bentonite, organic complexing in solution by fulvic acid and by humic acid, carbonate reactions, hydrolysis, and desorption of the cations from a clay-bound phase and from their metal hydroxides by the organic acids. The organic acids increase the solubility of Cu, Zn and Ni in the presence of clay. The Pb solubility is variable and can even decrease, particularly at acidic pH, with organic complexing likely due to colloidal coagulation. Both Zn and Ni are influenced by hydrolysis at basic pH. When carbonate was added to the metal-organic acid-clay mixtures, a further decrease in solubility was observed for Ni and, to a lesser extent, for Zn. The organic acids prove capable of remobilizing Cu, Pb, Zn and Ni from the solid phases examined. However, there is a general kinetic hindrance to the desorption particularly at basic pH. Copper desorption appears to be the most kinetically hindered. Conclusions pertinent to the geochemical dispersion of these metals are drawn.     

A resistivity section of Holy Island (Anglesey)

An account is given of the use of square array resistivity techniques to derive an apparent resistivity section down the length of Holy Island (Anglesey). Possible interpretations of this section in terms of the distribution of true resistivity with depth are given and discussed briefly in relation to proposed geological structure.     

From the Sun to the Earth: The 13 May 2005 Coronal Mass Ejection

We report the results of a multi-instrument, multi-technique, coordinated study of the solar eruptive event of 13 May 2005. We discuss the resultant Earth-directed (halo) coronal mass ejection (CME), and the effects on the terrestrial space environment and upper Earth atmosphere. The interplanetary CME (ICME) impacted the Earth’s magnetosphere and caused the most-intense geomagnetic storm of 2005 with a Disturbed Storm Time (Dst) index reaching ?263 nT at its peak. The terrestrial environment responded to the storm on a global scale. We have combined observations and measurements from coronal and interplanetary remote-sensing instruments, interplanetary and near-Earth in-situ measurements, remote-sensing observations and in-situ measurements of the terrestrial magnetosphere and ionosphere, along with coronal and heliospheric modelling. These analyses are used to trace the origin, development, propagation, terrestrial impact, and subsequent consequences of this event to obtain the most comprehensive view of a geo-effective solar eruption to date. This particular event is also part of a NASA-sponsored Living With a Star (LWS) study and an on-going US NSF-sponsored Solar, Heliospheric, and INterplanetary Environment (SHINE) community investigation.     

Trace Element Analysis of NIST SRM 614 and 616 Glass Reference Materials by Laser Ablation Microprobe-Inductively Coupled Plasma-Mass Spectrometry

Fifty elements in NIST SRM 614 and 616 glass reference materials were determined by laser ablation microprobe-inductively coupled plasma-mass spectrometry (LAM-ICP-MS). The values determined for NIST SRM 614 agreed well with the NIST-certified and information values (mean relative difference ± 3.6%), except for B, Sc and Sb. The values determined for NIST SRM 616 agreed with the NIST-certified and information values within a mean relative difference of ± 1.5%, except for B, Sc and Ga. In addition, at an 80 μm sampling scale, NIST SRM 614 and 616 glass discs were homogeneous for trace elements within the observed precisions of 5 and 15% (mean), respectively. Detection limits were in the range 0.01 - 0.3 μg g⁻¹ for elements of lower mass numbers (amu < 80) and 1 - 10 ng g⁻¹ for heavy elements (amu > 80). Detection at the sub ng g⁻¹ level is possible for most of the heavy elements by using an ablation pit size larger than 10 0 μm.     

Thermodynamically constrained averaging theory approach for modeling flow and transport phenomena in porous medium systems: 6. Two-fluid-phase flow

This work is the sixth in a series of papers on the thermodynamically constrained averaging theory (TCAT) approach for modeling flow and transport phenomena in multiscale porous medium systems. Building upon the general TCAT framework and the mathematical foundation presented in previous works, the limiting case of connected two-fluid-phase flow is considered. A constrained entropy inequality is developed based upon a set of primary restrictions. Formal approximations are introduced to deduce a general simplified entropy inequality (SEI). The SEI is used along with secondary restrictions and closure approximations consistent with the SEI to produce a general functional form of a two-phase-flow model. The general model is in turn simplified to yield a hierarchy of models by neglecting common curves and by neglecting both common curves and interfaces. The simplest case considered corresponds to a traditional two-phase-flow model. The more sophisticated models including interfaces and common curves are more physically realistic than traditional models. All models in the hierarchy are posed in terms of precisely defined variables that allow for a rigorous connection with the microscale. The explicit nature of the restrictions and approximations used in developing this hierarchy of models provides a clear means to both understand the limitations of traditional models and to build upon this work to produce more realistic models.