Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity profiles”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Method for Spatiotemporal Solar Power Profile Estimation for a Proposed U.S.–Caribbean–South America Super Grid under Hurricanes

Solar photovoltaic (PV) generation technology stands out as a scalable and cost-effective solution to enable the transition toward decarbonization. However, PV solar output, beyond the daily solar irradiance variability and unavailability during nights, is very sensitive to weather events like hurricanes. Hurricanes nucleate massive amounts of clouds around their centers, shading hundreds of kilometers in their path, reducing PV power output. This research proposes a spatiotemporal method, implemented in MATLAB R2023b coding, to estimate the shading effect of hurricanes over a wide distribution of PV solar plants connected to a high-voltage power infrastructure called the U.S.–Caribbean–South America super grid. The complete interconnection of the U.S., the Caribbean, and South America results in the lowest power valley levels, i.e., an overall percentual reduction in PV power output caused by hurricane shading. The simulations assess the impact of hurricanes in 10 synthetic trajectories spanning from Texas to Florida. The Caribbean would also experience lower power valleys with expanded interconnectivity schemes. The U.S.–Caribbean–South America super grid reduces Caribbean variability from 37.8% to 8.9% in the case study. The proposed spatiotemporal method for PV power profile estimation is a valuable tool for future solar power generation expansion, transmission planning, and system design considering the impact of hurricanes.

14 SOLAR ENERGY↗
DOE OSTI2020

Validation and Sensitivity Analysis of a 1-D Lake Model Across Global Lakes

Lakes have important influence on weather and climate from local to global scales. However, their prediction using numerical models is notoriously difficult because lakes are highly heterogeneous across the globe, but observations are sparse. In this study, we assessed the performance of a 1-D lake model in simulating the thermal structures of 58 lakes with diverse morphometric and geographic characteristics by following the phase 2a local lake protocol of the Inter-sectoral Impact Model Intercomparison Project (ISIMIP2a). After calibration, the root-mean-square errors (RMSE) were below 2 °C for 70% and 75% of the lakes for epilimnion and full-profile temperature simulations, with an average of 1.71 °C and 1.43 °C, respectively. The model performance mainly depended on lake shape rather than location, supporting the possibility of grouping model parameters by lake shape for global applications. Furthermore, through machine-learning based parameter sensitivity tests, we identified turbulent heat fluxes, wind-driven mixing and water transparency as the major processes controlling lake thermal and mixing regimes. Snow density was also important for modeling the ice phenology of high-latitude lakes. The relative influence of the key processes and the corresponding parameters mainly depended on lake latitude and depth. Turbulent heat fluxes showed a decreasing importance in affecting epilimnion temperature with increasing latitude. Wind-driven mixing was less influential to lake stratification for deeper lakes while the impact of light extinction, on the contrary, showed a positive correlation with depth. Our findings may guide improvements in 1-D lake model parameterizations to achieve higher fidelity in simulating global lake thermal dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Design and commissioning of resistive foil bolometer diagnostics on the ST40 tokamak

Following successful campaigns on the compact high field spherical tokamak ST40, resistive gold foil bolometers have now been installed to measure the radiated power profile. Positioned on the midplane, two bolometer cameras offer perspectives of the horizontal and vertical planes, while the third camera, situated above the midplane, provides a vertical view of the top divertor plates, including the X-point and the strike points. Further, these cameras use a commercial off-the-shelf data acquisition platform with modules specifically designed for bolometers in a Wheatstone bridge configuration, providing both phase sensitive detection and in situ calibration. This paper presents an overview of the design and functionality of all three cameras and the commissioning of the horizontal plane camera. The commissioning of the horizontal camera shows data from a range of plasma pulses, presenting a positive correlation between the relative profiles measured from the bolometer camera and an absolute extreme ultraviolet diode camera viewing along similar lines of sight.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Remote detection of radioactive material using a short-pulse CO 2 laser

