×

Low-storage, explicit Runge-Kutta schemes for the compressible Navier-Stokes equations. (English) Zbl 0986.76060

From the summary: We derive low-storage, explicit Runge-Kutta (ERK) schemes in the context of integrating compressible Navier-Stokes equations via direct numerical simulation. Optimization of ERK methods is done across the broad range of properties, such as linear and nonlinear stability, error control reliability, step change stability, and dissipation/dispersion accuracy, subject to varying degrees of memory economization. Sixteen ERK pairs are presented using from two to five registers of memory per equation and per grid point, and having accuracies from third- to fifth-order. Methods have been tested with the one-dimensional wave equation.

MSC:

76M20 Finite difference methods applied to problems in fluid mechanics
76M25 Other numerical methods (fluid mechanics) (MSC2010)
76N15 Gas dynamics (general theory)
65M20 Method of lines for initial value and initial-boundary value problems involving PDEs
Full Text: DOI

References:

[1] Abarbanel, S.; Gottlieb, D.; Carpenter, M. H., On the removal of boundary errors caused by Runge-Kutta integration of nonlinear partial differential equations, SIAM J. Sci. Comput., 17, 3, 777-782 (1996) · Zbl 0858.65096
[2] Alonso-Mallo, I., Explicit single-step methods with optimal order of convergence for partial differential equations, Appl. Numer. Math., 31, 2, 117-131 (1999) · Zbl 0941.65095
[3] Aves, M. A.; Griffiths, D. F.; Higham, D. J., Does error control suppress spuriosity?, SIAM J. Numer. Anal., 34, 2, 756-778 (1997) · Zbl 0873.65083
[4] Baker, T. S.; Dormand, J. R.; Gilmore, J. P.; Prince, P. J., Continuous approximation with embedded Runge-Kutta methods, Appl. Numer. Math., 22, 1-3, 51-62 (1996) · Zbl 0871.65077
[5] Blaisdell, G. A.; Mansour, N. N.; Reynolds, W. C., Compressibility effects on the growth and structure of homogeneous turbulent shear flow, J. Fluid Mech., 256, 443-485 (1993) · Zbl 0800.76186
[6] Bogacki, P.; Shampine, L. F., A (3,2) pair of Runge-Kutta formulas, Appl. Math. Lett., 2, 4, 321-325 (1989) · Zbl 0705.65055
[7] Bogacki, P.; Shampine, L. F., An efficient Runge-Kutta (4,5) pair, Comput. Math. Appl., 32, 6, 15-28 (1996) · Zbl 0857.65077
[8] Butcher, J. C., Coefficients for the study of Runge-Kutta integration processes, J. Austral. Math. Soc. Ser. B, 3, 2, 185-201 (1964) · Zbl 0223.65031
[9] Butcher, J. C., On Runge-Kutta processes of high order, J. Austral. Math. Soc. Ser. B, 4, 2, 179-194 (1964) · Zbl 0244.65046
[10] Butcher, J. C., The Numerical Analysis of Ordinary Differential Equations: Runge-Kutta and General Linear Methods (1987), Wiley: Wiley Chichester · Zbl 0616.65072
[11] Calvo, M.; Franco, J. M.; Montijano, J. I.; Rández, L., Explicit Runge-Kutta methods for initial problems with oscillating solutions, J. Comput. Appl. Math., 76, 1-2, 195-212 (1996) · Zbl 0871.65073
[12] Carpenter, M. H.; Kennedy, C. A., Fourth-order 2N-storage Runge-Kutta schemes, NASA TM-109112 (1994), National Aeronautics and Space Administration, Langley Research Center: National Aeronautics and Space Administration, Langley Research Center Hampton, VA
[13] Carpenter, M. H.; Kennedy, C. A., Third-order 2N-storage Runge-Kutta schemes with error control, NASA TM-109111 (1994), National Aeronautics and Space Administration, Langley Research Center: National Aeronautics and Space Administration, Langley Research Center Hampton, VA
[14] Carpenter, M. H.; Gottlieb, D.; Abarbanel, S.; Don, W.-S., The theoretical accuracy of Runge-Kutta time discretizations for the initial boundary value problem: a study of the boundary error, SIAM J. Sci. Comput., 16, 6, 1241-1252 (1995) · Zbl 0839.65098
[15] Conte, S. D.; Reeves, R. F., A Kutta third-order procedure for solving differential equations requiring minimum storage, SIAM J. Numer. Anal., 12, 4, 22-25 (1956)
[16] Cooper, G. J., A generalization of algebraic stability for Runge-Kutta methods, IMA J. Numer. Anal., 4, 4, 427-440 (1984) · Zbl 0557.65045
[17] Dahlquist, G. G.; Jeltsch, R., Generalized disks of contractivity for explicit and implicit Runge-Kutta methods, Report TRITA-NA-7906 (1979), Department of Numerical Analysis and Computer Science, Royal Inst. Tech.: Department of Numerical Analysis and Computer Science, Royal Inst. Tech. Stockholm
[18] Dekker, K.; Verwer, J. G., Stability of Runge-Kutta Methods for Stiff Nonlinear Differential Equations (1984), North-Holland: North-Holland Amsterdam · Zbl 0571.65057
[19] Dennis, J. E.; Schnabel, R. B., Numerical Methods for Unconstrained Optimization and Nonlinear Equations (1983), Prentice-Hall: Prentice-Hall Englewood Cliffs, NJ · Zbl 0579.65058
[20] Dormand, J. R.; Prince, P. J., A family of embedded Runge-Kutta formulae, J. Comput. Appl. Math., 6, 1, 19-26 (1980) · Zbl 0448.65045
[21] Dormand, J. R.; Prince, P. J., A reconsideration of some embedded Runge-Kutta formulae, J. Comput. Appl. Math., 15, 2, 203-211 (1986) · Zbl 0602.65046
[22] Dormand, J. R., Numerical Methods for Differential Equations (1996), CRC Press: CRC Press Boca Raton, FL · Zbl 0847.65046
[23] Enright, W. H.; Jackson, K. R.; Nørsett, S. P.; Thomsen, P. G., Interpolants for Runge-Kutta formulas, ACM Trans. Math. Software, 12, 3, 193-218 (1986) · Zbl 0617.65068
[24] Enright, W. H.; Pryce, J. D., Two FORTRAN packages for assessing initial value methods, ACM Trans. Math. Software, 13, 1, 1-27 (1987) · Zbl 0617.65069
[25] Enright, W. H.; Higham, D. J.; Owren, B.; Sharp, P. W., A survey of the explicit Runge-Kutta method, Technical Report, 291/94 (1994), Department of Computer Science, University of Toronto: Department of Computer Science, University of Toronto Toronto
[26] Erlebacher, G.; Hussaini, M. Y.; Shu, C.-W., Interaction of a shock with a longitudinal vortex, J. Fluid Mech., 337, 129-153 (1997) · Zbl 0889.76033
[27] Fehlberg, E., Low-order classical Runge-Kutta formulas with stepsize control and their application to some heat transfer problems, NASA Technical Report TR R-315 (1969), National Aeronautics and Space Administration, Marshall Space Flight Center: National Aeronautics and Space Administration, Marshall Space Flight Center Marshall, AL
[28] Fyfe, D. J., Economical evaluation of Runge-Kutta formulae, Math. Comp., 20, 392-398 (1966) · Zbl 0143.17207
[29] Gay, D. M., Computing optimal locally constrained steps, SIAM J. Sci. Statist. Comput., 2, 2, 186-197 (1981) · Zbl 0467.65027
[30] Gay, D. M., ALGORITHM 611. Subroutines for unconstrained minimization using a model/trust region approach, ACM Trans. Math. Software, 9, 4, 503-524 (1983) · Zbl 0519.65039
[31] Gill, S., A process for the step-by-step integration of differential equations in an automatic digital computing machine, Proc. Cambridge Phil. Soc., 47, 96-108 (1951) · Zbl 0042.13202
[32] Gottlieb, S.; Shu, C.-W., Total variation diminishing Runge-Kutta schemes, Math. Comp., 67, 221, 73-85 (1998) · Zbl 0897.65058
[33] Gran, I. R., Negative flame speed in an unsteady 2-D premixed flame: A computational study, (26th Symposium (International) on Combustion (1996), The Combustion Institute: The Combustion Institute Pittsburgh, PA), 323-329
[34] Griffiths, D. F.; Christie, I.; Mitchell, A. R., Analysis of error growth for explicit difference schemes in conduction-convection problems, Internat. J. Numer. Methods Engrg., 15, 7, 1075-1081 (1980) · Zbl 0432.76077
[35] Gustafsson, K., Control theoretic techniques for stepsize selection in Runge-Kutta methods, ACM Trans. Math. Software, 17, 4, 533-554 (1991) · Zbl 0900.65256
[36] Hairer, E.; Nørsett, S. P.; Wanner, G., Solving Ordinary Differential Equations I, Nonstiff Problems (1993), Springer-Verlag: Springer-Verlag Berlin · Zbl 0789.65048
[37] Hairer, E.; Wanner, G., Solving Ordinary Differential Equations II, Stiff and Differential-Algebraic Problems (1996), Springer-Verlag: Springer-Verlag Berlin · Zbl 0859.65067
[38] Hall, G., A new stepsize strategy for explicit Runge-Kutta codes, Adv. Comput. Math., 3, 4, 343-352 (1995) · Zbl 0833.65086
[39] Higham, D. J., Regular Runge-Kutta pairs, Appl. Numer. Math., 25, 2/3, 229-241 (1997) · Zbl 0890.65076
[40] Horváth, Z., Positivity of Runge-Kutta and diagonally split Runge-Kutta methods, Appl. Numer. Math., 28, 2-4, 309-326 (1998) · Zbl 0926.65073
[41] Hosea, M. E.; Shampine, L. F., Efficiency comparisons of methods for integrating ODEs, Comput. Math. Appl., 28, 6, 45-55 (1994) · Zbl 0807.65083
[42] van der Houwen, P. J., Explicit Runge-Kutta formulas with increased stability boundaries, Numer. Math., 20, 2, 149-164 (1972) · Zbl 0233.65039
[43] van der Houwen, P. J., Construction of Integration Formulas for Initial Value Problems (1977), North-Holland: North-Holland Amsterdam · Zbl 0359.65057
[44] van der Houwen, P. J.; Sommeijer, B. P., On the internal stability of explicit, \(m\)-stage Runge-Kutta methods for large \(m\)-values, Z. Angew. Math. Mech., 60, 10, 479-485 (1980) · Zbl 0455.65052
[45] van der Houwen, P. J.; Sommeijer, B. P., Explicit Runge-Kutta (-Nyström) methods with reduced phase errors for computing oscillating solutions, SIAM J. Numer. Anal., 24, 3, 595-616 (1987) · Zbl 0624.65058
[46] van der Houwen, P. J.; Sommeijer, B. P., Phase-lag analysis of implicit Runge-Kutta methods, SIAM J. Numer. Anal., 26, 1, 214-229 (1989) · Zbl 0669.65055
[47] Hu, F. Q.; Hussaini, M. Y.; Manthey, J. L., Low-dissipation and low-dispersion Runge-Kutta schemes for computational acoustics, J. Comput. Phys., 124, 1, 177-191 (1996) · Zbl 0849.76046
[48] Hundsdorfer, W. H., Numerical solution of advection-diffusion-reaction equations, CWI Report NM-N9603 (1996), Centrum voor Wiskunde en Informatica: Centrum voor Wiskunde en Informatica Amsterdam
[49] Kennedy, C. A.; Carpenter, M. H., A comparison of several new numerical methods for the simulation of compressible shear layers, Appl. Numer. Math., 14, 4, 397-433 (1994) · Zbl 0804.76062
[50] Kennedy, C. A.; Carpenter, M. H.; Lewis, R. M., Low-storage, explicit Runge-Kutta schemes for the compressible Navier-Stokes equations, ICASE Report No. 99-22, NASA/CR-1999-209349 (1999), Institute for Computer Applications in Science and Engineering: Institute for Computer Applications in Science and Engineering Hampton
[51] Koto, T., Explicit Runge-Kutta schemes for evolutionary problems in partial differential equations, Ann. Numer. Math., 1, 335-346 (1994) · Zbl 0836.65101
[52] Kraaijevanger, J. F.B. M.; Lenferink, H. W.J.; Spijker, M. N., Stepsize restrictions for stability in the numerical solution of ordinary and partial differential equations, J. Comput. Appl. Math., 20, 1, 67-81 (1987) · Zbl 0633.65065
[53] Kraaijevanger, J. F.B. M., Absolute monotonicity of polynomials occurring in the solution of initial value problems, Numer. Math., 48, 3, 303-322 (1986) · Zbl 0561.65055
[54] Kraaijevanger, J. F.B. M., Contractivity of Runge-Kutta methods, BIT, 31, 3, 482-528 (1991) · Zbl 0763.65059
[55] Kraaijevanger, J. F.B. M., Contractivity in the maximum norm for Runge-Kutta methods, (Cash, J. R.; Gladwell, I., Computational Ordinary Differential Equations (1992), Clarendon Press: Clarendon Press New York), 9-18 · Zbl 0774.65049
[56] Lambert, J. D., Numerical Methods for Ordinary Differential Systems. The Initial Value Problem (1991), Wiley: Wiley Chichester · Zbl 0745.65049
[57] Li, S. F., The nonlinear stability of explicit and diagonally implicit Runge-Kutta methods, Math. Numer. Sinica, 9, 4, 419-430 (1987), (in Chinese) · Zbl 0645.65037
[58] Lorenz, E. N., An N-cycle time differencing scheme for stepwise numerical integration, Monthly Weather Rev., 99, 8, 644-648 (1971)
[59] van der Marel, R. P., Stability radius of polynomials occurring in the numerical solution of initial value problems, BIT, 30, 3, 516-528 (1990) · Zbl 0711.65067
[60] Moré, J. J.; Garbow, B. S.; Hillstrom, K. E., User Guide for MINPACK-1, Argonne National Laboratories Report ANL-80-74 (1980), Argonne National Laboratories: Argonne National Laboratories Argonne, IL
[61] Müller, B., Linear stability condition for explicit Runge-Kutta methods to solve the compressible Navier-Stokes equations, Math. Methods Appl. Sci., 12, 2, 139-151 (1990) · Zbl 0687.76072
[62] Papakostas, S. N.; Papageorgiou, G., A family of fifth-order Runge-Kutta pairs, Math. Comp., 65, 215, 1165-1181 (1996) · Zbl 0938.65098
[63] Pathria, D., The correct formulation of intermediate boundary conditions for Runge-Kutta time integration of initial boundary value problems, SIAM J. Sci. Comput., 18, 5, 1255-1266 (1997) · Zbl 0897.65057
[64] Prince, P. J.; Dormand, J. R., High order embedded Runge-Kutta formulae, J. Comput. Appl. Math., 7, 1, 67-75 (1981) · Zbl 0449.65048
[65] Pruett, C. D.; Zang, T. A.; Chang, C. L.; Carpenter, M. H., Spatial direct numerical simulation of high-speed boundary-layer flows. I. Algorithmic considerations and validation, Theor. Comput. Fluid Dynamics, 7, 1, 49-76 (1995) · Zbl 0824.76055
[66] Pruett, C. D.; Zang, T. A., On simulation and analysis of instability and transition in high-speed boundary-layer flows, Comput. Systems Engrg., 6, 6, 563-575 (1995)
[67] Reddy, S. C.; Trefethen, L. N., Stability of the method of lines, Numer. Math., 62, 2, 235-267 (1992) · Zbl 0734.65077
[68] Sanz-Serna, J. M.; Verwer, J. G.; Hundsdorfer, W. H., Convergence and order reduction of Runge-Kutta schemes applied to evolutionary problems in partial differential equations, Numer. Math., 50, 4, 405-418 (1987) · Zbl 0589.65069
[69] Shampine, L. F., Storage reduction for Runge-Kutta codes, ACM Trans. Math. Software, 5, 3, 245-250 (1979) · Zbl 0413.65055
[70] Shampine, L. F., Some practical Runge-Kutta formulas, Math. Comp., 46, 173, 135-150 (1986) · Zbl 0594.65046
[71] Shampine, L. F., Numerical Solution of Ordinary Differential Equations (1994), Chapman-Hall: Chapman-Hall New York · Zbl 0826.65082
[72] Sharp, P. W., New low-order explicit Runge-Kutta pairs, Technical Report, 222/89 (1989), Department of Computer Science, University of Toronto: Department of Computer Science, University of Toronto Toronto
[73] Sharp, P. W., Numerical comparisons of some explicit Runge-Kutta pairs of orders 4-8, ACM Trans. Math. Software, 17, 3, 387-409 (1991) · Zbl 0900.65236
[74] Sharp, P. W.; Smart, E., Explicit Runge-Kutta pairs with one more derivative evaluation than the minimum, SIAM J. Sci. Comput., 14, 2, 338-348 (1993) · Zbl 0773.65051
[75] Sharp, P. W.; Verner, J. H., Completely imbedded Runge-Kutta pairs, SIAM J. Numer. Anal., 31, 4, 1169-1190 (1994) · Zbl 0809.65081
[76] Shu, C.-W.; Osher, S., Efficient implementation of essentially non-oscillatory schemes, J. Comput. Phys., 77, 2, 439-471 (1988) · Zbl 0653.65072
[77] Sowa, J., Stability of a Runge-Kutta method for the Navier-Stokes equation, BIT, 30, 3, 542-560 (1990) · Zbl 0697.76083
[78] Spijker, M. N., Contractivity of Runge-Kutta methods, (Dahlquist, G.; Jeltsch, R., Proc. of Numerical Methods for Solving Stiff Initial Value Problems (1981), Institut für Geometrie und Praktische Mathematik der RWTH Aachen, Bericht Nr. 9) · Zbl 0571.65063
[79] Stanescu, D.; Habashi, W. G., 2N-storage low dissipation and dispersion Runge-Kutta schemes for computational acoustics, J. Comput. Phys., 143, 2, 674-681 (1998) · Zbl 0952.76063
[80] Stein, O., Bifurcations of hyperbolic fixed points for explicit Runge-Kutta methods, IMA J. Numer. Anal., 17, 2, 151-175 (1997), (Errata 18 (2) (1998) 329) · Zbl 0872.65077
[81] Stuart, A. M., Nonlinear instability in dissipative finite-difference schemes, SIAM Rev., 31, 2, 191-220 (1989) · Zbl 0675.65099
[82] Stuart, A. M.; Humphries, A. R., Dynamical Systems and Numerical Analysis (1996), Cambridge University Press: Cambridge University Press Cambridge · Zbl 0869.65043
[83] Trouve, A.; Poinsot, T., The evolution equation for the flame surface density in turbulent premixed combustion, J. Fluid Mech., 278, 1-31 (1994) · Zbl 0825.76899
[84] Tsitouras, Ch.; Papkostas, S. N., Cheap error estimation for Runge-Kutta methods, SIAM J. Sci. Comput., 20, 6, 2067-2088 (1999) · Zbl 0935.65074
[85] Vadillo, F., On spurious fixed points of Runge-Kutta methods, J. Comput. Phys., 132, 1, 78-90 (1997) · Zbl 0877.65050
[86] Verner, J. H., Explicit Runge-Kutta methods with estimates of the local truncation error, SIAM J. Numer. Anal., 15, 4, 772-790 (1978) · Zbl 0403.65029
[87] Verner, J. H., Some Runge-Kutta formula pairs, SIAM J. Numer. Anal., 28, 2, 496-511 (1991) · Zbl 0726.65082
[88] Verner, J. H., High-order explicit Runge-Kutta pairs with low stage order, Appl. Numer. Math., 22, 1-3, 345-357 (1996) · Zbl 0868.65045
[89] Williamson, J. H., Low-storage Runge-Kutta schemes, J. Comput. Phys., 35, 48-56 (1980) · Zbl 0425.65038
[90] Wolfram, S., The Mathematica Book (1996), Cambridge University Press: Cambridge University Press Cambridge · Zbl 0878.65001
[91] Research, Wolfram, Mathematica 3.0 Standard Add-On Packages (1996), Cambridge University Press: Cambridge University Press Cambridge · Zbl 0861.68047
[92] A.A. Wray, Minimal storage time advancement schemes for spectral methods, unpublished, NASA Ames Research Center, Moffett Field, CA, 1986; A.A. Wray, Minimal storage time advancement schemes for spectral methods, unpublished, NASA Ames Research Center, Moffett Field, CA, 1986
[93] Zennaro, M., Contractivity of Runge-Kutta methods with respect to forcing terms, Appl. Numer. Math., 11, 4, 321-345 (1993) · Zbl 0774.65054
[94] Zingg, D. W.; Chisholm, T. T., Runge-Kutta methods for linear ordinary differential equations, Appl. Numer. Math., 31, 2, 227-238 (1999) · Zbl 0935.65077
[95] Zingg, D. W.; Lomax, H.; Jurgens, H., High-accuracy finite-difference schemes for linear wave propagation, SIAM J. Sci. Comput., 17, 2, 328-346 (1996) · Zbl 0877.65063
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.