SEARCH · Search NASA

Results for “Shock fitting”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS2006

Multi-Spacecraft Observations of Interplanetary Shock Accelerated Particle Events

We use simultaneous measurements from the Wind and ACE spacecraft to determine the spatial properties of both interplanetary (IP) shocks and the shock-associated energetic particle events. We combine plasma, magnetic field and energetic particle data from ACE and Wind for 124 energetic storm particle (ESP) events from 1998 to 2003 and examine the spatial and temporal variations of these events in the Earth's vicinity. We find that even though the two spacecraft were occasionally separated by more than 400 RE, the plasma, field, and energetic particle time-intensity profiles during the events were very similar. In addition, we find that the ion composition and energy spectra in individual IP shock events are identical at the two spacecraft locations. We also use the fitted shock velocity along the normal from ACE and estimate the shock transit time to Wind location. In general, there is poor agreement between the estimated transit time and the actual measured transit time. Hence, our assumptions that a) the IP shock at 1 AU propagates radially, and/or b) the IP shock is spherically symmetric at 1 AU are not valid. In this paper, we will also study, for the first time, the anisotropy measurements of low-energy IP shock-associated ions at both ACE and Wind. We will then compare these new anisotropy analyses with locally measured shock parameters and identify possible signatures of different shock acceleration processes as predicted by the first-order Fermi and shock-drift models.

Ho, G. C.↗

NASA NTRS1983

Psuedospectral calculation of shock turbulence interactions

A Chebyshev-Fourier discretization with shock fitting is used to solve the unsteady Euler equations. The method is applied to shock interactions with plane waves and with a simple model of homogeneous isotropic turbulence. The plane wave solutions are compared to linear theory.

Zang, T. A.↗

NASA NTRS1983

Pseudospectral calculation of shock turbulence interactions

Zang, T. A.↗

NASA NTRS1983

Transient flow analysis of the AEDC/HPDE MHD generator

A hybrid Lax-Wendroff/Method of Characteristics computer code has been developed for numerical simulation of flow transients associated with the operation of MHD generator facilities. The code employs the shock-fitting method, with an Eulerian formulation of the basic conservation equations and explicit tracking of shock waves. Pressure, temperature, and velocity are used as primary integration variables to simplify interfacing of the code with real-gas thermodynamic and transport property tables. Application of the code to the simulation of selected transients for the AEDC/HPDE MHD generator produced results that are in good agreement with experimental observations.

Wilson, D. R.↗

NASA NTRS1991

Spectral methods for the Euler equations - The blunt body problem revisited

The present use of the Chebyshev spectral collocation method, in conjunction with shock-fitting, to solve the blunt-body problem gives attention to the boundary and the shock-acceleration equations. The crux of these procedures is the use of the characteristic compatibility relations to compute the body pressure and shock velocity. It is shown that converged solutions are obtainable without artificial smoothing, and that spectral accuracy is achieved.

Kopriva, David A.↗

NASA NTRS1982

A computer program for the calculation of the flow field including boundary layer effects for mixed-compression inlets at angle of attack

A computer program was developed which is capable of calculating the flow field in the supersonic portion of a mixed compression aircraft inlet operating at angle of attack. The supersonic core flow is computed using a second-order three dimensional method-of-characteristics algorithm. The bow shock and the internal shock train are treated discretely using a three dimensional shock fitting procedure. The boundary layer flows are computed using a second-order implicit finite difference method. The shock wave-boundary layer interaction is computed using an integral formulation. The general structure of the computer program is discussed, and a brief description of each subroutine is given. All program input parameters are defined, and a brief discussion on interpretation of the output is provided. A number of sample cases, complete with data listings, are provided.

Vadyak, J.↗

NASA NTRS1995

Neptune inbound bow shock

Voyager 2 crossed the inbound or upstream Neptunian bow shock at 1430 spacecraft event time on August 24 in 1989 (Belcher et al., 1989). The plasma and magnetic field measurements allow us to study the solar wind interaction with the outermost gas giant. To fully utilize all of the spacecraft observations, an improved nonlinear least squares, 'Rankine-Hugoniot' magnetohydrodynamic shock-fitting technique has been developed (Szabo, 1994). This technique is applied to the Neptunian data set. We find that the upstream bow shock normal points nearly exactly toward the Sun consistent with any reasonable large-scale model of the bow shock for a near subsolar crossing. The shock was moving outward with a speed of 14 +/- 12 km/s. The shock can be characterized as a low beta, high Mach number, strong quasi-perpendicular shock. Finally, the shock microstructure features are resolved and found to scale well with theoretical expectations.

