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.

SourceResearch connectionFrom yearThrough year

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

At least 1,765 records · Page 98

All sources NASA NTRS DOE OSTI

DOE OSTI2021

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Characteristics of thermospheric infrasound predicted using ray tracing and weakly non-linear waveform analyses

Infrasonic signals refracted by thermal gradients in the rarefied upper atmosphere are modeled using a combination of ray tracing and weak shock theory to develop an understanding of thermospheric infrasound signals produced by energetic, transient sources. Canonical arrival structures in the form of u-wave signatures are identified for returns refracted at lower altitudes within the thermosphere, and possible multi-pathing produced by effective sound speed inflections are investigated to elucidate more complex arrival structures, which are found to be spatially localized. Furthermore, variability in the source characteristics is investigated and it is found that whereas some waveform phase information is lost due to finite amplitude effects, arrival characteristics are strongly dependent on the peak overpressure near the source. Variability in the propagation path is considered using archived atmospheric specifications and implies that despite uncertainties related to the dynamic and sparsely sampled nature of the atmosphere, thermospheric signatures might be useful in estimating the yield for explosive sources. Last, thermospheric arrivals from a failed rocket launch, as well as those from several large chemical explosions, are analyzed and it is found that qualitative trends match those predicted, and analyses here provide additional insight into such signatures.

42 ENGINEERING↗

NASA NTRS1992

Characteristics of hypervelocity impact craters on LDEF experiment S1003 and implications of small particle impacts on reflective surfaces

The Ion Beam Textured and Coated Surfaces Experiment (IBEX) was designated S1003 on the Long Duration Exposure Facility (LDEF) at a location of 98 degrees relative to the ram direction. Thirty-six diverse materials were exposed to the micrometeoroid (and some debris) environment for 5.8 years. Optical property measurements indicated no changes for almost all of the materials except S-13G, Kapton, and Kapton-coated surfaces, and these changes can be explained by other environmental effects. From the predicted micrometeoroid flux of NASA SP-8013, no changes in optical properties of the surfaces due to micrometeoroids were expected. However, there were hypervelocity impacts on the various diverse materials flown on IBEX. The characteristics of these craters were documented using scanning electron microscopy (SEM) and are presented. Interest in placing large solar concentrator/solar dynamic systems in space for power generation has again brought up a concern for maintaining the integrity of the optical properties of highly specular reflecting surfaces in the near-Earth space environment. It has been shown that highly reflective polished metals and thin film coatings degrade when exposed to simulated micrometeoroids in the laboratory. At LeRC, a shock tube was used to simulate the phenomenon of micrometeoroid optical properties of surfaces exposed to this impact were then evaluated. A calibrated sensor, 2000 A Al/stainless steel, was developed to not only detect the small size micrometeoroid environment, but also to evaluate the degradation of the optical properties of thin aluminum films in space. This sensor was flown on LDEF experiment S1003 and also on the OSO 3 and SERT 2 satellites that were launched in 1967 and 1970, respectively. No changes in the optical properties of the highly reflective surface sensor on SERT 2 were measured during 20 years in space. The results, as determined by the accuracy of the sensor, indicate that a highly reflective surface should lose less than 1 percent of its specular reflectance in near-Earth orbit during 20 years.

Mirtich, Michael J.↗

NASA NTRS1993

Uniform high-order spectral methods for one- and two-dimensional Euler equations

Uniform high order spectral methods to solve multi-dimensional Euler equations for gas dynamics are discussed. Uniform high order spectral approximations with spectral accuracy in smooth regions of solutions are constructed by introducing the idea of the Essentially Non-Oscillatory (ENO) polynomial interpolations into the spectral methods. The authors present numerical results for the inviscid Burgers' equation, and for the one-dimensional Euler equations including the interactions between a shock wave and density disturbance, Sod's and Lax's shock tube problems, and the blast wave problem. The interaction between a Mach 3 two-dimensional shock wave and a rotating vortex is simulated.

Cai, Wei↗

NASA NTRS1993

Chromospheric heating by acoustic shocks - A confrontation of GHRS observations of Alpha Tauri (K5 III) with ab initio calculations

