Lab on a chip for a low-carbon future.

Transitioning our society to a sustainable future, with low or net-zero carbon emissions to the atmosphere, will require a wide-spread transformation of energy and environmental technologies. In this perspective article, we describe how lab-on-a-chip (LoC) systems can help address this challenge by providing insight into the fundamental physical and geochemical processes underlying new technologies critical to this transition, and developing the new processes and materials required. We focus on six areas: (I) subsurface carbon sequestration, (II) subsurface hydrogen storage, (III) geothermal energy extraction, (IV) bioenergy, (V) recovering critical materials, and (VI) water filtration and remediation. We hope to engage the LoC community in the many opportunities within the transition ahead, and highlight the potential of LoC approaches to the broader community of researchers, industry experts, and policy makers working toward a low-carbon future.

Probing multiscale dissolution dynamics in natural rocks through microfluidics and compositional analysis.

Mineral dissolution significantly impacts many geological systems. Carbon released by diagenesis, carbon sequestration, and acid injection are examples where geochemical reactions, fluid flow, and solute transport are strongly coupled. The complexity in these systems involves interplay between various mechanisms that operate at timescales ranging from microseconds to years. Current experimental techniques characterize dissolution processes using static images that are acquired with long measurement times and/or low spatial resolution. These limitations prevent direct observation of how dissolution reactions progress within an intact rock with spatially heterogeneous mineralogy and morphology. We utilize microfluidic cells embedded with thin rock samples to visualize dissolution with significant temporal resolution (100 ms) in a large observation window (3 × 3 mm). We injected acidic fluid into eight shale samples ranging from 8 to 86 wt % carbonate. The pre- and postreaction microstructures are characterized at the scale of pores (0.1 to 1 µm) and fractures (1 to 1,000 µm). We observe that nonreactive particle exposure, fracture morphology, and loss of rock strength are strongly dependent on both the relative volume of reactive grains and their distribution. Time-resolved images of the rock unveil the spatiotemporal dynamics of dissolution, including two-phase flow effects in real time and illustrate the changes in the fracture interface across the range of compositions. Moreover, the dynamical data provide an approach for characterizing reactivity parameters of natural heterogeneous samples when porous media effects are not negligible. The platform and workflow provide real-time characterization of geochemical reactions and inform various subsurface engineering processes.

Module-Fluidics: Building Blocks for Spatio-Temporal Microenvironment Control.

Generating the desired solute concentration signal in micro-environments is vital to many applications ranging from micromixing to analyzing cellular response to a dynamic microenvironment. We propose a new modular design to generate targeted temporally varying concentration signals in microfluidic systems while minimizing perturbations to the flow field. The modularized design, here referred to as module-fluidics, similar in principle to interlocking toy bricks, is constructed from a combination of two building blocks and allows one to achieve versatility and flexibility in dynamically controlling input concentration. The building blocks are an oscillator and an integrator, and their combination enables the creation of controlled and complex concentration signals, with different user-defined time-scales. We show two basic connection patterns, in-series and in-parallel, to test the generation, integration, sampling and superposition of temporally-varying signals. All such signals can be fully characterized by analytic functions, in analogy with electric circuits, and allow one to perform design and optimization before fabrication. Such modularization offers a versatile and promising platform that allows one to create highly customizable time-dependent concentration inputs which can be targeted to the specific application of interest.

Chemical and Reactive Transport Processes Associated with Hydraulic Fracturing of Unconventional Oil/Gas Shales.

Hydraulic fracturing of unconventional oil/gas shales has changed the energy landscape of the U.S. Recovery of hydrocarbons from tight, hydraulically fractured shales is a highly inefficient process, with estimated recoveries of <25% for natural gas and <5% for oil. This review focuses on the complex chemical interactions of additives in hydraulic fracturing fluid (HFF) with minerals and organic matter in oil/gas shales. These interactions are intended to increase hydrocarbon recovery by increasing porosities and permeabilities of tight shales. However, fluid-shale interactions result in the dissolution of shale minerals and the release and transport of chemical components. They also result in mineral precipitation in the shale matrix, which can reduce permeability, porosity, and hydrocarbon recovery. Competition between mineral dissolution and mineral precipitation processes influences the amounts of oil and gas recovered. We review the temporal/spatial origins and distribution of unconventional oil/gas shales from mudstones and shales, followed by discussion of their global and U.S. distributions and compositional differences from different U.S. sedimentary basins. We discuss the major types of chemical additives in HFF with their intended purposes, including drilling muds. Fracture distribution, porosity, permeability, and the identity and molecular-level speciation of minerals and organic matter in oil/gas shales throughout the hydraulic fracturing process are discussed. Also discussed are analysis methods used in characterizing oil/gas shales before and after hydraulic fracturing, including permeametry and porosimetry measurements, X-ray diffraction/Rietveld refinement, X-ray computed tomography, scanning/transmission electron microscopy, and laboratory- and synchrotron-based imaging/spectroscopic methods. Reactive transport and spatial scaling are discussed in some detail in order to relate fundamental molecular-scale processes to fluid transport. Our review concludes with a discussion of potential environmental impacts of hydraulic fracturing and important knowledge gaps that must be bridged to achieve improved mechanistic understanding of fluid transport in oil/gas shales.

Dynamic Modeling of Fouling in Reverse Osmosis Membranes.