Szabo, Adam↗

NASA NTRS1985

A 3-dimensional Navier-Stokes-Euler code for blunt-body flow computations

The shock-layer flowfield is obtained with or without viscous and heat-conducting dissipations from the conservative laws of fluid dynamics equations using a shock-fitting implicity finite-difference technique. The governing equations are cast in curvilinear-orthogonal coordinates and transformed to the domain between the shock and the body. Another set of equations is used for the singular coordinate axis, which, together with a cone generator away from the stagnation point, encloses the computation domain. A time-dependent alternating direction implicit factorization technique is applied to integrate the equations with local-time increment until a steady solution is reached. The shock location is updated after the flowfield computation, but the wall conditions are implemented into the implicit procedure. Innovative procedures are introduced to define the initial flowfield, to treat both perfect and equilibrium gases, to advance the solution on a coarse-to-fine grid sequence, and to start viscous flow computations from their corresponding inviscid solutions. The results are obtained from a grid no greater than 28 by 18 by 7 and converged within 300 integration steps. They are of sufficient accuracy to start parabolized Navier-Stokes or Euler calculations beyond the nose region, to compare with flight and wind-tunnel data, and to evaluate conceptual designs of reentry spacecraft.

Li, C. P.↗

NASA NTRS1988

Integral equation solution for transonic and subsonic aerodynamics

Two methods are presented to solve for the subsonic and transonic flows around airfoils. The first method is based on the integral solution of the full-potential equation with a shock-capturing technique only or with shock capturing-shock fitting technique. In the second method, the integral soluton of the full potential equation is coupled with an embedded region of Euler equations around the shock location. The second method is a computationally efficient technique for flows with strong shocks where the entropy increase and vorticity production across the shock are not small. Several numerical examples are presented and compared with the experimental data and other computational results.

Kandil, Osama A.↗

NASA NTRS1993

Comparison of coordinate-invariant and coordinate-aligned upwinding for the Euler equations

A floating-shock fitting method for the Euler equations has been developed that uses one-sided spatial differences along and across streamlines. This method has been applied to unsteady shocked flow in a duct with a ramp, to supercritical flow over a circular cylinder, and to transonic flow over airfoils. Compared to methods using coordinate-aligned upwind differencing, the coordinate-invariant upwinding generally required fewer grid points to produce crisp shocks and shears. For transonic airfoils, coordinate-invariant upwind differencing advances the agreement between computations and experiment, and it reduces the grid dependency of the computed results for strong shocks. The computational expenditure is comparable to that of coordinate-aligned upwind methods for the Euler equations.

Hartwich, Peter M.↗

NASA NTRS1996

Investigation of shock-induced combustion past blunt projectiles

A numerical study is conducted to simulate shock-induced combustion in premixed hydrogen-air mixtures at various free-stream conditions and parameters. Two-dimensional axisymmetric, reacting viscous flow over blunt projectiles is computed to study shock-induced combustion at Mach 5.11 and Mach 6.46 in hydrogen-air mixture. A seven-species, seven reactions finite rate hydrogen-air chemical reaction mechanism is used combined with a finite-difference, shock-fitting method to solve the complete set of Navier-Stokes and species conservation equations. The study has allowed an improved understanding of the physics of shock-induced combustion over blunt projectiles and the numerical results can now be explained more readily with one-dimensional wave-interaction model.

Ahuja, J. K.↗

NASA NTRS1982

Three-dimensional flow calculations including boundary layer effects for supersonic inlets at angle of attack

An analysis is presented for calculating the steady three-dimensional flow field in supersonic mixed-compression inlets at incidence. A zonal modeling approach is employed to obtain the solution. The supersonic core flow is computed using a second-order pentahedral bicharacteristic algorithm. The bow shock wave and the reflected internal shock train are determined using a three-dimensional discrete shock fitting procedure. The boundary layer flow adjacent to both the centerbody and the cowl is computed using a second-order implicit finite difference method. The flow in a shock wave-boundary layer interaction region is computed using an integral formulation. The culmination of the present research effort is the development of a production-type computer program capable of analyzing flow in a variety of mixed-compression aircraft inlets. Numerical results and experimental correlations are presented to illustrate application of the analysis.