We compare ab initio calculations of semiforbidden C II line profiles near 2325 A with recently published observations of the inactive red giant Alpha Tau (K5 III) obtained using the GHRS on board the Hubble Space Telescope. Our one-dimensional, time-dependent calculations assume that the chromosphere is heated by stochastic acoustic shocks generated by photospheric convection. We calculate various models using results from traditional (mixing length) convection zone calculations as input to hydrodynamical models. The semiforbidden C II line profiles and ratios provide sensitive diagnostics of chromospheric velocity fields, electron densities, and temperatures. We identify major differences between observed and computed line profiles which are related to basic gas dynamics and which are probably not due to technical modeling restrictions. If the GHRS observations are representative of chromospheric conditions at all epochs, then one (or more) of our model assumptions must be incorrect. Several possibilities are examined. We predict time variability of semiforbidden C II lines for comparison with observations. Based upon data from the IUE archives, we argue that photospheric motions associated with supergranulation or global pulsation modes are unimportant in heating the chromosphere of Alpha Tau.

Judge, P. G.↗

NASA NTRS2004

Input/output properties of the lateral vestibular nucleus

This article is a review of work in three species, squirrel monkey, cat, and rat studying the inputs and outputs from the lateral vestibular nucleus (LVN). Different electrophysiological shock paradigms were used to determine the synaptic inputs derived from thick to thin diameter vestibular nerve afferents. Angular and linear mechanical stimulations were used to activate and study the combined and individual contribution of inner ear organs and neck afferents. The spatio-temporal properties of LVN neurons in the decerebrated rat were studied in response to dynamic acceleration inputs using sinusoidal linear translation in the horizontal head plane. Outputs were evaluated using antidromic identification techniques and identified LVN neurons were intracellularly injected with biocytin and their morphology studied.

Efferent Pathways/anatomy & histology/physiology↗

NASA NTRS2011

OVERFLOW Validation for Predicting Plume Impingement of Underexpanded Axisymmetric Jets onto Angled Flat Plates

This report documents how OVERFLOW, a computational fluid dynamics code, predicts plume impingement of underexpanded axisymmetric jets onto both perpendicular and inclined flat plates. The effects of the plume impinging on a range of plate inclinations varying from 90deg to 30deg are investigated and compared to the experimental results in Reference 1 and 2. The flow fields are extremely complex due to the interaction between the shock waves from the free jet and those deflected by the plate. Additionally, complex mixing effects create very intricate structures in the flow. The experimental data is very limited, so these validation studies will focus only on cold plume impingement on flat and inclined plates. This validation study will help quantify the error in the OVERFLOW simulation when applied to stage separation scenarios.

Lee, Henry C.↗

NASA NTRS1980

Physicochemical and dynamical processes in cometary ionospheres. I - The basic flow profile

The paper deals with the basic aspects of solar wind interaction with the ionosphere of a large comet and with detailed hydrodynamic models developed to describe the flow exterior to the tangential discontinuity surface that separates the outflowing purely cometary ions from the inflowing contaminated solar wind. The equation of the basic flow model and the results of a preliminary analysis of these equations are given, and the shape and structure of the shock layer are analyzed. The shock layer structure associated with fan-shaped comet emanating from the cometary nucleus in the general direction to the sun is identified.

Houpis, H. L. F.↗

DOE OSTI2023

Intricate structure of the plasma Rayleigh–Taylor instability in shock tubes

Spikes and bubbles grow on unstable interfaces that are accelerated in high-energy-density conditions. If a shock propagates ahead of the interface, the plasma can be heated to extreme conditions where conduction and radiation fluxes influence the hydrodynamics. For example, a National Ignition Facility experiment found reduced single-mode nonlinear mixed-width growth in conditions scaled from a supernova explosion. We present high-resolution two-dimensional radiation hydrodynamic simulations with the Flash code that quantitatively reproduce the experiment. Radiative fluxes are primarily responsible for ablating the spike and removing the mushroom caps. The ablated plasma increases the mixed mass and forms a low-density halo with spikes forming in both directions. This is considerably more complex than classical instability. The halo is sensitive to ablative physics, so radiographing it may aid in the verification of energy transport modeling. Although ablation changes the spike shape, it has little effect on the overall mixed width for these parameters. This is because ablation enhances the bubble velocity but it has the opposite effect on the spike. The radiation transport instead suppresses the growth via increasing the shocked foam density, thus decreasing the Atwood number. Furthermore, a terminal velocity model including the rarefaction expansion agrees with the experimental mixed-width growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRSDate not supplied

Summary of Shock Wave Turbulent Boundary Layer Interaction Experiments In a Circular Test Section

