SEARCH · Engineering Papers

Results for “shock capturing”

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.

SourceResearch connectionFrom yearThrough year

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

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Large perturbation flow field analysis and simulation for supersonic inlets

An analysis technique for simulation of supersonic mixed compression inlets with large flow field perturbations is presented. The approach is based upon a quasi-one-dimensional inviscid unsteady formulation which includes engineering models of unstart/restart, bleed, bypass, and geometry effects. Numerical solution of the governing time dependent equations of motion is accomplished through a shock capturing finite difference algorithm, of which five separate approaches are evaluated. Comparison with experimental supersonic wind tunnel data is presented to verify the present approach for a wide range of transient inlet flow conditions.

Varner, M. O.↗

NASA NTRS1985

Numerical solution of the unsteady Navier-Stokes equation

The construction and the analysis of nonoscillatory shock capturing methods for the approximation of hyperbolic conservation laws are discussed. These schemes share many desirable properties with total variation diminishing schemes, but TVD schemes have at most first-order accuracy, in the sense of truncation error, at extrema of the solution. In this paper a uniformly second-order approximation is constructed, which is nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time. This is achieved via a nonoscillatory piecewise linear reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell.

Osher, Stanley J.↗

NASA NTRS1986

Transonic flow solutions using a composite velocity procedure for potential, Euler and RNS equations

Solutions for transonic viscous and inviscid flows using a composite velocity procedure are presented. The velocity components of the compressible flow equations are written in terms of a multiplicative composite consisting of a viscous or rotational velocity and an inviscid, irrotational, potential-like function. This provides for an efficient solution procedure that is locally representative of both asymptotic inviscid and boundary layer theories. A modified conservative form of the axial momentum equation that is required to obtain rotational solutions in the inviscid region is presented and a combined conservation/nonconservation form is applied for evaluation of the reduced Navier-Stokes (RNS), Euler and potential equations. A variety of results is presented and the effects of the approximations on entropy production, shock capturing, and viscous interaction are discussed.

Gordnier, R. E.↗

NASA NTRS1987

Numerical simulation of shock wave diffraction by TVD schemes

An upwind total variation diminishing (TVD) scheme and a predictor-corrector symmetric TVD scheme were used to numerically simulate the blast wave diffraction on a stationary object. The objective is to help design an optimum configuration so that lateral motion is minimized and at the same time vortex shedding and flow separation are reduced during a blast wave encounter. Results are presented for a generic configuration for both a coarse grid and a fine grid to illustrate the global and local diffraction flow fields. Numerical experiments for the shock wave reflection on a wedge are also included to validate the current approach. Numerical study indicated that these TVD schemes are more stable and produced higher shock resolution than classical shock capturing methods such as the explicit MacCormack scheme.

Young, Victor Y. C.↗

NASA NTRS1986

A Petrov-Galerkin finite element method for the compressible Euler and Navier-Stokes equations

An overview of a new finite element method for the compressible Euler and Navier-Stokes equations is presented. The discretization is based on entropy variables. The method is developed within the framework of a Petrov-Galerkin formulation. Two perturbations are added to the weighting function; one is a generalization of the SUPG operator and the other is designed to enhance shock capturing capability. The treatment of boundary conditions and the consistent calculation of boundary fluxes are addressed. Results of numerical tests are presented which confirm the robustness and wide applicability of the method.

Franca, L. P.↗

NASA NTRS1987

An implicit flux-split algorithm to calculate hypersonic flowfields in chemical equilibrium

An implicit, finite-difference, shock-capturing algorithm that calculates inviscid, hypersonic flows in chemical equilibrium is presented. The flux vectors and flux Jacobians are differenced using a first-order, flux-split technique. The equilibrium composition of the gas is determined by minimizing the Gibbs free energy at every node point. The code is validated by comparing results over an axisymmetric hemisphere against previously published results. The algorithm is also applied to more practical configurations. The accuracy, stability, and versatility of the algorithm have been promising.

Palmer, Grant↗

NASA NTRS1987

Uniformly high order accurate essentially non-oscillatory schemes. III

