Asymptotics of eigenvalues of the nonlinear eigenvalue problem arising from the near mixed-mode crack-tip stress-strain field problems

In the present paper, approximate analytical and numerical solutions to nonlinear eigenvalue problems arising in nonlinear fracture mechanics in studying stress-strain fields near a crack tip under mixed-mode loading are presented. Asymptotic solutions are obtained by the perturbation method (the artificial small parameter method). The artificial small parameter is the difference between the eigenvalue corresponding to the nonlinear eigenvalue problem and the eigenvalue related to the linear “undisturbed” problem. It is shown that the perturbation technique is an effective method of solving nonlinear eigenvalue problems in nonlinear fracture mechanics. A comparison of numerical and asymptotic results for different values of the mixity parameter and hardening exponent shows good agreement. Thus, the perturbation theory technique for studying nonlinear eigenvalue problems is offered and applied to eigenvalue problems arising in fracture mechanics analysis in the case of mixed-mode loading.     

Approximation and eigenvalue extrapolation of Stokes eigenvalue problem by nonconforming finite element methods

In this paper we analyze the stream function-vorticity-pressure method for the Stokes eigenvalue problem. Further, we obtain full order convergence rate of the eigenvalue approximations for the Stokes eigenvalue problem based on asymptotic error expansions for two nonconforming finite elements, Q ₁^rot and EQ ₁^rot. Using the technique of eigenvalue error expansion, the technique of integral identities and the extrapolation method, we can improve the accuracy of the eigenvalue approximations. This project is supported in part by the National Natural Science Foundation of China (10471103) and is subsidized by the National Basic Research Program of China under the grant 2005CB321701.     

Eigenvalue problem for a coupled channel Schrödinger equation with application to the description of deformed nuclear systems

We discuss the accurate computation of the eigensolutions of systems of coupled channel Schrödinger equations as they appear in studies of real physical phenomena like fission, alpha decay and proton emission. A specific technique is used to compute the solution near the singularity in the origin, while on the rest of the interval the solution is propagated using a piecewise perturbation method. Such a piecewise perturbation method allows us to take large steps even for high energy-values. We consider systems with a deformed potential leading to an eigenvalue problem where the energies are given and the required eigenvalue is related to the adjustment of the potential, viz, the eigenvalue is the depth of the nuclear potential. A shooting technique is presented to determine this eigenvalue accurately.     

A Legendre spectral element method for eigenvalues in hydrodynamic stability

A Legendre polynomial-based spectral technique is developed to be applicable to solving eigenvalue problems which arise in linear and nonlinear stability questions in porous media, and other areas of Continuum Mechanics. The matrices produced in the corresponding generalised eigenvalue problem are sparse, reducing the computational and storage costs, where the superimposition of boundary conditions is not needed due to the structure of the method. Several eigenvalue problems are solved using both the Legendre polynomial-based and Chebyshev tau techniques. In each example, the Legendre polynomial-based spectral technique converges to the required accuracy utilising less polynomials than the Chebyshev tau method, and with much greater computational efficiency.     

Linear response eigenvalue problem solved by extended locally optimal preconditioned conjugate gradient methods

The locally optimal block preconditioned 4-d conjugate gradient method (LOBP4dCG) for the linear response eigenvalue problem was proposed by Bai and Li (2013) and later was extended to the generalized linear response eigenvalue problem by Bai and Li (2014). We put forward two improvements to the method: A shifting deflation technique and an idea of extending the search subspace. The deflation technique is able to deflate away converged eigenpairs from future computation, and the idea of extending the search subspace increases convergence rate per iterative step. The resulting algorithm is called the extended LOBP4dCG (ELOBP4dCG). Numerical results of the ELOBP4dCG strongly demonstrate the capability of deflation technique and effectiveness the search space extension for solving linear response eigenvalue problems arising from linear response analysis of two molecule systems.     

求解陀螺系统特征值问题的收缩二阶Lanczos方法

本文研究陀螺系统特征值问题的数值解法,利用反对称矩阵Lanczos算法,提出了求解陀螺系统特征值问题的二阶Lanczos方法.基于提出的陀螺系统特征值问题的非等价低秩收缩技术,给出了计算陀螺系统极端特征值的收缩二阶Lanczos方法.数值结果说明了算法的有效性.     

A Two-Level Method for Nonsymmetric Eigenvalue Problems

