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.

DOI: https://doi.org/10.5802/smai-jcm.28

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{c}ois 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}, pages = {1--56}, publisher = {Soci\'et\'e de Math\'ematiques Appliqu\'ees et Industrielles}, volume = {4}, year = {2018}, doi = {10.5802/smai-jcm.28}, mrnumber = {3774647}, zbl = {1416.76243}, language = {en}, url = {https://smai-jcm.centre-mersenne.org/articles/10.5802/smai-jcm.28/} }

François Bouchut; Yann Jobic; Roberto Natalini; René Occelli; Vincent Pavan. 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/articles/10.5802/smai-jcm.28/

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

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

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

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

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

[6] Link-wise artificial compressibility method, J. Comput. Phys., Volume 231 (2012), pp. 5109-5143 | Article | MR 2929935 | Zbl 1247.76063

[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 | Article | Numdam | MR 2781135 | Zbl 1290.35194

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

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

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

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

[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 | MR 1938602 | Zbl 1032.76064

[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 | MR 1957492 | Zbl 1173.82351

[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 | Article | MR 1705583 | Zbl 0957.82028

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

[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 | Zbl 1086.65091

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

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

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

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

[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 | Article | Zbl 0549.35055

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

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

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

[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 | Article

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

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

[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 | Article | MR 3514275 | Zbl 1373.76120

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

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

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

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

[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 | Article | MR 2775032 | Zbl 1364.76129

[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 | Article | MR 2576236 | Zbl 1329.76228

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

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

[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) | Article

[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 | Article | MR 1902782 | Zbl 1001.76081

[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 | Article | Zbl 0983.76054

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

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

[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 | Article | Zbl 1237.76102

[48] Positivity of flux vector splitting schemes, J. Comput. Phys., Volume 155 (1999), pp. 199-220 | Article | MR 1716481 | Zbl 0953.76064

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

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

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

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

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

[54] Compressible and incompressible fluids, Comm. Pure Appl. Math., Volume 35 (1982), pp. 629-651 | Article | MR 668409 | Zbl 0478.76091

[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 | Article

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

[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 | MR 2414850 | Zbl 1388.76207

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

[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 | Article

[60] The respiratory system in equations, Modeling, Simulation and Applications 7, Springer, Milan, 2013 | MR 3076523 | Zbl 1312.92005

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

[62] Consistent initial conditions for lattice Boltzmann simulations, Computers Fluids, Volume 35 (2006), pp. 855-862 | Article | MR 2239882 | Zbl 1177.76319

[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 | Article

[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 | Article | MR 3293458 | Zbl 1321.76053

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

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

[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 | Article

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

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

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

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

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

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

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

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

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

[77] Fast singular limits of hyperbolic PDEs, J. Diff. Eq., Volume 114 (1994), pp. 476-512 | Article | MR 1303036 | Zbl 0838.35071

[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 | Article | MR 2572429 | Zbl 1172.76035

[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 | MR 2576808 | Zbl 1352.76084

[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 | Article | Zbl 0910.76067

[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 | Zbl 0990.76001

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

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

[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 | Article

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

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