Engineering Papers⌕ Search

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.

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

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

The Origin of Inlet Buzz in a Mach 1.7 Low Boom Inlet Design

Supersonic inlets with external compression, having a good level performance at the critical operating point, exhibit a marked instability of the flow in some subcritical operation below a critical value of the capture mass flow ratio. This takes the form of severe oscillations of the shock system, commonly known as "buzz". The underlying purpose of this study is to indicate how Detached Eddy Simulation (DES) analysis of supersonic inlets will alter how we envision unsteady inlet aerodynamics, particularly inlet buzz. Presented in this paper is a discussion regarding the physical explanation underlying inlet buzz as indicated by DES analysis. It is the normal shock wave boundary layer separation along the spike surface which reduces the capture mass flow that is the controlling mechanism which determines the onset of inlet buzz, and it is the aerodynamic characteristics of a choked nozzle that provide the feedback mechanism that sustains the buzz cycle by imposing a fixed mean corrected inlet weight flow. Comparisons between the DES analysis of the Lockheed Martin Corporation (LMCO) N+2 inlet and schlieren photographs taken during the test of the Gulfstream Large Scale Low Boom (LSLB) inlet in the NASA 8x6 ft. Supersonic Wind Tunnel (SWT) show a strong similarity both in turbulent flow field structure and shock wave formation during the buzz cycle. This demonstrates the value of DES analysis for the design and understanding of supersonic inlets.

Mach 1.7↗
NASA NTRS1978

Lagrangian computation of inviscid compressible flows

A Lagrangian method is developed to solve the Euler equations of gas dynamics. The solution of the equations is obtained by a numerical computation with the well-known Flux-Corrected-Transport (FCT) numerical method. This procedure is modified so that the boundary treatment is accurate and relatively simple. Shock waves and other flow discontinuities are captured monotonically without any type of fitting procedures. The Lagrangian method is employed so that the problem of mesh generation is completely avoided. The method is applicable to all Mach numbers except the low subsonic range where compressibility effects are small. The method is applied to a one-dimensional Riemann problem (shock tube) and to a two-dimensional supersonic channel flow with reflecting shock waves.

Klopfer, G. H.↗
NASA NTRS1988

A computational analysis of under-expanded jets in the hypersonic regime

Underexpanded axisymmetric jets are studied numerically using a full Navier-Stokes solver. Emphasis has been given to supersonic and hypersonic jets in supersonic and hypersonic ambient flows, a phenomenon previously overlooked. It is demonstrated that the shear layers and shock patterns in a jet plume can be captured without complicated viscous/inviscid and subsonic/supersonic coupling schemes. In addition, a supersonic pressure relief effect has been identified for underexpanded jets in supersonic ambient flows. While it is well known that an underexpanded jet in a quiescent ambience (or subsonic ambience) contains multiple shock cells, the present study shows that because of the supersonic pressure relief effect, an underexpanded jet in a supersonic or hypersonic ambience contains only one major shock cell.

Hsu, Andrew T.↗
NASA NTRS1988

A computational analysis of under-expanded jets in the hypersonic regime

Underexpanded axisymmetric jets are studied numerically using a full Navier-Stokes solver. Emphasis has been given to supersonic and hypersonic jets in supersonic and hypersonic ambient flows, a phehomenon previously being overlooked for the most part. The present work demonstrates that the shear layers and shock patterns in a jet plume can be captured without complicated viscous/inviscid and subsonic/supersonic coupling schemes. In addition, a supersonic pressure relief effect has been identified for underexpanded jets in supersonic ambient flows. While it is well known that an underexpanded jet in a quiescent ambience (or subsonic ambience) contains multiple shock cells, the present study shows that because of the supersonic pressure relief effect, an underexpanded jet in a supersonic or hypersonic ambience contains only one major shock cell.

Hsu, Andrew T.↗
NASA NTRS1975

Numerical methods for the calculation of three-dimensional nozzle exhaust flow fields

Numerical codes developed for the calculation of three-dimensional nozzle exhaust flow fields associated with hypersonic airbreathing aircraft are described. Both codes employ reference plane grid networks with respect to three coordinate systems. Program CHAR3D is a characteristic code utilizing a new wave preserving network within the reference planes, while program BIGMAC is a finite difference code utilizing conservation variables and a one-sided difference algorithm. Secondary waves are numerically captured by both codes, while the underexpansion shock and plume boundary are treated discretely. The exhaust gas properties consist of hydrogen-air combustion product mixtures in local chemical equilibrium. Nozzle contours are treated by a newly developed geometry package based on dual cubic splines. Results are presented for simple configurations demonstrating two- and three-dimensional multiple wave interactions.

