Kinetic BGK numerical schemes for the approximation of incompressible Navier-Stokes equations are derived via classical discrete velocity vector BGK approximations, but applied to an inviscid compressible gas dynamics system with small Mach number parameter, according to the approach of Carfora and Natalini (2008). As the Mach number, the grid size and the timestep tend to zero, the low Mach number limit and the time-space convergence of the scheme are achieved simultaneously, and the numerical viscosity tends to the physical viscosity of the Navier-Stokes system. The method is analyzed and formulated as an explicit finite volume/difference flux vector splitting (FVS) scheme over a Cartesian mesh. It is close in spirit to lattice Boltzmann schemes, but it has several advantages. The first is that the scheme is expressed only in terms of momentum and mass compressible variables. It is therefore very easy to implement, and several types of boundary conditions are straightforward to apply. The second advantage is that the scheme satisfies a discrete entropy inequality, under a CFL condition of parabolic type and a subcharacteristic stability condition involving a cell Reynolds number that ensures that diffusion dominates advection at the level of the grid size. This ensures the robustness of the method, with explicit uniform bounds on the approximate solution. Moreover the scheme is proved to be second-order accurate in space if the parameters are well chosen, this is the case in particular for the Lax-Friedrichs scheme with Mach number proportional to the grid size. The scheme falls then into the class of artificial compressibility methods, the novelty being its exceptionally good theoretical properties. We show the efficiency of the method in terms of accuracy and robustness on a variety of classical two-dimensional benchmark tests. The method is finally applied in three dimensions to compute the permeability of a porous medium defined by a complex idealized Kelvin-like cell. Relations between our scheme and compressible low Mach number schemes are discussed.

Classification: 65M12, 65X68

Keywords: incompressible Navier-Stokes equations, vector BGK schemes, flux vector splitting, low Mach number limit, discrete entropy inequality, cell Reynolds number, lattice Boltzmann schemes

