With cedram.org   version française
Search for an article
Table of contents for this volume | Next article
François Bouchut; Yann Jobic; Roberto Natalini; René Occelli; Vincent Pavan
Second-order entropy satisfying BGK-FVS schemes for incompressible Navier-Stokes equations
SMAI-Journal of computational mathematics, 4 (2018), p. 1-56, doi: 10.5802/smai-jcm.28
Article PDF
Class. Math.: 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


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.


[1] T. Alazard, “Low Mach number limit of the full Navier-Stokes equations”, Arch. Rat. Mech. Anal. 180 (2006), p. 1-73
[2] F. Alexander, H. Chen, S. Chen & G. Doolen, “A lattice Boltzmann model for compressible fluids”, Phys. Rev. A 46 (1992), p. 1967-1970
[3] D. Aregba-Driollet & R. Natalini, “Discrete kinetic schemes for multidimensional systems of conservation laws”, SIAM J. Numer. Anal. 37 (2000), p. 1973-2004
[4] D. Aregba-Driollet, R. Natalini & S. Tang, “Explicit diffusive kinetic schemes for nonlinear degenerate parabolic systems”, Math. Comput. 73 (2004), p. 63-94
[5] B.F. Armaly, F. Durst, J.C.F. Pereira & B. Schonung, “Experimental and theoretical investigation of backward-facing step flow”, J. Fluid Mech. 127 (1983), p. 473-496
[6] P. Asinari, T. Ohwada, E. Chiavazzo & A.F. Di Rienzo, “Link-wise artificial compressibility method”, J. Comput. Phys. 231 (2012), p. 5109-5143
[7] E. Audusse, M.-O. Bristeau, B. Perthame & J. Sainte-Marie, “A multilayer Saint-Venant system with mass exchanges for shallow water flows. Derivation and numerical validation”, ESAIM: Math. Modelling Numer. Anal. 45 (2011), p. 169-200
[8] C. Levermore B. Boghosian, “A cellular automaton for Burgers’s equation”, Complex Systems 1 (1987), p. 17-30
[9] C. Bardos, F. Golse & C.D. Levermore, “Fluid dynamic limits of kinetic equations. I. Formal derivations”, J. Stat. Phys. 63 (1991), p. 323-344
[10] C. Bardos, F. Golse & C.D. Levermore, “Fluid dynamic limits of kinetic equations-II Convergence proofs for the Boltzmann-equation”, Comm. Pure Appl. Math. 46 (1993), p. 667-753
[11] C. Bardos, F. Golse & C.D. Levermore, “The acoustic limit for the Boltzmann equation”, Archive Rat. Mech. Anal. 153 (2000), p. 177-204
[12] J. Bernsdorf, Simulation of complex flows and multi-physics with the lattice-Boltzmann method, Ph. D. Thesis, Univ. Amsterdam, 2008
[13] F. Berthelin & F. Bouchut, “Kinetic invariant domains and relaxation limit from a BGK model to isentropic gas dynamics”, Asymptotic analysis 31 (2002), p. 153-176
[14] F. Berthelin & F. Bouchut, “Relaxation to isentropic gas dynamics for a BGK system with single kinetic entropy”, Meth. and Appl. of Analysis 9 (2002), p. 313-327
[15] F. Bouchut, “Construction of BGK models with a family of kinetic entropies for a given system of conservation laws”, J. Stat. Phys. 95 (1999), p. 113-170
[16] F. Bouchut, “Entropy satisfying flux vector splittings and kinetic BGK models”, Numer. Math. 94 (2003), p. 623-672
[17] F. Bouchut, 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] F. Bouchut, “A reduced stability condition for nonlinear relaxation to conservation laws”, J. Hyp. Diff. Eq. 1 (2004), p. 149-170
[19] F. Bouchut & H. Frid, “Finite difference schemes with cross derivatives correctors for multidimensional parabolic systems”, J. Hyp. Diff. Eq. 3 (2006), p. 27-52
[20] F. Bouchut, F.R. Guarguaglini & R. Natalini, “Diffusive BGK approximations for nonlinear multidimensional parabolic equations”, Indiana Univ. Math. J. 49 (2000), p. 723-749
[21] F. Bouchut, C. Klingenberg & K. Waagan, “A multiwave approximate Riemann solver for ideal MHD based on relaxation I - Theoretical framework”, Numer. Math. 108 (2007), p. 7-42
[22] F. Boyer & P. Fabrie, Mathematical tools for the study of the incompressible Navier-Stokes equations and related models, Applied Mathematical Sciences, Springer, New York, 2012
[23] Y. Brenier, “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. 50 (1983), p. 375-390
[24] Y. Brenier, “Averaged multivalued solutions for scalar conservation laws”, SIAM J. Numer. Anal. 21 (1984), p. 1013-1037
[25] D. Bresch, B. Desjardins, E. Grenier & C.-K. Lin, “Low Mach number limit of viscous polytropic flows: formal asymptotics in the periodic case”, Stud. Appl. Math. 109 (2002), p. 125-149
[26] D. Bresch & G. Métivier, “Anelastic limits for Euler-type systems”, Appl. Math. Res. Express. 2010 (2010), p. 119-141
[27] M. Breuer, J. Bernsdorf, T. Zeiser & F. Durst, “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 21 (2000), p. 186-196
[28] C.-H. Bruneau & P. Fabrie, “New efficient boundary conditions for incompressible Navier-Stokes equations: a well-posedness result”, Math. Model. Numer. Anal. 30 (1996), p. 815-840
[29] M. Carfora & R. Natalini, “A discrete kinetic approximation for the incompressible Navier-Stokes equations”, ESAIM: Math. Modelling Numer. Anal. 42 (2008), p. 93-112
[30] C. Chalons, M. Girardin & S. Kokh, “An all-regime Lagrange-Projection like scheme for the gas dynamics equations on unstructured meshes”, Comm. Comput. Phys. 20 (2016), p. 188-233
[31] G.Q. Chen, C.D. Levermore & T.-P. Liu, “Hyperbolic conservation laws with stiff relaxation terms and entropy”, Comm. Pure Appl. Math. 47 (1994), p. 787-830
[32] S. Chikatamarla & I. Karlin, “Lattices for the lattice Boltzmann method”, Phys. Rev. E 79 (2009)
[33] D.S. Clague, B.D. Kandhai, R. Zhang & P.M.A. Sloot, “Hydraulic permeability of (un)bounded fibrous media using the lattice Boltzmann method”, Phys. Rev. E 61 (2000), p. 616-625
[34] F. Cordier, P. Degond & A. Kumbaro, “An asymptotic-preserving all-speed scheme for the Euler and Navier-Stokes equations”, J. Comput. Phys. 231 (2012), p. 5685-5704
[35] H. Darcy, Les fontaines publiques de la ville de Dijon, Dalmont, Paris, 1856
[36] P. Degond & M. Tang, “All speed scheme for the low Mach number limit of the isentropic Euler equations”, Commun. Comput. Phys. 10 (2011), p. 1-31
[37] S. Dellacherie, “Analysis of Godunov type schemes applied to the compressible Euler system at low Mach number”, J. Comput. Phys. 229 (2010), p. 978-1016
[38] S. Dellacherie, “Construction and analysis of lattice Boltzmann methods applied to a 1D convection-diffusion equation”, Acta Applicandae Mathematicae 131 (2014), p. 69-140
[39] S. Dellacherie & P. Omnes, “On the Godunov scheme applied to the variable cross-section linear wave equation”, Springer Proc. Math. 4 (2011), p. 313-321
[40] S.M. Deshpande, “Kinetic theory based new upwind methods for inviscid compressible flows”, AIAA 24th Aerospace Science Meeting, Jan 6-9 1986, Nevada, USA, AIAA paper 86-0275, 1986
[41] S.M. Deshpande, “On the Maxwellian distribution, symmetric form and entropy conservation for the Euler equations”, NASA Langley Research Centre, Hampton, VA, 1986, NASA TP2613, 1986
[42] D. d’Humieres, Generalized lattice-Boltzmann equations, Rarefied Gas Dynamics: Theory and Simulations, vol. 159 of AIAA Progress in Astronautics and Astronautics, 1992, p. 450-458
[43] D. d’Humieres, I. Ginzburg, M. Krafczyk, P. Lallemand & L.S. Luo, “Multiple-relaxation-time lattice Boltzmann models in three dimensions”, Phil. Trans. Royal Soc. London series A-Math. Phys. Eng. Sci. 360 (2002), p. 437-451
[44] F. Dubois, “Décomposition de flux et discontinuité de contact stationnaire”, C. R. Acad. Sci. Paris Sér. I Math. 330 (2000), p. 847-850
[45] F. Dubois, “Equivalent partial differential equations of a lattice Boltzmann scheme”, Computers Math. Applic. 55 (2008), p. 1441-1449
[46] B. Elton, C. Levermore & G. Rodrigue, “Convergence of convective-diffusive lattice Boltzmann methods”, SIAM J. Numer. Anal. 32 (1995), p. 1327-1354
[47] E. Erturk, “Numerical solutions of 2-D steady incompressible flow over a backward-facing step, Part I: High Reynolds number solutions”, Computers Fluids 37 (2008), p. 633-655
[48] J. Gressier, P. Villedieu & J.-M. Moschetta, “Positivity of flux vector splitting schemes”, J. Comput. Phys. 155 (1999), p. 199-220
[49] F.J. Higuera & J. Jimenez, “Boltzmann approach to lattice gas simulations”, Europhys. Letters 9 (1989), p. 663-668
[50] T. Inamuro, M. Yoshino & F. Ogino, “Accuracy of the lattice Boltzmann method for small Knudsen number with finite Reynolds number”, Phys. Fluids 9 (1997)
[51] S. Jin, L. Pareschi & G. Toscani, “Diffusive relaxation schemes for multiscale discrete-velocity kinetic equations”, SIAM J. Numer. Anal. 35 (1998), p. 2405-2439
[52] S. Jin, L. Pareschi & G. Toscani, “Uniformly accurate diffusive relaxation schemes for multiscale transport equations”, SIAM J. Numer. Anal. 38 (2000), p. 913-936
[53] I.V. Karlin, A.N. Gorban, S. Succi & V. Boffi, “Maximum entropy principle for lattice kinetic equations”, Phys. Rev. Lett. 81 (1998), p. 6-9
[54] S. Klainerman & A. Majda, “Compressible and incompressible fluids”, Comm. Pure Appl. Math. 35 (1982), p. 629-651
[55] P. Kumar & F. Topin, “Micro-structural impact of different strut shapes and porosity on hydraulic properties of Kelvin-like metal foams”, Transport in Porous Media 105 (2014), p. 57-81
[56] P. Lallemand & L.-S. Luo, “Theory of the lattice Boltzmann method: dispersion, dissipation, isotropy, Galilean invariance, and stability”, Physical Review E 61 (2000), p. 6546-6562
[57] X.-S. Li & C.-W. Gu, “An all-speed Roe-type scheme and its asymptotic analysis of low Mach number behaviour”, J. Comput. Phys. 227 (2008), p. 5144-5159
[58] P.-L. Lions & N. Masmoudi, “Incompressible limit for a viscous compressible fluid”, J. Math. Pures Appl. 77 (1998), p. 585-627
[59] J.C. Mandal & S.M. Deshpande, 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, p. 384-392
[60] B. Maury, The respiratory system in equations, Modeling, Simulation and Applications 7, Springer, Milan, 2013
[61] G.R. McNamara & G. Zanetti, “Use of the Boltzmann-equation to simulate lattice-gas automata”, Phys. Rev. Letters 61 (1988), p. 2332-2335
[62] R. Mei, L.S. Luo, P. Lallemand & D. d’Humieres, “Consistent initial conditions for lattice Boltzmann simulations”, Computers Fluids 35 (2006), p. 855-862
[63] X. Nicolas, M. Medale, S. Glockner & S. Gounand, “Benchmark solution for a three-dimensional mixed-convection flow, part 1: reference solutions”, Numer. Heat Transfer B-Fundamentals 60 (2011), p. 325-345
[64] S. Noelle, G. Bispen, K.R. Arun, M. Lukacova-Medvidova & C.-D. Munz, “A weakly asymptotic preserving low Mach number scheme for the Euler equations of gas dynamics”, SIAM J. Sci. Comput. 36 (2014), p. B989-B1024
[65] T. Ohwada & P. Asinari, “Artificial compressibility method revisited: Asymptotic numerical method for incompressible Navier-Stokes equations”, J. Comput. Phys. 229 (2010), p. 1698-1723
[66] T. Ohwada, P. Asinari & D. Yabusaki, “Artificial compressibility method and lattice Boltzmann method: similarities and differences”, Computers Math. Applic. 61 (2011), p. 3461-3474
[67] M.A. Ol’shanskii & V.M. Staroverov, “On simulation of outflow boundary conditions in finite difference calculations for incompressible fluid”, Int. J. Numer. Meth. Fluids 33 (2000), p. 499-534
[68] H. Paillére, C. Viozat, A. Kumbaro & I. Toumi, “Comparison of low Mach number models for natural convection problems”, Heat and Mass Transfer 36 (2000), p. 567-573
[69] T.C. Papanastasiou, N. Malamataris & K. Ellwood, “A new outflow boundary-condition”, Int. J. Numer. Meth. Fluids 14 (1992), p. 587-608
[70] B. Perthame, “Boltzmann type schemes for gas dynamics and the entropy property”, SIAM J. Numer. Anal. 27 (1990), p. 1405-1421
[71] B. Perthame & E. Tadmor, “A kinetic-equation with kinetic entropy functions for scalar conservation laws”, Comm. Math. Phys. 136 (1991), p. 501-517
[72] A. Poux, S. Glockner, E. Ahusborde & M. Azaiez, “Open boundary conditions for the velocity-correction scheme of the Navier-Stokes equations”, Computers Fluids 70 (2012), p. 29-43
[73] D.I. Pullin, “Direct simulation methods for compressible inviscid ideal-gas flow”, J. Comput. Phys. 34 (1980), p. 231-244
[74] R.D. Reitz, “One-dimensional compressible gas dynamics calculations using the Boltzmann equations”, J. Comput. Phys. 42 (1981), p. 108-123
[75] P. Renard, A. Genty & F. Stauffer, “Laboratory determination of the full permeability tensor”, J. Geophysical Res. 106 (2001), p. 26443-26452
[76] L. Saint-Raymond, Hydrodynamic limits of the Boltzmann equation, Lecture Notes in Mathematics 1971, Springer, Berlin, 2009
[77] S. Schochet, “Fast singular limits of hyperbolic PDEs”, J. Diff. Eq. 114 (1994), p. 476-512
[78] A. Shah & L. Yuan, “Flux-difference splitting-based upwind compact schemes for the incompressible Navier-Stokes equations”, Int. J. Numer. Meth. Fluids 61 (2009), p. 552-568
[79] A. Shah, L. Yuan & A. Khan, “Upwind compact finite difference scheme for time-accurate solution of the incompressible Navier-Stokes equations”, Applied Math. Comput. 215 (2010), p. 3201-3213
[80] A. Sohankar, C. Norberg & L. Davidson, “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 26 (1998), p. 39-56
[81] S. Succi, 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] S. Succi, E. Foti & F. Higuera, “3-dimensional flows in complex geometries with the lattice Boltzmann method”, Europhys. Letters 10 (1989), p. 433-438
[83] S. Succi, I.V. Karlin & H. Chen, “Role of the H theorem in lattice Boltzmann hydrodynamic simulations”, Rev. Mod. Phys. 74 (2002), p. 1203-1220
[84] L. Talon, D. Bauer, N. Gland, S. Youssef, H. Auradou & I. Ginzburg, “Assessment of the two relaxation time lattice-Boltzmann scheme to simulate Stokes flow in porous media”, Water Resources Res. 48 (2012)
[85] E. Turkel, “Preconditioned methods for solving the incompressible and low speed compressible equations”, J. Comput. Phys. 72 (1987), p. 277-298
[86] W.A. Yong & L.S. Luo, “Nonexistence of H theorem for some lattice Boltzmann models”, J. Stat. Phys. 121 (2005), p. 91-103