Bio-inspired Murray materials for mass transfer and activity

NASA Astrophysics Data System (ADS)

Zheng, Xianfeng; Shen, Guofang; Wang, Chao; Li, Yu; Dunphy, Darren; Hasan, Tawfique; Brinker, C. Jeffrey; Su, Bao-Lian

2017-04-01

Both plants and animals possess analogous tissues containing hierarchical networks of pores, with pore size ratios that have evolved to maximize mass transport and rates of reactions. The underlying physical principles of this optimized hierarchical design are embodied in Murray's law. However, we are yet to realize the benefit of mimicking nature's Murray networks in synthetic materials due to the challenges in fabricating vascularized structures. Here we emulate optimum natural systems following Murray's law using a bottom-up approach. Such bio-inspired materials, whose pore sizes decrease across multiple scales and finally terminate in size-invariant units like plant stems, leaf veins and vascular and respiratory systems provide hierarchical branching and precise diameter ratios for connecting multi-scale pores from macro to micro levels. Our Murray material mimics enable highly enhanced mass exchange and transfer in liquid-solid, gas-solid and electrochemical reactions and exhibit enhanced performance in photocatalysis, gas sensing and as Li-ion battery electrodes.

Frequency Dependent Macro-dispersion Induced by Oscillatory Inputs and Spatial Heterogeneity

NASA Astrophysics Data System (ADS)

Rajabi, F.; Battiato, I.

2017-12-01

Elucidating flow and transport processes at the pore scale is the cornerstone of most hydrologic studies in the subsurface. This becomes even more imperative when the system is subject to a cyclic forcing. Such temporal variations with evolving heterogeneity of time scales spanning from days to years can influence transport phenomena at the pore level, e.g. yearly freeze/thaw in the thin active layer of soil above permafrost zone whose thickness increases throughout the thaw season. Moreover, understanding the interactions of different physical phenomena at the pore scale is key to predict the behavior at the continuum scale. Yet, the connection between periodic inputs at the pore scale and macrotransport is to a great extent unknown. In the spirit of homogenization technique, we derived a macrotime continuum-scale equation as well as expressions for the effective transport coefficients. The macrodispersion arises from contributions of molecular diffusion, spatial heterogeneity and time-dependent fluctuations. Moreover, we have quantified the solute spreading by effective dispersion in terms of dimensionless numbers (Pe, Da, and Strouhal), i.e. expressing the interplay of molecular diffusion, advection, reaction and signal frequency. Yet, as every macroscopic model, spatiotemporally averaged models can breakdown when certain criteria are violated. This makes the continuum scale equation a poor approximation for the processes at the pore scale. To this end, we also provide the conditions under which the space-time averaged equations accurately describe pore-scale processes. In addition, this study gives a robust evidence that transverse mixing can in fact benefit from fluctuating boundary forcing due to the interaction of temporal fluctuations and molecular diffusion. Furthermore, it provides a robust quantitative foundation for designing the desired systems since the interplay of geometry and external forcing has been directly connected to each other in terms of dimensionless (St) number. We compare our theoretical framework with data from an experiment performed on several micro-channels with different geometry and different frequencies of injection at the inlet. The proposed formulation is found to provide remarkably good predictions and correctly explain the experimental mixing dynamics.

Application of Mortar Coupling in Multiscale Modelling of Coupled Flow, Transport, and Biofilm Growth in Porous Media

NASA Astrophysics Data System (ADS)

Laleian, A.; Valocchi, A. J.; Werth, C. J.

2017-12-01

Multiscale models of reactive transport in porous media are capable of capturing complex pore-scale processes while leveraging the efficiency of continuum-scale models. In particular, porosity changes caused by biofilm development yield complex feedbacks between transport and reaction that are difficult to quantify at the continuum scale. Pore-scale models, needed to accurately resolve these dynamics, are often impractical for applications due to their computational cost. To address this challenge, we are developing a multiscale model of biofilm growth in which non-overlapping regions at pore and continuum spatial scales are coupled with a mortar method providing continuity at interfaces. We explore two decompositions of coupled pore-scale and continuum-scale regions to study biofilm growth in a transverse mixing zone. In the first decomposition, all reaction is confined to a pore-scale region extending the transverse mixing zone length. Only solute transport occurs in the surrounding continuum-scale regions. Relative to a fully pore-scale result, we find the multiscale model with this decomposition has a reduced run time and consistent result in terms of biofilm growth and solute utilization. In the second decomposition, reaction occurs in both an up-gradient pore-scale region and a down-gradient continuum-scale region. To quantify clogging, the continuum-scale model implements empirical relations between porosity and continuum-scale parameters, such as permeability and the transverse dispersion coefficient. Solutes are sufficiently mixed at the end of the pore-scale region, such that the initial reaction rate is accurately computed using averaged concentrations in the continuum-scale region. Relative to a fully pore-scale result, we find accuracy of biomass growth in the multiscale model with this decomposition improves as the interface between pore-scale and continuum-scale regions moves downgradient where transverse mixing is more fully developed. Also, this decomposition poses additional challenges with respect to mortar coupling. We explore these challenges and potential solutions. While recent work has demonstrated growing interest in multiscale models, further development is needed for their application to field-scale subsurface contaminant transport and remediation.

Changes in Pore Water Quality After Peatland Restoration: Assessment of a Large-Scale, Replicated Before-After-Control-Impact Study in Finland

NASA Astrophysics Data System (ADS)

Menberu, Meseret Walle; Marttila, Hannu; Tahvanainen, Teemu; Kotiaho, Janne S.; Hokkanen, Reijo; Kløve, Bjørn; Ronkanen, Anna-Kaisa

2017-10-01

Drainage is known to affect peatland natural hydrology and water quality, but peatland restoration is considered to ameliorate peatland degradation. Using a replicated BACIPS (Before-After-Control-Impact Paired Series) design, we investigated 24 peatlands, all drained for forestry and subsequently restored, and 19 pristine control boreal peatlands with high temporal and spatial resolution data on hydroclimate and pore water quality. In drained conditions, total nitrogen (Ntot), total phosphorus (Ptot), and dissolved organic carbon (DOC) in pore water were several-fold higher than observed at pristine control sites, highlighting the impacts of long-term drainage on pore water quality. In general, pore water DOC and Ntot decreased after restoration measures but still remained significantly higher than at pristine control sites, indicating long time lags in restoration effects. Different peatland classes and trophic levels (vegetation gradient) responded differently to restoration, primarily due to altered hydrology and varying acidity levels. Sites that were hydrologically overrestored (inundated) showed higher Ptot, Ntot, and DOC than well-restored or insufficiently restored sites, indicating the need to optimize natural-like hydrological regimes when restoring peatlands drained for forestry. Rich fens (median pH 6.2-6.6) showed lower pore water Ptot, Ntot, and DOC than intermediate and poor peats (pH 4.0-4.6) both before and after restoration. Nutrients and DOC in pore water increased in the first year postrestoration but decreased thereafter. The most important variables related to pore water quality were trophic level, peatland class, water table level, and soil and air temperature.

Micro-CT Pore Scale Study Of Flow In Porous Media: Effect Of Voxel Resolution

NASA Astrophysics Data System (ADS)

Shah, S.; Gray, F.; Crawshaw, J.; Boek, E.

2014-12-01

In the last few years, pore scale studies have become the key to understanding the complex fluid flow processes in the fields of groundwater remediation, hydrocarbon recovery and environmental issues related to carbon storage and capture. A pore scale study is often comprised of two key procedures: 3D pore scale imaging and numerical modelling techniques. The essence of a pore scale study is to test the physics implemented in a model of complicated fluid flow processes at one scale (microscopic) and then apply the model to solve the problems associated with water resources and oil recovery at other scales (macroscopic and field). However, the process of up-scaling from the pore scale to the macroscopic scale has encountered many challenges due to both pore scale imaging and modelling techniques. Due to the technical limitations in the imaging method, there is always a compromise between the spatial (voxel) resolution and the physical volume of the sample (field of view, FOV) to be scanned by the imaging methods, specifically X-ray micro-CT (XMT) in our case In this study, a careful analysis was done to understand the effect of voxel size, using XMT to image the 3D pore space of a variety of porous media from sandstones to carbonates scanned at different voxel resolution (4.5 μm, 6.2 μm, 8.3 μm and 10.2 μm) but keeping the scanned FOV constant for all the samples. We systematically segment the micro-CT images into three phases, the macro-pore phase, an intermediate phase (unresolved micro-pores + grains) and the grain phase and then study the effect of voxel size on the structure of the macro-pore and the intermediate phases and the fluid flow properties using lattice-Boltzmann (LB) and pore network (PN) modelling methods. We have also applied a numerical coarsening algorithm (up-scale method) to reduce the computational power and time required to accurately predict the flow properties using the LB and PN method.

Micro-computed tomography pore-scale study of flow in porous media: Effect of voxel resolution

NASA Astrophysics Data System (ADS)

Shah, S. M.; Gray, F.; Crawshaw, J. P.; Boek, E. S.

2016-09-01

A fundamental understanding of flow in porous media at the pore-scale is necessary to be able to upscale average displacement processes from core to reservoir scale. The study of fluid flow in porous media at the pore-scale consists of two key procedures: Imaging - reconstruction of three-dimensional (3D) pore space images; and modelling such as with single and two-phase flow simulations with Lattice-Boltzmann (LB) or Pore-Network (PN) Modelling. Here we analyse pore-scale results to predict petrophysical properties such as porosity, single-phase permeability and multi-phase properties at different length scales. The fundamental issue is to understand the image resolution dependency of transport properties, in order to up-scale the flow physics from pore to core scale. In this work, we use a high resolution micro-computed tomography (micro-CT) scanner to image and reconstruct three dimensional pore-scale images of five sandstones (Bentheimer, Berea, Clashach, Doddington and Stainton) and five complex carbonates (Ketton, Estaillades, Middle Eastern sample 3, Middle Eastern sample 5 and Indiana Limestone 1) at four different voxel resolutions (4.4 μm, 6.2 μm, 8.3 μm and 10.2 μm), scanning the same physical field of view. Implementing three phase segmentation (macro-pore phase, intermediate phase and grain phase) on pore-scale images helps to understand the importance of connected macro-porosity in the fluid flow for the samples studied. We then compute the petrophysical properties for all the samples using PN and LB simulations in order to study the influence of voxel resolution on petrophysical properties. We then introduce a numerical coarsening scheme which is used to coarsen a high voxel resolution image (4.4 μm) to lower resolutions (6.2 μm, 8.3 μm and 10.2 μm) and study the impact of coarsening data on macroscopic and multi-phase properties. Numerical coarsening of high resolution data is found to be superior to using a lower resolution scan because it avoids the problem of partial volume effects and reduces the scaling effect by preserving the pore-space properties influencing the transport properties. This is evidently compared in this study by predicting several pore network properties such as number of pores and throats, average pore and throat radius and coordination number for both scan based analysis and numerical coarsened data.

Role of medium heterogeneity and viscosity contrast in miscible flow regimes and mixing zone growth: A computational pore-scale approach

NASA Astrophysics Data System (ADS)

Afshari, Saied; Hejazi, S. Hossein; Kantzas, Apostolos

2018-05-01

Miscible displacement of fluids in porous media is often characterized by the scaling of the mixing zone length with displacement time. Depending on the viscosity contrast of fluids, the scaling law varies between the square root relationship, a sign for dispersive transport regime during stable displacement, and the linear relationship, which represents the viscous fingering regime during an unstable displacement. The presence of heterogeneities in a porous medium significantly affects the scaling behavior of the mixing length as it interacts with the viscosity contrast to control the mixing of fluids in the pore space. In this study, the dynamics of the flow and transport during both unit and adverse viscosity ratio miscible displacements are investigated in heterogeneous packings of circular grains using pore-scale numerical simulations. The pore-scale heterogeneity level is characterized by the variations of the grain diameter and velocity field. The growth of mixing length is employed to identify the nature of the miscible transport regime at different viscosity ratios and heterogeneity levels. It is shown that as the viscosity ratio increases to higher adverse values, the scaling law of mixing length gradually shifts from dispersive to fingering nature up to a certain viscosity ratio and remains almost the same afterwards. In heterogeneous media, the mixing length scaling law is observed to be generally governed by the variations of the velocity field rather than the grain size. Furthermore, the normalization of mixing length temporal plots with respect to the governing parameters of viscosity ratio, heterogeneity, medium length, and medium aspect ratio is performed. The results indicate that mixing length scales exponentially with log-viscosity ratio and grain size standard deviation while the impact of aspect ratio is insignificant. For stable flows, mixing length scales with the square root of medium length, whereas it changes linearly with length during unstable flows. This scaling procedure allows us to describe the temporal variation of mixing length using a generalized curve for various combinations of the flow conditions and porous medium properties.

Global sensitivity analysis of multiscale properties of porous materials

NASA Astrophysics Data System (ADS)

Um, Kimoon; Zhang, Xuan; Katsoulakis, Markos; Plechac, Petr; Tartakovsky, Daniel M.

2018-02-01

Ubiquitous uncertainty about pore geometry inevitably undermines the veracity of pore- and multi-scale simulations of transport phenomena in porous media. It raises two fundamental issues: sensitivity of effective material properties to pore-scale parameters and statistical parameterization of Darcy-scale models that accounts for pore-scale uncertainty. Homogenization-based maps of pore-scale parameters onto their Darcy-scale counterparts facilitate both sensitivity analysis (SA) and uncertainty quantification. We treat uncertain geometric characteristics of a hierarchical porous medium as random variables to conduct global SA and to derive probabilistic descriptors of effective diffusion coefficients and effective sorption rate. Our analysis is formulated in terms of solute transport diffusing through a fluid-filled pore space, while sorbing to the solid matrix. Yet it is sufficiently general to be applied to other multiscale porous media phenomena that are amenable to homogenization.

Coupling lattice Boltzmann and continuum equations for flow and reactive transport in porous media.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Coon, Ethan; Porter, Mark L.; Kang, Qinjun

2012-06-18

In spatially and temporally localized instances, capturing sub-reservoir scale information is necessary. Capturing sub-reservoir scale information everywhere is neither necessary, nor computationally possible. The lattice Boltzmann Method for solving pore-scale systems. At the pore-scale, LBM provides an extremely scalable, efficient way of solving Navier-Stokes equations on complex geometries. Coupling pore-scale and continuum scale systems via domain decomposition. By leveraging the interpolations implied by pore-scale and continuum scale discretizations, overlapping Schwartz domain decomposition is used to ensure continuity of pressure and flux. This approach is demonstrated on a fractured medium, in which Navier-Stokes equations are solved within the fracture while Darcy'smore » equation is solved away from the fracture Coupling reactive transport to pore-scale flow simulators allows hybrid approaches to be extended to solve multi-scale reactive transport.« less

The Backscattering Phase Function for a Sphere with a Two-Scale Relief of Rough Surface

NASA Astrophysics Data System (ADS)

Klass, E. V.

2017-12-01

The backscattering of light from spherical surfaces characterized by one and two-scale roughness reliefs has been investigated. The analysis is performed using the three-dimensional Monte-Carlo program POKS-RG (geometrical-optics approximation), which makes it possible to take into account the roughness of objects under study by introducing local geometries of different levels. The geometric module of the program is aimed at describing objects by equations of second-order surfaces. One-scale roughness is set as an ensemble of geometric figures (convex or concave halves of ellipsoids or cones). The two-scale roughness is modeled by convex halves of ellipsoids, with surface containing ellipsoidal pores. It is shown that a spherical surface with one-scale convex inhomogeneities has a flatter backscattering phase function than a surface with concave inhomogeneities (pores). For a sphere with two-scale roughness, the dependence of the backscattering intensity is found to be determined mostly by the lower-level inhomogeneities. The influence of roughness on the dependence of the backscattering from different spatial regions of spherical surface is analyzed.

Bio-inspired Murray materials for mass transfer and activity

PubMed Central

Zheng, Xianfeng; Shen, Guofang; Wang, Chao; Li, Yu; Dunphy, Darren; Hasan, Tawfique; Brinker, C. Jeffrey; Su, Bao-Lian

2017-01-01

Both plants and animals possess analogous tissues containing hierarchical networks of pores, with pore size ratios that have evolved to maximize mass transport and rates of reactions. The underlying physical principles of this optimized hierarchical design are embodied in Murray's law. However, we are yet to realize the benefit of mimicking nature's Murray networks in synthetic materials due to the challenges in fabricating vascularized structures. Here we emulate optimum natural systems following Murray's law using a bottom-up approach. Such bio-inspired materials, whose pore sizes decrease across multiple scales and finally terminate in size-invariant units like plant stems, leaf veins and vascular and respiratory systems provide hierarchical branching and precise diameter ratios for connecting multi-scale pores from macro to micro levels. Our Murray material mimics enable highly enhanced mass exchange and transfer in liquid–solid, gas–solid and electrochemical reactions and exhibit enhanced performance in photocatalysis, gas sensing and as Li-ion battery electrodes. PMID:28382972

Pore-level mechanics of foam generation and coalescence in the presence of oil.

PubMed

Almajid, Muhammad M; Kovscek, Anthony R

2016-07-01

The stability of foam in porous media is extremely important for realizing the advantages of foamed gas on gas mobility reduction. Foam texture (i.e., bubbles per volume of gas) achieved is dictated by foam generation and coalescence processes occurring at the pore-level. For foam injection to be widely applied during gas injection projects, we need to understand these pore-scale events that lead to foam stability/instability so that they are modeled accurately. Foam flow has been studied for decades, but most efforts focused on studying foam generation and coalescence in the absence of oil. Here, the extensive existing literature is reviewed and analyzed to identify open questions. Then, we use etched-silicon micromodels to observe foam generation and coalescence processes at the pore-level. Special emphasis is placed on foam coalescence in the presence of oil. For the first time, lamella pinch-off as described by Myers and Radke [40] is observed in porous media and documented. Additionally, a new mechanism coined "hindered generation" is found. Hindered generation refers to the role oil plays in preventing the successful formation of a lamella following snap-off near a pore throat. Copyright © 2015 Elsevier B.V. All rights reserved.

Multi-scale characterization of pore evolution in a combustion metamorphic complex, Hatrurim basin, Israel: Combining (ultra) small-angle neutron scattering and image analysis

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wang, Hsiu-Wen; Anovitz, Lawrence; Burg, Avihu

Backscattered scanning electron micrograph and ultra small- and small-angle neutron scattering data have been combined to provide statistically meaningful data on the pore/grain structure and pore evolution of combustion metamorphic complexes from the Hatrurim basin, Israel. Three processes, anti-sintering roughening, alteration of protolith (dehydration, decarbonation, and oxidation) and crystallization of high-temperature minerals, occurred simultaneously, leading to significant changes in observed pore/grain structures. Pore structures in the protoliths, and in lowand high-grade metamorphic rocks show surface (Ds) and mass (Dm) pore fractal geometries with gradual increases in both Ds and Dm values as a function of metamorphic grade. This suggests thatmore » increases in pore volume and formation of less branching pore networks are accompanied by a roughening of pore/grain interfaces. Additionally, pore evolution during combustion metamorphism is also characterized by reduced contributions from small-scale pores to the cumulative porosity in the high-grade rocks. At high temperatures, small-scale pores may be preferentially closed by the formation of high-temperature minerals, producing a rougher morphology with increasing temperature. Alternatively, large-scale pores may develop at the expense of small-scale pores. These observations (pore fractal geometry and cumulative porosity) indicate that the evolution of pore/grain structures is correlated with the growth of high-temperature phases and is a consequence of the energy balance between pore/grain surface energy and energy arising from heterogeneous phase contacts. The apparent pore volume density further suggests that the localized time/temperature development of the high-grade Hatrurim rocks is not simply an extension of that of the low-grade rocks. The former likely represents the "hot spots (burning foci)" in the overall metamorphic terrain while the latter may represent contact aureoles.« less

Persistent Homology fingerprinting of microstructural controls on larger-scale fluid flow in porous media

NASA Astrophysics Data System (ADS)

Moon, C.; Mitchell, S. A.; Callor, N.; Dewers, T. A.; Heath, J. E.; Yoon, H.; Conner, G. R.

2017-12-01

Traditional subsurface continuum multiphysics models include useful yet limiting geometrical assumptions: penny- or disc-shaped cracks, spherical or elliptical pores, bundles of capillary tubes, cubic law fracture permeability, etc. Each physics (flow, transport, mechanics) uses constitutive models with an increasing number of fit parameters that pertain to the microporous structure of the rock, but bear no inter-physics relationships or self-consistency. Recent advances in digital rock physics and pore-scale modeling link complex physics to detailed pore-level geometries, but measures for upscaling are somewhat unsatisfactory and come at a high computational cost. Continuum mechanics rely on a separation between small scale pore fluctuations and larger scale heterogeneity (and perhaps anisotropy), but this can break down (particularly for shales). Algebraic topology offers powerful mathematical tools for describing a local-to-global structure of shapes. Persistent homology, in particular, analyzes the dynamics of topological features and summarizes into numeric values. It offers a roadmap to both "fingerprint" topologies of pore structure and multiscale connectedness as well as links pore structure to physical behavior, thus potentially providing a means to relate the dependence of constitutive behaviors of pore structures in a self-consistent way. We present a persistence homology (PH) analysis framework of 3D image sets including a focused ion beam-scanning electron microscopy data set of the Selma Chalk. We extract structural characteristics of sampling volumes via persistence homology and fit a statistical model using the summarized values to estimate porosity, permeability, and connectivity—Lattice Boltzmann methods for single phase flow modeling are used to obtain the relationships. These PH methods allow for prediction of geophysical properties based on the geometry and connectivity in a computationally efficient way. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia LLC, a wholly owned subsidiary of Honeywell International Inc. for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

Generalized network modeling of capillary-dominated two-phase flow

NASA Astrophysics Data System (ADS)

Raeini, Ali Q.; Bijeljic, Branko; Blunt, Martin J.

2018-02-01

We present a generalized network model for simulating capillary-dominated two-phase flow through porous media at the pore scale. Three-dimensional images of the pore space are discretized using a generalized network—described in a companion paper [A. Q. Raeini, B. Bijeljic, and M. J. Blunt, Phys. Rev. E 96, 013312 (2017), 10.1103/PhysRevE.96.013312]—which comprises pores that are divided into smaller elements called half-throats and subsequently into corners. Half-throats define the connectivity of the network at the coarsest level, connecting each pore to half-throats of its neighboring pores from their narrower ends, while corners define the connectivity of pore crevices. The corners are discretized at different levels for accurate calculation of entry pressures, fluid volumes, and flow conductivities that are obtained using direct simulation of flow on the underlying image. This paper discusses the two-phase flow model that is used to compute the averaged flow properties of the generalized network, including relative permeability and capillary pressure. We validate the model using direct finite-volume two-phase flow simulations on synthetic geometries, and then present a comparison of the model predictions with a conventional pore-network model and experimental measurements of relative permeability in the literature.

Pore-scale simulation of CO2-water-rock interactions

NASA Astrophysics Data System (ADS)

Deng, H.; Molins, S.; Steefel, C. I.; DePaolo, D. J.

2017-12-01

In Geologic Carbon Storage (GCS) systems, the migration of scCO2 versus CO2-acidifed brine ultimately determines the extent of mineral trapping and caprock integrity, i.e. the long-term storage efficiency and security. While continuum scale multiphase reactive transport models are valuable for large scale investigations, they typically (over-)simplify pore-scale dynamics and cannot capture local heterogeneities that may be important. Therefore, pore-scale models are needed in order to provide mechanistic understanding of how fine scale structural variations and heterogeneous processes influence the transport and geochemistry in the context of multiphase flow, and to inform parameterization of continuum scale modeling. In this study, we investigate the interplay of different processes at pore scale (e.g. diffusion, reactions, and multiphase flow) through the coupling of a well-developed multiphase flow simulator with a sophisticated reactive transport code. The objectives are to understand where brine displaced by scCO2 will reside in a rough pore/fracture, and how the CO2-water-rock interactions may affect the redistribution of different phases. In addition, the coupled code will provide a platform for model testing in pore-scale multiphase reactive transport problems.

Pore-Scale Simulation and Sensitivity Analysis of Apparent Gas Permeability in Shale Matrix

PubMed Central

Zhang, Pengwei; Hu, Liming; Meegoda, Jay N.

2017-01-01

Extremely low permeability due to nano-scale pores is a distinctive feature of gas transport in a shale matrix. The permeability of shale depends on pore pressure, porosity, pore throat size and gas type. The pore network model is a practical way to explain the macro flow behavior of porous media from a microscopic point of view. In this research, gas flow in a shale matrix is simulated using a previously developed three-dimensional pore network model that includes typical bimodal pore size distribution, anisotropy and low connectivity of the pore structure in shale. The apparent gas permeability of shale matrix was calculated under different reservoir pressures corresponding to different gas exploitation stages. Results indicate that gas permeability is strongly related to reservoir gas pressure, and hence the apparent permeability is not a unique value during the shale gas exploitation, and simulations suggested that a constant permeability for continuum-scale simulation is not accurate. Hence, the reservoir pressures of different shale gas exploitations should be considered. In addition, a sensitivity analysis was also performed to determine the contributions to apparent permeability of a shale matrix from petro-physical properties of shale such as pore throat size and porosity. Finally, the impact of connectivity of nano-scale pores on shale gas flux was analyzed. These results would provide an insight into understanding nano/micro scale flows of shale gas in the shale matrix. PMID:28772465

Pore-Scale Simulation and Sensitivity Analysis of Apparent Gas Permeability in Shale Matrix.

PubMed

Zhang, Pengwei; Hu, Liming; Meegoda, Jay N

2017-01-25

Extremely low permeability due to nano-scale pores is a distinctive feature of gas transport in a shale matrix. The permeability of shale depends on pore pressure, porosity, pore throat size and gas type. The pore network model is a practical way to explain the macro flow behavior of porous media from a microscopic point of view. In this research, gas flow in a shale matrix is simulated using a previously developed three-dimensional pore network model that includes typical bimodal pore size distribution, anisotropy and low connectivity of the pore structure in shale. The apparent gas permeability of shale matrix was calculated under different reservoir pressures corresponding to different gas exploitation stages. Results indicate that gas permeability is strongly related to reservoir gas pressure, and hence the apparent permeability is not a unique value during the shale gas exploitation, and simulations suggested that a constant permeability for continuum-scale simulation is not accurate. Hence, the reservoir pressures of different shale gas exploitations should be considered. In addition, a sensitivity analysis was also performed to determine the contributions to apparent permeability of a shale matrix from petro-physical properties of shale such as pore throat size and porosity. Finally, the impact of connectivity of nano-scale pores on shale gas flux was analyzed. These results would provide an insight into understanding nano/micro scale flows of shale gas in the shale matrix.

Numerical simulation on hydromechanical coupling in porous media adopting three-dimensional pore-scale model.

PubMed

Liu, Jianjun; Song, Rui; Cui, Mengmeng

2014-01-01

A novel approach of simulating hydromechanical coupling in pore-scale models of porous media is presented in this paper. Parameters of the sandstone samples, such as the stress-strain curve, Poisson's ratio, and permeability under different pore pressure and confining pressure, are tested in laboratory scale. The micro-CT scanner is employed to scan the samples for three-dimensional images, as input to construct the model. Accordingly, four physical models possessing the same pore and rock matrix characteristics as the natural sandstones are developed. Based on the micro-CT images, the three-dimensional finite element models of both rock matrix and pore space are established by MIMICS and ICEM software platform. Navier-Stokes equation and elastic constitutive equation are used as the mathematical model for simulation. A hydromechanical coupling analysis in pore-scale finite element model of porous media is simulated by ANSYS and CFX software. Hereby, permeability of sandstone samples under different pore pressure and confining pressure has been predicted. The simulation results agree well with the benchmark data. Through reproducing its stress state underground, the prediction accuracy of the porous rock permeability in pore-scale simulation is promoted. Consequently, the effects of pore pressure and confining pressure on permeability are revealed from the microscopic view.

Numerical Simulation on Hydromechanical Coupling in Porous Media Adopting Three-Dimensional Pore-Scale Model

PubMed Central

Liu, Jianjun; Song, Rui; Cui, Mengmeng

2014-01-01

A novel approach of simulating hydromechanical coupling in pore-scale models of porous media is presented in this paper. Parameters of the sandstone samples, such as the stress-strain curve, Poisson's ratio, and permeability under different pore pressure and confining pressure, are tested in laboratory scale. The micro-CT scanner is employed to scan the samples for three-dimensional images, as input to construct the model. Accordingly, four physical models possessing the same pore and rock matrix characteristics as the natural sandstones are developed. Based on the micro-CT images, the three-dimensional finite element models of both rock matrix and pore space are established by MIMICS and ICEM software platform. Navier-Stokes equation and elastic constitutive equation are used as the mathematical model for simulation. A hydromechanical coupling analysis in pore-scale finite element model of porous media is simulated by ANSYS and CFX software. Hereby, permeability of sandstone samples under different pore pressure and confining pressure has been predicted. The simulation results agree well with the benchmark data. Through reproducing its stress state underground, the prediction accuracy of the porous rock permeability in pore-scale simulation is promoted. Consequently, the effects of pore pressure and confining pressure on permeability are revealed from the microscopic view. PMID:24955384

Effective pore-scale dispersion upscaling with a correlated continuous time random walk approach

NASA Astrophysics Data System (ADS)

Le Borgne, T.; Bolster, D.; Dentz, M.; de Anna, P.; Tartakovsky, A.

2011-12-01

We investigate the upscaling of dispersion from a pore-scale analysis of Lagrangian velocities. A key challenge in the upscaling procedure is to relate the temporal evolution of spreading to the pore-scale velocity field properties. We test the hypothesis that one can represent Lagrangian velocities at the pore scale as a Markov process in space. The resulting effective transport model is a continuous time random walk (CTRW) characterized by a correlated random time increment, here denoted as correlated CTRW. We consider a simplified sinusoidal wavy channel model as well as a more complex heterogeneous pore space. For both systems, the predictions of the correlated CTRW model, with parameters defined from the velocity field properties (both distribution and correlation), are found to be in good agreement with results from direct pore-scale simulations over preasymptotic and asymptotic times. In this framework, the nontrivial dependence of dispersion on the pore boundary fluctuations is shown to be related to the competition between distribution and correlation effects. In particular, explicit inclusion of spatial velocity correlation in the effective CTRW model is found to be important to represent incomplete mixing in the pore throats.

Quantitative analysis of nano-pore geomaterials and representative sampling for digital rock physics

NASA Astrophysics Data System (ADS)

Yoon, H.; Dewers, T. A.

2014-12-01

Geomaterials containing nano-pores (e.g., shales and carbonate rocks) have become increasingly important for emerging problems such as unconventional gas and oil resources, enhanced oil recovery, and geologic storage of CO2. Accurate prediction of coupled geophysical and chemical processes at the pore scale requires realistic representation of pore structure and topology. This is especially true for chalk materials, where pore networks are small and complex, and require characterization at sub-micron scale. In this work, we apply laser scanning confocal microscopy to characterize pore structures and microlithofacies at micron- and greater scales and dual focused ion beam-scanning electron microscopy (FIB-SEM) for 3D imaging of nanometer-to-micron scale microcracks and pore distributions. With imaging techniques advanced for nano-pore characterization, a problem of scale with FIB-SEM images is how to take nanometer scale information and apply it to the thin-section or larger scale. In this work, several texture characterization techniques including graph-based spectral segmentation, support vector machine, and principal component analysis are applied for segmentation clusters represented by 1-2 FIB-SEM samples per each cluster. Geometric and topological properties are analyzed and lattice-Boltzmann method (LBM) is used to obtain permeability at several different scales. Upscaling of permeability to the Darcy scale (e.g., the thin-section scale) with image dataset will be discussed with emphasis on understanding microfracture-matrix interaction, representative volume for FIB-SEM sampling, and multiphase flow and reactive transport. Funding from the DOE Basic Energy Sciences Geosciences Program is gratefully acknowledged. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

Dynamic Pore-Scale Imaging of Reactive Transport in Heterogeneous Carbonates at Reservior Conditions

NASA Astrophysics Data System (ADS)

Menke, Hannah; Bijeljic, Branko; Andrew, Matthew; Blunt, Martin

2014-05-01

Sequestering carbon in deep geologic formations is one way of reducing anthropogenic CO2 emissions. Carbon capture, Utilization, and Storage (CCUS) in carbonate reservoirs has the added benefit of mobilizing more oil for extraction, increasing oil reservoir yield, and generating revenue while also mitigating climate change. The magnitude, speed, and type of dissolution are dependent the intrinsic properties of the rock. Understanding how small changes in the pore structure affect dissolution is paramount for successful predictive modelling both on the pore-scale and for up-scaled reservoir simulations. We propose an experimental method whereby both 'Pink Beam' synchrotron radiation and a Micro-CT lab source are used in dynamic X-ray microtomography to investigate the pore structure changes in carbonate rocks of varying heterogeneity at high temperatures and pressures. Four carbonate rock types were studied, two relatively homogeneous carbonates, Ketton and Mt. Gambier, and two very heterogeneous carbonates, Estalliades and Portland Basebed. Each rock type was imaged under the same reservoir and flow conditions to gain insight into the impact of heterogeneity. A 4-mm carbonate core was injected with CO2-saturated brine at 10 MPa and 50oC for 2 hours. Depending on sample heterogeneity and X-ray source, tomographic images were taken at between 30-second and 20-minute time-resolutions and a 4-micron spatial resolution during injection. Changes in porosity, permeability, and structure were obtained by first binning and filtering the images, then binarizing them with watershed segmentation, and finally extracting a pore/throat network. Furthermore, pore-scale flow modelling was performed directly on the binarized image and used to track velocity distributions as the pore network evolved. Significant differences in dissolution type and magnitude were found for each rock type. The most homogeneous carbonate, Ketton, was seen to have predominately uniform dissolution with minor dissolution rate differences between the pores and pore throats. This was not true for the heterogeneous carbonates, Estalliades and Portland Basebed, which formed wormholes. Pore-scale modelling of flow directly on the voxels showed the differences in the evolution of complex flow fields with changes in dissolution regime. The PDFs of normalized velocity for uniform dissolution showed that the maximum pore velocity within the system decreased as dissolution occurred. This is due to dissolution enlarging pores and pore throats. However, in the wormholing regime, there was a large increase in maximum velocity once the wormhole broke through the length of the core and a preferential flow path was created. Additionally, this study serves as a unique benchmark for pore-scale reactive transport modelling directly on the binarized Micro-CT images. This dynamic pore-scale imaging method offers advantages in helping fully explain the dominant physical and chemical processes at the pore scale so that they may be up-scaled to the reservoir scale for increased accuracy in model prediction.

Micro- and nano-X-ray computed-tomography: A step forward in the characterization of the pore network of a leached cement paste

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bossa, Nathan, E-mail: bossanathan@gmail.com; INERIS, Parc Technologique Alata, BP2, 60550 Verneuil-en-Halatte; iCEINT, CNRS, Duke Univ. International Consortium for the Environmental Implications of Nanotechnology, Aix-en-Provence

2015-01-15

Pore structure of leached cement pastes (w/c = 0.5) was studied for the first time from micro-scale down to the nano-scale by combining micro- and nano-X-ray computed tomography (micro- and nano-CT). This allowed assessing the 3D heterogeneity of the pore network along the cement profile (from the core to the altered layer) of almost the entire range of cement pore size, i.e. from capillary to gel pores. We successfully quantified an increase of porosity in the altered layer at both resolutions. Porosity is increasing from 1.8 to 6.1% and from 18 to 58% at the micro-(voxel = 1.81 μm) andmore » nano-scale (voxel = 63.5 nm) respectively. The combination of both CT allowed to circumvent weaknesses inherent of both investigation scales. In addition the connectivity and the channel size of the pore network were also evaluated to obtain a complete 3D pore network characterization at both scales.« less

From Pore to Core: Do Engineered Nanoparticles Violate Upscaling Assumptions? A Microtomographic Investigation

NASA Astrophysics Data System (ADS)

Molnar, I. L.; O'Carroll, D. M.; Gerhard, J.; Willson, C. S.

2014-12-01

The recent success in using Synchrotron X-ray Computed Microtomography (SXCMT) for the quantification of nanoparticle concentrations within real, three-dimensional pore networks [1] has opened up new opportunities for collecting experimental data of pore-scale flow and transport processes. One opportunity is coupling SXCMT with nanoparticle/soil transport experiments to provide unique insights into how pore-scale processes influence transport at larger scales. Understanding these processes is a key step in accurately upscaling micron-scale phenomena to the continuum-scale. Upscaling phenomena from the micron-scale to the continuum-scale typically involves the assumption that the pore space is well mixed. Using this 'well mixed assumption' it is implicitly assumed that the distribution of nanoparticles within the pore does not affect its retention by soil grains. This assumption enables the use of volume-averaged parameters in calculating transport and retention rates. However, in some scenarios, the well mixed assumption will likely be violated by processes such as deposition and diffusion. These processes can alter the distribution of the nanoparticles in the pore space and impact retention behaviour, leading to discrepancies between theoretical predictions and experimental observations. This work investigates the well mixed assumption by employing SXCMT to experimentally examine pore-scale mixing of silver nanoparticles during transport through sand packed columns. Silver nanoparticles were flushed through three different sands to examine the impact of grain distribution and nanoparticle retention rates on mixing: uniform silica (low retention), well graded silica sand (low retention) and uniform iron oxide coated silica sand (high retention). The SXCMT data identified diffusion-limited retention as responsible for violations of the well mixed assumption. A mathematical description of the diffusion-limited retention process was created and compared to the experimental data at the pore and column-scale. The mathematical description accurately predicted trends observed within the SXCMT-datasets such as concentration gradients away from grain surfaces and also accurately predicted total retention of nanoparticles at the column scale. 1. ES&T 2014, 48, (2), 1114-1122.

Toward direct pore-scale modeling of three-phase displacements

NASA Astrophysics Data System (ADS)

Mohammadmoradi, Peyman; Kantzas, Apostolos

2017-12-01

A stable spreading film between water and gas can extract a significant amount of bypassed non-aqueous phase liquid (NAPL) through immiscible three-phase gas/water injection cycles. In this study, the pore-scale displacement mechanisms by which NAPL is mobilized are incorporated into a three-dimensional pore morphology-based model under water-wet and capillary equilibrium conditions. The approach is pixel-based and the sequence of invasions is determined by the fluids' connectivity and the threshold capillary pressure of the advancing interfaces. In addition to the determination of three-phase spatial saturation profiles, residuals, and capillary pressure curves, dynamic finite element simulations are utilized to predict the effective permeabilities of the rock microtomographic images as reasonable representations of the geological formations under study. All the influential features during immiscible fluid flow in pore-level domains including wetting and spreading films, saturation hysteresis, capillary trapping, connectivity, and interface development strategies are taken into account. The capabilities of the model are demonstrated by the successful prediction of saturation functions for Berea sandstone and the accurate reconstruction of three-phase fluid occupancies through a micromodel.

