Space–time discontinuous Galerkin finite element method for shallow water flows

A space–time discontinuous Galerkin (DG) finite element method is presented for the shallow water equations over varying bottom topography. The method results in nonlinear equations per element, which are solved locally by establishing the element communication with a numerical HLLC flux. To deal with spurious oscillations around discontinuities, we employ a dissipation operator only around discontinuities using Krivodonova's discontinuity detector. The numerical scheme is verified by comparing numerical and exact solutions, and validated against a laboratory experiment involving flow through a contraction. We conclude that the method is second order accurate in both space and time for linear polynomials.     

Relaxation approximation to bed-load sediment transport

In this work we propose and apply a numerical method based on finite volume relaxation approximation for computing the bed-load sediment transport in shallow water flows, in one and two space dimensions. The water flow is modeled by the well-known nonlinear shallow water equations which are coupled with a bed updating equation. Using a relaxation approximation, the nonlinear set of equations (and for two different formulations) is transformed to a semilinear diagonalizable problem with linear characteristic variables. A second order MUSCL-TVD method is used for the advection stage while an implicit–explicit Runge–Kutta scheme solves the relaxation stage. The main advantages of this approach are that neither Riemann problem solvers nor nonlinear iterations are required during the solution process. For the two different formulations, the applicability and effectiveness of the presented scheme is verified by comparing numerical results obtained for several benchmark test problems.     

Runge–Kutta convolution quadrature methods for well-posed equations with memory

Runge–Kutta based convolution quadrature methods for abstract, well-posed, linear, and homogeneous Volterra equations, non necessarily of sectorial type, are developed. A general representation of the numerical solution in terms of the continuous one is given. The error and stability analysis is based on this representation, which, for the particular case of the backward Euler method, also shows that the numerical solution inherits some interesting qualitative properties, such as positivity, of the exact solution. Numerical illustrations are provided.     

Construction of a multirate RODAS method for stiff ODEs

Multirate time stepping is a numerical technique for efficiently solving large-scale ordinary differential equations (ODEs) with widely different time scales localized over the components. This technique enables one to use large time steps for slowly varying components, and small steps for rapidly varying ones. Multirate methods found in the literature are normally of low order, one or two. Focusing on stiff ODEs, in this paper we discuss the construction of a multirate method based on the fourth-order RODAS method. Special attention is paid to the treatment of the refinement interfaces with regard to the choice of the interpolant and the occurrence of order reduction. For stiff, linear systems containing a stiff source term, we propose modifications for the treatment of the source term which overcome order reduction originating from such terms and which we can implement in our multirate method.     

Conserved quantities of some Hamiltonian wave equations after full discretization

