An official website of the United States government
Abstract Large discrepancies between well-mixed reaction rates and effective reactions rates estimated under fluid flow conditions have been a major issue for predicting reactive transport in porous media systems. In this study, we introduce a framework that accurately predicts effective reaction rates directly from pore structural features by combining 3D pore-scale numerical simulations with machine learning (ML). We first perform pore-scale reactive transport simulations with fluid–solid reactions in hundreds of porous media and calculate effective reaction rates from pore-scale concentration fields. We then train a Random Forests model with 11 pore structural features and effective reaction rates to quantify the importance of structural features in determining effective reaction rates. Based on the importance information, we train artificial neural networks with varying number of features and demonstrate that effective reaction rates can be accurately predicted with only three pore structural features, which are specific surface, pore sphericity, and coordination number. Finally, global sensitivity analyses using the ML model elucidates how the three structural features affect effective reaction rates. The proposed framework enables accurate predictions of effective reaction rates directly from a few measurable pore structural features, and the framework is readily applicable to a wide range of applications involving porous media flows.
more »
« less
Simulating fluid flow in complex porous materials by integrating the governing equations with deep-layered machines
Abstract Fluid flow in heterogeneous porous media arises in many systems, from biological tissues to composite materials, soil, wood, and paper. With advances in instrumentations, high-resolution images of porous media can be obtained and used directly in the simulation of fluid flow. The computations are, however, highly intensive. Although machine learning (ML) algorithms have been used for predicting flow properties of porous media, they lack a rigorous, physics-based foundation and rely on correlations. We introduce an ML approach that incorporates mass conservation and the Navier–Stokes equations in its learning process. By training the algorithm to relatively limited data obtained from the solutions of the equations over a time interval, we show that the approach provides highly accurate predictions for the flow properties of porous media at all other times and spatial locations, while reducing the computation time. We also show that when the network is used for a different porous medium, it again provides very accurate predictions.
more »
« less
- Award ID(s):
- 2000968
- PAR ID:
- 10389851
- Date Published:
- Journal Name:
- npj Computational Materials
- Volume:
- 7
- Issue:
- 1
- ISSN:
- 2057-3960
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We present a new computational fluid dynamics approach for simulating two‐phase flow in hybrid systems containing solid‐free regions and deformable porous matrices. Our approach is based on the derivation of a unique set of volume‐averaged partial differential equations that asymptotically approach the Navier‐Stokes Volume‐of‐Fluid equations in solid‐free regions and multiphase Biot Theory in porous regions. The resulting equations extend our recently developed Darcy‐Brinkman‐Biot framework to multiphase flow. Through careful consideration of interfacial dynamics (relative permeability and capillary effects) and extensive benchmarking, we show that the resulting model accurately captures the strong two‐way coupling that is often exhibited between multiple fluids and deformable porous media. Thus, it can be used to represent flow‐induced material deformation (swelling, compression) and failure (cracking, fracturing). The model's open‐source numerical implementation,hybridBiotInterFoam, effectively marks the extension of computational fluid mechanics into modeling multiscale multiphase flow in deformable porous systems. The versatility of the solver is illustrated through applications related to material failure in poroelastic coastal barriers and surface deformation due to fluid injection in poro‐visco‐plastic systems.more » « less
-
Abstract Wettability is one of the critical parameters affecting multiphase flow in porous media. The wettability is determined by the affinity of fluids to the rock surface, which varies due to factors such as mineral heterogeneity, roughness, ageing, and pore-space geometry. It is well known that wettability varies spatially in natural rocks, and it is still generally considered a constant parameter in pore-scale simulation studies. The accuracy of pore-scale simulation of multiphase flow in porous media is undermined by such inadequate wettability models. The advent of in situ visualization techniques, e.g. X-ray imaging and microtomography, enables us to characterize the spatial distribution of wetting more accurately. There are several approaches for such characterization. Most include the construction of a meshed surface of the interface surfaces in a segmented X-ray image and are known to have significant errors arising from insufficient resolution and surface-smoothing algorithms. This work presents a novel approach for spatial determination of wetting properties using local lattice-Boltzmann simulations. The scheme is computationally efficient as the segmented X-ray image is divided into subdomains before conducting the lattice-Boltzmann simulations, enabling fast simulations. To test the proposed method, it was applied to two synthetic cases with known wettability and three datasets of imaged fluid distributions. The wettability map was obtained for all samples using local lattice-Boltzmann calculations on trapped ganglia and optimization on surface affinity parameters. The results were quantitatively compared with a previously developed geometrical contact angle determination method. The two synthetic cases were used to validate the results of the developed workflow, as well as to compare the wettability results with the geometrical analysis method. It is shown that the developed workflow accurately characterizes the wetting state in the synthetic porous media with an acceptable uncertainty and is better to capture extreme wetting conditions. For the three datasets of imaged fluid distributions, our results show that the obtained contact angle distributions are consistent with the geometrical method. However, the obtained contact angle distributions tend to have a narrower span and are considered more realistic compared to the geometrical method. Finally, our results show the potential of the proposed scheme to efficiently obtain wettability maps of porous media using X-ray images of multiphase fluid distributions. The developed workflow can help for more accurate characterization of the wettability map in the porous media using limited experimental data, and hence more accurate digital rock analysis of multiphase flow in porous media.more » « less
-
McCartney, J.S.; Tomac, I. (Ed.)Multiphase flow patterns in porous media largely depend on the properties of the fluids and interfaces such as viscosity, surface tension, and contact angle. Microorganisms in soils change the fluid and interfacial properties, and thus can alter multiphase fluid flow in porous media. This study investigates the impact of motile bacterium Escherichia coli ( E. coli ) on fluid displacement patterns in a microfluidic chip. The fluid displacement is observed during the saturation and the desaturation processes of the microfluidic chip with and without E.coli suspension. Time-lapse photography results show that the presence of E.coli alters the displacement patterns during the wetting and drying process and changes the residual saturation of the chip. Although studies of the impacts of motility on interfacial properties remain elusive, these results bring the expectation to the manipulation of multiphase transport in porous media and the adaptive control of industrial and environmental flow processes using active particles.more » « less
-
This study aims to bridge length scales in immiscible multiphase flow simulation by connecting two published governing equations at the pore-scale and continuum-scale through a novel validation framework. We employ Niessner and Hassnaizadeh's [“A model for two-phase flow in porous media including fluid-fluid interfacial area,” Water Resour. Res. 44(8), W08439 (2008)] continuum-scale model for multiphase flow in porous media, combined with the geometric equation of state of McClure et al. [“Modeling geometric state for fluids in porous media: Evolution of the Euler characteristic,” Transp. Porous Med. 133(2), 229–250 (2020)]. Pore-scale fluid configurations simulated with the lattice-Boltzmann method are used to validate the continuum-scale results. We propose a mapping from the continuum-scale to pore-scale utilizing a generalized additive model to predict non-wetting phase Euler characteristics during imbibition, effectively bridging the continuum-to-pore length scale gap. Continuum-scale simulated measures of specific interfacial area, saturation, and capillary pressure are directly compared to up-scaled pore-scale simulation results. This research develops a numerical framework capable of capturing multiscale flow equations establishing a connection between pore-scale and continuum-scale simulations.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1038/s41524-021-00598-2
