skip to main content


Title: Tandem Molecular Dynamics and Continuum Studies of Shock‐Induced Pore Collapse in TATB
Abstract

All‐atom molecular dynamics (MD) and Eulerian continuum simulations, performed using the LAMMPS and SCIMITAR3D codes, respectively, were used to study thermo‐mechanical aspects of the shock‐induced collapse of an initially cylindrical 50 nm diameter pore in single crystals of 1,3,5‐triamino‐2,4,6‐trinitrobenzene (TATB). Three impact speeds, 0.5 km s−1, 1.0 km s−1and 2.0 km s−1, were used to generate the shocks. These impact conditions are expected to yield collapse mechanisms ranging from predominantly visco‐plastic to hydrodynamic. For the MD studies, three crystal orientations (relative to shock‐propagation direction) were examined that span the limiting cases with respect to the crystalline anisotropy in TATB. An isotropic constitutive model was used for the continuum simulations, thus crystal anisotropy is absent. The evolution of spatiotemporally resolved quantities during collapse is reported including local pressure, temperature, pore size and shape, and material flow. Multiple models for the melting curve and specific heat were studied. Within the isotropic elastic/perfectly plastic continuum framework and for the range of impact conditions studied, the specific heat and melting curve sub‐models are shown to have modest effects on the continuum hotspot predictions in the present inert calculation. Treating the MD predictions as ‘ground truth’, albeit with a classical rather than quantum‐like heat capacity, it is clear that – at a minimum – an extension of the constitutive model to account for crystal plasticity and anisotropic strength will be required for a high‐fidelity continuum description.

 
more » « less
NSF-PAR ID:
10133684
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Propellants, Explosives, Pyrotechnics
Volume:
45
Issue:
2
ISSN:
0721-3115
Page Range / eLocation ID:
p. 196-222
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Field investigations of the properties of heavily melted “rotten” Arcticsea ice were carried out on shorefast and drifting ice off the coast ofUtqiaġvik (formerly Barrow), Alaska, during the melt season. While noformal criteria exist to qualify when ice becomes rotten, the objectiveof this study was to sample melting ice at the point at which its structural andoptical properties are sufficiently advanced beyond the peak of the summerseason. Baseline data on the physical (temperature, salinity, density,microstructure) and optical (light scattering) properties of shorefast icewere recorded in May and June 2015. In July of both 2015 and 2017, smallboats were used to access drifting rotten ice within ∼32 km of Utqiaġvik. Measurements showed that pore space increased as icetemperature increased (−8 to 0 C), ice salinitydecreased (10 to 0 ppt), and bulk density decreased (0.9 to0.6 g cm−3). Changes in pore space were characterized with thin-sectionmicrophotography and X-ray micro-computed tomography in the laboratory. Theseanalyses yielded changes in average brine inclusion number density (whichdecreased from 32 to 0.01 mm−3), mean pore size (whichincreased from 80 µm to 3 mm), and total porosity (increased from0 % to > 45 %) and structural anisotropy (variable, withvalues of generally less than 0.7). Additionally, light-scattering coefficientsof the ice increased from approximately 0.06 to > 0.35 cm−1 as the ice melt progressed. Together, these findings indicate thatthe properties of Arctic sea ice at the end of melt season are significantlydistinct from those of often-studied summertime ice. If such rotten ice wereto become more prevalent in a warmer Arctic with longer melt seasons, thiscould have implications for the exchange of fluid and heat at the oceansurface.

     
    more » « less
  2. ABSTRACT

    When collapse of the iron core in a massive red or yellow supergiant does not lead to an energetic supernova, a significant fraction of the convective hydrogen envelope will fall in towards the black hole formed from the collapsing core. The random velocity field in the convective envelope results in finite specific angular momentum in each infalling shell. Using 3D hydrodynamical simulations, we follow the infall of this material to small radii, resolving the circularization radii of the flow. We show that infall of the convective envelope leads to nearly complete envelope ejection in a ≳1048 erg explosion with outflow speeds of ≳200 km s−1. The light curve of such an explosion would show a characteristic, red plateau as the ejecta cools and a hydrogen recombination front recedes through the expanding ejecta. Adopting supernova IIp scalings, the event would have a plateau luminosity of ≳1040 erg s−1 and a duration of several hundreds of days. These events would appear quite similar to luminous red novae with red or yellow supergiant progenitors; some luminous red novae may, in fact, be signposts of black hole formation. The mechanism studied here produces more energetic explosions than the weak shock generated from radiation of neutrino energy during the protoneutron star phase. Because we cannot simulate all the way to the horizon, our results are likely lower limits on the energy and luminosity of transients produced during the collapse of a red or yellow supergiant to form a black hole.

     
    more » « less
  3. Abstract A criterion to predict the onset of disordering under biaxial loading based on a critical potential energy per atom was studied. In contrast to previous theories for disordering, this criterion incorporates the effects of strain rate and strain state. The strain state (or stress state) is defined by the combination of strain (or stress) magnitudes and directions that are applied to each sample during the simulation. The validity of this criterion was studied using molecular dynamic (MD) simulations of Ag conducted over a wide range of biaxial strain rates, strain configurations, and crystal orientations with respect to the applied stress state. Biaxial strains were applied in two different planes, ( 11 2 ¯ ) and (001) in eight directions in each plane. Results showed that, when larger strain rates were applied, there was a transition from plastic deformation driven by the nucleation and propagation of dislocations to disordering and viscous flow. Although the critical strain rate to initiate disorder was found to vary in the range of ε ˙ = 1 × 10 11  s −1 to ε ˙ = 4 × 10 11  s −1 , a consistent minimum PE/atom of −2.7 eV was observed over a broad range of strain states and for both crystallographic orientations that were studied. This indicates that the critical PE/atom is a material property that can be used to predict the onset of disordering under biaxial loading. Further, the results showed that this criterion can be applied successfully even when non-uninform strain states arise in the crystal. 
    more » « less
  4. Abstract

    Impact metamorphic effects from quartz and feldspar and to a lesser extent olivine and pyroxene have been studied in detail. Comparatively, studies documenting shock effects in other minerals, such as double chain inosilicates, phyllosilicates, carbonates, and sulfates, are lacking. In this study, we investigate impact metamorphism recorded in crystalline basement rocks from the Steen River impact structure (SRIS), a 25 km diameter complex crater inNWAlberta, Canada. An array of advanced analytical techniques was used to characterize the breakdown of biotite in two distinct settings: along the margins of localized regions of shock melting and within granitic target rocks entrained as clasts in a breccia. In response to elevated temperature gradients along shock vein margins, biotite transformed at high pressure to an almandine‐Ca/Fe majorite‐rich garnet with a density of 4.2 g cm−3. The shock‐produced garnets are poikilitic, with oxide and silicate glass inclusions. Areas interstitial to garnets are vesiculated, in support of models for the formation of shock veins via oscillatory slip, with deformation continuing during pressure release. Biotite within granitic clasts entrained within the hot breccia matrix thermally decomposed at ambient pressure to produce a fine‐grained mineral assemblage of orthopyroxene + sanidine + titanomagnetite. These minerals are aligned to the (001) cleavage plane of the original crystal. In this and previous work, the transformation of an inosilicate (pargasite) and a phyllosilicate (biotite) to form garnet, an easily identifiable, robust mineral, has been documented. We contend that in deeply eroded astroblemes, high‐pressure minerals that form within or in the environs of shock veins may serve as one of the possibly few surviving indicators of impact metamorphism.

     
    more » « less
  5. Abstract

    We present the discovery and extensive follow-up of a remarkable fast-evolving optical transient, AT 2022aedm, detected by the Asteroid Terrestrial impact Last Alert Survey (ATLAS). In the ATLASoband, AT 2022aedm exhibited a rise time of 9 ± 1 days, reaching a luminous peak withMg≈ −22 mag. It faded by 2 mag in thegband during the next 15 days. These timescales are consistent with other rapidly evolving transients, though the luminosity is extreme. Most surprisingly, the host galaxy is a massive elliptical with negligible current star formation. Radio and X-ray observations rule out a relativistic AT 2018cow–like explosion. A spectrum in the first few days after explosion showed short-lived Heiiemission resembling young core-collapse supernovae, but obvious broad supernova features never developed; later spectra showed only a fast-cooling continuum and narrow, blueshifted absorption lines, possibly arising in a wind withv≈ 2700 km s−1. We identify two further transients in the literature (Dougie in particular, as well as AT 2020bot) that share similarities in their luminosities, timescales, color evolution, and largely featureless spectra and propose that these may constitute a new class of transients: luminous fast coolers. All three events occurred in passive galaxies at offsets of ∼4–10 kpc from the nucleus, posing a challenge for progenitor models involving massive stars or black holes. The light curves and spectra appear to be consistent with shock breakout emission, though this mechanism is usually associated with core-collapse supernovae. The encounter of a star with a stellar-mass black hole may provide a promising alternative explanation.

     
    more » « less