Highly Accurate Golden Section Search Algorithms and Fictitious Time Integration Method for Solving Nonlinear Eigenvalue Problems

This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.

Some Trace Formulas of a Eigenvalue Problem

In present paper, using some methods of approximation theory, the trace formulas for eigenvalues of a eigenvalue problem are calculated under the periodic condition and the decaying condition at x∞.

One Computational Method of the Eigenvalues of the Horizontal Vibration Problem of Beam

This paper considers one computational method of the eigenvalues approximate value of the horizontal vibration problem of beam. The proof of our main result is based on the variational formula. First of all, Cauchy inequality is used to obtain a basic inequality. Secondly, the functions of basis are made by Galerkin method, and the error estimates of eignevalues are obtained by Cauchy inequality. At last, the computational method of the approximate value of the eigenvalues turns out immediately, and acc...

ASYMPTOTIC BEHAVIOUR OF EIGENVALUES FOR THE DISCONTINUOUS BOUNDARY-VALUE PROBLEM WITH FUNCTIONAL-TRANSMISSION CONDITIONS

In this study, the boundary-value problem with eigenvalue parameter generated by the differential equation with discontinuous coefficients and boundary conditions which contains not only endpoints of the considered interval, but also point of discontinuity and linear functionals is investigated. So, the problem is not pure boundary-value. The authors single out a class of linear functionals and find simple algebraic conditions on coefficients, which garantee the existence of infinit number eigenvalues. Also the asymptotic formulas for eigenvalues are found.

THE EIGENVALUE PROBLEM FOR THE LAPLACIAN EQUATIONS

This article studies the Dirichlet eigenvalue problem for the Laplacian equations △u = -λu, x ∈Ω , u = 0, x ∈δΩ, where Ω belong to R^n is a smooth bounded convex domain. By using the method of appropriate barrier function combined with the maximum principle, authors obtain a sharp lower bound of the difference of the first two eigenvalues for the Dirichlet eigenvalue problem. This study improves the result of S.T. Yau et al.

Mode decomposition of nonlinear eigenvalue problems and application in flow stability

Direct numerical simulations are carried out with different disturbance forms introduced into the inlet of a flat plate boundary layer with the Mach number 4.5. According to the biorthogonal eigenfunction system of the linearized Navier-Stokes equations and the adjoint equations, the decomposition of the direct numerical simulation results into the discrete normal mode is easily realized. The decomposition coefficients can be solved by doing the inner product between the numerical results and the eigenfunctions of the adjoint equations. For the quadratic polynomial eigenvalue problem, the inner product operator is given in a simple form, and it is extended to an Nth-degree polynomial eigenvalue problem. The examples illustrate that the simplified mode decomposition is available to analyze direct numerical simulation results.

Two-level stabilized finite element method for Stokes eigenvalue problem

A two-level stabilized finite element method for the Stokes eigenvalue problem based on the local Gauss integration is considered. This method involves solving a Stokes eigenvalue problem on a coarse mesh with mesh size H and a Stokes problem on a fine mesh with mesh size h -- O（H2）, which can still maintain the asymptotically optimal accuracy. It provides an approximate solution with the convergence rate of the same order as the usual stabilized finite element solution, which involves solving a Stokes eigenvalue problem on a fine mesh with mesh size h. Hence, the two-level stabilized finite element method can save a large amount of computational time. Moreover, numerical tests confirm the theoretical results of the present method.

Eigenvalue Problem of Doubly Stochastic Hamiltonian Systems with Boundary Conditions

In this paper, we investigate the eigenvalue problem of forward-backward doubly stochastic dii^erential equations with boundary value conditions. We show that this problem can be represented as an eigenvalue problem of a bounded continuous compact operator. Hence using the famous Hilbert-Schmidt spectrum theory, we can characterize the eigenvalues exactly.

Solvability of Inverse Eigenvalue Problem for Dense Singular Symmetric Matrices

