11,984 research outputs found
On high-order pressure-robust space discretisations, their advantages for incompressible high Reynolds number generalised Beltrami flows and beyond
An improved understanding of the divergence-free constraint for the
incompressible Navier--Stokes equations leads to the observation that a
semi-norm and corresponding equivalence classes of forces are fundamental for
their nonlinear dynamics. The recent concept of {\em pressure-robustness}
allows to distinguish between space discretisations that discretise these
equivalence classes appropriately or not. This contribution compares the
accuracy of pressure-robust and non-pressure-robust space discretisations for
transient high Reynolds number flows, starting from the observation that in
generalised Beltrami flows the nonlinear convection term is balanced by a
strong pressure gradient. Then, pressure-robust methods are shown to outperform
comparable non-pressure-robust space discretisations. Indeed, pressure-robust
methods of formal order are comparably accurate than non-pressure-robust
methods of formal order on coarse meshes. Investigating the material
derivative of incompressible Euler flows, it is conjectured that strong
pressure gradients are typical for non-trivial high Reynolds number flows.
Connections to vortex-dominated flows are established. Thus,
pressure-robustness appears to be a prerequisite for accurate incompressible
flow solvers at high Reynolds numbers. The arguments are supported by numerical
analysis and numerical experiments.Comment: 43 pages, 18 figures, 2 table
Stabilized mixed finite element methods for linear elasticity on simplicial grids in
In this paper, we design two classes of stabilized mixed finite element
methods for linear elasticity on simplicial grids. In the first class of
elements, we use - and
- to approximate the stress
and displacement spaces, respectively, for , and employ a
stabilization technique in terms of the jump of the discrete displacement over
the faces of the triangulation under consideration; in the second class of
elements, we use - to
approximate the displacement space for , and adopt the
stabilization technique suggested by Brezzi, Fortin, and Marini. We establish
the discrete inf-sup conditions, and consequently present the a priori error
analysis for them. The main ingredient for the analysis is two special
interpolation operators, which can be constructed using a crucial
bubble function space of polynomials on each
element. The feature of these methods is the low number of global degrees of
freedom in the lowest order case. We present some numerical results to
demonstrate the theoretical estimates.Comment: 16 pages, 1 figur
Implicit High-Order Flux Reconstruction Solver for High-Speed Compressible Flows
The present paper addresses the development and implementation of the first
high-order Flux Reconstruction (FR) solver for high-speed flows within the
open-source COOLFluiD (Computational Object-Oriented Libraries for Fluid
Dynamics) platform. The resulting solver is fully implicit and able to simulate
compressible flow problems governed by either the Euler or the Navier-Stokes
equations in two and three dimensions. Furthermore, it can run in parallel on
multiple CPU-cores and is designed to handle unstructured grids consisting of
both straight and curved edged quadrilateral or hexahedral elements. While most
of the implementation relies on state-of-the-art FR algorithms, an improved and
more case-independent shock capturing scheme has been developed in order to
tackle the first viscous hypersonic simulations using the FR method. Extensive
verification of the FR solver has been performed through the use of
reproducible benchmark test cases with flow speeds ranging from subsonic to
hypersonic, up to Mach 17.6. The obtained results have been favorably compared
to those available in literature. Furthermore, so-called super-accuracy is
retrieved for certain cases when solving the Euler equations. The strengths of
the FR solver in terms of computational accuracy per degree of freedom are also
illustrated. Finally, the influence of the characterizing parameters of the FR
method as well as the the influence of the novel shock capturing scheme on the
accuracy of the developed solver is discussed
- …