Dynamic Inhibition of Calcite Dissolution in Flowing Acidic Pb2+ Solutions.

Reactions of mineral surfaces with dissolved metal ions at far-from-equilibrium conditions can deviate significantly from those in near-equilibrium systems due to steep concentration gradients, ion-surface interactions, and reactant transport effects that can lead to emergent behavior. We explored the effect of dissolved Pb2+ on the dissolution rate and topographic evolution of calcite (104) surfaces under far-from-equilibrium acidic conditions (pH 3.7) in a confined single-pass laminar-flow geometry. Operando measurements by digital holographic microscopy were conducted over a range of Pb2+ concentrations ([Pb2+] = 0 to 5 × 10-2 M) and flow velocities (v = 1.67-53.3 mm s-1). Calcite (104) surface dissolution rates decreased with increasing [Pb2+]. The inhibition of dissolution and the emergence of unique topographic features, including micropyramids, variable etch pit shapes, and larger scale topographic patterns, became increasingly apparent at [Pb2+] ≥ 5 × 10-3 M. A better understanding of such dynamic reactivity could be crucial for constructing accurate models of geochemical transport in aqueous carbonate systems.

The clogging effect of nanoscale zero valent iron corrosion in bicarbonate anaerobic water on porous media: A real-time pore-scale visualization.

The corrosion-induced permeability changes of nanoscale zero-valent iron (NZVI) are one of the crucial factors constraining the successful application of NZVI in the remediation of contaminated groundwater. It is of great significance to study the dynamic evolution of corrosion products of NZVI after NZVI is injected into porous media and its influence on pore plugging effect from the pore scale. Micro computed tomography (Micro-CT) imaging technology, mineralogical characterization and theoretical calculations were used to understand the details of NZVI corrosion plugging porous media at the pore scale. This study reveals the factors of NZVI corrosion plugging porous media, namely, gas production (H2) in the early and middle stages of corrosion (before 90 days) and solid phase changes (NZVI volume increase and migration) in the later stages (after 90 days). The permeability loss rate of the porous media was 66.8%, 87.3%, 79.4%, and 53.6% at the corrosion times of 30, 60, 90, and 120 days, respectively. After 90 d of corrosion, the particle size of NZVI increases by 7.9%, and the secondary minerals formed by corrosion are mainly Fe3O4/γ-Fe2O3 and FeOOH. In addition, this study also found that the migration of NZVI after 90 d was due to its corrosion reducing the magnetic attraction between particles, dissociating into smaller particles or agglomerates under the action of fluid dynamics, resulting in its redistribution in the porous medium and causing blockage. This study clarifies that NZVI corrosion plays a vital influence in affecting the permeability and clogging of porous media, providing valuable guidance for optimizing NZVI-based remediation technologies.

The effect of root hairs on root water uptake is determined by root-soil contact and root hair shrinkage.

The effect of root hairs on water uptake remains controversial. In particular, the key root hair and soil parameters that determine their importance have been elusive. We grew maize plants (Zea mays) in microcosms and scanned them using synchrotron-based X-ray computed microtomography. By means of image-based modelling, we investigated the parameters determining the effectiveness of root hairs in root water uptake. We explicitly accounted for rhizosphere features (e.g. root-soil contact and pore structure) and took root hair shrinkage of dehydrated root hairs into consideration. Our model suggests that > 85% of the variance in root water uptake is explained by the hair-induced increase in root-soil contact. In dry soil conditions, root hair shrinkage reduces the impact of hairs substantially. We conclude that the effectiveness of root hairs on root water uptake is determined by the hair-induced increase in root-soil contact and root hair shrinkage. Although the latter clearly reduces the effect of hairs on water uptake, our model still indicated facilitation of water uptake by root hairs at soil matric potentials from -1 to -0.1 MPa. Our findings provide new avenues towards a mechanistic understanding of the role of root hairs on water uptake.

Fluid transport in porous media based on differences in filter media morphology and biofilm growth in bioreactors.

