High order finite difference numerical methods for time-dependent convection-dominated problems

DOI10.1016/J.APNUM.2005.04.038MaRDI QIDQ2573277zbMATH OpenOpenAlexFDO

Authors Maria Antonietta Pirozzi, M. Morandi Cecchi

Publication date 7 November 2005

Published in Applied Numerical Mathematics (Search for Journal in Brave)

Full work available at URL https://doi.org/10.1016/j.apnum.2005.04.038

convection-diffusion equation Crank-Nicolson method error bounds finite difference method initial value problems multidimensional scheme splitting method stability viscous Burgers' equation

Mathematics Subject Classification ID

Hyperbolic conservation laws (35L65) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs (65M12) Error bounds for initial value and initial-boundary value problems involving PDEs (65M15)

Review text

The paper continues the investigation of specific implicit finite difference methods applied to initial value problems of convection-diffusion equations. \textit{M. A. Pirozzi} [High-accuracy finite difference methods for systems of conservation laws, II, Tech. Rep. No. 6, Instituto di Matematica, Fisica e Applicazioni, Universita di Napoli `Parthenope' (2002)] has constructed a parameter-dependent implicit discretisation of the linear convection-diffusion equation, which represents a parabolic PDE in case of one space dimension. The investigation of the corresponding modified equation allows to set the parameters such that high order schemes are achieved for the convection part in space. Furthermore, the parameters can be used to reduce dispersion and dissipative behaviour in the approximations. Based on a linearisation technique, the method can be used for solving nonlinear convection-diffusion equations, too. In the paper at hand, this class of implicit finite difference methods is generalised to two and three space dimensions. A splitting scheme arises, where just subsequent tridiagonal linear systems with the complexity of one space dimension have to be solved. Nevertheless, the parameters can be set such that the local truncation error in the convection part exhibits order two or three. Furthermore, the choices produce linear systems with strictly diagonal dominant matrices. The criterion of von Neumann does not restrict the size of time steps, which indicates good stability properties. In the main part of the paper, the authors present numerical simulations, which verify the resulting orders and indicate the magnitude of errors in amplitude and phase for corresponding approximations. However, the used equations are just the one- and two-dimensional linear convection-diffusion equation and the one-dimensional viscous Burgers' equation. Several choices for the parameters in the schemes and for the Courant number, which couples space and time step size, are used. Since only the convection part exhibits a high order approximation, the diffusion coefficient is reduced according to the space step size, where a hyperbolic partial differential equation arises in the limit. Thus the high order approximation is just given in the case of dominating convection, i.e., the diffusion part does not influence the complete error significantly. The results for the one-dimensional linear equation with a certain choice of parameters show that the resolution of amplitude and phase is as good as for the Crank-Nicolson mass operator scheme. Moreover, tests with discontinuous initial data determine sharp fronts correctly without oscillations or too much dissipation at later times.

Recommendations

Cites work

Cited in

(15)

This page was built for publication: High order finite difference numerical methods for time-dependent convection-dominated problems

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2573277)