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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

The solar wind interaction with Mars - Mariner 4, Mars 2, Mars 3, Mars 5, and Phobos 2 observations of bow shock position and shape

An aggregate Mars bow shock data set using Mariner 4, Mars 2, Mars 3, Mars 5, and Phobos 2 observations has been analyzed. The results support the earlier conclusion that the mean distance to the subsolar shock at Mars is nearly 1.5 planetary radii, from which gas dynamic models predict an obstacle altitude of 500 km. The Martian bow shock does not appear to vary significantly in shape or altitude with the phase of the solar cycle. The unusually distant dayside bow shock crossings reported by Mars 2 and 3 also appear in the Phobos 3 observations, suggesting that the dayside obstacle can on rare occasions reach altitudes over 1000 km. The Martian bow shock differs from that of Venus in that its mean altitude is greater, it lacks a strong solar cycle variation, and its location is far more variable, including the occurrence of strong bow shocks over the dayside hemisphere at distances at least as great as the orbit of Phobos 2, i.e., 2.8 Mars radii.

Slavin, J. A.↗
NASA NTRS2010

Dynamic (Vibration) Testing: Design-Certification of Aerospace System

Various types of dynamic testing of structures for certification purposes are described, including vibration, shock and acoustic testing. Modal testing is discussed as it frequently complements dynamic testing and is part of the structural verification/validation process leading up to design certification. Examples of dynamic and modal testing are presented as well as the common practices, procedures and standards employed.

Aggarwal, Pravin K.↗
DOE OSTI2021

Drop impact tester development for spent nuclear fuel vibration integrity study

The accumulated damage of SNF rod is defined as the integrated damage induced by the SNF assembly harmonic vibrations and SNF assembly & cask system contact interactions induced transient shocks, experienced by SNF rods during normal conditions of transportation (NCT). The CIRFT technology has provided means to resolve/understand the SNF rod harmonic vibration reliability issues associated with maximum axial clad tubing stress induced by SNF inertia related dynamic bending loading. However, the SNF rod integrity in association with the transient shock loads during NCT through dynamic contact interactions among SNF assembly cannot be investigated systematically through CIRFT technology alone. Thus, the proposed drop impact tester development become a viable tool for SNF system transient shock effect investigation. Furthermore, the transient shocks induced loading mode is normally perpendicular to the SNF rod axial orientation, i.e., in a transverse orientation to the SNF rod. To proceed with such loading mode investigation, a new device, “SNF Drop Impact Tester” that provides the transverse impact shock load onto a SNF rod through a free-drop projectile was developed in FY20 for SNF System Vibration Reliability Investigation; where the SNF system contact impact intensity can be calibrated accordingly with projectile weight and projectile head geometry. This newly developed research tool can further provide a detailed understanding about the effect of dynamic contact-interaction loadings; and combined harmonic vibration and transient shock loading modes on the fatigue damage evolution of the HBU SNF under NCT.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

The LISE package: solvers for static and time-dependent superfluid local density approximation equations in three dimensions

Nuclear implementation of the density functional theory (DFT) is at present the only microscopic framework applicable to the whole nuclear landscape. The extension of DFT to superfluid systems in the spirit of the Kohn-Sham approach, the superfluid local density approximation (SLDA) and its extension to time-dependent situations, time-dependent superfluid local density ap- proximation (TDSLDA), have been extensively used to describe various static and dynamical problems in nuclear physics, neutron star crust, and cold atom systems. In this paper, we present the codes that solve the static and time-dependent SLDA equations in three-dimensional coordinate space without any symmetry restriction. These codes are fully parallelized with the message passing interface (MPI) library and take advantage of graphic processing units (GPU) for accelerating execution. The dynamic codes have checkpoint/restart capabilities and for initial conditions one can use any generalized Slater determinant type of wave function. The code can describe a large number of physical problems: nuclear fission, collisions of heavy ions, the interaction of quantized vor- tices with nuclei in the nuclear star crust, excitation of superfluid fermion systems by time dependent external fields, quantum shock waves, domain wall generation and propagation, the dynamics of the Anderson-Bogoliubov-Higgs mode, dynamics of fragmented condensates, vortex rings dynamics, generation and dynamics of quantized vortices, their crossing and recombinations and the incipient phases of quantum turbulence.

Jin, Shi↗
DOE OSTI2026

Mechanism of the wurtzite to rocksalt phase transformation in cadmium sulfide single crystals shock compressed along the 𝑐-axis to elastic impact stresses of ∼ 5 GPa

