Zero-finder methods derived from Obreshkov’s techniques

In this paper two families of zero-finding iterative methods for solving nonlinear equations f(x)=0 are presented. The key idea to derive them is to solve an initial value problem applying Obreshkov-like techniques. More explicitly, Obreshkov’s methods have been used to numerically solve an initial value problem that involves the inverse of the function f that defines the equation. Carrying out this procedure, several methods with different orders of local convergence have been obtained. An analysis of the efficiency of these methods is given. Finally we introduce the concept of extrapolated computational order of convergence with the aim of numerically test the given methods. A procedure for the implementation of an iterative method with an adaptive multi-precision arithmetic is also presented.     

Numerical transformation methods: Blasius problem and its variants

Blasius problem is the simplest nonlinear boundary-layer problem. We hope that any approach developed for this epitome can be extended to more difficult hydrodynamics problems. With this motivation we review the so called Töpfer transformation, which allows us to find a non-iterative numerical solution of the Blasius problem by solving a related initial value problem and applying a scaling transformation. The applicability of a non-iterative transformation method to the Blasius problem is a consequence of its partial invariance with respect to a scaling group. Several problems in boundary-layer theory lack this kind of invariance and cannot be solved by non-iterative transformation methods. To overcome this drawback, we can modify the problem under study by introducing a numerical parameter, and require the invariance of the modified problem with respect to an extended scaling group involving this parameter. Then we apply initial value methods to the most recent developments involving variants and extensions of the Blasius problem.     

On multistep interval methods for solving the initial value problem

In this paper we shortly complete our previous considerations on interval versions of Adams multistep methods [M. Jankowska, A. Marciniak, Implicit interval multistep methods for solving the initial value problem, Comput. Meth. Sci. Technol. 8(1) (2002) 17–30; M. Jankowska, A. Marciniak, On explicit interval methods of Adams–Bashforth type, Comput. Meth. Sci. Technol. 8(2) (2002) 46–57; A. Marciniak, Implicit interval methods for solving the initial value problem, Numerical Algorithms 37 (2004) 241–251]. It appears that there exist two families of implicit interval methods of this kind. More considerations are dealt with two new kinds of interval multistep methods based on conventional well-known Nyström and Milne–Simpson methods. For these new interval methods we prove that the exact solution of the initial value problem belongs to the intervals obtained. Moreover, we present some estimations of the widths of interval solutions. Some conclusions bring this paper to the end.     

Symmetric multistep Obrechkoff methods with zero phase-lag for periodic initial value problems of second order differential equations

In this paper, symmetric multistep Obrechkoff methods of orders 8 and 12, involving a parameter p to solve a special class of second order initial value problems in which the first order derivative does not appear explicitly, are discussed. It is shown that the methods have zero phase-lag when p is chosen as 2π times the frequency of the given initial value problem.     

Zero-finder methods derived using Runge-Kutta techniques

In this paper some families of zero-finding iterative methods for nonlinear equations are presented. The key idea to derive them is to solve an initial value problem applying Runge-Kutta techniques. More explicitly, these methods are used to solve the problem that consists in a differential equation in what appears the inverse function of the one which zero will be computed and the condition given by the value attained by it at the initial approximation. Carrying out this procedure several families of different orders of local convergence are obtained. Furthermore, the efficiency of these families are computed and two new families using like-Newton’s methods that improve the most efficient one are also given.     

On efficient numerical methods for an initial-boundary value problem with nonlocal boundary conditions

Many physical phenomena are modeled by nonclassical hyperbolic boundary value problems with nonlocal boundary conditions. In this paper, the problem of solving the one-dimensional wave equation subject to given initial and non-local boundary conditions is considered. These non-local conditions arise mainly when the data on the boundary cannot be measured directly. Several finite difference methods with low order have been proposed in other papers for the numerical solution of this one dimensional non-classic boundary value problem. Here, we derive a new family of efficient three-level algorithms with higher order to solve the wave equation and also use a Simpson formula with higher order to approximate the integral conditions. Additionally, the fourth-order formula is also adapted to nonlinear equations, in particular to the well-known nonlinear Klein–Gordon equations which many physical problems can be modeled with. Numerical results are presented and are compared with some existing methods showing the efficiency of the new algorithms.     

