Search

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. [Copyright &y& Elsevier]

Numerical simulation of fluid structure interaction (FSI) problems is one of the most challenging topics in computational fluid dynamics. In this paper, coupling edge-based smoothed finite element method (ES-FEM) and smoothed particle hydrodynamics (SPH) method (ES-FEM-SPH) is proposed for solving FSI problems, where the edge-based smoothed finite element method is used to model the movement and deformation of structures, and the smoothed particle hydrodynamics is used to model the fluid flow. In ES-FEM, the gradient smoothing technique is applied over the smoothing domain and it can effectively overcome the "overly-stiff" effect in conventional FEM model. Some correction algorithms including density correction, kernel gradient correction and particle shift technique are integrated into the SPH method to improve computational stability and accuracy. A virtual particle coupling scheme is used to implement the coupling of ES-FEM and SPH with complex geometry interface. As ES-FEM is more accurate than conventional FEM, and it is expected that this ES-FEM-SPH coupling approach should be superior than existing FEM-SPH coupling approaches. A number of test examples with FSI are investigated with the presented ES-FEM-SPH, and compared with results from other approaches including FEM-SPH. From the obtained numerical results, we can conclude that the ES-FEM-SPH coupling approach is effective to simulate FSI problems. • Coupling edge-based smoothed finite element method (ES-FEM) with SPH is proposed for solving FSI problems. • The coupling strategy is based on the automatically generated virtual particles. • The gradient smoothing technique in ES-FEM overcomes the "overly-stiff" effect in conventional FEM. • The SPH method with integrated algorithms is effective and robust in modeling fluid flows. • ES-FEM-SPH is demonstrated to be effective in modeling diversified FSI problems. [ABSTRACT FROM AUTHOR]