A discontinuous skeletal method for the viscosity-dependent Stokes problem

Abstract

International audienceWe devise and analyze arbitrary-order nonconforming methods for the discretization of the viscosity-dependent Stokes equations on simplicial meshes. We keep track explicitly of the viscosity and aim at pressure-robust schemes that can deal with the practically relevant case of body forces with large curl-free part in a way that the discrete velocity error is not spoiled by large pressures. The method is inspired from the recent Hybrid High-Order (HHO) methods for linear elasticity. After elimination of the auxiliary variables by static condensation, the linear system to be solved involves only discrete face-based velocities, which are polynomials of degree $k \ge 0$, and cell-wise constant pressures. Our main result is a pressure-independent energy-error estimate on the velocity of order $(k+1)$. The main ingredient to achieve pressure-independence is the use of a divergence-preserving velocity reconstruction operator in the discretization of the body forces. We also prove an $L^2$-pressure estimate of order $(k+1)$ and an $L^2$-velocity estimate of order $(k+2)$, the latter under elliptic regularity. The local mass and momentum conservation properties of the discretization are also established. Finally, two-and three-dimensional numerical results are presented to support the analysis

On high-order pressure-robust space discretisations, their advantages for incompressible high Reynolds number generalised Beltrami flows and beyond

“…Fortunately, pressure-robust space discretisations behave in a robust manner when confronted with strong pressure gradients, and many different ways to construct such schemes have been found recently. To name only a few, inf-sup stable H 1 -conforming and divergence-free mixed methods [65], inf-sup stable H(div)-conforming DG methods [19,45] and inf-sup stable H 1 -conforming and nonconforming finite element methods (FEM), finite volume (FVM) methods, and Hybrid High Order methods (HHO) with appropriately modified velocity test functions [46,47,23,42,52,48] are pressure-robust. Moreover, also in the context of isogeometric analysis various pressure-robust discretisations have been developed [15,29,30].…”

A Hybrid High‐Order method for finite elastoplastic deformations within a logarithmic strain framework

“…The discrete problem (35) expresses the principle of virtual work at the global level, and adapting the ideas introduced in the work of Cockburn et al 39 (see also the works of Abbas et al 1 and Botti et al 33 ), it is possible to infer a local principle of virtual work in terms of face-based discrete tractions that comply with the law of action and reaction.…”

“…To evaluate the influence of the stabilization parameter 0 , we compare the total number of Newton's iterations needed to solve the nonlinear problem (35) versus the magnitude of the stabilization parameter 0 . Newton's iterations are stopped under the relative criterion…”

Hybrid Discretization Methods with Adaptive Yield Surface Detection for Bingham Pipe Flows

“…For earlier work on agglomeration-based coarsening, we refer, e.g., to [3,4] and references therein. For simplicity, we focus on Bingham pipe flows, but we mention that incompressible flows can be approximated by hybrid discretization methods, as shown, e.g., in [14,37] for HDG methods and in [19] for HHO methods.…”