Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic inference”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Probing the lattice structure of dynamically compressed and released single crystal iron through the alpha to epsilon phase transition

In this work, experiments using broadband Laue x-ray diffraction (XRD) were used to examine the lattice structure of dynamically compressed [100]-oriented single crystal iron samples at the Dynamic Compression Sector at the Advanced Photon Source. These experiments used 1 μm thick iron single crystal samples sandwiched between a polyimide ablator and a polycarbonate window. A 100 J, 10 ns duration laser pulse incident on the polyimide ablator was used to shock compress the iron samples to initial stresses greater than 25 GPa, exceeding the ~13 GPa alpha (body-centered-cubic or bcc structure) to epsilon (hexagonal-close-packed or hcp structure) phase transition stress. XRD measurements were performed at various times relative to the shock wave entering the iron sample: early times, <~150 ps while the initial shock waves propagated through the iron; intermediate times, after the iron equilibrated with the ablator and window reaching a plateau stress state (12 or 17 GPa) lasting several nanoseconds; and late times, during uniaxial strain release. The early time measurements show that in <~150 ps, the high-pressure hcp phase is relaxed with a c/a ratio of 1.61, contrary to previous laser shock experiments where a c/a ratio of 1.7 was inferred. In the plateau stress state and partially released states, XRD measurements showed that the hcp structure retained a c/a ratio of 1.61 with no observable changes in the microstructure. Upon stress release at ~1 GPa/ns release rate, the reverse phase transition (hcp to bcc) to the original single crystal orientation (implying a transformation memory effect) was observed to reach completion somewhere between 13 and 11 GPa, indicating little stress hysteresis under rapid uniaxial strain release. A similar memory effect for the reverse hcp to bcc transformation has been previously observed under hydrostatic compression. However, the bcc/hcp orientation relationships differ somewhat between dynamic and static compression experiments, implying that the transformation pathway under uniaxial dynamic strain differs from the Burgers mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multiscale Reactive Model for 1,3,5-Triamino-2,4,6-trinitrobenzene Inferred by Reactive MD Simulations and Unsupervised Learning

When high-energy-density materials are subjected to thermal or mechanical insults at extreme conditions (shock loading), a coupled response between the thermo-mechanical and chemical behaviors is systematically induced. Herein we develop a reaction model for the fast chemistry of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) at the mesoscopic scale, where the chemical behavior is determined by underlying microscopic reactive simulations. The slow carbon cluster formation is not discussed in the present work. All-atom reactive molecular dynamics (MD) simulations are performed with the ReaxFF potential, and a reduced-order chemical kinetics model for TATB is fitted to isothermal and adiabatic simulations of single crystal chemical decomposition. Unsupervised machine learning techniques based on non-negative matrix factorization are applied to MD trajectories to model the decomposition kinetics of TATB in terms of a four-component model. The associated heats of reaction are fit to the temperature evolution from adiabatic decomposition trajectories. Using a chemical species analysis, we show that non-negative matrix factorization captures the main chemical decomposition steps of TATB and provides an accurate estimation of their evolution with temperature. The final analytical formulation, coupled to a diffusion term, is incorporated into a continuum formalism, and simulation results are compared one-to-one against MD simulations of 1D reaction propagation along different crystallographic directions and with different initial temperatures. A good agreement is found for both the temporal and spatial evolution of the temperature field.

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

Direct Observation of the Dynamics of Single-Cell Metabolic Activity during Microbial Diauxic Growth

Understanding how genetic information is realized as the behavior of individual cells is a long-term goal of biology but represents a significant technological challenge. In clonal microbial populations, variation in gene regulation is often interpreted as metabolic heterogeneity. This follows the central dogma of biology, in which information flows from DNA to RNA to protein and ultimately manifests as activity. At present, DNA and RNA can be characterized in single cells, but the abundance and activity of proteins cannot. Inferences about metabolic activity usually therefore rely on the assumption that transcription reflects activity. By tracking the atoms from which they build their biomass, we make direct observations of growth rate and substrate specialization in individual cells throughout a period of growth in a changing environment. This approach allows the flow of information from DNA to be constrained from the distal end of the regulatory cascade and will become an essential tool in the rapidly advancing field of single-cell metabolism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Using "AI Poincare" to analyze non-linear integrable optics