Dash, S. M.↗
NASA NTRS1980

Hyperbolic/parabolic development for the GIM-STAR code

Flow fields in supersonic inlet configurations were computed using the eliptic GIM code on the STAR computer. Spillage flow under the lower cowl was calculated to be 33% of the incoming stream. The shock/boundary layer interaction on the upper propulsive surface was computed including separation. All shocks produced by the flow system were captured. Linearized block implicit (LBI) schemes were examined to determine their application to the GIM code. Pure explicit methods have stability limitations and fully implicit schemes are inherently inefficient; however, LBI schemes show promise as an effective compromise. A quasiparabolic version of the GIM code was developed using elastical parabolized Navier-Stokes methods combined with quasitime relaxation. This scheme is referred to as quasiparabolic although it applies equally well to hyperbolic supersonic inviscid flows. Second order windward differences are used in the marching coordinate and either explicit or linear block implicit time relaxation can be incorporated.

Spradley, L. W.↗
NASA NTRS1993

Implicit multigrid Euler solutions with symmetric Total-Variation-Diminishing dissipation

A symmetric Total-Variation-Diminishing (TVD) formulation of the numerical dissipation terms has been incorporated into a diagonalized alternating direction implicit multigrid algorithm to solve the Euler equations of inviscid compressible flow. The new treatment of the dissipation makes is possible to capture both very strong and very weak shocks, virtually without oscillation for the steady flows of interest here. In addition, the TVD constraint fixes one of the two previously arbitrary constants in the formulation of the dissipation, and results in both converged solutions and convergence rates which are relatively insensitive to the choice of the remaining dissipation parameter.

Caughey, David A.↗
NASA NTRS2014

A Wolf-Rayet-Like Progenitor of SN 2013cu from Spectral Observations of a Stellar Wind

The explosive fate of massive Wolf-Rayet stars (WRSs) is a key open question in stellar physics. An appealing option is that hydrogen- deficient WRSs are the progenitors of some hydrogen-poor supernova explosions of types IIb, Ib and Ic. A blue object, having luminosity and colours consistent with those of some WRSs, has recently been identified in pre-explosion images at the location of a supernova of type Ib, but has not yet been conclusively determined to have been the progenitor. Similar work has so far only resulted in non-detections. Comparison of early photometric observations of type Ic supernovae with theoretical models suggests that the progenitor stars had radii of less than 10(exp 12) centimetres, as expected for some WRSs. The signature of WRSs, their emission line spectra, cannot be probed by such studies. Here we report the detection of strong emission lines in a spectrum of type IIb supernova 2013cu (iPTF13ast) obtained approximately 15.5 hours after explosion (by 'flash spectroscopy', which captures the effects of the supernova explosion shock breakout flash on material surrounding the progenitor star).We identify Wolf-Rayet-like wind signatures, suggesting a progenitor of the WN(h) subclass (those WRSs with winds dominated by helium and nitrogen, with traces of hydrogen). The extent of this dense wind may indicate increased mass loss from the progenitor shortly before its explosion, consistent with recent theoretical predictions.

observations↗
NASA NTRS1985

Nuclear reactions in shock wave front during supernova events