Intercomparison of 3D pore-scale flow and solute transport simulation methods

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yang, Xiaofan; Mehmani, Yashar; Perkins, William A.

2016-09-01

Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include methods that 1) explicitly model the three-dimensional geometry of pore spaces and 2) those that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of class 1, based on direct numerical simulation using computational fluid dynamics (CFD) codes, against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of class 1 based on the immersed-boundary method (IMB),more » lattice Boltzmann method (LBM), smoothed particle hydrodynamics (SPH), as well as a model of class 2 (a pore-network model or PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and nonreactive solute transport, and intercompare the model results with previously reported experimental observations. Experimental observations are limited to measured pore-scale velocities, so solute transport comparisons are made only among the various models. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations).« less

Intercomparison of 3D pore-scale flow and solute transport simulation methods

DOE PAGES

Mehmani, Yashar; Schoenherr, Martin; Pasquali, Andrea; ...

2015-09-28

Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include 1) methods that explicitly model the three-dimensional geometry of pore spaces and 2) methods that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of the first type, using computational fluid dynamics (CFD) codes employing a standard finite volume method (FVM), against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of the first type based onmore » the lattice Boltzmann method (LBM) and smoothed particle hydrodynamics (SPH), as well as a model of the second type, a pore-network model (PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (FVM-based CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and (for capable codes) nonreactive solute transport, and intercompare the model results. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations). Generally good agreement was achieved among the various approaches, but some differences were observed depending on the model context. The intercomparison work was challenging because of variable capabilities of the codes, and inspired some code enhancements to allow consistent comparison of flow and transport simulations across the full suite of methods. This paper provides support for confidence in a variety of pore-scale modeling methods and motivates further development and application of pore-scale simulation methods.« less

Intercomparison of 3D pore-scale flow and solute transport simulation methods

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yang, Xiaofan; Mehmani, Yashar; Perkins, William A.

2016-09-01

Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include 1) methods that explicitly model the three-dimensional geometry of pore spaces and 2) methods that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of the first type, using computational fluid dynamics (CFD) codes employing a standard finite volume method (FVM), against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of the first type based onmore » the lattice Boltzmann method (LBM) and smoothed particle hydrodynamics (SPH), as well as a model of the second type, a pore-network model (PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (FVM-based CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and (for capable codes) nonreactive solute transport, and intercompare the model results. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations). Generally good agreement was achieved among the various approaches, but some differences were observed depending on the model context. The intercomparison work was challenging because of variable capabilities of the codes, and inspired some code enhancements to allow consistent comparison of flow and transport simulations across the full suite of methods. This study provides support for confidence in a variety of pore-scale modeling methods and motivates further development and application of pore-scale simulation methods.« less

Final Report for DE-SC0008059

DOE Office of Scientific and Technical Information (OSTI.GOV)

Karnik, Rohit

Graphene, a one-atom thick material comprising hexagonally bonded carbon atoms, is one of the strongest materials known and is impermeable to even helium gas. The potential for creation of tunable nanometer-scale pores in graphene, combined with mechanical strength, chemical resistance, and atomic thickness, make it a promising material for improving selectivity, permeability, and energy efficiency in a diverse range of membrane separations. However, fundamental understanding of creation of porous graphene, its mass transport properties, and the relationship between porosity and membrane performance is lacking. This project studied methods to create pores in graphene by ion irradiation and etching, and investigatedmore » the effect of the porosity of graphene on mass transport properties at the membrane level as well as across individual pores. Detailed characterization of pore size distributions in graphene using high-resolution scanning transmission electron microscopy was performed to elucidate the relationship between different pore creation processes and the resulting porosity. New methods to create tunable nanometer-scale pores in graphene at high density over large areas were developed. Graphene membranes were fabricated comprising single-layer porous graphene on appropriate support membranes. Mass transport of ions, molecules, water, and gases across the membranes was investigated, and predictive models were developed to relate the graphene porosity and support membrane structure to the mass transport properties of the membranes. The effect of graphene porosity on permeance of the membranes to water, salts, and small molecules was experimentally quantified. Single-pore measurements were performed and a model was developed to understand ion transport across graphene pores, which led to the experimental observation of heterogeneity in transport behaviors and spontaneous fluctuations of ion transport. Functional single-layer nanoporous graphene membranes comprising a high density of pores over macro-scale areas with molecular selectivity were demonstrated for the first time. This project led to three patent applications and 9 journal articles. The fundamental developments in understanding nanoporous graphene membranes have laid the foundations to develop energy-efficient nanoporous graphene membranes for various applications in water desalination, oil and natural gas separations, biotechnology, and chemical processing.« less

Pore-scale dynamics of enzyme adsorption, swelling and reactive dissolution determine sugar yield in hemicellulose hydrolysis for biofuel production

NASA Astrophysics Data System (ADS)

Dutta, Sajal Kanti; Chakraborty, Saikat

2016-12-01

Hemicelluloses are the earth’s second most abundant structural polymers, found in lignocellulosic biomass. Efficient enzymatic depolymerization of xylans by cleaving their β-(1 → 4)-glycosidic bonds to produce soluble sugars is instrumental to the cost-effective production of liquid biofuels. Here we show that the multi-scale two-phase process of enzymatic hydrolysis of amorphous hemicelluloses is dominated by its smallest scale-the pores. In the crucial first five hours, two to fourfold swelling of the xylan particles allow the enzymes to enter the pores and undergo rapid non-equilibrium adsorption on the pore surface before they hydrolyze the solid polymers, albeit non-competitively inhibited by the products xylose and xylobiose. Rapid pore-scale reactive dissolution increases the solid carbohydrate’s porosity to 80-90%. This tightly coupled experimental and theoretical study quantifies the complex temporal dynamics of the transport and reaction processes coupled across scales and phases to show that this unique pore-scale phenomenon can be exploited to accelerate the depolymerization of hemicelluloses to monomeric sugars in the first 5-6 h. We find that an ‘optimal substrate loading’ of 5 mg/ml (above which substrate inhibition sets in) accelerates non-equilibrium enzyme adsorption and solid hemicellulose depolymerization at the pore-scale, which contributes three-quarters of the soluble sugars produced for bio-alcohol fermentation.

Dynamic Pore-Scale Imaging of Reactive Transport in Heterogeneous Carbonates at Reservoir Conditions Across Multiple Dissolution Regimes

NASA Astrophysics Data System (ADS)

Menke, H. P.; Bijeljic, B.; Andrew, M. G.; Blunt, M. J.

2014-12-01

Sequestering carbon in deep geologic formations is one way of reducing anthropogenic CO2 emissions. When supercritical CO2 mixes with brine in a reservoir, the acid generated has the potential to dissolve the surrounding pore structure. However, the magnitude and type of dissolution are condition dependent. Understanding how small changes in the pore structure, chemistry, and flow properties affect dissolution is paramount for successful predictive modelling. Both 'Pink Beam' synchrotron radiation and a Micro-CT lab source are used in dynamic X-ray microtomography to investigate the pore structure changes during supercritical CO2 injection in carbonate rocks of varying heterogeneity at high temperatures and pressures and various flow-rates. Three carbonate rock types were studied, one with a homogeneous pore structure and two heterogeneous carbonates. All samples are practically pure calcium carbonate, but have widely varying rock structures. Flow-rate was varied in three successive experiments by over an order of magnitude whlie keeping all other experimental conditions constant. A 4-mm carbonate core was injected with CO2-saturated brine at 10 MPa and 50oC. Tomographic images were taken at 30-second to 20-minute time-resolutions during a 2 to 4-hour injection period. A pore network was extracted using a topological analysis of the pore space and pore-scale flow modelling was performed directly on the binarized images with connected pathways and used to track the altering velocity distributions. Significant differences in dissolution type and magnitude were found for each rock type and flowrate. At the highest flow-rates, the homogeneous carbonate was seen to have predominately uniform dissolution with minor dissolution rate differences between the pores and pore throats. Alternatively, the heterogeneous carbonates which formed wormholes at high flow rates. At low flow rates the homogeneous rock developed wormholes, while the heterogeneous samples showed evidence of compact dissolution. This study serves as a unique benchmark for pore-scale reactive transport modelling directly on the binarized Micro-CT images. Dynamic pore-scale imaging methods offer advantages in helping explain the dominant processes at the pore scale so that they may be up-scaled for accurate model prediction.

Nanometer-Scale Pore Characteristics of Lacustrine Shale, Songliao Basin, NE China

PubMed Central

Wang, Min; Yang, Jinxiu; Wang, Zhiwei; Lu, Shuangfang

2015-01-01

In shale, liquid hydrocarbons are accumulated mainly in nanometer-scale pores or fractures, so the pore types and PSDs (pore size distributions) play a major role in the shale oil occurrence (free or absorbed state), amount of oil, and flow features. The pore types and PSDs of marine shale have been well studied; however, research on lacustrine shale is rare, especially for shale in the oil generation window, although lacustrine shale is deposited widely around the world. To investigate the relationship between nanometer-scale pores and oil occurrence in the lacustrine shale, 10 lacustrine shale core samples from Songliao Basin, NE China were analyzed. Analyses of these samples included geochemical measurements, SEM (scanning electron microscope) observations, low pressure CO2 and N2 adsorption, and high-pressure mercury injection experiments. Analysis results indicate that: (1) Pore types in the lacustrine shale include inter-matrix pores, intergranular pores, organic matter pores, and dissolution pores, and these pores are dominated by mesopores and micropores; (2) There is no apparent correlation between pore volumes and clay content, however, a weak negative correlation is present between total pore volume and carbonate content; (3) Pores in lacustrine shale are well developed when the organic matter maturity (Ro) is >1.0% and the pore volume is positively correlated with the TOC (total organic carbon) content. The statistical results suggest that oil in lacustrine shale mainly occurs in pores with diameters larger than 40 nm. However, more research is needed to determine whether this minimum pore diameter for oil occurrence in lacustrine shale is widely applicable. PMID:26285123

Nanometer-Scale Pore Characteristics of Lacustrine Shale, Songliao Basin, NE China.

PubMed

Wang, Min; Yang, Jinxiu; Wang, Zhiwei; Lu, Shuangfang

2015-01-01

In shale, liquid hydrocarbons are accumulated mainly in nanometer-scale pores or fractures, so the pore types and PSDs (pore size distributions) play a major role in the shale oil occurrence (free or absorbed state), amount of oil, and flow features. The pore types and PSDs of marine shale have been well studied; however, research on lacustrine shale is rare, especially for shale in the oil generation window, although lacustrine shale is deposited widely around the world. To investigate the relationship between nanometer-scale pores and oil occurrence in the lacustrine shale, 10 lacustrine shale core samples from Songliao Basin, NE China were analyzed. Analyses of these samples included geochemical measurements, SEM (scanning electron microscope) observations, low pressure CO2 and N2 adsorption, and high-pressure mercury injection experiments. Analysis results indicate that: (1) Pore types in the lacustrine shale include inter-matrix pores, intergranular pores, organic matter pores, and dissolution pores, and these pores are dominated by mesopores and micropores; (2) There is no apparent correlation between pore volumes and clay content, however, a weak negative correlation is present between total pore volume and carbonate content; (3) Pores in lacustrine shale are well developed when the organic matter maturity (Ro) is >1.0% and the pore volume is positively correlated with the TOC (total organic carbon) content. The statistical results suggest that oil in lacustrine shale mainly occurs in pores with diameters larger than 40 nm. However, more research is needed to determine whether this minimum pore diameter for oil occurrence in lacustrine shale is widely applicable.

Using Multiscale Modeling to Study Coupled Flow, Transport, Reaction and Biofilm Growth Processes in Porous Media

NASA Astrophysics Data System (ADS)

Valocchi, A. J.; Laleian, A.; Werth, C. J.

2017-12-01

Perturbation of natural subsurface systems by fluid inputs may induce geochemical or microbiological reactions that change porosity and permeability, leading to complex coupled feedbacks between reaction and transport processes. Some examples are precipitation/dissolution processes associated with carbon capture and storage and biofilm growth associated with contaminant transport and remediation. We study biofilm growth due to mixing controlled reaction of multiple substrates. As biofilms grow, pore clogging occurs which alters pore-scale flow paths thus changing the mixing and reaction. These interactions are challenging to quantify using conventional continuum-scale porosity-permeability relations. Pore-scale models can accurately resolve coupled reaction, biofilm growth and transport processes, but modeling at this scale is not feasible for practical applications. There are two approaches to address this challenge. Results from pore-scale models in generic pore structures can be used to develop empirical relations between porosity and continuum-scale parameters, such as permeability and dispersion coefficients. The other approach is to develop a multiscale model of biofilm growth in which non-overlapping regions at pore and continuum spatial scales are coupled by a suitable method that ensures continuity of flux across the interface. Thus, regions of high reactivity where flow alteration occurs are resolved at the pore scale for accuracy while regions of low reactivity are resolved at the continuum scale for efficiency. This approach thus avoids the need for empirical upscaling relations in regions with strong feedbacks between reaction and porosity change. We explore and compare these approaches for several two-dimensional cases.

Nonaqueous Phase Liquid Dissolution in Porous Media: Multi-Scale Effects of Multi-Component Dissolution Kinetics on Cleanup Time

DOE Office of Scientific and Technical Information (OSTI.GOV)

McNab, W; Ezzedine, S; Detwiler, R

2007-02-26

Industrial organic solvents such as trichloroethylene (TCE) and tetrachloroethylene (PCE) constitute a principal class of groundwater contaminants. Cleanup of groundwater plume source areas associated with these compounds is problematic, in part, because the compounds often exist in the subsurface as dense nonaqueous phase liquids (DNAPLs). Ganglia (or 'blobs') of DNAPL serve as persistent sources of contaminants that are difficult to locate and remediate (e.g. Fenwick and Blunt, 1998). Current understanding of the physical and chemical processes associated with dissolution of DNAPLs in the subsurface is incomplete and yet is critical for evaluating long-term behavior of contaminant migration, groundwater cleanup, andmore » the efficacy of source area cleanup technologies. As such, a goal of this project has been to contribute to this critical understanding by investigating the multi-phase, multi-component physics of DNAPL dissolution using state-of-the-art experimental and computational techniques. Through this research, we have explored efficient and accurate conceptual and numerical models for source area contaminant transport that can be used to better inform the modeling of source area contaminants, including those at the LLNL Superfund sites, to re-evaluate existing remediation technologies, and to inspire or develop new remediation strategies. The problem of DNAPL dissolution in natural porous media must be viewed in the context of several scales (Khachikian and Harmon, 2000), including the microscopic level at which capillary forces, viscous forces, and gravity/buoyancy forces are manifested at the scale of individual pores (Wilson and Conrad, 1984; Chatzis et al., 1988), the mesoscale where dissolution rates are strongly influenced by the local hydrodynamics, and the field-scale. Historically, the physico-chemical processes associated with DNAPL dissolution have been addressed through the use of lumped mass transfer coefficients which attempt to quantify the dissolution rate in response to local dissolved-phase concentrations distributed across the source area using a volume-averaging approach (Figure 1). The fundamental problem with the lumped mass transfer parameter is that its value is typically derived empirically through column-scale experiments that combine the effects of pore-scale flow, diffusion, and pore-scale geometry in a manner that does not provide a robust theoretical basis for upscaling. In our view, upscaling processes from the pore-scale to the field-scale requires new computational approaches (Held and Celia, 2001) that are directly linked to experimental studies of dissolution at the pore scale. As such, our investigation has been multi-pronged, combining theory, experiments, numerical modeling, new data analysis approaches, and a synthesis of previous studies (e.g. Glass et al, 2001; Keller et al., 2002) aimed at quantifying how the mechanisms controlling dissolution at the pore-scale control the long-term dissolution of source areas at larger scales.« less

Estimating Pore Properties from NMR Relaxation Time Measurements in Heterogeneous Media

NASA Astrophysics Data System (ADS)

Grunewald, E.; Knight, R.

2008-12-01

The link between pore geometry and the nuclear magnetic resonance (NMR) relaxation time T2 is well- established for simple systems but is poorly understood for complex media with heterogeneous pores. Conventional interpretation of NMR relaxation data employs a model of isolated pores in which each hydrogen proton samples only one pore type, and the T2-distribution is directly scaled to estimate a pore-size distribution. During an actual NMR measurement, however, each proton diffuses through a finite volume of the pore network, and so may sample multiple pore types encountered within this diffusion cell. For cases in which heterogeneous pores are strongly coupled by diffusion, the meaning of the T2- distribution is not well understood and further research is required to determine how such measurements should be interpreted. In this study we directly investigate the implications of pore coupling in two groups of laboratory NMR experiments. We conduct two suites of experiments, in which samples are synthesized to exhibit a range of pore coupling strengths using two independent approaches: (a) varying the scale of the diffusion cell and (b) varying the scale over which heterogeneous pores are encountered. In the first set of experiments, we vary the scale of the diffusion cell in silica gels which have a bimodal pore-size distribution comprised of intragrannular micropores and much larger intergrannular pores. The untreated gel exhibits strong pore coupling with a single broad peak observed in the T2-distribution. By treating the gel with varied amounts of paramagnetic iron surface coatings, we decrease the surface relaxation time, T2S, and effectively decrease both the size of the diffusion cell and the degree of pore coupling. As more iron is coated to the grain surfaces, we observe a separation of the broad T2-distribution into two peaks that more accurately represent the true bimodal pore-size distribution. In the second set of experiments, we vary the scale over which heterogeneous pores are encountered in bimodal grain packs of pure quartz (long T2S) and hematite (short T2S). The scale of heterogeneity is varied by changing the mean grain size and relative mineral concentrations. When the mean grain size is small and the mineral concentrations are comparable, the T2-distribution is roughly monomodal indicating strong pore coupling. As the grain size is increased or the mineral concentrations are made increasingly uneven, the T2- distribution develops a bimodal character, more representative of the actual distribution of pore types. Numerical simulations of measurements in both experiment groups allow us to more closely investigate how the relaxing magnetization evolves in both time and space. Collectively, these experiments provide important insights into the effects of pore coupling on NMR measurements in heterogeneous systems and contribute to our ultimate goal of improving the interpretation of these data in complex near-surface sediments.

Application of a pore-scale reactive transport model to a natural analog for reaction-induced pore alterations

DOE PAGES

Yoon, Hongkyu; Major, Jonathan; Dewers, Thomas; ...

2017-01-05

Dissolved CO 2 in the subsurface resulting from geological CO 2 storage may react with minerals in fractured rocks, confined aquifers, or faults, resulting in mineral precipitation and dissolution. The overall rate of reaction can be affected by coupled processes including hydrodynamics, transport, and reactions at the (sub) pore-scale. In this work pore-scale modeling of coupled fluid flow, reactive transport, and heterogeneous reactions at the mineral surface is applied to account for permeability alterations caused by precipitation-induced pore-blocking. This paper is motivated by observations of CO 2 seeps from a natural CO 2 sequestration analog, Crystal Geyser, Utah. Observations alongmore » the surface exposure of the Little Grand Wash fault indicate the lateral migration of CO 2 seep sites (i.e., alteration zones) of 10–50 m width with spacing on the order of ~100 m over time. Sandstone permeability in alteration zones is reduced by 3–4 orders of magnitude by carbonate cementation compared to unaltered zones. One granular porous medium and one fracture network systems are used to conceptually represent permeable porous media and locations of conduits controlled by fault-segment intersections and/or topography, respectively. Simulation cases accounted for a range of reaction regimes characterized by the Damköhler (Da) and Peclet (Pe) numbers. Pore-scale simulation results demonstrate that combinations of transport (Pe), geochemical conditions (Da), solution chemistry, and pore and fracture configurations contributed to match key patterns observed in the field of how calcite precipitation alters flow paths by pore plugging. This comparison of simulation results with field observations reveals mechanistic explanations of the lateral migration and enhances our understanding of subsurface processes associated with the CO 2 injection. In addition, permeability and porosity relations are constructed from pore-scale simulations which account for a range of reaction regimes characterized by the Da and Pe numbers. Finally, the functional relationships obtained from pore-scale simulations can be used in a continuum scale model that may account for large-scale phenomena mimicking lateral migration of surface CO 2 seeps.« less

Investigation of Biogrout processes by numerical analysis at pore scale

NASA Astrophysics Data System (ADS)

Bergwerff, Luke; van Paassen, Leon A.; Picioreanu, Cristian; van Loosdrecht, Mark C. M.

2013-04-01

Biogrout is a soil improving process that aims to improve the strength of sandy soils. The process is based on microbially induced calcite precipitation (MICP). In this study the main process is based on denitrification facilitated by bacteria indigenous to the soil using substrates, which can be derived from pretreated waste streams containing calcium salts of fatty acids and calcium nitrate, making it a cost effective and environmentally friendly process. The goal of this research is to improve the understanding of the process by numerical analysis so that it may be improved and applied properly for varying applications, such as borehole stabilization, liquefaction prevention, levee fortification and mitigation of beach erosion. During the denitrification process there are many phases present in the pore space including a liquid phase containing solutes, crystals, bacteria forming biofilms and gas bubbles. Due to the amount of phases and their dynamic changes (multiphase flow with (non-linear) reactive transport), there are many interactions making the process very complex. To understand this complexity in the system, the interactions between these phases are studied in a reductionist approach, increasing the complexity of the system by one phase at a time. The model will initially include flow, solute transport, crystal nucleation and growth in 2D at pore scale. The flow will be described by Navier-Stokes equations. Initial study and simulations has revealed that describing crystal growth for this application on a fixed grid can introduce significant fundamental errors. Therefore a level set method will be employed to better describe the interface of developing crystals in between sand grains. Afterwards the model will be expanded to 3D to provide more realistic flow, nucleation and clogging behaviour at pore scale. Next biofilms and lastly gas bubbles may be added to the model. From the results of these pore scale models the behaviour of the system may be studied and eventually observations may be extrapolated to a larger continuum scale.

Pore-scale and Continuum Simulations of Solute Transport Micromodel Benchmark Experiments

DOE Office of Scientific and Technical Information (OSTI.GOV)

Oostrom, Martinus; Mehmani, Yashar; Romero Gomez, Pedro DJ

Four sets of micromodel nonreactive solute transport experiments were conducted with flow velocity, grain diameter, pore-aspect ratio, and flow focusing heterogeneity as the variables. The data sets were offered to pore-scale modeling groups to test their simulators. Each set consisted of two learning experiments, for which all results was made available, and a challenge experiment, for which only the experimental description and base input parameters were provided. The experimental results showed a nonlinear dependence of the dispersion coefficient on the Peclet number, a negligible effect of the pore-aspect ratio on transverse mixing, and considerably enhanced mixing due to flow focusing.more » Five pore-scale models and one continuum-scale model were used to simulate the experiments. Of the pore-scale models, two used a pore-network (PN) method, two others are based on a lattice-Boltzmann (LB) approach, and one employed a computational fluid dynamics (CFD) technique. The learning experiments were used by the PN models to modify the standard perfect mixing approach in pore bodies into approaches to simulate the observed incomplete mixing. The LB and CFD models used these experiments to appropriately discretize the grid representations. The continuum model use published non-linear relations between transverse dispersion coefficients and Peclet numbers to compute the required dispersivity input values. Comparisons between experimental and numerical results for the four challenge experiments show that all pore-scale models were all able to satisfactorily simulate the experiments. The continuum model underestimated the required dispersivity values and, resulting in less dispersion. The PN models were able to complete the simulations in a few minutes, whereas the direct models needed up to several days on supercomputers to resolve the more complex problems.« less

Pore-scale dynamics of enzyme adsorption, swelling and reactive dissolution determine sugar yield in hemicellulose hydrolysis for biofuel production

PubMed Central

Dutta, Sajal Kanti; Chakraborty, Saikat

2016-01-01

Hemicelluloses are the earth’s second most abundant structural polymers, found in lignocellulosic biomass. Efficient enzymatic depolymerization of xylans by cleaving their β-(1 → 4)-glycosidic bonds to produce soluble sugars is instrumental to the cost-effective production of liquid biofuels. Here we show that the multi-scale two-phase process of enzymatic hydrolysis of amorphous hemicelluloses is dominated by its smallest scale–the pores. In the crucial first five hours, two to fourfold swelling of the xylan particles allow the enzymes to enter the pores and undergo rapid non-equilibrium adsorption on the pore surface before they hydrolyze the solid polymers, albeit non-competitively inhibited by the products xylose and xylobiose. Rapid pore-scale reactive dissolution increases the solid carbohydrate’s porosity to 80–90%. This tightly coupled experimental and theoretical study quantifies the complex temporal dynamics of the transport and reaction processes coupled across scales and phases to show that this unique pore-scale phenomenon can be exploited to accelerate the depolymerization of hemicelluloses to monomeric sugars in the first 5–6 h. We find that an ‘optimal substrate loading’ of 5 mg/ml (above which substrate inhibition sets in) accelerates non-equilibrium enzyme adsorption and solid hemicellulose depolymerization at the pore-scale, which contributes three-quarters of the soluble sugars produced for bio-alcohol fermentation. PMID:27905534

Pore and Continuum Scale Study of the Effect of Subgrid Transport Heterogeneity on Redox Reaction Rates

DOE Office of Scientific and Technical Information (OSTI.GOV)

Liu, Yuanyuan; Liu, Chongxuan; Zhang, Changyong

2015-08-01

A micromodel system with a pore structure for heterogeneous flow and transport was used to investigate the effect of subgrid transport heterogeneity on redox reaction rates. Hematite reductive dissolution by injecting a reduced form of flavin mononucleotide (FMNH2) at variable flow rates was used as an example to probe the variations of redox reaction rates in different subgrid transport domains. Experiments, pore-scale simulations, and macroscopic modeling were performed to measure and simulate in-situ hematite reduction and to evaluate the scaling behavior of the redox reaction rates from the pore to macroscopic scales. The results indicated that the measured pore-scale ratesmore » of hematite reduction were consistent with the predictions from a pore scale reactive transport model. A general trend is that hematite reduction followed reductant transport pathways, starting from the advection-dominated pores toward the interior of diffusion-dominated domains. Two types of diffusion domains were considered in the micromodel: a micropore diffusion domain, which locates inside solid grains or aggregates where reactant transport is limited by diffusion; and a macropore diffusion domain, which locates at wedged, dead-end pore spaces created by the grain-grain contacts. The rate of hematite reduction in the advection-dominated domain was faster than those in the diffusion-controlled domains, and the rate in the macropore diffusion domain was faster than that in the micropore domain. The reduction rates in the advection and macropore diffusion domains increased with increasing flow rate, but were affected by different mechanisms. The rate increase in the advection domain was controlled by the mass action effect as a faster flow supplied more reactants, and the rate increase in the macropore domain was more affected by the rate of mass exchange with the advection domain, which increased with increasing flow rate. The hematite reduction rate in the micropore domain was, however, not affected by the flow rate because molecular diffusion limits reductant supply to the micropore domain interior. Domain-based macroscopic models were evaluated to scale redox reaction rates from the pore to macroscopic scales. A single domain model, which ignores subgrid transport heterogeneity deviated significantly from the pore-scale results. Further analysis revealed that the rate expression for hematite reduction was not scalable from the pore to porous media using the single domain model. A three-domain model, which effectively considers subgrid reactive diffusion in the micropore and macropore domains, significantly improved model description. Overall this study revealed the importance of subgrid transport heterogeneity in the manifestation of redox reaction rates in porous media and in scaling reactions from the pore to porous media. The research also supported that the domain-based scaling approach can be used to directly scale redox reactions in porous media with subgrid transport heterogeneity.« less

Multiscale pore structure and constitutive models of fine-grained rocks

NASA Astrophysics Data System (ADS)

Heath, J. E.; Dewers, T. A.; Shields, E. A.; Yoon, H.; Milliken, K. L.

2017-12-01

A foundational concept of continuum poromechanics is the representative elementary volume or REV: an amount of material large enough that pore- or grain-scale fluctuations in relevant properties are dissipated to a definable mean, but smaller than length scales of heterogeneity. We determine 2D-equivalent representative elementary areas (REAs) of pore areal fraction of three major types of mudrocks by applying multi-beam scanning electron microscopy (mSEM) to obtain terapixel image mosaics. Image analysis obtains pore areal fraction and pore size and shape as a function of progressively larger measurement areas. Using backscattering imaging and mSEM data, pores are identified by the components within which they occur, such as in organics or the clastic matrix. We correlate pore areal fraction with nano-indentation, micropillar compression, and axysimmetic testing at multiple length scales on a terrigenous-argillaceous mudrock sample. The combined data set is used to: investigate representative elementary volumes (and areas for the 2D images); determine if scale separation occurs; and determine if transport and mechanical properties at a given length scale can be statistically defined. Clear scale separation occurs between REAs and observable heterogeneity in two of the samples. A highly-laminated sample exhibits fine-scale heterogeneity and an overlapping in scales, in which case typical continuum assumptions on statistical variability may break down. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia LLC, a wholly owned subsidiary of Honeywell International Inc. for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

Pore-scale simulation of calcium carbonate precipitation and dissolution under highly supersaturated conditions in a microfludic pore network

NASA Astrophysics Data System (ADS)

Yoon, H.; Dewers, T. A.; Valocchi, A. J.; Werth, C. J.

2011-12-01

Dissolved CO2 during geological CO2 storage may react with minerals in fractured rocks or confined aquifers and cause mineral precipitation. The overall rate of reaction can be affected by coupled processes among hydrodynamics, transport, and reactions at pore-scale. Pore-scale models of coupled fluid flow, reactive transport, and CaCO3 precipitation and dissolution are applied to account for transient experimental results of CaCO3 precipitation and dissolution under highly supersaturated conditions in a microfluidic pore network (i.e., micromodel). Pore-scale experiments in the micromodel are used as a basis for understanding coupled physics of systems perturbed by geological CO2 injection. In the micromodel, precipitation is induced by transverse mixing along the centerline in pore bodies. Overall, the pore-scale model qualitatively captured the governing physics of reactions such as precipitate morphology, precipitation rate, and maximum precipitation area in first few pore spaces. In particular, we found that proper estimation of the effective diffusion coefficient and the reactive surface area is necessary to adequately simulate precipitation and dissolution rates. As the model domain increases, the effect of flow patterns affected by precipitation on the overall reaction rate also increases. The model is also applied to account for the effect of different reaction rate laws on mineral precipitation and dissolution at pore-scale. Reaction rate laws tested include the linear rate law, nonlinear power law, and newly-developed rate law based on in-situ measurements at nano scale in the literature. Progress on novel methods for upscaling pore-scale models for reactive transport are discussed, and are being applied to mineral precipitation patterns observed in natural analogues. H.Y. and T. D. were supported as part of the Center for Frontiers of Subsurface Energy Security, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Award Number DE-SC0001114. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

Fluid displacement fronts in porous media: pore scale interfacial jumps, pressure bursts and acoustic emissions

NASA Astrophysics Data System (ADS)

Moebius, Franziska; Or, Dani

2014-05-01

The macroscopically smooth and regular motion of fluid fronts in porous media is composed of numerous rapid pore-scale interfacial jumps and pressure bursts that involve intense interfacial energy release in the form of acoustic emissions. The characteristics of these pore scale events affect residual phase entrapment and transport properties behind the front. We present experimental studies using acoustic emission technique (AE), rapid imaging, and liquid pressure measurements to characterize these processes during drainage and imbibition in simple porous media. Imbibition and drainage produce different AE signatures (AE amplitudes obey a power law). For rapid drainage, AE signals persist long after cessation of front motion reflecting fluid redistribution and interfacial relaxation. Imaging revealed that the velocity of interfacial jumps often exceeds front velocity by more than 50 fold and is highly inertial component (Re>1000). Pore invasion volumes reduced deduced from pressure fluctuations waiting times (for constant withdrawal rates) show remarkable agreement with geometrically-deduced pore volumes. Discrepancies between invaded volumes and geometrical pores increase with increasing capillary numbers due to constraints on evacuation opportunity times and simultaneous invasion events. A mechanistic model for interfacial motions in a pore-throat network was developed to investigate interfacial dynamics focusing on the role of inertia. Results suggest that while pore scale dynamics were sensitive to variations in pore geometry and boundary conditions, inertia exerted only a minor effect on phase entrapment. The study on pore scale invasion events paints a complex picture of rapid and inertial motions and provides new insights on mechanisms at displacement fronts that are essential for improved macroscopic description of multiphase flows in porous media.

Parametric study of thin film evaporation from nanoporous membranes

NASA Astrophysics Data System (ADS)

Wilke, Kyle L.; Barabadi, Banafsheh; Lu, Zhengmao; Zhang, TieJun; Wang, Evelyn N.

2017-10-01

The performance and lifetime of advanced electronics are often dictated by the ability to dissipate heat generated within the device. Thin film evaporation from nanoporous membranes is a promising thermal management approach, which reduces the thermal transport distance across the liquid film while also providing passive capillary pumping of liquid to the evaporating interface. In this work, we investigated the dependence of thin film evaporation from nanoporous membranes on a variety of geometric parameters. Anodic aluminum oxide membranes were used as experimental templates, where pore radii of 28-75 nm, porosities of 0.1-0.35, and meniscus locations down to 1 μm within the pore were tested. We demonstrated different heat transfer regimes and observed more than an order of magnitude increase in dissipated heat flux by operating in the pore-level evaporation regime. The pore diameter had little effect on pore-level evaporation performance due to the negligible conduction resistance from the pore wall to the evaporating interface. The dissipated heat flux scaled with porosity as the evaporative area increased. Furthermore, moving the meniscus as little as 1 μm into the pore decreased the dissipated heat flux by more than a factor of two due to the added resistance to vapor escaping the pore. The experimental results elucidate thin film evaporation from nanopores and confirm findings of recent modeling efforts. This work also provides guidance for the design of future thin film evaporation devices for advanced thermal management. Furthermore, evaporation from nanopores is relevant to water purification, chemical separations, microfluidics, and natural processes such as transpiration.

VISUALIZATION AND SIMULATION OF NON-AQUEOUS PHASE LIQUIDS SOLUBILIZATION IN PORE NETWORKS

EPA Science Inventory

The design of in-situ remediation of contaminated soils is mostly based on a description at the macroscopic scale using a averaged quantities. These cannot address issues at the pore and pore network scales. In this paper, visualization experiments and numerical simulations in ...

Improved transfer of graphene for gated Schottky-junction, vertical, organic, field-effect transistors.

PubMed

Lemaitre, Maxime G; Donoghue, Evan P; McCarthy, Mitchell A; Liu, Bo; Tongay, Sefaattin; Gila, Brent; Kumar, Purushottam; Singh, Rajiv K; Appleton, Bill R; Rinzler, Andrew G

2012-10-23

An improved process for graphene transfer was used to demonstrate high performance graphene enabled vertical organic field effect transistors (G-VFETs). The process reduces disorder and eliminates the polymeric residue that typically plagues transferred films. The method also allows for purposely creating pores in the graphene of a controlled areal density. Transconductance observed in G-VFETs fabricated with a continuous (pore-free) graphene source electrode is attributed to modulation of the contact barrier height between the graphene and organic semiconductor due to a gate field induced Fermi level shift in the low density of electronic-states graphene electrode. Pores introduced in the graphene source electrode are shown to boost the G-VFET performance, which scales with the areal pore density taking advantage of both barrier height lowering and tunnel barrier thinning. Devices with areal pore densities of 20% exhibit on/off ratios and output current densities exceeding 10(6) and 200 mA/cm(2), respectively, at drain voltages below 5 V.

Understanding fluid transport through the multiscale pore network of a natural shale

NASA Astrophysics Data System (ADS)

Davy, Catherine A.; Nguyen Kim, Thang; Song, Yang; Troadec, David; Blanchenet, Anne-Marie; Adler, Pierre M.

2017-06-01

The pore structure of a natural shale is obtained by three imaging means. Micro-tomography results are extended to provide the spatial arrangement of the minerals and pores present at a voxel size of 700 nm (the macroscopic scale). FIB/SEM provides a 3D representation of the porous clay matrix on the so-called mesoscopic scale (10-20 nm); a connected pore network, devoid of cracks, is obtained for two samples out of five, while the pore network is connected through cracks for two other samples out of five. Transmission Electron Microscopy (TEM) is used to visualize the pore space with a typical pixel size of less than 1 nm and a porosity ranging from 0.12 to 0.25. On this scale, in the absence of 3D images, the pore structure is reconstructed by using a classical technique, which is based on truncated Gaussian fields. Permeability calculations are performed with the Lattice Boltzmann Method on the nanoscale, on the mesoscale, and on the combination of the two. Upscaling is finally done (by a finite volume approach) on the bigger macroscopic scale. Calculations show that, in the absence of cracks, the contribution of the nanoscale pore structure on the overall permeability is similar to that of the mesoscale. Complementarily, the macroscopic permeability is measured on a centimetric sample with a neutral fluid (ethanol). The upscaled permeability on the macroscopic scale is in good agreement with the experimental results.

Quantifying uncertainty and computational complexity for pore-scale simulations

NASA Astrophysics Data System (ADS)

Chen, C.; Yuan, Z.; Wang, P.; Yang, X.; Zhenyan, L.

2016-12-01

Pore-scale simulation is an essential tool to understand the complex physical process in many environmental problems, from multi-phase flow in the subsurface to fuel cells. However, in practice, factors such as sample heterogeneity, data sparsity and in general, our insufficient knowledge of the underlying process, render many simulation parameters and hence the prediction results uncertain. Meanwhile, most pore-scale simulations (in particular, direct numerical simulation) incur high computational cost due to finely-resolved spatio-temporal scales, which further limits our data/samples collection. To address those challenges, we propose a novel framework based on the general polynomial chaos (gPC) and build a surrogate model representing the essential features of the underlying system. To be specific, we apply the novel framework to analyze the uncertainties of the system behavior based on a series of pore-scale numerical experiments, such as flow and reactive transport in 2D heterogeneous porous media and 3D packed beds. Comparing with recent pore-scale uncertainty quantification studies using Monte Carlo techniques, our new framework requires fewer number of realizations and hence considerably reduce the overall computational cost, while maintaining the desired accuracy.

Upscaling of reaction rates in reactive transport using pore-scale reactive transport model

NASA Astrophysics Data System (ADS)

Yoon, H.; Dewers, T. A.; Arnold, B. W.; Major, J. R.; Eichhubl, P.; Srinivasan, S.

2013-12-01

Dissolved CO2 during geological CO2 storage may react with minerals in fractured rocks, confined aquifers, or faults, resulting in mineral precipitation and dissolution. The overall rate of reaction can be affected by coupled processes among hydrodynamics, transport, and reactions at the (sub) pore-scale. In this research pore-scale modeling of coupled fluid flow, reactive transport, and heterogeneous reaction at the mineral surface is applied to account for permeability alterations caused by precipitation-induced pore-blocking. This work is motivated by the observed CO2 seeps from a natural analog to geologic CO2 sequestration at Crystal Geyser, Utah. A key observation is the lateral migration of CO2 seep sites at a scale of ~ 100 meters over time. A pore-scale model provides fundamental mechanistic explanations of how calcite precipitation alters flow paths by pore plugging under different geochemical compositions and pore configurations. In addition, response function of reaction rates will be constructed from pore-scale simulations which account for a range of reaction regimes characterized by the Damkohler and Peclet numbers. Newly developed response functions will be used in a continuum scale model that may account for large-scale phenomena mimicking lateral migration of surface CO2 seeps. Comparison of field observations and simulations results will provide mechanistic explanations of the lateral migration and enhance our understanding of subsurface processes associated with the CO2 injection. This work is supported as part of the Center for Frontiers of Subsurface Energy Security, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Award Number DE-SC0001114. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

Measurements of pore-scale flow through apertures

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chojnicki, Kirsten

Pore-scale aperture effects on flow in pore networks was studied in the laboratory to provide a parameterization for use in transport models. Four cases were considered: regular and irregular pillar/pore alignment with and without an aperture. The velocity field of each case was measured and simulated, providing quantitatively comparable results. Two aperture effect parameterizations were considered: permeability and transmission. Permeability values varied by an order of magnitude between the cases with and without apertures. However, transmission did not correlate with permeability. Despite having much greater permeability the regular aperture case permitted less transmission than the regular case. Moreover, both irregularmore » cases had greater transmission than the regular cases, a difference not supported by the permeabilities. Overall, these findings suggest that pore-scale aperture effects on flow though a pore-network may not be adequately captured by properties such as permeability for applications that are interested in determining particle transport volume and timing.« less

Reply to 'Comments on upscaling geochemical reaction rates usingpore-scale network modeling' by Peter C. Lichtner and Qinjun Kang

DOE Office of Scientific and Technical Information (OSTI.GOV)

Li, Li; Peters, Catherine A.; Celia, Michael A.

2006-05-03

Our paper "Upscaling geochemical reaction rates usingpore-scale network modeling" presents a novel application of pore-scalenetwork modeling to upscale mineral dissolution and precipitationreaction rates from the pore scale to the continuum scale, anddemonstrates the methodology by analyzing the scaling behavior ofanorthite and kaolinite reaction kinetics under conditions related to CO2sequestration. We conclude that under highly acidic conditions relevantto CO2 sequestration, the traditional continuum-based methodology may notcapture the spatial variation in concentrations from pore to pore, andscaling tools may be important in correctly modeling reactive transportprocesses in such systems. This work addresses the important butdifficult question of scaling mineral dissolution and precipitationreactionmore » kinetics, which is often ignored in fields such as geochemistry,water resources, and contaminant hydrology. Although scaling of physicalprocesses has been studied for almost three decades, very few studieshave examined the scaling issues related to chemical processes, despitetheir importance in governing the transport and fate of contaminants insubsurface systems.« less

Characterization and Quantification of the Pore Structures of the Shale Oil Reservoir Formations in Multiscale

NASA Astrophysics Data System (ADS)

Liu, K.; Ostadhassan, M.

2016-12-01

Due to the fast development of hydraulic fracturing and horizontal drilling, shale formations now are one important resource of energy in North America. Characterizing the pore structure of these shale formations is of critical importance in understanding the original oil/gas in place and also the flow properties of the rock matrix. Pore with different properties such as pore size and pore shape can impact the physical, mechanical and chemical properties including strength, elastic modulus, permeability and conductivity. Nowadays, image analysis has been a robust method to quantify the pore information from the porous medium.SEM has been one of the most useful tools to study the pore microstructures due to its high depth of focus which can provide detailed topographical information about the surface. The suitable difference between solid matrix and pores due to the different gray level pixels can be used to study the pore structures.In this paper, we characterized and quantified the pore structures of rock samples from Middle Bakken Formation which is a typical unconventional reservoir in North America. High resolution SEM images of five samples we chose based on the gamma logs were derived after sample preparation. After determining the threshold of the images, we extracted the pore spaces. Then we analyzed the pore structures properties such as pore size distributions and pore shape distributions of the five samples and compared based on their mineral compositions. After that, we analyzed their heterogeneity and isotropy properties which have been identified as an important factor affecting reservoir productivity. Finally, we studied the impact of scale effect on the pore structures characterization.

Developing an Effective Model for Shale Gas Flow in Nano-scale Pore Clusters based on FIB-SEM Images

NASA Astrophysics Data System (ADS)

Jiang, W. B.; Lin, M.; Yi, Z. X.; Li, H. S.

2016-12-01

Nano-scale pores existed in the form of clusters are the controlling void space in shale gas reservoir. Gas transport in nanopores which has a significant influence on shale gas' recoverability displays multiple transport regimes, including viscous, slippage flow and Knudsen diffusion. In addition, it is also influenced by pore space characteristics. For convenience and efficiency consideration, it is necessary to develop an upscaling model from nano pore to pore cluster scale. Existing models are more like framework functions that provide a format, because the parameters that represent pore space characteristics are underdetermined and may have multiple possibilities. Therefore, it is urgent to make them clear and obtained a model that is closer to reality. FIB-SEM imaging technology is able to acquire three dimensional images with nanometer resolution that nano pores can be visible. Based on the images of two shale samples, we used a high-precision pore network extraction algorithm to generate equivalent pore networks and simulate multiple regime (non-Darcy) flow in it. Several structural parameters can be obtained through pore network modelling. It is found that although the throat-radius distributions are very close, throat flux-radius distributions of different samples can be divided into two categories. The variation of tortuosity with pressure and the overall trend of throat-flux distribution changes with pressure are disclosed. A deeper understanding of shale gas flow in nano-scale pore clusters is obtained. After all, an upscaling model that connects absolute permeability, apparent permeability and other characteristic parameters is proposed, and the best parameter scheme considering throat number-radius distribution and flowing porosity for this model is selected out of three schemes based on pore scale results, and it can avoid multiple-solution problem and is useful in reservoir modelling and experiment result analysis, etc. This work is supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB10020302), the National Natural Science Foundation of China (41574129), the Key Instrument Developing Project of the CAS (ZDYZ2012-1-08-02), the 973 Program (2014CB239004).

Impact of spatially correlated pore-scale heterogeneity on drying porous media

NASA Astrophysics Data System (ADS)

Borgman, Oshri; Fantinel, Paolo; Lühder, Wieland; Goehring, Lucas; Holtzman, Ran

2017-07-01

We study the effect of spatially-correlated heterogeneity on isothermal drying of porous media. We combine a minimal pore-scale model with microfluidic experiments with the same pore geometry. Our simulated drying behavior compares favorably with experiments, considering the large sensitivity of the emergent behavior to the uncertainty associated with even small manufacturing errors. We show that increasing the correlation length in particle sizes promotes preferential drying of clusters of large pores, prolonging liquid connectivity and surface wetness and thus higher drying rates for longer periods. Our findings improve our quantitative understanding of how pore-scale heterogeneity impacts drying, which plays a role in a wide range of processes ranging from fuel cells to curing of paints and cements to global budgets of energy, water and solutes in soils.

Pore-scale simulation of microbial growth using a genome-scale metabolic model: Implications for Darcy-scale reactive transport

NASA Astrophysics Data System (ADS)

Tartakovsky, G. D.; Tartakovsky, A. M.; Scheibe, T. D.; Fang, Y.; Mahadevan, R.; Lovley, D. R.

2013-09-01

Recent advances in microbiology have enabled the quantitative simulation of microbial metabolism and growth based on genome-scale characterization of metabolic pathways and fluxes. We have incorporated a genome-scale metabolic model of the iron-reducing bacteria Geobacter sulfurreducens into a pore-scale simulation of microbial growth based on coupling of iron reduction to oxidation of a soluble electron donor (acetate). In our model, fluid flow and solute transport is governed by a combination of the Navier-Stokes and advection-diffusion-reaction equations. Microbial growth occurs only on the surface of soil grains where solid-phase mineral iron oxides are available. Mass fluxes of chemical species associated with microbial growth are described by the genome-scale microbial model, implemented using a constraint-based metabolic model, and provide the Robin-type boundary condition for the advection-diffusion equation at soil grain surfaces. Conventional models of microbially-mediated subsurface reactions use a lumped reaction model that does not consider individual microbial reaction pathways, and describe reactions rates using empirically-derived rate formulations such as the Monod-type kinetics. We have used our pore-scale model to explore the relationship between genome-scale metabolic models and Monod-type formulations, and to assess the manifestation of pore-scale variability (microenvironments) in terms of apparent Darcy-scale microbial reaction rates. The genome-scale model predicted lower biomass yield, and different stoichiometry for iron consumption, in comparison to prior Monod formulations based on energetics considerations. We were able to fit an equivalent Monod model, by modifying the reaction stoichiometry and biomass yield coefficient, that could effectively match results of the genome-scale simulation of microbial behaviors under excess nutrient conditions, but predictions of the fitted Monod model deviated from those of the genome-scale model under conditions in which one or more nutrients were limiting. The fitted Monod kinetic model was also applied at the Darcy scale; that is, to simulate average reaction processes at the scale of the entire pore-scale model domain. As we expected, even under excess nutrient conditions for which the Monod and genome-scale models predicted equal reaction rates at the pore scale, the Monod model over-predicted the rates of biomass growth and iron and acetate utilization when applied at the Darcy scale. This discrepancy is caused by an inherent assumption of perfect mixing over the Darcy-scale domain, which is clearly violated in the pore-scale models. These results help to explain the need to modify the flux constraint parameters in order to match observations in previous applications of the genome-scale model at larger scales. These results also motivate further investigation of quantitative multi-scale relationships between fundamental behavior at the pore scale (where genome-scale models are appropriately applied) and observed behavior at larger scales (where predictions of reactive transport phenomena are needed).

Pore-scale simulation of microbial growth using a genome-scale metabolic model: Implications for Darcy-scale reactive transport

NASA Astrophysics Data System (ADS)

Scheibe, T. D.; Tartakovsky, G.; Tartakovsky, A. M.; Fang, Y.; Mahadevan, R.; Lovley, D. R.

2012-12-01

Recent advances in microbiology have enabled the quantitative simulation of microbial metabolism and growth based on genome-scale characterization of metabolic pathways and fluxes. We have incorporated a genome-scale metabolic model of the iron-reducing bacteria Geobacter sulfurreducens into a pore-scale simulation of microbial growth based on coupling of iron reduction to oxidation of a soluble electron donor (acetate). In our model, fluid flow and solute transport is governed by a combination of the Navier-Stokes and advection-diffusion-reaction equations. Microbial growth occurs only on the surface of soil grains where solid-phase mineral iron oxides are available. Mass fluxes of chemical species associated with microbial growth are described by the genome-scale microbial model, implemented using a constraint-based metabolic model, and provide the Robin-type boundary condition for the advection-diffusion equation at soil grain surfaces. Conventional models of microbially-mediated subsurface reactions use a lumped reaction model that does not consider individual microbial reaction pathways, and describe reactions rates using empirically-derived rate formulations such as the Monod-type kinetics. We have used our pore-scale model to explore the relationship between genome-scale metabolic models and Monod-type formulations, and to assess the manifestation of pore-scale variability (microenvironments) in terms of apparent Darcy-scale microbial reaction rates. The genome-scale model predicted lower biomass yield, and different stoichiometry for iron consumption, in comparison to prior Monod formulations based on energetics considerations. We were able to fit an equivalent Monod model, by modifying the reaction stoichiometry and biomass yield coefficient, that could effectively match results of the genome-scale simulation of microbial behaviors under excess nutrient conditions, but predictions of the fitted Monod model deviated from those of the genome-scale model under conditions in which one or more nutrients were limiting. The fitted Monod kinetic model was also applied at the Darcy scale; that is, to simulate average reaction processes at the scale of the entire pore-scale model domain. As we expected, even under excess nutrient conditions for which the Monod and genome-scale models predicted equal reaction rates at the pore scale, the Monod model over-predicted the rates of biomass growth and iron and acetate utilization when applied at the Darcy scale. This discrepancy is caused by an inherent assumption of perfect mixing over the Darcy-scale domain, which is clearly violated in the pore-scale models. These results help to explain the need to modify the flux constraint parameters in order to match observations in previous applications of the genome-scale model at larger scales. These results also motivate further investigation of quantitative multi-scale relationships between fundamental behavior at the pore scale (where genome-scale models are appropriately applied) and observed behavior at larger scales (where predictions of reactive transport phenomena are needed).

Pore-scale simulation of microbial growth using a genome-scale metabolic model: Implications for Darcy-scale reactive transport

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tartakovsky, Guzel D.; Tartakovsky, Alexandre M.; Scheibe, Timothy D.

2013-09-07

Recent advances in microbiology have enabled the quantitative simulation of microbial metabolism and growth based on genome-scale characterization of metabolic pathways and fluxes. We have incorporated a genome-scale metabolic model of the iron-reducing bacteria Geobacter sulfurreducens into a pore-scale simulation of microbial growth based on coupling of iron reduction to oxidation of a soluble electron donor (acetate). In our model, fluid flow and solute transport is governed by a combination of the Navier-Stokes and advection-diffusion-reaction equations. Microbial growth occurs only on the surface of soil grains where solid-phase mineral iron oxides are available. Mass fluxes of chemical species associated withmore » microbial growth are described by the genome-scale microbial model, implemented using a constraint-based metabolic model, and provide the Robin-type boundary condition for the advection-diffusion equation at soil grain surfaces. Conventional models of microbially-mediated subsurface reactions use a lumped reaction model that does not consider individual microbial reaction pathways, and describe reactions rates using empirically-derived rate formulations such as the Monod-type kinetics. We have used our pore-scale model to explore the relationship between genome-scale metabolic models and Monod-type formulations, and to assess the manifestation of pore-scale variability (microenvironments) in terms of apparent Darcy-scale microbial reaction rates. The genome-scale model predicted lower biomass yield, and different stoichiometry for iron consumption, in comparisonto prior Monod formulations based on energetics considerations. We were able to fit an equivalent Monod model, by modifying the reaction stoichiometry and biomass yield coefficient, that could effectively match results of the genome-scale simulation of microbial behaviors under excess nutrient conditions, but predictions of the fitted Monod model deviated from those of the genome-scale model under conditions in which one or more nutrients were limiting. The fitted Monod kinetic model was also applied at the Darcy scale; that is, to simulate average reaction processes at the scale of the entire pore-scale model domain. As we expected, even under excess nutrient conditions for which the Monod and genome-scale models predicted equal reaction rates at the pore scale, the Monod model over-predicted the rates of biomass growth and iron and acetate utilization when applied at the Darcy scale. This discrepancy is caused by an inherent assumption of perfect mixing over the Darcy-scale domain, which is clearly violated in the pore-scale models. These results help to explain the need to modify the flux constraint parameters in order to match observations in previous applications of the genome-scale model at larger scales. These results also motivate further investigation of quantitative multi-scale relationships between fundamental behavior at the pore scale (where genome-scale models are appropriately applied) and observed behavior at larger scales (where predictions of reactive transport phenomena are needed).« less

Pore-scale and continuum simulations of solute transport micromodel benchmark experiments

DOE PAGES

Oostrom, M.; Mehmani, Y.; Romero-Gomez, P.; ...

2014-06-18

Four sets of nonreactive solute transport experiments were conducted with micromodels. Three experiments with one variable, i.e., flow velocity, grain diameter, pore-aspect ratio, and flow-focusing heterogeneity were in each set. The data sets were offered to pore-scale modeling groups to test their numerical simulators. Each set consisted of two learning experiments, for which our results were made available, and one challenge experiment, for which only the experimental description and base input parameters were provided. The experimental results showed a nonlinear dependence of the transverse dispersion coefficient on the Peclet number, a negligible effect of the pore-aspect ratio on transverse mixing,more » and considerably enhanced mixing due to flow focusing. Five pore-scale models and one continuum-scale model were used to simulate the experiments. Of the pore-scale models, two used a pore-network (PN) method, two others are based on a lattice Boltzmann (LB) approach, and one used a computational fluid dynamics (CFD) technique. Furthermore, we used the learning experiments, by the PN models, to modify the standard perfect mixing approach in pore bodies into approaches to simulate the observed incomplete mixing. The LB and CFD models used the learning experiments to appropriately discretize the spatial grid representations. For the continuum modeling, the required dispersivity input values were estimated based on published nonlinear relations between transverse dispersion coefficients and Peclet number. Comparisons between experimental and numerical results for the four challenge experiments show that all pore-scale models were all able to satisfactorily simulate the experiments. The continuum model underestimated the required dispersivity values, resulting in reduced dispersion. The PN models were able to complete the simulations in a few minutes, whereas the direct models, which account for the micromodel geometry and underlying flow and transport physics, needed up to several days on supercomputers to resolve the more complex problems.« less

Effects of pore-scale physics on uranium geochemistry in Hanford sediments

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hu, Qinhong; Ewing, Robert P.

Overall, this work examines a key scientific issue, mass transfer limitations at the pore-scale, using both new instruments with high spatial resolution, and new conceptual and modeling paradigms. The complementary laboratory and numerical approaches connect pore-scale physics to macroscopic measurements, providing a previously elusive scale integration. This Exploratory research project produced five peer-reviewed journal publications and eleven scientific presentations. This work provides new scientific understanding, allowing the DOE to better incorporate coupled physical and chemical processes into decision making for environmental remediation and long-term stewardship.

Pore characteristics and their emergent effect on water adsorption and transport in clays using small-angle neutron scattering with contrast variation

NASA Astrophysics Data System (ADS)

Ding, M.; Hartl, M.; Wang, Y.; Hjelm, R.

2013-12-01

In nuclear waste management, clays are canonical materials in the construction of engineered barriers. They are also naturally occurring reactive minerals which play an important role in retention and colloidal facilitated reactive transport in subsurface systems. Knowledge of total and accessible porosity in clays is crucial in determining fluids transport behavior in clays. It will provide fundamental insight on the performance efficiency of specific clays as a barrier material and their role in regulating radionuclide transport in subsurface environments. The aim of the present work is to experimentally investigate the change in pore characteristics of clays as function of moisture content, and to determine their pore character in relation to their water retention capacity. Recent developments in small-angle neutron scattering (SANS) techniques allow quantitative measurement of pore morphology and size distribution of various materials in their pristine state under various sample environments (exposure to solution, high temperature, and so on). Furthermore, due to dramatic different neutron scattering properties of hydrogen and deuterium, one can readily use contrast variation, which is the isotopic labeling with various ratios of H and D (e.g. mixture of H2O/D2O) to highlight or suppress features of the sample. This is particularly useful in the study of complex pore system such as clays. In this study, we have characterized the pore structures for a number of clays including clay minerals and field samples which are relevant to high-level waste systems under various sample environments (e.g., humidity, temperature and pressure) using SANS. Our results suggest that different clays show unique pore features under various sample environments. To distinguish between accessible/non-accessible pores and the nature of pore filling (e.g. the quantity of H2O adsorbed by clays, and the distribution of H2O in relation to pore character) to water, clays were exposed for various periods to a specific humidity (e.g., relative humidity: RH=100%, RH=75%). The humidity is controlled by using saturated aqueous solutions, consisting of specific H2O/D2O mixtures. Our results have shown distinct variations in water adsorption and moisture diffusivity among clays. Our results allow us to obtain on the pore scale porosity changes due to water movement in clays. As emergent transport property, nano- to micro-scale structural characterization is crucial in providing insights into pore-scale transport processes, which are pertinent to upscale continuum model development involving flow and transport at low water content, flow and phase behavior under confinement, and low-permeability media.

Passive advection-dispersion in networks of pipes: Effect of connectivity and relationship to permeability

NASA Astrophysics Data System (ADS)

Bernabé, Y.; Wang, Y.; Qi, T.; Li, M.

2016-02-01

The main purpose of this work is to investigate the relationship between passive advection-dispersion and permeability in porous materials presumed to be statistically homogeneous at scales larger than the pore scale but smaller than the reservoir scale. We simulated fluid flow through pipe network realizations with different pipe radius distributions and different levels of connectivity. The flow simulations used periodic boundary conditions, allowing monitoring of the advective motion of solute particles in a large periodic array of identical network realizations. In order to simulate dispersion, we assumed that the solute particles obeyed Taylor dispersion in individual pipes. When a particle entered a pipe, a residence time consistent with local Taylor dispersion was randomly assigned to it. When exiting the pipe, the particle randomly proceeded into one of the pipes connected to the original one according to probabilities proportional to the outgoing volumetric flow in each pipe. For each simulation we tracked the motion of at least 6000 solute particles. The mean fluid velocity was 10-3 ms-1, and the distance traveled was on the order of 10 m. Macroscopic dispersion was quantified using the method of moments. Despite differences arising from using different types of lattices (simple cubic, body-centered cubic, and face-centered cubic), a number of general observations were made. Longitudinal dispersion was at least 1 order of magnitude greater than transverse dispersion, and both strongly increased with decreasing pore connectivity and/or pore size variability. In conditions of variable hydraulic radius and fixed pore connectivity and pore size variability, the simulated dispersivities increased as power laws of the hydraulic radius and, consequently, of permeability, in agreement with previously published experimental results. Based on these observations, we were able to resolve some of the complexity of the relationship between dispersivity and permeability.

Pore-Scale Modeling of Pore Structure Effects on P-Wave Scattering Attenuation in Dry Rocks

PubMed Central

Li, Tianyang; Qiu, Hao; Wang, Feifei

2015-01-01

Underground rocks usually have complex pore system with a variety of pore types and a wide range of pore size. The effects of pore structure on elastic wave attenuation cannot be neglected. We investigated the pore structure effects on P-wave scattering attenuation in dry rocks by pore-scale modeling based on the wave theory and the similarity principle. Our modeling results indicate that pore size, pore shape (such as aspect ratio), and pore density are important factors influencing P-wave scattering attenuation in porous rocks, and can explain the variation of scattering attenuation at the same porosity. From the perspective of scattering attenuation, porous rocks can safely suit to the long wavelength assumption when the ratio of wavelength to pore size is larger than 15. Under the long wavelength condition, the scattering attenuation coefficient increases as a power function as the pore density increases, and it increases exponentially with the increase in aspect ratio. For a certain porosity, rocks with smaller aspect ratio and/or larger pore size have stronger scattering attenuation. When the pore aspect ratio is larger than 0.5, the variation of scattering attenuation at the same porosity is dominantly caused by pore size and almost independent of the pore aspect ratio. These results lay a foundation for pore structure inversion from elastic wave responses in porous rocks. PMID:25961729

Microfluidic Experiments Studying Pore Scale Interactions of Microbes and Geochemistry

NASA Astrophysics Data System (ADS)

Chen, M.; Kocar, B. D.

2016-12-01

Understanding how physical phenomena, chemical reactions, and microbial behavior interact at the pore-scale is crucial to understanding larger scale trends in groundwater chemistry. Recent studies illustrate the utility of microfluidic devices for illuminating pore-scale physical-biogeochemical processes and their control(s) on the cycling of iron, uranium, and other important elements 1-3. These experimental systems are ideal for examining geochemical reactions mediated by microbes, which include processes governed by complex biological phenomenon (e.g. biofilm formation, etc.)4. We present results of microfluidic experiments using a model metal reducing bacteria and varying pore geometries, exploring the limitations of the microorganisms' ability to access tight pore spaces, and examining coupled biogeochemical-physical controls on the cycling of redox sensitive metals. Experimental results will provide an enhanced understanding of coupled physical-biogeochemical processes transpiring at the pore-scale, and will constrain and compliment continuum models used to predict and describe the subsurface cycling of redox-sensitive elements5. 1. Vrionis, H. A. et al. Microbiological and geochemical heterogeneity in an in situ uranium bioremediation field site. Appl. Environ. Microbiol. 71, 6308-6318 (2005). 2. Pearce, C. I. et al. Pore-scale characterization of biogeochemical controls on iron and uranium speciation under flow conditions. Environ. Sci. Technol. 46, 7992-8000 (2012). 3. Zhang, C., Liu, C. & Shi, Z. Micromodel investigation of transport effect on the kinetics of reductive dissolution of hematite. Environ. Sci. Technol. 47, 4131-4139 (2013). 4. Ginn, T. R. et al. Processes in microbial transport in the natural subsurface. Adv. Water Resour. 25, 1017-1042 (2002). 5. Scheibe, T. D. et al. Coupling a genome-scale metabolic model with a reactive transport model to describe in situ uranium bioremediation. Microb. Biotechnol. 2, 274-286 (2009).

Dispersion upscaling from a pore scale characterization of Lagrangian velocities

NASA Astrophysics Data System (ADS)

Turuban, Régis; de Anna, Pietro; Jiménez-Martínez, Joaquín; Tabuteau, Hervé; Méheust, Yves; Le Borgne, Tanguy

2013-04-01

Mixing and reactive transport are primarily controlled by the interplay between diffusion, advection and reaction at pore scale. Yet, how the distribution and spatial correlation of the velocity field at pore scale impact these processes is still an open question. Here we present an experimental investigation of the distribution and correlation of pore scale velocities and its relation with upscaled dispersion. We use a quasi two-dimensional (2D) horizontal set up, consisting of two glass plates filled with cylinders representing the grains of the porous medium : the cell is built by soft lithography technique, wich allows for full control of the system geometry. The local velocity field is quantified from particle tracking velocimetry using microspheres that are advected with the pore scale flow. Their displacement is purely advective, as the particle size is chosen large enough to avoid diffusion. We thus obtain particle trajectories as well as lagrangian velocities in the entire system. The measured velocity field shows the existence of a network of preferential flow paths in channels with high velocities, as well as very low velocity in stagnation zones, with a non Gaussian distribution. Lagrangian velocities are long range correlated in time, which implies a non-fickian scaling of the longitudinal variance of particle positions. To upscale this process we develop an effective transport model, based on correlated continous time random walk, which is entirely parametrized by the pore scale velocity distribution and correlation. The model predictions are compared with conservative tracer test data for different Peclet numbers. Furthermore, we investigate the impact of different pore geometries on the distribution and correlation of Lagrangian velocities and we discuss the link between these properties and the effective dispersion behavior.

Application Of Positron Beams For The Characterization Of Nano-scale Pores In Thin Films

NASA Astrophysics Data System (ADS)

Hirata, K.; Ito, K.; Kobayashi, Y.; Suzuki, R.; Ohdaira, T.; Eijt, S. W. H.; Schut, H.; van Veen, A.

2003-08-01

We applied three positron annihilation techniques, positron 3γ-annihilation spectroscopy, positron annihilation lifetime spectroscopy, and angular correlation of annihilation radiation, to the characterization of nano-scale pores in thin films by combining them with variable-energy positron beams. Characterization of pores in thin films is an important part of the research on various thin films of industrial importance. The results of our recent studies on pore characterization of thin films by positron beams will be reported here.

Evaluating the hydraulic and transport properties of peat soil using pore network modeling and X-ray micro computed tomography

NASA Astrophysics Data System (ADS)

Gharedaghloo, Behrad; Price, Jonathan S.; Rezanezhad, Fereidoun; Quinton, William L.

2018-06-01

Micro-scale properties of peat pore space and their influence on hydraulic and transport properties of peat soils have been given little attention so far. Characterizing the variation of these properties in a peat profile can increase our knowledge on the processes controlling contaminant transport through peatlands. As opposed to the common macro-scale (or bulk) representation of groundwater flow and transport processes, a pore network model (PNM) simulates flow and transport processes within individual pores. Here, a pore network modeling code capable of simulating advective and diffusive transport processes through a 3D unstructured pore network was developed; its predictive performance was evaluated by comparing its results to empirical values and to the results of computational fluid dynamics (CFD) simulations. This is the first time that peat pore networks have been extracted from X-ray micro-computed tomography (μCT) images of peat deposits and peat pore characteristics evaluated in a 3D approach. Water flow and solute transport were modeled in the unstructured pore networks mapped directly from μCT images. The modeling results were processed to determine the bulk properties of peat deposits. Results portray the commonly observed decrease in hydraulic conductivity with depth, which was attributed to the reduction of pore radius and increase in pore tortuosity. The increase in pore tortuosity with depth was associated with more decomposed peat soil and decreasing pore coordination number with depth, which extended the flow path of fluid particles. Results also revealed that hydraulic conductivity is isotropic locally, but becomes anisotropic after upscaling to core-scale; this suggests the anisotropy of peat hydraulic conductivity observed in core-scale and field-scale is due to the strong heterogeneity in the vertical dimension that is imposed by the layered structure of peat soils. Transport simulations revealed that for a given solute, the effective diffusion coefficient decreases with depth due to the corresponding increase of diffusional tortuosity. Longitudinal dispersivity of peat also was computed by analyzing advective-dominant transport simulations that showed peat dispersivity is similar to the empirical values reported in the same peat soil; it is not sensitive to soil depth and does not vary much along the soil profile.

Flow regimes during immiscible displacement

DOE Office of Scientific and Technical Information (OSTI.GOV)

Armstrong, Ryan T.; Mcclure, James; Berrill, Mark A.

Fractional ow of immiscible phases occurs at the pore scale where grain surfaces and phases interfaces obstruct phase mobility. However, the larger scale behavior is described by a saturation-dependent phenomenological relationship called relative permeability. As a consequence, pore-scale parameters, such as phase topology and/ or geometry, and details of the flow regime cannot be directly related to Darcy-scale flow parameters. It is well understood that relative permeability is not a unique relationship of wetting-phase saturation and rather depends on the experimental conditions at which it is measured. Herein we use fast X-ray microcomputed tomography to image pore-scale phase arrangements duringmore » fractional flow and then forward simulate the flow regimes using the lattice-Boltzmann method to better understand the underlying pore-scale flow regimes and their influence on Darcy-scale parameters. We find that relative permeability is highly dependent on capillary number and that the Corey model fits the observed trends. At the pore scale, while phase topologies are continuously changing on the scale of individual pores, the Euler characteristic of the nonwetting phase (NWP) averaged over a sufficiently large field of view can describe the bulk topological characteristics; the Euler characteristic decreases with increasing capillary number resulting in an increase in relative permeability. Lastly, we quantify the fraction of NWP that flows through disconnected ganglion dynamics and demonstrate that this can be a significant fraction of the NWP flux for intermediate wetting-phase saturation. Furthermore, rate dependencies occur in our homogenous sample (without capillary end effect) and the underlying cause is attributed to ganglion flow that can significantly influence phase topology during the fractional flow of immiscible phases.« less

Flow regimes during immiscible displacement

DOE PAGES

Armstrong, Ryan T.; Mcclure, James; Berrill, Mark A.; ...

2017-02-01

Fractional ow of immiscible phases occurs at the pore scale where grain surfaces and phases interfaces obstruct phase mobility. However, the larger scale behavior is described by a saturation-dependent phenomenological relationship called relative permeability. As a consequence, pore-scale parameters, such as phase topology and/ or geometry, and details of the flow regime cannot be directly related to Darcy-scale flow parameters. It is well understood that relative permeability is not a unique relationship of wetting-phase saturation and rather depends on the experimental conditions at which it is measured. Herein we use fast X-ray microcomputed tomography to image pore-scale phase arrangements duringmore » fractional flow and then forward simulate the flow regimes using the lattice-Boltzmann method to better understand the underlying pore-scale flow regimes and their influence on Darcy-scale parameters. We find that relative permeability is highly dependent on capillary number and that the Corey model fits the observed trends. At the pore scale, while phase topologies are continuously changing on the scale of individual pores, the Euler characteristic of the nonwetting phase (NWP) averaged over a sufficiently large field of view can describe the bulk topological characteristics; the Euler characteristic decreases with increasing capillary number resulting in an increase in relative permeability. Lastly, we quantify the fraction of NWP that flows through disconnected ganglion dynamics and demonstrate that this can be a significant fraction of the NWP flux for intermediate wetting-phase saturation. Furthermore, rate dependencies occur in our homogenous sample (without capillary end effect) and the underlying cause is attributed to ganglion flow that can significantly influence phase topology during the fractional flow of immiscible phases.« less

Tracer counterpermeation analysis of diffusivity in finite-length nanopores with and without single-file dynamics

DOE PAGES

Ackerman, David M.; Evans, James W.

2017-01-19

Here, we perform a tracer counterpermeation (TCP) analysis for a stochastic model of diffusive transport through a narrow linear pore where passing of species within the pore is inhibited or even excluded (single-file diffusion). TCP involves differently labeled but otherwise identical particles from two decoupled infinite reservoirs adsorbing into opposite ends of the pore, and desorbing from either end. In addition to transient behavior, we assess steady-state concentration profiles, spatial correlations, particle number fluctuations, and diffusion fluxes through the pore. From the profiles and fluxes, we determine a generalized tracer diffusion coefficient D tr(x), at various positions x within themore » pore. D tr(x) has a plateau value in the pore center scaling inversely with the pore length, but it is enhanced near the pore openings. The latter feature reflects the effect of fluctuations in adsorption and desorption, and it is also associated with a nontrivial scaling of the concentration profiles near the pore openings.« less

Tracer counterpermeation analysis of diffusivity in finite-length nanopores with and without single-file dynamics

NASA Astrophysics Data System (ADS)

Ackerman, David M.; Evans, James W.

2017-01-01

We perform a tracer counterpermeation (TCP) analysis for a stochastic model of diffusive transport through a narrow linear pore where passing of species within the pore is inhibited or even excluded (single-file diffusion). TCP involves differently labeled but otherwise identical particles from two decoupled infinite reservoirs adsorbing into opposite ends of the pore, and desorbing from either end. In addition to transient behavior, we assess steady-state concentration profiles, spatial correlations, particle number fluctuations, and diffusion fluxes through the pore. From the profiles and fluxes, we determine a generalized tracer diffusion coefficient Dtr(x ) , at various positions x within the pore. Dtr(x ) has a plateau value in the pore center scaling inversely with the pore length, but it is enhanced near the pore openings. The latter feature reflects the effect of fluctuations in adsorption and desorption, and it is also associated with a nontrivial scaling of the concentration profiles near the pore openings.

An Image-based Micro-continuum Pore-scale Model for Gas Transport in Organic-rich Shale

NASA Astrophysics Data System (ADS)

Guo, B.; Tchelepi, H.

2017-12-01

Gas production from unconventional source rocks, such as ultra-tight shales, has increased significantly over the past decade. However, due to the extremely small pores ( 1-100 nm) and the strong material heterogeneity, gas flow in shale is still not well understood and poses challenges for predictive field-scale simulations. In recent years, digital rock analysis has been applied to understand shale gas transport at the pore-scale. An issue with rock images (e.g. FIB-SEM, nano-/micro-CT images) is the so-called "cutoff length", i.e., pores and heterogeneities below the resolution cannot be resolved, which leads to two length scales (resolved features and unresolved sub-resolution features) that are challenging for flow simulations. Here we develop a micro-continuum model, modified from the classic Darcy-Brinkman-Stokes framework, that can naturally couple the resolved pores and the unresolved nano-porous regions. In the resolved pores, gas flow is modeled with Stokes equation. In the unresolved regions where the pore sizes are below the image resolution, we develop an apparent permeability model considering non-Darcy flow at the nanoscale including slip flow, Knudsen diffusion, adsorption/desorption, surface diffusion, and real gas effect. The end result is a micro-continuum pore-scale model that can simulate gas transport in 3D reconstructed shale images. The model has been implemented in the open-source simulation platform OpenFOAM. In this paper, we present case studies to demonstrate the applicability of the model, where we use 3D segmented FIB-SEM and nano-CT shale images that include four material constituents: organic matter, clay, granular mineral, and pore. In addition to the pore structure and the distribution of the material constituents, we populate the model with experimental measurements (e.g. size distribution of the sub-resolution pores from nitrogen adsorption) and parameters from the literature and identify the relative importance of different physics on gas production. Overall, the micro-continuum model provides a novel tool for digital rock analysis of organic-rich shale.

Evidence for Enhanced Matrix Diffusion in Geological Environment

NASA Astrophysics Data System (ADS)

Sato, Kiminori; Fujimoto, Koichiro; Nakata, Masataka; Shikazono, Naotatsu

2013-01-01

Molecular diffusion in rock matrix, called as matrix diffusion, has been appreciated as a static process for elemental migration in geological environment that has been acknowledged in the context of geological disposal of radioactive waste. However, incomprehensible enhancement of matrix diffusion has been reported at a number of field test sites. Here, the matrix diffusion of saline water at Horonobe, Hokkaido, Japan is highlighted directly probing angstrom-scale pores on a field scale up to 1 km by positron--positronium annihilation spectroscopy. The first application of positron--positronium annihilation spectroscopy to field-scale geophysical research reveals the slight variation of angstrom-scale pores influenced by saline water diffusion with complete accuracy. We found widely interconnected 3 Å pores, which offer the pathway of saline water diffusion with the highly enhanced effective matrix diffusion coefficient of 4× 10-6 cm2 s-1. The present findings provide unambiguous evidence that the angstrom-scale pores enhance effective matrix diffusion on a field scale in geological environment.

Persistent Homology to describe Solid and Fluid Structures during Multiphase Flow

NASA Astrophysics Data System (ADS)

Herring, A. L.; Robins, V.; Liu, Z.; Armstrong, R. T.; Sheppard, A.

2017-12-01

The question of how to accurately and effectively characterize essential fluid and solid distributions and structures is a long-standing topic within the field of porous media and fluid transport. For multiphase flow applications, considerable research effort has been made to describe fluid distributions under a range of conditions; including quantification of saturation levels, fluid-fluid pressure differences and interfacial areas, and fluid connectivity. Recent research has effectively used topological metrics to describe pore space and fluid connectivity, with researchers demonstrating links between pore-scale nonwetting phase topology to fluid mobilization and displacement mechanisms, relative permeability, fluid flow regimes, and thermodynamic models of multiphase flow. While topology is clearly a powerful tool to describe fluid distribution, topological metrics by definition provide information only on the connectivity of a phase, not its geometry (shape or size). Physical flow characteristics, e.g. the permeability of a fluid phase within a porous medium, are dependent on the connectivity of the pore space or fluid phase as well as the size of connections. Persistent homology is a technique which provides a direct link between topology and geometry via measurement of topological features and their persistence from the signed Euclidean distance transform of a segmented digital image (Figure 1). We apply persistent homology analysis to measure the occurrence and size of pore-scale topological features in a variety of sandstones, for both the dry state and the nonwetting phase fluid during two-phase fluid flow (drainage and imbibition) experiments, visualized with 3D X-ray microtomography. The results provide key insights into the dominant topological features and length scales of a media which control relevant field-scale engineering properties such as fluid trapping, absolute permeability, and relative permeability.

Validation of model predictions of pore-scale fluid distributions during two-phase flow

NASA Astrophysics Data System (ADS)

Bultreys, Tom; Lin, Qingyang; Gao, Ying; Raeini, Ali Q.; AlRatrout, Ahmed; Bijeljic, Branko; Blunt, Martin J.

2018-05-01

Pore-scale two-phase flow modeling is an important technology to study a rock's relative permeability behavior. To investigate if these models are predictive, the calculated pore-scale fluid distributions which determine the relative permeability need to be validated. In this work, we introduce a methodology to quantitatively compare models to experimental fluid distributions in flow experiments visualized with microcomputed tomography. First, we analyzed five repeated drainage-imbibition experiments on a single sample. In these experiments, the exact fluid distributions were not fully repeatable on a pore-by-pore basis, while the global properties of the fluid distribution were. Then two fractional flow experiments were used to validate a quasistatic pore network model. The model correctly predicted the fluid present in more than 75% of pores and throats in drainage and imbibition. To quantify what this means for the relevant global properties of the fluid distribution, we compare the main flow paths and the connectivity across the different pore sizes in the modeled and experimental fluid distributions. These essential topology characteristics matched well for drainage simulations, but not for imbibition. This suggests that the pore-filling rules in the network model we used need to be improved to make reliable predictions of imbibition. The presented analysis illustrates the potential of our methodology to systematically and robustly test two-phase flow models to aid in model development and calibration.

Pore-Scale X-ray Micro-CT Imaging and Analysis of Oil Shales

NASA Astrophysics Data System (ADS)

Saif, T.

2015-12-01

The pore structure and the connectivity of the pore space during the pyrolysis of oil shales are important characteristics which determine hydrocarbon flow behaviour and ultimate recovery. We study the effect of temperature on the evolution of pore space and subsequent permeability on five oil shale samples: (1) Vernal Utah United States, (2) El Lajjun Al Karak Jordan, (3) Gladstone Queensland Australia (4) Fushun China and (5) Kimmerdige United Kingdom. Oil Shale cores of 5mm in diameter were pyrolized at 300, 400 and 500 °C. 3D imaging of 5mm diameter core samples was performed at 1μm voxel resolution using X-ray micro computed tomography (CT) and the evolution of the pore structures were characterized. The experimental results indicate that the thermal decomposition of kerogen at high temperatures is a major factor causing micro-scale changes in the internal structure of oil shales. At the early stage of pyrolysis, micron-scale heterogeneous pores were formed and with a further increase in temperature, the pores expanded and became interconnected by fractures. Permeability for each oil shale sample at each temperature was computed by simulation directly on the image voxels and by pore network extraction and simulation. Future work will investigate different samples and pursue insitu micro-CT imaging of oil shale pyrolysis to characterize the time evolution of the pore space.

Pore-Scale Determination of Gas Relative Permeability in Hydrate-Bearing Sediments Using X-Ray Computed Micro-Tomography and Lattice Boltzmann Method

NASA Astrophysics Data System (ADS)

Chen, Xiongyu; Verma, Rahul; Espinoza, D. Nicolas; Prodanović, Maša.

2018-01-01

This work uses X-ray computed micro-tomography (μCT) to monitor xenon hydrate growth in a sandpack under the excess gas condition. The μCT images give pore-scale hydrate distribution and pore habit in space and time. We use the lattice Boltzmann method to calculate gas relative permeability (krg) as a function of hydrate saturation (Shyd) in the pore structure of the experimental hydrate-bearing sand retrieved from μCT data. The results suggest the krg - Shyd data fit well a new model krg = (1-Shyd)·exp(-4.95·Shyd) rather than the simple Corey model. In addition, we calculate krg-Shyd curves using digital models of hydrate-bearing sand based on idealized grain-attaching, coarse pore-filling, and dispersed pore-filling hydrate habits. Our pore-scale measurements and modeling show that the krg-Shyd curves are similar regardless of whether hydrate crystals develop grain-attaching or coarse pore-filling habits. The dispersed pore filling habit exhibits much lower gas relative permeability than the other two, but it is not observed in the experiment and not compatible with Ostwald ripening mechanisms. We find that a single grain-shape factor can be used in the Carman-Kozeny equation to calculate krg-Shyd data with known porosity and average grain diameter, suggesting it is a useful model for hydrate-bearing sand.

Soil Pore Characteristics, an Underappreciated Regulatory Factor in GHGs Emission and C Stabilization

NASA Astrophysics Data System (ADS)

Toosi, E. R.; Yu, J.; Doane, T. A.; Guber, A.; Rivers, M. L.; Marsh, T. L.; Ali, K.; Kravchenko, A. N.

2015-12-01

Enduring challenges in understanding soil organic matter (SOM) stability and emission of greenhouse gases (GHGs) from soil stem from complexities of soil processes, many of which occur at micro-scales. The goal of this study is to evaluate the interactive effects soil pore characteristics, soil moisture levels, inherent SOM levels and properties, and substrate quality, on GHGs emission, and accelerated decomposition of native SOM following addition of fresh substrate i.e. priming. Our core hypothesis is that soil pore characteristics play a major role as a mediator in (i) the decomposition of organic matter regardless of its source (i.e. litter vs. native SOM) or substrate quality, as well as in (ii) GHGs emissions. Samples with prevalence of small (<10 μm) vs. large (>30 μm) pores were prepared from soils with similar properties but under long-term contrasting management. The samples were incubated (110 d) at low and optimum soil moisture conditions after addition of high quality (13C-soybean) and low quality (13C-corn) substrate. Headspace gas was analyzed for 13C-CO2 and GHGs on a regularly basis (day 1, 3, 7, 14, 24, 36, 48, 60, 72, 90, and 110). Selected samples were scanned at the early stage of decomposition (7, 14, 24 d) at 2-6 μm resolutions using X-ray computed μ tomography in order to: (1) quantify soil pore characteristics; (2) visualize and quantify distribution of soil moisture within samples of different pore characteristics; and (3) to visualize and measure losses of decomposing plant residue. Initial findings indicate that, consistent with our hypotheses, pore characteristics influenced GHGs emission, and intensity and pattern of plant residue decomposition. The importance of pores was highly pronounced in presence of added plant residue where greater N2O emission occurred in samples with dominant large pores, in contrast to CO2. Further findings will be discussed upon completion of the study and analysis of the results.

Multi-Scale Multi-Physics Modeling of Matrix Transport Properties in Fractured Shale Reservoirs

NASA Astrophysics Data System (ADS)

Mehmani, A.; Prodanovic, M.

2014-12-01

Understanding the shale matrix flow behavior is imperative in successful reservoir development for hydrocarbon production and carbon storage. Without a predictive model, significant uncertainties in flowback from the formation, the communication between the fracture and matrix as well as proper fracturing practice will ensue. Informed by SEM images, we develop deterministic network models that couple pores from multiple scales and their respective fluid physics. The models are used to investigate sorption hysteresis as an affordable way of inferring the nanoscale pore structure in core scale. In addition, restricted diffusion as a function of pore shape, pore-throat size ratios and network connectivity is computed to make correct interpretation of the 2D NMR maps possible. Our novel pore network models have the ability to match sorption hysteresis measurements without any tuning parameters. The results clarify a common misconception of linking type 3 nitrogen hysteresis curves to only the shale pore shape and show promising sensitivty for nanopore structre inference in core scale. The results on restricted diffusion shed light on the importance of including shape factors in 2D NMR interpretations. A priori "weighting factors" as a function of pore-throat and throat-length ratio are presented and the effect of network connectivity on diffusion is quantitatively assessed. We are currently working on verifying our models with experimental data gathered from the Eagleford formation.

Pore-scale spectral induced polarization (SIP) signaturesassociated with FeS biomineral transformations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Slater, Lee; Ntarlagiannis, Dimitrios; Personna, Yves R.

2007-10-01

The authors measured Spectral Induced Polarization (SIP) signatures in sand columns during (1) FeS biomineralization produced by sulfate reducing bacteria (D. vulgaris) under anaerboci conditions, and (2) subsequent biomineral dissolution upon return to an aerobic state. The low-frequency (0.1-10 Hz peak) relaxations produced during biomineralization can be modeled with a Cole-Cole formulation, from which the evolution of the polarization magnitude and relaxation length scale can be estimated. They find that the modeled time constant is consistent with the polarizable elements being biomineral encrused pores. Evolution of the model parameters is consistent with FeS surface area increases and pore-size reduction duringmore » biomineral growth, and subsequent biomineral dissolution (FeS surface area decreases and pore expansion) upon return to the aerobic state. They conclude that SIP signatures are diagnostic of pore-scale geometrical changes associated with FeS biomineralization by sulfate reducing bacteria.« less

Water permeability in hydrate-bearing sediments: A pore-scale study

NASA Astrophysics Data System (ADS)

Dai, Sheng; Seol, Yongkoo

2014-06-01

Permeability is a critical parameter governing methane flux and fluid flow in hydrate-bearing sediments; however, limited valid data are available due to experimental challenges. Here we investigate the relationship between apparent water permeability (k') and hydrate saturation (Sh), accounting for hydrate pore-scale growth habit and meso-scale heterogeneity. Results from capillary tube models rely on cross-sectional tube shapes and hydrate pore habits, thus are appropriate only for sediments with uniform hydrate distribution and known hydrate pore character. Given our pore network modeling results showing that accumulating hydrate in sediments decreases sediment porosity and increases hydraulic tortuosity, we propose a modified Kozeny-Carman model to characterize water permeability in hydrate-bearing sediments. This model agrees well with experimental results and can be easily implemented in reservoir simulators with no empirical variables other than Sh. Results are also relevant to flow through other natural sediments that undergo diagenesis, salt precipitation, or bio-clogging.

Micro/Nano-pore Network Analysis of Gas Flow in Shale Matrix

PubMed Central

Zhang, Pengwei; Hu, Liming; Meegoda, Jay N.; Gao, Shengyan

2015-01-01

The gas flow in shale matrix is of great research interests for optimized shale gas extraction. The gas flow in the nano-scale pore may fall in flow regimes such as viscous flow, slip flow and Knudsen diffusion. A 3-dimensional nano-scale pore network model was developed to simulate dynamic gas flow, and to describe the transient properties of flow regimes. The proposed pore network model accounts for the various size distributions and low connectivity of shale pores. The pore size, pore throat size and coordination number obey normal distribution, and the average values can be obtained from shale reservoir data. The gas flow regimes were simulated using an extracted pore network backbone. The numerical results show that apparent permeability is strongly dependent on pore pressure in the reservoir and pore throat size, which is overestimated by low-pressure laboratory tests. With the decrease of reservoir pressure, viscous flow is weakening, then slip flow and Knudsen diffusion are gradually becoming dominant flow regimes. The fingering phenomenon can be predicted by micro/nano-pore network for gas flow, which provides an effective way to capture heterogeneity of shale gas reservoir. PMID:26310236

Micro/Nano-pore Network Analysis of Gas Flow in Shale Matrix.

PubMed

Zhang, Pengwei; Hu, Liming; Meegoda, Jay N; Gao, Shengyan

2015-08-27

The gas flow in shale matrix is of great research interests for optimized shale gas extraction. The gas flow in the nano-scale pore may fall in flow regimes such as viscous flow, slip flow and Knudsen diffusion. A 3-dimensional nano-scale pore network model was developed to simulate dynamic gas flow, and to describe the transient properties of flow regimes. The proposed pore network model accounts for the various size distributions and low connectivity of shale pores. The pore size, pore throat size and coordination number obey normal distribution, and the average values can be obtained from shale reservoir data. The gas flow regimes were simulated using an extracted pore network backbone. The numerical results show that apparent permeability is strongly dependent on pore pressure in the reservoir and pore throat size, which is overestimated by low-pressure laboratory tests. With the decrease of reservoir pressure, viscous flow is weakening, then slip flow and Knudsen diffusion are gradually becoming dominant flow regimes. The fingering phenomenon can be predicted by micro/nano-pore network for gas flow, which provides an effective way to capture heterogeneity of shale gas reservoir.

Pore-scale observation and 3D simulation of wettability effects on supercritical CO2 - brine immiscible displacement in drainage

NASA Astrophysics Data System (ADS)

Hu, R.; Wan, J.; Chen, Y.

2016-12-01

Wettability is a factor controlling the fluid-fluid displacement pattern in porous media and significantly affects the flow and transport of supercritical (sc) CO2 in geologic carbon sequestration. Using a high-pressure micromodel-microscopy system, we performed drainage experiments of scCO2 invasion into brine-saturated water-wet and intermediate-wet micromodels; we visualized the scCO2 invasion morphology at pore-scale under reservoir conditions. We also performed pore-scale numerical simulations of the Navier-Stokes equations to obtain 3D details of fluid-fluid displacement processes. Simulation results are qualitatively consistent with the experiments, showing wider scCO2 fingering, higher percentage of scCO2 and more compact displacement pattern in intermediate-wet micromodel. Through quantitative analysis based on pore-scale simulation, we found that the reduced wettability reduces the displacement front velocity, promotes the pore-filling events in the longitudinal direction, delays the breakthrough time of invading fluid, and then increases the displacement efficiency. Simulated results also show that the fluid-fluid interface area follows a unified power-law relation with scCO2 saturation, and show smaller interface area in intermediate-wet case which suppresses the mass transfer between the phases. These pore-scale results provide insights for the wettability effects on CO2 - brine immiscible displacement in geologic carbon sequestration.

Investigation of Coupled model of Pore network and Continuum in shale gas

NASA Astrophysics Data System (ADS)

Cao, G.; Lin, M.

2016-12-01

Flow in shale spanning over many scales, makes the majority of conventional treatment methods disabled. For effectively simulating, a coupled model of pore-scale and continuum-scale was proposed in this paper. Based on the SEM image, we decompose organic-rich-shale into two subdomains: kerogen and inorganic matrix. In kerogen, the nanoscale pore-network is the main storage space and migration pathway so that the molecular phenomena (slip and diffusive transport) is significant. Whereas, inorganic matrix, with relatively large pores and micro fractures, the flow is approximate to Darcy. We use pore-scale network models (PNM) to represent kerogen and continuum-scale models (FVM or FEM) to represent matrix. Finite element mortars are employed to couple pore- and continuum-scale models by enforcing continuity of pressures and fluxes at shared boundary interfaces. In our method, the process in the coupled model is described by pressure square equation, and uses Dirichlet boundary conditions. We discuss several problems: the optimal element number of mortar faces, two categories boundary faces of pore network, the difference between 2D and 3D models, and the difference between continuum models FVM and FEM in mortars. We conclude that: (1) too coarse mesh in mortars will decrease the accuracy, while too fine mesh will lead to an ill-condition even singular system, the optimal element number is depended on boundary pores and nodes number. (2) pore network models are adjacent to two different mortar faces (PNM to PNM, PNM to continuum model), incidental repeated mortar nodes must be deleted. (3) 3D models can be replaced by 2D models under certain condition. (4) FVM is more convenient than FEM, for its simplicity in assigning interface nodes pressure and calculating interface fluxes. This work is supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB10020302), the 973 Program (2014CB239004), the Key Instrument Developing Project of the CAS (ZDYZ2012-1-08-02), the National Natural Science Foundation of China (41574129).

High-resolution face verification using pore-scale facial features.

PubMed

Li, Dong; Zhou, Huiling; Lam, Kin-Man

2015-08-01

Face recognition methods, which usually represent face images using holistic or local facial features, rely heavily on alignment. Their performances also suffer a severe degradation under variations in expressions or poses, especially when there is one gallery per subject only. With the easy access to high-resolution (HR) face images nowadays, some HR face databases have recently been developed. However, few studies have tackled the use of HR information for face recognition or verification. In this paper, we propose a pose-invariant face-verification method, which is robust to alignment errors, using the HR information based on pore-scale facial features. A new keypoint descriptor, namely, pore-Principal Component Analysis (PCA)-Scale Invariant Feature Transform (PPCASIFT)-adapted from PCA-SIFT-is devised for the extraction of a compact set of distinctive pore-scale facial features. Having matched the pore-scale features of two-face regions, an effective robust-fitting scheme is proposed for the face-verification task. Experiments show that, with one frontal-view gallery only per subject, our proposed method outperforms a number of standard verification methods, and can achieve excellent accuracy even the faces are under large variations in expression and pose.

Investigating the relative permeability behavior of microporosity-rich carbonates and tight sandstones with multiscale pore network models

NASA Astrophysics Data System (ADS)

Bultreys, Tom; Stappen, Jeroen Van; Kock, Tim De; Boever, Wesley De; Boone, Marijn A.; Hoorebeke, Luc Van; Cnudde, Veerle

2016-11-01

The relative permeability behavior of rocks with wide ranges of pore sizes is in many cases still poorly understood and is difficult to model at the pore scale. In this work, we investigate the capillary pressure and relative permeability behavior of three outcrop carbonates and two tight reservoir sandstones with wide, multimodal pore size distributions. To examine how the drainage and imbibition properties of these complex rock types are influenced by the connectivity of macropores to each other and to zones with unresolved small-scale porosity, we apply a previously presented microcomputed-tomography-based multiscale pore network model to these samples. The sensitivity to the properties of the small-scale porosity is studied by performing simulations with different artificial sphere-packing-based networks as a proxy for these pores. Finally, the mixed-wet water-flooding behavior of the samples is investigated, assuming different wettability distributions for the microporosity and macroporosity. While this work is not an attempt to perform predictive modeling, it seeks to qualitatively explain the behavior of the investigated samples and illustrates some of the most recent developments in multiscale pore network modeling.

Pore space analysis of NAPL distribution in sand-clay media

USGS Publications Warehouse

Matmon, D.; Hayden, N.J.

2003-01-01

This paper introduces a conceptual model of clays and non-aqueous phase liquids (NAPLs) at the pore scale that has been developed from a mathematical unit cell model, and direct micromodel observation and measurement of clay-containing porous media. The mathematical model uses a unit cell concept with uniform spherical grains for simulating the sand in the sand-clay matrix (???10% clay). Micromodels made with glass slides and including different clay-containing porous media were used to investigate the two clays (kaolinite and montmorillonite) and NAPL distribution within the pore space. The results were used to understand the distribution of NAPL advancing into initially saturated sand and sand-clay media, and provided a detailed analysis of the pore-scale geometry, pore size distribution, NAPL entry pressures, and the effect of clay on this geometry. Interesting NAPL saturation profiles were observed as a result of the complexity of the pore space geometry with the different packing angles and the presence of clays. The unit cell approach has applications for enhancing the mechanistic understanding and conceptualization, both visually and mathematically, of pore-scale processes such as NAPL and clay distribution. ?? 2003 Elsevier Science Ltd. All rights reserved.

Computational study of 3-D hot-spot initiation in shocked insensitive high-explosive

NASA Astrophysics Data System (ADS)

Najjar, F. M.; Howard, W. M.; Fried, L. E.; Manaa, M. R.; Nichols, A., III; Levesque, G.

2012-03-01

High-explosive (HE) material consists of large-sized grains with micron-sized embedded impurities and pores. Under various mechanical/thermal insults, these pores collapse generating hightemperature regions leading to ignition. A hydrodynamic study has been performed to investigate the mechanisms of pore collapse and hot spot initiation in TATB crystals, employing a multiphysics code, ALE3D, coupled to the chemistry module, Cheetah. This computational study includes reactive dynamics. Two-dimensional high-resolution large-scale meso-scale simulations have been performed. The parameter space is systematically studied by considering various shock strengths, pore diameters and multiple pore configurations. Preliminary 3-D simulations are undertaken to quantify the 3-D dynamics.

Towards a rigorous mesoscale modeling of reactive flow and transport in an evolving porous medium and its applications to soil science

NASA Astrophysics Data System (ADS)

Ray, Nadja; Rupp, Andreas; Knabner, Peter

2016-04-01

Soil is arguably the most prominent example of a natural porous medium that is composed of a porous matrix and a pore space. Within this framework and in terms of soil's heterogeneity, we first consider transport and fluid flow at the pore scale. From there, we develop a mechanistic model and upscale it mathematically to transfer our model from the small scale to that of the mesoscale (laboratory scale). The mathematical framework of (periodic) homogenization (in principal) rigorously facilitates such processes by exactly computing the effective coefficients/parameters by means of the pore geometry and processes. In our model, various small-scale soil processes may be taken into account: molecular diffusion, convection, drift emerging from electric forces, and homogeneous reactions of chemical species in a solvent. Additionally, our model may consider heterogeneous reactions at the porous matrix, thus altering both the porosity and the matrix. Moreover, our model may additionally address biophysical processes, such as the growth of biofilms and how this affects the shape of the pore space. Both of the latter processes result in an intrinsically variable soil structure in space and time. Upscaling such models under the assumption of a locally periodic setting must be performed meticulously to preserve information regarding the complex coupling of processes in the evolving heterogeneous medium. Generally, a micro-macro model emerges that is then comprised of several levels of couplings: Macroscopic equations that describe the transport and fluid flow at the scale of the porous medium (mesoscale) include averaged time- and space-dependent coefficient functions. These functions may be explicitly computed by means of auxiliary cell problems (microscale). Finally, the pore space in which the cell problems are defined is time- and space dependent and its geometry inherits information from the transport equation's solutions. Numerical computations using mixed finite elements and potentially random initial data, e.g. that of porosity, complement our theoretical results. Our investigations contribute to the theoretical understanding of the link between soil formation and soil functions. This general framework may be applied to various problems in soil science for a range of scales, such as the formation and turnover of microaggregates or soil remediation.

A New Dual-Pore Formation Factor Model: A Percolation Network Study and Comparison to Experimental Data

NASA Astrophysics Data System (ADS)

Tang, Y. B.; Li, M.; Bernabe, Y.

2014-12-01

We modeled the electrical transport behavior of dual-pore carbonate rocks in this paper. Based on experimental data of a carbonate reservoir in China, we simply considered the low porosity samples equivalent to the matrix (micro-pore system) of the high porosity samples. For modeling the bimodal porous media, we considered that the matrix is homogeneous and interconnected. The connectivity and the pore size distribution of macro-pore system are varied randomly. Both pore systems are supposed to act electrically in parallel, connected at the nodes, where the fluid exchange takes place, an approach previously used by Bauer et al. (2012). Then, the effect of the properties of matrix, the pore size distribution and connectivity of macro-pore system on petrophysical properties of carbonates can be investigated. We simulated electrical current through networks in three-dimensional simple cubic (SC) and body-center cubic (BCC) with different coordination numbers and different pipe radius distributions of macro-pore system. Based on the simulation results, we found that the formation factor obeys a "universal" scaling relationship (i.e. independent of lattice type), 1/F∝eγz, where γ is a function of the normalized standard deviation of the pore radius distribution of macro-pore system and z is the coordination number of macro-pore system. This relationship is different from the classic "universal power law" in percolation theory. A formation factor model was inferred on the basis of the scaling relationship mentioned above and several scale-invariant quantities (such as hydraulic radius rH and throat length l of macro-pore). Several methods were developed to estimate corresponding parameters of the new model with conventional core analyses. It was satisfactorily tested against experimental data, including some published experimental data. Furthermore, the relationship between water saturation and resistivity in dual-pore carbonates was discussed based on the new model.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yoon, Hongkyu; Major, Jonathan; Dewers, Thomas

Dissolved CO 2 in the subsurface resulting from geological CO 2 storage may react with minerals in fractured rocks, confined aquifers, or faults, resulting in mineral precipitation and dissolution. The overall rate of reaction can be affected by coupled processes including hydrodynamics, transport, and reactions at the (sub) pore-scale. In this work pore-scale modeling of coupled fluid flow, reactive transport, and heterogeneous reactions at the mineral surface is applied to account for permeability alterations caused by precipitation-induced pore-blocking. This paper is motivated by observations of CO 2 seeps from a natural CO 2 sequestration analog, Crystal Geyser, Utah. Observations alongmore » the surface exposure of the Little Grand Wash fault indicate the lateral migration of CO 2 seep sites (i.e., alteration zones) of 10–50 m width with spacing on the order of ~100 m over time. Sandstone permeability in alteration zones is reduced by 3–4 orders of magnitude by carbonate cementation compared to unaltered zones. One granular porous medium and one fracture network systems are used to conceptually represent permeable porous media and locations of conduits controlled by fault-segment intersections and/or topography, respectively. Simulation cases accounted for a range of reaction regimes characterized by the Damköhler (Da) and Peclet (Pe) numbers. Pore-scale simulation results demonstrate that combinations of transport (Pe), geochemical conditions (Da), solution chemistry, and pore and fracture configurations contributed to match key patterns observed in the field of how calcite precipitation alters flow paths by pore plugging. This comparison of simulation results with field observations reveals mechanistic explanations of the lateral migration and enhances our understanding of subsurface processes associated with the CO 2 injection. In addition, permeability and porosity relations are constructed from pore-scale simulations which account for a range of reaction regimes characterized by the Da and Pe numbers. Finally, the functional relationships obtained from pore-scale simulations can be used in a continuum scale model that may account for large-scale phenomena mimicking lateral migration of surface CO 2 seeps.« less

New insights on the complex dynamics of two-phase flow in porous media under intermediate-wet conditions.

PubMed

Rabbani, Harris Sajjad; Joekar-Niasar, Vahid; Pak, Tannaz; Shokri, Nima

2017-07-04

Multiphase flow in porous media is important in a number of environmental and industrial applications such as soil remediation, CO 2 sequestration, and enhanced oil recovery. Wetting properties control flow of immiscible fluids in porous media and fluids distribution in the pore space. In contrast to the strong and weak wet conditions, pore-scale physics of immiscible displacement under intermediate-wet conditions is less understood. This study reports the results of a series of two-dimensional high-resolution direct numerical simulations with the aim of understanding the pore-scale dynamics of two-phase immiscible fluid flow under intermediate-wet conditions. Our results show that for intermediate-wet porous media, pore geometry has a strong influence on interface dynamics, leading to co-existence of concave and convex interfaces. Intermediate wettability leads to various interfacial movements which are not identified under imbibition or drainage conditions. These pore-scale events significantly influence macro-scale flow behaviour causing the counter-intuitive decline in recovery of the defending fluid from weak imbibition to intermediate-wet conditions.

Pore Scale Investigation of Wettability Alteration Through Chemically-Tuned Waterflooding in Oil-Wet Carbonate Rocks Using X-Ray Micro-Ct Imaging

NASA Astrophysics Data System (ADS)

Tawfik, M. S.; Karpyn, Z.

2017-12-01

Carbonate reservoirs host more than half of the remaining oil reserves worldwide. Due to their complex pore structure and intermediate to oil-wet nature, it is challenging to produce the remaining oil from these formations. For two decades, chemically tuned waterflooding (CTWF) has gained the attention of many researchers. Experimental, numerical, and field studies suggest that changes in ion composition of injected brine can increase oil recovery in carbonate reservoirs via wettability alteration. However, previous studies explaining the improvement in oil recovery by wettability alteration deduce wettability based on indirect measurements, including sessile drop contact angle measurements on polished rocks, relative permeability, chromatographic separation of SCN- and potential determining ions (PDIs), etc. CTWF literature offers no direct measurement of wettability alteration at the pore scale. This study proposes a direct pore-scale measurement of changes in interfacial curvatures before and after CTWF. Micro-coreflood experiments are performed to investigate the effect of injection brine salinity, ion composition and temperature on rock wettability at the pore scale. X-ray micro-CT scanning is used to obtain 3D image sets to calculate in-situ contact angle distributions. The study also aims to find a correlation between the magnitude of improvement in oil recovery at the macro-scale and the corresponding contact angle distribution at the pore-scale at different experimental conditions. Hence, macro-scale coreflood experiments are performed using the same conditions as the micro-corefloods. Macro-scale coreflood experiments have shown that brines with higher concentration of Ca2+, Mg2+ and SO42- ions have higher recoveries compared to standard seawater. This translates to wettability alteration into a more intermediate-wet state. This study enhances the understanding of the pore-scale physico-chemical mechanisms controlling wettability alteration via CTWF, which helps tune existing CTWF models, and therefore results in more well-informed candidate reservoir selection and the development of a workflow to determine the optimum injection brine properties for a given crude oil-brine-rock system.

Multi-scale modeling of multi-component reactive transport in geothermal aquifers

NASA Astrophysics Data System (ADS)

Nick, Hamidreza M.; Raoof, Amir; Wolf, Karl-Heinz; Bruhn, David

2014-05-01

In deep geothermal systems heat and chemical stresses can cause physical alterations, which may have a significant effect on flow and reaction rates. As a consequence it will lead to changes in permeability and porosity of the formations due to mineral precipitation and dissolution. Large-scale modeling of reactive transport in such systems is still challenging. A large area of uncertainty is the way in which the pore-scale information controlling the flow and reaction will behave at a larger scale. A possible choice is to use constitutive relationships relating, for example the permeability and porosity evolutions to the change in the pore geometry. While determining such relationships through laboratory experiments may be limited, pore-network modeling provides an alternative solution. In this work, we introduce a new workflow in which a hybrid Finite-Element Finite-Volume method [1,2] and a pore network modeling approach [3] are employed. Using the pore-scale model, relevant constitutive relations are developed. These relations are then embedded in the continuum-scale model. This approach enables us to study non-isothermal reactive transport in porous media while accounting for micro-scale features under realistic conditions. The performance and applicability of the proposed model is explored for different flow and reaction regimes. References: 1. Matthäi, S.K., et al.: Simulation of solute transport through fractured rock: a higher-order accurate finite-element finite-volume method permitting large time steps. Transport in porous media 83.2 (2010): 289-318. 2. Nick, H.M., et al.: Reactive dispersive contaminant transport in coastal aquifers: Numerical simulation of a reactive Henry problem. Journal of contaminant hydrology 145 (2012), 90-104. 3. Raoof A., et al.: PoreFlow: A Complex pore-network model for simulation of reactive transport in variably saturated porous media, Computers & Geosciences, 61, (2013), 160-174.

Concentration of Immunoglobulins in Microfiltration Permeates of Skim Milk: Impact of Transmembrane Pressure and Temperature on the IgG Transmission Using Different Ceramic Membrane Types and Pore Sizes.

PubMed

Heidebrecht, Hans-Jürgen; Toro-Sierra, José; Kulozik, Ulrich

2018-06-28

The use of bioactive bovine milk immunoglobulins (Ig) has been found to be an alternative treatment for certain human gastrointestinal diseases. Some methodologies have been developed with bovine colostrum. These are considered in laboratory scale and are bound to high cost and limited availability of the raw material. The main challenge remains in obtaining high amounts of active IgG from an available source as mature cow milk by the means of industrial processes. Microfiltration (MF) was chosen as a process variant, which enables a gentle and effective concentration of the Ig fractions (ca. 0.06% in raw milk) while reducing casein and lactose at the same time. Different microfiltration membranes (ceramic standard and gradient), pore sizes (0.14⁻0.8 µm), transmembrane pressures (0.5⁻2.5 bar), and temperatures (10, 50 °C) were investigated. The transmission of immunoglobulin G (IgG) and casein during the filtration of raw skim milk (<0.1% fat) was evaluated during batch filtration using a single channel pilot plant. The transmission levels of IgG (~160 kDa) were measured to be at the same level as the reference major whey protein β-Lg (~18 kDa) at all evaluated pore sizes and process parameters despite the large difference in molecular mass of both fractions. Ceramic gradient membranes with a pore sizes of 0.14 µm showed IgG-transmission rates between 45% to 65% while reducing the casein fraction below 1% in the permeates. Contrary to the expectations, a lower pore size of 0.14 µm yielded fluxes up to 35% higher than 0.2 µm MF membranes. It was found that low transmembrane pressures benefit the Ig transmission. Upscaling the presented results to a continuous MF membrane process offers new possibilities for the production of immunoglobulin enriched supplements with well-known processing equipment for large scale milk protein fractionation.

Toward multiscale modelings of grain-fluid systems

NASA Astrophysics Data System (ADS)

Chareyre, Bruno; Yuan, Chao; Montella, Eduard P.; Salager, Simon

2017-06-01

Computationally efficient methods have been developed for simulating partially saturated granular materials in the pendular regime. In contrast, one hardly avoid expensive direct resolutions of 2-phase fluid dynamics problem for mixed pendular-funicular situations or even saturated regimes. Following previous developments for single-phase flow, a pore-network approach of the coupling problems is described. The geometry and movements of phases and interfaces are described on the basis of a tetrahedrization of the pore space, introducing elementary objects such as bridge, meniscus, pore body and pore throat, together with local rules of evolution. As firmly established local rules are still missing on some aspects (entry capillary pressure and pore-scale pressure-saturation relations, forces on the grains, or kinetics of transfers in mixed situations) a multi-scale numerical framework is introduced, enhancing the pore-network approach with the help of direct simulations. Small subsets of a granular system are extracted, in which multiphase scenario are solved using the Lattice-Boltzman method (LBM). In turns, a global problem is assembled and solved at the network scale, as illustrated by a simulated primary drainage.

Study into the correlation of dominant pore throat size and SIP relaxation frequency

NASA Astrophysics Data System (ADS)

Kruschwitz, Sabine; Prinz, Carsten; Zimathies, Annett

2016-12-01

There is currently a debate within the SIP community about the characteristic textural length scale controlling relaxation time of consolidated porous media. One idea is that the relaxation time is dominated by the pore throat size distribution or more specifically the modal pore throat size as determined in mercury intrusion capillary pressure tests. Recently new studies on inverting pore size distributions from SIP data were published implying that the relaxation mechanisms and controlling length scale are well understood. In contrast new analytical model studies based on the Marshall-Madden membrane polarization theory suggested that two relaxation processes might compete: the one along the short narrow pore (the throat) with one across the wider pore in case the narrow pores become relatively long. This paper presents a first systematically focused study into the relationship of pore throat sizes and SIP relaxation times. The generality of predicted trends is investigated across a wide range of materials differing considerably in chemical composition, specific surface and pore space characteristics. Three different groups of relaxation behaviors can be clearly distinguished. The different behaviors are related to clay content and type, carbonate content, size of the grains and the wide pores in the samples.

Mineral transformation controls speciation and pore-fluid transmission of contaminants in waste-weathered Hanford sediments

NASA Astrophysics Data System (ADS)

Perdrial, Nicolas; Thompson, Aaron; O'Day, Peggy A.; Steefel, Carl I.; Chorover, Jon

2014-09-01

Portions of the Hanford Site (WA, USA) vadose zone were subjected to weathering by caustic solutions during documented releases of high level radioactive waste (containing Sr, Cs and I) from leaking underground storage tanks. Previous studies have shown that waste-sediment interactions can promote variable incorporation of contaminants into neo-formed mineral products (including feldspathoids and zeolites), but processes regulating the subsequent contaminant release from these phases into infiltrating background pore waters remain poorly known. In this paper, reactive transport experiments were conducted with Hanford sediments previously weathered for one year in simulated hyper-alkaline waste solutions containing high or low 88Sr, 127I, and 133Cs concentrations, with or without CO2(aq). These waste-weathered sediments were leached in flow-through column experiments with simulated background pore water (characteristic of meteoric recharge) to measure contaminant release from solids formed during waste-sediment interaction. Contaminant sorption-desorption kinetics and mineral transformation reactions were both monitored using continuous-flow and wet-dry cycling regimes for ca. 300 pore volumes. Less than 20% of contaminant 133Cs and 88Sr mass and less than 40% 127I mass were released over the course of the experiment. To elucidate molecular processes limiting contaminant release, reacted sediments were studied with micro- (TEM and XRD) and molecular- (Sr K-edge EXAFS) scale methods. Contaminant dynamics in column experiments were principally controlled by rapid dissolution of labile solids and competitive exchange reactions. In initially feldspathoidic systems, time-dependent changes in the local zeolitic bonding environment observed with X-ray diffraction and EXAFS are responsible for limiting contaminant release. Linear combination fits and shell-by-shell analysis of Sr K-edge EXAFS data revealed modification in Sr-Si/Al distances within the zeolite cage. Wet-dry cycling did not affect significantly molecular-scale transformations relative to continuous-flow controls. Results indicate that contaminants bound to the solid phase in distinct micro- and molecular-scale coordinative environments can generate similar macro-scale release behaviors, highlighting the need for multi-scale interrogations to constrain mechanisms of reactive transport. Data also indicate that weathering-induced change in ion exchange selectivity coefficients should be incorporated in simulations of contaminant release from caustic high-level radioactive waste impacted sediments.

The effects of burial diagenesis on multiscale porosity in the St. Peter Sandstone: An imaging, small-angle, and ultra-small-angle neutron scattering analysis

DOE PAGES

Anovitz, Lawrence M.; Freiburg, Jared T.; Wasbrough, Matthew; ...

2017-11-06

To examine the effects of burial diagenesis on heirarchical pore structures in sandstone and compare those with the effects of overgrowth formation, we obtained samples of St. Peter Sandstone from drill cores obtained in the Illinois and Michigan Basins. The multiscale pore structure of rocks in sedimentary reservoirs and the mineralogy associated with those pores are critical factors for estimating reservoir properties, including fluid mass in place, permeability, and capillary pressures, as well as geochemical interactions between the rock and the fluid. The combination of small- and ultra-small-angle neutron scattering with backscattered electron or X ray-computed tomographic imaging, or both,more » provided a means by which pore structures were quantified at scales ranging from aproximately 1 nm to 1 cm—seven orders of magnitude. Larger scale (>10 µm) porosity showed the expected logarithmic decrease in porosity with depth, although there was significant variation in each sample group. However, small- and ultra-small-angle neutron scattering data showed that the proportion of small-scale porosity increased with depth. Porosity distributions were not continuous, but consisted of a series of log normal-like distributions at several distinct scales within these rocks. Fractal dimensions at larger scales decreased (surfaces smoothed) with increasing depth, and those at smaller scales increased (surfaces roughened) and pores become more isolated (higher lacunarity). Furthermore, data suggest that changes in pore-size distributions are controlled by both physical (compaction) and chemical effects (precipitation, cementation, dissolution).« less

The effects of burial diagenesis on multiscale porosity in the St. Peter Sandstone: An imaging, small-angle, and ultra-small-angle neutron scattering analysis

DOE Office of Scientific and Technical Information (OSTI.GOV)

Anovitz, Lawrence M.; Freiburg, Jared T.; Wasbrough, Matthew

To examine the effects of burial diagenesis on heirarchical pore structures in sandstone and compare those with the effects of overgrowth formation, we obtained samples of St. Peter Sandstone from drill cores obtained in the Illinois and Michigan Basins. The multiscale pore structure of rocks in sedimentary reservoirs and the mineralogy associated with those pores are critical factors for estimating reservoir properties, including fluid mass in place, permeability, and capillary pressures, as well as geochemical interactions between the rock and the fluid. The combination of small- and ultra-small-angle neutron scattering with backscattered electron or X ray-computed tomographic imaging, or both,more » provided a means by which pore structures were quantified at scales ranging from aproximately 1 nm to 1 cm—seven orders of magnitude. Larger scale (>10 µm) porosity showed the expected logarithmic decrease in porosity with depth, although there was significant variation in each sample group. However, small- and ultra-small-angle neutron scattering data showed that the proportion of small-scale porosity increased with depth. Porosity distributions were not continuous, but consisted of a series of log normal-like distributions at several distinct scales within these rocks. Fractal dimensions at larger scales decreased (surfaces smoothed) with increasing depth, and those at smaller scales increased (surfaces roughened) and pores become more isolated (higher lacunarity). Furthermore, data suggest that changes in pore-size distributions are controlled by both physical (compaction) and chemical effects (precipitation, cementation, dissolution).« less

Rapid Polymer Transport in a Single Nanometer-Scale Pore

NASA Astrophysics Data System (ADS)

Kasianowicz, J. J.

1998-03-01

Protein ion channels are nanometer-scale pores that control the transport of ions and polymers across cell membranes. We compared the ability of charged and nonelectrolyte linear polymers to partition into a single channel reconstituted into a planar lipid bilayer membrane. The entry of each polymer (e.g. monodisperse length single-stranded homopolymeric RNA1 or poly(ethylene glycol)2,3) into the pore caused characteristic transient decreases in the channel's ionic conductance. The ionic current blockades yield detailed information about the physical properties of the polymers and the pore. The biological and technological significance of the results will be discussed.

Machine learning framework for analysis of transport through complex networks in porous, granular media: A focus on permeability

NASA Astrophysics Data System (ADS)

van der Linden, Joost H.; Narsilio, Guillermo A.; Tordesillas, Antoinette

2016-08-01

We present a data-driven framework to study the relationship between fluid flow at the macroscale and the internal pore structure, across the micro- and mesoscales, in porous, granular media. Sphere packings with varying particle size distribution and confining pressure are generated using the discrete element method. For each sample, a finite element analysis of the fluid flow is performed to compute the permeability. We construct a pore network and a particle contact network to quantify the connectivity of the pores and particles across the mesoscopic spatial scales. Machine learning techniques for feature selection are employed to identify sets of microstructural properties and multiscale complex network features that optimally characterize permeability. We find a linear correlation (in log-log scale) between permeability and the average closeness centrality of the weighted pore network. With the pore network links weighted by the local conductance, the average closeness centrality represents a multiscale measure of efficiency of flow through the pore network in terms of the mean geodesic distance (or shortest path) between all pore bodies in the pore network. Specifically, this study objectively quantifies a hypothesized link between high permeability and efficient shortest paths that thread through relatively large pore bodies connected to each other by high conductance pore throats, embodying connectivity and pore structure.

Molecular mechanism for lipid flip-flops.

PubMed

Gurtovenko, Andrey A; Vattulainen, Ilpo

2007-12-06

Transmembrane lipid translocation (flip-flop) processes are involved in a variety of properties and functions of cell membranes, such as membrane asymmetry and programmed cell death. Yet, flip-flops are one of the least understood dynamical processes in membranes. In this work, we elucidate the molecular mechanism of pore-mediated transmembrane lipid translocation (flip-flop) acquired from extensive atomistic molecular dynamics simulations. On the basis of 50 successful flip-flop events resolved in atomic detail, we demonstrate that lipid flip-flops may spontaneously occur in protein-free phospholipid membranes under physiological conditions through transient water pores on a time scale of tens of nanoseconds. While the formation of a water pore is induced here by a transmembrane ion density gradient, the particular way by which the pore is formed is irrelevant for the reported flip-flop mechanism: the appearance of a transient pore (defect) in the membrane inevitably leads to diffusive translocation of lipids through the pore, which is driven by thermal fluctuations. Our findings strongly support the idea that the formation of membrane defects in terms of water pores is the rate-limiting step in the process of transmembrane lipid flip-flop, which, on average, requires several hours. The findings are consistent with available experimental and computational data and provide a view to interpret experimental observations. For example, the simulation results provide a molecular-level explanation in terms of pores for the experimentally observed fact that the exposure of lipid membranes to electric field pulses considerably reduces the time required for lipid flip-flops.

Pore-scale investigation on the response of heterotrophic respiration to moisture conditions in heterogeneous soils

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yan, Zhifeng; Liu, Chongxuan; Todd-Brown, Katherine E.

The relationship between microbial respiration rate and soil moisture content is an important property for understanding and predicting soil organic carbon degradation, CO 2 production and emission, and their subsequent effects on climate change. This paper reports a pore-scale modeling study to investigate the response of heterotrophic respiration to moisture conditions in soils and to evaluate various factors that affect this response. X-ray computed tomography was used to derive soil pore structures, which were then used for pore-scale model investigation. The pore-scale results were then averaged to calculate the effective respiration rates as a function of water content in soils.more » The calculated effective respiration rate first increases and then decreases with increasing soil water content, showing a maximum respiration rate at water saturation degree of 0.75 that is consistent with field and laboratory observations. The relationship between the respiration rate and moisture content is affected by various factors, including pore-scale organic carbon bioavailability, the rate of oxygen delivery, soil pore structure and physical heterogeneity, soil clay content, and microbial drought resistivity. Simulations also illustrates that a larger fraction of CO 2 produced from microbial respiration can be accumulated inside soil cores under higher saturation conditions, implying that CO 2 flux measured on the top of soil cores may underestimate or overestimate true soil respiration rates under dynamic moisture conditions. Overall, this study provides mechanistic insights into the soil respiration response to the change in moisture conditions, and reveals a complex relationship between heterotrophic microbial respiration rate and moisture content in soils that is affected by various hydrological, geochemical, and biophysical factors.« less

Capillary filling rules and displacement mechanisms for spontaneous imbibition of CO2 for carbon storage and EOR using micro-model experiments and pore scale simulation

NASA Astrophysics Data System (ADS)

Chapman, E.; Yang, J.; Crawshaw, J.; Boek, E. S.

2012-04-01

In the 1980s, Lenormand et al. carried out their pioneering work on displacement mechanisms of fluids in etched networks [1]. Here we further examine displacement mechanisms in relation to capillary filling rules for spontaneous imbibition. Understanding the role of spontaneous imbibition in fluid displacement is essential for refining pore network models. Generally, pore network models use simple capillary filling rules and here we examine the validity of these rules for spontaneous imbibition. Improvement of pore network models is vital for the process of 'up-scaling' to the field scale for both enhanced oil recovery (EOR) and carbon sequestration. In this work, we present our experimental microfluidic research into the displacement of both supercritical CO2/deionised water (DI) systems and analogous n-decane/air - where supercritical CO2 and n-decane are the respective wetting fluids - controlled by imbibition at the pore scale. We conducted our experiments in etched PMMA and silicon/glass micro-fluidic hydrophobic chips. We first investigate displacement in single etched pore junctions, followed by displacement in complex network designs representing actual rock thin sections, i.e. Berea sandstone and Sucrosic dolomite. The n-decane/air experiments were conducted under ambient conditions, whereas the supercritical CO2/DI water experiments were conducted under high temperature and pressure in order to replicate reservoir conditions. Fluid displacement in all experiments was captured via a high speed video microscope. The direction and type of displacement the imbibing fluid takes when it enters a junction is dependent on the number of possible channels in which the wetting fluid can imbibe, i.e. I1, I2 and I3 [1]. Depending on the experiment conducted, the micro-models were initially filled with either DI water or air before the wetting fluid was injected. We found that the imbibition of the wetting fluid through a single pore is primarily controlled by the geometry of the pore body rather than the downstream pore throat sizes, contrary to the established capillary filling rules as used in current pore network models. Our experimental observations are confirmed by detailed lattice-Boltzmann pore scale computer simulations of fluid displacement in the same geometries. This suggests that capillary filling rules for imbibition as used in pore network models may need to be revised. [1] G. Lenormand, C. Zarcone and A. Sarr, J. Fluid Mech. 135 , 337-353 (1983).

Lattice Boltzmann-Based Approaches for Pore-Scale Reactive Transport

DOE PAGES

Yoon, Hongkyu; Kang, Qinjun; Valocchi, Albert J.

2015-07-29

Here an important geoscience and environmental applications such as geologic carbon storage, environmental remediation, and unconventional oil and gas recovery are best understood in the context of reactive flow and multicomponent transport in the subsurface environment. The coupling of chemical and microbiological reactions with hydrological and mechanical processes can lead to complex behaviors across an enormous range of spatial and temporal scales. These coupled responses are also strongly influenced by the heterogeneity and anisotropy of the geologic formations. Reactive transport processes can change the pore morphology at the pore scale, thereby leading to nonlinear interactions with advective and diffusive transport,more » which can strongly influence larger-scale properties such as permeability and dispersion.« less

Pore scale study of multiphase multicomponent reactive transport during CO 2 dissolution trapping

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chen, Li; Wang, Mengyi; Kang, Qinjun

Solubility trapping is crucial for permanent CO 2 sequestration in deep saline aquifers. For the first time, a pore-scale numerical method is developed to investigate coupled scCO 2-water two-phase flow, multicomponent (CO 2(aq), H +, HCO 3 –, CO 3 2 – and OH –) mass transport, heterogeneous interfacial dissolution reaction, and homogeneous dissociation reactions. Pore-scale details of evolutions of multiphase distributions and concentration fields are presented and discussed. Time evolutions of several variables including averaged CO 2(aq) concentration, scCO 2 saturation, and pH value are analyzed. Specific interfacial length, an important variable which cannot be determined but is requiredmore » by continuum models, is investigated in detail. Mass transport coefficient or efficient dissolution rate is also evaluated. The pore-scale results show strong non-equilibrium characteristics during solubility trapping due to non-uniform distributions of multiphase as well as slow mass transport process. Complicated coupling mechanisms between multiphase flow, mass transport and chemical reactions are also revealed. Lastly, effects of wettability are also studied. The pore-scale studies provide deep understanding of non-linear non-equilibrium multiple physicochemical processes during CO 2 solubility trapping processes, and also allow to quantitatively predict some important empirical relationships, such as saturation-interfacial surface area, for continuum models.« less

Pore scale study of multiphase multicomponent reactive transport during CO 2 dissolution trapping

DOE PAGES

Chen, Li; Wang, Mengyi; Kang, Qinjun; ...

2018-04-26

Solubility trapping is crucial for permanent CO 2 sequestration in deep saline aquifers. For the first time, a pore-scale numerical method is developed to investigate coupled scCO 2-water two-phase flow, multicomponent (CO 2(aq), H +, HCO 3 –, CO 3 2 – and OH –) mass transport, heterogeneous interfacial dissolution reaction, and homogeneous dissociation reactions. Pore-scale details of evolutions of multiphase distributions and concentration fields are presented and discussed. Time evolutions of several variables including averaged CO 2(aq) concentration, scCO 2 saturation, and pH value are analyzed. Specific interfacial length, an important variable which cannot be determined but is requiredmore » by continuum models, is investigated in detail. Mass transport coefficient or efficient dissolution rate is also evaluated. The pore-scale results show strong non-equilibrium characteristics during solubility trapping due to non-uniform distributions of multiphase as well as slow mass transport process. Complicated coupling mechanisms between multiphase flow, mass transport and chemical reactions are also revealed. Lastly, effects of wettability are also studied. The pore-scale studies provide deep understanding of non-linear non-equilibrium multiple physicochemical processes during CO 2 solubility trapping processes, and also allow to quantitatively predict some important empirical relationships, such as saturation-interfacial surface area, for continuum models.« less

Pore scale study of multiphase multicomponent reactive transport during CO2 dissolution trapping

NASA Astrophysics Data System (ADS)

Chen, Li; Wang, Mengyi; Kang, Qinjun; Tao, Wenquan

2018-06-01

Solubility trapping is crucial for permanent CO2 sequestration in deep saline aquifers. For the first time, a pore-scale numerical method is developed to investigate coupled scCO2-water two-phase flow, multicomponent (CO2(aq), H+, HCO3-, CO32- and OH-) mass transport, heterogeneous interfacial dissolution reaction, and homogeneous dissociation reactions. Pore-scale details of evolutions of multiphase distributions and concentration fields are presented and discussed. Time evolutions of several variables including averaged CO2(aq) concentration, scCO2 saturation, and pH value are analyzed. Specific interfacial length, an important variable which cannot be determined but is required by continuum models, is investigated in detail. Mass transport coefficient or efficient dissolution rate is also evaluated. The pore-scale results show strong non-equilibrium characteristics during solubility trapping due to non-uniform distributions of multiphase as well as slow mass transport process. Complicated coupling mechanisms between multiphase flow, mass transport and chemical reactions are also revealed. Finally, effects of wettability are also studied. The pore-scale studies provide deep understanding of non-linear non-equilibrium multiple physicochemical processes during CO2 solubility trapping processes, and also allow to quantitatively predict some important empirical relationships, such as saturation-interfacial surface area, for continuum models.

Porous materials for thermal management under extreme conditions.

PubMed

Clyne, T W; Golosnoy, I O; Tan, J C; Markaki, A E

2006-01-15

A brief analysis is presented of how heat transfer takes place in porous materials of various types. The emphasis is on materials able to withstand extremes of temperature, gas pressure, irradiation, etc. i.e. metals and ceramics, rather than polymers. A primary aim is commonly to maximize either the thermal resistance (i.e. provide insulation) or the rate of thermal equilibration between the material and a fluid passing through it (i.e. to facilitate heat exchange). The main structural characteristics concern porosity (void content), anisotropy, pore connectivity and scale. The effect of scale is complex, since the permeability decreases as the structure is refined, but the interfacial area for fluid-solid heat exchange is, thereby, raised. The durability of the pore structure may also be an issue, with a possible disadvantage of finer scale structures being poor microstructural stability under service conditions. Finally, good mechanical properties may be required, since the development of thermal gradients, high fluid fluxes, etc. can generate substantial levels of stress. There are, thus, some complex interplays between service conditions, pore architecture/scale, fluid permeation characteristics, convective heat flow, thermal conduction and radiative heat transfer. Such interplays are illustrated with reference to three examples: (i) a thermal barrier coating in a gas turbine engine; (ii) a Space Shuttle tile; and (iii) a Stirling engine heat exchanger. Highly porous, permeable materials are often made by bonding fibres together into a network structure and much of the analysis presented here is oriented towards such materials.

Enhancing Tensile Response of Sn Using Cu at Nano Length Scale and High Temperature Extrusion

DTIC Science & Technology

2009-02-01

temperature extruded Sn-1.1Cu 664 samples suggesting the presence of lenticular pores. This aspect ratio of pores was only 1.7 for high temperature Sn...resulting in filling the voids or breaking the lenticular pores into small pores besides higher atomic diffusion rates [8...relatively round pores were observed for hot extruded Sn-Cu samples that helps to increase the strength. The lenticular pores (higher aspect ratio) in

Shifts in pore connectivity from precipitation versus groundwater rewetting increases soil carbon loss after drought

DOE Office of Scientific and Technical Information (OSTI.GOV)

Smith, A. Peyton; Bond-Lamberty, Ben; Benscoter, Brian W.

Droughts and other extreme precipitation events are predicted to increase in intensity, duration and extent, with uncertain implications for terrestrial carbon (C) sequestration. Soil wetting from above (precipitation) results in a characteristically different pattern of pore-filling than wetting from below (groundwater), with larger, well-connected pores filling before finer pore spaces, unlike groundwater rise in which capillary forces saturate the finest pores first. Here we demonstrate that pore-scale wetting patterns interact with antecedent soil moisture conditions to alter pore-, core- and field-scale C dynamics. Drought legacy and wetting direction are perhaps more important determinants of short-term C mineralization than current soilmore » moisture content in these soils. Our results highlight that microbial access to C is not solely limited by physical protection, but also by drought or wetting-induced shifts in hydrologic connectivity. We argue that models should treat soil moisture within a three-dimensional framework emphasizing hydrologic conduits for C and resource diffusion.« less

Macroscopic modeling of heat and water vapor transfer with phase change in dry snow based on an upscaling method: Influence of air convection

NASA Astrophysics Data System (ADS)

Calonne, N.; Geindreau, C.; Flin, F.

2015-12-01

At the microscopic scale, i.e., pore scale, dry snow metamorphism is mainly driven by the heat and water vapor transfer and the sublimation-deposition process at the ice-air interface. Up to now, the description of these phenomena at the macroscopic scale, i.e., snow layer scale, in the snowpack models has been proposed in a phenomenological way. Here we used an upscaling method, namely, the homogenization of multiple-scale expansions, to derive theoretically the macroscopic equivalent modeling of heat and vapor transfer through a snow layer from the physics at the pore scale. The physical phenomena under consideration are steady state air flow, heat transfer by conduction and convection, water vapor transfer by diffusion and convection, and phase change (sublimation and deposition). We derived three different macroscopic models depending on the intensity of the air flow considered at the pore scale, i.e., on the order of magnitude of the pore Reynolds number and the Péclet numbers: (A) pure diffusion, (B) diffusion and moderate convection (Darcy's law), and (C) strong convection (nonlinear flow). The formulation of the models includes the exact expression of the macroscopic properties (effective thermal conductivity, effective vapor diffusion coefficient, and intrinsic permeability) and of the macroscopic source terms of heat and vapor arising from the phase change at the pore scale. Such definitions can be used to compute macroscopic snow properties from 3-D descriptions of snow microstructures. Finally, we illustrated the precision and the robustness of the proposed macroscopic models through 2-D numerical simulations.

Pore-Scale Investigation on Stress-Dependent Characteristics of Granular Packs and Their Impact on Multiphase Fluid Distribution

NASA Astrophysics Data System (ADS)

Torrealba, V.; Karpyn, Z.; Yoon, H.; Hart, D. B.; Klise, K. A.

2013-12-01

The pore-scale dynamics that govern multiphase flow under variable stress conditions are not well understood. This lack of fundamental understanding limits our ability to quantitatively predict multiphase flow and fluid distributions in natural geologic systems. In this research, we focus on pore-scale, single and multiphase flow properties that impact displacement mechanisms and residual trapping of non-wetting phase under varying stress conditions. X-ray micro-tomography is used to image pore structures and distribution of wetting and non-wetting fluids in water-wet synthetic granular packs, under dynamic load. Micro-tomography images are also used to determine structural features such as medial axis, surface area, and pore body and throat distribution; while the corresponding transport properties are determined from Lattice-Boltzmann simulations performed on lattice replicas of the imaged specimens. Results are used to investigate how inter-granular deformation mechanisms affect fluid displacement and residual trapping at the pore-scale. This will improve our understanding of the dynamic interaction of mechanical deformation and fluid flow during enhanced oil recovery and geologic CO2 sequestration. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

Pores-scale hydrodynamics in a progressively bio-clogged three-dimensional porous medium: 3D particle tracking experiments and stochastic transport modelling

NASA Astrophysics Data System (ADS)

Morales, V. L.; Carrel, M.; Dentz, M.; Derlon, N.; Morgenroth, E.; Holzner, M.

2017-12-01

Biofilms are ubiquitous bacterial communities growing in various porous media including soils, trickling and sand filters and are relevant for applications such as the degradation of pollutants for bioremediation, waste water or drinking water production purposes. By their development, biofilms dynamically change the structure of porous media, increasing the heterogeneity of the pore network and the non-Fickian or anomalous dispersion. In this work, we use an experimental approach to investigate the influence of biofilm growth on pore scale hydrodynamics and transport processes and propose a correlated continuous time random walk model capturing these observations. We perform three-dimensional particle tracking velocimetry at four different time points from 0 to 48 hours of biofilm growth. The biofilm growth notably impacts pore-scale hydrodynamics, as shown by strong increase of the average velocity and in tailing of Lagrangian velocity probability density functions. Additionally, the spatial correlation length of the flow increases substantially. This points at the formation of preferential flow pathways and stagnation zones, which ultimately leads to an increase of anomalous transport in the porous media considered, characterized by non-Fickian scaling of mean-squared displacements and non-Gaussian distributions of the displacement probability density functions. A gamma distribution provides a remarkable approximation of the bulk and the high tail of the Lagrangian pore-scale velocity magnitude, indicating a transition from a parallel pore arrangement towards a more serial one. Finally, a correlated continuous time random walk based on a stochastic relation velocity model accurately reproduces the observations and could be used to predict transport beyond the time scales accessible to the experiment.

Investigating Darcy-scale assumptions by means of a multiphysics algorithm

NASA Astrophysics Data System (ADS)

Tomin, Pavel; Lunati, Ivan

2016-09-01

Multiphysics (or hybrid) algorithms, which couple Darcy and pore-scale descriptions of flow through porous media in a single numerical framework, are usually employed to decrease the computational cost of full pore-scale simulations or to increase the accuracy of pure Darcy-scale simulations when a simple macroscopic description breaks down. Despite the massive increase in available computational power, the application of these techniques remains limited to core-size problems and upscaling remains crucial for practical large-scale applications. In this context, the Hybrid Multiscale Finite Volume (HMsFV) method, which constructs the macroscopic (Darcy-scale) problem directly by numerical averaging of pore-scale flow, offers not only a flexible framework to efficiently deal with multiphysics problems, but also a tool to investigate the assumptions used to derive macroscopic models and to better understand the relationship between pore-scale quantities and the corresponding macroscale variables. Indeed, by direct comparison of the multiphysics solution with a reference pore-scale simulation, we can assess the validity of the closure assumptions inherent to the multiphysics algorithm and infer the consequences for macroscopic models at the Darcy scale. We show that the definition of the scale ratio based on the geometric properties of the porous medium is well justified only for single-phase flow, whereas in case of unstable multiphase flow the nonlinear interplay between different forces creates complex fluid patterns characterized by new spatial scales, which emerge dynamically and weaken the scale-separation assumption. In general, the multiphysics solution proves very robust even when the characteristic size of the fluid-distribution patterns is comparable with the observation length, provided that all relevant physical processes affecting the fluid distribution are considered. This suggests that macroscopic constitutive relationships (e.g., the relative permeability) should account for the fact that they depend not only on the saturation but also on the actual characteristics of the fluid distribution.

Topological characteristics underpin intermittency and anomalous transport behavior in soil-like porous media

NASA Astrophysics Data System (ADS)

Holzner, M.; Morales, V.; Willmann, M.; Jerjen, I.; Kaufmann, R.; Dentz, M.

2016-12-01

Continuum models of porous media are based on the validity of the Darcy equation for fluid and Fick's law for scalar fluxes on a representative elementary volume. Fluctuations of pore-scale flow and scalar transport are averaged out and represented in terms of effective parameters such as hydrodynamic dispersion. However, the intermittent behavior of pore-scale flow impacts on the nature of particle and scalar transport, and it determines the way dissolved substances mix and react. The understanding of the origin of these processes is of both fundamental and practical importance in applications ranging from reactive transport in groundwater flow to diffusion in fuel cells or biological systems. A central issue in porous medium flow is therefore to relate intermittent behavior of Lagrangian velocity at pore scale imposed by the complex pore network geometry to transport properties at larger scales. Lagrangian measurements in porous systems are nonetheless scarce and most experimental techniques do not provide access to all three velocity components. In this contribution we report 3D measurements of Lagrangian velocity in soil-like porous media. We complement these measurements with detailed X-ray scans of the pore network. We find sharp velocity transitions close to pore throats, and low flow variability in the pore bodies, which gives rise to stretched exponential Lagrangian velocity and acceleration distributions characterized by a sharp peak at low velocity and a superlinear evolution of particle dispersion. We demonstrate that porosity and pore size distribution alone cannot explain the observed features of the flow. Rather, anomalous transport is better interpreted in terms of how pores of various geometries are interconnected. We reproduce the main observations using a continuous-time random walk (CTRW) model revealing the main features that control the system and showing the potential of this simple model to capture transport in complex geometries.

Visualizing and Quantifying Pore Scale Fluid Flow Processes With X-ray Microtomography

NASA Astrophysics Data System (ADS)

Wildenschild, D.; Hopmans, J. W.; Vaz, C. M.; Rivers, M. L.

2001-05-01

When using mathematical models based on Darcy's law it is often necessary to simplify geometry, physics or both and the capillary bundle-of-tubes approach neglects a fundamentally important characteristic of porous solids, namely interconnectedness of the pore space. New approaches to pore-scale modeling that arrange capillary tubes in two- or three-dimensional pore space have been and are still under development: Network models generally represent the pore space by spheres while the pore throats are usually represented by cylinders or conical shapes. Lattice Boltzmann approaches numerically solve the Navier-Stokes equations in a realistic microscopically disordered geometry, which offers the ability to study the microphysical basis of macroscopic flow without the need for a simplified geometry or physics. In addition to these developments in numerical modeling techniques, new theories have proposed that interfacial area should be considered as a primary variable in modeling of a multi-phase flow system. In the wake of this progress emerges an increasing need for new ways of evaluating pore-scale models, and for techniques that can resolve and quantify phase interfaces in porous media. The mechanisms operating at the pore-scale cannot be measured with traditional experimental techniques, however x-ray computerized microtomography (CMT) provides non-invasive observation of, for instance, changing fluid phase content and distribution on the pore scale. Interfacial areas have thus far been measured indirectly, but with the advances in high-resolution imaging using CMT it is possible to track interfacial area and curvature as a function of phase saturation or capillary pressure. We present results obtained at the synchrotron-based microtomography facility (GSECARS, sector 13) at the Advanced Photon Source at Argonne National Laboratory. Cylindrical sand samples of either 6 or 1.5 mm diameter were scanned at different stages of drainage and for varying boundary conditions. A significant difference in fluid saturation and phase distribution was observed for different drainage conditions, clearly showing preferential flow and a dependence on the applied flow rate. For the 1.5 mm sample individual pores and water/air interfaces could be resolved and quantified using image analysis techniques. Use of the Advanced Photon Source was supported by the U.S. Department of Energy, Basic Energy Sciences, Office of Science, under Contract No. W-31-109-Eng-38.

Understanding fluid transport through the multiscale pore network of a natural shale

NASA Astrophysics Data System (ADS)

Davy, Catherine; Adler, Pierre; Song, Yang; Nguyen, Thang Kim; Troadec, David; Dhénin, Jean-Francois

2017-04-01

Natural shales have a complex pore structure, which is only partly understood today. In the present contribution, a combination of different techniques is used to get information on three different scales. On each scale, the relevant flow equation is solved and provides input for the flow equation of the next higher scale. More precisely, micro-CT, FIB/SEM (Focused Ion Beam/Scanning Electron Microscopy) and TEM (Transmission Electron Microscopy) provide a full representative 3D pore space on the macroscopic scale, the mesoscale and the nanoscale. The corresponding typical voxel sizes are 0.7 μm, 10 nm and 1 nm, respectively. The porosity on the micro-CT images is 0.5 %, and it is not connected. One can distinguish between the pores, the porous clay matrix and non porous minerals; the volume percentages of these last two phases are 0.6 and 0.395, respectively. Samples of the porous clay matrix were analyzed by FIB/SEM which yields 3D information. They have a porosity ranging from 2 to 6 %. In some of them, the pore space is connected. Finally, TEM provides 2D images with a porosity of about 10 to 25 %. These information were used in the following way to estimate the macroscopic permeability which has been measured independently and found equal to 6 x10-20 m2. At the nanoscopic scale analyzed by 2D TEM, in the absence of 3D images, the pore structure is reconstructed by using a technique based on truncated Gaussian fields. Then, the Stokes equations are solved by using a 3D Lattice Boltzmann method. The resulting velocity field is averaged and this provides the permeability K_n. The permeability of the nanoscale structure varies between 0.7x 10-20 and 1.8x10-19 m2. As expected, the material is anisotropic. At the mesoscale, percolation of the FIB/SEM pore volume occurs only along a single direction. The Stokes equations are again solved by the same method and the mesoscopic permeability Km varies between 3.3 10-20 and 1.20 10-18 m2, depending on the nature of the percolating volume. The influence of the nanoscale porosity on the mesoscopic permeability is also studied. Two examples show that despite the scale ratio between the mesoscopic and nanoscopic pores, the nanoscopic pore structure cannot be neglected to estimate the permeability of the pore clay matrix. Finally, the sample provided by micro-CT is considered as a porous medium composed of three phases with permeabilities 0 (for the non porous minerals), 1 (for the porous clay matrix) and infinity (for the macroscopic pores). The overall permeability Kmacro is obtained by solving the Darcy's equation with a variable local permeability with spatially periodic boundary conditions. Kmacro is found of the order of 0.4 and the medium is relatively isotropic on this scale. This estimation of Kmacro is in agreement with the measured value.

Finite-difference method Stokes solver (FDMSS) for 3D pore geometries: Software development, validation and case studies

NASA Astrophysics Data System (ADS)

Gerke, Kirill M.; Vasilyev, Roman V.; Khirevich, Siarhei; Collins, Daniel; Karsanina, Marina V.; Sizonenko, Timofey O.; Korost, Dmitry V.; Lamontagne, Sébastien; Mallants, Dirk

2018-05-01

Permeability is one of the fundamental properties of porous media and is required for large-scale Darcian fluid flow and mass transport models. Whilst permeability can be measured directly at a range of scales, there are increasing opportunities to evaluate permeability from pore-scale fluid flow simulations. We introduce the free software Finite-Difference Method Stokes Solver (FDMSS) that solves Stokes equation using a finite-difference method (FDM) directly on voxelized 3D pore geometries (i.e. without meshing). Based on explicit convergence studies, validation on sphere packings with analytically known permeabilities, and comparison against lattice-Boltzmann and other published FDM studies, we conclude that FDMSS provides a computationally efficient and accurate basis for single-phase pore-scale flow simulations. By implementing an efficient parallelization and code optimization scheme, permeability inferences can now be made from 3D images of up to 109 voxels using modern desktop computers. Case studies demonstrate the broad applicability of the FDMSS software for both natural and artificial porous media.

Hydraulic conductivity of variably saturated porous media: Film and corner flow in angular pore space

NASA Astrophysics Data System (ADS)

Tuller, Markus; Or, Dani

2001-05-01

Many models for hydraulic conductivity of partially saturated porous media rely on oversimplified representation of the pore space as a bundle of cylindrical capillaries and disregard flow in liquid films. Recent progress in modeling liquid behavior in angular pores of partially saturated porous media offers an alternative framework. We assume that equilibrium liquid-vapor interfaces provide well-defined and stable boundaries for slow laminar film and corner flow regimes in pore space comprised of angular pores connected to slit-shaped spaces. Knowledge of liquid configuration in the assumed geometry facilitates calculation of average liquid velocities in films and corners and enables derivation of pore-scale hydraulic conductivity as a function of matric potential. The pore-scale model is statistically upscaled to represent hydraulic conductivity for a sample of porous medium. Model parameters for the analytical sample-scale expressions are estimated from measured liquid retention data and other measurable medium properties. Model calculations illustrate the important role of film flow, whose contribution dominates capillary flow (in full pores and corners) at relatively high matric potentials (approximately -100 to -300 J kg-1, or -1 to 3 bars). The crossover region between film and capillary flow is marked by a significant change in the slope of the hydraulic conductivity function as often observed in measurements. Model predictions are compared with the widely applied van Genuchten-Mualem model and yield reasonable agreement with measured retention and hydraulic conductivity data over a wide range of soil textural classes.

One Hundred Ways to be Non-Fickian - A Rigorous Multi-Variate Statistical Analysis of Pore-Scale Transport

NASA Astrophysics Data System (ADS)

Most, Sebastian; Nowak, Wolfgang; Bijeljic, Branko

2015-04-01

Fickian transport in groundwater flow is the exception rather than the rule. Transport in porous media is frequently simulated via particle methods (i.e. particle tracking random walk (PTRW) or continuous time random walk (CTRW)). These methods formulate transport as a stochastic process of particle position increments. At the pore scale, geometry and micro-heterogeneities prohibit the commonly made assumption of independent and normally distributed increments to represent dispersion. Many recent particle methods seek to loosen this assumption. Hence, it is important to get a better understanding of the processes at pore scale. For our analysis we track the positions of 10.000 particles migrating through the pore space over time. The data we use come from micro CT scans of a homogeneous sandstone and encompass about 10 grain sizes. Based on those images we discretize the pore structure and simulate flow at the pore scale based on the Navier-Stokes equation. This flow field realistically describes flow inside the pore space and we do not need to add artificial dispersion during the transport simulation. Next, we use particle tracking random walk and simulate pore-scale transport. Finally, we use the obtained particle trajectories to do a multivariate statistical analysis of the particle motion at the pore scale. Our analysis is based on copulas. Every multivariate joint distribution is a combination of its univariate marginal distributions. The copula represents the dependence structure of those univariate marginals and is therefore useful to observe correlation and non-Gaussian interactions (i.e. non-Fickian transport). The first goal of this analysis is to better understand the validity regions of commonly made assumptions. We are investigating three different transport distances: 1) The distance where the statistical dependence between particle increments can be modelled as an order-one Markov process. This would be the Markovian distance for the process, where the validity of yet-unexplored non-Gaussian-but-Markovian random walks start. 2) The distance where bivariate statistical dependence simplifies to a multi-Gaussian dependence based on simple linear correlation (validity of correlated PTRW/CTRW). 3) The distance of complete statistical independence (validity of classical PTRW/CTRW). The second objective is to reveal characteristic dependencies influencing transport the most. Those dependencies can be very complex. Copulas are highly capable of representing linear dependence as well as non-linear dependence. With that tool we are able to detect persistent characteristics dominating transport even across different scales. The results derived from our experimental data set suggest that there are many more non-Fickian aspects of pore-scale transport than the univariate statistics of longitudinal displacements. Non-Fickianity can also be found in transverse displacements, and in the relations between increments at different time steps. Also, the found dependence is non-linear (i.e. beyond simple correlation) and persists over long distances. Thus, our results strongly support the further refinement of techniques like correlated PTRW or correlated CTRW towards non-linear statistical relations.

