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 1,315 records · Page 73

Modeling and Flight Performance of Supersonic Disk-Gap-Band Parachutes in Slender Body Wakes

NASA's ASPIRE (Advanced Supersonic Parachute Inflation Research and Experiments) project is investigating the supersonic deployment and inflation of full-scale parachutes. To achieve Mars-relevant conditions, the parachutes are deployed at high altitudes over Earth on a sounding rocket platform. During the flight test, Disk-Gap-Band parachutes of 21.5 meter diameter are deployed behind a slender payload 1/6th the diameter of the blunt Mars2020 capsule. Due to the differences in leading body geometry between the test flight and a parachute deployment at Mars, high fidelity numerical simulations of slender and blunt bodywakes, and of rigid parachutes behind them, were used to understand differences and similarities in the flow and the effect on parachute drag. The slender body wake is thinner, closes earlier, and presents a smaller wake deficit. Thus, a parachute deployed in the wake of a slender body is more likely to see a higher dynamic pressure than a parachute deployed behind a blunt body. In the presence of a parachute, the interaction of the unsteady wake with the parachute bow shock is stronger behind the blunt body. Simulations yield highly unsteady forces on the parachute, which was modeled as a rigid body. The mean parachute force behind a slender body is between 3 and 12 percent higher than behind a blunt body, depending on the angle of the parachute with the flow. As the angle of incidence increases, more of the parachute moves out of the leading body wakes, decreasing the sensitivity to leading body shape. To compare the flow past parachutes in Earth's and Mars' atmospheres, simulations were also performed in CO2. At the Mach number considered (1.75), the shock standoff distance ahead of the parachute, post-shock jump conditions, and the resulting parachute forces were found to be very similar in both air and CO2, indicating that a high altitude test is a good proxy for a Mars descent. The results of these numerical simulations and available data on past flight and wind tunnel tests of supersonic Disk-Gap-Band parachutes behind slender bodies were used to generate a parachute drag model for ASPIRE, which in turn was used to help design the flight test. The first flight test occurred in October 2017. The parachute was successfully deployed at Mach 1.77 and an altitude of 42 kilometers. Test instrumentation provided the atmospheric conditions, test vehicle trajectory, and the loads on the parachute along with detailed high-resolution imagery of the inflation process. Reconstruction of the flight test indicated that the measured forces on the parachute were within the model's bounds, although the model over-predicted the parachute force during the first few seconds. The parachute forces during the long subsonic period were well-predicted by the ASPIRE drag model.

Muppidi, Suman↗

Investigation of Transitional Shock-Wave/Boundary Layer Interactions Using Direct Numerical Simulations

Interaction of a transitional boundary layer with a normal-shock is investigated using unstructured tetrahedral meshes under the numerical framework of the space-time conservation element, solution element (CESE) method. The computations mimic recent experimental efforts at the University of Tennessee Space Institute, where a Mach 2.0 flow interacts with a tall cylinder attached to a flat plate. The location of the cylinder with respect to the flat plate leading edge determines if the incoming boundary layer is laminar, transitional or fully turbulent. Four representative flow conditions exemplifying laminar and transitional boundary layers are analyzed by direct numerical simulations. Similar to what was observed in the experiments for the case of transitional interaction, the computations reveal an intermittent upstream influence (UI) shock that repeatedly travels upstream from the lambda-foot toward the leading edge before vanishing. Through detailed unsteady flow analysis obtained using Fourier analysis and dynamic mode decomposition techniques, the presence of disturbances with similar frequencies as those measured in experiments were identified in the flow along with locations that appear to influence the dynamics of the flow.

Venkatachari, Balaji Shankar↗

Design of Hypervelocity Flow Generator (HFG) and Its Diagnostics

