A Stable and High Accurate Numerical Simulator for Convection-Diffusion Equation

Year 2022, Volume: 8 Issue: 1, 45 - 53, 10.03.2022

Osman Ünal

https://doi.org/10.28979/jarnas.985946

Abstract

This study presents a simulator to obtain numerical solution of convection-diffusion equation. It includes explicit, fully implicit and semi-implicit time discretization techniques. Although the explicit method is efficient for some simple conduction problems, it needs stability criteria. Otherwise, it is unstable especially for high Courant number and large time step size. When using the explicit method, the time step size is chosen with care to get well-posed numerical solution. This requirement sets a serious constraint for the explicit method because choosing small time step size causes the simulation time to become quite long. Even though the implicit method is stable for large time step size and high Courant number, it leads to numerical dispersion like the explicit method. On the other hand, second-order accurate semi-implicit method significantly reduces numerical dispersion. In addition to time discretization techniques, this simulator contains several space discretization methods such as first-order upstream and UMIST (University of Manchester Institute of Science and Technology) techniques. The proposed numerical simulator is suitable for easily using the different combinations of time and space discretization methods. Secondly, this study is to propose the use of semi-implicit time discretization technique with UMIST space discretization method to minimize numerical dispersion and suppress unphysical oscillation. In spite of the fact that the UMIST method suppresses unphysical oscillation, it causes a small and undesired oscillation at flood front for very large Courant number. Third objective of this study is to propose minor modification on the UMIST method to eliminate this unphysical oscillation.

Keywords

Finite difference method, high-resolution scheme, explicit method, fully implicit method, semi-implicit method

References

• Crank, J., & Nicolson, P. (1947, January). A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type. In Mathematical Proceedings of the Cambridge Philosophical Society (Vol. 43, No. 1, pp. 50-67). Cambridge University Press.
• Ertekin, T., Abou-Kassem, J. H., & King, G. R. (2001). Basic applied reservoir simulation (Vol. 7). Richardson, TX: Society of Petroleum Engineers.
• Harten, A. (1984). On a class of high resolution total-variation-stable finite-difference schemes. SIAM Journal on Numerical Analysis, 21(1), 1-23.
• Kamalyar, K., Kharrat, R., & Nikbakht, M. (2014). Numerical Aspects of the Convection–Dispersion Equation. Petroleum science and technology, 32(14), 1729-1762.
• Kurganov, A., & Tadmor, E. (2000). New high-resolution central schemes for nonlinear conservation laws and convection–diffusion equations. Journal of Computational Physics, 160(1), 241-282.
• Leonard, B. P. (1979). A stable and accurate convective modelling procedure based on quadratic upstream interpolation. Computer methods in applied mechanics and engineering, 19(1), 59-98.
• Morton, K. W. (2019). Numerical solution of convection-diffusion problems. CRC Press.
• Peaceman, D. W. (2000). Fundamentals of numerical reservoir simulation. Elsevier.
• Peng, Y., Liu, C., & Shi, L. (2013, August). Soution of Convection-Diffusion Equations. In International Conference on Information Computing and Applications (pp. 546-555). Springer, Berlin, Heidelberg.
• Sarra, S. A. (2003). The method of characteristics with applications to conservation laws. Journal of Online mathematics and its Applications, 3, 1-16.
• Sweby, P. K. (1984). High resolution schemes using flux limiters for hyperbolic conservation laws. SIAM journal on numerical analysis, 21(5), 995-1011.
• Versteeg, H. K., & Malalasekera, W. (2007). An introduction to computational fluid dynamics: the finite volume method. Pearson education.
• Wolcott, D. S., Kazemi, H., & Dean, R. H. (1996, October). A practical method for minimizing the grid orientation effect in reservoir simulation. In SPE annual technical conference and exhibition. OnePetro.