To elucidate the effect of pore structure on bioclogging and seepage flow in bioreactors, we used X-ray computed tomography (X-CT) to investigate the changes in seepage flow of porous media in zeolite, gravel and ceramsite bioreactors with biofilm growth by injecting a non-ionic contrast medium iohexol. Based on the X-CT images using a ball-and-stick model, the highest average pore radius (R‾) and the average pore throat radius (r‾) in the ceramsite column were found under the initial conditions, which facilitated its permeability. The pore and throat of the gravel column were small and homogeneous relatively. Biofilm growth decreased the pore and pore throat in the columns. The total throat area of zeolite, gravel and ceramsite columns declined by 74%, 73% and 79% respectively. The zeolite column had the highest average pore throat, which contributed to its maximum conductivity subsequently after biofilm growth. Further, the fractal dimensions of the pore structure increased with biofilm growth, especially in the zeolite and ceramsite columns. The heterogeneity of the porous media was reinforced by the biofilm growth in the zeolite and ceramsite columns due to their higher heterogeneity initially. We also observed that an increase in heterogeneity of porous medium amplified the preferential flow and flow heterogeneities, especially in the zeolite and ceramsite columns.

Direct characterization of solute transport in unsaturated porous media using fast X-ray synchrotron microtomography.

Solute transport in unsaturated porous materials is a complex process, which exhibits some distinct features differentiating it from transport under saturated conditions. These features emerge mostly due to the different transport time scales at different regions of the flow network, which can be classified into flowing and stagnant regions, predominantly controlled by advection and diffusion, respectively. Under unsaturated conditions, the solute breakthrough curves show early arrivals and very long tails, and this type of transport is usually referred to as non-Fickian. This study directly characterizes transport through an unsaturated porous medium in three spatial dimensions at the resolution of 3.25 µm and the time resolution of 6 s. Using advanced high-speed, high-spatial resolution, synchrotron-based X-ray computed microtomography (sCT) we obtained detailed information on solute transport through a glass bead packing at different saturations. A large experimental dataset (>50 TB) was produced, while imaging the evolution of the solute concentration with time at any given point within the field of view. We show that the fluids' topology has a critical signature on the non-Fickian transport, which yet needs to be included in the Darcy-scale solute transport models. The three-dimensional (3D) results show that the fully mixing assumption at the pore scale is not valid, and even after injection of several pore volumes the concentration field at the pore scale is not uniform. Additionally, results demonstrate that dispersivity is changing with saturation, being twofold larger at the saturation of 0.52 compared to that at the fully saturated domain.

Pore-Scale Heterogeneities Improve the Degradation of a Self-Inhibiting Substrate: Insights from Reactive Transport Modeling.

In situ bioremediation is a common remediation strategy for many groundwater contaminants. It was traditionally believed that (in the absence of mixing-limitations) a better in situ bioremediation is obtained in a more homogeneous medium where the even distribution of both substrate and bacteria facilitates the access of a larger portion of the bacterial community to a higher amount of substrate. Such conclusions were driven with the typical assumption of disregarding substrate inhibitory effects on the metabolic activity of enzymes at high concentration levels. To investigate the influence of pore matrix heterogeneities on substrate inhibition, we use a numerical approach to solve reactive transport processes in the presence of pore-scale heterogeneities. To this end, a rigorous reactive pore network model is developed and used to model the reactive transport of a self-inhibiting substrate under both transient and steady-state conditions through media with various, spatially correlated, pore-size distributions. For the first time, we explore on the basis of a pore-scale model approach the link between pore-size heterogeneities and substrate inhibition. Our results show that for a self-inhibiting substrate, (1) pore-scale heterogeneities can consistently promote degradation rates at toxic levels, (2) the effect reverses when the concentrations fall to levels essential for microbial growth, and (3) an engineered combination of homogeneous and heterogeneous media can increase the overall efficiency of bioremediation.

Pore-Scale Investigation of Dynamic Immiscible Displacement in Layered Media using Synchrotron X-ray Microtomography.

