Bulletin of the American Physical Society
70th Annual Meeting of the APS Division of Fluid Dynamics
Volume 62, Number 14
Sunday–Tuesday, November 19–21, 2017; Denver, Colorado
Session F35: Porous Media Flows: PoreScale IPorous

Hide Abstracts 
Chair: Mohammad Nomeli, Massachusetts Institute of Technology Room: 301 
Monday, November 20, 2017 8:00AM  8:13AM 
F35.00001: Prediction of velocity distribution in simple 3D porous media Mohammad AlAdwani, Pietro de Anna, Ruben Juanes Fluid flow and particle transport through porous media are determined by the geometry of the host medium itself. Despite the fundamental importance of the velocity distribution in controlling earlytime and latetime transport properties (e.g., early breakthrough and superdiffusive spreading), direct relations linking velocity distribution with the statistics of pore structure in 3D porous media have not been established yet. High velocities are controlled by the formation of channels, while low velocities are dominated by stagnation zones. Recent studies have proposed phenomenological models for the distribution of high velocities including stretched exponential and powerexponential distributions but without an underlying mechanistic or statistical physics theory. Here, we investigate the relationship between the structure of the host medium and the resulting fluid flow in random dense spherical packs. We simulate flow at low Reynolds numbers by solving the Stokes equations and imposing a noslip boundary condition at the boundary of each sphere. We show that the distribution of low velocities in 3D porous media is described by a Gamma distribution, which is robust to variations in the geometry of the porous media. [Preview Abstract] 
Monday, November 20, 2017 8:13AM  8:26AM 
F35.00002: Fast parallel algorithms for flow in porous media Romir Moza, George Biros In porous media flows, the variation in permeability over a large range of length scales gives rise to fluctuations in velocity which in turn lead to phenomena of interest, such as anomalous diffusion. However, due to the complex 3D structure of the medium, numerical simulation of flows in porous media remains a challenge and there is a need to develop efficient parallel algorithms for fast computation. We model porous media with a set of arbitrary shapes that are homeomorphic to a sphere. The stokes equations are often valid for these flows, which we solve using the boundary integral method. This reduces the dimensionality and number of unknowns for the problem. We evaluate farfield kernel summations using a parallel fast multipole method. Since the separation between spheres can be arbitrarily small, it is a challenge to resolve the flow near the surface. We use special singular and near singular quadratures to tackle this challenge and achieve high accuracy close to the spheres. We study hydrodynamic properties of the flow, such as dispersion, using the correlated continuous time random walk model. Our solver is highly scalable, solving for velocity fields in arbitrary geometries of several thousand spheres in a few minutes. [Preview Abstract] 
Monday, November 20, 2017 8:26AM  8:39AM 
F35.00003: Effect of particle size distribution on permeability in the randomly packed porous media Bojan Markicevic An answer of how porous medium heterogeneity influences the medium permeability is still inconclusive, where both increase and decrease in the permeability value are reported. A numerical procedure is used to generate a randomly packed porous material consisting of spherical particles. Six different particle size distributions are used including \textit{mono}, \textit{bi} and \textit{three}disperse particles, as well as uniform, normal and \textit{log}normal particle size distribution with the maximum to minimum particle size ratio ranging from three to eight for different distributions. In all six cases, the average particle size is kept the same. For all media generated, the stochastic homogeneity is checked from distribution of three coordinates of particle centers, where uniform distribution of $x$, $y$ and $z$ positions is found. The medium surface area remains essentially constant except for \textit{bi}modal distribution in which medium area decreases, while no changes in the porosity are observed (around 0.36). The fluid flow is solved in such domain, and after checking for the pressure axial linearity, the permeability is calculated from the Darcy law. The permeability comparison reveals that the permeability of the \textit{mono}disperse medium is smallest, and the permeability of all \textit{poly}disperse samples is less than ten percent higher. For \textit{bi}modal particles, the permeability is for a quarter higher compared to the other media which can be explained by volumetric contribution of larger particles and larger passages for fluid flow to take place. [Preview Abstract] 
Monday, November 20, 2017 8:39AM  8:52AM 
F35.00004: Pore size distribution effect on rarefied gas transport in porous media Takuma Hori, Yuta Yoshimoto, Shu Takagi, Ikuya Kinefuchi Gas transport phenomena in porous media are known to strongly influence the performance of devices such as gas separation membranes and fuel cells. Knudsen diffusion is a dominant flow regime in these devices since they have nanoscale pores. Many experiments have shown that these porous media have complex structures and pore size distributions; thus, the diffusion coefficient in these media cannot be easily assessed. Previous studies have reported that the characteristic pore diameter of porous media can be defined in light of the pore size distribution; however, tortuosity factor, which is necessary for the evaluation of diffusion coefficient, is still unknown without gas transport measurements or simulations. Thus, the relation between pore size distributions and tortuosity factors is required to obtain the gas transport properties. We perform numerical simulations to prove the relation between them. Porous media are numerically constructed while satisfying given pore size distributions. Then, the meansquare displacement simulation is performed to obtain the tortuosity factors of the constructed porous media.. [Preview Abstract] 
Monday, November 20, 2017 8:52AM  9:05AM 
F35.00005: Pressure drop for inertial flows in elastic porous media Martin Pauthenet, Alessandro Bottaro, Yohan Davit, Michel Quintard The effect of the porosity and of the elastic properties of anisotropic solid skeletons saturated by a fluid is studied for flows displaying unsteady inertial effects. Insight is achieved by direct numerical simulations of the NavierStokes equations for model porous media, with inclusions which can oscillate with respect to their reference positions because of the presence of a restoring elastic force modeled by a spring. The numerical technique is based on the immersed boundary method, to easily allow for the displacement of pores of arbitrary shapes and dimensions. Solid contacts are anelastic. The parameters examined include the local Reynolds number, $Re_d$, based on the mean velocity through the reference unit cell and the characteristic size of the inclusions, the direction of the macroscopic forcing pressure gradient, the reduced frequency, $f^*$, ratio of the flow frequency to the natural frequency of the springmass system, and the reduced mass, $m^*$, ratio of the solid to the fluid density. Results demonstrate the effect of these parameters, and permit to determine the filtration laws useful for the subsequent macroscopic modeling of these flows through the volume averaged NavierStokes equations. [Preview Abstract] 
Monday, November 20, 2017 9:05AM  9:18AM 
F35.00006: ReducedOrder Direct Numerical Simulation of Solute Transport in Porous Media Yashar Mehmani, Hamdi Tchelepi Porescale models are an important tool for analyzing fluid dynamics in porous materials (e.g., rocks, soils, fuel cells). Current direct numerical simulation (DNS) techniques, while very accurate, are computationally prohibitive for sample sizes that are statistically representative of the porous structure. Reducedorder approaches such as porenetwork models (PNM) aim to approximate the porespace geometry and physics to remedy this problem. Predictions from current techniques, however, have not always been successful. This work focuses on singlephase transport of a passive solute under advectiondominated regimes and delineates the minimum set of approximations that consistently produce accurate PNM predictions. Novel network extraction (discretization) and particle simulation techniques are developed and compared to highfidelity DNS simulations for a wide range of micromodel heterogeneities and a single sphere pack. Moreover, common modeling assumptions in the literature are analyzed and shown that they can lead to firstorder errors under advectiondominated regimes. This work has implications for optimizing material design and operations in manufactured (electrodes) and natural (rocks) porous media pertaining to energy systems. [Preview Abstract] 
Monday, November 20, 2017 9:18AM  9:31AM 
F35.00007: The Mechanism behind Erosive Bursts in Porous Media Robin Jaeger, Miller Mendoza, Hans Herrmann We implemented a new model based on the Lattice Boltzmann method to simulate erosion and deposition in suspension flows through porous media. Using this model we show that the cause of erosive bursts in filtration experiments is the reopening of clogged pores when the pressure difference between two opposite sites of the pore surpasses a certain threshold. We perform numerical simulations and find excellent agreement to experimental results when comparing shape and size distribution of pressure loss jumps, which are the direct result of erosive bursts. Furthermore, we find that erosive bursts only occur for pressure gradient thresholds within the range of two critical values, independent on how the flow is driven. We believe that our findings provide a better understanding of sudden sand production in oil wells and breakthrough in filtration. [Preview Abstract] 
Monday, November 20, 2017 9:31AM  9:44AM 
F35.00008: Lattice Boltzmann simulations of immiscible displacement process with large viscosity ratios Parthib Rao, Laura Schaefer Immiscible displacement is a key physical mechanism involved in enhanced oil recovery and carbon sequestration processes. This multiphase flow phenomenon involves a complex interplay of viscous, capillary, inertial and wettability effects. The lattice Boltzmann (LB) method is an accurate and efficient technique for modeling and simulating multiphase/multicomponent flows especially in complex flow configurations and media. In this presentation we present numerical simulation results of displacement process in thin long channels. The results are based on a new psuedopotential multicomponent LB model with multiple relaxation time collision (MRT) model and explicit forcing scheme. We demonstrate that the proposed model is capable of accurately simulating the displacement process involving fluids with a wider range of viscosity ratios (>100) and which also leads to viscosityindependent interfacial tension and reduction of some important numerical artifacts. [Preview Abstract] 
Monday, November 20, 2017 9:44AM  9:57AM 
F35.00009: Analysis of Nonequilibrium Capillary PressureSaturation Relation using Direct Numerical Simulations with VolumeOfFluid (VOF) Method Santosh Konangi, Nikhil K Palakurthi, Nikolaos Karadimitriou, An Fu, Ken Comer, Urmila Ghia In traditional twophase flow models of porous media, capillary pressure (P$_{\mathrm{c}})$ and saturation (S$^{\mathrm{w}})$ are hysteretically related, i.e. different P$_{\mathrm{c}}$S$^{\mathrm{w}}$ curves are obtained for drainage and imbibition. Extended twophase flow theories hypothesize that inclusion of specific interfacial area (a$^{\mathrm{wn}})$ will result in a unique relation between capillary pressure, saturation and interfacial area (P$_{\mathrm{c}}$S$^{\mathrm{w}}$a$^{\mathrm{wn}})$. Several studies have confirmed the reduction of hysteresis in the P$_{\mathrm{c}}$S$^{\mathrm{w}}$a$^{\mathrm{wn}}$ relation under quasistatic conditions. However, the uniqueness of the P$_{\mathrm{c}}$S$^{\mathrm{w}}$a$^{\mathrm{wn}}$ relation under transient conditions is not clear. We investigate role of specific interfacial area under dynamic conditions using porescale direct numerical simulations (DNS) on two micromodels with volumeoffluid (VOF) method. From the DNS data, PcS$^{\mathrm{w}}$ curves are estimated for drainage and imbibition; validity of different macroscopic capillary pressure definitions is evaluated. The quasistatic and dynamic P$_{\mathrm{c}}$S$^{\mathrm{w}}$a$^{\mathrm{wn}}$ relations are examined for uniqueness. [Preview Abstract] 
Follow Us 
Engage
Become an APS Member 
My APS
Renew Membership 
Information for 
About APSThe American Physical Society (APS) is a nonprofit membership organization working to advance the knowledge of physics. 
© 2024 American Physical Society
 All rights reserved  Terms of Use
 Contact Us
Headquarters
1 Physics Ellipse, College Park, MD 207403844
(301) 2093200
Editorial Office
100 Motor Pkwy, Suite 110, Hauppauge, NY 11788
(631) 5914000
Office of Public Affairs
529 14th St NW, Suite 1050, Washington, D.C. 200452001
(202) 6628700