Calculating the effective permeability of sandstone with multiscale lattice Boltzmann/finite element simulations

The lattice Boltzmann (LB) method is an efficient technique for simulating fluid flow through individual pores of complex porous media. The ease with which the LB method handles complex boundary conditions, combined with the algorithm’s inherent parallelism, makes it an elegant approach to solving flow problems at the sub-continuum scale. However, the realities of current computational resources can limit the size and resolution of these simulations. A major research focus is developing methodologies for upscaling microscale techniques for use in macroscale problems of engineering interest. In this paper, we propose a hybrid, multiscale framework for simulating diffusion through porous media. We use the finite element (FE) method to solve the continuum boundary-value problem at the macroscale. Each finite element is treated as a sub-cell and assigned permeabilities calculated from subcontinuum simulations using the LB method. This framework allows us to efficiently find a macroscale solution while still maintaining information about microscale heterogeneities. As input to these simulations, we use synchrotron-computed 3D microtomographic images of a sandstone, with sample resolution of 3.34 μm. We discuss the predictive ability of these simulations, as well as implementation issues. We also quantify the lower limit of the continuum (Darcy) scale, as well as identify the optimal representative elementary volume for the hybrid LB–FE simulations.     

40Ar/39Ar and K–Ar geochronology of magmatic and hydrothermal events in a classic low-suphidation epithermal bonanza deposit: El Peñon, northern Chile

The epithermal El Peñon gold–silver deposit consists of quartz–adularia veins emplaced within a late Upper Paleocene rhyolitic dome complex, located in the Paleocene–Lower Eocene Au–Ag belt of northern Chile. Detailed K–Ar and ⁴⁰Ar/³⁹Ar geochronology on volcano–plutonic rocks and hydrothermal minerals were carried out to constrain magmatic and hydrothermal events. The Paleocene to Lower Eocene magmatism in the El Peñon area is confined to a rhomb-shaped basin, which was controlled by N–S trending normal faults and both NE- and NW-trending transtensional fault systems. The earliest products of the basin-filling sequences comprise of Middle to Upper Paleocene (～59–55 Ma) welded rhyolitic ignimbrites and andesitic to dacitic lavas, with occasional dacitic dome complexes. Later, rhyolitic and dacitic dome complexes (～55–52 Ma) represent the waning stages of volcanism during the latest Upper Paleocene and the earliest Eocene. Lower Eocene porphyry intrusives (～48–43 Ma) mark the end of the magmatism in the basin and a change to a compressive tectonomagmatic regime. ⁴⁰Ar/³⁹Ar geochronology of hydrothermal adularia from the El Peñon deposit yields ages between 51.0±0.6 and 53.1±0.5 Ma. These results suggest that mineralization occurred slightly after the emplacement of the El Peñon rhyolitic dome at 54.5±0.6 Ma (⁴⁰Ar/³⁹Ar age) and was closely tied to later dacitic–rhyodacitic bodies of 52 to 53 Ma (K–Ar ages), probably as short-lived pulses related to single volcanic events.     

Globally increased pelagic carbonate production during the Mid-Brunhes dissolution interval and the CO2 paradox of MIS 11

The Mid-Brunhes dissolution interval (MBDI) represents a period of global carbonate dissolution, lasting several hundred thousand years, centred around Marine Isotope Stage (MIS) 11. Here we report the effects of dissolution in ODP core 982, taken from 1134 m in the North Atlantic. Paradoxically, records of atmospheric CO₂ from Antarctic ice-cores reveal no long term trend over the last 400 kyr and suggest that CO₂ during MIS 11 was no higher than during the present interglacial. We suggest that a global increase in pelagic carbonate production during this period, possibly related to the proliferation of the Gephyrocapsa coccolithophore, could have altered marine carbonate chemistry in such a way as to drive increased dissolution under the constraints of steady state. An increase in the production of carbonate in surface waters would cause a drawdown of global carbonate saturation and increase dissolution at the seafloor. In order to reconcile the record of atmospheric CO₂ variability we suggest that an increase in the flux of organic matter from the surface to deep ocean, associated with either a net increase in primary production or the enhanced ballasting effect provided by an increased flux of CaCO₃, could have countered the effect of increased calcification on CO₂.     

Structural evolution of zeolite levyne under hydrostatic and non-hydrostatic pressure: geometric modelling