Understanding the dynamics of immiscible fluid in a porous media is critical in many chemical and environmental engineering processes. However, the geological heterogeneity effect on multiphase flow behavior remains unclear. Here, the dynamics of immiscible fluid displacement and entrapment were experimentally demonstrated at pore-level using time-lapse synchrotron X-ray microtomography. A drainage-imbibition experiment was designed using an unconsolidated layered sand pack that comprised coarse sand and fine sand zones. There were significant differences between the two zones, with regard to the temporal variations in fluid saturation and morphological evolution of nonwetting fluid (oil) during imbibition. Highly connected oil clusters in the coarse zone broke up into many small fragments, whereas the cluster in the fine zone remained connected while spanning multiple pores. To further understand the impacts of pore size and connectivity on multiphase fluid dynamics, a new approach that tracks the temporal variation of immiscible fluid in individual pores was conducted. The surface area at the oil-water interface increased during imbibition, which is expected to facilitate mass transfer and surface interactions. Understanding immiscible fluid displacement in layered porous media at the pore-level could lead to more effective environmental remediation.

Study on Oil Recovery Mechanism of Polymer-Surfactant Flooding Using X-ray Microtomography and Integral Geometry.

Understanding pore-scale morphology and distribution of remaining oil in pore space are of great importance to carry out in-depth tapping of oil potential. Taking two water-wet cores from a typical clastic reservoir in China as an example, X-ray CT imaging is conducted at different experimental stages of water flooding and polymer-surfactant (P-S) flooding by using a high-resolution X-ray microtomography. Based on X-ray micro-CT image processing, 3D visualization of rock microstructure and fluid distribution at the pore scale is achieved. The integral geometry newly developed is further introduced to characterize pore-scale morphology and distribution of remaining oil in pore space. The underlying mechanism of oil recovery by P-S flooding is further explored. The results show that the average diameter of oil droplets gradually decreases, and the topological connectivity becomes worse after water flooding and P-S flooding. Due to the synergistic effect of "1 + 1 > 2" between the strong sweep efficiency of surfactant and the enlarged swept volume of the polymer, oil droplets with a diameter larger than 124.58 µm can be gradually stripped out by the polymer-surfactant system, causing a more scattered distribution of oil droplets in pore spaces of the cores. The network-like oil clusters are still dominant when water flooding is continued to 98% of water cut, but the dominant pore-scale oil morphology has evolved from network-like to porous-type and isolated-type after P-S flooding, which can provide strong support for further oil recovery in the later stage of chemical flooding.

Phase Saturation Control on Mixing-Driven Reactions in 3D Porous Media.

Transported chemical reactions in unsaturated porous media are relevant to environmental and industrial applications. Continuum scale models are based on equivalent parameters derived from analogy with saturated conditions and cannot appropriately account for incomplete mixing. It is also unclear how the third dimension controls mixing and reactions. We obtain three-dimensional (3D) images by magnetic resonance imaging using an immiscible nonwetting liquid as a second phase and a fast irreversible bimolecular reaction. We study the impact of phase saturation on the dynamics of mixing and the reaction front. We quantify the temporally resolved effective reaction rate and describe it using the lamellar theory of mixing, which explains faster than Fickian (t0.5) rate of product formation by accounting for the deformation of the mixing interface between the two reacting fluids. For a given Péclet, although stretching and folding of the reactive front enhance as saturation decreases, enhancing the product formation, the product formation is larger as saturation increases. After breakthrough, the extinction of the reaction takes longer as saturation decreases because of the larger nonmixed volume behind the front. These results are the basis for a general model to better predict reactive transport in unsaturated porous media not achievable by the current continuum paradigm.

The effects of long-term fertilizations on soil hydraulic properties vary with scales.

