Multistep methods for differential algebraic equations

Multistep methods for the differential/algebraic equations (DAEs) in the form of

The application of the concept of reaction layer to the study of electrode processes coupled with the second-order chemical reactions at microelectrodes under steady-state conditions

Under steady-state conditions, the current equations of the second-order EC, ECE and DISPI reactions at microdisk, microspherical and microring electrodes are derived with the aid of the concept of the reaction layer. The conditions under which these equations would be valid are also discussed. Using these equations, methods to determine the kinetic parameters for the second-order EC, ECE and DISPI reactions are presented. The reduction of 2,6-diphenyl-pyrylium cation and oxidation of triphenylamine were investigated as examples of the second-order EC and ECE reactions.     

D-CONVERGENCE AND STABILITY OF A CLASS OF LINEAR MULTISTEP METHODS FOR NONLINEAR DDES

1. IntroductionConsider the following nonlinear delay problem{:;\f>>:v{t(tf,?,,<,>3,<'~">>,:: 5:,3Ti:,,,,, [l:::;where y: R - C",T > 0 is a delay term, f: [t.,T] x CN x CN - CN and W(t):[to -- T, tol - CN denotes a given initial function. Thoroughout this paper 9 the problem(1.1) is supposed to have a unique solution y(t), which satisfies11 y(')(t) 115 Mi, t e [to ~ T,T]here norm 11. 11 is defined by 11 x II'=< xgx > (Vx E C"), and Mi > 0 are someconstants.Definition 1.1.[1] The clas…     

CHARACTERIZATIONS OF SYMMETRIC MULTISTEP RUNGE-KUTTA METHODS

Some characterizations for symmetric multistep Runge-Kutta(RK) methods are obtained. Symmetric two-step RK methods with one and two-stages are presented. Numerical examples show that symmetry of multistep RK methods alone is not sufficient for long time integration for reversible Hamiltonian systems. This is an important difference between one-step and multistep symmetric RK methods.     

Application of ‘Operational Quadrature Methods’ in Time Domain Boundary Element Methods

The usual time domain Boundary Element Method (BEM) contains fundamentalsolutions which are convoluted with time-dependent boundary data andintegrated over the boundary surface. Here, a new approach for theevaluation of the convolution integrals, the so-called OperationalQuadrature Methods developed by Lubich, is presented. In thisformulation, the convolution integral is numerically approximated by aquadrature formula whose weights are determined using the Laplacetransform of the fundamental solution and a linear multisep method. Tostudy the behaviour of the method, the numerical convolution of afundamental solution with a unit step function is compared with theanalytical result. Then, a time domain Boundary Element formulationapplying the Operational Quadrature Methods is derived. For thisformulation only the fundamental solutions in Laplace domain arenecessary. The properties of the new formulation are studied with anumerical example.     

求解非线性中立型延迟微分方程一类线性多步方法的收敛性

本文致力于带有Lagrang插值的一类线性多步法求解非线性中立型延迟微分方程的误差分析.证明了一个p′阶的线性多步方法配上一个q阶的Lagrang插值导致一个minf[p′,q 1]阶的E-(或EB-)收敛的非线性中立型延迟微分方程数值方法.     

A multistep flux-corrected transport scheme

A multistep flux-corrected transport (MFCT) scheme is developed to achieve conservative and monotonic tracer transports for multistep dynamical cores. MFCT extends Zalesak two-time level scheme to any multistep time-differencing schemes by including multiple high-order fluxes in the antidiffusive flux, while computing the two-time level low-order monotone solution. The multistep time-differencing scheme used in this study is the third-order Adams–Bashforth (AB3) scheme implemented in a finite-volume icosahedral shallow-water model. The accuracy of AB3 MFCT is quantified by the shape-preserving advection experiments in non-divergent flow, as well as a cosine bell whose shape changes during advection in shear flow. AB3 MFCT has been shown to be insensitive to time step size. This make AB3 MFCT an attractive transport scheme for explicit high resolution model applications with small time step. MFCT is tested in shallow-water model simulations to demonstrate that the use of MFCT maintains positive-definite tracer transport, while at the same time conserving both fluid mass and tracer mass within round-off errors in the AB3 dynamic core.     

New open modified Newton Cotes type formulae as multilayer symplectic integrators

New modified open Newton Cotes integrators are introduced in this paper. For the new proposed integrators the connection between these new algorithms, differential methods and symplectic integrators is studied. Much research has been done on one step symplectic integrators and several of them have obtained based on symplectic geometry. However, the research on multistep symplectic integrators is very poor. Zhu et al. [1] studied the well known open Newton Cotes differential methods and they presented them as multilayer symplectic integrators. Chiou and Wu [2] studied the development of multistep symplectic integrators based on the open Newton Cotes integration methods. In this paper we introduce a new open modified numerical method of Newton Cotes type and we present it as symplectic multilayer structure. The new obtained symplectic schemes are applied for the solution of Hamilton’s equations of motion which are linear in position and momentum. An important remark is that the Hamiltonian energy of the system remains almost constant as integration proceeds. We have applied also efficiently the new proposed method to a nonlinear orbital problem and an almost periodic orbital problem.     

Differentiation of Rhizoma Curcumas Longae and Radix Curcumae by a Multistep Infrared Macro-Fingerprint Method

A multistep infrared macro-fingerprint method was applied to identify two Chinese herbal drugs, Rhizoma Curcumas Longae (RCL) and Radix Curcumae (RC). Fourier transform infrared (FT-IR) spectra of the two were similar to each other and consistent with the 11 peaks of the spectrum of starch. RCL had a characteristic absorption peak at approximately 1514 cm^?1 that correlated to the strong peak near 1509 cm^?1 of curcumin. Between 900 cm^?1–1700 cm^?1 of the second derivative infrared (SD-IR) spectra, with higher resolution, RCL, and curcumin had 10 common peaks. In the FT-IR and SD-IR spectra of the ethanol extract, the spectra of the RCL extract and curcumin were similar, but RC was different. According to the fingerprint characteristics of the infrared spectra for RC and its extracts, the strongest peak at 1055 cm^?1; the C-O absorption peaks at 1124 cm^?1, 1106 cm^?1, and 996 cm^?1; and the strong methylene peaks at 2925 cm^?1 and 2853 cm^?1 suggest that RC contains more saccharides. In the range of 1350 cm^?1–1700 cm^?1, RCL and RC had similar two-dimensional infrared (2D-IR) correlation spectra. Both of them had three autopeaks, but the autopeaks were located at 1458 cm^?1, 1560 cm^?1, and 1641 cm^?1 for RCL and 1458 cm^?1, 1560 cm^?1, and 1669 cm^?1 for RC, suggested that the aromatic components of the two were not identical. The average correlation for the 18 RCL and 18 RC samples were 0.9906 and 0.9878, respectively, and this method achieves a good classification of the sample type.     

一类反应扩散方程组的隐-显多步有限元方法及其分析

考虑一类非线性反应扩散方程组，提出了隐-显多步有限元格式逼近，证明了格式最优的L^2模误差估计．