Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock dynamics”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Effects of H2O Vapor on Vibrational Relaxation in Expanding and Contracting Flows

As opposed to previous explanations based on the effects of anharmonicity of simple diatomic molecules, traces of water vapor are suggested to be the most likely cause of the anomalously fast vibrational relaxation of such gases observed in supersonic and hypersonic nozzles. The mechanism is the strong V-VR coupling with H2O molecules that dramatically facilitates the collisional transfer of vibrational energy. Slight moisture content is thus a real world aspect of gas dynamics that must be considered in characterizations of shock tubes, reflected shock tunnels, and expansion tubes.

Meador, Willard E.↗
NASA NTRS1997

TURBO-AE: An Aeroelastic Code for Propulsion Applications

NASA's Advanced Subsonic Technology (AST) program is developing new technologies to increase the fuel efficiency of commercial aircraft engines, improve the safety of engine operation, and reduce engine emissions and noise. With the development of new designs for ducted fans, compressors, and turbines to achieve these goals, a basic aeroelastic requirement is that there should be no flutter or high resonant blade stresses in the operating regime. To verify the aeroelastic soundness of these designs, we need an accurate prediction and analysis code. Such a two-dimensional viscous propulsion aeroelastic code, named TURBO-AE, is being developed at the NASA Lewis Research Center. The TURBO-AE aeroelastic code is based on a three-dimensional unsteady aerodynamic Euler/Navier-Stokes turbomachinery code TURBO, developed under a grant from NASA Lewis. TURBO-AE can model viscous flow effects that play an important role in certain aeroelastic problems, such as flutter with flow separation (or stall flutter) and flutter in the presence of shock and boundary-layer interaction. The structural dynamics representation of the blade in the TURBO-AE code is based on a normal mode representation. A finite element analysis code, such as NASTRAN, is used to calculate in-vacuum vibration modes and the associated natural frequency. A work-per-cycle approach is used to determine aeroelastic (flutter) stability. With this approach, the motion of the blade is prescribed to be a harmonic vibration in a specified in vacuum normal mode. The aerodynamic forces acting on the vibrating blade and the work done by these forces on the vibrating blade during a cycle of vibration are calculated. If positive work is being done on the blade by the aerodynamic forces, the blade is dynamically unstable, since it will extract energy from the flow, leading to an increase in the amplitude of the blade's oscillation. Initial calculations have been done for a configuration representative of the Energy Efficient Engine fan rotor. The accompanying figure shows the work-per-cycle after each cycle of vibration. It can be seen that the work-per-cycle does not vary much after the fourth cycle. The negative sign of the converged work-per-cycle shows that the fan blade is dynamically stable and will not flutter. TURBO-AE will provide a useful aeroelastic prediction/analysis capability for engine manufacturers. It will reduce design cycle times by allowing new blade designs to be verified for aeroelastic soundness before blades are built and tested. With this prediction capability, it will be possible to build thinner, lighter, and faster rotating blades without encountering aeroelastic problems like stall flutter and high-cycle fatigue due to forced vibrations.

Bakhle, Milind A.↗
NASA NTRS2016

Investigation of Non-Equilibrium Radiation for Earth Entry

