Search NASA⌕ Search

SEARCH · Search NASA

Results for “COLLAPSE”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Protostellar Collapse with a Shock

We reexamine both numerically and analytically the collapse of the singular isothermal sphere in the context of low-mass star formation. We consider the case where the onset of collapse is initiated by some arbitrary process which is accompanied by a central output of either heat or kinetic energy. We find two classes of numerical solutions describing this manner of collapse. The first approaches in time the expansion wave solution of Shu, while the second class is characterized by an ever-decreasing central accretion rate and the presence of an outwardly propagating weak shock. The collapse solution which represents the dividing case between these two classes is determined analytically by a similarity analysis. This solution shares with the expansion wave solution the properties that the gas remains stationary with an r(exp -2) density profile at large radius and that, at small radius, the gas free-falls onto a nascent core at a constant rate which depends only on the isothermal sound speed. This accretion rate is a factor of approx. 0.1 that predicted by the expansion wave solution. This reduction is due in part to the presence of a weak shock which propagates outward at 1.26 times the sound speed. Gas in the postshock region first moves out subsonically but is then decelerated and begins to collapse. The existence of two classes of numerical collapse solutions is explained in terms of the instability to radial perturbations of the analytic solution. Collapse occurring in the manner described by some of our solutions would eventually unbind a finite-sized core. However, this does not constitute a violation of the instability properties of the singular isothermal sphere which is unstable both to collapse and to expansion. To emphasize this, we consider a purely expanding solution for isothermal spheres. This solution is found to be self-similar and results in a uniform density core in the central regions of the gas. Our solutions may be relevant to the 'luminosity' problem of protostellar cores since the predicted central accretion rates are significantly reduced relative to that of the expansion wave solution. Furthermore, our calculations indicate that star-forming cloud cores are not very tightly bound and that modest disturbances can easily result in both termination of infall and dispersal of unaccreted material.

Tsai, John C.↗
NASA NTRS1995

Protostellar Collapse with a Shock

We reexamine both numerically and analytically the collapse of the singular isothermal sphere in the context of low-mass star formation. We consider the case where the onset of collapse is initiated by some arbitrary process which is accompanied by a central output of either heat or kinetic energy. We find two classes of numerical solutions describing this manner of collapse. The first approaches in time the expansion wave solution of Shu, while the second class is characterized by an ever-decreasing central accretion rate and the presence of an outwardly propagating weak shock. The collapse solution which represents the dividing case between these two classes is determined analytically by a similarity analysis. This solution shares with the expansion wave solution the properties that the gas remains stationary with an r(sup -2) density profile at large radius and that, at small radius, the gas free-falls onto a nascent core at a constant rate which depends only on the isothermal sound speed. This accretion rate is a factor of approx. 0.1 that predicted by the expansion wave solution. This reduction is due in part to the presence of a weak shock which propagates outward at 1.26 times the sound speed. Gas in the postshock region first moves out subsonically but is then decelerated and begins to collapse. The existence of two classes of numerical collapse solutions is explained in terms of the instability to radial perturbations of the analytic solution. Collapse occurring in the manner described by some of our solutions would eventually unbind a finite-sized core. However, this does not constitute a violation of the instability properties of the singular isothermal sphere which is unstable both to collapse and to expansion. To emphasize this, we consider a purely expanding solution for isothermal spheres. This solution is found to be self-similar and results in a uniform density core in the central regions of the gas. Our solutions may be relevant to the 'luminosity' problem of protostellar cores since the predicted central accretion rates are significantly reduced relative to that of the expansion wave solution. Furthermore, our calculations indicate that star-forming cloud cores are not very tightly bound and that modest disturbances can easily result in both termination of infall and dispersal of unaccreted material.

Tsai, John C.↗
DOE OSTI2022

Phase Contrast X‐Ray Imaging of the Collapse of an Engineered Void in Single‐Crystal HMX

