Related papers: Gamra: Simple Meshes for Complex Earthquakes
Iterative geostatistical seismic inversion integrates seismic and well data to infer the spatial distribution of subsurface elastic properties. These methods provide limited assessment to the spatial uncertainty of the inverted elastic…
Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport…
Fluid discontinuities, such as shock fronts and vortex sheets, can reflect waves and become unstable to corrugation. Analytical calculations of these phenomena are tractable in the simplest cases only, while their numerical simulations are…
A marker-and-cell finite difference method is developed for solving the two dimensional and three dimensional linear elasticity in the displacement-stress formulation on staggered grids. The method employs a staggered grid arrangement,…
We present a multi-resolution approach for constructing model-based simulations of hydraulic fracturing, wherein flow through porous media is coupled with fluid-driven fracture. The approach consists of a hybrid scheme that couples a…
Computational earthquake sequence models provide generative estimates of the time, location, and size of synthetic seismic events that can be compared with observed earthquake histories and assessed as rupture forecasts. Here we describe a…
The ground structure can substantially influence seismic ground motion underscoring the need to develop a ground structure model with sufficient reliability in terms of ground motion estimation for earthquake damage mitigation. While many…
The immersed boundary method is a mathematical formulation and numerical method for solving fluid-structure interaction problems. For many biological problems, such as models that include the cell membrane, the immersed structure is a…
A new immersed finite element (IFE) method is developed for second-order elliptic problems with discontinuous diffusion coefficient. The IFE space is constructed based on the rotated Q1 nonconforming finite elements with the integral-value…
Dislocations are the main carriers of the permanent deformation of crystals. For simulations of engineering applications, continuum models where material microstructures are represented by continuous density distributions of dislocations…
We present a variational quantum algorithm for structural mechanical problems, specifically addressing crack opening simulations that traditionally require extensive computational resources. Our approach provides an alternative solution for…
Efficient simulation of plasmas in various contexts often involves the use of meshes that conform to the intrinsic geometry of the system under consideration. We present here a description of a new magnetohydrodynamic code, Gamera (Grid…
We present a unified framework to tie overlapping meshes in solid mechanics applications. This framework is a combination of the X-FEM method and the mortar method, which uses Lagrange multipliers to fulfill the tying constraints. As known,…
Within the framework of the displacement-based Virtual Element Method (VEM) for plane elasticity a significant problem is represented by an accurate evaluation of the stress field. In particular, in the classical VEM formulation, a suitable…
This paper presents two approaches to mathematical modelling of a synthetic seismic pulse, and a comparison between them. First, a new analytical model is developed in two-dimensional Cartesian coordinates. Combined with an initial…
We develop a method for solving elliptic partial differential equations on surfaces described by CAD patches that may have gaps/overlaps. The method is based on hybridization using a three-dimensional mesh that covers the gap/overlap…
Simulations of fluid flow in naturally fractured rocks have implications for several subsurface applications, including energy storage and extraction, and waste storage. We are interested in flow in discrete fracture networks, which…
We consider numerical simulation of the isotropic elastic wave equations arising from seismic applications with non-trivial land topography. The more flexible finite element method is applied to the shallow region of the simulation domain…
We present a finite element discretisation to model the interaction between a poroelastic structure and an elastic medium. The consolidation problem considers fully coupled deformations across an interface, ensuring continuity of…
We develop a new finite element method for solving planar elasticity problems involving of heterogeneous materials with a mesh not necessarily aligning with the interface of the materials. This method is based on the `broken'…