Hybrid High-Order formulations with turbulence modelling capabilities for incompressible flow problems

We propose a Hybrid High-Order (HHO) formulation of the incompressible Navier--Stokes equations, that is well suited to be employed for the simulation of turbulent flows. The spatial discretization relies on hybrid velocity and pressure spaces and the temporal discretization is based on Explicit Singly Diagonal Implicit Runge-Kutta (ESDIRK) methods. The formulation possesses some attractive features that can be fruitfully exploited when high-fidelity computations are required, namely: pressure-robustness, conservation of mass enforced cell-by-cell up to machine precision, robustness in the inviscid limit, implicit high-order accurate time stepping with local time step adaptation, reduced memory footprint thanks to static condensation of both velocity and pressure, possibility to exploit inherited $p$-multilevel solution strategies to improve performance of iterative solvers. After demonstrating the relevant properties of the scheme in practice, performing challenging 2D and 3D test cases, we consider the simulation of the Taylor--Green Vortex flow problem at Reynolds 1600.