Fifth order finite volume WENO in general orthogonally - curvilinear coordinates

High order reconstruction in the finite volume (FV) approach is achieved by a more fundamental form of the fifth order WENO reconstruction in the framework of orthogonally-curvilinear coordinates, for solving the hyperbolic conservation equations. The derivation employs a piecewise parabolic polynomial approximation to the zone averaged values to reconstruct the right, middle, and left interface values. The grid dependent linear weights of the WENO are recovered by inverting a Vandermode-like linear system of equations with spatially varying coefficients. A scheme for calculating the linear weights, optimal weights, and smoothness indicator on a regularly- and irregularly-spaced grid in orthogonally-curvilinear coordinates is proposed. A grid independent relation for evaluating the smoothness indicator is derived from the basic definition. Finally, the procedures for the source term integration and extension to multi-dimensions are proposed. Analytical values of the linear and optimal weights, and also the weights required for the source term integration and flux averaging, are provided for a regularly-spaced grid in Cartesian, cylindrical, and spherical coordinates. Conventional fifth order WENO reconstruction for the regularly-spaced grids in the Cartesian coordinates can be fully recovered in the case of limiting curvature. The fifth order finite volume WENO-C (orthogonally-curvilinear version of WENO) reconstruction scheme is tested for several 1D and 2D benchmark test cases involving smooth and discontinuous flows in cylindrical and spherical coordinates.

[1]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

[2]  E. Johnsen,et al.  High-order schemes for the Euler equations in cylindrical/spherical coordinates , 2017, 1701.04834.

[3]  Jun Zhu,et al.  A new third order finite volume weighted essentially non-oscillatory scheme on tetrahedral meshes , 2017, J. Comput. Phys..

[4]  Chi-Wang Shu.  High-order Finite Difference and Finite Volume WENO Schemes and Discontinuous Galerkin Methods for CFD , 2003 .

[5]  Dinshaw S. Balsara.  Divergence-free reconstruction of magnetic fields and WENO schemes for magnetohydrodynamics , 2009, J. Comput. Phys..

[6]  Eric Johnsen,et al.  Implementation of WENO schemes in compressible multicomponent flow problems , 2005, J. Comput. Phys..

[7]  Phillip Colella,et al.  A limiter for PPM that preserves accuracy at smooth extrema , 2008, J. Comput. Phys..

[8]  Dinshaw S. Balsara,et al.  An efficient class of WENO schemes with adaptive order , 2016, J. Comput. Phys..

[9]  Eleuterio F. Toro,et al.  WENO schemes based on upwind and centred TVD fluxes , 2005 .

[10]  Andrea Mignone,et al.  High-order conservative reconstruction schemes for finite volume methods in cylindrical and spherical coordinates , 2014, J. Comput. Phys..

[11]  Michael Dumbser,et al.  Efficient implementation of high order unstructured WENO schemes for cavitating flows , 2013 .

[12]  V. Rusanov,et al.  The calculation of the interaction of non-stationary shock waves and obstacles , 1962 .

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

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

[15]  Mengping Zhang,et al.  On the Order of Accuracy and Numerical Performance of Two Classes of Finite Volume WENO Schemes , 2011 .

[16]  J. M. Powers,et al.  Mapped weighted essentially non-oscillatory schemes: Achieving optimal order near critical points , 2005 .

[17]  Chi-Wang Shu,et al.  High Order Weighted Essentially Nonoscillatory Schemes for Convection Dominated Problems , 2009, SIAM Rev..

[18]  Z. Wang,et al.  Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids , 2002 .

[19]  Chi-Wang Shu,et al.  A technique of treating negative weights in WENO schemes , 2000 .

[20]  Phillip Colella,et al.  A HIGH-ORDER FINITE-VOLUME METHOD FOR CONSERVATION LAWS ON LOCALLY REFINED GRIDS , 2011 .

[21]  G. Sod.  A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws , 1978 .

[22]  J. Blondin,et al.  The piecewise-parabolic method in curvilinear coordinates , 1993 .

[23]  Chi-Wang Shu,et al.  High order finite difference and finite volume WENO schemes and discontinuous Galerkin methods for CFD , 2001 .

[24]  Gabriella Puppo,et al.  A third order central WENO scheme for 2D conservation laws , 2000 .

[26]  U. Ziegler,et al.  A semi-discrete central scheme for magnetohydrodynamics on orthogonal-curvilinear grids , 2011, J. Comput. Phys..

