Generalized Eigenvalue Problem Research Articles

A computationally efficient [formula omitted]-spectral form for partial dispersion analyses within the wave finite element framework

This paper addresses the computation of frequency-dependent dispersion curves (i.e., k(ω)) and wave modes within the framework of the Wave Finite Element Method (WFEM) and in the context of high-dimensional periodic unit cell models. Numerous applications, ranging from phononics to vibroacoustics, now rely on dispersion analyses or wave expansion over a subset of eigensolutions –complex wavenumbers and Bloch waves– resulting from the resolution of an eigenvalue problem with a T-palindromic quadratic structure (T-PQEP). To exploit the structure of finite element models, various structure-preserving linearizations such as the Zhong-Williams and the (S+S−1)-transform have already been developed to achieve partial wave resolution of large T-PQEP, primarily targeting the dominating (least decaying) waves. In this paper we derive an alternative linearization of the T-PQEP for the k(ω) problem, which leads to enhanced targeting of the eigenvalues around the unit circle and reduces the inaccuracies induced by root multiplicity. A specific form of the problem is then proposed as an optimal compromise between ease of implementation, numerical stability, convergence and accuracy enhancement. The performance of our proposed linearization is compared against existing ones across various iterative eigensolvers, since the generalized eigenvalue problems involve complex non-hermitian matrices, which are not extensively included in eigensolvers. Results indicate that the proposed linearization should be favoured for the WFEM, as it provides numerical enhancements in dispersion and wave vectors computation for large eigenvalue problems, as well as for further wave expansion applications.

Fluid-structure coupling analysis in liquid-filled containers using scaled boundary finite element method

In this study, a semi-analytical model is developed to investigate the fluid-structure coupling characteristics of liquid sloshing in an elastic rectangular container subjected to horizontal external excitation based on the scaled boundary finite element method (SBFEM) for the first time. The fluid inside the container is assumed to be incompressible, inviscid and irrotational, with the hydrodynamic pressure chosen as independent nodal variable in the governing equations. The container walls are considered as cantilever beams. The coupled fluid-structure system is initially divided into the structural domain and fluid domain, after which the SBFEM is employed to obtain the governing equations for each sub-domain. In the framework of the SBFEM, only the boundary of each sub-domain, rather than the entire computational domain, needs to be meshed and discretized. This method reduces the spatial dimension of the problem by one and offers an efficient approach to model the computational domain, while allowing for analytical formulations to be derived for the internal of the domain, resulting in an accurate description of the field variables. The fundamental equation of the entire coupled fluid-structure system is then assembled by performing the equilibrium condition and compatibility condition to ensure the balance of interaction forces at the interface between container walls and the liquid. The free vibrations analysis of the fluid-structure coupling system is solved by utilizing the generalized eigenvalue problem, and the transient dynamic response is determined using the synchronous solution algorithm in conjunction with the implicit-implicit scheme of the Newmark method. To validate the excellent accuracy and stability of the proposed formulation, several numerical examples are presented to investigate the free vibration and transient dynamic characteristics for the fluid-structure coupling problem. The obtained results show good agreement with reference solutions available in the literature. Additionally, the effects of geometrical and material parameters on the system responses are examined and discussed.

Analysis of eigenvalue condition numbers for a class of randomized numerical methods for singular matrix pencils

The numerical solution of the generalized eigenvalue problem for a singular matrix pencil is challenging due to the discontinuity of its eigenvalues. Classically, such problems are addressed by first extracting the regular part through the staircase form and then applying a standard solver, such as the QZ algorithm, to that regular part. Recently, several novel approaches have been proposed to transform the singular pencil into a regular pencil by relatively simple randomized modifications. In this work, we analyze three such methods by Hochstenbach, Mehl, and Plestenjak that modify, project, or augment the pencil using random matrices. All three methods rely on the normal rank and do not alter the finite eigenvalues of the original pencil. We show that the eigenvalue condition numbers of the transformed pencils are unlikely to be much larger than the δ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\delta $$\\end{document}-weak eigenvalue condition numbers, introduced by Lotz and Noferini, of the original pencil. This not only indicates favorable numerical stability but also reconfirms that these condition numbers are a reliable criterion for detecting simple finite eigenvalues. We also provide evidence that, from a numerical stability perspective, the use of complex instead of real random matrices is preferable even for real singular matrix pencils and real eigenvalues. As a side result, we provide sharp left tail bounds for a product of two independent random variables distributed with the generalized beta distribution of the first kind or Kumaraswamy distribution.

