Related papers: Integral equation method for microseismic wavefiel…
Three numerical algorithms are proposed to solve the time-dependent elastodynamic equations in elastic solids. All algorithms are based on approximating the solution of the equations, which can be written as a matrix exponential. By…
The geometry of mesoscopic inhomogeneities plays an important role in determining the macroscopic propagation behaviors of elastic waves in a heterogeneous medium. Nonequiaxed inhomogeneities can lead to anisotropic wave velocity and…
We provide an analytical formulation to model the propagation of elastic waves in a homogeneous half-space supporting an array of thin plates. The technique provides the displacement field obtained from the interaction between an incident…
This paper deals with the numerical modeling of transient mechanical waves in linear viscoelastic solids. Dissipation mechanisms are described using the generalized Zener model. No time convolutions are required thanks to the introduction…
Starting from the hydrostatic Boussinesq equations, we derive a time-averaged `hydrostatic wave equation' that describes the propagation of inertia-gravity internal waves through quasi-geostrophic flow. The derivation uses a…
This paper proposes a frequency-time hybrid solver for the time-dependent wave equation in two-dimensional interior spatial domains. The approach relies on four main elements, namely, 1) A multiple scattering strategy that decomposes a…
We present a novel immersed boundary method that implements acoustic perturbation theory to model an acoustically levitated droplet. Instead of resolving sound waves numerically, our hybrid method solves acoustic scattering…
An approximate solution of the Dyson equation related to a stochastic Helmholtz equation, which describes the acoustic dynamics of a three-dimensional isotropic random medium with elastic tensor fluctuating in space, is obtained in the…
Full-Waveform Inversion (FWI) is a nonlinear iterative seismic imaging technique that, by reducing the misfit between recorded and predicted seismic waveforms, can produce detailed estimates of subsurface geophysical properties.…
We present an adjoint-based optimization method to invert for stress and frictional parameters used in earthquake modeling. The forward problem is linear elastodynamics with nonlinear rate-and-state frictional faults. The misfit functional…
This paper deals with the spectral element modeling of seismic wave propagation at the global scale. Two aspects relevant to low-frequency studies are particularly emphasized. First, the method is generalized beyond the Cowling…
We present a stable finite element method for incompressible nonlinear elasticity based on a four-field mixed formulation involving the displacement, displacement gradient, first Piola--Kirchhoff stress and pressure. Unlike existing…
This paper presents a novel approach for implementing frequency-dependent hydrodynamic coefficients in Morison's equation, which is widely used in hydrodynamics modeling. Accurate hydrodynamic predictions using Morison's equation…
In this paper we develop a multiple scattering model for elastic waves in random anisotropic media. It relies on a kinetic approach of wave propagation phenomena pertaining to the situation whereby the wavelength is comparable to the…
We present a novel method to investigate the dynamics of a single semiflexible polymer, subject to anisotropic friction in a viscous fluid. In contrast to previous approaches, we do not rely on a discrete bead-rod model, but introduce a…
Derivatives of the displacement tensor with respect to the independent model parameters of the subsurface, also called Fr\'echet derivatives (or sensitivity kernels), are a key ingredient for seismic full-waveform inversion with a…
We extend the theory of complete Bernstein functions to matrix-valued functions and apply it to analyze Green's function of an anisotropic multi-dimension\-al linear viscoelastic problem. Green's function is given by the superposition of…
The spectral decomposition of the elastic wave operator in a layered isotropic half-space is derived by means of standard functional analytic methods. Particular attention is paid to the coupled $P$-$SV$ waves. The problem is formulated…
We developed a 3-D elastic Boundary Element Method (BEM) computer code, called AstroSeis, to model seismic wavefields in a body with an arbitrary shape, such as an asteroid. Besides the AstroSeis can handle arbitrary surface topography, it…
Boundary integral methods are attractive for solving homogeneous linear constant coefficient elliptic partial differential equations on complex geometries, since they can offer accurate solutions with a computational cost that is linear or…