This is an exploratory study on the high-pressure (HP) structural evolution of a zeolitic framework (with LEV topology) on the basis of geometric modelling and previously published accurate unit-cell constants measured by means of single-crystal X-ray diffraction. The geometric simulations for 11 P values from 0 to 5 GPa gives more insight into the HP-behaviour of levyne, showing that the anomalous elastic behaviour of this zeolite observed under hydrostatic conditions at low P (P<1 GPa) is due to a double change in the compressional mechanism. Since the geometric simulation is not restricted to using the experimentally determined cell parameters, simulations of uniaxial compression along the [001] direction and of compression in the (001) plane have been performed, shedding more light on the compression mechanisms under non-hydrostatic regimes, which are difficult to access experimentally. The mechanisms associated with compressions along different axes provide insight into the hydrostatic compression mechanisms leading to the anomalous elastic behaviour.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Design, analysis, and inference for studies comparing thematic accuracy of classified remotely sensed data: a special case of map comparison

Assessing thematic map accuracy is a special type of map comparison that is frequently applied to remote sensing classification problems. For map comparisons in the accuracy assessment setting, one map represents the classified output and the other map represents the true or “reference” condition. Several articles in this special issue describe state-of-the-art map comparison analysis tools that could serve to quantify accuracy of a single map. However, accuracy assessment objectives generally extend beyond describing accuracy of a single map to comparing accuracy of several maps. Consequently, interest focuses on comparing map comparison measures when these measures are used to represent accuracy. The virtual workshop emphasizes the analysis component of map comparisons, but it is also important to examine the underlying study designs generating the data input into these analyses. The study designs for accuracy comparisons implemented in remote sensing practice often investigate only a single test site, thus limiting our ability to generalize the results of these accuracy comparisons. Map accuracy comparison studies can be designed to provide stronger generalizations by incorporating experimental design principles such as replication and blocking, and identifying an experimental unit appropriate for the application. It is also important to recognize the role of statistical hypothesis testing and inference for different objectives that motivate map accuracy comparisons. Deciding which of two maps to use for a particular site can be addressed by enumerative inference and does not require hypothesis testing. For the objective of a more general comparison of classification procedures, analytic inference is appropriate and hypothesis testing plays a more prominent role.     

Contents of Volume 26

Table of Contents

Contents of Volume 26     

A Wind Tunnel Model for Quantifying Fluxes in the Urban Boundary Layer

Transport of pollution and heatout of streets into the boundary layer above is not currently understood and so fluxes cannot be quantified. Scalar concentration within the street is determined by the flux out of it and so quantifying fluxes for turbulent flow over a rough urban surface is essential. We have developed a naphthalene sublimation technique to measure transfer from a two-dimensional street canyon in a wind tunnel for the case of flow perpendicular to the street. The street was coated with naphthalene, which sublimes at room temperature, so that the vapour represented the scalar source. The transfer velocity w_T relates the flux out of the canyon to the concentration within it and is shown to be linearly related to windspeed above the street. The dimensionless transfer coefficient w_T/U represents the ventilation efficiency of the canyon (here, w_T is a transfer velocity,U is the wind speed at the boundary-layer top). Observed values are between 1.5 and 2.7 ×10^-3 and, for the case where H/W0 (ratio of buildingheight to street width), values are in the same range as estimates of transfer from a flat plate, giving confidence that the technique yields accurate values for street canyon scalar transfer. w_T/U varies with aspect ratio (H/W), reaching a maximum in the wake interference regime (0.3 < H/W < 0.65). However, when upstream roughness is increased, the maximum in w_T/U reduces, suggesting that street ventilation is less sensitive to H/W when the flow is in equilibrium with the urban surface. The results suggest that using naphthalene sublimation with wind-tunnel models of urban surfaces can provide a direct measure of area-averaged scalar fluxes.     

Preface

Preface

Preface     

A comparison of empirical orthogonal decomposition methods in baroclinic flows

The relative merits of three contrasting empirical orthogonal decomposition methods in common use (namely, Proper Orthogonal Decomposition, Biorthogonal Decomposition and Multivariate Singular Systems Analysis) are considered as applied to baroclinic flow data. The regimes analysed are a steady, drifting wave, a modulated amplitude vacillating wave flow and a neighbouring multi-mode state which exhibits intermittency. The results are used to make a qualitative comparison of the methods in terms of convergence properties, variance capture and eigenfunction structure. The feasibility of using the resulting empirical orthogonal functions to transform partial differential equations to ordinary differential equations by Galerkin projection is mentioned.