Given a list of real numbers ∧={λ1,…, λn}, we determine the conditions under which ∧will form the spectrum of a dense n × n singular symmetric matrix. Based on a solvability lemma, an algorithm to compute the elements of the matrix is derived for a given list ∧ and dependency parameters. Explicit computations are performed for n≤5 and r≤4 to illustrate the result.

THE TWO-LEVEL STABILIZED FINITE ELEMENT METHOD BASED ON MULTISCALE ENRICHMENT FOR THE STOKES EIGENVALUE PROBLEM

In this paper,we first propose a new stabilized finite element method for the Stokes eigenvalue problem.This new method is based on multiscale enrichment,and is derived from the Stokes eigenvalue problem itself.The convergence of this new stabilized method is proved and the optimal priori error estimates for the eigenfunctions and eigenvalues are also obtained.Moreover,we combine this new stabilized finite element method with the two-level method to give a new two-level stabilized finite element method for the Stokes eigenvalue problem.Furthermore,we have proved a priori error estimates for this new two-level stabilized method.Finally,numerical examples confirm our theoretical analysis and validate the high effectiveness of the new methods.

Physics-constrained neural network for solving discontinuous interface K-eigenvalue problem with application to reactor physics

Machine learning-based modeling of reactor physics problems has attracted increasing interest in recent years.Despite some progress in one-dimensional problems,there is still a paucity of benchmark studies that are easy to solve using traditional numerical methods albeit still challenging using neural networks for a wide range of practical problems.We present two networks,namely the Generalized Inverse Power Method Neural Network(GIPMNN)and Physics-Constrained GIPMNN(PC-GIPIMNN)to solve K-eigenvalue problems in neutron diffusion theory.GIPMNN follows the main idea of the inverse power method and determines the lowest eigenvalue using an iterative method.The PC-GIPMNN additionally enforces conservative interface conditions for the neutron flux.Meanwhile,Deep Ritz Method(DRM)directly solves the smallest eigenvalue by minimizing the eigenvalue in Rayleigh quotient form.A comprehensive study was conducted using GIPMNN,PC-GIPMNN,and DRM to solve problems of complex spatial geometry with variant material domains from the fleld of nuclear reactor physics.The methods were compared with the standard flnite element method.The applicability and accuracy of the methods are reported and indicate that PC-GIPMNN outperforms GIPMNN and DRM.

THE SUFFICIENT CONDITIONS FOR SOLVABILITY OF ALGEBRAIC INVERSE EIGENVALUE PROBLEMS

Applying constructed homotopy and its properties,we gel some sufficient conditions for the solvability of algebraic inverse eigenvalue problems,which are better than that of the paper [4] in some cases. Inverse eigenvalue problems,solvability,sufficient conditions.

THE UNSOLVABILITY OF GENERALIZED INVERSE EIGENVALUE PROBLEMS ALMOST EVERYWHERE

In this paper the unsolvability of generalized inverse eigenvalue problems almost everywhere is discussed.We first give the definitions for the unsolvability of generalized inverse eigenvalue problems almost everywhere.Then adopting the method used in [14],we present some sufficient conditions such that the generalized inverse eigenvalue problems are unsohable almost everywhere.

Uniqueness theorem,theorem of reciprocity,and eigenvalue problems in linear theory of porous piezoelectricity

The uniqueness theorem and the theorem of reciprocity in the linearized porous piezoelectricity are established under the assumption of positive definiteness of elastic and electric fields. General theorems in the linear theory of porous piezoelectric materials are proved for the quasi-static electric field approximation. The uniqueness theorem is also proved without using the positive definiteness of the elastic field. An eigenvalue problem associated with free vibrations of a porous piezoelectric body is stud- ied using the abstract formulation. Some properties of operators are also proved. The problem of frequency shift due to small disturbances, based on an abstract formulation, is studied using a variational and operator approach. A perturbation analysis of a special ease is also given.

A BIHARMONIC EIGENVALUE PROBLEM AND ITS APPLICATION