Continuous extensions to Nyström methods for second order initial value problems

In this work we consider interpolants for Nyström methods, i.e., methods for solving second order initial value problems. We give a short introduction to the theory behind the discrete methods, and extend some of the work to continuous, explicit Nyström methods. Interpolants for continuous, explicit Runge-Kutta methods have been intensively studied by several authors, but there has not been much effort devoted to continuous Nyström methods. We therefore extend some of the work by Owren.     

Towards a new framework for evaluating systemic problem structuring methods

Operational researchers and social scientists often make significant claims for the value of systemic problem structuring and other participative methods. However, when they present evidence to support these claims, it is usually based on single case studies of intervention. There have been very few attempts at evaluating across methods and across interventions undertaken by different people. This is because, in any local intervention, contextual factors, the skills of the researcher and the purposes being pursued by stakeholders affect the perceived success or failure of a method. The use of standard criteria for comparing methods is therefore made problematic by the need to consider what is unique in each intervention. So, is it possible to develop a single evaluation approach that can support both locally meaningful evaluations and longer-term comparisons between methods? This paper outlines a methodological framework for the evaluation of systemic problem structuring methods that seeks to do just this.     

热传导型半导体瞬态问题的迎风有限体积元方法

半导体瞬态问题的数学模型是由四个方程组成的非线性偏微分方程组的初边值问题所决定．其中电子浓度和空穴浓度方程往往是对流占优扩散问题,普通的方法已不适用,为此本文用迎风格式处理对流项部分,提出一种全离散迎风有限体积元方法,并进行收敛性分析,在最一般的情况下得到了一阶精度L2模误差估计结果．     

B-CONVERGENCE PROPERTIES OF GENERAL LINEAR METHODS

In this paper, for general linear methods applied to strictly dissipative initial value problem in Hilbert spaces, we prove that algebraic stability implies B-convergence, which extends and improves the existing results on Runge-Kutta methods. Specializing our results for the case of multi-step Runge-Kutta methods, a series of B-convergence results are obtained.     

Construction of Solutions and L^1-error Estimates of Viscous Methods for Scalar Conservation Laws with Boundary

This paper is concerned with an initial boundary value problem for strictly convex conservation laws whose weak entropy solution is in the piecewise smooth solution class consisting of finitely many discontinuities. By the structure of the weak entropy solution of the corresponding initial value problem and the boundary entropy condition developed by Bardos-Leroux Nedelec, we give a construction method to the weak entropy solution of the initial boundary value problem. Compared with the initial value problem, the weak entropy solution of the initial boundary value problem includes the following new interaction type： an expansion wave collides with the boundary and the boundary reflects a new shock wave which is tangent to the boundary. According to the structure and some global estimates of the weak entropy solution, we derive the global L^1-error estimate for viscous methods to this initial boundary value problem by using the matching travelling wave solutions method. If the inviscid solution includes the interaction that an expansion wave collides with the boundary and the boundary reflects a new shock wave which is tangent to the boundary, or the inviscid solution includes some shock wave which is tangent to the boundary, then the error of the viscosity solution to the inviscid solution is bounded by O（ε^1/2） in L^1-norm; otherwise, as in the initial value problem, the L^1-error bound is O（ε｜ In ε｜）.     

A boundary value approach to the numerical solution of initial value problems by multistep methos †

A boundary value appraoch to the numerical solution of initial value problems by means of linear multistep methods is presented. This theory is based on the study of linear difference equations when their general solution is computed by imposing boundary conditions. All the main stability and convergence properties of the obtained methods are investigated abd compared to those of the classical multistep methods. Then, as an example, new itegration formulas, called extended trapezoidal rules, are derived. For any order they have the same stability properties (in the sense of the definitions given in this paper) of the trapezoidal rule, which is the first method in this class. Some numerical examples are presented to confirm the theoretical expectations and to allow us to trust a future code based on boundary value methods.     

