6,247 research outputs found
Spectral/hp element methods: recent developments, applications, and perspectives
The spectral/hp element method combines the geometric flexibility of the
classical h-type finite element technique with the desirable numerical
properties of spectral methods, employing high-degree piecewise polynomial
basis functions on coarse finite element-type meshes. The spatial approximation
is based upon orthogonal polynomials, such as Legendre or Chebychev
polynomials, modified to accommodate C0-continuous expansions. Computationally
and theoretically, by increasing the polynomial order p, high-precision
solutions and fast convergence can be obtained and, in particular, under
certain regularity assumptions an exponential reduction in approximation error
between numerical and exact solutions can be achieved. This method has now been
applied in many simulation studies of both fundamental and practical
engineering flows. This paper briefly describes the formulation of the
spectral/hp element method and provides an overview of its application to
computational fluid dynamics. In particular, it focuses on the use the
spectral/hp element method in transitional flows and ocean engineering.
Finally, some of the major challenges to be overcome in order to use the
spectral/hp element method in more complex science and engineering applications
are discussed
A matrix-free high-order discontinuous Galerkin compressible Navier-Stokes solver: A performance comparison of compressible and incompressible formulations for turbulent incompressible flows
Both compressible and incompressible Navier-Stokes solvers can be used and
are used to solve incompressible turbulent flow problems. In the compressible
case, the Mach number is then considered as a solver parameter that is set to a
small value, , in order to mimic incompressible flows.
This strategy is widely used for high-order discontinuous Galerkin
discretizations of the compressible Navier-Stokes equations. The present work
raises the question regarding the computational efficiency of compressible DG
solvers as compared to a genuinely incompressible formulation. Our
contributions to the state-of-the-art are twofold: Firstly, we present a
high-performance discontinuous Galerkin solver for the compressible
Navier-Stokes equations based on a highly efficient matrix-free implementation
that targets modern cache-based multicore architectures. The performance
results presented in this work focus on the node-level performance and our
results suggest that there is great potential for further performance
improvements for current state-of-the-art discontinuous Galerkin
implementations of the compressible Navier-Stokes equations. Secondly, this
compressible Navier-Stokes solver is put into perspective by comparing it to an
incompressible DG solver that uses the same matrix-free implementation. We
discuss algorithmic differences between both solution strategies and present an
in-depth numerical investigation of the performance. The considered benchmark
test cases are the three-dimensional Taylor-Green vortex problem as a
representative of transitional flows and the turbulent channel flow problem as
a representative of wall-bounded turbulent flows
On the stability of projection methods for the incompressible Navier-Stokes equations based on high-order discontinuous Galerkin discretizations
The present paper deals with the numerical solution of the incompressible
Navier-Stokes equations using high-order discontinuous Galerkin (DG) methods
for discretization in space. For DG methods applied to the dual splitting
projection method, instabilities have recently been reported that occur for
coarse spatial resolutions and small time step sizes. By means of numerical
investigation we give evidence that these instabilities are related to the
discontinuous Galerkin formulation of the velocity divergence term and the
pressure gradient term that couple velocity and pressure. Integration by parts
of these terms with a suitable definition of boundary conditions is required in
order to obtain a stable and robust method. Since the intermediate velocity
field does not fulfill the boundary conditions prescribed for the velocity, a
consistent boundary condition is derived from the convective step of the dual
splitting scheme to ensure high-order accuracy with respect to the temporal
discretization. This new formulation is stable in the limit of small time steps
for both equal-order and mixed-order polynomial approximations. Although the
dual splitting scheme itself includes inf-sup stabilizing contributions, we
demonstrate that spurious pressure oscillations appear for equal-order
polynomials and small time steps highlighting the necessity to consider inf-sup
stability explicitly.Comment: 31 page
Finite element formulation of general boundary conditions for incompressible flows
We study the finite element formulation of general boundary conditions for
incompressible flow problems. Distinguishing between the contributions from the
inviscid and viscid parts of the equations, we use Nitsche's method to develop
a discrete weighted weak formulation valid for all values of the viscosity
parameter, including the limit case of the Euler equations. In order to control
the discrete kinetic energy, additional consistent terms are introduced. We
treat the limit case as a (degenerate) system of hyperbolic equations, using a
balanced spectral decomposition of the flux Jacobian matrix, in analogy with
compressible flows. Then, following the theory of Friedrich's systems, the
natural characteristic boundary condition is generalized to the considered
physical boundary conditions. Several numerical experiments, including standard
benchmarks for viscous flows as well as inviscid flows are presented
High-order DG solvers for under-resolved turbulent incompressible flows: A comparison of and (div) methods
The accurate numerical simulation of turbulent incompressible flows is a
challenging topic in computational fluid dynamics. For discretisation methods
to be robust in the under-resolved regime, mass conservation as well as energy
stability are key ingredients to obtain robust and accurate discretisations.
Recently, two approaches have been proposed in the context of high-order
discontinuous Galerkin (DG) discretisations that address these aspects
differently. On the one hand, standard -based DG discretisations enforce
mass conservation and energy stability weakly by the use of additional
stabilisation terms. On the other hand, pointwise divergence-free
-conforming approaches ensure exact mass conservation
and energy stability by the use of tailored finite element function spaces. The
present work raises the question whether and to which extent these two
approaches are equivalent when applied to under-resolved turbulent flows. This
comparative study highlights similarities and differences of these two
approaches. The numerical results emphasise that both discretisation strategies
are promising for under-resolved simulations of turbulent flows due to their
inherent dissipation mechanisms.Comment: 24 pages, 13 figure
- …