For Earth re-entry at velocities between 8 and 11.5 km/s, the accuracy of NASA's computational uid dynamic and radiative simulations of non-equilibrium shock layer radiation is assessed through comparisons with measurements. These measurements were obtained in the NASA Ames Research Center's Electric Arc Shock Tube (EAST) facility. The experiments were aimed at measuring the spatially and spectrally resolved radiance at relevant entry conditions for both an approximate Earth atmosphere (79% N2 : 21% O2 by mole) as well as a more accurate composition featuring the trace species Ar and CO2 (78.08% N2 : 20.95% O2 : 0.04% CO2 : 0.93% Ar by mole). The experiments were configured to target a wide range of conditions, of which shots from 8 to 11.5 km/s at 0.2 Torr (26.7 Pa) are examined in this paper. The non-equilibrium component was chosen to be the focus of this study as it can account for a significant percentage of the emitted radiation for Earth re-entry, and more importantly, non-equilibrium has traditionally been assigned a large uncertainty for vehicle design. The main goals of this study are to present the shock tube data in the form of a non-equilibrium metric, evaluate the level of agreement between the experiment and simulations, identify key discrepancies and to examine critical aspects of modeling non-equilibrium radiating flows. Radiance pro les integrated over discreet wavelength regions, ranging from the Vacuum Ultra Violet (VUV) through to the Near Infra-Red (NIR), were compared in order to maximize both the spectral coverage and the number of experiments that could be used in the analysis. A previously defined non-equilibrium metric has been used to allow comparisons with several shots and reveal trends in the data. Overall, LAURA/HARA is shown to under-predict EAST by as much as 40% and over-predict by as much as 12% depending on the shock speed. DPLR/NEQAIR is shown to under-predict EAST by as much as 50% and over-predict by as much as 20% depending on the shock speed. The one standard deviation scatter in the EAST results was calculated to be 31%. An estimate for the upper bound of the absolute error in wall-directed heat flux was calculated. Below 9 km/s, where the relative difference is large, the absolute error in radiative heat flux due to non-equilibrium models is estimated to be less then 1 W/sq cm. At the highest shock speed of 11 km/s, the error in non-equilibrium is estimated to be less than 20 W/sq cm.

Shock Tube↗
NASA NTRS2019

Stars at High Spatial Resolution

We summarize some of the compelling new scientific opportunities for understanding stars and stellar systems that can be enabled by sub-milliarcsec (sub-mas) angular resolution, UV/Optical spectral imaging observations, which can reveal the details of the many dynamic processes (e.g., evolving magnetic fields, accretion, convection, shocks, pulsations, winds, and jets) that affect stellar formation, structure, and evolution. These observations can only be provided by long baseline interferometers or sparse aperture telescopes in space, since the aperture diameters required are in excess of 500 m – a regime in which monolithic or segmented designs are not and will not be feasible - and since they require observations at wavelengths (UV) not accessible from the ground. Such observational capabilities would enable tremendous gains in our understanding of the individual stars and stellar systems that are the building blocks of our Universe and which serve as the hosts for life throughout the Cosmos.

Kenneth G Carpenter↗
NASA NTRS1996

Computation of Thermally Perfect Properties of Oblique Shock Waves

A set of compressible flow relations describing flow properties across oblique shock waves, derived for a thermally perfect, calorically imperfect gas, is applied within the existing thermally perfect gas (TPG) computer code. The relations are based upon a value of cp expressed as a polynomial function of temperature. The updated code produces tables of compressible flow properties of oblique shock waves, as well as the original properties of normal shock waves and basic isentropic flow, in a format similar to the tables for normal shock waves found in NACA Rep. 1135. The code results are validated in both the calorically perfect and the calorically imperfect, thermally perfect temperature regimes through comparisons with the theoretical methods of NACA Rep. 1135, and with a state-of-the-art computational fluid dynamics code. The advantages of the TPG code for oblique shock wave calculations, as well as for the properties of isentropic flow and normal shock waves, are its ease of use, and its applicability to any type of gas (monatomic, diatomic, triatomic, polyatomic, or any specified mixture thereof).

Tatum, Kenneth E.↗
NASA NTRS2012

Quasi 1D Modeling of Mixed Compression Supersonic Inlets

The AeroServoElasticity task under the NASA Supersonics Project is developing dynamic models of the propulsion system and the vehicle in order to conduct research for integrated vehicle dynamic performance. As part of this effort, a nonlinear quasi 1-dimensional model of the 2-dimensional bifurcated mixed compression supersonic inlet is being developed. The model utilizes computational fluid dynamics for both the supersonic and subsonic diffusers. The oblique shocks are modeled utilizing compressible flow equations. This model also implements variable geometry required to control the normal shock position. The model is flexible and can also be utilized to simulate other mixed compression supersonic inlet designs. The model was validated both in time and in the frequency domain against the legacy LArge Perturbation INlet code, which has been previously verified using test data. This legacy code written in FORTRAN is quite extensive and complex in terms of the amount of software and number of subroutines. Further, the legacy code is not suitable for closed loop feedback controls design, and the simulation environment is not amenable to systems integration. Therefore, a solution is to develop an innovative, more simplified, mixed compression inlet model with the same steady state and dynamic performance as the legacy code that also can be used for controls design. The new nonlinear dynamic model is implemented in MATLAB Simulink. This environment allows easier development of linear models for controls design for shock positioning. The new model is also well suited for integration with a propulsion system model to study inlet/propulsion system performance, and integration with an aero-servo-elastic system model to study integrated vehicle ride quality, vehicle stability, and efficiency.

