The spectral methods for parabolic Volterra integro-differential equations

In this paper we study the numerical solutions to parabolic Volterra integro-differential equations in one-dimensional bounded and unbounded spatial domains. In a bounded domain, the given parabolic Volterra integro-differential equation is converted to two equivalent equations. Then, a Legendre-collocation method is used to solve them and finally a linear algebraic system is obtained. For an unbounded case, we use the algebraic mapping to transfer the problem on a bounded domain and then apply the same presented approach for the bounded domain. In both cases, some numerical examples are presented to illustrate the efficiency and accuracy of the proposed method.     

Improved radial basis function methods for multi-dimensional option pricing

In this paper, we have derived a radial basis function (RBF) based method for the pricing of financial contracts by solving the Black–Scholes partial differential equation. As an example of a financial contract that can be priced with this method we have chosen the multi-dimensional European basket call option. We have shown numerically that our scheme is second-order accurate in time and spectrally accurate in space for constant shape parameter. For other non-optimal choices of shape parameter values, the resulting convergence rate is algebraic. We propose an adapted node point placement that improves the accuracy compared with a uniform distribution. Compared with an adaptive finite difference method, the RBF method is 20–40 times faster in one and two space dimensions and has approximately the same memory requirements.     

Stability estimates based on numerical ranges with an application to a spectral method

This paper is concerned with the stability of numerical processes for solving initial value problems. We present a stability result which is related to a well-known theorem by von Neumann, but the requirements to be satisfied are less severe and easier to verify.As an illustration we consider a simple convection-diffusion equation. For the spatial discretization we use a spectral collocation method (based on so-called Legendre-Gauss-Lobatto points). We show that the fully discretized numerical process is stable, provided that the temporal step size is bounded by a constant depending only on the convection-diffusion equation, the number of collocation points and the time-stepping method under consideration.This research has been supported by the Netherlands Organization for Scientific Research (N.W.O.).     

GUARANTEED BOUNDS FOR THE SOLUTION OF THE WAVE EQUATION

Abstract

Numerical solution of the wave equation in the form of close lower and upper bounds provides a secure a posteriori error estimate that can be used for efficient accuracy control. The method considered in this paper uses some monotone properties of the differential operator in the wave equation to construct bounds for the solution in the form of trigonometric polynomials of x. Aspects of the numerical implementation, the accuracy of the computed bounds and some numerical examples are discussed.     

Least-squares spectral element method for non-linear hyperbolic differential equations

The least-squares spectral element method has been applied to the one-dimensional inviscid Burgers equation which allows for discontinuous solutions. In order to achieve high order accuracy both in space and in time a space–time formulation has been applied. The Burgers equation has been discretized in three different ways: a non-conservative formulation, a conservative system with two variables and two equations: one first order linear PDE and one linearized algebraic equation, and finally a variant on this conservative formulation applied to a direct minimization with a QR-decomposition at elemental level. For all three formulations an h/p-convergence study has been performed and the results are discussed in this paper.     

A qualocation method for Burgers’ equation

In this paper, a qualocation method for the one-dimensional Burgers’ equation is proposed. A semidiscrete scheme along with optimal error estimates is discussed. Results of a numerical experiment performed support the theoretical results.     

Orthogonal spline collocation methods for Schr\"{o}dinger-type equations in one space variable

Summary. We examine the use of orthogonal spline collocation for the semi-discreti\-za\-tion of the cubic Schr\"{o}dinger equation and the two-dimensional parabolic equation of Tappert. In each case, an optimal order estimate of the error in the semidiscrete approximation is derived. For the cubic Schr\"{o}dinger equation, we present the results of numerical experiments in which the integration in time is performed using a routine from a software library. Received February 14, 1992 / Revised version received December 29, 1992     

Vibration analysis of Kirchhoff plates by the Morley element method

Vibration analysis of Kirchhoff plates is of great importance in many engineering fields. The semi-discrete and the fully discrete Morley element methods are proposed to solve such a problem, which are effective even when the region of interest is irregular. The rigorous error estimates in the energy norm for both methods are established. Some reasonable approaches to choosing the initial functions are given to keep the good convergence rate of the fully discrete method. A number of numerical results are provided to illustrate the computational performance of the method in this paper.     

A Gautschi-type method for oscillatory second-order differential equations