Experimental evidence of the role of pores on movement and distribution of bacteria in soil

NASA Astrophysics Data System (ADS)

Kravchenko, Alexandra N.; Rose, Joan B.; Marsh, Terence L.; Guber, Andrey K.

2014-05-01

It has been generally recognized that micro-scale heterogeneity in soil environments can have a substantial effect on movement, fate, and survival of soil microorganisms. However, only recently the development of tools for micro-scale soil analyses, including X-ray computed micro-tomography (μ-CT), enabled quantitative analyses of these effects. The long-term goal of our work is to explore how differences in micro-scale characteristics of pore structures influence movement, spatial distribution patterns, and activities of soil microorganisms. Using X-ray μ-CT we found that differences in land use and management practices lead to development of contrasting patterns in pore size-distributions within intact soil aggregates. Then our experiments with Escherichia coli added to intact soil aggregates demonstrated that the differences in pore structures can lead to substantial differences in bacteria redistribution and movement within the aggregates. Specifically, we observed more uniform E.coli redistribution in aggregates with homogeneously spread pores, while heterogeneous pore structures resulted in heterogeneous E.coli patterns. Water flow driven by capillary forces through intact aggregate pores appeared to be the main contributor to the movement patterns of the introduced bacteria. Influence of pore structure on E.coli distribution within the aggregates further continued after the aggregates were subjected to saturated water flow. E. coli's resumed movement with saturated water flow and subsequent redistribution within the soil matrix was influenced by porosity, abundance of medium and large pores, pore tortuosity, and flow rates, indicating that greater flow accompanied by less convoluted pores facilitated E. coli transport within the intra-aggregate space. We also found that intra-aggregate heterogeneity of pore structures can have an effect on spatial distribution patterns of indigenous microbial populations. Preliminary analysis showed that in aggregates from an organic agricultural system with cover crops, characterized by greater intra-aggregate pore heterogeneity, bacteria of Actinobacteria and Firmicutes groups were more abundant in presence of large as compared to small pores. In contrast, no differences were observed in the aggregates from conventionally managed soil, overall characterized by homogeneous intra-aggregate pore patterns. Further research efforts are being directed towards quantification of the pore structure effects on activities and community composition of soil microorganisms.

