Relations between Transfer Matrices and Numerical Stability Analysis to Avoid the Ømega d Problem

The transfer matrix method is usually employed to study problems described by $N$ equations of matrix Sturm--Liouville (MSL) kind. In some cases a numerical degradation (the so called $\Omega d$ problem) appears thus impairing the performance of the method. We present here a procedure that can overcome this problem in the case of multilayer systems having piecewise constant coefficients. This is performed by studying the relations between the associated transfer matrix $\boldsymbol{(T)}$ and other transfer matrix variants. In this way it was possible to obtain the matrices which can overcome the $\Omega d$ problem in the general case and then in problems which are particular cases of the general one. In this framework different strategies are put forward to solve different boundary condition problems by means of these numerically stable matrices. Numerical and analytic examples are presented to show that these stable variants are more adequate than other matrix methods to overcome the $\Omega d$ problem. ...

[1]  Eng Leong Tan,et al.  Matrix Algorithms for Modeling Acoustic Waves in Piezoelectric Multilayers , 2007, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[2]  G. Bonnet Orthotropic elastic media having a closed form expression of the Green tensor , 2009 .

[3]  O. Poncelet,et al.  Love waves in two-dimensional phononic crystals with depth-dependent properties , 2013 .

[4]  Karl Meerbergen,et al.  The Quadratic Eigenvalue Problem , 2001, SIAM Rev..

[5]  A. Shuvalov A sextic formalism for three–dimensional elastodynamics of cylindrically anisotropic radially inhomogeneous materials , 2003, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[6]  J. W. Dunkin,et al.  Computation of modal solutions in layered, elastic media at high frequencies , 1965 .

[7]  A. Shcherbakov,et al.  Rigorous S-matrix approach to the modeling of the optical properties of OLEDs , 2011 .

[8]  Federico García Moliner,et al.  Transfer matrix, green function and related techniques: tools for the study of multilayer heterostructures , 2004 .

[9]  Eng Leong Tan,et al.  Stiffness matrix method with improved efficiency for elastic wave propagation in layered anisotropic media , 2005 .

[10]  F. Szmulowicz Kronig - Penney model: a new solution , 1997 .

[11]  Joel N. Schulman,et al.  Wave Mechanics Applied to Semiconductor Heterostructures , 1991 .

[12]  R. Pérez-Álvarez,et al.  Long wave polar modes in semiconductor heterostructures , 1998 .

[13]  W. Penney,et al.  Quantum Mechanics of Electrons in Crystal Lattices , 1931 .

[14]  Resolvent Method for Calculating Dispersion Spectra of the Shear Waves in the Phononic Plates and Waveguides , 2014 .

[15]  F. Szmulowicz New Kronig–Penney equation emphasizing the band edge conditions , 2008 .

[16]  Feruza A. Amirkulova,et al.  Stable methods to solve the impedance matrix for radially inhomogeneous cylindrically anisotropic structures , 2012, 1212.5840.

[17]  S. Rokhlin,et al.  Stable recursive algorithm for elastic wave propagation in layered anisotropic media: stiffness matrix method. , 2002, The Journal of the Acoustical Society of America.

[18]  E. L. Tan Generalized eigenproblem of hybrid matrix for Floquet wave propagation in one-dimensional phononic crystals with solids and fluids. , 2010, Ultrasonics.

[19]  Eng Leong Tan,et al.  A concise and efficient scattering matrix formalism for stable analysis of elastic wave propagation in multilayered anisotropic solids. , 2003, Ultrasonics.

[20]  Weiqiu Chen,et al.  Guided wave propagation in multilayered piezoelectric structures , 2009 .

[21]  Yu Jiangong,et al.  On dispersion relations of waves in multilayered magneto-electro-elastic plates , 2012 .

[22]  A. Norris,et al.  Wave impedance matrices for cylindrically anisotropic radially inhomogeneous elastic solids , 2010, 1003.5713.

[23]  Witold Hurewicz,et al.  Lectures on Ordinary Differential Equations , 1959 .

[24]  Eng Leong Tan,et al.  Hybrid compliance-stiffness matrix method for stable analysis of elastic wave propagation in multilayered anisotropic media. , 2006, The Journal of the Acoustical Society of America.

[25]  L. Wang,et al.  A compliance/stiffness matrix formulation of general Green's function and effective permittivity for piezoelectric multilayers , 2004, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[26]  M.J.S. Lowe,et al.  Matrix techniques for modeling ultrasonic waves in multilayered media , 1995, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[27]  Min-zhong Wang,et al.  Three-dimensional Green’s functions for infinite anisotropic piezoelectric media , 2007 .

[28]  L. Leija,et al.  Dispersion curves of shear horizontal wave surface velocities in multilayer piezoelectric systems , 2010 .

[29]  S. Rokhlin,et al.  Ultrasonic wave interaction with a thin anisotropic layer between two anisotropic solids: Exact and asymptotic‐boundary‐condition methods , 1992 .

[30]  Ernian Pan,et al.  Wave propagation in magneto-electro-elastic multilayered plates , 2007 .

[31]  Nicholas J. Higham,et al.  INVERSE PROBLEMS NEWSLETTER , 1991 .

[32]  Mei Han An,et al.  accuracy and stability of numerical algorithms , 1991 .