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 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Computation of three-dimensional inviscid flow over hypersonic missile configurations using the GIM code

A three-dimensional computational technique was used to obtain flowfield solutions to the Euler equations over selected hypersonic missile configurations. The General Interpolants Method (GIM) computer code was used with interpolation functions in an algebraic approach to generate a discrete computational grid for each configuration. The spatial marching version of the GIM code, which treats the parabolized Navier-Stokes (PNS) equations or the Euler equations with a shock capturing, 'MacCormack-like' scheme, was used to advance the solution hyperbolically over each configuration. The inviscid flowfield solutions over the two three-dimensional missile configurations, calculated using the GIM hyperbolic scheme, are presented here. The flow field over a wing/body configuration at zero degree angle of attack is presented. Flow over the fuselage of a tactical missile, termed the TAME 10, at both zero degree and 7.5 degree angles of attack is presented. In addition, an inviscid, two-dimensional analysis of an inlet configuration designed to mount on the TAME 10 is included. Contour maps of velocity and pressure are included for each configuration. Comparison of calculation and data show good agreement.

Xiques, K. E.↗

NASA NTRS1982

Nonlinear Green's function method for unsteady transonic flows

Advantages to employing Green's function in describing unsteady three-dimensional transonic flows are explored. The development of the function for application to linear subsonic and supersonic unsteady aerodynamics is reviewed. It is shown that unique solutions are possible for external flows, with all functional expressions being defined in Prandtl-Glauert space. The development of methods of using the Green's function for transonic flows is traced, noting the necessity of including the effects of significant nonlinear terms. The steady-state problem is considered to demonstrate the shock-capturing ability of the method and the usefulness of the function in the incompressible, subsonic, transonic, and supersonic areas of potential unsteady three-dimensional flows around complex configurations. Computational time is asserted to be an order of magnitude less than with finite difference methods.

Tseng, K.↗

NASA NTRS1984

Spectral methods for compressible flow problems

Recent results concerning numerical simulation of shock waves using spectral methods are reviewed. Shock fitting techniques were discussed as well as shock capturing techniques with finite difference artificial viscosity. Also the notion of the information contained in the numerical results obtained by spectral methods and show how this information is recovered was discussed.

Gottlieb, D.↗

NASA NTRS1984

Implicit treatment of the unsteady full potential equation in conservation form

An implicit, conservative treatment for the unsteady full potential equation in two-dimensions is presented. The method employs a local time linearization for density, and introduces flux biasing concepts based on sonic conditions for the generation of artificial viscosity to capture shocks without any overshoots. The boundary condition is treated implicitly using a splitting procedure consistent with the approximate factorization scheme. This allows for extremely large Courant numbers, even for nonorthogonal grid at the body. The method has application not only to unsteady problems, but also to generate the starting blunt body solution for a supersonic full potential marching code. Results are presented for flows over cylinders, spheres and airfoils. Comparisons are made with available Euler and full potential results, and are in excellent agreement.

Shankar, V.↗

NASA NTRS1983

Global MHD model of the earth's magnetosphere

A global MHD model of the earth's magnetosphere is defined. An introduction to numerical methods for solving the MHD equations is given with emphasis on the shock-capturing technique. Finally, results concerning the shape of the magnetosphere and the plasma flows inside the magnetosphere are presented.

Wu, C. C.↗

NASA NTRS1984

Interactive phenomena in supersonic jet mixing problems. I Phenomenology and numerical modeling techniques

The interactive phenomena that occur in supersonic jet mixing flowfields, and numerical modeling techniques developed to analyze such phenomena are discussed. A spatial marching procedure based on solving the parabolized Navier-Stokes jet mixing equations is presented. This procedure combines shock-capturing methodology for the analysis of supersonic mixing regions with pressure-split methodology for the analysis of subsonic mixing regions. The two regions are coupled at viscous sonic lines utilizing a viscous-characteristic coupling procedure. Specialized techniques for the treatment of jet boundary growth, strong discontinuties (Mach disks), and small embedded subsonic zones (behind Mach disks) are presented. Turbulent processes are represented by two-equation turbulence model formulations. In Part II of this article, numerical studies are presented for a variety of supersonic jet interactive phenomena.