Abstract Imaging the collapse of a single void that creates a hot spot initiation site in an otherwise defect‐free explosive is challenging given the spatial and temporal scales involved in explosive systems. This work presents our attempt to examine a single hot spot mode (void collapse) in single‐crystal octahydro‐l,3,5,7‐tetranitro‐l,3,5,7‐tetrazocine (HMX) embedded in Sylgard. The hot spot heating mechanisms involved with pore collapse include adiabatic heating, jetting, and viscoplastic dissipation. Quantifying the dynamics of a pore collapse is a crucial step to understanding which mechanisms dominate during ignition events. Our experiments were conducted with a single‐stage, light‐gas gun at Argonne National Laboratory's Advanced Photon Source, applying the phase contrast imaging technique while collecting high‐speed video. The details of HMX single crystal production, defect (pore) engineering, and sample construction, along with experimental results are presented here. These results demonstrate that detailed collapse dynamics can be obtained from homogeneous, single‐crystal explosives with this approach. Qualitative comparisons are made with simulation data which show good agreement in the transition between a quasi‐symmetric pore collapse and an asymmetric collapse with jetting across the pore as measured with normalized pore area and pore circularity.

42 ENGINEERING↗
DOE OSTI2023

Energy localization efficiency in 1,3,5-trinitro-2,4,6-triaminobenzene pore collapse mechanisms

Atomistic and continuum scale modeling efforts have shown that the shock-induced collapse of porosity can occur via a wide range of mechanisms dependent on pore morphology, the shockwave pressure, and material properties. The mechanisms that occur under weaker shocks tend to be more efficient at localizing thermal energy but do not result in high, absolute temperatures or spatially large localizations compared to mechanisms found under strong shock conditions. However, the energetic material 1,3,5-trinitro-2,4,6-triaminobenzene (TATB) undergoes a wide range of collapse mechanisms that are not typical of similar materials, leaving the collapse mechanisms and the resultant energy localization from the collapse, i.e., hotspots, relatively uncharacterized. Therefore, we present the pore collapse simulations of cylindrical pores in TATB for a wide range of pore sizes and shock strengths that trigger viscoplastic collapses that occur almost entirely perpendicular to the shock direction for weak shocks and hydrodynamic-like collapses for strong shocks that do not break the strong hydrogen bonds of the TATB basal planes. The resulting hotspot temperature fields from these mechanisms follow trends that differ considerably from other energetic materials; hence, we compare them under normalized temperature values to assess the relative efficiency of each mechanism to localize energy. The local intra-molecular strain energy of the hotspots is also assessed to better understand the physical mechanisms behind the phenomena that lead to a latent potential energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Dynamic Rupture Simulations of Caldera Collapse Earthquakes: Effects of Wave Radiation, Magma Viscosity, and Evidence of Complex Nucleation at Kilauea 2018

All instrumented basaltic caldera collapses have generated M w > 5 very long period earthquakes. However, previous studies of source dynamics have been limited to lumped models treating the caldera block as rigid, leaving open questions related to how ruptures initiate and propagate around the ring fault, and the seismic expressions of those dynamics. Here, we present the first 3D numerical model capturing the nucleation and propagation of ring fault rupture, the mechanical coupling to the underlying viscoelastic magma, and the associated seismic wavefield. We demonstrate that seismic radiation, neglected in previous models, acts as a damping mechanism reducing coseismic slip by up to half, with effects most pronounced for large magma chamber volume/ring fault radius or highly compliant crust/compressible magma. Viscosity of basaltic magma has negligible effect on collapse dynamics. In contrast, viscosity of silicic magma significantly reduces ring fault slip. We use the model to simulate the 2018 Kilauea caldera collapse. Three stages of collapse, characterized by ring fault rupture initiation and propagation, deceleration of the downward-moving caldera block and magma column, and post-collapse resonant oscillations, in addition to chamber pressurization, are identified in simulated and observed (unfiltered) near-field seismograms. A detailed comparison of simulated and observed displacement waveforms corresponding to collapse earthquakes with hypocenters at various azimuths of the ring fault reveals a complex nucleation phase for earthquakes initiated on the northwest. Our numerical simulation framework will enhance future efforts to reconcile seismic and geodetic observations of caldera collapse with conceptual models of ring fault and magma chamber dynamics.

