Related papers: Gamra: Simple Meshes for Complex Earthquakes
Accurate and high-resolution Earth system model (ESM) simulations are essential to assess the ecological and socio-economic impacts of anthropogenic climate change, but are computationally too expensive to be run at sufficiently high…
Driven by the growing interest in numerical simulations of dislocation-interface interactions in general crystalline materials with elastic anisotropy, we develop algorithms for the integration of interface tractions needed to couple…
Noise is one of the primary sources of interference in seismic exploration. Many authors have proposed various methods to remove noise from seismic data; however, in the face of strong noise conditions, satisfactory results are often not…
Numerical simulations are essential for evaluating the performance and safety of geological engineered systems such as geologic carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately…
Immersed boundary methods have attracted substantial interest in the last decades due to their potential for computations involving complex geometries. Often these cannot be efficiently discretized using boundary-fitted finite elements.…
We present a novel isogeometric method, namely the Immersed Boundary-Conformal Method (IBCM), that features a layer of discretization conformal to the boundary while employing a simple background mesh for the remaining domain. In this…
The finite element method (FEM) is among the most commonly used numerical methods for solving engineering problems. Due to its computational cost, various ideas have been introduced to reduce computation times, such as domain decomposition,…
Within this work, we develop a phase-field description for simulating fractures in incompressible materials. Standard formulations are subject to volume-locking when the solid is (nearly) incompressible. We propose an approach that builds…
In the constraint programming framework, state-of-the-art static and dynamic decomposition techniques are hard to apply to problems with complete initial constraint graphs. For such problems, we propose a hybrid approach of these techniques…
A non-iterative waveform sensing approach is proposed toward (i) geometric reconstruction of penetrable fractures, and (ii) quantitative identification of their heterogeneous contact condition by seismic i.e. elastic waves. To this end, the…
Dynamic response evaluation in structural engineering is the process of determining the response of a structure, such as member forces, node displacements, etc when subjected to dynamic loads such as earthquakes, wind, or impact. This is an…
The mapped bases or Fake Nodes Approach (FNA), introduced in [10], allows to change the set of nodes without the need of resampling the function. Such scheme has been successfully applied in preventing the appearance of the Gibbs phenomenon…
We introduce a high-performance simulation framework that permits the semi-independent, task-based solution of sets of partial differential equations, typically manifesting as updates to a collection of `patches' in space-time. A hybrid…
Nearly all aspects of earthquake rupture are controlled by the friction along the fault that progressively increases with tectonic forcing, but in general cannot be directly measured. We show that fault friction can be determined at any…
Faults are complex systems embedded in an evolving medium fractured by seismic ruptures. This off-fault damage zone is shown to be thermo-hydro-mechano-chemically coupled to the main fault plane by a growing number of studies. Yet,…
The dissolution of solids has created spectacular geomorphologies ranging from centimeter-scale cave scallops to the kilometer-scale "stone forests" of China and Madagascar. Mathematically, dissolution processes are modeled by a Stefan…
Coupling of rupture processes in solids with waves also propagating in fluids is a prominent phenomenon arising during tectonic earthquakes. It is executed here in a single `monolithic' model which can asymptotically capture both damageable…
In this paper, we present an immersed weak Galerkin method for solving second-order elliptic interface problems. The proposed method does not require the meshes to be aligned with the interface. Consequently, uniform Cartesian meshes can be…
We develop a novel nonlocal model of dislocations based on the framework of peridynamics. By embedding interior discontinuities into the nonlocal constitutive law, the displacement jump in the Volterra dislocation model is reproduced,…
The immersed boundary method is a numerical and mathematical formulation for solving fluid-structure interaction problems. It relies on solving fluid equations on an Eulerian fluid grid and interpolating the resulting velocity back onto…