SEARCH · Search NASA

Results for “finite element simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

739 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Realization of two-dimensional discrete time crystals with anisotropic Heisenberg coupling

A discrete time crystal (DTC) is an out-of-equilibrium phase of matter that spontaneously breaks discrete time-translation symmetry. Previous studies have been limited to a set of models with Ising-like couplings - and mostly only in one dimension - thus precluding our understanding of the existence (or not) of DTCs in models with more realistic interactions. In this work, by combining the latest generation of IBM quantum processors with state-of-the-art tensor network methods, we demonstrate the existence of a DTC in a two-dimensional system governed by anisotropic Heisenberg interactions. We uncover a rich phase diagram encompassing spin-glass, ergodic, and time-crystalline phases, and identify the interplay of initialization, interaction anisotropy, and driving protocols in stabilizing the DTC phase. By extending the study of Floquet matter beyond simplified models, we lay the groundwork for exploring how driven systems bridge the gap between quantum coherence and emergent non-equilibrium thermodynamics.

Phase transitions and critical phenomena↗

NASA NTRS2026

Comparison Testing and Analysis of Braided Cords Constructed from Varying Kevlar Fiber Types

Material advancements move quickly, resulting in material obsolescence issues that require parachute engineers to understand how the material properties change with a material upgrade and if there are expected impacts to overall performance of the parachute. The objective of this test was to measure the elongation of different Kevlar types at both the cord-level and at the yarn-level to determine if a change in Kevlar type impacted parachute performance. An Instron machine with a pin-to-pin connection was used to measure load and displacement time histories for various Kevlar cord and yarn samples in a cyclic load profile. Twenty-four cord samples and twenty-four yarn samples were tested; each sample was cycled five times to the relevant maximum load (either to breaking strength or to the 30% peak load expected in flight). In addition to the Instron measurements, a photogrammetric technique was also implemented to measure elongation of the samples using button targets. Modulus test data shows approximately a 10% difference in elongation between Kevlar 29AP and Kevlar 129, and another 10% difference in elongation between Kevlar 29 and Kevlar 29AP. Although this is a measurable change in modulus, the effect on parachute performance is likely small, particularly for non-primary structure elements of relatively small length. This test campaign provides an example of the type of testing and analysis that parachute engineers could do to understand how material upgrades affect performance of their systems.

Parachutes↗

NASA NTRSDate not supplied

A GPU-based Approach for Turbomachinery Application

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Gas Turbine↗

DOE OSTI2023

Delineating the Impact of Transition‐Metal Crossover on Solid‐Electrolyte Interphase Formation with Ion Mass Spectrometry

Abstract Lithium‐metal batteries (LMB) employing cobalt‐free layered‐oxide cathodes are a sustainable path forward to achieving high energy densities, but these cathodes exhibit substantial transition‐metal dissolution during high‐voltage cycling. While transition‐metal crossover is recognized to disrupt solid‐electrolyte interphase (SEI) formation on graphite anodes, experimental evidence is necessary to demonstrate this for lithium‐metal anodes. In this work, advanced high‐resolution 3D chemical analysis is conducted with time‐of‐flight secondary‐ion mass spectrometry (TOF‐SIMS) to establish spatial correlations between the transition metals and electrolyte decomposition products found on cycled lithium‐metal anodes. Insights into the localization of various chemistries linked to crucial processes that define LMB performance, such as lithium deposition, SEI growth, and transition‐metal deposition are deduced from a precise elemental and spatial analysis of the SEI. Heterogenous transition‐metal deposition is found to perpetuate both heterogeneous SEI growth and lithium deposition on lithium‐metal anodes. These correlations are confirmed across various lithium‐metal anodes that are cycled with different cobalt‐free cathodes and electrolytes. An advanced electrolyte that is stable to higher voltages is shown to minimize transition‐metal crossover and its effects on lithium‐metal anodes. Overall, these results highlight the importance of maintaining uniform SEI coverage on lithium‐metal anodes, which is disrupted by transition‐metal crossover during operation at high voltages.

Chemistry↗

NASA NTRSDate not supplied

LES-Based Analysis of Film Cooling in a High-Pressure Turbine Rotor Including Effects of Purge and Tip Flow

