Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic shock loading”

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 199 records · Page 11

30 cm Drop Modeling

This report documents modeling and analysis of the 30 cm cask drop scenario for spent nuclear fuel (SNF) casks. The 30 cm drop is an important hypothetical dynamic loading scenario for SNF because 10 CFR 71.71 defines testing requirements for normal conditions of transport (NCT) of SNF packages, and free drop (of 30 cm for the heaviest packages) tests are required. Other NCT conditions specified in 10 CFR 71.71 include vibration normally incident to transport, and DOE’s multimodal transportation test (MMTT) and related analysis concluded that mechanical shock and vibration incident to transportation is negligible for SNF structural integrity because the loads are so small. The SNF cask 30 cm drop has a lot more potential to cause damage to SNF because the kinetic energy and momentum of the system is stopped abruptly during the cask impact event. Test data and model predictions both agree that the SNF assemblies are subjected to significantly higher loads in a cask drop event than normal transportation shock and vibration, but currently there is no reason to expect gross structural damage of the fuel assembly or widespread SNF cladding failure as a result of a 30 cm cask drop. The ultimate goal of this work is to close the knowledge gap of stress profiles on SNF rods and fuel assembly components during NCT, and this report describes the progress toward reaching that goal as well as a plan for finishing this work over the next year or two. Key remaining tasks include refining the fuel assembly model, calculating a selection of specific fuel assembly cases to populate the damage model, and finally running the damage model to predict the full general response range of 30 cm drop loads on SNF.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

'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.↗

Anisotropic strength behavior of single-crystal TATB

High-rate strength behavior plays an important role in the shock initiation of high explosives, with plastic deformation serving to localize heat into hot spots and as a mechanochemical means to enhance reactivity. Recent simulations predict that detonation-like shocks produce highly reactive nanoscale shear bands in the layered crystalline explosive TATB (1,3,5-triamino-2,4,6-trinitrobenzene), but the thresholds leading to this response are poorly understood. We utilize molecular dynamics to simulate the high-rate compressive stress–strain response of TATB, with a focus on understanding flow behavior. Here, the dependence of strength on pressure and loading axis (crystal orientation) is explored. The deformation mechanisms fall broadly into two categories, with compression along crystal layers activating a buckling/twinning mode and compression normal to the layers producing nanoscale shear bands. Despite the complexity of the underlying mechanisms, the crystal exhibits relatively straightforward stress–strain curves. Most of the crystal orientations studied show rapid strain softening following the onset of yielding, which settles to a steady flow state. Trajectories are analyzed using five metrics for local states and structural order, but most of these metrics yield similar distributions for these deformation mechanisms. On the other hand, a recently proposed measure of intramolecular strain energy is found to most cleanly distinguish between these mechanisms, while also providing a plausible connection with mechanochemically accelerated decomposition kinetics. Localization of intramolecular strain energy is found to depend strongly on crystal orientation and pressure.

36 MATERIALS SCIENCE↗

Ultrafast x-ray diffraction of high-pressure phases in dynamically compressed TiO 2

Here, we investigate the high-pressure polymorphism of Ti⁢O 2 under laser-shock compression from 54(5) to 137(7) GPa using in situ femtosecond x-ray diffraction. Our results provide experimental evidence of the 𝑃⁢𝑐⁢𝑎⁢2 1 -type distorted fluorite structure formed from polycrystalline Ti⁢O 2 dynamically compressed to 54(5) GPa. Upon higher compression, we observe the direct formation of the ninefold coordinated F⁡e 2 ⁢P-type phase at 68(4) and 78(3) GPa in polycrystalline and [001]-oriented Ti⁢O 2 , respectively. This represents an unprecedented 100 GPa reduction in the shock synthesis pressure of the F⁡e 2⁢ P-type structure relative to quasihydrostatic loading conditions. On pressure release, the F⁡e 2⁢ P-type phase transforms to the α-Pb⁢O 2 structure and, at later times, reverts to rutile. Thus, the rutile →F⁡e 2⁢ P and α-Pb⁢O 2 →rutile transformations are both observed to occur on nanosecond timescales. Our results highlight the unique ability of high-strain-rate uniaxial compression to synthesize novel high-pressure phases and also indicate the importance of in situ atomic-level probes in developing pressure-temperature phase diagrams.

Crystal structures↗

Analytic model of dislocation density evolution in fcc polycrystals accounting for dislocation generation, storage, and dynamic recovery mechanisms

