Modal method for the 2D wave propagation in heterogeneous anisotropic media

A multimodal method based on a generalization of the admittance matrix is used to analyze wave propagation in heterogeneous two-dimensional anisotropic media. The heterogeneity of the medium can be due to the presence of anisotropic inclusions with arbitrary shapes, to a succession of anisotropic media with complex interfaces between them, or both. Using a modal expansion of the wave field, the problem is reduced to a system of two sets of first-order differential equations for the modal components of the field, similar to the system obtained in the rigorous coupled wave analysis. The system is solved numerically, using the admittance matrix, which leads to a stable numerical method, the basic properties of which are discussed. The convergence of the method is discussed, considering arrays of anisotropic inclusions with complex shapes, which tend to show that Li's rules are not concerned within our approach. The method is validated by comparison with a subwavelength layered structure presenting an effective anisotropy at the wave scale.