Simulation of hydraulic fracturing processes in rocks by coupling the lattice Boltzmann model and the Park-Paulino-Roesler potential-based cohesive zone model

Author(s):  
Luis A. Mejía Camones ◽  
Eurípedes do A. Vargas ◽  
Raquel Q. Velloso ◽  
Glaucio H. Paulino
2021 ◽  
Author(s):  
Chang Huang ◽  
Shengli Chen

Abstract The difficulty of hydraulic fracturing in organic-rich shale caused by the increased ductility has not been well interpreted quantitatively, although it is well perceived that the increased shale ductility can impede the propagation of hydraulic fractures and enhance the healing of created fractures upon injection shutdown. This study aims to quantitatively study the impacts of increased ductility on the stimulated reservoir volume (SRV) using an advanced XFEM-based simulator. To achieve this goal, a modified cohesive zone model has been integrated into an in-house fully coupled poroelastic XFEM framework. The study continues by extending the functionality of the numerical framework to simulating multiple interacting fractures. The utilization of the object-oriented programming paradigm in the development of the framework makes it an easy extension to include the multi-fracture network by creating more instances of crack segments. A main hydraulic fracture with an arbitrary number of intersected branches can thus be modeled. A series of parametric studies will be conducted to investigate the impacts of increased ductility on the induced SRV by varying four involved material parameters individually. The modified cohesive zone model, which is essentially a traction-separation law (TSL), is characterized by four parameters: the initial tensile strength Tini, ultimate tensile strength Tkrg, the critical separation Dc, and the final crack separation Dmax. It can flexibly model different crack opening scenarios and simulate more realistically the increased shale ductility. The fully coupled poroelastic XFEM framework has been comprehensively verified against the latest semi-analytical solutions on the four well-known propagation regimes. The numerical results show that the shape of TSL does affect the main hydraulic fracture growth as well as the evolvement of the fracture network, given the same cohesive crack energy and tensile strength. It infers that ductility is not only controlled by cohesive crack energy and tensile strength, which further indicates the necessity of the newly proposed cohesive zone model. The magnitude of the initial tensile strength, controlling when the cohesive crack starts propagating, is found to have the greatest impacts on the fracture length, and SRV, among all four TSL parameters. The novelty of this study is two-fold. First, the newly modified cohesive zone model can more realistically represent the increased shale ductility. Second, the advanced XFEM framework that enables the simulation of a fracture network can study the impacts of increased ductility on the whole SRV but not a single crack.


2020 ◽  
Vol 10 (19) ◽  
pp. 6640
Author(s):  
Zhonghua Shi ◽  
Zhenhang Kang ◽  
Qiang Xie ◽  
Yuan Tian ◽  
Yueqing Zhao ◽  
...  

An effective deicing system is needed to be designed to conveniently remove ice from the surfaces of structures. In this paper, an ultrasonic deicing system for different configurations was estimated and verified based on finite element simulations. The research focused on deicing efficiency factor (DEF) discussions, prediction, and validations. Firstly, seven different configurations of Lead zirconate titanate (PZT) disk actuators with the same volume but different radius and thickness were adopted to conduct harmonic analysis. The effects of PZT shape on shear stresses and optimal frequencies were obtained. Simultaneously, the average shear stresses at the ice/substrate interface and total energy density needed for deicing were calculated. Then, a coefficient named deicing efficiency factor (DEF) was proposed to estimate deicing efficiency. Based on these results, the optimized configuration and deicing frequency are given. Furthermore, four different icing cases for the optimize configuration were studied to further verify the rationality of DEF. The effects of shear stress distributions on deicing efficiency were also analyzed. At same time, a cohesive zone model (CZM) was introduced to describe interface behavior of the plate and ice layer. Standard-explicit co-simulation was utilized to model the wave propagation and ice layer delamination process. Finally, the deicing experiments were carried out to validate the feasibility and correctness of the deicing system.


2021 ◽  
Vol 11 (1) ◽  
pp. 456
Author(s):  
Yanglong Zhong ◽  
Liang Gao ◽  
Xiaopei Cai ◽  
Bolun An ◽  
Zhihan Zhang ◽  
...  

The interface crack of a slab track is a fracture of mixed-mode that experiences a complex loading–unloading–reloading process. A reasonable simulation of the interaction between the layers of slab tracks is the key to studying the interface crack. However, the existing models of interface disease of slab track have problems, such as the stress oscillation of the crack tip and self-repairing, which do not simulate the mixed mode of interface cracks accurately. Aiming at these shortcomings, we propose an improved cohesive zone model combined with an unloading/reloading relationship based on the original Park–Paulino–Roesler (PPR) model in this paper. It is shown that the improved model guaranteed the consistency of the cohesive constitutive model and described the mixed-mode fracture better. This conclusion is based on the assessment of work-of-separation and the simulation of the mixed-mode bending test. Through the test of loading, unloading, and reloading, we observed that the improved unloading/reloading relationship effectively eliminated the issue of self-repairing and preserved all essential features. The proposed model provides a tool for the study of interface cracking mechanism of ballastless tracks and theoretical guidance for the monitoring, maintenance, and repair of layer defects, such as interfacial cracks and slab arches.


2021 ◽  
Vol 33 (3) ◽  
pp. 036607
Author(s):  
Y. Feng ◽  
J. Miranda-Fuentes ◽  
J. Jacob ◽  
P. Sagaut

Sign in / Sign up

Export Citation Format

Share Document