Development of a sixth-order two-dimensional convection-diffusion scheme via Cole-Hopf transformation
Journal
Computer Methods in Applied Mechanics and Engineering
Journal Volume
191
Journal Issue
27-28
Pages
2979-2995
Date Issued
2002
Author(s)
Abstract
In this paper, we develop a two-dimensional finite-difference scheme for solving the time-dependent convection- diffusion equation. The numerical method exploits Cole-Hopf equation to transform the nonlinear scalar transport equation into the linear heat conduction equation. Within the semi-discretization context, the time derivative term in the transformed parabolic equation is approximated by a second-order accurate time-stepping scheme, resulting in an inhomogeneous Helmholtz equation. We apply the alternating direction implicit scheme of Polezhaev to solve the Helmholtz equation. As the key to success in the present simulation, we develop a Helmholtz scheme with sixth-order spatial accuracy. As is standard practice, we validated the code against test problems which were amenable to exact solutions. Results show excellent agreement for the one-dimensional test problems and good agreement with the analytical solution for the two-dimensional problem. © 2002 Elsevier Science B.V. All rights reserved.
Type
journal article