Cadmium sulfide (CdS), which exhibits a wurzite (WZ) to rocksalt (RS) phase transformation at elevated stresses, is an ideal system to address the role of deformation on phase transformation mechanisms and kinetics. CdS has been shown to have very different elastic-inelastic behavior for shock propagation along the WZ 𝑐- and 𝑎-axes, and also exhibits significant differences in the time scale for the phase transformation when shocked along these different axes. As an important first step in examining the role of deformation in shock-induced phase transformations, here we present in situ, time-resolved x-ray diffraction (XRD) measurements in single-event, shock wave experiments on CdS single crystals shocked along the WZ 𝑐-axis to elastic impact stresses near or above ∼ 5 GPa, where a marked increase in the kinetics of this transformation has been reported. The XRD measurements are compared with forward diffraction simulations to evaluate different proposed phase transformation mechanisms and are shown to be consistent with a shearing mechanism that results in the RS [001] direction being aligned with the original WZ 𝑐-axis and the RS [110] direction being aligned with the original WZ 𝑎-axis. In order to reproduce all observed RS diffraction spots the forward diffraction simulations required three different RS crystalline domains, each with a [110] direction aligned with one of the three original WZ 𝑎-axes, and a significant mosaicity, particularly about the RS [001] direction (the direction of shock propagation). In conclusion, this large mosaicity is likely due to variations in the possible shearing mechanisms that connect the WZ and RS structures, which can create crystallites with up to ∼ ±10° of relative rotations about the RS [001] direction, as well as further rotation of these crystallites as they grow, merge, and coalesce.

Crystal structure↗
NASA NTRS1994

Validation of the RPLUS3D Code for Supersonic Inlet Applications Involving Three-Dimensional Shock Wave-Boundary Layer Interactions

A three-dimensional computational fluid dynamics code, RPLUS3D, which was developed for the reactive propulsive flows of ramjets and scramjets, was validated for glancing shock wave-boundary layer interactions. Both laminar and turbulent flows were studied. A supersonic flow over a wedge mounted on a flat plate was numerically simulated. For the laminar case, the static pressure distribution, velocity vectors, and particle traces on the flat plate were obtained. For turbulent flow, both the Baldwin-Lomax and Chien two-equation turbulent models were used. The static pressure distributions, pitot pressure, and yaw angle profiles were computed. In addition, the velocity vectors and particle traces on the flat plate were also obtained from the computed solution. Overall, the computed results for both laminar and turbulent cases compared very well with the experimentally obtained data.

Kapoor, Kamlesh↗
DOE OSTI2021

A fourth-order Lagrangian discontinuous Galerkin method using a hierarchical orthogonal basis on curvilinear grids

The existing high-order Lagrangian discontinuous Galerkin (DG) hydrodynamic methods are restricted to using quadratic meshes with quadratic polynomials (P2), which in turn, yield up to third-order accuracy. Additionally, these existing DG hydrodynamic schemes, when extended to work with cubic meshes and cubic polynomials (P3), can be unstable on strong-shock problems. Therefore, this paper presents a new Lagrangian discontinuous Galerkin (DG) hydrodynamic method to simulate compressible material dynamics ( e.g., gasses, fluids, and solids) with strong-shocks using cubic meshes and cubic polynomials, and delivers up to fourth-order accuracy on smooth flows. The stability on shock problems is achieved using new hierarchical orthogonal basis functions and a new subcell mesh stabilization (SMS) scheme for cubic meshes. The accuracy and robustness of the new high-order accurate Lagrangian DG hydrodynamic method is demonstrated by simulating a diverse suite of challenging test problems covering gas and solid dynamic problems on curvilinear meshes.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1986

The clumpy circumstellar medium around young supernova remnants

Each of the youngest supernova remnants known in the Milky Way, Cas A, Kepler's SNR, and Tycho's SNR, shows a different morphological structure caused by different conditions in the progenitor stars and their surroundings. In all three cases, however, the observed shells have a thickness of about 1/4 the radius, a sharp outer edge, and significant brightness irregularities. These features require that the circumstellar medium be highly clumped. To investigate the phenomenon, models of the expansion have been constructed using a one-dimensional spherical hydrodynamic code. As a supernova shock moves down the external density gradient of the star, material behind the shock begins to go into free expansion. Then as surrounding material is encountered a reverse shock moving back into the ejectum will be formed. Until the expansion has swept up about eight times the ejected mass when the situation can be considered as a point explosion in its surroundings, the dynamics are controlled by conditions between the shocks. The region is also where the synchrotron radio emission from relativistic electrons trapped in magnetic fields arises. Initial particles and fields are accelerated and amplified by eddy motion at the interface between the ejected and swept-up material and at the boundaries of clumps. Polarimetry shows that these SNR have a net radial orientation of their magnetic fields apparently from stretching by Rayleigh-Taylor instabilities at the contact surfaces. Without clumps the observed shell is much too narrow and steep on the inside.

