skip to main content


Title: A Framework for Simulation of Multiple Elastic Scattering in Two Dimensions
Consider the elastic scattering of a time-harmonic wave by multiple well-separated rigid particles with smooth boundaries in two dimensions. Instead of using the complex Green's tensor of the elastic wave equation, we utilize the Helmholtz decomposition to convert the boundary value problem of the elastic wave equation into a coupled boundary value problem of the Helmholtz equation. Based on single, double, and combined layer potentials with the simpler Green's function of the Helmholtz equation, we present three different boundary integral equations for the coupled boundary value problem. The well-posedness of the new integral equations is established. Computationally, a scattering matrix based method is proposed to evaluate the elastic wave for arbitrarily shaped particles. The method uses the local expansion for the incident wave and the multipole expansion for the scattered wave. The linear system of algebraic equations is solved by GMRES with fast multipole method (FMM) acceleration. Numerical results show that the method is fast and highly accurate for solving elastic scattering problems with multiple particles.  more » « less
Award ID(s):
1912704
NSF-PAR ID:
10182372
Author(s) / Creator(s):
;
Date Published:
Journal Name:
SIAM journal on scientific computing
Volume:
41
ISSN:
1064-8275
Page Range / eLocation ID:
A3276-A3299
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Consider the scattering of a time-harmonic acoustic plane wave by a bounded elastic obstacle which is immersed in a homogeneous acoustic medium. This paper is concerned with an inverse acoustic-elastic interaction problem, which is to determine the location and shape of the elastic obstacle by using either the phased or phaseless far-field data. By introducing the Helmholtz decomposition, the model problem is reduced to a coupled boundary value problem of the Helmholtz equations. The jump relations are studied for the second derivatives of the single-layer potential in order to deduce the corresponding boundary integral equations. The well-posedness is discussed for the solution of the coupled boundary integral equations. An efficient and high order Nyström-type discretization method is proposed for the integral system. A numerical method of nonlinear integral equations is developed for the inverse problem. For the case of phaseless data, we show that the modulus of the far-field pattern is invariant under a translation of the obstacle. To break the translation invariance, an elastic reference ball technique is introduced. We prove that the inverse problem with phaseless far-field pattern has a unique solution under certain conditions. In addition, a numerical method of the reference ball technique based nonlinear integral equations is proposed for the phaseless inverse problem. Numerical experiments are presented to demonstrate the effectiveness and robustness of the proposed methods. 
    more » « less
  2. We present two (a decoupled and a coupled) integral-equation-based methods for the Morse-Ingard equations subject to Neumann boundary conditions on the exterior domain. Both methods are based on second-kind integral equation (SKIE) formulations. The coupled method is well-conditioned and can achieve high accuracy. The decoupled method has lower computational cost and more flexibility in dealing with the boundary layer; however, it is prone to the ill-conditioning of the decoupling transform and cannot achieve as high accuracy as the coupled method. We show numerical examples using a Nyström method based on quadrature-by-expansion (QBX) with fast-multipole acceleration. We demonstrate the accuracy and efficiency of the solvers in both two and three dimensions with complex geometry. 
    more » « less
  3. Consider the scattering of a time-harmonic elastic plane wave by a periodic rigid surface. The elastic wave propagation is governed by the two-dimensional Navier equation. Based on a Dirichlet-to-Neumann (DtN) map, a transparent boundary condition (TBC) is introduced to reduce the scattering problem into a boundary value problem in a bounded domain. By using the finite element method, the discrete problem is considered, where the TBC is replaced by the truncated DtN map. A new duality argument is developed to derive the a posteriori error estimate, which contains both the finite element approximation error and the DtN truncation error. An a posteriori error estimate based adaptive finite element algorithm is developed to solve the elastic surface scattering problem. Numerical experiments are presented to demonstrate the effectiveness of the proposed method. 
    more » « less
  4. There is scientific and industrial interest in understanding how geologic faults respond to transient sources of fluid. Natural and artificial sources can elevate pore fluid pressure on the fault frictional interface, which may induce slip. We consider a simple boundary value problem to provide an elementary model of the physical process and to provide a benchmark for numerical solution procedures. We examine the slip of a fault that is an interface of two elastic half-spaces. Injection is modelled as a line source at constant pressure and fluid pressure is assumed to diffuse along the interface. The resulting problem is an integro-differential equation governing fault slip, which has a single dimensionless parameter. The expansion of slip is self-similar and the rupture front propagates at a factor $\lambda$ of the diffusive length scale $\sqrt {\alpha t}$ . We identify two asymptotic regimes corresponding to $\lambda$ being small or large and perform a perturbation expansion in each limit. For large $\lambda$ , in the regime of a so-called critically stressed fault, a boundary layer emerges on the diffusive length scale, which lags far behind the rupture front. We demonstrate higher-order matched asymptotics for the integro-differential equation, and in doing so, we derive a multipole expansion to capture successive orders of influence on the outer problem for fault slip for a driving force that is small relative to the crack dimensions. Asymptotic expansions are compared with accurate numerical solutions to the full problem, which are tabulated to high precision. 
    more » « less
  5. In some applications, it is reasonable to assume that geodesics (rays) have a consistent orientation so that a time-harmonic elastic wave equation may be viewed as an evolution equation in one of the spatial directions. With such applications in mind, motivated by our recent work [Hadamard- Babich ansatz for point-source elastic wave equations in variable media at high frequencies, Multiscale Model Simul. 19/1 (2021) 46–86], we propose a new truncated Hadamard-Babich ansatz based globally valid asymptotic method, dubbed the fast Huygens sweeping method, for computing Green’s functions of frequency-domain point-source elastic wave equations in inhomogeneous media in the high-frequency asymptotic regime and in the presence of caustics. The first novelty of the fast Huygens sweeping method is that the Huygens-Kirchhoff secondary-source principle is used to integrate many locally valid asymptotic solutions to yield a globally valid asymptotic solution so that caustics can be treated automatically. This yields uniformly accurate solutions both near the source and away from it. The second novelty is that a butterfly algorithm is adapted to accelerate matrix-vector products induced by the Huygens-Kirchhoff integral. The new method enjoys the following desired features: (1) it treats caustics automatically; (2) precomputed asymptotic ingredients can be used to construct Green’s functions of elastic wave equations for many different point sources and for arbitrary frequencies; (3) given a specified number of points per wavelength, it constructs Green’s functions in nearly optimal complexity O(N logN) in terms of the total number of mesh points N, where the prefactor of the complexity depends only on the specified accuracy and is independent of the frequency parameter. Three-dimensional numerical examples are presented to demonstrate the performance and accuracy of the new method. 
    more » « less