Determination of pore-scale hydrate phase equilibria in sediments using lab-on-a-chip technology.

PubMed

Almenningen, Stian; Flatlandsmo, Josef; Kovscek, Anthony R; Ersland, Geir; Fernø, Martin A

2017-11-21

We present an experimental protocol for fast determination of hydrate stability in porous media for a range of pressure and temperature (P, T) conditions. Using a lab-on-a-chip approach, we gain direct optical access to dynamic pore-scale hydrate formation and dissociation events to study the hydrate phase equilibria in sediments. Optical pore-scale observations of phase behavior reproduce the theoretical hydrate stability line with methane gas and distilled water, and demonstrate the accuracy of the new method. The procedure is applicable for any kind of hydrate transitions in sediments, and may be used to map gas hydrate stability zones in nature.

EVALUATION AND ANALYSIS OF MICROSCALE FLOW AND TRANSPORT DURING REMEDIATION

EPA Science Inventory

The design of in-situ remediation is currently based on a description at the macroscopic scale. Phenomena at the pore and pore-network scales are typically lumped in terms of averaged quantities, using empirical or ad hoc expressions. These models cannot address fundamental rem...

Imaging of forced-imbibition in carbonate rocks using synchrotron X-ray micro-tomography

NASA Astrophysics Data System (ADS)