Soil structure is an indicator of soil quality and its alterations following cropping system conversion or fertilization change evolve slowly. How such alterations vary with scale remains elusive. We investigated this based on the Rothamsted long-term wheat experiment (since 1843) in the UK. Triplicate cores 7 cm high and 10 cm in diameter were taken from plots that have been under different fertilizations or returned to natural woodland for more than one century for imaging using X-ray computed tomography with the voxel size being 40 µm. We then broke each core and sampled three aggregates from it to scan with the voxel size being 1.5 µm. For each core and aggregate sample, we calculated its pore size distribution, permeability and tortuosity. The results showed that the fertilization change >170 years ago reshaped the soil structure differently between the core scale and the aggregate scale. Macro-porosity of the pores (>40 µm) in the cores unfertilized or fertilized with inorganic fertilizers was low and the pores were poorly connected in the top 10 cm of soil, compared to those given farmyard manure or in the woodland. In all treatments, the pores in the core images were hydraulically anisotropic with their permeability in the horizontal direction being higher than that in the vertical direction, whereas the aggregates were comparatively isotropic. The fertilization affected image porosity and permeability at core scale more significantly than at aggregate scale, and the aggregates fertilized with farmyard manure and in the woodland were more permeable than the aggregates in other treatments. It was also found that compared to no-fertilization or fertilization with complete fertilizers, fertilizing without phosphorus over the past 20 years increased the porosity and permeability of the aggregates but not of the cores. Fertilization with inorganic fertilizers increased the tortuosity of the macropores in the cores but not of the intra-aggregate pores, compared to no-fertilization. Porosity-permeability relationship for aggregates unfertilized or fertilized with inorganic fertilisers followed a power law with R 2 > 0.8. In contrast, the permeability of the aggregates in farmyard manure and in the woodland trended differently as their porosity increased. The results also revealed that the transport ability of the aggregates and cores responded differently to carbon in that with soil carbon increasing, the permeability of the aggregates increased asymptotically while the permeability of the cores, especially its horizontal component, increased exponentially.

Innovations in laboratory-based dynamic micro-CT to accelerate in situ research.

In the past few years, dynamic computed tomography (CT) approaches or uninterrupted acquisitions of deforming materials have rapidly emerged as an essential technique to understand material evolution, facilitating in situ investigations ranging from mechanical deformation to fluid flow in porous materials and beyond. Developments at synchrotron facilities have led this effort, pointing to the future of the technique. In the laboratory, recent developments at TESCAN XRE have made it possible to image, reconstruct and inspect dynamic processes in the laboratory with a temporal resolution below 10 s, meaning that an entire acquisition from 0 to 360° is completed within 10 s. The aim of this study is to explore the challenges and innovations that have led to the ability to perform high speed, dynamic acquisitions. A unique horizontally rotating gantry based micro-CT system was developed to facilitate complex in situ experiments. In doing so, the sample stays fixed while source and detector are uninterruptedly rotating around a vertical axis. In this work, the dynamic CT method with this rotating gantry based system will be described by two application examples: (1) deformation and collapse of a delicate beer foam and (2) in situ baking process of pastry. For the pastry baking process, an oven was needed to reach baking temperature. In a conventional micro-CT system, where the sample rotates, it is not so obvious to rotate an oven with sensor and heating cables. On the other hand, the delicate foam of a collapsing beer head is able to rotate, but because of the tangential convection during fast rotation (<10 s), it could influence the bubble detachment and liquid drainage and thus also the foam degradation. To investigate both processes, a horizontally rotating gantry based micro-CT is required. For both examples it was possible to quantify the key parameters such as pore size and distribution to better understand the rise and fall of porous foams. These examples will highlight the recent progress in adapting micro-CT workflows to accommodate uninterrupted imaging of dynamic events and point to opportunities for future continued development. LAY DESCRIPTION: Micro-CT allows the nondestructive visualisation of internal structures and is being used routinely in the field of Material Science, Geoscience, Life Science and more. Because of its nondestructive aspect, micro-CT is optimal to take repetitive scans of the same sample over time. The combination of taking different scans over time is so called time-resolved CT. By doing so, crucial insights can be obtained on how materials form, deform and perform over time or under certain external conditions. TESCAN XRE have made it possible to image, reconstruct and inspect dynamic processes in the laboratory with a temporal resolution below 10 s. The dynamic CT method will be described through the lens of two application examples: (1) deformation and collapse of a delicate beer foam and (2) in situ baking process of pastry. These examples will highlight the recent progress in adapting micro-CT workflows to accommodate imaging of dynamic events and point to opportunities for future continued development.

Dynamic fluid connectivity during steady-state multiphase flow in a sandstone.

The current conceptual picture of steady-state multiphase Darcy flow in porous media is that the fluid phases organize into separate flow pathways with stable interfaces. Here we demonstrate a previously unobserved type of steady-state flow behavior, which we term "dynamic connectivity," using fast pore-scale X-ray imaging. We image the flow of N2 and brine through a permeable sandstone at subsurface reservoir conditions, and low capillary numbers, and at constant fluid saturation. At any instant, the network of pores filled with the nonwetting phase is not necessarily connected. Flow occurs along pathways that periodically reconnect, like cars controlled by traffic lights. This behavior is consistent with an energy balance, where some of the energy of the injected fluids is sporadically converted to create new interfaces.

