skip to main content


Title: Comprehensive comparison of pore-scale models for multiphase flow in porous media
Multiphase flows in porous media are important in many natural and industrial processes. Pore-scale models for multiphase flows have seen rapid development in recent years and are becoming increasingly useful as predictive tools in both academic and industrial applications. However, quantitative comparisons between different pore-scale models, and between these models and experimental data, are lacking. Here, we perform an objective comparison of a variety of state-of-the-art pore-scale models, including lattice Boltzmann, stochastic rotation dynamics, volume-of-fluid, level-set, phase-field, and pore-network models. As the basis for this comparison, we use a dataset from recent microfluidic experiments with precisely controlled pore geometry and wettability conditions, which offers an unprecedented benchmarking opportunity. We compare the results of the 14 participating teams both qualitatively and quantitatively using several standard metrics, such as fractal dimension, finger width, and displacement efficiency. We find that no single method excels across all conditions and that thin films and corner flow present substantial modeling and computational challenges.  more » « less
Award ID(s):
1803989
NSF-PAR ID:
10106702
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; « less
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Volume:
116
Issue:
28
ISSN:
0027-8424
Page Range / eLocation ID:
13799 to 13806
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    The Amundsen Sea sector of Antarctica has long been considered the most vulnerable part of the West Antarctic Ice Sheet (WAIS) because of the great water depth at the grounding line and the absence of substantial ice shelves. Glaciers in this configuration are thought to be susceptible to rapid or runaway retreat. Ice flowing into the Amundsen Sea Embayment is undergoing the most rapid changes of any sector of the Antarctic Ice Sheet outside the Antarctic Peninsula, including changes caused by substantial grounding-line retreat over recent decades, as observed from satellite data. Recent models suggest that a threshold leading to the collapse of WAIS in this sector may have been already crossed and that much of the ice sheet could be lost even under relatively moderate greenhouse gas emission scenarios. Drill cores from the Amundsen Sea provide tests of several key questions about controls on ice sheet stability. The cores offer a direct record of glacial history offshore from a drainage basin that receives ice exclusively from the WAIS, which allows clear comparisons between the WAIS history and low-latitude climate records. Today, warm Circumpolar Deep Water (CDW) is impinging onto the Amundsen Sea shelf and causing melting of the underside of the WAIS in most places. Reconstructions of past CDW intrusions can assess the ties between warm water upwelling and large-scale changes in past grounding-line positions. Carrying out these reconstructions offshore from the drainage basin that currently has the most substantial negative mass balance of ice anywhere in Antarctica is thus of prime interest to future predictions. The scientific objectives for this expedition are built on hypotheses about WAIS dynamics and related paleoenvironmental and paleoclimatic conditions. The main objectives are 1. To test the hypothesis that WAIS collapses occurred during the Neogene and Quaternary and, if so, when and under which environmental conditions; 2. To obtain ice-proximal records of ice sheet dynamics in the Amundsen Sea that correlate with global records of ice-volume changes and proxy records for atmospheric and ocean temperatures; 3. To study the stability of a marine-based WAIS margin and how warm deep-water incursions control its position on the shelf; 4. To find evidence for earliest major grounded WAIS advances onto the middle and outer shelf; 5. To test the hypothesis that the first major WAIS growth was related to the uplift of the Marie Byrd Land dome. International Ocean Discovery Program (IODP) Expedition 379 completed two very successful drill sites on the continental rise of the Amundsen Sea. Site U1532 is located on a large sediment drift, now called Resolution Drift, and penetrated to 794 m with 90% recovery. We collected almost-continuous cores from the Pleistocene through the Pliocene and into the late Miocene. At Site U1533, we drilled 383 m (70% recovery) into the more condensed sequence at the lower flank of the same sediment drift. The cores of both sites contain unique records that will enable study of the cyclicity of ice sheet advance and retreat processes as well as bottom-water circulation and water mass changes. In particular, Site U1532 revealed a sequence of Pliocene sediments with an excellent paleomagnetic record for high-resolution climate change studies of the previously sparsely sampled Pacific sector of the West Antarctic margin. Despite the drilling success at these sites, the overall expedition experienced three unexpected difficulties that affected many of the scientific objectives: 1. The extensive sea ice on the continental shelf prevented us from drilling any of the proposed shelf sites. 2. The drill sites on the continental rise were in the path of numerous icebergs of various sizes that frequently forced us to pause drilling or leave the hole entirely as they approached the ship. The overall downtime caused by approaching icebergs was 50% of our time spent on site. 3. An unfortunate injury to a member of the ship's crew cut the expedition short by one week. Recovery of core on the continental rise at Sites U1532 and U1533 cannot be used to precisely indicate the position of ice or retreat of the ice sheet on the shelf. However, these sediments contained in the cores offer a range of clues about past WAIS extent and retreat. At Sites U1532 and U1533, coarse-grained sediments interpreted to be ice-rafted debris (IRD) were identified throughout all recovered time periods. A dominant feature of the cores is recorded by lithofacies cyclicity, which is interpreted to represent relatively warmer periods variably characterized by higher microfossil abundance, greater bioturbation, and higher counts of IRD alternating with colder periods characterized by dominantly gray laminated terrigenous muds. Initial comparison of these cycles to published records from the region suggests that the units interpreted as records of warmer time intervals in the core tie to interglacial periods and the units interpreted as deposits of colder periods tie to glacial periods. The cores from the two drill sites recovered sediments of purely terrigenous origin intercalated or mixed with pelagic or hemipelagic deposits. In particular, Site U1533, which is located near a deep-sea channel originating from the continental slope, contains graded sands and gravel transported downslope from the shelf to the abyssal plain. The channel is likely the path of such sediments transported downslope by turbidity currents or other sediment-gravity flows. The association of lithologic facies at both sites predominantly reflects the interplay of downslope and contouritic sediment supply with occasional input of more pelagic sediment. Despite the lack of cores from the shelf, our records from the continental rise reveal the timing of glacial advances across the shelf and thus the existence of a continent-wide ice sheet in West Antarctica at least during longer time periods since the late Miocene. Cores from both sites contain abundant coarse-grained sediments and clasts of plutonic origin transported either by downslope processes or by ice rafting. If detailed provenance studies confirm our preliminary assessment that the origin of these samples is from the plutonic bedrock of Marie Byrd Land, their thermochronological record will potentially reveal timing and rates of denudation and erosion linked to crustal uplift. The chronostratigraphy of both sites enables the generation of a seismic sequence stratigraphy not only for the Amundsen Sea rise but also for the western Amundsen Sea along the Marie Byrd Land margin through a connecting network of seismic lines. 
    more » « less
  2. Summary

    Industrial symbiosis (IS), as a subfield of industrial ecology, is concerned with cooperation among industrial firms in managing resources, particularly by‐products, such that the waste of one firm becomes the input of another. This “closed‐loop” pattern also lies at the heart of the concept of the circular economy (CE). Both concepts are typically considered at scales ranging from industrial parks to global supply chains, but rarely at the scale of a single facility housing multiple firms. Moreover, both concepts have requirements and implications that reach beyond the tracking and coordination of material flows. These additional requirements include the cultivation of new business models and political support, as well as engagement of local communities and cultures. Thus, both IS and the CE face the challenge of bridging the gap between the technical and sociocultural aspects of industrial development and adapting to the facility scale. We address this challenge by proposing a framework for the development of facility‐scale industrial symbiosis (Facility‐IS). Facility‐IS accounts for the presence of symbiotic interfirm relationships—including, but not limited to, material and knowledge exchange—within a single physical facility. The conditions for planning, facilitating, and growing Facility‐IS are outlined, as well as three approaches to implement such projects. Overall, this research forms a Facility‐IS framework rooted in existing best practices for understanding the CE at the local scale, tracking business sustainability factors and assessing viable strategies that facilitate IS.

     
    more » « less
  3. Abstract

    Geophysical fluid‐granular flows, such as pyroclastic currents and debris flows, owe much of their runout and hazard behavior to the occurrence and time‐variant decay of a flow‐internal fluid pore pressure. However, modeling the effects of fluid pore pressure to forecast hazards is challenging because a unified method in Earth Sciences to quantitatively determine the permeability of these natural mixtures is currently missing. Here we combine experiments on fluidization and defluidization of pyroclastic materials, eolian sediments, and glass beads mixtures with numerical multiphase simulations to compare previous attempts to compute the permeability of complex natural particle‐fluid mixtures. In analogy to particle‐engineering studies on simple gas‐particle mixtures, we demonstrate that the effective length‐scale in the characterization of the fluid‐particle interaction of complex natural mixtures is the product of the Sauter mean diameter and the particle sphericity. Its use in the Kozeny‐Carman equation allows accurate prediction of mixture permeability, and we suggest the routine calculation of the Sauter mean from grain size distributions of the deposits of geophysical mass flows in Earth Sciences. We also show, through defluidization experiments, that the duration of gas retention in natural mixtures is well described when using the Sauter mean as the effective particle size. Further, we show through multiphase simulations that initial bed expansion extends the pore pressure diffusion timescale up to nine times. These results can be applied to small‐to‐large volume dense pyroclastic currents where the ranges of Sauter mean diameter predict gas retention for long duration and to debris flows and snow avalanches.

     
    more » « less
  4. Dense granular heap flows are common in nature, such as during avalanches and landslides, as well as in industrial flows. In granular heap flows, rapid flow is localized near the free surface with the thickness of the rapidly flowing layer dependent on the overall flow rate. In the region deep beneath the surface, exponentially decaying creeping flow dominates with characteristic decay length depending only on the geometry and not the overall flow rate. Existing continuum models for dense granular flow based upon local constitutive equations are not able to simultaneously predict both of these experimentally observed features – failing to even predict the existence of creeping flow beneath the surface. In this work, we apply a scale-dependent continuum approach – the non-local granular fluidity model – to steady, dense granular flows on a heap between two smooth, frictional side walls. We show that the model captures the salient features of both the flow-rate-dependent, rapidly flowing surface layer and the flow-rate-independent, slowly creeping bulk under steady flow conditions. 
    more » « less
  5. Complex and active fluids find broad applications in flows through porous materials. Nontrivial rheology can couple to porous microstructure leading to surprising flow patterns and associated transport properties in geophysical, biological, and industrial systems. Viscoelastic instabilities are highly sensitive to pore geometry and can give rise to chaotic velocity fluctuations. A number of recent studies have begun to untangle how the pore-scale geometry influences the sample-scale flow topology and the resulting dispersive transport properties of these complex systems. Beyond classical rheological properties, active colloids and swimming cells exhibit a range of unique properties, including reduced effective viscosity, collective motion, and random walks, that present novel challenges to understanding their mechanics and transport in porous media flows. This review article aims to provide a brief overview of essential, fundamental concepts followed by an in-depth summary of recent developments in this rapidly evolving field. The chosen topics are motivated by applications, and new opportunities for discovery are highlighted.

     
    more » « less