On spectral methods for Volterra-type integro-differential equations

This paper considers the spectral methods for a Volterra-type integro-differential equation. Firstly, the Volterra-type integro-differential equation is equivalently restated as two integral equations of the second kind. Secondly, a Legendre-collocation method is used to solve them. Then the error analysis is conducted based on the _L∞$L_{\infty }$-norm. In addition, numerical results are presented to confirm our analysis.     

High-order adaptive finite element-singly implicit Runge-Kutta methods for parabolic differential equations

We describe an adaptive mesh refinement finite element method-of-lines procedure for solving one-dimensional parabolic partial differential equations. Solutions are calculated using Galerkin's method with a piecewise hierarchical polynomial basis in space and singly implicit Runge-Kutta (SIRK) methods in time. A modified SIRK formulation eliminates a linear systems solution that is required by the traditional SIRK formulation and leads to a new reduced-order interpolation formula. Stability and temporal error estimation techniques allow acceptance of approximate solutions at intermediate stages, yielding increased efficiency when solving partial differential equations. A priori energy estimates of the local discretization error are obtained for a nonlinear scalar problem. A posteriori estimates of local spatial discretization errors, obtained by order variation, are used with the a priori error estimates to control the adaptive mesh refinement strategy. Computational results suggest convergence of the a posteriori error estimate to the exact discretization error and verify the utility of the adaptive technique.This research was partially supported by the U.S. Air Force Office of Scientific Research, Air Force Systems Command, USAF, under Grant Number AFOSR-90-0194; the U.S. Army Research Office under Contract Number DAAL 03-91-G-0215; by the National Science Foundation under Grant Number CDA-8805910; and by a grant from the Committee on Research, Tulane University.     

High order explicit methods for parabolic equations

This paper discusses explicit embedded integration methods with large stability domains of order 3 and 4. The high order produces accurate results, the large stability domains allow some reasonable stiffness, the explicitness enables the method to treat very large problems, often space discretization of parabolic PDEs, and the embedded formulas permit an efficient stepsize control. The construction of these methods is achieved in two steps: firstly we compute stability polynomials of a given order with optimal stability domains, i.e., possessing a Chebyshev alternation; secondly we realize a corresponding explicit Runge-Kutta method with the help of the theory of composition methods. This work was supported by the Russian Fund Fundamental Researches and the Swiss National Science Foundation 20-43.314.95.     

Preconditioners for spectral element methods for elliptic and parabolic problems

In this paper we propose preconditioners for spectral element methods for elliptic and parabolic problems. These preconditioners are constructed using separation of variables and are easy to invert. Moreover they are spectrally equivalent to the quadratic forms which they are used to approximate.     

High-order finite element methods for time-fractional partial differential equations

The aim of this paper is to develop high-order methods for solving time-fractional partial differential equations. The proposed high-order method is based on high-order finite element method for space and finite difference method for time. Optimal convergence rate O((Δt)^2−α+N^−r) is proved for the (r−1)th-order finite element method (r≥2).     

Superconvergence of numerical solutions to weakly singular Volterra integro-differential equations

Summary We discuss the application of a class of spline collocation methods to first-order Volterra integro-differential equations (VIDEs) which contain a weakly singular kernel (t–s)^– with 0<<1. It will be shown that superconvergence properties may be obtained by using appropriate collocation parameters and graded meshes. The grading exponents of graded meshes used are not greater thanm (the polynomial degree) which is independent of . This is in contrast to the theories of spline collocation methods for Volterra (or Fredholm) integral equation of the second kind. Numerical examples are given to illustrate the theoretical results.     

Stability and accuracy of time-extrapolated ADI-FDTD methods for solving wave equations

Some recent work on the ADI-FDTD method for solving Maxwell's equations in 3-D have brought out the importance of extrapolation methods for the time stepping of wave equations. Such extrapolation methods have previously been used for the solution of ODEs. The present context (of wave equations) brings up two main questions which have not been addressed previously: (1) when will extrapolation in time of an unconditionally stable scheme for a wave equation again feature unconditional stability, and (2) how will the accuracy and computational efficiency depend on how frequently in time the extrapolations are carried out. We analyze these issues here.     

Locally linearized fractional step methods for nonlinear parabolic problems

This work deals with the efficient numerical solution of a class of nonlinear time-dependent reaction-diffusion equations. Via the method of lines approach, we first perform the spatial discretization of the original problem by applying a mimetic finite difference scheme. The system of ordinary differential equations arising from that process is then integrated in time with a linearly implicit fractional step method. For that purpose, we locally decompose the discrete nonlinear diffusion operator using suitable Taylor expansions and a domain decomposition splitting technique. The totally discrete scheme considers implicit time integrations for the linear terms while explicitly handling the nonlinear ones. As a result, the original problem is reduced to the solution of several linear systems per time step which can be trivially decomposed into a set of uncoupled parallelizable linear subsystems. The convergence of the proposed methods is illustrated by numerical experiments.     

Finite element methods for optimal control problems governed by integral equations and integro-differential equations

