A CFD-DEM Eulerian-Lagrangian solver for particle-laden viscoelastic flows (for oral presentation)

Abstract

The ability to simulate the behavior of dense suspensions using computationally-efficient Eulerian-Lagrangian techniques requires accurate particulate-phase drag models that are valid for a wide range of material parameters. The present work aims at developing appropriate drag models for moderately-dense suspensions, in which the continuous phase also has viscoelastic characteristics. To this end, we parametrize the suspension properties through the Deborah number and the particle volume fraction, and compute the evolution in the drag coefficient of spheres translating through a viscoelastic fluid that is described by the Oldroyd-B model. To calculate the drag coefficient, we resort to 3D direct numerical simulations (DNS) of unconfined viscoelastic creeping flows (Re < 0.1) past random arrays of stationary spheres, over a wide range of Deborah numbers (De < 5), volume fractions (φ < 20%) and particle configurations. From these calculations we obtain a closure law F(De, φ) for the drag force in the viscoelastic fluid (with fixed retardation ratio 𝛽𝛽 = 0.5), which is on average within 4.7% of the DNS results. Subsequently, this closure law was incorporated into a CFD-DEM Eulerian-Lagrangian solver to handle particle-laden viscoelastic flow calculations, and two case studies were simulated to assess the accuracy and robustness of our numerical approach. These tests consisted of simulating the settling process in Newtonian and viscoelastic fluids within eccentric annular pipes and rectangular channels; configurations commonly employed in hydraulic fracturing operations. The numerical results obtained were found to be in good agreement with experimental data available for suspensions in Newtonian matrix fluids. For the case of viscoelastic fluids, the resulting particle distribution is presented for different elasticity numbers (i.e., El = De/Re) and particle volume fractions, and the results provide insight into the pronounced effects of viscoelastic matrix fluids in hydraulic fracturing operations