Abstract
We study the numerical anisotropy existent in compact difference schemes as applied to hyperbolic partial differential equations, and propose an approach to reduce this error and to improve the stability restrictions based on a previous analysis applied to explicit schemes. A prefactorization of compact schemes is applied to avoid the inversion of a large matrix when calculating the derivatives at the next time level, and a predictor–corrector time marching scheme is used to update the solution in time. A reduction of the isotropy error is attained for large wave numbers and, most notably, the stability restrictions associated with MacCormack time marching schemes are considerably improved. Compared to conventional compact schemes of similar order of accuracy, the multidimensional schemes employ larger stencils which would presumably demand more processing time, but we show that the new stability restrictions render the multidimensional schemes to be in fact more efficient, while maintaining the same dispersion and dissipation characteristics of the one dimensional schemes
Similar content being viewed by others
References
Vichnevetsky, R., Bowles, J.B.: Fourier Analysis of Numerical Approximations of Hyperbolic Equations. SIAM Studies in Applied Mathematics, Philadelphia (1982)
Trefethen, L.N.: Group velocity in finite difference schemes. SIAM Rev. 24, 113 (1982)
Zingg, D.W., Lomax, H.: Finite difference schemes on regular triangular grids. J. Comput. Phys. 108, 306–313 (1993)
Tam, C.K.W., Webb, J.C.: Radiation boundary condition and anisotropy correction for finite difference solutions of the Helmholtz equation. J. Comput. Phys. 113, 122–133 (1994)
Lin, R.K., Sheu, T.W.H.: Application of dispersion-relation-preserving theory to develop a two-dimensional convection-diffusion scheme. J. Comput. Phys. 208, 493–526 (2005)
Tam, C.K.W., Webb, J.C.: Dispersion-relation-preserving finite difference schemes for computational aeroacoustics. J. Comput. Phys. 107, 262–281 (1993)
Sun, G., Trueman, C.W.: Suppression of numerical anisotropy and dispersion with optimized finite-difference time-domain methods. IEEE Trans. Antennas Propag. 53, 4121–4128 (2005)
Shen, G., Cangellaris, A.C.: A new FDTD stencil for reduced numerical anisotropy in the computer modeling of wave phenomena. Int. J. RF Microw. Comput.-Aided Eng. 17, 447–454 (2007)
Sescu, A., Afjeh, A.A., Hixon, R.: Optimized difference schemes for multidimensional hyperbolic PDEs. Electron. J. Differ. Equ. Conf. 17, 213–225 (2009)
Sescu, A., Hixon, R., Afjeh, A.A.: Multidimensional optimization of finite difference schemes for computational aeroacoustics. J. Comput. Phys. 227, 4563–4588 (2008)
Sescu, A., Afjeh, A.A., Hixon, R., Sescu, C.: Conditionally stable multidimensional schemes for advective equations. J. Sci. Comput. 42, 96–117 (2009)
Sescu, A., Hixon, R., Sescu, C., Abdollah, A.A.: Stability Investigation of Multidimensional Optimized Spatial Stencils, AIAA-Paper 2009-0005 (2009)
Hixon, R.: Prefactored small-stencil compact schemes. J. Comput. Phys. 165, 522–541 (2000)
MacCormack, R.W.: The Effect of Viscosity in Hypervelocity Impact Cratering. AIAA Paper AIAA-69-354 (1969)
Hirsch, C.: Numerical Computation of Internal and External Flows. Volume 1: Fundamentals of Numerical Discretization. Wiley, ISBN-10: 0471923850 (2001)
Gottlieb, D., Turkel, E.: Dissipative two-four methods for time-dependent problems. Math. Comput. 30, 703–723 (1976)
Bayliss, A., Parikh, P., Maestrello, L., Turkel, E.: A Fourth-Order Scheme for the Unsteady Compressible Navier–Stokes Equations, NASA CR-177994, ICASE Rep., No. 85-44 (1984)
Hixon, R., Turkel, E.: Compact implicit MacCormack-type schemes with high accuracy. J. Comput. Phys. 158, 51–70 (2000)
Kennedy, C.A., Carpenter, M.H.: Comparison of Several Numerical Methods for Simulation of Compressible Shear Layers, Technical Report NASA-97-TP3484 (1997)
Wendroff, B.: The stability of MacCormack’s method for the scalar advection equation. Appl. Math. Lett. 4, 89–91 (1991)
Hong, H.: The exact region of stability for MacCormack scheme. Computing 56, 371–383 (1996)
Acknowledgments
The authors would like to express many thanks to the reviewers for their very constructive comments. AS thanks to Dr. A. Afjeh, from the University of Toledo, for many fruitful discussions and for the substantial support.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Sescu, A., Hixon, R. Numerical Anisotropy Study of a Class of Compact Schemes. J Sci Comput 61, 327–342 (2014). https://doi.org/10.1007/s10915-014-9826-0
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10915-014-9826-0