In this article, we focus on the eigenvalue problem of the following linear biharmonic equation in R^N：△^2u-αu＋λg（x）u=0 with u ∈H^2（R^N）,u≠0,N≥5Note that there are two parameters α and λ in it, which is different from the usual eigenvalue problems. Here, we consider λ as an eigenvalue and seek Ior a sulble range of parameter α, which ensures that problem （＊） has a maximal eigenvalue. As the loss of strong maximum principle for our problem, we can only get the existence of non-trivial solutions, not positive solutions, in this article. As an application, by using these results, we studied also the existence of non-trivial solutions for an asymptotically linear biharmonic equation in R^N.

THE BI-SELF-CONJUGATE AND NONNEGATIVE DEFINITE SOLUTIONS TO THE INVERSE EIGENVALUE PROBLEM OF QUATERNION MATRICES

The main aim of this paper is to discuss the following two problems: Problem I: Given X ∈ Hn×m (the set of all n×m quaternion matrices), A = diag(λ1,…, λm) EEEEE Hm×m, find A ∈ BSHn×n≥such that AX = X(?), where BSHn×n≥ denotes the set of all n×n quaternion matrices which are bi-self-conjugate and nonnegative definite. Problem Ⅱ2= Given B ∈ Hn×m, find B ∈ SE such that ||B-B||Q = minAE∈=sE ||B-A||Q, where SE is the solution set of problem I , || ·||Q is the quaternion matrix norm. The necessary and sufficient conditions for SE being nonempty are obtained. The general form of elements in SE and the expression of the unique solution B of problem Ⅱ are given.

NETWORK ANALYSIS OF EIGENVALUE PROBLEM FOR MULTILAYER DIELECTRIC WAVEGUIDE CONSISTING OF ARBITRARY NUMBER OF LAYERS

In this paper,the eigenvalue problem of a multilayer dielectric waveguide consisting of arbitrarynumber of layers is solved by the microwave network method.A general program with the function of com-puter graphics has been developed for analyzing the dispersion characteristics and the electromagnetic fielddistributions of an N layer dielectric waveguide.As an example of practical applications,the procedure ofmode conversion and mode separation in dielectric branching waveguides is vividly demonstrated throughanalyzing the field distributions of asymmetric multilayer dielectric structures and the general rules of modeconversion are discussed.

A NEW MATRIX PERTURBATION METHOD FOR ANALYTICAL SOLUTION OF THE COMPLEX MODAL EIGENVALUE PROBLEM OF VISCOUSLY DAMPED LINEAR VIBRATION SYSTEMS

A new matrix perturbation analysis method is presented for efficient approximate solution of the complex modal quadratic generalized eigenvalue problem of viscously damped linear vibration systems. First, the damping matrix is decomposed into the sum of a proportional-and a nonproportional-damping parts, and the solutions of the real modal eigenproblem with the proportional dampings are determined, which are a set of initial approximate solutions of the complex modal eigenproblem. Second, by taking the nonproportional-damping part as a small modification to the proportional one and using the matrix perturbation analysis method, a set of approximate solutions of the complex modal eigenvalue problem can be obtained analytically. The result is quite simple. The new method is applicable to the systems with viscous dampings-which do not deviate far away from the proportional-damping case. It is particularly important that the solution technique be also effective to the systems with heavy, but not over, dampings. The solution formulas of complex modal eigenvlaues and eigenvectors are derived up to second-order perturbation terms. The effectiveness of the perturbation algorithm is illustrated by an exemplar numerical problem with heavy dampings. In addition, the practicability of approximately estimating the complex modal eigenvalues, under the proportional-damping hypothesis, of damped vibration systems is discussed by several numerical examples.

A new algorithm for an inverse eigenvalue problem on Jacobi matrices

In this paper, an inverse problem on Jacobi matrices presented by Shieh in 2004 is studied. Shieh＇s result is improved and a new and stable algorithm to reconstruct its solution is given. The numerical examples is also given.

Multiple Solutions for the Eigenvalue Problem of Nonlinear Fractional Differential Equations

In this article, multiple solutions for the eigenvalue problem of nonlinear fractional differential equation is considered. We obtain the existence and multiplicity results of positive solutions by using some fixed point theorems.