A series of experiments were performed at Mach 2.5 in a 17 cm diameter circular test section to characterize an impinging/reflected shock wave turbulent boundary layer interaction generated by a cone-cylinder centerbody. The cone-cylinder centerbody generates a conical shock wave that interacts with the naturally occurring boundary layer developing on the test section wall. Three different cone angles were used in the experiment to study unseparated, incipiently separated, and separated interactions. When the cone-cylinder centerbody is positioned on the centerline, a flowfield which is two-dimensional in the mean is generated. Three dimensional interactions were also created by offsetting the cone-cylinder centerbody from the test section centerline. The results are intended to provide benchmark quality datasets for computational fluid dynamics (CFD) validation without the pitfalls inherent in rectangular configurations where corner effects prohibit a truly two-dimensional flow in the mean. The experimental measurements included surface flow visualization, wall static pressure, flowfield Pitot tube pressure, constant-voltage anemometry (CVA) normal hot-wire, and particle image velocimetry (PIV) measurements. The hot-wire measurements were used to calculate mean mass flux and total temperature profiles, mass flux and total temperature turbulence intensities, and the mass flux-total temperature correlation. The PIV measurements provide three-dimensional mean velocity measurements. Agreement between the pressure, hot-wire, and PIV measurements is established in the undisturbed upstream flowfield.

Supersonic↗

NASA NTRSDate not supplied

Summary of Shock Wave Turbulent Boundary Layer Interaction Experiments in a Circular Test Section

compressible flow↗

NASA NTRS1989

Dynamic grid refinement for partial differential equations on parallel computers

The fast adaptive composite grid method (FAC) is an algorithm that uses various levels of uniform grids to provide adaptive resolution and fast solution of PDEs. An asynchronous version of FAC, called AFAC, that completely eliminates the bottleneck to parallelism is presented. This paper describes the advantage that this algorithm has in adaptive refinement for moving singularities on multiprocessor computers. This work is applicable to the parallel solution of two- and three-dimensional shock tracking problems.

Mccormick, S.↗

DOE OSTI2024

Recent developments in mathematical aspects of relativistic fluids

Abstract We review some recent developments in mathematical aspects of relativistic fluids. The goal is to provide a quick entry point to some research topics of current interest that is accessible to graduate students and researchers from adjacent fields, as well as to researches working on broader aspects of relativistic fluid dynamics interested in its mathematical formalism. Instead of complete proofs, which can be found in the published literature, here we focus on the proofs’ main ideas and key concepts. After an introduction to the relativistic Euler equations, we cover the following topics: a new wave-transport formulation of the relativistic Euler equations tailored to applications; the problem of shock formation for relativistic Euler; rough (i.e., low-regularity) solutions to the relativistic Euler equations; the relativistic Euler equations with a physical vacuum boundary; relativistic fluids with viscosity. We finish with a discussion of open problems and future directions of research.

Disconzi, Marcelo (ORCID:0000000234497778)↗

DOE OSTI2022

Flame acceleration and transition to detonation in a pre-/main-chamber combustion system

Numerical simulations are performed to study the mechanism of deflagration to detonation transition (DDT) in a pre-/main-chamber combustion system with a stoichiometric ethylene–oxygen mixture. A Godunov algorithm, fifth-order in space, and third-order in time, is used to solve the fully compressible Navier–Stokes equations on a dynamically adapting mesh. A single-step, calibrated chemical diffusive model described by Arrhenius kinetics is used for energy release and conservation between the fuel and the product. The two-dimensional simulation shows that a laminar flame grows in the pre-chamber and then develops into a jet flame as it passes through the orifice. A strong shock forms immediately ahead of the flame, reflecting off the walls and interacting with the flame front. The shock–flame interactions are crucial for the development of flame instabilities, which trigger the subsequent flame development. The DDT arises due to a shock-focusing mechanism, where multiple shocks collide at the flame front. A chemical explosive mode analysis (CEMA) criterion is developed to study the DDT ignition mode. Preliminary one-dimensional computations for a laminar propagating flame, a fast flame deflagration, and a Chapman–Jouguet detonation are conducted to demonstrate the validity of CEMA on the chemical-diffusive model, as well as to determine the proper conditioning value for CEMA diagnostic. The two-dimensional analysis with CEMA indicates that the DDT initiated by the shock-focusing mechanism can form a strong thermal expansion region at the flame front that features large positive eigenvalues for the chemical explosive mode and dominance of the local autoignition mode. Thus, the CEMA criterion proposed in this study provides a robust diagnostic for identifying autoignition-supported DDT, of which the emergence of excessive local autoignition mode is found to be a precursor. The effect of grid size, initial temperature, and orifice size are then evaluated, and results show that although the close-chamber DDT is highly stochastic, the detonation initiation mechanism remains robust.

42 ENGINEERING↗

NASA NTRS1993

Development of a tethered satellite force transducer

