Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shock”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Laboratory-Based Thermal Shock Investigation of Heat Flux Sensors for the Mars 2020 Backshell

In 2012 during the entry, descent, and landing of the Mars Science Laboratory (MSL), the MSL Entry, Descent, and Landing Instrumentation (MEDLI) sensor suite was collecting in-flight heatshield pressure and temperature data. The data collected by the MEDLI instruments has since been used for reconstruction of vehicle aerodynamics, atmospheric conditions, aerothermal heating, and Thermal Protection System (TPS) performance as well as material response model validation and refinement. The Mars Entry, Descent, and Landing Instrumentation 2 (MEDLI2) sensor suite for the Mars 2020 heatshield and backshell is being designed to expand on the measurements and knowledge gained from MEDLI. Similar to MEDLI, MEDLI2 will measure the pressure and temperature of the heatshield. MEDLI2 will additionally measure the temperature, pressure, total heat flux, and radiative heat flux on the backshell.Since the backshell instrumentation is new to MEDLI2, Do No Harm (DNH) testing was conducted on instrumented backshell TPS (SLA-561V) panels. The panels consisted of four pressure port holes, one Mars Entry Atmospheric Data System (MEADS) pressure port plug, one MEDLI2 Integrated Sensor Plug (MISP) thermal plug, and one heat flux sensor. DNH testing was conducted to ensure the performance of the TPS was not degraded due to sensor integration and to characterize any TPS performance changes. The testing consisted of environmental testing— vibration, shock, thermal vacuum (TVAC) cycling— and bounding aerothermal (arc jet) testing. During arc jet testing, the heat flux sensors embedded in the SLA-561V panels exhibited an unexpected temporary reduction in the heat flux sensor temperature and response. After review of the test results, it was determined that this unexpected response was confined to the two heat flux sensors that experienced the greatest thermal shock condition. This condition consisted of a liquid nitrogen (LN2) bath that induced temperatures of approximately -190°C, and then a transition (thermal shock) to an arc jet test at a heat rate of approximately 21 W/cm2. Both heat flux sensors that were exposed to this thermal shock experienced a blister in the thermal coating during the arc jet test.Two heat flux sensor thermal shock test series were performed to investigate the cause of the blistering and subsequent energy release. In these tests, the heat flux sensor was first cold soaked in either a dry ice or LN2 bath to induce temperatures of approximately -78°C or -190°C, respectively. Then the sensors were thermally shocked using two propane torches with a heat rate of either approximately 8 W/cm2 or 21 W/cm2. The key findings indicated that there is a correlation between thermal shock and the blistering observed in the DNH test series, and that the cause appeared to be rooted in the heat flux sensor epoxy that encapsulates the sensor thermopile.Since the heat flux sensors are required to measure heat fluxes up to 15 W/cm2 during the Mars 2020 entry, a third test series was designed to determine if blistering is an issue at this maximum expected flight heat flux. Results from all three thermal shock test series and a discussion about whether or not blistering of the heat flux sensor thermal coating could be an issue for the Mars 2020 mission will be presented.

Miller, R. A.↗
NASA NTRS2020

X-ray Evolution of the Nova V959 Mon Suggests a Delayed Ejection and a Non-Radiative Shock

X-ray observations of shocked gas in novae can provide a useful probe of the dynamics of the ejecta. Here we report on X-ray observations of the nova V959 Mon, which was also detected in GeV gamma-rays with the Fermi satellite. We find that the X-ray spectra are consistent with a two-temperature plasma model with non-solar abundances. We interpret the X-rays as due to shock interaction between the slow equatorial torus and the fast polar outflow that were inferred from radio observations ofV959 Mon. We further propose that the hotter component, responsible for most of the flux, is from the reverse shock driven into the fast outflow. We find a systematic drop in the column density of the absorber between days 60 and 140, consistent with the expectations for such a picture. We present intriguing evidence for a delay of around 40 d in the expulsion of the ejecta from the central binary. Moreover, we infer a relatively small (a few times 10−6Msun) ejecta mass ahead of the shock, considerably lower than the mass of 104K gas inferred from radio observations. Finally, we infer that the dominant X-ray shock was likely notradiative at the time of our observations, and that the shock power was considerably higher than the observed X-ray luminosity. It is unclear why high X-ray luminosity, closer to the inferred shock power, is never seen in novae at early times, when the shock is expected to have high enough density to be radiative.

