Abstract:In this article we are interested for the numerical computation of spectra of non-self adjoint quadratic operators, in two and three spatial dimensions. Indeed, in the multidimensional case very few results are known on the location of the eignevalues. This leads to solve nonlinear eigenvalue problems. In introduction we begin with a review of theoretical results and numerical results obtained for the one dimensional case. Then we present the numerical methods developed to compute the spectra (finite difference discretization) for the two and three dimensional cases. The numerical results obtained are presented and analyzed. One difficulty here is that we have to compute eigenvalues of strongly non-self-adjoint operators which are unstable. This work is in continuity of a previous work in one spatial dimension.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to calculate the spectra of non - self - adjoint quadratic operators in two - dimensional and three - dimensional spaces. Specifically, the paper focuses on the eigenvalue location problem of non - self - adjoint quadratic operators in the multi - dimensional case (two - dimensional and three - dimensional). Since theoretical results are very limited in these high - dimensional cases, numerical methods need to be developed to solve nonlinear eigenvalue problems.
### Main Research Contents
1. **Background and Motivation**:
- In one - dimensional space, there have been some theoretical and numerical results regarding the spectra of non - self - adjoint operators.
- But in two - dimensional and three - dimensional spaces, very little is known about the eigenvalue location of non - self - adjoint quadratic operators.
- This situation has led to the need for nonlinear eigenvalue problems, that is, to solve problems of the form \( L(\lambda)u(x)=0 \), where \( L(\lambda) \) is a family of operators depending on the complex parameter \( \lambda \).
2. **Research Objectives**:
- The goal of the paper is to extend previous work in one - dimensional space and study the spectra of non - self - adjoint quadratic operators in higher - dimensional (two - dimensional and three - dimensional) spaces.
- In particular, study the operator \( L_P(\lambda)=\Delta+(P(x)\lambda)^2 \), where \( P(x) \) is an elliptic positive polynomial, defined on the Hilbert space \( L^2(\mathbb{R}^n) \).
3. **Numerical Methods**:
- Use the finite - difference method to discretize the problem in order to solve cases with homogeneous Dirichlet boundary conditions and periodic boundary conditions.
- For homogeneous Dirichlet boundary conditions, use five - point and seven - point finite - difference schemes to handle two - dimensional and three - dimensional problems respectively.
- For periodic boundary conditions, consider cases of different coefficients \( a_j(x_j) \), including constant coefficients and sine function coefficients.
4. **Numerical Results and Discussion**:
- Present and analyze the numerical results obtained in two - dimensional and three - dimensional spaces.
- The results show that under certain conditions, the distribution of eigenvalues is consistent with the theoretical results in one - dimensional space.
- Discuss the influence of the parameter \( c \), the grid resolution \( N \) and the size of the computational domain \( L \) on the eigenvalue distribution.
### Conclusions
- The paper successfully calculates the spectra of non - self - adjoint quadratic operators in two - dimensional and three - dimensional spaces by numerical methods.
- The numerical results show that the theoretical results in one - dimensional space still hold in two - dimensional and three - dimensional spaces.
- Propose directions for future research, especially in terms of improving computational efficiency and accuracy.
### Formula Summary
The formulas involved in the paper mainly include:
- Operator \( L(\lambda)=H_0+\lambda H_1+\lambda^2I \)
- Laplace operator \( \Delta u(x,y)=\frac{\partial^2 u}{\partial x^2}+\frac{\partial^2 u}{\partial y^2} \)
- Finite - difference approximation \( \frac{\partial^2 u}{\partial x^2}\approx\frac{u_{i - 1,j}-2u_{i,j}+u_{i + 1,j}}{h^2} \)
These formulas are used to describe and solve nonlinear eigenvalue problems.