In this work, high order splitting methods have been used for calculating the numerical solutions of the Burgers' equation in one space dimension with periodic and Dirichlet boundary conditions. However, splitting methods with real coefficients of order higher than two necessarily have negative coefficients and can not be used for time-irreversible systems, such as Burgers equations, due to the time-irreversibility of the Laplacian operator. Therefore, the splitting methods with complex coefficients and extrapolation methods with real and positive coefficients have been employed. If we consider the system as the perturbation of an exactly solvable problem(or can be easily approximated numerically), it is possible to employ highly efficient methods to approximate Burgers' equation. The numerical results show that the methods with complex time steps having one set of coefficients real and positive, say a i ∈ R + and b i ∈ C + , and high order extrapolation methods derived from a lower order splitting method produce very accurate solutions of the Burgers' equation.
Abstract:The purpose of this work is to introduce a new kind of finite difference formulation inspired from Fourier analysis, for reaction-diffusion equations. Compared to classical schemes, the proposed scheme is much more accurate and has interesting stability properties. Convergence properties and stability of the scheme are discussed. Numerical examples are provided to show better performance of the method, compared with other existing methods in the literature.
Purpose
The purpose of this study is to propose a non-classical method to obtain efficient and accurate numerical solutions of the advection–diffusion–reaction equations.
Design/methodology/approach
Unlike conventional numerical methods, this study proposes a numerical scheme using outer Newton iteration applied to a time-dependent PDE. The linearized time dependent PDE is discretized by trapezoidal rule, which is second order in time, and by spline-based finite difference method of fourth order in space.
Findings
Using the proposed technique, even when relatively large time step sizes are used in computations, the efficiency of the proposed procedure is very clear for the numerical examples in comparison with the existing classical methods.
Originality/value
This study, unlike these classical methods, proposes an alternative approach based on linearizing the nonlinear problem at first, and then discretizing it by an appropriate scheme. This technique helps to avoid considering the convergence issues of Newton iteration applied to nonlinear algebraic system containing many unknowns at each time step if an implicit method is used in time discretization. The linearized PDE can be solved by implicit time integrator, which enables the use of large time step size.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.