Ground facilities for hypersonic research are a key element for successful development of aerodynamically proven hypervelocity vehicles. Design concepts and diagnostics of a hypervelocity flow generator (HFG) were made as a test platform for hypersonic/hypervelocity spacecraft models at the NASA Langley Research Center. The HFG is a hypersonic flow field generator using optically heated gas which is blown into an 80 cu m vacuum chamber. The vacuum chamber is kept at a stable vacuum pressure with a combination of three large vacuum pumps, while the HFG is in the test mode. The HFG provides a relatively small test section with approximately a 20 cm window. This facility was designed to generate 2.45 km/sec of flow speed, and potentially generate a continuous flow with the nozzle and vacuum system. The window in test section provides a direct view of the shock wave around a model in order to measure temperature, pressure, and density profile within the shock layer. One of the key test goals under this project is to understand why the emission spectra from the standing shockwave plasma predicted by Lora-Loran codes are significantly different from the measured emission spectra from the Flight Investigation Reentry Environment (FIRE II) Flight. The correct estimate of the thermal loading on the leading edge of hypersonic vehicles greatly affects the aerodynamic design, the material selection for the vehicle, and the cooling requirement and can be obtained by the precise modeling of emission spectra from shockwave plasma. However, the estimation of thermal loading is not an easy task due to complex non-equilibrium radiative process within high temperature shock layers that still falls into a category of cold plasma. Direct flight experiments are the most desirable, but not a cost-effective approach. Analysis by computational fluid dynamics (CFD) offers many test flexibilities. However, the CFD codes must be fully tested and validated with experimental data before the codes are effectively used for practical design. Large discrepancies between experiments and codes appear in hypersonic/hypervelocity flow regimes at high altitudes of 60 km ~ 90 km. This HFG facility offers some important parameters for CFD code validation, such as collision cross-sections, relaxation times, reaction rate coefficients and transportation coefficients. The HFG test facility is based on the ejection flow of high temperature gas heated over to 3500 K through a nozzle. The tungsten gas chamber of the HFG is heated up to a desired temperature by a 60 kW optical power beam source. This system consists of an optical power source, a thermal chamber, an expansion nozzle, a test section, and an 80 cu m vacuum tank. 60 kW optical input power is obtained from the 150 kW Vortek arc lamp system (by Vortek Industries, Vancouver, Canada). This optical beam is focused to heat the gas chamber within which a flow media is heated. The maximum achievable temperature of the flow medium reached approximately 3500 K or even higher but is limited by the melting point of the chamber material used. The exhaust velocity through the nozzle was determined by the stagnation temperature and the molecular weight of the working medium at the test section. To provide design parameters, a NASA Chemical Equilibrium with Application (CEA) computer program is used for the simulation of aerothermal data. This CEA program can calculate chemical equilibrium and properties of complex mixtures using shock tube parameters. For nitrogen gas at 2666 K stagnation temperature, the maximum achievable velocity at the test section is approximately 2.45 km/sec which is within the range of the thermal velocity of 8000 K shock layer. Based on the calculation through the CEA program, the design parameters of a HFG were determined and implemented for the test section that includes an expansion nozzle. The installed test section of the experimental facility can sustain a condition of a re-entry vehicle from the Space at an altitude of 60-90 km. After installation of the HFG, the system was fully tested and its operational parameters were measured. An 80 cu m vacuum chamber of HFG was set at 1 torr level to keep a stable downstream condition. The pumping time to reach the minimum vacuum pressure (~ 1 torr) at the test chamber from atmospheric pressure was approximately an hour using Kinney (MBV-14000/MB -1600/KT-300) and Stokes (Model 1772 and 412) pumps. Such a setup condition allowed a continued stable operation of the HFG experiment with flowrates through 1-, 3-, or 5-mm diameters of nozzle throat. The flow characteristics of the HFG for various operating conditions were performed using a focused Schlieren method. At 0.8 torr chamber pressure, a barrel shock was observed at the test section. With a sphere obstacle of 12 mm diameter in the flow, a bow shock (~ 2-mm thickness) was observed by a focused Schlieren visualization method.

