Optimised prefactored compact schemes for linear wave propagation phenomena

详细信息    查看全文

• 作者：A. Rona^a ; ^{ar45@leicester.ac.uk" class="auth_mail" title="E-mail the corresponding author} ; I. Spisso^b ; ^{i.spisso@cineca.it" class="auth_mail" title="E-mail the corresponding author} ; E. Hall^a ; ^{eh171@leicester.ac.uk" class="auth_mail" title="E-mail the corresponding author} ; M. Bernardini^c ; ^{matteo.bernardini@uniroma1.it" class="auth_mail" title="E-mail the corresponding author} ; S. Pirozzoli^c ; ^{sergio.pirozzoli@uniroma1.it" class="auth_mail" title="E-mail the corresponding author}
• 关键词：Computational aeroacoustics ; Prefactored compact schemes ; Optimised algorithms ; Numerical dispersion ; Numerical dissipation ; Reduced numerical error
• 刊名：Journal of Computational Physics
• 出版年：2017
• 出版时间：1 January 2017
• 年：2017
• 卷：328
• 期：Complete
• 页码：66-85
• 全文大小：1787 K

文摘

A family of space- and time-optimised prefactored compact schemes are developed that minimise the computational cost for given levels of numerical error in wave propagation phenomena, with special reference to aerodynamic sound. This work extends the approach of Pirozzoli [1] to the MacCormack type prefactored compact high-order schemes developed by Hixon [2], in which their shorter Padé stencil from the prefactorisation leads to a simpler enforcement of numerical boundary conditions. An explicit low-storage multi-step Runge–Kutta integration advances the states in time. Theoretical predictions for spatial and temporal error bounds are derived for the cost-optimised schemes and compared against benchmark schemes of current use in computational aeroacoustic applications in terms of computational cost for a given relative numerical error value. One- and two-dimensional test cases are presented to examine the effectiveness of the cost-optimised schemes for practical flow computations. An effectiveness up to about 50% higher than the standard schemes is verified for the linear one-dimensional advection solver, which is a popular baseline solver kernel for computational physics problems. A substantial error reduction for a given cost is also obtained in the more complex case of a two-dimensional acoustic pulse propagation, provided the optimised schemes are made to operate close to their nominal design points.