Condition of extremum for eigenvalues of elliptic boundary-value problems

In this paper, we consider problems of eigenvalue optimization for elliptic boundary-value problems. The coefficients of the higher derivatives are determined by the internal characteristics of the medium and play the role of control. The necessary conditions of the first and second order for problems of the first eigenvalue maximization are presented. In the case where the maximum is reached on a simple eigenvalue, the second-order condition is formulated as completeness condition for a system of functions in Banach space. If the maximum is reached on a double eigenvalue, the necessary condition is presented in the form of linear dependence for a system of functions. In both cases, the system is comprised of the eigenfunctions of the initial-boundary value problem. As an example, we consider the problem of maximization of the first eigenvalue of a buckling column that lies on an elastic foundation.     

On the method of complementary functions for nonlinear boundary-value problems

We show that, like the method of adjoints, the method of complementary functions can be effectively used to solve nonlinear boundary-value problems.This work was supported by the Alexander von Humboldt Foundation. The author is thankful to Prof. G. Hämmerlin for providing the facilities and to Miss J. Gumberger for performing numerical tests. The author is also indebted to Dr. S. M. Roberts for his suggestions on the first draft of this paper.     

On Warner's algorithm for the solution of boundary-value problems for ordinary differential equations

The paper discusses the solution of boundary-value problems for ordinary differential equations by Warner's algorithm. This shooting algorithm requires that only the original system of differential equations is solved once in each iteration, while the initial conditions for a new iteration are evaluated from a matrix equation. Numerical analysis performed shows that the algorithm converges even for very bad starting values of the unknown initial conditions and that the number of iterations is small and weakly dependent on the starting point. Based on this algorithm, a general subroutine can be realized for the solution of a large class of boundary-value problems.     

On computing eigenvalues of second-order linear pencils

^** Email: mhannaby{at}yahoo.com^*** Email: zahraa26{at}yahoo.com In this paper, we use sinc techniques to compute the eigenvaluesof a second-order operator pencil of the form Q – P approximately.Here Q and P are self-adjoint differential operators of thesecond and first order, respectively. Also the eigenparameterappears in the boundary conditions linearly.     

Parallel initial-value algorithms for singularly perturbed boundary-value problems

Initial-value methods for linear and semilinear singularly perturbed boundary-value problems are examined with a view to designing and implementing algorithms on parallel architectures. Practical experiments on a CRAY Y-MP 8/432 multiprocessor have been performed, showing the reliability and performance of several proposed parallel schemes.This work was supported by CNR, Rome, Italy (Progetto Finalizzato Sistemi Informatici e Calcolo Parallelo, Sottoprogetto 1).The authors wish to thank Dr. A. Papini, who carried out most of the computations reported in this work.     

Fully discrete approximations of parabolic boundary-value problems with nonsmooth boundary data

We study a numerical scheme for the approximation of parabolic boundary-value problems with nonsmooth boundary data. This fully discrete scheme requires no boundary constraints on the approximating elements. Our principal result is the derivation of optimal convergence estimates in L_p[0,T; L₂()] norms for boundary data in L_p[0, T; L₂()], 1p . For the same algorithms, we also show that the convergence remains optimal even in higher norms. The techniques employed are based on the theory of analytic semigroups combined with singular integrals.This paper was written in 1990, when the author was in the Department of Mathematical Sciences, University of Cincinnati. A preliminary version of this research was presented at the SIAM Annual Meeting in July 1989.     

On boundary-value problems for second order perturbed differential inclusions

This article studies the existence of solutions to boundary-value problems for second order multi-valued perturbed differential inclusions under the mixed Lipschitz and Carathéodory conditions. The existence of extremal solutions is also obtained under certain monotonicity conditions and the weaker nonconvexity conditions for multi-valued functions.     

Initial-value methods for second-order singularly perturbed boundary-value problems

In this paper, we present a numerical method for solving linear and nonlinear second-order singularly perturbed boundary-value-problems. For linear problems, the method comes from the well-known WKB method. The required approximate solution is obtained by solving the reduced problem and one or two suitable initial-value problems, directly deduced from the given problem. For nonlinear problems, the quasilinearization method is applied. Numerical results are given showing the accuracy and feasibility of the proposed method.This work was supported in part by the Consiglio Nazionale delle Ricerche (Contract No. 86.02108.01 and Progetto Finalizzatto Sistemi Informatia e Calcolo Paralello, Sottoprogetto 1), and in part by the Ministero della Pubblica Istruzione, Rome, Italy.     

On the distribution of real eigenvalues in linear viscoelastic oscillators

In this paper, a linear viscoelastic system is considered where the viscoelastic force depends on the past history of motion via a convolution integral over an exponentially decaying kernel function. The free‐motion equation of this nonviscous system yields a nonlinear eigenvalue problem that has a certain number of real eigenvalues corresponding to the nonoscillatory nature. The quality of the current numerical methods for deriving those eigenvalues is directly related to damping properties of the viscoelastic system. The main contribution of this paper is to explore the structure of the set of nonviscous eigenvalues of the system while the damping coefficient matrices are rank deficient and the damping level is changing. This problem will be investigated in the cases of low and high levels of damping, and a theorem that summarizes the possible distribution of real eigenvalues will be proved. Moreover, upper and lower bounds are provided for some of the eigenvalues regarding the damping properties of the system. Some physically realistic examples are provided, which give us insight into the behavior of the real eigenvalues while the damping level is changing.     

