An extended range of stable-symmetric-conservative Flux Reconstruction correction functions

The Flux Reconstruction (FR) approach offers an efficient route to achieving high-order accuracy on unstructured grids. Additionally, FR offers a flexible framework for defining a range of numerical schemes in terms of so-called FR correction functions. Recently, a one-parameter family of FR correction functions were identified that lead to stable schemes for 1D linear advection problems. In this study we develop a procedure for identifying an extended range of stable, symmetric, and conservative FR correction functions. The procedure is applied to identify ranges of such correction functions for various orders of accuracy. Numerical experiments are undertaken, and the results found to be in agreement with the theoretical findings

An Analysis of Solution Point Coordinates for Flux Reconstruction Schemes on Tetrahedral Elements

The flux reconstruction (FR) approach offers an efficient route to high-order accuracy on unstructured grids. In this work we study the effect of solution point placement on the stability and accuracy of FR schemes on tetrahedral grids. To accomplish this we generate a large number of solution point candidates that satisfy various criteria at polynomial orders ℘=3,4,5℘=3,4,5 . We then proceed to assess their properties by using them to solve the non-linear Euler equations on both structured and unstructured meshes. The results demonstrate that the location of the solution points is important in terms of both the stability and accuracy. Across a range of cases it is possible to outperform the solution points of Shunn and Ham for specific problems. However, there appears to be a degree of problem-dependence with regards to the optimal point set, and hence overall it is concluded that the Shunn and Ham points offer a good compromise in terms of practical utility

Towards green aviation with Python at petascale

Accurate simulation of unsteady turbulent flow is critical for improved design of greener aircraft that are quieter and more fuel-efficient. We demonstrate application of PyFR, a Python based computational fluid dynamics solver, to petascale simulation of such flow problems. Rationale behind algorithmic choices, which offer increased levels of accuracy and enable sustained computation at up to 58% of peak DP-FLOP/s on unstruc- tured grids, will be discussed in the context of modern hardware. A range of software innovations will also be detailed, including use of runtime code generation, which enables PyFR to efficiently target multiple platforms, including heterogeneous systems, via a single implemen- tation. Finally, results will be presented from a full- scale simulation of flow over a low-pressure turbine blade cascade, along with weak/strong scaling statistics from the Piz Daint and Titan supercomputers, and performance data demonstrating sustained computation at up to 13.7 DP-PFLOP/s

Tutorial on Hybridizable Discontinous Galerkin (HDG) for second-order elliptic problems

The HDG is a new class of discontinuous Galerkin (DG) methods that shares favorable properties with classical mixed methods such as the well known Raviart-Thomas methods. In particular, HDG provides optimal convergence of both the primal and the dual variables of the mixed formulation. This property enables the construction of superconvergent solutions, contrary to other popular DG methods. In addition, its reduced computational cost, compared to other DG methods, has made HDG an attractive alternative for solving problems governed by partial differential equations. A tutorial on HDG for the numerical solution of second-order elliptic problems is presented. Particular emphasis is placed on providing all the necessary details for the implementation of HDG methods.Peer ReviewedPreprin

A survey of free software for the design, analysis, modelling, and simulation of an unmanned aerial vehicle

The objective of this paper is to analyze free software for the design, analysis, modelling, and simulation of an unmanned aerial vehicle (UAV). Free software is the best choice when the reduction of production costs is necessary; nevertheless, the quality of free software may vary. This paper probably does not include all of the free software, but tries to describe or mention at least the most interesting programs. The first part of this paper summarizes the essential knowledge about UAVs, including the fundamentals of flight mechanics and aerodynamics, and the structure of a UAV system. The second section generally explains the modelling and simulation of a UAV. In the main section, more than 50 free programs for the design, analysis, modelling, and simulation of a UAV are described. Although the selection of the free software has been focused on small subsonic UAVs, the software can also be used for other categories of aircraft in some cases; e.g. for MAVs and large gliders. The applications with an historical importance are also included. Finally, the results of the analysis are evaluated and discussed—a block diagram of the free software is presented, possible connections between the programs are outlined, and future improvements of the free software are suggested. © 2015, CIMNE, Barcelona, Spain.Internal Grant Agency of Tomas Bata University in Zlin [IGA/FAI/2015/001, IGA/FAI/2014/006

Inline vector compression for computational physics

A novel inline data compression method is presented for single-precision vectors in three dimensions. The primary application of the method is for accelerating computational physics calculations where the throughput is bound by memory bandwidth. The scheme employs spherical polar coordinates, angle quantisation, and a bespoke floating-point representation of the magnitude to achieve a fixed compression ratio of 1.5. The anisotropy of this method is considered, along with companding and fractional splitting techniques to improve the efficiency of the representation. We evaluate the scheme numerically within the context of high-order computational fluid dynamics. For both the isentropic convecting vortex and the Taylor–Green vortex test cases, the results are found to be comparable to those without compression. Performance is evaluated for a vector addition kernel on an NVIDIA Titan V GPU; it is demonstrated that a speedup of 1.5 can be achieved

On Fourier analysis of polynomial multigrid for arbitrary multi-stage cycles

The Fourier analysis of the \emph{p}-multigrid acceleration technique is considered for a dual-time scheme applied to the advection-diffusion equation with various cycle configurations. It is found that improved convergence can be achieved through \emph{V}-cycle asymmetry where additional prolongation smoothing is applied. Experiments conducted on the artificial compressibility formulation of the Navier--Stokes equations found that these analytic findings could be observed numerically in the pressure residual, whereas velocity terms---which are more hyperbolic in character---benefited primarily from increased pseudo-time steps