Asymptotics of solutions to the spectral elasticity problem for a spatial body with a thin coupler

We construct asymptotics for the eigenvalues and vector eigenfunctions of the elasticity problem for an anisotropic body with a thin coupler (of diameter h) attached to its surface. In the spectrum we select two series of eigenvalues with stable asymptotics. The first series is formed by eigenvalues O(h ²) corresponding to the transverse oscillations of the rod with rigidly fixed ends, while the second is generated by the longitudinal oscillations and twisting of the rod, as well as eigenoscillations of the body without the coupler. We check the convergence theorem for the first series and derive the error estimates for both series.     

Perturbation Analysis for the Eigenvalue Problem of a Formal Product of Matrices

We study the perturbation theory for the eigenvalue problem of a formal matrix product A ₁ ^{s 1} ··· A _p ^{s p}, where all A _k are square and s _k {–1, 1}. We generalize the classical perturbation results for matrices and matrix pencils to perturbation results for generalized deflating subspaces and eigenvalues of such formal matrix products. As an application we then extend the structured perturbation theory for the eigenvalue problem of Hamiltonian matrices to Hamiltonian/skew-Hamiltonian pencils.     

Eigenvalue asymptotics for randomly perturbed non-selfadjoint operators

We consider quite general h-pseudodifferential operators on R ⁿ with small random perturbations and show that in the limit h → 0 the eigenvalues are distributed according to a Weyl law with a probabality that tends to 1. The first author has previously obtained a similar result in dimension 1. Our class of perturbations is different.     

Analysis of Two Linearized Problems Modeling Viscous Two‐Layer Flows

Two problems that appear in the linearization of certain free boundary value problems of the hydrodynamics of two viscous fluids are studied in the strip‐like domain Π = {x = (x₁, x₂) ∈ ℝ² : x₁ ∈ ℝ¹, (0 < x₂ < h_*) ∨ (h_* < x₂ < 1)}. The first problem arises in the linearization of a two‐layer flow down a geometrically perturbed inclined plane. The second one appears after the linearization of a two‐layer flow in a geometrically perturbed inclined channel with one moving (smooth) wall. For this purpose the unknown flow domain was mapped onto the double strip Π. The arising linear elliptic problems contain additional unknown functions in the boundary conditions. The paper is devoted to the investigation of these boundary problems by studying the asymptotics of the eigenvalues of corresponding operator pencils. It can be proved that the boundary value problems are uniquely solvable in weighted Sobolev spaces with exponential weight. The study of the full (nonlinear) free boundary value problems will be the topic of a forthcoming paper.     

A degenerate kernel method for eigenvalue problems of compact integral operators

We consider the approximation of eigenfunctions of a compact integral operator with a smooth kernel by a degenerate kernel method. By interpolating the kernel of the integral operator in both the variables, we prove that the error bounds for eigenvalues and for the distance between the spectral subspaces are of the orders h ^2r and h ^r respectively. By iterating the eigenfunctions we show that the error bounds for eigenfunctions are of the orders h ^2r . We give the numerical results.      

On spectra of a certain class of quadratic operator pencils with one-dimensional linear part

We consider a class of quadratic operator pencils that occur in many problems of physics. The part of such a pencil linear with respect to the spectral parameter describes viscous friction in problems of small vibrations of strings and beams. Patterns in the location of eigenvalues of such pencils are established. If viscous friction (damping) is pointwise, then the operator in the linear part of the pencil is one-dimensional. For this case, rules in the location of purely imaginary eigenvalues are found. Published in Ukrains’kyi Matematychnyi Zhurnal, Vol. 59, No. 5, pp. 702–716, May, 2007.     

On the vibrations of a plate with a concentrated mass and very small thickness

We consider the vibrations of an elastic plate that contains a small region whose size depends on a small parameter ε. The density is of order O(ε^–m) in the small region, the concentrated mass, and it is of order O(1) outside; m is a positive parameter. The thickness plate h being fixed, we describe the asymptotic behaviour, as ε→O, of the eigenvalues and eigenfunctions of the corresponding spectral problem, depending on the value of m: Low‐ and high‐frequency vibrations are studied for m>2. We also consider the case where the thickness plate h depends on ε; then, different values of m are singled out. Copyright © 2003 John Wiley & Sons, Ltd.     

On linear vibrational systems with one dimensional damping II

We are interested in the quadratic eigenvalue problem of damped oscillations where the damping matrix has dimension one. This describes systems with one point damper. A generic example is a linearn-mass oscillator fixed on one end and damped on the other end. We prove that in this case the system parameters (mass and spring constants) are uniquely (up to a multiplicative constant) determined by any given set of the eigenvalues in the left half plane. We also design an effective construction of the system parameters from the spectral data. We next propose an efficient method for solving the Ljapunov equation generated by arbitrary stiffness and mass matrices and a one dimensional damping matrix. The method is particularly efficient if the Ljapunov equation has to be solved many times where only the damping dyadic is varied. In particular, the method finds an optimal position of a damper in some 60n ³ operations. We apply this method to our generic example and show, at least numerically, that the damping is optimal (in the sense that the solution of a corresponding Ljapunov equation has a minimal trace) if all eigenvalues are brought together. We include some perturbation results concerning the damping factor as the varying parameter. The results are hoped to be of some help in studying damping matrices of the rank much smaller than the dimension of the problem.     

