CFD-DEM modeling of fracture initiation with polymer injection in granular media--颗粒学报PARTICUOLOGY

Kazidenov, D., & Amanbek, Y. (2025). CFD-DEM modeling of fracture initiation with polymer injection in granular media. Particuology, 97, 58-68. https://doi.org/10.1016/j.partic.2024.12.001

CFD-DEM modeling of fracture initiation with polymer injection in granular media

Daniyar Kazidenov, Yerlan Amanbek^*

Department of Mathematics, Nazarbayev University, Kabanbay batyr 53, Astana, 01000, Kazakhstan

10.1016/j.partic.2024.12.001

Volume 97, February 2025, Pages 58-68

Received 12 October 2024, Revised 26 November 2024, Accepted 3 December 2024, Available online 20 December 2024, Version of Record 8 January 2025.

E-mail: yerlan.amanbek@nu.edu.kz

Highlights

• Power-law drag force-based CFD-DEM is used to study fracture initiation induced by shear-thinning polymer solutions.

• Numerical fracture initiation model is successfully validated against a real laboratory experiment.

• Impact of various factors, including fluid and solid material properties, on fracture initiation is analyzed.

• A fracture initiation model based on dimensionless parameters is proposed.

Abstract

We numerically study the mechanisms and conditions for fracture initiation in weakly cohesive granular media induced by non-Newtonian polymer solutions. A coupled computational fluid dynamics–discrete element method (CFD-DEM) approach is utilized to model fluid flow in a porous medium. The flow behavior of polymer solutions and the drag force acting on particles are calculated using a power-law model. The adequacy of the numerical model is confirmed by comparing the results with a laboratory experiment. The numerical results are consistent with the experimental data presenting similar trends in dimensionless parameters that incorporate fluid flow rate, rheology, peak pressure, and confining stress. The results show that fluid flow rate, rheology, and solid material characteristics strongly influence fracture initiation behavior. Injection of a more viscous guar-based solution results in wider fractures induced by grain displacement, whereas a less viscous XG-based solution creates more linear fractures dominated by infiltration. The ratio of peak pressures between two fluids is higher in the rigid material than in the softer material. Finally, the dimensionless parameters 1/Π₁ and τ₂, which account for fluid and solid material properties accordingly, are effective indicators in determining fracture initiation induced by shear-thinning fluids. Our numerical results show that fracture initiation occurs above 1/Π₁ = 0.06 and τ₂ = 2 ⋅ 10⁻⁷.

Graphical abstract

Keywords

Fracture initiation; Fracture propagation; Non-Newtonian fluid; Power-law model

文章导读