HHO Methods for the Incompressible Navier-Stokes and the Incompressible Euler Equations

• Published in: Journal of scientific computing 2022-07, Vol.92 (1), p.28, Article 28
• Main Authors: Botti, Lorenzo, Massa, Francesco Carlo
• Format: Article
• Language: English
• Subjects: Algorithms; Approximation; Computational Mathematics and Numerical Analysis; Condensation; Divergence; Energy; Error analysis; Euler-Lagrange equation; Fluid flow; Kinetic energy; Mathematical and Computational Engineering; Mathematical and Computational Physics; Mathematics; Mathematics and Statistics; Methods; Navier-Stokes equations; Polynomials; Reynolds number; Riemann solver; Robustness (mathematics); Theoretical; Velocity; Velocity coupling; Velocity distribution
• ISSN: 0885-7474; 1573-7691
• DOI: 10.1007/s10915-022-01864-1

We propose two Hybrid High-Order (HHO) methods for the incompressible Navier-Stokes equations and investigate their robustness with respect to the Reynolds number. While both methods rely on a HHO formulation of the viscous term, the pressure-velocity coupling is fundamentally different, up to the point that the two approaches can be considered antithetical. The first method is kinetic energy preserving, meaning that the skew-symmetric discretization of the convective term is guaranteed not to alter the kinetic energy balance. The approximated velocity fields exactly satisfy the divergence free constraint and continuity of the normal component of the velocity is weakly enforced on the mesh skeleton, leading to H-div conformity. The second scheme relies on Godunov fluxes for pressure-velocity coupling: a Harten, Lax and van Leer approximated Riemann Solver designed for cell centered formulations is adapted to hybrid face centered formulations. The resulting numerical scheme is robust up to the inviscid limit, meaning that it can be applied for seeking approximate solutions of the incompressible Euler equations. The schemes are numerically validated performing steady and unsteady two dimensional test cases and evaluating the convergence rates on h -refined mesh sequences. In addition to standard benchmark flow problems, specifically conceived test cases are conducted for studying the error behaviour when approaching the inviscid limit.