Elastic surface waves in crystals – Part 2: Cross-check of two full-wave numerical modeling methods

Abstract

Nathalie Favretto-Cristini. Elastic surface waves in crystals -part 2: cross-check of two full-wave numerical modeling methods. Ultrasonics, Elsevier, 2011, 51 (8) AbstractWe obtain the full-wave solution for the wave propagation at the surface of anisotropic media using two spectral numerical modeling algorithms. The simulations focus on media of cubic and hexagonal symmetries, for which the physics has been reviewed and clarified in a companion paper. Even in the case of homogeneous media, the solution requires the use of numerical methods because the analytical Green's function cannot be obtained in the whole space. The algorithms proposed here allow for a general material variability and the description of arbitrary crystal symmetry at each grid point of the numerical mesh. They are based on high-order spectral approximations of the wave field for computing the spatial derivatives. We test the algorithms by comparison to the analytical solution and obtain the wave field at different faces (stress-free surfaces) of apatite, zinc and copper. Finally, we perform simulations in heterogeneous media, where no analytical solution exists in general, showing that the modeling algorithms can handle large impedance variations at the interface.

The perfectly matched layer (PML) for hyperbolic wave propagation problems: A review

“…A typical example where the geometric stability condition is violated is in advective acoustics or linearised Euler equations with non-vanishing mean flow [70,105,84,1,71,63,33]. Other examples often include certain crystals and anisotropic elastic media [75,35,41,15]. See for example the anisotropic material AM2 displayed in Figure 2 d).…”

“…The following theorem and corollary were proven in [39] Theorem 4 Consider the PML equation in the Laplace space (70) with constant damping d ξ ≥ 0, and source terms F(x, y, z, s), f , subject to homogeneous initial data, the boundary conditions (71), and the jump condition (72) at discontinuities in material parameters. Let the energy norms E p (s, d ξ ) > 0, E f (s, d ξ ) > 0 in the Laplace space be defined in (75). For any Re{s} = a > 0 we have…”

A stable discontinuous Galerkin method for linear elastodynamics in 3D geometrically complex media using physics based numerical fluxes

“…This can have a huge impact in high performance computing applications, since we can easily design efficient communication avoiding parallel algorithms. We will present numerical simulations verifying accuracy and stability of the method, using community developed benchmark problems [16,17,1,2]. We will also present simulation on real geologically constrained complex geometries.…”

Assessment of Risks Induced by Countermining Unexploded Large-Charge Historical Ordnance in a Shallow Water Environment—Part II: Modeling of Seismo-Acoustic Wave Propagation

“…Very efficient perfectly matched layers (PML) can be used to limit the size of the studied domain (e.g., [17]), and thus, to reduce the required computational resources that may otherwise become prohibitive for large size of 3-D domains and high-frequency (HF) simulations. In addition, the effect of wave attenuation can be accurately taken into account, and it has been shown that the behavior of surface and interface waves is accurately modeled (e.g., [18]), which is particularly important in this study. More specifically, we use the SPECFEM software package (https://geodynamics.org/cig/software/ specfem2d/).…”