High-order flux correction/finite difference schemes for strand grids

A novel high-order method combining unstructured flux correction along body surfaces and high-order finite differences normal to surfaces is formulated for unsteady viscous flows on strand grids. The flux correction algorithm is applied in each unstructured layer of the strand grid, and the layers are then coupled together via a source term containing derivatives in the strand direction. Strand-direction derivatives are approximated to high-order via summation-by-parts operators for first derivatives and second derivatives with variable coefficients. We show how this procedure allows for the proper truncation error canceling properties required for the flux correction scheme. The resulting scheme possesses third-order design accuracy, but often exhibits fourth-order accuracy when higher-order derivatives are employed in the strand direction, especially for highly viscous flows. We prove discrete conservation for the new scheme and time stability in the absence of the flux correction terms. Results in two dimensions are presented that demonstrate improvements in accuracy with minimal computational and algorithmic overhead over traditional second-order algorithms.

[1]  A. Jameson,et al.  Finite volume solution of the two-dimensional Euler equations on a regular triangular mesh , 1985 .

[2]  A. Wissink,et al.  Strand Grid Solution Procedures for Sharp Corners , 2014 .

[3]  Andrew M. Wissink,et al.  Application of strand meshes to complex aerodynamic flow fields , 2010, J. Comput. Phys..

[4]  Andrew M. Wissink,et al.  Validation of the Strand Grid Approach , 2009 .

[5]  Marcel Vinokur,et al.  Spectral difference method for unstructured grids I: Basic formulation , 2006, J. Comput. Phys..

[6]  Erik Schnetter,et al.  Optimized High-Order Derivative and Dissipation Operators Satisfying Summation by Parts, and Applications in Three-dimensional Multi-block Evolutions , 2005, J. Sci. Comput..

[7]  P. Frederickson,et al.  Higher order solution of the Euler equations on unstructured grids using quadratic reconstruction , 1990 .

[8]  Magnus Svärd,et al.  A stable high-order finite difference scheme for the compressible Navier-Stokes equations, far-field boundary conditions , 2007, J. Comput. Phys..

[9]  P. Roache.  Code Verification by the Method of Manufactured Solutions , 2002 .

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

[11]  HarrisRob,et al.  Efficient quadrature-free high-order spectral volume method on unstructured grids , 2008 .

[12]  H. Kreiss,et al.  Finite Element and Finite Difference Methods for Hyperbolic Partial Differential Equations , 1974 .

[13]  C. Norberg.  Fluctuating lift on a circular cylinder: review and new measurements , 2003 .

[14]  D. Mavriplis.  Multigrid solution of the two-dimensional Euler equations on unstructured triangular meshes , 1987 .

[15]  D. Gottlieb,et al.  The Stability of Numerical Boundary Treatments for Compact High-Order Finite-Difference Schemes , 1993 .

[16]  Antony Jameson,et al.  Advances in Bringing High-order Methods to Practical Applications in Computational Fluid Dynamics , 2011 .

[17]  Andrew M. Wissink,et al.  Application of Strand Meshes to Complex Aerodynamic Flow Fields , 2010 .

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

[19]  Magnus Svärd,et al.  A stable high-order finite difference scheme for the compressible Navier-Stokes equations: No-slip wall boundary conditions , 2008, J. Comput. Phys..

[20]  J. L. Steger,et al.  A chimera grid scheme , 2011 .

[21]  Dimitri J. Mavriplis,et al.  AGGLOMERATION MULTIGRID FOR THE THREE-DIMENSIONAL EULER EQUATIONS , 1994 .

[22]  Magnus Svärd,et al.  Well-Posed Boundary Conditions for the Navier-Stokes Equations , 2005, SIAM J. Numer. Anal..

[23]  D. Brandt,et al.  Multi-level adaptive solutions to boundary-value problems math comptr , 1977 .

[24]  Jan Nordström,et al.  High-order finite difference methods, multidimensional linear problems, and curvilinear coordinates , 2001 .

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

[26]  Magnus Svärd,et al.  On Coordinate Transformations for Summation-by-Parts Operators , 2004, J. Sci. Comput..

[27]  W. K. Anderson,et al.  Petrov-Galerkin and discontinuous-Galerkin methods for time-domain and frequency-domain electromagnetic simulations , 2011, J. Comput. Phys..

[28]  James D. Baeder,et al.  Implicit Hole Cutting - A New Approach to Overset Grid Connectivity , 2003 .

[29]  Aaron Katz,et al.  An Efficient Correction Method to Obtain a Formally Third-Order Accurate Flow Solver for Node-Centered Unstructured Grids , 2012, J. Sci. Comput..

[30]  Zhi Jian Wang,et al.  Efficient quadrature-free high-order spectral volume method on unstructured grids: Theory and 2D implementation , 2008, J. Comput. Phys..

[31]  J. Benek,et al.  A flexible grid embedding technique with application to the Euler equations , 1983 .

[32]  Dimitri J. Mavriplis,et al.  Revisiting the Least-squares Procedure for Gradient Reconstruction on Unstructured Meshes , 2003 .

[33]  Jan Nordström,et al.  Finite volume methods, unstructured meshes and strict stability for hyperbolic problems , 2003 .

[34]  Ken Mattsson,et al.  Summation by Parts Operators for Finite Difference Approximations of Second-Derivatives with Variable Coefficients , 2012, J. Sci. Comput..

[35]  Andrew M. Wissink,et al.  Parallel domain connectivity algorithm for unsteady flow computations using overlapping and adaptive grids , 2010, J. Comput. Phys..

[36]  D. Gottlieb,et al.  Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: methodology and application to high-order compact schemes , 1994 .

[37]  A. Katz,et al.  High-Order Flux Correction for Viscous Flows on Arbitrary Unstructured Grids , 2013 .

[38]  Aaron Katz,et al.  Strand Grid Solution Procedures for Sharp Corners , 2013 .

[39]  Yen Liu,et al.  Quadratic reconstruction finite volume schemes on 3D arbitrary unstructured polyhedral grids , 1999 .

[40]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

[41]  D. Kopriva.  A Conservative Staggered-Grid Chebyshev Multidomain Method for Compressible Flows. II. A Semi-Structured Method , 1996 .

[42]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[43]  Andrew M. Wissink,et al.  On Strand Grids for Complex Flows , 2007 .

[44]  Amir Nejat,et al.  On Obtaining High-Order Finite-Volume Solutions to the Euler Equations on Unstructured Meshes , 2007 .