2,054 research outputs found
Implementation of a low-mach number modification for high-order finite-volume schemes for arbitrary hybrid unstructured meshes
An implementation of a novel low-mach number treatment for high-order finite-volume schemes using arbitrary hybrid unstructured meshes is presented in this paper. Low-Mach order modifications for Godunov type finite-volume schemes have been implemented successfully for structured and unstructured meshes, however the methods break down for hybrid mesh topologies containing multiple element types. The modification is applied to the UCNS3D finite-volume framework for compressible flow configurations, which have been shown as very capable of handling any type of grid topology. The numerical methods under consideration are the Monotonic Upstream-Centered Scheme for Conservation Laws (MUSCL) and the Weighted Essentially Non-Oscillatory (WENO) schemes for two-dimensional mixed-element type unstructured meshes. In the present study the HLLC Approximate Riemann Solver is used with an explicit TVD Runge-Kutta 3rd-order method due to its excellent scalability. These schemes (up to 5th-order) are applied to well established two-dimensional and three-dimensional test cases. The challenges that occur when applying these methods to low-mach flow configurations is thoroughly analysed and possible improvements and further test cases are suggested
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)
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
Addressing the challenges of implementation of high-order finite volume schemes for atmospheric dynamics of unstructured meshes
The solution of the non-hydrostatic compressible Euler equations using Weighted Essentially Non-Oscillatory (WENO) schemes in two and three-dimensional unstructured meshes, is presented. Their key characteristics are their simplicity; accuracy; robustness; non-oscillatory properties; versatility in handling any type of grid topology; computational and parallel efficiency. Their defining characteristic is a non-linear combination of a series of high-order reconstruction polynomials arising from a series of reconstruction stencils. In the present study an explicit TVD Runge-Kutta 3rd -order method is employed due to its lower computational resources requirement compared to implicit type time advancement methods. The WENO schemes (up to 5th -order) are applied to the two dimensional and three dimensional test cases: a 2D rising
Heterogeneous Computing on Mixed Unstructured Grids with PyFR
PyFR is an open-source high-order accurate computational fluid dynamics
solver for mixed unstructured grids that can target a range of hardware
platforms from a single codebase. In this paper we demonstrate the ability of
PyFR to perform high-order accurate unsteady simulations of flow on mixed
unstructured grids using heterogeneous multi-node hardware. Specifically, after
benchmarking single-node performance for various platforms, PyFR v0.2.2 is used
to undertake simulations of unsteady flow over a circular cylinder at Reynolds
number 3 900 using a mixed unstructured grid of prismatic and tetrahedral
elements on a desktop workstation containing an Intel Xeon E5-2697 v2 CPU, an
NVIDIA Tesla K40c GPU, and an AMD FirePro W9100 GPU. Both the performance and
accuracy of PyFR are assessed. PyFR v0.2.2 is freely available under a 3-Clause
New Style BSD license (see www.pyfr.org).Comment: 21 pages, 9 figures, 6 table
A High-Order Unifying Discontinuous Formulation for the Navier-Stokes Equations on 3D Mixed Grids
This is the published version. Copyright 2011 © EDP SciencesThe newly developed unifying discontinuous formulation named the correction procedure via reconstruction (CPR) for conservation laws is extended to solve the Navier-Stokes equations for 3D mixed grids. In the current development, tetrahedrons and triangular prisms are considered. The CPR method can unify several popular high order methods including the discontinuous Galerkin and the spectral volume methods into a more efficient differential form. By selecting the solution points to coincide with the flux points, solution reconstruction can be completely avoided. Accuracy studies confirmed that the optimal order of accuracy can be achieved with the method. Several benchmark test cases are computed by solving the Euler and compressible Navier-Stokes equations to demonstrate its performance
- …