On an eigenvalue for the Laplace operator in a disk with Dirichlet boundary condition on a small part of the boundary in a critical case

A boundary-value problem of finding eigenvalues is considered for the negative Laplace operator in a disk with Neumann boundary condition on almost all the circle except for a small arc of vanishing length, where the Dirichlet boundary condition is imposed. A complete asymptotic expansion with respect to a parameter (the length of the small arc) is constructed for an eigenvalue of this problem that converges to a double eigenvalue of the Neumann problem.     

Steklov-Neumann eigenproblems and nonlinear elliptic equations with nonlinear boundary conditions

We study the solvability of nonlinear second order elliptic partial differential equations with nonlinear boundary conditions. We introduce the notion of “eigenvalue-lines” in the plane; these eigenvalue-lines join each Steklov eigenvalue to the first eigenvalue of the Neumann problem with homogeneous boundary condition. We prove existence results when the nonlinearities involved asymptotically stay, in some sense, below the first eigenvalue-lines or in a quadrilateral region (depicted in Fig. 1) enclosed by two consecutive eigenvalue-lines. As a special case we derive the so-called nonresonance results below the first Steklov eigenvalue as well as between two consecutive Steklov eigenvalues. The case in which the eigenvalue-lines join each Neumann eigenvalue to the first Steklov eigenvalue is also considered. Our method of proof is variational and relies mainly on minimax methods in critical point theory.     

Laguerre collocation solutions to boundary layer type problems

In this paper a Laguerre collocation type method based on usual Laguerre functions is designed in order to solve high order nonlinear boundary value problems as well as eigenvalue problems, on semi-infinite domain. The method is first applied to Falkner–Skan boundary value problem. The solution along with its first two derivatives are computed inside the boundary layer on a fine grid which cluster towards the fixed boundary. Then the method is used to solve a generalized eigenvalue problem which arise in the study of the stability of the Ekman boundary layer. The method provides reliable numerical approximations, is robust and easy implementable. It introduces the boundary condition at infinity without any truncation of the domain. A particular attention is payed to the treatment of boundary conditions at origin. The dependence of the set of solutions to Falkner–Skan problem on the parameter embedded in the system is reproduced correctly. For Ekman eigenvalue problem, the critical Reynolds number which assure the linear stability is computed and compared with existing results. The leftmost part of the spectrum is validated using QZ as well as some Jacobi–Davidson type methods.     

Determination of the number of an eigenvalue of a singular nonlinear self-adjoint spectral problem for a linear Hamiltonian system of differential equations

We suggest a method for determining the number of an eigenvalue of a self-adjoint spectral problem nonlinear with respect to the spectral parameter, for some class of Hamiltonian systems of ordinary differential equations on the half-line. The standard boundary conditions are posed at zero, and the solution boundedness condition is posed at infinity. We assume that the matrix of the system is monotone with respect to the spectral parameter. The number of an eigenvalue is determined by the properties of the corresponding nontrivially solvable homogeneous boundary value problem. For the considered class of systems, it becomes possible to compute the numbers of eigenvalues lying in a given range of the spectral parameter without finding the eigenvalues themselves.     

Existence of an extremal ground state energy of a nanostructured quantum dot

This paper is concerned with two rearrangement optimization problems. These problems are motivated by two eigenvalue problems which depend nonlinearly on the eigenvalues. We consider a rational and a quadratic eigenvalue problem with Dirichlet’s boundary condition and investigate two related optimization problems where the goal function is the corresponding first eigenvalue. The first eigenvalue in the rational eigenvalue problem represents the ground state energy of a nanostructured quantum dot. In both the problems, the admissible set is a rearrangement class of a given function.     

On the numerical treatment of the eigenparameter dependent boundary conditions

In this paper, we consider the numerical treatment of singular eigenvalue problems supplied with eigenparameter dependent boundary conditions using spectral methods. On the one hand, such boundary conditions hinder the construction of test and trial space functions which could incorporate them and thus providing well-conditioned Galerkin discretization matrices. On the other hand, they can generate surprising behavior of the eigenvectors hardly detected by analytic methods. These singular problems are often indirectly approximated by regular ones. We argue that spectral collocation as well as tau method offer remedies for the first two issues and provide direct and efficient treatment to such problems. On a finite domain, we consider the so-called Petterson-König’s rod eigenvalue problem and on the half line, we take into account the Charney’s baroclinic stability problem and the Fourier eigenvalue problem. One boundary condition in these problems depends on the eigenparameter and additionally, this also could depend on some physical parameters. The Chebyshev collocation based on both, square and rectangular differentiation and a Chebyshev tau method are used to discretize the first problem. All these schemes cast the problems into singular algebraic generalized eigenvalue ones which are solved by the QZ and/or Arnoldi algorithms as well as by some target oriented Jacobi-Davidson methods. Thus, the spurious eigenvalues are completely eliminated. The accuracy of square Chebyshev collocation is roughly estimated and its order of approximation with respect to the eigenvalue of interest is determined. For the problems defined on the half line, we make use of the Laguerre-Gauss-Radau collocation. The method proved to be reliable, accurate, and stable with respect to the order of approximation and the scaling parameter.     