Intuitionistic fuzzy generalized eigenvalue proximal support vector machine

Generalized eigenvalue proximal support vector machine (GEPSVM) has attracted widespread attention due to its simple architecture, rapid execution, and commendable performance. GEPSVM gives equal significance to all samples, thereby diminishing its robustness and efficacy when confronted with real-world datasets containing noise and outliers. In order to reduce the impact of noises and outliers, we propose a novel intuitionistic fuzzy generalized eigenvalue proximal support vector machine (IF-GEPSVM). The proposed IF-GEPSVM assigns the intuitionistic fuzzy score to each training sample based on its location and surroundings in the high-dimensional feature space by using a kernel function. The solution of the IF-GEPSVM optimization problem is obtained by solving a generalized eigenvalue problem. Further, we propose an intuitionistic fuzzy improved generalized eigenvalue proximal support vector machine (IF-IGEPSVM) by solving the standard eigenvalue decomposition resulting in simpler optimization problems with less computation cost which leads to an efficient intuitionistic fuzzy-based model. We conduct a comprehensive evaluation of the proposed IF-GEPSVM and IF-IGEPSVM models on UCI and KEEL benchmark datasets. Moreover, to evaluate the robustness of the proposed IF-GEPSVM and IF-IGEPSVM models, label noise is introduced into some UCI and KEEL datasets. The experimental findings showcase the superior generalization performance of the proposed IF-GEPSVM and IF-IGEPSVM models when compared to the existing baseline models, both with and without label noise. Our experimental results, supported by rigorous statistical analyses, confirm the superior generalization abilities of the proposed IF-GEPSVM and IF-IGEPSVM models over the baseline models. Furthermore, we implement the proposed IF-GEPSVM and IF-IGEPSVM models on the USPS recognition dataset, yielding promising results that underscore the models’ effectiveness in practical and real-world applications. The source code of the proposed IF-GEPSVM and IF-IGEPSVM models are available at https://github.com/mtanveer1/IF-GEPSVM.

Study of circulating liquid fuel in a 1D critical system with thermal feedback

This research focuses on the description and modeling of a one-dimensional molten salt reactor (MSR), in the presence of thermal feedback. Following the example of previous works, a simple one-dimensional system is proposed, describing a molten salt reactor with a main neutron-multiplying zone called core and a recirculation loop where the salt cools down. Specific attention is paid to the precursors’ drift by modifying the neutron balance equation. Liquid nuclear fuels are characterized by a high volumetric expansion coefficient in comparison to customary solid fuels. Therefore, a strong coupling between neutronics and thermal-hydraulics is expected. As a consequence, a highly negative density coefficient characterizes the thermal feedback on the neutron reactivity. The precursor equation is here inverted analytically and combined with the neutron balance equation to obtain a generalized eigenvalue problem with the neutron flux distribution as the unknown. The balance equations are derived by finite volume integration over a discretized mesh, and the coupling between the two physical models is treated by Picard iterations. The numerical solution is finally extended to time-dependent calculations and compared to an analytical work for a one-dimensional circulating fuel reactor already existing in the literature.

Finding eigenvectors with a quantum variational algorithm

This paper presents a hybrid variational quantum algorithm that finds a random eigenvector of a unitary matrix with a known quantum circuit. The algorithm is based on the SWAP test on trial states generated by a parametrized quantum circuit. The eigenvector is described by a compact set of classical parameters that can be used to reproduce the found approximation to the eigenstate on demand. This variational eigenvector finder can be adapted to solve the generalized eigenvalue problem, to find the eigenvectors of normal matrices and to perform quantum principal component analysis on unknown input mixed states. These algorithms can all be run with low-depth quantum circuits, suitable for an efficient implementation on noisy intermediate-scale quantum computers and, with some restrictions, on linear optical systems. In full-scale quantum computers, where there might be optimization problems due to barren plateaus in larger systems, the proposed algorithms can be used as a primitive to boost known quantum algorithms. Limitations and potential applications are discussed.

Vibration of embedded restrained composite tube shafts with nonlocal and strain gradient effects

