Artificial Compressibility Methods for the Incompressible Navier–Stokes Equations Using Lowest-Order Face-Based Schemes on Polytopal Meshes (original) (raw)

Abstract

We investigate artificial compressibility (AC) techniques for the time discretization of the incompressible Navier–Stokes equations. The space discretization is based on a lowest-order face-based scheme supporting polytopal meshes, namely discrete velocities are attached to the mesh faces and cells, whereas discrete pressures are attached to the mesh cells. This face-based scheme can be embedded into the framework of hybrid mixed mimetic schemes and gradient schemes, and has close links to the lowest-order version of hybrid high-order methods devised for the steady incompressible Navier–Stokes equations. The AC time-stepping uncouples at each time step the velocity update from the pressure update. The performances of this approach are compared against those of the more traditional monolithic approach which maintains the velocity-pressure coupling at each time step. We consider both first-order and second-order time schemes and either an implicit or an explicit treatment of the nonlinear convection term. We investigate numerically the CFL stability restriction resulting from an explicit treatment, both on Cartesian and polytopal meshes. Finally, numerical tests on large 3D polytopal meshes highlight the efficiency of the AC approach and the benefits of using second-order schemes whenever accurate discrete solutions are to be attained.

Award Identifier / Grant number: 2017/0526

Funding statement: The PhD fellowship of R. Milani was partially supported by EDF R&D and ANRT (No. 2017/0526).

References

[1] J. Aghili, S. Boyaval and D. A. Di Pietro, Hybridization of mixed high-order methods on general meshes and application to the Stokes equations, Comput. Methods Appl. Math. 15 (2015), no. 2, 111–134.10.1515/cmam-2015-0004Search in Google Scholar

[2] P. R. Amestoy, I. S. Duff, J.-Y. L’Excellent and J. Koster, A fully asynchronous multifrontal solver using distributed dynamic scheduling, SIAM J. Matrix Anal. Appl. 23 (2001), no. 1, 15–41.10.1137/S0895479899358194Search in Google Scholar

[3] F. Archambeau, N. Méchitoua and M. Sakiz, Code Saturne: A finite volume code for the computation of turbulent incompressible flows – industrial applications, Int. J. Finite Vol. 1 (2004), no. 1, 1–62.Search in Google Scholar

[4] M. Arioli, Generalized Golub–Kahan bidiagonalization and stopping criteria, SIAM J. Matrix Anal. Appl. 34 (2013), no. 2, 571–592.10.1137/120866543Search in Google Scholar

[5] M. Arioli, C. Kruse, R. Ulrich and N. Tardieu, An iterative generalized Golub–Kahan algorithm for problems in structural mechanics, Technical Report, 2018.Search in Google Scholar

[6] D. N. Arnold, R. S. Falk and R. Winther, Preconditioning in H ⁢ ( div ) and applications, Math. Comp. 66 (1997), no. 219, 957–984.10.1090/S0025-5718-97-00826-0Search in Google Scholar

[7] F. Bassi and S. Rebay, A high-order accurate discontinuous finite element method for the numerical solution of the compressible Navier–Stokes equations, J. Comput. Phys. 131 (1997), no. 2, 267–279.10.1006/jcph.1996.5572Search in Google Scholar

[8] L. Beirão da Veiga, V. Gyrya, K. Lipnikov and G. Manzini, Mimetic finite difference method for the Stokes problem on polygonal meshes, J. Comput. Phys. 228 (2009), no. 19, 7215–7232.10.1016/j.jcp.2009.06.034Search in Google Scholar

[9] L. Beirão da Veiga, C. Lovadina and G. Vacca, Divergence free virtual elements for the Stokes problem on polygonal meshes, ESAIM Math. Model. Numer. Anal. 51 (2017), no. 2, 509–535.10.1051/m2an/2016032Search in Google Scholar

[10] L. Beirão da Veiga, C. Lovadina and G. Vacca, Virtual elements for the Navier–Stokes problem on polygonal meshes, SIAM J. Numer. Anal. 56 (2018), no. 3, 1210–1242.10.1137/17M1132811Search in Google Scholar