Detection of radioactive material at distances greater than the radiated particle range is an important goal with applications in areas such as national defense and disaster response. Here, we demonstrate avalanche-breakdown-based remote detection of a 3.6 mCi α-particle source at a standoff distance of 10 m, using 70 ps, long-wave infrared (λ = 9.2 µm) CO 2 laser pulses. This is ∼10 times longer than our previous results using a mid-IR laser. The primary detection method is direct backscatter from microplasmas generated in the laser focal volume. The backscatter signal is amplified as it propagates back through the CO 2 laser chain, enhancing sensitivity by >100 times. Here we also characterize breakdown plasmas with fluorescence imaging, and present a simple model to estimate backscattered signals as a function of the seed density profile in the laser focal volume. All of this is achieved with a relatively long-drive laser focal geometry (f/200) that is readily scalable to >100 m.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

How wind speed shear and directional veer affect the power production of a megawatt-scale operational wind turbine

Most megawatt-scale wind turbines align themselves into the wind as defined by the wind speed at or near the center of the rotor (hub height). However, both wind speed and wind direction can change with height across the area swept by the turbine blades. A turbine aligned to hub-height winds might experience suboptimal or superoptimal power production, depending on the changes in the vertical profile of wind, also known as shear. Using observed winds and power production over 6 months at a site in the high plains of North America, we quantify the sensitivity of a wind turbine's power production to wind speed shear and directional veer as well as atmospheric stability. We measure shear using metrics such as α (the log-law wind shear exponent), β bulk (a measure of bulk rotor-disk-layer veer), β total (a measure of total rotor-disk-layer veer), and rotor-equivalent wind speed (REWS; a measure of actual momentum encountered by the turbine by accounting for shear). We also consider the REWS with the inclusion of directional veer, REWS θ , although statistically significant differences in power production do not occur between REWS and REWS θ at our site. When REWS differs from the hub-height wind speed (as measured by either the lidar or a transfer function-corrected nacelle anemometer), the turbine power generation also differs from the mean power curve in a statistically significant way. This change in power can be more than 70 kW or up to 5 % of the rated power for a single 1.5 MW utility-scale turbine. Over a theoretical 100-turbine wind farm, these changes could lead to instantaneous power prediction gains or losses equivalent to the addition or loss of multiple utility-scale turbines. At this site, REWS is the most useful metric for segregating the turbine's power curve into high and low cases of power production when compared to the other shear or stability metrics. Therefore, REWS enables improved forecasts of power production.

17 WIND ENERGY↗
DOE OSTI2020

Changing the rotational direction of a wind turbine under veering inflow: a parameter study

All current-day wind-turbine blades rotate in clockwise direction as seen from an upstream perspective. The choice of the rotational direction impacts the wake if the wind profile changes direction with height. Here, we investigate the respective wakes for veering and backing winds in both hemispheres by means of large-eddy simulations. We quantify the sensitivity of the wake to the strength of the wind veer, the wind speed, and the rotational frequency of the rotor in the Northern Hemisphere. A veering wind in combination with counterclockwise-rotating blades results in a larger streamwise velocity output, a larger spanwise wake width, and a larger wake deflection angle at the same downwind distance in comparison to a clockwise-rotating turbine in the Northern Hemisphere. In the Southern Hemisphere, the same wake characteristics occur if the turbine rotates counterclockwise. These downwind differences in the wake result from the amplification or weakening or reversion of the spanwise wind component due to the effect of the superimposed vortex of the rotor rotation on the inflow's shear. An increase in the directional shear or the rotational frequency of the rotor under veering wind conditions increases the difference in the spanwise wake width and the wake deflection angle between clockwise- and counterclockwise-rotating actuators, whereas the wind speed lacks a significant impact.

17 WIND ENERGY↗
DOE OSTI2024

Ultra-Long Distance BOTDA Sensor System Employing Hybrid Amplification and Advanced Noise Reduction Techniques