Secondary Flows↗

DOE OSTI2026

Machine Learning for Multipactor Susceptibility Prediction in Planar RF Gaps

Multipactor discharge is a nonlinear electron avalanche that limits the performance of high-power radio-frequency (RF) and vacuum electronic devices. Predicting multipactor susceptibility traditionally relies on Monte Carlo or particle-in-cell (PIC) simulations, which become computationally expensive for large parametric studies. In this work, we present a supervised machine-learning (ML) framework for prediction of multipactor susceptibility in a two-surface planar geometry. The models are trained using high-fidelity PIC simulation generated susceptibility data and learn the relationship between operational parameters, geometry, and material-dependent secondary electron emission properties. The proposed approach enables rapid reconstruction of susceptibility charts while preserving the physical structure of multipactor growth regions.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Resonant Metasurfaces with Van Der Waals Hyperbolic Nanoantennas and Extreme Light Confinement

This work reports on a metasurface based on optical nanoantennas made of van der Waals material hexagonal boron nitride. The optical nanoantenna made of hyperbolic material was shown to support strong localized resonant modes stemming from the propagating high-k waves in the hyperbolic material. An analytical approach was used to determine the mode profile and type of cuboid nanoantenna resonances. An electric quadrupolar mode was demonstrated to be associated with a resonant magnetic response of the nanoantenna, which resembles the induction of resonant magnetic modes in high-refractive-index nanoantennas. The analytical model accurately predicts the modes of cuboid nanoantennas due to the strong boundary reflections of the high-k waves, a capability that does not extend to plasmonic or high-refractive-index nanoantennas, where the imperfect reflection and leakage of the mode from the cavity complicate the analysis. In the reported metasurface, excitations of the multipolar resonant modes are accompanied by directional scattering and a decrease in the metasurface reflectance to zero, which is manifested as the resonant Kerker effect. Van der Waals nanoantennas are envisioned to support localized resonances and can become an important functional element of metasurfaces and transdimensional photonic components. By designing efficient subwavelength scatterers with high-quality-factor resonances, this work demonstrates that this type of nanoantenna made of naturally occurring hyperbolic material is a viable substitute for plasmonic and all-dielectric nanoantennas in developing ultra-compact photonic components.

Chemistry↗

DOE OSTI2026

Cation-Diffusive Carbon Interlayers Stabilize Na Metal and Double the Current in Na-S Redox-Flow Batteries for Grid-Scale Energy Storage

The sodium-sulfur nonaqueous redox-flow batteries (Na-S NARFBs) using earth-abundant elements are highly attractive due to the low material cost and improved energy density for grid-scale energy storage. However, the low current performance, poor Na0/Na+ redox kinetics, and Na dendrite growth pose severe challenges. We introduce cation-diffusive layers (CDLs): thin and Na+ affinitive interlayers at the Na anode that direct Na+ transport and stabilize Na deposition. Benchmarking three archetypal materials—carbon paper (CP), glass microfiber paper (GF), and foam—across Na-Na and Na-Cu, and Na-S cells identifies CP as the optimum. CP reduces symmetric cell overpotential by more than 70%, achieves 98% Na plating-stripping efficiency, and doubles the Na-S cell current density from 0.5 to 1.0 mA cm−2 without sacrificing capacity or efficiency. Ex situ electrochemical and SEM/XPS analysis, combined with molecular dynamics (MD) studies, reveal that electron-rich carbon fibers disperse supporting salt aggregates, enrich near-surface Na+ density, and create ion transport pathways for fast Na0/Na+ exchange while mitigating membrane degradation. Because of the ion-centric mechanism, CDLs can be generalized to other metal-anode designs. Further, this work establishes CDL design rules—cationic affinity and appropriate micro/nanostructure—as a simple, scalable route to high-current, durable metal-anode flow batteries.

Wu, Wenda [ORNL] (ORCID:000900033307687X)↗

DOE OSTI2026

Spinel high-entropy oxides (FeNiCrMnZnX) 3 O 4 (X = Al, mg) as anode materials for high-performance lithium-ion batteries

