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,387 records · Page 77

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Multi-frame x-ray radiography and image tracking for quantification of expansion in laser-driven tin ejecta microjets

One regime of experimental particle-laden flow study involves ejecta microjets—often defined as a stream of micrometer-scale particles generated through shock interaction with a non-uniform surface and generally travel above 1 km/s. In order to capture the change in characteristics as a function of propagation time, we apply a multi-frame x-ray radiography platform to observe and track the jet transport dynamics. A synchrotron x-ray source allows us to perform quantitative analyses and comparisons between the eight images captured by the imaging system. Observation of a single jet through time allows the use of a cross correlation algorithm to independently track various regions within the jet and quantify the jet expansion over time using normalized area and normalized areal density values. Through a comparison with the calculated values of ballistic transport, these findings show less expansion than expected for ballistically transporting particles. In conclusion, this work combines multi-frame synchrotron radiography with image tracking to establish a foundation for future studies on jet transport and particle interaction dynamics.

Sun, Yuchen [Lawrence Livermore National Laborator↗

NASA NTRS1996

Dynamics of the Impacts

The SL9 impacts are best known by their plumes, several of which were imaged towering 3000 km above Jupiter's limb. The heat released when they fen back on the atmosphere produced the famous infrared main events. The reentry shock was observed directly in CO emission and indirectly through thermally glowing dust. The relatively low color temperature of the latter implies that most of the energy released on reentry was radiated. If so, the R impact released some 3 x 10(exp 26) ergs, and the larger L, K, and G impacts as much as 3 x 10(exp 27) ergs. These estimates agree well with impact energies deduced from tidal disruption computations. The duration of the infrared events measures time aloft and hence provides a second measure of plume height. These indicate that plume height was only weakly dependent on impact energy. Evidently all the plumes were launched at roughly 10-13 km/s. Using a semi-analytic model for the deceleration, disintegration, and destruction of intruding bodies by an ever-vigilant atmosphere, we find that similar plume heights is a direct consequence of smaller impactors exploding at higher altitudes, in such a way that the different explosions were geometrically similar. We then compare the predictions of our model to the Venerian cratering record, which provides a good statistical test of impact deceleration by a thick atmosphere. Chemistry should have provided an independent measure of explosion altitude: abundant shock-generated CS, CS2 and HCN indicates a source in dry jovian air, above the putative water clouds. However, the Galileo Probe results seem to imply that we should expect no more. Observed water and S2 are consistent with a somewhat oxidized gas (presumably the comet itself), but the absence of SO2 and CO2 shows that conditions were neither too oxidizing nor the shocks too hot.

Zahnle, Kevin↗

NASA NTRS2020

Wind Tunnel Investigation of the Supersonic Stage Separation Aerodynamics of a Generic 0.0175-Scale Bimese Two-Stage-to-Orbit Reusable Launch Vehicle Configuration

A wind tunnel investigation was conducted of the supersonic stage separation aerodynamics of a generic two-stage-to-orbit bimese wingbody configuration in the NASA Langley Research Center Unitary Plan Wind Tunnel. Proximity and isolated model testing was conducted at Mach numbers of 2.3, 3.0, and 4.5 and a unit Reynolds number of 2.0 million per foot using 0.0175-scale models of the Langley Glide-Back Booster concept designated as the orbiter and booster in belly-to-belly and back-to-belly configurations. Longitudinal forces and moments were obtained on both models and surface static pressure measurements were obtained on the orbiter model at 328 relative proximity locations and at relative angles of attack of 0 degrees and 5 degrees. The test results supported a larger effort to develop and validate experimental and computational tools applicable to the design and simulation of stage separation and abort procedures for reusable launch vehicles composed of multiple bodies, including winged bodies. An initial proof-of-concept experiment featuring low-cost uninstrumented models was conducted to verify an emerging automated model control system and new support system hardware, and to identify potential model and support system blockage and unsteady aerodynamics/model dynamics prior to committing to higher-fidelity instrumented models. This investigation led to upgrades in the facility stage separation hardware, calibration and testing techniques and capabilities, and data analysis and documentation methodologies that have been extended to the more recent NASA Constellation and Space Launch System crew and cargo launch vehicle programs. A virtual diagnostics interface methodology was used to facilitate the design of the stage separation support hardware, to position the models in the test section, and to define the experimental test space. Advances in the facility automated model positioning system established a foundation for the development of a continuous-sweep data acquisition technique that is responsible for significant productivity improvements to the current NASA Space Launch System test program. The automated model positioning capability was leveraged to conduct a companion statistically-designed stage separation experiment requiring randomization of the relative proximity positions of the orbiter and booster models. The respective zones of influence and interference effects of the orbiter and booster were identified from three-dimensional scatter plots, contour and influence maps, and two-dimensional plotting methods. The highly-nonlinear, shock-dominated aerodynamic characteristics of the orbiter and booster in the Unitary Plan Wind Tunnel exhibited good agreement with independent test data obtained in a NASA Marshall Space Flight Center wind tunnel and with computational fluid dynamics predictions using a compressible, three-dimensional flow solver and an inviscid, unstructured Cartesian method.

Erickson, Gary E.↗

NASA NTRSDate not supplied

Comparison of OVERFLOW Computational and Experimental Results for a Blunt Mars Entry Vehicle Concept during Supersonic Retropropulsion

Simulations of unsteady supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle were performed using the OVERFLOW Computational Fluid Dynamics (CFD) solver. High-fidelity flow solver techniques, including Detached Eddy Simulation (DES) turbulence modeling and Adaptive Mesh Refinement (AMR), were employed to obtain improved realism in CFD predictions. Simulation conditions and geometry configurations were designed to match specific runs in the Descent System Study (DSS) wind tunnel testing (WTT) campaign. The accuracy of each simulation is assessed by direct comparison to experimental data. Comparisons of computational predictions of the SRP flowfield and bow shock shape to experimental schlieren imaging show reasonable prediction of mean shock shape, with approximately 10% similarity in shock standoff distance for selected conditions, as well as similarity in local, time-varying fluctuations of the shock-plume interaction. Comparisons of discrete measurements of surface pressure coefficient (Cp) indicate CFD accuracy within approximately 10% of the experiment across the majority of the model heatshield, with larger variations at some of the heatshield edge locations with stronger flow unsteadiness. Simulated unsteadiness of these chaotic flows, which were highly dynamic and multi-modal, was shown to be within 20-40% of experimentally-measured pressure standard deviation (SD) for the majority of the sampled locations.

Supersonic Retropropulsion↗

NASA NTRSDate not supplied

Comparison of OVERFLOW Computational and Experimental Results for a Blunt Mars Entry Vehicle Concept during Supersonic Retropropulsion

Supersonic Retropropulsion↗

DOE OSTI2021

Novel Mechanisms for the Formation of Dislocation Cell Patterns in BCC Metal

In this study, we present the first simulation results of the formation of dislocation cell wall microstructures in tantalum subjected to shock loading. Dislocation patterns and cell wall formation are important to understanding the mechanical properties of the materials in which they spontaneously arise, and yet the processing and self-assembly mechanisms leading to their formation are poorly understood. By employing transmission electron microscopy and discrete dislocation dynamics, we propose a new mechanism involving coplanar dislocations and pseudo-dipole mixed dislocation arrays that is essential to the pattern formation process. Our large-scale 3D DDD simulations demonstrate the self-organization of dislocation networks into cell walls in deformed BCC metal (tantalum) persisting at strain ε = 20%. The simulation analysis captures several crucial aspects of how the dislocation cell pattern affects metal plasticity, as observed in experiments. Although experimental evidence is inconclusive regarding whether cell wall formation takes place at the shock front, after the shock, during release, or when the sample has had enough time to relax post-recovery, our simulations indicate cell wall formation occurs after the shock and before release. The extended Taylor hardening composite model effectively considers the non-uniform dislocation density when cell walls form and accurately describes the corresponding flow stress.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Numerical relativity in spherical coordinates: A new dynamical spacetime and general relativistic MHD evolution framework for the Einstein Toolkit

We present SphericalNR, a new framework for the publicly available Einstein Toolkit that numerically solves the Einstein field equations coupled to the equations of general relativistic magnetohydrodynamic (GRMHD) in a 3+1 split of spacetime in spherical coordinates without symmetry assumptions. The spacetime evolution is performed using reference-metric versions of either the Baumgarte-Shapiro-Shibata-Nakamura equations or the fully covariant and conformal Z4 system with constraint damping. We have developed a reference-metric version of the Valencia formulation of GRMHD with a vector potential method, guaranteeing the absence of magnetic monopoles during the evolution. In our framework, every dynamical field (both spacetime and matter) is evolved using its components in an orthonormal basis with respect to the spherical reference metric. Furthermore, all geometric information about the spherical coordinate system is encoded in source terms appearing in the evolution equations. This allows for the straightforward extension of Cartesian high-resolution shock-capturing finite volume codes to use spherical coordinates with our framework. To this end, we have adapted GRHydro, a Cartesian finite volume GRMHD code already available in the Einstein Toolkit, to use spherical coordinates. We present the full evolution equations of the framework, as well as details of its implementation in the Einstein Toolkit. Finally, we validate SphericalNR by demonstrating it passes a variety of challenging code tests in static and dynamical spacetimes.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2019

WRLES: Wave Resolving Large-Eddy Simulation Code, Theory and Usage

A computational fluid dynamics code has been developed for large-eddy simulations (LES) of turbulent flow. The code uses high-order of accuracy and high-resolution numerical methods to minimize solution error and maximize the resolution of the turbulent structures. Spatial discretization is performed using explicit central differencing. The central differencing schemes in the code include 2nd- to 12th-order standard central difference methods as well as 7-, 9-, 11- and 13-point dispersion relation preserving schemes. Solution filtering and high-order shock capturing are included for stability. Time discretization is performed using multistage Runge-Kutta methods that are up to 4th order accurate. Several options are available to model turbulence including: Baldwin-Lomax and Spalart-Allmaras Reynolds-averaged Navier-Stokes turbulence models, and Smagorinsky, Dynamic Smagorinsky and Vreman sub-grid scale models for LES. This report presents the theory behind the numerical and physical models used in the code and provides a user's manual to the operation of the code.

High-Order↗

NASA NTRS1990

New developments in adaptive methods for computational fluid dynamics

New developments in a posteriori error estimates, smart algorithms, and h- and h-p adaptive finite element methods are discussed in the context of two- and three-dimensional compressible and incompressible flow simulations. Applications to rotor-stator interaction, rotorcraft aerodynamics, shock and viscous boundary layer interaction and fluid-structure interaction problems are discussed.

Oden, J. T.↗

NASA NTRS1990

Solar radio astronomy at low frequencies

The characteristics of solar radio emissions at decametric to kilometric wavelengths are reviewed. Special attention is given to the radiation of the quiet sun at several metric and decametric wavelengths and to nonthermal radiation from the active sun, including radio bursts of type III (electron beams), type-III bursts from behind the sun, storms of type III bursts, the flare-associated radio bursts, type II bursts (shock waves), and shock-associated bursts. It is pointed out that almost no observations have been made so far of solar radiation between about 20 MHz and about 2 MHz. Below about 2 MHz, dynamic spectra of flux densities of solar burst have been recorded in space and observations were made of the directions of centroids and characteristic sizes of the emitting sources.

Dulk, George A.↗

DOE OSTI2025

Pressure-temperature equation of state of Al 2 ⁢O 3 up to 14 Mbar and 40 kK

Sapphire (Al 2 ⁢O 3 ), known for its remarkable incompressibility at ambient conditions, plays a pivotal role in both static and dynamic compression research. Accurately characterizing its equation of state (EoS) is essential for these applications. Here, we present a complete Hugoniot of Al 2 ⁢ O 3 as locus of experimentally assessed, high-precision, pressure, density and temperature states up to 14 Mbar and 43 kK. The Hugoniot is established with single shock experiments using magnetically launched hyper velocity flyers on the Z Accelerator at Sandia National Laboratories. We explore principal Hugoniot states at very high shock 𝑇 and 𝑝 in the solid phase, tracking the solid-liquid boundary and culminating at 2.4-fold compression, where data provides a direct constraint on the liquid phase. Corresponding shock release data probe thermodynamic states complementary to the Hugoniot and place additional constraints on tabular EoS models. Our findings indicate a significant deviation from existing tabular EoS models for Al 2 ⁢ O 3 dictating a comprehensive overhaul. We develop two advanced EoSs for Al 2 ⁢ O 3 the SESAME 97412 model, featuring an extensive phase diagram that includes three solid phases and the liquid phase, and the updated LEOS 2200m2 model. EoS development is assisted with Quantum Molecular Dynamics simulations. Our experimental data allows for stringent testing of our EoSs. Both models accurately capture the Hugoniot of Al 2 ⁢O 3 up to the highest pressures and temperatures. Rigorous experimental determination of extreme pressures and temperatures, paired with sophisticated models, advances the frontier of EoS development beyond 1 terapascal.

Kalita, Patricia [Sandia National Laboratories (SN↗

DOE OSTI2024

Large-Eddy Simulation Study of Flow and Combustion Dynamics in a Full-Scale Hydrogen–Air Rotating Detonation Combustor-Stator Integrated System

In the present work, a first-of-its-kind three-dimensional (3D) large-eddy simulation (LES) study is conducted to numerically investigate the combustion dynamics as well as aero-thermal phenomena in a full-scale nonpremixed hydrogen–air rotating detonation engine (RDE) (with a diverging-shaped lower-end wall), when integrated with nozzle guide vanes (NGV) acting as the turbine stator. The wall-modeled LES framework incorporates hydrogen–air detailed chemical kinetics and adaptive mesh refinement (AMR). A comparative analysis is carried out for two operating conditions with different fuel/air mass flow rates but global equivalence ratio of unity, and considering RDE configurations without and with stator. The LES model is validated against available experimental data for the low mass flux condition with respect to detonation wave speed/height, wave dynamics, and axial static pressure distribution. Numerical results indicate significant deflagrative combustion occurring in the fill region near the inner wall due to formation of recirculation zones in the injection near-field driven by the backward facing step. The leading detonation wave is found to be trailed by an azimuthal reflected-shock combustion (ARSC) wave, consistent with experimental observations, which consumes unburned vitiated reactants that leak through the main detonation wave. The main detonation wave characteristics, such as detonation wave speed/height and combustion efficiency, do not change appreciably with the presence of NGV. A novel combustion diagnostic technique based on chemical explosive mode analysis (CEMA) is employed to quantify the fraction of heat release occurring in the detonative mode versus deflagrative mode for the simulated conditions. The exit flow is found to be nearly fully subsonic and supersonic for the low and high mass flux conditions, respectively. Further analysis of the exit flow profiles shows that the presence of NGV renders the flow more axial and significantly impacts the exit Mach number and total pressure, while the total temperature shows negligible change. In addition, the low mass flux operating point, despite exhibiting more deflagrative losses within the combustor, yields overall lower pressure drop from plenum to exhaust, which is mainly attributed to lower pressure drop across the injectors. Lastly, the rotating detonation engine-nozzle guide vanes (RDE-NGV) configuration exhibits higher total pressure loss compared to rotating detonation engine (RDE) without stator across both the mass flux conditions. In conclusion, this study extends the state-of-the-art in numerical modeling of pressure gain combustion (PGC) systems by demonstrating high-fidelity 3D reacting LES of full-scale RDE-NGV systems relevant to RDE-turbine integration for stationary power generation.

Combustion↗

DOE OSTI2020

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1985

Plasma heating at collisionless shocks due to the kinetic cross-field streaming instability

Heating at collisionless shocks due to the kinetic cross-field streaming instability, which is the finite beta (ratio of plasma to magnetic pressure) extension of the modified two stream instability, is studied. Heating rates are derived from quasi-linear theory and compared with results from particle simulations to show that electron heating relative to ion heating and heating parallel to the magnetic field relative to perpendicular heating for both the electrons and ions increase with beta. The simulations suggest that electron dynamics determine the saturation level of the instability, which is manifested by the formation of a flattop electron distribution parallel to the magnetic field. As a result, both the saturation levels of the fluctuations and the heating rates decrease sharply with beta. Applications of these results to plasma heating in simulations of shocks and the earth's bow shock are described.

Winske, D.↗

DOE OSTI2024

Atomic cluster expansion potential for large scale simulations of hydrocarbons under shock compression

We present an Atomic Cluster Expansion (ACE) machine learned potential developed for high-fidelity atomistic simulations of hydrocarbons, targeting pressures and temperatures near and above supercritical fluid regimes for molecular fluids. A diverse set of stoichiometries were covered in training, including 1:0 (pure carbon), 1:4 (methane), and 1:1 (benzene), and rich bonding environments sampled at supercritical temperatures, hydrogen rich, reactive mixtures where metastable stoichiometries arise, including 1:2 (ethylene) and 1:3 (ethane). A high-fidelity training database was constructed by performing large-scale quantum molecular dynamic simulations [density functional theory (DFT) MD] of diamond, graphite, methane, and benzene. A novel approach to selecting structures from DFT MD is also presented, which allows for the rapid selection of unique DFT MD frames from complex trajectories. Comparisons to DFT and experimental data demonstrate that the presented ACE potential accurately reproduces isotherms, carbon melting curves, radial distribution functions, and shock Hugoniots for carbon and hydrocarbon systems for pressures up to 100 GPa and temperatures up to 6000 K for hydrocarbon systems and up to 9000 K for pure carbon systems. This work delivers a potential that can be used for accurate, large-scale simulations of shocked hydrocarbons and demonstrates a methodology for fitting and validating machine learning interatomic potentials to complex molecular environments, which can be applied to energetic materials in future works.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

New dynamic diamond anvil cell for time-resolved radial x-ray diffraction

The dynamic diamond anvil cell (dDAC) is a recently developed experimental platform that has shown promise for studying the behavior of materials at strain rates ranging from intermediate to quasi-static and shock compression regimes. Combining dDAC with time-resolved x-ray diffraction (XRD) in the radial geometry (i.e., with incident x-rays perpendicular to the axis of compression) enables the study of material properties such as strength, texture evolution, and deformation mechanisms. This work describes a radial XRD dDAC setup at beamline P02.2 (Extreme Conditions Beamline) at DESY’s PETRA III synchrotron. Time-resolved radial XRD data are collected for titanium, zirconium, and zircon samples, demonstrating the ability to study the strength and texture of materials at compression rates above 300 GPa/s. In addition, the simultaneous optical imaging of the DAC sample chamber is demonstrated. The ability to conduct simultaneous radial XRD and optical imaging provides the opportunity to characterize plastic strain and deviatoric strain rates in the DAC at intermediate rates, exploring the strength and deformation mechanisms of materials in this regime.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

New dynamic diamond anvil cell for time-resolved radial x-ray diffraction

47 OTHER INSTRUMENTATION↗

NASA NTRS1973

Theoretical and experimental studies of reentry plasmas

A viscous shock-layer analysis was developed and used to calculate nonequilibrium-flow species distributions in the plasma layer of the RAM vehicle. The theoretical electron-density results obtained are in good agreement with those measured in flight. A circular-aperture flush-mounted antenna was used to obtain a comparison between theoretical and experimental antenna admittance in the presence of ionized boundary layers of low collision frequency. The electron-temperature and electron-density distributions in the boundary layer were independently measured. The antenna admittance was measured using a four-probe microwave reflectometer and these measured values were found to be in good agreement with those predicted. Measurements were also performed with another type of circular-aperture antenna and good agreement was obtained between the calculations and the experimental results. A theoretical analysis has been completed which permits calculation of the nonequilibrium, viscous shock-layer flow field for a sphere-cone body. Results are presented for two different bodies at several different altitudes illustrating the influences of bluntness and chemical nonequilibrium on several gas dynamic parameters of interest. Plane-wave transmission coefficients were calculated for an approximate space-shuttle body using a typical trajectory.

Dunn, M. G.↗