Approximation du spectre d'un operateur lineaire; application aux operateurs differentiels elliptiques non autoadjoints

To be computed, the eigenvalues of a closed linear operatorT in a Banach space are usually approximated by the eigenvalues ofT h, a linear operator approximatingT in a finite dimensional space (for example, finite difference method, Galerkin method),h is a parameter which tends to 0. This approximation is studied in [2]; stability ofT h implies the continuity of the spectrum ofT h, whenh tends to 0.We present here a new kind of sufficient condition. For that purpose, we disconnect the continuity of the spectrum ofT h into lower and upper semicontinuities. And we give two different criteria for these semi-continuities. Applications to the approximation of nonselfadjoint elliptic operators by finite difference schemes, are given.     

A note on harmonic Ritz values and their reciprocals

This note summarizes an investigation of harmonic Ritz values to approximate the interior eigenvalues of a real symmetric matrix A while avoiding the explicit use of the inverse A^?1. We consider a bounded functional ψ that yields the reciprocals of the harmonic Ritz values of a symmetric matrix A. The crucial observation is that with an appropriate residual s, many results from Rayleigh quotient and Rayleigh–Ritz theory naturally extend. The same is true for the generalization to matrix pencils (A, B) when B is symmetric positive definite. These observations have an application in the computation of eigenvalues in the interior of the spectrum of a large sparse matrix. The minimum and maximum of ψ correspond to the eigenpairs just to the left and right of zero (or a chosen shift). As a spectral transformation, this distinguishes ψ from the original harmonic approach where an interior eigenvalue remains at the interior of the transformed spectrum. As a consequence, ψ is a very attractive vehicle for a matrix‐free, optimization‐based eigensolver. Instead of computing the smallest/largest eigenvalues by minimizing/maximizing the Rayleigh quotient, one can compute interior eigenvalues as the minimum/maximum of ψ. Copyright © 2009 John Wiley & Sons, Ltd.     

Saddle points and overdetermined problems for the Helmholtz equation

In a seminal 1971 paper, James Serrin showed that the only open, smoothly bounded domain in ⁿ on which the positive Dirichlet eigenfunction of the Laplacian has constant (nonzero) normal derivative on the boundary, is then-dimensional ball. The positivity of the eigenfunction is crucial to his proof. To date it is an open conjecture that the same result is true for Dirichlet eigenvalues other than the least. We show that for simply connected, plane domains, the absence of saddle points is a condition sufficient to validate this conjecture. This condition is also sufficient to prove Schiffer's conjecture: the only simply connected planar domain, on the boundary of which a nonconstant Neumann eigenfunction of the Laplacian can take constant value, is the disc.     

Eigenvalue bifurcation in multiparameter families of non-self-adjoint operator matrices

We consider two-point non-self-adjoint boundary eigenvalue problems for linear matrix differential operators. The coefficient matrices in the differential expressions and the matrix boundary conditions are assumed to depend analytically on the complex spectral parameter λ and on the vector of real physical parameters p. We study perturbations of semi-simple multiple eigenvalues as well as perturbations of non-derogatory eigenvalues under small variations of p. Explicit formulae describing the bifurcation of the eigenvalues are derived. Application to the problem of excitation of unstable modes in rotating continua such as spherically symmetric MHD α ²-dynamo and circular string demonstrates the efficiency and applicability of the approach.     

Determination of differential pencils with spectral parameter dependent boundary conditions from interior spectral data

Second‐order differential pencils L(p,q,h₀,h₁,H₀,H₁) on a finite interval with spectral parameter dependent boundary conditions are considered. We prove the following: (i) a set of values of eigenfunctions at the mid‐point of the interval [0,π] and one full spectrum suffice to determine differential pencils L(p,q,h₀,h₁,H₀,H₁); and (ii) some information on eigenfunctions at some an internal point and parts of two spectra suffice to determine differential pencils L(p,q,h₀,h₁,H₀,H₁). Copyright © 2013 The Authors. Mathematical Methods in the Applied Sciences published by John Wiley & Sons, Ltd.     

Homotopy analysis method for multiple solutions of the fractional Sturm-Liouville problems

In this paper, Homotopy Analysis Method (HAM) is applied to numerically approximate the eigenvalues of the fractional Sturm-Liouville problems. The eigenvalues are not unique. These multiple solutions, i.e., eigenvalues, can be calculated by starting the HAM algorithm with one and the same initial guess and linear operator L\mathcal{L}. It can be seen in this paper that the auxiliary parameter (^h/_2p),\hbar, which controls the convergence of the HAM approximate series solutions, has another important application. This important application is predicting and calculating multiple solutions.     

