Effective high-order energy stable flux reconstruction methods for first-order hyperbolic linear and nonlinear systems

Abstract The paper devises high-order energy stable flux reconstruction method towards first-order hyperbolic linear advection-diffusion and nonlinear Navier-Stokes problems (HFR). The underlying idea is to take the advantages of high-order differential framework from flux reconstruction and accurate gradient solution from first-order hyperbolic system. By means of the additional gradient equation, the linear advection-diffusion system is recast as hyperbolic formulation, which is equivalent to original equation when pseudo time derivatives are dropped, and then diffusive fluxes are directly solved by flux reconstruction. The resulting scheme has compact stencils and high accuracy without increasing computation cost. Mathematical analysis is presented to confirm that HFR method is consistent and stable for all orders of accuracy, at least for linear problems. In addition, numerical results for linear and nonlinear problems demonstrate that HFR method can achieve the desired accuracy for both primitive and gradient variable, reduce the absolute error of gradient variable, and allow a larger maximum CFL number compared to the existing methods, such as BR2 and direct flux reconstruction (DFR). It indicates that HFR method is potential in accuracy and efficiency for solving steady viscous problems.

[1]  Antony Jameson,et al.  Erratum to: On the Stability of the Flux Reconstruction Schemes on Quadrilateral Elements for the Linear Advection Equation , 2016, J. Sci. Comput..

[2]  Haijin Wang,et al.  Local discontinuous Galerkin methods with implicit-explicit time-marching for time-dependent incompressible fluid flow , 2017, Math. Comput..

[3]  David C. Del Rey Fernández,et al.  Multidimensional Summation-by-Parts Operators: General Theory and Application to Simplex Elements , 2015, SIAM J. Sci. Comput..

[4]  John H. Kolias,et al.  A CONSERVATIVE STAGGERED-GRID CHEBYSHEV MULTIDOMAIN METHOD FOR COMPRESSIBLE FLOWS , 1995 .

[5]  David L. Darmofal,et al.  p-Multigrid solution of high-order discontinuous Galerkin discretizations of the compressible Navier-Stokes equations , 2005 .

[6]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..

[7]  Alexandros Syrakos,et al.  A critical analysis of some popular methods for the discretisation of the gradient operator in finite volume methods , 2016 .

[8]  Christopher J. Roy,et al.  Review of code and solution verification procedures for computational simulation , 2005 .

[9]  Hiroaki Nishikawa,et al.  Alternative Formulations for First-, Second-, and Third-Order Hyperbolic Navier-Stokes Schemes , 2015 .

[10]  H. Nishikawa,et al.  Explicit Hyperbolic Reconstructed Discontinuous Galerkin Methods for Time-Dependent Problems , 2018, 2018 Fluid Dynamics Conference.

[11]  A. Jameson,et al.  A Direct Flux Reconstruction Scheme for Advection–Diffusion Problems on Triangular Grids , 2017, J. Sci. Comput..

[12]  H. T. Huynh,et al.  A Reconstruction Approach to High -Order Schemes Including Discontinuous Galerkin for Diffusion , 2009 .

[13]  Philipp Öffner,et al.  Extended skew-symmetric form for summation-by-parts operators and varying Jacobians , 2017, J. Comput. Phys..

[14]  Antony Jameson,et al.  On consistency and rate of convergence of Flux Reconstruction for time-dependent problems , 2017, J. Comput. Phys..

[15]  Hong Luo,et al.  Reconstructed discontinuous Galerkin methods for linear advection-diffusion equations based on first-order hyperbolic system , 2018, J. Comput. Phys..

[16]  Hiroaki Nishikawa,et al.  First-, second-, and third-order finite-volume schemes for diffusion , 2014, J. Comput. Phys..

[17]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[18]  Zachary Grant,et al.  Optimal explicit strong stability preserving Runge-Kutta methods with high linear order and optimal nonlinear order , 2014, Math. Comput..

[19]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[20]  Hiroaki Nishikawa,et al.  A first-order system approach for diffusion equation. I: Second-order residual-distribution schemes , 2007, J. Comput. Phys..

[21]  J. Romero,et al.  A Simplified Formulation of the Flux Reconstruction Method , 2016, J. Sci. Comput..

[22]  Antony Jameson,et al.  Energy stable flux reconstruction schemes for advection-diffusion problems on triangles , 2013, J. Comput. Phys..

[23]  Hiroaki Nishikawa,et al.  Divergence formulation of source term , 2012, J. Comput. Phys..

[24]  M. Carpenter,et al.  Fourth-order 2N-storage Runge-Kutta schemes , 1994 .

[25]  Magnus Svärd,et al.  On the convergence rates of energy-stable finite-difference schemes , 2019, J. Comput. Phys..

[26]  Antony Jameson,et al.  Insights from von Neumann analysis of high-order flux reconstruction schemes , 2011, J. Comput. Phys..

[27]  Hiroaki Nishikawa.  New-Generation Hyperbolic Navier-Stokes Schemes: O(1=h) Speed-Up and Accurate Viscous/Heat Fluxes , 2011 .