A boundary element method for the Dirichlet eigenvalue problem of the Laplace operator

The solution of eigenvalue problems for partial differential operators by using boundary integral equation methods usually involves some Newton potentials which may be resolved by using a multiple reciprocity approach. Here we propose an alternative approach which is in some sense equivalent to the above. Instead of a linear eigenvalue problem for the partial differential operator we consider a nonlinear eigenvalue problem for an associated boundary integral operator. This nonlinear eigenvalue problem can be solved by using some appropriate iterative scheme, here we will consider a Newton scheme. We will discuss the convergence and the boundary element discretization of this algorithm, and give some numerical results.     

Propagation of TM waves in a layer with arbitrary nonlinearity

A boundary value problem for Maxwell’s equations describing propagation of TM waves in a nonlinear dielectric layer with arbitrary nonlinearity is considered. The layer is located between two linear semi-infinite media. The problem is reduced to a nonlinear boundary eigenvalue problem for a system of second-order nonlinear ordinary differential equations. A dispersion equation for the eigenvalues of the problem (propagation constants) is derived. For a given nonlinearity function, the dispersion equation can be studied both analytically and numerically. A sufficient condition for the existence of at least one eigenvalue is formulated.     

A new type of full multigrid method for the elasticity eigenvalue problem

This paper introduces a new type of full multigrid method for the elasticity eigenvalue problem. The main idea is to avoid solving large scale elasticity eigenvalue problem directly by transforming the solution of the elasticity eigenvalue problem into a series of solutions of linear boundary value problems defined on a multilevel finite element space sequence and some small scale elasticity eigenvalue problems defined on the coarsest correction space. The involved linear boundary value problems will be solved by performing some multigrid iterations. Besides, some efficient techniques such as parallel computing and adaptive mesh refinement can also be absorbed in our algorithm. The efficiency and validity of the multigrid methods are verified by several numerical experiments.     

Solvability and stability of the inverse Sturm–Liouville problem with analytical functions in the boundary condition

The paper deals with the Sturm–Liouville eigenvalue problem with the Dirichlet boundary condition at one end of the interval and with the boundary condition containing entire functions of the spectral parameter at the other end. We study the inverse problem, which consists in recovering the potential from a part of the spectrum. This inverse problem generalizes partial inverse problems on finite intervals and on graphs and also the inverse transmission eigenvalue problem. We obtain sufficient conditions for global solvability of the studied inverse problem, which prove its local solvability and stability. In addition, application of our main results to the partial inverse Sturm–Liouville problem on the star-shaped graph is provided.     

Principal Eigenvalue Minimization for an Elliptic Problem with Indefinite Weight and Robin Boundary Conditions

This paper focuses on the study of a linear eigenvalue problem with indefinite weight and Robin type boundary conditions. We investigate the minimization of the positive principal eigenvalue under the constraint that the absolute value of the weight is bounded and the total weight is a fixed negative constant. Biologically, this minimization problem is motivated by the question of determining the optimal spatial arrangement of favorable and unfavorable regions for a species to survive. For rectangular domains with Neumann boundary condition, it is known that there exists a threshold value such that if the total weight is below this threshold value then the optimal favorable region is like a section of a disk at one of the four corners; otherwise, the optimal favorable region is a strip attached to the shorter side of the rectangle. Here, we investigate the same problem with mixed Robin-Neumann type boundary conditions and study how this boundary condition affects the optimal spatial arrangement.     

A mixed method of subspace iteration for Dirichlet eigenvalue problems

A full multigrid scheme was used in computing some eigenvalues of the Laplace eigenvalue problem with the Dirichlet boundary condition. We get a system of algebraic equations with an aid of finite difference method and apply subspace iteration method to the system to compute first some eigenvalues. The result shows that this is very effective in calculating some eigenvalues of this problem.     

A full multigrid method for nonlinear eigenvalue problems

We introduce a type of full multigrid method for the nonlinear eigenvalue problem. The main idea is to transform the solution of the nonlinear eigenvalue problem into a series of solutions of the corresponding linear boundary value problems on the sequence of finite element spaces and nonlinear eigenvalue problems on the coarsest finite element space. The linearized boundary value problems are solved by some multigrid iterations. Besides the multigrid iteration, all other efficient iteration methods for solving boundary value problems can serve as the linear problem solver. We prove that the computational work of this new scheme is truly optimal, the same as solving the linear corresponding boundary value problem. In this case, this type of iteration scheme certainly improves the overfull efficiency of solving nonlinear eigenvalue problems. Some numerical experiments are presented to validate the efficiency of the new method.     

On eigenvalue problems with Robin type boundary conditions having indefinite coefficients