A two-level discretization method for eigenvalue problems is studied. Compared to the standard Galerkin finite element discretization technique performed on a fine grid this method discretizes the eigenvalue problem on a coarse grid and obtains an improved eigenvector (eigenvalue) approximation by solving only a linear problem on the fine grid (or two linear problems for the case of eigenvalue approximation of nonsymmetric problems). The improved solution has the asymptotic accuracy of the Galerkin discretization solution. The link between the method and the iterated Galerkin method is established. Error estimates for the general nonsymmetric case are derived.     

The Unsolvability of Multiplicative Inverse Eigenvalue Problems Almost Everywhere

The idea and technique used in [7] are applied to the multiplicative inverse eigenvalue problems as well. Some sufficient and necessary conditions that the multiplicative inverse eigenvalue problems be unsolvable almost everywhere are given. The results are similar to those of [7], but the proofs are more complicated.     

Logarithmic potential <Emphasis Type="Italic">β</Emphasis>-ensembles and Feynman graphs

We present a diagrammatic technique for calculating the free energy of the matrix eigenvalue model (the model with an arbitrary power of the Vandermonde determinant) to all orders of the 1/N expansion in the case when the limiting eigenvalue distribution spans an arbitrary (but fixed) number of disjoint intervals (curves) and when logarithmic terms are present. This diagrammatic technique is corrected and refined as compared to our first paper with B. Eynard of 2006.     

Perturbation method for solving the nonlinear eigenvalue problem arising from fatigue crack growth problem in a damaged medium

An analytical solution of the nonlinear eigenvalue problem arising from the fatigue crack growth problem in a damaged medium in coupled formulation is obtained. The perturbation technique for solving the nonlinear eigenvalue problem is used. The method allows to find the analytical formula expressing the eigenvalue as the function of parameters of the damage evolution law. It is shown that the eigenvalues of the nonlinear eigenvalue problem are fully determined by the exponents of the damage evolution law. In the paper the third-order (four-term) asymptotic expansions of the angular functions determining the stress and continuity fields in the neighborhood of the crack tip are given. The asymptotic expansions of the angular functions permit to find the closed-form solution for the problem considered.     

Approximation and eigenvalue extrapolation of biharmonic eigenvalue problem by nonconforming finite element methods

In this paper, we analyze the biharmonic eigenvalue problem by two nonconforming finite elements, Q and E Q. We obtain full order convergence rate of the eigenvalue approximations for the biharmonic eigenvalue problem based on asymptotic error expansions for these two nonconforming finite elements. Using the technique of eigenvalue error expansion, the technique of integral identities, and the extrapolation method, we can improve the accuracy of the eigenvalue approximations. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Locating a nearest matrix with an eigenvalue of prespecified algebraic multiplicity

The Wilkinson distance of a matrix A is the two-norm of the smallest perturbation E so that A + E has a multiple eigenvalue. Malyshev derived a singular value optimization characterization for the Wilkinson distance. In this work we generalize the definition of the Wilkinson distance as the two-norm of the smallest perturbation so that the perturbed matrix has an eigenvalue of prespecified algebraic multiplicity. We provide a singular value characterization for this generalized Wilkinson distance. Then we outline a numerical technique to solve the derived singular value optimization problems. In particular the numerical technique is applicable to Malyshev’s formula to compute the Wilkinson distance as well as to retrieve a nearest matrix with a multiple eigenvalue.     

Spectral Galerkin approximation and rigorous error analysis for the Steklov eigenvalue problem in circular domain

In this paper, we present spectral Galerkin approximation and rigorous error analysis for the Steklov eigenvalue problem in a circular domain. First of all, we use the polar coordinate transformation and technique of separation of variables to reduce the problem to a sequence of equivalent 1‐dimensional eigenvalue problems that can be solved individually in parallel. Then, we derive the pole conditions and introduce weighted Sobolev space according to pole conditions. Together with the approximate properties of orthogonal polynomials, we prove the error estimates of approximate eigenvalues for each 1‐dimensional eigenvalue problem. Finally, we provide some numerical experiments to validate the theoretical results and algorithms.     

THE SHIFTED-INVERSE POWER WEAK GALERKIN METHOD FOR EIGENVALUE PROBLEMS

