High-order flux reconstruction method for the hyperbolic formulation of the incompressible Navier-Stokes equations on unstructured grids

A high-order Flux reconstruction implementation of the hyperbolic formulation for the incompressible Navier-Stokes equation is presented. The governing equations employ Chorin's classical artificial compressibility (AC) formulation cast in hyperbolic form. Instead of splitting the second-order conservation law into two equations, one for the solution and another for the gradient, the Navier-Stokes equation is cast into a first-order hyperbolic system of equations. Including the gradients in the AC iterative process results in a significant improvement in accuracy for the pressure, velocity, and its gradients. Furthermore, this treatment allows for taking larger time-steps since the hyperbolic formulation eliminates the restriction due to diffusion. Tests using the method of manufactured solutions show that solving the conventional form of the Navier-Stokes equation lowers the order of accuracy for gradients, while the hyperbolic method is shown to provide equal orders of accuracy for both the velocity and its gradients which may be beneficial in several applications. Two- and three-dimensional benchmark tests demonstrate the superior accuracy and computational efficiency of the developed solver in comparison to the conventional method and other published works. This study shows that the developed high-order hyperbolic solver for incompressible flows is attractive due to its accuracy, stability and efficiency in solving diffusion dominated problems.

[2]  E. Erturk,et al.  Numerical solutions of 2‐D steady incompressible driven cavity flow at high Reynolds numbers , 2004, ArXiv.

[3]  Tobias Martin,et al.  Efficient WENO library for OpenFOAM , 2020, SoftwareX.

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

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

[6]  Yi Liu,et al.  Hyperbolic advection-diffusion schemes for high-Reynolds-number boundary-layer problems , 2018, J. Comput. Phys..

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

[8]  Carl Ollivier-Gooch,et al.  Obtaining and Verifying High-Order Unstructured Finite Volume Solutions to the Euler Equations , 2009 .

[9]  Hong Luo,et al.  Reconstructed discontinuous Galerkin methods for compressible flows based on a new hyperbolic Navier-Stokes system , 2020, J. Comput. Phys..

[10]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

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

[12]  Alessandro Colombo,et al.  Artificial compressibility Godunov fluxes for variable density incompressible flows , 2017, Computers & Fluids.

[13]  S. Zaghi,et al.  OFF, Open source Finite volume Fluid dynamics code: A free, high-order solver based on parallel, modular, object-oriented Fortran API , 2014, Comput. Phys. Commun..

[14]  Chunlei Liang,et al.  A high-order solver for unsteady incompressible Navier-Stokes equations using the flux reconstruction method on unstructured grids with implicit dual time stepping , 2015, J. Comput. Phys..

[15]  Freddie D. Witherden,et al.  On the development and implementation of high-order flux reconstruction schemes for computational fluid dynamics , 2015 .

[16]  Freddie D. Witherden,et al.  Heterogeneous Computing on Mixed Unstructured Grids with PyFR , 2014, ArXiv.

[17]  Chunhua Sheng,et al.  A new method towards high-order weno schemes on structured and unstructured grids , 2020 .

[18]  Fermín Navarrina,et al.  A new shock-capturing technique based on Moving Least Squares for higher-order numerical schemes on unstructured grids , 2010 .

[19]  Chao Yan,et al.  Effective high-order energy stable flux reconstruction methods for first-order hyperbolic linear and nonlinear systems , 2020, J. Comput. Phys..

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

[21]  Dinshaw S. Balsara,et al.  An efficient class of WENO schemes with adaptive order for unstructured meshes , 2020, J. Comput. Phys..

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

[23]  H. Ahn.  Hyperbolic cell-centered finite volume method for steady incompressible Navier-Stokes equations on unstructured grids , 2020 .

[24]  S. Osher,et al.  Uniformly high order accuracy essentially non-oscillatory schemes III , 1987 .

[25]  A. Chorin.  A Numerical Method for Solving Incompressible Viscous Flow Problems , 1997 .

[26]  Panagiotis Tsoutsanis,et al.  Stencil selection algorithms for WENO schemes on unstructured meshes , 2019, J. Comput. Phys. X.

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

