期刊
ENGINEERING WITH COMPUTERS
卷 -, 期 -, 页码 -出版社
SPRINGER
DOI: 10.1007/s00366-022-01626-5
关键词
Random arrays of spheres; Drag coefficient; Viscoelastic fluids; Oldroyd-B model; Eulerian-Lagrangian formulation; Discrete particle method
资金
- FEDER funds through the COMPETE 2020 Programme
- FCT (Portuguese Foundation for Science and Technology) under MIT Portugal program [UID-B/05256/2020, UID-P/05256/2020, MIT-EXPL/TDI/0038/2019, POCI-01-0145-FEDER-016665]
- University of Minho cluster [NORTE-07-0162-FEDER-000086]
- Minho Advanced Computing Center (MACC) [CPCA_A2_6052_2020]
- PRACE - Partnership for Advanced Computing in Europe [icei-prace-2020-0009]
- Fundação para a Ciência e a Tecnologia [MIT-EXPL/TDI/0038/2019] Funding Source: FCT
Accurately resolving the coupled momentum transfer between the liquid and solid phases of complex fluids is essential in multiphase transport processes. This study uses direct numerical simulations and closure drag laws to investigate the flow characteristics of viscoelastic fluids through static random arrays of spherical particles.
Accurately resolving the coupled momentum transfer between the liquid and solid phases of complex fluids is a fundamental problem in multiphase transport processes, such as hydraulic fracture operations. Specifically we need to characterize the dependence of the normalized average fluid-particle force < F > on the volume fraction of the dispersed solid phase and on the rheology of the complex fluid matrix, parameterized through the Weissenberg number Wi measuring the relative magnitude of elastic to viscous stresses in the fluid. Here we use direct numerical simulations (DNS) to study the creeping flow (Re << 1) of viscoelastic fluids through static random arrays of monodisperse spherical particles using a finite volume Navier-Stokes/ Cauchy momentum solver. The numerical study consists of N = 150 different systems, in which the normalized average fluid-particle force < F > is obtained as a function of the volume fraction phi (0 < phi <= 0.2) of the dispersed solid phase and the Weissenberg number Wi (0 <= Wi <= 4). From these predictions a closure law < F(phi, Wi)> for the drag force is derived for the quasi-linear Oldroyd-B viscoelastic fluid model (with fixed retardation ratio beta = 0.5) which is, on average, within 5.7% of the DNS results. In addition, a flow solver able to couple Eulerian and Lagrangian phases (in which the particulate phase is modeled by the discrete particle method (DPM)) is developed, which incorporates the viscoelastic nature of the continuum phase and the closed-form drag law. Two case studies were simulated using this solver, to assess the accuracy and robustness of the newly developed approach for handling particle-laden viscoelastic flow configurations with O(10(5) - 10(6)) rigid spheres that are representative of hydraulic fracture operations. Three-dimensional settling processes in a Newtonian fluid and in a quasi-linear Oldroyd-B viscoelastic fluid are both investigated using a rectangular channel and an annular pipe domain. Good agreement is obtained for the particle distribution measured in a Newtonian fluid, when comparing numerical results with experimental data. For the cases in which the continuous fluid phase is viscoelastic we compute the evolution in the velocity fields and predicted particle distributions are presented at different elasticity numbers 0 <= El <= 30 (where El = Wi/Re) and for different suspension particle volume fractions.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据