This study dives into the applicability of using automated discovery of conserved quantities in dynamical systems relevant to accelerator physics. Specifically, we explore the performance of AI Poincaré in analyzing numerical trajectory data obtained using the McMillan system of non-linear integrable optics. A comprehensive evaluation of the algorithm's performance is conducted through diverse methodologies. These include the analysis of the estimated number of conserved quantities embedded in a dataset and the deviation of interpolated points on the inferred manifold with respect to points in actually in the dataset. the investigation identifies an optimal range of perturbation distances where the underlying manifold extraction algorithm inside AI Poincaré exhibits optimal performance. Additionally, an improved neural network architecture is proposed based on the observed results. Finally, we apply the algorithm to preliminary experimental data from the Integrable Optics Test Accelerator at Fermilab to successfully infer the number of conserved quantities even in the presence of fast decoherence of the measured signal.

Osmanov, Lazare [Free U. Tbilisi]↗
DOE OSTI2025

A segmented approach to modeling building height: Delineating high-rise and low-rise buildings for enhanced height estimation

Understanding building height is imperative to the overall study of energy efficiency, population distribution, urban morphologies, emergency response, among others. Currently, existing approaches for modeling building height at scale are hindered by two pervasive issues. First, there is no consistent approach to quantify what a high-rise building is at a macro scale, leaving researchers unable to accurately compare results across geographies and domains. Second, high-rise buildings represent a small fraction of the built environment, implying data imbalance challenges that negatively affect current approaches. This is a problem of practical relevance since information on high-rise buildings is important for studies on urban heat islands, population dynamics, and pollution dispersion. Here, we introduce a novel approach to map building height which first identifies two distinct distributions within the built environment, with one being composed of low-rise buildings and one composed of high-rise buildings. We then develop an ensemble scheme where discrete specialist models are trained for each subset of low-rise buildings and high-rise buildings to infer building height from morphology features. For experiments mapping heights of 4.85 million buildings in Japan, we show an increase of 34 % in accuracy within 3m error when compared to the current state-of-the-art when modeling high-rise buildings, which based on KNN experimentation we define as any building > 12m . Our findings show that such an ensemble framework outperforms the current state-of-the-art approaches, which is especially relevant in relation to inferring height for high-rise buildings, a prominent issue of existing approaches for mapping the built environment.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Active learning with multifidelity modeling for efficient rare event simulation

Here, while multifidelity modeling provides a cost-effective way to conduct uncertainty quantification with computationally expensive models, much greater efficiency can be achieved by adaptively deciding the number of required high-fidelity (HF) simulations, depending on the type and complexity of the problem and the desired accuracy in the results. We propose a framework for active learning with multifidelity modeling emphasizing the efficient estimation of rare events. Our framework works by fusing a low-fidelity (LF) prediction with an HF-inferred correction, filtering the corrected LF prediction to decide whether to call the high-fidelity model, and for enhanced subsequent accuracy, adapting the correction for the LF prediction after every HF model call. The framework does not make any assumptions as to the LF model type or its correlations with the HF model. In addition, for improved robustness when estimating smaller failure probabilities, we propose using dynamic active learning functions that decide when to call the HF model. We demonstrate our framework using several academic case studies (including some high-dimensional problems) and two finite element model case studies: estimating Navier-Stokes velocities using the Stokes approximation and estimating stresses in a transversely isotropic model subjected to displacements via a coarsely meshed isotropic model. Across these case studies, not only did the proposed framework estimate the failure probabilities accurately, but compared with either Monte Carlo or a standard variance reduction method, it also required only a small fraction of the calls to the HF model.

42 ENGINEERING↗
DOE OSTI2022

The field of human building interaction for convergent research and innovation for intelligent built environments