In the present continuation of development and analysis efforts towards essentially nonoscillatory shock-capturing approximations of hyperbolic conservation laws, the Godunov (1959) scheme and its MUSCL second-order-accurate extension are generalized to an arbitrary order of accuracy by means of a hierarchy of uniformly high-order-accurate schemes. The design method, which involves the solution's essentially nonoscillatory piecewise polynomial reconstruction on the basis of its cell averages, followed by time-evolution through an approximate solution of the resulting initial value problem and its averaging over each cell, employs an adaptive stencil of grid points to yield schemes that are highly nonlinear.

Harten, Ami↗

NASA NTRS1987

Modification of the Osher upwind scheme for use in three dimensions

The Osher upwind scheme, a shock-capturing algorithm based on an approximate Riemann solver, which has been developed for the solution of hyperbolic conservation laws, is briefly described, and the need for modification of the scheme for use in three dimensions is demonstrated. Two forms of a scheme modification for use on a three-dimensional arbitrary grid are then presented.

Hessenius, Kristin A.↗

NASA NTRS1988

An upwind parabolized Navier-Stokes code for real gas flows

A real gas, upwind, parabolized Navier-Stokes (PNS) code has been developed to compute the two-dimensional/axisymmetric hypersonic flow of equilibrium air around various body shapes. The new code is an extension of the upwind (perfect gas) PNS code of Lawrence, Tannehill and Chaussee. The upwind algorithm is based on Roe's flux-difference splitting scheme which has been modified to account for real gas effects. Simplified curve fits are used to obtain the thermodynamic and transport properties of equilibrium air. The new code has been validated by computing the hypersonic laminar flow of air over a flat plate, a wedge, a ramp, and a cone. The results of these computations are compared with the results from a conventional centrally-differenced, real gas, PNS code and the agreement is excellent, except in the vicinity of shock waves where the present code exhibits far superior shock capturing capabilities.

Tannehill, John C.↗

NASA NTRS1988

An upwind parabolized Navier-Stokes code for chemically reacting flows

A new upwind, parabolized Navier-Stokes (PNS) code has been developed to compute the hypersonic, viscous, chemically reacting flow around two-dimensional or axisymmetric bodies. The new code is an extension of the upwind (perfect gas) PNS code of Lawrence et al. (1986). The upwind algorithm is based on Roe's flux-difference splitting scheme which has been modified to account for real gas effects. The algorithm solves the gas dynamic and species continuity equations in a 'loosely' coupled manner. The new code has been validated by computing the laminar flow (at free stream Mach number 25) of chemically reacting air over a wedge and a cone. The results of these computations are compared with the results from a centrally-differenced, fully coupled, nonequilibrium PNS code. The agreement is excellent, except in the vicinity of the shock wave where the present code exhibits superior shock capturing capabilities.

Tannehill, John C.↗

NASA NTRS1988

Artificial dissipation models for hypersonic external flow

Four artificial dissipation models which augment central difference schemes were examined for hypersonic external flows. The models were a first and third order dissipation model, a directionally scaled first and third order dissipation model, a flux limited dissipation model, and a flux difference split dissipation model. Each model was implemented in the lower-upper symmetric-Gauss-Seidel (LU-SGS) algorithm to solve the full Navier-Stokes equations. The latter two models can be regarded as total variation diminishing (TVD) schemes. Test results for model problems showed that the flux limited dissipation model was robust enough to predict a high speed blunt body flow with strong shock and expansion waves. The flux difference split dissipation model was capable of shock capturing with higher resolution, but was less robust. First and third order dissipation models turned out to be neither accurate nor robust enough for high Mach number flow computations.

Yoon, Seokkwan↗

NASA NTRS1988

An improved flux-split algorithm applied to hypersonic flows in chemical equilibrium

An explicit, finite-difference, shock-capturing numerical algorithm is presented and applied to hypersonic flows assumed to be in thermochemical equilibrium. Real-gas chemistry is either loosely coupled to the gasdynamics by way of a Gibbs free energy minimization package or fully coupled using species mass conservation equations with finite-rate chemical reactions. A scheme is developed that maintains stability in the explicit, finite-rate formulation while allowing relatively high time steps. The codes use flux vector splitting to difference the inviscid fluxes and employ real-gas corrections to viscosity and thermal conductivity. Numerical results are compared against existing ballistic range and flight data. Flows about complex geometries are also computed.

Palmer, Grant↗

