×

Completed repeated Richardson extrapolation for compressible fluid flows. (English) Zbl 1443.76160

Summary: Richardson extrapolation is a powerful approach for reducing spatial discretization errors and increasing, in this way, the accuracy of the computed solution obtained by use of many numerical methods for solving different scientific and engineering problems. This approach has been used in a variety of computational fluid dynamics problems to reduce numerical errors, but its use has been restricted mainly to the computation of incompressible fluid flows and on grids with coincident nodes. In this work we present a completed repeated Richardson extrapolation (CRRE) procedure for a more generic type of grid not necessarily with coincident nodes, and test it on compressible fluid flows. Three tests are performed for one-dimensional and quasi-one-dimensional Euler equations: (i) Rayleigh flow, (ii) isentropic flow, and (iii) adiabatic flow through a nozzle. The last test involves a normal shock wave. To build a simple solver, these problems are solved by a first-order upwind-type finite difference method as the base scheme. The normal shock wave problem is also solved with a high-order weighted essentially nonoscillatory (WENO) scheme to compare it with the CRRE procedure. The procedure we propose can increase the achieved accuracy and significantly decrease the magnitude of the spatial error in all three tests. Its performance is best demonstrated in the Rayleigh flow test, where the spatial discretization error is reduced by seven orders of magnitude and the achieved accuracy is increased from 0.998 to 6.62 on a grid with 10,240 nodes. Similar performance is observed for isentropic flow, for which the spatial discretization error is reduced by nine orders of magnitude and the achieved accuracy is increased from 0.996 to 6.73 on a grid with 10,240 nodes. Finally, in adiabatic flow with a normal shock wave, the procedure can reduce the spatial discretization error both upstream and downstream of the shock. However, the more expensive high-order WENO scheme results in errors of lower magnitude upstream of the shock and a sharper shock transition for this shocked test case.

MSC:

76M20 Finite difference methods applied to problems in fluid mechanics
65B15 Euler-Maclaurin formula in numerical analysis
65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Hirsch, C., Numerical Computation of Iternal and External flows (2007), Elsevier
[2] Kaps, P.; Poon, S. W.H.; Bui, T. D., Rosenbrock methods for stiff ODEs: a comparison of Richardson extrapolation and embedding technique, Computing, 34, 17-40 (1985) · Zbl 0554.65054
[3] 4th International Symposium on Computational Wind Engineering (CWE2006).
[4] Richardson, L. F., The approximate arithmetical solution by finite differences of physical problems involving differential equations, with an application to the stresses in a masonry dam, Philos. Trans. R. Soc. Lond. Ser. A, Contain. Pap. Math. Phys. Character, 210, 307-357 (1910) · JFM 42.0873.02
[5] Richardson, L. F.; Gaunt, J. A., The deferred approach to the limit. Part I. Single lattice. Part II. Interpenetrating lattices, Philos. Trans. R. Soc. Lond. Ser. A Contain. Pap. Math. Phys. Character, 226, 299-361 (1927) · JFM 53.0432.02
[6] Roache, P. J.; Knupp, P. M., Completed Richardson extrapolation, Commun. Numer. Methods Eng., 9, 365-374 (1993) · Zbl 0782.65049
[7] Marchi, C.; Giacomini, F. F.; Santiago, C. D., Repeated Richardson extrapolation to reduce the field discretization error in computational fluid dynamics, Numer. Heat Transf. Part B, 70, 340-353 (2016)
[8] Martn-Vaquero, J.; Kleefeld, B., Extrapolated stabilized explicit Runge-Kutta methods, J. Comput. Phys., 326, 141-155 (2016) · Zbl 1422.65233
[9] Amore, P.; Boyd, J. P.; Fernndez, F. M.; Rsler, B., High order eigenvalues for the Helmholtz equation in complicated non-tensor domains through Richardson extrapolation of second order finite differences, J. Comput. Phys., 312, 252-271 (2016) · Zbl 1352.65458
[10] Dang, A.; Kehtarnavaz, H.; Coats, D., The use of Richardson extrapolation in PNs solutions of rocket nozzle flow, Proceedings of the 25th Joint Propulsion Conference (1989), American Institute of Aeronautics and Astronautics: American Institute of Aeronautics and Astronautics Monterey, Califrnia
[11] Marchi, C. H.; Martins, M. A.; Novak, L. A.; Araki, L. K.; Pinto, M. A.V.; de F. T. Gonalves, S.; Moro, D. F.; Freitas, I.d. S., Polynomial interpolation with repeated Richardson extrapolation to reduce discretization error in CFD, Appl. Math. Model., 40, 8872-8885 (2016) · Zbl 1480.65032
[12] Shu, C. W., Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws, (Quarteroni, A., Advanced Numerical Approximation of Nonlinear Hyperbolic Equations (1998), Springer: Springer Berlin), 325-432 · Zbl 0927.65111
[13] Gottlieb, S.; Gottlieb, D.; Shu, C. W., Recovering high-order accuracy in WENO computations of steady-state hyperbolic systems, J. Sci. Comput., 28, 307-318 (2006) · Zbl 1158.76365
[14] Marchi, C. H.; Novak, L. A.; Santiago, C. D.; Vargas, A. P.d. S., Highly accurate numerical solutions with repeated Richardson extrapolation for 2D Laplace equation, Appl. Math. Model., 37, 7386-7397 (2013) · Zbl 1426.65154
[15] Marchi, C. H.; Germer, E. M., Effect of ten CFD numerical schemes on repeated Richardson extrapolation (RRE), J. Appl. Comput. Math., 2, 1-9 (2013)
[16] Borges, R.; Carmona, M.; Costa, B.; Don, W. S., An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws, J. Comput. Phys., 227, 3191-3211 (2008) · Zbl 1136.65076
[17] Anderson, J. D., Modern Compressible Flow: With Historical Perspective (2003), McGraw-Hill: McGraw-Hill New York
[18] Tan, S.; Shu, C. W., Inverse Lax-Wendroff procedure for numerical boundary conditions of conservation laws, J. Comput. Phys., 229, 8144-8166 (2010) · Zbl 1198.65174
[19] Tan, S.; Wang, C.; Shu, C.-W.; Ning, J., Efficient implementation of high order inverse Lax-Wendroff boundary treatment for conservation laws, J. Comput. Phys., 231, 2510-2527 (2012) · Zbl 1430.65005
[20] Standard for Verification and Validation in Computational Fluid Dynamics and Heat Transfer, American Society of Mechanical Engineers, New York, 2009. https://www.asme.org/codes-standards/find-codes-standards/v-v-20-standard-verification-validation-computational-fluid-dynamics-heat-transfer#ASME-digital-books.
[21] Marchi, C. H.; Silva, A. F.C.d., Unidimensional numerical solution error estimation for convergent apparent order, Numer. Heat Transf., 42, 167-188 (2002)
[22] Roache, P. J., Fundamentals of Verification and Validation (2009), Hermosa Publishers: Hermosa Publishers New Mexico
[23] Don, W.-S.; Borges, R., Accuracy of the weighted essentially non-oscillatory conservative finite difference schemes, J. Comput. Phys., 250, 347-372 (2013) · Zbl 1349.65285
[24] Zhu, J.; Shu, C. W., Numerical study on the convergence to steady state solutions of a new class of high order WENO schemes, J. Comput. Phys., 349, 80-96 (2017) · Zbl 1380.76074
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.