The new unique isotopic anomalous coponent of Xe(XeX) was found in the carbonaceous chondrites. It is enriched in light shielded isotopes (124Xe and 126Xe) and in heavy nonshielded isotopes (134Xe and 136Xe. All characteristics of Xe-X can be explained by a model of nucleosynthesis of the Xe isotopes in shock wave front passed through the He envelope during supernova events. The light isotopes are created by p process and the heavy isotopes are created by n process (slow r process). They were captured with high temperature carbon grains condensing by supernova shock waves.

Lavrukhina, A. K.↗
NASA NTRS1996

Effects of Shocks on Emission from Central Engines of Active Galactic Nuclei. I

In this paper we show that perturbations of the accretion flow within the central engines of some active galactic nuclei (AGNS) are likely to form shock waves in the accreting plasma. Such shocks, which may be either collisional or collisionless, can contribute to the observed high-energy temporal and spectral variability. Our rationale is the following: Observations show that the continuum emission probably originates in an optically thin, hot plasma in the AGN central engine. The flux and spectrum from this hot plasma varies significantly over light crossing timescales. Several authors have suggested that macroscopic perturbations contained within this plasma are the sources of this variability. In order to produce the observed emission the perturbations must be radiatively coupled with the optically thin hot matter and must also move with high velocities. We suggest that shocks, which can be very effective in randomizing the bulk motion of the perturbations, are responsible for this coupling. Shocks should form in the central engine, because the temperatures and magnetic fields are probably reduced below their virial values by radiative dissipation. Perturbations moving at Keplerian speeds, or strong non-linear excitations, result in supersonic and super-Alfvenic velocities leading to shock waves within the hot plasma. We show that even a perturbation smaller than the emitting region can form a shock that significantly modifies the continuum emission in an AGN, and that the spectral and temporal variability from such a shock generally resembles those of radio-quiet AGNS. As an example, the shock inducing perturbation in our model is a small main-sequence star, the capturing and eventual accretion of which are known to be a plausible process. We argue that shocks in the central engine may also provide a natural triggering mechanism for the "cold" component of Guilbert & Rees two-phase medium and an efficient mecha- nism for angular momentum transfer. Current and future missions, such as ASCA, XTE, XMM, AXAF, and ASTRO-E may determine the importance of shock-related emission from the central engines of AGNS.

Sivron, R.↗
NASA NTRS1985

Molecular processes in a high temperature shock layer

The development of techniques for the calculation of electron capture widths, electronic wave functions, cross sections and rates needed for the description of the dissociative recombination (DR) of molecular ions with electrons were described. The cross sections and rates were calculated by using harmonic oscillator wave functions for the ion and a delta function approximation for the continuum vibrational wave function in the repulsive dissociative channel. In order to obtain DR cross sections of quantitative accuracy, a computer program which solves the one dimensional nuclear motion wave equation was revised to calculate the cross sections and rates. The program and the new results are described. Included is a discussion of large windows found in the dissociative recombination cross sections from excited ion vibrational levels. These windows have not been previously reported in the literature. The magnitude of the DR cross sections for several dissociative routes are sensitive to the location of the crossing of the neutral and ion potential curves. Studies of the effects of basis set and CI wave function size on vertical excitation energies are described. Preliminary studies on N2 and O2 using large scale wave functions are also reported.

Guberman, S. L.↗
NASA NTRS2008

Experimental Plans for Subsystems of a Shock Wave Driven Gas Core Reactor

This Contractor Report proposes a number of plans for experiments on subsystems of a shock wave driven pulsed magnetic induction gas core reactor (PMI-GCR, or PMD-GCR pulsed magnet driven gas core reactor). Computer models of shock generation and collision in a large-scale PMI-GCR shock tube have been performed. Based upon the simulation results a number of issues arose that can only be addressed adequately by capturing experimental data on high pressure (approx.1 atmosphere or greater) partial plasma shock wave effects in large bore shock tubes ( 10 cm radius). There are three main subsystems that are of immediate interest (for appraisal of the concept viability). These are (1) the shock generation in a high pressure gas using either a plasma thruster or pulsed high magnetic field, (2) collision of MHD or gas dynamic shocks, their interaction time, and collision pile-up region thickness, and (3) magnetic flux compression power generation (not included here).

Kazeminezhad, F.↗
NASA NTRSDate not supplied

How Is the Solar Wind Shocked? A Distribution-Function Perspective From the MMS Mission

In this talk, we will discuss how plasma distribution functions from the Magnetospheric Multiscale (MMS) mission transformed the picture on how the solar wind is shocked. For quasi-parallel shocks with typical Alfven Mach numbers, we will examine how ion deceleration, reflection, and thermalization take place, and illustrate how the new picture expands on the framework envisioned by Schwartz and Burgess in 1991. For quasi-perpendicular shocks with the range of Alfven Mach numbers (~1-60) captured by MMS, we will examine examples showing the processes of ion and electron thermalization go beyond the existing paradigm and call for future investigation. The knowledge advanced by MMS measurements provides a ground to understand planetary shocks in our solar system as well as astrophysical shocks arising from stellar wind interaction with exoplanets and supernova remnants.

Li-Jen Chen↗
NASA NTRS1986

Application of NCOREL to aircraft configurations

The NCOREL computer program and methodology (Grossman, 1979) modified for computation of flows over complex geometries is described. The newly introduced features include a more flexible grid-generation package, capable of generating grids for realistic wing-body cross sections, and a more general numerical method, developed for the capture of highly three-dimensional embedded, oblique shocks. In addition, an inlet capability is included, by means of which the cross section is allowed to change discontinuously with the addition of the inlet geometry. The new cross section is regridded and the potentials and their derivatives are interpolated for the marching to continue. Hence, a discontinuity in the grid or mappings is taken into account at the inlet station. The results of the flow-field computation for a realistic fighter configuration are discussed.

Siclari, M. J.↗
NASA NTRS1990

A flux-split solution procedure for unsteady flow calculations

The solution of reduced Navier Stokes (RNS) equations is considered using a flux-split procedure. Unsteady flow in a two dimensional engine inlet is computed. The problems of unstart and restart are investigated. A sparse matrix direct solver combined with domain decomposition strategy is used to compute the unsteady flow field at each instant of time. Strong shock-boundary layer interaction, time varying shocks and time varying recirculation regions are efficiently captured.

Pordal, H. S.↗
NASA NTRS1992

The geocentric particulate distribution: Cometary, asteroidal, or space debris?

Definition of the Low Earth Orbit (LEO) particulate environment has been refined considerably with the analysis of data from NASA's Long Duration Exposure Facility (LDEF). Measurements of the impact rates from particulates ranging from sub-micron to millimetres in dimension and, especially, information on their directionality has permitted new scrunity of the sources of the particulates. Modelling of the dynamics of both bound (Earth orbital) and unbound (hyperbolic interplanetary) particulates intercepting LDEF's faces leads to the conclusion that the source is dominantly interplanetary for particle dimensions of greater than some 5 microns diameter; however the anisotropy below this dimension demands lower velocities and is compatible with an orbital component. Characteristics of the LDEF interplanetary component are compatible with familiar meteoroid sources and deep space measurements. Understanding of the orbital component which exceeds the interplanetary flux by a factor of 4 is less clear; although the very small particulates in orbit have been associated with space debris (Lawrance and Brownlee, 1986) this data conflicts with other measurements (McDonnell, Carey and Dixon, 1984) at the same epoch. By analysis of trajectories approaching the Earth and its atmosphere, we have shown that a significant contribution could be captured by aerocapture, i.e., atmospheric drag, from either asteroidal or cometary sources; such enhancement is unlikely however to provide the temporal and spatial fluctuations observed by the LDEF Interplanetary Dust Experiment (Mullholland et al. 1992). A further new mechanism is also examined, that of aerofragmentation capture, where an atmospheric grazing trajectory, which would not normally lead to capture, leads to fragmentation by thermal or mechanical shock; the microparticulates thus created can be injected in large numbers, but only into short-lifetime orbits. The concentration in one particular orbit plane, could explain the temporal fluctuations seen on LDEF; space debris could also explain the phenomenon.

Mcdonnell, J. A. M.↗
NASA NTRS1993

A new procedure for dynamic adaption of three-dimensional unstructured grids

A new procedure is presented for the simultaneous coarsening and refinement of three-dimensional unstructured tetrahedral meshes. This algorithm allows for localized grid adaption that is used to capture aerodynamic flow features such as vortices and shock waves in helicopter flowfield simulations. The mesh-adaption algorithm is implemented in the C programming language and uses a data structure consisting of a series of dynamically-allocated linked lists. These lists allow the mesh connectivity to be rapidly reconstructed when individual mesh points are added and/or deleted. The algorithm allows the mesh to change in an anisotropic manner in order to efficiently resolve directional flow features. The procedure has been successfully implemented on a single processor of a Cray Y-MP computer. Two sample cases are presented involving three-dimensional transonic flow. Computed results show good agreement with conventional structured-grid solutions for the Euler equations.

Biswas, Rupak↗
NASA NTRS1993

Numerical analysis of base flowfield at high altitude for a four-engine clustered nozzle configuration

The objective of this study is to benchmark a four-engine clustered nozzle base flowfield with a computational fluid dynamics (CFD) model. The CFD model is a pressure based, viscous flow formulation. An adaptive upwind scheme is employed for the spatial discretization. The upwind scheme is based on second and fourth order central differencing with adaptive artificial dissipation. Qualitative base flow features such as the reverse jet, wall jet, recompression shock, and plume-plume impingement have been captured. The computed quantitative flow properties such as the radial base pressure distribution, model centerline Mach number and static pressure variation, and base pressure characteristic curve agreed reasonably well with those of the measurement. Parametric study on the effect of grid resolution, turbulence model, inlet boundary condition and difference scheme on convective terms has been performed. The results showed that grid resolution and turbulence model are two primary factors that influence the accuracy of the base flowfield prediction.

Wang, Ten-See↗