To address the high cost, cobalt dependency, and resource constraints typical of conventional high-entropy oxide (HEO) anodes, this study reports the successful synthesis of two Co-free, six-component spinel-type HEOs(FeNiCrMnZnAl) 3 O 4 (HEO-Al) and (FeNiCrMnZnMg) 3 O 4 (HEO-Mg), via a sol-gel method. The distinct effects of Al 3+ and Mg 2+ incorporation on the electrochemical performance and lithium storage kinetics were systematically investigated. XRD, Raman, and TEM characterizations confirm that both materials possess a pure spinel phase, uniform particle size, and homogeneous elemental distribution. Notably, electrochemical evaluations reveal that HEO-Al delivers a superior reversible capacity of 480.7 mAh g −1 after 100 cycles at 0.1 A g −1 , and maintains 354.5 mAh g −1 after 1000 long-term cycles at 1 A g −1 , significantly outperforming HEO-Mg. Kinetic analysis indicates that HEO-Al exhibits lower charge transfer resistance, a higher Li + diffusion coefficient, and a pseudocapacitive contribution of up to 82%. Furthermore, these findings demonstrate that Al substitution effectively optimizes the structural stability and lithium storage kinetics of Co-free HEOs, providing a viable strategy for designing low-cost, highly stable HEO anode systems.

Anode materials↗

DOE OSTI2026

Evaluation of high-temperature properties of binary, eutectic Ni-X alloys (X = B, Ca, Ce, La, Y, Zr)

We investigated the high-temperature behavior of binary eutectic nickel-base alloys with low-solubility elements, towards developing alloys with new combinations of structural and functional (e.g. thermal or corrosion) properties. Based on thermodynamic screening, we cast binary near-eutectic Ni-B, Ni-Ca, Ni-Ce, Ni-La, Ni-Y, and Ni-Zr alloys and characterized their microstructure and high-temperature properties. All alloys (except Ni-Ca) contain a fine eutectic microstructure, containing submicron alternating lamellae of γ-Ni and continuous intermetallic phases. The eutectic microstructure did not significantly coarsen at 700°C (except for Ni-B), but measurably coarsened in all alloys at 900°C, most quickly in Ni-B, followed by Ni-Ce, Ni-La, Ni-Y, and Ni-Zr. Coarsening in these alloys occurs by fault migration; diffusion within the intermetallic, and its interfacial energy, likely influence coarsening kinetics. At room-temperature, all alloys possess negligible ductility due to the continuous, brittle intermetallic. At 700 and 900°C, the Ni-Zr alloy showed the best combination of strength and ductility. Oxidation experiments at 900°C on Ni-Zr and Ni-Ce revealed heavy internal oxidation of the continuous intermetallic phase. Overall, the Ni-Zr alloy showed good high-temperature coarsening resistance and strength, though ductility and oxidation resistance require improvement for use in structural applications. Possible routes for improving alloy properties via additive manufacturing are discussed.

Ekaputra, Clement [ORNL] (ORCID:0000000307166479)↗

NASA NTRS2026

The Canadian Galactic Emission Mapper: A New 8-10 GHz Telescope To Map the Polarization of the Northern Sky

The Canadian Galactic Emission Mapper (CGEM) radio telescope is mapping polarized Galactic foregrounds from 8-10 GHz to aid in the search for B-modes in the Cosmic Microwave Background (CMB). Here we describe the design and early on-sky characterization of the CGEM optics. CGEM employs an on-axis, hat feed optical design that is a body of revolution (BOR) and exhibits excellent intrinsic polarization purity. We describe how we have optimized the optics with a novel framework that can directly minimize simulated intensity to polarized intensity (T → P) leakage in angular power spectrum space. The optimized optics exhibit simulated T → P leakage that is orders of magnitude below anticipated CℓBB for r = 10 −3 when scaled from 8-10 GHz to the CMB observing window near 95 GHz. We go on to describe the mechanical design of the optics, including a novel secondary mirror support made from Astroquartz composite that has low loss, low dielectric constant, and preserves the BOR symmetry of the optics. We then showcase the early on-sky performance of the optics with observations of the Sun and satellites, which probe the beam to 40dB down from the peak. A beam model for the deployed CGEM optics based on electromagnetic simulations is in excellent agreement with the on-sky data.

cosmic microwave background↗