Here, an analytic model of the evolution of dislocation density in fcc polycrystals is described. The evolution equations approximately account for most known dislocation storage, dynamic recovery, and dislocation generation mechanisms in fcc polycrystals. Specifically, the model incorporates network (forest) and grain boundary storage, mobile-network and mobile–mobile annihilation, screw–screw annihilation via athermal and thermal single cross-slip, generation by double cross-slip (Koehler mechanism, including dipole formation), Frank-Read sources, grain boundary nucleation, and mobile–immobile dislocation nucleation due to shock loading. Single cross-slip is assumed to proceed through the Friedel–Escaig (FE) mechanism; the corresponding activation energy is calculated using a modified FE model. The activation energy for double cross-slip is calculated for the first time by extending the FE model. The exact evolution equations are integro-differential equations, and as such are difficult to implement in a code; hence, the evolution equations are simplified by making several approximations. Preliminary results on copper are presented, including comparisons to experimental data.

36 MATERIALS SCIENCE↗

Sierra/SD – Verification Test Manual – 5.22

Verification and validation (V&V) of scientific computing programs are important at Sandia National Labs due to the expanding role of computational simulation in managing the United States nuclear stockpile. The complexities of structural response calculations used to analyze physical problems, the varieties of codes applied to the calculations, and the importance of accurate predictions when assessing field conditions demand confidence in the consistency and accuracy of computer codes. Confidence in the accuracy of the predictions arising from computer simulations must ultimately be gained through verification and validation. The Sierra salinas structural dynamics analysis code, Sierra/SD, is used at the DOE Laboratories, and in several DOD projects. The roles of Sierra/SD in the qualification of weapon systems and components for normal and hostile environments throughout the Stockpile-to-Target Sequence include to, • Redesign weapon components. • Certify weapon components and systems for target environments such as hypersonic vehicles. • Certify that components will survive the thermal mechanical shock loads associated with hostile environments. • Evaluate current stockpile issues, including issues associated with uncertainty quantification. • Address many other problems that are encountered in stockpile management. The Sierra/SD verification plan is described, and an evolving set of key verification tests are described in detail. The verification tests ensure the correctness of the mathematics and numerical algorithms associated with functionality describing engineering phenomena. Development is in accordance with a set of tailored Software Quality Engineering (SQE) practices. SQE practices guide the overall verification and validation effort.

97 MATHEMATICS AND COMPUTING↗

Fluid dynamics of the unsteady two phase processes leading to DDT in granular solid propellants

Deflagration to Detonation (DDT) was predicted to occur in porous beds of high-energy solid propellants by solving the unsteady fluid mechanical convective heat transfer from hot gas products, obtained from the rapid burning at high pressures, provides the impetus to develop a narrow combustion zone and a resulting strong shock. A parametric study clearly indicates that DDT occurs only when a combination of the solids loading fraction, the burning rate constants, the propellant chemical energy, and the particle size provide for critical energy and gas release to support a detonation wave. Predictions for the run-up length to detonation as a function of these parameters are presented.

Krier, H.↗

Landing-Impact Characteristics of Load-Alleviating Struts on a Model of a Winged Space Vehicle

A brief experimental investigation was made of the landing-impact characteristics of a 1/9-scale dynamic model of a winged space vehicle. The landing tests were made by catapulting a free model onto a hard; surface runway and onto water. The model had a conical fuselage and a flat - plate wing with a basic delta planform and 75 deg sweepback of the leading edge. The use of yielding-metal shock absorbers and various landing-gear arrangements was investigated during landing impact. The basic landing gear consisted of a dual rubber-tired nose wheel and twin main skids aft of the center of gravity near the wing tips. landing motion and acceleration data were obtained over a range of landing attitudes, gross weights, and initial sinking speeds. Brief tests were made with an alternate nose-wheel location. An all-skid configuration also was briefly evaluated for hard-surface and water landings. The landing gear employing yielding struts for impact-energy absorption during hard-surface landings resulted in accelerations of approximately 5 1/2 g near the nose gear over a range of landing parameters. Replacing the nose wheel and tire with a skid did not significantly change the accelerations. Landings in smooth water with rigid struts and adequate planing area at the nose skid resulted in a maximum landing acceleration of approximately 4g.

Blanchard, Ulysse J.↗

Multiscale Reactive Model for 1,3,5-Triamino-2,4,6-trinitrobenzene Inferred by Reactive MD Simulations and Unsupervised Learning

When high-energy-density materials are subjected to thermal or mechanical insults at extreme conditions (shock loading), a coupled response between the thermo-mechanical and chemical behaviors is systematically induced. Herein we develop a reaction model for the fast chemistry of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) at the mesoscopic scale, where the chemical behavior is determined by underlying microscopic reactive simulations. The slow carbon cluster formation is not discussed in the present work. All-atom reactive molecular dynamics (MD) simulations are performed with the ReaxFF potential, and a reduced-order chemical kinetics model for TATB is fitted to isothermal and adiabatic simulations of single crystal chemical decomposition. Unsupervised machine learning techniques based on non-negative matrix factorization are applied to MD trajectories to model the decomposition kinetics of TATB in terms of a four-component model. The associated heats of reaction are fit to the temperature evolution from adiabatic decomposition trajectories. Using a chemical species analysis, we show that non-negative matrix factorization captures the main chemical decomposition steps of TATB and provides an accurate estimation of their evolution with temperature. The final analytical formulation, coupled to a diffusion term, is incorporated into a continuum formalism, and simulation results are compared one-to-one against MD simulations of 1D reaction propagation along different crystallographic directions and with different initial temperatures. A good agreement is found for both the temporal and spatial evolution of the temperature field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stress field measurements using quantitative schlieren

