scholarly journals Identifying linear absolute instabilities from differential eigenvalue problems using sensitivity analysis

2019 ◽  
Vol 870 ◽  
pp. 941-969 ◽  
Author(s):  
L. S. de B. Alves ◽  
S. C. Hirata ◽  
M. Schuabb ◽  
A. Barletta
Keyword(s):  
Sensitivity Analysis ◽  
Dispersion Relation ◽  
Eigenvalue Problem ◽  
Group Velocity ◽  
Algebraic System ◽  
Eigenvalue Problems ◽  
Computational Cost ◽  
Root Finding ◽  
Alternative Methodology ◽  
Zero Group

Identifying the convective/absolute instability nature of a local base flow requires an analysis of its linear impulse response. One must find the appropriate singularity in the eigenvalue problem with complex frequencies and wavenumbers and prove causality. One way to do so is to show that the appropriate integration contour of this response, a steepest decent path through the relevant singularity, exists. Due to the inherent difficulties of such a proof, one often verifies instead whether this singularity satisfies the collision criterion. In other words, one must show that the branches involved in the formation of this singularity come from distinct halves of the complex wavenumber plane. However, this graphical search is computationally intensive in a single plane and essentially prohibitive in two planes. A significant computational cost reduction can be achieved when root finding procedures are applied instead of graphical ones to search for singularities. They focus on locating these points, with causality being verified graphically a posteriori for a small parametric sample size. The use of root-finding procedures require auxiliary equations, often derived by applying the zero group velocity conditions to the dispersion relation. This relation, in turn, is derived by applying matrix forming to the differential eigenvalue problem and taking the determinant of the resulting system of algebraic equations. Taking the derivative of the dispersion relation with respect to the wavenumbers generates the auxiliary equations. If the algebraic system is decoupled, this derivation is straightforward. However, its computational cost is often prohibitive when the algebraic system is coupled. Other methods exist, but often they can also be too costly and/or not reliable for two wavenumber plane searches. This paper describes an alternative methodology based on sensitivity analysis and adjoints that allow the zero group velocity conditions to be applied directly to the differential eigenvalue problem. In doing so, the direct and auxiliary differential eigenvalue problems can be solved simultaneously using standard shooting methods to directly locate singularities. Auxiliary dispersion relations no longer have to be derived, although it is shown that they are the algebraic equivalent of the auxiliary differential eigenvalue problems obtained by this alternative methodology. Using the latter dramatically reduces computational costs. The search for arbitrary singularities is then not only accelerated in single wavenumber planes but it also becomes viable in two wavenumber planes. Finally, the new method also allows group velocity calculations, greatly facilitating the verification of causality. Several test cases are presented to illustrate the capabilities of this new method.

scholarly journals Sensitivity analysis of waveguide eigenvalue problems

2011 ◽  
Vol 9 ◽  
pp. 85-89 ◽  
Author(s):  
N. Burschäpers ◽  
S. Fiege ◽  
R. Schuhmann ◽  
A. Walther
Keyword(s):  
Sensitivity Analysis ◽  
Eigenvalue Problem ◽  
Eigenvalue Problems ◽  
Wave Number ◽  
Design Parameters ◽  
Dielectric Waveguides ◽  
Integration Technique ◽  
Numerical Examples ◽  
Numerical Effort ◽  
Complex Solutions

Abstract. We analyze the sensitivity of dielectric waveguides with respect to design parameters such as permittivity values or simple geometric dependencies. Based on a discretization using the Finite Integration Technique the eigenvalue problem for the wave number is shown to be non-Hermitian with possibly complex solutions even in the lossless case. Nevertheless, the sensitivity can be obtained with negligible numerical effort. Numerical examples demonstrate the validity of the approach.

Transparent boundary conditions for discretized non-Hermitian eigenvalue problems

2021 ◽  
pp. 2150138
Author(s):  
Jonathan Heinz ◽  
Miroslav Kolesik
Keyword(s):  
Quantum Mechanics ◽  
Eigenvalue Problem ◽  
Boundary Conditions ◽  
Eigenvalue Problems ◽  
Complex Scaling ◽  
Transparent Boundary Conditions ◽  
Drastic Reduction ◽  
Optical Cavities ◽  
Energy Dependent ◽  
Lossy Waveguides

A method is presented for transparent, energy-dependent boundary conditions for open, non-Hermitian systems, and is illustrated on an example of Stark resonances in a single-particle quantum system. The approach provides an alternative to external complex scaling, and is applicable when asymptotic solutions can be characterized at large distances from the origin. Its main benefit consists in a drastic reduction of the dimesnionality of the underlying eigenvalue problem. Besides application to quantum mechanics, the method can be used in other contexts such as in systems involving unstable optical cavities and lossy waveguides.

scholarly journals Adjoint-based parametric sensitivity analysis for swirling M-flames

2018 ◽  
Vol 859 ◽  
pp. 516-542 ◽  
Author(s):  
Calum S. Skene ◽  
Peter J. Schmid
Keyword(s):  
Sensitivity Analysis ◽  
Frequency Response ◽  
Stokes Equations ◽  
Numerical Study ◽  
Computational Cost ◽  
Base Flow ◽  
Sensitivity Analyses ◽  
Perturbation Expansions ◽  
Mean Flow ◽  
Adjoint Solutions