[11] M. Benzi, M. A. Olshanskii and Z. Wang, Modified augmented Lagrangian preconditioners for the incompressible Navier–Stokes equations, Internat. J. Numer. Methods Fluids 66 (2011), no. 4, 486–508.10.1002/fld.2267Search in Google Scholar

[12] J. Bonelle, Compatible Discrete Operator schemes on polyhedral meshes for elliptic and Stokes equations, PhD thesis, Université Paris-Est, 2014.Search in Google Scholar

[13] J. Bonelle, D. A. Di Pietro and A. Ern, Low-order reconstruction operators on polyhedral meshes: Application to compatible discrete operator schemes, Comput. Aided Geom. Design 35/36 (2015), 27–41.10.1016/j.cagd.2015.03.015Search in Google Scholar

[14] J. Bonelle and A. Ern, Analysis of compatible discrete operator schemes for elliptic problems on polyhedral meshes, ESAIM Math. Model. Numer. Anal. 48 (2014), no. 2, 553–581.10.1051/m2an/2013104Search in Google Scholar

[15] J. Bonelle and A. Ern, Analysis of compatible discrete operator schemes for the Stokes equations on polyhedral meshes, IMA J. Numer. Anal. 35 (2015), no. 4, 1672–1697.10.1093/imanum/dru051Search in Google Scholar

[16] J. Bonelle, A. Ern and R. Milani, Compatible discrete operator schemes for the steady incompressible Stokes and Navier–Stokes equations, Finite Volumes for Complex Applications IX - Methods, Theoretical Aspects, Examples, Springer Proc. Math. Stat. 323, Springer, Cham (2020), 93–101.10.1007/978-3-030-43651-3_6Search in Google Scholar

[17] L. Botti, D. A. Di Pietro and J. Droniou, A hybrid high-order method for the incompressible Navier–Stokes equations based on Temam’s device, J. Comput. Phys. 376 (2019), 786–816.10.1016/j.jcp.2018.10.014Search in Google Scholar

[18] A. L. Bowers, S. Le Borne and L. G. Rebholz, Error analysis and iterative solvers for Navier–Stokes projection methods with standard and sparse grad-div stabilization, Comput. Methods Appl. Mech. Engrg. 275 (2014), 1–19.10.1016/j.cma.2014.02.021Search in Google Scholar

[19] F. Boyer, S. Krell and F. Nabet, Inf-Sup stability of the discrete duality finite volume method for the 2D Stokes problem, Math. Comp. 84 (2015), no. 296, 2705–2742.10.1090/mcom/2956Search in Google Scholar

[20] C. Cancès and P. Omnes, Finite Volumes for Complex Applications VIII – Methods and Theoretical Aspects, Springer Proc. Math. Stat. 199, Springer, Cham, 2017.10.1007/978-3-319-57397-7Search in Google Scholar

[21] A. Cangiani, V. Gyrya and G. Manzini, The nonconforming virtual element method for the Stokes equations, SIAM J. Numer. Anal. 54 (2016), no. 6, 3411–3435.10.1137/15M1049531Search in Google Scholar

[22] P. Cantin, Approximation of scalar and vector transport problems on polyhedral meshes, PhD thesis, Université Paris-Est, 2016.Search in Google Scholar

[23] P. Cantin and A. Ern, An edge-based scheme on polyhedral meshes for vector advection-reaction equations, ESAIM Math. Model. Numer. Anal. 51 (2017), no. 5, 1561–1581.10.1051/m2an/2016075Search in Google Scholar

[24] S. Charnyi, T. Heister, M. A. Olshanskii and L. G. Rebholz, On conservation laws of Navier–Stokes Galerkin discretizations, J. Comput. Phys. 337 (2017), 289–308.10.1016/j.jcp.2017.02.039Search in Google Scholar

[25] A. J. Chorin, Numerical solution of the Navier–Stokes equations, Math. Comp. 22 (1968), 745–762.10.1090/S0025-5718-1968-0242392-2Search in Google Scholar

[26] B. Cockburn, G. Kanschat, D. Schötzau and C. Schwab, Local discontinuous Galerkin methods for the Stokes system, SIAM J. Numer. Anal. 40 (2002), no. 1, 319–343.10.1137/S0036142900380121Search in Google Scholar