Review of Microfluidic Devices and Imaging Techniques for Fluid Flow Study in Porous Geomaterials.

Understanding transport phenomena and governing mechanisms of different physical and chemical processes in porous media has been a critical research area for decades. Correlating fluid flow behaviour at the micro-scale with macro-scale parameters, such as relative permeability and capillary pressure, is key to understanding the processes governing subsurface systems, and this in turn allows us to improve the accuracy of modelling and simulations of transport phenomena at a large scale. Over the last two decades, there have been significant developments in our understanding of pore-scale processes and modelling of complex underground systems. Microfluidic devices (micromodels) and imaging techniques, as facilitators to link experimental observations to simulation, have greatly contributed to these achievements. Although several reviews exist covering separately advances in one of these two areas, we present here a detailed review integrating recent advances and applications in both micromodels and imaging techniques. This includes a comprehensive analysis of critical aspects of fabrication techniques of micromodels, and the most recent advances such as embedding fibre optic sensors in micromodels for research applications. To complete the analysis of visualization techniques, we have thoroughly reviewed the most applicable imaging techniques in the area of geoscience and geo-energy. Moreover, the integration of microfluidic devices and imaging techniques was highlighted as appropriate. In this review, we focus particularly on four prominent yet very wide application areas, namely "fluid flow in porous media", "flow in heterogeneous rocks and fractures", "reactive transport, solute and colloid transport", and finally "porous media characterization". In summary, this review provides an in-depth analysis of micromodels and imaging techniques that can help to guide future research in the in-situ visualization of fluid flow in porous media.

On the challenges of greyscale-based quantifications using X-ray computed microtomography.

For X-ray computed microtomography (µ-CT) images of porous rocks where the grains and pores are not fully resolved, the greyscale values of each voxel can be used for quantitative calculations. This study addresses the challenges that arise with greyscale-based quantifications by conducting experiments designed to investigate the sources of error/uncertainty. We conduct greyscale-based calculations of porosity, concentration and diffusivity from various µ-CT experiments using a Bentheimer sandstone sample. The dry sandstone is imaged overtime to test the variation of greyscale values over sequential scans due to instrumentation stability. The sandstone is then imaged in a dry and contrast-agent saturated state at low resolution to determine a porosity map, which is compared to a porosity map derived from segmented high-resolution data. Then the linearity of the relationship between the concentration of a contrast agent and its corresponding attenuation coefficient is tested by imaging various solutions of known concentration. Lastly, a diffusion experiment is imaged at low resolution under dynamic conditions to determine local diffusivity values for the sandstone, which is compared to values derived from direct pore-scale simulations using high-resolution data. Overall, we identify the main errors associated with greyscale-based quantification and provide practical suggestions to alleviate these issues. LAY DESCRIPTION: X-ray computed microtomography (CT) imaging has become an important way to study the pore space of a porous medium. Using segmented images, we can build 3D pore space models for porous media and characterize the morphology and/or run simulations on the models. So, image segmentation is a critical image processing step. However, for low resolution images where image segmentation is not possible, grayscales are directly used for quantifications such as porosity and concentration calculations. Although these types of calculations have been widely accepted and used, the uncertainties and errors associated with grayscale-based quantifications are not fully discussed. Here we specifically design experiments with X-ray CT imaging to address the challenges that arise in grayscale-based quantifications. For instance, in order to validate porosity calculation results from low resolution images (with the help of high attenuating tracer), high resolution images are also acquired, which serve as a benchmark. The errors associated with concentration calculation using grayscale values are also discussed. In addition, numerical simulations using grayscale values are performed on a diffusion experiment images with X-ray CT. The problems that arise in dynamic imaging and the subsequent numerical simulations are discussed. The experiments, calculations and discussions provide a more comprehensive understanding on grayscale-based quantifications and aid in designing better X-ray CT experiments.

Droplet fragmentation: 3D imaging of a previously unidentified pore-scale process during multiphase flow in porous media.

