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 solve...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.展开更多
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 th...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.展开更多
The dynamic analysis of damped structural system by using finite element method leads to nonlinear eigenvalue problem(NEP)(particularly,quadratic eigenvalue problem).In general,the parameters of NEP are considered as ...The dynamic analysis of damped structural system by using finite element method leads to nonlinear eigenvalue problem(NEP)(particularly,quadratic eigenvalue problem).In general,the parameters of NEP are considered as exact values.But in actual practice because of different errors and incomplete information,the parameters may have uncertain or vague values and such uncertain values may be considered in terms of fuzzy numbers.This article proposes an efficient fuzzy-affine approach to solve fully fuzzy nonlinear eigenvalue problems(FNEPs)where involved parameters are fuzzy numbers viz.triangular and trapezoidal.Based on the parametric form,fuzzy numbers have been transformed into family of standard intervals.Further due to the presence of interval overestimation problem in standard interval arithmetic,affine arithmetic based approach has been implemented.In the proposed method,the FNEP has been linearized into a generalized eigenvalue problem and further solved by using the fuzzy-affine approach.Several application problems of structures and also general NEPs with fuzzy parameters are investigated based on the proposed procedure.Lastly,fuzzy eigenvalue bounds are illustrated with fuzzy plots with respect to its membership function.Few comparisons are also demonstrated to show the reliability and efficacy of the present approach.展开更多
We study the smooth LU decomposition of a given analytic functional A-matrix A(A) and its block-analogue. Sufficient conditions for the existence of such matrix decompositions are given, some differentiability about...We study the smooth LU decomposition of a given analytic functional A-matrix A(A) and its block-analogue. Sufficient conditions for the existence of such matrix decompositions are given, some differentiability about certain elements arising from them are proved, and several explicit expressions for derivatives of the specified elements are provided. By using these smooth LU decompositions, we propose two numerical methods for computing multiple nonlinear eigenvalues of A(A), and establish their locally quadratic convergence properties. Several numerical examples are provided to show the feasibility and effectiveness of these new methods.展开更多
In this paper we consider the existence of infinitely many positive solutions for second order nonlinear eigenvalue problem with singular coefficient function. By the use of Krasnosel'skii fixed point theorem of c...In this paper we consider the existence of infinitely many positive solutions for second order nonlinear eigenvalue problem with singular coefficient function. By the use of Krasnosel'skii fixed point theorem of cone expansion-compression type we give several sufficient conditions.展开更多
Nonlinear rank-one modification of the symmetric eigenvalue problem arises from eigenvibrations of mechanical structures with elastically attached loads and calculation of the propagation modes in optical fiber. In th...Nonlinear rank-one modification of the symmetric eigenvalue problem arises from eigenvibrations of mechanical structures with elastically attached loads and calculation of the propagation modes in optical fiber. In this paper, we first study the existence and uniqueness of eigenvalues, and then investigate three numerical algorithms, namely Picard iteration, nonlinear Rayleigh quotient iteration and successive linear approximation method (SLAM). The global convergence of the SLAM is proven under some mild assumptions. Numerical examples illustrate that the SLAM is the most robust method.展开更多
The transmission eigenvalue problem is an eigenvalue problem that arises in the scatter- ing of time-harmonic waves by an inhomogeneous medium of compact support. Based on a fourth order formulation, the transmission ...The transmission eigenvalue problem is an eigenvalue problem that arises in the scatter- ing of time-harmonic waves by an inhomogeneous medium of compact support. Based on a fourth order formulation, the transmission eigenvalue problem is discretized by the Mor- ley element. For the resulting quadratic eigenvalue problem, a recursive integral method is used to compute real and complex eigenvalues in prescribed regions in the complex plane. Numerical examples are presented to demonstrate the effectiveness of the proposed method.展开更多
In this paper we propose a uniformly convergent numerical method for discretizing singularly perturbed nonlinear eigenvalue problems under constraints with applications in Bose-Einstein condensation and quantum chemis...In this paper we propose a uniformly convergent numerical method for discretizing singularly perturbed nonlinear eigenvalue problems under constraints with applications in Bose-Einstein condensation and quantum chemistry.We begin with the time-independent Gross-Pitaevskii equation and show how to reformulate it into a singularly perturbed nonlinear eigenvalue problem under a constraint.Matched asymptotic approximations for the problem are presented to locate the positions and characterize the widths of boundary layers and/or interior layers in the solution.A uniformly convergent numerical method is proposed by using the normalized gradient flow and piecewise uniform mesh techniques based on the asymptotic approximations for the problem.Extensive numerical results are reported to demonstrate the effectiveness of our numerical method for the problems.Finally,the method is applied to compute ground and excited states of Bose-Einstein condensation in the semiclassical regime and some conclusive findings are reported.展开更多
In this paper, we investigate a fourth-order three-point nonlinear eigenvalue problem. By the Krasnosel’skii’s fixed-point theorem in cone, some new results on the existence of positive solutions are obtained.
In this paper we study the existence of nontrivial solutions of a class of asymptotically linear elliptic resonant problems at higher eigenvalues with the nonlinear term which may be un- bounded by making use of the M...In this paper we study the existence of nontrivial solutions of a class of asymptotically linear elliptic resonant problems at higher eigenvalues with the nonlinear term which may be un- bounded by making use of the Morse theory for a C^2-function at both isolated critical point and infinity.展开更多
The focus of this paper is on determination of the dynamic parameters of structural systems with viscoelastic (VE) dampers described by Maxwell rheological models. Such parameters could be obtained after solving the a...The focus of this paper is on determination of the dynamic parameters of structural systems with viscoelastic (VE) dampers described by Maxwell rheological models. Such parameters could be obtained after solving the appropriately defined nonlinear eigenvalue problem for frames with VE dampers. The solution to the nonlinear eigenvalue problem is obtained by equating to zero the determinant of the considered system of equations. Apart from complex conjugate eigenvalues, the real ones occurred when dampers that are described by the classic Maxwell model, are also determined.展开更多
We propose a finite element method to compute the band structures of dispersive photonic crystals in 3D.The nonlinear Maxwell’s eigenvalue problem is formulated as the eigenvalue problem of a holomorphic operator fun...We propose a finite element method to compute the band structures of dispersive photonic crystals in 3D.The nonlinear Maxwell’s eigenvalue problem is formulated as the eigenvalue problem of a holomorphic operator function.The N´ed´elec edge elements are employed to discretize the operators,where the divergence free condition for the electric field is realized by a mixed form using a Lagrange multiplier.The convergence of the eigenvalues is proved using the abstract approximation theory for holomorphic operator functions with the regular approximation of the edge elements.The spectral indicator method is then applied to compute the discrete eigenvalues.Numerical examples are presented demonstrating the effectiveness of the proposed method.展开更多
In this work the one-band effective Hamiltonian governing the electronic states of a quantum dot/ring in a homogenous magnetic field is used to derive a pair/quadruple of nonlinear eigenvalue problems corresponding to...In this work the one-band effective Hamiltonian governing the electronic states of a quantum dot/ring in a homogenous magnetic field is used to derive a pair/quadruple of nonlinear eigenvalue problems corresponding to different spin orientations and in case of rotational symmetry additionally to quantum number±ℓ.We show,that each of those pair/quadruple of nonlinear problems allows for the minmax characterization of its eigenvalues under certain conditions,which are satisfied for our examples and the common InAs/GaAs heterojunction.Exploiting the minmax property we devise efficient iterative projection methods simultaneously handling the pair/quadruple of nonlinear problems and thereby saving up to 40%of the computational time as compared to the nonlinear Arnoldi method applied to each of the problems separately.展开更多
In this paper we apply a cone theoretic fixed-point theorem proved by Krasnosel'skii to obtain sufficient conditions for the existence of positive solutions to some boundary value problems for a class of functional d...In this paper we apply a cone theoretic fixed-point theorem proved by Krasnosel'skii to obtain sufficient conditions for the existence of positive solutions to some boundary value problems for a class of functional difference equations. We consider the case that the nonlinear term satisfies asvrnntoticallv linear growth.展开更多
基金the National Science and Tech-nology Council,Taiwan for their financial support(Grant Number NSTC 111-2221-E-019-048).
文摘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.
基金supported by the National Natural Science Foundation of China(Nos.1133200711202147+2 种基金and 9216111)the Specialized Research Fund for the Doctoral Program of Higher Education(No.20120032120007)the Open Fund from State Key Laboratory of Aerodynamics(Nos.SKLA201201 and SKLA201301)
文摘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.
文摘The dynamic analysis of damped structural system by using finite element method leads to nonlinear eigenvalue problem(NEP)(particularly,quadratic eigenvalue problem).In general,the parameters of NEP are considered as exact values.But in actual practice because of different errors and incomplete information,the parameters may have uncertain or vague values and such uncertain values may be considered in terms of fuzzy numbers.This article proposes an efficient fuzzy-affine approach to solve fully fuzzy nonlinear eigenvalue problems(FNEPs)where involved parameters are fuzzy numbers viz.triangular and trapezoidal.Based on the parametric form,fuzzy numbers have been transformed into family of standard intervals.Further due to the presence of interval overestimation problem in standard interval arithmetic,affine arithmetic based approach has been implemented.In the proposed method,the FNEP has been linearized into a generalized eigenvalue problem and further solved by using the fuzzy-affine approach.Several application problems of structures and also general NEPs with fuzzy parameters are investigated based on the proposed procedure.Lastly,fuzzy eigenvalue bounds are illustrated with fuzzy plots with respect to its membership function.Few comparisons are also demonstrated to show the reliability and efficacy of the present approach.
基金supported by the National Basic Research Program(No.2005CB321702)the China Outstanding Young Scientist F0undation(No.10525102)the National Natural Science Foundation (No.10471146),P.R.China
文摘We study the smooth LU decomposition of a given analytic functional A-matrix A(A) and its block-analogue. Sufficient conditions for the existence of such matrix decompositions are given, some differentiability about certain elements arising from them are proved, and several explicit expressions for derivatives of the specified elements are provided. By using these smooth LU decompositions, we propose two numerical methods for computing multiple nonlinear eigenvalues of A(A), and establish their locally quadratic convergence properties. Several numerical examples are provided to show the feasibility and effectiveness of these new methods.
文摘In this paper we consider the existence of infinitely many positive solutions for second order nonlinear eigenvalue problem with singular coefficient function. By the use of Krasnosel'skii fixed point theorem of cone expansion-compression type we give several sufficient conditions.
基金supported in part by NSF grants DMS-0611548 and OCI-0749217 and DOE grant DE-FC02-06ER25794supported in part by NSF of China under the contract number 10871049 and Shanghai Down project 200601.
文摘Nonlinear rank-one modification of the symmetric eigenvalue problem arises from eigenvibrations of mechanical structures with elastically attached loads and calculation of the propagation modes in optical fiber. In this paper, we first study the existence and uniqueness of eigenvalues, and then investigate three numerical algorithms, namely Picard iteration, nonlinear Rayleigh quotient iteration and successive linear approximation method (SLAM). The global convergence of the SLAM is proven under some mild assumptions. Numerical examples illustrate that the SLAM is the most robust method.
文摘The transmission eigenvalue problem is an eigenvalue problem that arises in the scatter- ing of time-harmonic waves by an inhomogeneous medium of compact support. Based on a fourth order formulation, the transmission eigenvalue problem is discretized by the Mor- ley element. For the resulting quadratic eigenvalue problem, a recursive integral method is used to compute real and complex eigenvalues in prescribed regions in the complex plane. Numerical examples are presented to demonstrate the effectiveness of the proposed method.
基金Singapore Ministry of Education grant No.R-146-000-083-112 and would like to thank Professor Tao Tang for very helpful discussion on the subject.
文摘In this paper we propose a uniformly convergent numerical method for discretizing singularly perturbed nonlinear eigenvalue problems under constraints with applications in Bose-Einstein condensation and quantum chemistry.We begin with the time-independent Gross-Pitaevskii equation and show how to reformulate it into a singularly perturbed nonlinear eigenvalue problem under a constraint.Matched asymptotic approximations for the problem are presented to locate the positions and characterize the widths of boundary layers and/or interior layers in the solution.A uniformly convergent numerical method is proposed by using the normalized gradient flow and piecewise uniform mesh techniques based on the asymptotic approximations for the problem.Extensive numerical results are reported to demonstrate the effectiveness of our numerical method for the problems.Finally,the method is applied to compute ground and excited states of Bose-Einstein condensation in the semiclassical regime and some conclusive findings are reported.
基金the National Natural Science Foundation of China(11261053)The Natural Science Foundation of Gansu Province(1308RJZA125)
文摘In this paper, we investigate a fourth-order three-point nonlinear eigenvalue problem. By the Krasnosel’skii’s fixed-point theorem in cone, some new results on the existence of positive solutions are obtained.
文摘In this paper we study the existence of nontrivial solutions of a class of asymptotically linear elliptic resonant problems at higher eigenvalues with the nonlinear term which may be un- bounded by making use of the Morse theory for a C^2-function at both isolated critical point and infinity.
基金the financial support received from the Poznan University of Technology(Grant No.DS 11-088/12)in connection with this work.
文摘The focus of this paper is on determination of the dynamic parameters of structural systems with viscoelastic (VE) dampers described by Maxwell rheological models. Such parameters could be obtained after solving the appropriately defined nonlinear eigenvalue problem for frames with VE dampers. The solution to the nonlinear eigenvalue problem is obtained by equating to zero the determinant of the considered system of equations. Apart from complex conjugate eigenvalues, the real ones occurred when dampers that are described by the classic Maxwell model, are also determined.
基金China Postdoctoral Science Foundation Grant 2019M650460the NSF grant DMS-2011148.The research of J.Sun is supported partially by the Simons Foundation Grant 711922.
文摘We propose a finite element method to compute the band structures of dispersive photonic crystals in 3D.The nonlinear Maxwell’s eigenvalue problem is formulated as the eigenvalue problem of a holomorphic operator function.The N´ed´elec edge elements are employed to discretize the operators,where the divergence free condition for the electric field is realized by a mixed form using a Lagrange multiplier.The convergence of the eigenvalues is proved using the abstract approximation theory for holomorphic operator functions with the regular approximation of the edge elements.The spectral indicator method is then applied to compute the discrete eigenvalues.Numerical examples are presented demonstrating the effectiveness of the proposed method.
基金We would like to thank Oleksandr Voskoboynikov for his comments on the physical relevance of the model under consideration.We also thank the anonymous referees for their comments helping us to improve this manuscript.
文摘In this work the one-band effective Hamiltonian governing the electronic states of a quantum dot/ring in a homogenous magnetic field is used to derive a pair/quadruple of nonlinear eigenvalue problems corresponding to different spin orientations and in case of rotational symmetry additionally to quantum number±ℓ.We show,that each of those pair/quadruple of nonlinear problems allows for the minmax characterization of its eigenvalues under certain conditions,which are satisfied for our examples and the common InAs/GaAs heterojunction.Exploiting the minmax property we devise efficient iterative projection methods simultaneously handling the pair/quadruple of nonlinear problems and thereby saving up to 40%of the computational time as compared to the nonlinear Arnoldi method applied to each of the problems separately.
文摘In this paper we apply a cone theoretic fixed-point theorem proved by Krasnosel'skii to obtain sufficient conditions for the existence of positive solutions to some boundary value problems for a class of functional difference equations. We consider the case that the nonlinear term satisfies asvrnntoticallv linear growth.