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

DOI10.1137/140993442zbMATH Open1322.34097arXiv1503.09038OpenAlexW1692423489MaRDI QIDQ5261727FDOQ5261727

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

Publication date: 7 July 2015

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

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.

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

Recommendations

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

zbMATH Keywords

numerical stabilitymatrix Sturm-Liouville problemtransfer matrixquadratic eigenvalues\(\Omega d\) problem

Mathematics Subject Classification ID

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

Cites Work

• Accuracy and Stability of Numerical Algorithms
• The quadratic eigenvalue problem
• Wave propagation in magneto-electro-elastic multilayered plates
• Wave Impedance Matrices for Cylindrically Anisotropic Radially Inhomogeneous Elastic Solids
• On dispersion relations of waves in multilayered magneto-electro-elastic plates
• Title not available (Why is that?)
• Resolvent Method for Calculating Dispersion Spectra of the Shear Waves in the Phononic Plates and Waveguides
• A sextic formalism for three–dimensional elastodynamics of cylindrically anisotropic radially inhomogeneous materials
• Three-dimensional Green's functions for infinite anisotropic piezoelectric media
• Orthotropic elastic media having a closed form expression of the Green tensor
• Title not available (Why is that?)

Cited In (5)

• Title not available (Why is that?)
• Unitarity and symmetries of the multicomponent scattering matrix
• Quantum thermostatted disordered systems and sensitivity under compression
• Transfer function matrix approach to decoupling problem with stability
• Photonic properties of 1D multilayered structures based on quasiperiodic Rudin-Shapiro sequence

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)