En
期刊文章
文章精选
当期目录
过刊浏览
高被引
高下载
特邀文章
专刊合集
待刊文章
Volume 105
您当前的位置：首页 > 期刊文章 > 过刊浏览 > Volumes 96-107 (2025) > Volume 105
CFD-DEM modeling of fluid-driven fracture induced by temperature-dependent polymer injection
Daniyar Kazidenov a, Sagyn Omirbekov b, Yerlan Amanbek a *
a Department of Mathematics, Nazarbayev University, Astana, Kazakhstan
b National Laboratory Astana, Nazarbayev University, Astana, Kazakhstan
10.1016/j.partic.2025.08.002
Volume 105, October 2025, Pages 259-276
Received 18 June 2025, Revised 29 July 2025, Accepted 2 August 2025, Available online 18 August 2025, Version of Record 30 August 2025.
E-mail: yerlan.amanbek@nu.edu.kz
Highlights

• A coupled CFD-DEM model is used to simulate fracturing in porous media with temperature-dependent non-Newtonian polymer flow.

• A fluid-particle heat transfer model with temperature-dependent power-law rheology is implemented and verified analytically.

• Effects of polymer concentration and thermal conditions on fracture initiation and propagation behavior are analyzed.

• A fracture initiation criterion using dimensionless parameters predicts fracturing under varying thermal conditions.


Abstract

This study numerically investigates fracture initiation and propagation during polymer-based solution injection under varying thermal conditions. A coupled computational fluid dynamics and discrete element method (CFD-DEM) framework is used to model non-Newtonian fluid flow through a granular medium. The rheology of shear-thinning fluids and fluid-particle heat transfer are modeled with temperature-dependent power-law parameters. The current model is validated by comparing fracture propagation behavior and peak pressures against the similar numerical study. The adequacy of the fluid-particle heat transfer model is confirmed by comparing the results with an analytical approach. The simulation results show that polymer concentration significantly influences fracturing behavior. Less concentrated, lower-viscosity fluids are more likely to create linear fracture paths with enhanced fluid infiltration. In contrast, fluids with higher polymer concentrations and viscosities tend to produce wider fractures characterized by greater particle displacement. An increase in the fluid temperature injected into the cooler medium leads to a reduction of fracture size for the 0.4 % (w/w) XG solution, while the 0.6 % (w/w) XG solution tends to form more linear fracture tips. At sufficiently elevated medium temperatures, the injection of cooler fluids prevents fracture initiation for both concentrations. Lower-viscosity cases, dominated by infiltration, reflect broader thermal transitions in particle temperature distribution, whereas higher-viscosity cases, characterized by particle displacement, exhibit narrower transition regions along fracture boundaries. A fracture initiation criterion for shear-thinning fluids is proposed based on the dimensionless parameters Π1 and τ2. Fracture occurs when Π1 > 73 and τ2 > 3.58 × 10−9. The 0.4 % solution exhibits lower thermal sensitivity with relatively minimal variations in the dimensionless parameters, while the 0.6 % solution shows a greater response to temperature changes, reflected in broader variations of these parameters.

Graphical abstract
Keywords
Fracture initiation; Non-Newtonian fluid; Temperature-dependent rheology; Thermal fracturing; CFD-DEM modeling; Polymer flooding
文章导读