5,268 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
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)
Hybrid Spectral Difference/Embedded Finite Volume Method for Conservation Laws
A novel hybrid spectral difference/embedded finite volume method is
introduced in order to apply a discontinuous high-order method for large scale
engineering applications involving discontinuities in the flows with complex
geometries. In the proposed hybrid approach, the finite volume (FV) element,
consisting of structured FV subcells, is embedded in the base hexahedral
element containing discontinuity, and an FV based high-order shock-capturing
scheme is employed to overcome the Gibbs phenomena. Thus, a discontinuity is
captured at the resolution of FV subcells within an embedded FV element. In the
smooth flow region, the SD element is used in the base hexahedral element.
Then, the governing equations are solved by the SD method. The SD method is
chosen for its low numerical dissipation and computational efficiency
preserving high-order accurate solutions. The coupling between the SD element
and the FV element is achieved by the globally conserved mortar method. In this
paper, the 5th-order WENO scheme with the characteristic decomposition is
employed as the shock-capturing scheme in the embedded FV element, and the
5th-order SD method is used in the smooth flow field.
The order of accuracy study and various 1D and 2D test cases are carried out,
which involve the discontinuities and vortex flows. Overall, it is shown that
the proposed hybrid method results in comparable or better simulation results
compared with the standalone WENO scheme when the same number of solution DOF
is considered in both SD and FV elements.Comment: 27 pages, 17 figures, 2 tables, Accepted for publication in the
Journal of Computational Physics, April 201
GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
We present two new Lagrangian methods for hydrodynamics, in a systematic
comparison with moving-mesh, SPH, and stationary (non-moving) grid methods. The
new methods are designed to simultaneously capture advantages of both
smoothed-particle hydrodynamics (SPH) and grid-based/adaptive mesh refinement
(AMR) schemes. They are based on a kernel discretization of the volume coupled
to a high-order matrix gradient estimator and a Riemann solver acting over the
volume 'overlap.' We implement and test a parallel, second-order version of the
method with self-gravity & cosmological integration, in the code GIZMO: this
maintains exact mass, energy and momentum conservation; exhibits superior
angular momentum conservation compared to all other methods we study; does not
require 'artificial diffusion' terms; and allows the fluid elements to move
with the flow so resolution is automatically adaptive. We consider a large
suite of test problems, and find that on all problems the new methods appear
competitive with moving-mesh schemes, with some advantages (particularly in
angular momentum conservation), at the cost of enhanced noise. The new methods
have many advantages vs. SPH: proper convergence, good capturing of
fluid-mixing instabilities, dramatically reduced 'particle noise' & numerical
viscosity, more accurate sub-sonic flow evolution, & sharp shock-capturing.
Advantages vs. non-moving meshes include: automatic adaptivity, dramatically
reduced advection errors & numerical overmixing, velocity-independent errors,
accurate coupling to gravity, good angular momentum conservation and
elimination of 'grid alignment' effects. We can, for example, follow hundreds
of orbits of gaseous disks, while AMR and SPH methods break down in a few
orbits. However, fixed meshes minimize 'grid noise.' These differences are
important for a range of astrophysical problems.Comment: 57 pages, 33 figures. MNRAS. A public version of the GIZMO code,
user's guide, test problem setups, and movies are available at
http://www.tapir.caltech.edu/~phopkins/Site/GIZMO.htm
ECHO: an Eulerian Conservative High Order scheme for general relativistic magnetohydrodynamics and magnetodynamics
We present a new numerical code, ECHO, based on an Eulerian Conservative High
Order scheme for time dependent three-dimensional general relativistic
magnetohydrodynamics (GRMHD) and magnetodynamics (GRMD). ECHO is aimed at
providing a shock-capturing conservative method able to work at an arbitrary
level of formal accuracy (for smooth flows), where the other existing GRMHD and
GRMD schemes yield an overall second order at most. Moreover, our goal is to
present a general framework, based on the 3+1 Eulerian formalism, allowing for
different sets of equations, different algorithms, and working in a generic
space-time metric, so that ECHO may be easily coupled to any solver for
Einstein's equations. Various high order reconstruction methods are implemented
and a two-wave approximate Riemann solver is used. The induction equation is
treated by adopting the Upwind Constrained Transport (UCT) procedures,
appropriate to preserve the divergence-free condition of the magnetic field in
shock-capturing methods. The limiting case of magnetodynamics (also known as
force-free degenerate electrodynamics) is implemented by simply replacing the
fluid velocity with the electromagnetic drift velocity and by neglecting the
matter contribution to the stress tensor. ECHO is particularly accurate,
efficient, versatile, and robust. It has been tested against several
astrophysical applications, including a novel test on the propagation of large
amplitude circularly polarized Alfven waves. In particular, we show that
reconstruction based on a Monotonicity Preserving filter applied to a fixed
5-point stencil gives highly accurate results for smooth solutions, both in
flat and curved metric (up to the nominal fifth order), while at the same time
providing sharp profiles in tests involving discontinuities.Comment: 20 pages, revised version submitted to A&
- …