Thomas Nelson↗
NASA NTRS2020

Magnetospheric Multiscale Observations of Earth’s Oblique Bow Shock Reformation by Foreshock Ultra-Low Frequency Waves

Collisionless shocks can be nonstationary with periodic reformation shown in many simulation results, but direct observations are still tenuous and difficult to conclusively interpret. In this study, using Magnetospheric Multiscale (MMS) observations, we report direct observational evidence of Earth’s oblique bow shock reformation driven by the foreshock ultra-low frequency (ULF)waves. When the four MMS spacecraft were in a string-of-pearls formation roughly along the bow shock normal, they observed that when each period of foreshock ULF waves encountered the bow shock, a new shock ramp formed. Meanwhile, in the magnetosheath, the old bow shock’s remnants were observed periodically convecting downstream. We propose that the reformation mechanism of the oblique bow shock is the variation of the upstream conditions by the periodic ULF waves as they encounter the bow shock. We also examine the nature of reflected ions during the reformation process.

Terry Z Liu↗
NASA NTRSDate not supplied

Zircon (U-Th)/He Impact Crater Thermochronometry and the Effects of Shock Microstructures on Helium Diffusion Kinetics

Absolute age determination of impact cratering events remains difficult and often controversial; a challenge that has resulted in < 50% of known impact craters regarded as accurately and precisely dated. Besides conventional 40Ar/39Ar and U-Pb methods, zircon (U-Th)/He (ZHe) dating of impactites has been applied to large- to medium-sized impact structures. ZHe dates can be fully reset in minutes at 1000°C, which is commonly reached in central sections of the melt sheet, whereas resetting of ZHe at <300°C, which might be encountered near the crater margins or persist in post-impact hydrothermal systems, may take >103-6 years. There is a critical need to quantify the effects of shock-induced microstructures and impact metamorphism on helium diffusion kinetics in well-characterized, variably shocked zircon to further establish the reliability of (U-Th)/He for dating impacts. For this purpose, we investigated suevite and impact melt samples from two impact structures, the Chicxulub multi-ring basin and the Ries complex crater, which enables us to compare zircon helium diffusion kinetics from impact structures with differing sizes, ages, and hydrothermal system longevities. Shock microstructures were characterized by backscattered-electron (BSE) imaging prior to diffusion step-heating fractional release experiments using light-bulb furnace with prograde and retrograde incremental 10°C steps from 250°C to 600°C. Afterward, we characterize the diffusion domain sizes and their interconnectivity within the shocked zircon grains using electron backscatter diffraction (EBSD). We find that zircon with few shock microstructures exhibit no significant deviation from helium diffusion kinetics of undamaged zircon. In contrast, zircon grains with planar deformation features (PDFs) and granular textures classified by BSE and EBSD are characterized by a dramatic decrease in helium retentivity, similar to radiation damaged grains, due to a reduction in the effective domain size and the introduction of interconnected fast diffusion pathways created by shock microstructures. A subset of grains were dated by ZHe after the external morphology of the grains was determined by BSE imaging. The euhedral grains yielded a weighted mean age within the uncertainty of the accepted impact ages, whereas the grains with PDFs or granular textures exhibited younger ages. Thus, these diffusion experiments and ZHe dates suggest that the dramatic decrease in domain size likely renders shocked grains more susceptible to impact-induced hydrothermal resetting and subsequent overprinting. Hence, characterization of shock microstructures is critical for determining accurate impact ages using ZHe methods especially when applied to previously unconstrained craters. The thermochronometer also offers the opportunity to determine the magnitude and duration of post-impact hydrothermal circulation.

zircon↗
NASA NTRS2022