58 GEOSCIENCES↗
DOE OSTI2023

Modeling shock-induced void collapse in single-crystal Ta systems at the mesoscales

Understanding the role of microstructural heterogeneities on the shock wave propagation and defect evolution behavior is essential to predicting the dynamic response of metals. Heterogeneities, such as voids, provide challenges to understanding the wave propagation behavior as the shock-void interaction can collapse the void and result in large plastic strains and significant heating (hot spot formation) in the metal. Accurate modeling of this phenomenon requires predicting the void collapse mechanisms and the related heat generation and dissipation mechanisms (hotspot formation) that determine the microstructure evolution. While molecular dynamics (MD) simulations can model the void collapse behavior, the time/length scale capabilities pose a challenge to connect with continuum models or the experimental scales. Here, this study presents the capability of the newly developed quasi-coarse-grained dynamics (QCGD) method that extends the MD simulations to larger system sizes and longer times to model this phenomenon. This study uses QCGD simulations to investigate the mechanisms of shock wave interactions with pre-existing voids in single-crystal Ta microstructures for variations in shock pressures, void size, and loading orientations. For a given orientation, the plasticity contributions and rates of void collapse are observed to vary with shock pressures and void size. The larger void sizes and higher pressures result in increased temperatures (hot spots) and faster void collapse rates and unravel the variations in the plasticity contributions. In addition, QCGD simulations investigate the post-collapse microstructure evolution as a release wave travels through the hot spot region. The simulations reveal that the reduction in temperatures due to heat dissipation initiates the dynamic recrystallization behavior in the hotspot regions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Exemption of lattice collapse in Ni–MnO 2 birnessite regulated by the structural water mobility

Lattice collapse and associated mechanical fracture frequently occur in Li-intercalated metal oxide cathodes at the deep charge state upon Li-ion removal, governed by chemical compositions and the resulting electron density of the oxygen atoms. However, similar lattice collapse for metal oxide electrodes in aqueous storage and its mitigation have not been well studied. In this work, we reported the lattice collapse of MnO 2 layered birnessite during the aqueous de-sodiation process at high voltage due to the structural water motion, as evidenced by in situ XRD. Unlike non-aqueous Li-intercalated electrodes, Ni-dopants mitigated the lattice collapse of birnessite at deep charge states. Moreover, density functional theory (DFT) calculations showed that Ni doping induces charge depletion of lattice oxygen due to its higher electronegativity than Mn. This charge reduction, in turn, yields a significant decrease in electrostatic repulsion between the oxygens belonging to lattice and structural water. Classical molecular dynamics simulations based on atomic charges obtained from DFT elucidate that Ni-doping facilitates immobilization of structural water in (Ni)MnO 2 and prevents lattice collapse upon Na-ion removal. (Ni)MnO 2 exempted from lattice collapse shows an improved storage capacity relative to MnO 2 while maintaining similar cycling stability.

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

Three-dimensional Hydrodynamic Simulations of Convective Nuclear Burning in Massive Stars Near Iron Core Collapse

