High-order compact boundary value method for the solution of unsteady convection–diffusion problems
Introduction
Consider the two-dimensional unsteady convection–diffusion equationwith initial conditionand the boundary conditionswhere is a transported (advected and diffused) scalar variable, and are arbitrary constants which show the speed of convection and the diffusion coefficients, i.e. and are positive constants.
For a positive integer n let denote the step size of spatial derivatives and for step size of temporal derivative. So we define
This equation may be seen in computational hydraulics and fluid dynamics to model convection–diffusion of quantities such as mass, heat, energy, vorticity, etc. The applications of Eq. (1.1) have been stated in [7], [8], [18].
Various finite difference schemes have been proposed to solve convection–diffusion problems approximately. For solving steady convection–diffusion equations, some compact finite difference approximations have been developed in the last two decay [12], [20], [16]. The authors of [17] proposed a nine-point high-order compact implicit scheme for Eq. (1.1) which is third-order accurate in space and second-order accurate in time, and has a large zone of stability. Spotz and Carey [21] presented an extension of higher order compact difference techniques for steady-state [20] to the time-dependent problems. These methods are conditionally stable and have fourth-order accuracy in space and second order or lower in time component. Cecchi and Pirozzi [5] presented a family of unconditionally stable finite difference methods for Eq. (1.1) and the possible maximum order of accuracy is . Other new and high-order formulas for solving Eq. (1.1) have been developed by some authors [14], [15], [23]. Kalita et al. [14] proposed some schemes for solving Eq. (1.1) which are based on high-order compact scheme and weighted time discretization. Their schemes are second or lower order accurate in time and fourth-order accurate in space. Karaa and Zhang [15] presented a high-order alternating direction implicit (ADI) method for solving Eq. (1.1) which is fourth and second order in space and time respectively and is unconditionally stable. More research works can be found in Refs. [10], [23], [25], [26], [27]. In this paper, we introduce an unconditionally stable scheme for Eq. (1.1) which is fourth-order accurate in both space and time variables. We apply a compact finite difference scheme of fourth order for discretizing spatial derivatives and boundary value methods for the time integration of Eq. (1.1). Boundary value methods (BVMs) are unconditionally stable and are high-accuracy schemes for solving ordinary differential equations based on the linear multi-step formulas [2], [3], [4]. Unlike Runge–Kutta or other initial value methods (IVM), BVMs achieve the advantage of both good stability and high-order accuracy [1], [2], [3], [4], [9], [11].
The outline of this paper is as follows: in Section 2, we present the fourth-order compact finite difference approximation of one and two-dimensional steady-state convection–diffusion equations. In Section 3 we briefly introduce the boundary value methods and present a fourth-order boundary value method for solving IVPs. In Section 4, we apply the fourth-order compact difference scheme for discretizing spatial derivatives and the fourth-order boundary value method for the resulted linear system of ordinary differential equations. Also in this section we propose a scheme for the solution of resulted large linear system of equations. In Section 5 we describe the numerical experiments of solving one and two-dimensional unsteady convection–diffusion equations with the method of this article for four problems and compare the numerical results with some well-known numerical methods. The concluding remarks in drawn in Section 6. Some references are given for possible extension of the approach [28], [29], [30].
Section snippets
Compact finite difference schemes
The basic approach for high-order compact difference methods is to introduce the standard compact difference approximations to the differential equations and then by repeated differentiation and associated compact differencing, a new high-order compact scheme will be developed that incorporates the effect of the leading truncation error terms in the standard method [20]. Recently due to the high order, compactness and high resolution, we have seen increasing population for high-order compact
The boundary value methods
Boundary value methods are the recent classes of initial value solvers which can be interpreted as a generalization of the linear multi-step methods (LMMs) [2], [4]. In comparison with the other initial value solvers, BVMs additional to having high-order accuracy, are unconditionally stable methods.
Consider the following initial value problem
After discretizing the (3.1) with a k-step linear multi-step method, the k-conditions are needed to obtain the discrete
Compact boundary value methods
In this section we will combine high-order compact difference schemes presented in Section 2 with boundary value methods to obtain a high-order method for solving the two-dimensional unsteady convection–diffusion Eq. (1.1).
At first we discretize partial differential Eq. (1.1) in space with previous high-order methods to obtain a system of ordinary differential equations with unknown function at each spatial grid point. Then we will apply the boundary value methods for solving the resulted
Numerical experiments
In this section we present the numerical results of the new method on several problems. We tested the accuracy and stability of the method presented in this paper by performing the mentioned method for different values of , . We performed our computations in Matlab 7 software on a Pentium IV, 2800 MHz CPU machine with 1 GB of memory. Also we calculated the computational orders of the method presented in this article (denoted by C-order) with the in which and are
Concluding remarks
Some high-order methods for solving unsteady two-dimensional convection–diffusion Eq. (1.1) have been introduced in Refs. [10], [14], [15], [17], [23]. In this paper we proposed another class of high-order compact schemes for solving the mentioned equation. We combined a high-order compact finite difference scheme to approximate the spatial derivatives and the fourth-order boundary value method for the time integration of resulted linear system of ordinary differential equations. The proposed
Acknowledgment
The authors thank the anonymous referees for comments and suggestions.
References (30)
- et al.
Stability properties of some BVM methods
Appl. Numer. Math.
(1993) - et al.
Boundary value methods: the third way between linear multistep and Runge–Kutta methods
Comput. Math. Appl.
(1998) - et al.
Unsymmetric and symmetric meshless schemes for the unsteady convection–diffusion equation
Comput. Methods Appl. Mech. Eng.
(2006) Weighted finite difference techniques for the one-dimensional advection-diffusion equation
Appl. Math. Comput.
(2004)Numerical solution of the three-dimensional advection-diffusion equation
Appl. Math. Comput.
(2004)Finite difference procedures for solving a problem arising in modeling and design of certain optoelectronic devices
Math. Comput. Simulat.
(2006)- et al.
Stability of some boundary value methods for IVPs
Appl. Numer. Math.
(1995) - et al.
High order ADI method for solving unsteady convection–diffusion problems
J. Comput. Phys.
(2004) High-order compact difference schemes for time-dependent Maxwell equations
J. Comput. Phys.
(1999)- et al.
A fourth-order compact ADI method for solving two-dimensional unsteady convection–diffusion problems
J. Comput. Appl. Math.
(2007)
Multigrid method and fourth order compact difference scheme for 2D Poisson equation with unequal meshsize discretization
J. Comput. Phys.
Determination of a control parameter in the two-dimensional diffusion equation
Appl. Numer. Math.
Fourth-order techniques for identifying a control parameter in the parabolic equations
Inter. J. Eng. Sci.
The dual reciprocity boundary element method (DRBEM) for two-dimensional sine-Gordon equation
Comput. Methods Appl. Mech. Eng.
Stability of some boundary value methods for the solution of initial value problems
BIT
Cited by (63)
Study and application of pollutant diffusion based on SPH method
2023, Engineering Analysis with Boundary ElementsA new simultaneously compact finite difference scheme for high-dimensional time-dependent PDEs
2023, Mathematics and Computers in SimulationDahlquist's barriers and much beyond
2023, Journal of Computational PhysicsA new conservative fourth-order accurate difference scheme for the nonlinear Schrödinger equation with wave operator
2022, Applied Numerical MathematicsEfficient and accurate numerical methods for the multidimensional convection–diffusion equations
2019, Mathematics and Computers in SimulationCitation Excerpt :Therefore, it is more efficient than non-compact scheme. There are growing interest in the numerical solution for CD equations in recent years [3,7–9,18,20,21,24,25,27], ranging from time-splitting method, ADI method, high order compact (HOC) method [3,7,20] to high order compact exponential finite difference methods [25]. Sanyasiraju et al. considered exponential HOC method for one-dimensional CD equation [24,25].
An accurate meshless formulation for the simulation of linear and fully nonlinear advection diffusion reaction problems
2018, Advances in Engineering SoftwareCitation Excerpt :This technique was based on the use of the modified equivalent partial differential equation. In [12] new classes of high-order accurate methods were proposed for solving two-dimensional unsteady convection diffusion equations. These techniques were based on the method of lines approach.