Hypervelocity flow↗

Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX

A multiscale modeling strategy is used to quantify factors governing the temperature rise in hot spots formed by pore collapse from supported and unsupported shock waves in the high explosive HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Two physical aspects are examined in detail, namely the melting temperature and liquid shear viscosity. All-atom molecular dynamics simulations of phase coexistence are used to predict the pressure-dependent melting temperature up to 5 GPa. Equilibrium simulations and the Green–Kubo formalism are used to obtain the temperature- and pressure-dependent liquid shear viscosity. Starting from a simplified continuum-based grain-scale model of HMX, in this study we systematically increase the complexity of treatments for the solid–liquid phase transition and liquid shear viscosity in simulations of pore collapse. Using a realistic pressure-dependent melting temperature completely suppresses melting for supported shocks, which is otherwise predicted when treating it as a constant determined at atmospheric pressure. Alternatively, melt pools form around collapsed pores when the pressure (and melting temperature) are reduced during the release stage of unsupported shocks. Capturing the pressure dependence of the shear viscosity increases the peak temperature of melt pools by hundreds of Kelvin through viscous work. The complicated interplay of the solid-phase plastic work, solid–liquid phase transition, and liquid-phase viscous work identified here motivate taking a systematic approach to building increasingly complex grain-scale models.

36 MATERIALS SCIENCE↗

Surface Properties of Double-Fin Generated Shock-Wave/Boundary-Layer Interactions

Shock-Wave/Boundary-Layer Interaction (SBLI) is ubiquitous in high-speed air vehicles. SBLI surface flowfield properties of crossing shocks induced by symmetric Double-Fins (DFs) at fin angles of 8° and 10° at Mach 2 are experimentally investigated; a 10° Single-Fin (SF) is also explored as a comparison. Surface visualization is used to capture distinct flowfield structures such as separation and upstream influence. Steady and unsteady Pressure-Sensitive Paints (PSPs) are used to obtain global surface pressure fields. Characteristics of DF-SBLI surface topology are discussed by comparing the corresponding SF-SBLI. Here, the effect of fin angle on DF-SBLI is characterized, and unsteady dynamics of the pressure field are examined. Spectral Proper Orthogonal Decomposition (SPOD) is performed on unsteady PSP, and traveling surface pressure waves and their frequency-dependent behavior are discussed. Dispersion relation of surface pressure waves is revealed by performing two-dimensional space–time Fourier transform on unsteady PSP, and the results are compared with the phase velocity distributions obtained from SPOD. Accompanying group velocity behavior is examined for distinct SBLI regions wherein wave propagation is found dependent on its frequency. Steady surface pressures along the centerline are compared with previous experimental data. Fin-tip spacing is found to be an essential parameter for scaling of upstream influence length.

33 ADVANCED PROPULSION SYSTEMS↗

Towards Understanding the Fluid Dynamic Phenomenon of Interest to Rocket Base Heating: A Review

The significance of the base heating problem for rockets during ascent is due to the complex interaction between the rocket nozzle plumes and the external-flow which can change the flow characteristics in the base region dramatically. At lower altitudes the external-flow merges with the plume-flow, without the formation of a large separated flow region, and the cooler external-flow promotes convective cooling of the base wall. Under these conditions the majority of the base heating is due to radiative heating from the shock heated plume gases. At higher altitudes, however, the process of base heating is not so straightforward. The plume and the base flow expands dramatically and separated flow regions occur in the base area. Hot exhaust gases from the rocket nozzle will be entrained into the separated flow regions and produce a convective component to the base wall heating. Further, if the rocket exhaust-gas contains soot, the soot can increase the emission from the gas and dramatically increase the wall absorption coefficient for radiative heating if it is deposited on the walls . In addition, if the rocket exhaust gas is fuel rich, the fuel can bum in the separated flow regions and further increase the base heating. The base burning phenomenon, and the increased base heating caused by it at higher altitudes, have been observed for the Space Shuttle and Saturn Rocket. Under these conditions, the total heating is significantly higher than the heating without separated flow in the base region, and the increase in heating is directly attributable to the fluid dynamic complexity of the base region. Realistic simulation of the base heating requires that the calculated flow environment reproduce the fluid dynamic flow features accurately. Thus, it will be necessary to introduce into the CFD codes the capability for the flow to respond to the complex vehicle geometry, the effect of turbulence, the ability to accurately reproduce the plume shock/shear layer structures and capture such phenomenon as shock induced base separation and base-burning phenomenon. A survey of experimental, theoretical and computational work that details the fluid dynamics of the base flow environment will be presented in the proposed paper. CFD simulations of rocket base flows using standard CFD codes such as OVERFLOW or GASP will be explored to capture these phenomenon accurately. Merits and limitations of these codes for base flow environment predictions will be explored.

