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

Abstract

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.

Authors:
 [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). E-mail: liumb@ntu.edu.sg
  2. Center for Advanced Modeling and Simulation, Idaho National Laboratory, P.O. Box 1625, MS 2211, Idaho Falls, ID 83415-2211 (United States)
Publication Date:
OSTI Identifier:
20991561
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 222; Journal 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)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; COMPUTERIZED SIMULATION; FRACTURES; GASES; GEOMETRY; HYDRODYNAMICS; LIQUIDS; MATHEMATICAL MODELS; PARTICLE INTERACTIONS; SOLIDS; THIN FILMS; WEIGHTING FUNCTIONS

Citation Formats

Liu Moubin, Meakin, Paul, and Huang Hai. Dissipative particle dynamics simulation of fluid motion through an unsaturated fracture and fracture junction. United States: N. p., 2007. Web.
Liu Moubin, Meakin, Paul, & Huang Hai. Dissipative particle dynamics simulation of fluid motion through an unsaturated fracture and fracture junction. United States.
Liu Moubin, Meakin, Paul, and Huang Hai. Thu . "Dissipative particle dynamics simulation of fluid motion through an unsaturated fracture and fracture junction". United States. doi:.
@article{osti_20991561,
title = {Dissipative particle dynamics simulation of fluid motion through an unsaturated fracture and fracture junction},
author = {Liu Moubin and Meakin, Paul and Huang Hai},
abstractNote = {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.},
doi = {},
journal = {Journal of Computational Physics},
number = 1,
volume = 222,
place = {United States},
year = {Thu Mar 01 00:00:00 EST 2007},
month = {Thu Mar 01 00:00:00 EST 2007}
}
  • 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. Ourmore » 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.« less
  • Multiphase fluid flow through porous media involves complex fluid dynamics, and it is difficult to model such complex behavior, on the pore scale, using grid-based continuum models. In this paper, the application of dissipative particle dynamics (DPD), a relatively new mesoscale method, to the simulation of pore-scale multiphase fluid flows under a variety of flow conditions is described. We demonstrate that the conventional DPD method using purely repulsive conservative (nondissipative) particle-particle interactions is capable of modeling single-phase flow fields in saturated porous media. In order to simulate unsaturated multiphase flow through porous media, we applied a modified model for themore » conservative particle-particle interactions that combines short-range repulsive and long-range attractive interactions. This form for the conservative particle-particle interactions allows the behavior of multiphase systems consisting of gases, liquids, and solids to be simulated. We also demonstrated that the flow of both wetting and nonwetting fluids through porous media can be simulated by controlling the ratios between the fluid-fluid and fluid-solid (fluid-wall) interparticle interaction strengths.« less
  • Complex fluid behavior in unsaturated fracture and fracture networks, such as film flow, the migration, fragmentation, and coalescence of droplets, and rivulet flow with or without meandering or pulsation, has been widely observed in laboratory experiments. In this study, a modified two-dimensional volume of fluid (VOF) method was used to simulate liquid motion in partially saturated fracture apertures under a variety of flow conditions. This modeling approach systematically incorporates the effects of inertial forces, viscosity, gravity acting on the fluid densities, fracture wall wetting, and the pressure drop across curved fluid-fluid interfaces due to surface tension. This allows us tomore » obtain a better understanding of the fundamental physics governing unsaturated fluid flow in fracture apertures. The VOF method is able to handle the complex dynamics of fluid-fluid interfaces and free surfaces in unsaturated fractures by using a fixed Eulerian grid. Fragmentation and coalescence of the fluids are automatically handled without resorting to complex adaptive mesh refinement or interface repairing algorithms. The wetting of fracture walls was modeled by imposing contact angles near the contact lines (contact points in two-dimensional simulations), and different contact angles were automatically chosen depending on whether the liquid interface is advancing, receding, or essentially stationary. The qualitative agreements between the numerical simulations and complex multiphase fluid dynamics reported in laboratory experiments clearly demonstrate the potential value of the VOF method for the mechanistically based modeling of immiscible liquid motion in unsaturated fracture networks.« less
  • We rederive the equations of motion of dissipative relativistic fluid dynamics from kinetic theory. In contrast with the derivation of Israel and Stewart, which considered the second moment of the Boltzmann equation to obtain equations of motion for the dissipative currents, we directly use the latter's definition. Although the equations of motion obtained via the two approaches are formally identical, the coefficients are different. We show that, for the one-dimensional scaling expansion, our method is in better agreement with the solution obtained from the Boltzmann equation.