Relations between transfer matrices and numerical stability analysis to avoid the d problem

DOI10.1137/140993442MaRDI QIDQ5261727zbMATH OpenOpenAlexFDO

Authors R. Pérez-Álvarez, René Pernas-Salomón, V. R. Velasco

Publication date 7 July 2015

Published in SIAM Journal on Applied Mathematics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1503.09038

zbMATH Keywords

\(\Omega d\) problemmatrix Sturm-Liouville problemnumerical stabilityquadratic eigenvaluestransfer matrix

Mathematics Subject Classification ID

Implicit ordinary differential equations, differential-algebraic equations (34A09) Linear ordinary differential equations and systems (34A30) Sturm-Liouville theory (34B24) Qualitative investigation and simulation of ordinary differential equation models (34C60) Numerical approximation of eigenvalues and of other parts of the spectrum of ordinary differential operators (34L16) Waves and radiation in optics and electromagnetic theory (78A40)

Abstract: 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 $Omegad$ 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 and other transfer matrix variants. In this way it was possible to obtain the matrices which can overcome the $Omegad$ 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 $Omegad$ problem. Due to the ubiquity of the MSL system, these results can be applied to the study of many elementary excitations in multilayer structures.

Recommendations

• scientific article; zbMATH DE number 2156328
• Transfer matrices in symmetry problems
• scientific article; zbMATH DE number 92720
• The Inverse Method for Transfer Matrices
• scientific article; zbMATH DE number 194990

Cites work

• A sextic formalism for three–dimensional elastodynamics of cylindrically anisotropic radially inhomogeneous materials
• Accuracy and Stability of Numerical Algorithms
• scientific article; zbMATH DE number 3135277 (Why is no real title available?)
• scientific article; zbMATH DE number 3756878 (Why is no real title available?)
• On dispersion relations of waves in multilayered magneto-electro-elastic plates
• Orthotropic elastic media having a closed form expression of the Green tensor
• Resolvent method for calculating dispersion spectra of the shear waves in the phononic plates and waveguides
• The quadratic eigenvalue problem
• Three-dimensional Green's functions for infinite anisotropic piezoelectric media
• Wave impedance matrices for cylindrically anisotropic radially inhomogeneous elastic solids
• Wave propagation in magneto-electro-elastic multilayered plates

Cited in

(5)

• Transfer function matrix approach to decoupling problem with stability
• Quantum thermostatted disordered systems and sensitivity under compression
• Unitarity and symmetries of the multicomponent scattering matrix
• scientific article; zbMATH DE number 2156328 (Why is no real title available?)
• Photonic properties of 1D multilayered structures based on quasiperiodic Rudin-Shapiro sequence

Describes a project that uses

Uses Software

• mctoolbox

This page was built for publication: Relations between transfer matrices and numerical stability analysis to avoid the \(\Omega d\) problem

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5261727)