An Investigation Into Some Important Aspects of Droplet Breakup/Vaporization Behavior Caused By a Gas Flow (Both Nonreacting & Detonative Combustion) in a Shock Tube

As a part of the rotating detonation engine (RDE) technology enablement project at NASA Glenn Research Center (GRC), an effort was undertaken to extend our current computational capabilities of OpenNCC in some important ways with the implementation of a modeling approach to account for the droplet breakup caused by a shock-induced gas motion, & a vaporization model valid over a wide range of pressure conditions encountered in multiphase detonation. With the modified code, a study was undertaken to investigate the individual droplet behavior followed by the passage of a shock front. The study is carried out by tracking a sparse group of droplets to gain some understanding of shock induced droplet behavior under various shock strengths & fuel injector conditions. The study also looks into the effect of randomization involved in determining the droplet breakup outcomes. Over a wide range of sub-critical conditions examined, larger droplets are observed to undergo significant changes in droplet behavior following their breakup. However, smaller droplets (10 µm or less ) remain unaffected by any shock induced breakup. In a follow-on work, we investigated the impact of shock and droplet interaction in a detonation study involving both gaseous as well as gas/liquid (droplet clouds) fuel/air stoichiometric mixtures in a simple 3D shock-tube configuration. The droplet clouds are made up of different initial droplet sizes of either 6, 10, or 30 µm. We also investigated the individual droplet behavior followed by the passage of a detonation front. The results represent conditions that lead to both overdriven and C-J (ChapmanJouguet) detonations. Under both test conditions, most of the droplet vaporization is completed within a short distance (duration) behind the detonation front & well within the region of complete combustion observed in a corresponding equivalent gas-phase fuel/air mixture. The impact of the shock-induced droplet breakup is found to be significant in the calculations involving the 30-µm droplets. Subsequent to the breakup, the drop sizes vary from 1 to 10 µm. Another factor that contributed to the observed rapid vaporization is the result of vaporization taking place under supercritical conditions. The overall detonation properties of various droplet clouds (made up of different initial sizes) are similar to those observed in a corresponding gaseous fuel/air mixture. In the calculation involving a gaseous fuel, the calculated C-J detonation velocity is 1822 m/s involving Jet-A/air and φ = 1. In the overdriven detonation, it is 2044 m/s. In the calculations involving droplet clouds, the corresponding detonation velocities are lower. The impact of increased droplet size is primarily seen in a higher reduction in the detonation velocity.

shock/droplet interaction↗
NASA NTRS2022

The Extended Field-Aligned Suprathermal Proton Beam and Long-Lasting Trapped Energetic Particle Population Observed Upstream of a Transient Interplanetary Shock

The properties of the suprathermal particle distributions observed upstream of interplanetary shocks depend not only on the properties of the shocks but also on the transport conditions encountered by the particles as they propagate away from the shocks. The confinement of particles in close proximity to the shocks, as well as particle scattering processes during propagation to the spacecraft, lead to the common observation of upstream diffuse particle distributions. We present observations of a rare extended anisotropic low-energy(≤30 keV)proton beam together with a trapped ≥500 keV proton population observed in association with the arrival of an oblique interplanetary shock at the Advanced Composition Explorer, the Interplanetary Monitoring Platform-8, and the Wind spacecraft on 2001 January 31. Continuous injection of particles by the traveling shock into a smooth radial magnetic field region formed in the tail of a modest high-speed solar wind stream produced an extended foreshock region of energetic particles. The absence of enhanced magnetic field fluctuations upstream of the shock results in the observation of a prolonged anisotropic field-aligned beam of ≤30 keV protons as well as a population of higher-energy(≥500 keV)protons with small pitch-angle cosine(μ∼0)extending far from the shock.

D. Lario↗
NASA NTRSDate not supplied

Implicit Large-Eddy Simulation of Axisymmetric and Asymmetric Turbulent Shock Wave Boundary Layer Interaction at Mach 2.5