Hamiltonian PDEs have some invariant quantities, which would be good to conserve with the numerical integration. In this paper, we concentrate on the nonlinear wave and Schrödinger equations. Under hypotheses of regularity and periodicity, we study how a symmetric space discretization makes that the space discretized system also has some invariants or `nearly' invariants which well approximate the continuous ones. We conjecture some facts which would explain the good numerical approximation of them after time integration when using symplectic Runge-Kutta methods or symmetric linear multistep methods for second-order systems.     

A note on the Gautschi-type method for oscillatory second-order differential equations

The Gautschi-type method has been proposed by Hochbruck and Lubich for oscillatory second-order differential equations. They conjecture that this method allows for a uniform error bound independent of the size of the system. The conjecture is proved in this note.     

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.     

Dimension splitting for evolution equations

In this paper, we are concerned with splitting methods for the time integration of abstract evolution equations. We introduce an analytic framework which allows us to prove optimal convergence orders for various splitting methods, including the Lie and Peaceman–Rachford splittings. Our setting is applicable for a wide variety of linear equations and their dimension splittings. In particular, we analyze parabolic problems with Dirichlet boundary conditions, as well as degenerate equations on bounded domains. We further illustrate our theoretical results with a set of numerical experiments. This work was supported by the Austrian Science Fund under grant M961-N13.     

A two-grid method for finite volume element approximations of second-order nonlinear hyperbolic equations

The two-grid method is studied for solving a two-dimensional second-order nonlinear hyperbolic equation using finite volume element method. The method is based on two different finite element spaces defined on one coarse grid with grid size H and one fine grid with grid size h, respectively. The nonsymmetric and nonlinear iterations are only executed on the coarse grid and the fine grid solution can be obtained in a single symmetric and linear step. It is proved that the coarse grid can be much coarser than the fine grid. A prior error estimate in the H¹-norm is proved to be O(h+H³|lnH|) for the two-grid semidiscrete finite volume element method. With these proposed techniques, solving such a large class of second-order nonlinear hyperbolic equations will not be much more difficult than solving one single linearized equation. Finally, a numerical example is presented to validate the usefulness and efficiency of the method.     

Efficient computation of characteristic roots of delay differential equations using LMS methods

We aim at the efficient computation of the rightmost, stability-determining characteristic roots of a system of delay differential equations. The approach we use is based on the discretization of the time integration operator by a linear multistep (LMS) method. The size of the resulting algebraic eigenvalue problem is inversely proportional to the steplength. We summarize theoretical results on the location and numerical preservation of roots. Furthermore, we select nonstandard LMS methods, which are better suited for our purpose. We present a new procedure that aims at computing efficiently and accurately all roots in any right half-plane. The performance of the new procedure is demonstrated for small- and large-scale systems of delay differential equations.     

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.      

Maximum norm error estimates of efficient difference schemes for second-order wave equations

The three-level explicit scheme is efficient for numerical approximation of the second-order wave equations. By employing a fourth-order accurate scheme to approximate the solution at first time level, it is shown that the discrete solution is conditionally convergent in the maximum norm with the convergence order of two. Since the asymptotic expansion of the difference solution consists of odd powers of the mesh parameters (time step and spacings), an unusual Richardson extrapolation formula is needed in promoting the second-order solution to fourth-order accuracy. Extensions of our technique to the classical ADI scheme also yield the maximum norm error estimate of the discrete solution and its extrapolation. Numerical experiments are presented to support our theoretical results.     

Non-separable splines and numerical computation of evolution equations by the Galerkin methods

We construct new non-separable splines and we apply the spline sampling approximation to the computation of numerical solutions of evolution equations. The non-separable splines are basis functions which give a fine sampling approximation which enables us to compute numerical solutions by means of the method of lines combined with the Galerkin method. To demonstrate our approach we compute numerical solutions of the Burgers equation and the Kadomtsev–Petviashvili equation.     

Septic spline solutions of sixth-order boundary value problems

Septic spline is used for the numerical solution of the sixth-order linear, special case boundary value problem. End conditions for the definition of septic spline are derived, consistent with the sixth-order boundary value problem. The algorithm developed approximates the solution and their higher-order derivatives. The method has also been proved to be second-order convergent. Three examples are considered for the numerical illustrations of the method developed. The method developed in this paper is also compared with that developed in [M. El-Gamel, J.R. Cannon, J. Latour, A.I. Zayed, Sinc-Galerkin method for solving linear sixth order boundary-value problems, Mathematics of Computation 73, 247 (2003) 1325–1343], as well and is observed to be better.     

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.     

A class of explicit multistep exponential integrators for semilinear problems

A class of explicit multistep exponential methods for abstract semilinear equations is introduced and analyzed. It is shown that the k-step method achieves order k, for appropriate starting values, which can be computed by auxiliary routines or by one strategy proposed in the paper. Together with some implementation issues, numerical illustrations are also provided.     

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.     

A dual graph-norm refinement indicator for finite volume approximations of the Euler equations

Summary. We generalise and apply a refinement indicator of the type originally designed by Mackenzie, Süli and Warnecke in [15] and [16] for linear Friedrichs systems to the Euler equations of inviscid, compressible fluid flow. The Euler equations are symmetrized by means of entropy variables and locally linearized about a constant state to obtain a symmetric hyperbolic system to which an a posteriori error analysis of the type introduced in [15] can be applied. We discuss the details of the implementation of the refinement indicator into the DLR--Code which is based on a finite volume method of box type on an unstructured grid and present numerical results. Received May 15, 1995 / Revised version received April 17, 1996     

A second order backward difference method with variable steps for a parabolic problem

The numerical solution of a parabolic problem is studied. The equation is discretized in time by means of a second order two step backward difference method with variable time step. A stability result is proved by the energy method under certain restrictions on the ratios of successive time steps. Error estimates are derived and applications are given to homogenous equations with initial data of low regularity.     

A numerical verification method for a periodic solution of a delay differential equation

We describe a numerical method with guaranteed accuracy to enclose a periodic solution for a system of delay differential equations. Using a certain system of equations corresponding to the original system, we derive sufficient conditions for the existence of the solution, the satisfaction of which can be verified computationally. We describe the verification procedure in detail and give a numerical example.