On the p version of the finite element method in the presence of numerical integration

We analyze the error in the p version of the finite element method when the effect of the quadrature error is taken into account. We extend some results by Banerjee and Suri [Math. Comput. 59 , 1–20 (1992)] on the H¹-norm error to the case of the error in the L₂ norm. We investigate three sources of quadrature error that can occur: the error due to the numerical integration of the right-hand side, that due to nonconstant coefficients, and that due to the presence of mapped elements. Presented are various theoretical and computational examples regarding the sharpness of our results. In addition, we make a note on the use of numerical quadrature in conjunction with p-adaptive procedures and on the necessity of overintegration in the h version with linear elements, when the L₂ norm is of interest. © 1993 John Wiley & Sons, Inc.     

A finite-difference approximation for the eigenvalues of the clamped plate

Summary A finite-difference scheme is given for the eigenproblem of the clamped plate. The discrete eigenvalues and eigenvectors are shown to converge to the continuous eigenvalues and eigenvectors likeO(h ²) andO (h ² logh 1/2) respectively.This work, supported by the U.S. Department of the Navy under Contract N 00017-62-C-0604.     

Eigenvalues of Two Parameter Polynomial Operator Pencils of Waveguide Type

The spectral structure of two parameter unbounded operator pencils of waveguide type is studied. Theorems on discreteness of the spectrum for a fixed parameter are proved. Variational principles for real eigenvalues in some parts of the root zones are established. In the case of n = 1 (quadratic pencils) domains containing the spectrum are described (see Fig. 1–3). Conditions in the definition of the pencils of waveguide type arise naturally from physical problems and each of them has a physical meaning. In particular a connection between the energetic stability condition and a perturbation problem for the coefficients is given.     

Oscillation theorems for symplectic difference systems

We consider symplectic difference systems involving a spectral parameter, together with the Dirichlet boundary conditions. The main result of the paper is a discrete version of the so-called oscillation theorem which relates the number of finite eigenvalues less than a given number to the number of focal points of the principal solution of the symplectic system. In two recent papers the same problem was treated and an essential ingredient was to establish the concept of the multiplicity of a focal point. But there was still a rather restrictive condition needed, which is eliminated here by using the concept of finite eigenvalues (or zeros) from the theory of matrix pencils.     

Asymptotic and spectral analysis of non‐selfadjoint operators generated by a filament model with a critical value of a boundary parameter

We consider a class of non‐selfadjoint operators generated by the equation and the boundary conditions, which govern small vibrations of an ideal filament with non‐conservative boundary conditions at one end and a heavy load at the other end. The filament has a non‐constant density and is subject to a viscous damping with a non‐constant damping coefficient. The boundary conditions contain two arbitrary complex parameters. In our previous paper (Mathematical Methods in the Applied Sciences 2001; 24 (15) : 1139–1169), we have derived the asymptotic approximations for the eigenvalues and eigenfunctions of the aforementioned non‐selfadjoint operators when the boundary parameters were arbitrary complex numbers except for one specific value of one of the parameters. We call this value the critical value of the boundary parameter. It has been shown (in Mathematical Methods in the Applied Sciences 2001; 24 (15) : 1139–1169) that the entire set of the eigenvalues is located in a strip parallel to the real axis. The latter property is crucial for the proof of the fact that the set of the root vectors of the operator forms a Riesz basis in the state space of the system. In the present paper, we derive the asymptotics of the spectrum exactly in the case of the critical value of the boundary parameter. We show that in this case, the asymptotics of the eigenvalues is totally different, i.e. both the imaginary and real parts of eigenvalues tend to ∞as the number of an eigenvalue increases. We will show in our next paper, that as an indirect consequence of such a behaviour of the eigenvalues, the set of the root vectors of the corresponding operator is not uniformly minimal (let alone the Riesz basis property). Copyright © 2003 John Wiley & Sons, Ltd.     

Resonance phenomena in a semi-infinite string with a periodic end

We study the propagation of linear waves, generated by a compactly supported time-harmonic force distribution, in a semi-infinite string under the assumption that the material properties depend p-period-ically on the space variable outside a sufficiently large interval [0, a]. The spectrum of the self-adjoint extension A of the spatial part of the differential operator consists of a finite or countable number of bands and a (possibly empty) discrete set of eigenvalues located in the gaps of the continuous spectrum. We show that resonances of order t or t^½, respectively, occur if either ω² is an eigenvalue of A or (i) ω² is a boundary point of the continuous spectrum of A and (ii) the corresponding time-independent homogeneous problem has a non-trivial solution which is p-periodic or p-semiperiodic for x > a (‘standing wave’).