Correction of finite element estimates for Sturm-Liouville eigenvalues

Summary It is shown that a simple asymptotic correction technique of Paine, de Hoog and Anderssen reduces the error in the estimate of thekth eigenvalue of a regular Sturm-Liouville problem obtained by the finite element method, with linear hat functions and mesh lengthh, fromO(k ⁴ h ²) toO(k h ²). The result still holds when the matrix elements are evaluated by Simpson's rule, but if the trapezoidal rule is used the error isO(k ² h ²). Numerical results demonstrate the usefulness of the correction even for low values ofk.     

Small potential corrections for the discrete eigenvalues of the Sturm-Liouville problem

Summary The inverse Sturm-Liouville problem is the problem of finding a good approximation of a potential functionq such that the eigenvalue problem (^*)–y +qy=y holds on (0, ) fory(0)=y()=0 and a set of given eigenvalues . Since this problem has to be solved numerically by discretization and since the higher discrete eigenvalues strongly deviate from the corresponding Sturm-Liouville eigenvalues , asymptotic corrections for the 's serve to get better estimates forq. Let _k (1kn) be the first eigenvalues of (^*), let _k be the corresponding discrete eigenvalues obtained by the finite element method for (^*) and let _{k k} for the special caseq=0. Then, starting from an asymptotic correction technique proposed by Paine, de Hoog and Anderssen, new estimates for the errors of the corrected discrete eigenvalues are obtained and confirm and improve the knownO(kh ²)(h:=/(n+1)) behaviour. The estimates are based on new Sobolev inequalities and on Fourier analysis and it is shown that for 4+c ₂ k(n+1)/2, wherec ₁ andc ₂ are constants depending onq which tend to 0 for vanishingq.     

Nonlinear eigenvalue approximation

Summary For each in some domainD in the complex plane, letF() be a linear, compact operator on a Banach spaceX and letF be holomorphic in . Assuming that there is a so thatI–F() is not one-to-one, we examine two local methods for approximating the nonlinear eigenvalue . In the Newton method the smallest eigenvalue of the operator pencil [I–F(),F()] is used as increment. We show that under suitable hypotheses the sequence of Newton iterates is locally, quadratically convergent. Second, suppose 0 is an eigenvalue of the operator pencil [I–F(),I] with algebraic multiplicitym. For fixed leth() denote the arithmetic mean of them eigenvalues of the pencil [I–F(),I] which are closest to 0. Thenh is holomorphic in a neighborhood of andh()=0. Under suitable hypotheses the classical Muller's method applied toh converges locally with order approximately 1.84.     

B-Convergence of Lobatto IIIC formulas

Summary A completion ofB-convergence results of Lobatto IIIC schemes is presented. In particular, it is shown that Lobatto IIIC schemes with more than two stages areB-convergent when applied to IVPs with a negative one-sided Lipschitz constantm; they are notB-convergent, however, for IVPs with a non-negativem.     

Equilibria of Runge-Kutta methods

Summary It is known that certain Runge-Kutta methods share the property that, in a constant-step implementation, if a solution trajectory converges to a bounded limit then it must be a fixed point of the underlying differential system. Such methods are calledregular. In the present paper we provide a recursive test to check whether given method is regular. Moreover, by examining solution trajectories of linear equations, we prove that the order of ans-stage regular method may not exceed 2[(s+2)/2] and that the maximal order of regular Runge-Kutta method with an irreducible stability function is 4.     

Rosenbrock methods for Differential Algebraic Equations

Summary This paper deals with the numerical solution of Differential/Algebraic Equations (DAE) of index one. It begins with the development of a general theory on the Taylor expansion for the exact solutions of these problems, which extends the well-known theory of Butcher for first order ordinary differential equations to DAE's of index one. As an application, we obtain Butcher-type results for Rosenbrock methods applied to DAE's of index one, we characterize numerical methods as applications of certain sets of trees. We derive convergent embedded methods of order 4(3) which require 4 or 5 evaluations of the functions, 1 evaluation of the Jacobian and 1 LU factorization per step.     

General framework,stability and error analysis for numerical stiff boundary value methods

Summary This paper provides a general framework, called theoretical multiple shooting, within which various numerical methods for stiff boundary value ordinary differential problems can be analyzed. A global stability and error analysis is given, allowing (as much as possible) the specificities of an actual numerical method to come in only locally. We demonstrate the use of our results for both one-sided and symmetric difference schemes. The class of problems treated includes some with internal (e.g. turning point) layers.     

On error behaviour of partitioned linearly implicit runge-kutta methods for stiff and differential algebraic systems

This paper studies partitioned linearly implicit Runge-Kutta methods as applied to approximate the smooth solution of a perturbed problem with stepsizes larger than the stiffness parameter. Conditions are supplied for construction of methods of arbitrary order. The local and global error are analyzed and the limiting case 0 considered yielding a partitioned linearly implicit Runge-Kutta method for differential-algebraic equations of index one. Finally, some numerical experiments demonstrate our theoretical results.     

The stability function for multistep collocation methods

Summary C-polynomials for rational approximation to the exponential function was introduced by Nørsett [7] to study stability properties of one-step methods. For one-step collocation methods theC-polynomial has a very simple form. In this paper we studyC-polynomials for multistep collocation methods and obtain results that generalize those in the one-step case, and provide a way to analyze linear stability of such methods.     