This paper proposes and analyzes a new weak Galerkin method for the eigenvalue problem by using the shifted-inverse power technique. A high order lower bound can be obtained at a relatively low cost via the proposed method. The error estimates for both eigenvalue and eigenfunction are provided and asymptotic lower bounds are shown as well under some conditions. Numerical examples are presented to validate the theoretical analysis.     

The groupoid algebra of an eigenvalue pattern

The eigenvalue pattern of a *-homomorphism between two matrix algebras over commutative C-algebras is a generalization of the Gelfand map in the commutative case. We give a systematic formulation of abstract eigenvalue pattern and extend the classical results by using a technique involving the groupoid algebras of eigenvalue patterns. In the case with matrix algebras over the one-dimensional circle, we characterize all the *-homomorphisms up to unitary equivalence by their eigenvalue patterns. Moreover, this technique has an application to recent classification theorems of C-algebras proved by the present author.

     

Topological multiplicity of the maximum eigenvalue of graph 1-Laplacian

In this paper, we reveal the connection between the independent number of a graph and the topological multiplicity of the maximal eigenvalue of the corresponding graph 1-Laplacian. The pseudo independent number of a graph is introduced, which provides a better lower estimate of the topological multiplicity of the maximum eigenvalue. The technique of our proof is based on the localization property of the eigenvector for graph 1-Laplacian, the Krasnoselski genus, and its relation to the topological join.     

TWO-GRID DISCRETIZATION SCHEMES OF THE NONCONFORMING FEM FOR EIGENVALUE PROBLEMS

This paper extends the two-grid discretization scheme of the conforming finite elements proposed by Xu and Zhou （Math. Comput., 70 （2001）, pp.17-25） to the nonconforming finite elements for eigenvalue problems. In particular, two two-grid discretization schemes based on Rayleigh quotient technique are proposed. By using these new schemes, the solution of an eigenvalue problem on a fine mesh is reduced to that on a much coarser mesh together with the solution of a linear algebraic system on the fine mesh. The resulting solution still maintains an asymptotically optimal accuracy. Comparing with the two-grid discretization scheme of the conforming finite elements, the main advantages of our new schemes are twofold when the mesh size is small enough. First, the lower bounds of the exact eigenvalues in our two-grid discretization schemes can be obtained. Second, the first eigenvalue given by the new schemes has much better accuracy than that obtained by solving the eigenvalue problems on the fine mesh directly.     

A dynamical method of DAEs for the smallest eigenvalue problem

This article gives a new method based on the dynamical system of differential-algebraic equations for the smallest eigenvalue problem of a symmetric matrix. First, the smallest eigenvalue problem is converted into an equivalent constrained optimization problem. Second, from the Karush–Kuhn–Tucker conditions for this special equality-constrained problem, a special continuous dynamical system of differential-algebraic equations is obtained. Lastly, based on the implicit Euler method and an analogous trust-region technique, we obtain a prediction-correction method to compute a steady-state solution of this special system of differential-algebraic equations, and consequently obtain the smallest eigenvalue of the original problem. We also analyze the local superlinear property for this new method, and present the promising numerical results, in comparison with other methods.     

隐式重启精化全局Lanczos方法

The numerical methods for solving large symmetric eigenvalue problems are considered in this paper.Based on the global Lanczos process,a global Lanczos method for solving large symmetric eigenvalue problems is presented.In order to accelerate the convergence of the F-Ritz vectors,the refined global Lanczos method is developed.Combining the implicitly restarted strategy with the deflation technique,an implicitly restarted and refined global Lanczos method for computing some eigenvalues of large symmetric matrices is proposed.Numerical results show that the proposed methods are efficient.     

A multilevel correction method for Stokes eigenvalue problems and its applications

In this paper, a new multilevel correction scheme is proposed to solve Stokes eigenvalue problems by the finite element method. This new scheme contains a series of correction steps, and the accuracy of eigenpair approximation can be improved after each step. In each correction step, we only need to solve a Stokes problem on the corresponding fine finite element space and a Stokes eigenvalue problem on the coarsest finite element space. This correction scheme can improve the efficiency of solving Stokes eigenvalue problems by the finite element method. As applications of this multilevel correction method, a multigrid method and an adaptive finite element technique are introduced for Stokes eigenvalue problems. Some numerical results are given to validate our schemes. Copyright © 2013 John Wiley & Sons, Ltd.