High-speed inlet flows feature turbulent shock wave boundary layer interactions at supersonic Mach numbers. A canonical Mach 2.5 shock wave boundary layer interaction experiment at the NASA Glenn Research Center was designed for axisymmetric and asymmetric interactions depending on the placement of the shock generator. Four different wall-resolved implicit large-eddy simulations of the experiment were carried out at one-tenth of the Reynolds number of the experiment: (1) Two axisymmetric interactions with shock generator radii of 9.2% and 14.7% of the test section diameter and (2) two asymmetric interactions with 14.7% shock generator radius and shock generator centerline offsets of one-third and one-fourth of the test section radius. For the two axisymmetric interactions, various sub-domains upstream of the interaction, within the interaction, and downstream of the interaction were analyzed with the proper orthogonal decomposition. For all three regions, the dominant modes are streamwise coherent structures with Strouhal numbers based on separation length of 0.15 and below. Fourier spectra of the wall-pressure coefficient for the interaction region have peaks at Strouhal numbers based on separation length of 0.15 and 0.04. For the two asymmetric interactions, the streamwise extent of the separation varies considerably in the azimuthal direction. For both shock generator offsets, the interaction is strongest for the azimuthal location that is closest to the shock generator. Compared to the axisymmetric interaction, the streamwise extent of the separated flow region is noticeably reduced and a substantial crossflow is observed. Fourier analyses of the unsteady flow indicate low-frequency content f or the upstream part of the interaction which is more two-dimensional. Away from the upstream region, with increasing sweep angle and crossflow, the low-frequency content is diminished.

Computational Fluid Dynamics↗
DOE OSTI2015

Kinetic Simulations of Particle Acceleration at Shocks

Collisionless shocks are mediated by collective electromagnetic interactions and are sources of non-thermal particles and emission. The full particle-in-cell approach and a hybrid approach are sketched, simulations of collisionless shocks are shown using a multicolor presentation. Results for SN 1006, a case involving ion acceleration and B field amplification where the shock is parallel, are shown. Electron acceleration takes place in planetary bow shocks and galaxy clusters. It is concluded that acceleration at shocks can be efficient: >15%; CRs amplify B field via streaming instability; ion DSA is efficient at parallel, strong shocks; ions are injected via reflection and shock drift acceleration; and electron DSA is efficient at oblique shocks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Shock waves in a Z-pinch and the formation of high energy density plasma

A Z-pinch liner, imploding onto a target plasma, evolves in a step-wise manner, producing a stable, magneto-inertial, high-energy-density plasma compression. The typical configuration is a cylindrical, high-atomic-number liner imploding onto a low-atomic-number target. The parameters for a terawatt-class machine (e.g., Zebra at the University of Nevada, Reno, Nevada Terawatt Facility) have been simulated. The 2-1/2 D MHD code, MACH2, was used to study this configuration. The requirements are for an initial radius of a few mm for stable implosion; the material densities properly distributed, so that the target is effectively heated initially by shock heating and finally by adiabatic compression; and the liner's thickness adjusted to promote radial current transport and subsequent current amplification in the target. Since the shock velocity is smaller in the liner, than in the target, a stable-shock forms at the interface, allowing the central load to accelerate magnetically and inertially, producing a magneto-inertial implosion and high-energy density plasma. Comparing the implosion dynamics of a low-Z target with those of a high-Z target demonstrates the role of shock waves in terms of compression and heating. In the case of a high-Z target, the shock wave does not play a significant heating role. The shock waves carry current and transport the magnetic field, producing a high density on-axis, at relatively low temperature. Whereas, in the case of a low-Z target, the fast moving shock wave preheats the target during the initial implosion phase, and the later adiabatic compression further heats the target to very high energy density. As a result, the compression ratio required for heating the low-Z plasma to very high energy densities is greatly reduced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Interaction of a weak shock wave with a discontinuous heavy-gas cylinder

