×

Hp-adaptive time integration based on the BDF for viscous flows. (English) Zbl 1349.76217

Summary: This paper presents a procedure based on the Backward Differentiation Formulas of order 1 to 5 to obtain efficient time integration of the incompressible Navier-Stokes equations. The adaptive algorithm performs both stepsize and order selections to control respectively the solution accuracy and the computational efficiency of the time integration process. The stepsize selection (h-adaptivity) is based on a local error estimate and an error controller to guarantee that the numerical solution accuracy is within a user prescribed tolerance. The order selection (p-adaptivity) relies on the idea that low-accuracy solutions can be computed efficiently by low order time integrators while accurate solutions require high order time integrators to keep computational time low. The selection is based on a stability test that detects growing numerical noise and deems a method of order p stable if there is no method of lower order that delivers the same solution accuracy for a larger stepsize. Hence, it guarantees both that (1) the used method of integration operates inside of its stability region and (2) the time integration procedure is computationally efficient. The proposed time integration procedure also features a time-step rejection and quarantine mechanisms, a modified Newton method with a predictor and dense output techniques to compute solution at off-step points.

MSC:

76M10 Finite element methods applied to problems in fluid mechanics
65M60 Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs
65L06 Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations
65M20 Method of lines for initial value and initial-boundary value problems involving PDEs
76N15 Gas dynamics (general theory)

Software:

