Source author record

Adrian J. Lew

Adrian J. Lew appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Computational Geometry math.MG physics.geo-ph

Catalog footprint

What is connected

7works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2020arXiv

Spring-Slider and Finite Element Modeling of Microseismic Events and Fault Slip during Hydraulic Fracturing

Hydraulic fracturing increases reservoir permeability by opening fractures and triggering slip on natural fractures and faults. While seismic slip of small faults or fault patches is detectable as microseismic events, the role of aseismic slip is poorly understood. From a modeling standpoint, geomechanical analysis using the Coulomb criterion can determine if faults slip but not whether slip is seismic or aseismic. Here we propose a computational methodology to predict fault slip, and whether slip is seismic or aseismic, using rate-and-state friction. To avoid computational costs associated with resolving small faults, we use the spring-slider idealization that treats faults as points. Interaction between faults is neglected. The method is applied to study fault slip from a hydraulic fracture that grows past a fault, without intersecting it. We represent the hydraulic fracture stressing using an asymptotic expansion of stresses around the tip of a tensile crack. We investigate the effect of fault length, orientation, and distance from the hydraulic fracture. For velocity-weakening faults with stiffness smaller than a critical stiffness, slip is seismic, whereas faults with stiffness greater than the critical stiffness slip aseismically. Furthermore, we compare the spring-slider idealization with a finite element analysis that resolves spatially variable slip. The spring-slider idealization provides reasonably accurate predictions of moment and even moment-rate history, especially for faults having stiffness close to or larger than the critical stiffness. Differences appear for large faults where rupture propagation is important, though differences might still be negligible for many applications.

preprint2015arXiv

Computing stress intensity factors for curvilinear cracks

The use of the interaction integral to compute stress intensity factors around a crack tip requires selecting an auxiliary field and a material variation field. We formulate a family of these fields accounting for the curvilinear nature of cracks that, in conjunction with a discrete formulation of the interaction integral, yield optimally convergent stress intensity factors. We formulate three pairs of auxiliary and material variation fields chosen to yield a simple expression of the interaction integral for different classes of problems. The formulation accounts for crack face tractions and body forces. Distinct features of the fields are their ease of construction and implementation. The resulting stress intensity factors are observed converging at a rate that doubles the one of the stress field. We provide a sketch of the theoretical justification for the observed convergence rates, and discuss issues such as quadratures and domain approximations needed to attain such convergent behavior. Through two representative examples, a circular arc crack and a loaded power function crack, we illustrate the convergence rates of the computed stress intensity factors. The numerical results also show the independence of the method on the size of the domain of integration.

preprint2015arXiv

Universal Meshes for the Simulation of Brittle Fracture and Moving Boundary Problems

Universal meshes have recently appeared in the literature as a compu- tationally efficient and robust paradigm for the generation of conforming simpli- cial meshes for domains with evolving boundaries. The main idea behind a univer- sal mesh is to immerse the moving boundary in a background mesh (the universal mesh), and to produce a mesh that conforms to the moving boundary at any given time by adjusting a few of elements of the background mesh. In this manuscript we present the application of universal meshes to the simulation of brittle fracturing. To this extent, we provide a high level description of a crack propagation algorithm and showcase its capabilities. Alongside universal meshes for the simulation of brit- tle fracture, we provide other examples for which universal meshes prove to be a powerful tool, namely fluid flow past moving obstacles. Lastly, we conclude the manuscript with some remarks on the current state of universal meshes and future directions.

preprint2014arXiv

High-Order Finite Element Methods for Moving Boundary Problems with Prescribed Boundary Evolution