Linearly implicit extrapolation methods for differential-algebraic systems

Summary We consider the numerical solution of implicit differential equations in which the solution derivative appears multiplied by a solution-dependent singular matrix. We study extrapolation methods based on two linearly implicit Euler discretizations. Their error behaviour is explained by perturbed asymptotic expansions.     

Discretization of multiparameter eigenvalue problems

Summary Although multiparameter eigenvalue problems, as for example Mathieu's differential equation, have been known for a long time, so far no work has been done on the numerical treatment of these problems. So in this paper we extend the spectral theory for one parameter (cf. [7, II, VII]) to multiparameter eigenvalue problmes, formulate in the framework of discrete approximation a convergent numerical treatment, establish algebraic bifurcation equations for the intersection points of the eigenvalue curves and illustrate this with some numerical examples.     

Natural Runge-Kutta and projection methods

Summary Recently the author defined the class of natural Runge-Kutta methods and observed that it includes all the collocation methods. The present paper is devoted to a complete characterization of this class and it is shown that it coincides with the class of the projection methods in some polynomial spaces.This work was supported by the Italian Ministero della Pubblica Istruzione, funds 40%     

On the order conditions of Runge-Kutta methods with higher derivatives

Summary Runge-Kutta methods have been generalized to procedures with higher derivatives of the right side ofy=f(t,y) e.g. by Fehlberg 1964 and Kastlunger and Wanner 1972. In the present work some sufficient conditions for the order of consistence are derived for these methods using partially the degree of the corresponding numerical integration formulas. In particular, methods of Gauß, Radau, and Lobatto type are generalized to methods with higher derivatives and their maximum order property is proved. The applied technique was developed by Crouzeix 1975 for classical Runge-Kutta methods. Examples of simple explicit and semi-implicit methods are given up to order 7 and 6 respectively.     

Energy conserving,arbitrary order numerical solutions of theN-body problem

Summary In this paper an energy conserving modification of the well-known extrapolation methods for solving theN-body problem is presented. The method is compared with the most commonly used extrapolation methods. It appears that the presented modification yields a better approximation of the exact solution. Computer examples are given.     

The application of Rosenbrock-Wanner type methods with stepsize control in differential-algebraic equations

Summary Two Rosenbrock-Wanner type methods for the numerical treatment of differential-algebraic equations are presented. Both methods possess a stepsize control and an index-1 monitor. The first method DAE34 is of order (3)4 and uses a full semi-implicit Rosenbrock-Wanner scheme. The second method RKF4DA is derived from the Runge-Kutta-Fehlberg 4(5)-pair, where a semi-implicit Rosenbrock-Wanner method is embedded, in order to solve the nonlinear equations. The performance of both methods is discussed in artificial test problems and in technical applications.     

Partitioned adaptive Runge-Kutta methods for the solution of nonstiff and stiff systems

Summary For the numerical solution of initial value problems of ordinary differential equations partitioned adaptive Runge-Kutta methods are studied. These methods consist of an adaptive Runge-Kutta methods for the treatment of a stiff system and a corresponding explicit Runge-Kutta method for a nonstiff system. First we modify the theory of Butcher series for partitioned adaptive Runge-Kutta methods. We show that for any explicit Runge-Kutta method there exists a translation invariant partitoned adaptive Runge-Kutta method of the same order. Secondly we derive a special translaton invariant partitioned adaptive Runge-Kutta method of order 3. An automatic stiffness detection and a stepsize control basing on Richardson-extrapolation are performed. Extensive tests and comparisons with the partitioned RKF4RW-algorithm from Rentrop [16] and the partitioned algorithm LSODA from Hindmarsh [9] and Petzold [15] show that the partitoned adaptive Runge-Kutta algorithm works reliable and gives good numericals results. Furthermore these tests show that the automatic stiffness detection in this algorithm is effective.     

A relative backward perturbation theorem for the eigenvalue problem

Summary A-posteriori bound is given for the computed eigenpair ( ), of the eigenvalue problemAx=x, which is shown to be more realistic than the available one. A simple expression is further presented for calculating the backward error.     

One-step and extrapolation methods for differential-algebraic systems

Summary The paper analyzes one-step methods for differential-algebraic equations (DAE) in terms of convergence order. In view of extrapolation methods, certain perturbed asymptotic expansions are shown to hold. For the special DAE extrapolation solver based on the semi-implicit Euler discretization, the perturbed order pattern of the extrapolation tableau is derived in detail. The theoretical results lead to modifications of the known code. The efficiency of the modifications is illustrated by numerical comparisons over critical examples mainly from chemical combustion.     

An algebraic characterization ofB-convergent Runge-Kutta methods

Summary In the analysis of discretization methods for stiff intial value problems, stability questions have received most part of the attention in the past.B-stability and the equivalent criterion algebraic stability are well known concepts for Runge-Kutta methods applied to dissipative problems. However, for the derivation ofB-convergence results — error bounds which are not affected by stiffness — it is not sufficient in many cases to requireB-stability alone. In this paper, necessary and sufficient conditions forB-convergence are determined.This paper was written while J. Schneid was visiting the Centre for Mathematics and Computer Science with an Erwin-Schrödinger stipend from the Fonds zur Förderung der wissenschaftlichen Forschung