×

A note on the accuracy of the generalized-\(\alpha\) scheme for the incompressible Navier-Stokes equations. (English) Zbl 07863084

Summary: We investigate the temporal accuracy of two generalized-\(\alpha\) schemes for the incompressible Navier-Stokes equations. In a widely-adopted approach, the pressure is collocated at the time step \(t_{n+1}\) while the remainder of the Navier-Stokes equations is discretized following the generalized-\(\alpha\) scheme. That scheme has been claimed to be second-order accurate in time. We developed a suite of numerical code using inf-sup stable higher-order non-uniform rational B-spline (NURBS) elements for spatial discretization. In doing so, we are able to achieve high spatial accuracy and to investigate asymptotic temporal convergence behavior. Numerical evidence suggests that only first-order accuracy is achieved, at least for the pressure, in this aforesaid temporal discretization approach. On the other hand, evaluating the pressure at the intermediate time step \(t_{n+\alpha_f}\) recovers second-order accuracy, and the numerical implementation is simplified. We recommend this second approach as the generalized-\(\alpha\) scheme of choice when integrating the incompressible Navier-Stokes equations.
{© 2020 John Wiley & Sons Ltd}

MSC:

76Mxx Basic methods in fluid mechanics
74Sxx Numerical and other methods in solid mechanics
65Nxx Numerical methods for partial differential equations, boundary value problems

Software:

GitHub

References:

[1] ChungJ, HulbertGM. A time integration algorithm for structural dynamics with improved numerical dissipation: the generalized‐α method. J Appl Mech. 1993;60:371‐375. · Zbl 0775.73337
[2] HilberHM, HughesTJR. Collocation, dissipation and ’overshoot’ for time integration schemes in structural dynamics. Earthq Eng Struct Dyn. 1978;6:99‐117.
[3] JansenKE, WhitingCH, HulbertGM. A generalized‐α method for integrating the filtered Navier‐Stokes equations with a stabilized finite element method. Comput Methods Appl Mech Eng. 2000;190:305‐319. · Zbl 0973.76048
[4] KadapaC, DettmerWG, PerićD. On the advantages of using the first‐order generalised‐alpha scheme for structural dynamic problems. Comput Struct. 2017;193:226‐238.
[5] svSolver GitHub repository; 2019. https://github.com/SimVascular/svSolver Accessed August 08, 2019.
[6] BazilevsY, CaloVM, CottrellJA, HughesTJR, RealiA, ScovazziG. Variational multiscale residual‐based turbulence modeling for large eddy simulation of incompressible flows. Comput Methods Appl Mech Eng. 2007;197:173‐201. · Zbl 1169.76352
[7] BazilevsY, CaloVM, ZhangY, HughesTJR. Isogeometric fluid‐structure interaction analysis with applications to arterial blood flow. Comput Mech. 2006;38:310‐322. · Zbl 1161.74020
[8] BazilevsY, TakizawaK, TezduyarTE. Computational Fluid‐Structure Interaction: Methods and Applications. Hoboke, NJ: John Wiley & Sons; 2013. · Zbl 1286.74001
[9] CasqueroH, Bona‐CasasC, ToshniwalD, HughesTJR, GomezH, ZhangJY. The divergence‐conforming immersed boundary method: application to vesicle and capsule dynamics; 2020. arXiv preprint arXiv:2001.08244.
[10] CasqueroH, LiuL, Bona‐CasasC, ZhangY, GomezH. A hybrid variational‐collocation immersed method for fluid‐structure interaction using unstructured T‐splines. Int J Numer Methods Eng. 2016;105(11):855‐880. · Zbl 07868662
[11] FigueroaCA. A Coupled‐Momentum Method to Model Blow Flow and Vessel Deformation in Human Arteries: Applications in Disease Research and Simulation‐Based Medical Planning [PhD thesis]. Stanford University; 2006.
[12] GamnitzerP, GravemeierV, WallWA. Time‐dependent subgrid scales in residual‐based large eddy simulation of turbulent channel flow. Comput Methods Appl Mech Eng. 2010;199:819‐827. · Zbl 1406.76026
[13] JoshiV, JaimanRK. A hybrid variational Allen‐Cahn/ALE scheme for the coupled analysis of two‐phase fluid‐structure interaction. Int J Numer Methods Eng. 2019;117(4):405‐429. · Zbl 07865252
[14] KangS, ChoiHG, YooJY. Investigation of fluid‐structure interactions using a velocity‐linked P2/P1 finite element method and the generalized‐α method. Int J Numer Methods Eng. 2012;90(12):1529‐1548. · Zbl 1246.74059
[15] PetersEL, EvansJA. A divergence‐conforming hybridized discontinuous Galerkin method for the incompressible Reynolds‐averaged Navier‐Stokes equations. Int J Numer Methods Fluids. 2019;91(3):112‐133.
[16] Vignon‐ClementelIE, FigueroaCA, JansenKE, TaylorCA. Outflow boundary conditions for three‐dimensional finite element modeling of blood flow and pressure in arteries. Comput Methods Appl Mech Eng. 2006;195:3776‐3796. · Zbl 1175.76098
[17] YanJ, YanW, LinS, WagnerGJ. A fully coupled finite element formulation for liquid‐solid‐gas thermo‐fluid flow with melting and solidification. Comput Methods Appl Mech Eng. 2018;336:444‐470. · Zbl 1440.74266
[18] ZhuQ, XuandF, XuS, HsuMC, YanJ. An immersogeometric formulation for free‐surface flows with application to marine engineering problems. Comput Methods Appl Mech Eng. 2020;361:112748. · Zbl 1442.76019
[19] KimHJ. Three‐Dimensional Finite Element Modeling of Blood Flow in the Coronary Arteries [PhD thesis]. Stanford University; 2009.
[20] LiuJ, YangW, DongM, MarsdenAL. The nested block preconditioning technique for the incompressible Navier‐Stokes equations with emphasis on hemodynamic simulations. Comput Methods Appl Mech Eng. 2020;367:113122. · Zbl 1442.76045
[21] Vignon‐ClementelIE. A Coupled Multidomain Method for Computational Modeling of Blood Flow [PhD thesis]. Stanford University; 2006.
[22] MoghadamME, Vignon‐ClementelIE, FigliolaR, MarsdenAL. Modeling of congenital hearts alliance (MOCHA) investigators. a modular numerical method for implicit 0D/3D coupling in cardiovascular finite element simulations. J Comput Phys. 2013;244:63‐79. · Zbl 1377.76041
[23] CodinaR, PrincipeJ, GuaschO, BadiaS. Time dependent subscales in the stabilized finite element approximation of incompressible flow problems. Comput Methods Appl Mech Eng. 2007;196:2413‐2430. · Zbl 1173.76335
[24] JodlbauerD, LangerU, WickT. Parallel block‐preconditioned monolithic solvers for fluid‐structure interaction problems. Int J Numer Methods Eng. 2019;117(6):623‐643. · Zbl 07865264
[25] SaksonoPH, DettmerWG, PerićD. An adaptive remeshing strategy for flows with moving boundariesand fluid‐structure interaction. Int J Numer Methods Eng. 2007;71:1009‐1050. · Zbl 1194.76140
[26] RannacherR. Methods for Numerical Flow Simulation. Basel: Birkhäuser; 2008:275‐332.
[27] GravemeierV, KronbichlerM, GeeMW, WallWA. An algebraic variational multiscale‐multigrid method for large‐eddy simulation: generalized‐α time integration, Fourier analysis and application to turbulent flow past a square‐section cylinder. Comput Mech. 2011;47:217‐233. · Zbl 1398.76105
[28] KadapaC, DettmerWG, PerićD. Accurate iteration‐free mixed‐stabilised formulation for laminar incompressible Navier‐Stokes: applications to fluid‐structure interaction; 2020. arXiv preprint arXiv:2001.04925.
[29] LiuJ, MarsdenAL. A unified continuum and variational multiscale formulation for fluids, solids, and fluid‐structure interaction. Comput Methods Appl Mech Eng. 2018;337:549‐597. · Zbl 1440.74013
[30] RossiS, AbboudN, ScovazziG. Implicit finite incompressible elastodynamics with linear finite elements: a stabilized method in rate form. Comput Methods Appl Mech Eng. 2016;311:208‐249. · Zbl 1439.74464
[31] GroppWD, KaushikDK, KeyesDE, SmithBF. Performance modeling and tuning of an unstructured mesh CFD application. Paper presented at: Proceedings of the 2000 ACM/IEEE Conference on Supercomputing SC’00; Dallas, TX, 2000.
[32] HsuMC, BazilevsY, CaloVM, TezduyarTE, HughesTJR. Improving stability of stabilized and multiscale formulations in flow simulations at small time steps. Comput Methods Appl Mech Eng. 2010;19:828‐840. · Zbl 1406.76028
[33] HughesTJR, CottrellJA, BazilevsY. Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. Comput Methods Appl Mech Eng. 2005;194:4135‐4195. · Zbl 1151.74419
[34] DettmerW, PerićD. An analysis of the time integration algorithms for the finite element solutions of incompressible Navier‐Stokes equations based on a stabilised formulation. Comput Methods Appl Mech Eng. 2003;192:1177‐1226. · Zbl 1091.76521
[35] BazilevsY, MichlerC, CaloVM, HughesTJR. Isogeometric variational multiscale modeling of wall‐bounded turbulent flows with weakly enforced boundary conditions on unstretched meshes. Comput Methods Appl Mech Eng. 2010;199:780‐790. · Zbl 1406.76023
[36] LiuJ, YangW, LanIS, MarsdenAL. Fluid‐structure interaction modeling of blood flow in the pulmonary arteries using the unified continuum and variational multiscale formulation. Mech Res Commun. 2020;107:103556.
[37] EthierCR, SteinmanDA. Exact fully 3D Navier‐Stokes solutions for benchmarking. Int J Numer Methods Fluids. 1994;19:369‐375. · Zbl 0814.76031
[38] WomersleyJR. Method for the calculation of velocity, rate of flow and viscous drag in arteries when the pressure gradient is known. J Physiol. 1955;127:553‐563.
[39] TakacsT, JüttlerB. Existence of stiffness matrix integrals for singularly parameterized domains in isogeometric analysis. Comput Methods Appl Mech Eng. 2011;200:3568‐3582. · Zbl 1239.65014
[40] GreshoPM, SaniRL. Incompressible Flow and the Finite Element Method. Volume 1: Advection‐Diffusion and Isothermal Laminar Flow. Hoboken, NJ: John Wiley & Sons Ltd; 1998. · Zbl 0941.76002
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.