Trigonometric–fitted hybrid four–step methods of sixth order for solving

A two–stage, explicit, hybrid four–step method of sixth order for the solution of the special second order initial value problem is presented here. The new method is trigonometric fitted, thus it uses variable coefficients. Numerical tests illustrate the superiority of our proposal over similar methods found in the relevant literature on a set of standard problems.     

A Survey of Algebraic Multilevel Iteration (AMLI) Methods

For the solution by preconditioned conjugate gradient methods of symmetric positive definite equations as arising in boundary value problems we consider preconditioning methods of AMLI type. Particular attention is devoted to providing methods of optimal order of computational complexity which in addition promise to be robust, i.e. with a convergence rate which is bounded above independently of size of discretization parameter h, jumps in problem coefficients, and shape of finite elements or, equivalently, anisotropy of problem coefficients. In addition, the computational cost per iteration step must have optimal order.New results on upper bounds of one of the important parameters in the methods, the Cauchy—Bunyakowski—Schwarz constant are given and an algebraic method how to improve its value is presented.This revised version was published online in October 2005 with corrections to the Cover Date.     

On the order of deferred correction

New deferred correction methods for the numerical solution of initial value problems in ordinary differential equations have recently been introduced by Dutt, Greengard and Rokhlin. A convergence proof is presented for these methods, based on the abstract Stetter-Lindberg-Skeel framework and Spijker-type norms. It is shown that p corrections of an order-r one-step solver yield order-r(p+1) accuracy.     

Energy-preserving trigonometrically fitted continuous stage Runge-Kutta-Nyström methods for oscillatory Hamiltonian systems

Numerical Algorithms - Recently, continuous-stage Runge-Kutta-Nyström (CSRKN) methods for solving numerically second-order initial value problem $q^{\prime \prime }= f(q)$ have been proposed...     

Ninth-order,explicit, two-step methods for second-order inhomogeneous linear IVPs

In the present paper, the second-order linear inhomogeneous initial value problems (LI-IVP) with constant coefficients is of our interest. Two-step hybrid methods (ie, of Numerov-type) are considered for addressing this problem. Thus, a special set of order conditions is given and solved for derivation a low cost new method. Actually, we manage to save two stages per step. This is crucial as shown in various numerical tests where our new proposal outperforms standard methods in the relevant literature.     

New Bundle Methods for Solving Lagrangian Relaxation Dual Problems

Bundle methods have been used frequently to solve nonsmooth optimization problems. In these methods, subgradient directions from past iterations are accumulated in a bundle, and a trial direction is obtained by performing quadratic programming based on the information contained in the bundle. A line search is then performed along the trial direction, generating a serious step if the function value is improved by or a null step otherwise. Bundle methods have been used to maximize the nonsmooth dual function in Lagrangian relaxation for integer optimization problems, where the subgradients are obtained by minimizing the performance index of the relaxed problem. This paper improves bundle methods by making good use of near-minimum solutions that are obtained while solving the relaxed problem. The bundle information is thus enriched, leading to better search directions and less number of null steps. Furthermore, a simplified bundle method is developed, where a fuzzy rule is used to combine linearly directions from near-minimum solutions, replacing quadratic programming and line search. When the simplified bundle method is specialized to an important class of problems where the relaxed problem can be solved by using dynamic programming, fuzzy dynamic programming is developed to obtain efficiently near-optimal solutions and their weights for the linear combination. This method is then applied to job shop scheduling problems, leading to better performance than previously reported in the literature.     

P-stable exponentially-fitted Obrechkoff methods of arbitrary order for second-order differential equations

We consider the construction of P-stable exponentially-fitted symmetric two-step Obrechkoff methods for solving second order differential equations related to an initial value problem. Our approach is based on two ideas: for the exponential fitting, we follow the ideas of Ixaru and Vanden Berghe; for the P-stability we introduce exponentially-fitted Padé approximants to the exponential function. By combining these two ideas, we are able to construct P-stable methods of arbitrary (even) order.