@article{SMAI-JCM_2018__4__1_0, author = {Fran\c cois Bouchut and Yann Jobic and Roberto Natalini and Ren\'e Occelli and Vincent Pavan}, title = {Second-order entropy satisfying BGK-FVS schemes for incompressible Navier-Stokes equations}, journal = {The SMAI journal of computational mathematics}, publisher = {Soci\'et\'e de Math\'ematiques Appliqu\'ees et Industrielles}, volume = {4}, year = {2018}, pages = {1-56}, doi = {10.5802/smai-jcm.28}, language = {en}, url={smai-jcm.centre-mersenne.org/item/SMAI-JCM_2018__4__1_0/} }

Bouchut, François; Jobic, Yann; Natalini, Roberto; Occelli, René; Pavan, Vincent. Second-order entropy satisfying BGK-FVS schemes for incompressible Navier-Stokes equations. The SMAI journal of computational mathematics, Volume 4 (2018) , pp. 1-56. doi : 10.5802/smai-jcm.28. https://smai-jcm.centre-mersenne.org/item/SMAI-JCM_2018__4__1_0/

[1] Low Mach number limit of the full Navier-Stokes equations, Arch. Rat. Mech. Anal., Volume 180 (2006), pp. 1-73

[2] A lattice Boltzmann model for compressible fluids, Phys. Rev. A, Volume 46 (1992), pp. 1967-1970

[3] Discrete kinetic schemes for multidimensional systems of conservation laws, SIAM J. Numer. Anal., Volume 37 (2000), pp. 1973-2004

[4] Explicit diffusive kinetic schemes for nonlinear degenerate parabolic systems, Math. Comput., Volume 73 (2004), pp. 63-94

[5] Experimental and theoretical investigation of backward-facing step flow, J. Fluid Mech., Volume 127 (1983), pp. 473-496

[6] Link-wise artificial compressibility method, J. Comput. Phys., Volume 231 (2012), pp. 5109-5143

[7] A multilayer Saint-Venant system with mass exchanges for shallow water flows. Derivation and numerical validation, ESAIM: Math. Modelling Numer. Anal., Volume 45 (2011), pp. 169-200

[8] A cellular automaton for Burgers’s equation, Complex Systems, Volume 1 (1987), pp. 17-30

[9] Fluid dynamic limits of kinetic equations. I. Formal derivations, J. Stat. Phys., Volume 63 (1991), pp. 323-344

[10] Fluid dynamic limits of kinetic equations-II Convergence proofs for the Boltzmann-equation, Comm. Pure Appl. Math., Volume 46 (1993), pp. 667-753

[11] The acoustic limit for the Boltzmann equation, Archive Rat. Mech. Anal., Volume 153 (2000), pp. 177-204

[12] Simulation of complex flows and multi-physics with the lattice-Boltzmann method (2008), 114, 76 pages (Ph. D. Thesis)

[13] Kinetic invariant domains and relaxation limit from a BGK model to isentropic gas dynamics, Asymptotic analysis, Volume 31 (2002), pp. 153-176

[14] Relaxation to isentropic gas dynamics for a BGK system with single kinetic entropy, Meth. and Appl. of Analysis, Volume 9 (2002), pp. 313-327

[15] Construction of BGK models with a family of kinetic entropies for a given system of conservation laws, J. Stat. Phys., Volume 95 (1999), pp. 113-170

[16] Entropy satisfying flux vector splittings and kinetic BGK models, Numer. Math., Volume 94 (2003), pp. 623-672

[17] Nonlinear stability of finite volume methods for hyperbolic conservation laws, and well-balanced schemes for sources, Frontiers in Mathematics, Birkhäuser Verlag, Basel, Boston, Berlin, 2004

[18] A reduced stability condition for nonlinear relaxation to conservation laws, J. Hyp. Diff. Eq., Volume 1 (2004), pp. 149-170

[19] Finite difference schemes with cross derivatives correctors for multidimensional parabolic systems, J. Hyp. Diff. Eq., Volume 3 (2006), pp. 27-52

[20] Diffusive BGK approximations for nonlinear multidimensional parabolic equations, Indiana Univ. Math. J., Volume 49 (2000), pp. 723-749

[21] A multiwave approximate Riemann solver for ideal MHD based on relaxation I - Theoretical framework, Numer. Math., Volume 108 (2007), pp. 7-42

[22] Mathematical tools for the study of the incompressible Navier-Stokes equations and related models, Applied Mathematical Sciences, Springer, New York, 2012

[23] Résolution d’équations d’évolution quasilinéaires en dimension N d’espace à l’aide d’équations linéaires en dimension N+1, J. Diff. Eq., Volume 50 (1983), pp. 375-390

[24] Averaged multivalued solutions for scalar conservation laws, SIAM J. Numer. Anal., Volume 21 (1984), pp. 1013-1037

[25] Low Mach number limit of viscous polytropic flows: formal asymptotics in the periodic case, Stud. Appl. Math., Volume 109 (2002), pp. 125-149

[26] Anelastic limits for Euler-type systems, Appl. Math. Res. Express., Volume 2010 (2010), pp. 119-141

[27] Accurate computations of the laminar flow past a square cylinder based on two different methods: lattice-Boltzmann and finite-volume, Intern. J. Heat Fluid Flow, Volume 21 (2000), pp. 186-196

[28] New efficient boundary conditions for incompressible Navier-Stokes equations: a well-posedness result, Math. Model. Numer. Anal., Volume 30 (1996), pp. 815-840

[29] A discrete kinetic approximation for the incompressible Navier-Stokes equations, ESAIM: Math. Modelling Numer. Anal., Volume 42 (2008), pp. 93-112

[30] An all-regime Lagrange-Projection like scheme for the gas dynamics equations on unstructured meshes, Comm. Comput. Phys., Volume 20 (2016), pp. 188-233

[31] Hyperbolic conservation laws with stiff relaxation terms and entropy, Comm. Pure Appl. Math., Volume 47 (1994), pp. 787-830

[32] Lattices for the lattice Boltzmann method, Phys. Rev. E, Volume 79 (2009), 046701 pages

[33] Hydraulic permeability of (un)bounded fibrous media using the lattice Boltzmann method, Phys. Rev. E, Volume 61 (2000), pp. 616-625

[34] An asymptotic-preserving all-speed scheme for the Euler and Navier-Stokes equations, J. Comput. Phys., Volume 231 (2012), pp. 5685-5704

[35] Les fontaines publiques de la ville de Dijon, Dalmont, Paris, 1856

[36] All speed scheme for the low Mach number limit of the isentropic Euler equations, Commun. Comput. Phys., Volume 10 (2011), pp. 1-31

[37] Analysis of Godunov type schemes applied to the compressible Euler system at low Mach number, J. Comput. Phys., Volume 229 (2010), pp. 978-1016

[38] Construction and analysis of lattice Boltzmann methods applied to a 1D convection-diffusion equation, Acta Applicandae Mathematicae, Volume 131 (2014), pp. 69-140

[39] On the Godunov scheme applied to the variable cross-section linear wave equation, Springer Proc. Math., Volume 4 (2011), pp. 313-321

[40] Kinetic theory based new upwind methods for inviscid compressible flows (1986) (AIAA 24th Aerospace Science Meeting, Jan 6-9 1986, Nevada, USA, AIAA paper 86-0275)

[41] On the Maxwellian distribution, symmetric form and entropy conservation for the Euler equations (1986) (NASA Langley Research Centre, Hampton, VA, 1986, NASA TP2613)

[42] Generalized lattice-Boltzmann equations, Rarefied Gas Dynamics: Theory and Simulations, vol. 159 of AIAA Progress in Astronautics and Astronautics, 1992, pp. 450-458

[43] Multiple-relaxation-time lattice Boltzmann models in three dimensions, Phil. Trans. Royal Soc. London series A-Math. Phys. Eng. Sci., Volume 360 (2002), pp. 437-451

[44] Décomposition de flux et discontinuité de contact stationnaire, C. R. Acad. Sci. Paris Sér. I Math., Volume 330 (2000), pp. 847-850

[45] Equivalent partial differential equations of a lattice Boltzmann scheme, Computers Math. Applic., Volume 55 (2008), pp. 1441-1449

[46] Convergence of convective-diffusive lattice Boltzmann methods, SIAM J. Numer. Anal., Volume 32 (1995), pp. 1327-1354

[47] Numerical solutions of 2-D steady incompressible flow over a backward-facing step, Part I: High Reynolds number solutions, Computers Fluids, Volume 37 (2008), pp. 633-655

[48] Positivity of flux vector splitting schemes, J. Comput. Phys., Volume 155 (1999), pp. 199-220

[49] Boltzmann approach to lattice gas simulations, Europhys. Letters, Volume 9 (1989), pp. 663-668

[50] Accuracy of the lattice Boltzmann method for small Knudsen number with finite Reynolds number, Phys. Fluids, Volume 9 (1997), 3535 pages

[51] Diffusive relaxation schemes for multiscale discrete-velocity kinetic equations, SIAM J. Numer. Anal., Volume 35 (1998), pp. 2405-2439

[52] Uniformly accurate diffusive relaxation schemes for multiscale transport equations, SIAM J. Numer. Anal., Volume 38 (2000), pp. 913-936

[53] Maximum entropy principle for lattice kinetic equations, Phys. Rev. Lett., Volume 81 (1998), pp. 6-9

[54] Compressible and incompressible fluids, Comm. Pure Appl. Math., Volume 35 (1982), pp. 629-651

[55] Micro-structural impact of different strut shapes and porosity on hydraulic properties of Kelvin-like metal foams, Transport in Porous Media, Volume 105 (2014), pp. 57-81

[56] Theory of the lattice Boltzmann method: dispersion, dissipation, isotropy, Galilean invariance, and stability, Physical Review E, Volume 61 (2000), pp. 6546-6562

[57] An all-speed Roe-type scheme and its asymptotic analysis of low Mach number behaviour, J. Comput. Phys., Volume 227 (2008), pp. 5144-5159

[58] Incompressible limit for a viscous compressible fluid, J. Math. Pures Appl., Volume 77 (1998), pp. 585-627

[59] Higher order accurate kinetic flux vector splitting method for Euler equations, Nonlinear Hyperbolic Equations – Theory, Computation Methods, and Applications (Notes on Numerical Fluid Mechanics Series 24), 1989, pp. 384-392

[60] The respiratory system in equations, Modeling, Simulation and Applications 7, Springer, Milan, 2013

[61] Use of the Boltzmann-equation to simulate lattice-gas automata, Phys. Rev. Letters, Volume 61 (1988), pp. 2332-2335

[62] Consistent initial conditions for lattice Boltzmann simulations, Computers Fluids, Volume 35 (2006), pp. 855-862

[63] Benchmark solution for a three-dimensional mixed-convection flow, part 1: reference solutions, Numer. Heat Transfer B-Fundamentals, Volume 60 (2011), pp. 325-345

[64] A weakly asymptotic preserving low Mach number scheme for the Euler equations of gas dynamics, SIAM J. Sci. Comput., Volume 36 (2014), p. B989-B1024

[65] Artificial compressibility method revisited: Asymptotic numerical method for incompressible Navier-Stokes equations, J. Comput. Phys., Volume 229 (2010), pp. 1698-1723

[66] Artificial compressibility method and lattice Boltzmann method: similarities and differences, Computers Math. Applic., Volume 61 (2011), pp. 3461-3474

[67] On simulation of outflow boundary conditions in finite difference calculations for incompressible fluid, Int. J. Numer. Meth. Fluids, Volume 33 (2000), pp. 499-534

[68] Comparison of low Mach number models for natural convection problems, Heat and Mass Transfer, Volume 36 (2000), pp. 567-573

[69] A new outflow boundary-condition, Int. J. Numer. Meth. Fluids, Volume 14 (1992), pp. 587-608

[70] Boltzmann type schemes for gas dynamics and the entropy property, SIAM J. Numer. Anal., Volume 27 (1990), pp. 1405-1421

[71] A kinetic-equation with kinetic entropy functions for scalar conservation laws, Comm. Math. Phys., Volume 136 (1991), pp. 501-517

[72] Open boundary conditions for the velocity-correction scheme of the Navier-Stokes equations, Computers Fluids, Volume 70 (2012), pp. 29-43

[73] Direct simulation methods for compressible inviscid ideal-gas flow, J. Comput. Phys., Volume 34 (1980), pp. 231-244

[74] One-dimensional compressible gas dynamics calculations using the Boltzmann equations, J. Comput. Phys., Volume 42 (1981), pp. 108-123

[75] Laboratory determination of the full permeability tensor, J. Geophysical Res., Volume 106 (2001), pp. 26443-26452

[76] Hydrodynamic limits of the Boltzmann equation, Lecture Notes in Mathematics 1971, Springer, Berlin, 2009

[77] Fast singular limits of hyperbolic PDEs, J. Diff. Eq., Volume 114 (1994), pp. 476-512

[78] Flux-difference splitting-based upwind compact schemes for the incompressible Navier-Stokes equations, Int. J. Numer. Meth. Fluids, Volume 61 (2009), pp. 552-568

[79] Upwind compact finite difference scheme for time-accurate solution of the incompressible Navier-Stokes equations, Applied Math. Comput., Volume 215 (2010), pp. 3201-3213

[80] Low-Reynolds-number flow around a square cylinder at incidence: study of blockage, onset of vortex shedding and outlet boundary condition, Intern. J. Numer. Meth. Fluids, Volume 26 (1998), pp. 39-56

[81] The lattice Boltzmann equation for fluid dynamics and beyond, Numerical Mathematics and Scientific Computation, Oxford Science Publications, the Clarendon Press, Oxford University Press, New York, 2001

[82] 3-dimensional flows in complex geometries with the lattice Boltzmann method, Europhys. Letters, Volume 10 (1989), pp. 433-438

[83] Role of the H theorem in lattice Boltzmann hydrodynamic simulations, Rev. Mod. Phys., Volume 74 (2002), pp. 1203-1220

[84] Assessment of the two relaxation time lattice-Boltzmann scheme to simulate Stokes flow in porous media, Water Resources Res., Volume 48 (2012), W04526 pages

[85] Preconditioned methods for solving the incompressible and low speed compressible equations, J. Comput. Phys., Volume 72 (1987), pp. 277-298

[86] Nonexistence of H theorem for some lattice Boltzmann models, J. Stat. Phys., Volume 121 (2005), pp. 91-103