Lorenzo Botti

Lorenzo Botti, Daniele A. Di Pietro, Jérôme Droniou

In this work, we develop and analyse a novel Hybrid High-Order discretisation of the Brinkman problem. The method hinges on hybrid discrete velocity unknowns at faces and elements and on discontinuous pressures. Based on the discrete unknowns, we reconstruct inside each element a Stokes velocity one degree higher than face unknowns, and a Darcy velocity in the Raviart-Thomas-Nédélec space. These reconstructed velocities are respectively used to formulate the discrete versions of the Stokes and Darcy terms in the momentum equation, along with suitably designed penalty contributions. The proposed construction is tailored to yield optimal error estimates that are robust throughout the entire spectrum of local (Stokes- or Darcy-dominated) regimes, as identified by a dimensionless number which can be interpreted as a friction coefficient. The singular limit corresponding to the Darcy equation is also fully supported by the method. Numerical examples corroborate the theoretical results. This paper also contains two contributions whose interest goes beyond the specific method and application treated in this work: an investigation of the dependence of the constant in the second Korn inequality on star-shaped domains and its application to the study of the approximation properties of the strain projector in general Sobolev seminorms.

Matteo Franciolini, Lorenzo Botti, Alessandro Colombo et al.

In recent years several research efforts focused on the development of high-order discontinuous Galerkin (dG) methods for scale resolving simulations of turbulent flows. Nevertheless, in the context of incompressible flow computations, the computational expense of solving large scale equation systems characterized by indefinite Jacobian matrices has often prevented from dealing with industrially-relevant computations. In this work we seek to improve the efficiency of Rosenbrock-type linearly-implicit Runge-Kutta methods by devising robust, scalable and memory-lean solution strategies. In particular, we introduce memory saving p-multigrid preconditioners coupling matrix-free and matrix-based Krylov iterative smoothers. The p-multigrid preconditioner relies on cheap block-diagonal smoother's preconditioners on the fine space to reduce assembly costs and memory allocation, and ensures an adequate resolution of the coarsest space of the multigrid iteration using Additive Schwarz precondioned smoothers to obtain satisfactory convergence rates and optimal parallel efficiency of the method. Extensive numerical validation is performed. The Rosenbrock formulation is applied to test cases of growing complexity: the laminar unsteady flow around a two-dimensional cylinder at Re=200 and around a sphere at Re=300, the transitional flow problem of the ERCOFTAC T3L test case suite with different levels of free-stream turbulence. As proof of concept, the numerical solution of the Boeing Rudimentary Landing Gear test case at Re=10^6 is reported. A good agreement of the solutions with experimental data is documented, as well as strong memory savings and execution time gains with respect to state-of-the art solution strategies.

Lorenzo Botti, Daniele Di Pietro, Jérôme Droniou

In this work we propose a novel Hybrid High-Order method for the incompressible Navier--Stokes equations based on a formulation of the convective term including Temam's device for stability. The proposed method has several advantageous features: it supports arbitrary approximation orders on general meshes including polyhedral elements and non-matching interfaces; it is inf-sup stable; it is locally conservative; it supports both the weak and strong enforcement of velocity boundary conditions; it is amenable to efficient computer implementations where a large subset of the unknowns is eliminated by solving local problems inside each element. Particular care is devoted to the design of the convective trilinear form, which mimicks at the discrete level the non-dissipation property of the continuous one. The possibility to add a convective stabilisation term is also contemplated, and a formulation covering various classical options is discussed. The proposed method is theoretically analysed, and an energy error estimate in $h^{k+1}$ (with $h$ denoting the meshsize) is proved under the usual data smallness assumption. A thorough numerical validation on two and three-dimensional test cases is provided both to confirm the theoretical convergence rates and to assess the method in more physical configurations (including, in particular, the well-known two- and three-dimensional lid-driven cavity problems).