Fluid–particle flow simulations using two-way-coupled mesoscale SPH–DEM and validation

• Published in: International journal of multiphase flow Vol. 59; pp. 121 - 134
• Main Authors: Robinson, Martin, Ramaioli, Marco, Luding, Stefan
• Format: Journal Article
• Language: English
• Published: Kidlington Elsevier Ltd 01-02-2014; Elsevier
• Subjects: Buoyancy-driven instability; Computational fluid dynamics; Computational methods in fluid dynamics; DEM; Discrete element method; Discrete particle model; Exact sciences and technology; Finite element method; Flexibility; Fluid dynamics; Fluid flow; Fluids; Fluid–particle flow; Fundamental areas of phenomenology (including applications); Hydrodynamic stability; Mathematical analysis; Multiphase and particle-laden flows; Navier-Stokes equations; Nonhomogeneous flows; PARDEM; Physics; Rayleigh–Taylor instability; Sedimentation; SPH; Discrete particle model; Fluid–particle flow; PARDEM; SPH; Rayleigh–Taylor instability; DEM; Sedimentation; Digital simulation; Smoothed particle hydrodynamics method; Discrete element method; Particle suspension; Rayleigh-Taylor instability; Two-phase flow; Modelling; Fluid-particle flow
• ISSN: 0301-9322; 1879-3533
• DOI: 10.1016/j.ijmultiphaseflow.2013.11.003

[Display omitted]
•We present a simulation method for fluid–particle flows using coupled SPH and DEM.•Purely particle-based method, no prescribed mesh.•Suitable for problems such as free surface flow, complex and/or moving geometries.•Successfully reproduces expected behaviour in 3D sedimentation test cases.•Limitations for too fine fluid resolutions or high porosity gradients are discussed.
First, a meshless simulation method is presented for multiphase fluid–particle flows with a two-way coupled Smoothed Particle Hydrodynamics (SPH) for the fluid and the Discrete Element Method (DEM) for the solid phase. The unresolved fluid model, based on the locally averaged Navier Stokes equations, is expected to be considerably faster than fully resolved models. Furthermore, in contrast to similar mesh-based Discrete Particle Models (DPMs), our purely particle-based method enjoys the flexibility that comes from the lack of a prescribed mesh. It is suitable for problems such as free surface flow or flow around complex, moving and/or intermeshed geometries and is applicable to both dilute and dense particle flows.
Second, a comprehensive validation procedure for fluid–particle simulations is presented and applied here to the SPH–DEM method, using simulations of single and multiple particle sedimentation in a 3D fluid column and comparison with analytical models. Millimetre-sized particles are used along with three different test fluids: air, water and a water–glycerol solution. The velocity evolution for a single particle compares well (less than 1% error) with the analytical solution as long as the fluid resolution is coarser than two times the particle diameter. Two more complex multiple particle sedimentation problems (sedimentation of a homogeneous porous block and an inhomogeneous Rayleigh Taylor Instability) are also reproduced well for porosities 0.6⩽∊⩽1.0, although care should be taken in the presence of high porosity gradients.
Overall the SPH–DEM method successfully reproduces quantitatively the expected behaviour in these test cases, and promises to be a flexible and accurate tool for other, realistic fluid–particle system simulations (for which other problem-relevant test cases have to be added for validation).