[28]  B. Strand.  Summation by parts for finite difference approximations for d/dx , 1994 .

[29]  Hiroaki Nishikawa,et al.  A first-order system approach for diffusion equation. II: Unification of advection and diffusion , 2010, J. Comput. Phys..

[30]  H. T. Huynh,et al.  High-Order Methods for Computational Fluid Dynamics: A Brief Review of Compact Differential Formulations on Unstructured Grids , 2013 .

[31]  Chi-Wang Shu,et al.  The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems , 1998 .

[32]  J. Hesthaven,et al.  Nodal high-order methods on unstructured grids , 2002 .

[33]  Philip L. Roe,et al.  Towards High-Order Fluctuation-Splitting Schemes for Navier-Stokes Equations , 2005 .

[34]  David C. Del Rey Fernández,et al.  Review of summation-by-parts operators with simultaneous approximation terms for the numerical solution of partial differential equations , 2014 .

[35]  Jan Nordström,et al.  Summation-by-Parts operators with minimal dispersion error for coarse grid flow calculations , 2017, J. Comput. Phys..

[36]  Brian C. Vermeire,et al.  On the behaviour of fully-discrete flux reconstruction schemes , 2017 .

[37]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[38]  P. Moin,et al.  DIRECT NUMERICAL SIMULATION: A Tool in Turbulence Research , 1998 .

[39]  Jesse Chan,et al.  Skew-Symmetric Entropy Stable Modal Discontinuous Galerkin Formulations , 2019, Journal of Scientific Computing.

[40]  Zhi Jian Wang,et al.  A unifying lifting collocation penalty formulation including the discontinuous Galerkin, spectral volume/difference methods for conservation laws on mixed grids , 2009, J. Comput. Phys..

[41]  Hiroaki Nishikawa,et al.  First, second, and third order finite-volume schemes for advection-diffusion , 2013, J. Comput. Phys..

[42]  H. T. Huynh,et al.  A Flux Reconstruction Approach to High-Order Schemes Including Discontinuous Galerkin Methods , 2007 .

[43]  Antony Jameson,et al.  High-Order Flux Reconstruction Schemes with Minimal Dispersion and Dissipation , 2015, J. Sci. Comput..

[44]  Yi Liu,et al.  Hyperbolic Navier-Stokes Method for High-Reynolds-Number Boundary Layer Flows , 2017 .

[45]  Antony Jameson,et al.  A Proof of the Stability of the Spectral Difference Method for All Orders of Accuracy , 2010, J. Sci. Comput..

[46]  Spencer J. Sherwin,et al.  Spatial eigensolution analysis of energy-stable flux reconstruction schemes and influence of the numerical flux on accuracy and robustness , 2018, J. Comput. Phys..

[47]  Alireza Mazaheri,et al.  Efficient high-order discontinuous Galerkin schemes with first-order hyperbolic advection-diffusion system approach , 2016, J. Comput. Phys..

[48]  Hiroaki Nishikawa,et al.  First, Second, and Third Order Finite-Volume Schemes for Navier-Stokes Equations , 2014 .

[49]  Antony Jameson,et al.  A New Class of High-Order Energy Stable Flux Reconstruction Schemes for Triangular Elements , 2012, J. Sci. Comput..

[50]  Antony Jameson,et al.  Energy Stable Flux Reconstruction Schemes for Advection–Diffusion Problems on Tetrahedra , 2013, Journal of Scientific Computing.

[51]  Zhi J. Wang,et al.  The direct discontinuous Galerkin (DDG) viscous flux scheme for the high order spectral volume method , 2010 .

[52]  Haijin Wang,et al.  Third order implicit-explicit Runge-Kutta local discontinuous Galerkin methods with suitable boundary treatment for convection-diffusion problems with Dirichlet boundary conditions , 2018, J. Comput. Appl. Math..

[53]  Haiyang Gao,et al.  A High-Order Lifting Collocation Penalty Formulation for the Navier-Stokes Equations on 2-D Mixed Grids , 2009 .

[54]  Magnus Svärd,et al.  Review of summation-by-parts schemes for initial-boundary-value problems , 2013, J. Comput. Phys..

[55]  Boris Diskin,et al.  A Critical Study of Agglomerated Multigrid Methods for Diffusion , 2009 .

[56]  Antony Jameson,et al.  A New Class of High-Order Energy Stable Flux Reconstruction Schemes , 2011, J. Sci. Comput..

[57]  A. Ern,et al.  Mathematical Aspects of Discontinuous Galerkin Methods , 2011 .

[58]  Yi Liu,et al.  Third-Order Inviscid and Second-Order Hyperbolic Navier-Stokes Solvers for Three-Dimensional Unsteady Inviscid and Viscous Flows , 2017 .

[59]  S. Rebay,et al.  A High-Order Accurate Discontinuous Finite Element Method for the Numerical Solution of the Compressible Navier-Stokes Equations , 1997 .