Search for: All records

Abstract Dike propagation is an intrinsically multiphase problem, where deformation and fluid flow are intricately coupled in a fracture process. Here we perform the first fully coupled simulations of dike propagation in two dimensions, accounting for depressurization of a circular magma chamber, dynamic fluid flow, fracture formation, and elastic deformation. Despite the complexity of the governing equations, we observe that the lengthening is well explained by a simple model, whereis the dike length,is time, andandare constants. We compare the model to seismic data from eight dikes in Iceland and Ethiopia, and, in spite of the assumption of plane strain, we find good agreement between the data and the model. In addition, we derive an approximate model for the depressurization of the chamber with the dike length. These models may help forecast the growth of lateral dikes and magma chamber depressurization. 

Summary Methods to compute the stress intensity factors along a three‐dimensional (3D) crack front often display a tenuous rate of convergence under mesh refinement or, worse, do not converge, particularly when applied on unstructured meshes. In this work, we propose an alternative formulation of the interaction integral functional and a method to compute stress intensity factors along the crack front which can be shown to converge. The novelty of our method is the decoupling of the two discretizations: the bulk mesh for the finite element solution and the mesh along the crack front for the numerical stress intensity factors, and hence we term it the multiple mesh interaction integral (MMII) method. Through analysis of the convergence of the functional and method, we find scalings of these two mesh sizes to guarantee convergence of the computed stress intensity factors in a variety of norms, including maximum pointwise error and total variation. We demonstrate the MMII on four examples: a semiinfinite straight crack with the asymptotic displacement fields, the same geometry with a nonuniform stress intensity factor along the crack front, a spherical cap crack in a cylinder under tension, and the elliptical crack under far‐field tension and shear. 

Summary We present an algorithm to simulate curvilinear hydraulic fractures in plane strain and axisymmetry. We restrict our attention to sharp fractures propagating in an isotropic, linear elastic medium and driven by the injection of a laminar, Newtonian fluid governed by lubrication theory, and we require the existence of a finite lag region between the fluid front and the crack tip. The key novelty of our approach is in how we discretize the evolving crack and fluid domains: we utilize universal meshes (UMs), a technique to create conforming triangulations of a problem domain by only perturbing nodes of a universal background mesh in the vicinity of the boundary. In this way, we construct meshes, which conform to the crack and to the fluid front. This allows us to build standard piecewise linear finite element spaces and to monolithically solve the quasistatic hydraulic fracture problem for the displacement field in the rock and the pressure in the fluid. We demonstrate the performance of our algorithms through three examples: a convergence study in plane strain, a comparison with experiments in axisymmetry, and a novel case of a fracture in a narrow pay zone. 

In this work, we present and prove results underlying a method which uses functionals derived from the interaction integral to approximate the stress intensity factors along a three-dimensional crack front. We first prove that the functionals possess a pair of important properties. The functionals are well-defined and continuous for square-integrable tensor fields, such as the gradient of a finite element solution. Furthermore, the stress intensity factors are representatives of such functionals in a space of functions over the crack front. Our second result is an error estimate for the numerical stress intensity factors computed via our method. The latter property of the functionals provides a recipe for numerical stress intensity factors; we apply the functionals to the gradient of a finite element approximation for a specific set of crack front variations, and we calculate the stress intensity factors by inverting the mass matrix for those variations.