The interaction between a cylindrical inhomogeneity and a weak planar shock wave is investigated experimentally and numerically, and special attention is given to the wave patterns and vortex dynamics in this scenario. A soap-film technique is realized to generate a well-controlled discontinuous cylinder (SF{sub 6} surrounded by air) with no supports or wires in the shock-tube experiment. The symmetric evolving interfaces and few disturbance waves are observed in a high-speed schlieren photography. Numerical simulations are also carried out for a detailed analysis. The refracted shock wave inside the cylinder is perturbed by the diffracted shock waves and divided into three branches. When these shock branches collide, the shock focusing occurs. A nonlinear model is then proposed to elucidate effects of the wave patterns on the evolution of the cylinder. A distinct vortex pair is gradually developing during the shock-cylinder interaction. The numerical results show that a low pressure region appears at the vortex core. Subsequently, the ambient fluid is entrained into the vortices which are expanding at the same time. Based on the relation between the vortex motion and the circulation, several theoretical models of circulation in the literature are then checked by the experimental and numerical results. Most of these theoretical circulation models provide a reasonably good prediction of the vortex motion in the present configuration.

42 ENGINEERING↗
DOE OSTI2021

Insight into the chemistry of TNT during shock compression through ultrafast absorption spectroscopies

Thin films of trinitrotoluene (TNT) were shock compressed using the ultrafast laser shock apparatus at Los Alamos National Laboratory. Visible (VIS) and mid-infrared (MIR) transient absorption spectroscopies were simultaneously performed to probe for electronic and vibrational changes during shock compression of TNT. Additionally, three shock pressures (16 GPa, 33 GPa, and 45 GPa) were selected to observe no reaction, incipient reaction, and strongly developed reactions for TNT within the experimental time scale of <250 ps. Negligible absorption changes in MIR or VIS absorptions were observed at 16 GPa. At 33 GPa, MIR absorptions in the 3000 cm -1 –4000 cm -1 range were observed to increase during the shock and continue to increase during the rarefaction, in contrast to the VIS absorption measurements, which increased during the shock and almost fully recovered during rarefaction. At 45 GPa, both VIS and MIR absorptions were strong and irreversible. The intense and spectrally broad MIR absorptions were attributed to short lived intermediates with strong, spectrally broad absorptions that dominate the spectral response. The MIR and VIS absorption changes observed at 33 GPa and 45 GPa were credited to shock induced chemistry, most likely including the formation of a very broad hydrogenic stretch feature. The results from these experiments are consistent with the chemical mechanisms that include O–H or N–H formation such as CH 3 oxidation or C–N homolysis.

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

Perpendicular Subcritical Shock Structure in a Collisional Plasma Experiment

We present a study of perpendicular subcritical shocks in a collisional laboratory plasma. Shocks are produced by placing obstacles into the supermagnetosonic outflow from an inverse wire array z pinch. We demonstrate the existence of subcritical shocks in this regime and find that secondary shocks form in the downstream. Detailed measurements of the subcritical shock structure confirm the absence of a hydrodynamic jump. We calculate the classical (Spitzer) resistive diffusion length and show that it is approximately equal to the shock width. As a result, we measure little heating across the shock (<10 % of the ion kinetic energy) which is consistent with an absence of viscous dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Electron Energy Partition across Interplanetary Shocks. III. Analysis

An analysis of model fit results of 15,210 electron velocity distribution functions (VDFs), observed within ±2 hr of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 au, is presented as the third and final part on electron VDFs near IP shocks. The core electrons and protons dominate in the magnitude and change in the partialto-total thermal pressure ratio, with the core electrons often gaining as much or more than the protons. Only a moderate positive correlation is observed between the electron temperature and the kinetic energy change across the shock, while weaker, if any, correlations were found with any other macroscopic shock parameter. No VDF parameter correlated with the shock normal angle. The electron VDF evolves from a narrowly peaked core with flaring suprathermal tails in the upstream to either a slightly hotter core with steeper tails or much hotter flattop core with even steeper tails downstream of the weaker and strongest shocks, respectively. Both quasi-static and fluctuating fields are examined as possible mechanisms modifying the VDF, but neither is sufficient alone. For instance, flattop VDFs can be generated by nonlinear ion acoustic wave stochastic acceleration (i.e., inelastic collisions), while other work suggested they result from the combination of quasi-static and fluctuating fields. This three-part study shows that not only are these systems not thermodynamic in nature; even kinetic models may require modification to include things like inelastic collision operators to properly model electron VDF evolution across shocks or in the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2020