On the distribution of Laplacian eigenvalues of trees

For a tree T$T$ with n$n$ vertices, we apply an algorithm due to Jacobs and Trevisan (2011) to study how the number of small Laplacian eigenvalues behaves when the tree is transformed by a transformation defined by Mohar (2007). This allows us to obtain a new bound for the number of eigenvalues that are smaller than 2. We also report our progress towards a conjecture on the number of eigenvalues that are smaller than the average degree.     

On eigenvalues of spherical fractal drums

We will give better estimates of a kind of nonisotropic fractal drum.     

Numerical solution of second-order nonlinear singularly perturbed boundary-value problems by initial-value methods

Nonlinear singularly perturbed boundary-value problems are considered, with one or two boundary layers but no turning points. The theory of differential inequalities is used to obtain a numerical procedure for quasilinear and semilinear problems. The required solution is approximated by combining the solutions of suitable auxiliary initial-value problems easily deduced from the given problem. From the numerical results, the method seems accurate and solutions to problems with extremely thin layers can be obtained at reasonable cost.This work was supported by CNR, Rome, Italy (Progetto Finalizzato Sistemi Informatici e Calcolo Parallelo, Sottoprogetto 1).     

Differential quadrature method (DQM) and Boubaker Polynomials Expansion Scheme (BPES) for efficient computation of the eigenvalues of fourth-order Sturm-Liouville problems

The differential quadrature method (DQM) and the Boubaker Polynomials Expansion Scheme (BPES) are applied in order to compute the eigenvalues of some regular fourth-order Sturm-Liouville problems. Generally, these problems include fourth-order ordinary differential equations together with four boundary conditions which are specified at two boundary points. These problems concern mainly applied-physics models like the steady-state Euler-Bernoulli beam equation and mechanicals non-linear systems identification. The approach of directly substituting the boundary conditions into the discrete governing equations is used in order to implement these boundary conditions within DQM calculations. It is demonstrated through numerical examples that accurate results for the first kth eigenvalues of the problem, where k = 1, 2, 3, … , can be obtained by using minimally 2(k + 4) mesh points in the computational domain. The results of this work are then compared with some relevant studies.     

On the problem of spurious eigenvalues in the approximation of linear elliptic problems in mixed form

In the approximation of linear elliptic operators in mixed form, it is well known that the so-called inf-sup and ellipticity in the kernel properties are sufficient (and, in a sense to be made precise, necessary) in order to have good approximation properties and optimal error bounds. One might think, in the spirit of Mercier-Osborn-Rappaz-Raviart and in consideration of the good behavior of commonly used mixed elements (like Raviart-Thomas or Brezzi-Douglas-Marini elements), that these conditions are also sufficient to ensure good convergence properties for eigenvalues. In this paper we show that this is not the case. In particular we present examples of mixed finite element approximations that satisfy the above properties but exhibit spurious eigenvalues. Such bad behavior is proved analytically and demonstrated in numerical experiments. We also present additional assumptions (fulfilled by the commonly used mixed methods already mentioned) which guarantee optimal error bounds for eigenvalue approximations as well.

     

A modified Rayleigh conjecture for static problems

A modified Rayleigh conjecture (MRC) in scattering theory was proposed and justified by the author [A.G. Ramm, Modified Rayleigh conjecture and applications, J. Phys. A 35 (2002) L357–L361]. The MRC allows one to develop efficient numerical algorithms for solving boundary-value problems. It gives an error estimate for solutions. In this paper the MRC is formulated and proved for static problems.     

On principal eigenvalues for boundary value problems with indefinite weight and Robin boundary conditions

We investigate the existence of principal eigenvalues (i.e., eigenvalues corresponding to positive eigenfunctions) for the boundary value problem on ; on , where is a bounded region in , is an indefinite weight function and may be positive, negative or zero.

     

On a two-point boundary-value problem in geophysics

We present some results on the existence and uniqueness of solutions to a two-point boundary-value problem that models gyre flows in rotating spherical coordinates.     

Sextic spline method for the solution of a system of obstacle problems

We use sextic spline function to develop numerical method for the solution of system of second-order boundary-value problems associated with obstacle, unilateral, and contact problems. We show that the approximate solutions obtained by the present method are better than those produced by other collocation, finite difference and spline methods. A numerical example is given to illustrate practical usefulness of our method.     

On solving the second four-element Riemann-type boundary-value problem on a disk

The paper is devoted to solving one of the main four-element Riemann-type boundary-value problems in classes of piecewise-bianalytic functions with unit circumference as jump line. We prove a theorem on reduction of solving the stated problem to solving two vector-matrix Riemann problems with respect to piecewise-analytic vector-functions. Using it, we obtain an algorithm for solving the problem and conditions under which one can get a constructive and explicit solution in terms of Cauchy-type integrals. __________ Translated from Lietuvos Matematikos Rinkinys, Vol. 46, No. 3, pp. 377–385, July–September, 2006.