[28]  R. Mittal.  A Fourier–Chebyshev spectral collocation method for simulating flow past spheres and spheroids , 1999 .

[29]  X. Nogueira,et al.  Moving Kriging reconstruction for high-order finite volume computation of compressible flows , 2013 .

[30]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[31]  D. Hill,et al.  Unstructured-Grid Third-Order Finite Volume Discretization Using a Multistep Quadratic Data-Reconstruction Method , 2010 .

[32]  Panagiotis Tsoutsanis,et al.  WENO schemes on unstructured meshes using a relaxed a posteriori MOOD limiting approach , 2020 .

[33]  Freddie D. Witherden,et al.  Locally adaptive pseudo-time stepping for high-order Flux Reconstruction , 2019, J. Comput. Phys..

[34]  Freddie D. Witherden,et al.  An extended range of stable-symmetric-conservative Flux Reconstruction correction functions , 2015 .

[35]  Jingyang Guo,et al.  A RBF-WENO finite volume method for hyperbolic conservation laws with the monotone polynomial interpolation method , 2017 .

[36]  F. Navarrina,et al.  High‐order finite volume schemes on unstructured grids using moving least‐squares reconstruction. Application to shallow water dynamics , 2006 .

[37]  G. Taylor.  Stability of a Viscous Liquid Contained between Two Rotating Cylinders , 1923 .

[38]  Hiroaki Nishikawa,et al.  A hyperbolic Poisson solver for tetrahedral grids , 2020, J. Comput. Phys..

[39]  Antony Jameson,et al.  Symmetric quadrature rules for simplexes based on sphere close packed lattice arrangements , 2014, J. Comput. Appl. Math..

[40]  Hiroaki Nishikawa,et al.  First order hyperbolic approach for Anisotropic Diffusion equation , 2019, J. Comput. Phys..

[41]  Weiwei Zhang,et al.  A high-order finite volume method on unstructured grids using RBF reconstruction , 2016, Comput. Math. Appl..

[42]  Frank E. Ham,et al.  Symmetric quadrature rules for tetrahedra based on a cubic close-packed lattice arrangement , 2012, J. Comput. Appl. Math..

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

[44]  Fermín Navarrina,et al.  Finite volume solvers and movingleast-squares approximations for thecompressible Navier-Stokes equations onunstructured grids , 2007 .

[45]  Timothy J. Barth,et al.  The design and application of upwind schemes on unstructured meshes , 1989 .

[46]  Hiroaki Nishikawa,et al.  Dimensional scaling and numerical similarity in hyperbolic method for diffusion , 2018, J. Comput. Phys..

[47]  Valerio D’Alessandro,et al.  Assessment of a high-order discontinuous Galerkin method for incompressible three-dimensional Navier–Stokes equations: Benchmark results for the flow past a sphere up to Re = 500 , 2013 .

[48]  Freddie D. Witherden,et al.  A high-order cross-platform incompressible Navier-Stokes solver via artificial compressibility with application to a turbulent jet , 2018, Comput. Phys. Commun..

[49]  O. Botella,et al.  BENCHMARK SPECTRAL RESULTS ON THE LID-DRIVEN CAVITY FLOW , 1998 .

[50]  Hiroaki Nishikawa.  On hyperbolic method for diffusion with discontinuous coefficients , 2018, J. Comput. Phys..

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

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

[53]  William W. Willmarth,et al.  Some experimental results on sphere and disk drag , 1971 .

[54]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[55]  Freddie D. Witherden,et al.  On the utility of GPU accelerated high-order methods for unsteady flow simulations: A comparison with industry-standard tools , 2017, J. Comput. Phys..

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

[57]  ChorinAlexandre Joel.  A Numerical Method for Solving Incompressible Viscous Flow Problems , 1997 .

[58]  Panagiotis Tsoutsanis,et al.  Improvement of the computational performance of a parallel unstructured WENO finite volume CFD code for Implicit Large Eddy Simulation , 2018, Computers & Fluids.

[59]  Brian C. Vermeire,et al.  Optimal embedded pair Runge-Kutta schemes for pseudo-time stepping , 2020, J. Comput. Phys..

[60]  Patrick Knupp,et al.  Code Verification by the Method of Manufactured Solutions , 2000 .