Explicit multistep method for the numerical solution of stiff differential equations

An explicit multistep method of variable order for integrating stiff systems with high accuracy and low computational costs is examined. To stabilize the computational scheme, componentwise estimates are used for the eigenvalues of the Jacobian matrix having the greatest moduli. These estimates are obtained at preliminary stages of the integration step. Examples are given to demonstrate that, for certain stiff problems, the method proposed is as efficient as the best implicit methods.     

The combined Laplace transform and new homotopy perturbation methods for stiff systems of ODEs

In this paper, we propose new technique for solving stiff system of ordinary differential equations. This algorithm is based on Laplace transform and homotopy perturbation methods. The new technique is applied to solving two mathematical models of stiff problem. We show that the present approach is relatively easy, efficient and highly accurate.     

Period two implies chaos for a class of ODEs

We extend a result of J. Andres and K. Pastor concerning scalar time-periodic first order ordinary differential equations without uniqueness, by proving that the existence of just one subharmonic implies the existence of large sets of subharmonics of all given orders. Since these periodic solutions must coexist with complicated dynamics, we might paraphrase T. Y. Li and J. A. Yorke by loosely saying that in this setting even period two implies chaos. Similar results are obtained for a class of differential inclusions.

     

A review of theoretical and numerical analysis for nonlinear stiff Volterra functional differential equations

In this review, we present the recent work of the author in comparison with various related results obtained by other authors in literature. We first recall the stability, contractivity and asymptotic stability results of the true solution to nonlinear stiff Volterra functional differential equations (VFDEs), then a series of stability, contractivity, asymptotic stability and B-convergence results of Runge-Kutta methods for VFDEs is presented in detail. This work provides a unified theoretical foundation for the theoretical and numerical analysis of nonlinear stiff problems in delay differential equations (DDEs), integro-differential equations (IDEs), delayintegro-differential equations (DIDEs) and VFDEs of other type which appear in practice.      

线性微分方程的微分算子级数解法

介绍了微分算子级数法及其求解线性常微分方程通解、特解的原理、方法和实例.这个方法和其它解法的差别,在于不借助其它学科知识的启示,直接通过方程中微分算子的运算求出方程的特解或通解.     

An almost L‐stable BDF‐type method for the numerical solution of stiff ODEs arising from the method of lines

A new BDF‐type scheme is proposed for the numerical integration of the system of ordinary differential equations that arises in the Method of Lines solution of time‐dependent partial differential equations. This system is usually stiff, so it is desirable for the numerical method to solve it to have good properties concerning stability. The method proposed in this article is almost L‐stable and of algebraic order three. Numerical experiments illustrate the performance of the new method on different stiff systems of ODEs after discretizing in the space variable some PDE problems. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

二阶微分方程第二边值问题模糊状态的数值求解

本讨论二阶微分方程的第二边值问题具有模糊不确定性时，运用模糊仿真原理和差分方法，求其边值问题的数值解法。     

Adaptivity and computational complexity in the numerical solution of ODEs

In this paper we analyze the problem of adaptivity for one-step numerical methods for solving ODEs, both IVPs and BVPs, with a view to generating grids of minimal computational cost for which the local error is below a prescribed tolerance (optimal grids). The grids are generated by introducing an auxiliary independent variable τ and finding a grid deformation map, t=Θ(τ), that maps an equidistant grid {τ_j} to a non-equidistant grid in the original independent variable, {t_j}. An optimal deformation map Θ is determined by a variational approach. Finally, we investigate the cost of the solution procedure and compare it to the cost of using equidistant grids. We show that if the principal error function is non-constant, an adaptive method is always more efficient than a non-adaptive method.     

Numerical Solution of Boundary Value Problems for Selfadjoint Differential Equations of 2nth Order

The paper is devoted to solving boundary value problems for self-adjoint linear differential equations of 2nth order in the case that the corresponding differential operator is self-adjoint and positive semidefinite. The method proposed consists in transforming the original problem to solving several initial value problems for certain systems of first order ODEs. Even if this approach may be used for quite general linear boundary value problems, the new algorithms described here exploit the special properties of the boundary value problems treated in the paper. As a consequence, we obtain algorithms that are much more effective than similar ones used in the general case. Moreover, it is shown that the algorithms studied here are numerically stable.     

