Projektive Newton-Verfahren bei elliptischen Randwertaufgaben

Summary In this paper we investigate projective Newton methods for nonlinear elliptic boundary value problems. These methods yield approximations by solving the linear equations of the Newton method by a projection method, e.g. the Ritz method. Using subspaces of finite elements or polynominals we obtain error estimates and optimal convergence theorems (inH ₁ andL ₂).

     

A choice of forcing terms in inexact Newton method

Inexact Newton method is one of the effective tools for solving systems of nonlinear equations. In each iteration step of the method, a forcing term, which is used to control the accuracy when solving the Newton equations, is required. The choice of the forcing terms is of great importance due to their strong influence on the behavior of the inexact Newton method, including its convergence, efficiency, and even robustness. To improve the efficiency and robustness of the inexact Newton method, a new strategy to determine the forcing terms is given in this paper. With the new forcing terms, the inexact Newton method is locally Q-superlinearly convergent. Numerical results are presented to support the effectiveness of the new forcing terms.     

Construction of Newton-like iteration methods for solving nonlinear equations

In this paper, we present a simple, and yet powerful and easily applicable scheme in constructing the Newton-like iteration formulae for the computation of the solutions of nonlinear equations. The new scheme is based on the homotopy analysis method applied to equations in general form equivalent to the nonlinear equations. It provides a tool to develop new Newton-like iteration methods or to improve the existing iteration methods which contains the well-known Newton iteration formula in logic; those all improve the Newton method. The orders of convergence and corresponding error equations of the obtained iteration formulae are derived analytically or with the help of Maple. Some numerical tests are given to support the theory developed in this paper.     

Numerische Quadratur bei Projektionsverfahren

Summary In this paper we investigate the influence of the numerical quadrature in projection methods. In particular we derive conditions for the order of the quadrature formulas in finite element methods under which the order of convergence is not perturbed. It seems that this question has been discussed only for the Ritz method. There is an essential difference between this method on one side and the Galerkin and least squares methods on the other side. The methods using numerical integration are only in the latter case still projection methods. The resulting conditions for the quadrature formulas are often much weaker than those for the Ritz method. Numerical examples using cubic splines and polynomials show that the conditions derived are realistic. These examples also allow the comparison of some projection methods.

     

An improved approximate Newton method for implicit Runge-Kutta formulas

Implicit Runge-Kutta (IRK) methods (such as the s-stage Radau IIA method with s=3,5, or 7) for solving stiff ordinary differential equation systems have excellent stability properties and high solution accuracy orders, but their high computing costs in solving their nonlinear stage equations have seriously limited their applications to large scale problems. To reduce such a cost, several approximate Newton algorithms were developed, including a commonly used one called the simplified Newton method. In this paper, a new approximate Jacobian matrix and two new test rules for controlling the updating of approximate Jacobian matrices are proposed, yielding an improved approximate Newton method. Theoretical and numerical analysis show that the improved approximate Newton method can significantly improve the convergence and performance of the simplified Newton method.     

Variational iteration technique for solving nonlinear equations

In this paper, we use the variational iteration technique to suggest some new iterative methods for solving nonlinear equations f(x)=0. We also discuss the convergence criteria of these new iterative methods. Comparison with other similar methods is also given. These new methods can be considered as an alternative to the Newton method. We also give several examples to illustrate the efficiency of these methods. This technique can be used to suggest a wide class of new iterative methods for solving system of nonlinear equations.     

Schwarz methods for inequalities with contraction operators

We prove the convergence of some multiplicative and additive Schwarz methods for inequalities which contain contraction operators. The problem is stated in a reflexive Banach space and it generalizes the well-known fixed-point problem in the Hilbert spaces. Error estimation theorems are given for three multiplicative algorithms and two additive algorithms. We show that these algorithms are in fact Schwarz methods if the subspaces are associated with a decomposition of the domain. Also, for the one- and two-level methods in the finite element spaces, we write the convergence rates as functions of the overlapping and mesh parameters. They are similar with the convergence rates of these methods for linear problems. Besides the direct use of the five algorithms for the inequalities with contraction operators, we can use the above results to obtain the convergence rate of the Schwarz method for other types of inequalities or nonlinear equations. In this way, we prove the convergence and estimate the error of the one- and two-level Schwarz methods for some inequalities in Hilbert spaces which are not of the variational type, and also, for the Navier–Stokes problem. Finally, we give conditions of existence and uniqueness of the solution for all problems we consider. We point out that these conditions and the convergence conditions of the proposed algorithms are of the same type.     

Generalized Newton’s method based on graphical derivatives

