THE STABILITY OF A VARIATIONAL METHOD FOR THE SOLUTION OF A LINEAR HYPERBOLIC INITIAL BOUNDARY VALUE PROBLEM

Abstract

The numerical stability of a variational method that is used to obtain the solution of a one space dimension wave equation with initial and boundary conditions is analyzed. The phase speed and group velocity of the numerical solution are also investigated with respect to that of the exact solution.     

A finite volume scheme for nonlinear parabolic equations derived from one-dimensional local dirichlet problems

In this paper, we propose a new method to compute the numerical flux of a finite volume scheme, used for the approximation of the solution of the nonlinear partial differential equation u_t+div(qf(u))−ΔΦ(u)=0 in a 1D, 2D or 3D domain. The function Φ is supposed to be strictly increasing, but some values s such that Φ′(s)=0 can exist. The method is based on the solution, at each interface between two control volumes, of the nonlinear elliptic two point boundary value problem (qf(υ)+(Φ(υ))′)′=0 with Dirichlet boundary conditions given by the values of the discrete approximation in both control volumes. We prove the existence of a solution to this two point boundary value problem. We show that the expression for the numerical flux can be yielded without referring to this solution. Furthermore, we prove that the so designed finite volume scheme has the expected stability properties and that its solution converges to the weak solution of the continuous problem. Numerical results show the increase of accuracy due to the use of this scheme, compared to some other schemes.     

Approximation with frequency filter for backward parabolic equations

Summary. We propose a numerical method for the initial (and boundary) value problem for the equation of the form where A is an unbounded, selfadjoint operator with negative spectrum. Roundoff errors in the numerical solution of such problem may generate a parasite term growing very quickly with time. To eliminate this parasite term, we apply a special finite difference equation with r free parameters. Similar ideas may be useful also for another numerically difficult differential problems. Received October 6, 1997 / revised version received November 26, 1998 / Published online October 16, 2000     

Marching schemes for inverse acoustic scattering problems

Summary For the numerical solution of inverse Helmholtz problems the boundary value problem for a Helmholtz equation with spatially variable wave number has to be solved repeatedly. For large wave numbers this is a challenge. In the paper we reformulate the inverse problem as an initial value problem, and describe a marching scheme for the numerical computation that needs only n² log n operations on an n × n grid. We derive stability and error estimates for the marching scheme. We show that the marching solution is close to the low-pass filtered true solution. We present numerical examples that demonstrate the efficacy of the marching scheme.     

Convergence of the Lambert-McLeod trajectory solver and of the celf method

Summary A trajectory problem is an initial value problemd y/dt=f(y),y(0)= where the interest lies in obtaining the curve traced by the solution (the trajectory), rather than in finding the actual correspondanc between values of the parametert and points on that curve. We prove the convergence of the Lambert-McLeod scheme for the numerical integration of trajectory problems. We also study the CELF method, an explicit procedure for the integration in time of semidiscretizations of PDEs which has some useful conservation properties. The proofs rely on the concept of restricted stability introduced by Stetter. In order to show the convergence of the methods, an idea of Strang is also employed, whereby the numerical solution is compared with a suitable perturbation of the theoretical solution, rather than with the theoretical solution itself.     

Numerical analysis for a conservative difference scheme to solve the Schrödinger-Boussinesq equation

In this paper, we present a finite difference scheme for the solution of an initial-boundary value problem of the Schrödinger-Boussinesq equation. The scheme is fully implicit and conserves two invariable quantities of the system. We investigate the existence of the solution for the scheme, give computational process for the numerical solution and prove convergence of iteration method by which a nonlinear algebra system for unknown Vⁿ⁺¹ is solved. On the basis of a priori estimates for a numerical solution, the uniqueness, convergence and stability for the difference solution is discussed. Numerical experiments verify the accuracy of our method.     

Error estimate for finite volume scheme

We study the convergence of a finite volume scheme for the linear advection equation with a Lipschitz divergence-free speed in R ^d . We prove a h ^1/2-error estimate in the L ^∞(0,t;L ¹)-norm for BV data. This result was expected from numerical experiments and is optimal.     

Runge–Kutta methods and viscous wave equations

We study the numerical time integration of a class of viscous wave equations by means of Runge–Kutta methods. The viscous wave equation is an extension of the standard second-order wave equation including advection–diffusion terms differentiated in time. The viscous wave equation can be very stiff so that for time integration traditional explicit methods are no longer efficient. A-Stable Runge–Kutta methods are then very good candidates for time integration, in particular diagonally implicit ones. Special attention is paid to the question how the A-Stability property can be translated to this non-standard class of viscous wave equations.      

Computation and analysis for a constrained entropy optimization problem in finance

In [T. Coleman, C. He, Y. Li, Calibrating volatility function bounds for an uncertain volatility model, Journal of Computational Finance (2006) (submitted for publication)], an entropy minimization formulation has been proposed to calibrate an uncertain volatility option pricing model (UVM) from market bid and ask prices. To avoid potential infeasibility due to numerical error, a quadratic penalty function approach is applied. In this paper, we show that the solution to the quadratic penalty problem can be obtained by minimizing an objective function which can be evaluated via solving a Hamilton–Jacobian–Bellman (HJB) equation. We prove that the implicit finite difference solution of this HJB equation converges to its viscosity solution. In addition, we provide computational examples illustrating accuracy of calibration.     

Error growth analysis via stability regions for discretizations of initial value problems