Kopasakis, George↗
NASA NTRS1984

Detailed studies of the dynamics and energetics of coronal bullets

Coronal bullets are small ejecta of cool, dense plasma observed to accelerate through the solar atmosphere from 20 to 450 km/s. The NRL Dynamic Flux Tube Model has been used to simulate the evolving physical properties of these dynamic events. The present calculations utilize an adaptive-gridding technique to resolve the fine structure within and around the bullets. In this work, an identification was made of a component of shocked plasma which piles up ahead of the bullet and eventually dominates both the dynamics and heating of the original bullet mass. The observational consequences of this shocked component are discussed in terms of the available HRTS EUV data, and suggestions are made for optimizing future observations of this phenomenon. An investigation has also been conducted of the structure of the bullet material visible in EUV spectral lines and the observable characteristics of the EUV-emitting plasma. Finally, the most likely mechanisms for accelerating the bullets, as well as favorable sites of origin are evaluated.

Karpen, J. T.↗
NASA NTRS2017

Observations of Deep Ionospheric F-Region Density Depletions with FPMU Instrumentation and their Relationship with the Global Dynamics of the June 22-23, 2015 Geomagnetic Storm

The magnetic storm that commenced on June 22-23, 2015 was one of the largest storms in our current solar cycle, resulting from an active region on the Sun that produced numerous coronal mass ejections (CMEs) and associated interplanetary shock waves. On June 22 at 18:36 UT the magnetosphere was impacted by the shock wave on the magnetosphere. Observations from several spacecraft observed the dynamic response of the magnetosphere and ionosphere. MMS observatories in the near earth tail These low altitude measurements are correlated in the magnetosphere with particle flux dropouts measured by MMS We follow the timing of this storm in the ionosphere with the density depletions throughout the ISS orbits, DMSP drift velocities, and enhanced AMPERE Birkland currents. Together these observations and simulation results will be assembled to provide each region's context to the global dynamics and time evolution of the storm. The models during these event support and flesh out the puzzle of the global dynamics.

Coffey, Victoria↗
NASA NTRS2006

Shock Structure Analysis and Aerodynamics in a Weakly Ionized Gas Flow

The structure of a shock wave propagating through a weakly ionized gas is analyzed using an electrofluid dynamics model composed of classical conservation laws and Gauss Law. A viscosity model is included to correctly model the spatial scale of the shock structure, and quasi-neutrality is not assumed. A detailed analysis of the structure of a shock wave propagating in a weakly ionized gas is presented, together with a discussion of the physics underlying the key features of the shock structure. A model for the flow behind a shock wave propagating through a weakly ionized gas is developed and used to analyze the effect of the ionization on the aerodynamics and performance of a two-dimensional hypersonic lifting body.

Saeks, R.↗
NASA NTRS1985

Acceleration of cosmic rays in supernova-remnants

It is commonly accepted that supernova-explosions are the dominant source of cosmic rays up to an energy of 10 to the 14th power eV/nucleon. Moreover, these high energy particles provide a major contribution to the energy density of the interstellar medium (ISM) and should therefore be included in calculations of interstellar dynamic phenomena. For the following the first order Fermi mechanism in shock waves are considered to be the main acceleration mechanism. The influence of this process is twofold; first, if the process is efficient (and in fact this is the cas) it will modify the dynamics and evolution of a supernova-remnant (SNR), and secondly, the existence of a significant high energy component changes the overall picture of the ISM. The complexity of the underlying physics prevented detailed investigations of the full non-linear selfconsistent problem. For example, in the context of the energy balance of the ISM it has not been investigated how much energy of a SN-explosion can be transfered to cosmic rays in a time-dependent selfconsistent model. Nevertheless, a lot of progress was made on many aspects of the acceleration mechanism.