This paper concerns developing a numerical method of the Newton type to solve systems of nonlinear equations described by nonsmooth continuous functions. We propose and justify a new generalized Newton algorithm based on graphical derivatives, which have never been used to derive a Newton-type method for solving nonsmooth equations. Based on advanced techniques of variational analysis and generalized differentiation, we establish the well-posedness of the algorithm, its local superlinear convergence, and its global convergence of the Kantorovich type. Our convergence results hold with no semismoothness and Lipschitzian assumptions, which is illustrated by examples. The algorithm and main results obtained in the paper are compared with well-recognized semismooth and B-differentiable versions of Newton’s method for nonsmooth Lipschitzian equations.     

Convergence and optimization of the parallel method of simultaneous directions for the solution of elliptic problems

For the solution of elliptic problems, fractional step methods and in particular alternating directions (ADI) methods are iterative methods where fractional steps are sequential. Therefore, they only accept parallelization at low level. In [T. Lu, P. Neittaanmäki, X.C. Tai, A parallel splitting-up method for partial differential equations and its applications to Navier–Stokes equations, RAIRO Modél. Math. Anal. Numér. 26 (6) (1992) 673–708], Lu et al. proposed a method where the fractional steps can be performed in parallel. We can thus speak of parallel fractional step (PFS) methods and, in particular, simultaneous directions (SDI) methods. In this paper, we perform a detailed analysis of the convergence and optimization of PFS and SDI methods, complementing what was done in [T. Lu, P. Neittaanmäki, X.C. Tai, A parallel splitting-up method for partial differential equations and its applications to Navier–Stokes equations, RAIRO Modél. Math. Anal. Numér. 26 (6) (1992) 673–708]. We describe the behavior of the method and we specify the good choice of the parameters. We also study the efficiency of the parallelization. Some 2D, 3D and high-dimensional tests confirm our results.     

Über ein Verfahren der Ordnung1 + \sqrt 2 zur Nullstellenbestimmung

Summary A new iterative method for solving nonlinear equations is presented which is shown to converge locally withR-order of convergence at least under suitable differentiability assumptions. The method needs as many function evaluations per step as the classical Newton method.

     

Convergence of a generalized Newton and an inexact generalized Newton algorithms for solving nonlinear equations with nondifferentiable terms

In this paper, we consider two versions of the Newton-type method for solving a nonlinear equations with nondifferentiable terms, which uses as iteration matrices, any matrix from B-differential of semismooth terms. Local and global convergence theorems for the generalized Newton and inexact generalized Newton method are proved. Linear convergence of the algorithms is obtained under very mild assumptions. The superlinear convergence holds under some conditions imposed on both terms of equation. Some numerical results indicate that both algorithms works quite well in practice.      

A new fitted operator finite difference method to solve systems of evolutionary reaction-diffusion equations

Abstract

In recent years, fitted operator finite difference methods (FOFDMs) have been developed for numerous types of singularly perturbed ordinary differential equations. The construction of most of these methods differed though the final outcome remained similar. The most crucial aspect was how the difference operator was designed to approximate the differential operator in question. Very often the approaches for constructing these operators had limited scope in the sense that it was difficult to extend them to solve even simple one-dimensional singularly perturbed partial differential equations. However, in some of our most recent work, we have successfully designed a class of FOFDMs and extended them to solve singularly perturbed time-dependent partial differential equations. In this paper, we design and analyze a robust FOFDM to solve a system of coupled singularly perturbed parabolic reaction-diffusion equations. We use the backward Euler method for the semi-discretization in time. An FOFDM is then developed to solve the resulting set of boundary value problems. The proposed method is analyzed for convergence. Our method is uniformly convergent with order one and two, respectively, in time and space, with respect to the perturbation parameters. Some numerical experiments supporting the theoretical investigations are also presented.     

Zur Konvergenz des SSOR-Verfahrens für nichtlineare Gleichungssysteme

Zusammenfassung In der vorliegenden Arbeit wird für die Konvergenz des symmetrischen Relaxationsverfahrens (SSOR-Verfahren) bei Anwendung auf eine Klassenichtlinearer Gleichungssysteme die globale Konvergenz bewiesen. Diese Gleichungssysteme treten z.B. bei der Diskretisierung nichtlinearer partieller Differentialgleichungen auf.

---

Convergence of the SSOR-method for nonlinear systems of simultaneous equations

Summary In this paper we prove the convergence of the symmetric successive overrelaxation method if it is applied to certain nonlinear systems of simultaneous equations. These equations are obtained, for example, by discretizing nonlinear elliptic partial differential equations.

Herrn Helmut Brakhage, Kaiserslautern, anläßlich seines sechzigsten Geburtstages am 8. 7. 1986 gewidmet     

Lösung von gewöhnlichen Differentialgleichungen mit nichtlinearen Splines