A linear numerical study is conducted to quantify the effect of swirl on the response behaviour of premixed lean flames to general harmonic excitation in the inlet, upstream of combustion. This study considers axisymmetric M-flames and is based on the linearised compressible Navier–Stokes equations augmented by a simple one-step irreversible chemical reaction. Optimal frequency response gains for both axisymmetric and non-axisymmetric perturbations are computed via a direct–adjoint methodology and singular value decompositions. The high-dimensional parameter space, containing perturbation and base-flow parameters, is explored by taking advantage of generic sensitivity information gained from the adjoint solutions. This information is then tailored to specific parametric sensitivities by first-order perturbation expansions of the singular triplets about the respective parameters. Valuable flow information, at a negligible computational cost, is gained by simple weighted scalar products between direct and adjoint solutions. We find that for non-swirling flows, a mode with azimuthal wavenumber $m=2$ is the most efficiently driven structure. The structural mechanism underlying the optimal gains is shown to be the Orr mechanism for $m=0$ and a blend of Orr and other mechanisms, such as lift-up, for other azimuthal wavenumbers. Further to this, velocity and pressure perturbations are shown to make up the optimal input and output showing that the thermoacoustic mechanism is crucial in large energy amplifications. For $m=0$ these velocity perturbations are mainly longitudinal, but for higher wavenumbers azimuthal velocity fluctuations become prominent, especially in the non-swirling case. Sensitivity analyses are carried out with respect to the Mach number, Reynolds number and swirl number, and the accuracy of parametric gradients of the frequency response curve is assessed. The sensitivity analysis reveals that increases in Reynolds and Mach numbers yield higher gains, through a decrease in temperature diffusion. A rise in mean-flow swirl is shown to diminish the gain, with increased damping for higher azimuthal wavenumbers. This leads to a reordering of the most effectively amplified mode, with the axisymmetric ($m=0$) mode becoming the dominant structure at moderate swirl numbers.

Dispersion relations of dust lattice waves in two-dimensional honeycomb configuration

2013 ◽  
Vol 79 (5) ◽  
pp. 629-633
Author(s):  
B. FAROKHI
Keyword(s):  
Dispersion Relation ◽  
Group Velocity ◽  
Hexagonal Lattice ◽  
Dispersion Relations ◽  
Two Dimensional ◽  
The Third ◽  
Lattice Waves ◽  
The Waves

AbstractThe linear dust lattice waves propagating in a two-dimensional honeycomb configuration is investigated. The interaction between particles is considered up to distance 2a, i.e. the third-neighbor interactions. Longitudinal and transverse (in-plane) dispersion relations are derived for waves in arbitrary directions. The study of dispersion relations with more neighbor interactions shows that in some cases the results change physically. Also, the dispersion relation in the different direction displays anisotropy of the group velocity in the lattice. The results are compared with dispersion relations of the waves in the hexagonal lattice.

scholarly journals Linearization schemes for Hermite matrix polynomials

2015 ◽  
Author(s):  
Nikta Shayanfar ◽  
Heike Fassbender
Keyword(s):  
Eigenvalue Problem ◽  
Eigenvalue Problems ◽  
Matrix Polynomial ◽  
Multiple Input Multiple Output ◽  
Matrix Pencil ◽  
Least Square ◽  
Coupled Systems ◽  
Matrix Polynomials ◽  
Polynomial Eigenvalue Problem ◽  
The Matrix

The polynomial eigenvalue problem is to find the eigenpair of $(\lambda,x) \in \mathbb{C}\bigcup \{\infty\} \times \mathbb{C}^n \backslash \{0\}$ that satisfies $P(\lambda)x=0$, where $P(\lambda)=\sum_{i=0}^s P_i \lambda ^i$ is an $n\times n$ so-called matrix polynomial of degree $s$, where the coefficients $P_i, i=0,\cdots,s$, are $n\times n$ constant matrices, and $P_s$ is supposed to be nonzero. These eigenvalue problems arise from a variety of physical applications including acoustic structural coupled systems, fluid mechanics, multiple input multiple output systems in control theory, signal processing, and constrained least square problems. Most numerical approaches to solving such eigenvalue problems proceed by linearizing the matrix polynomial into a matrix pencil of larger size. Such methods convert the eigenvalue problem into a well-studied linear eigenvalue problem, and meanwhile, exploit and preserve the structure and properties of the original eigenvalue problem. The linearizations have been extensively studied with respect to the basis that the matrix polynomial is expressed in. If the matrix polynomial is expressed in a special basis, then it is desirable that its linearization be also expressed in the same basis. The reason is due to the fact that changing the given basis ought to be avoided \cite{H1}. The authors in \cite{ACL} have constructed linearization for different bases such as degree-graded ones (including monomial, Newton and Pochhammer basis), Bernstein and Lagrange basis. This contribution is concerned with polynomial eigenvalue problems in which the matrix polynomial is expressed in Hermite basis. In fact, Hermite basis is used for presenting matrix polynomials designed for matching a series of points and function derivatives at the prescribed nodes. In the literature, the linearizations of matrix polynomials of degree $s$, expressed in Hermite basis, consist of matrix pencils with $s+2$ blocks of size $n \times n$. In other words, additional eigenvalues at infinity had to be introduced, see e.g. \cite{CSAG}. In this research, we try to overcome this difficulty by reducing the size of linearization. The reduction scheme presented will gradually reduce the linearization to its minimal size making use of ideas from \cite{VMM1}. More precisely, for $n \times n$ matrix polynomials of degree $s$, we present linearizations of smaller size, consisting of $s+1$ and $s$ blocks of $n \times n$ matrices. The structure of the eigenvectors is also discussed.

Sign in / Sign up

Export Citation Format

Share Document