This paper deals with numerical methods for the solution of linear initial value problems. Two main theorems are presented on the stability of these methods. Both theorems give conditions guaranteeing a mild error growth, for one-step methods characterized by a rational function ϕ(z). The conditions are related to the stability regionS={z:z∈ℂ with |ϕ(z)|≤1}, and can be viewed as variants to the resolvent condition occurring in the reputed Kreiss matrix theorem. Stability estimates are presented in terms of the number of time stepsn and the dimensions of the space. The first theorem gives a stability estimate which implies that errors in the numerical process cannot grow faster than linearly withs orn. It improves previous results in the literature where various restrictions were imposed onS and ϕ(z), including ϕ′(z)≠0 forz∈σS andS be bounded. The new theorem is not subject to any of these restrictions. The second theorem gives a sharper stability result under additional assumptions regarding the differential equation. This result implies that errors cannot grow faster thann ^β, with fixed β<1. The theory is illustrated in the numerical solution of an initial-boundary value problem for a partial differential equation, where the error growth is measured in the maximum norm.     

A second-order accurate numerical method for a fractional wave equation

We study a generalized Crank–Nicolson scheme for the time discretization of a fractional wave equation, in combination with a space discretization by linear finite elements. The scheme uses a non-uniform grid in time to compensate for the singular behaviour of the exact solution at t = 0. With appropriate assumptions on the data and assuming that the spatial domain is convex or smooth, we show that the error is of order k ² + h ², where k and h are the parameters for the time and space meshes, respectively.     

A GLOBAL-LOCAL FINITE ELEMENT METHOD IN SPACE-TIME FOR A HYPERBOLIC PROBLEM

ABSTRACT

A finite element method for solving the wave equation with couples boundary conditions is presented. In this approach finite elements are applied globally with respect to space and simultaneously but locally with respect to time. This gives rise to a single-step method in time. The method is a practical and economic one and the numerical results obtained compare favourably with the available analytic solution.     

A normal compliance contact problem in viscoelasticity: An a posteriori error analysis and computational experiments

In this work, the numerical approximation of a viscoelastic contact problem is studied. The classical Kelvin-Voigt constitutive law is employed, and contact is assumed with a deformable obstacle and modelled using the normal compliance condition. The variational formulation leads to a nonlinear parabolic variational equation. An existence and uniqueness result is recalled. Then, a fully discrete scheme is introduced, by using the finite element method to approximate the spatial variable and the implicit Euler scheme to discretize time derivatives. A priori error estimates recently proved for this problem are recalled. Then, an a posteriori error analysis is provided, extending some preliminary results obtained in the study of the heat equation and other parabolic equations. Upper and lower error bounds are proved. Finally, some numerical experiments are presented to demonstrate the accuracy and the numerical behaviour of the error estimates.     

High-order one-way model nesting in dispersive non-uniform media

Global-regional model interaction is considered for two-dimensional linear time dependent waves in a dispersive non-uniform medium with a continuously varying wave speed. The setup, which is sometimes called ‘one-way nesting,’ arises in Numerical Weather Prediction (NWP) as well as in other fields concerning waves in very large domains. The Carpenter scheme for this type of problem is revisited, in the context of the dispersive wave equation with a variable wave speed. The original Carpenter scheme is based on the Sommerfeld radiation operator, and thus is associated with low-order accuracy. By replacing the Sommerfeld operator with the high-order Hagstrom-Warburton absorbing operator, a modified Carpenter open boundary condition emerges which possesses high-order accuracy. This is demonstrated via a numerical example in a wave guide with a wave speed which varies linearly in the cross section.     

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     

Solving the generalized regularized long wave equation on the basis of a reproducing kernel space

On the basis of a reproducing kernel space, an iterative algorithm for solving the generalized regularized long wave equation is presented. The analytical solution in the reproducing kernel space is shown in a series form and the approximate solution u_n is constructed by truncating the series to n terms. The convergence of u_n to the analytical solution is also proved. Results obtained by the proposed method imply that it can be considered as a simple and accurate method for solving such evolution equations.     

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.     

A posteriori error analysis for numerical approximations of Friedrichs systems

The global error of numerical approximations for symmetric positive systems in the sense of Friedrichs is decomposed into a locally created part and a propagating component. Residual-based two-sided local a posteriori error bounds are derived for the locally created part of the global error. These suggest taking the -norm as well as weaker, dual norms of the computable residual as local error indicators. The dual graph norm of the residual is further bounded from above and below in terms of the norm of where h is the local mesh size. The theoretical results are illustrated by a series of numerical experiments. Received January 10, 1997 / Revised version received March 5, 1998     

Boundary controllability of a linear semi-discrete 1-D wave equation derived from a mixed finite element method

In this article one discusses the controllability of a semi-discrete system obtained by discretizing in space the linear 1-D wave equation with a boundary control at one extremity. It is known that the semi-discrete models obtained with finite difference or the classical finite element method are not uniformly controllable as the discretization parameter h goes to zero (see [8]). Here we introduce a new semi-discrete model based on a mixed finite element method with two different basis functions for the position and velocity. We show that the controls obtained with these semi-discrete systems can be chosen uniformly bounded in L²(0,T) and in such a way that they converge to the HUM control of the continuous wave equation, i.e. the minimal L²-norm control. We illustrate the mathematical results with several numerical experiments. Supported by Grant BFM 2002-03345 of MCYT (Spain) and the TMR projects of the EU ``Homogenization and Multiple Scales" and ``New materials, adaptive systems and their nonlinearities: modelling, control and numerical simulations". Partially Supported by Grant BFM 2002-03345 of MCYT (Spain), Grant 17 of Egide-Brancusi Program and Grant 80/2005 of CNCSIS (Romania).