Brillouin Optical Time Domain Analysis (BOTDA) sensor systems play a pivotal role in distributed sensing, which enables precise measurements of strain and temperature across extensive fiber lengths. This research offers a comprehensive strategy to extend the sensing range of BOTDA systems beyond tens of kilometers while maintaining high spatial resolutions. Such enhanced sensing is realized through the integration of distributed Raman amplification, inline amplification using erbium-doped fiber amplifiers (EDFA), and advanced noise reduction techniques. Optimizing Raman pump and probe wave profiles ensures the robustness of Brillouin scattering signals, effectively countering attenuation-induced losses. Additionally, enhancing weak Brillouin signals by strategically placing EDFAs at optimal fiber locations keeps sensor sensitivity high. Leveraging inherent redundancy in measured data as a function of frequency and fiber distance, the non-local means (NLM) filter removes noise while preserving essential physical information. In summary, this research has shown a holistic exploration of extending BOTDA’s distance sensing capabilities up to 150 km with spatial resolutions of 8 meters.

Bhatta, Hari↗
DOE OSTI2022

Microstructure-Sensitive modeling of surface roughness and notch effects on extreme value fatigue response

Fatigue life knockdown due to surface roughness induces local stress concentrations. Physical fatigue experiments are crucial to qualify materials for use but are limited in their capability to examine the mechanisms underlying the fatigue crack formation driving force due to these stress concentrations. Microstructure-sensitive computational models can aid in this regard and are employed in this work to examine the effects of surface roughness in Al 7075-T6. Digital microstructure models are subjected to crystal plasticity finite element method (CPFEM) simulations to examine the combined effects of microstructure and grain scale asperities on the driving force for fatigue crack formation. Following elastic–plastic shakedown, mesoscale volume-averaged fatigue indicator parameters (FIPs) are computed within fatigue damage process zones of grains. An increase in the intensity of realistic surface roughness profiles corresponds accordingly to larger FIPs but it is difficult to precisely interpret the mechanisms and effects of this intensification. We thus parametrically investigate surface asperities that couple with microstructure and lack of constraint on slip at the free surface to describe microstructure-sensitive surface roughness knockdown effects on fatigue resistance. The depth of a single valley (i.e., notch) due to surface roughness is found to be more detrimental to fatigue resistance than the radius (i.e., notch width, not to be confused with notch acuity). Finally, we also find that different crystallographic textures, grain morphologies, and number of sampled grains have a limited effect on the depth of influence of these notches.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An Experimental Kinetics Study of Isopropanol Pyrolysis and Oxidation behind Reflected Shock Waves

Isopropanol has potential as a future bio-derived fuel and is a promising substitute for ethanol in gasoline blends. Even so, little has been done in terms of high-temperature chemical kinetic speciation studies of this molecule. To this end, experiments were conducted in a shock tube using simultaneous CO and H2O laser absorption measurements. Water and CO formation during isopropanol pyrolysis was also examined at temperatures between 1127 and 2162 K at an average pressure of 1.42 atm. Species profiles were collected at temperatures between 1332 and 1728 K and at an average pressure of 1.26 atm for equivalence ratios of 0.5, 1.0, and 2.0 in highly diluted mixtures of 20% helium and 79.5% argon. Species profiles were also compared to four modern C3 alcohol mechanisms, including the impact of recent rate constant measurements. The Li et al. (2019) and Saggese et al. (2021) models both best predict CO and water production under pyrolysis conditions, while the AramcoMech 3.0 and Capriolo and Konnov models better predict the oxidation experimental profiles. Additionally, previous studies have collected ignition delay time (τign) data for isopropanol but are limited to low pressures in highly dilute mixtures. Therefore, real fuel–air experiments were conducted in a heated shock tube with isopropanol for stoichiometric and lean conditions at 10 and 25 atm between 942 and 1428 K. Comparisons to previous experimental results highlight the need for real fuel–air experiments and proper interpretation of shock-tube data. The AramcoMech 3.0 model over predicts τign values, while the Li et al. model severely under predicts τign. The models by Capriolo and Konnov and Saggese et al. show good agreement with experimental τign values. A sensitivity analysis using these two models highlights the underlying chemistry for isopropanol combustion at 25 atm. Additionally, modifying the Li et al. model with a recently measured reaction rate shows improvement in the model’s ability to predict CO and water profiles during dilute oxidation. Finally, a regression analysis was performed to quantify τign results from this study.

