- Journal Home
- Volume 36 - 2024
- Volume 35 - 2024
- Volume 34 - 2023
- Volume 33 - 2023
- Volume 32 - 2022
- Volume 31 - 2022
- Volume 30 - 2021
- Volume 29 - 2021
- Volume 28 - 2020
- Volume 27 - 2020
- Volume 26 - 2019
- Volume 25 - 2019
- Volume 24 - 2018
- Volume 23 - 2018
- Volume 22 - 2017
- Volume 21 - 2017
- Volume 20 - 2016
- Volume 19 - 2016
- Volume 18 - 2015
- Volume 17 - 2015
- Volume 16 - 2014
- Volume 15 - 2014
- Volume 14 - 2013
- Volume 13 - 2013
- Volume 12 - 2012
- Volume 11 - 2012
- Volume 10 - 2011
- Volume 9 - 2011
- Volume 8 - 2010
- Volume 7 - 2010
- Volume 6 - 2009
- Volume 5 - 2009
- Volume 4 - 2008
- Volume 3 - 2008
- Volume 2 - 2007
- Volume 1 - 2006
Commun. Comput. Phys., 20 (2016), pp. 486-511.
Published online: 2018-04
Cited by
- BibTex
- RIS
- TXT
Implicit time integration schemes are popular because their relaxed stability constraints can result in better computational efficiency. For time-accurate unsteady simulations, it has been well recognized that the inherent dispersion and dissipation errors of implicit Runge-Kutta schemes will reduce the computational accuracy for large time steps. Yet for steady state simulations using the time-dependent governing equations, these errors are often overlooked because the intermediate solutions are of less interest. Based on the model equation dy/dt = (µ+iλ)y of scalar convection diffusion systems, this study examines the stability limits, dispersion and dissipation errors of four diagonally implicit Runge-Kutta-type schemes on the complex (µ+iλ)∆t plane. Through numerical experiments, it is shown that, as the time steps increase, the A-stable implicit schemes may not always have reduced CPU time and the computations may not always remain stable, due to the inherent dispersion and dissipation errors of the implicit Runge-Kutta schemes. The dissipation errors may decelerate the convergence rate, and the dispersion errors may cause large oscillations of the numerical solutions. These errors, especially those of high wavenumber components, grow at large time steps. They lead to difficulty in the convergence of the numerical computations, and result in increasing CPU time or even unstable computations as the time step increases. It is concluded that an optimal implicit time integration scheme for steady state simulations should have high dissipation and low dispersion.
}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.OA-2016-0032}, url = {http://global-sci.org/intro/article_detail/cicp/11161.html} }Implicit time integration schemes are popular because their relaxed stability constraints can result in better computational efficiency. For time-accurate unsteady simulations, it has been well recognized that the inherent dispersion and dissipation errors of implicit Runge-Kutta schemes will reduce the computational accuracy for large time steps. Yet for steady state simulations using the time-dependent governing equations, these errors are often overlooked because the intermediate solutions are of less interest. Based on the model equation dy/dt = (µ+iλ)y of scalar convection diffusion systems, this study examines the stability limits, dispersion and dissipation errors of four diagonally implicit Runge-Kutta-type schemes on the complex (µ+iλ)∆t plane. Through numerical experiments, it is shown that, as the time steps increase, the A-stable implicit schemes may not always have reduced CPU time and the computations may not always remain stable, due to the inherent dispersion and dissipation errors of the implicit Runge-Kutta schemes. The dissipation errors may decelerate the convergence rate, and the dispersion errors may cause large oscillations of the numerical solutions. These errors, especially those of high wavenumber components, grow at large time steps. They lead to difficulty in the convergence of the numerical computations, and result in increasing CPU time or even unstable computations as the time step increases. It is concluded that an optimal implicit time integration scheme for steady state simulations should have high dissipation and low dispersion.