Related papers: Solving seismic wave propagation in elastic media …
We propose a computationally efficient and systematically convergent approach for elastodynamics simulations. We recast the second-order dynamical equation of elastodynamics into an equivalent first-order system of coupled equations, so as…
Exponential time integrators have been applied successfully in several physics-related differential equations. However, their application in hyperbolic systems with absorbing boundaries, like the ones arising in seismic imaging, still lacks…
We present a one-step algorithm that solves the Maxwell equations for systems with spatially varying permittivity and permeability by the Chebyshev method. We demonstrate that this algorithm may be orders of magnitude more efficient than…
We review some recent developments in numerical algorithms to solve the time-dependent Maxwell equations for systems with spatially varying permittivity and permeability. We show that the Suzuki product-formula approach can be used to…
In geophysics, wave propagation in elastic media is a crucial subject. In this context, seismology has made significant progress as a result of numerous advances, among these stands out the advancement of numerical methods such as the…
We propose a computationally efficient technique for extrapolating seismic waves in an arbitrary isotropic elastic medium. The method is based on factorizing the full elastic wave equation into a product of pseudo-differential operators.…
This paper discusses new simulation algorithms for stochastic chemical kinetics that exploit the linearity of the chemical master equation and its matrix exponential exact solution. These algorithms make use of various approximations of the…
The development of finite/fixed-time stable optimization algorithms typically involves study of specific problem instances. The lack of a unified framework hinders understanding of more sophisticated algorithms, e.g., primal-dual gradient…
Elastic wave propagation is studied in a heterogeneous 2-D medium consisting of an elastic matrix containing randomly distributed circular elastic inclusions. The aim of this study is to determine the effective wavenumbers when the incident…
We analyse and compare several algorithms to compute numerically periodic solutions of high-dimensional dynamical systems and investigate their Floquet stability without building the monodromy matrix. The solution and its perturbation are…
An FFT-based algorithm is developed to simulate the propagation of elastic waves in heterogeneous $d$-dimensional rectangular shape domains. The method allows one to prescribe the displacement as a function of time in a subregion of the…
A framework for exponential time discretization of the multilayer rotating shallow water equations is developed in combination with a mimetic discretization in space. The method is based on a combination of existing exponential time…
Simulating the propagation of elastic waves in multi-layered media has many applications. A common approach is to use matrix methods where the elastic wave-field within each material layer is represented by a sum of partial-waves along with…
In this paper we study quantum simulation algorithms on the elastic wave equations using the Schr\"odingerisation method. The Schr\"odingerisation method transforms any linear PDEs into a system of Schr\"odinger-type PDEs -with unitary…
Approximated numerical techniques, for the solution of the elastic wave scattering problem over semi-infinite domains are reviewed. The approximations involve the representation of the half-space by a boundary condition described in terms…
Matrix evolution equations occur in many applications, such as dynamical Lyapunov/Sylvester systems or Riccati equations in optimization and stochastic control, machine learning or data assimilation. In many such problems, the dominant…
We study step-wise time approximations of non-linear hyperbolic initial value problems. The technique used here is a generalization of the minimizing movements method, using two time-scales: one for velocity, the other (potentially much…
We present a practical algorithm to approximate the exponential of skew-Hermitian matrices up to round-off error based on an efficient computation of Chebyshev polynomials of matrices and the corresponding error analysis. It is based on…
Consider the elastic scattering of a time-harmonic wave by multiple well separated rigid particles in two dimensions. To avoid using the complex Green's tensor of the elastic wave equation, we utilize the Helmholtz decomposition to convert…
In contrast to the traditional elastodynamic equations, a more comprehensive formulation of one dimensional (1D) elastodynamic equations is given for inhomogeneous media by using the coordinate transformation method. These modified…