Nonspherical structure in massive stars at the point of iron core collapse can have a qualitative impact on the properties of the ensuing core-collapse supernova explosions and the multimessenger signals they produce. Strong perturbations can aid successful explosions by strengthening turbulence in the postshock region. Here we report on a set of 4π 3D hydrodynamic simulations of O- and Si-shell burning in massive star models of varied initial masses using MESA and the FLASH simulation framework. We evolve four separate 3D models for roughly the final 10 minutes prior to and including iron core collapse. We consider initial 1D MESA models with masses of 14, 20, and 25 M ⊙ to survey a range of O/Si-shell density and compositional configurations. Further, we characterize the convective shells in our 3D models and compare them to the corresponding 1D models. In general, we find that the angle-average convective speeds in our 3D simulations near collapse are three to four times larger than the convective speeds predicted by MESA at the same epoch for our chosen mixing length parameter of α MLT = 1.5. In three of our simulations, we observe significant power in the spherical harmonic decomposition of the radial velocity field at harmonic indices of ℓ = 1–3 near collapse. Our results suggest that large-scale modes are common in massive stars near collapse and should be considered a key aspect of presupernova progenitor models.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2013

Estimating the Collapse Pressure of an Inflatable Aerodynamic Decelerator

The collapse pressure of an inflatable membrane is the minimum differential pressure which will sustain a specific desired shape under an applied load. In this paper, we present a method for estimating the collapse pressure of a tension-cone inflatable aerodynamic decelerator (IAD) that is subject to a static aerodynamic load. The IAD surface is modeled as an elastic membrane. For a given aerodynamic load and sufficiently high torus differential pressure, the IAD assumes a stable axisymmetric equilibrium shape. When the torus pressure is reduced sufficiently, the symmetric equilibrium state becomes unstable and we define this instance to be the critical pressure Pcr. In this paper, we will compare our predicted critical torus pressure with the corresponding observed torus collapse pressure (OTCP) for fifteen tests that were conducted by the third author and his collaborators at the NASA Glenn Research Center 10x10 Supersonic Wind Tunnel in April 2008. One of the difficulties with these types of comparisons is establishing the instance of torus collapse and determining the OTCP from quantities measured during the experiment. In many cases, torus collapse is gradual and the OTCP is not well-defined. However, in eight of the fifteen wind tunnel tests where the OTCP is well-defined, we find that the average of the relative differences (Pcr - OTCP/Pcr) was 8.9%. For completeness, we will also discuss the seven tests where the observed torus collapse pressure is not well-defined.

Baginski, Frank E.↗
NASA NTRS2017

Atmospheric Collapse on Early Mars: The Role of CO2 Clouds

The abundance of evidence that liquid water flowed on the surface early in Mars' history strongly implies that the early Martian atmosphere was significantly more massive than it is today. While it seems clear that the total CO2 inventory was likely substantially larger in the past, the fundamental question about the physical state of that CO2 is not completely understood. Because the temperature at which CO2 condenses increases with surface pressure, surface CO2 ice is more likely to form and persist as the atmospheric mass increases. For the atmosphere to remain stable against collapse, there must be enough energy, distributed planet wide, to stave off the formation of permanent CO2 caps that leads to atmospheric collapse. The presence of a "faint young sun" that was likely about 25 percent less luminous 3.8 billion years ago than the sun today makes this even more difficult. Several physical processes play a role in the ultimate stability of a CO2 atmosphere. The system is regulated by the energy balance between solar insolation, the radiative effects of the atmosphere and its constituents, atmospheric heat transport, heat exchange between the surface and the atmosphere, and latent heating/cooling. Specific considerations in this balance for a given orbital obliquity/eccentricity and atmospheric mass are the albedo of the caps, the dust content of the atmosphere, and the presence of water and/or CO2 clouds. Forget et al. show that, for Mars' current obliquity (in a circular orbit), CO2 atmospheres ranging in surface pressure from 500 hectopascals to 3000 hectopascals would have been stable against collapsing into permanent surface ice reservoirs. Soto et al. examined a similar range in initial surface pressure to investigate atmospheric collapse and to compute collapse rates. CO2 clouds and their radiative effects were included in Forget et al. but they were not included in Soto et al. Here we focus on how CO2 clouds affect the stability of the atmosphere against collapse.

Kahre, M. A.↗
DOE OSTI2021