Using X-ray computed microtomography, we have visualized and quantified the in situ structure of a trapped nonwetting phase (oil) in a highly heterogeneous carbonate rock after injecting a wetting phase (brine) at low and high capillary numbers. We imaged the process of capillary desaturation in 3D and demonstrated its impacts on the trapped nonwetting phase cluster size distribution. We have identified a previously unidentified pore-scale event during capillary desaturation. This pore-scale event, described as droplet fragmentation of the nonwetting phase, occurs in larger pores. It increases volumetric production of the nonwetting phase after capillary trapping and enlarges the fluid-fluid interface, which can enhance mass transfer between the phases. Droplet fragmentation therefore has implications for a range of multiphase flow processes in natural and engineered porous media with complex heterogeneous pore spaces.

Transport properties and size exclusion effects in wide-pore superficially porous particles.

The effects of hydrodynamic radius on the transport of solute molecules in packed beds of wide-pore superficially porous particles (SPP) are studied using pore-scale simulation. The free molecular diffusion rate varies with radius through the Stokes-Einstein relation. Lattice Boltzmann and Langevin methods are used to model fluid motion and the transport of an ensemble of solute molecules in the fluid, providing statistics on solute concentration, flux, molecule age and residence time, as a function of depth in the SPP. Intraparticle effective diffusion and bed dispersion coefficients are calculated and correlated with the hydrodynamic radius and accessible porosity. The relative importance of convection and diffusion are found to depend on the molecule (tracer) size through the diffusion rate, and convection effects are more significant for larger, slower-diffusing molecules. When larger molecules are utilized, the intraparticle concentration is reduced in proportion to the local particle porosity, leading to a natural definition of the accessible porosity used in size exclusion chromatography (SEC). Although the pore shape is complex, the SEC constant K can be calculated directly from simulation. Simulation demonstrates that the effective diffusion coefficient is elevated near the particle hull, which is largely open to interstitial flow, and decreases with depth into the particle. All molecules studied here have transport access to the entire particle depth, although the accessible volume at a given depth depends on their size. The first passage time into the particle is well predicted by the diffusion rate, but residence time is influenced by convection, shortening the average visit duration. These results are of interest in "perfusion" chromatography where convection is thought to increase separation efficiency for large biomolecules.

Dynamic Pore-Scale Model of Drainage in Granular Porous Media: The Pore-Unit Assembly Method.

Dynamics of drainage is analyzed for packings of spheres, using numerical experiments. For this purpose, a dynamic pore-scale model was developed to simulate water flow during drainage. The pore space inside a packing of spheres was extracted using regular triangulation, resulting in an assembly of grain-based tetrahedra. Then, pore units were constructed by identifying and merging tetrahedra that belong to the same pore, resulting in an assembly of pore units. Each pore unit was approximated by a volume-equivalent regular shape (e.g., cube and octahedron), for which a local capillary pressure-saturation relationship was obtained. To simulate unsaturated flow, a pore-scale version of IMPES (implicit pressure solver and explicit saturation update) was employed in order to calculate pressure and saturation distributions as a function of time for the assembly of pore units. To test the dynamic model, it was used on a packing of spheres to reproduce the corresponding measured quasi-static capillary pressure-saturation curve for a sand packing. Calculations were done for a packing of spheres with the same grain size distribution and porosity as the sand. We obtained good agreement, which confirmed the ability of the dynamic code to accurately describe drainage under low flow rates. Simulations of dynamic drainage revealed that drainage occurred in the form of finger-like infiltration of air into the pore space, caused by heterogeneities in the pore structure. During the finger-like infiltration, the pressure difference between air and water was found to be significantly higher than the capillary pressure. Furthermore, we tested the effects of the averaging, boundary conditions, domain size, and viscosity on the dynamic flow behavior. Finally, the dynamic coefficient was determined and compared to experimental data.

Water flow across the interface of contrasting materials: Pressure discontinuity and its implications.