Dorfi, E. A.↗
NASA NTRS2002

Turbomachinery Course

In this viewgraph presentation, information and diagrams are provided on rocket engine turbopumps. These turbomachines are highly complex and have several unique features: (1) They are generally very high power density machines; (2) They experience high fluid dynamic loads; (3) They are exposed to severe thermal shocks in terms of rapid starts and stops and extremely high heat transfer coefficients; (4) They have stringent suction performance requirements to minimize tank weight; (5) Their working fluids significantly impact the design: oxidizers are generally explosive, they afford almost no lubrication for bearings and seals, some fuels can degrade material properties, cryogenics result in severe thermal gradients; (6) Their life requirements are short relative to other turbomachines in that there are hundreds of cycles and a few hours of operation for reusable systems. The design of rocket engine turbomachines is a systems engineering challenge because multiple engineering disciplines must be integrated to deal with issues pertaining to stress, structural dynamics, hydrodynamics, aerodynamics, thermodynamics, and materials and process selection.

Stinson, Henry↗
NASA NTRS2003

Finite Element Method for Capturing Ultra-relativistic Shocks

While finite element methods are used extensively by researchers solving computational fluid dynamics in fields other than astrophysics, their use in astrophysical fluid simulations has been predominantly overlooked. Current simulations using other methods such as finite difference and finite volume (based on finite difference) have shown remarkable results, but these methods are limited by their fundamental properties in aspects that are important for simulations with complex geometries and widely varying spatial and temporal scale differences. We have explored the use of finite element methods for astrophysical fluids in order to establish the validity of using such methods in astrophysical environments. We present our numerical technique applied to solving ultra-relativistic (Lorentz Factor Gamma >> 1) shocks which are prevalent in astrophysical studies including relativistic jets and gamma-ray burst studies. We show our finite element formulation applied to simulations where the Lorentz factor ranges up to 2236 and demonstrate its stability in solving ultra-relativistic flows. Our numerical method is based on the Flowfield Dependent Variation (FDV) Method, unique in that numerical diffusion is derived from physical parameters rather than traditional artificial viscosity methods. Numerical instabilities account for most of the difficulties when capturing shocks in this regime. Our method results in stable solutions and accurate results as compared with other methods.

Richardson, G. A.↗
NASA NTRS1993

Numerical study of shock-induced combustion in methane-air mixtures

The shock-induced combustion of methane-air mixtures in hypersonic flows is investigated using a new reaction mechanism consisting of 19 reacting species and 52 elementary reactions. This reduced model is derived from a full kinetic mechanism via the Detailed Reduction technique. Zero-dimensional computations of several shock-tube experiments are presented first. The reaction mechanism is then combined with a fully implicit Navier-Stokes computational fluid dynamics (CFD) code to conduct numerical simulations of two-dimensional and axisymmetric shock-induced combustion experiments of stoichiometric methane-air mixtures at a Mach number of M = 6.61. Applications to the ram accelerator concept are also presented.

Yungster, Shaye↗
NASA NTRS2013

Understanding the Flow Physics of Shock Boundary-Layer Interactions Using CFD and Numerical Analyses

Computational fluid dynamic (CFD) analyses of the University of Michigan (UM) Shock/Boundary-Layer Interaction (SBLI) experiments were performed as an extension of the CFD SBLI Workshop held at the 48th AIAA Aerospace Sciences Meeting in 2010. In particular, the UM Mach 2.75 Glass Tunnel with a semi-spanning 7.75deg wedge was analyzed in attempts to explore key physics pertinent to SBLI's, including thermodynamic and viscous boundary conditions as well as turbulence modeling. Most of the analyses were 3D CFD simulations using the OVERFLOW flow solver, with additional quasi-1D simulations performed with an in house MATLAB code interfacing with the NIST REFPROP code to explore perfect verses non-ideal air. A fundamental exploration pertaining to the effects of particle image velocimetry (PIV) on post-processing data is also shown. Results from the CFD simulations showed an improvement in agreement with experimental data with key contributions including adding a laminar zone upstream of the wedge and the necessity of mimicking PIV particle lag for comparisons. Results from the quasi-1D simulation showed that there was little difference between perfect and non-ideal air for the configuration presented.