Human-Building Interaction (HBI) is a convergent field that represents the growing complexities of the dynamic interplay between human experience and intelligence within built environments. This paper provides core definitions, research dimensions, and an overall vision for the future of HBI as developed through consensus among 25 interdisciplinary experts in a series of facilitated workshops. Three primary areas contribute to and require attention in HBI research: humans (human experiences, performance, and well-being), buildings (building design and operations), and technologies (sensing, inference, and awareness). Three critical interdisciplinary research domains intersect these areas: control systems and decision making, trust and collaboration, and modeling and simulation. Finally, at the core, it is vital for HBI research to center on and support equity, privacy, and sustainability. Compelling research questions are posed for each primary area, research domain, and core principle. State-of-the-art methods used in HBI studies are discussed, and examples of original research are offered to illustrate opportunities for the advancement of HBI research.

42 ENGINEERING↗
DOE OSTI2022

Functional Data Analysis for Extracting the Intrinsic Dimensionality of Spectra: Application to Chemical Homogeneity in the Open Cluster M67

High-resolution spectroscopic surveys of the Milky Way have entered the Big Data regime and have opened avenues for solving outstanding questions in Galactic archeology. However, exploiting their full potential is limited by complex systematics, whose characterization has not received much attention in modern spectroscopic analyses. In this work, we present a novel method to disentangle the component of spectral data space intrinsic to the stars from that due to systematics. Using functional principal component analysis on a sample of 18,933 giant spectra from APOGEE, we find that the intrinsic structure above the level of observational uncertainties requires ≈10 functional principal components (FPCs). Our FPCs can reduce the dimensionality of spectra, remove systematics, and impute masked wavelengths, thereby enabling accurate studies of stellar populations. To demonstrate the applicability of our FPCs, we use them to infer stellar parameters and abundances of 28 giants in the open cluster M67. We employ Sequential Neural Likelihood, a simulation-based Bayesian inference method that learns likelihood functions using neural density estimators, to incorporate non-Gaussian effects in spectral likelihoods. By hierarchically combining the inferred abundances, we limit the spread of the following elements in M67: Fe ≲ 0.02 dex; C ≲ 0.03 dex; O, Mg, Si, Ni ≲ 0.04 dex; Ca ≲ 0.05 dex; N, Al ≲ 0.07 dex (at 68% confidence). Our constraints suggest a lack of self-pollution by core-collapse supernovae in M67, which has promising implications for the future of chemical tagging to understand the star formation history and dynamical evolution of the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A Dataset of CFD Simulated Industrial Furnace Images for Conditional Automatic Generation with GANs

The steel industry is constantly looking for ways to automate processes and improve efficiency. A standard practice in industry is to simulate how complex systems will operate before they are actually used. Some complex systems, including steel industry processes such as blast furnaces, require complex physics-based simulations utilizing computational fluid dynamics (CFD). These CFD physics-based simulations are very accurate but can take significant time and computational resources to process, resulting in challenges for the implementation of the models in real-world operational environments. In recent years, deep learning (DL) has been considered as a substitute for these CFD models. DL models can be trained on validated CFD simulation data and then used for industrial process inference. Previous DL-based solutions have made great contributions for industrial automation but are currently missing the additional visualization component that CFD simulations also provide. In this paper, we propose a dataset for simple DL generative approaches that can help to address this issue. The dataset and methodology under development to approach this prediction are discussed in this work.

Calix, Ricardo↗
DOE OSTI2024

Defect-Limited Mobility and Defect Thermochemistry in Mixed A-Cation Tin Perovskites: (CH3NH3)1-xCsxSnBr3