[27] B. Cockburn, N. C. Nguyen and J. Peraire, A comparison of HDG methods for Stokes flow, J. Sci. Comput. 45 (2010), no. 1–3, 215–237.10.1007/s10915-010-9359-0Search in Google Scholar

[28] V. DeCaria, W. Layton and M. McLaughlin, A conservative, second order, unconditionally stable artificial compression method, Comput. Methods Appl. Mech. Engrg. 325 (2017), 733–747.10.1016/j.cma.2017.07.033Search in Google Scholar

[29] S. Delcourte, Développement de méthodes de volumes finis pour la mécanique des fluides, PhD thesis, Université Paul Sabatier, 2007.Search in Google Scholar

[30] D. A. Di Pietro, J. Droniou and A. Ern, A discontinuous-skeletal method for advection-diffusion-reaction on general meshes, SIAM J. Numer. Anal. 53 (2015), no. 5, 2135–2157.10.1137/140993971Search in Google Scholar

[31] D. A. Di Pietro and A. Ern, A hybrid high-order locking-free method for linear elasticity on general meshes, Comput. Methods Appl. Mech. Engrg. 283 (2015), 1–21.10.1016/j.cma.2014.09.009Search in Google Scholar

[32] D. A. Di Pietro, A. Ern, A. Linke and F. Schieweck, A discontinuous skeletal method for the viscosity-dependent Stokes problem, Comput. Methods Appl. Mech. Engrg. 306 (2016), 175–195.10.1016/j.cma.2016.03.033Search in Google Scholar

[33] D. A. Di Pietro and S. Krell, A hybrid high-order method for the steady incompressible Navier–Stokes problem, J. Sci. Comput. 74 (2018), no. 3, 1677–1705.10.1007/s10915-017-0512-xSearch in Google Scholar

[34] D. A. Di Pietro and S. Lemaire, An extension of the Crouzeix–Raviart space to general meshes with application to quasi-incompressible linear elasticity and Stokes flow, Math. Comp. 84 (2015), no. 291, 1–31.10.1090/S0025-5718-2014-02861-5Search in Google Scholar

[35] J. Droniou and R. Eymard, Study of the mixed finite volume method for Stokes and Navier–Stokes equations, Numer. Methods Partial Differential Equations 25 (2009), no. 1, 137–171.10.1002/num.20333Search in Google Scholar

[36] J. Droniou and R. Eymard, Benchmark: Two hybrid mimetic mixed schemes for the lid-driven cavity, Finite Volumes for Complex Applications VIII – Methods and Theoretical Aspects, Springer Proc. Math. Stat. 199, Springer, Cham (2017), 107–124.10.1007/978-3-319-57397-7_8Search in Google Scholar

[37] J. Droniou, R. Eymard and P. Feron, Gradient schemes for Stokes problem, IMA J. Numer. Anal. 36 (2016), no. 4, 1636–1669.10.1093/imanum/drv061Search in Google Scholar

[38] J. Droniou, R. Eymard, T. Gallouet and R. Herbin, Gradient schemes: A generic framework for the discretisation of linear, nonlinear and nonlocal elliptic and parabolic equations, Math. Models Methods Appl. Sci. 23 (2013), no. 13, 2395–2432.10.1142/S0218202513500358Search in Google Scholar

[39] A. Ern and J.-L. Guermond, Finite Elements. III, Texts Appl. Math. 74, Springer, Cham, 2021.10.1007/978-3-030-57348-5Search in Google Scholar

[40] R. Eymard, P. Feron and C. Guichard, Family of convergent numerical schemes for the incompressible Navier–Stokes equations, Math. Comput. Simulation 144 (2018), 196–218.10.1016/j.matcom.2017.08.003Search in Google Scholar

[41] R. Eymard, T. Gallouët and R. Herbin, Discretization of heterogeneous and anisotropic diffusion problems on general nonconforming meshes SUSHI: A scheme using stabilization and hybrid interfaces, IMA J. Numer. Anal. 30 (2010), no. 4, 1009–1043.10.1093/imanum/drn084Search in Google Scholar