RODAS; DIFSUB; PARDISO; DASSL
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Hairer, E.; Nørsett, S.; Wanner, G., Solving Ordinary Differential Equations, Volume 1: Nonstiff Problems (2002), Springer: Springer Germany
[2] Hairer, E.; Wanner, G., Solving Ordinary Differential Equations, Volume 2: Stiff and Differential-Algebraic Problems (2002), Springer: Springer Germany
[3] Asher, U. M.; Petzold, L. R., Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations (1998), Society for Industrial and Applied Mathematics (SIAM): Society for Industrial and Applied Mathematics (SIAM) Philadelphia, Pennsylvania
[4] Gresho, P. M.; Sani, R. L., Incompressible Flow and the Finite-Element Method: Volume 1, Advection-Diffusion (1998), Wiley: Wiley England · Zbl 0941.76002
[5] Jansen, K.; Whiting, C.; Hulbert, G., A generalized-\(α\) method for integrating the filtered Navier-Stokes equations with a stabilized finite-element method, Comput. Methods Appl. Mech. Eng., 190, 3, 305-319 (2000) · Zbl 0973.76048
[6] Jothiprasad, G.; Mavriplis, D. J.; Caughey, D. A., Higher-order time integration schemes for the unsteady Navier-Stokes equations on unstructured meshes, J. Comput. Phys., 191, 542-566 (2003) · Zbl 1134.76428
[7] Montlaur, A.; Fernandez-Mendez, S.; Huerta, A., High-order time integration for unsteady incompressible flows, Int. J. Numer. Methods Fluids, 70, 603-626 (2012) · Zbl 1412.65085
[8] Hay, A.; Yu, K. R.; Etienne, S.; Garon, A.; Pelletier, D., High-order temporal accuracy for 3D finite-element ALE flow simulations, Comput. Fluids, 100, 204-217 (2014) · Zbl 1391.76335
[9] Sanderse, B., Energy-conserving Runge-Kutta methods for the incompressible Navier-Stokes equations, J. Comput. Phys., 233, 100-131 (2013) · Zbl 1286.76034
[10] Etienne, S.; Garon, A.; Pelletier, D., Perspective on the geometric conservation law and finite element methods for ALE simulations of incompressible flows, J. Comput. Phys., 228, 2313-2333 (2009) · Zbl 1275.76159
[11] Gear, C. W., Numerical Initial Value Problems in Ordinary Differential Equations (1971), Prentice Hall · Zbl 0217.21701
[12] Malidi, A.; Dufour, S.; N’dri, D., A study of time integration schemes for the numerical modelling of free surface flows, Int. J. Numer. Methods Fluids, 48, 1123-1147 (2005) · Zbl 1071.76030
[13] Kay, D.; Gresho, P.; Griffiths, D.; Silverster, D., Adaptive time-stepping for incompressible flow part II: Navier-Stokes equations, SIAM J. Sci. Comput., 32, 1, 111-128 (2010) · Zbl 1410.76293
[14] John, V.; Rang, J., Adaptive time step control for the incompressible Navier-Stokes equations, Comput. Methods Appl. Mech. Eng., 199, 514-524 (2010) · Zbl 1227.76048
[15] Gear, C. W., The automatic integration of ordinary differential equations, Commun. ACM, 14, 176-179 (1971) · Zbl 0217.21701
[16] Gear, C. W., DIFSUB for solution of ordinary differential equations, Algorithm 407, Commun. ACM, 14, 185-190 (1971)
[17] Skelboe, S., The control of order and steplength for backward differentiation methods, BIT Numer. Math., 17, 91-107 (1977) · Zbl 0361.65062
[18] Brenan, K. E.; Campbell, S. L.; Petzold, L. R., Numerical Solution of Initial-Value Problems in Differential-Algebraic Equations (1989), Elsevier: Elsevier New York · Zbl 0699.65057
[19] Shampine, L. F.; Gordon, M. K., Computer Solution of Ordinary Differential Equations (1975), W.H. Freeman: W.H. Freeman San Fransisco, California · Zbl 0347.65001
[20] Hindmarsh, A. C.; Byrne, G. D., A polyalgorithm for the numerical solution of ordinary differential equations, ACM Trans. Math. Softw., 1, 71-96 (1975) · Zbl 0311.65049
[21] Hindmarsh, A. C., ODE solvers for use with the method of lines, (Vichnevetsky, R.; Stepleman, R., Advances in Computer Methods for Partial Differential Equations - IV (1981), New Brunswick: New Brunswick New Jersey), 312-316
[22] Petzold, L. R., A description of DASSL: a differential/algebraic system solver (1982), Applied Mathematics Division, Sandia National Laboratories: Applied Mathematics Division, Sandia National Laboratories Livermore, California, Tech. Rep. SAND82-8637
[23] Jackson, K. R.; Sacks-Davis, R., An alternative implementation of variable step-size multistep formulas for stiff ODEs, ACM Trans. Math. Softw., 6, 295-318 (1980) · Zbl 0434.65046
[24] Dettmer, W.; Perić, D., An analysis of the time integration algorithms for the finite element solutions of incompressible Navier-Stokes equations based on the stabilized formulation, Comput. Methods Appl. Mech. Eng., 192, 1177-1226 (2003) · Zbl 1091.76521
[25] Grigorieff, R. D., Stability of multistep-methods on variable grids, Numer. Math., 42, 3, 359-377 (1983) · Zbl 0554.65051
[26] Krogh, F. T., A test for instability in the numerical solution of ordinary differential equations, J. ACM, 14, 2, 351-354 (1967) · Zbl 0189.16502
[27] Gustafsson, K.; Lundh, M.; Söderlind, G., A PI stepsize control for the numerical solution of ordinary differential equations, BIT Numer. Math., 28, 270-287 (1988) · Zbl 0645.65039
[28] Gustafsson, K., Control-theoretic techniques for stepsize selection in implicit Runge-Kutta methods, ACM Trans. Math. Softw., 20, 496-517 (1994) · Zbl 0888.65096
[29] Pelletier, D., Adaptive finite element computations of complex flows, Int. J. Numer. Methods Fluids, 31, 189-202 (1999) · Zbl 0982.76063
[30] Braza, M.; Chassaing, P.; Minh, H. H., Numerical study and physical analysis of the pressure and velocity fields in the near wake of a circular cylinder, J. Fluid Mech., 165, 79-130 (1986) · Zbl 0596.76047
[31] Schenk, O.; Gärtner, K., Solving unsymmetric sparse systems of linear equations with PARDISO, Future Gener. Comput. Syst., 20, 3, 475-487 (2004) · Zbl 1062.65035
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.