Singh, K.; Menke, H. P.; Andrew, M. G.; Lin, Q.; Saif, T.; Al-Khulaifi, Y.; Reynolds, C. A.; Bijeljic, B.; Rau, C.; Blunt, M. J.

2016-12-01

We have investigated the pore-scale behavior of brine-oil systems and oil trapping during forced-imbibition in a water-wet carbonate rock in a capillary-dominated flow regime at reservoir pressure conditions. To capture the dynamics of the brine-oil front progression and snap-off process, real-time tomograms with a time resolution of 38 s (24 s for imaging and 14 s for recording the data) and a spatial resolution of 3.28 µm were acquired at Diamond Light Source (UK). The data were first analyzed at global scale (complete imaged rock) for overall front behavior. From the saturation profiles, we obtain the location of the tail of the desaturation front that progresses with a velocity of 13 µm/min. This velocity is smaller than average flow velocity 16.88 µm/min, which explains why it needs slightly more than 1 pore volume of brine injection to reach the residual saturation of oil in a water-wet rock. The data were further analyzed at local scale to investigate the pore-scale mechanisms of oil trapping during brine flooding. We isolated various trapping events which resulted in the creation of discrete oil ganglia occupying one to several pore bodies. We perform pore-scale curvature analysis of brine-oil interfaces to obtain local capillary pressure that will be related to the shape and the size of throats in which ganglia were trapped.