Zusammenfassung In dieser Arbeit werden nichtlineare Splines zur Lösung von Anfangswertaufgaben bei gewöhnlichen Differentialgleichungen herangezogen. In der Nähe von Singularitäten besitzen z.B. verallgemeinerte rationale Splines mit variablen Exponenten gute Approximationseigenschaften. Bei Polynomsplines können Konvergenzaussagen hergeleitet werden, indem Äquivalenz dieser Verfahren mit gewissen linearen Mehrschrittverfahren gezeigt wird. In dieser Arbeit behandeln wir den nichtlinearen Fall, indem wir die lokalen Fehler in den Knoten direkt verfolgen. Einige numerische Beispiele zeigen die Güte dieser Verfahren insbesondere bei solchen Lösungen, die sehr steil anwachsen oder sogar im betrachteten Intervall singulär werden.

---

Solution of ordinary differential equations with nonlinear splines

Summary We consider the technique of using nonlinear splines to solve the initial value problem of ordinary differential equations. It is known, for example, that generalized rational splines with variable exponents yield good approximations to the exact solution in the neighborhood of a singularity. In the case of polynomial splines, convergence results may be derived by demonstrating the equivalence of the method to linear multistep methods. This sort of analysis has been done by many authors. In this paper we treat the nonlinear case and are able to prove convergence by directly estimating the local errors at interior knots. Some computational examples are given which illustrate the power of the method near a singularity.

     

A cascadic multigrid algorithm for semilinear elliptic problems

Summary. We propose a cascadic multigrid algorithm for a semilinear elliptic problem. The nonlinear equations arising from linear finite element discretizations are solved by Newton's method. Given an approximate solution on the coarsest grid on each finer grid we perform exactly one Newton step taking the approximate solution from the previous grid as initial guess. The Newton systems are solved iteratively by an appropriate smoothing method. We prove that the algorithm yields an approximate solution within the discretization error on the finest grid provided that the start approximation is sufficiently accurate and that the initial grid size is sufficiently small. Moreover, we show that the method has multigrid complexity. Received February 12, 1998 / Revised version received July 22, 1999 / Published online June 8, 2000     

Dirac delta methods for Helmholtz transmission problems

In this paper we use a boundary integral method with single layer potentials to solve a class of Helmholtz transmission problems in the plane. We propose and analyze a novel and very simple quadrature method to solve numerically the equivalent system of integral equations which provides an approximation of the solution of the original problem with linear convergence (quadratic in some special cases). Furthermore, we also investigate a modified quadrature approximation based on the ideas of qualocation methods. This new scheme is again extremely simple to implement and has order three in weak norms.      

On an Aitken–Newton type method

We study the solving of nonlinear equations by an iterative method of Aitken type, which has the interpolation nodes controlled by the Newton method. We obtain a local convergence result which shows that the q-convergence order of this method is 6 and its efficiency index is $\sqrt[5]{6},$ which is higher than the efficiency index of the Aitken or Newton methods. Monotone sequences are obtained for initial approximations farther from the solution, if they satisfy the Fourier condition and the nonlinear mapping satisfies monotony and convexity assumptions on the domain.     

Newton methods for a class of nonlinear hypersingular integral equations

Different iterative schemes based on collocation methods have been well studied and widely applied to the numerical solution of nonlinear hypersingular integral equations (Capobianco et al. 2005). In this paper we apply Newton’s method and its modified version to solve the equations obtained by applying a collocation method to a nonlinear hypersingular integral equation of Prandtl’s type. The corresponding convergence results are derived in suitable Sobolev spaces. Some numerical tests are also presented to validate the theoretical results.     

Fourier analysis of schwarz alternating methods for piecewise hermite bicubic orthogonal spline collocation

The rates of convergence of two Schwarz alternating methods are analyzed for the iterative solution of a discrete problem which arises when orthogonal spline collocation with piecewise Hermite bicubics is applied to the Dirichlet problem for Poisson's equation on a rectangle. In the first method, the rectangle is divided into two overlapping subrectangles, while three overlapping subrectangles are used in the second method. Fourier analysis is used to obtain explicit formulas for the convergence factors by which theH ¹-norm of the errors is reduced in one iteration of the Schwarz methods. It is shown numerically that while these factors depend on the size of overlap, they are independent of the partition stepsize. Results of numerical experiments are presented which confirm the established rates of convergence of the Schwarz methods.This research was supported in part by funds from the National Science Foundation grant CCR-9103451.     

A COMPLETE ANALYSIS FOR SOME A POSTERIORI ERROR ESTIMATES WITH THE FINITE ELEMENT METHOD OF LINES FOR A NONLINEAR PARABOLIC EQUATION

ABSTRACT

A posteriori error estimates for semidiscrete finite element methods for a nonlinear parabolic initial-boundary value problem are considered. The error estimates are obtained by solving local parabolic or elliptic equations for corrections to the solution on each element. The convergence results improve previous results where unnecessary assumptions are imposed on the approximate solution and the elliptic projection of the exact solution.