Related papers: Gamra: Simple Meshes for Complex Earthquakes
We present an algorithm to solve the dispersive depth-averaged Serre-Green-Naghdi (SGN) equations using patch-based adaptive mesh refinement. These equations require adding additional higher derivative terms to the nonlinear shallow water…
Recently, the phase field method has been increasingly used for brittle fractures in soft materials like polymers, elastomers, and biological tissues. When considering finite deformations to account for the highly deformable nature of soft…
We propose a model order reduction approach to speed up the computation of seismograms, i.e. the solution of the seismic wave equation evaluated at a receiver location, for different model parameters. Our approach achieves a reduction of…
We present a computational framework for solving the equations of inviscid gas dynamics using structured grids with embedded geometries. The novelty of the proposed approach is the use of high-order discontinuous Galerkin (dG) schemes and a…
The modeling of cracks has been an intensely researched topic for decades - both from the mechanical as well as from the mathematics point of view. As far as the modeling of sharp cracks/interfaces is concerned, the resulting free boundary…
In the fast developing countries it is hard to trace new buildings construction or old structures destruction and, as a result, to keep the up-to-date cadastre maps. Moreover, due to the complexity of urban regions or inconsistency of data…
We consider variational problems that model the bending behavior of curves that are constrained to belong to given hypersurfaces. Finite element discretizations of corresponding functionals are justified rigorously via Gamma-convergence.…
SUMMARY Geophysical imaging using the inversion procedure is a powerful tool for the exploration of the Earth's subsurface. However, the interpretation of inverted images can sometimes be difficult, due to the inherent limitations of…
Travel time tomography is used to infer the underlying three-dimensional wavespeed structure of the Earth by fitting seismic travel time data collected at surface stations. Data interpolation and denoising techniques are important…
The $\delta f$ particle-in-cell algorithm has been a useful tool in studying the physics of plasmas, particularly turbulent magnetized plasmas in the context of gyrokinetics. The reduction in noise due to not having to resolve the full…
Our understanding of earthquakes is based on the theory of plate tectonics. Earthquake dynamics is the study of the interactions of plates (solid disjoint parts of the lithosphere) which produce seismic activity. Over the last about fifty…
We design the conforming virtual element method for the numerical approximation of the two dimensional elastodynamics problem. We prove stability and convergence of the semi-discrete approximation and derive optimal error estimates under…
Numerical and analytical methods are developed for the investigation of contact sets in electrostatic-elastic deflections modeling micro-electro mechanical systems. The model for the membrane deflection is a fourth-order semi-linear partial…
Accurate interpolation of seismic data is crucial for improving the quality of imaging and interpretation. In recent years, deep learning models such as U-Net and generative adversarial networks have been widely applied to seismic data…
We present a collision model for phase-resolved Direct Numerical Simulations of sediment transport that couple the fluid and particles by the Immersed Boundary Method. Typically, a contact model for these types of simulations comprises a…
The study of flow in fractured porous media is a key ingredient for many geoscience applications, such as reservoir management and geothermal energy production. Modelling and simulation of these highly heterogeneous and geometrically…
Reliable earthquake detection and seismic phase classification is often challenging especially in the circumstances of low magnitude events or poor signal-to-noise ratio. With improved seismometers and better global coverage, a sharp…
The Gradient Scheme framework provides a unified analysis setting for many different families of numerical methods for diffusion equations. We show in this paper that the Gradient Scheme framework can be adapted to elasticity equations, and…
Integrated Nested Laplace Approximation provides a fast and effective method for marginal inference on Bayesian hierarchical models. This methodology has been implemented in the R-INLA package which permits INLA to be used from within R…
We present results from a comparative analysis of empirical and synthetic shapes for isoseismals of low intensity (I = IV-VI on the MCS scale) for six Italian earthquakes of ML = 4.5-6. Our modelling of isoseismals is based on a plane…