Quantitative schlieren analysis is extended here to optically transparent solids in quasi-static and dynamic experiments to measure stress distributions. The quasi-static experiments in polymethyl methacrylate (PMMA) compared refraction angles and stress gradients calculated from schlieren images to the analytical Flamant solution of a line load on a half-space. The quantitative schlieren measurements of the stress field in the thin sample with a load compared well to the analytical solution. The analysis method was then extended to explosive induced shock waves in PMMA. The explosive induced response of PMMA was experimentally studied using high-speed schlieren to visualize the shock propagation in conjunction with Photon Doppler Velocimetry (PDV) to record surface velocity histories. The stress state estimated from the schlieren images was compared to the stress calculated from the PDV measurements. High-speed imaging limitations caused the shock wave to not be fully resolved in the images, but was resolved in the PDV measurement. The stress state behind the shock calculated from the high-speed images followed a similar trend to the stress calculated from the PDV measurements.

Physics↗

A platform for planar dynamic compression of crystalline hydrogen toward the terapascal regime

Here we describe a method for laser-driven planar compression of crystalline hydrogen that starts with a sample of solid para-hydrogen (even-valued rotational quantum number j) having an entropy of 0.06 k B /molecule at 10 K and 2 atm, with Boltzmann constant k B . Starting with this low-entropy state, the sample is compressed using a small initial shock (<0.2 GPa), followed by a pressure ramp that approaches isentropic loading as the sample is taken to hundreds of GPa. Planar loading allows for quantitative compression measurements; the objective of our low-entropy compression is to keep the sample cold enough to characterize crystalline hydrogen toward the terapascal range.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

X-ray imaging of liquid-liquid Mg cavitation

The coupling of dynamic loading platforms with synchrotron light sources have enabled experiments examining void formation in optically opaque materials. Much like spallation in a solid, cavitation occurs when rarefaction waves collide within a liquid. The resulting deformation leads to a complex flow field with numerous voids. It has been suggested that compaction of such a structure may be a significant source of ejecta (i.e. shallow bubble collapse). Experiments imaging supported and unsupported shocks interacting with a free surface in shock melted Mg with synchrotron radiation are presented. Here, the results are used to examine the structure of a cavitation plane that forms in the unsupported case and provide an estimate of the porosity within this region.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Unsteady aerodynamic analysis of space shuttle vehicles. Part 1: Summary report

An analysis of the unsteady aerodynamics of space shuttle vehicles was performed. The results show that slender wing theory can be modified to give the potential flow static and dynamic characteristics over a large Mach number range from M = 0 to M 1. A semi-empirical analytic approximation is derived for the loads induced by the leading edge vortex; and it is shown that the developed analytic technique gives good prediction of experimentally determined steady and unsteady delta wing aerodynamics, including the effects of leading edge roundness. At supersonic speeds, attached leading edge flow is established and shock-induced flow separation effects become of concern. Analysis of experimental results for a variety of boost configurations led to a definition of the main features of the flow interference effects between orbiter (delta wing) and booster. The effects of control deflection on the unsteady aerodynamics of the delta-wing orbiter were also evaluated.

Ericsson, L. E.↗

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Unsteady Analyses and Risk Assessment - Part 2

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Design Support of the feasibility of operating conditions and procedures is critical in such cases due to the possibility of startup/shutdown transients, moving shock structures, unsteady shock-boundary layer interactions and engine and diffuser unstart modes that can result in catastrophic failure. Analyses of such systems is difficult due to resolution requirements needed to accurately capture moving shock structures, shock-boundary layer interactions, two-phase flow regimes and engine unstart modes. In a companion paper, we will demonstrate with the use of CFD, steady analyses advanced capability to evaluate supersonic diffuser and steam ejector performance in the sub-scale A-3 facility. In this paper we will address transient issues with the operation of the facility especially at startup and shutdown, and assess risks related to afterburning due to the interaction of a fuel rich plume with oxygen that is a by-product of the steam ejectors. The primary areas that will be addressed in this paper are: (1) analyses of unstart modes due to flow transients especially during startup/ignition, (2) engine safety during the shutdown process (3) interaction of steam ejectors with the primary plume i.e. flow transients as well as probability of afterburning. In this abstract we discuss unsteady analyses of the engine shutdown process. However, the final paper will include analyses of a staged startup, drawdown of the engine test cell pressure, and risk assessment of potential afterburning in the facility. Unsteady simulations have been carried out to study the engine shutdown process in the facility and understand the physics behind the interactions between the steam ejectors, the test cell and the supersonic diffuser. As a first approximation, to understand the dominant unsteady mechanisms in the engine test cell and the supersonic diffuser, the turning duct in the facility was removed. As the engine loses power a rarefaction wave travels downstream that disrupts the shock cell structure in the supersonic diffuser. Flow from the test cell is seen to expand into the supersonic diffuser section and re-pressurizes the area around the nozzle along with a upstream traveling compression wave that emanates from near the first stage ejectors. Flow from the first stage ejector expands to the center of the duct and a new shock train is formed between the first and second stage ejectors. Both stage ejectors keep the facility pressurized and prevent any large amplitude pressure fluctuations from affecting the engine nozzle. The resultant pressure loads the nozzle experiences in the shutdown process are small.