NASA NTRS1988

A comparison of ENO and TVD schemes

The numerical performance of a second-order upwind-based TVD scheme is compared with that of a uniform second-order ENO scheme on shock capturing. The cases considered include flows with Mach numbers of 2.9, 5.0, and 10.0. For cases with Mach numbers of 5.0 and 10.0, the computed ENO results are inferior to the corresponding TVD results.

Chang, Shih-Hung↗

NASA NTRS1988

The development of flux-split algorithms for flows with non-equilibrium thermodynamics and chemical reactions

A finite-volume method for the numerical computation of flows with nonequilibrium thermodynamics and chemistry is presented. A thermodynamic model is described which simplifies the coupling between the chemistry and thermodynamics and also results in the retention of the homogeneity property of the Euler equations (including all the species continuity and vibrational energy conservation equations). Flux-splitting procedures are developed for the fully coupled equations involving fluid dynamics, chemical production and thermodynamic relaxation processes. New forms of flux-vector split and flux-difference split algorithms are embodied in a fully coupled, implicit, large-block structure, including all the species conservation and energy production equations. Several numerical examples are presented, including high-temperature shock tube and nozzle flows. The methodology is compared to other existing techniques, including spectral and central-differenced procedures, and favorable comparisons are shown regarding accuracy, shock-capturing and convergence rates.

Grossman, B.↗

NASA NTRS1988

Universal limiter for transient interpolation modeling of the advective transport equations: The ULTIMATE conservative difference scheme

A fresh approach is taken to the embarrassingly difficult problem of adequately modeling simple pure advection. An explicit conservative control-volume formation makes use of a universal limiter for transient interpolation modeling of the advective transport equations. This ULTIMATE conservative difference scheme is applied to unsteady, one-dimensional scalar pure advection at constant velocity, using three critical test profiles: an isolated sine-squared wave, a discontinuous step, and a semi-ellipse. The goal, of course, is to devise a single robust scheme which achieves sharp monotonic resolution of the step without corrupting the other profiles. The semi-ellipse is particularly challenging because of its combination of sudden and gradual changes in gradient. The ULTIMATE strategy can be applied to explicit conservation schemes of any order of accuracy. Second-order schemes are unsatisfactory, showing steepening and clipping typical of currently popular so-called high resolution shock-capturing of TVD schemes. The ULTIMATE third-order upwind scheme is highly satisfactory for most flows of practical importance. Higher order methods give predictably better step resolution, although even-order schemes generate a (monotonic) waviness in the difficult semi-ellipse simulation. Little is to be gained above ULTIMATE fifth-order upwinding which gives results close to the ultimate for which one might hope.

Leonard, B. P.↗

NASA NTRS1989

Transonic flow solutions using a composite velocity procedure for potential, Euler and RNS equations

Gordnier, R. E.↗

NASA NTRS1989

Integral equation solution of the full potential equation for transonic flows

An integral equation method for solving the full potential equation has been developed for arbitrary configurations in two- or three-dimensional transonic flows. This method is capable of capturing shocks using Murman-Cole type of finite difference scheme and is capable of predicting accurate and force-free wake shape as well. A rectangular grid combined with a technique of local grid refinement greatly improved the computational efficiency.

Chu, Li-Chuan↗

NASA NTRS1989

Comparisons of TVD schemes applied to the Navier-Stokes equations

In this study, the following total variation diminishing (TVD) schemes for solving the Navier-Stokes equations have been tested: the Chakravarthy and Szema (1985) upwind biased TVD scheme, the Harten's upwind TVD scheme described by Yee et al. (1983), and the Yee's (1985) symmetric TVD scheme. The schemes have been compared using three test cases. The first case was the one-dimensional shock tube problem which tested the shock-capturing abilities of the schemes. Chakravarthy's and Harten's schemes gave similar results which were found to be more accurate than the results from Yee's scheme. The second case was a compressible boundary layer which tested the schemes's abilities to solve fiscous flows. In this case, the three schemes yielded almost identical results. Finally, the shock/boundary-layer interaction case studied experimentally by Hakkinen et al. (1959) was computed. Here, Chakravarthy's and Yee's schemes compared most favorably with the published data, with Yee's scheme giving slightly better results.

Buelow, Philip E.↗