skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Smoothed Particle Hydrodynamics and its applications for multiphase flow and reactive transport in porous media

Abstract

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian method based on a meshless discretization of partial differential equations. In this review, we present SPH discretization of the Navier-Stokes and Advection-Diffusion-Reaction equations, implementation of various boundary conditions, and time integration of the SPH equations, and we discuss applications of the SPH method for modeling pore-scale multiphase flows and reactive transport in porous and fractured media.

Authors:
; ; ; ; ;
Publication Date:
Research Org.:
Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1334861
Report Number(s):
PNNL-SA-107017
Journal ID: ISSN 1420-0597; KJ0401000
DOE Contract Number:
AC05-76RL01830
Resource Type:
Journal Article
Resource Relation:
Journal Name: Computational Geosciences; Journal Volume: 20; Journal Issue: 4
Country of Publication:
United States
Language:
English

Citation Formats

Tartakovsky, Alexandre M., Trask, Nathaniel, Pan, K., Jones, Bruce D., Pan, Wenxiao, and Williams, John R. Smoothed Particle Hydrodynamics and its applications for multiphase flow and reactive transport in porous media. United States: N. p., 2016. Web. doi:10.1007/s10596-015-9468-9.
Tartakovsky, Alexandre M., Trask, Nathaniel, Pan, K., Jones, Bruce D., Pan, Wenxiao, & Williams, John R. Smoothed Particle Hydrodynamics and its applications for multiphase flow and reactive transport in porous media. United States. doi:10.1007/s10596-015-9468-9.
Tartakovsky, Alexandre M., Trask, Nathaniel, Pan, K., Jones, Bruce D., Pan, Wenxiao, and Williams, John R. 2016. "Smoothed Particle Hydrodynamics and its applications for multiphase flow and reactive transport in porous media". United States. doi:10.1007/s10596-015-9468-9.
@article{osti_1334861,
title = {Smoothed Particle Hydrodynamics and its applications for multiphase flow and reactive transport in porous media},
author = {Tartakovsky, Alexandre M. and Trask, Nathaniel and Pan, K. and Jones, Bruce D. and Pan, Wenxiao and Williams, John R.},
abstractNote = {Smoothed Particle Hydrodynamics (SPH) is a Lagrangian method based on a meshless discretization of partial differential equations. In this review, we present SPH discretization of the Navier-Stokes and Advection-Diffusion-Reaction equations, implementation of various boundary conditions, and time integration of the SPH equations, and we discuss applications of the SPH method for modeling pore-scale multiphase flows and reactive transport in porous and fractured media.},
doi = {10.1007/s10596-015-9468-9},
journal = {Computational Geosciences},
number = 4,
volume = 20,
place = {United States},
year = 2016,
month = 3
}
  • A numerical model based on smoothed particle hydrodynamics (SPH) was used to simulate reactive transport and mineral precipitation in porous and fractured porous media. The model was used to study effects of the Damkohler and Peclet numbers and pore-scale heterogeneity on reactive transport and the character of mineral precipitation, and to estimate effective reaction coefficients and mass transfer coefficients. The changes in porosity, fluid and solute fluxes and transport parameters resulting from mineral precipitation were also investigated. The simulation results show that the SPH, Lagrangian particle method, is an effective tool for studying pore scale flow and transport. The particlemore » nature of SPH allows complex physical processes such as diffusion, reaction and mineral precipitation to be modeled with relative ease.« less
  • The development of a framework to support smoothed particle hydrodynamics (SPH) simulations of fluid flow and transport in porous media is described. The framework is built using the Common Component Architecture (CCA) toolkit and supports SPH simulations using a variety of different SPH models and setup formats. The SPH simulation code is decomposed into independent components that represent self-contained units of functionality. Different physics models can be developed within the framework by re-implementing key components but no modification of other components is required. The model for defining components and developing abstract interfaces for them that support a high degree ofmore » modularity and minimal dependencies between components is discussed in detail.« less
  • We have conducted a series of high-resolution numerical experiments using the Pair-Wise Force Smoothed Particle Hydrodynamics (PF-SPH) multiphase flow model. First, we derived analytical expressions relating parameters in the PF-SPH model to the surface tension and static contact angle. Next, we used the model to study viscous fingering, capillary fingering, and stable displacement of immiscible fluids in porous media for a wide range of capillary numbers and viscosity ratios. We demonstrated that the steady state saturation profiles and the boundaries of viscous fingering, capillary fingering, and stable displacement regions compare favorably with micromodel laboratory experimental results. For displacing fluid withmore » low viscosity, we observed that the displacement pattern changes from viscous fingering to stable displacement with increasing injection rate. When a high viscosity fluid is injected, transition behavior from capillary fingering to stable displacement occurred as the flow rate was increased. These observation also agree with the results of the micromodel laboratory experiments.« less
  • A numerical model based on smoothed particle hydrodynamics (SPH) was used to simulate pore-scale liquid and gas flow in synthetic two-dimensional porous media consisting of non-overlapping grains. The model was used to study effects of pore scale heterogeneity and anisotropy on relationship between the average saturation and the Bond number. The effect of the wetting fluid properties on drainage was also investigated. It is shown that pore-scale heterogeneity and anisotropy can cause saturation/Bond number and entry (bubbling) pressures to be dependent on the flow direction suggesting that these properties should be described by tensor rather than scalar quantities.
  • A meso-scale stochastic Lagrangian particle model was developed and used to simulate conservative and reactive transport in porous media. In the stochastic model, the fluid flow in a porous continuum is governed by a combination of a Langevin equation and continuity equation. Pore-scale velocity fluctuations, the source of hydrodynamic dispersion, are represented by the white noise. A smoothed particle hydrodynamics method was used to solve the governing equations. Changes in the properties of the fluid particles (e.g., the solute concentration) are governed by the advection-diffusion equation. The separate treatment of advective and diffusive mixing in the stochastic transport model ismore » more realistic than the classical advection-dispersion theory, which uses a single effective diffusion coefficient (the dispersion coefficient) to describe both types of mixing leading to over-prediction of mixing induced effective reaction rates. The stochastic model predicts much lower reaction product concentrations in mixing induced reactions. In addition, the dispersion theory predicts more stable fronts (with a higher effective fractal dimension) than the stochastic model during the growth of Rayleigh-Taylor instabilities.« less