of the discretization is then minimized. This idea is called global discretisation where attention is no longer paid to The three-dimensional, time-dependent convection-diffusion equation (CDE) is considered. An exponential transformation is used the individual spatial and temporal derivatives but to the to collectively transform the CDE. The idea of global discretization whole equation. With this idea, four finite difference is used, where attention is paid to the whole transformed CDE, but schemes are established for both the CDE and the transnot to the individual spatial and temporal derivatives in the equation. formed CDE. The modified PDEs, following Warming and Four finite difference schemes for both CDE and transformed CDE Hyett [10], are obtained, which indicate that the schemes are established. The modified partial differential equations of these schemes are obtained, which indicate that the trunction errors of are of either second or fourth order if the time step is the schemes can be of second and fourth order, depending on the properly prescribed. A series of analytical solutions to both prescription of the time step length. Some characteristic physical the Navier-Stokes and Burgers equations are chosen as parameters, i.e., local Reynolds number, local Strouhal number, benchmark cases against which the new method is assessed. and viscous diffusive length, are introduced into the schemes and
The paper is organized as follows. Sections 2 and 3 demonthe viscous diffusive length is found to be a significant parameter in relating temporal discretisation with spatial discretisation. A se-strate the major steps in establishing the finite difference ries of benchmark analytical solutions of Navier-Stokes and Burgers schemes. Section 4 provides the formulation of the schemes equations, as well as the numerical solutions using the well-known for the CDE and the schemes for the transformed CDE discretisation schemes, are used to investigate the properties of the are given in the Appendix. Section 5 presents the analytical derived schemes. The high-order schemes achieve higher resolusolutions and the numerical algorithms chosen to benchtions over the conventional schemes without decreasing much the sparsity of the matrix structures. Grid refinement studies reveal mark the schemes. The numerical experiments are dethat the inverse exponential transformation of the finite difference scribed in Section 6. The results are discussed in Section schemes tends to destroy some resolution of the schemes, espe-7 and conclusions drawn from the work close out the cially for large local Reynolds number.