Highly-Accurate Filter-Based Artificial-Dissipation Schemes for Stiff Unsteady Fluid Systems

It is well known that the unmodified application of central difference schemes to the Euler equations produces numerically unstable results because such schemes do not naturally damp the high-frequency modes involved in odd-even decoupling. This shortcoming is usually overcome by adding artificial-dissipation terms, thereby producing stable schemes at the price of potential loss of solution accuracy. For unsteady fluid dynamics, solution filtering schemes have been proposed as a more accurate alternative to artificial dissipation especially when explicit physical-time integration is utilized. However, to solve computationally stiff problems efficiently, it is necessary to use a dual-time stepping approach, to which the application of solution filtering is not straightforward. Restricting the solution filtering only to the physical-time level does not guarantee numerical stability, as errors can accumulate in pseudo-time, causing divergence, while including it at the pseudo-time level introduces inconsistencies that lead to convergence problems. In the present work, Shapiroand Purser-type explicit solution filters are used to derive a new class of filter-equivalent artificial-dissipation operators that can be applied in pseudo time to produce stable, convergent, low-dissipation solutions. These novel artificial-dissipation operators are shown to be indistinguishable from their corresponding filtering procedure for explicit, single-time schemes and are just as effective within a dual-time framework. In addition, the filter-based artificial-dissipation schemes are formulated for use with local preconditioning methods and applied to stiff problems such as low-Mach unsteady flows.

[1]  Thomas H. Pulliam,et al.  Artificial Dissipation Models for the Euler Equations , 1985 .

[2]  A. Prothero,et al.  On the stability and accuracy of one-step methods for solving stiff systems of ordinary differential equations , 1974 .

[3]  U. Shaked,et al.  H∞ nonlinear filtering , 1996 .

[4]  J. Butcher Numerical Methods for Ordinary Differential Equations: Butcher/Numerical Methods , 2005 .

[5]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

[6]  S. Venkateswaran,et al.  Analysis of preconditioning methods for the euler and navier-stokes equations , 1999 .

[7]  Charles L. Merkle,et al.  Comparison of Artificial Dissipation and Filtering Schemes for Time-Accurate Simulations , 2015 .

[8]  Li Wang,et al.  Implicit solution of the unsteady Euler equations for high-order accurate discontinuous Galerkin discretizations , 2006, J. Comput. Phys..

[9]  Ralph Shapiro,et al.  Smoothing, filtering, and boundary effects , 1970 .

[10]  M. Carpenter,et al.  Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations , 2003 .

[11]  Ding Li,et al.  Influence of Large-Scale Pressure Changes on Preconditioned Solutions at Low Speeds , 2004 .

[12]  J. Neumann,et al.  Numerical Integration of the Barotropic Vorticity Equation , 1950 .

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

[14]  S. Venkateswaran,et al.  Efficiency and accuracy issues in contemporary CFD algorithms , 2000 .

[15]  W. Habashi,et al.  2N-Storage Low Dissipation and Dispersion Runge-Kutta Schemes for Computational Acoustics , 1998 .

[16]  V. Vatsa,et al.  Pressure updating methods for the steady-state fluid equations , 1995 .

[17]  M.Y. Hussaini,et al.  Low-Dissipation and Low-Dispersion Runge-Kutta Schemes for Computational Acoustics , 1994 .

[18]  J. Butcher Numerical methods for ordinary differential equations , 2003 .

[19]  R. Purser The Filtering of Meteorological Fields , 1987 .

[20]  Hester Bijl,et al.  Implicit Time Integration Schemes for the Unsteady Compressible Navier–Stokes Equations: Laminar Flow , 2002 .

[21]  Nathan L. Mundis,et al.  The Role of Dispersion and Dissipation on Stabilization Strategies for Time-Accurate Simulations , 2016 .

[22]  C. Merkle,et al.  The application of preconditioning to reacting flow computations , 1995 .

[23]  R. Shapiro The Use of Linear Filtering as a Parameterization of Atmospheric Diffusion , 1971 .

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

[25]  Andrew M. Wissink,et al.  Feature-driven Cartesian adaptive mesh refinement for vortex-dominated flows , 2011, J. Comput. Phys..

[26]  D. Gaitonde,et al.  Pade-Type Higher-Order Boundary Filters for the Navier-Stokes Equations , 2000 .

[27]  E. Turkel,et al.  Preconditioned methods for solving the incompressible low speed compressible equations , 1987 .

[28]  Nathan L. Mundis,et al.  Optimal Runge-Kutta Schemes for High-order Spatial and Temporal Discretizations , 2015 .

[29]  R. D. Richtmyer,et al.  A Method for the Numerical Calculation of Hydrodynamic Shocks , 1950 .