Friedlander, David J.↗
NASA NTRS1980

'GAIM' - Gas-addition, impedance-matched arc driver

A conceptual view for a GAIM energy/driver system to maximize shock-tube performance through efficient interfacing of the energy source with the gas dynamics of the arc driver is presented. Electrical and arc-chamber requirements are evaluated utilizing two new computer codes. One code calculates the shock wave generated for a selected time rate and magnitude of arc-energy input; the other computes the values of external circuit elements required to produce the selected energy input, with the driver represented as the load element of the electrical discharge circuit. Results indicate that the energy-storage capability and the driver arrangement needed to produce the highest shock Mach number can be achieved by means of driver gas addition and by impedance matching (GAIM). Design criteria are presented for arc energy requirements necessary to produce given shock-wave speeds. Shock velocities as high as the 70 km/sec required for simulating Jovian entry now seem possible in shock-tube operation. Practical implementation of a GAIM system is discussed.

Dannenberg, R. E.↗
NASA NTRS1996

An MHD Study of the Interaction Between the Solar Wind and the Interstellar Medium

The overall objective of this research program is to obtain a better understanding of the interaction between the solar wind and the interstellar medium through the use of numerical solutions of the time-dependent magnetohydrodynamic (MHD) equations. The simulated results have been compared with observations where possible and with the results from previous analytic and numerical studies. The primary accomplishment of this project has been the development of codes for 2-D models in both spherical and cylindrical coordinates and the application of the codes to the solar wind/interstellar medium interaction. Computations have been carried out for both a relatively simple gas-dynamic interaction and a flow-aligned interstellar magnetic field. The results have been shown to compare favorably with models that use more approximations and to modify and extend the previous results as would be expected. The simulations have also been used along with a data analysis study to provide a quantitative estimate of the distance to the termination and bow shocks. Some of the specific topics that have been studied are: (1) gas dynamic models of the solar wind/interstellar medium interaction, (2) termination shock response to large-scale solar wind fluctuations, and (3) distances to the termination shock and heliopause. The main results from each of these studies are summarized. The results were published in three papers which are included as attachments.

Steinolfson, R. S.↗
NASA NTRS1979

Position and shape of the Venus bow shock - Pioneer Venus Orbiter observations

Magnetometer data from the Pioneer Venus Orbiter is used to examine the position and shape of this planet's bow shock. Utilizing crossings identified on 86 occasions during the first 65 orbits a mean shock surface is defined for sun-Venus-satellite angles of 60-110 deg. Both the shock shape and variance in location are found to be very similar to the terrestrial case for the range in SVS angle considered. However, while the spread in shock positions at the earth is due predominantly to the magnetopause location varying in response to solar wind dynamic pressure, ionopause altitude variations can have little effect on total obstacle radius. Thus, the Cytherean shock is sometimes observed much closer to or farther from the planet than previously predicted by gasdynamic theory applied to the deflection of flow about a blunt body which acts neither as source nor sink for any portion of the flow.

Slavin, J. A.↗
NASA NTRS1987

Increase of ion kinetic temperature across a collisionless shock. II - A simulation study

The present paper is a continuation of the preceding article by Lee et al., (1986) in which it is suggested that the nonadiabatic motion of the directly transmitted ions in a quasi-perpendicular shock wave can result in an increase of the ion kinetic temperature transverse to the ambient magnetic field in the downstream. A series of computer simulations based on a hybrid code have been carried out to examine the dynamics of the transmitted ions in both the subcritical and supercritical shock waves. It is found that, in both cases, the directly transmitted ions can contribute to the heating process. In the case of a resistiveless supercritical shock, the reflected and transmitted ions can be equally important; whereas for a subcritical shock, the transmitted ions are primarily responsible for the ion heating.

Lee, L. C.↗