SEARCH · Search NASA

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Unsteady transonic airfoil computation using the integral solution of full-potential equation

The shock-capturing integral-equation scheme developed by Kandil and Hu (1987) for the analysis of steady transonic flow over airfoils is extended to the unsteady case. The full potential formulation of the governing equations is reviewed; the solution method is outlined; and results for a NACA 0012 airfoil in forced pitching oscillation at Mach 0.755 are presented in extensive graphs and briefly characterized. The present technique is shown to require significantly less computation time than finite-difference or finite-volume methods, and to give shock-motion predictions in good agreement with those of an implicit finite-volume Euler solver; the surface-pressure peaks are slightly underpredicted.

Kandil, O. A.↗

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.↗

NASA NTRS1988

Second order accurate finite difference approximations for the transonic small disturbance equation and the full potential equation

New shock-capturing finite difference approximations for solving two scalar conservation law nonlinear partial differential equations describing inviscid, isentropic, compressible flows of aerodynamics at transonic speeds are presented. A global linear stability theorem is applied to these schemes in order to derive a necessary and sufficient condition for the finite element method. A technique is proposed to render the described approximations total variation-stable by applying the flux limiters to the nonlinear terms of the difference equation dimension by dimension. An entropy theorem applying to the approximations is proved, and an implicit, forward Euler-type time discretization of the approximation is presented. Results of some numerical experiments using the approximations are reported.

Mostrel, M. M.↗

NASA NTRS1989

The development of an explicit thermochemical nonequilibrium algorithm and its application to compute three dimensional AFE flowfields

This study presents a three-dimensional explicit, finite-difference, shock-capturing numerical algorithm applied to viscous hypersonic flows in thermochemical nonequilibrium. The algorithm employs a two-temperature physical model. Equations governing the finite-rate chemical reactions are fully-coupled to the gas dynamic equations using a novel coupling technique. The new coupling method maintains stability in the explicit, finite-rate formulation while allowing relatively large global time steps. The code uses flux-vector accuracy. Comparisons with experimental data and other numerical computations verify the accuracy of the present method. The code is used to compute the three-dimensional flowfield over the Aeroassist Flight Experiment (AFE) vehicle at one of its trajectory points.

Palmer, Grant↗

NASA NTRS1989

Numerical study of unsteady viscous hypersonic blunt body flows with an impinging shock

A complex two-dimensional, unsteady, viscous hypersonic shock wave interaction is numerically simulated by a high-resolution, second-order fully implicit shock-capturing scheme. The physical model consists of a nonstationary oblique shock impinging on the bow shock of a blunt body. Studies indicate that the unsteady flow patterns are slightly different from their steady counterparts. However, for the sample cases investigated the peak surface pressures for the unsteady flows seem to occur at very different impingement locations than for the steady flow cases.

Klopfer, G. H.↗

NASA NTRS1989

Hypersonic blunt body computations including real gas effects

Various second-order explicit and implicit TVD shock-capturing methods, a generalization of Roe's approximate Riemann solver, and a generalized flux-vector splitting scheme are used to study two-dimensional hypersonic real-gas flows. Special attention is given to the identification of 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 numbers, for cases involving steady-state inviscid blunt flows. Blunt body calculations at Mach numbers of greater than 15 are performed to treat real-gas effects, and impinging shock results are obtained to test the treatment of slip surfaces and complex structures. Even with the addition of improvements, the convergence rate of algorithms in the hypersonic flow regime is found to be generally slower for a real gas than for a perfect gas.

Montagne, J.-L.↗

NASA NTRS1990

Integral solution of unsteady full-potential equation for a transonic pitching airfoil

The unsteady full-potential equation formulation in a moving frame of reference has been has been developed and used to solve unsteady transonic flow problems. An unsteady integral-equation shock-capturing (IE-SC) scheme has been developed. The resulting unsteady IE-SC scheme is applied to a NACA 0012 airfoil undergoing a pitching oscillation. The numerical results are compared with those of an implicit, approximately factored, finite-volume Euler scheme. The present scheme is efficient in terms of the number of iterations as compared to the other existing schemes, which use finite-difference or finite-volume methods.

Kandil, Osama A.↗

NASA NTRS1990

Enhanced thermochemical nonequilibrium computations of flow around the aeroassist flight experiment vehicle

A three-dimensional explicit, finite-rate, shock-capturing numerical algorithm is used to calculate thermochemical nonequilibrium flowfields about the Aeroassist Flight Experiment vehicle at one of its flight trajectory points. The full Navier-Stokes equations and an eleven species chemical model with the latest reaction rates are incorporated into the code. Results are compared against experimental data and other numerical solutions. The effects of changes to the physical model on the stagnation line and base region flow and surface quantities on the forebody are investigated.

Palmer, Grant↗

NASA NTRS1991

High-Resolution Numerical Simulation Of Shock Waves

NASA technical memorandum compares results of upwind and symmetric shock-capturing methods in numerical simulation of gas-dynamic flows. Both methods find shocks as sharp variations in fluid properties over few grid points. Methods differ in type of artificial viscosity introduced to stabilize computations. Symmetric method shown as accurate as upwind method, but with fewer and simpler time steps in transient case or fewer iteration steps in steady-state case.

Yee, H. C.↗