Hybrid organic-inorganic semiconductors crystallizing in the perovskite structure present a significant opportunity for realizing defect-tolerant semiconductors. In this work, we examine the solid solution, (CH3NH3)1-xCsxSnBr3, and identify the thermochemistry dictating the intrinsic carrier concentrations and how local structural distortions influences this electronic behavior. This family of compounds exhibits the expected systematic trend in decreasing optical gap with the cesium to methylammonium A-site mixing ratio, x, in the visible region. However, the carrier mobility, as determined from time-resolved microwave conductivity measurements, trends opposite to that expected from first-principles calculations combined with Boltzmann scattering theory calculations of the carrier mobility. We propose that this is a result of increasing carrier scattering with x in (CH3NH3)1-xCsxSnBr3 due to a significant increase in the carrier density with x. By examining the dependence of the carrier density as a function of x, we infer the compositional-dependence of the average enthalpy and nonconfigurational entropy per defect. While diffraction reveals a cubic aristotypic perovskite structure as a function of x at room temperature, the pair distribution functions obtained from synchrotron X-ray total scattering from these materials are better described by symmetry-adapted displacement modes of the Pm3m crystal structure, which we attribute to a large degree of anharmonic dynamics of the atom positions. This analysis shows that CH3NH3+-rich compositions retain mostly linear Sn-Br-Sn bonding environments through the displacements, while Cs+-rich compositions lead to more significantly bent Sn-Br-Sn environments. We propose that these bent bonding arrangements in Cs-rich compositions yield a higher propensity for defect formation. This also provides a rationale for carrier trapping that gives rise to anomalous microwave transients. Together, these results provide insight into the structure-dynamics-properties relationships in this highly anharmonic system with high amplitude atomic motions and low defect formation energies.

defect thermochemistry↗
DOE OSTI2026

Informing Robust Functional Relationship Benchmarks: An Evaluation of the Temperature Sensitivity of Ecosystem Respiration Across the Arctic-Boreal Region

During land model development, simulated carbon dynamics are often benchmarked against observational data sets to evaluate model performance. Functional relationship benchmarks are the relationship between a driving variable (e.g., temperature) and a response variable (e.g., ecosystem respiration) and are a promising tool for assessing model performance by evaluating modeled sensitivities to changing environmental conditions. However, observed functional relationships can be influenced by choices made during data collection and throughout the benchmarking process, impacting the inferred skill of land models. To avoid misrepresenting a model's true performance, it is necessary to systematically evaluate best practices when constructing functional relationship benchmarks. We developed a set of guidelines for constructing functional relationship benchmarks, considering the choice of data set, number of daily observations, temporal extent, and temporal resolution across Alaska and Canada over a 20-year period from 2001 to 2020. The temperature sensitivity of ecosystem respiration from observations, evaluated through an apparent Q 10 , is highly variable both spatially and as a result of the data processing approach applied in the benchmark formation. When benchmarking 13 models from the Warming Permafrost Model Intercomparison Project (WrPMIP), the range in inferred model skill is substantially impacted by the choices applied in constructing functional relationship benchmarks. The inferred performance of a given model is most sensitive to the number of daily observations and temporal extent, followed by choice of benchmark data set and temporal averaging. Results from this analysis can guide the development of consistent and robust functional relationships for future model evaluation studies.

