Magnetic Resonance Imaging of Flow-Distributed Oscillations

Multiscale Porosity Microfluidics to Study Bacterial Transport in Heterogeneous Chemical Landscapes

Microfluidic models are proving to be powerful systems to study fundamental processes in porous media, due to their ability to replicate topologically complex environments while allowing detailed, quantitative observations at the pore scale. Yet, while porous media such as living tissues, geological substrates, or industrial systems typically display a porosity that spans multiple scales, most microfluidic models to date are limited to a single porosity or a small range of pore sizes. Here, a novel microfluidic system with multiscale porosity is presented. By embedding polyacrylamide (PAAm) hydrogel structures through in-situ photopolymerization in a landscape of microfabricated polydimethylsiloxane (PDMS) pillars with varying spacing, micromodels with porosity spanning several orders of magnitude, from nanometers to millimeters are created. Experiments conducted at different porosity patterns demonstrate the potential of this approach to characterize fundamental and ubiquitous biological and geochemical transport processes in porous media. Accounting for multiscale porosity allows studies of the resulting heterogeneous fluid flow and concentration fields of transported chemicals, as well as the biological behaviors associated with this heterogeneity, such as bacterial chemotaxis. This approach brings laboratory studies of transport in porous media a step closer to their natural counterparts in the environment, industry, and medicine.

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 (t^0.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.

MRI of chemical reactions and processes

As magnetic resonance imaging (MRI) can spatially resolve a wealth of molecular information available from nuclear magnetic resonance (NMR), it is able to non-invasively visualise the composition, properties and reactions of a broad range of spatially-heterogeneous molecular systems. Hence, MRI is increasingly finding applications in the study of chemical reactions and processes in a diverse range of environments and technologies. This article will explain the basic principles of MRI and how it can be used to visualise chemical composition and molecular properties, providing an overview of the variety of information available. Examples are drawn from the disciplines of chemistry, chemical engineering, environmental science, physics, electrochemistry and materials science. The review introduces a range of techniques used to produce image contrast, along with the chemical and molecular insight accessible through them. Methods for mapping the distribution of chemical species, using chemical shift imaging or spatially-resolved spectroscopy, are reviewed, as well as methods for visualising physical state, temperature, current density, flow velocities and molecular diffusion. Strategies for imaging materials with low signal intensity, such as those containing gases or low sensitivity nuclei, using compressed sensing, para-hydrogen or polarisation transfer, are discussed. Systems are presented which encapsulate the diversity of chemical and physical parameters observable by MRI, including one- and two-phase flow in porous media, chemical pattern formation, phase transformations and hydrodynamic (fingering) instabilities. Lastly, the emerging area of electrochemical MRI is discussed, with studies presented on the visualisation of electrochemical deposition and dissolution processes during corrosion and the operation of batteries, supercapacitors and fuel cells.

Eigenvector centrality for geometric and topological characterization of porous media

Solving flow and transport through complex geometries such as porous media is computationally difficult. Such calculations usually involve the solution of a system of discretized differential equations, which could lead to extreme computational cost depending on the size of the domain and the accuracy of the model. Geometric simplifications like pore networks, where the pores are represented by nodes and the pore throats by edges connecting pores, have been proposed. These models, despite their ability to preserve the connectivity of the medium, have difficulties capturing preferential paths (high velocity) and stagnation zones (low velocity), as they do not consider the specific relations between nodes. Nonetheless, network theory approaches, where a complex network is a graph, can help to simplify and better understand fluid dynamics and transport in porous media. Here we present an alternative method to address these issues based on eigenvector centrality, which has been corrected to overcome the centralization problem and modified to introduce a bias in the centrality distribution along a particular direction to address the flow and transport anisotropy in porous media. We compare the model predictions with millifluidic transport experiments, which shows that, albeit simple, this technique is computationally efficient and has potential for predicting preferential paths and stagnation zones for flow and transport in porous media. We propose to use the eigenvector centrality probability distribution to compute the entropy as an indicator of the "mixing capacity" of the system.

Helical Turing patterns in the Lengyel-Epstein model in thin cylindrical layers

The formation of Turing patterns was investigated in thin cylindrical layers using the Lengyel-Epstein model of the chlorine dioxide-iodine-malonic acid reaction. The influence of the width of the layer W and the diameter D of the inner cylinder on the pattern with intrinsic wavelength l were determined in simulations with initial random noise perturbations to the uniform state for W < l/2 and D ∼ l or lower. We show that the geometric constraints of the reaction domain may result in the formation of helical Turing patterns with parameters that give stripes (b = 0.2) or spots (b = 0.37) in two dimensions. For b = 0.2, the helices were composed of lamellae and defects were likely as the diameter of the cylinder increased. With b = 0.37, the helices consisted of semi-cylinders and the orientation of stripes on the outer surface (and hence winding number) increased with increasing diameter until a new stripe appeared.

Chemical patterns in translating vortices: inter- and intra-cellular mixing effects

Stationary chemical patterns-flow distributed oscillations (FDOs)-are obtained when the Belousov-Zhabotinsky reaction is coupled with translating vortex flow in a Vortex Flow Reactor. For certain conditions, the FDOs are unstable with the observation of disappearing bands or complex patterns. The transitions between modes of pattern formation are reproduced in a modified Oregonator model consisting of two-zone cells connected in series. We show that increasing inter-cellular mixing of the outer zones results in a transition from FDO to absolute instabilities (AI) and increasing intra-cellular mixing between the core and outer zones can drive the reverse transition between modes (AI to FDO).

Distinguishing similar patterns with different underlying instabilities: effect of advection on systems with Hopf, Turing-Hopf, and wave instabilities

Systems with the same local dynamics but different types of diffusive instabilities may show the same type of patterns. In this paper, we show that under the influence of advective flow the scenario of patterns that is formed at different velocities change; therefore, we propose the use of advective flow as a tool to uncover the underlying instabilities of a reaction-diffusion system.

Magnetic resonance imaging of chemistry

Magnetic resonance imaging (MRI) has long been recognized as one of the most important tools in medical diagnosis and research. However, MRI is also well placed to image chemical reactions and processes, determine the concentration of chemical species, and look at how chemistry couples with environmental factors, such as flow and heterogeneous media. This tutorial review will explain how magnetic resonance imaging works, reviewing its application in chemistry and its ability to directly visualise chemical processes. It will give information on what resolution and contrast are possible, and what chemical and physical parameters can be measured. It will provide examples of the use of MRI to study chemical systems, its application in chemical engineering and the identification of contrast agents for non-clinical applications. A number of studies are presented including investigation of chemical conversion and selectivity in fixed-bed reactors, temperature probes for catalyst pellets, ion mobility during tablet dissolution, solvent dynamics and ion transport in Nafion polymers and the formation of chemical waves and patterns.

Inward propagating chemical waves in Taylor vortices

Advection-reaction-diffusion (ARD) waves in the Belousov-Zhabotinsky reaction in steady Taylor-Couette vortices have been visualized using magnetic-resonance imaging and simulated using an adapted Oregonator model. We show how propagating wave behavior depends on the ratio of advective, chemical and diffusive time scales. In simulations, inward propagating spiral flamelets are observed at high Damköhler number (Da). At low Da, the reaction distributes itself over several vortices and then propagates inwards as contracting ring pulses--also observed experimentally.

Temperature Changes Visualization during Chemical Wave Propagation

Magnetic resonance imaging was used for two-dimensional temperature visualization of chemical waves propagation in the autocatalytic exothermal reaction of thiosulfate oxidation by chlorite. The technique presented is based on the temperature dependence of the water chemical shift. Temperature maps were acquired by employing the TurboFLASH imaging method. The results obtained allow one to judge about directions of buoyancy flows. Two types of convection critical modes in a vertical tube during the wave propagation were detected.

Magnetic resonance imaging of chemical waves in porous media

Magnetic resonance imaging (MRI) provides a powerful tool for the investigation of chemical structures in optically opaque porous media, in which chemical concentration gradients can be visualized, and diffusion and flow properties are simultaneously determined. In this paper we give an overview of the MRI technique and review theory and experiments on the formation of chemical waves in a tubular packed bed reactor upon the addition of a nonlinear chemical reaction. MR images are presented of reaction-diffusion waves propagating in the three-dimensional (3D) network of channels in the reactor, and the 3D structure of stationary concentration patterns formed via the flow-distributed oscillation mechanism is demonstrated to reflect the local hydrodynamics in the packed bed. Possible future directions regarding the influence of heterogeneities on transport and reaction are discussed.