Electron Energy Partition Across Interplanetary Shocks. III. Analysis

An analysis of model fit results of 15,210 electron velocity distribution functions (VDFs), observed within ±2 hr of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 au, is presented as the third and final part on electron VDFs near IP shocks. The core electrons and protons dominate in the magnitude and change in the partialto-total thermal pressure ratio, with the core electrons often gaining as much or more than the protons. Only a moderate positive correlation is observed between the electron temperature and the kinetic energy change across the shock, while weaker, if any, correlations were found with any other macroscopic shock parameter. No VDF parameter correlated with the shock normal angle. The electron VDF evolves from a narrowly peaked core with flaring suprathermal tails in the upstream to either a slightly hotter core with steeper tails or much hotter flattop core with even steeper tails downstream of the weaker and strongest shocks, respectively. Both quasi-static and fluctuating fields are examined as possible mechanisms modifying the VDF, but neither is sufficient alone. For instance, flattop VDFs can be generated by nonlinear ion acoustic wave stochastic acceleration (i.e., inelastic collisions), while other work suggested they result from the combination of quasi-static and fluctuating fields. This three-part study shows that not only are these systems not thermodynamic in nature; even kinetic models may require modification to include things like inelastic collision operators to properly model electron VDF evolution across shocks or in the solar wind.

Solar wind↗
NASA NTRS2023

Geoeffectiveness of Interplanetary Shocks Controlled By Impact Angles: Past Research, Recent Advancements, and Future Work

Interplanetary shocks are disturbances commonly observed in the solar wind. IP shock impacts can cause a myriad of space weather effects in the Earth’s magnetopause, inner magnetosphere, ionosphere, thermosphere, and ground magnetic field. The shock impact angle, measured as the angle the shock normal vector performs with the Sun-Earth line, has been shown to be a very important parameter that controls shock geoeffectivess. An extensive review provided by Oliveira and Samsonov (2018) summarized all the work known at the time with respect to shock impact angles and geomagnetic activity; however, this topic has had some progress since Oliveira and Samsonov (2018) and the main goal of this mini review is to summarize all achievements to date in the topic to the knowledge of the author. Finally, this mini review also brings a few suggestions and ideas for future research in the area of IP shock impact angle geoeffectiveness.

interplanetary shocks↗
DOE OSTI2025

Thermal shock resistance of lightweight cements developed for geothermal conditions

Cements are a critical component in well construction, as they act to prevent well fluid and gas escape, prevent corrosion of the casing, and strengthen the wellbore to prevent deformation. Under the high temperature/pressure conditions common in geothermal systems, the injection of cold water for energy production is expected to induce cyclic damage to the borehole cement through the rapid temperature fluctuations. These “thermal shocks” are expected to cause casing shrinkage, annulus formation, and cement tensile stresses. To understand the effect of cold water injection on the wellbore environment, a set of rock-cement-steel samples were created to simulate the structure of a geothermal well. Lightweight thermally-insulating cement blends were tested under thermal shock conditions in this study. In each test, the samples were pressurized to an effective pressure of ~3.5 MPa and placed at high temperatures. Thermal shocks were performed by injecting cold water (~10-15 °C) through the samples at a constant rate while keeping the samples at high temperatures until the sample temperature stopped decreasing and deformation ceased. Eight thermal shock tests were conducted with each sample – two at 100 °C and six at 200 °C. Post-tests analysis was then conducted by cutting open each sample to examine the damage in each component of the simulated wellbore. Experimental results suggest that all samples experienced similar degrees of axial and lateral contraction during cold water injection, but for the most part this contraction is recoverable when injection halts. Post-test analysis revealed that fly ash cenosphere pre-treatment had the best effect on improving thermal shock resistance in the cement blends. Thermomechanical modeling of likely stress paths experienced by the cements during heating/cooling cycles shows that elasto-plastic cement constitutive behavior results in most plastic strain occurring during the initial heating steps, with mostly elastic strain occurring during the thermal shock cycles. In conclusion, this agrees with experimental evidence, suggesting that cement damage from shocking occurs via other mechanisms such as chemical alteration, corrosion, and fatigue.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A space-time smooth artificial viscosity method with wavelet noise indicator and shock collision scheme, Part 2: The 2-D case