09 BIOMASS FUELS↗
DOE OSTI2024

A magnetic‐directed micro‐particle with near‐ IR light triggered guest‐release property

Abstract In this study, guest carriers, composed of magnetic‐directed hydrogel particles with near‐infrared (NIR)‐light‐triggered guest‐release properties, are prepared via the one‐pot method. The gel particles contain photostable NIR‐responsive polypyrrole nanoparticles with temperature‐sensitive poly(N‐isopropyl acrylamide), the combination of which can induce the accelerated release of encapsulated drugs by 337% upon 5 minutes of NIR light irradiation based on the release profile of NRMG‐1 at 10–15 min. Additionally, with magnetic iron oxide particles crosslinking the polyvinyl alcohol frames, these gel particles can form a physical barrier that decreases the release rate by 20% under an external magnetic field over 24 h. The magnetic iron oxide particles also enable the gel particles to reach their target destination without extra drug losses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Low-mode nonuniformity in direct-drive ICF implosions due to laser smoothing techniques employed on OMEGA

For successful laser-direct-drive inertial confinement fusion implosions, the laser irradiation must be highly uniform over the target surface. On OMEGA, multiple laser beams are used to illuminate targets quasi-uniformly. High-mode-number nonuniformities due to laser speckle on each individual beam are reduced by splitting each beam into two orthogonal polarizations (i.e., polarization smoothing, or PS) and a range of wavelengths (i.e., smoothing by spectral dispersion) that are dispersed at the target plane. However, cross-beam energy transfer (CBET) is sensitive to both the polarizations and wavelengths of the interacting beams, so the interplay between CBET and the laser-smoothing schemes results in unique intensity variation across each beam profile, which is a systematic source of low-mode drive nonuniformity on OMEGA. Here, we model these effects and find that the predicted ℓ = 1 mode in the laser-absorption distribution is consistent with the systematic core-flow direction that has been determined from the OMEGA implosion database. We also observe good agreement with the measured core-flow directions for two specific sets of implosions (one with PS, the other without PS) when we also account for the measured beam mispointing and the beam power imbalance.

Crossed beam scattering↗
DOE OSTI2024

Astrophysical constraints from synchrotron emission on very massive decaying dark matter

If the cosmological dark matter (DM) couples to Standard Model (SM) fields, it can decay promptly to SM states in a highly energetic hard process, which subsequently showers and hadronizes to give stable particles including e ± , γ , p ± , and ν ν ¯ at lower energy. If the DM particle is very heavy, the high-energy e ± , due to the Klein-Nishina cross section suppression, preferentially lose energy via synchrotron emission which, in turn, can be of unusually high energies. Here, we present previously unexplored bounds on heavy decaying DM up to the Planck scale, by studying the synchrotron emission from the e ± produced in the ambient Galactic magnetic field. In particular, we explore the sensitivity of the resulting constraints on the DM decay width to (i) different SM decay channels, to (ii) the Galactic magnetic field configurations, and to (iii) various different DM density profiles proposed in the literature. We find that constraints from the synchrotron component complement and improve on constraints from very high-energy cosmic-ray and gamma-ray observatories targeting the prompt emission when the DM is sufficiently massive, most significantly for masses in excess of 10 12 GeV . Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2021

Laser-induced fluorescence thermometry of supercritical CO 2 flows inside a micro-channel

This work demonstrates a thermometric technique using laser-induced fluorescence (LIF) in supercritical carbon dioxide flows in a micro-channel. Rhodamine 6G was used as a temperature-sensitive fluorescent dye. The flow conditions were at a pressure of 7.9 MPa and temperature in the range of 23°–90°C. 2D spatial distributions and time-resolved temperature profiles were obtained at this high pressure. Measured LIF signals showed close relations to the temperatures obtained from resistance temperature detectors.