Torsional vibration response of a circular nanoshaft, which is restrained by the means of elastic springs at both ends, is a matter of great concern in the field of nano-/micromechanics. Hence, the complexities arising from the deformable boundary conditions present a formidable obstacle to the attainment of closed-form solutions. In this study, a general method is presented to calculate the torsional vibration frequencies of functionally graded porous tube nanoshafts under both deformable and rigid boundary conditions. Classical continuum theory, upgraded with nonlocal strain gradient elasticity theory, is employed to reformulate the partial differential equation of the nanoshaft. First, torsional vibration equation based on the nonlocal strain gradient theory is derived for functionally graded porous nanoshaft embedded in an elastic media via Hamilton’s principle. The ordinary differential equation is found by discretizing the partial differential equation with the separation of variables method. Then, Fourier sine series is used as the rotation function. The necessary Stokes' transformation is applied to establish the general eigenvalue problem including the different parameters. For the first time in the literature, a solution that can analyze the torsional vibration frequencies of functionally graded porous tube shafts embedded in an elastic media under general (elastic and rigid) boundary conditions on the basis of nonlocal strain gradient theory is presented in this study. The results obtained show that while the increase in the material length scale parameter, elastic media and spring stiffnesses increase the frequencies of nanoshafts, the increase in the nonlocal parameter and functionally grading index values decreases the frequencies of nanoshafts. The detailed effects of these parameters are discussed in the article.

R I biorthogonal polynomials of Hahn type

ABSTRACT A finite family of R I polynomials is introduced and studied. It consists in a set of polynomials of 3 F 2 form whose biorthogonality to an ensemble of rational functions is spelled out. These polynomials are shown to satisfy two generalized eigenvalue problems: in addition to their recurrence relation of R I type, they are also found to obey a difference equation. Underscoring this bispectrality is a triplet of operators with tridiagonal actions. The algebra associated to these operators is provided.

New vector transport operators extending a Riemannian CG algorithm to generalized Stiefel manifold with low-rank applications

This paper proposes two innovative vector transport operators, leveraging the Cayley transform, for the generalized Stiefel manifold embedded with a non-standard metric. Specifically, it introduces the differentiated retraction and an approximation of the Cayley transform to the differentiated matrix exponential. These vector transports are demonstrated to satisfy the Ring–Wirth non-expansive condition under non-standard metrics, and one of them is also isometric. Building upon the novel vector transport operators, we extend the modified Polak–Ribière–Polyak (PRP) conjugate gradient method to the generalized Stiefel manifold. Under a non-monotone line search condition, we prove our algorithm globally converges to a stationary point. The efficiency of the proposed vector transport operators is empirically validated through numerical experiments involving generalized eigenvalue problems and canonical correlation analysis.

Conical-Shaped Shells of Non-Uniform Thickness Vibration Analysis Using Higher-Order Shear Deformation Theory

Conical-Shaped Shells of Non-Uniform Thickness Vibration Analysis Using Higher-Order Shear Deformation Theory

The pseudospectra of black holes in AdS

We study the stability of quasinormal modes (QNMs) in electrically charged black brane spacetimes that asymptote to AdS by means of the pseudospectrum. Methodologically, we adopt ingoing Eddington-Finkelstein coordinates to cast QNMs in terms of a generalised eigenvalue problem involving a non-selfadjoint operator; this simplifies the computation significantly in comparison with previous results in the literature. Our analysis reveals spectral instability for (neutral) scalar as well as gravitoelectric perturbations. This indicates that the equilibration process of perturbed black branes is sensitive to external perturbations. Particular attention is given on the hydrodynamic modes, which are found to be the least unstable. In contrast with computations in hyperboloidal coordinates, we find that the pseudospectral contour lines cross to the upper half plane. This indicates the existence of pseudo-resonances as well as the possibility of transient instabilities. We also investigate the asymptotic structure of pseudospectral contour levels and we find remarkable universality across all sectors, persistent in the extremal limit.

On the magnetohydrodynamic stability of channel flow of Navier–Stokes–Voigt fluid

