(Received ?; revised ?; accepted ?. -To be entered by editorial office)During improved oil recovery, gas may be introduced into a porous reservoir filled with surfactant solution in order to form foam. A model for the evolution of the resulting foam front known as 'pressure-driven growth' is analysed. An asymptotic solution of this model for long times is derived that shows that foam can propagate indefinitely into the reservoir without gravity override. Moreover 'pressure-driven growth' is shown to correspond to a special case of the more general 'viscous froth' model. In particular, it is a singular limit of the viscous froth, corresponding to the elimination of a surface tension term, permitting sharp corners and kinks in the predicted shape of the front. Sharp corners tend to develop from concave regions of the front. The principal solution of interest has a convex front, however, so that although this solution itself has no sharp corners (except for some kinks that develop spuriously owing to errors in a numerical scheme), it is found nevertheless to exhibit milder singularities in front curvature, as the long-time asymptotic analytical solution makes clear. Numerical schemes for the evolving front shape which perform robustly (avoiding the development of spurious kinks) are also developed. Generalisations of this solution to geologically heterogeneous reservoirs should exhibit concavities and/or sharp corner singularities as an inherent part of their evolution: propagation of fronts containing such 'inherent' singularities can be readily incorporated into these numerical schemes.

Wrobel, M., Mishuris, G. (2015). Hydraulic fracture revisited: Particle velocity based simulation. International Journal of Engineering Science, 94, 23-58.We develop a new effective mathematical formulation and resulting universal computational algorithm capable of tackling various HF models in the framework of a unified approach. The scheme is not limited to any particular elasticity operator or crack propagation regime. Its basic assumptions are: (i) proper choice of independent and dependent variables (with the direct utilization of a new one ? the reduced particle velocity), (ii) tracing the fracture front by use of the Stefan condition (speed equation), which can be integrated in closed form and provides an explicit relation between the crack propagation speed and the coefficients in the asymptotic expansion of the crack opening, (iii) proper regularization techniques, (iv) improved temporal approximation, (v) modular algorithm architecture. The application of the new dependent variable, the reduced particle velocity, instead of the usual fluid flow rate, facilitates the computation of the crack propagation speed from the local relation based on the speed equation. This way, we avoid numerical evaluation of the undetermined limit of the product of fracture aperture and pressure gradient at the crack tip (or alternatively the limit resulting from ratio of the fluid flow rate and the crack opening), which always poses a considerable computational challenge. As a result, the position of the crack front is accurately determined from an explicit formula derived from the speed equation. This approach leads to a robust numerical scheme. Its performance is demonstrated using classical examples of 1D models for hydraulic fracturing: PKN and KGD models under various fracture propagation regimes. Solution accuracy is verified against dedicated analytical benchmarks and other solutions available in the literature. The scheme can be directly extended to more general 2D and 3D cases.publishersversionPeer reviewe

Piccolroaz, A; Mishuris, G; Movchan AB. Symmetric and skew-symmetric weight functions in 2D perturbation models for semi-infinite interfacial cracks. Journal of the Mechanics and Physics of Solids. 2009, 57(9), 1657-1682In this paper we address the vector problem of a 2D half-plane interfacial crack loaded by a general asymmetric distribution of forces acting on its faces. It is shown that the general integral formula for the evaluation of stress intensity factors, as well as high-order terms, requires both symmetric and skew-symmetric weight function matrices. The symmetric weight function matrix is obtained via the solution of a Wiener-Hopf functional equation, whereas the derivation of the skew-symmetric weight function matrix requires the construction of the corresponding full field singular solution. The weight function matrices are then used in the perturbation analysis of a crack advancing quasi-statically along the interface between two dissimilar media. A general and rigorous asymptotic procedure is developed to compute the perturbations of stress intensity factors as well as high-order terms.Peer reviewe

No abstract

The paper is concerned with the problem of a semi-infinite crack at the interface between two dissimilar elastic half-spaces, loaded by a general asymmetrical system of forces distributed along the crack faces. On the basis of the weight function approach and the fundamental reciprocal identity (Betti formula), we formulate the elasticity problem in terms of singular integral equations relating the applied loading and the resulting crack opening. Such formulation is fundamental in the theory of elasticity and extensively used to solve several problems in linear elastic fracture mechanics (for instance various classic crack problems in homogeneous and heterogeneous media).Peer reviewe

The problem of hydraulic fracture propagation is considered by using its recently suggested modified formulation in terms of the particle velocity, the opening in the proper degree, appropriate spatial coordinates and e-regularization. We show that the formulation may serve for significant increasing the efficiency of numerical tracing the fracture propagation. Its advantages are illustrated by re-visiting the Nordgren problem. It is shown that the modified formulation facilitates (i) possibility to have various stiffness of differential equations resulting after spatial discretization, (ii) obtaining highly accurate and stable numerical results with moderate computational effort, and (iii) sensitivity analysis. The exposition is extensively illustrated by numerical examples.?Peer reviewe

Spontaneous brittle fracture is studied based on the model of a body, recently introduced by two of the authors, where only the prospective crack path is specified as a discrete set of alternating initially stretched and compressed bonds. In such a structure, a bridged crack destroying initially stretched bonds may propagate under a certain level of the internal energy without external sources. The general analytical solution with the crack speed-energy relation is presented in terms of the crack-related dynamic Green's function. For anisotropic chains and lattices considered earlier in quasi-statics, the dynamic problems are examined and discussed in detail. The crack speed is found to grow unboundedly as the energy approaches its upper limit. The steady-state sub-and supersonic regimes found analytically are confirmed by numerical simulations. In addition, irregular growth, clustering and crack speed oscillation modes are detected at a lower bound of the internal energy. It is observed, in numerical simulations, that the spontaneous fracture can occur in the form of a pure bridged, partially bridged or fully open crack depending on the internal energy level.

A novel hydraulic fracture (HF) formulation is introduced which accounts for the hydraulically induced shear stress at the crack faces. It utilizes a general form of the elasticity operator alongside a revised fracture propagation condition based on the critical value of the energy release rate. It is shown that the revised formulation describes the underlying physics of HF in a more accurate way and is in agreement with the asymptotic behaviour of the linear elastic fracture mechanics. A number of numerical simulations by means of the universal HF algorithm previously developed in Wrobel & Mishuris (2015) are performed in order to: i) compare the modified HF formulation with its classic counterpart and ii) investigate the peculiarities of the former. Computational advantages of the revised HF model are demonstrated. Asymptotic estimations of the main solution elements are provided for the cases of small and large toughness. The modified formulation opens new ways to analyse the physical phenomenon of HF and also improves the reliability and efficiency of its numerical simulations.