Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shock-capturing schemes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

114 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

On various treatments of potential equations at shocks

The potential equation that governs transonic inviscid flowfields is elliptic in subsonic regions and hyperbolic in supersonic regions. The transition from supersonic to subsonic flow may occur either continuously or across a surface of discontinuity known as a shock. Because of the changing type of governing equation and the presence of the surfaces of discontinuity, special treatments of the governing equations are needed at shocks to preserve mass conservation and to ensure satisfaction of the governing equation. In the so-called shock-capturing techniques, shocks are captured by smearing the discontinuity over several mesh cells in the computational domain. Most existing conservative schemes conserve mass flux throughout the flowfield, including shock boundaries; however, they also introduce zero-order errors in the approximation to the governing equation. It is shown that the zero-order errors, which do not diminish as the mesh spacing approaches zero, can cause discrepancies in the prediction of shock strength and location. These parameters are extremely sensitive features of transonic inviscid flowfield calculations. It is also shown that the zero-order errors can be avoided by developing higher-order schemes that properly model the governing equation at shocks.

Chen, L. T.↗
NASA NTRS1989

A class of high resolution explicit and implicit shock-capturing methods

An attempt is made to give a unified and generalized formulation of a class of high resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock wave computations. Included is a systematic review of the basic design principle of the various related numerical methods. Special emphasis is on the construction of the basis nonlinear, spatially second and third order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and the flux vector splitting approaches. Generalization of these methods to efficiently include equilibrium real gases and large systems of nonequilibrium flows are discussed. Some issues concerning the applicability of these methods that were designed for homogeneous hyperbolic conservation laws to problems containing stiff source terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for 1-, 2- and 3-dimensional gas dynamics problems.

Yee, H. C.↗
NASA NTRS1992

Numerical simulation of entry flow over blunt swept-wing planes

Viscous, heat-conducting flow with chemical and vibrational relaxation processes of the constituent gases surrounding winged spacecraft is considered in the continuum regime. The Navier-Stokes equations are appended by additional vibrational energy and species rate equations and supplemented by the equations of state and the phenomenological laws based on mixture rules or collisional cross sections. Numerical convective flux can be obtained from several forms of one-dimensional Riemann solver, with or without entropy correction. High-order accuracy is obtained from two types of reconstructive interpolation. A number of explicit and implicit numerical schemes have been implemented as a means to yield converged solutions. Both shock-fitting, finite-difference and shock-capturing, finite-volume techniques have been tested for configurations such as a sphere, double ellipsoid, blunt-edge delta wing, a European Hermes vehicle, and the U.S. Shuttle Orbiter. The shock-fitting code provides excellent results only for simple configurations, whereas the shock-capturing code leads to overall satisfying solutions for complex geometries.

Li, C. P.↗
NASA NTRS1998

Low Dissipative High Order Shock-Capturing Methods using Characteristic-Based Filters

An approach which closely maintains the non-dissipative nature of classical fourth or higher- order spatial differencing away from shock waves and steep gradient regions while being capable of accurately capturing discontinuities, steep gradient and fine scale turbulent structures in a stable and efficient manner is described. The approach is a generalization of the method of Gustafsson and Olsson and the artificial compression method (ACM) of Harten. Spatially non-dissipative fourth or higher-order compact and non-compact spatial differencings are used as the base schemes. Instead of applying a scalar filter as in Gustafsson and Olsson, an ACM like term is used to signal the appropriate amount of second or third-order TVD or ENO types of characteristic based numerical dissipation. This term acts as a characteristic filter to minimize numerical dissipation for the overall scheme. For time-accurate computations, time discretizations with low dissipation are used. Numerical experiments on 2-D vortical flows, vortex-shock interactions and compressible spatially and temporally evolving mixing layers showed that the proposed schemes have the desired property with only a 10% increase in operations count over standard second-order TVD schemes. Aside from the ability to accurately capture shock-turbulence interaction flows, this approach is also capable of accurately preserving vortex convection. Higher accuracy is achieved with fewer grid points when compared to that of standard second-order TVD or ENO schemes. To demonstrate the applicability of these schemes in sustaining turbulence where shock waves are absent, a simulation of 3-D compressible turbulent channel flow in a small domain is conducted.

Yee, H. C.↗
NASA NTRS1998

Low Dissipative High Order Shock-Capturing Methods Using Characteristic-Based Filters

An approach which closely maintains the non-dissipative nature of classical fourth or higher- order spatial differencing away from shock waves and steep gradient regions while being capable of accurately capturing discontinuities, steep gradient and fine scale turbulent structures in a stable and efficient manner is described. The approach is a generalization of the method of Gustafsson and Oisson and the artificial compression method (ACM) of Harten. Spatially non-dissipative fourth or higher-order compact and non-compact spatial differencings are used as the base schemes. Instead of applying a scalar filter as in Gustafsson and Olsson, an ACM like term is used to signal the appropriate amount of second or third-order TVD or ENO types of characteristic based numerical dissipation. This term acts as a characteristic filter to minimize numerical dissipation for the overall scheme. For time-accurate computations, time discretizations with low dissipation are used. Numerical experiments on 2-D vortical flows, vortex-shock interactions and compressible spatially and temporally evolving mixing layers showed that the proposed schemes have the desired property with only a 10% increase in operations count over standard second-order TVD schemes. Aside from the ability to accurately capture shock-turbulence interaction flows, this approach is also capable of accurately preserving vortex convection. Higher accuracy is achieved with fewer grid points when compared to that of standard second-order TVD or ENO schemes. To demonstrate the applicability of these schemes in sustaining turbulence where shock waves are absent, a simulation of 3-D compressible turbulent channel flow in a small domain is conducted.

Yee, H. C.↗
NASA NTRSDate not supplied

An Improved Approach to the Predictability & Reliability of the Onset of Turbulence With Shocks

The construction of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and for (b) obtaining correct propagation speed of discontinuities in the presence of stiff source terms share one important ingredient – minimization of numerical dissipation while maintaining numerical stability. The dual requirements to achieve both numerical stability and minimal numerical dissipation are often conflicting since existing shock capturing schemes were designed mainly to be robust for rapidly developed turbulence-free flows and for shock waves without stiff source term. For the past two decades, Yee and collaborators have focused on an improved understanding of the nonlinear behavior of different high order shock-capturing methods. It was found that even very high order methods without proper nonlinear stability and numerical dissipation control can either numerically smear the onset of turbulence due to excess numerical dissipation, or induce (onset) numerical turbulence that is not physical turbulence due to lack of proper numerical dissipation to improve nonlinear stability for long time integration. Our approach is to combine (I) and (II) below for obtaining the physically correct onset of turbulence with shocks, including problems with stiff source terms: (I) Nonlinear dynamics is utilized to complement the traditional linearized stability theory (Yee & Sweby, Yee et al., Griffiths et al., Lafon & Yee, Yee, Wang et al., Kotov et al. 1990- 2015) in order to (i) Minimize numerically induced false transition to turbulence, (ii) Minimize numerical instability due to long time integration of turbulent flows, (iii) Minimize numerically induced standing wave solutions, and (iv) Minimize wrong propagation of speed of discontinuities due to the presence of stiff source terms. (II) Our recently developed physical preserving (structural preserving) high order methods with improved nonlinear stability & accuracy that are essential in minimizing spurious numerics are used.

HECC↗