This is the second part to our companion paper. Herein, we generalize to two space dimensions the C-method developed in for adding localized, space-time smooth artificial viscosity to nonlinear systems of conservation laws that propagate shock waves, rarefaction waves, and contact discontinuities. For gas dynamics, the C-method couples the Euler equations to scalar reaction-diffusion equations, which we call C-equations, whose solutions serve as space-time smooth artificial viscosity indicators for shocks and contacts. We develop a high-order numerical algorithm for gas dynamics in 2-D which can accurately simulate the Rayleigh-Taylor (RT) instability with Kelvin-Helmholtz (KH) roll-up of the contact discontinuity, as well as shock collision and bounce-back. Solutions to our C-equations not only indicate the location of the shocks and contacts, but also track the geometry of the evolving fronts. This allows us to implement both directionally isotropic and anisotropic artificial viscosity schemes, the latter adding diffusion only in directions tangential to the evolving front. We additionally produce a novel shock collision indicator function, which naturally activates during shock collision, and then smoothly deactivates. Moreover, we implement a high-frequency 2-D wavelet-based noise detector together with an efficient and localized noise removal algorithm. To test the methodology, we use a highly simplified WENO-based discretization scheme. We provide numerical results for some classical 2-D test problems, including the RT problem, the Noh problem, a circular explosion problem from the Liska & Wendroff [13] review paper, the Sedov blast wave problem, the double Mach 10 reflection test, and a shock-wall collision problem. In particular, we show that our artificial viscosity method can eliminate the wall-heating phenomenon for the Noh problem, and thereby produce an accurate, non-oscillatory solution, even though our simplified WENO-type scheme fails to run for this problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Dynamic crushing of metal lattice metamaterials: Shock mode diagrams and transition to topology-independent compaction regime

Additively manufactured lattice metamaterials offer design versatility in strength and energy absorption and provide an additional degree of freedom through the selection of the lattice topology. Under quasistatic loading, the unit cell structure can strongly affect the stiffness, yield, and post-yield behavior, but whether and to what degree the effect of lattice topology persists into dynamic loading scenarios, up to the compaction shock regime, has not been established. LLNL’ s ALE3D hydrocode was used to perform a computational investigation of dynamic loading in multiple lattice types, including the gyroid, octet, Schwarz D, and rhombic dodecahedron, under impact velocities from 0.25 to 2.25 km/s. Shock Hugoniots for each lattice topology are generated and compared, suggesting that above a critical velocity, distinctions between architectures may not persevere and compacted lattices behave similarly. Here, to investigate the transition between topology-dependent quasistatic compression and the topology-independent regime above the critical velocity, a one-dimensional elastic-linear hardening plasticity-densified solid (E-LHP-DS) shock model for lattice materials was developed that relies upon confined compression to link the quasistatic and shock mechanics. Unlike similar works, the model does not assume rigid behavior prior to yield or locking behavior at densification, allowing a richer exploration of lattice mechanics. With only six parameters, the analytical model simultaneously fit quasistatic confined compression simulations for relative densities 0.1 $≤ \bar{ρ} ≤$ 0.9 and predicted dynamic compaction behavior to traverse several distinct shock modes, each defined by a critical impact speed (equivalently, critical stresses). Comparing the numerical results to the one-dimensional E-LHP-DS shock model predictions suggests that the topology-independence under strong shocks is linked to the onset of densification, which can be predicted based on quasistatic confined compression results.

Cellular material↗