Dickel, J. R.↗
NASA NTRS1994

Cosmic bullets as particle accelerators and radio sources

We have simulated in two dimensions the dynamical evolution of dense gas clouds(`cosmic bullets') moving supersonically through a uniform low-density medium. The diffusive shock acceleration of relativistic protons (cosmic rays) and their dynamical feedback on the background flow are included by the two-fluid model for this process. The acceleration of relativistic electrons is approximated by a test-particle model, and a passive magnetic field is followed by a simple advection scheme. Strong bow shocks, with Mach numbers similar to that of a bullet's motion, are the most important particle accelerators in the flow, while tail shocks and shocks inside the bullets do not play generally significant roles in this regard. For our simulation parameters, approximately greater than 10% of the initial bullet kinetic energy is converted to a combination of internal energy of gas and cosmic-ray protons by the time the bullets begin to be disrupted. Characteristically, the cosmic rays gain several percent of the available kinetic energy. Bullet destruction on timescales only a little larger than the ram pressure bullet crushing time begins in response to Kelvin-Helmholtz and especially to Rayleigh-Taylor instabilities along the forward bullet surface. For dense bullets this happens before the bullet is stopped by ram pressure. According to our simple model for synchrotron emission from relativistic electrons accelerated and transported within the flows, that emission increases rapidly as the bullet begins to fragment, when it is strongly dominated by field enhancement in sheared flows. Synchrotron emission from the acceleration region within the bow shock is, by contrast, much weaker.

Jones, T. W.↗
DOE OSTI2020

Nonideal mixing effects in warm dense matter studied with first-principles computer simulations

Here, we study nonideal mixing effects in the regime of warm dense matter (WDM) by computing the shock Hugoniot curves of BN, MgO, and MgSiO 3 . First, we derive these curves from the equations of state (EOS) of the fully interacting systems, which were obtained using a combination of path integral Monte Carlo calculations at high temperature and density functional molecular dynamics simulations at lower temperatures. We then use the ideal mixing approximation at constant pressure and temperature to rederive these Hugoniot curves from the EOS tables of the individual elements. We find that the linear mixing approximation works remarkably well at temperatures above ~2 × 10 5 K, where the shock compression ratio exceeds ~3.2. The shape of the Hugoniot curve of each compound is well reproduced. Regions of increased shock compression, which emerge because of the ionization of L and K shell electrons, are well represented, and the maximum compression ratio of the Hugoniot curves is reproduced with high precision. Some deviations are seen near the onset of the L shell ionization regime, where ionization equilibrium in the fully interacting system cannot be well reproduced by the ideal mixing approximation. This approximation also breaks down at lower temperatures, where chemical bonds play an increasingly important role. However, the results imply that the equilibrium properties of binary and ternary mixtures in the regime of WDM can be derived from the EOS tables of the individual elements. This significantly simplifies the characterization of binary and ternary mixtures in the WDM and plasma phases, which otherwise requires large numbers of more computationally expensive first-principles computer simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Unraveling the implications of finite specimen size on the interpretation of dynamic experiments for polycrystalline aluminum through direct numerical simulations

Normal and Pressure-shear plate impact (NPI and PSPI) experiments are popular experimental techniques for studying the mean-field macroscopic behavior of polycrystalline metals under high-rate dynamic loading. However, since both configurations rely upon geometry for subjecting the specimen to high strain rates, these experiments often involve a limited specimen size. Moreover, because of the inherent heterogeneities present within polycrystalline metals, it is difficult to ascertain if the size of the specimen and/or regions where measurements are made are sufficiently large for making representative inferences about the mean-field macroscopic properties from single-point velocity measurements. In the present study, we quantify the expected measurement variability on observable point measurements in NPI and PSPI experiments by carrying out direct numerical simulations (DNS) of statistically representative polycrystalline microstructures subjected to dynamic compression and compression-shear loading. In particular, we consider the role of specific material heterogeneities (e.g. the grain-to-grain difference in size, crystallographic orientation) on dispersion in the normal and transverse particle velocity records and on local fluctuations in key state variables (e.g. velocity, accumulated plastic strain) by incorporating these effects directly into a representative synthetic microstructure geometry and crystalline description of pure polycrystalline aluminum. The form of the present study is a large parametric investigation, consisting of ten ensembles of one hundred simulations. Each of the thousand simulations reflects a randomly realized synthetic microstructure in one of five cases of decreasing average grain size for the two loading configurations. Our analysis of the DNS results demonstrates that for both of these experimental configurations, the grain size directly correlates with the coefficient of variation (CV) in simulated point measurements, showing a convergent decrease in CV to zero (i.e. particle velocity record approaches the mean-field value) with decreasing grain size. Remarkably, the magnitude of variations in the particle velocity record is shown to be largest where the deviatoric stresses are most significant. In the case of NPI, this occurs at the elastic and plastic wavefront, whereas, in the case of PSPI, the magnitude of fluctuations are approximately constant throughout the experimental window time. The reasoning for the scatter in particle velocity due to the heterogeneous microstructure is demonstrated to be dependent on the mechanisms for accommodating deformation and on the interaction of reflection waves generated at sites of heterogeneities occurring at the scale of grains. Lastly, we develop a power-law description for the magnitude of scattering versus characteristic length, which provides a statistical framework for assessing the required number of grains per characteristic specimen dimension for minimizing scatter within these two experimental configurations (NPI, PSPI).

36 MATERIALS SCIENCE↗
NASA NTRS2009

CFD Analysis of Nozzle Jet Plume Effects on Sonic Boom Signature

An axisymmetric full Navier-Stokes computational fluid dynamics study is conducted to examine nozzle exhaust jet plume effects on the sonic boom signature of a supersonic aircraft. A simplified axisymmetric nozzle geometry, representative of the nozzle on the NASA Dryden NF-15B Lift and Nozzle Change Effects on Tail Shock research airplane, is considered. The computational fluid dynamics code is validated using available wind-tunnel sonic boom experimental data. The effects of grid size, spatial order of accuracy, grid type, and flow viscosity on the accuracy of the predicted sonic boom pressure signature are quantified. Grid lines parallel to the Mach wave direction are found to give the best results. Second-order accurate upwind methods are required as a minimum for accurate sonic boom simulations. The highly underexpanded nozzle flow is found to provide significantly more reduction in the tail shock strength in the sonic boom N-wave pressure signature than perfectly expanded and overexpanded nozzle flows. A tail shock train in the sonic boom signature is observed for the highly underexpanded nozzle flow. Axisymmetric computational fluid dynamics simulations show the flow physics inside the F-15 nozzle to be nonisentropic and complex. Although the one-dimensional isentropic nozzle plume results look reasonable, they fail to capture the sonic boom shock train in the highly underexpanded nozzle flow.

Bui, Trong T.↗
NASA NTRS1983

Interactive calculation procedures for mixed compression inlets

The proper design of engine nacelle installations for supersonic aircraft depends on a sophisticated understanding of the interactions between the boundary layers and the bounding external flows. The successful operation of mixed external-internal compression inlets depends significantly on the ability to closely control the operation of the internal compression portion of the inlet. This portion of the inlet is one where compression is achieved by multiple reflection of oblique shock waves and weak compression waves in a converging internal flow passage. However weak these shocks and waves may seem gas-dynamically, they are of sufficient strength to separate a laminar boundary layer and generally even strong enough for separation or incipient separation of the turbulent boundary layers. An understanding was developed of the viscous-inviscid interactions and of the shock wave boundary layer interactions and reflections.

Reshotko, Eli↗
NASA NTRS1993

Theoretical studies in interstellar cloud chemistry

This final report represents the completion of the three tasks under the purchase order no. SCPDE5620,1,2F. Chemical composition of gravitationally contracting, but otherwise quiescent, interstellar clouds and of interstellar clouds traversed by high velocity shocks, were modeled in a comprehensive manner that represents a significant progress in modeling these objects. The evolutionary chemical modeling, done under this NASA contract, represents a notable advance over the 'classical' fixed condition equilibrium models because the evolutionary models consider not only the chemical processes but also the dynamical processes by which the dark interstellar clouds may have assumed their present state. The shock calculations, being reported here, are important because they extend the limited chemical composition derivable from dynamical calculations for the total density and temperature structures behind the shock front. In order to be tractable, the dynamical calculations must severely simplify the chemistry. The present shock calculations take the shock profiles from the dynamical calculations and derive chemical composition in a comprehensive manner. The results of the present modeling study are still to be analyzed with reference to astronomical observational data and other contemporary model predictions. As far as humanly possible, this analysis will be continued with CRE's (Creative Research Enterprises's) IR&D resources, until a sponsor is found.

Chiu, Y. T.↗
NASA NTRS1991

Generation of Pc 1 waves by the ion temperature anisotropy associated with fast shocks caused by sudden impulses

The high correlation of Pc 1 events with magnetospheric compressions is known. A mechanism is proposed which leads to the generation of Pc 1 waves. The interaction of a dynamic pressure pulse with the earth's bow shock leads to the formation of a weak fast-mode shock propagating into the magnetoshealth. The shock wave can pass right through a tangential discontinuity (magnetopause) and into the magnetosphere, without disturbing either of the structures. In a quasiperpendicular geometry, the shock wave exhibits anisotropic heating. This anisotropy drives unstable ion-cyclotron waves which can contribute to the generation of the Pc 1 waves which are detected. The viability of the mechanism is demonstrated with simulations. This mechanism could explain the peak in the occurrence of observed Pc 1 waves in the postnoon sector where a field-aligned discontinuity in the solar wind would most often be parallel to the magnetopause surface due to the average Parker-spiral magnetic-field configuration.

Mandt, M. E.↗
DOE OSTI2021

Experimental Observations of Laser-Driven Tin Ejecta Microjet Interactions

The study of high-velocity particle-laden flow interactions is of importance for the understanding of a wide range of natural phenomena, ranging from planetary formation to cloud interactions. Experimental observations of particle dynamics are sparse given the difficulty of generating high-velocity flows of many particles. Ejecta microjets are micron-scale jets formed by strong shocks interacting with imprinted surfaces to generate particle plumes traveling at several kilometers per second. As such, the interaction of two ejecta microjets provides a novel experimental methodology to study interacting particle streams. In this Letter, we report the first time sequences of x-ray radiography images of two interacting tin ejecta microjets taken on a platform designed for the OMEGA Extended Performance (OMEGA EP) laser. We observe that the microjets pass through each other unattenuated for the case of 11.7±3.2 GPa shock pressures and jet velocities of 2.2±0.5 km/s but show strong interaction dynamics for 116.0±6.1 GPa shock pressures and jet velocities of 6.5±0.5 km/s. Furthermore, we find that radiation-hydrodynamic simulations of the experiments are able to capture many aspects of the collisional behavior, such as the attenuation of jet velocity in the direction of propagation, but are unable to match the full spread of the strongly interacting cloud.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The Potential Energy Hotspot: Effects of Impact Velocity, Defect Geometry, and Crystallographic Orientation

In energetic materials, the localization of energy into “hotspots” is known to dictate the initiation of chemical reactions and detonation. Recent all-atom simulations have shown that more energy is localized as internal potential energy (PE) than can be inferred from the kinetic energy (KE) alone. The mechanisms associated with pore collapse and hotspot formation are known to depend on pore geometry and dynamic material response such as plasticity. Therefore, we use molecular dynamics (MD) simulations to characterize shock-induced pore collapse and the subsequent formation of hotspots in 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a highly anisotropic molecular crystal, for various defect shapes, shock strengths, and crystallographic orientations. We find that the localization of energy as PE is consistently larger than the KE in cases with significant plastic deformation. Additionally, an analysis of MD trajectories reveals the underlying molecular- and crystal-level processes that govern the effect of orientation and pore shape on PE localization. We find that the regions of highest PE relate to the areas of maximum plastic deformation, while KE is maximized at the point of impact. Comparisons against octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) reveal less energy localization in TATB, which could be a contributing factor to the latter’s insensitivity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Imaging velocity interferometer system for any reflector (VISAR) diagnostics for high energy density sciences

Two variants of optical imaging velocimetry, specifically the one-dimensional streaked line-imaging and the two-dimensional time-resolved area-imaging versions of the Velocity Interferometer System for Any Reflector (VISAR), have become important diagnostics in high energy density sciences, including inertial confinement fusion and dynamic compression of condensed matter. Here, we give a brief review of the historical development of these techniques, then describe the current implementations at major high energy density (HED) facilities worldwide, including the OMEGA Laser Facility and the National Ignition Facility. We illustrate the versatility and power of these techniques by reviewing diverse applications of imaging VISARs for gas-gun and laser-driven dynamic compression experiments for materials science, shock physics, condensed matter physics, chemical physics, plasma physics, planetary science and astronomy, as well as a broad range of HED experiments and laser-driven inertial confinement fusion research.

47 OTHER INSTRUMENTATION↗