690 research outputs found
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
Hybridized Formulations of Flux Reconstruction Schemes for Advection-Diffusion Problems
We present the hybridization of flux reconstruction methods for
advection-diffusion problems. Hybridization introduces a new variable into the
problem so that it can be reduced via static condensation. This allows the
solution of implicit discretizations to be done more efficiently. We derive an
energy statement from a stability analysis considering a range of correction
functions on hybridized and embedded flux reconstruction schemes. Then, we
establish connections to standard formulations. We devise a post-processing
scheme that leverages existing flux reconstruction operators to enhance
accuracy for diffusion-dominated problems. Results show that the implicit
convergence of these methods for advection-diffusion problems can result in
performance benefits of over an order of magnitude. In addition, we observe
that the superconvergence property of hybridized methods can be extended to the
family of FR schemes for a range of correction functions
PyFR: An Open Source Framework for Solving Advection-Diffusion Type Problems on Streaming Architectures using the Flux Reconstruction Approach
High-order numerical methods for unstructured grids combine the superior
accuracy of high-order spectral or finite difference methods with the geometric
flexibility of low-order finite volume or finite element schemes. The Flux
Reconstruction (FR) approach unifies various high-order schemes for
unstructured grids within a single framework. Additionally, the FR approach
exhibits a significant degree of element locality, and is thus able to run
efficiently on modern streaming architectures, such as Graphical Processing
Units (GPUs). The aforementioned properties of FR mean it offers a promising
route to performing affordable, and hence industrially relevant,
scale-resolving simulations of hitherto intractable unsteady flows within the
vicinity of real-world engineering geometries. In this paper we present PyFR,
an open-source Python based framework for solving advection-diffusion type
problems on streaming architectures using the FR approach. The framework is
designed to solve a range of governing systems on mixed unstructured grids
containing various element types. It is also designed to target a range of
hardware platforms via use of an in-built domain specific language based on the
Mako templating engine. The current release of PyFR is able to solve the
compressible Euler and Navier-Stokes equations on grids of quadrilateral and
triangular elements in two dimensions, and hexahedral elements in three
dimensions, targeting clusters of CPUs, and NVIDIA GPUs. Results are presented
for various benchmark flow problems, single-node performance is discussed, and
scalability of the code is demonstrated on up to 104 NVIDIA M2090 GPUs. The
software is freely available under a 3-Clause New Style BSD license (see
www.pyfr.org)
High-Order Methods for Computational Fluid Dynamics: A Brief Review of Compact Differential Formulations on Unstructured Grids
Popular high-order schemes with compact stencils for Computational Fluid Dynamics (CFD) include Discontinuous Galerkin (DG), Spectral Difference (SD), and Spectral Volume (SV) methods. The recently proposed Flux Reconstruction (FR) approach or Correction Procedure using Reconstruction (CPR) is based on a differential formulation and provides a unifying framework for these high-order schemes. Here we present a brief review of recent developments for the FR/CPR schemes as well as some pacing items
Dealiasing techniques for high-order spectral element methods on regular and irregular grids
High-order methods are becoming increasingly attractive in both academia and industry, especially in the context of computational fluid dynamics. However, before they can be more widely adopted, issues such as lack of robustness in terms of numerical stability need to be addressed, particularly when treating industrial-type problems where challenging geometries and a wide range of physical scales, typically due to high Reynolds numbers, need to be taken into account. One source of instability is aliasing effects which arise from the nonlinearity of the underlying problem. In this work we detail two dealiasing strategies based on the concept of consistent integration. The first uses a localised approach, which is useful when the nonlinearities only arise in parts of the problem. The second is based on the more traditional approach of using a higher quadrature. The main goal of both dealiasing techniques is to improve the robustness of high order spectral element methods, thereby reducing aliasing-driven instabilities. We demonstrate how these two strategies can be effectively applied to both continuous and discontinuous discretisations, where, in the latter, both volumetric and interface approximations must be considered. We show the key features of each dealiasing technique applied to the scalar conservation law with numerical examples and we highlight the main differences in terms of implementation between continuous and discontinuous spatial discretisations
Flux form Semi-Lagrangian methods for parabolic problems
A semi-Lagrangian method for parabolic problems is proposed, that extends
previous work by the authors to achieve a fully conservative, flux-form
discretization of linear and nonlinear diffusion equations. A basic consistency
and convergence analysis are proposed. Numerical examples validate the proposed
method and display its potential for consistent semi-Lagrangian discretization
of advection--diffusion and nonlinear parabolic problems
Aspects of Unstructured Grids and Finite-Volume Solvers for the Euler and Navier-Stokes Equations
One of the major achievements in engineering science has been the development of computer algorithms for solving nonlinear differential equations such as the Navier-Stokes equations. In the past, limited computer resources have motivated the development of efficient numerical schemes in computational fluid dynamics (CFD) utilizing structured meshes. The use of structured meshes greatly simplifies the implementation of CFD algorithms on conventional computers. Unstructured grids on the other hand offer an alternative to modeling complex geometries. Unstructured meshes have irregular connectivity and usually contain combinations of triangles, quadrilaterals, tetrahedra, and hexahedra. The generation and use of unstructured grids poses new challenges in CFD. The purpose of this note is to present recent developments in the unstructured grid generation and flow solution technology
Issures in Discontinuous High-Order Methods: Broadband Wave Computation and Viscous Boundary Layer Resolution
A new discontinuous formulation named Correction Procedure via Reconstruction (CPR) was developed for conservation laws. CPR is an efficient nodal differential formulation unifying the discontinuous Galerkin (DG), spectral volume (SV) and spectral difference (SD) methods, is easy to implement. In this thesis, we deal with two issues: the efficient computation of broadband waves, and the proper resolution of a viscous boundary layer with the high-order CPR method.
A hybrid discontinuous space including polynomial and Fourier bases is employed in the CPR formulation in order to compute broad-band waves. The polynomial bases are used to achieve a certain order of accuracy, while the Fourier bases are able to exactly resolve waves at a certain frequency. Free-parameters introduced in the Fourier bases are optimized in order to minimize both dispersion and dissipation errors by mimicking the dispersion-relation-preserving (DRP) method for a one-dimensional wave problem. For the one-dimensional wave problem, the dispersion and dissipation properties and the optimization procedure are investigated through a wave propagation analysis. The optimization procedure is verified with a wave propagation analysis and several numerical tests. The two-dimensional wave behavior is investigated through a wave propagation analysis and the wave propagation properties are verified with a numerical test of the two-dimensional acoustic wave equation.
In order to understand the mesh size requirement to resolve a viscous boundary layer using CPR method, grid resolution studies are performed. . It is well known that the mesh size, which is defined from non-dimensional wall distance y^+=1, gives accepted results to simulate viscous boundary layer problem for 2nd order finite volume method. For high-order CPR formulation, what grid size y^+ is required to match the results from the 2nd order finite volume method with y^+=1. 1D and 2D burger\u27s equation are used as the viscous boundary layer model problem. Skin friction is used as the indicator of accuracy for the resolution of a boundary layer.
Keywords: (Correction Procedure via Reconstruction), A Hybrid Discontinuous Space, Wave Propagation Analysis, Grid Resolution Study, Method of Manufactured Solution
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
On the Properties of Energy Stable Flux Reconstruction Schemes for Implicit Large Eddy Simulation
We begin by investigating the stability, order of accuracy, and dispersion and dissipation characteristics of the extended range of energy stable flux reconstruction (E-ESFR) schemes in the context of implicit large eddy simulation (ILES). We proceed to demonstrate that subsets of the E-ESFR schemes are more stable than collocation nodal discontinuous Galerkin methods recovered with the flux reconstruction approach (FRDG) for marginally-resolved ILES simulations of the Taylor-Green vortex. These schemes are shown to have reduced dissipation and dispersion errors relative to FRDG schemes of the same polynomial degree and, simultaneously, have increased CourantFriedrichs-Lewy (CFL) limits. Finally, we simulate turbulent flow over an SD7003 aerofoil using two of the most stable E-ESFR schemes identified by the aforementioned Taylor-Green vortex experiments. Results demonstrate that subsets of E-ESFR schemes appear more stable than the commonly used FRDG method, have increased CFL limits, and are suitable for ILES of complex turbulent flows on unstructured grids
- …