Vadyak, J.↗

NASA NTRS1988

Numerical simulation of hypersonic flow over an aeroassist flight experiment vehicle

Solutions of wind-tunnel and entry-flight flow around the vehicle are obtained from the Navier-Stokes equations coupled with the chemical species continuity equations if needed. The time-iterative method employs several techniques: shock fitting, chemistry-split ADI and an algebraic grid in conformal spherical-polar space. Sensitivities of the results to numerical parameters and to frozen, equilibrium and finite rate reactions are investigated in the forebody computation. Quantitative results are obtained for the shock layer and the near wake for the entire vehicle corresponding to both ground test and flight conditions. Complex flow characteristics are analyzed on the basis of the complete flowfield over the aerobrake and simplified afterbodies. The method is stable and cost effective, and has yielded shock locations and wall pressure distributions which are in good agreement with wind-tunnel data.

Li, C. P.↗

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 NTRS1993

Spectral solution of the viscous blunt-body problem

The viscous blunt-body problem is solved with a shock-fitted Chebyshev spectral method. No explicit artificial viscosity or filtering is needed to obtain smooth, converged solutions. The method is applied to two problems. First, results for the flow over a right circular cylinder in the Mach number range of 5.5-6.0 are compared with experimental data. Second, a solution for a Mach 25 flow over a hyperbolic cone is compared with a viscous shock-layer calculation.

Kopriva, David A.↗

NASA NTRS1988

Analysis and Optimization of Truncated Scarf Nozzles Subject to External Flow Conditions

Results of a calculation of an optimized truncated scarfed nozzle were compared. The truncated scarfed nozzle was designed for an exit Mach number of 6.0, i.e., the Mach number at the last nozzle characteristic is 6.0, with an external flow Mach number of 5.0. The nozzle was designed by the Rao method for optimum thrust nozzles modified for 2-D flow and truncated scarfed nozzle applications. This design was analyzed using a shock-fitting method for 2-D supersonic flows. Excellent agreement was achieved between the design and analysis. Truncation of the lower nozzle wall (cowl) revealed that there is an optimum length for truncating the cowl without degrading the nozzle performance. Truncation of the nozzle cowl past this optimal length should be analyzed in trade-off studies for thrust loss versus gross vehicle weight. Plots of the oblique shock wave equations were also identified which will allow computation of slip line angle, dynamic pressure coefficient, or ambient Mach number for various specific heat ratios.

Shyne, Rickey J.↗

NASA NTRS1982

RAXJET: A computer program for predicting transonic, axisymmetric flow over nozzle afterbodies with supersonic jet exhausts

A viscous-inviscid interaction method to calculate the subsonic and transonic flow over nozzle afterbodies with supersonic jet exhausts was developed. The method iteratively combines a relaxation solution of the full potential equation for the inviscid external flow, a shock capturing-shock fitting inviscid jet solution, an integral boundary layer solution, a control volume method for treating separated flows, and an overlaid mixing layer solution. A computer program called RAXJET which incorporates the method, illustrates the predictive capabilities of the method by comparison with experimental data is described, a user's guide to the computer program is provided. The method accurately predicts afterbody pressures, drag, and flow field properties for attached and separated flows for which no shock induced separation occurs.

Wilmoth, R. G.↗

NASA NTRS1977

Supersonic flow about circular cones at large angles of attack - A floating discontinuity approach

The technique of floating shock fitting is adapted to the computation of the inviscid flowfield about circular cones in a supersonic free stream at angles of attack that exceed the cone half-angle. In those regions in which the governing conical equations are mixed elliptic-hyperbolic, the fully hyperbolic form is obtained by the addition of the temporal derivative. The resulting equations are applicable over the complete range of free-stream Mach numbers, angles of attack and cone half-angles for which the bow shock is attached. An explicit finite-difference algorithm is used to obtain the solution by an unsteady relaxation approach. The bow shock, embedded crossflow shock, and vortical singularity in the leeward symmetry plane are all treated as floating discontinuities in a fixed computational mesh. The method yields excellent results for the bow and embedded shocks, however, the solution in the leeward symmetry plane exhibits viscous-like effects and does not appear to adequately predict the behavior of the vortical singularity.

Daywitt, J.↗