Park, Suhyeon↗
DOE OSTI2022

Geothermal Fault Zone and Fluid Imaging through Joint Airborne ZTEM and Ground MT Data Inversion Analysis

This project has aimed to achieve detailed electrical resistivity resolution at geothermal reservoir scales by combining airborne natural electromagnetic (EM) field surveying (ZTEM) with ground magnetotelluric (MT) measurements to approximate an airborne MT geophysical method. MT alone is relatively expensive and may have permitting challenges in sensitive areas. Airborne ZTEM field data contains only the magnetic field, requires a background assumption, and has been limited to relatively high frequencies, thus suffering uniqueness problems. Based on proto-type 2D simulations, ZTEM ambiguities may be reduced through formal incorporation with possibly sparse ground MT soundings, which we pursued in full 3D for this project. The methodology was tested at the high-temperature Roosevelt Hot Springs geothermal system, Utah, which was considered advantageous given the near total exposure of crystalline reservoir rocks across the project area. ZTEM and ground MT survey data were acquired in 2017, subcontracted to outside parties with which we have worked in the past. These included 80 remote-referenced tensor MT soundings over the Mineral Mountains and adjacent Roosevelt Hot Spring producing geothermal system. These MT stations abut later coverage of a similar number of MT stations taken for the Utah FORGE project providing excellent total data aperture to re-solve structure beneath both project areas better than either set alone. The airborne ZTEM survey covered 704 line kilometers in E-W flight lines with a 250 m line spacing. Although this survey was timed during a maintenance-related shutdown of power production at the Roosevelt Hot Springs, other noise sources difficult to identify but including two high-voltage state-scale transmission lines compromised the ZTEM survey badly leading to unusable responses. Thus, with DOE management concurrence, the project proceeded to emphasize inversion and interpretation of the joint SubTER-FORGE MT data sets with regard to the Roosevelt Hot Springs reservoir recharge and to deep heat sources for both it and the Utah FORGE EGS project area. We also investigated the joint ZTEM-MT sampling concept with data sets from the Eleven Mile Canyon prospect area donated by the U.S. Navy (A. Sabin, PoC). Inversion of the SubTER-FORGE MT data using the HexMT 3D finite element algorithm reveals a large, low-resistivity anomaly extending sub-vertically through the depth range of the crust beneath the western Mineral Mountains. The steep conductive zone connects in the lower crust to a more tabular conductor characteristic of much of the Great Basin that generally is ascribed to current mafic magmatic underplating, hybridization and fluid release. The location of the resolved anomaly relative to the recent (0.5-0.8 Ma) eruptive centers of the Mineral Mountains implicates it as remnants of the magma body which fed these centers. This structure appears to be currently feeding heat and fluids upward into the Roosevelt Hot Springs hydrothermal system, as well as heat laterally to the FORGE project area. Separate and joint inversion models were carried out for the donated Eleven Mile Canyon MT-ZTEM data set to demonstrate concept. ZTEM only inversion showed two main alteration zones in the western portion of the project area known from geological mapping. Joint inversion including an E-W profile of MT soundings sharpened these features considerably. It also resolved in much greater detail the graben related normal faulting structure of the central project area which lies at depths exceeding the sensitivity of ZTEM alone. The sparse number of MT da-ta relative to the ZTEM required upweighting the former by a factor of several, but an exact procedure awaits future research. Our final impression is that sparse MT data can improve resolution of the subsurface over that of ZTEM alone. However, well sampled MT data are to be preferred and offer the simplicity of interpreting just one data type, and possess the superior resolution capability coming with the electric field everywhere, and from their high bandwidth.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Thermal modeling of three lakes within the continuous permafrost zone in Alaska using the LAKE 2.0 model

