Construction of Explicit and Implicit Symmetric TVD Schemes and Their Applications

A one-parameter family of second-order explicit and implicit total variation diminishing (TVD) schemes is reformulated so that a simplier and wider group of limiters is included. The resulting scheme can be viewed as a symmetrical algorithm with a variety of numerical dissipation terms that are designed for weak solutions of hyperbolic problems. This is a generalization of recent works of Roe and Davis to a wider class of symmetric schemes other than Lax-Wendroff. The main properties of the present class of schemes are that they can be implicit, and, when steady-state calculations are sought, the numerical solution is independent of the time step. Numerical experiments with two-dimensional unsteady and steady-state airfoil calculations show that the proposed symmetric TVD schemes are quite robust and accurate

Construction of Explicit and Implicit Symmetric TVD Schemes and Their Applications

A one-parameter family of second-order explicit and implicit total variation diminishing (TVD) schemes is reformulated so that a simplier and wider group of limiters is included. The resulting scheme can be viewed as a symmetrical algorithm with a variety of numerical dissipation terms that are designed for weak solutions of hyperbolic problems. This is a generalization of recent works of Roe and Davis to a wider class of symmetric schemes other than Lax-Wendroff. The main properties of the present class of schemes are that they can be implicit, and, when steady-state calculations are sought, the numerical solution is independent of the time step. Numerical experiments with two-dimensional unsteady and steady-state airfoil calculations show that the proposed symmetric TVD schemes are quite robust and accurate

Non-Linear Filtering and Limiting in High Order Methods for Ideal and Non-Ideal MHD

The adaptive nonlinear filtering and limiting in spatially high order schemes (Yee et al. J. Comput. Phys. 150, 199–238, (1999), Sjogreen and Yee, J. Scient. Comput. 20, 211–255, (2004)) for the compressible Euler and Navier–Stokes equations have been recently extended to the ideal and non-ideal magnetohydrodynamics (MHD) equations, (Sjogreen and Yee, (2003), Proceedings of the 16th AIAA/CFD conference, June 23–26, Orlando F1; Yee and Sjogreen (2003), Proceedings of the International Conference on High Performance Scientific Computing, March, 10–14, Honai, Vietnam; Yee and Sjogreen (2003), RIACS Technical Report TR03. 10, July, NASA Ames Research Center; Yee and Sjogreen (2004), Proceedings of the ICCF03, July 12–16, Toronto, Canada). The numerical dissipation control in these adaptive filter schemes consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free from numerical dissipation contamination. The numerical dissipation considered consists of high order linear dissipation for the suppression of high frequency oscillation and the nonlinear dissipative portion of high-resolution shock-capturing methods for discontinuity capturing. The applicable nonlinear dissipative portion of high-resolution shock-capturing methods is very general. The objective of this paper is to investigate the performance of three commonly used types of discontinuity capturing nonlinear numerical dissipation for both the ideal and non-ideal MHD

Grid convergence of high order methods for multiscale complex unsteady viscous compressible flows

Grid convergence of several high order methods for the computation of rapidly developing complex unsteady viscous compressible flows with a wide range of physical scales is studied. The recently developed adaptive numerical dissipation control high order methods referred to as the ACM and wavelet filter schemes are compared with a fifth-order weighted ENO (WENO) scheme. The two 2-D compressible full Navier–Stokes models considered do not possess known analytical and experimental data. Fine grid solutions from a standard second-order TVD scheme and a MUSCL scheme with limiters are used as reference solutions. The first model is a 2-D viscous analog of a shock tube problem which involves complex shock/shear/boundary-layer interactions. The second model is a supersonic reactive flow concerning fuel breakup. The fuel mixing involves circular hydrogen bubbles in air interacting with a planar moving shock wave. Both models contain fine scale structures and are stiff in the sense that even though the unsteadiness of the flows are rapidly developing, extreme grid refinement and time step restrictions are needed to resolve all the flow scales as well as the chemical reaction scales. Our computations were all made on uniform grids, and our conclusions cannot be directly carried over to, for example, curvilinear grids

Simulation of Richtmyer–Meshkov instability by sixth-order filter methods

Simulation of a 2-D Richtmyer–Meshkov instability (RMI), including inviscid, viscous and magnetic field effects was conducted comparing recently developed sixthorder filter schemes with various standard shock-capturing methods. The suppression of the inviscid gas dynamics RMI in the presence of a magnetic field was investigated by Samtaney and Wheatley et al. Numerical results illustrated here exhibit behavior similar to the work of Samtaney. Due to the different amounts and different types of numerical dissipation contained in each scheme, the structures and the growth of eddies for the chaotic-like inviscid gas dynamics RMI case are highly grid size and scheme dependent, even with many levels of refinement. The failure of grid refinement for all studied numerical methods extends to the viscous gas dynamics case for high Reynolds number. For lower Reynolds number, grid convergence has been achieved by all studied methods. To achieve similar resolution, standard shock-capturing methods require more grid points than filter schemes and yet the CPU times using the same grid for all studied methods are comparable

Accuracy consideration by DRP schemes for DNS and LES of compressible flow computations