The SEDS (Small Expendable Deployer System)/Delta II is a demonstration flight of a tethered satellite deployer system. The satellite or end mass payload, which is deployed and tethered to the Delta II rocket, is instrumented to determine its dynamics during deployment. A three-component force transducer or strain-gage balance was developed to measure tether tension on the end-mass side of the tether. The transducer was designed to measure tensions up to 20N in each direction and resolve tensions as low as 0.002N. In addition, the transducer was required to withstand the shock, vibration, and temperature excursions of the mission. The development of the transducer from design to testing and integration is discussed.

Rhew, Ray D.↗

NASA NTRS1994

Advanced ballistic range technology

The research conducted supported two facilities at NASA Ames Research Center: the Hypervelocity Free-Flight Aerodynamic Facility and the 16-Inch Shock Tunnel. During the grant period, a computerized film-reading system was developed, and five- and six-degree-of-freedom parameter-identification routines were written and successfully implemented. Studies of flow separation were conducted, and methods to extract phase shift information from finite-fringe interferograms were developed. Methods for constructing optical images from Computational Fluid Dynamics solutions were also developed, and these methods were used for one-to-one comparisons of experiment and computations.

Yates, Leslie A.↗

NASA NTRS2001

Nature of Fluctuations on Directional Discontinuities Inside a Solar Ejection: Wind and IMP 8 Observations

A solar ejection passed the Wind spacecraft between December 23 and 26, 1996. On closer examination, we find a sequence of ejecta material, as identified by abnormally low proton temperatures, separated by plasmas with typical solar wind temperatures at 1 AU. Large and abrupt changes in field and plasma properties occurred near the separation boundaries of these regions. At the one boundary we examine here, a series of directional discontinuities was observed. We argue that Alfvenic fluctuations in the immediate vicinity of these discontinuities distort minimum variance normals, introducing uncertainty into the identification of the discontinuities as either rotational or tangential. Carrying out a series of tests on plasma and field data including minimum variance, velocity and magnetic field correlations, and jump conditions, we conclude that the discontinuities are tangential. Furthermore, we find waves superposed on these tangential discontinuities (TDs). The presence of discontinuities allows the existence of both surface waves and ducted body waves. Both probably form in the solar atmosphere where many transverse nonuniformities exist and where theoretically they have been expected. We add to prior speculation that waves on discontinuities may in fact be a common occurrence. In the solar wind, these waves can attain large amplitudes and low frequencies. We argue that such waves can generate dynamical changes at TDs through advection or forced reconnection. The dynamics might so extensively alter the internal structure that the discontinuity would no longer be identified as tangential. Such processes could help explain why the occurrence frequency of TDs observed throughout the solar wind falls off with increasing heliocentric distance. The presence of waves may also alter the nature of the interactions of TDs with the Earth's bow shock in so-called hot flow anomalies.

Vasquez, Bernard J.↗

NASA NTRS2008

Transient Three-Dimensional Startup Side Load Analysis of a Regeneratively Cooled Nozzle

The objective of this effort is to develop a computational methodology to capture the startup side load physics and to anchor the computed aerodynamic side loads with the available data from a regeneratively cooled, high-aspect-ratio nozzle, hot-fired at sea level. The computational methodology is based on an unstructured-grid, pressure-based, reacting flow computational fluid dynamics and heat transfer formulation, a transient 5 s inlet history based on an engine system simulation, and a wall temperature distribution to reflect the effect of regenerative cooling. To understand the effect of regenerative wall cooling, two transient computations were performed using the boundary conditions of adiabatic and cooled walls, respectively. The results show that three types of shock evolution are responsible for side loads: generation of combustion wave; transitions among free-shock separation, restricted-shock separation, and simultaneous free-shock and restricted shock separations; along with the pulsation of shocks across the lip, although the combustion wave is commonly eliminated with the sparklers during actual test. The test measured two side load events: a secondary and lower side load, followed by a primary and peak side load. Results from both wall boundary conditions captured the free-shock separation to restricted-shock separation transition with computed side loads matching the measured secondary side load. For the primary side load, the cooled wall transient produced restricted-shock pulsation across the nozzle lip with peak side load matching that of the test, while the adiabatic wall transient captured shock transitions and free-shock pulsation across the lip with computed peak side load 50% lower than that of the measurement. The computed dominant pulsation frequency of the cooled wall nozzle agrees with that of a separate test, while that of the adiabatic wall nozzle is more than 50% lower than that of the measurement. The computed teepee-like formation and the tangential motion of the shocks during lip pulsation also qualitatively agree with those of test observations. Moreover, a third transient computation was performed with a proportionately shortened 1 s sequence, and lower side loads were obtained with the higher ramp rate.

Wang, Ten-See↗