Venkatapathy, E.↗

3D probabilistic fracture mechanics / computational fluid dynamics simulation of a reactor pressure vessel under transient conditions

Reactor pressure vessels (RPVs) are safety-critical light-water-reactor components that, under irradiation, experience long-term material degradation in the form of embrittlement. This can increase their susceptibility to fracture under thermal-shock conditions, which could occur during off-normal transients such as loss-of-coolant accidents (LOCAs). During a LOCA, the most severe conditions for the RPV occur when emergency core cooling water is injected through the cold legs into the water-and-steam-filled RPV. The rapid cooling of the downcomer and internal RPV surface causes decreased temperature and elevated thermally driven tensile stresses in the RPV wall. This, combined with long-term material embrittlement, may cause fracture initiation at pre-existing flaws, challenging the integrity of the RPV. Assessing RPV integrity during transients with a large spatial variation in the coolant temperature requires a modeling approach that considers the effects of spatially varying coolant temperature on the fracture probability of a population of flaws distributed throughout the RPV, accounting for spatially varying embrittlement. Here, the present study addresses this need by demonstrating first-of-its-kind coupling of a high-fidelity 3D computational fluid dynamics code with 3D probabilistic fracture mechanics This was accomplished using representative models of a pressurized-water reactor subjected to small- and medium-break LOCA conditions, both of which can result in large spatial temperature variations. While the observed impact of accounting for 3D effects was minimal under the small-break LOCA this study indicates a significant increase in the probability of fracture initiation under the medium-break LOCA when 3D effects are considered, relative to a spatially uniform cooling scenario.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Probing the lattice structure of dynamically compressed and released single crystal iron through the alpha to epsilon phase transition

In this work, experiments using broadband Laue x-ray diffraction (XRD) were used to examine the lattice structure of dynamically compressed [100]-oriented single crystal iron samples at the Dynamic Compression Sector at the Advanced Photon Source. These experiments used 1 μm thick iron single crystal samples sandwiched between a polyimide ablator and a polycarbonate window. A 100 J, 10 ns duration laser pulse incident on the polyimide ablator was used to shock compress the iron samples to initial stresses greater than 25 GPa, exceeding the ~13 GPa alpha (body-centered-cubic or bcc structure) to epsilon (hexagonal-close-packed or hcp structure) phase transition stress. XRD measurements were performed at various times relative to the shock wave entering the iron sample: early times, <~150 ps while the initial shock waves propagated through the iron; intermediate times, after the iron equilibrated with the ablator and window reaching a plateau stress state (12 or 17 GPa) lasting several nanoseconds; and late times, during uniaxial strain release. The early time measurements show that in <~150 ps, the high-pressure hcp phase is relaxed with a c/a ratio of 1.61, contrary to previous laser shock experiments where a c/a ratio of 1.7 was inferred. In the plateau stress state and partially released states, XRD measurements showed that the hcp structure retained a c/a ratio of 1.61 with no observable changes in the microstructure. Upon stress release at ~1 GPa/ns release rate, the reverse phase transition (hcp to bcc) to the original single crystal orientation (implying a transformation memory effect) was observed to reach completion somewhere between 13 and 11 GPa, indicating little stress hysteresis under rapid uniaxial strain release. A similar memory effect for the reverse hcp to bcc transformation has been previously observed under hydrostatic compression. However, the bcc/hcp orientation relationships differ somewhat between dynamic and static compression experiments, implying that the transformation pathway under uniaxial dynamic strain differs from the Burgers mechanism.