[42] P. E. Farrell, L. Mitchell and F. Wechsung, An augmented Lagrangian preconditioner for the 3D stationary incompressible Navier–Stokes equations at high Reynolds number, SIAM J. Sci. Comput. 41 (2019), no. 5, A3073–A3096.10.1137/18M1219370Search in Google Scholar

[43] P. Feron, Gradient Schemes for some elliptic and parabolic, linear and non-linear problems, PhD thesis, Université Paris-Est, 2016.Search in Google Scholar

[44] J. Fořt, J. Fürst, J. Halama, R. Herbin and F. Hubert, Finite Volumes for Complex Applications VI – Problems & Perspectives, Springer Proc. Math. Stat. 4, Springer, Berlin, 2011.10.1007/978-3-642-20671-9Search in Google Scholar

[45] G. N. Gatica, M. Munar and F. A. Sequeira, A mixed virtual element method for the Navier–Stokes equations, Math. Models Methods Appl. Sci. 28 (2018), no. 14, 2719–2762.10.1142/S0218202518500598Search in Google Scholar

[46] T. Goudon, S. Krell and G. Lissoni, DDFV method for Navier–Stokes problem with outflow boundary conditions, Numer. Math. 142 (2019), no. 1, 55–102.10.1007/s00211-018-1014-ySearch in Google Scholar

[47] J.-L. Guermond and P. Minev, High-order time stepping for the incompressible Navier–Stokes equations, SIAM J. Sci. Comput. 37 (2015), no. 6, A2656–A2681.10.1137/140975231Search in Google Scholar

[48] J.-L. Guermond and P. Minev, High-order time stepping for the Navier–Stokes equations with minimal computational complexity, J. Comput. Appl. Math. 310 (2017), 92–103.10.1016/j.cam.2016.04.033Search in Google Scholar

[49] J.-L. Guermond and P. Minev, High-order adaptive time stepping for the incompressible Navier–Stokes equations, SIAM J. Sci. Comput. 41 (2019), no. 2, A770–A788.10.1137/18M1209301Search in Google Scholar

[50] J. L. Guermond, P. Minev and J. Shen, An overview of projection methods for incompressible flows, Comput. Methods Appl. Mech. Engrg. 195 (2006), no. 44–47, 6011–6045.10.1016/j.cma.2005.10.010Search in Google Scholar

[51] P. Hansbo and M. G. Larson, Discontinuous Galerkin methods for incompressible and nearly incompressible elasticity by Nitsche’s method, Comput. Methods Appl. Mech. Engrg. 191 (2002), no. 17–18, 1895–1908.10.1016/S0045-7825(01)00358-9Search in Google Scholar

[52] T. Heister and G. Rapin, Efficient augmented Lagrangian-type preconditioning for the Oseen problem using grad-div stabilization, Internat. J. Numer. Methods Fluids 71 (2013), no. 1, 118–134.10.1002/fld.3654Search in Google Scholar

[53] J. G. Heywood and R. Rannacher, Finite element approximation of the nonstationary Navier–Stokes problem. I. Regularity of solutions and second-order error estimates for spatial discretization, SIAM J. Numer. Anal. 19 (1982), no. 2, 275–311.10.1137/0719018Search in Google Scholar

[54] J. G. Heywood and R. Rannacher, Finite element approximation of the nonstationary Navier–Stokes problem. III. Smoothing property and higher order error estimates for spatial discretization, SIAM J. Numer. Anal. 25 (1988), no. 3, 489–512.10.1137/0725032Search in Google Scholar

[55] Y. Jeon, E.-J. Park and D. Sheen, A hybridized finite element method for the Stokes problem, Comput. Math. Appl. 68 (2014), no. 12, 2222–2232.10.1016/j.camwa.2014.08.005Search in Google Scholar

[56] M. Jobelin, C. Lapuerta, J.-C. Latché, P. Angot and B. Piar, A finite element penalty-projection method for incompressible flows, J. Comput. Phys. 217 (2006), no. 2, 502–518.10.1016/j.jcp.2006.01.019Search in Google Scholar

[57] S. Krell and G. Manzini, The discrete duality finite volume method for Stokes equations on three-dimensional polyhedral meshes, SIAM J. Numer. Anal. 50 (2012), no. 2, 808–837.10.1137/110831593Search in Google Scholar

