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

Title: Dissipative particle dynamics simulation of fluid motion through an unsaturated fracture and fracture junction

Journal Article · · Journal of Computational Physics
OSTI ID:20991561
 [1];  [2];  [2]
  1. Center for Advanced Modeling and Simulation, Idaho National Laboratory, P.O. Box 1625, MS 2211, Idaho Falls, ID 83415-2211 (United States) and College of Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798 (Singapore)
  2. Center for Advanced Modeling and Simulation, Idaho National Laboratory, P.O. Box 1625, MS 2211, Idaho Falls, ID 83415-2211 (United States)

Multiphase fluid motion in unsaturated fractures and fracture networks involves complicated fluid dynamics, which is difficult to model using grid-based continuum methods. In this paper, the application of dissipative particle dynamics (DPD), a relatively new mesoscale method to simulate fluid motion in unsaturated fractures is described. Unlike the conventional DPD method that employs a purely repulsive conservative (non-dissipative) particle-particle interaction to simulate the behavior of gases, we used conservative particle-particle interactions that combine short-range repulsive and long-range attractive interactions. This new conservative particle-particle interaction allows the behavior of multiphase systems consisting of gases, liquids and solids to be simulated. Our simulation results demonstrate that, for a fracture with flat parallel walls, the DPD method with the new interaction potential function is able to reproduce the hydrodynamic behavior of fully saturated flow, and various unsaturated flow modes including thin film flow, wetting and non-wetting flow. During simulations of flow through a fracture junction, the fracture junction can be fully or partially saturated depending on the wetting property of the fluid, the injection rate and the geometry of the fracture junction. Flow mode switching from a fully saturated flow to a thin film flow can also be observed in the fracture junction.

OSTI ID:
20991561
Journal Information:
Journal of Computational Physics, Vol. 222, Issue 1; Other Information: DOI: 10.1016/j.jcp.2006.07.017; PII: S0021-9991(06)00320-2; Copyright (c) 2006 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9991
Country of Publication:
United States
Language:
English