Search NASA⌕ Search

SEARCH · Search NASA

Results for “Observing System 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.

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

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The aluminium-26 distribution in a cosmological simulation of a Milky Way-type Galaxy

Context . The 1.8 MeV γ -rays corresponding to the decay of the radioactive isotope 26 Al (with a half-life of 0.72 Myr ) have been observed by the SPI detector on the INTEGRAL spacecraft and extensively used as a tracer of star formation and current nucleosynthetic activity in the Milky Way Galaxy. Further information is encoded in the observation related to the higher 26 Al content found in regions of the Galaxy with the highest line-of-sight (LoS) velocity relative to an observer located in the Solar System. However, this feature remains unexplained. Aims . We ran a cosmological “zoom-in” chemodynamical simulation of a Milky Way-type galaxy, including the production and decays of radioactive nuclei in a fully self-consistent way. We then analyzed the results to follow the evolution of 26 Al throughout the lifetime of the simulated galaxy to provide a new method for interpreting the 26 Al observations. Methods . We included the massive star sources of 26 Al in the Galaxy and its radioactive decay into a state-of-the-art galactic chemical evolution model, coupled with cosmological growth and hydrodynamics. This approach allowed us to follow the spatial and temporal evolution of the 26 Al content in the simulated galaxy. Results . Our results are in agreement with the observations with respect to the fact that gas particles in the simulation with relatively higher 26 Al content also have the highest LoS velocities. On the other hand, gas particles with relatively lower 26 Al content (i.e., not bright enough to be observed) generally display the lowest LoS velocities. However, this result is not conclusive because the overall rotational velocity of our simulated galaxy is higher than that observed for cold CO gas in the Milky Way Galaxy. Furthermore, we found no significant correlation between gas temperature, rotational velocity, and 26 Al content at any given radius. We also found the presence of transient 26 Al-rich spots at low LoS velocities and we show that one such spot had been captured by the INTEGRAL/SPI data. Based on our model, we present a prediction for the detection of 1.8 MeV γ -rays by the future COSI mission. We find that according to our model, the new instrument will be able to observe similar 26 Al-emission patterns to those seen by INTEGRAL/SPI.

Wehmeyer, B. (ORCID:0009000254149292)↗
DOE OSTI2024

Medium-assisted enhancement of χ c 1 ( 3872 ) production from small to large colliding systems

Studies of exotic hadrons such as the χ c 1 ( 3872 ) state provide crucial insights into the fundamental force governing the strong interaction dynamics, with an emerging frontier to investigate their production in high energy collisions where a partonic medium is present. The latest experimental measurements from the Large Hadron Collider show an intriguing evolution pattern of the χ c 1 ( 3872 ) -to- ψ ( 2 S ) yield ratio from proton-proton collisions with increasing multiplicities toward proton-lead and lead-lead collisions. Here we propose a mechanism of medium-assisted enhancement for the χ c 1 ( 3872 ) production, which competes with the more conventional absorption-induced suppression and results in a nonmonotonic trend from small to large colliding systems. Realistic simulations from this model offer a quantitative description of all available data. Predictions are made for the centrality dependence of this observable in PbPb collisions as well as for its system-size dependence from OO and ArAr to XeXe and PbPb collisions. In both cases, a nonmonotonic behavior emerges as the imprint of the competition between enhancement and suppression and can be readily tested by future data. Published by the American Physical Society 2024

Physics↗
DOE OSTI2026

Connecting Current and Future Dual Active Galactic Nucleus Searches to LISA and Pulsar Timing Array Gravitational-wave Detections

Abstract Dual active galactic nuclei (DAGN) mark the observable stage of massive black hole (MBH) pairing during galaxy mergers and are the progenitors of the MBH binaries that generate low-frequency gravitational waves. Using the large-volume ASTRID cosmological simulation, we construct mock DAGN catalogs tailored to the selection functions of current (COSMOS-Web and DESI) and forthcoming (Nancy Grace Roman Space Telescope (Roman) and Lynx X-ray Observatory (Lynx)) surveys, enabling direct comparisons between simulations and observations. With realistic observational selections, ASTRID reproduces the observed dual fractions, projected separations, and host-galaxy properties across redshifts. We predict a substantial population of small-separation (<5 kpc) duals that remain inaccessible to current surveys, demonstrating that the apparent paucity of subkiloparsec systems in COSMOS-Web is primarily a consequence of observational selection rather than an intrinsic absence. Following each simulated dual to coalescence, we show that DAGN are effective tracers of MBH mergers: ∼30%–70% merge within ≲1 Gyr, and 20%–60% of these mergers produce gravitational-wave signals detectable by the Laser Interferometer Space Antenna (LISA). Duals observable with Roman and Lynx are the progenitors of ∼10%–50% of low-redshift LISA sources and contribute ∼30% of the PTA-band stochastic gravitational-wave background. We further identify massive green-valley galaxies hosting moderate-luminosity active galactic nuclei (AGN), together with massive star-forming galaxies containing bright quasars at z > 1, as the environments most likely to host imminent MBH binaries. These results establish a unified cosmological framework connecting DAGN demographics, MBH binary evolution, and gravitational-wave sources, while identifying high-priority targets for coordinated electromagnetic and multimessenger observations in the coming decade.