[58] O. A. Ladyzhenskaya, The Mathematical Theory of Viscous Incompressible Flow, Math. Appl. 2, Gordon and Breach, New York, 1969.Search in Google Scholar

[59] W. Layton and M. McLaughlin, Doubly-adaptive artificial compression methods for incompressible flow, J. Numer. Math. 28 (2020), no. 3, 179–196.10.1515/jnma-2019-0015Search in Google Scholar

[60] C. Lehrenfeld and J. Schöberl, High order exactly divergence-free hybrid discontinuous Galerkin methods for unsteady incompressible flows, Comput. Methods Appl. Mech. Engrg. 307 (2016), 339–361.10.1016/j.cma.2016.04.025Search in Google Scholar

[61] R. Milani, Compatible discrete operator schemes for the unsteady incompressible Navier–Stokes equations, PhD thesis, Université Paris-Est, 2020.Search in Google Scholar

[62] L. Mu, X. Wang and X. Ye, A modified weak Galerkin finite element method for the Stokes equations, J. Comput. Appl. Math. 275 (2015), 79–90.10.1016/j.cam.2014.08.006Search in Google Scholar

[63] N. C. Nguyen, J. Peraire and B. Cockburn, An implicit high-order hybridizable discontinuous Galerkin method for the incompressible Navier–Stokes equations, J. Comput. Phys. 230 (2011), no. 4, 1147–1170.10.1016/j.jcp.2010.10.032Search in Google Scholar

[64] Y. Notay and P. S. Vassilevski, Recursive Krylov-based multigrid cycles, Numer. Linear Algebra Appl. 15 (2008), no. 5, 473–487.10.1002/nla.542Search in Google Scholar

[65] T. Ohwada and P. Asinari, Artificial compressibility method revisited: asymptotic numerical method for incompressible Navier–Stokes equations, J. Comput. Phys. 229 (2010), no. 5, 1698–1723.10.1016/j.jcp.2009.11.003Search in Google Scholar

[66] M. Olshanskii, G. Lube, T. Heister and J. Löwe, Grad-div stabilization and subgrid pressure models for the incompressible Navier–Stokes equations, Comput. Methods Appl. Mech. Engrg. 198 (2009), no. 49–52, 3975–3988.10.1016/j.cma.2009.09.005Search in Google Scholar

[67] A. Prohl, Projection and Quasi-Compressibility Methods for Solving the Incompressible Navier–Stokes Equations, Adv. Numer. Math., B. G. Teubner, Stuttgart, 1997.10.1007/978-3-663-11171-9Search in Google Scholar

[68] J. Shen, On error estimates of the penalty method for unsteady Navier–Stokes equations, SIAM J. Numer. Anal. 32 (1995), no. 2, 386–403.10.1137/0732016Search in Google Scholar

[69] G. I. Taylor and A. E. Green, Mechanism of the production of small eddies from large ones, Proc. Roy. Soc. Lond. A 158 (1937), no. 895, 499–521.10.1098/rspa.1937.0036Search in Google Scholar

[70] R. Témam, Sur l’approximation de la solution des équations de Navier–Stokes par la méthode des pas fractionnaires. I, Arch. Rational Mech. Anal. 32 (1969), 135–153.10.1007/BF00247678Search in Google Scholar

[71] N. Vladimirova, B. Kuznetsov and N. N. Yanenko, Numerical calculation of the symmetrical flow of viscous incompressible liquid around a plate, Some Problems in Computational and Applied Mathematics, Nauka, Moscow (1966).Search in Google Scholar

[72] N. N. Yanenko, The Method of Fractional Steps, Springer, New York, 1971.10.1007/978-3-642-65108-3Search in Google Scholar

[73] L. Yang, S. Badia and R. Codina, A pseudo-compressible variational multiscale solver for turbulent incompressible flows, Comput. Mech. 58 (2016), no. 6, 1051–1069.10.1007/s00466-016-1332-9Search in Google Scholar

Received: 2021-05-02

Revised: 2021-08-03

Accepted: 2021-08-10

Published Online: 2021-08-27

Published in Print: 2022-01-01

© 2021 Walter de Gruyter GmbH, Berlin/Boston