Surface waves Anisotropy Modeling

a b s t r a c t We 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 symme-tries, 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 mate-rial 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 deriv-atives. We test the algorithms by comparison to the analytical solution and obtain the wave field at dif-ferent 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 algo-rithms can handle large impedance variations at the interface. 2011 Elsevier B.V. All rights reserved. 1