×

Application of the rosenbrock methods to the solution of unsteady 3D incompressible Navier-Stokes equations. (English) Zbl 1411.76059

Summary: We consider the Rosenbrock methods, namely a family of methods for differential algebraic equations, for the solution of the unsteady three-dimensional Navier-Stokes equations. These multistage schemes are attractive for non-linear problems because they achieve high order in time, ensuring stability properties and linearizing the system to be solved at each timestep. Moreover, as they provide inexpensive ways to estimate the local truncation error, adaptive timestep strategies can be easily devised. In this work, we test the Rosenbrock methods for the solution of three-dimensional unsteady incompressible flows. We derive the correct essential boundary conditions to impose at each stage in order to retain the convergence order of the schemes. Then, we consider two benchmark tests: a flow problem with imposed oscillatory pressure gradient whose analytical solution is known and the classical flow past a cylinder. In the latter case, we especially focus on the accuracy in the approximation of the drag and lift coefficients. In both benchmarks we test the performance of a time adaptivity scheme.

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
65L80 Numerical methods for differential-algebraic equations
76D05 Navier-Stokes equations for incompressible viscous fluids
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Atkinson, K.; Han, W.; Stewart, D. E., Numerical solution of ordinary differential equations, 108 (2011), John Wiley & Sons
[2] Bassi, F.; Botti, L.; Colombo, A.; Crivellini, A.; Ghidoni, A.; Massa, F., On the development of an implicit high-order discontinuous Galerkin method for DNS and implicit LES of turbulent flows, Eur J Mech B/Fluids, 55, 367-379 (2016) · Zbl 1408.76360
[3] Bassi, F.; Botti, L.; Colombo, A.; Ghidoni, A.; Massa, F., Linearly implicit Rosenbrock-type Runge-Kutta schemes applied to the Discontinuous Galerkin solution of compressible and incompressible unsteady flows, Comput Fluids, 118, 305-320 (2015) · Zbl 1390.76833
[4] Bayraktar, E.; Mierka, O.; Turek, S., Benchmark computations of 3D laminar flow around a cylinder with CFX, OpenFOAM and FeatFlow, Int J Comput Sci Eng, 7, 3, 253-266 (2012)
[5] Bertagna L., Deparis S., Formaggia L., Forti D., Veneziani A. The LifeV library: beyond the proof of concept. Submitted: arXiv:171006596; Bertagna L., Deparis S., Formaggia L., Forti D., Veneziani A. The LifeV library: beyond the proof of concept. Submitted: arXiv:171006596
[6] Blom, D. S.; Birken, P.; Bijl, H.; Kessels, F.; Meister, A.; van Zuijlen, A. H., A comparison of Rosenbrock and ESDIRK methods combined with iterative solvers for unsteady compressible flows, Adv Comput Math, 42, 6, 1401-1426 (2016) · Zbl 1388.76166
[7] Cellier, F. E.; Kofman, E., Continuous system simulation (2006), Springer Science & Business Media · Zbl 1112.93004
[8] Dedè, L., Optimal flow control for Navier-Stokes equations: drag minimization, Int J Numer Methods Fluids, 55, 4, 347-366 (2007) · Zbl 1388.76074
[9] Gustafsson, K.; Lundh, M.; Söderlind, G., Api stepsize control for the numerical solution of ordinary differential equations, BIT Numer Math, 28, 2, 270-287 (1988) · Zbl 0645.65039
[10] Guyon, E., Physical hydrodynamics (2001), Oxford University Press · Zbl 1012.76003
[11] Hairer, E.; Wanner, G., Solving ordinary differential equations. II, Springer Series in Computational Mathematics, 14 (2010), Springer-Verlag, Berlin · Zbl 1192.65097
[12] Ilchmann, A.; Reis, T., Surveys in differential-algebraic equations II (2014), Springer
[13] John, V., On the efficiency of linearization schemes and coupled multigrid methods in the simulation of a 3D flow around a cylinder, Int J Numer Methods Fluids, 50, 7, 845-862 (2006) · Zbl 1086.76039
[14] John, V.; Matthies, G.; Rang, J., A comparison of time-discretization/linearization approaches for the incompressible Navier-Stokes equations, Comput Methods Appl Mech Eng, 195, 44, 5995-6010 (2006) · Zbl 1124.76041
[15] John, V.; Rang, J., Adaptive time step control for the incompressible Navier-Stokes equations, Comput Methods Appl Mech Eng, 199, 9, 514-524 (2010) · Zbl 1227.76048
[16] Lang, J.; Teleaga, D., Towards a fully space-time adaptive FEM for magnetoquasistatics, IEEE Trans Magn, 44, 6, 1238-1241 (2008)
[17] Liu, X.; Xia, Y.; Luo, H.; Xuan, L., A comparative study of Rosenbrock-type and implicit Runge-Kutta time integration for discontinuous Galerkin method for unsteady 3D compressible Navier-Stokes equations, Commun Comput Phys, 20, 4, 1016-1044 (2016) · Zbl 1373.76096
[18] Massa, F.; Noventa, G.; Lorini, M.; Bassi, F.; Ghidoni, A., High-order linearly implicit two-step peer schemes for the discontinuous Galerkin solution of the incompressible Navier-Stokes equations, Comput Fluids, 162, 55-71 (2018) · Zbl 1390.76338
[19] Nørsett, S. P.; Wolfbrandt, A., Order conditions for Rosenbrock type methods, Numerische Mathematik, 32, 1, 1-15 (1979) · Zbl 0471.65044
[20] 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
[21] Quarteroni, A., Numerical models for differential problems, 8 (2014), Springer-Verlag
[22] Quarteroni, A.; Sacco, R.; Saleri, F., Numerical mathematics, 37 (2010), Springer Science & Business Media · Zbl 0913.65002
[23] Rang, J., Improved traditional Rosenbrock-Wanner methods for stiff ODEs and DAEs, J Comput Appl Math, 286, 128-144 (2015) · Zbl 1326.65085
[24] Rang, J., The Prothero and Robinson example: Convergence studies for Runge-Kutta and Rosenbrock-Wanner methods, Appl Numer Math, 108, 37-56 (2016) · Zbl 1346.65038
[25] Rang, J.; Angermann, L., New Rosenbrock W-methods of order 3 for partial differential algebraic equations of index 1, BIT Numer Math, 45, 4, 761-787 (2005) · Zbl 1093.65097
[26] Roche, M., Rosenbrock methods for differential algebraic equations, Numerishe Mathematik, 52, 1, 45-63 (1987) · Zbl 0613.65076
[27] Rosenbrock, H. H., Some general implicit processes for the numerical solution of differential equations, Comput J, 5, 4, 329-330 (1963) · Zbl 0112.07805
[28] Sanderse, B.; Koren, B., Accuracy analysis of explicit Runge-Kutta methods applied to the incompressible Navier-Stokes equations, J Comput Phys, 231, 8, 3041-3063 (2012) · Zbl 1402.65116
[29] Schäref, M.; Turek, S., Benchmark computations of laminar flow around a cylinder, Flow Simul High-Perform Comput II, 127, 553-563 (1996)
[30] Segal, A.; ur Rehman, M.; Vuik, C., Preconditioners for incompressible Navier-Stokes solvers, Numer Math, 3, 3, 245-275 (2010) · Zbl 1240.65098
[31] Shirokoff, D.; Rosales, R. R., An efficient method for the incompressible Navier-Stokes equations on irregular domains with no-slip boundary conditions, high order up to the boundary, J Comput Phys, 230, 23, 8619-8646 (2011) · Zbl 1426.76489
[32] Söderlind, G., Automatic control and adaptive time-stepping, Numer Algorithms, 31, 1-4, 281-310 (2002) · Zbl 1012.65080
[33] Söderlind, G., Digital filters in adaptive time-stepping, ACM Trans Math Softw (TOMS), 29, 1, 1-26 (2003) · Zbl 1097.93516
[34] Steinebach G. Order-reduction of ROW-methods for DAEs and method of lines applications; 1995. Preprint-Nr. 1741, FB Mathematik, TH Darmstadt.; Steinebach G. Order-reduction of ROW-methods for DAEs and method of lines applications; 1995. Preprint-Nr. 1741, FB Mathematik, TH Darmstadt.
[35] Verwer, J. G.; Spee, E. J.; Blom, J. G.; Hundsdorfer, W., A second-order Rosenbrock method applied to photochemical dispersion problems, SIAM J Sci Comput, 29, 4, 1456-1480 (1999) · Zbl 0928.65116
[36] Womersley, J. R., Method for the calculation of velocity, rate of flow and viscous drag in arteries when the pressure gradient is known, J Physiol (Lond), 127, 3, 553-563 (1955)
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.