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

Title: Pairwise Force Smoothed Particle Hydrodynamics model for multiphase flow: Surface tension and contact line dynamics

Abstract

We present a novel formulation of the Pairwise Force Smoothed Particle Hydrodynamics Model (PF-SPH) and use it to simulate two- and three-phase flows in bounded domains. In the PF-SPH model, the Navier-Stokes equations are discretized with the Smoothed Particle Hydrodynamics (SPH) method and the Young-Laplace boundary condition at the fluid-fluid interface and the Young boundary condition at the fluid-fluid-solid interface are replaced with pairwise forces added into the Navier-Stokes equations. We derive a relationship between the parameters in the pairwise forces and the surface tension and static contact angle. Next, we demonstrate the accuracy of the model under static and dynamic conditions. Finally, to demonstrate the capabilities and robustness of the model we use it to simulate flow of three fluids in a porous material.

Authors:
;
Publication Date:
Research Org.:
Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1233759
Report Number(s):
PNNL-SA-107016
Journal ID: ISSN 0021-9991; KJ0401000
DOE Contract Number:
AC05-76RL01830
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 305
Country of Publication:
United States
Language:
English
Subject:
multiphase flow, smoothed particle hydrodynamics, multiscale modeling, contact angle

Citation Formats

Tartakovsky, Alexandre M., and Panchenko, Alexander. Pairwise Force Smoothed Particle Hydrodynamics model for multiphase flow: Surface tension and contact line dynamics. United States: N. p., 2016. Web. doi:10.1016/j.jcp.2015.08.037.
Tartakovsky, Alexandre M., & Panchenko, Alexander. Pairwise Force Smoothed Particle Hydrodynamics model for multiphase flow: Surface tension and contact line dynamics. United States. doi:10.1016/j.jcp.2015.08.037.
Tartakovsky, Alexandre M., and Panchenko, Alexander. 2016. "Pairwise Force Smoothed Particle Hydrodynamics model for multiphase flow: Surface tension and contact line dynamics". United States. doi:10.1016/j.jcp.2015.08.037.
@article{osti_1233759,
title = {Pairwise Force Smoothed Particle Hydrodynamics model for multiphase flow: Surface tension and contact line dynamics},
author = {Tartakovsky, Alexandre M. and Panchenko, Alexander},
abstractNote = {We present a novel formulation of the Pairwise Force Smoothed Particle Hydrodynamics Model (PF-SPH) and use it to simulate two- and three-phase flows in bounded domains. In the PF-SPH model, the Navier-Stokes equations are discretized with the Smoothed Particle Hydrodynamics (SPH) method and the Young-Laplace boundary condition at the fluid-fluid interface and the Young boundary condition at the fluid-fluid-solid interface are replaced with pairwise forces added into the Navier-Stokes equations. We derive a relationship between the parameters in the pairwise forces and the surface tension and static contact angle. Next, we demonstrate the accuracy of the model under static and dynamic conditions. Finally, to demonstrate the capabilities and robustness of the model we use it to simulate flow of three fluids in a porous material.},
doi = {10.1016/j.jcp.2015.08.037},
journal = {Journal of Computational Physics},
number = ,
volume = 305,
place = {United States},
year = 2016,
month = 1
}
  • Cited by 4
  • A two-dimensional numerical model based on smoothed particle hydrodynamics (SPH) was used to simulate unsaturated (multiphase) flow through fracture junctions. A combination of standard SPH equations with pair-wise fluid-fluid and fluid-solid particle-particle interactions allowed surface tension and three-phase contact dynamics to be simulated. The validity of the model was verified by calculating the surface tension in four different ways: 1. from small amplitude oscillations of fluid drops; 2. from the dependence of the capillary pressure on drop radius; 3. from capillary rise simulations; and 4. from the behavior of a fluid drop confined between parallel walls under the influence ofmore » gravity. All four simulations led to consistent values for the surface tension.. The dependency of receding and advancing contact angles on droplet velocity was studied. Incorporation of surface tension and fluid-solid interactions allowed unsaturated flow through fracture junctions to be realistically simulated, and the simulation results compare well with the laboratory experiments of Dragila, and Weisbrod.« less
  • A two-dimensional numerical model based on smoothed particle hydrodynamics (SPH) was used to simulate unsaturated (multiphase) flow through fracture junctions. A combination of standard SPH equations with pairwise fluid-fluid and fluid-solid particle-particle interactions allowed surface tension and three-phase contact dynamics to be simulated. The model was validated by calculating the surface tension in four different ways: (i) from small-amplitude oscillations of fluid drops, (ii) from the dependence of the capillary pressure on drop radius, (iii) from capillary rise simulations, and (iv) from the behavior of a fluid drop confined between parallel walls under the influence of gravity. All four simulationsmore » led to consistent values for the surface tension. The dependence of receding and advancing contact angles on droplet velocity was studied. Incorporation of surface tension and fluid-solid interactions allowed unsaturated flow through fracture junctions to be realistically simulated, and the simulation results compare well with the laboratory experiments of Dragila and Weisbrod.« less
  • 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.