一类强阻尼非线性波动方程的广义解

用Galerkin方法研究了一类非线性波动方程的初边值问题,证明其在一定条件下强解的存在性.     

Spectral methods for some singularly perturbed third order ordinary differential equations

Spectral methods with interface point are presented to deal with some singularly perturbed third order boundary value problems of reaction-diffusion and convection-diffusion types. First, linear equations are considered and then non-linear equations. To solve non-linear equations, Newton’s method of quasi-linearization is applied. The problem is reduced to two systems of ordinary differential equations. And, then, each system is solved using spectral collocation methods. Our numerical experiments show that the proposed methods are produce highly accurate solutions in little computer time when compared with the other methods available in the literature.      

常微分方程初值问题解的一个存在唯一性结果

提出并证明了常微分方程初值问题解的一个存在唯一性结果.     

Novel numerical techniques based on Fokas transforms,for the solution of initial boundary value problems

The unified transform method of A. S. Fokas has led to important new developments, regarding the analysis and solution of various types of linear and nonlinear PDE problems. In this work we use these developments and obtain the solution of time-dependent problems in a straightforward manner and with such high accuracy that cannot be reached within reasonable time by use of the existing numerical methods. More specifically, an integral representation of the solution is obtained by use of the A. S. Fokas approach, which provides the value of the solution at any point, without requiring the solution of linear systems or any other calculation at intermediate time levels and without raising any stability problems. For instance, the solution of the initial boundary value problem with the non-homogeneous heat equation is obtained with accuracy 10⁻¹⁵, while the well-established Crank–Nicholson scheme requires 2048 time steps in order to reach a 10⁻⁸ accuracy.     

Numerical methods based on multipoint Hermite interpolating polynomials for solving the Cauchy problem for stiff systems of ordinary differential equations

Families of A-, L-, and L(δ)-stable methods are constructed for solving the Cauchy problem for a system of ordinary differential equations (ODEs). The L(δ)-stability of a method with a parameter δ ∈ (0, 1) is defined. The methods are based on the representation of the right-hand sides of an ODE system at the step h in terms of two-or three-point Hermite interpolating polynomials. Comparative results are reported for some test problems. The multipoint Hermite interpolating polynomials are used to derive formulas for evaluating definite integrals. Error estimates are given.     

An analysis of the order of Runge-Kutta methods that use an iterative scheme to compute their internal stage values

We employ B-series to analyze the order of Runge-Kutta (RK) methods that use simple iteration, modified Newton iteration or full Newton iteration to compute their internal stage values. Our assumptions about the initial guess for the internal stage values cover a wide range of practical schemes. Moreover, the analytical techniques developed in this paper can be applied in many other similar contexts.This work was supported in part by the Natural Sciences and Engineering Research Council of Canada, the Information Technology Research Centre of Ontario, and the Norwegian Research Council for Science and Humanities.     

Discretized Multisplitting AOR Waveform Relaxation Algorithms for Initial Value Problem of Systems of ODEs

1.IntroductionInthedevelopmentofnewelectricalcircuits,thesimulationofthebehaviourofthecircuithasbecomeanessentialtoolforelectricalengineers.Fromthelayoutofthecircuitanonlinearsystemofordinarydifferentialequationsisgeneratedwhichdescribesthedynamicalbehaviourofthecircuit.Inthesimulationofverylargescaleintegrated(VLSI)circuitsthedimensionofthesystemofODEscanbecomeverylarge.Moreoversincethesystemisstiff,solvingthesesystemsisaverycomputionallyintensivetaskandtheuseofsupercomputersbecomesin-evit…     

关于四阶变系数线性微分方程解的不稳定性

本文利用李雅普诺夫第二方法[1]给出了至少有一个特征根具有正实部的四阶变系数线性微分方程解的不稳定性的充分条件。     

变系数线性微分方程初值问题数值解的小波方法

通过利用小波尺度函数的正交性并结合配点法 ,本文给出了一种求解变系数线性微分方程初值问题数值解的小波算法 .在一定的假设条件下 ,对算法的收敛性进行了理论分析 .最后 ,我们还给出了一个具体的数值计算例子 .