NASA NTRS1991

Overview Of Methods For Computation Of Shocks

NASA technical memorandum provides systematic overview of class of conservative finite-difference shock-capturing numerical-integration methods for solution of hyperbolic conservation laws. Unified and generalized formulation for one class of methods presented.

Yee, H. C.↗

NASA NTRS1991

Numerical simulation of fluid flow around a scramaccelerator projectile

Numerical simulations of the fluid motion and temperature distribution around a 'scramaccelerator' projectile are obtained for Mach numbers in the 5-10 range. A finite element method is used to solve the equations of motion for inviscid and viscous two-dimensional or axisymmetric compressible flow. The time-dependent equations are solved explicitly, using bilinear isoparametric quadrilateral elements, mass lumping, and a shock-capturing Petrov-Galerkin formulation. Computed results indicate that maintaining on-design performance for controlling and stabilizing oblique detonation waves is critically dependent on projectile shape and Mach number.

Pepper, Darrell W.↗

NASA NTRS1991

High-order ENO methods for the unsteady compressible Navier-Stokes equations

The adaptive stencil concepts of ENO (Essentially Non-Oscillatory) methods are applied to the laminar Navier-Stokes equations to yield a high-order, time-accurate algorithm with a shock-capturing capability. The method targets problems in the areas of nonlinear acoustics, compressible transition, and turbulence which, due to the presence of shocks or complex geometries, are not easily solved by spectral methods. The present approach has been implemented and tested for the full three-dimensional Navier-Stokes equations in a transformed curvilinear coordinate system. Validation results are presented for a variety of problems which verify the method's accuracy properties and shock capturing capabilities, as well as demonstrate its use as a direct simulation tool.

Atkins, H. L.↗

NASA NTRS1991

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 remove 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 viscous flow over a blunt biconic.

Turkel, E.↗

NASA NTRS1991

Flowfield computations and comparison with Shuttle aerodynamic data

An in-house developed flow solver, E3D, has been applied to investigate the flow field around the Shuttle Orbiter over an angle-of-attack range of 0 deg to 60 deg at Mach 3.5 and 10.0. The 3D Euler equations are integrated by means of a time-marching finite-volume shock-capturing method, based on cell-centered and upwind evaluation of the cell face fluxes with van Leer or Roe's flux-splitting techniques. A comparison of CFD, wind tunnel, and aeronautical data base is presented. Solutions are in very good agreement with available experimental data. Through the comparisons, it is demonstrated that the E3D code is capable of evaluating total aerodynamic performance.

Ma, E. C.↗

NASA NTRS1989

Hypersonic 3-D flow past winged bodies

The hypersonic flow past winged bodies is calculated using the 3D Navier-Stokes equations in both a space- and time-marching finite difference code. The bodies are a blunt ogive-cylinder with a delta wing planform and an elliptical body referred to as the allbody configuration. Laminar flow solutions for the ogive-cylinder-wing body and for the allbody configuration are presented. It is shown that the present technique has an excellent shock-capturing capability for all speed regimes. An axial zonal capability is incorporated to make the procedure more versatile. The grid system is created offline using an efficient hyperbolic grid generator which can handle the rapid variations in the body cross sections.

Chaussee, Denny S.↗

NASA NTRS1991

Multi-Dimensional ENO Schemes for General Geometries

A class of ENO schemes is presented for the numerical solution of multidimensional hyperbolic systems of conservation laws in structured and unstructured grids. This is a class of shock-capturing schemes which are designed to compute cell-averages to high order accuracy. The ENO scheme is composed of a piecewise-polynomial reconstruction of the solution form its given cell-averages, approximate evolution of the resulting initial value problem, and averaging of this approximate solution over each cell. The reconstruction algorithm is based on an adaptive selection of stencil for each cell so as to avoid spurious oscillations near discontinuities while achieving high order of accuracy away from them.

Harten, Ami↗

NASA NTRS1993

Conservative high-order-accurate finite-difference methods for curvilinear grids

Two fourth-order-accurate finite-difference methods for numerically solving hyperbolic systems of conservation equations on smooth curvilinear grids are presented. The first method uses the differential form of the conservation equations; the second method uses the integral form of the conservation equations. Modifications to these schemes, which are required near boundaries to maintain overall high-order accuracy, are discussed. An analysis that demonstrates the stability of the modified schemes is also provided. Modifications to one of the schemes to make it total variation diminishing (TVD) are also discussed. Results that demonstrate the high-order accuracy of both schemes are included in the paper. In particular, a Ringleb-flow computation demonstrates the high-order accuracy and the stability of the boundary and near-boundary procedures. A second computation of supersonic flow over a cylinder demonstrates the shock-capturing capability of the TVD methodology. An important contribution of this paper is the dear demonstration that higher order accuracy leads to increased computational efficiency.

Rai, Man M.↗

NASA NTRS1990

A comparison of optimization-based approaches for solving the aerodynamic design problem

Three optimization-based methods for solving aerodynamic design problems are compared. The Euler equations for one-dimensional duct flow was used as a model problem, and the three methods are compared for efficiency, robustness, and implementation difficulty. The smoothness of the design problem with respect to different shock-capturing finite difference schemes, and in the presence of grid refinement, is investigated.

Frank, Paul D.↗