DOE OSTI2026

Capturing Secondary Kinetic Instabilities in Three‐Dimensional Dayside Reconnection Using an Improved Gradient‐Based Closure

Magnetic reconnection is a highly dynamic process that excites a wide variety of kinetic waves and instabilities. Transverse current sheet instabilities such as the lower-hybrid drift and secondary drift-kink instabilities in particular have been shown by kinetic simulations to modify the reconnection and introduce significant turbulence and mixing to the reconnection layer. Past studies using the ten-moment fluid model to capture important kinetic physics such as the electron inertia and full representation of the pressure tensor proved advantageous to a two-fluid representation of reconnection, but the model struggled when using a local relaxation closure for the heat flux to replicate the current sheet instabilities and subsequent mixing seen in kinetic simulations. This work uses the Gkeyll software framework to perform simulations of asymmetric reconnection based on the 16 October 2015 MMS crossing of a diffusion region, the Burch event. An improved gradient-based heat flux closure is implemented, showing significant improvement in secondary kinetic instabilities that grow in the current sheet. These instabilities generate turbulence which leads to growth of secondary magnetic islands and flux ropes.

Bradshaw, K. [Princeton University, NJ (United Sta↗

DOE OSTI2026

String-Breaking Dynamics in Quantum Adiabatic and Diabatic Processes

Confinement prohibits isolation of color charges, e.g., quarks, in nature via a process called string breaking : the separation of two charges results in an increase in the energy of a color flux, visualized as a string, connecting those charges. Eventually, creating additional charges is energetically favored, hence breaking the string. Such a phenomenon can be probed in simpler models, including quantum spin chains, enabling enhanced understanding of string-breaking dynamics. A challenging task is to understand how string breaking occurs as time elapses, in an out-of-equilibrium setting. This work establishes the phenomenology of dynamical string breaking induced by a gradual increase of string tension over time. It, thus, goes beyond instantaneous quench processes and enables tracking the real-time evolution of strings in a more controlled setting. We focus on domain-wall confinement in a family of quantum Ising chains. Our results indicate that, for sufficiently short strings and slow evolution, string breaking can be described by the transition dynamics of a two-state quantum system akin to a Landau-Zener process. For longer strings, a more intricate spatiotemporal pattern emerges: the string breaks by forming a superposition of bubbles (domains of flipped spins of varying sizes), which involve highly excited states. We finally demonstrate that string breaking driven only by quantum fluctuations can be realized in the presence of sufficiently long-ranged interactions. This work holds immediate relevance for studying string breaking in quantum-simulation experiments.

Ising model↗

DOE OSTI2025

Bi‐Continuous W‐Rich Refractory High Entropy Alloy‐Cu Composite: Toward Material Innovation of Nuclear Reactor Coolant System

Abstract Refractory high‐entropy alloys (RHEAs) are considered promising candidate materials for next‐generation nuclear reactors due to their superior mechanical strength, irradiation resistance, and thermal stability at high temperatures. However, the significant positive heat of mixing between refractory alloying elements and Cu, commonly used in cooling systems, poses challenges in forming composite structures. This study addresses the issue using a liquid metal dealloying (LMD) process. A precursor alloy (WTaVTi) with a directional dendrite‐interdendrite structure is fabricated and reacted with molten Cu at 1200 °C for 96 h. This approach produced a RHEA‐Cu composite with a stable interface between RHEA (W 31.5 Ta 30.9 V 21.4 Ti 14.3 ) and Cu, featuring a spontaneously formed W‐rich interlayer that enhances interfacial bonding. The composite showed excellent irradiation resistance, with 30% less swelling under α‐ion irradiation than pure W. It also exhibited low thermal conductivity at room temperature, but reached ≈120 W m −1 ·K −1 at ≈650 °C, surpassing pure W. This temperature‐dependent rise in κ, with a positive gradient of +0.075 W m −1 ·K − 2 , is attributed to decreasing diffuse mismatch at elevated temperatures. The large‐scale reaction and stable microstructure achieved through LMD process highlight its industrial potential. This work offers a strategy for developing high‐performance materials by combining RHEA's radiation resistance with Cu's thermal conductivity for extreme environments.

Chemistry↗

NASA NTRSDate not supplied

Preliminary Assessment of the Impact on the V-Band Oxygen Channels From Satellite Communication Uplinks

We calculate the percentage of time that an ATMS-like instrument [1] will be illuminated by the uplink beam of one of the proposed V-band communication system and estimate the damage resulting from such exposure. Using a combination of openly available information and educated guesses about the location and characteristics of the up/down link terminals, we constructed the ground segment of a hypothetical high-speed communication network. The space segment of the network was constructed from the orbital data of the existing Starlink constellation [2] of 6223 communication satellites (comsats) which is used as strawman to represent any other possible constellation of communication satellites. It is shown that without a very delicate balance of frequency allocations (science vs telecommunications), coupled with extremely steep and deep bandpass-defining filters, and strict adherence to the agreed limits (i.e. no out-of-band transmissions) the deployment of the telecommunication network leads to almost-complete loss of some important geophysical data. For the analysis we use the spectral characteristics of the ATMS instrument with the ephemeris for the NOAA-21 satellite [3]. The analysis is conducted for the USA and the simulation covers 8 consecutive days in July 2024. Effective and accurate microwave remote sensing of the atmosphere depends on the availability of interference-free spectrum windows at frequencies which are prescribed by physical processes [e.g. 4]. The family of resonant lines of the oxygen molecule near 60 GHz provides a unique opportunity to sample the vertical distribution of temperature and density from space, and it has been exploited for weather and climate studies from polar-orbiting satellites since 1978 (MSU on TIROS-N [5]). It remains a staple in the payloads operated by Russia, China, USA, Japan, France, India, UK, Ukraine [6] which are built around a common blueprint: a few wide-band (hundreds of MHz) channels around 50 GHz to sample the atmosphere and the surface while several more channels with high spectral resolution (few MHz) sample the individual resonant lines. Accurate retrieval of the environmental parameters depends upon the data provided by both sets of channels, and the their location in frequency space is not arbitrary and cannot be altered at will [7, 8]. The introduction of 5G technology in 2019 has driven telecommunication companies to request more bandwidth to be dedicated to their devices. This additional bandwidth is only available in spectral regions traditionally reserved for environmental and astrophysical research, such as the V-band between 50 and 60 GHz for up/downlink between satellites in low-earth orbits and terminals connected to fiberoptics network for distribution to high-speed local internet services. The power broadcast by the uplink communication leg is many orders of magnitude greater than the natural thermal signal emitted from the Earth scene. If the ground antenna were to perfectly align with the passive instrument’s antenna, the spaceborne receiver would suffer permanent, irreparable damage. While a direct boresight-to-boresight conjunction is extremely unlikely (even with a large constellation of satellites the fraction of the celestial sphere occupied by the satellites remains minuscule) the finite size of the ground station’s antenna beam in the sky suggests that the ATMS will be in the near background (as seen from the ground station) of one of the communication satellites and will be illuminated by either the main lobe or the near sidelobes of the uplink antenna more often than it is desirable. For our analysis we first calculate the position of the ATMS with respect to each of the ground stations at a resolution of 0.2 sec, then calculate the position of each of the comsats which are at least 25 deg above the station’s local horizon; finally we calculate the angle between the line-of-sight of the ATMS and the line-of-sight of the comsat. We assume that the gain pattern of the ground station is circularly symmetric; the angle-off-station-boresight then provides an attenuation of the uplink power which we use to assess the likely effect upon the passive instrument’s operations. We assume that each ground station can communicate with all the comsats in its field of view; this implies that, on average, a ground station can engage with 46 comsats simultaneously. The analysis is repeated for the case when the uplink broadcast within the ATMS passive channels (in-band scenario) and for the case when the uplink is limited to frequencies adjacent to the ATMS channels (out-of-band scenario). The antenna of the ground station is modelled as having a HPBW (Half-Power Beam Width) of 0.16 deg and EIRP (Equivalent Isotropic Radiated Power) of 70 dBW. We account for the geometric dissipation of the signal caused by the satellite orbital altitude, the attenuation induced by atmospheric gasses at 51 GHz and the mismatch between the circular polarization of the ground-based transmitting antenna and the linear polarization of the satellite-borne receiving antenna. The damages on ATMS are estimated from bench-level measurement conducted at the ATMS’ manufacturer facilities [unpublished].

passive microwave↗

DOE OSTI2025

Microstructure, mechanical, and thermal properties of compositionally complex (Hf,Zr,Nb,Ti)B 2 ‒LaB 6 ceramics

Novel compositionally complex borides, (Hf,Zr,Nb,Ti)B 2 and (Hf,Zr,Nb,Ti)B 2 ‒LaB 6 , were fabricated using spark plasma sintering process. (Hf,Zr,Nb,Ti)B 2 ‒LaB 6 exhibits a dual-phase microstructure, in which (Hf,Zr,Nb,Ti)B 2 is a primary phase with the hexagonal structure and LaB 6 is a secondary phase with a cubic structure. The mechanical properties of both (Hf,Zr,Nb,Ti)B 2 and (Hf,Zr,Nb,Ti)B 2 ‒LaB 6 are comparable, with a combination of high hardness and moderate fracture toughness. Thermal diffusivity and conductivity of (Hf,Zr,Nb,Ti)B 2 are much lower than the individual transition metal borides but are significantly increased by the addition of LaB 6 . Herein, it is implied that the thermal properties of boride ceramics can be controlled through the appropriate design of principal metal element compositions.

borides↗

NASA NTRS2026

Isotopic Signature of Organic Molecules from Beyond the Solar System: An Enriched Methane D/H Ratio in the Interstellar Object 3I/ATLAS

Interstellar objects are interlopers from other planetary systems, and their volatile compositions provide a glimpse into planet formation around their host star. We present near-infrared spectra of the coma of interstellar object 3I/ATLAS measured with the James Webb Space Telescope. Our results demonstrate an unexpectedly high D/H = (3.31 ± 0.34)% for methane and represent an exceedingly rare detection of deuterated organic molecules in an interstellar object. This D/H ratio is a factor of 14 ± 2 higher than that measured in comet 67P/Churyumov-Gerasimenko by the Rosetta spacecraft, the only other comet for which CH3D has been detected, yet the ratio of deuteration in methane compared with water is consistent for both comets within 1.2σ. The D/H ratio in methane is observationally unconstrained in extrasolar sources to date, but the enriched ratio in 3I/ATLAS is most similar to those measured in other organic molecules toward primitive environments. The high D/H ratios of water and methane in 3I/ATLAS are a natural consequence of formation in a high D/H elemental ratio environment as a result of locally cold conditions in the protoplanetary disk and prior interstellar cloud. Thus, 3I/ATLAS formed in an environment very different from that in which our Sun and planets originated.

Interstellar Object 3I/ATLAS↗

DOE OSTI2026

Impacts of Multidimensional Progenitor Perturbations on Core-collapse Supernova Explosions

Numerical studies of core-collapse supernovae have demonstrated the importance of nonradial motions in precollapse progenitors on the explosion outcome. We use the Chimera neutrino radiation hydrodynamics code running seven two-dimensional simulations of 15 M⊙ progenitors with different progenitor structures introduced by different one- and two-dimensional precollapse stellar evolution environments to examine the impacts of stellar structure and nonspherical motion in the precollapse progenitor on the development of explosions. We compare the explosion evolution of these models in terms of shock dynamics, diagnostic energy, neutrino heating, accretion, explosion geometry, nuclear abundances, and turbulent convection. We also analyze how stochastic variation impacts our simulations. Contrary to results reported in prior studies examining the impacts of multidimensional progenitors, we observe similar shock revival times and explosion development in our simulations despite differences in initial compositions and structures. We find no discernible impact from the accretion of nonradial perturbations from a multi-D progenitor onto the stalled shock in the revival and strength of explosion, as fully developed neutrino-driven convection behind the stalled shock is similar for all our models. For models with physically sourced noise in the iron core, a strong oscillation of the shock occurs after bounce and deflects infall laterally, and accelerates the saturation of the lateral turbulent kinetic energy. An examination of model stochasticity shows that any prior expected impacts on explosive outcome due to convection-related perturbations lie below the detectable threshold of numerical variation.

Chen, Chien-Hui [North Carolina State University]↗