Conservative Eulerian numerical methods for geometrical shock dynamics (GSD) and detonation shock dynamics (DSD)
Explore the source record for details and available documents.
97 MATHEMATICS AND COMPUTING↗
SEARCH · Search NASA
Results for “geometrical 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.
Generating Three-Axis Shock Response on Resonant Plates without Additional Dynamic Fixturing through Geometric Optimization
Explore the source record for details and available documents.
Alvis, Tyler Harrison [Sandia National Laboratorie↗
TPSAS-NF1676L-31264-DND
Nonlinear unsteady aerodynamic models (TSD, Euler/Navier-Stokes (RANS), LES) are complex, expensive solutions. CFD used by industry mostly for steady/static loads cases. Limit Cycle Oscillation (LCO) is an equilibrium point of the nonlinear aeroelastic system induced by nonlinear unsteady aerodynamics (shock dynamics, separated flow) and/or nonlinear structural dynamics (geometric and material nonlinearities). While Euler solutions can generate LCOs, viscosity yields improved accuracy (shock strength, position). Most common aeroelastic CFD codes currently consist of nonlinear unsteady aerodynamics with linear structural dynamics. Van der Pol system is a simple example of LCOs.
Walter A Silva↗
Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability
Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Coronal Magnetic Field Measurement from EUV Images Made by the Solar Dynamics Observatory
By measuring the geometrical properties of the coronal mass ejection (CME) flux rope and the leading shock observed on 2010 June 13 by the Solar Dynamics Observatory (SDO) mission's Atmospheric Imaging Assembly we determine the Alfven speed and the magnetic field strength in the inner corona at a heliocentric distance of approx. 1.4 Rs The basic measurements are the shock standoff distance (Delta R) ahead of the CME flux rope, the radius of curvature of the flux rope (R(sub c)), and the shock speed. We first derive the Alfvenic Mach number (M) using the relationship, Delta R/R(sub c) = 0.81[(gamma−1) M(exp 2) + 2] / [(gamma +1)(M2 − 1)], where gamma is the only parameter that needed to be assumed. For gamma = 4/3, the Mach number declined from 3.7 to 1.5 indicating shock weakening within the field of view of the imager. The shock formation coincided with the appearance of a type II radio burst at a frequency of approx. 300 MHz (harmonic component), providing an independent confirmation of the shock. The shock compression ratio derived from the radio dynamic spectrum was found to be consistent with that derived from the theory of fast-mode MHD shocks. From the measured shock speed and the derived Mach number, we found the Alfven speed to increase from approx 140 km/s to 460 km/s over the distance range 1.2-1.5 Rs. By deriving the upstream plasma density from the emission frequency of the associated type II radio burst, we determined the coronal magnetic field to be in the range 1.3-1.5 G. The derived magnetic field values are consistent with other estimates in a similar distance range. This work demonstrates that the EUV imagers, in the presence of radio dynamic spectra, can be used as coronal magnetometers
magnetic fields↗
Understanding Microstructural Effects on Dynamic Performance Towards the Development of Shock Metamaterials
With the recent advances in additive manufacturing (AM), long-range periodic lattice assemblies are being developed for vibration and shock mitigation components in aerospace and military applications with unique geometric and topological structures. There has been extensive work in understanding the static properties associated with varying topology of these lattice architectures, but there is almost no understanding of microstructural affects in such structures under high-strain rate dynamic loading conditions. Here we report the shock behavior of lattices with varying intrinsic grain structures achieved by post process annealing. High resolution 316L stainless steel lattices were 3D printed by a laser-powder bed fusion machine and characterized by computed tomography. Subsequent annealing resulted in stress-relieved and recrystallized lattices. Overall the lattices had strong cubic texture aligning with the x-, y- and z-directions of the build with a preference outside the build direction (z). The recrystallized sample had more equiaxed polygonal grains and a layer of BCC ferrite at the surface of the structure approximately 1 grain thick. Upon dynamic compression the as-deposited lattice showed steady compaction behavior while the heat-treated lattices exhibit negative velocity behavior indicative of failure. We attribute this to the stiffer BCC ferrite in the annealed lattices becoming damaged and fragmenting during compression.
36 MATERIALS SCIENCE↗
Photogrammetric Measurements of CEV Airbag Landing Attenuation Systems
High-speed photogrammetric measurements are being used to assess the impact dynamics of the Orion Crew Exploration Vehicle (CEV) for ground landing contingency upon return to earth. Test articles representative of the Orion capsule are dropped at the NASA Langley Landing and Impact Research (LandIR) Facility onto a sand/clay mixture representative of a dry lakebed from elevations as high as 62 feet (18.9 meters). Two different types of test articles have been evaluated: (1) half-scale metal shell models utilized to establish baseline impact dynamics and soil characterization, and (2) geometric full-scale drop models with shock-absorbing airbags which are being evaluated for their ability to cushion the impact of the Orion CEV with the earth s surface. This paper describes the application of the photogrammetric measurement technique and provides drop model trajectory and impact data that indicate the performance of the photogrammetric measurement system.
Barrows, Danny A.↗
High-speed x-ray phase contrast imaging and digital image correlation analysis of microscale shock response of an additively manufactured energetic material simulant
The performance of energetic materials subjected to dynamic loading significantly depends on their micro- and meso-scale structural morphology. The geometric versatility offered by additive manufacturing opens new pathways to tailor the performance of these materials. Additively manufactured energetic materials (AMEMs) have a wide range of structural characteristics with a hierarchy of length scales and process-inherent heterogeneities, which are hitherto difficult to precisely control. It is important to understand how these features affect AMEMs’ response under dynamic/shock loading. Therefore, temporally and spatially resolved measurements of both macroscopic behavior and micro- and meso-level processes influencing macroscopic behavior are required. In this paper, we analyze the shock compression response of an AMEM simulant loaded under several impact conditions and orientations. Furthermore, x-ray phase contrast imaging (PCI) is used to track features across the observed shock front and determine the linear shock velocity vs particle velocity equation of state, as well as to quantify the interior deformation fields via digital image correlation (DIC) analyses. Photon Doppler velocimetry is simultaneously used to measure the particle velocities of the specimens, which are consistent with those obtained from x-ray PCI. The DIC analyses provide an assessment of the average strain fields inside the material, showing that the average axial strain depends on the loading intensity and reaches as high as 0.23 for impact velocities up to 1.5 km/s. The overall results demonstrate the utility of x-ray PCI for probing “in-material” equation of state and interior strains associated with dynamic shock compression behavior of the AMEM simulant.
3D printing↗
Research study of space plasma boundary processes
Representation of the Earth's bow shock and magnetopause and their geometrically determined macrostructure was investigated. Computer graphic depictions of the global distributions of bow shock structures and elementary animation of the dynamics of those distributions in the changing solar wind were developed. The shock-foreshock boundary and subcritical bow shocks as observed by ISEE 1 and 2 are discussed.
Greenstadt, E. W.↗
Influence of Shock Wave on the Flutter Behavior of Fan Blades Investigated
Modern fan designs have blades with forward sweep; a lean, thin cross section; and a wide chord to improve performance and reduce noise. These geometric features coupled with the presence of a shock wave can lead to flutter instability. Flutter is a self-excited dynamic instability arising because of fluid-structure interaction, which causes the energy from the surrounding fluid to be extracted by the vibrating structure. An in-flight occurrence of flutter could be catastrophic and is a significant design issue for rotor blades in gas turbines. Understanding the flutter behavior and the influence of flow features on flutter will lead to a better and safer design. An aeroelastic analysis code, TURBO, has been developed and validated for flutter calculations at the NASA Glenn Research Center. The code has been used to understand the occurrence of flutter in a forward-swept fan design. The forward-swept fan, which consists of 22 inserted blades, encountered flutter during wind tunnel tests at part speed conditions.
Srivastava, Rakesh↗
User's Manual for Aerofcn: a FORTRAN Program to Compute Aerodynamic Parameters
The computer program AeroFcn is discussed. AeroFcn is a utility program that computes the following aerodynamic parameters: geopotential altitude, Mach number, true velocity, dynamic pressure, calibrated airspeed, equivalent airspeed, impact pressure, total pressure, total temperature, Reynolds number, speed of sound, static density, static pressure, static temperature, coefficient of dynamic viscosity, kinematic viscosity, geometric altitude, and specific energy for a standard- or a modified standard-day atmosphere using compressible flow and normal shock relations. Any two parameters that define a unique flight condition are selected, and their values are entered interactively. The remaining parameters are computed, and the solutions are stored in an output file. Multiple cases can be run, and the multiple case solutions can be stored in another output file for plotting. Parameter units, the output format, and primary constants in the atmospheric and aerodynamic equations can also be changed.
Conley, Joseph L.↗
Reflection and Refraction of Acoustic Waves by a Shock Wave
The presence of sound waves in one or the other of the fluid regions on either side of a shock wave is made apparent, in the region under superpressure, by acoustic waves (reflected or refracted according to whether the incident waves lie in the region of superpressure or of subpressure) and by thermal waves. The characteristics of these waves are calculated for a plane, progressive, and uniform incident wave. In the case of refraction, the refracted acoustic wave can, according to the incidence, be plane, progressive, and uniform or take the form of an 'accompanying wave' which remains attached to the front of the shock while sliding parallel to it. In all cases, geometrical constructions permit determination of the kinematic characteristics of the reflected or refractive acoustic waves. The dynamic relationships show that the amplitude of the reflected wave is always less than that of the incident wave. The amplitude of the refracted wave, whatever its type, may in certain cases be greater than that of the incident wave.
Brillouin, J.↗
Influences of shock imprinting on mix in a 3D-printed porous media
Mixing of materials in porous media can cause a significant impact on fusion yield as previously demonstrated by the National Ignition Facility (NIF) MARBLE Campaign. Initially, the reactants are separated, with deuterium in the lattice struts and a tritium gas fill in the voids. Lattice parameters such as the strut thickness and relative pitch, provide a control for the mix parameters in the experiment. Los Alamos National Laboratory’s (LANL) BOSQUE project looks to better understand how the mix of the reactants and shell materials impact the fusion burn and resultant yield on various laser platforms. xRAGE’s Eulerian hydrodynamics and adaptive mesh refinement (AMR) provide the unique ability to study the impacts of multiscale features of complex lattice structures. This modeling provides the ability to measure shock front variations as the wave progress’ through a given media. Initial conditions of the lattice are essential to accurately model mix and burn measured by experiment. By varying the initial orientation and densities of these lattice regions the early time dynamics of how the shock is launched into the system is changed and advocates for the study of resulting effects. In this work, we will study the sensitivities of shock effects in varying 3D printed geometric systems and how these shocks alter the structure and mix in the lattice. We will discuss both preliminary experimental results and simulations to help plan and constrain future experiments where we will study the impact of different lattice geometries and lattice bulk densities. This work concludes with the relative impacts of lattice geometries on shock speeds at different bulk densities and the resultant mix widths due to those shock interactions. We see agreement with theory at the higher end of our bulk density study, however, as we approach lower bulk density systems the dynamics of these interactions begin to change.
2 Photon Polymerization↗
An infrasound source analysis of the OSIRIS-REx sample return capsule hypersonic re-entry
The OSIRIS-REx sample return capsule's hypersonic re-entry into the atmosphere is a rare opportunity to test a variety of sonic boom source models since the projectile dimensions are well characterized. While the as-flown flight path is unknown, the predicted flight path enables a rough approximation of the source Mach number and location. Six infrasound microphones deployed in the boom carpet along the predicted flight path recorded impulsive signals from the OSIRIS-REx re-entry. Using a suite of atmosphere profiles and the geometric acoustics approximation, we estimate locations with uncertainty estimates along the flight path from which the signals were emitted. Acoustic overpressure and signal duration predictions from Whitham's far field theory, Carlson's simplified sonic boom prediction method, and a drag-dominated hypersonic model are analyzed with uncertainty estimates from the location estimate. While the Carlson simplified sonic boom prediction method could be accurate, our preference is for the drag-dominated source model. Using this source model with an inviscid Burgers's equation solver for propagation, we obtained an excellent match to the recorded data. In conclusion, these results will help better inform future sample return capsule re-entry observation campaigns as well as contribute to a better understanding of high altitude infrasonic sources.
58 GEOSCIENCES↗
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↗
Investigation on the explosion dynamics of large-format lithium-ion pouch cells
Explosion is the most extreme case of thermal runaway of lithium-ion (Li-ion) batteries. In this study, explosion dynamics of large-format Li-ion cells are investigated experimentally and numerically. Overcharge-to-explosion tests are conducted on 40 Ah Li-ion cells with Li[Ni 0.8 Co 0.1 Mn 0.1 ]O 2 cathode. Based on the explosion physics, shockwave and detonation models are used to characterize the shock effect of the cell explosion and evaluate the explosion equivalent. Von Neumann peak is observed on the pressure curves, and the shockwave velocity is supersonic at this time; the unwrinkled spherical flame phenomenon observed in the experiment indicates that it is detonation. Additionally, a geometric model is established based on the real testing scenario, and the explosion behavior is numerically studied. The characteristics of the explosion dynamics process are interpreted and the propagation mechanism of the shockwave are revealed; the deflagration to detonation transition (DDT) phenomenon in this process is caused by the formation of “hot spots”, and the explosion of the cells eventually turns into stable combustion. Here this study fills the gap in the research on thermal runaway of Li-ion cells, especially in the extreme cases of fire and explosion, and provide useful guidance for battery safety.
42 ENGINEERING↗
Carbon-carbon interactions in warm dense titanium carbide
Encouraged by experimental reports of diamond precipitation in carbon-containing materials under warm dense matter conditions and in shock compression experiments, we examine the behavior of carbon in TiC using density functional theory-molecular dynamics simulations. Two polymorphs of TiC are considered, the ambient-pressure B1 ($Fm$$\overline{3}$$m$) structure in which C-C interactions are prevented by geometric constraints, and a high-pressure Cmcm structure in which C atoms condense into zigzag chains. Chemistry-inspired bonding analyses confirm the covalently bound nature of these chains and further illuminate important interatomic interactions in both structures. Upon melting of B1 TiC, new short-range C-C interactions develop in the liquid, while the pre-existing C-C interactions in Cmcm TiC persist in the liquid. Here, the resulting carbon networks in the melt provide a promising environment for eventual diamond nucleation.
Ab initio molecular dynamics↗
Prospects for drag prediction using computational fluid dynamics
An account is given of CFD methods deemed likely to yield aerodynamic drag prediction improvements. Improvements in drag prediction are noted to largely proceed from improvements in viscous-layer resolution. Geometric modeling and grid generation are acknowledged to be major pacing items, and may be more important than turbulence modeling. The use of upwind differencing has led to improved solution accuracy due to the nature of the artificial viscosity in cases involving shocks and nonphysical dissipation across viscous layers. Multiblock algorithms will allow more complex geometries to be treated.
Roberts, Thomas W.↗