A mathematical multiscale model of bone remodeling, accounting for pore space-specific mechanosensation.

PubMed

Pastrama, Maria-Ioana; Scheiner, Stefan; Pivonka, Peter; Hellmich, Christian

2018-02-01

While bone tissue is a hierarchically organized material, mathematical formulations of bone remodeling are often defined on the level of a millimeter-sized representative volume element (RVE), "smeared" over all types of bone microstructures seen at lower observation scales. Thus, there is no explicit consideration of the fact that the biological cells and biochemical factors driving bone remodeling are actually located in differently sized pore spaces: active osteoblasts and osteoclasts can be found in the vascular pores, whereas the lacunar pores host osteocytes - bone cells originating from former osteoblasts which were then "buried" in newly deposited extracellular bone matrix. We here propose a mathematical description which considers size and shape of the pore spaces where the biological and biochemical events take place. In particular, a previously published systems biology formulation, accounting for biochemical regulatory mechanisms such as the rank-rankl-opg pathway, is cast into a multiscale framework coupled to a poromicromechanical model. The latter gives access to the vascular and lacunar pore pressures arising from macroscopic loading. Extensive experimental data on the biological consequences of this loading strongly suggest that the aforementioned pore pressures, together with the loading frequency, are essential drivers of bone remodeling. The novel approach presented here allows for satisfactory simulation of the evolution of bone tissue under various loading conditions, and for different species; including scenarios such as mechanical dis- and overuse of murine and human bone, or in osteocyte-free bone. Copyright © 2017 Elsevier Inc. All rights reserved.

Non-scaling behavior of electroosmotic flow in voltage-gated nanopores

DOE PAGES

Lian, Cheng; Gallegos, Alejandro; Liu, Honglai; ...

2016-11-17

Ionic transport through nanopores is of fundamental importance for the design and development of nanofiltration membranes and novel electrochemical devices including supercapacitors, fuel cells and batteries. Recent experiments have shown an unusual variation of electrical conductance with the pore size and the electrolyte parameters that defies conventional scaling relations. Here ionic transport through voltage-gated nanopores was studied by using the classical density functional theory for ion distributions in combination with the Navier–Stokes equation for the electroosmotic flow. We also identified a significant influence of the gating potential on the scaling behavior of the conductance with changes in the pore sizemore » and the salt concentration. Finally, for ion transport in narrow pores with a high gating voltage, the conductivity shows an oscillatory dependence on the pore size owing to the strong overlap of electric double layers.« less

Reversible Self-Actuated Thermo-Responsive Pore Membrane

PubMed Central

Park, Younggeun; Gutierrez, Maria Paz; Lee, Luke P.

2016-01-01

Smart membranes, which can selectively control the transfer of light, air, humidity and temperature, are important to achieve indoor climate regulation. Even though reversible self-actuation of smart membranes is desirable in large-scale, reversible self-regulation remains challenging. Specifically, reversible 100% opening/closing of pore actuation showing accurate responsiveness, reproducibility and structural flexibility, including uniform structure assembly, is currently very difficult. Here, we report a reversible, thermo-responsive self-activated pore membrane that achieves opening and closing of pores. The reversible, self-actuated thermo-responsive pore membrane was fabricated with hybrid materials of poly (N-isopropylacrylamide), (PNIPAM) within polytetrafluoroethylene (PTFE) to form a multi-dimensional pore array. Using Multiphysics simulation of heat transfer and structural mechanics based on finite element analysis, we demonstrated that pore opening and closing dynamics can be self-activated at environmentally relevant temperatures. Temperature cycle characterizations of the pore structure revealed 100% opening ratio at T = 40 °C and 0% opening ratio at T = 20 °C. The flexibility of the membrane showed an accurate temperature-responsive function at a maximum bending angle of 45°. Addressing the importance of self-regulation, this reversible self-actuated thermo-responsive pore membrane will advance the development of future large-scale smart membranes needed for sustainable indoor climate control. PMID:27991563

Reversible Self-Actuated Thermo-Responsive Pore Membrane

NASA Astrophysics Data System (ADS)

Park, Younggeun; Gutierrez, Maria Paz; Lee, Luke P.

2016-12-01

Smart membranes, which can selectively control the transfer of light, air, humidity and temperature, are important to achieve indoor climate regulation. Even though reversible self-actuation of smart membranes is desirable in large-scale, reversible self-regulation remains challenging. Specifically, reversible 100% opening/closing of pore actuation showing accurate responsiveness, reproducibility and structural flexibility, including uniform structure assembly, is currently very difficult. Here, we report a reversible, thermo-responsive self-activated pore membrane that achieves opening and closing of pores. The reversible, self-actuated thermo-responsive pore membrane was fabricated with hybrid materials of poly (N-isopropylacrylamide), (PNIPAM) within polytetrafluoroethylene (PTFE) to form a multi-dimensional pore array. Using Multiphysics simulation of heat transfer and structural mechanics based on finite element analysis, we demonstrated that pore opening and closing dynamics can be self-activated at environmentally relevant temperatures. Temperature cycle characterizations of the pore structure revealed 100% opening ratio at T = 40 °C and 0% opening ratio at T = 20 °C. The flexibility of the membrane showed an accurate temperature-responsive function at a maximum bending angle of 45°. Addressing the importance of self-regulation, this reversible self-actuated thermo-responsive pore membrane will advance the development of future large-scale smart membranes needed for sustainable indoor climate control.

Tunable acoustic absorbers with periodical micro-perforations having varying pore shapes

NASA Astrophysics Data System (ADS)

Ren, Shuwei; Liu, Xuewei; Gong, Junqing; Tang, Yufan; Xin, Fengxian; Huang, Lixi; Lu, Tian Jian

2017-11-01

Circular pores with sub-millimeter diameters have been widely used to construct micro-perforated panels (MPPs), the acoustical performance of which can be predicted well using the Maa theory (MAA D.-Y., J. Acoust. Soc. Am., 104 (1998) 2861). We present a tunable MPP absorber with periodically arranged cylindrical pores, with their cross-sectional shapes systematically altered around the circle while maintaining their cross-sectional areas unchanged. Numerical analyses based on the viscous-thermal coupled acoustical equations are utilized to investigate the tunable acoustic performance of the proposed absorbers and to reveal the underlying physical mechanisms. We demonstrate that pore morphology significantly affects the sound absorbption of MPPs by modifying the velocity field (and hence viscous dissipation) in the pores. Pore shapes featured as meso-scale circular pores accompanied with micro-scale bulges along the boundaries can lead to perfect sound absorption at relatively low frequencies. This work not only enriches the classical Maa theory on MPPs having circular perforations, but it also opens a new avenue for designing subwavelength acoustic metamaterials of superior sound absorption in target frequency ranges.

Pore-scale Modeling of CO2 Local Trapping in Heterogeneous Porous Media with Inter-granular Cements

NASA Astrophysics Data System (ADS)

Wang, D.; Li, Y.

2017-12-01

Based on pore-scale modeling of CO2/brine multiphase flow in heterogeneous porous media with inter-granular cements, we numerically analyze the effects of cement-modified pore structure on CO2 local trapping. Results indicate: 1) small pore throat is the main reason for causing CO2 local trapping in front of low-porosity layers (namely dense layers) formed by inter-granular cements; 2) in the case of the same pore throat size, the smaller particle size can increase the number of flow paths for CO2 plume and equivalently enhances local permeability, which may counteract the impediment of high capillary pressure on CO2 migration to some extent and consequently disables CO2 local capillary trapping; 3) the isolated pores by inter-granular cements can lead to dramatic reduction of CO2 saturation inside the dense layers, whereas the change of connectivity of some pores due to the cements can increase CO2 accumulation in front of the dense layers by lowering the displacement area of CO2 plume.

Dynamic three-dimensional pore-scale imaging of reaction in a carbonate at reservoir conditions.

PubMed

Menke, Hannah P; Bijeljic, Branko; Andrew, Matthew G; Blunt, Martin J

2015-04-07

Quantifying CO2 transport and average effective reaction rates in the subsurface is essential to assess the risks associated with underground carbon capture and storage. We use X-ray microtomography to investigate dynamic pore structure evolution in situ at temperatures and pressures representative of underground reservoirs and aquifers. A 4 mm diameter Ketton carbonate core is injected with CO2-saturated brine at 50 °C and 10 MPa while tomographic images are taken at 15 min intervals with a 3.8 μm spatial resolution over a period of 2(1/2) h. An approximate doubling of porosity with only a 3.6% increase in surface area to volume ratio is measured from the images. Pore-scale direct simulation and network modeling on the images quantify an order of magnitude increase in permeability and an appreciable alteration of the velocity field. We study the uniform reaction regime, with dissolution throughout the core. However, at the pore scale, we see variations in the degree of dissolution with an overall reaction rate which is approximately 14 times lower than estimated from batch measurements. This work implies that in heterogeneous rocks, pore-scale transport of reactants limits dissolution and can reduce the average effective reaction rate by an order of magnitude.

Polarity control at interfaces: Quantifying pseudo-solvent effects in nano-confined systems

DOE PAGES

Singappuli-Arachchige, Dilini; Manzano, J. Sebastian; Sherman, Lindy M.; ...

2016-08-02

Surface functionalization controls local environments and induces solvent-like effects at liquid–solid interfaces. We explored structure–property relationships between organic groups bound to pore surfaces of mesoporous silica nanoparticles and Stokes shifts of the adsorbed solvatochromic dye Prodan. Correlating shifts of the dye on the surfaces with its shifts in solvents resulted in a local polarity scale for functionalized pores. The scale was validated by studying the effects of pore polarity on quenching of Nile Red fluorescence and on the vibronic band structure of pyrene. Measurements were done in aqueous suspensions of porous particles, proving that the dielectric properties in the poresmore » are different from the bulk solvent. The precise control of pore polarity was used to enhance the catalytic activity of TEMPO in the aerobic oxidation of furfuryl alcohol in water. Furthermore, an inverse relationship was found between pore polarity and activity of TEMPO in the pores, demonstrating that controlling the local polarity around an active site allows modulating the activity of nanoconfined catalysts.« less

Condensation in Nanoporous Packed Beds.

PubMed

Ally, Javed; Molla, Shahnawaz; Mostowfi, Farshid

2016-05-10

In materials with tiny, nanometer-scale pores, liquid condensation is shifted from the bulk saturation pressure observed at larger scales. This effect is called capillary condensation and can block pores, which has major consequences in hydrocarbon production, as well as in fuel cells, catalysis, and powder adhesion. In this study, high pressure nanofluidic condensation studies are performed using propane and carbon dioxide in a colloidal crystal packed bed. Direct visualization allows the extent of condensation to be observed, as well as inference of the pore geometry from Bragg diffraction. We show experimentally that capillary condensation depends on pore geometry and wettability because these factors determine the shape of the menisci that coalesce when pore filling occurs, contrary to the typical assumption that all pore structures can be modeled as cylindrical and perfectly wetting. We also observe capillary condensation at higher pressures than has been done previously, which is important because many applications involving this phenomenon occur well above atmospheric pressure, and there is little, if any, experimental validation of capillary condensation at such pressures, particularly with direct visualization.

Comparison of Pore-Network and Lattice Boltzmann Models for Pore-Scale Modeling of Geological Storage of CO2 in Natural Reservoir Rocks

NASA Astrophysics Data System (ADS)

Kohanpur, A. H.; Chen, Y.; Valocchi, A. J.; Tudek, J.; Crandall, D.

2016-12-01

CO2-brine flow in deep natural rocks is the focus of attention in geological storage of CO2. Understanding rock/flow properties at pore-scale is a vital component in field-scale modeling and prediction of fate of injected CO2. There are many challenges in working at the pore scale, such as size and selection of representative elementary volume (REV), particularly for material with complex geometry and heterogeneity, and the high computational costs. These issues factor into trade-offs that need to be made in choosing and applying pore-scale models. On one hand, pore-network modeling (PNM) simplifies the geometry and flow equations but can provide characteristic curves on fairly large samples. On the other hand, the lattice Boltzmann method (LBM) solves Navier-Stokes equations on the real geometry but is limited to small samples due to its high computational costs. Thus, both methods have some advantages but also face some challenges, which warrants a more detailed comparison and evaluation. In this study, we used industrial and micro-CT scans of actual reservoir rock samples to characterize pore structure at different resolutions. We ran LBM models directly on the characterized geometry and PNM on the equivalent 3D extracted network to determine single/two-phase flow properties during drainage and imbibition processes. Specifically, connectivity, absolute permeability, relative permeability curve, capillary pressure curve, and interface location are compared between models. We also did simulations on several subsamples from different locations including different domain sizes and orientations to encompass analysis of heterogeneity and isotropy. This work is primarily supported as part of the Center for Geologic Storage of CO2, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science and partially supported by the International Institute for Carbon-Neutral Energy Research (WPI-I2CNER) based at Kyushu University, Japan.

Comparison of Pore-scale CO2-water-glass System Wettability and Conventional Wettability Measurement on a Flat Plate for Geological CO2 Sequestration

NASA Astrophysics Data System (ADS)

Jafari, M.; Cao, S. C.; Jung, J.

2017-12-01

Goelogical CO2 sequestration (GCS) has been recently introduced as an effective method to mitigate carbon dioxide emission. CO2 from main producer sources is collected and then is injected underground formations layers to be stored for thousands to millions years. A safe and economical storage project depends on having an insight of trapping mechanisms, fluids dynamics, and interaction of fluids-rocks. Among different forces governing fluids mobility and distribution in GCS condition, capillary pressure is of importance, which, in turn, wettability (measured by contact angel (CA)) is the most controversial parameters affecting it. To explore the sources of discrepancy in the literature for CA measurement, we conducted a series of conventional captive bubble test on glass plates under high pressure condition. By introducing a shape factor, we concluded that surface imperfection can distort the results in such tests. Since the conventional methods of measuring the CA is affected by gravity and scale effect, we introduced a different technique to measure pore-scale CA inside a transparent glass microchip. Our method has the ability to consider pore sizes and simulate static and dynamics CA during dewetting and imbibition. Glass plates shows a water-wet behavior (CA 30° - 45°) by a conventional experiment consistent with literature. However, CA of miniature bubbles inside of the micromodel can have a weaker water-wet behavior (CA 55° - 69°). In a more realistic pore-scale condition, water- CO2 interface covers whole width of a pore throats. Under this condition, the receding CA, which is used for injectability and capillary breakthrough pressure, increases with decreasing pores size. On the other hand, advancing CA, which is important for residual or capillary trapping, does not show a correlation with throat sizes. Static CA measured in the pores during dewetting is lower than static CA on flat plate, but it is much higher when measured during imbibition implying weaker water-wet behavior. Pore-scale CA, which realistically represents rocks wettability behavior, shows weaker water-wet behavior than conventional measurement methods, which must be considered for safety of geological storage.

Importance of filter’s microstructure in dynamic filtration modeling of gasoline particulate filters (GPFs): Inhomogeneous porosity and pore size distribution

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gong, Jian; Stewart, Mark L.; Zelenyuk, Alla

The state-of-the-art multiscale modeling of GPFs including channel scale, wall scale, and pore scale is described. The microstructures of two GPFs were experimentally characterized. The pore size distributions of the GPFs were determined by mercury porosimetry. The porosity was measured by X-ray computed tomography (CT) and found to be inhomogeneous across the substrate wall. The significance of pore size distribution with respect to filtration performance was analyzed. The predictions of filtration efficiency were improved by including the pore size distribution in the filtration model. A dynamic heterogeneous multiscale filtration (HMF) model was utilized to simulate particulate filtration on a singlemore » channel particulate filter with realistic particulate emissions from a spark-ignition direct-injection (SIDI) gasoline engine. The dynamic evolution of filter’s microstructure and macroscopic filtration characteristics including mass- and number-based filtration efficiencies and pressure drop were predicted and discussed. The microstructure of the GPF substrate including inhomogeneous porosity and pore size distribution is found to significantly influence local particulate deposition inside the substrate and macroscopic filtration performance and is recommended to be resolved in the filtration model to simulate and evaluate the filtration performance of GPFs.« less

Importance of filter’s microstructure in dynamic filtration modeling of gasoline particulate filters (GPFs): Inhomogeneous porosity and pore size distribution

DOE PAGES

Gong, Jian; Stewart, Mark L.; Zelenyuk, Alla; ...

2018-01-03

The state-of-the-art multiscale modeling of gasoline particulate filter (GPF) including channel scale, wall scale, and pore scale is described. The microstructures of two GPFs were experimentally characterized. The pore size distributions of the GPFs were determined by mercury porosimetry. The porosity was measured by X-ray computed tomography (CT) and found to be inhomogeneous across the substrate wall. The significance of pore size distribution with respect to filtration performance was analyzed. The predictions of filtration efficiency were improved by including the pore size distribution in the filtration model. A dynamic heterogeneous multiscale filtration (HMF) model was utilized to simulate particulate filtrationmore » on a single channel particulate filter with realistic particulate emissions from a spark-ignition direct-injection (SIDI) gasoline engine. The dynamic evolution of filter’s microstructure and macroscopic filtration characteristics including mass- and number-based filtration efficiencies and pressure drop were predicted and discussed. In conclusion, the microstructure of the GPF substrate including inhomogeneous porosity and pore size distribution is found to significantly influence local particulate deposition inside the substrate and macroscopic filtration performance and is recommended to be resolved in the filtration model to simulate and evaluate the filtration performance of GPFs.« less

Importance of filter’s microstructure in dynamic filtration modeling of gasoline particulate filters (GPFs): Inhomogeneous porosity and pore size distribution

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gong, Jian; Stewart, Mark L.; Zelenyuk, Alla

The state-of-the-art multiscale modeling of gasoline particulate filter (GPF) including channel scale, wall scale, and pore scale is described. The microstructures of two GPFs were experimentally characterized. The pore size distributions of the GPFs were determined by mercury porosimetry. The porosity was measured by X-ray computed tomography (CT) and found to be inhomogeneous across the substrate wall. The significance of pore size distribution with respect to filtration performance was analyzed. The predictions of filtration efficiency were improved by including the pore size distribution in the filtration model. A dynamic heterogeneous multiscale filtration (HMF) model was utilized to simulate particulate filtrationmore » on a single channel particulate filter with realistic particulate emissions from a spark-ignition direct-injection (SIDI) gasoline engine. The dynamic evolution of filter’s microstructure and macroscopic filtration characteristics including mass- and number-based filtration efficiencies and pressure drop were predicted and discussed. In conclusion, the microstructure of the GPF substrate including inhomogeneous porosity and pore size distribution is found to significantly influence local particulate deposition inside the substrate and macroscopic filtration performance and is recommended to be resolved in the filtration model to simulate and evaluate the filtration performance of GPFs.« less

High Fidelity Computational Analysis of CO2 Trapping at Pore Scales

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kumar, Vinod

2013-07-13

With an alarming rise in carbon dioxide (CO2) emission from anthropogenic sources, CO2 sequestration has become an attractive choice to mitigate the emission. Some popular storage media for CO{sub 2} are oil reservoirs, deep coal-bed, and deep oceanic-beds. These have been used for the long term CO{sub 2} storage. Due to special lowering viscosity and surface tension property of CO{sub 2}, it has been widely used for enhanced oil recovery. The sites for CO{sub 2} sequestration or enhanced oil recovery mostly consist of porous rocks. Lack of knowledge of molecular mobility under confinement and molecule-surface interactions between CO2 and naturalmore » porous media results in generally governed by unpredictable absorption kinetics and total absorption capacity for injected fluids, and therefore, constitutes barriers to the deployment of this technology. Therefore, it is important to understand the flow dynamics of CO{sub 2} through the porous microstructures at the finest scale (pore-scale) to accurately predict the storage potential and long-term dynamics of the sequestered CO{sub 2}. This report discusses about pore-network flow modeling approach using variational method and analyzes simulated results this method simulations at pore-scales for idealized network and using Berea Sandstone CT scanned images. Variational method provides a promising way to study the kinetic behavior and storage potential at the pore scale in the presence of other phases. The current study validates variational solutions for single and two-phase Newtonian and single phase non-Newtonian flow through angular pores for special geometries whose analytical and/or empirical solutions are known. The hydraulic conductance for single phase flow through a triangular duct was also validated against empirical results derived from lubricant theory.« less

Connecting the molecular scale to the continuum scale for diffusion processes in smectite-rich porous media.

PubMed

Bourg, Ian C; Sposito, Garrison

2010-03-15

In this paper, we address the manner in which the continuum-scale diffusive properties of smectite-rich porous media arise from their molecular- and pore-scale features. Our starting point is a successful model of the continuum-scale apparent diffusion coefficient for water tracers and cations, which decomposes it as a sum of pore-scale terms describing diffusion in macropore and interlayer "compartments." We then apply molecular dynamics (MD) simulations to determine molecular-scale diffusion coefficients D(interlayer) of water tracers and representative cations (Na(+), Cs(+), Sr(2+)) in Na-smectite interlayers. We find that a remarkably simple expression relates D(interlayer) to the pore-scale parameter δ(nanopore) ≤ 1, a constrictivity factor that accounts for the lower mobility in interlayers as compared to macropores: δ(nanopore) = D(interlayer)/D(0), where D(0) is the diffusion coefficient in bulk liquid water. Using this scaling expression, we can accurately predict the apparent diffusion coefficients of tracers H(2)0, Na(+), Sr(2+), and Cs(+) in compacted Na-smectite-rich materials.

Experimental Quantification of Pore-Scale Flow Phenomena in 2D Heterogeneous Porous Micromodels: Multiphase Flow Towards Coupled Solid-Liquid Interactions

NASA Astrophysics Data System (ADS)

Li, Y.; Kazemifar, F.; Blois, G.; Christensen, K. T.

2017-12-01

Geological sequestration of CO2 within saline aquifers is a viable technology for reducing CO2 emissions. Central to this goal is accurately predicting both the fidelity of candidate sites pre-injection of CO2 and its post-injection migration. Moreover, local fluid pressure buildup may cause activation of small pre-existing unidentified faults, leading to micro-seismic events, which could prove disastrous for societal acceptance of CCS, and possibly compromise seal integrity. Recent evidence shows that large-scale events are coupled with pore-scale phenomena, which necessitates the representation of pore-scale stress, strain, and multiphase flow processes in large-scale modeling. To this end, the pore-scale flow of water and liquid/supercritical CO2 is investigated under reservoir-relevant conditions, over a range of wettability conditions in 2D heterogeneous micromodels that reflect the complexity of a real sandstone. High-speed fluorescent microscopy, complemented by a fast differential pressure transmitter, allows for simultaneous measurement of the flow field within and the instantaneous pressure drop across the micromodels. A flexible micromodel is also designed and fabricated, to be used in conjunction with the micro-PIV technique, enabling the quantification of coupled solid-liquid interactions.

Quantitative high-resolution mapping of phenanthrene sorption to black carbon particles.

PubMed

Obst, Martin; Grathwohl, Peter; Kappler, Andreas; Eibl, Oliver; Peranio, Nicola; Gocht, Tilman

2011-09-01

Sorption of hydrophobic organic contaminants such as polycyclic aromatic hydrocarbons (PAHs) to black carbon (BC) particles has been the focus of numerous studies. Conclusions on sorption mechanisms of PAH on BC were mostly derived from studies of sorption isotherms and sorption kinetics, which are based on batch experiments. However, mechanistic modeling approaches consider processes at the subparticle scale, some including transport within the pore-space or different spatial pore-domains. Direct evidence based on analytical techniques operating at the submicrometer scale for the location of sorption sites and the adsorbed species is lacking. In this work, we identified, quantified, and mapped the sorption of PAHs on different BC particles (activated carbon, charcoal and diesel soot) on a 25-100 nm scale using scanning transmission X-ray microscopy (STXM). In addition, we visualized the pore structure of the particles by transmission electron microscopy (TEM) on the 1-10 nm-scale. The combination of the chemical information from STXM with the physical information from TEM revealed that phenanthrene accumulates in the interconnected pore-system along primary "cracks" in the particles, confirming an adsorption mechanism.

2D and 3D characterization of pore defects in die cast AM60

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yang, Zhuofei; CanmetMATERIALS, 183 Longwood Road South, Hamilton L8P 0A5, Ontario Canada; Maurey, Alexandre

2016-04-15

The widespread application of die castings can be hampered due to the potential of large scale porosity to act as nucleation sites for fracture and fatigue. It is therefore important to develop robust approaches to the characterization of porosity providing parameters that can be linked to the material's mechanical properties. We have tackled this problem in a study of the AM60 die cast Mg alloy, using samples extracted from a prototype shock tower. A quantitative characterization of porosity has been undertaken, analyzing porosity in both 2D (using classical metallographic methods) and in 3D (using X-ray computed tomography (XCT)). Metallographic characterizationmore » results show that shrinkage pores and small gas pores can be distinguished based on their distinct geometrical features. Shrinkage pores are irregular with multiple arms, resulting in a form factor less than 0.4. In contrast, gas pores are generally more circular in shape yielding form factors larger than 0.6. XCT provides deeper insight into the shape of pores, although this understanding is limited by the resolution obtainable by laboratory based XCT. It also shows how 2D sectioning can produce artefacts as single complex pores are sectioned into multiple small pores. - Highlights: • Mg (e.g. AM60) die castings may contain large scale porosity that act as nucleation sites for fracture and fatigue • Quantitative characterization of porosity metallography (2D) and X-ray tomography (3D) is used • Shrinkage pores and small gas pores can be distinguished based on their distinct geometrical features. • Shrinkage pores are irregular giving a form factor < 0.4; gas pores are rounder with form factors > 0.6 • XCT enables pore visualization, although limited by the resolution obtainable by laboratory based XCT.« less

Monitoring of stainless-steel slag carbonation using X-ray computed microtomography.

PubMed

Boone, Marijn A; Nielsen, Peter; De Kock, Tim; Boone, Matthieu N; Quaghebeur, Mieke; Cnudde, Veerle

2014-01-01

Steel production is one of the largest contributors to industrial CO2 emissions. This industry also generates large amounts of solid byproducts, such as slag and sludge. In this study, fine grained stainless-steel slag (SSS) is valorized to produce compacts with high compressive strength without the use of a hydraulic binder. This carbonation process is investigated on a pore-scale level to identify how the mineral phases in the SSS react with CO2, where carbonates are formed, and what the impact of these changes is on the pore network of the carbonated SSS compact. In addition to conventional research techniques, high-resolution X-ray computed tomography (HRXCT) is applied to visualize and quantify the changes in situ during the carbonation process. The results show that carbonates mainly precipitate at grain contacts and in capillary pores and this precipitation has little effect on the connectivity of the pore space. This paper also demonstrates the use of a custom-designed polymer reaction cell that allows in situ HRXCT analysis of the carbonation process. This shows the distribution and influence of water and CO2 in the pore network on the carbonate precipitation and, thus, the influence on the compressive strength development of the waste material.

Monitoring CO2 invasion processes at the pore scale using geological labs on chip.

PubMed

Morais, S; Liu, N; Diouf, A; Bernard, D; Lecoutre, C; Garrabos, Y; Marre, S

2016-09-21

In order to investigate at the pore scale the mechanisms involved during CO2 injection in a water saturated pore network, a series of displacement experiments is reported using high pressure micromodels (geological labs on chip - GLoCs) working under real geological conditions (25 < T (°C) < 75 and 4.5 < p (MPa) < 8). The experiments were focused on the influence of three experimental parameters: (i) the p, T conditions, (ii) the injection flow rates and (iii) the pore network characteristics. By using on-chip optical characterization and imaging approaches, the CO2 saturation curves as a function of either time or the number of pore volume injected were determined. Three main mechanisms were observed during CO2 injection, namely, invasion, percolation and drying, which are discussed in this paper. Interestingly, besides conventional mechanisms, two counterintuitive situations were observed during the invasion and drying processes.

Quantification of soil structure based on Minkowski functions

NASA Astrophysics Data System (ADS)

Vogel, H.-J.; Weller, U.; Schlüter, S.

2010-10-01

The structure of soils and other geologic media is a complex three-dimensional object. Most of the physical material properties including mechanical and hydraulic characteristics are immediately linked to the structure given by the pore space and its spatial distribution. It is an old dream and still a formidable challenge to relate structural features of porous media to their functional properties. Using tomographic techniques, soil structure can be directly observed at a range of spatial scales. In this paper we present a scale-invariant concept to quantify complex structures based on a limited set of meaningful morphological functions. They are based on d+1 Minkowski functionals as defined for d-dimensional bodies. These basic quantities are determined as a function of pore size or aggregate size obtained by filter procedures using mathematical morphology. The resulting Minkowski functions provide valuable information on the size of pores and aggregates, the pore surface area and the pore topology having the potential to be linked to physical properties. The theoretical background and the related algorithms are presented and the approach is demonstrated for the pore structure of an arable soil and the pore structure of a sand both obtained by X-ray micro-tomography. We also analyze the fundamental problem of limited resolution which is critical for any attempt to quantify structural features at any scale using samples of different size recorded at different resolutions. The results demonstrate that objects smaller than 5 voxels are critical for quantitative analysis.

Campaign 2 Level 2 Milestone Review 2009: Milestone # 3131 Grain Scale Simulation of Pore Collapse

DOE Office of Scientific and Technical Information (OSTI.GOV)

Schwartz, A J

2009-09-28

The milestone reviewed on Sept. 16, 2009 was 'High-fidelity simulation of shock initiation of high explosives at the grain scale using coupled hydrodynamics, thermal transport and chemistry'. It is the opinion of the committee that the team has satisfied the milestone. A detailed description of how the goals were met is provided. The milestone leveraged capabilities from ASC Physics and Engineering Materials program combined with experimental input from Campaign 2. A combined experimental-multiscale simulation approach was used to create and validate the various TATB model components. At the lowest length scale, quantum chemical calculations were used to determine equations ofmore » state, thermal transport properties and reaction rates for TATB as it is decomposing. High-pressure experiments conducted in diamond anvil cells, gas guns and the Z machine were used to validate the EOS, thermal conductivity, specific heat and predictions of water formation. The predicted reaction networks and chemical kinetic equations were implemented in Cheetah and validated against the lower length scale data. Cheetah was then used within the ASC code ALE3D for high-resolution, thermo-mechanically coupled simulations of pore collapse at the micron size scale to predict conditions for detonation initiation.« less

Accounting for microbial habitats in modeling soil organic matter dynamics

NASA Astrophysics Data System (ADS)

Chenu, Claire; Garnier, Patricia; Nunan, Naoise; Pot, Valérie; Raynaud, Xavier; Vieublé, Laure; Otten, Wilfred; Falconer, Ruth; Monga, Olivier

2017-04-01

The extreme heterogeneity of soils constituents, architecture and inhabitants at the microscopic scale is increasingly recognized. Microbial communities exist and are active in a complex 3-D physical framework of mineral and organic particles defining pores of various sizes, more or less inter-connected. This results in a frequent spatial disconnection between soil carbon, energy sources and the decomposer organisms and a variety of microhabitats that are more or less suitable for microbial growth and activity. However, current biogeochemical models account for C dynamics at the macroscale (cm, m) and consider time- and spatially averaged relationships between microbial activity and soil characteristics. Different modelling approaches have intended to account for this microscale heterogeneity, based either on considering aggregates as surrogates for microbial habitats, or pores. Innovative modelling approaches are based on an explicit representation of soil structure at the fine scale, i.e. at µm to mm scales: pore architecture and their saturation with water, localization of organic resources and of microorganisms. Three recent models are presented here, that describe the heterotrophic activity of either bacteria or fungi and are based upon different strategies to represent the complex soil pore system (Mosaic, LBios and µFun). These models allow to hierarchize factors of microbial activity in soil's heterogeneous architecture. Present limits of these approaches and challenges are presented, regarding the extensive information required on soils at the microscale and to up-scale microbial functioning from the pore to the core scale.

Determining Representative Elementary Volume For Multiple Petrophysical Parameters using a Convex Hull Analysis of Digital Rock Data

NASA Astrophysics Data System (ADS)

Shah, S.; Gray, F.; Yang, J.; Crawshaw, J.; Boek, E.

2016-12-01