We introduce a framework for the design of finite element methods for two-dimensional moving boundary problems with prescribed boundary evolution that have arbitrarily high order of accuracy, both in space and in time. At the core of our approach is the use of a universal mesh: a stationary background mesh containing the domain of interest for all times that adapts to the geometry of the immersed domain by adjusting a small number of mesh elements in the neighborhood of the moving boundary. The resulting method maintains an exact representation of the (prescribed) moving boundary at the discrete level, or an approximation of the appropriate order, yet is immune to large distortions of the mesh under large deformations of the domain. The framework is general, making it possible to achieve any desired order of accuracy in space and time by selecting a preferred and suitable finite-element space on the universal mesh for the problem at hand, and a preferred and suitable time integrator for ordinary differential equations. We illustrate our approach by constructing a particular class of methods, and apply them to a prescribed-boundary variant of the Stefan problem. We present numerical evidence for the order of accuracy of our schemes in one and two dimensions.

preprint2014arXiv

Supercloseness of Orthogonal Projections onto Nearby Finite Element Spaces

We derive upper bounds on the difference between the orthogonal projections of a smooth function $u$ onto two finite element spaces that are nearby, in the sense that the support of every shape function belonging to one but not both of the spaces is contained in a common region whose measure tends to zero under mesh refinement. The bounds apply, in particular, to the setting in which the two finite element spaces consist of continuous functions that are elementwise polynomials over shape-regular, quasi-uniform meshes that coincide except on a region of measure $O(h^γ)$, where $γ$ is a nonnegative scalar and $h$ is the mesh spacing. The projector may be, for example, the orthogonal projector with respect to the $L^2$- or $H^1$-inner product. In these and other circumstances, the bounds are superconvergent under a few mild regularity assumptions. That is, under mesh refinement, the two projections differ in norm by an amount that decays to zero at a faster rate than the amounts by which each projection differs from $u$. We present numerical examples to illustrate these superconvergent estimates and verify the necessity of the regularity assumptions on $u$.

preprint2013arXiv

Analysis of a method to parameterize planar curves immersed in triangulations

We prove that a planar $C^2$-regular boundary $Γ$ can always be parameterized with its closest point projection $π$ over a certain collection of edges $Γ_h$ in an ambient triangulation, by making simple assumptions on the background mesh. For $Γ_h$, we select the edges that have both vertices on one side of $Γ$ and belong to a triangle that has a vertex on the other side. By imposing restrictions on the size of triangles near the curve and by requesting that certain angles in the mesh be strictly acute, we prove that $π:Γ_h\rightarrowΓ$ is a homeomorphism, that it is $C^1$ on each edge in $Γ_h$ and provide bounds for the Jacobian of the parameterization. The assumptions on the background mesh are both easy to satisfy in practice and conveniently verified in computer implementations. The parameterization analyzed here was previously proposed by the authors and applied to the construction of high-order curved finite elements on a class of planar piecewise $C^2$-curves.

preprint2012arXiv

Universal Meshes: A new paradigm for computing with nonconforming triangulations

We describe a method for discretizing planar C2-regular domains immersed in non-conforming triangulations. The method consists in constructing mappings from triangles in a background mesh to curvilinear ones that conform exactly to the immersed domain. Constructing such a map relies on a novel way of parameterizing the immersed boundary over a collection of nearby edges with its closest point projection. By interpolating the mappings to curvilinear triangles at select points, we recover isoparametric mappings for the immersed domain defined over the background mesh. Indeed, interpolating the constructed mappings just at the vertices of the background mesh yields a fast meshing algorithm that involves only perturbing a few vertices near the boundary. For the discretization of a curved domain to be robust, we have to impose restrictions on the background mesh. Conversely, these restrictions define a family of domains that can be discretized with a given background mesh. We then say that the background mesh is a universal mesh for such a family of domains. The notion of universal meshes is particularly useful in free/moving boundary problems because the same background mesh can serve as the universal mesh for the evolving domain for time intervals that are independent of the time step. Hence it facilitates a framework for finite element calculations over evolving domains while using a fixed background mesh. Furthermore, since the evolving geometry can be approximated with any desired order, numerical solutions can be computed with high-order accuracy. We demonstrate these ideas with various numerical examples.