36 MATERIALS SCIENCE↗

Novel Hot-Spot Ignition Designs for Inertial Confinement Fusion with Liquid-Deuterium-Tritium Spheres

A new class of ignition designs is proposed for inertial confinement fusion (ICF) experiments. These designs are based on the hot-spot-ignition approach, but instead of conventional targets that comprise of spherical shells with thin frozen deuterium{tritium (DT) layers, homogeneous-density (solid or liquid) DT spheres are used where the lower-density central region and higher-density shell are created dynamically by appropriately shaping the laser pulse. These offer several advantages, including simplicity in target production and lower sensitivity to both laser imprint and physics uncertainty in shock interaction with the ice-vapor interface. The design evolution starts by launching an ~1-Mbar shock into a homogeneous DT sphere. After bouncing from the center, the reflected shock reaches the outer surface of the sphere and the shocked material starts to expand outward until its pressure drops below the ablation pressure. At this point, an adjustment shock is launched inward by supporting ablation pressure. This shock compresses the ablator and fuel, forming a shell. The shell is then accelerated and compressed by appropriately shaping the drive laser pulse, similar to the conventional thin-shell, hot-spot designs. Here, this paper demonstrates the feasibility of the new concept using hydrodynamic simulations and discusses the advantages and disadvantages of the concept compared with more-traditional ICF designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study of the solar corona using radio and space observations

Optical, radio and X-ray evidence of violent mass motions in the corona has existed for some years but only recently have the form, nature, frequency and implication of the transients become obvious. The observed properties of coronal transients concentrating on the white-light and radio manifestations. The possible mechanisims involved in the radio bursts are discussed. The estimates of various forms of energy are reviewed. It appears that the magnetic energy transported from the Sun by the transient exceeds that of any other form, and that magnetic forces dominate in the dynamics of the motions. The conversion of magnetic energy into mechanical energy, by expansion of the fields, provides a possible driving force for the coronal and interplanetary shock waves.

Dulk, G. A.↗

Radio and white-light observations of coronal transients

Optical, radio and X-ray evidence of violent mass motions in the corona has existed for some years but only recently have the form, nature, frequency and implication of the transients become obvious. In this paper the observed properties of coronal transients are reviewed, with concentration on the white-light and radio manifestations. The classification according to speeds seems to be meaningful, with the slow transients having thermal emissions at radio wavelengths and the fast ones nonthermal. The possible mechanisms involved in the radio bursts are then discussed and estimates of various forms of energy are reviewed. It appears that the magnetic energy transported from the sun by the transient exceeds that of any other form, and that magnetic forces dominate in the dynamics of the motions. The conversion of magnetic energy into mechanical energy, by expansion of the field, provides a possible driving force for the coronal and interplanetary shock waves.

Dulk, G. A.↗

Fluid mechanics mechanisms in the stall process of airfoils for helicopters

Phenomena that control the flow during the stall portion of a dynamic stall cycle are analyzed, and their effect on blade motion is outlined. Four mechanisms by which dynamic stall may be initiated are identified: (1) bursting of the separation bubble, (2) flow reversal in the turbulent boundary layer on the airfoil upper surface, (3) shock wave-boundary layer interaction behind the airfoil crest, and (4) acoustic wave propagation below the airfoil. The fluid mechanics that contribute to the identified flow phenomena are summarized, and the usefulness of a model that incorporates the required fluid mechanics mechanisms is discussed.

Young, W. H., Jr.↗

Magnetohydrodynamic interaction of high-speed streams

Numerical solutions of a magnetohydrodynamic model are carried out to describe the nonlinear interaction of corotating high-speed streams near the solar equatorial plane. Two problems are studied. The first problem is to simulate the evolution of an idealized high-speed stream. Numerical solutions are obtained to represent the variations of flow velocity, magnetic field, plasma density, temperature, and conduction heat flux in the interaction region. They demonstrate that the dynamical interaction and heat conduction process are responsible for the thermal structure of a high-speed stream. The second problem deals with the formation of corotating shock waves near the leading edge of a broad stream resulting from the merging of characteristic curves. Corotating shocks do not necessarily occur in pairs; a reverse shock can be formed without a forward shock nearby.

Whang, Y. C.↗

Unsteady turbulent shear flow in shock tube discontinuities

A pressure-ruptured shock tube and an arc driven shock tube, have been used to study the evolution of turbulent fluctuations at contact surfaces with N2O4-2NO2 mixtures and at ionizing shock fronts in argon. The study has focused on point density diagnostics derived from crossed light beam correlations and electric probes. Turbulent bursts are found for which dynamical and spectral analyses suggest a particle-like evolution of fluctuation segments with a unique and characteristic frequency, independent of flow history and overall flow conditions.

Johnson, J. A., III↗

The solar wind interaction with Venus

The relation between Venus and the solar wind is analyzed. The effects of the intrinsic field, neutral atmosphere, and ionosphere of Venus on the solar wind are examined. The solar wind interaction phenomena is studied; consideration is given to the free-stream solar wind, bow shock, magnetosheath, boundary layer, ionospheric features, neutral atmosphere features, and wake and magnetotail. Further research on the boundary layer and tail formation, global models, and high-dynamic pressure interaction is proposed.

Luhmann, J. G.↗

The evolution of the system of globular clusters

The physical processes affecting the evolution of the system of globular clusters are reviewed. Aspects of internal structure evolution discussed include evaporation, core collapse, expansion, tidal shocks, and interaction with other massive components of the Galaxy such as black holes, molecular clouds, and other clusters. The evolution of orbital parameters is examined, including dynamical friction, scattering by massive subcomponents, and interactions with nearby systems. In addition, an overview of the system evolution is presented.

Ostriker, Jeremiah P.↗

Real-time measurement of materials properties at high temperatures by laser produced plasmas

Determination of elemental composition and thermophysical properties of materials at high temperatures, as visualized in the context of containerless materials processing in a microgravity environment, presents a variety of unusual requirements owing to the thermal hazards and interferences from electromagnetic control fields. In addition, such information is intended for process control applications and thus the measurements must be real time in nature. A new technique is described which was developed for real time, in-situ determination of the elemental composition of molten metallic alloys such as specialty steel. The technique is based on time-resolved spectroscopy of a laser produced plasma (LPP) plume resulting from the interaction of a giant laser pulse with a material target. The sensitivity and precision were demonstrated to be comparable to, or better than, the conventional methods of analysis which are applicable only to post-mortem specimens sampled from a molten metal pool. The LPP technique can be applied widely to other materials composition analysis applications. The LPP technique is extremely information rich and therefore provides opportunities for extracting other physical properties in addition to the materials composition. The case in point is that it is possible to determine thermophysical properties of the target materials at high temperatures by monitoring generation and transport of acoustic pulses as well as a number of other fluid-dynamic processes triggered by the LPP event. By manipulation of the scaling properties of the laser-matter interaction, many different kinds of flow events, ranging from shock waves to surface waves to flow induced instabilities, can be generated in a controllable manner. Time-resolved detection of these events can lead to such thermophysical quantities as volume and shear viscosities, thermal conductivity, specific heat, mass density, and others.

Kim, Yong W.↗