The linear stability of a pressure-driven channel flow of an electrically conducting Navier–Stokes–Voigt type of the viscoelastic fluid subject to a transverse magnetic field is investigated. The validity of Squire's theorem is proved, and the generalized eigenvalue problem for two-dimensional modes is obtained by adopting the Galerkin method, which is subsequently solved using the QZ-algorithm. Although the base flow retains its Newtonian fluid characteristics, the noticeable influence of the Kelvin–Voigt parameter in conjunction with the Hartmann number on the stability of fluid flow is perceived. Instability is exclusively identified within a specific range of the Kelvin–Voigt parameter, markedly affected by the Hartmann number. Furthermore, closed neutral stability curves arise, indicating the requirement of two values of the Reynolds number to completely assess the linear stability criteria, in contrast to the typical single value observed in the case of Newtonian fluids.

Tensor networks for solving the time-independent Boltzmann neutron transport equation

Tensor network techniques, known for their low-rank approximation ability that breaks the curse of dimensionality, are emerging as a foundation of new mathematical methods for ultra-fast numerical solutions of high-dimensional Partial Differential Equations (PDEs). Here, we present a mixed Tensor Train (TT)/Quantized Tensor Train (QTT) approach for the numerical solution of time-independent Boltzmann Neutron Transport equations (BNTEs) in Cartesian geometry. Discretizing a realistic three-dimensional (3D) BNTE by (i) diamond differencing in space, (ii) multigroup-in-energy, and (iii) discrete ordinates collocation in angle leads to large generalized eigenvalue problems that generally require a matrix-free approach and large computer clusters. Starting from this discretization, we construct a TT representation of the PDE fields and discrete operators, followed by a QTT representation of the TT cores. We then solve the tensorized generalized eigenvalue problem using a fixed-point scheme with tensor network optimization techniques. We validate our approach by applying the method to two examples of 3D neutron transport problems, currently solved by the Los Alamos National Laboratory PARallel TIme-dependent SN (PARTISN) solver.1 We demonstrate that our TT/QTT method, executed on a standard desktop computer, leads to large compression. This allows for the storage of terrabyte-sized neutron angular flux eigenvectors in megabytes. Additionally, we create megabyte-sized full access TT representations of yottabyte-sized transport matrix operators. By leveraging the TT operators and solution methods, we obtain a 7500 times speedup when compared to the PARTISN solution time with an error of less than 10−5.

Stability of a photosurfactant-laden viscous liquid thread under illumination

This paper investigates the effects of a light-actuated photosurfactant on the canonical problem of the linear stability of a viscous thread surrounded by a dynamically passive fluid. A model consisting of the Navier–Stokes equations and a set of molar concentration equations is presented that capture light-induced switching between two stable surfactant isomer states, trans and cis. These two states display significantly different interfacial properties, allowing for some external control of the stability behaviour of the thread via incident light. Normal modes are used to generate a generalized eigenvalue problem for the growth rate which is solved with a hybrid analytical and numerical method. The results are validated with appropriate analytical solutions of increasing complexity, beginning with a solution to a clean interface, then analytical solutions for one insoluble surfactant, one soluble surfactant and a special case of two photosurfactants with a spatially uniform undisturbed state. Presenting each of these cases allows for a holistic discussion of the effect of surfactants in general on the stability of a liquid thread. Finally, the numerical solutions in the presence of two photosurfactants that display radially non-uniform undisturbed states are presented, and details of the impact of the illumination on the linear stability of the thread are discussed.

Non-Bloch band theory of generalized eigenvalue problems

Non-Bloch band theory of generalized eigenvalue problems

An application of stability charts to prediction of buckling instability in tapered columns via Galerkin’s method

This work tackles the mathematical modeling of buckling problem to obtain their critical loads in tapered columns subjected to concentrated and axial distributed loads. The governing model is a general eigenvalue problem that has no exact solution due to some new terms included. A semi-analytical technique satisfying the boundary conditions is proposed for the solution procedure. The minimum residual Galerkin’s method is suggested due to its effectiveness as a semi-analytical tool for the buckling problem to obtain the buckling shape modes by using admissible periodic functions. The study investigates the buckling instability and the responses of tapered columns with different periodic trial shape functions as approximations to the exact solutions. Based on the eigenvalue problem, Galerkin’s method is employed to obtain the transition curves to represent the critical loads. The stability charts (Ince–Strutt diagrams) among the parameters of the problem are proposed to explain the elastic stability of different tapered columns subjected to different shapes of cross sections and distributed weights. Consequently, the influences of the included parameters on the critical buckling loads are discussed. Among the different tapered columns presented, some parameters in the proposed distributions have a big influence on the critical buckling load and the creation of the instability regions in the chart for the clamped-clamped boundary conditions. The results are verified using the analytical solutions for some specific known problems.