Dash, S. M.↗

NASA NTRS1985

Uniformly high-order accurate non-oscillatory schemes, 1

The construction and the analysis of nonoscillatory shock capturing methods for the approximation of hyperbolic conservation laws was begun. These schemes share many desirable properties with total variation diminishing schemes (TVD), but TVD schemes have at most first order accuracy, in the sense of truncation error, at extreme of the solution. A uniformly second order approximation was constucted, 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.

Harten, A.↗

NASA NTRS1984

Metric-discontinuous zonal grid calculations using the Osher scheme

Computations on zonal grids - in particular, grids with metric discontinuities resulting from the interspersion of highly clustered regions with coarse regions - are possible using a fully conservative form of the Osher upwind scheme. These zonal grids can result from an abrupt clustering of points near solution discontinuities or near other flow features that require improved resolution. The zonal approach is shown to capture shocks with almost 'shock-fitting' quality but with minimal effort. Results for inviscid flow, including quasi-one-dimensional nozzle flow, supersonic flow over a cylinder, and blast-wave diffraction by a ramp, are presented. These calculations demonstrate the powerful capabilities of the Osher scheme used in conjunction with zonal grids in simulating flow fields with complex shock patterns.

Rai, M. M.↗

NASA NTRS1985

Adaptive mesh solution for supersonic conical flow in a rectilinear inlet

Solutions for the inviscid and viscous supersonic conical flow in a complete rectilinear inlet consisting of four planes intersecting at arbitrary wedge and sweep angles are obtained. To compute the flow on a specially constructed mesh, a three-dimensional flow solver ARC3D is used. It is shown that a single-pass mesh re-adaptation procedure can be used to obtain improved shock capture without the necessity of modifying the flow code, provided that the flow solver works in curvilinear coordinates and is restartable. Results computed by the method show good agreement with experimental measurements and previous calculations.

Kerlick, G. D.↗

NASA NTRS1985

Large-scale computations in fluid mechanics; Proceedings of the Fifteenth Summer Seminar on Applied Mathematics, University of California, La Jolla, CA, June 27-July 8, 1983. Parts 1 & 2

Papers are presented on such topics as the use of semi-Lagrangian advective schemes in meteorological modeling; computation with high-resolution upwind schemes for hyperbolic equations; dynamics of flame propagation in a turbulent field; a modified finite element method for solving the incompressible Navier-Stokes equations; computational fusion magnetohydrodynamics; and a nonoscillatory shock capturing scheme using flux-limited dissipation. Consideration is also given to the use of spectral techniques in numerical weather prediction; numerical methods for the incorporation of mountains in atmospheric models; techniques for the numerical simulation of large-scale eddies in geophysical fluid dynamics; high-resolution TVD schemes using flux limiters; upwind-difference methods for aerodynamic problems governed by the Euler equations; and an MHD model of the earth's magnetosphere.

Engquist, B. E.↗

NASA NTRS1985

An MHD model of the earth's magnetosphere

It is pointed out that the earth's magnetosphere arises from the interaction of the solar wind with the earth's geomagnetic field. A global magnetohydrodynamics (MHD) model of the earth's magnetosphere has drawn much attention in recent years. In this model, MHD equations are used to describe the solar wind interaction with the magnetosphere. In the present paper, some numerical aspects of the model are considered. Attention is given to the ideal MHD equations, an equation of state for the plasma, the model as an initial- and boundary-value problem, the shock capturing technique, computational requirements and techniques for global MHD modeling, a three-dimensional mesh system employed in the global MHD model, and some computational results.

Wu, C. C.↗

NASA NTRS1985

Application of TVD schemes for the Euler equations of gas dynamics

Highly accurate and yet stable shock-capturing finite difference schemes have been designed for the computation of the Euler equations of gas dynamics. Four different principles for the construction of high resolution total variation diminishing (TVD) schemes are available, including hybrid schemes, a second-order extension of Godunov's scheme by van Leer (1979), the modified flux approach of Harten (1983, 1984), and the numerical fluctuation approach of Roe (1985). The present paper has the objective to review the class of second-order TVD schemes via the modified flux approach. Attention is given to first-order TVD schemes, a second-order accurate explicit TVD scheme, the global order of accuracy of the second-order TVD scheme, extensions to systems and two-dimensional conservation laws, numerical experiments with a second-order explicit TVD scheme, implicit TVD schemes, and second-order implicit TVD schemes.

