6,696 research outputs found
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
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
Non-modal analysis of spectral element methods: Towards accurate and robust large-eddy simulations
We introduce a \textit{non-modal} analysis technique that characterizes the
diffusion properties of spectral element methods for linear
convection-diffusion systems. While strictly speaking only valid for linear
problems, the analysis is devised so that it can give critical insights on two
questions: (i) Why do spectral element methods suffer from stability issues in
under-resolved computations of nonlinear problems? And, (ii) why do they
successfully predict under-resolved turbulent flows even without a
subgrid-scale model? The answer to these two questions can in turn provide
crucial guidelines to construct more robust and accurate schemes for complex
under-resolved flows, commonly found in industrial applications. For
illustration purposes, this analysis technique is applied to the hybridized
discontinuous Galerkin methods as representatives of spectral element methods.
The effect of the polynomial order, the upwinding parameter and the P\'eclet
number on the so-called \textit{short-term diffusion} of the scheme are
investigated. From a purely non-modal analysis point of view, polynomial orders
between and with standard upwinding are well suited for under-resolved
turbulence simulations. For lower polynomial orders, diffusion is introduced in
scales that are much larger than the grid resolution. For higher polynomial
orders, as well as for strong under/over-upwinding, robustness issues can be
expected. The non-modal analysis results are then tested against under-resolved
turbulence simulations of the Burgers, Euler and Navier-Stokes equations. While
devised in the linear setting, our non-modal analysis succeeds to predict the
behavior of the scheme in the nonlinear problems considered
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
A CutFEM method for two-phase flow problems
In this article, we present a cut finite element method for two-phase
Navier-Stokes flows. The main feature of the method is the formulation of a
unified continuous interior penalty stabilisation approach for, on the one
hand, stabilising advection and the pressure-velocity coupling and, on the
other hand, stabilising the cut region. The accuracy of the algorithm is
enhanced by the development of extended fictitious domains to guarantee a well
defined velocity from previous time steps in the current geometry. Finally, the
robustness of the moving-interface algorithm is further improved by the
introduction of a curvature smoothing technique that reduces spurious
velocities. The algorithm is shown to perform remarkably well for low capillary
number flows, and is a first step towards flexible and robust CutFEM algorithms
for the simulation of microfluidic devices
Simulation of flows with violent free surface motion and moving objects using unstructured grids
This is the peer reviewed version of the following article: [Löhner, R. , Yang, C. and Oñate, E. (2007), Simulation of flows with violent free surface motion and moving objects using unstructured grids. Int. J. Numer. Meth. Fluids, 53: 1315-1338. doi:10.1002/fld.1244], which has been published in final form at https://doi.org/10.1002/fld.1244. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving.A volume of fluid (VOF) technique has been developed and coupled with an incompressible Euler/NavierâStokes solver operating on adaptive, unstructured grids to simulate the interactions of extreme waves and three-dimensional structures. The present implementation follows the classic VOF implementation for the liquidâgas system, considering only the liquid phase. Extrapolation algorithms are used to obtain velocities and pressure in the gas region near the free surface. The VOF technique is validated against the classic dam-break problem, as well as series of 2D sloshing experiments and results from SPH calculations. These and a series of other examples demonstrate that the ability of the present approach to simulate violent free surface flows with strong nonlinear behaviour.Peer ReviewedPostprint (author's final draft
- âŠ