Poe, Jeralyn [Northern Arizona University, Flagsta↗
DOE OSTI2020

Measurement of Kinetic-Scale Current Filamentation Dynamics and Associated Magnetic Fields in Interpenetrating Plasmas

Here, we present the first local, quantitative measurements of ion current filamentation and magnetic field amplification in interpenetrating plasmas, characterizing the dynamics of the ion Weibel instability. The interaction of a pair of laser-generated, counterpropagating, collisionless, supersonic plasma flows is probed using optical Thomson scattering (TS). Analysis of the TS ion-feature revealed anticorrelated modulations in the density of the two ion streams at the spatial scale of the ion skin depth c/ω pi = 120 μm, and a correlated modulation in the plasma current. The inferred current profile implies a magnetic field amplitude ~30 ± 6 T , corresponding to ~1% of the flow kinetic energy, indicating that magnetic trapping is the dominant saturation mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

FPGA Acceleration of GCN in Light of the Symmetry of Graph Adjacency Matrix

Graph Convolutional Neural Networks (GCNs) are widely used to process large-scale graph data. Different from deep neural networks (DNNs), GCNs are sparse, irregular, and unstructured, posing unique challenges to hardware acceleration with regular processing elements (PEs). In particular, the adjacency matrix of a GCN is extremely sparse, leading to frequent but irregular memory access, low spatial/temporal data locality and poor data reuse. Furthermore, a realistic graph usually consists of unstructured data (e.g., unbalanced distributions), creating significantly different processing times and imbalanced workload for each node in GCN acceleration. To overcome these challenges, we propose an end-to-end hardware-software co-design to accelerate GCNs on resource-constrained FPGAs with the features including: (1) A custom dataflow that leverages symmetry along the diagonal of the adjacency matrix to accelerate feature aggregation for undirected graphs. We utilize either the upper or the lower triangular matrix of the adjacency matrix to perform aggregation in GCN to improve data reuse. (2) Unified compute cores for both aggregation and transform phases, with full support to the symmetry-based dataflow. These cores can be dynamically reconfigured to the systolic mode for transformation or as individual accumulators for aggregation in GCN processing. (3) Preprocessing of the graph in software to rearrange the edges and features to match the custom dataflow. This step improves the regularity in memory access and data reuse in the aggregation phase. Moreover, we quantize the GCN precision from FP32 to INT8 to reduce the memory footprint without losing the inference accuracy. We implement our accelerator design in Intel Stratix10 MX FPGA board with HBM2, and demonstrate 1.3x-110.5x improvement in end-to-end GCN latency as compared to the state-of the-art FPGA implementations, on the graph datasets of Cora, Pubmed, Citeseer and Reddit.

Nair, Gopikrishnan R.↗
DOE OSTI2023

Mapping Glacier Basal Sliding Applying Machine Learning

During the RESOLVE project (“High-resolution imaging in subsurface geophysics: development of a multi-instrument platform for interdisciplinary research”), continuous surface displacement and seismic array observations were obtained on Glacier d’Argentière in the French Alps for 35 days in May 2018. The data set is used to perform a detailed study of targeted processes within the highly dynamic cryospheric environment. In particular, the physical processes controlling glacial basal motion are poorly understood and remain challenging to observe directly. Especially in the Alpine region for temperate based glaciers where the ice rapidly responds to changing climatic conditions and thus, processes are strongly intermittent in time and heterogeneous in space. Spatially dense seismic and Global Positioning System (GPS) measurements are analyzed applying machine learning to gain insight into the processes controlling glacial motions of Glacier d’Argentière. Using multiple bandpass-filtered copies of the continuous seismic waveforms, we compute energy-based features, develop a matched field beamforming catalog and include meteorological observations. Features describing the data are analyzed with a gradient boosting decision tree model to directly estimate the GPS displacements from the seismic noise. We posit that features of the seismic noise provide direct access to the dominant parameters that drive displacement on the highly variable and unsteady surface of the glacier. The machine learning model infers daily fluctuations and longer term trends. The results show on-ice displacement rates are strongly modulated by activity at the base of the glacier. The techniques presented provide a new approach to study glacial basal sliding and discover its full complexity.

58 GEOSCIENCES↗
DOE OSTI2022

Characterization of self-magnetic pinch (SMP) radiographic diode performance on RITS-6 at Sandia National Laboratories. I. Diode dynamics, DC heating to extend radiation pulse

Radiographic diodes focus on an intense electron beam to a small spot size to minimize the source area of energetic photons for radiographic interrogation. The self-magnetic pinch (SMP) diode has been developed as such a source and operated as a load for the six-cavity radiographic integrated test stand (RITS-6) inductive voltage adder driver. While experiments support the generally accepted conclusion that a 1:1 aspect diode (cathode diameter equals anode–cathode gap) delivers optimum SMP performance, such experiments also show that reducing the cathode diameter, while reducing spot size, also results in reduced radiation dose, by as much as 50%, and degraded shot reproducibility. Analysis of the effective electron impingement angle on the anode converter with time made possible by a newly developed dose-rate array diagnostic indicates that fast-developing oscillations of the angle are correlated with early termination of the radiation pulse on many of the smaller-diameter SMP shots. This behavior as a function of relative cathode size persists through experiments with output voltages and currents up to 11.5 MV and 225 kA, respectively, and with spot sizes below approximately few millimeters. Since simulations to date have not predicted such oscillatory behavior, considerable discussion of the angle behavior of SMP shots is made to lend credence to the inference. We report there is clear anecdotal evidence that DC heating of the SMP diode region leads to stabilization of this oscillatory behavior. This is the first of two papers on the performance of the SMP diode on the RITS-6 accelerator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Particle control via cryopumping and its impact on the edge plasma profiles of Alcator C-Mod

At the high n e proposed for high-field fusion reactors, it is uncertain whether ionization, as opposed to plasma transport, will be most influential in determining n e at the pedestal and separatrix. A database of Alcator C-Mod discharges is analyzed to evaluate the impact of source modification via cryopumping. The database contains similarly-shaped H-modes at fixed I P = 0.8 MA and B t = 5.4 T, spanning a large range in P net and ionization. Measurements from an edge Thomson scattering system are combined with those from a midplane-viewing Ly α camera to evaluate changes to n e and T e in response to changes in ionization rates, S ion ∙ $n^{sep}_e$ and $T^{ped}_e$ are found to be most sensitive to changes to $S^{sep}_{ion}$, as opposed to $n^{ped}_e$ and $T^{sep}_e$. Dimensionless quantities, namely α MHD and v*, are found to regulate attainable pedestal values. Select discharges at different values of P net and in different pumping configurations are analyzed further using SOLPS-ITER. It is determined that changes to plasma transport coefficients are required to self-consistently model both plasma and neutral edge dynamics. Pumping is found to modify the poloidal distribution of atomic neutral density, n 0 , along the separatrix, increasing n 0 at the active X-point. Opaqueness to neutrals from high n e in the divertor is found to play a role in mediating neutral penetration lengths and hence, the poloidal distribution of neutrals along the separatrix. Pumped discharges thus require a larger particle diffusion coefficient than that inferred purely from 1D experimental profiles at the outer midplane.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Learned adaptive properties for mitigation of weight perturbations in embedded spiking networks

Recent years have seen an increased importance of neural network inference in edge-based scenarios, which impose size and power constraints requiring novel computing devices. These same edge scenarios may require operating over long periods of time, or exposure to extreme environments, resulting in a drift of neural network weights that cause degraded performance. In searching for ways to develop neural network approaches that perform robustly under these conditions, we propose a biologically-inspired mechanism for the dynamic adaptation of within-neuron parameters that is guided by a global context signal carrying information about perturbations and variability in incoming stimuli. Specifically, we demonstrate that adaptive voltage thresholds or neuronal time constants, when informed by a global context signal, can enable network-level mechanisms to recover from perturbed synaptic weights. Consistent with prior literature, the context-modulated approach is effective for recurrent, but not feedforward networks, by modulating network level dynamics. We demonstrate this approach successfully recovers performance in image classification tasks and spatiotemporal tracking tasks under idealized and Gaussian noise as well as for realistic perturbations from a memristive device when exposed to ionizing radiation. Finally, we discuss how this approach enables the design of robust and energy-efficient neuromorphic systems that perform well, even in resource-constrained scenarios with extreme environments such as edge processing.

context modulation↗
DOE OSTI2020

Constraining the Milky Way Mass Profile with Phase-space Distribution of Satellite Galaxies

We estimate the Milky Way (MW) halo properties using satellite kinematic data including the latest measurements from Gaia DR2. With a simulation-based 6D phase-space distribution function (DF) of satellite kinematics, we can infer halo properties efficiently and without bias, and handle the selection function and measurement errors rigorously in the Bayesian framework. Applying our DF from the EAGLE simulation to 28 satellites, we obtain an MW halo mass of $M={1.23}_{-0.18}^{+0.21}\times {10}^{12}{M}_{\odot }$ and a concentration of $c={9.4}_{-2.1}^{+2.8}$ with the prior based on the M–c relation. The inferred mass profile is consistent with previous measurements but with better precision and reliability due to the improved methodology and data. Potential improvement is illustrated by combining satellite data and stellar rotation curves. Using our EAGLE DF and best-fit MW potential, we provide much more precise estimates of the kinematics for those satellites with uncertain measurements. Compared to the EAGLE DF, which matches the observed satellite kinematics very well, the DF from the semi-analytical model based on the dark-matter-only simulation Millennium II (SAM-MII) over-represents satellites with small radii and velocities. We attribute this difference to less disruption of satellites with small pericenter distances in the SAM-MII simulation. Finally, by varying the disruption rate of such satellites in this simulation, we estimate a ~5% scatter in the inferred MW halo mass among hydrodynamics-based simulations.

79 ASTRONOMY AND ASTROPHYSICS↗