During reverse osmosis (RO) membrane filtration, performance is dramatically affected by fouling, which concurrently decreases the permeate flux while increasing the energy required to operate the system. Comprehensive design and optimization of RO systems are best served by an understanding of the coupling between membrane shape, local flow field, and fouling; however, current studies focus exclusively on simplified steady-state models that ignore the dynamic coupling between fluid flow, solute transport, and foulant accumulation. We developed a customized solver (SUMs: Stanford University Membrane Solver) under the open source finite volume simulator OpenFOAM to solve transient Navier-Stokes, advection-diffusion, and adsorption-desorption equations for foulant accumulation. We implemented two permeate flux reduction models at the membrane boundary: the resistance-in-series (RIS) model and the effective-pressure-drop (EPD) model. The two models were validated against filtration experiments by comparing the equilibrium flux, pressure drop, and fouling pattern on the membrane. Both models not only predict macroscopic quantities (e.g., permeate flux and pressure drop) but also the fouling pattern developed on the membrane, with a good match with experimental results. Furthermore, the models capture the temporal evolution of foulant accumulation and its coupling with flux reduction.

Universal scaling-law for flow resistance over canopies with complex morphology.

Flow resistance caused by vegetation is a key parameter to properly assess flood management and river restoration. However, quantifying the friction factor or any of its alternative metrics, e.g. the drag coefficient, in canopies with complex geometry has proven elusive. We explore the effect of canopy morphology on vegetated channels flow structure and resistance by treating the canopy as a porous medium characterized by an effective permeability, a property that describes the ease with which water can flow through the canopy layer. We employ a two-domain model for flow over and within the canopy, which couples the log-law in the free layer to the Darcy-Brinkman equation in the vegetated layer. We validate the model analytical solutions for the average velocity profile within and above the canopy, the volumetric discharge and the friction factor against data collected across a wide range of canopy morphologies encountered in riverine systems. Results indicate agreement between model predictions and data for both simple and complex plant morphologies. For low submergence canopies, we find a universal scaling law that relates friction factor with canopy permeability and a rescaled bulk Reynolds number. This provides a valuable tool to assess habitats sustainability associated with hydro-dynamical conditions.

An analysis platform for multiscale hydrogeologic modeling with emphasis on hybrid multiscale methods.

One of the most significant challenges faced by hydrogeologic modelers is the disparity between the spatial and temporal scales at which fundamental flow, transport, and reaction processes can best be understood and quantified (e.g., microscopic to pore scales and seconds to days) and at which practical model predictions are needed (e.g., plume to aquifer scales and years to centuries). While the multiscale nature of hydrogeologic problems is widely recognized, technological limitations in computation and characterization restrict most practical modeling efforts to fairly coarse representations of heterogeneous properties and processes. For some modern problems, the necessary level of simplification is such that model parameters may lose physical meaning and model predictive ability is questionable for any conditions other than those to which the model was calibrated. Recently, there has been broad interest across a wide range of scientific and engineering disciplines in simulation approaches that more rigorously account for the multiscale nature of systems of interest. In this article, we review a number of such approaches and propose a classification scheme for defining different types of multiscale simulation methods and those classes of problems to which they are most applicable. Our classification scheme is presented in terms of a flowchart (Multiscale Analysis Platform), and defines several different motifs of multiscale simulation. Within each motif, the member methods are reviewed and example applications are discussed. We focus attention on hybrid multiscale methods, in which two or more models with different physics described at fundamentally different scales are directly coupled within a single simulation. Very recently these methods have begun to be applied to groundwater flow and transport simulations, and we discuss these applications in the context of our classification scheme. As computational and characterization capabilities continue to improve, we envision that hybrid multiscale modeling will become more common and also a viable alternative to conventional single-scale models in the near future.

Effective medium theory for drag-reducing micro-patterned surfaces in turbulent flows.

Many studies in the last decade have revealed that patterns at the microscale can reduce skin drag. Yet, the mechanisms and parameters that control drag reduction, e.g. Reynolds number and pattern geometry, are still unclear. We propose an effective medium representation of the micro-features, that treats the latter as a porous medium, and provides a framework to model turbulent flow over patterned surfaces. Our key result is a closed-form expression for the skin friction coefficient in terms of frictional Reynolds (or Kármán) number in turbulent regime, the viscosity ratio between the fluid in and above the features, and their geometrical properties. We apply the proposed model to turbulent flows over superhydrophobic ridged surfaces. The model predictions agree with laboratory experiments for Reynolds numbers ranging from 3000 to 10000.

Flow-induced shear instabilities of cohesive granulates.

In this work we use a multiscale framework to calculate the fluidization threshold of three-dimensional cohesive granulates under shear forces exerted by a creeping flow. A continuum model of flow through porous media provides an analytical expression for the average drag force on a single grain. The balance equation for the forces and a force propagation model are then used to investigate the effects of porosity and packing structure on the stability of the pile. We obtain a closed-form expression for the instability threshold of a regular packing of monodisperse frictionless cohesive spherical grains in a planar fracture. Our result quantifies the compound effect of structural (packing orientation and porosity) and dynamical properties of the system on its stability.

Elastic response of carbon nanotube forests to aerodynamic stresses.

The ability to determine static and (hydro)dynamic properties of carbon nanotubes (CNTs) is crucial for many applications. While their static properties (e.g., solubility and wettability) are fairly well understood, their mechanical responses (e.g., deflection under shear) to ambient fluid flow are to a large extent unknown. We analyze the elastic response of single-walled CNT forests, attached to the bottom wall of a channel, to the aerodynamic loading exerted by both laminar and turbulent flows. Our analysis yields analytical expressions for velocity distributions, the drag coefficient, and bending profiles of individual CNTs. This enables us to determine flexural rigidity of CNTs in wind-tunnel experiments. The model predictions agree with laboratory experiments for a large range of channel velocities.