Ahuja, Vineet↗

Understanding the phase transformation mechanisms that affect the dynamic response of Fe-based microstructures at the atomic scales

Large-scale molecular dynamics (MD) simulations were carried out to investigate the shock-induced evolution of microstructure in Fe-based systems comprising single-crystal and layered Cu/Fe alloys with a distribution of interfaces. The shock compression of pure single-crystal Fe oriented along [110] above a threshold pressure results in a BCC (α)→HCP (ε) phase transformation behavior that generates a distribution of ε phase variants in the phase transformed region of the microstructure behind the shock front. The propagation of the release wave through a phase transformed ε phase causes a reverse ε→α phase transformation and renders a distribution of twins for the [110] oriented Fe that serve as void nucleation sites during spall failure. The simulations reveal that the α→ε→α transformation-induced twinning for shock loading along the [110] direction is due to a dominant ε phase variant formed during compression that rotates on the arrival of the release wave followed by a reverse phase transformation to twins in the α phase. The modifications in the evolution of the ε phase variants and twins in Fe behavior are also studied for Cu–Fe layered microstructures due to the shock wave interactions with the Cu/Fe interfaces using a newly constructed Cu–Fe alloy potential. Here, the MD simulations suggest that interfaces affect the observed variants during shock compression and, hence, distributions of twins during shock release that affects the void nucleation stresses in the Fe phase of Cu/Fe microstructures.

36 MATERIALS SCIENCE↗

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dynamic and Flight Tests on Rubber-Cord and Oleo-Rubber-Disk Landing Gears for an F6C-4 Airplane

The investigation described in this report was conducted for the purpose of comparing an oleo-rubber-disk and a rubber-cord landing gear, built for use on an F6C-4 airplane. The investigation consisted of drop tests under various loading conditions and flight tests on an F6C-4 airplane. In the drop tests the total work done on each gear and the work done on each of the shock-absorbing units were determined. For both drop tests and flight tests the maximum loads and accelerations were determined. The comparative results showed that the oleo gear was slightly superior in reducing the ordinary landing shocks, that it had a greater capacity for work, and that it was very superior in the reduction of the rebound. The results further showed that for drops comparable to very severe landings, the rubber-cord gear was potentially more effective as a shock-reducing mechanism. However, due to the construction of this chassis, which limited the maximum elongation of the cords, this gear was incapable of withstanding as severe tests as the oleo gear. The action of the oleo gear was greatly inferior to the action of an ideal gear. The maximum accelerations encountered during the flight tests for severe landings were 3.64g for the rubber-cord gear and 2.27g for the oleo gear. These were less than those experienced in free drops of 7 inches on either gear.

Peck, William C↗

Analytical and experimental investigations of the oblique detonation wave engine concept

Wave combustors, which include the oblique detonation wave engine (ODWE), are attractive propulsion concepts for hypersonic flight. These engines utilize oblique shock or detonation waves to rapidly mix, ignite, and combust the air-fuel mixture in thin zones in the combustion chamber. Benefits of these combustion systems include shorter and lighter engines which require less cooling and can provide thrust at higher Mach numbers than conventional scramjets. The wave combustor's ability to operate at lower combustor inlet pressures may allow the vehicle to operate at lower dynamic pressures which could lessen the heating loads on the airframe. The research program at NASA-Ames includes analytical studies of the ODWE combustor using Computational Fluid Dynamics (CFD) codes which fully couple finite rate chemistry with fluid dynamics. In addition, experimental proof-of-concept studies are being performed in an arc heated hypersonic wind tunnel. Several fuel injection design were studied analytically and experimentally. In-stream strut fuel injectors were chosen to provide good mixing with minimal stagnation pressure losses. Measurements of flow field properties behind the oblique wave are compared to analytical predictions.

Menees, Gene P.↗