Design sensitivity analysis of modal frequencies of elastic structures submerged in an infinite fluid domain

SummaryThis paper presents a numerical approach for the design sensitivity analysis of modal frequencies of elastic structures submerged in an unbounded heavy fluid domain. Because the feedback of sound pressure onto the fluid‐loaded structures has to be taken into account in this case, a fully coupled scheme which combines the structural finite element method (FEM) and the acoustic boundary element method (BEM) is employed in the numerical simulation. The resulting nonlinear eigenvalue problem created by the coupled finite element and boundary element (FE‐BE) method is converted into a generalized eigenvalue problem (GEP) of reduced dimension by using a contour integral method. The design sensitivity formulation for non‐repeated eigenvalues is derived based on an adjoint method which uses both the left and right eigenvectors, and a small GEP is formed to calculate the gradients of repeated eigenvalues. The Burton‐Miller formulation is employed to shift the fictitious eigenfrequencies and their derivatives, and the parameter setting which is beneficial to the filtering of fictitious eigenfrequencies is investigated. Numerical examples are given to verify the accuracy and effectiveness of the developed method.

Effects of Rotation and Magnetic Field on Rayliegh Benard Convection

In this paper, a numerical method based on the Chebyshev tau method is applied to analyze the effects of rotation and magnetic fields on Rayleigh-Bénard convection. The rotation and magnetic fields are assumed to be parallel to the vertical direction. The perturbation equations and boundary conditions are analyzed using normal mode analysis. The equations are then converted into a non-dimensional form and transformed into a generalized eigenvalue problem of the form AX=RBX, where R represents the eigenvalue corresponding to the Rayleigh number. The MATLAB software package is utilized to determine the relationship between the Rayleigh number and the Taylor number (rate of rotation), as well as the relationship between the Rayleigh number and the magnetic parameter (strength of the magnetic field) for different boundary conditions (free-free, rigid-rigid, or one free and the other rigid). The numerical and graphical results are presented and found to be in full agreement with the results obtained from previous analytical and numerical studies of the problem.

Effect of Oldroyd-B fluid on thermal convection stability in a bidispersive porous media with relatively large macropores

This study investigates the effect of Oldroyd-B fluids on the thermal convection instability in a bi-disperse porous medium with relatively large macropores. The modified momentum equation of Oldroyd-B fluid in the macro- and micropores is first derived, assuming local thermal equilibrium. The Orr–Sommerfeld (O–S) equation is derived using the new modified incompressible Navier–Stokes equations of Oldroyd-B fluid. The Chebyshev collocation method is employed to solve the generalized eigenvalue problem, and the neutral stability curves and temporal growth rates are obtained. The results demonstrate that the relaxation time enhances the convective instability, while the retardation time suppresses it. Furthermore, the effect of relaxation time is generally small for smaller wave number, but the effect of relaxation time is significant for large wave number, and this effect becomes accentuated as the wave number increases. The retardation time does not have any influence when the wavenumber is small, and the influence of retardation time becomes small with the large relaxation time. Moreover, the convection instability is enhanced with the higher micropermeability and a higher ratio of momentum transfer coefficient between the macro- and microphases will suppress the instability.

A family of Nonoverlapping Spectral Additive Schwarz methods (NOSAS) and their economic versions

NOSAS (Nonoverlapping Spectral Additive Schwarz) methods are nonoverlapping iterative domain decomposition methods for efficiently solving large sparse linear systems arising from elliptic problems with heterogeneous coefficients. NOSAS methods are of additive Schwarz types with local Dirichlet solver on each nonoverlapping subdomain Ωi, plus a coarse problem on the interface of the subdomains. The coarse problem has global and local interaction components which are associated respectively to low-frequency and high-frequency modes obtained locally from a generalized eigenvalue problem defined on each subdomain’s interface. In this paper, we propose several variants of the NOSAS methods by selecting different generalized eigenvalue problems and different coarse space extensions. Additionally, we propose an economic variant to reduce the complexity of these generalized eigenvalue problems. The analysis of all the variants is established, and numerical tests are provided.