Stability properties of numerical methods for solving delay differential equations

Stability properties of numerical methods for delay differential equations are considered. Some suitable definitions for the stability of the numerical methods are included and Runge-Kutta type methods satisfying these properties are tested on a numerical example.     

时滞积分微分方程的Rosenbrock方法

本文研究时滞积分微分方程的数值方法.通过改造现有常及离散型延迟微分方程的数值方法,并匹配以适当数值求积公式,构造了求解时滞积分微分方程的Rosenbrock方法,导出了其稳定性准则.数值例子阐明了所获方法的计算有效性.     

Multi-modal method for chatter stability prediction and control in milling of thin-walled workpiece

Chatter stability in milling can be predicted by analytical methods or numerical methods. The system should be considered as multi-modal in milling of thin-walled workpiece. This paper proposes a numerical difference method based on Adams-Bashforth scheme. Moreover, multi-modal scheme of numerical methods is proposed. Analytical methods and numerical methods are verified by performing a series of milling trials. The experimental results are consistent with the predicted critical stability boundaries. Moreover, a new method for analyzing the computational time of analytical methods and numerical methods, which is based on the time complexity modeling, is presented. Computational time can be expressed as exact mathematical expression. By using the expression, the rate of increase of computational time can be derived.     

Convergence of general linear methods on differential-algebraic systems of index 3

Many numerical methods used to solve ordinary differential equations or differential-algebraic equations can be written as general linear methods. The purpose of this paper is to extend the known convergence results for Runge-Kutta and linear multistep methods to a large class of new promising numerical schemes. The theoretical results are illustrated by some numerical experiments.     

Numerical Stability and Oscillations of Runge-Kutta Methods for Differential Equations with Piecewise Constant Arguments of Advanced Type

For differential equations with piecewise constant arguments of advanced type, numerical stability and oscillations of Runge-Kutta methods are investigated. The necessary and sufficient conditions under which the numerical stability region contains the analytic stability region are given. The conditions of oscillations for the Runge-Kutta methods are obtained also. We prove that the Runge-Kutta methods preserve the oscillations of the analytic solution. Moreover, the relationship between stability and oscillations is discussed. Several numerical examples which confirm the results of our analysis are presented.     

Comparing the contractivity properties of semi-implicit methods

Contractivity is a desirable property of numerical integration methods for stiff systems of ordinary differential equations. In this paper, numerical parameters are used to allow a direct and quantitative comparison of the contractivity properties of various methods for non-linear stiff problems. Results are provided for popular Rosenbrock methods and some more recently developed semi-implicit methods.     

Numerical modeling of unsteady flow in steam turbine stage

This work deals with numerical solution of unsteady flow in turbine stage. We use models of compressible single-phase flow of air and two-phase flow of wet steam. Presented numerical methods are based on different stator-rotor matching algorithms, as well as different numerical schemes. Numerical results achieved by both methods and flow models are discussed.     

Numerical Determination of Two Characteristics of a Compressible Ion Exchanger

Two mathematical models of ion exchange with a compressible ion exchanger are considered. The inverse problems are analyzed, uniqueness of solutions is proved, and numerical methods are proposed. The efficiency of the numerical methods is investigated in computer experiments.     

On the stability of Volterra integral equations with a lagging argument

Stability properties of numerical methods for Volterra integral equations with lagging argument are considered. Some suitable definitions for the stability of the numerical methods are included, and plots of the stability regions for particular cases are included.     

Numerical solution of the finite horizon stochastic linear quadratic control problem

The treatment of the stochastic linear quadratic optimal control problem with finite time horizon requires the solution of stochastic differential Riccati equations. We propose efficient numerical methods, which exploit the particular structure and can be applied for large‐scale systems. They are based on numerical methods for ordinary differential equations such as Rosenbrock methods, backward differentiation formulas, and splitting methods. The performance of our approach is tested in numerical experiments.     

The Numerical Simulation of Space-Time Variable Fractional Order Diffusion Equation

Many physical processes appear to exhibit fractional order behavior that may vary with time or space. The continuum of order in the fractional calculus allows the order of the fractional operator to be considered as a variable. Numerical methods and analysis of stability and convergence of numerical scheme for the variable fractional order partial differential equations are quite limited and difficult to derive. This motivates us to develop efficient numerical methods as well as stability and convergence of the implicit numerical methods for the space-time variable fractional order diffusion equation on a finite domain. It is worth mentioning that here we use the Coimbra-definition variable time fractional derivative which is more efficient from the numerical standpoint and is preferable for modeling dynamical systems. An implicit Euler approximation is proposed and then the stability and convergence of the numerical scheme are investigated. Finally, numerical examples are provided to show that the implicit Euler approximation is computationally efficient.     

Synthesis of Numerical Methods for Pareto Set Approximation Based on a Universal Procedure

Computational Mathematics and Mathematical Physics - A universal computational procedure is used to construct numerical methods for Pareto set approximation. The numerical methods are developed...     

Approximate analytical solutions of a class of boundary layer equations over nonlinear stretching surface

Third order nonlinear ordinary differential equations, subject to appropriate boundary conditions arising in fluid dynamics, are solved using three different methods viz., the Dirichlet series, method of stretching of variables, and asymptotic function method. Similarity transformations are used to convert the governing partial differential equations into nonlinear ordinary differential equations. The numerical results obtained from the above methods for various problems are given in terms of skin friction. Our study revealed that the results obtained from these methods agree well with those of direct numerical simulation of ordinary differential equations. Also, these methods have advantages over pure numerical methods in obtaining derived quantities such as velocity profile accurately for various values of the parameters at a stretch.     