Summary. We study a numerical method for second-order differential equations in which high-frequency oscillations are generated by a linear part. For example, semilinear wave equations are of this type. The numerical scheme is based on the requirement that it solves linear problems with constant inhomogeneity exactly. We prove that the method admits second-order error bounds which are independent of the product of the step size with the frequencies. Our analysis also provides new insight into the m ollified impulse method of García-Archilla, Sanz-Serna, and Skeel. We include results of numerical experiments with the sine-Gordon equation. Received January 21, 1998 / Published online: June 29, 1999     

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.     

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.     

The Laguerre spectral method for solving Neumann boundary value problems

In this paper, we propose a Laguerre spectral method for solving Neumann boundary value problems. This approach differs from the classical spectral method in that the homogeneous boundary condition is satisfied exactly. Moreover, a tridiagonal matrix is employed, instead of the full stiffness matrix encountered in the classical variational formulation of such problems. For analyzing the numerical errors, some basic results on Laguerre approximations are established. The convergence is proved. The numerical results demonstrate the efficiency of this approach.     

A stabilized semi-implicit Galerkin scheme for Navier-Stokes equations

By introducing a time relaxation term for the time derivative of higher frequency components, we proposed a stabilized semi-implicit Galerkin scheme for evolutionary Navier-Stokes equations in this paper. Analysis shows that such a scheme has weaker stability conditions than that of a classical semi-implicit Galerkin scheme and, when a suitable relaxation parameter σ is chosen, it generates an approximate solution with the same accuracy as the classical one. That means the proposed scheme might use a larger time step to generate a bounded approximate solution. Thus it is more suitable for long time simulations.     

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.     

Error analysis and applications of the Fourier-Galerkin Runge-Kutta schemes for high-order stiff PDEs

An integrating factor mixed with Runge-Kutta technique is a time integration method that can be efficiently combined with spatial spectral approximations to provide a very high resolution to the smooth solutions of some linear and nonlinear partial differential equations. In this paper, the novel hybrid Fourier-Galerkin Runge-Kutta scheme, with the aid of an integrating factor, is proposed to solve nonlinear high-order stiff PDEs. Error analysis and properties of the scheme are provided. Application to the approximate solution of the nonlinear stiff Korteweg-de Vries (the 3rd order PDE, dispersive equation), Kuramoto-Sivashinsky (the 4th order PDE, dissipative equation) and Kawahara (the 5th order PDE) equations are presented. Comparisons are made between this proposed scheme and the competing method given by Kassam and Trefethen. It is found that for KdV, KS and Kawahara equations, the proposed method is the best.     

Orthogonal spline collocation methods for biharmonic problems

Orthogonal spline collocation methods are formulated and analyzed for the solution of certain biharmonic problems in the unit square. Particular attention is given to the Dirichlet biharmonic problem which is solved using capacitance matrix techniques. Received November 11, 1996     

A time-splitting spectral method for coupled Gross–Pitaevskii equations with applications to rotating Bose–Einstein condensates

We propose a time-splitting spectral method for the coupled Gross–Pitaevskii equations, which describe the dynamics of rotating two-component Bose–Einstein condensates at a very low temperature. The new numerical method is explicit, unconditionally stable, time reversible, time transverse invariant, and of spectral accuracy in space and second-order accuracy in time. Moreover, it conserves the position densities in the discretized level. Numerical applications on studying the generation of topological modes and the vortex lattice dynamics for the rotating two-component Bose–Einstein condensates are presented in detail.     

Spectral methods for initial boundary value problems with nonsmooth data

Summary. In this paper we consider hyperbolic initial boundary value problems with nonsmooth data. We show that if we extend the time domain to minus infinity, replace the initial condition by a growth condition at minus infinity and then solve the problem using a filtered version of the data by the Galerkin-Collocation method using Laguerre polynomials in time and Legendre polynomials in space, then we can recover pointwise values with spectral accuracy, provided that the actual solution is piecewise smooth. For this we have to perform a local smoothing of the computed solution. Received August 1, 1995 / Revised version received August 19, 1997     

Numerical solution of the 3D time dependent Schrödinger equation in spherical coordinates: Spectral basis and effects of split-operator technique

We study a numerical solution of the multi-dimensional time dependent Schrödinger equation using a split-operator technique for time stepping and a spectral approximation in the spatial coordinates. We are particularly interested in systems with near spherical symmetries. One expects these problems to be most efficiently computed in spherical coordinates as a coarse grain discretization should be sufficient in the angular directions. However, in this coordinate system the standard Fourier basis does not provide a good basis set in the radial direction. Here, we suggest an alternative basis set based on Chebyshev polynomials and a variable transformation.