Several dispersion relation-preserving (DRP) spatially central discretizations are considered as the base scheme in the framework of the Yee & Sjögreen low dissipative nonlinear filter approach. In addition, the nonlinear filter of Yee & Sjögreen with shock-capturing and long time integration capabilities is used to replace the standard DRP linear filter for both smooth flows and flows containing discontinuities. DRP schemes for computational aeroacoustics (CAA) focus on dispersion error consideration for long time lin- ear wave propagation rather than the formal order of accuracy of the scheme. The resulting DRP schemes usually have wider grid stencils and increased CPU operations count compared with standard central schemes of the same formal order of accuracy. For discontinuous initial data and long time wave propa- gation of smooth acoustic waves, various space and time DRP linear filter are needed. For acoustic waves interacting with shocks and turbulence induced noise, DRP schemes with linear filters alone usually are not capable of simulating such flows. The investigation presented in this paper is focused on the pos- sible gain in efficiency and accuracy by spatial DRP schemes over standard central schemes having the same grid stencil width for general direct numerical simulations (DNS) and large eddy simulations (LES) of compressible flows. Representative test cases for both smooth flows and problems containing discontinuities for 3D DNS of compressible gas dynamics are included

High Order Filter Methods for Wide Range of Compressible Flow Speeds

This paper extends the accuracy of the high order nonlinear filter finite difference method of Yee and Sjogreen [Development of Low Dissipative High Order Filter Schemes for Multiscale Navier-Stokes/MHD Systems, J. Comput. Phys., 225 (2007) 910–934] and Sjogreen and Yee [Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for Shock-Turbulence Computation, RIACS Technical Report TR01.01, NASA Ames research center (Oct 2000); Also J. Scient. Comput., 20 (2004) 211–255] for compressible turbulence with strong shocks to a wider range of flow speeds without having to tune the key filter parameter. Such a filter method consists of two steps: a full time step using a spatially high-order non-dissipative base scheme, followed by a post-processing filter step. The postprocessing filter step consists of the products of wavelet-based flow sensors and nonlinear numerical dissipations. For low speed turbulent flows and long time integration of smooth flows, the existing flow sensor relies on tuning the amount of shock-dissipation in order to obtain highly accurate turbulent numerical solutions. The improvement proposed here is to solve the conservative skew-symmetric form of the governing equations in conjunction with an added flow speed and shock strength indicator to minimize the tuning of the key filter parameter. Test cases illustrate the improved accuracy by the proposed ideas without tuning the key filter parameter of the nonlinear filter step

Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for High Order Methods

The recently developed essentially fourth-order or higher low dissipative shockcapturing scheme of Yee, Sandham, and Djomehri [25] aimed at minimizing numerical dissipations for high speed compressible viscous flows containing shocks, shears and turbulence. To detect non-smooth behavior and control the amount of numerical dissipation to be added, Yee et al. employed an artificial compression method (ACM) of Harten [4] but utilize it in an entirely different context than Harten originally intended. The ACM sensor consists of two tuning parameters and is highly physical problem dependent. To minimize the tuning of parameters and physical problem dependence, new sensors with improved detection properties are proposed. The new sensors are derived from utilizing appropriate non-orthogonal wavelet basis functions and they can be used to completely switch off the extra numerical dissipation outside shock layers. The non-dissipative spatial base scheme of arbitrarily high order of accuracy can be maintained without compromising its stability at all parts of the domain where the solution is smooth. Two types of redundant non-orthogonal wavelet basis functions are considered. One is the B-spline wavelet (Mallat and Zhong [14]) used by Gerritsen and Olsson [3] in an adaptive mesh refinement method, to determine regions where refinement should be done. The other is the modification of the multiresolution method of Harten [5] by converting it to a new, redundant, non-orthogonal wavelet. The wavelet sensor is then obtained by computing the estimated Lipschitz exponent of a chosen physical quantity (or vector) to be sensed on a chosen wavelet basis function. Both wavelet sensors can be viewed as dual purpose adaptive methods leading to dynamic numerical dissipation control and improved grid adaptation indicators. Consequently, they are useful not only for shock-turbulence computations but also for computational aeroacoustics and numerical combustion. In addition, these sensors are scheme independent and can be stand-alone options for numerical algorithms other than the Yee et al. scheme

High order numerical simulation of sound generated by the Kirchhoff vortex

An improved high order finite difference method for low Mach number computational aeroacoustics (CAA) is described. The improvements involve the conditioning of the Euler equations to minimize numerical cancellation errors, and the use of a stable non-dissipative sixth-order central spatial interior scheme and a third-order boundary scheme. Spurious high frequency oscillations are damped by a third-order characteristic-based filter. The objective of this paper is to apply these improvements in the simulation of sound generated by the Kirchhoff vortex

Development of low dissipative high order filter schemes for multiscale Navier–Stokes/MHD systems

Recent progress in the development of a class of low dissipative high order (fourth-order or higher) filter schemes for multiscale Navier–Stokes, and ideal and non-ideal magnetohydrodynamics (MHD) systems is described. The four main features of this class of schemes are: (a) multiresolution wavelet decomposition of the computed flow data as sensors for adaptive numerical dissipative control, (b) multistep filter to accommodate efficient application of different numerical dissipation models and different spatial high order base schemes, (c) a unique idea in solving the ideal conservative MHD system (a non-strictly hyperbolic conservation law) without having to deal with an incomplete eigensystem set while at the same time ensuring that correct shock speeds and locations are computed, and (d) minimization of the divergence of the magnetic field numerical error. By design, the flow sensors, different choice of high order base schemes and numerical dissipation models are stand-alone modules. A whole class of low dissipative high order schemes can be derived at ease, making the resulting computer software very flexible with widely applicable. Performance of multiscale and multiphysics test cases are illustrated with many levels of grid refinement and comparison with commonly used schemes in the literature