Nektar++: Design and implementation of an implicit, spectral/hp element, compressible flow solver using a Jacobian-free Newton Krylov approach

File Description SizeFormat 
2002.04222v1.pdfAccepted version3.5 MBAdobe PDFView/Open
Title: Nektar++: Design and implementation of an implicit, spectral/hp element, compressible flow solver using a Jacobian-free Newton Krylov approach
Authors: Yan, Z-G
Pan, Y
Castiglioni, G
Hillewaert, K
Peiró, J
Moxey, D
Sherwin, SJ
Item Type: Journal Article
Abstract: At high Reynolds numbers the use of explicit in time compressible flow simulations with spectral/ element discretization can become significantly limited by time step. To alleviate this limitation we extend the capability of the spectral/ element open-source software framework, Nektar++, to include an implicit discontinuous Galerkin compressible flow solver. The integration in time is carried out by a singly diagonally implicit Runge–Kutta method. The non-linear system arising from the implicit time integration is iteratively solved by the Jacobian-free Newton Krylov (JFNK) method. A favorable feature of the JFNK approach is its extensive use of the explicit operators available from the previous explicit in time implementation. The functionalities of different building blocks of the implicit solver are analyzed from the point of view of software design and placed in appropriate hierarchical levels in the C++ libraries. In the detailed implementation, the contributions of different parts of the solver to computational cost, memory consumption and programming complexity are also analyzed. A combination of analytical and numerical methods is adopted to simplify the programming complexity in forming the preconditioning matrix. The solver is verified and tested using cases such as manufactured compressible Poiseuille flow, Taylor–Green vortex, turbulent flow over a circular cylinder at and shock wave boundary-layer interaction. The results show that the implicit solver can speed-up the simulations while maintaining good simulation accuracy.
Issue Date: 1-Jan-2021
Date of Acceptance: 14-Mar-2020
URI: http://hdl.handle.net/10044/1/80588
DOI: 10.1016/j.camwa.2020.03.009
ISSN: 0898-1221
Publisher: Elsevier BV
Start Page: 351
End Page: 372
Journal / Book Title: Computers & Mathematics with Applications
Volume: 81
Copyright Statement: © 2020 Elsevier Ltd. All rights reserved. This manuscript is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International Licence http://creativecommons.org/licenses/by-nc-nd/4.0/
Keywords: Science & Technology
Physical Sciences
Mathematics, Applied
Mathematics
Nektar plus
Implicit time integration
Spectral/hp element
Discontinuous Galerkin
Jacobian-free Newton Krylov
DISCONTINUOUS GALERKIN METHOD
NAVIER-STOKES EQUATIONS
TENSOR-PRODUCT PRECONDITIONERS
SIMULATION
SCHEMES
STRATEGY
physics.flu-dyn
physics.flu-dyn
physics.comp-ph
76N99
Numerical & Computational Mathematics
01 Mathematical Sciences
08 Information and Computing Sciences
15 Commerce, Management, Tourism and Services
Publication Status: Published
Online Publication Date: 2020-04-21
Appears in Collections:Faculty of Engineering
Aeronautics