Core-collapse contamination in photometric samples of Type Ia Supernovae

This is an exciting time for cosmology with type Ia supernovae (SNe Ia). The recentlyconcluded Dark Energy Survey SN programme (DES-SN) has obtained the largest anddeepest high-redshift cosmological SN Ia sample, and the Vera Rubin Observatory isexpected to observe at least one order of magnitude more SNe Ia in the next decade.In both these experiments, only a limited fraction (.10 per cent) of the SNe can bespectroscopically classified. This leaves us with large ‘photometric’ SN samples, withthe potential for significant contamination by core-collapse SNe that may bias SN Iacosmological measurements. This thesis demonstrates how this contamination can bemodelled and accounted for in current and future cosmological analyses.First, I present state-of-the-art simulations of the SN universe. These are designedto accurately model the population of SNe Ia, peculiar SNe Ia and core-collapse SNe, aswell as their host galaxies. To improve the diversity and quality of the simulated corecollapseSNe, I build a new library of core-collapse SN templates using spectroscopicand photometric (optical and near-ultraviolet) data of 67 core-collapse SNe from theliterature. I account for our incomplete knowledge of core-collapse SN properties bygenerating a set of SN simulations (rather than a single one), each exploring differentmodelling choices and template libraries. I then characterise selection effects in theDES-SN survey and incorporate them in the simulations, thus obtaining a series ofDES-like simulated SN samples that can be compared to the observed DES-SN data.The agreement between the simulations and data is excellent across many observed SNproperties, including Hubble residuals. These simulations are the first to reproduce theobserved photometric SN and host galaxy properties in high-redshift surveys with no fine-tuning of the input parameters.I use my simulation framework to train and test the performance of SuperNNova,a photometric SN classifier based on recurrent neural networks. I explore differenttraining and validation strategies and show that, across all the DES-SN simulationstested, SuperNNova reduces core-collapse SN contamination to 0.8–3.5 per cent. Ithen show that biases due to contamination on the equation-of-state of dark energy,w, are < 0:008 when using our reference SuperNNova model. This compares to anexpected statistical uncertainty on w from the DES-SN sample of 0:039, thus showing that contamination is not a limiting systematic for the cosmological analysis of theDES-SN sample.The results presented in this thesis are the foundation of the DES SN Ia cosmologicalanalysis; they also provide important implications for the future of SN cosmology,as they demonstrate that contamination is not expected to significantly degrade thecosmological figure of merit of the Rubin SN Ia analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Embedding condition describing shearing, dissipative collapse

Highlights: • Perturbative approach to the time-dependent Karmarkar condition in shearing collapse. • Analysis of the (in)stability of the collapsing fluid using the TOV equation. • Study of the temperature profiles using a causal heat transport equation. In this work we employ the time-dependent Karmarkar condition to model dissipative gravitational collapse in the presence of shear which leads to a final static configuration. We believe that this is a first attempt at using the time-dependent Karmarkar condition describing a radiating star with shear in which the collapsing core dissipates energy in the form of a radial heat flux. The final static configuration which is described by an embedding class I spacetime was obtained by Hansraj and Moodly who showed that this solution satisfies all the physical requirements for a static compact object. As the collapse proceeds and the heat flux tends to zero the dynamical core evolves into a static compact core with the exterior spacetime being described by the vacuum Schwarzschild solution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Storage-Induced Collapse of Lignin Macromolecular Structure and Its Impacts on the Biorefinery