Numerical solution of the three‐dimensional parabolic equation with an integral condition

Developement of numerical methods for obtaining approximate solutions to the three dimensional diffusion equation with an integral condition will be carried out. The numerical techniques discussed are based on the fully explicit (1,7) finite difference technique and the fully implicit (7,1) finite difference method and the (7,7) Crank‐Nicolson type finite difference formula. The new developed methods are tested on a problem. Truncation error analysis and numerical examples are used to illustrate the accuracy of the new algorithms. The results of numerical testing show that the numerical methods based on the finite difference techniques discussed in the present article produce good results. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 193–202, 2002; DOI 10.1002/num.1040     

Inverse linear multistep methods for the numerical solution of initial value problems of second order differential equations

In this paper inverse linear multistep methods for the numerical solution of second order differential equations are presented. Local accuracy and stability of the methods are defined and discussed. The methods are applicable to a class of special second order initial value problems, not explicitly involving the first derivative. The methods are not convergent, but yield good numerical results if applied to problems they are designed for. Numerical results are presented for both the linear and nonlinear initial value problems.     

Delay-dependent stability analysis of multistep methods for delay differential equations

This paper deals with the delay-dependent stability of numerical methods for delay differential equations. First, a stability criterion of Runge-Kutta methods is extended to the case of general linear methods. Then, linear multistep methods are considered and a class of r（0）-stable methods are found. Later, some examples of r（0）-stable multistep multistage methods are given. Finally, numerical experiments are presented to confirm the theoretical results.     

Stability analysis of parabolic partial differential equations with piecewise continuous arguments

This article deals with the analytic and numerical stability of numerical methods for a parabolic partial differential equation with piecewise continuous arguments of alternately retarded and advanced type. First, application of the theory of separation of variables in matrix form and the Fourier method, the necessary and sufficient condition under which the analytic solution is asymptotically stable is derived. Then, the θ‐methods are applied to solve the corresponding initial value problem, the sufficient conditions for the asymptotic stability of numerical methods are obtained. Finally, several numerical examples are presented to support the theoretical results. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 531–545, 2017     

Finite-difference methods for solving the reaction-diffusion equations of a simple isothermal chemical system

Numerical methods are proposed for the numerical solution of a system of reaction-diffusion equations, which model chemical wave propagation. The reaction terms in this system of partial differential equations contain nonlinear expressions. Nevertheless, it is seen that the numerical solution is obtained by solving a linear algebraic system at each time step, as opposed to solving a nonlinear algebraic system, which is often required when integrating nonlinear partial differential equations. The development of each numerical method is made in the light of experience gained in solving the system of ordinary differential equations, which model the well-stirred analogue of the chemical system. The first-order numerical methods proposed for the solution of this initialvalue problem are characterized to be implicit. However, in each case it is seen that the numerical solution is obtained explicitly. In a series of numerical experiments, in which the ordinary differential equations are solved first of all, it is seen that the proposed methods have superior stability properties to those of the well-known, first-order, Euler method to which they are compared. Incorporating the proposed methods into the numerical solution of the partial differential equations is seen to lead to two economical and reliable methods, one sequential and one parallel, for solving the travelling-wave problem. © 1994 John Wiley & Sons, Inc.     

Wavelets collocation methods for the numerical solution of elliptic BV problems

Based on collocation with Haar and Legendre wavelets, two efficient and new numerical methods are being proposed for the numerical solution of elliptic partial differential equations having oscillatory and non-oscillatory behavior. The present methods are developed in two stages. In the initial stage, they are developed for Haar wavelets. In order to obtain higher accuracy, Haar wavelets are replaced by Legendre wavelets at the second stage. A comparative analysis of the performance of Haar wavelets collocation method and Legendre wavelets collocation method is carried out. In addition to this, comparative studies of performance of Legendre wavelets collocation method and quadratic spline collocation method, and meshless methods and Sinc–Galerkin method are also done. The analysis indicates that there is a higher accuracy obtained by Legendre wavelets decomposition, which is in the form of a multi-resolution analysis of the function. The solution is first found on the coarse grid points, and then it is refined by obtaining higher accuracy with help of increasing the level of wavelets. The accurate implementation of the classical numerical methods on Neumann’s boundary conditions has been found to involve some difficulty. It has been shown here that the present methods can be easily implemented on Neumann’s boundary conditions and the results obtained are accurate; the present methods, thus, have a clear advantage over the classical numerical methods. A distinct feature of the proposed methods is their simple applicability for a variety of boundary conditions. Numerical order of convergence of the proposed methods is calculated. The results of numerical tests show better accuracy of the proposed method based on Legendre wavelets for a variety of benchmark problems.     

Numerical treatment of delay differential equations by Hermite Interpolation

Summary A class of numerical methods for the treatment of delay differential equations is developed. These methods are based on the wellknown Runge-Kutta-Fehlberg methods. The retarded argument is approximated by an appropriate multipoint Hermite Interpolation. The inherent jump discontinuities in the various derivatives of the solution are considered automatically.Problems with piecewise continuous right-hand side and initial function are treated too. Real-life problems are used for the numerical test and a comparison with other methods published in literature.