On central-difference and upwind schemes

A class of numerical dissipation models for central-difference schemes constructed with second- and fourth-difference terms is considered. The notion of matrix dissipation associated with upwind schemes is used to establish improved shock capturing capability for these models. In addition, conditions are given that guarantee that such dissipation models produce a Total Variation Diminishing (TVD) scheme. Appropriate switches for this type of model to ensure satisfaction of the TVD property are presented. Significant improvements in the accuracy of a central-difference scheme are demonstrated by computing both inviscid and viscous transonic airfoil flows

Multigrid for hypersonic viscous two- and three-dimensional flows

The use of a multigrid method with central differencing to solve the Navier-Stokes equations for hypersonic flows is considered. The time dependent form of the equations is integrated with an explicit Runge-Kutta scheme accelerated by local time stepping and implicit residual smoothing. Variable coefficients are developed for the implicit process that removes the diffusion limit on the time step, producing significant improvement in convergence. A numerical dissipation formulation that provides good shock capturing capability for hypersonic flows is presented. This formulation is shown to be a crucial aspect of the multigrid method. Solutions are given for two-dimensional viscous flow over a NACA 0012 airfoil and three-dimensional flow over a blunt biconic

Hypersonic blunt body computations including real gas effects

The recently developed second-order explicit and implicit total variation diminishing (TVD) shock-capturing methods of the Harten and Yee, Yee, and van Leer types in conjunction with a generalized Roe's approximate Riemann solver of Vinokur and the generalized flux-vector splittings of Vinokur and Montagne for two-dimensional hypersonic real gas flows are studied. A previous study on one-dimensional unsteady problems indicated that these schemes produce good shock-capturing capability and that the state equation does not have a large effect on the general behavior of these methods for a wide range of flow conditions for equilibrium air. The objective of this paper is to investigate the applicability and shock resolution of these schemes for two-dimensional steady-state hypersonic blunt body flows. The main contribution of this paper is to identify some of the elements and parameters which can affect the convergence rate for high Mach numbers or real gases but have negligible effect for low Mach number cases for steady-state inviscid blunt body flows

An Adaptive Characteristic-wise Reconstruction WENOZ scheme for Gas Dynamic Euler Equations

Due to its excellent shock-capturing capability and high resolution, the WENO scheme family has been widely used in varieties of compressive flow simulation. However, for problems containing strong shocks and contact discontinuities, such as the Lax shock tube problem, the WENO scheme still produces numerical oscillations. To avoid such numerical oscillations, the characteristic-wise construction method should be applied. Compared to component-wise reconstruction, characteristic-wise reconstruction leads to much more computational cost and thus is not suite for large scale simulation such as direct numeric simulation of turbulence. In this paper, an adaptive characteristic-wise reconstruction WENO scheme, i.e. the AdaWENO scheme, is proposed to improve the computational efficiency of the characteristic-wise reconstruction method. The new scheme performs characteristic-wise reconstruction near discontinuities while switching to component-wise reconstruction for smooth regions. Meanwhile, a new calculation strategy for the WENO smoothness indicators is implemented to reduce over-all computational cost. Several one dimensional and two dimensional numerical tests are performed to validate and evaluate the AdaWENO scheme. Numerical results show that AdaWENO maintains essentially non-oscillatory flow field near discontinuities as the characteristic-wise reconstruction method. Besieds, compared to the component-wise reconstruction, AdaWENO is about 40\% faster which indicates its excellent efficiency

Nonlinear spectral-like schemes for hybrid schemes

In spectral-like resolution-WENO hybrid schemes, if the switch function takes more grid points as discontinuity points, the WENO scheme is often turned on, and the numerical solutions may be too dissipative. Conversely, if the switch function takes less grid points as discontinuity points, the hybrid schemes usually are found to produce oscillatory solutions or just to be unstable. Even if the switch function takes less grid points as discontinuity points, the final hybrid scheme is inclined to be more stable, provided the spectral-like resolution scheme in the hybrid scheme has moderate shock-capturing capability. Following this idea, we propose nonlinear spectral-like schemes named weighted group velocity control (WGVC) schemes. These schemes show not only high-resolution for short waves but also moderate shock capturing capability. Then a new class of hybrid schemes is designed in which the WGVC scheme is used in smooth regions and the WENO scheme is used to capture discontinuities. These hybrid schemes show good resolution for small-scales structures and fine shock-capturing capabilities while the switch function takes less grid points as discontinuity points. The seven-order WGVC-WENO scheme has also been applied successfully to the direct numerical simulation of oblique shock wave-turbulent boundary layer interaction

A time-accurate high-resolution TVD scheme for solving the Navier-Stokes equations

A total variation diminishing (TVD) scheme has been developed and incorporated into an existing time-accurate high-resolution Navier-Stokes code. The accuracy and the robustness of the resulting solution procedure have been assessed by performing many calculations in four different areas: shock tube flows, regular shock reflection, supersonic boundary layer, and shock boundary layer interactions. These numerical results compare well with corresponding exact solutions or experimental data

Computational analysis of hypersonic airbreathing aircraft flow fields

The general problem of calculating the flow fields associated with hypersonic airbreathing aircraft is presented. Unique aspects of hypersonic aircraft aerodynamics are introduced and their demands on computational fluid dynamics are outlined. Example calculations associated with inlet/forebody integration and hypersonic nozzle design are presented to illustrate the nature of the problems considered