A two-stage fourth order time-accurate discretization for Lax-Wendroff type flow solvers II. High order numerical boundary conditions

Abstract This paper serves to treat boundary conditions numerically with high order accuracy in order to suit the two-stage fourth-order finite volume schemes for hyperbolic problems developed in J. Li and Z. Du (2016) [17] . As such, it is significant when capturing small scale structures near physical boundaries. Different from previous contributions in literature, the current approach constructs a fourth order accurate approximation to boundary conditions by only using the Jacobian matrix of the flux function (characteristic information) instead of its successive differentiation of governing equations leading to tensors of high ranks in the inverse Lax–Wendroff method. Technically, data in several ghost cells are constructed with interpolation so that the interior scheme can be implemented over boundary cells, and theoretical boundary condition has to be modified properly at intermediate stages so as to make the two-stage scheme over boundary cells fully consistent with that over interior cells. This is nonintuitive and highlights the fact that theoretical boundary conditions are only prescribed for continuous partial differential equations (PDEs), while they must be approximated in a consistent way (even though they could be exactly valued) when the PDEs are discretized . Several numerical examples are provided to illustrate the performance of the current approach when dealing with general boundary conditions.

[1]  Randall J. LeVeque,et al.  A High-Resolution Rotated Grid Method for Conservation Laws with Embedded Geometries , 2005, SIAM J. Sci. Comput..

[2]  Heinz-Otto Kreiss,et al.  Difference Approximations of the Neumann Problem for the Second Order Wave Equation , 2004, SIAM J. Numer. Anal..

[3]  Lilia Krivodonova,et al.  High-order accurate implementation of solid wall boundary conditions in curved geometries , 2006 .

[4]  B. Sjögreen,et al.  A Cartesian embedded boundary method for hyperbolic conservation laws , 2006 .

[5]  Shiyi Chen,et al.  Experimental investigation of freely falling thin disks. Part 2. Transition of three-dimensional motion from zigzag to spiral , 2013, Journal of Fluid Mechanics.

[6]  Joseph Falcovitz,et al.  Generalized Riemann Problems in Computational Fluid Dynamics , 2003 .

[7]  Chi-Wang Shu,et al.  Inverse Lax-Wendroff procedure for numerical boundary conditions of conservation laws , 2010, J. Comput. Phys..

[8]  Randall J. LeVeque,et al.  H-Box Methods for the Approximation of Hyperbolic Conservation Laws on Irregular Grids , 2003, SIAM J. Numer. Anal..

[9]  David I. Gottlieb,et al.  The Theoretical Accuracy of Runge-Kutta Time Discretizations for the Initial Boundary Value Problem: A Study of the Boundary Error , 1995, SIAM J. Sci. Comput..

[10]  Cheng Wang,et al.  Efficient implementation of high order inverse Lax-Wendroff boundary treatment for conservation laws , 2012, J. Comput. Phys..

[11]  H. Kreiss,et al.  Initial-Boundary Value Problems and the Navier-Stokes Equations , 2004 .

[12]  Jiequan Li,et al.  A Hermite WENO reconstruction for fourth order temporal accurate schemes based on the GRP solver for hyperbolic conservation laws , 2017, J. Comput. Phys..

[13]  Xi Chen,et al.  Interactions between second mode and low-frequency waves in a hypersonic boundary layer , 2017, Journal of Fluid Mechanics.

[14]  Kun Xu,et al.  An efficient and accurate two-stage fourth-order gas-kinetic scheme for the Euler and Navier-Stokes equations , 2016, J. Comput. Phys..

[15]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[16]  Jiequan Li,et al.  Hyperbolic balance laws: Riemann invariants and the generalized Riemann problem , 2007, Numerische Mathematik.

[17]  R. Lindsay,et al.  Elements of gasdynamics , 1957 .

[18]  Shiyi Chen,et al.  Experimental investigation of freely falling thin disks. Part 1. The flow structures and Reynolds number effects on the zigzag motion , 2013, Journal of Fluid Mechanics.

[19]  Rolf Jeltsch,et al.  A Higher-Order Boundary Treatment for Cartesian-Grid Methods , 1998 .

[20]  Chi-Wang Shu,et al.  Stability analysis of the inverse Lax-Wendroff boundary treatment for high order upwind-biased finite difference schemes , 2016, J. Comput. Appl. Math..

[21]  D. Pathria,et al.  The Correct Formulation of Intermediate Boundary Conditions for Runge-Kutta Time Integration of Initial Boundary Value Problems , 1997, SIAM J. Sci. Comput..

[22]  Jiequan Li,et al.  A Two-Stage Fourth Order Time-Accurate Discretization for Lax-Wendroff Type Flow Solvers I. Hyperbolic Conservation Laws , 2015, SIAM J. Sci. Comput..

[23]  HEINZ-OTTO KREISS,et al.  A Second Order Accurate Embedded Boundary Method for the Wave Equation with Dirichlet Data , 2005, SIAM J. Sci. Comput..

[24]  J. Falcovitz,et al.  A Two-Dimensional Conservation Laws Scheme for Compressible Flows with Moving Boundaries , 1997 .

[25]  Ling,et al.  NUMERICAL BOUNDARY CONDITIONS FOR THE FAST SWEEPING HIGH ORDER WENO METHODS FOR SOLVING THE EIKONAL EQUATION , 2008 .

[26]  Gerald Warnecke,et al.  A direct Eulerian GRP scheme for compressible fluid flows , 2006, J. Comput. Phys..