Abstract. Lakes in the Arctic are important reservoirs of heat with much lower albedo in summer and greater absorption of solar radiation than surrounding tundra vegetation. In the winter, lakes that do not freeze to their bed have a mean annual bed temperature >0 ∘C in an otherwise frozen landscape. Under climate warming scenarios, we expect Arctic lakes to accelerate thawing of underlying permafrost due to warming water temperatures in the summer and winter. Previous studies of Arctic lakes have focused on ice cover and thickness, the ice decay process, catchment hydrology, lake water balance, and eddy covariance measurements, but little work has been done in the Arctic to model lake heat balance. We applied the LAKE 2.0 model to simulate water temperatures in three Arctic lakes in northern Alaska over several years and tested the sensitivity of the model to several perturbations of input meteorological variables (precipitation, shortwave radiation, and air temperature) and several model parameters (water vertical resolution, sediment vertical resolution, depth of soil column, and temporal resolution). The LAKE 2.0 model is a one-dimensional model that explicitly solves vertical profiles of water state variables on a grid. We used a combination of meteorological data from local and remote weather stations, as well as data derived from remote sensing, to drive the model. We validated modeled water temperatures with data of observed lake water temperatures at several depths over several years for each lake. Our validation of the LAKE 2.0 model is a necessary step toward modeling changes in Arctic lake ice regimes, lake heat balance, and thermal interactions with permafrost. The sensitivity analysis shows us that lake water temperature is not highly sensitive to small changes in air temperature or precipitation, while changes in shortwave radiation and large changes in precipitation produced larger effects. Snow depth and lake ice strongly affect water temperatures during the frozen season, which dominates the annual thermal regime of Arctic lakes. These findings suggest that reductions in lake ice thickness and duration could lead to more heat storage by lakes and enhanced permafrost degradation.

Clark, Jason A. (ORCID:0000000248105261)↗
DOE OSTI2021

Migration Kinetics of Surface Ions in Oxygen-Deficient Perovskite During Topotactic Transitions

Oxygen diffusivity and surface exchange kinetics underpin the ionic, electronic, and catalytic functionalities of complex multivalent oxides. Towards understanding and controlling the kinetics of oxygen transport in emerging technologies, it is highly desirable to reveal the underlying lattice dynamics and ionic activities related to oxygen variation. In this study, the evolution of oxygen content is identified in real-time during the progress of a topotactic phase transition in La 0.7 Sr 0.3 MnO 3-δ epitaxial thin films, both at the surface and throughout the bulk. Using polarized neutron reflectometry, a quantitative depth profile of the oxygen content gradient is achieved, which, alongside atomic-resolution scanning transmission electron microscopy, uniquely reveals the formation of a novel structural phase near the surface. Surface-sensitive X-ray spectroscopies further confirm a significant change of the electronic structure accompanying the transition. The anisotropic features of this novel phase enable a distinct oxygen diffusion pathway in contrast to conventional observation of oxygen motion at moderate temperatures. The results provide insights furthering the design of solid oxygen ion conductors within the framework of topotactic phase transitions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Response of Atmospheric River Width and Intensity to Aquaplanet Warming: A Detection Algorithm‐ and Background Moisture‐Independent Approach

The width of an atmospheric river (AR) is an important parameter when evaluating its impact. Although previous research suggests ARs will widen with global warming, a precise response has been muddled by the large sensitivity of width to a diverse set of AR detection techniques (ARDTs). Here, we propose a methodology that removes the influence of the ARDT by modeling AR‐integrated vapor transport (IVT) profiles as idealized exponential curves with free parameters given by background IVT, intensity above background IVT, and profile width. Notably, our definition for AR profile width does not include any explicit numerical thresholds, relative or absolute, for IVT. We apply our approach to a series of idealized aquaplanet experiments, first with a baseline sea surface temperature (SST) distribution, and then with +2K, +4K, and +6K uniform warming, so as to determine the contributions of each free parameter to AR width. We also apply our approach to high‐resolution atmosphere‐only models forced with SSTs modified to preserve historical variability but following projected warming over 2016–2050. Our results show that contributions to impacts‐relevant AR widening comes primarily from enhancements in background IVT and AR intensity, as opposed to from dynamic width changes.

54 ENVIRONMENTAL SCIENCES↗
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↗