×

zbMATH — the first resource for mathematics

PyFR: an open source framework for solving advection-diffusion type problems on streaming architectures using the flux reconstruction approach. (English) Zbl 1348.65005
Summary: 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 http://www.pyfr.org).

MSC:
65-04 Software, source code, etc. for problems pertaining to numerical analysis
65M60 Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs
65Y05 Parallel numerical computation
PDF BibTeX XML Cite
Full Text: DOI arXiv
References:
[1] Reed, W. H.; Hill, T. R., Triangular mesh methods for the neutron transport equation. technical report LA-UR-73-479, (1973), Los Alamos Scientific Laboratory
[2] Kopriva, David A.; Kolias, John H., A conservative staggered-grid Chebyshev multidomain method for compressible flows, J. Comput. Phys., 125, 1, 244-261, (1996) · Zbl 0847.76069
[3] Sun, Yuzhi; Wang, Zhi Jian; Liu, Yen, High-order multidomain spectral difference method for the Navier-Stokes equations on unstructured hexahedral grids, Commun. Comput. Phys., 2, 2, 310-333, (2007) · Zbl 1164.76360
[4] Huynh, H. T., A flux reconstruction approach to high-order schemes including discontinuous Galerkin methods, AIAA paper, 4079, 2007, (2007)
[5] Hesthaven, Jan S.; Warburton, Tim, Nodal discontinuous Galerkin methods: algorithms, analysis, and applications, vol. 54, (2008), Springer Verlag New York · Zbl 1134.65068
[6] Vincent, P. E.; Castonguay, P.; Jameson, A., A new class of high-order energy stable flux reconstruction schemes, J. Sci. Comput., 47, 1, 50-72, (2011) · Zbl 1433.76094
[7] Castonguay, P.; Vincent, P. E.; Jameson, A., A new class of high-order energy stable flux reconstruction schemes for triangular elements, J. Sci. Comput., (2011) · Zbl 1433.76094
[8] Patrice Castonguay, P.E. Vincent, Antony Jameson, Application of high-order energy stable flux reconstruction schemes to the Euler equations, in: 49th AIAA Aerospace Sciences Meeting, volume 686, 2011.
[9] Jameson, A.; Vincent, P. E.; Castonguay, P., On the non-linear stability of flux reconstruction schemes, J. Sci. Comput., 50, 2, 434-445, (2011) · Zbl 1457.65116
[10] Vincent, P. E.; Castonguay, P.; Jameson, A., Insights from von Neumann analysis of high-order flux reconstruction schemes, J. Comput. Phys., 230, 22, 8134-8154, (2011) · Zbl 1343.65117
[11] Castonguay, Patrice; Vincent, Peter E.; Jameson, Antony, A new class of high-order energy stable flux reconstruction schemes for triangular elements, J. Sci. Comput., 51, 1, 224-256, (2012) · Zbl 1457.65101
[12] Williams, D. M.; Castonguay, P.; Vincent, P. E.; Jameson, A., Energy stable flux reconstruction schemes for advection-diffusion problems on triangles, J. Comput. Phys., (2013) · Zbl 1349.65528
[13] Castonguay, P.; Williams, D. M.; Vincent, P. E.; Jameson, A., Energy stable flux reconstruction schemes for advection-diffusion problems, Comput. Methods Appl. Mech. Engrg., (2013) · Zbl 1286.65119
[14] Williams, D. M.; Jameson, A., Energy stable flux reconstruction schemes for advection-diffusion problems on tetrahedra, J. Sci. Comput., 1-39, (2013)
[15] Kopriva, David A., A staggered-grid multidomain spectral method for the compressible Navier-Stokes equations, J. Comput. Phys., 143, 1, 125-158, (1998) · Zbl 0921.76121
[16] SymPy Development Team. Sympy: Python library for symbolic mathematics, 2013.
[17] Michael Bayer, Mako: Templates for python, 2013.
[18] Klöckner, Andreas; Pinto, Nicolas; Lee, Yunsup; Catanzaro, Bryan; Ivanov, Paul; Fasih, Ahmed, Pycuda and pyopencl: A scripting-based approach to gpu run-time code generation, Parallel Comput., 38, 3, 157-174, (2012)
[19] Lisandro Dalcin, mpi4py: Mpi for python, 2013.
[20] Toro, Eleuterio F., Riemann solvers and numerical methods for fluid dynamics: a practical introduction, (2009), Springer · Zbl 1227.76006
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.