Lignin plays a vital role in the economics of biorefineries, serving as a source of process energy and a feedstock for sustainable fuels and chemical production. While understanding lignin’s chemical composition is crucial, emerging evidence suggests that a more comprehensive understanding of its macromolecular structure is critical to explaining its complex behavior in the biorefinery. This study investigated the partial collapse of the lignin network in corn stover feedstock after harvest and storage as a result of the microbial digestion of hemicellulose. Fluorescence microscopy was used to detect the collapse of lignin in terms of lignin’s inter-molecular interaction and the re-orientation of lignin’s chromophores, by the changes in lignin’s fluorescence lifetime, anisotropy, and the number of effective emitters. With minimal sample perturbation, our in-situ microscopic results revealed lignin's coil-globule transition phenomena, which was only previously predicted by molecular dynamics modeling extracted lignin in solvent. This collapse of lignin macromolecular structure was confirmed by results from NMR, IR, Raman, and powder X-ray diffraction. We also investigated the impact of this storage-induced collapse on the downstream biorefinery processes. Our study revealed that the two major approaches for lignin valorization in the lignin-first biorefinery model, namely monomer extraction and milled wood lignin extraction, were negatively impacted by the lignin collapse. As changes during storage are a source of feedstock variability, our study highlights the importance of understanding the effect of feedstock handling on biorefinery operations and economics.

09 BIOMASS FUELS↗
DOE OSTI2022

Dynamics of an inertially collapsing gas bubble between two parallel, rigid walls

The collapse of cavitation bubbles in channel flows can give rise to structural damage along neighbouring walls. Although the collapse of a bubble near a single wall has been studied extensively, less is known about bubble collapse between two walls, e.g. as in a channel. We conduct highly resolved, direct simulations of the Navier–Stokes equations to investigate the bubble dynamics and pressures produced by the collapse of a bubble between two parallel rigid walls. We examine the dependence of the dynamics and pressures on the initial bubble location, confinement and driving pressure. For a fixed initial stand-off distance, as the channel width increases the bubble volume, migration distance and re-entrant jet speed approach their single-wall counterparts. We obtain an expression for the minimum channel width at which the confinement does not affect the bubble dynamics depending on the driving pressure difference and initial stand-off distance. For a fixed channel width, varying stand-off distance reduced the maximum wall pressures in the channel relative to the single wall; the trend was consistent for three different driving pressures. Two different jetting behaviours are seen when the bubble is centred in the channel, depending on the channel width. Under significant confinement, wall-parallel re-entrant jets impinge upon each other and further intensify the collapse of the vortex ring.

Rodriguez, Jr, Mauro (ORCID:0000000305450265)↗
DOE OSTI2025

Accurately simulating core-collapse self-interacting dark matter halos

The properties of satellite halos provide a promising probe for dark matter (DM) physics. Observations have motivated current efforts to explain surprisingly compact DM halos. If DM is not collisionless, but has strong self-interactions, halos can undergo gravothermal collapse, leading to higher densities in the central region of the halo. However, it is challenging to model this collapse phase from first principles. To improve on this, we sought to better understand the numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations in the collapse phase. Especially, our aim was to better understand the evolution of satellite halos. To do so, we ran SIDM N-body simulations of a low-mass halo in isolation and within an external gravitational potential. The simulation set-up was motivated by the perturber of the stellar stream GD-1. We find that the halo evolution is very sensitive to energy conservation errors, and a SIDM kernel size that is too large can artificially speed up the collapse. Moreover, we demonstrate that the King model can describe the density profile at small radii for the late stages that we have simulated. Furthermore, for our most highly resolved simulation (N = 5 × 10 7 ) we have made the data public. It can serve as a benchmark. Overall, we find that the current numerical methods do not suffer from convergence problems in the late collapse phase and provide guidance on how to choose numerical parameters, for example that the energy conservation error is better kept well below 1%. This allows simulations to be run of halos that become concentrated enough to explain observations of GD-1-like stellar streams or strong gravitational lensing systems.

dark matter↗
DOE OSTI2022

Sudden collapse of a pressure profile generated by off-axis heating in a linear magnetized plasma