Yee, H. C.↗

NASA NTRS1986

Nonlinear Supersonic Full Potential Analysis

Supersonic Implicit Marching Program (SIMP) applies numerical method, based on conservative form of full potential equation, to problem of threedimensional supersonic flows with embedded subsonic regions. Conservative formulation of problem provides ability to capture shocks and to assess accurately impact of sweep, thickness, and lift for conditions where linear theory unsatisfactory. Technique uses characteristic signal-propagation theory to control density biasing for treatment of shocks (including embedded shocks) and mixed elliptic/hyperbolic crossflow. SIMP written in FORTRAN 77.

Shankar, V.↗

NASA NTRS1985

Shock fitting applied to the prediction of high-speed rotor noise

A shock fitting method applied to the transonic small disturbance (TSD) potential equation is described. This method is then applied to a simple, two dimensional (2-D) rotating disturbance which is analogous to a shock radiating from the tip of a rotor blade in high speed hover. A comparison is made between the results of this method and the more standard shock capturing method. This comparison makes it clear that the effect of the results on the acoustic signature of the 2-D model is significant, and similar results can be expected when the method is extended to the three dimensional (3-D) case.

Rutherford, J. W.↗

NASA NTRS1986

An LU implicity scheme for high speed inlet analysis

A numerical method is developed to analyze the inviscid flowfield of a high speed inlet by the solution of the Euler equations. The lower-upper implicit scheme in conjunction with adaptive dissipation proves to be an efficient and robust nonoscillatory shock capturing technique for high Mach number flows as well as for transonic flows.

Yoon, S.↗

NASA NTRS1986

Uniformly high order accurate essentially non-oscillatory schemes 3

In this paper (a third in a series) the construction and the analysis of essentially non-oscillatory shock capturing methods for the approximation of hyperbolic conservation laws are presented. Also presented is a hierarchy of high order accurate schemes which generalizes Godunov's scheme and its second order accurate MUSCL extension to arbitrary order of accuracy. The design involves an essentially non-oscillatory piecewise polynomial 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. The reconstruction algorithm is derived from a new interpolation technique that when applied to piecewise smooth data gives high-order accuracy whenever the function is smooth but avoids a Gibbs phenomenon at discontinuities. Unlike standard finite difference methods this procedure uses an adaptive stencil of grid points and consequently the resulting schemes are highly nonlinear.

Harten, A.↗

NASA NTRS1986

Thermodynamic evaluation of transonic compressor rotors using the finite volume approach

The finite volume explicit time marching method was refined and improved. Previously, extension had been made to the finite volume method to improve the accuracy of the calculation of total pressure in inviscid flow, extend the method to allow the calculation of laminar and turbulent boundary layers in internal flows, and improve the shock capturing properties of the method by introducing a Mach number dependent interpolation scheme for the pressure used in the calculating the density. The current work extends these developments by using the new pressure interpolation scheme in two dimensional viscous calculations, including a more complete description of the viscous stresses, introducing a criteria for the transverse upwind differencing which is a function of the ratio of transverse and streamwise mass fluxes, and allowing the calculation of internal flow where boundary layers are present on both walls of the duct. The manner in which the viscous stresses are evaluated in the nonorthogonal, nonuniform grid is detailed. The convergence is investigated and results for calculations of laminar flow in a converging duct are presented. Results for calculations of transonic flow in a converging-diverging nozzle are presented and the results are compared with Sajben's measurements and calculations by others.

Nicholson, S.↗

NASA NTRS1985

MHD modelling of the earth's magnetosphere

Attention is given to the shock capturing technique, nonuniform grid system, and multiple time scale problem aspects of a global MHD model of the earth's magnetosphere that uses numerical methods. Because of the fast Alfven waves near the earth, the model is more difficult to solve than problems of hydrodynamic flow past bodies. The Rusanov (1962) scheme, which is a variant of the Lax (1954) method that maintains viscosity at a minimum value except where a large one is required, is used. The magnetosphere shape obtained by these means is noted.

Wu, C. C.↗