Water flow along or across the interfaces of contrasting materials is ubiquitous in hydrology and how to solve them in macroscopic models derived from volumetric average of the pore-scale processes remains elusive. While the change in the average velocity and pressure at water-sediment interface has been well established for channel flow over porous beds, whether a volumetric average alerts the pressure continuity when water flows across the interface of two porous materials is poorly understood despite its imperative implications in hydrological modelling. The primary purpose of this paper is to provide evidences via pore-scale simulations that volumetrically averaging the pore-scale processes indeed yields a discontinuous pressure when water flows across a material interface. We simulated two columns numerically reconstructed by filling them with stratified media: One is an idealised two-layer system and the other one is a 3D column filled by fine glass beads over coarse glass beads with their pore geometry acquired using x-ray computed tomography. The pore-scale simulation is to mimic the column experiment by driving fluid to flow through the void space under an externally imposed pressure gradient. Once fluid flow reaches steady state, its velocity and pressure in all voxels are sampled and they are then spatially averaged over each section perpendicular to the average flow direction. The results show that the average pressure drops abruptly at the material interface no matter which direction the fluid flows. Compared with the effective permeability estimated from the homogenization methods well established in the literature, the emerged discontinuous pressure at the interface reduces the combined ability of the two strata to conduct water. It is also found that under certain circumstances fluid flow is direction-dependant, moving faster when flowing in the fine-coarse direction than in the coarse-to-fine direction under the same pressure gradient. Although significant efforts are needed to incorporate these findings into practical models, we do elicit the emergence of discontinuous pressure at material interface due to volumetric average as well as its consequent implications in modelling of flow in heterogeneous and stratified media.

Modeling multidimensional and multispecies biofilms in porous media.

Modeling multidimensional and multispecies biofilm in porous media at the pore scale is challenging due to the need to simultaneously track the microbial community in the biofilms and the interfaces between the biofilms and the fluid. Therefore, researchers usually assume that the model has only one dimension in space or has only one microbial species. This work uses bioremediation of U(VI)-contaminated groundwater as the context to develop a two-dimensional and multispecies biofilm model. The model simulates the transverse mixing zone in which U(VI) is mixed with propionate, a nutrient externally supplied to stimulate the growth of microorganisms. The model considers multiple interactions among fluid flow, transport and reaction of chemical species, and growth of biofilm. The biofilm consists of two types of active biomass (syntrophs and dissimilatory metal reducing bacteria [DMBR]) and inert biomass. The two types of active biomass collaboratively remove U(VI). The model outputs biomass distribution, chemical species concentrations, and fluid flow at the pore scale to fundamentally study the multiple interactions. The model also outputs the contaminant removal rate that can be potentially used for up-scaling studies. The simulated results are generally consistent with experimental observations from other studies in trend. The trend can be explained by the multiple interactions based on thermodynamics and microbial kinetics. Biotechnol. Bioeng. 2017;114: 1679-1687. © 2017 Wiley Periodicals, Inc.

Monitoring oil displacement processes with k-t accelerated spin echo SPI.

Magnetic resonance imaging (MRI) is a robust tool to monitor oil displacement processes in porous media. Conventional MRI measurement times can be lengthy, which hinders monitoring time-dependent displacements. Knowledge of the oil and water microscopic distribution is important because their pore scale behavior reflects the oil trapping mechanisms. The oil and water pore scale distribution is reflected in the magnetic resonance T2 signal lifetime distribution. In this work, a pure phase-encoding MRI technique, spin echo SPI (SE-SPI), was employed to monitor oil displacement during water flooding and polymer flooding. A k-t acceleration method, with low-rank matrix completion, was employed to improve the temporal resolution of the SE-SPI MRI measurements. Comparison to conventional SE-SPI T2 mapping measurements revealed that the k-t accelerated measurement was more sensitive and provided higher-quality results. It was demonstrated that the k-t acceleration decreased the average measurement time from 66.7 to 20.3 min in this work. A perfluorinated oil, containing no (1) H, and H2 O brine were employed to distinguish oil and water phases in model flooding experiments. High-quality 1D water saturation profiles were acquired from the k-t accelerated SE-SPI measurements. Spatially and temporally resolved T2 distributions were extracted from the profile data. The shift in the (1) H T2 distribution of water in the pore space to longer lifetimes during water flooding and polymer flooding is consistent with increased water content in the pore space. Copyright © 2015 John Wiley & Sons, Ltd.