Here, the features of an unexpected, large event that arises spontaneously during a basic heat transport experiment are presented. It consists of the sudden collapse of the radial plasma pressure profile, akin to disruption events observed in toroidal magnetic confinement devices. The experiment is performed on the Large Plasma Device at the University of California, Los Angeles (UCLA). It uses a LaB6 thermionic emitter of annular shape to induce off-axis heating of a cold, afterglow plasma, in a linear magnetic geometry. The temporal evolution consists of three regimes. During an early, quiescent period, classical heat transport along and across the magnetic field arises from Coulomb collisions. After significant pressure gradients develop, drift-Alfvén waves become unstable. Upon reaching large amplitude, they trigger avalanche events that flatten the outer part of the heated region, which, in turn, quenches the instability. Due to the sustained heating, the pressure profile rebuilds and the process repeats, leading to a relatively long, second regime that displays multiple avalanches, but suddenly, the annular pressure profile is observed to collapse. After this collapse, the system enters a third regime with large fluctuations. Before the collapse, a rapid, runaway heating environment arises whose time evolution exhibits a self-similar dependence on the applied voltage. The time evolution, morphology, and scaling of the collapse event are presented, and an examination is made of the underlying mechanisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Mechanics of pore array collapse and interaction in shock-compressed polymethyl methacrylate (PMMA)

Recent studies on dynamic pore collapse have revealed significant development of shear localization, which can lead to material failure in porous structures and hot spot generation in energetic materials. These findings have dramatically improved the understanding of failure mechanisms during pore collapse but also prompt further investigation of realistic porous materials. In particular, porous media consist of many pores and porous networks. Even in low-porosity materials, pores can form in close proximity during the manufacturing process, leading to the critical question of pore–pore interaction during collapse under dynamic loading conditions. This study investigates, via plate impact experiments coupled with high-speed internal digital image correlation and shadowgraphy techniques, the collapse of two pores in shock-compressed PMMA at stresses between 0.4 and 1 GPa. The results of these experiments provide new insights into shear localization in pore collapse, in addition to distinct interactions between pores. Shadowgraphy measurements reveal novel, direct visualization of shear band development and crack evolution from pore surfaces. Spacing between adiabatic shear bands is measured over a range of impact stresses and is predicted accurately by the Grady–Kipp model. Pore interactions are found to effect a transition in the impact stress threshold at which different failure mechanisms initiate and are also found to possibly influence preferential sites for shear cracking. Throughout the study, numerical and theoretical models are leveraged to understand shear localization behavior. The role of baroclinicity and wave interactions between the pores is used to elucidate interaction mechanisms between pores.

Lawlor, Barry P. [California Institute of Technolo↗
DOE OSTI2024

Numerical challenges in modeling gravothermal collapse in Self-Interacting Dark Matter halos

When dark matter has a large cross section for self scattering, halos can undergo a process known as gravothermal core collapse, where the inner core rapidly increases in density and temperature. To date, several methods have been used to implement Self-Interacting Dark Matter (SIDM) in N-body codes, but there has been no systematic study of these different methods or their accuracy in the core-collapse phase. In this paper, we compare three different numerical implementations of SIDM, including the standard methods from the GIZMO and Arepo codes, by simulating idealized dwarf halos undergoing significant dark matter self interactions (σ/m = 50 cm 2 /g). When simulating these halos, we also vary the mass resolution, time-stepping criteria, and gravitational force-softening scheme. The various SIDM methods lead to distinct differences in a halo's evolution during the core-collapse phase, as each results in spurious scattering rate differences and energy gains/losses. The use of adaptive force softening for gravity can lead to numerical heating that artificially accelerates core collapse, while an insufficiently small simulation time step can cause core evolution to stall or completely reverse. Additionally, particle numbers must be large enough to ensure that the simulated halos are not sensitive to noise in the initial conditions. Even for the highest-resolution simulations tested in this study (10 6 particles per halo), we find that variations of order 10% in collapse time are still present. Here, the results of this work underscore the sensitivity of SIDM modeling on the choice of numerical implementation and motivate a careful study of how these results generalize to halos in a cosmological context.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