Abstract

In numerical simulation of hydraulic fracture propagation, tangent component of the fluid velocity generally considered to be neglected near the crack front. Then Reynolds transport theorem yields that the limit of the particle velocity coincides with the vector of the front propagation speed. We use this fact in combination with the Poiseuille-type equation, which implies that the particle velocity is always collinear to pressure gradient. We show that this specific feature of the hydraulic fracture problem may serve to simplify tracing the front propagation. The latter may be traced without explicit evaluation of the normal to the front, which is needed in conventional applications of the theory of propagating interfaces. Numerical experiments confirm that, despite huge errors in pressure and even greater errors in its gradient, the propagation speed, statistically averaged over a distance of a mesh size, is found quite accurate. We conclude that suggested method may simplify numerical simulation of hydraulic fractures driven by Newtonian and non-Newtonian fluids.

Highlights

  • Numerical modeling is of significance for hydraulic fracture (HF) design [13]

  • This leads to the necessity of the normal to the fracture front evaluation to find the velocity vector entering the discretized continuity equation

  • The numerical results show that the method avoiding explicit evaluation of the normal may be used in practical calculations

Read more

Summary

Introduction

The mathematical formulation of the problem, used for modeling, includes the continuity equation and the movement equation of the Poiseuille type [14, 15] Combined, they contain second spatial derivatives of the pressure. Errors in the pressure gradient, which has non-integrable singularity, are naturally even greater; and they further aggravate when evaluating the second derivatives of the pressure For this reason, conventional schemes of HF modeling avoid evaluation of the pressure in tip elements. This leads to the necessity of the normal to the fracture front evaluation to find the velocity vector entering the discretized continuity equation In its turn, this leads to the need for using conventional methods of the theory of propagating interfaces [2, 16].

Problem Formulation
Numerical results for 1D bench-mark problem
Numerical Results for Axisymmetric Benchmark Problem
Numerical Results for the Benchmark Problem with Stress contrast
Conclusions
Full Text
Paper version not known

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.