[27]  B. Fryxell,et al.  FLASH: An Adaptive Mesh Hydrodynamics Code for Modeling Astrophysical Thermonuclear Flashes , 2000 .

[28]  Kun Xu,et al.  A high-order multidimensional gas-kinetic scheme for hydrodynamic equations , 2013 .

[29]  Nikolaus A. Adams,et al.  A family of high-order targeted ENO schemes for compressible-fluid simulations , 2016, J. Comput. Phys..

[30]  Wai-Sun Don,et al.  An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws , 2008, J. Comput. Phys..

[31]  Harold L. Atkins,et al.  A Finite-Volume High-Order ENO Scheme for Two-Dimensional Hyperbolic Systems , 1993 .

[32]  Eleuterio F. Toro,et al.  Finite-volume WENO schemes for three-dimensional conservation laws , 2004 .

[33]  Luka Sopta,et al.  Extension of ENO and WENO schemes to one-dimensional sediment transport equations , 2004 .

[34]  E. Müller,et al.  A conservative second-order difference scheme for curvilinear coordinates. I: Assignment of variables on a staggered grid , 1989 .

[35]  Nail K. Yamaleev,et al.  A systematic methodology for constructing high-order energy stable WENO schemes , 2009, J. Comput. Phys..

[36]  Gabriella Puppo,et al.  Cool WENO schemes , 2017, Computers & Fluids.

[37]  Christiane Helzel,et al.  Improved Accuracy of High-Order WENO Finite Volume Methods on Cartesian Grids , 2014, Journal of Scientific Computing.

[38]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..

[39]  M. Dumbser,et al.  High-Order Unstructured Lagrangian One-Step WENO Finite Volume Schemes for Non-Conservative Hyperbolic Systems: Applications to Compressible Multi-Phase Flows , 2013, 1304.4816.

[40]  Michael Dumbser,et al.  Efficient, high accuracy ADER-WENO schemes for hydrodynamics and divergence-free magnetohydrodynamics , 2008, Journal of Computational Physics.

[41]  S. Falle.  Self-similar jets , 1991 .

[42]  Manuel J. Castro,et al.  Third‐ and fourth‐order well‐balanced schemes for the shallow water equations based on the CWENO reconstruction , 2018, International Journal for Numerical Methods in Fluids.

[43]  Kun Xu,et al.  A gas-kinetic BGK scheme for the Navier-Stokes equations and its connection with artificial dissipation and Godunov method , 2001 .

[44]  Ahmadreza Pishevar,et al.  A New fourth order central WENO method for 3D hyperbolic conservation laws , 2012, Appl. Math. Comput..

[45]  G. Russo,et al.  Central WENO schemes for hyperbolic systems of conservation laws , 1999 .

[46]  Zhi J. Wang,et al.  Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids. Basic Formulation , 2002 .

[47]  Jun Zhu,et al.  New Finite Volume Weighted Essentially Nonoscillatory Schemes on Triangular Meshes , 2018, SIAM J. Sci. Comput..

[48]  Guang-Shan Jiang,et al.  A High-Order WENO Finite Difference Scheme for the Equations of Ideal Magnetohydrodynamics , 1999 .

[49]  Wai-Sun Don,et al.  Accuracy of the weighted essentially non-oscillatory conservative finite difference schemes , 2013, J. Comput. Phys..

[50]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[51]  M. Semplice,et al.  Adaptive Mesh Refinement for Hyperbolic Systems Based on Third-Order Compact WENO Reconstruction , 2014, Journal of Scientific Computing.

[52]  Gabriella Puppo,et al.  CWENO: Uniformly accurate reconstructions for balance laws , 2016, Math. Comput..

[53]  P. Lax,et al.  On Upstream Differencing and Godunov-Type Schemes for Hyperbolic Conservation Laws , 1983 .

[54]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

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

[56]  R. W. MacCormack,et al.  A Numerical Method for Solving the Equations of Compressible Viscous Flow , 1981 .

[57]  Chi-Wang Shu,et al.  Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy , 2000 .

[58]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[59]  Xu-Dong Liu,et al.  Solution of Two-Dimensional Riemann Problems of Gas Dynamics by Positive Schemes , 1998, SIAM J. Sci. Comput..

[60]  Xiangmin Jiao,et al.  WLS-ENO: Weighted-least-squares based essentially non-oscillatory schemes for finite volume methods on unstructured meshes , 2016, J. Comput. Phys..