Chen, Nianyi [Max-Planck-Institut für Astrophysik;↗
DOE OSTI2025

The Hunga Volcanic Eruption Atmospheric Impacts Report

On 15 January 2022 a highly explosive eruption of the Hunga volcano occurred in the Kingdom of Tonga in the South Pacific Ocean (175°24’ W, 20°33’ S). The Volcanic Explosivity Index (VEI) 6 eruption originated from a shallow submarine vent, making it distinct from large subaerial eruptions of recent decades (e.g., 1982 El Chichón, 1991 Mt. Pinatubo). In particular, seawater enhanced explosivity and dampened sulfur dioxide (SO 2 ) emissions. The eruption was the culmination of ~1 month of precursory activity; however, the timing and size of the eruption were unexpected, partly due to the challenges of monitoring submarine volcanoes. The stratospheric hydration caused by the eruption was unprecedented in magnitude, altitude, and duration in the satellite record. This Executive Summary reflects the current assessment of the Hunga eruption and its impact on the climate system. We report key observations of the eruption and its aftermath, as well as simulations of its impact by global chemistry-climate models. The Hunga eruption had an unprecedented impact on the stratosphere and mesosphere due to the plume height and large water content, which increased the global stratospheric water vapour burden by 10%. Most of this water has remained in the atmosphere into 2025. However, Hunga’s net impact on surface climate was small compared to that of earlier large-magnitude volcanic eruptions, due to limited sulfate aerosol loading in the stratosphere and the high altitude of the water vapour injection.

58 GEOSCIENCES↗
DOE OSTI2026

Ab Initio Polariton Transport Dynamics with the Classical Path Approximation

We present an ab initio framework for simulating polariton transport dynamics based on the classical path approximation (CPA). The quantum dynamics of polariton transport involves simulating many electronic degrees of freedom, making a fully ab initio dynamics simulation computationally expensive. We demonstrate that the CPA, which removes the need for excited-state nuclear gradients, is well-suited for polaritonic systems because collective light–matter coupling leads to vanishing excited-state forces. Benchmark comparisons between CPA and full evaluation of the excited-state forces show excellent agreement for polariton transport results in model light–matter systems such as polariton group velocities and mean-squared displacements. Ab initio simulations of polariton transport using CPA reproduce key physical trends that are observed in experiments with BODIPY molecules. Our work establishes the CPA as a highly efficient tool for ab initio investigations of transport and energy flow in hybrid light–matter systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

SIDM Concerto: Compilation and Data Release of Self-interacting Dark Matter Zoom-in Simulations

We present SIDM Concerto: 14 cosmological zoom-in simulations in cold dark matter (CDM) and self-interacting dark matter (SIDM) models based on the Symphony and Milky Way-est suites. SIDM Concerto includes one Large Magellanic Cloud– (LMC-) mass system (host mass ∼10 11 M ⊙ ), two Milky Way (MW) analogs (∼10 12 M ⊙ ), two group-mass hosts (∼10 13 M ⊙ ), and one low-mass cluster (∼10 14 M ⊙ ). Each host contains ≈2 × 10 7 particles and is run in CDM and one or more strong, velocity-dependent SIDM models. Our analysis of SIDM (sub)halo populations over seven subhalo mass decades reveals that (1) the fraction of core-collapsed isolated halos and subhalos peaks at a maximum circular velocity corresponding to the transition of the SIDM cross section from a v −4 to v 0 scaling; (2) SIDM subhalo mass functions are suppressed by ≈50% relative to CDM in LMC, MW, and group-mass hosts but are consistent with CDM in the low-mass cluster host; (3) subhalos’ inner density profile slopes, which are more diverse in SIDM than in CDM, are sensitive to both the amplitude and shape of the SIDM cross section. Our simulations provide a benchmark for testing SIDM predictions with astrophysical observations of field and satellite galaxies, strong lensing systems, and stellar streams. Data products are publicly available at doi:10.5281/zenodo.14933624.

dark matter↗
DOE OSTI2025

Warm-phase microphysical evolution in large-eddy simulations of tropical cumulus congestus: evaluating drop size distribution evolution using polarimetry retrievals, in situ measurements, and a thermal-based framework

Owing to uncertainties in convective microphysics processes, improving parameterizations in Earth system models (ESMs) can benefit from observationally constrained cases suitable for scaling between cloud-resolving models and ESMs. We propose a benchmark large-eddy simulation (LES) cumulus congestus case study from the NASA Cloud, Aerosol, and Monsoon Processes Philippines Experiment (CAMP 2 Ex) for evaluating and improving ESMs in single-column model (SCM) mode. We seek observational constraints using novel polarimetric retrievals and in situ cloud microphysics measurements. Simulations using bulk and bin microphysics initialized with observed aerosol profiles are compared to cloud-top retrievals of cloud droplet effective radius (R eff ), effective variance (ν eff ), and number concentration (N d ) from the airborne Research Scanning Polarimeter (RSP). Both schemes reproduce characteristics of cloud-top N d and R eff that increase and decrease with altitude, respectively. Cloud-top N d is low-biased relative to RSP retrievals in both schemes, potentially due to limitations in both simulations and retrieval assumptions. Cloud-top R eff is low-biased in the bulk scheme but reasonably reproduced by the bin scheme. Profiles of N d and R eff are sensitive to the collision–coalescence process and the vertical variation in aerosol size distribution. Comparison of simulated and in situ droplet size distributions (DSDs) shows that, to first order, integrated moments are always sensitive to sizes < ~ 30 µm and can also be sensitive to larger sizes if the DSDs are sufficiently broad, with implications for the assumed maximum observed size retrieved by the RSP. The bin scheme captures the observed extended tail of the DSD, while the bulk scheme is unable to due to parametric constraints. Differences in expected relationships between in situ measurements of cloud cores and cloud-top retrievals by RSP demonstrate difficulty in constraining well the case presented herein. Finally, a thermal-tracking framework demonstrates that the dilution of N d throughout a thermal's lifetime is heavily determined by collision–coalescence and the height-varying aerosol distribution and that, in the absence of these, the impact of entrainment on diluting N d is largely offset by secondary activation. Implications for evaluating warm-phase convective microphysics schemes in ESMs and translating results for use on global, space-based polarimetry platforms are discussed.

Stanford, McKenna Wallace [Columbia Univ., New Yor↗
DOE OSTI2024

Phonons reveal coupled cholesterol-lipid dynamics in ternary membranes

Experimental studies of collective dynamics in lipid bilayers have been challenging due to the energy resolution required to observe these low-energy phonon-like modes. However, inelastic X-ray scattering (IXS) measurements — a technique for probing vibrations in soft and biological materials — are now possible with sub-meV resolution, permitting direct observation of low energy, phonon-like modes in lipid membranes. Here, in this study, IXS measurements with sub-meV energy resolution reveal a low-energy optic-like phonon mode at roughly 3 meV in the liquid-ordered (L o ) and liquid-disordered (L d ) phases of a ternary lipid mixture. This mode is only observed experimentally at momentum transfers greater than 5 nm -1 in the L o system. A similar gapped mode is also observed in all-atom molecular dynamics (MD) simulations of the same mixture, indicating that the simulations accurately represnt the fast, collective dynamics in the L o phase. It’s optical nature and the Q range of the gap together suggest that the observed mode is due to the coupled motion of cholesterol-lipid pairs, separated by several hydrocarbon chains within the membrane plane. Analysis of the simulations provide molecular insight into the origin of the mode in transient, nanoscale substructures of hexagonally packed hydrocarbon chains. This nanoscale hexagonal packing was previously reported based on molecular dynamics simulations and later by NMR measurements. Here, however, the integration of IXS and MD simulations identifies a new signature of the L o substructure in the collective lipid dynamics, thanks to the recent confluence of IXS sensitivity and MD simulation capabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Optimization-Based Data-Driven Approach for Detecting Fault Location in Power Systems

In grids with large penetration of converterinterfaced resources (CIRs), measurements of voltage, current, and line parameters can fluctuate significantly during fault conditions. These fluctuations, combined with complex network topologies and extensive system branching, make accurate fault location challenging. Faults, such as short circuits, can cause prolonged outages with serious socio-economic impacts, highlighting the need for rapid fault identification to minimize downtime. However, current fault detection methods—such as relays and digital fault recorders—often relay information too slowly, impeding swift corrective action. Given the limited availability of high-resolution phasor measurement units, this paper introduces an optimization-based observer to estimate fault locations, grid line parameters, and voltages using local CIR measurements. To preserve the confidentiality of CIRs and enhance estimation accuracy, this study uses a black-box model of CIRs. This bottom-up, event-driven approach can enhances protection and control systems through optimized and real-time fault detection. Simulation results show that the optimization-based data-driven observer can accurately detect fault locations and estimate grid states and parameters, providing valuable insights for utilities and operators in grid applications.

Subedi, Sunil [ORNL] (ORCID:000000034069090X)↗
DOE OSTI2026

Observational Data for Next-Generation Climate Model Evaluation: Requirements, Considerations, and Best Practices

Climate model simulations are an important source of information about our planet’s climate system and also enable informed decision-making under different future scenarios. As a new archive of results from the next generation of climate models is anticipated to become available with the Coupled Model Intercomparison Project phase 7 (CMIP7), the need to develop efficient and robust methods to evaluate models is paramount. Observations are an integral part of model evaluation, providing a means to quantify and understand the degree to which climate models can faithfully reproduce Earth system processes. Such analysis is critical for constraining climate projections, identifying areas of focus for model development, and assisting analysts in deciphering the utility of models for specific applications. Observations of Earth system come from a diversity of sources, span different space–time domains, and are produced by different communities, and each dataset features different data structures and formats, metadata standards, and its own unique uncertainties. Uncertainties in an observational dataset may stem from gaps in temporal and spatial coverage, instrumentation errors, or assumptions in retrieval and processing methods. How then does one ensure that observational data are ready for use and utilized in the most appropriate way for robust, rapid, and routine climate model evaluation? The CMIP7 Model Benchmarking Task Team with input from the broader climate modeling, model evaluation, and observational data communities present a vision and considerations for best practices toward the optimal and appropriate use of observational data to support next-generation climate model evaluation.

Climate models↗
DOE OSTI2026

Understanding the Influence of Chain Architecture on the Transport Quantities of Polymer Electrolytes with Covalently Bonded Anions

Here, we use a combination of experiments and coarse-grained molecular dynamics simulations to elucidate the structure–property relationships in polymer electrolytes obtained by the copolymerization of poly(vinyl ethylene carbonate─lithium styrene bis(trifluoromethanesulfonyl)imide) or p(VEC-LiSTFSI). Experiments show that the conductivity reduces with increasing anion (i.e., STFSI) fraction on the chain, and the cation transference number (t + ) is found to be dependent on the anion fraction. Furthermore, a significant fraction of unpolymerized VEC monomers are observed. Since it is inherently difficult to experimentally control the chain architecture and the amount of unpolymerized VEC in these systems, we perform coarse-grained molecular dynamics simulations on model polymer systems with different chain architectures to mimic the plausible experimental systems. Specifically, we look at the differences in transference numbers arising from (i) a random copolymer of VEC and STFSI monomers; (ii) a blend of VEC-STFSI copolymer with VEC monomers; and (iii) a ternary blend of the VEC homopolymer, STFSI homopolymers, and VEC monomers. The ternary blend model demonstrates the closest resemblance with the experimental transference numbers and diffusivities. The lithium diffusivity obtained from the coarse-grained models with VEC monomers (plasticizers) is about 1.5 times that of the model without VEC monomers, showing that the plasticizing effect of VEC monomers is modest. We rationalize the experimental observations based on aggregate and cluster analyses obtained from molecular simulations. This work reveals that polymer electrolyte chain architecture and plasticizers can critically influence the transport properties, and these parameters should be considered when designing single ion conducting polymeric electrolytes.

cluster distribution↗
DOE OSTI2025

Hydrodynamic Predictions for the Next Outburst of T Coronae Borealis: It Will Be the Brightest Classical or Recurrent Nova Ever Observed in X-Rays*

T Coronae Borealis (TCrB) is a recurrent nova with recorded outbursts in 1866 and 1946 and possible outbursts in 1217 and 1787. It is predicted to explode again in 2025 or 2026 based on multiple observational studies. The system consists of a massive ( M wd ≳ 1.35 M ⊙ ) white dwarf (WD) and a red giant (M3–M4 III). We have performed 1D hydrodynamic simulations with NOVA to predict the behavior of the next outburst. These simulations consist of a range of mass accretion rates onto ∼1.35 M ⊙ WDs, designed to bound the conditions necessary to achieve ignition of an explosion after an ≈80 yr interoutburst period. We have used both carbon–oxygen and oxygen–neon initial compositions, in order to include the possible ejecta abundances to be measured in the observations of the next outburst. As the WD in the TCrB system is observed to be massive, theoretical predictions reported here imply that the WD is growing in mass as a consequence of the thermonuclear runaway. Therefore, the secular evolution of the WD may allow it to approach the Chandrasekhar limit and either explode as a Type Ia supernova or undergo accretion-induced collapse, depending on its underlying composition. We have followed the evolution of just the WD, after removing the ejected matter from the surface layers. Our intent is to illuminate the mystery of the unique, second maximum in the two well-observed outbursts and we have found conditions that bracket the predictions.

Cataclysmic variable stars↗
DOE OSTI2025

Observable measurement-induced transitions

One of the main postulates of quantum mechanics is that measurements destroy quantum coherence (wave function collapse). Recently it was discovered that in a many-body system dilute local measurements still preserve some coherence across the entire system. As the measurement density is increased, a phase transition occurs that is characterized by the disentanglement of different parts of the system. Unfortunately, this transition is impossible to observe experimentally for macroscopic systems because it requires an exponentially costly full tomography of the many-body wave function or a comparison with the simulation on an oracle classical computer. In this work we report the discovery of another measurement-induced phase transition that can be observed experimentally if quantum dynamics can be reversed. On one side of this phase transition the quantum information encoded in some part of the Hilbert space is fully recovered after the time inversion. On the other side, all quantum information is corrupted. This transition also manifests itself as the change in the behavior of the probability to observe the same measurement outcome in the process that consists of identical blocks repeated many times. In each block the unitary evolution is followed by the measurement. On one side of the transition the probability decreases exponentially with the number of repetitions, on the other it tends to a constant as the number of repetitions is increased. We confirm the existence of the proposed phase transition through numerical simulations of realistic quantum circuits and analytical calculations using an effective random-matrix theory model.

Measurement-induced phase transitions↗
DOE OSTI2026

Importance of Spatially Continuous Urban Surface Properties in Urban‐Resolving Earth System Modeling

Accurate representation of urban properties and processes at higher resolutions in global modeling systems is essential for advancing our ability to capture the complexities of urban systems and informing effective resilience strategies. However, the prescription of coarse global-scale urban properties in most state-of-the-art Earth system models (ESMs) is limiting their potential for capturing urban signals as they advance toward kilometer-scale simulation capabilities. To bridge this gap in inadequate urban property representation and to advance urban-resolving Earth system modeling, this work integrates the newly-developed global 1 km-resolution facet-level urban surface property data set, U-Surf, into the land component of Community Earth System Model (CESM)—Community Terrestrial System Model (CTSM). The land-only CTSM simulations are validated against satellite measurements, ground-based urban weather stations, flux tower observations, and reanalysis data. Results demonstrate that the enhanced urban properties allow improved simulations of urban meteorology and surface energy fluxes compared to the default coarse-resolution categorical urban canopy parameters. Spatial scaling analysis reveals regime-dependent information loss during resolution aggregation, as well as substantial scale-dependent variations in urban surface energy flux representation. Furthermore, these findings have critical implications for coupled Earth system modeling when including the effect of land-atmosphere interaction. This work establishes a foundation for future urban-resolving kilometer-scale ESM development, which will enable systematic intra- and inter-city comparisons that inform urban adaptation strategies across diverse global urban environments.

Cheng, Yifan [University of Illinois Urbana-Champa↗
DOE OSTI2025

Effect of Rocky Mountains and Tibetan Plateau 1998 Spring Land Temperature on N. American and East Asian Summer Precipitation Anomalies

This work follows up on the GEWEX/LS4P Phase I (LS4P-I) experiments, a community effort highlighting the spring land surface temperature anomalies in the Tibetan Plateau (TP) as a useful source for subseasonal to seasonal (S2S) prediction of summer precipitation in global hot spot regions, particularly in East Asia and North America. This paper extends the investigation to both the US Rocky Mountain (RM) region and the TP, considering the 1998 summer drought/flood event in North America/East Asia, respectively, as a case study. A previously developed initialization method for land surface temperature/subsurface temperature (LST/SUBT) is used in the NCEP Global Forecast System, coupled with a land model, SSiB2 (GFS/SSiB2), to produce observed RM cold May temperature anomaly. Forward simulation yields June precipitation anomalies at five remote locations. Likewise, the TP warm May temperature anomaly also produces June precipitation anomalies at these five locations. The effects of RM (cold) and TP (warm) temperature anomalies are consistent in the US South Coastal regions and the south Yangtze River Basin, yielding 49% (42%) of observed drought and 34% (44%) of observed flood, respectively. These LST/SUBT effects in RM and TP induce a global large-scale wave train linking North America with the TP, affecting the subtropical westerly jet and thereby modulating summer precipitation. Global SST effect is examined for comparison but does not yield statistically significant June precipitation anomalies in GFS/SSiB2. Furthermore, this study adds to evidence that high-mountain LST effects in the RM and TP are first-order sources of S2S precipitation predictability in summer months.

Nayak, Hara Prasad [University of California, Los ↗
DOE OSTI2024

Foreword: Message From the MTS Journal Editor

Here, it is my pleasure to present this general issue to the MTS community. This general issue features manuscripts on a variety of topics of interest to the MTS community including articles, reviews, and commentaries: technology gaps for monitoring birds and marine mammals at offshore wind facilities; adaptive tracking of the Barents Sea polar front using an autonomous underwater vehicle; advanced simulation and environmental impact assessment of combustion in maritime energy systems; performance of moored real time ocean observations during Cyclones in the Bay of Bengal; recent advances in assessing environmental effects of marine renewable energy; and the new blue economy and future for marine renewable energy.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Investigating SU(3) with Nf=8 fundamental fermions at strong renormalized coupling

Lattice simulations have observed a novel strong coupling symmetric mass generation (SMG) phase for the SU(3) gauge system with $N_f=8$ fundamental fermions (represented by two sets of staggered fields) at very large renormalized coupling ($g^2_{GF} \gtrsim 25$). The results of Phys.Rev.D 106 (2022) 014513 suggest that the SMG phase is separated from the weak coupling, conformal phase by a continuous phase transition, implying that the SMG phase exists in the continuum limit. To scrutinize these findings, we are generating a set of large volume zero temperature ensembles using nHYP improved staggered fermions with additional Pauli-Villars fields to tame gauge field fluctuations. We consider the low-lying meson spectrum and verify the existence of the SMG phase. Based on a finite size scaling analysis we predict that the phase transition between the strong and weak coupling phases is likely governed by a merged fixed point that is ultraviolet in the strong coupling but infrared in the weak coupling side. This finding suggests that the SU(3) 8-flavor system sits at the opening of the conformal window

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

A Methodology to Evaluate the Grid Reliability Impact of Oscillations Induced by Large Loads

The rapid growth of hyperscale AI data centers is bringing renewed attention to the reliability risk that sustained forced oscillations pose to bulk power systems, with cyclic computational workloads emerging as a new forcing source. Unlike the broadband, stochastic disturbances from traditional industrial loads such as arc furnaces, AI training and inference facilities can inject large active power swings concentrated at specific frequencies over extended durations - characteristics that existing grid planning practices do not account for. While the North American Electric Reliability Corporation (NERC) has recognized this gap and called for system-level studies of large load interconnections, no standardized methodology exists to screen, simulate, and quantify these risks at the planning stage. This report presents the Risk Assessment Tool for Large Load-induced Events (RATLLE), a Python-based, publicly available script suite developed at the Pacific Northwest National Laboratory to evaluate bulk power system reliability risks from data center-induced oscillations. RATLLE implements a three-module workflow: a screening module that identifies vulnerable interconnection locations and excitable system modes; a simulation module that models cyclic data center load behavior using a commercial positive sequence simulation platform; and an analysis module that computes risk metrics and generates interactive visualization dashboards. The risk metrics, formulated around simulation observables, map oscillation impacts to a three-stage severity scale spanning latent equipment fatigue through imminent cascading failure. The methodology is demonstrated on two Western Electricity Coordinating Council (WECC) system models: a publicly available 240-bus reduced representation and a detailed 2031 Heavy Winter planning case. Case studies illustrate that even modest 50 MW forced oscillations at resonant frequencies can produce wide-area power swings, N-1 security constraint violations, and cascading generator trips through protection actions - outcomes that would not occur under normal operating conditions without oscillations present. The results underscore the need for standardized oscillation impact assessment in large load interconnection studies and provide a reproducible, extensible framework for utilities to adopt or customize within their existing planning workflows.

Biswas, Shuchismita↗