Advances in 3D pore-scale imaging and computational methods have allowed an exceptionally detailed quantitative and qualitative analysis of the fluid flow in complex porous media. A fundamental problem in pore-scale imaging and modelling is how to represent and model the range of scales encountered in porous media, starting from the smallest pore spaces. In this study, a novel method is presented for determining the representative elementary volume (REV) of a rock for several parameters simultaneously. We calculate the two main macroscopic petrophysical parameters, porosity and single-phase permeability, using micro CT imaging and Lattice Boltzmann (LB) simulations for 14 different porous media, including sandpacks, sandstones and carbonates. The concept of the `Convex Hull' is then applied to calculate the REV for both parameters simultaneously using a plot of the area of the convex hull as a function of the sub-volume, capturing the different scales of heterogeneity from the pore-scale imaging. The results also show that the area of the convex hull (for well-chosen parameters such as the log of the permeability and the porosity) decays exponentially with sub-sample size suggesting a computationally efficient way to determine the system size needed to calculate the parameters to high accuracy (small convex hull area). Finally we propose using a characteristic length such as the pore size to choose an efficient absolute voxel size for the numerical rock.

A multi-scale Lattice Boltzmann model for simulating solute transport in 3D X-ray micro-tomography images of aggregated porous materials

NASA Astrophysics Data System (ADS)

Zhang, Xiaoxian; Crawford, John W.; Flavel, Richard J.; Young, Iain M.

2016-10-01

The Lattice Boltzmann (LB) model and X-ray computed tomography (CT) have been increasingly used in combination over the past decade to simulate water flow and chemical transport at pore scale in porous materials. Because of its limitation in resolution and the hierarchical structure of most natural soils, the X-ray CT tomography can only identify pores that are greater than its resolution and treats other pores as solid. As a result, the so-called solid phase in X-ray images may in reality be a grey phase, containing substantial connected pores capable of conducing fluids and solute. Although modified LB models have been developed to simulate fluid flow in such media, models for solute transport are relatively limited. In this paper, we propose a LB model for simulating solute transport in binary soil images containing permeable solid phase. The model is based on the single-relaxation time approach and uses a modified partial bounce-back method to describe the resistance caused by the permeable solid phase to chemical transport. We derive the relationship between the diffusion coefficient and the parameter introduced in the partial bounce-back method, and test the model against analytical solution for movement of a pulse of tracer. We also validate it against classical finite volume method for solute diffusion in a simple 2D image, and then apply the model to a soil image acquired using X-ray tomography at resolution of 30 μm in attempts to analyse how the ability of the solid phase to diffuse solute at micron-scale affects the behaviour of the solute at macro-scale after a volumetric average. Based on the simulated results, we discuss briefly the danger in interpreting experimental results using the continuum model without fully understanding the pore-scale processes, as well as the potential of using pore-scale modelling and tomography to help improve the continuum models.

Adaptive Multiscale Modeling of Geochemical Impacts on Fracture Evolution

NASA Astrophysics Data System (ADS)

Molins, S.; Trebotich, D.; Steefel, C. I.; Deng, H.

2016-12-01

Understanding fracture evolution is essential for many subsurface energy applications, including subsurface storage, shale gas production, fracking, CO2 sequestration, and geothermal energy extraction. Geochemical processes in particular play a significant role in the evolution of fractures through dissolution-driven widening, fines migration, and/or fracture sealing due to precipitation. One obstacle to understanding and exploiting geochemical fracture evolution is that it is a multiscale process. However, current geochemical modeling of fractures cannot capture this multi-scale nature of geochemical and mechanical impacts on fracture evolution, and is limited to either a continuum or pore-scale representation. Conventional continuum-scale models treat fractures as preferential flow paths, with their permeability evolving as a function (often, a cubic law) of the fracture aperture. This approach has the limitation that it oversimplifies flow within the fracture in its omission of pore scale effects while also assuming well-mixed conditions. More recently, pore-scale models along with advanced characterization techniques have allowed for accurate simulations of flow and reactive transport within the pore space (Molins et al., 2014, 2015). However, these models, even with high performance computing, are currently limited in their ability to treat tractable domain sizes (Steefel et al., 2013). Thus, there is a critical need to develop an adaptive modeling capability that can account for separate properties and processes, emergent and otherwise, in the fracture and the rock matrix at different spatial scales. Here we present an adaptive modeling capability that treats geochemical impacts on fracture evolution within a single multiscale framework. Model development makes use of the high performance simulation capability, Chombo-Crunch, leveraged by high resolution characterization and experiments. The modeling framework is based on the adaptive capability in Chombo which not only enables mesh refinement, but also refinement of the model-pore scale or continuum Darcy scale-in a dynamic way such that the appropriate model is used only when and where it is needed. Explicit flux matching provides coupling betwen the scales.

Assessing the utility of FIB-SEM images for shale digital rock physics

NASA Astrophysics Data System (ADS)

Kelly, Shaina; El-Sobky, Hesham; Torres-Verdín, Carlos; Balhoff, Matthew T.

2016-09-01

Shales and other unconventional or low permeability (tight) reservoirs house vast quantities of hydrocarbons, often demonstrate considerable water uptake, and are potential repositories for fluid sequestration. The pore-scale topology and fluid transport mechanisms within these nanoporous sedimentary rocks remain to be fully understood. Image-informed pore-scale models are useful tools for studying porous media: a debated question in shale pore-scale petrophysics is whether there is a representative elementary volume (REV) for shale models? Furthermore, if an REV exists, how does it differ among petrophysical properties? We obtain three dimensional (3D) models of the topology of microscale shale volumes from image analysis of focused ion beam-scanning electron microscope (FIB-SEM) image stacks and investigate the utility of these models as a potential REV for shale. The scope of data used in this work includes multiple local groups of neighboring FIB-SEM images of different microscale sizes, corresponding core-scale (milli- and centimeters) laboratory data, and, for comparison, series of two-dimensional (2D) cross sections from broad ion beam SEM images (BIB-SEM), which capture a larger microscale field of view than the FIB-SEM images; this array of data is larger than the majority of investigations with FIB-SEM-derived microscale models of shale. Properties such as porosity, organic matter content, and pore connectivity are extracted from each model. Assessments of permeability with single phase, pressure-driven flow simulations are performed in the connected pore space of the models using the lattice-Boltzmann method. Calculated petrophysical properties are compared to those of neighboring FIB-SEM images and to core-scale measurements of the sample associated with the FIB-SEM sites. Results indicate that FIB-SEM images below ∼5000 μm3 volume (the largest volume analyzed) are not a suitable REV for shale permeability and pore-scale networks; i.e. field of view is compromised at the expense of detailed, but often unconnected, nanopore morphology. Further, we find that it is necessary to acquire several local FIB-SEM or BIB-SEM images and correlate their extracted geometric properties to improve the likelihood of achieving representative values of porosity and organic matter volume. Our work indicates that FIB-SEM images of microscale volumes of shale are a qualitative tool for petrophysical and transport analysis. Finally, we offer alternatives for quantitative pore-scale assessments of shale.

Characterization of nanoporous shales with gas sorption

NASA Astrophysics Data System (ADS)

Joewondo, N.; Prasad, M.

2017-12-01

The understanding of the fluid flow in porous media requires the knowledge of the pore system involved. Fluid flow in fine grained shales falls under different regime than transport regime in conventional reservoir due to the different average pore sizes in the two materials; the average pore diameter of conventional sandstones is on the micrometer scale, while of shales can be as small as several nanometers. Mercury intrusion porosimetry is normally used to characterize the pores of conventional reservoir, however with increasingly small pores, the injection pressure required to imbibe the pores becomes infinitely large due to surface tension. Characterization of pores can be expressed by a pore size distribution (PSD) plot, which reflects distribution of pore volume or surface area with respect to pore size. For the case of nanoporous materials, the surface area, which serves as the interface between the rock matrix and fluid, becomes increasingly large and important. Physisorption of gas has been extensively studied as a method of nanoporous solid characterization (particularly for the application of catalysis, metal organic frameworks, etc). The PSD is obtained by matching the experimental result to the calculated theoretical result (using Density Functional Theory (DFT), a quantum mechanics based modelling method for molecular scale interactions). We present the challenges and experimental result of Nitrogen and CO2 gas sorption on shales with various mineralogy and the interpreted PSD obtained by DFT method. Our result shows significant surface area contributed by the nanopores of shales, hence the importance of surface area measurements for the characterization of shales.

Hydromechanics of Reservoir Induced Seismicity

NASA Astrophysics Data System (ADS)

Dura-Gomez, Inmaculada

Data from five reservoirs were analyzed to investigate the various factors and possible pore pressure thresholds associated with Reservoir Induced Seismicity (RIS). Data was obtained from the following reservoirs: Koyna and Warna Reservoirs in India, Itoiz Reservoir in the western Pyrenees, Spain, and Jocassee and Monticello Reservoirs in South Carolina, U.S.A. Koyna Reservoir is one out of four reservoirs in the world where M≥6.0 induced earthquakes have occurred, whereas Warna Reservoir accounts for one out of ten cases with 5.0≤M≤5.9 induced earthquakes. Induced seismicity in the Koyna-Warna region is associated with annual filling cycles in the two reservoirs, large water level changes (30 to 45 m) and the presence of regional scale fractures. The Koyna-Warna case includes 19 M≥5.0 earthquakes at non-repeating hypocenters. The calculation of excess pore pressures associated with these earthquakes suggests values >300 kPa or >600 kPa, before or after 1993 respectively. The need for larger pore pressures from 1993 suggests that M≥5 earthquakes were induced on stronger faults in the region. The exceedance of the previous water level maxima (stress memory) is the most important, although not determining factor in inducing these M≥5.0 earthquakes. Itoiz Reservoir is one of twenty nine reservoirs with 4.0≤M≤4.9 induced earthquakes. The analysis of the RIS associated with the Itoiz Reservoir impoundment, between January 2004 and the end of 2008, shows that that pore pressures diffuse away from Itoiz Reservoir through the carbonate megabreccia systems of the Early to Middle Eocene Hecho Group, and a series of near-vertical thrust faults above the gently dipping Gavarnie thrust. Excess diffused pore pressures destabilize saturated critically stressed seismogenic fractures where RIS takes place. In particular, M≥3.0 earthquakes in the region are associated with excess pore pressures of the order of 100 to 200 kPa. Jocassee and Monticello Reservoirs in South Carolina are among the best studied reservoir induced seismicity cases in the world, and have been associated with M<4 and M<3 earthquakes respectively. The analysis of some of these earthquakes emphasizes the contribution of the diffused pore pressures to the observed seismicity. In the case of Jocassee Reservoir, 2.0≤M≤3.0 earthquakes occurred in a homogenous, non-fractured crystalline rock and are associated with excess diffused pore pressures of the order of 600 to 700 kPa. Earthquakes associated with the impoundment of Monticello Reservoir occurred in a region of very complex geology, with many pre-existing local scale fractures. The calculation of the excess diffused pore pressures associated with 2.0≤M≤3.0 earthquakes yielded values of the order of 100 to 300 kPa. Synthesis of these data show that RIS occurs when excess pore pressures (which occur primarily by diffusion) reach threshold values needed to induce RIS. The occurrence of RIS and its magnitude are controlled by the filling history, availability of fluid filled saturated fractures and their hydrogeological properties.

Morphological analysis of pore size and connectivity in a thick mixed-culture biofilm.

PubMed

Rosenthal, Alex F; Griffin, James S; Wagner, Michael; Packman, Aaron I; Balogun, Oluwaseyi; Wells, George F

2018-05-19

Morphological parameters are commonly used to predict transport and metabolic kinetics in biofilms. Yet, quantification of biofilm morphology remains challenging due to imaging technology limitations and lack of robust analytical approaches. We present a novel set of imaging and image analysis techniques to estimate internal porosity, pore size distributions, and pore network connectivity to a depth of 1 mm at a resolution of 10 µm in a biofilm exhibiting both heterotrophic and nitrifying activity. Optical coherence tomography (OCT) scans revealed an extensive pore network with diameters as large as 110 µm directly connected to the biofilm surface and surrounding fluid. Thin section fluorescence in situ hybridization microscopy revealed ammonia oxidizing bacteria (AOB) distributed through the entire thickness of the biofilm. AOB were particularly concentrated in the biofilm around internal pores. Areal porosity values estimated from OCT scans were consistently lower than those estimated from multiphoton laser scanning microscopy, though the two imaging modalities showed a statistically significant correlation (r = 0.49, p<0.0001). Estimates of areal porosity were moderately sensitive to grey level threshold selection, though several automated thresholding algorithms yielded similar values to those obtained by manually thresholding performed by a panel of environmental engineering researchers (±25% relative error). These findings advance our ability to quantitatively describe the geometry of biofilm internal pore networks at length scales relevant to engineered biofilm reactors and suggest that internal pore structures provide crucial habitat for nitrifier growth. This article is protected by copyright. All rights reserved. This article is protected by copyright. All rights reserved.

Flow through internal elastic lamina affects shear stress on smooth muscle cells (3D simulations).

PubMed

Tada, Shigeru; Tarbell, John M

2002-02-01

We describe a three-dimensional numerical simulation of interstitial flow through the medial layer of an artery accounting for the complex entrance condition associated with fenestral pores in the internal elastic lamina (IEL) to investigate the fluid mechanical environment around the smooth muscle cells (SMCs) right beneath the IEL. The IEL was modeled as an impermeable barrier to water flow except for the fenestral pores, which were assumed to be uniformly distributed over the IEL. The medial layer was modeled as a heterogeneous medium composed of a periodic array of cylindrical SMCs embedded in a continuous porous medium representing the interstitial proteoglycan and collagen matrix. Depending on the distance between the IEL bottom surface and the upstream end of the proximal layer of SMCs, the local shear stress on SMCs right beneath the fenestral pore could be more than 10 times higher than that on the cells far removed from the IEL under the conditions that the fenestral pore diameter and area fraction of pores were kept constant at 1.4 microm and 0.05, respectively. Thus these proximal SMCs may experience shear stress levels that are even higher than endothelial cells exposed to normal blood flow (order of 10 dyn/cm(2)). Furthermore, entrance flow through fenestral pores alters considerably the interstitial flow field in the medial layer over a spatial length scale of the order of the fenestral pore diameter. Thus the spatial gradient of shear stress on the most superficial SMC is noticeably higher than computed for endothelial cell surfaces.

A multi-scale network method for two-phase flow in porous media

DOE Office of Scientific and Technical Information (OSTI.GOV)

Khayrat, Karim, E-mail: khayratk@ifd.mavt.ethz.ch; Jenny, Patrick

Pore-network models of porous media are useful in the study of pore-scale flow in porous media. In order to extract macroscopic properties from flow simulations in pore-networks, it is crucial the networks are large enough to be considered representative elementary volumes. However, existing two-phase network flow solvers are limited to relatively small domains. For this purpose, a multi-scale pore-network (MSPN) method, which takes into account flow-rate effects and can simulate larger domains compared to existing methods, was developed. In our solution algorithm, a large pore network is partitioned into several smaller sub-networks. The algorithm to advance the fluid interfaces withinmore » each subnetwork consists of three steps. First, a global pressure problem on the network is solved approximately using the multiscale finite volume (MSFV) method. Next, the fluxes across the subnetworks are computed. Lastly, using fluxes as boundary conditions, a dynamic two-phase flow solver is used to advance the solution in time. Simulation results of drainage scenarios at different capillary numbers and unfavourable viscosity ratios are presented and used to validate the MSPN method against solutions obtained by an existing dynamic network flow solver.« less

Methane and CO2 Adsorption and Transport in Carbon-based Systems from Experiments and Molecular Simulation

NASA Astrophysics Data System (ADS)

Wilcox, Jennifer; Firouzi, Mahnaz; Rupp, Erik; Haghapanah, Reza; Wang, Beibei

2013-04-01

Carbon capture and sequestration is one strategy that could potentially mitigate gigatons of CO2 emissions per year; however, technical obstacles have thus far hindered wide-scale deployment of this strategy. To design efficient and reliable strategies for either carbon capture or sequestration at the full-scale, one needs to understand the chemical and physical properties of CO2 and its interaction with its local surroundings at the molecular-scale. To investigate the chemical and physical properties of CO2 and its local surroundings at the molecular-scale, surface characterization studies are carried out alongside theoretical model efforts. Experimental investigation of CO2 interactions with organic-based porous materials ranging in complexity from functionalized graphene and activated carbon to various-rank coal and gas shale samples to create a set of realistic models that take into account both surface and pore heterogeneity. Integration of theory and experiments takes place to allow for the relevant physics at the molecular-level to be revealed. Determining adsorption and transport phenomena of CO2 (and mixtures, including H2O, and CH4) within the model pore systems can be used to understand the complex pore matrices of carbon-based sorbents, coal, and the organic components of gas shale that are crucial to determining their carbon capture or sequestration potential. Non-equilibrium molecular dynamics (NEMD) simulations of pure carbon dioxide, methane, helium and their mixtures have been carried out in carbon slit pores to investigate gas slippage and Klinkenberg effects in the organic matrices of coal and gas shale rocks. NEMD techniques are ideally suited for the experimental situation in which an external driving force, such as a chemical potential or pressure gradient, are applied on the system. Simulations have been conducted to determine the effect of pore size and exposure to an external potential on the velocity profile and slip-stick boundary conditions. The simulations indicate that molecule-wall collisions influence the velocity profile, which deviates significantly from the Navier-Stokes hydrodynamic prediction for micro and mesopores. Also, the shape of the velocity profile is found to be independent of the applied pressure gradient in micropores. The results indicate that the velocity profile is uniform for pore sizes less than 2 nm (micropores). As pore sizes increase to 10 nm, parabolic profiles are observed due to the reduced interaction of gas molecules with the pore walls. Interestingly, in small pores unlike in large pores, the gas velocity at the walls is non-zero and predicted gas transport is somewhat enhanced as the gas flow transitions from a parabolic velocity profile to plug-flow. In addition, a 3-D pore network, representative of porous carbon-based materials, has been generated atomistically using the Voronoi tessellation method. Simulations have been carried out to determine the effect of the pore structure and modeled viscosity on permeability and Klinkenberg parameters. The use of the bulk-phase viscosity for estimating the permeability of CO2 in units of Darcy in a 3-D micropore network is not an appropriate assumption as it significantly underestimates the CO2 permeability given that CO2 is an adsorbing gas with strong pore wall interactions. On the other hand, since the transport properties of CH4 are less influenced by the pore walls compared with CO2, the use of the bulk-phase CH4 viscosity estimates are a reasonable assumption. The application of this work is to advance our understanding of gas transport and to provide insight into mechanisms of gas-surface interactions in the complex natural systems such as gas shale so that we can make accurate capacity estimates in addition to assisting in enhancing natural gas recovery from these systems. These results will potentially have important implications on CO2 adsorption and transport in carbon-based materials and geologic formations and may provide an understanding of the limitations of the use of bulk-phase fluid viscosities to model transport properties for nanoconfined fluids.

Investigation of structural-scale levels of spall fracture induced by a nanosecond relativistic high-current electron beam in ultrafine-grained Ti-Al-V-Mo alloy

NASA Astrophysics Data System (ADS)

Dudarev, E. F.; Markov, A. B.; Bakach, G. P.; Maletkina, T. Yu.; Belov, N. N.; Tabachenko, A. N.; Skosirskii, A. B.; Habibullin, M. V.; Yakovlev, E. V.

2017-12-01

The results of an experimental and theoretical study of shock-wave processes and spall fracture in an ultrafine-grained and coarse-grained (α + β) Ti-Al-V-Mo alloy under the action of a nanosecond relativistic high-current electron beam are reported. Mathematical modeling is performed to show that when an electron beam with a power density of 1.65 × 1010 W/cm2 impacts this alloy, a shock wave with a compression amplitude of 13 GPa appears and its reflection gives rise to a tensile wave. Its amplitude increases with decreasing target thickness. The calculated increase in the thickness of the spalled layer at the rear surface of the target corresponds to the experimental data. It is established experimentally that plastic deformation precedes the spall fracture sequentially at three structural-scale levels. At the beginning pores are formed and merge, then microcracks are formed at different angles to the back surface of the target between the pores, and then a macrocrack is formed. As a result, the macrocrack surface is not smooth but exhibits pits of ductile fracture.

Direct pore-scale reactive transport modelling of dynamic wettability changes induced by surface complexation

NASA Astrophysics Data System (ADS)

Maes, Julien; Geiger, Sebastian

2018-01-01

Laboratory experiments have shown that oil production from sandstone and carbonate reservoirs by waterflooding could be significantly increased by manipulating the composition of the injected water (e.g. by lowering the ionic strength). Recent studies suggest that a change of wettability induced by a change in surface charge is likely to be one of the driving mechanism of the so-called low-salinity effect. In this case, the potential increase of oil recovery during waterflooding at low ionic strength would be strongly impacted by the inter-relations between flow, transport and chemical reaction at the pore-scale. Hence, a new numerical model that includes two-phase flow, solute reactive transport and wettability alteration is implemented based on the Direct Numerical Simulation of the Navier-Stokes equations and surface complexation modelling. Our model is first used to match experimental results of oil droplet detachment from clay patches. We then study the effect of wettability change on the pore-scale displacement for simple 2D calcite micro-models and evaluate the impact of several parameters such as water composition and injected velocity. Finally, we repeat the simulation experiments on a larger and more complex pore geometry representing a carbonate rock. Our simulations highlight two different effects of low-salinity on oil production from carbonate rocks: a smaller number of oil clusters left in the pores after invasion, and a greater number of pores invaded.

Coupled Electro-Hydrodynamic Effects of Electro-Osmosis from Pore Scale to Darcy Scale

NASA Astrophysics Data System (ADS)

Schotting, R.; Joekar-Niasar, V.; Leijnse, A.

2011-12-01

Electro-osmosis is "movement of a fluid under the effect of an electric field in a porous medium". This phenomenon has many applications in civil engineering (slope stabilization, dewatering), environmental engineering (soil remediation, sludge dewatering), chemical engineering (micro- or nano- mixers), medical engineering (drug delivery), etc. The key factor in electro-osmosis is the competition between the electrochemical and hydrodynamic forces as well as the coupling between the solid surface and the electrolyte properties. The objective of this research is to understand the influence of pore-scale heterogeneities of surface properties on the Darcy-scale behavior. We develop novel analytical solutions for the flow and transport of electrolyte including electro-hydrodynamic forces in a single micro-channel. We propose the complete analytical solution for monovalent electrolyte at full range overlapping double layers, and nonlinear electric field, including the Donan effect in transport of ions. These pore-scale formulations are numerically upscaled to obtain the Darcy-scale behavior. Our results show the contribution of electro-osmotic, chemical-osmotic and hydrodynamic components of the flow equation on pressure field evolution and multi-directional flow field at Darcy scale.

Using synchrotron-based X-ray micro-computed tomography and high-performance pore-scale simulation to evaluate hydraulic properties in biochar-amended soils

NASA Astrophysics Data System (ADS)

Zhou, H.; Yu, X.; Chen, C.; Zeng, L.; Lu, S.; Wu, L.

2016-12-01

In this research, we combined synchrotron-based X-ray micro-computed tomography (SR-mCT), with three-dimensional lattice Bolzmann (LB) method, to quantify how the change in pore space architecture affected macroscopic hydraulic of two clayey soils amended with biochar. SR-mCT was used to characterize pore structures of the soils before and after biochar addition. The high-resolution soil pore structures were then directly used as internal boundary conditions for three-dimensional water flow simulations with the LB method, which was accelerated by graphics processing unit (GPU) parallel computing. It was shown that, due to the changes in soil pore geometry, the application of biochar increased the soil permeability by at least 1 order of magnitude, and decreased the tortuosity by 20-30%. This work was the first physics based modeling study on the effect of biochar amendment on soil hydraulic properties. The developed theories and techniques have promising potential in understanding the mechanisms of water and nutrient transport in soil at the pore scale.

The effect of pore-scale geometry and wettability on two-phase relative permeabilities within elementary cells

NASA Astrophysics Data System (ADS)

Bianchi Janetti, Emanuela; Riva, Monica; Guadagnini, Alberto

2017-04-01

We study the relative role of the complex pore space geometry and wettability of the solid matrix on the quantification of relative permeabilities characterizing steady state immiscible two-phase flow in porous media. We do so by considering elementary cells, which are typically employed in upscaling frameworks based on, e.g., homogenization or volume averaging. In this context one typically relies on the solution of pore-scale physics at a scale which is much smaller than that of an investigated porous system. Pressure-driven two-phase flow following simultaneous co-current injection of water and oil is numerically solved for a suite of regular and stochastically generated two-dimensional explicit elementary cells with fixed porosity and sharing main topological/morphological features. We show that relative permeabilities of the randomly generated elementary cells are significantly influenced by the formation of preferential percolation paths (principal pathways), giving rise to a strongly nonuniform distribution of fluid fluxes. These pathways are a result of the spatially variable resistance that the random pore structures exert on the fluid. The overall effect on relative permeabilities of the diverse organization of principal pathways, as driven by a given random realization at the scale of the unit cell, is significantly larger than that of the wettability of the host rock. In contrast to what can be observed for the random cells analyzed, relative permeabilities of regular cells display a clear trend with contact angle at the investigated scale. Our findings suggest the need to perform systematic upscaling studies in a stochastic context, to propagate the effects of uncertain pore space geometries to a probabilistic description of relative permeability curves at the continuum scale.

Characterization of Mixed Wettability at Different Scales and its Impact on Oil Recovery Efficiency

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sharma, Mukul M.; Hirasaki, George J.

The objectives of this project was to: (1) quantify the pore scale mechanisms that determine the wettability state of a reservoir, (2) study the effect of crude oil, brine and mineral compositions in the establishment of mixed wet states, (3) clarify the effect of mixed - wettability on oil displacement efficiency in waterfloods, (4) develop a new tracer technique to measure wettability, fluid distributions, residual saturation's and relative permeabilities, and (5) develop methods for properly incorporating wettability in up-scaling from pore to core to reservoir scales.

Plasma flows and magnetic field interplay during the formation of a pore

NASA Astrophysics Data System (ADS)

Ermolli, I.; Cristaldi, A.; Giorgi, F.; Giannattasio, F.; Stangalini, M.; Romano, P.; Tritschler, A.; Zuccarello, F.

2017-04-01

Aims: Recent simulations of solar magneto-convection have offered new levels of understanding of the interplay between plasma motions and magnetic fields in evolving active regions. We aim at verifying some aspects of the formation of magnetic regions derived from recent numerical studies in observational data. Methods: We studied the formation of a pore in the active region (AR) NOAA 11462. We analysed data obtained with the Interferometric Bidimensional Spectrometer (IBIS) at the Dunn Solar Telescope on April 17, 2012, consisting of full Stokes measurements of the Fe I 617.3 nm lines. Furthermore, we analysed SDO/HMI observations in the continuum and vector magnetograms derived from the Fe I 617.3 nm line data taken from April 15 to 19, 2012. We estimated the magnetic field strength and vector components and the line-of-sight (LOS) and horizontal motions in the photospheric region hosting the pore formation. We discuss our results in light of other observational studies and recent advances of numerical simulations. Results: The pore formation occurs in less than 1 h in the leading region of the AR. We observe that the evolution of the flux patch in the leading part of the AR is faster (<12 h) than the evolution (20-30 h) of the more diffuse and smaller scale flux patches in the trailing region. During the pore formation, the ratio between magnetic and dark area decreases from 5 to 2. We observe strong downflows at the forming pore boundary and diverging proper motions of plasma in the vicinity of the evolving feature that are directed towards the forming pore. The average values and trends of the various quantities estimated in the AR are in agreement with results of former observational studies of steady pores and with their modelled counterparts, as seen in recent numerical simulations of a rising-tube process. The agreement with the outcomes of the numerical studies holds for both the signatures of the flux emergence process (e.g. appearance of small-scale mixed polarity patterns and elongated granules) and the evolution of the region. The processes driving the formation of the pore are identified with the emergence of a magnetic flux concentration and the subsequent reorganization of the emerged flux, by the combined effect of velocity and magnetic field, in and around the evolving structure. Movies associated to Figs. 1 and 4 are available at http://www.aanda.org

Bead-Based Microfluidic Sediment Analogues: Fabrication and Colloid Transport.

PubMed

Guo, Yang; Huang, Jingwei; Xiao, Feng; Yin, Xiaolong; Chun, Jaehun; Um, Wooyong; Neeves, Keith B; Wu, Ning

2016-09-13

Mobile colloids can act as carriers for low-solubility contaminants in the environment. However, the dominant mechanism for this colloid-facilitated transport of chemicals is unclear. Therefore, we developed a bead-based microfluidic platform of sediment analogues and measured both single and population transport of model colloids. The porous medium is assembled through a bead-by-bead injection method. This approach has the versatility to build both electrostatically homogeneous and heterogeneous media at the pore scale. A T-junction at the exit also allowed for encapsulation and enumeration of colloids effluent at single particle resolution to give population dynamics. Tortuosity calculated from pore-scale trajectory analysis and its comparison with lattice Boltzmann simulations revealed that transport of colloids was influenced by the size exclusion effect. The porous media packed by positively and negatively charged beads into two layers showed distinctive colloidal particle retention and significant remobilization and re-adsorption of particles during water flushing. We demonstrated the potential of our method to fabricate porous media with surface heterogeneities at the pore scale. With both single and population dynamics measurement, our platform has the potential to connect pore-scale and macroscale colloid transport on a lab scale and to quantify the impact of grain surface heterogeneities that are natural in the subsurface environment.

Multiscale Pore Throat Network Reconstruction of Tight Porous Media Constrained by Mercury Intrusion Capillary Pressure and Nuclear Magnetic Resonance Measurements

NASA Astrophysics Data System (ADS)

Xu, R.; Prodanovic, M.

2017-12-01

Due to the low porosity and permeability of tight porous media, hydrocarbon productivity strongly depends on the pore structure. Effective characterization of pore/throat sizes and reconstruction of their connectivity in tight porous media remains challenging. Having a representative pore throat network, however, is valuable for calculation of other petrophysical properties such as permeability, which is time-consuming and costly to obtain by experimental measurements. Due to a wide range of length scales encountered, a combination of experimental methods is usually required to obtain a comprehensive picture of the pore-body and pore-throat size distributions. In this work, we combine mercury intrusion capillary pressure (MICP) and nuclear magnetic resonance (NMR) measurements by percolation theory to derive pore-body size distribution, following the work by Daigle et al. (2015). However, in their work, the actual pore-throat sizes and the distribution of coordination numbers are not well-defined. To compensate for that, we build a 3D unstructured two-scale pore throat network model initialized by the measured porosity and the calculated pore-body size distributions, with a tunable pore-throat size and coordination number distribution, which we further determine by matching the capillary pressure vs. saturation curve from MICP measurement, based on the fact that the mercury intrusion process is controlled by both the pore/throat size distributions and the connectivity of the pore system. We validate our model by characterizing several core samples from tight Middle East carbonate, and use the network model to predict the apparent permeability of the samples under single phase fluid flow condition. Results show that the permeability we get is in reasonable agreement with the Coreval experimental measurements. The pore throat network we get can be used to further calculate relative permeability curves and simulate multiphase flow behavior, which will provide valuable insights into the production optimization and enhanced oil recovery design.

Inertial effects during irreversible meniscus reconfiguration in angular pores

NASA Astrophysics Data System (ADS)

Ferrari, Andrea; Lunati, Ivan

2014-12-01

In porous media, the dynamics of the invading front between two immiscible fluids is often characterized by abrupt reconfigurations caused by local instabilities of the interface. As a prototype of these phenomena we consider the dynamics of a meniscus in a corner as it can be encountered in angular pores. We investigate this process in detail by means of direct numerical simulations that solve the Navier-Stokes equations in the pore space and employ the Volume of Fluid method (VOF) to track the evolution of the interface. We show that for a quasi-static displacement, the numerically calculated surface energy agrees well with the analytical solutions that we have derived for pores with circular and square cross sections. However, the spontaneous reconfigurations are irreversible and cannot be controlled by the injection rate: they are characterized by the amount of surface energy that is spontaneously released and transformed into kinetic energy. The resulting local velocities can be orders of magnitude larger than the injection velocity and they induce damped oscillations of the interface that possess their own time scales and depend only on fluid properties and pore geometry. In complex media (we consider a network of cubic pores) reconfigurations are so frequent and oscillations last long enough that increasing inertial effects leads to a different fluid distribution by influencing the selection of the next pore to be invaded. This calls into question simple pore-filling rules based only on capillary forces. Also, we demonstrate that inertial effects during irreversible reconfigurations can influence the work done by the external forces that is related to the pressure drop in Darcy's law. This suggests that these phenomena have to be considered when upscaling multiphase flow because local oscillations of the menisci affect macroscopic quantities and modify the constitutive relationships to be used in macro-scale models. These results can be extrapolated to other interface instabilities that are at the origin of fast pore-scale events, such as Haines jumps, snap-off and coalescence.

Pore-Fractal Structure in Porous Carbons Made from Corn and Wheat

NASA Astrophysics Data System (ADS)

Kapoor, Y. M.; Schmidt, P. W.; Rice, Randall D.; Shulse, Laural; Voss, D. J.; Venkatraman, A.; Fan, L. T.; Walawender, W. P.; Rieker, T. P.

1998-03-01

Small-angle X-ray scattering has been used in a study of the pore structure of some porous and activated carbons on length scales between about 5 and 10^4 ÅThe carbons were obtained by pyrolysis and activation of wheat and American corn (maize). The scattering data showed that in each carbon there are at least two of the following four types of pores: (1) pores with diameters of at least 10^4 Åpores with smooth or fractal surfaces and diameters of at least 5 x 10^3 Åpore-fractals with diameters of no more than about 10^3 Åand (4) pores with diameters no larger than 100 ÅThe relation between the pore structure and the procedure used to obtain the carbon and will be discussed.

Hybrid discrete-continuum modeling for transport, biofilm development and solid restructuring including electrostatic effects

NASA Astrophysics Data System (ADS)

Prechtel, Alexander; Ray, Nadja; Rupp, Andreas

2017-04-01

We want to present an approach for the mathematical, mechanistic modeling and numerical treatment of processes leading to the formation, stability, and turnover of soil micro-aggregates. This aims at deterministic aggregation models including detailed mechanistic pore-scale descriptions to account for the interplay of geochemistry and microbiology, and the link to soil functions as, e.g., the porosity. We therefore consider processes at the pore scale and the mesoscale (laboratory scale). At the pore scale transport by diffusion, advection, and drift emerging from electric forces can be taken into account, in addition to homogeneous and heterogeneous reactions of species. In the context of soil micro-aggregates the growth of biofilms or other glueing substances as EPS (extracellular polymeric substances) is important and affects the structure of the pore space in space and time. This model is upscaled mathematically in the framework of (periodic) homogenization to transfer it to the mesoscale resulting in effective coefficients/parameters there. This micro-macro model thus couples macroscopic equations that describe the transport and fluid flow at the scale of the porous medium (mesoscale) with averaged time- and space-dependent coefficient functions. These functions may be explicitly computed by means of auxiliary cell problems (microscale). Finally, the pore space in which the cell problems are defined is time and space dependent and its geometry inherits information from the transport equation's solutions. The microscale problems rely on versatile combinations of cellular automata and discontiuous Galerkin methods while on the mesoscale mixed finite elements are used. The numerical simulations allow to study the interplay between these processes.

Multiscale pore networks and their effect on deformation and transport property alteration associated with hydraulic fracturing

NASA Astrophysics Data System (ADS)

Daigle, Hugh; Hayman, Nicholas; Jiang, Han; Tian, Xiao; Jiang, Chunbi

2017-04-01

Multiple lines of evidence indicate that, during a hydraulic fracture stimulation, the permeability of the unfractured matrix far from the main, induced tensile fracture increases by one to two orders of magnitude. This permeability enhancement is associated with pervasive shear failure in a large region surrounding the main induced fracture. We have performed low-pressure gas sorption, mercury intrusion, and nuclear magnetic resonance measurements along with high-resolution scanning electron microscope imaging on several preserved and unpreserved shale samples from North American basins before and after inducing failure in confined compressive strength tests. We have observed that the pore structure in intact samples exhibits multiscale behavior, with sub-micron-scale pores in organic matter connected in isolated, micron-scale clusters which themselves are connected to each other through a network of microcracks. The organic-hosted pore networks are poorly connected due to a significant number of dead-end pores within the organic matter. Following shear failure, we often observe an increase in pore volume in the sub-micron range, which appears to be related to the formation of microcracks that propagate along grain boundaries and other planes of mechanical strength contrast. This is consistent with other experimental and field evidence. In some cases these microcracks cross or terminate in organic matter, intersecting the organic-hosted pores. The induced microcrack networks typically have low connectivity and do not appreciably increase the connectivity of the overall pore network. However, in other cases the shear deformation results in an overall pore volume decrease; samples which exhibit this behavior tend to have more clay minerals. Our interpretation of these phenomena is as follows. As organic matter is converted to hydrocarbons, organic-hosted pores develop, and the hydrocarbons contained in these pores are overpressured. The disconnected nature of these clusters of organic-hosted pores prevents the overpressure from dissipating, resulting in localized overpressure at the micron scale. When the rock is subjected to a hydraulic fracture stimulation, the rock surrounding the main induced fracture experiences shear deformation. Those parts of the rock that contain overpressured fluids in the organic-hosted pores will be more likely to experience dilatancy in the form of brittle deformation; the portions of the rock lacking in organic-hosted pores will tend to experience compactive shear failure since the effective normal stresses are larger. The microcrack networks that propagate into the regions of organic-hosted porosity allow the hydrocarbons resident in those pores to migrate to the main induced tensile fractures. The disconnected nature of the microcrack networks causes only a slight increase in permeability, which is consistent with other observations. Our work illustrates how multiscale pore networks in shale interact with in situ stresses to affect the bulk shale rheology.

Flow and Transport in Complex Microporous Carbonates as a Consequence of Separation of Scales

NASA Astrophysics Data System (ADS)

Bijeljic, B.; Raeini, A. Q.; Lin, Q.; Blunt, M. J.

2017-12-01

Some of the most important examples of flow and transport in complex pore structures are found in subsurface applications such as contaminant hydrology, carbon storage and enhanced oil recovery. Carbonate rock structures contain most of the world's oil reserves, considerable amount of water reserves, and potentially hold a storage capacity for carbon dioxide. However, this type of pore space is difficult to represent due to complexities associated with a wide range of pore sizes and variation in connectivity which poses a considerable challenge for quantitative predictions of transport across multiple scales.A new concept unifying X-ray tomography experiment and direct numerical simulation has been developed that relies on full description flow and solute transport at the pore scale. Differential imaging method (Lin et al. 2016) provides rich information in microporous space, while advective and diffusive mass transport are simulated on micro-CT images of pore-space: Navier-Stokes equations are solved for flow in the image voxels comprising the pore space, streamline-based simulation is used to account for advection, and diffusion is superimposed by random walk.Quantitative validation has been done on analytical solutions for diffusion and by comparing the model predictions versus the experimental NMR measurements in the dual porosity beadpack. Furthermore, we discriminate signatures of multi-scale transport behaviour for a range of carbonate rock (Figure 1), dependent on the heterogeneity of the inter- and intra-grain pore space, heterogeneity in the flow field, and the mass transfer characteristics of the porous media. Finally, we demonstrate the predictive capabilities of the model through an analysis that includes a number of probability density functions flow and transport (PDFs) measures of non-Fickian transport on the micro-CT i935mages. In complex porous media separation of scales exists, leading to flow and transport signatures that need to be described by multiple functions with distinct flow field and transport characteristics. Reference: Lin, Q., Al-Khulaifi Y., Blunt, M.J. and Bijeljic B. (2016). Advances in Water Resources, 96, 306-322, doi:10.1016/j.advwatres.2016.08.002.

Pore-scale simulation of wettability and interfacial tension effects on flooding process for enhanced oil recovery.

PubMed

Zhao, Jin; Wen, Dongsheng

2017-08-27

For enhanced oil recovery (EOR) applications, the oil/water flow characteristics during the flooding process was numerically investigated with the volume-of-fluid method at the pore scale. A two-dimensional pore throat-body connecting structure was established, and four scenarios were simulated in this paper. For oil-saturated pores, the wettability effect on the flooding process was studied; for oil-unsaturated pores, three effects were modelled to investigate the oil/water phase flow behaviors, namely the wettability effect, the interfacial tension (IFT) effect, and the combined wettability/IFT effect. The results show that oil saturated pores with the water-wet state can lead to 25-40% more oil recovery than with the oil-wet state, and the remaining oil mainly stays in the near wall region of the pore bodies for oil-wet saturated pores. For oil-unsaturated pores, the wettability effects on the flooding process can help oil to detach from the pore walls. By decreasing the oil/water interfacial tension and altering the wettability from oil-wet to water-wet state, the remaining oil recovery rate can be enhanced successfully. The wettability-IFT combined effect shows better EOR potential compared with decreasing the interfacial tension alone under the oil-wet condition. The simulation results in this work are consistent with previous experimental and molecular dynamics simulation conclusions. The combination effect of the IFT reducation and wettability alteration can become an important recovery mechanism in future studies for nanoparticles, surfactant, and nanoparticle-surfactant hybrid flooding process.

Pore-scale simulation of wettability and interfacial tension effects on flooding process for enhanced oil recovery

PubMed Central

Zhao, Jin

2017-01-01

For enhanced oil recovery (EOR) applications, the oil/water flow characteristics during the flooding process was numerically investigated with the volume-of-fluid method at the pore scale. A two-dimensional pore throat-body connecting structure was established, and four scenarios were simulated in this paper. For oil-saturated pores, the wettability effect on the flooding process was studied; for oil-unsaturated pores, three effects were modelled to investigate the oil/water phase flow behaviors, namely the wettability effect, the interfacial tension (IFT) effect, and the combined wettability/IFT effect. The results show that oil saturated pores with the water-wet state can lead to 25–40% more oil recovery than with the oil-wet state, and the remaining oil mainly stays in the near wall region of the pore bodies for oil-wet saturated pores. For oil-unsaturated pores, the wettability effects on the flooding process can help oil to detach from the pore walls. By decreasing the oil/water interfacial tension and altering the wettability from oil-wet to water-wet state, the remaining oil recovery rate can be enhanced successfully. The wettability-IFT combined effect shows better EOR potential compared with decreasing the interfacial tension alone under the oil-wet condition. The simulation results in this work are consistent with previous experimental and molecular dynamics simulation conclusions. The combination effect of the IFT reducation and wettability alteration can become an important recovery mechanism in future studies for nanoparticles, surfactant, and nanoparticle–surfactant hybrid flooding process. PMID:29308190

Evolution of the pore structure during the early stages of the alkali-activation reaction: An in situ small-angle neutron scattering investigation

DOE PAGES

White, Claire E.; Olds, Daniel P.; Hartl, Monika; ...

2017-02-01

The long-term durability of cement-based materials is influenced by the pore structure and associated permeability at the sub-micrometre length scale. With the emergence of new types of sustainable cements in recent decades, there is a pressing need to be able to predict the durability of these new materials, and therefore nondestructive experimental techniques capable of characterizing the evolution of the pore structure are increasingly crucial for investigating cement durability. Here, small-angle neutron scattering is used to analyze the evolution of the pore structure in alkali-activated materials over the initial 24 h of reaction in order to assess the characteristic poremore » sizes that emerge during these short time scales. By using a unified fitting approach for data modeling, information on the pore size and surface roughness is obtained for a variety of precursor chemistries and morphologies (metakaolin- and slag-based pastes). Furthermore, the impact of activator chemistry is elucidatedviathe analysis of pastes synthesized using hydroxide- and silicate-based activators. It is found that the main aspect influencing the size of pores that are accessible using small-angle neutron scattering analysis (approximately 10–500 Å in diameter) is the availability of free silica in the activating solution, which leads to a more refined pore structure with smaller average pore size. Furthermore, as the reaction progresses the gel pores visible using this scattering technique are seen to increase in size.« less

Pore size distribution and accessible pore size distribution in bituminous coals

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sakurovs, Richard; He, Lilin; Melnichenko, Yuri B

2012-01-01

The porosity and pore size distribution of coals determine many of their properties, from gas release to their behavior on carbonization, and yet most methods of determining pore size distribution can only examine a restricted size range. Even then, only accessible pores can be investigated with these methods. Small-angle neutron scattering (SANS) and ultra small-angle neutron scattering (USANS) are increasingly used to characterize the size distribution of all of the pores non-destructively. Here we have used USANS/SANS to examine 24 well-characterized bituminous and subbituminous coals: three from the eastern US, two from Poland, one from New Zealand and the restmore » from the Sydney and Bowen Basins in Eastern Australia, and determined the relationships of the scattering intensity corresponding to different pore sizes with other coal properties. The range of pore radii examinable with these techniques is 2.5 nm to 7 {micro}m. We confirm that there is a wide range of pore sizes in coal. The pore size distribution was found to be strongly affected by both rank and type (expressed as either hydrogen or vitrinite content) in the size range 250 nm to 7 {micro}m and 5 to 10 nm, but weakly in intermediate regions. The results suggest that different mechanisms control coal porosity on different scales. Contrast-matching USANS and SANS were also used to determine the size distribution of the fraction of the pores in these coals that are inaccessible to deuterated methane, CD{sub 4}, at ambient temperature. In some coals most of the small ({approx} 10 nm) pores were found to be inaccessible to CD{sub 4} on the time scale of the measurement ({approx} 30 min - 16 h). This inaccessibility suggests that in these coals a considerable fraction of inherent methane may be trapped for extended periods of time, thus reducing the effectiveness of methane release from (or sorption by) these coals. Although the number of small pores was less in higher rank coals, the fraction of total pores that was inaccessible was not rank dependent. In the Australian coals, at the 10 nm to 50 nm size scales the pores in inertinites appeared to be completely accessible to CD{sub 4}, whereas the pores in the vitrinite were about 75% inaccessible. Unlike the results for total porosity that showed no regional effects on relationships between porosity and coal properties, clear regional differences in the relationships between fraction of closed porosity and coal properties were found. The 10 to 50 nm-sized pores of inertinites of the US and Polish coals examined appeared less accessible to methane than those of the inertinites of Australian coals. This difference in pore accessibility in inertinites may explain why empirical relationships between fluidity and coking properties developed using Carboniferous coals do not apply to Australian coals.« less

4D XMT of Reaction in Carbonates: Reactive Transport Dynamics at Multiples Scales

NASA Astrophysics Data System (ADS)

Menke, H. P.; Reynolds, C. A.; Andrew, M. G.; Nunes, J. P. P.; Bijeljic, B.; Blunt, M. J.

2016-12-01

Upscaling pore scale rock-fluid interaction processes for predictive modelling poses a challenge to underground carbon storage. We have completed experiments and flow modelling to investigate the impact of pore-space heterogeneity and scale on the dissolution of two limestones at both the mm and cm scales. Two samples were reacted with reservoir condition CO2-saturated brine at both scales and scanned dynamically as dissolution took place. First, 1-cm long 4-mm diameter micro cores were scanned during reactive flow at a 4-μm resolution between 4 and 40 times using 4D X-ray micro-tomography over the course of 1.5 hours using a laboratory μ-CT. Second, 3.8-cm diameter, 8-cm long macro cores were reacted at the same conditions inside a reservoir condition flow rig and imaged using a medical CT scanner. Each sample was imaged 10 times over the course of 1.5 hours at a 250 x 250 x 500-μm resolution. The reacted macro cores were then scanned inside a μ-CT at a 27-μm resolution to assess the alteration in pore-scale reaction-induced heterogeneity. It was found that both limestones showed channel formation at the pore-scale and progressive high porosity pathway dissolution at the core-scale with the more heterogeneous rock having dissolution progressing along direction of flow more quickly. Additionally, upon analysis of the high-resolution macro core images it was found that the dissolution pathways contained a distinct microstructure that was not visible at the resolution of the medical CT, where the reactive fluid had not completely dissolved the internal pore-structure. Flow was modelled in connected pathways, the flow streamlines were traced and streamline density for each voxel was calculated. It was found that the streamline density was highest in the most well-connected pathways and that density increased with increasing heterogeneity as the number of connected pathways decreased and flow was consolidated along fewer pathways. This work represents the first study of scale dependency using reservoir condition 4D X-ray tomography and provides insight into the mechanisms that control local reaction rates at multiple scales.

FINAL REPORT: Mechanistically-Base Field Scale Models of Uranium Biogeochemistry from Upscaling Pore-Scale Experiments and Models

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wood, Brian D.

2013-11-04

Biogeochemical reactive transport processes in the subsurface environment are important to many contemporary environmental issues of significance to DOE. Quantification of risks and impacts associated with environmental management options, and design of remediation systems where needed, require that we have at our disposal reliable predictive tools (usually in the form of numerical simulation models). However, it is well known that even the most sophisticated reactive transport models available today have poor predictive power, particularly when applied at the field scale. Although the lack of predictive ability is associated in part with our inability to characterize the subsurface and limitations inmore » computational power, significant advances have been made in both of these areas in recent decades and can be expected to continue. In this research, we examined the upscaling (pore to Darcy and Darcy to field) the problem of bioremediation via biofilms in porous media. The principle idea was to start with a conceptual description of the bioremediation process at the pore scale, and apply upscaling methods to formally develop the appropriate upscaled model at the so-called Darcy scale. The purpose was to determine (1) what forms the upscaled models would take, and (2) how one might parameterize such upscaled models for applications to bioremediation in the field. We were able to effectively upscale the bioremediation process to explain how the pore-scale phenomena were linked to the field scale. The end product of this research was to produce a set of upscaled models that could be used to help predict field-scale bioremediation. These models were mechanistic, in the sense that they directly incorporated pore-scale information, but upscaled so that only the essential features of the process were needed to predict the effective parameters that appear in the model. In this way, a direct link between the microscale and the field scale was made, but the upscaling process helped inform potential users of the model what kinds of information would be needed to accurately characterize the system.« less

Unifying Pore Network Modeling, Continuous Time Random Walk Theory and Experiment - Accomplishments and Future Directions

NASA Astrophysics Data System (ADS)

Bijeljic, B.

2008-05-01

This talk will describe and highlight the advantages offered by a methodology that unifies pore network modeling, CTRW theory and experiment in description of solute dispersion in porous media. Solute transport in a porous medium is characterized by the interplay of advection and diffusion (described by Peclet number, Pe) that cause spreading of solute particles. This spreading is traditionally described by dispersion coefficients, D, defined by σ 2 = 2Dt, where σ 2 is the variance of the solute position and t is the time. Using a pore-scale network model based on particle tracking, the rich Peclet- number dependence of dispersion coefficient is predicted from first principles and is shown to compare well with experimental data for restricted diffusion, transition, power-law and mechanical dispersion regimes in the asymptotic limit. In the asymptotic limit D is constant and can be used in an averaged advection-dispersion equation. However, it is highly important to recognize that, until the velocity field is fully sampled, the particle transport is non-Gaussian and D possesses temporal or spatial variation. Furthermore, temporal probability density functions (PDF) of tracer particles are studied in pore networks and an excellent agreement for the spectrum of transition times for particles from pore to pore is obtained between network model results and CTRW theory. Based on the truncated power-law interpretation of PDF-s, the physical origin of the power-law scaling of dispersion coefficient vs. Peclet number has been explained for unconsolidated porous media, sands and a number of sandstones, arriving at the same conclusion from numerical network modelling, analytic CTRW theory and experiment. Future directions for further applications of the methodology presented are discussed in relation to the scale- dependent solute dispersion and reactive transport. Significance of pre-asymptotic dispersion in porous media is addressed from pore-scale upwards and the impact of heterogeneity is discussed. The length traveled by solute plumes before Gaussian behaviour is reached increases with an increase in heterogeneity and/or Pe. This opens up the question on the nature of dispersion in natural systems where the heterogeneities at the larger scales will profoundly increase the range of velocities in the aquifer, thus considerably delaying the asymptotic approach to Gaussian behaviour. As a consequence, the asymptotic behaviour might not be reached at the field scale.

Multiscale modelling of dual-porosity porous media; a computational pore-scale study for flow and solute transport

NASA Astrophysics Data System (ADS)

de Vries, Enno T.; Raoof, Amir; van Genuchten, Martinus Th.

2017-07-01

Many environmental and agricultural applications involve the transport of water and dissolved constituents through aggregated soil profiles, or porous media that are structured, fractured or macroporous in other ways. During the past several decades, various process-based macroscopic models have been used to simulate contaminant transport in such media. Many of these models consider advective-dispersive transport through relatively large inter-aggregate pore domains, while exchange with the smaller intra-aggregate pores is assumed to be controlled by diffusion. Exchange of solute between the two domains is often represented using a first-order mass transfer coefficient, which is commonly obtained by fitting to observed data. This study aims to understand and quantify the solute exchange term by applying a dual-porosity pore-scale network model to relatively large domains, and analysing the pore-scale results in terms of the classical dual-porosity (mobile-immobile) transport formulation. We examined the effects of key parameters (notably aggregate porosity and aggregate permeability) on the main dual-porosity model parameters, i.e., the mobile water fraction (ϕm) and the mass transfer coefficient (α). Results were obtained for a wide range of aggregate porosities (between 0.082 and 0.700). The effect of aggregate permeability was explored by varying pore throat sizes within the aggregates. Solute breakthrough curves (BTCs) obtained with the pore-scale network model at several locations along the domain were analysed using analytical solutions of the dual-porosity model to obtain estimates of ϕm and α. An increase in aggregate porosity was found to decrease ϕm and increase α, leading to considerable tailing in the BTCs. Changes in the aggregate pore throat size affected the relative flow velocity between the intra- and inter-aggregate domains. Higher flow velocities within the aggregates caused a change in the transport regime from diffusion dominated to more advection dominated. This change increased the exchange rate of solutes between the mobile and immobile domains, with a related increase in the value of the mass transfer coefficient and less tailing in the BTCs.

Restore the change process of Longmaxi shale pore structure during the diagenesis by the potassium and the magnesium

NASA Astrophysics Data System (ADS)

Zhou, W.

2016-12-01

The pore structure of Longmaxi shale was changing during the diagenetic process, mainly caused by the illitization and serpentinzation. The evolution of shale pore structure mainly relates to the element migration. Based on the result of electron microprobe analyser (EMPA), it is possible to find the distribution of element in shale directly and to distinguish the destroyed primary pore structure as element will remain in the migration way. The migration of potassium in Longmaxi shale mainly happened during early diagenesis phase to middle diagenesis phase (Geothermal temperature: 60°-140°). During the illitization, potassium mainly came from potassium feldspar, migrated though the connected pore structure and reacted with smectite. Illite and illite/smectite in Longmaxi shale distribute continuously in 10micron-level flocculent formation, which means that primary connective pore structure in Longmaxi shale has a same scale. The concentration of potassium has an obvious gradient that potassium content in middle of flocculation of Illite/smectite is about 6.8% and 4.8% in the boundary parts (Fig.). In addition, as SiO2 was generated during the illitization, which makes Longmaxi shale very compacted. The migration of magnesium in Longmaxi shale happened during low temperature serpentinization (Geothermal temperature: 140°-350°). Magnesium mainly came from dolomite and migrated in primary pores. According to the result of EMPA, it can be recognized that the migration path of magnesium is much simpler than potassium, which is caused as serpentinization do not have much reaction with clay minerals around (Fig.). Serpentine jams the primary pores of Longmaxi shale too. As reaction temperature of serpentinization is higher than illitization and serpentine is inserts in illite/smectite, the formation process of Longmaxi shale pore structure can be mainly divided into two phases: geothermal temperature˜140° and˜140°.

Spatial and temporal distribution of pore gas concentrations during mainstream large-scale trough composting in China.

PubMed

Zeng, Jianfei; Shen, Xiuli; Sun, Xiaoxi; Liu, Ning; Han, Lujia; Huang, Guangqun

2018-05-01

With the advantages of high treatment capacity and low operational cost, large-scale trough composting has become one of the mainstream composting patterns in composting plants in China. This study measured concentrations of O 2 , CO 2 , CH 4 and NH 3 on-site to investigate the spatial and temporal distribution of pore gas concentrations during mainstream large-scale trough composting in China. The results showed that the temperature in the center of the pile was obviously higher than that in the side of the pile. Pore O 2 concentration rapidly decreased and maintained <5% (in volume) for 38 days or more in both the center and side of the pile and effective O 2 diffusion occurred at most in every two contiguous layers. Pore CO 2 and CH 4 concentrations at each measurement point were positively correlated (0.436 ≤ r ≤ 0.570, P < 0.01) and the concentrations in the side of the pile were obviously lower than those in the center. The top layer exhibited highest pore O 2 concentration and lowest CO 2 and CH 4 concentrations, and the bottom layer was on the contrary. No significant differences in pore NH 3 concentrations between different layers or between different measurement points in the same layer were found. Therefore, mixing the center and the side of the pile when mechanical turning and adjusting the height of the pile according to the physical properties of bulking agents are suggested to optimize the oxygen distribution and promote the composting process during large-scale trough composting when the pile was naturally aerated, which will contribute to improving the current undesirable atmosphere environment in China. Copyright © 2018 Elsevier Ltd. All rights reserved.

A Generalized Hybrid Multiscale Modeling Approach for Flow and Reactive Transport in Porous Media

NASA Astrophysics Data System (ADS)

Yang, X.; Meng, X.; Tang, Y. H.; Guo, Z.; Karniadakis, G. E.

2017-12-01

Using emerging understanding of biological and environmental processes at fundamental scales to advance predictions of the larger system behavior requires the development of multiscale approaches, and there is strong interest in coupling models at different scales together in a hybrid multiscale simulation framework. A limited number of hybrid multiscale simulation methods have been developed for subsurface applications, mostly using application-specific approaches for model coupling. The proposed generalized hybrid multiscale approach is designed with minimal intrusiveness to the at-scale simulators (pre-selected) and provides a set of lightweight C++ scripts to manage a complex multiscale workflow utilizing a concurrent coupling approach. The workflow includes at-scale simulators (using the lattice-Boltzmann method, LBM, at the pore and Darcy scale, respectively), scripts for boundary treatment (coupling and kriging), and a multiscale universal interface (MUI) for data exchange. The current study aims to apply the generalized hybrid multiscale modeling approach to couple pore- and Darcy-scale models for flow and mixing-controlled reaction with precipitation/dissolution in heterogeneous porous media. The model domain is packed heterogeneously that the mixing front geometry is more complex and not known a priori. To address those challenges, the generalized hybrid multiscale modeling approach is further developed to 1) adaptively define the locations of pore-scale subdomains, 2) provide a suite of physical boundary coupling schemes and 3) consider the dynamic change of the pore structures due to mineral precipitation/dissolution. The results are validated and evaluated by comparing with single-scale simulations in terms of velocities, reactive concentrations and computing cost.

Laboratory measurements of shock propagation through spherical cavities in an optically accessible polymer.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chojnicki, Kirsten; Cooper, Marcia A.; Guo, Shuyue

Pore-scale aperture effects on flow in pore networks was studied in the laboratory to provide a parameterization for use in transport models. Four cases were considered: regular and irregular pillar/pore alignment with and without an aperture. The velocity field of each case was measured and simulated, providing quantitatively comparable results. Two aperture effect parameterizations were considered: permeability and transmission. Permeability values varied by an order of magnitude between the cases with and without apertures. However, transmission did not correlate with permeability. Despite having much greater permeability the regular aperture case permitted less transmission than the regular case. Moreover, both irregularmore » cases had greater transmission than the regular cases, a difference not supported by the permeabilities. Overall, these findings suggest that pore-scale aperture effects on flow though a pore-network may not be adequately captured by properties such as permeability for applications that are interested in determining particle transport volume and timing.« less

Transport upscaling from pore- to Darcy-scale: Incorporating pore-scale Berea sandstone Lagrangian velocity statistics into a Darcy-scale transport CTRW model

NASA Astrophysics Data System (ADS)

Puyguiraud, Alexandre; Dentz, Marco; Gouze, Philippe

2017-04-01

For the past several years a lot of attention has been given to pore-scale flow in order to understand and model transport, mixing and reaction in porous media. Nevertheless we believe that an accurate study of spatial and temporal evolution of velocities could bring important additional information for the upscaling from pore to higher scales. To gather these pieces of information, we perform Stokes flow simulations on pore-scale digitized images of a Berea sandstone core. First, micro-tomography (XRMT) imaging and segmentation processes allow us to obtain 3D black and white images of the sample [1]. Then we used an OpenFoam solver to perform the Stokes flow simulations mentioned above, which gives us the velocities at the interfaces of a cubic mesh. Subsequently, we use a particle streamline reconstruction technique which uses the Eulerian velocity field previously obtained. This technique, based on a modified Pollock algorithm [2], enables us to make particle tracking simulations on the digitized sample. In order to build a stochastic pore-scale transport model, we analyze the Lagrangian velocity series in two different ways. First we investigate the velocity evolution by sampling isochronically (t-Lagrangian), and by studying its statistical properties in terms of one- and two-points statistics. Intermittent patterns can be observed. These are due to the persistance of low velocities over a characteristic space length. Other results are investigated, such as correlation functions and velocity PDFs, which permit us to study more deeply this persistence in the velocities and to compute the correlation times. However, with the second approach, doing these same analysis in space by computing the velocities equidistantly, enables us to remove the intermittency shown in the temporal evolution and to model these velocity series as a Markov process. This renders the stochastic particle dynamics into a CTRW [3]. [1] Gjetvaj, F., A. Russian, P. Gouze, and M. Dentz (2015), Dual control of flow field heterogeneity and immobile porosity on non-Fickian transport in Berea sandstone, Water Resour. Res., 51, 8273-8293, doi:10.1002/2015WR017645. [2] Mostaghimi, P., Bijeljic, B., Blunt, M. (2012). Simulation of Flow and Dispersion on Pore-Space Images. Society of Petroleum Engineers. doi:10.2118/135261-PA. [3] Dentz, M., P. K. Kang, A. Comolli, T. Le Borgne, and D. R. Lester, Continuous time random walks for the evolution of Lagrangian velocities, Phys. Rev. Fluids, 2016. Keywords: Porescale, particle tracking, transport, Lagrangian velocity, ergodicity, Markovianity, continuous time random walks, upscaling.

Evaluation of Micro- and Macro-Scale Petrophysical Characteristics of Lower Cretaceous Sandstone with Flow Modeling in µ-CT Imaged Geometry

NASA Astrophysics Data System (ADS)

Katsman, R.; Haruzi, P.; Waldmann, N.; Halisch, M.

2017-12-01

In this study petrophysical characteristics of rock samples from 3 successive outcrop layers of Hatira Formation Lower Cretaceous Sandstone in northen Israel were evaluated at micro- and macro-scales. The study was carried out by two complementary methods: using conventional experimental measurements of porosity, pore size distribution and permeability; and using a 3D µCT imaging and modeling of signle-phase flow in the real micro-scale sample geometry. The workfow included µ-CT scanning, image processing, image segmentation, and image analyses of pore network, followed by fluid flow simulations at a pore-scale. Upscaling the results of the micro-scale flow simulations yielded a macroscopic permeabilty tensor. Comparison of the upscaled and the experimentally measured rock properties demonstrated a reasonable agreement. In addition, geometrical (pore size distribution, surface area and tortuosity) and topological (Euler characteristic) characteristics of the grains and of the pore network were evaluated at a micro-scale. Statistical analyses of the samples for estimation of anisotropy and inhomogeneity of the porous media were conducted and the results agree with anisotropy and inhomogeneity of the upscaled permeabilty tensor. Isotropic pore orientation of the primary inter-granular porosity was identified in all three samples, whereas the characteristics of the secondary porosity were affected by precipitated cement and clay matrix within the primary pore network. Results of this study provide micro- and macro-scale characteristics of the Lower Cretaceous sandstone that is used in different places over the world as a reservoir for petroleum production and sequestration.

On the predictivity of pore-scale simulations: Estimating uncertainties with multilevel Monte Carlo

NASA Astrophysics Data System (ADS)

Icardi, Matteo; Boccardo, Gianluca; Tempone, Raúl

2016-09-01

A fast method with tunable accuracy is proposed to estimate errors and uncertainties in pore-scale and Digital Rock Physics (DRP) problems. The overall predictivity of these studies can be, in fact, hindered by many factors including sample heterogeneity, computational and imaging limitations, model inadequacy and not perfectly known physical parameters. The typical objective of pore-scale studies is the estimation of macroscopic effective parameters such as permeability, effective diffusivity and hydrodynamic dispersion. However, these are often non-deterministic quantities (i.e., results obtained for specific pore-scale sample and setup are not totally reproducible by another ;equivalent; sample and setup). The stochastic nature can arise due to the multi-scale heterogeneity, the computational and experimental limitations in considering large samples, and the complexity of the physical models. These approximations, in fact, introduce an error that, being dependent on a large number of complex factors, can be modeled as random. We propose a general simulation tool, based on multilevel Monte Carlo, that can reduce drastically the computational cost needed for computing accurate statistics of effective parameters and other quantities of interest, under any of these random errors. This is, to our knowledge, the first attempt to include Uncertainty Quantification (UQ) in pore-scale physics and simulation. The method can also provide estimates of the discretization error and it is tested on three-dimensional transport problems in heterogeneous materials, where the sampling procedure is done by generation algorithms able to reproduce realistic consolidated and unconsolidated random sphere and ellipsoid packings and arrangements. A totally automatic workflow is developed in an open-source code [1], that include rigid body physics and random packing algorithms, unstructured mesh discretization, finite volume solvers, extrapolation and post-processing techniques. The proposed method can be efficiently used in many porous media applications for problems such as stochastic homogenization/upscaling, propagation of uncertainty from microscopic fluid and rock properties to macro-scale parameters, robust estimation of Representative Elementary Volume size for arbitrary physics.

Nanoscale Pore Imaging and Pore Scale Fluid Flow Modeling in Chalk

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tomutsa, Liviu; Silin, Dmitriy

2004-08-19

For many rocks of high economic interest such as chalk, diatomite, tight gas sands or coal, nanometer scale resolution is needed to resolve the 3D-pore structure, which controls the flow and trapping of fluids in the rocks. Such resolutions cannot be achieved with existing tomographic technologies. A new 3D imaging method, based on serial sectioning and using the Focused Ion Beam (FIB) technology has been developed. FIB allows for the milling of layers as thin as 10 nanometers by using accelerated Ga+ ions to sputter atoms from the sample surface. After each milling step, as a new surface is exposed,more » a 2D image of this surface is generated. Next, the 2D images are stacked to reconstruct the 3D pore or grain structure. Resolutions as high as 10 nm are achievable using such a technique. A new robust method of pore-scale fluid flow modeling has been developed and applied to sandstone and chalk samples. The method uses direct morphological analysis of the pore space to characterize the petrophysical properties of diverse formations. Not only petrophysical properties (porosity, permeability, relative permeability and capillary pressures) can be computed but also flow processes, such as those encountered in various IOR approaches, can be simulated. Petrophysical properties computed with the new method using the new FIB data will be presented. Present study is a part of the development of an Electronic Core Laboratory at LBNL/UCB.« less

Effects of pore-scale dispersion, degree of heterogeneity, sampling size, and source volume on the concentration moments of conservative solutes in heterogeneous formations

Treesearch

Daniele Tonina; Alberto Bellin

2008-01-01

Pore-scale dispersion (PSD), aquifer heterogeneity, sampling volume, and source size influence solute concentrations of conservative tracers transported in heterogeneous porous formations. In this work, we developed a new set of analytical solutions for the concentration ensemble mean, variance, and coefficient of variation (CV), which consider the effects of all these...

Pore-pressure sensitivities to dynamic strains: observations in active tectonic regions

USGS Publications Warehouse

Barbour, Andrew J.

2015-01-01

Triggered seismicity arising from dynamic stresses is often explained by the Mohr-Coulomb failure criterion, where elevated pore pressures reduce the effective strength of faults in fluid-saturated rock. The seismic response of a fluid-rock system naturally depends on its hydro-mechanical properties, but accurately assessing how pore-fluid pressure responds to applied stress over large scales in situ remains a challenging task; hence, spatial variations in response are not well understood, especially around active faults. Here I analyze previously unutilized records of dynamic strain and pore-pressure from regional and teleseismic earthquakes at Plate Boundary Observatory (PBO) stations from 2006 through 2012 to investigate variations in response along the Pacific/North American tectonic plate boundary. I find robust scaling-response coefficients between excess pore pressure and dynamic strain at each station that are spatially correlated: around the San Andreas and San Jacinto fault systems, the response is lowest in regions of the crust undergoing the highest rates of secular shear strain. PBO stations in the Parkfield instrument cluster are at comparable distances to the San Andreas fault (SAF), and spatial variations there follow patterns in dextral creep rates along the fault, with the highest response in the actively creeping section, which is consistent with a narrowing zone of strain accumulation seen in geodetic velocity profiles. At stations in the San Juan Bautista (SJB) and Anza instrument clusters, the response depends non-linearly on the inverse fault-perpendicular distance, with the response decreasing towards the fault; the SJB cluster is at the northern transition from creeping-to-locked behavior along the SAF, where creep rates are at moderate to low levels, and the Anza cluster is around the San Jacinto fault, where to date there have been no statistically significant creep rates observed at the surface. These results suggest that the strength of the pore pressure response in fluid-saturated rock near active faults is controlled by shear strain accumulation associated with tectonic loading, which implies a strong feedback between fault strength and permeability: dynamic triggering susceptibilities may vary in space and also in time.

The effects of wettability and trapping on relationships between interfacial area, capillary pressure and saturation in porous media: A pore-scale network modeling approach

NASA Astrophysics Data System (ADS)

Raeesi, Behrooz; Piri, Mohammad

2009-10-01

SummaryWe use a three-dimensional mixed-wet random pore-scale network model to investigate the impact of wettability and trapping on the relationship between interfacial area, capillary pressure and saturation in two-phase drainage and imbibition processes. The model is a three-dimensional network of interconnected pores and throats of various geometrical shapes. It allows multiple phases to be present in each capillary element in wetting and spreading layers, as well as occupying the center of the pore space. Two different random networks that represent the pore space in Berea and a Saudi Arabia reservoir sandstone are used in this study. We allow the wettability of the rock surfaces contacted by oil to alter after primary drainage. The model takes into account both contact angle and trapping hystereses. We model primary oil drainage and water flooding for mixed-wet conditions, and secondary oil injection for a water-wet system. The total interfacial area for pores and throats are calculated when the system is at capillary equilibrium. They include contributions from the arc menisci (AMs) between the bulk and corner fluids, and from the main terminal menisci (MTMs) between different bulk fluids. We investigate hysteresis in these relationships by performing water injection into systems of varying wettability and initial water saturation. We show that trapping and contact angle hystereses significantly affect the interfacial area. In a strongly water-wet system, a sharp increase is observed at the beginning of water flood, which shifts the area to a higher level than primary drainage. As we change the wettability of the system from strongly water-wet to strongly oil-wet, the trapped oil saturation decreases significantly. Starting water flood from intermediate water saturations, greater than the irreducible water saturation, can also affect the non-wetting phase entrapment, resulting in different interfacial area behaviors. This can increase the interfacial area significantly in oil-wet systems. A qualitative comparison of our results with the experimental data available in literature for glass beads shows, with some expected differences, an encouraging agreement. Also, our results agree well with those generated by the previously developed models.

Percolation properties of 3-D multiscale pore networks: how connectivity controls soil filtration processes

NASA Astrophysics Data System (ADS)

Perrier, E. M. A.; Bird, N. R. A.; Rieutord, T. B.

2010-04-01

Quantifying the connectivity of pore networks is a key issue not only for modelling fluid flow and solute transport in porous media but also for assessing the ability of soil ecosystems to filter bacteria, viruses and any type of living microorganisms as well inert particles which pose a contamination risk. Straining is the main mechanical component of filtration processes: it is due to size effects, when a given soil retains a conveyed entity larger than the pores through which it is attempting to pass. We postulate that the range of sizes of entities which can be trapped inside soils has to be associated with the large range of scales involved in natural soil structures and that information on the pore size distribution has to be complemented by information on a Critical Filtration Size (CFS) delimiting the transition between percolating and non percolating regimes in multiscale pore networks. We show that the mass fractal dimensions which are classically used in soil science to quantify scaling laws in observed pore size distributions can also be used to build 3-D multiscale models of pore networks exhibiting such a critical transition. We extend to the 3-D case a new theoretical approach recently developed to address the connectivity of 2-D fractal networks (Bird and Perrier, 2009). Theoretical arguments based on renormalisation functions provide insight into multi-scale connectivity and a first estimation of CFS. Numerical experiments on 3-D prefractal media confirm the qualitative theory. These results open the way towards a new methodology to estimate soil filtration efficiency from the construction of soil structural models to be calibrated on available multiscale data.

Percolation properties of 3-D multiscale pore networks: how connectivity controls soil filtration processes

NASA Astrophysics Data System (ADS)

Perrier, E. M. A.; Bird, N. R. A.; Rieutord, T. B.

2010-10-01

Quantifying the connectivity of pore networks is a key issue not only for modelling fluid flow and solute transport in porous media but also for assessing the ability of soil ecosystems to filter bacteria, viruses and any type of living microorganisms as well inert particles which pose a contamination risk. Straining is the main mechanical component of filtration processes: it is due to size effects, when a given soil retains a conveyed entity larger than the pores through which it is attempting to pass. We postulate that the range of sizes of entities which can be trapped inside soils has to be associated with the large range of scales involved in natural soil structures and that information on the pore size distribution has to be complemented by information on a critical filtration size (CFS) delimiting the transition between percolating and non percolating regimes in multiscale pore networks. We show that the mass fractal dimensions which are classically used in soil science to quantify scaling laws in observed pore size distributions can also be used to build 3-D multiscale models of pore networks exhibiting such a critical transition. We extend to the 3-D case a new theoretical approach recently developed to address the connectivity of 2-D fractal networks (Bird and Perrier, 2009). Theoretical arguments based on renormalisation functions provide insight into multi-scale connectivity and a first estimation of CFS. Numerical experiments on 3-D prefractal media confirm the qualitative theory. These results open the way towards a new methodology to estimate soil filtration efficiency from the construction of soil structural models to be calibrated on available multiscale data.

Pore-scale investigation on stress-dependent characteristics of granular packs and the impact of pore deformation on fluid distribution

DOE PAGES

Yoon, Hongkyu; Klise, Katherine A.; Torrealba, Victor A.; ...

2015-05-25

Understanding the effect of changing stress conditions on multiphase flow in porous media is of fundamental importance for many subsurface activities including enhanced oil recovery, water drawdown from aquifers, soil confinement, and geologic carbon storage. Geomechanical properties of complex porous systems are dynamically linked to flow conditions, but their feedback relationship is often oversimplified due to the difficulty of representing pore-scale stress deformation and multiphase flow characteristics in high fidelity. In this work, we performed pore-scale experiments of single- and multiphase flow through bead packs at different confining pressure conditions to elucidate compaction-dependent characteristics of granular packs and their impactmore » on fluid flow. A series of drainage and imbibition cycles were conducted on a water-wet, soda-lime glass bead pack under varying confining stress conditions. Simultaneously, X-ray micro-CT was used to visualize and quantify the degree of deformation and fluid distribution corresponding with each stress condition and injection cycle. Micro-CT images were segmented using a gradient-based method to identify fluids (e.g., oil and water), and solid phase redistribution throughout the different experimental stages. Changes in porosity, tortuosity, and specific surface area were quantified as a function of applied confining pressure. Results demonstrate varying degrees of sensitivity of these properties to confining pressure, which suggests that caution must be taken when considering scalability of these properties for practical modeling purposes. Changes in capillary number with confining pressure are attributed to the increase in pore velocity as a result of pore contraction. Furthermore, this increase in pore velocity was found to have a marginal impact on average phase trapping at different confining pressures.« less

Multi-Scale Porous Ultra High Temperature Ceramics

DTIC Science & Technology

2015-01-08

different techniques: replica, particle stabilized foams, ice templating (freeze casting) and partial sintering. The pore morphology (closed-bubble...the porosity, pore size, shape and morphology . X-Ray Tomography was used to study their 3D microstructure. The 3D microstructures captured with...four different techniques: replica, particle stabilized foams, ice templating (freeze casting) and partial sintering. The pore morphology (closed-bubble

Quantification of Soil Pore Structure Based on Minkowski-Functions

NASA Astrophysics Data System (ADS)

Vogel, H.; Weller, U.; Schlüter, S.

2009-05-01

The porous structure in soils and other geologic media is typically a complex 3-dimensional object. Most of the physical material properties including mechanical and hydraulic characteristics are immediately linked to this structure which can be directly observed using non-invasive techniques as e.g. X-ray tomography. It is an old dream and still a formidable challenge to related structural features of porous media to their physical properties. In this contribution we present a scale-invariant concept to quantify pore structure based on a limited set of meaningful morphological functions. They are based on d+1 Minkowski functionals as defined for d-dimensional bodies. These basic quantities are determined as a function of pore size obtained by filter procedures using mathematical morphology. The resulting Minkowski functions provide valuable information on pore size, pore surface area and pore topology having the potential to be linked to physical properties. The theoretical background and the related algorithms are presented and the approach is demonstrated for the structure of an arable topsoil obtained by X-ray micro tomography. We also discuss the fundamental problem of limited resolution which is critical for any attempt to quantify structural features at any scale.

Statistical scaling of pore-scale Lagrangian velocities in natural porous media.

PubMed

Siena, M; Guadagnini, A; Riva, M; Bijeljic, B; Pereira Nunes, J P; Blunt, M J

2014-08-01

We investigate the scaling behavior of sample statistics of pore-scale Lagrangian velocities in two different rock samples, Bentheimer sandstone and Estaillades limestone. The samples are imaged using x-ray computer tomography with micron-scale resolution. The scaling analysis relies on the study of the way qth-order sample structure functions (statistical moments of order q of absolute increments) of Lagrangian velocities depend on separation distances, or lags, traveled along the mean flow direction. In the sandstone block, sample structure functions of all orders exhibit a power-law scaling within a clearly identifiable intermediate range of lags. Sample structure functions associated with the limestone block display two diverse power-law regimes, which we infer to be related to two overlapping spatially correlated structures. In both rocks and for all orders q, we observe linear relationships between logarithmic structure functions of successive orders at all lags (a phenomenon that is typically known as extended power scaling, or extended self-similarity). The scaling behavior of Lagrangian velocities is compared with the one exhibited by porosity and specific surface area, which constitute two key pore-scale geometric observables. The statistical scaling of the local velocity field reflects the behavior of these geometric observables, with the occurrence of power-law-scaling regimes within the same range of lags for sample structure functions of Lagrangian velocity, porosity, and specific surface area.

Membrane properties change in fine-pore aeration diffusers: full-scale variations of transfer efficiency and headloss.

PubMed

Rosso, Diego; Libra, Judy A; Wiehe, Wolfgang; Stenstrom, Michael K

2008-05-01

Fine-pore diffusers are the most common aeration system in municipal wastewater treatment. Punched polymeric membranes are often used in fine-pore aeration due to their advantageous initial performance. These membranes are subject to fouling and scaling, resulting in increased headloss and reduced oxygen transfer efficiency, both contributing to increased plant energy costs. This paper describes and discusses the change in material properties for polymeric fine-pore diffusers, comparing new and used membranes. Three different diffuser technologies were tested and sample diffusers from two wastewater treatment facilities were analysed. The polymeric membranes analysed in this paper were composed of ethylene-propylene-diene monomer (EPDM), polyurethane, and silicon. Transfer efficiency is usually lower with longer times in operation, as older, dilated orifices produce larger bubbles, which are unfavourable to mass transfer. At the same time, headloss increases with time in operation, since membranes increase in rigidity and hardness, and fouling and scaling phenomena occur at the orifice opening. Change in polymer properties and laboratory test results correlate with the decrease in oxygen transfer efficiency.

Development of a New Analog Test System Capable of Modeling Tectonic Deformation Incorporating the Effects of Pore Fluid Pressure

NASA Astrophysics Data System (ADS)

Zhang, M.; Nakajima, H.; Takeda, M.; Aung, T. T.

2005-12-01

Understanding and predicting the tectonic deformation within geologic strata has been a very important research subject in many fields such as structural geology and petroleum geology. In recent years, such research has also become a fundamental necessity for the assessment of active fault migration, site selection for geological disposal of radioactive nuclear waste and exploration for methane hydrate. Although analog modeling techniques have played an important role in the elucidation of the tectonic deformation mechanisms, traditional approaches have typically used dry materials and ignored the effects of pore fluid pressure. In order for analog models to properly depict the tectonic deformation of the targeted, large-prototype system within a small laboratory-scale configuration, physical properties of the models, including geometry, force, and time, must be correctly scaled. Model materials representing brittle rock behavior require an internal friction identical to the prototype rock and virtually zero cohesion. Granular materials such as sand, glass beads, or steel beads of dry condition have been preferably used for this reason in addition to their availability and ease of handling. Modeling protocols for dry granular materials have been well established but such model tests cannot account for the pore fluid effects. Although the concept of effective stress has long been recognized and the role of pore-fluid pressure in tectonic deformation processes is evident, there have been few analog model studies that consider the effects of pore fluid movement. Some new applications require a thorough understanding of the coupled deformation and fluid flow processes within the strata. Taking the field of waste management as an example, deep geological disposal of radioactive waste has been thought to be an appropriate methodology for the safe isolation of the wastes from the human environment until the toxicity of the wastes decays to non-hazardous levels. For the deep geological disposal concept, besides containing the wastes with engineering methods such as the glassification of the radioactive wastes, the geological formation itself is expected to serve as a natural barrier that retards migration of radionuclides. To evaluate the long-term safety of deep geological disposal, a better understanding of the fate and transport of radionuclides in a geologically heterogeneous environment is necessary. To meet such requirements, a new analog test sandbox model system was developed. This model system allows the pore fluid flows to be controlled during the model tests and permits the study of flow and transport phenomena in the deformed heterogeneous model. One- or two-dimensional fluid flow is controlled using a side-wall piston. Deformation processes can be observed through a transparent front panel, and pore fluid movement can be also visualized using a color tracer. In this study, the scaling requirements for analog modeling, including pore water pressure, are discussed based on the theory of dimensional analysis, supplemented by data from a series of laboratory shear tests, and a detailed description of the model system. Preliminary experimental results are presented.

Observations of wave-induced pore pressure gradients and bed level response on a surf zone sandbar

NASA Astrophysics Data System (ADS)

Anderson, Dylan; Cox, Dan; Mieras, Ryan; Puleo, Jack A.; Hsu, Tian-Jian

2017-06-01

Horizontal and vertical pressure gradients may be important physical mechanisms contributing to onshore sediment transport beneath steep, near-breaking waves in the surf zone. A barred beach was constructed in a large-scale laboratory wave flume with a fixed profile containing a mobile sediment layer on the crest of the sandbar. Horizontal and vertical pore pressure gradients were obtained by finite differences of measurements from an array of pressure transducers buried within the upper several centimeters of the bed. Colocated observations of erosion depth were made during asymmetric wave trials with wave heights between 0.10 and 0.98 m, consistently resulting in onshore sheet flow sediment transport. The pore pressure gradient vector within the bed exhibited temporal rotations during each wave cycle, directed predominantly upward under the trough and then rapidly rotating onshore and downward as the wavefront passed. The magnitude of the pore pressure gradient during each phase of rotation was correlated with local wave steepness and relative depth. Momentary bed failures as deep as 20 grain diameters were coincident with sharp increases in the onshore-directed pore pressure gradients, but occurred at horizontal pressure gradients less than theoretical critical values for initiation of the motion for compact beds. An expression combining the effects of both horizontal and vertical pore pressure gradients with bed shear stress and soil stability is used to determine that failure of the bed is initiated at nonnegligible values of both forces.