A Robin type boundary condition with a sign-changing coefficient is treated. First, the associated linear elliptic eigenvalue problem is studied, where the existence of a principal eigenvalue is discussed by the use of a variational approach. Second, the associated semilinear elliptic boundary value problem of logistic type is studied and the one parameter-dependent structure of positive solutions is investigated, where results obtained are due to the construction of suitable super- and subsolutions by using the principal positive eigenfunctions of the linear eigenvalue problem.     

Two-Parameter Asymptotics in a Boundary-Value Problem for the Laplacian

We study a model boundary-value problem for the Laplacian in the unit disk with closely-spaced and periodic alternation of the type of boundary condition for the case in which the Dirichlet problem is the limit one. We study and justify the two-parameter asymptotics of an eigenvalue of the perturbed problem converging to a simple eigenvalue of the limit problem.     

高阶椭圆算子Dirichlet本征值的下界

设Ω是R^m(m≥2）中的一个有界区域,其边界足够光滑,考察2p(p≥1）阶椭圆算子（－1）^p ∑│α│＝│β│＝pa^α(Aαβa^β)在Dirichlet边界条件下的本征值问题,给出了其本征值的一个下界,该下界除与维数m有关外仅依赖于区域Ω的体积。     

$p(x)$-Laplacian方程Robin边界条件下的特征值问题

本文研究了Robin边界条件下$p(x)$-Laplacian方程特征值问题. 利用变指数Sobolev空间理论, 我们用Luxemburg范数来定义Rayleigh商, 并给出该Rayleigh商的最小值点对应的Euler-Lagrange方程. 根据Ljusternik-Schnirelman原理, 我们证明了Robin边值问题存在无穷多特征值序列, 其中最小的特征值存在且是严格大于零的, 并且与最小的特征值相对应的特征函数不变号.     

Nonlinear Eigenvalue Problems in the Stability Analysis of Morphogen Gradients

This paper is concerned with several eigenvalue problems in the linear stability analysis of steady state morphogen gradients for several models of Drosophila wing imaginal discs including one not previously considered. These problems share several common difficulties including the following: (a) The steady state solution which appears in the coefficients of the relevant differential equations of the stability analysis is only known qualitatively and numerically. (b) Though the governing differential equations are linear, the eigenvalue parameter appears nonlinearly after reduction to a problem for one unknown. (c) The eigenvalues are determined not only as solutions of a homogeneous boundary value problem with homogeneous Dirichlet boundary conditions, but also by an alternative auxiliary condition to one of the Dirichlet conditions allowed by a boundary condition of the original problem. Regarding the stability of the steady state morphogen gradients, we prove that the eigenvalues must all be positive and hence the steady state morphogen gradients are asymptotically stable. The other principal finding is a novel result pertaining to the smallest (positive) eigenvalue that determines the slowest decay rate of transients and the time needed to reach steady state. Here we prove that the smallest eigenvalue does not come from the nonlinear Dirichlet eigenvalue problem but from the complementary auxiliary condition requiring only to find the smallest zero of a rational function. Keeping in mind that even the steady state solution needed for the stability analysis is only known numerically, not having to solve the nonlinear Dirichlet eigenvalue problem is both an attractive theoretical outcome and a significant computational simplification.     

The ice‐fishing problem: the fundamental sloshing frequency versus geometry of holes

We study an eignevalue problem with a spectral parameter in a boundary condition. This problem for the Laplace equation is relevant to sloshing frequencies that describe free oscillations of an inviscid, incompressible, heavy fluid in a half‐space covered by a rigid dock with some apertures (an ice sheet with fishing holes). The dependence of the fundamental eigenvalue on holes' geometry is investigated. We give conditions on a plane region guaranteeing that the fundamental eigenvalue corresponding to this region is larger than the fundamental eigenvalue corresponding to a single circular hole. Examples of regions satisfying these conditions and having the same area as the unit disk are given. New results are also obtained for the problem with a single circular hole. On the other hand, we construct regions for which the fundamental eigenfrequency is larger than the similar frequency for the circular hole of the same area and even as large as one wishes. In the latter examples, the hole regions are either not connected or bounded by a rather complicated curves. Copyright © 2004 John Wiley & Sons, Ltd.     

An accurate a posteriori error estimator for the Steklov eigenvalue problem and its applications

In this paper, a type of accurate a posteriori error estimator is proposed for the Steklov eigenvalue problem based on the complementary approach, which provides an asymptotic exact estimate for the approximate eigenpair. Besides, we design a type of cascadic adaptive finite element method for the Steklov eigenvalue problem based on the proposed a posteriori error estimator. In this new cascadic adaptive scheme, instead of solving the Steklov eigenvalue problem in each adaptive space directly, we only need to do some smoothing steps for linearized boundary value problems on a series of adaptive spaces and solve some Steklov eigenvalue problems on a low dimensional space. Furthermore, the proposed a posteriori error estimator provides the way to refine mesh and control the number of smoothing steps for the cascadic adaptive method. Some numerical examples are presented to validate the efficiency of the algorithm in this paper.