Hostname: page-component-586b7cd67f-g8jcs Total loading time: 0 Render date: 2024-11-22T23:30:47.915Z Has data issue: false hasContentIssue false
  • English
  • Français

On the Stability and CPU Time of the Implicit Runge-Kutta Schemes for Steady State Simulations

Part of: Basic methods in fluid mechanics Partial differential equations, initial value and time-dependent initial-boundary value problems

Published online by Cambridge University Press:  21 July 2016

Yongle Du  and
John A. Ekaterinaris
Yongle Du*
Affiliation:
Embry-Riddle Aeronautical University, Daytona Beach, FL 32114, USA
John A. Ekaterinaris*
Affiliation:
Embry-Riddle Aeronautical University, Daytona Beach, FL 32114, USA
*
*Corresponding author. Email addresses:[email protected] (Y. Du), [email protected] (J. A. Ekaterinaris)
*Corresponding author. Email addresses:[email protected] (Y. Du), [email protected] (J. A. Ekaterinaris)
Get access

Abstract

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.

Keywords

Runge-Kuttaimplicit time integrationdispersiondissipationnumerical stabilitysteady state simulation

MSC classification

Secondary: 65M22: Solution of discretized equations 65M12: Stability and convergence of numerical methods 76M20: Finite difference methods 76M25: Other numerical methods
Type
Research Article
Information
Communications in Computational Physics , Volume 20 , Issue 2 , August 2016 , pp. 486 - 511
Copyright
Copyright © Global-Science Press 2016 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1] Cranka, J., Nicolson, P., A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type, Mathematical Proceedings of the Cambridge Philosophical Society 43 (1) (1947) 50–67.Google Scholar
[2] Lax, P., Wendroff, B., Systems of conservation laws, Communications on Pure and Applied Mathematics 13 (2) (1960) 217–237.CrossRefGoogle Scholar
[3] Donea, J., A Taylor-Galerkin method for convective transport problems, International Journal for Numerical Methods in Engineering 20 (1) (1984) 101–119.Google Scholar
[4] Butcher, J., Numerical methods for ordinary differential equations, 2nd Edition, John Wiley & Sons Ltd., New York, 2008.Google Scholar
[5] Butcher, J., A history of Runge-Kutta methods, Applied Numerical Mathematics 20 (1996) 247–260.Google Scholar
[6] Alexander, R., Diagonally implicit Runge-Kutta methods for stiff O.D.E.’S, Journal of Computational Physics 14 (1977) 1006–1021.Google Scholar
[7] Hu, F. Q., Hussaini, M. Y., Manthey, J. L., Low-dissipation and low-dispersion Runge-Kutta schemes for computational acoustics, Journal of Computational Physics 124 (1996) 177–191.Google Scholar
[8] Bogey, C., Bailly, C., A family of low dispersive and low dissipative explicit schemes for flow and noise computations, Journal of Computational Physics 194 (2004) 194–214.Google Scholar
[9] Berland, J., Bogey, C., Bailly, C., Low-dissipation and low-dispersion fourth-order Runge-Kutta algorithm, Computers & Fluids 35 (2006) 1459–1463.CrossRefGoogle Scholar
[10] Najafi-Yazdi, A., Mongeau, L., A low-dispersion and low-dissipation implicit Runge-Kutta scheme, Journal of Computational Physics 233 (2013) 315–323.Google Scholar
[11] Nazari, F., Mohammadian, A., Charron, M., Zadra, A., Optimal high-order diagonally-implicit Runge-Kutta schemes for nonlinear diffusive systems on atmospheric boundary layer, Journal of Computational Physics 271 (2014) 118–130.Google Scholar
[12] Nazari, F., Mohammadian, A., Charron, M., High-order low-dissipation low-dispersion diagonally implicit Runge-Kutta schemes, Journal of Computational Physics 286 (2015) 38–48.Google Scholar
[13] Du, Y., Ekaterinaris, J. A., Optimization of the dinagonally implicit Runge-Kutta methods for convection diffusion equations, under review, AIAA Journal (2016).Google Scholar
[14] Jameson, A., Baker, T. J., Solution of the Euler equations for complex configurations, in: AIAA paper 83-1929-CP, 1983.Google Scholar
[15] Swanson, R. C., Turkel, E., Artificial dissipation and central difference schemes for the Euler and Navier-Stokes equations, in: AIAA paper 87-1107, 1987.Google Scholar
[16] Hairer, E., Wanner, G., Solving ordinary differential equations II: stiff and differential-algebraic problems, 2nd Edition, Springer, Sydney, Australia, 2004.Google Scholar
[17] Prothero, A., Robinson, A., On the stability and accuracy of one-step methods for solving stiff systems of ordinary differential equations, Mathematics of Computation 28 (1974) 145–162.Google Scholar
[18] Gottlieb, S., Shu, C., Total variation diminishing Runge-Kutta schemes, Mathematics of Computation 67 (221) (1998) 73–85.Google Scholar
[19] Gottlieb, S., Shu, C., Tadmor, E., Strong stability-preserving high-order time discretization methods, SIAM Review 43 (1) (2001) 89–112.Google Scholar
[20] Ruuth, S. J., Spiteri, R. J., High-order strong stability-preserving Runge-Kutta methods with downwind-biased spatial discretiztions, SIAM Journal on Numerical Analysis 42 (3) (2004) 974996.Google Scholar
[21] Gottlieb, S., On high order strong stability-preserving Runge-Kutta and multistep time discretizations, Journal of Scientific Computing 25 (112) (2005) 89–112.Google Scholar
[22] Ketcheson, D. I., Gottlieb, S., MacDonald, C. B., High-order strong stability-preserving Runge-Kutta methods with downwind-biased spatial discretiztions, SIAM Journal on Numerical Analysis 49 (6) (2011) 26182639.Google Scholar
[23] Tam, C. K. W., Webb, J. C., Dispersion-relation-preserving finite difference schemes for computational aeroacoustics, Journal of Computational Physics 107 (1) (1993) 262–281.Google Scholar
[24] Shu, C., Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws, Tech. Rep. NASA/CR-97-206253, NASA (2000).Google Scholar
[25] Henrick, A. K., Aslam, T. D., Powers, J. M., Mapped weighted essentially non-oscillatory schemes: achieving optimal order near critical points, Journal of Computational Physics 207 (2005) 542–567.Google Scholar
[26] Burgers, J. M., A mathematical model illustrating the theory of turbulence, Advances in Applied Mechanics 1 (1948) 171–199.Google Scholar
[27] Burns, J., Balogh, A., Gilliam, D. S., Shubov, V. I., Numerical stationary solutions for a viscous Burgers's equation, Journal of Mathematical Systems, Estimation, and Control 8 (2) (1998) 1–16.Google Scholar
[28] Yoon, S., Jameson, A., Lower-upper implicit schemes with multiple grids for the Euler equations, AIAA Journal 25 (7) (1987) 929–935.Google Scholar
[29] Yoon, S., Jameson, A., A lower-upper symmetric-Gaussian Seidel method for the Euler and Navier-Stokes equations, AIAA Journal 26 (9) (1988) 1025–1026.CrossRefGoogle Scholar
[30] White, F. M., Viscous Fluid Flow, 3rd Edition, McGraw Hill Inc., New York, 2005.Google Scholar