In this paper, we analyze finite-element Galerkin discretizations for a class of constrained optimal control problems that are governed by Fredholm integral or integro-differential equations. The analysis focuses on the derivation of a priori error estimates and a posteriori error estimators for the approximation schemes.Grants, communicated-by lines, or other notes about the article will be placed here between rules. Such notes are optional.     

Implicit-explicit multistep methods for quasilinear parabolic equations

Summary. Efficient combinations of implicit and explicit multistep methods for nonlinear parabolic equations were recently studied in [1]. In this note we present a refined analysis to allow more general nonlinearities. The abstract theory is applied to a quasilinear parabolic equation. Received March 10, 1997 / Revised version received March 2, 1998     

Superconvergence of collocation methods for a class of weakly singular Volterra integral equations

We discuss the application of spline collocation methods to a certain class of weakly singular Volterra integral equations. It will be shown that, by a special choice of the collocation parameters, superconvergence properties can be obtained if the exact solution satisfies certain conditions. This is in contrast with the theory of collocation methods for Abel type equations. Several numerical examples are given which illustrate the theoretical results.     

Contractivity of domain decomposition splitting methods for nonlinear parabolic problems

This work deals with the efficient numerical solution of nonlinear parabolic problems posed on a two-dimensional domain Ω. We consider a suitable decomposition of domain Ω and we construct a subordinate smooth partition of unity that we use to rewrite the original equation. Then, the combination of standard spatial discretizations with certain splitting time integrators gives rise to unconditionally contractive schemes. The efficiency of the resulting algorithms stems from the fact that the calculations required at each internal stage can be performed in parallel.     

Weighted finite difference methods for a class of space fractional partial differential equations with variable coefficients

A class of weighted finite difference methods (WFDMs) for solving a class of initial-boundary value problems of space fractional partial differential equations with variable coefficients is presented. Their stability and convergence properties are considered. It is proven that the WFDMs are unconditionally-stable for , and conditionally-stable for , here r is the weighting parameter subjected to 0≤r≤1. Some convergence results are given. These methods, problems and results generalize the corresponding methods, problems and results given in [7], [8] and [10]. Some numerical examples are provided to show the effectiveness of the methods with different weighting parameters.     

Collocation and iterated collocation methods for a class of weakly singular Volterra integral equations

We discuss the convergence properties of spline collocation and iterated collocation methods for a weakly singular Volterra integral equation associated with certain heat conduction problems. This work completes the previous studies of numerical methods for this type of equations with noncompact kernel. In particular, a global convergence result is obtained and it is shown that discrete superconvergence can be achieved with the iterated collocation if the exact solution belongs to some appropriate spaces. Some numerical examples illustrate the theoretical results.     

Solving parabolic and hyperbolic equations by the generalized finite difference method

Classical finite difference schemes are in wide use today for approximately solving partial differential equations of mathematical physics. An evolution of the method of finite differences has been the development of generalized finite difference (GFD) method, that can be applied to irregular grids of points.     

Unconditionally stable explicit methods for parabolic equations

Summary This paper discussesrational Runge-Kutta methods for stiff differential equations of high dimensions. These methods are explicit and in addition do not require the computation or storage of the Jacobian. A stability analysis (based onn-dimensional linear equations) is given. A second orderA ₀-stable method with embedded error control is constructed and numerical results of stiff problems originating from linear and nonlinear parabolic equations are presented.     

Error and extrapolation of a compact LOD method for parabolic differential equations

This paper is concerned with a compact locally one-dimensional (LOD) finite difference method for solving two-dimensional nonhomogeneous parabolic differential equations. An explicit error estimate for the finite difference solution is given in the discrete infinity norm. It is shown that the method has the accuracy of the second-order in time and the fourth-order in space with respect to the discrete infinity norm. A Richardson extrapolation algorithm is developed to make the final computed solution fourth-order accurate in both time and space when the time step equals the spatial mesh size. Numerical results demonstrate the accuracy and the high efficiency of the extrapolation algorithm.     

Multi-parameter error resolution for the collocation method of Volterra integral equations

In this paper, a multi-parameter error resolution technique is introduced and applied to the collocation method for Volterra integral equations. By using this technique, an approximation of higher accuracy is obtained by using a multi-processor in parallel. Additionally, a correction scheme for approximation of higher accuracy and a global superconvergence result are presented.     

Quadrature-difference methods for solving linear and nonlinear singular integro-differential equations

Here we propose and justify quadrature-difference methods for solving different kinds (linear, nonlinear and multidimensional) of periodic singular integro-differential equations.     

Spectral element methods for parabolic problems

A spectral element method for solving parabolic initial boundary value problems on smooth domains using parallel computers is presented in this paper. The space domain is divided into a number of shape regular quadrilaterals of size h and the time step k   is proportional to h²$h^2$. At each time step we minimize a functional which is the sum of the squares of the residuals in the partial differential equation, initial condition and boundary condition in different Sobolev norms and a term which measures the jump in the function and its derivatives across inter-element boundaries in certain Sobolev norms. The Sobolev spaces used are of different orders in space and time. We can define a preconditioner for the minimization problem which allows the problem to decouple. Error estimates are obtained for both the h and p versions of this method.