Search NASA⌕ Search

SEARCH · Search NASA

Results for “Distributed Estimation”

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 1,171 records · Page 65

Characterization of the Nucleus, Morphology, and Activity of Interstellar Comet 2I/Borisov by Optical and Near-infrared GROWTH, Apache Point, IRTF, ZTF, and Keck Observations

We present visible and near-infrared (NIR) photometric and spectroscopic observations of interstellar object (ISO) 2I/Borisov taken from 2019 September 10 to 2019 December 20 using the GROWTH, the Apache Point Observatory Astrophysical Research Consortium 3.5 m, and the NASA Infrared Telescope Facility 3.0 m combined with pre- and postdiscovery observations of 2I obtained by the Zwicky Transient Facility from 2019 March 17 to 2019 May 5. Comparison with imaging of distant solar system comets shows an object very similar to mildly active solar system comets with an outgassing rate of ∼10{sup 27} mol s{sup −1}. The photometry, taken in filters spanning the visible and NIR range, shows a gradual brightening trend of ∼0.03 mag day{sup −1} since 2019 September 10 UTC for a reddish object becoming neutral in the NIR. The light curve from recent and prediscovery data reveals a brightness trend suggesting the recent onset of significant H{sub 2}O sublimation with the comet being active with super volatiles such as CO at heliocentric distances >6 au consistent with its extended morphology. Using the advanced capability to significantly reduce the scattered light from the coma enabled by high-resolution NIR images from Keck adaptive optics taken on 2019 October 4, we estimate a diameter for 2I’s nucleus of ≲1.4 km. We use the size estimates of 1I/’Oumuamua and 2I/Borisov to roughly estimate the slope of the ISO size distribution, resulting in a slope of ∼3.4 ± 1.2, similar to solar system comets and bodies produced from collisional equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗

Velocity segregation and systematic biases in velocity dispersion estimates with the SPT-GMOS spectroscopic survey

The velocity distribution of galaxies in clusters is not universal; rather, galaxies are segregated according to their spectral type and relative luminosity. We examine the velocity distributions of different populations of galaxies within 89 Sunyaev Zel'dovich (SZ) selected galaxy clusters spanning $ 0.28 < z < 1.08$. Our sample is primarily draw from the SPT-GMOS spectroscopic survey, supplemented by additional published spectroscopy, resulting in a final spectroscopic sample of 4148 galaxy spectra---2868 cluster members. The velocity dispersion of star-forming cluster galaxies is $17\pm4$% greater than that of passive cluster galaxies, and the velocity dispersion of bright ($m < m^{*}-0.5$) cluster galaxies is $11\pm4$% lower than the velocity dispersion of our total member population. We find good agreement with simulations regarding the shape of the relationship between the measured velocity dispersion and the fraction of passive vs. star-forming galaxies used to measure it, but we find a small offset between this relationship as measured in data and simulations in which suggests that our dispersions are systematically low by as much as 3\% relative to simulations. We argue that this offset could be interpreted as a measurement of the effective velocity bias that describes the ratio of our observed velocity dispersions and the intrinsic velocity dispersion of dark matter particles in a published simulation result. Here, by measuring velocity bias in this way suggests that large spectroscopic surveys can improve dispersion-based mass-observable scaling relations for cosmology even in the face of velocity biases, by quantifying and ultimately calibrating them out.

79 ASTRONOMY AND ASTROPHYSICS↗

Cold Molecular Gas Along the Merger Sequence in Local Luminous Infrared Galaxies

We present an initial result from the {sup 12}CO (J = 1–0) survey of 79 galaxies in 62 local luminous and ultraluminous infrared galaxy (LIRG and ULIRG) systems obtained using the 45 m telescope at the Nobeyama Radio Observatory. This is a systematic {sup 12}CO (J = 1–0) survey of the Great Observatories All-sky LIRGs Survey (GOALS) sample. The molecular gas mass of the sample is in the range 2.2×10{sup 8}--7.0×10{sup 9} M{sub ⊙} within the central several kiloparsecs subtended by the 15{sup ′′} beam. A method to estimate the size of a CO gas distribution is introduced, which is combined with the total CO flux in the literature. This method is applied to part of our sample, and we find that the median CO radius is 1–4 kpc. From the early stage to the late stage of mergers, we find that the CO size decreases while the median value of the molecular gas mass in the central several-kiloparsec region is constant. Our results statistically support a scenario where molecular gas inflows toward the central region from the outer disk to replenish gas consumed by starburst, and that such a process is common in merging LIRGs.

79 ASTRONOMY AND ASTROPHYSICS↗

Filaments of the Slime Mold Cosmic Web and How They Affect Galaxy Evolution

Abstract We present a novel approach for identifying cosmic web filaments within theDisPerSEstructure identification framework, using cosmic density field estimates from the Monte Carlo Physarum Machine (MCPM), inspired by the slime mold organism. We apply our method to the IllustrisTNG (TNG100) cosmological simulations and investigate the impact of filaments on galaxies. The MCPM density field is superior to the Delaunay tessellation field estimator in tracing the true underlying matter distribution and allows filaments to be identified with higher fidelity, finding more low-prominence/diffuse filaments. Using our new filament catalogs, we find that ≳90% of galaxies are located within ∼1.5 Mpc of a filamentary spine, with little change in the median star formation activity with distance to the nearest filament. Instead, we uncover a differential effect of the local filament line density, Σ fil (MCPM)—the total MCPM overdensity per unit length along a filament segment—on galaxy formation: most galaxies are quenched and gas-poor near high-line density filaments atz≤ 1. At earlier times, the filamentary environment appears to have no effect on galactic gas supply and quenching. Atz= 0, quenching in log ( M * / M ⊙ ) ≳ 10.5 galaxies is mainly driven by mass, while lower-mass galaxies are significantly affected by the filament line density. Satellites are far more susceptible to filaments than centrals. The local environments of massive halos are not sufficient to account for the effect of filament line density on gas removal and quenching. Our new approach holds great promise for observationally identifying filaments from galaxy surveys such as SDSS and DESI.

Astronomy & Astrophysics↗

A derecho climatology (2004–2021) in the United States based on machine learning identification of bow echoes

Due to their persistent widespread severe winds, derechos pose significant threats to human safety and property, with impacts comparable to many tornadoes and hurricanes. Yet, automated detection of derechos remains challenging due to the absence of spatiotemporally continuous observations and the complex criteria employed to define the phenomenon. This study presents an objective derecho detection approach capable of automatically identifying derechos through both observations and model results. The approach is grounded in a physically based definition of derechos and integrates three algorithms: (1) the Python Flexible Object Tracker (PyFLEXTRKR) algorithm to track mesoscale convective systems (MCSs), (2) a semantic segmentation convolutional neural network to identify bow echoes, and (3) a comprehensive classification algorithm to detect derechos within MCS life cycles and distinguish derecho-producing from non-derecho-producing MCSs. Using this approach, we developed a novel high-resolution (4 km and hourly) observational dataset of derechos and accompanying derecho-producing MCSs over the United States east of the Rocky Mountains from 2004 to 2021. The dataset consists of two subsets based on different gust speed data sources and is analyzed to document the climatology of derechos in the United States. On average, 12–15 derechos are identified per year, aligning with previous estimations (∼6–21 events annually). The spatial distribution and seasonal variation patterns are consistent with prior studies, showing peak occurrences in the Great Plains and the Midwest during the warm season. Additionally, during the study period, derechos account for approximately 3.1 % of measured damaging gusts (≥25.93 m s−1) over the eastern United States. The dataset is publicly available at https://doi.org/10.5281/zenodo.14835362 (Li et al., 2025).

54 ENVIRONMENTAL SCIENCES↗

Using neural network ensembles to separate ocean biogeochemical and physical drivers of phytoplankton biogeography in Earth system models

Abstract. Earth system models (ESMs) are useful tools for predicting and understanding past and future aspects of the climate system. However, the biological and physical parameters used in ESMs can have wide variations in their estimates. Even small changes in these parameters can yield unexpected results without a clear explanation of how a particular outcome was reached. The standard method for estimating ESM sensitivity is to compare spatiotemporal distributions of variables from different runs of a single ESM. However, a potential pitfall of this method is that ESM output could match observational patterns because of compensating errors. For example, if a model predicts overly weak upwelling and low nutrient concentrations, it might compensate for this by allowing phytoplankton to have a high sensitivity to nutrients. Recently, we demonstrated that neural network ensembles (NNEs) are capable of extracting relationships between predictor and target variables within ocean biogeochemical models. Being able to view the relationships between variables, along with spatiotemporal distributions, allows for a more mechanistically based examination of ESM outputs. Here, we investigated whether we could apply NNEs to help us determine why different ESMs produce different spatiotemporal distributions of phytoplankton biomass. We tested this using three cases. The first and second case used different runs of the same ESM, except that the physical circulations differed between them in the first case, while the biological equations differed between them in the second. Our results indicated that the NNEs were capable of extracting the relationships between variables for different runs of a single ESM, allowing us to distinguish between differences due to changes in circulation (which do not change relationships) from changes in biogeochemical formulation (which do change relationships). In the third case, we applied NNEs to two different ESMs. The results of the third case highlighted the capability of NNEs to contrast the apparent relationships of different ESMs and some of the challenges it presents. Although applied specifically to the ocean components of an ESM, our study demonstrates that Earth system modelers can use NNEs to separate the contributions of different components of ESMs. Specifically, this allows modelers to compare the apparent relationships across different ESMs and observational datasets.

54 ENVIRONMENTAL SCIENCES↗

Influences of lidar scanning parameters on wind turbine wake retrievals in complex terrain

Abstract. Scanning lidars enable the collection of spatially distributed measurements of turbine wakes and the estimation of wake properties such as magnitude, extent, and trajectory. Lidar-based characterizations, however, may be subject to distortions due to the observational system. Distortions can arise from the resolution of the measurement points across the wake, the projection of the winds onto the beam, averaging along the beam probe volume, and intervening evolution of the flow over the scan duration. Using a large-eddy simulation and simulated measurements with a virtual lidar model, we assess how scanning lidar systems may influence the properties of the retrieved wake using a case study from the Perdigão campaign. We consider three lidars performing range-height indicator sweeps in complex terrain, based on the deployments of lidars from the Danish Technical University (DTU) and German Aerospace Center (DLR) at the Perdigão site. The unwaked flow, measured by the DTU lidar, is well-captured by the lidar, even without combining data into a multi-lidar retrieval. The two DLR lidars measure a waked transect from different downwind vantage points. In the region of the wake, the observation system reacts to the smaller spatial and temporal variations of the winds, allowing more significant observation distortions to arise. While the measurements largely capture the wake structure and trajectory over its 4–5 D extent, limited spatial resolution of measurement points and volume averaging lead to a quicker loss of the two lobes in the near wake, smearing of the vertical bounds of the wake (< 30 m), wake center displacements up to 10 m, and dampening of the maximum velocity deficit by up to a third. The virtual lidar tool, coupled with simulations, provides a means for assessing measurement capabilities in advance of measurement campaigns.

17 WIND ENERGY↗

Sensor Error Estimation for Reactor Coolant System with Generalization Error Equations

We consider a simplified version of the sensor error estimation problem in pressure sensors in the primary coolant system of a nuclear power plant under the condition that a majority of sensors are error-free. We present a sensor fusion method to estimate sensor errors in the primary coolant system of a power reactor by using measurements collected by multiple sensors. The underlying relationship between the measured and estimated parameters is inferred using measurements collected under controlled no-error conditions of the plant. We propose a machine learning approach to train fusers that provide an estimate of the sensor measurement by using measurements from other sensors, and the difference between the actual sensor measurement and its estimate provides the sensor estimate. We present generalization error equations for these methods that characterize the performance of their sensor error estimates on future measurements using the distribution-free machine learning formulation.

Rao, Nageswara S.↗

Method for enhancing a computer to estimate an uncertainty of an onset of a signal of interest in time-series noisy data

A computer-implemented method of enhancing a computer to estimate an uncertainty of an onset of a signal of interest in time-series noisy data. A first mathematical model of first time series data that contains only noise is calculated. A second mathematical model of second time series data that contains the noise and an onset of a signal of interest in the second time series data is calculated. A difference is evaluated between a first combination, being the first mathematical model and the second mathematical model, and a second combination, being the first time series data and the second time series data, wherein evaluating is performed using a generalized entropy metric. A specific time when an onset of the signal of interest occurs is estimated from the difference. An “a posteriori” distribution is derived for an uncertainty of the specific time at which the onset occurs.

Stracuzzi, David John↗

Weak Gravitational Lensing of Low Surface Brightness Galaxies in the DES Year 3 Catalog

We present galaxy-galaxy lensing measurements of a sample of low surface brightness galaxies (LSBGs) drawn from the Dark Energy Survey Year 3 (Y3) data. LSBGs are diffuse galaxies with a surface brightness dimmer than the ambient night sky. Given their faintness, the use of standard observational techniques proves challenging. Weak gravitational lensing probes both the baryonic and dark matter content of galaxies, rendering it a powerful technique to estimate LSBG masses. The LSBG lens sample consists of 23, 790 total extended galaxies separated into red and blue color types at \(g-i\ge 0.60\) and \(g-i\le 0.60\), respectively. We use the Y3 \metacal\ shape catalog as the source sample, with a number density of {\color{red}\(5.59 \mathrm{gal}/\mathrm{arcmin}^2\)}. We measure the tangential shear around the lens galaxies across angular scales of \(0.25-400\) arcmin and find a signal-to-noise of 6.67 for red galaxies, 2.17 for blue galaxies, and 5.30 for the total sample. W e fit a model built from two NFW profiles corresponding to the LSBG dark matter subhalo and host halo to the red LSBG shear measurements with an MCMC. We estimate the host halo mass at \(7.3 ^{+2.0}_{-1.6}\times 10^{12} M_\mathrm{\odot}\). We place a 95\% upper bound on the subhalo mass at \(2.8\times 10^{11} M_\mathrm{\odot}\). We utilize the lens sample photometry to obtain an estimate of the red LSBG stellar mass distribution. We compare the ratio between the stellar mass and the subhalo mass to the parameterized, satellite-specific SHMR \citep{Moster2010}. This work represents the first example of an attempted constraint on the masses of LSBGs using weak gravitational lensing.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of a Loop for Fission Gas Release and Monitoring at the PULSTAR Reactor

Modern nuclear fuel performance analysis tools aim at establishing a predictive approach for examining macroscopic fuel behavior due to fuel gaseous swelling and the associated pellet cladding mechanical interaction (PCMI). In general, the macroscopic behavior is highly sensitive to understanding microscopically driven phenomena, which in this case is the release of fission gas (i.e., Kr and Xe) from the fuel. Previous studies identified the uncertainties associated with the intra-granular fission gas diffusion coefficient as most limiting for the accuracy of fission gas release calculations. In particular, it was concluded that, in order to improve fission gas release as well as gaseous swelling predictions using fuel performance codes, a better characterization of the diffusion coefficient is needed. This work explores the implementation of an on-line fission gas release and measurement loop at the North Carolina State University (NCSU) PULSTAR reactor. A unique characteristic of this loop is the ability to control the temperature of the fuel samples during irradiation to establish a more accurate correlation between the irradiation conditions and an observable such as the intragranular gas diffusion coefficient. Analysis was performed of an 'in-pile' fission gas release experiment that is currently under implementation at the NCSU PULSTAR reactor. Using MCNP6 calculations and COMSOL simulations the fission rate density and temperature distributions in UO{sub 2} microspheres were estimated. To explore the feasibility of detecting fission gas isotopes, estimates of R/B for {sup 138}Xe were obtained using the equivalent sphere model and combined with MCNP6 predicted B estimates to obtain R. Assuming detection in a typical gamma-ray spectrometry HPGe system, it was conclude that reasonable count rates will be achieved to allow extraction of release information. At this stage, accurate assessment is underway of the performance of various loop components to set the final design criteria in support of loop fabrication. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Predicting Biological Cleanliness: An Empirical Bayes Approach for Spacecraft Bioburden Accounting

To comply with the international planetary protection policy set forth by the Committee on Space Research and NASA Agency level requirements, spacecraft destined to biologically sensitive planetary bodies have to minimize terrestrial biological contamination. Analysis, testing and inspection are the standard forward verification activities that are used to demonstrate compliance with the biological contamination requirements. For testing of spacecraft surface areas, a swab or wipe sample is collected from surfaces prior to last access and subsequently processed in the lab using NASA Approved Planetary Protection Methods for Culture Based Assays. Raw data resulting from this assay is then statistically treated employing a mathematical paradigm stemming from the 1970’s Viking Lander Project to generate the bioburden density and total microbial bioburden present. This standard approach arbitrarily accounts for error and provides an upper conservative bound as it reports the maximum number of spores estimated to be present on flight hardware surfaces. A bioburden density estimate factors in the following variables: the observed bioburden count, representative volume processed, sampling efficiencies. Notably, to account for error in the approach, a 0 observed count is arbitrarily changed to a count of 1 for each hardware grouping. The data generated by spacecraft bioburden verification campaigns in the past have resulted in <80% of wipes and <90% of swabs containing a bioburden count of 0. As such, having a robust and well documented statistical approach for dealing with the probability of low incident rates is necessary to be able to estimate spacecraft bioburden. Being able to statistically describe the bioburden distribution and associated confidence level is a gamechanger for the development of bioburden allocations during mission design and will allow for tighter management of risk throughout spacecraft build. Thus, Empirical Bayes statistical approach was evaluated to estimate the microbial bioburden on spacecraft to mitigate the aforementioned mathematical concerns and provide a probabilistic bioburden distribution of the flight hardware surface. For application of this approach to performing bioburden calculations, a range of non-informative prior assumptions on hardware surfaces are explored for Bayesian analyses while informative priors using posterior distributions from prior assays are utilized for Empirical Bayes analyses. Several non-informative priors are currently under investigation to assess fitness including use of these priors to serve as a foundation to build off of NASA specification values or a basis of risk to account for unknowns during the integration and testing process. Informative priors under consideration are generated using sampled bioburden values from hardware originating within like processing environments (e.g. vendor cleaning process or similar assembly process), temporal spacecraft status events as a prediction for hardware cleanliness of future samples, and heritage system bioburden actuals to predict allocation for subsequent missions. Informative priors and probabilistic bioburden distributions are then validated using data sets from the Mars Exploration Rover, Mars Science Laboratory, and InSight missions. Using Empirical Bayes approach to generate a probabilistic bioburden distribution as demonstrated through mission use cases provides a valid approach for use in the end-to-end requirements verification process.

97 MATHEMATICS AND COMPUTING↗

Bayesian Framework for Bioburden Density Estimation in Planetary Protection

To comply with the international planetary protection policy set forth by the Committee on Space Research and NASA Agency level requirements, spacecraft destined to biologically sensitive planetary bodies have to minimize terrestrial biological contamination. Analysis, testing and inspection are the standard forward verification activities that are used to demonstrate compliance with the biological contamination requirements. For testing of spacecraft surface areas, a swab or wipe sample is collected from surfaces prior to last access and subsequently processed in the lab using NASA Approved Planetary Protection Methods for Culture Based Assays. Raw data resulting from this assay is then statistically treated employing a mathematical paradigm stemming from the 1970’s Viking Lander Project to generate the bioburden density and total microbial bioburden present. This standard approach arbitrarily accounts for error and provides an upper conservative bound as it reports the maximum number of spores estimated to be present on flight hardware surfaces. A bioburden density estimate factors in the following variables: the observed bioburden count, representative volume processed, sampling efficiencies. Notably, to account for error in the approach, a 0 observed count is arbitrarily changed to a count of 1 for each hardware grouping. The data generated by spacecraft bioburden verification campaigns in the past have resulted in <80% of wipes and <90% of swabs containing a bioburden count of 0. As such, having a robust and well documented statistical approach for dealing with the probability of low incident rates is necessary to be able to estimate spacecraft bioburden. Being able to statistically describe the bioburden distribution and associated confidence level is a gamechanger for the development of bioburden allocations during mission design and will allow for tighter management of risk throughout spacecraft build. Thus, Empirical Bayes statistical approach was evaluated to estimate the microbial bioburden on spacecraft to mitigate the aforementioned mathematical concerns and provide a probabilistic bioburden distribution of the flight hardware surface. For application of this approach to performing bioburden calculations, a range of non-informative prior assumptions on hardware surfaces are explored for Bayesian analyses while informative priors using posterior distributions from prior assays are utilized for Empirical Bayes analyses. Several non-informative priors are currently under investigation to assess fitness including use of these priors to serve as a foundation to build off of NASA specification values or a basis of risk to account for unknowns during the integration and testing process. Informative priors under consideration are generated using sampled bioburden values from hardware originating within like processing environments (e.g. vendor cleaning process or similar assembly process), temporal spacecraft status events as a prediction for hardware cleanliness of future samples, and heritage system bioburden actuals to predict allocation for subsequent missions. Informative priors and probabilistic bioburden distributions are then validated using data sets from the Mars Exploration Rover, Mars Science Laboratory, and InSight missions. Using Empirical Bayes approach to generate a probabilistic bioburden distribution as demonstrated through mission use cases provides a valid approach for use in the end-to-end requirements verification process.

97 - MATHEMATICS AND COMPUTING↗

Bayesian operator inference for data-driven reduced-order modeling

This work proposes a Bayesian inference method for the reduced-order modeling of time-dependent systems. Informed by the structure of the governing equations, the task of learning a reduced-order model from data is posed as a Bayesian inverse problem with Gaussian prior and likelihood. The resulting posterior distribution characterizes the operators defining the reduced-order model, hence the predictions subsequently issued by the reduced-order model are endowed with uncertainty. The statistical moments of these predictions are estimated via a Monte Carlo sampling of the posterior distribution. Since the reduced models are fast to solve, this sampling is computationally efficient. Furthermore, the proposed Bayesian framework provides a statistical interpretation of the regularization term that is present in the deterministic operator inference problem, and the empirical Bayes approach of maximum marginal likelihood suggests a selection algorithm for the regularization hyperparameters. The proposed method is demonstrated on two examples: the compressible Euler equations with noise-corrupted observations, and a single-injector combustion process.

97 MATHEMATICS AND COMPUTING↗

Net-baryon multiplicity distribution consistent with lattice QCD

We determine a net-baryon multiplicity distribution which reproduces all cumulants measured so far by lattice QCD. We present the dependence on the volume and temperature of this distribution. We find that for temperatures and volumes encountered in heavy-ion reactions, the multiplicity distribution is very close to the Skellam distribution, making the experimental determination of it rather challenging. We further provide estimates for the statistics required to measure cumulants of the net-baryon and net-proton distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

OE0000878-Master thesis Observability Analysis of Power Distribution Systems with Distributed Energy Resources Using Correlational Measurements

In this thesis, an observability analysis framework is proposed for power distribution systems. The framework analyzes the sufficiency of the available measurements for monitoring the system. In the areas of the network where sufficient sensors have been deployed, the complex voltage at the nodes may be estimated. The framework also provides a metric that quantifies the accuracy of these voltage estimates. Due to the cost of sensors and the size of distribution systems, it is frequently the case that the available measurements are insufficient for complete observability of the system. In this thesis, the author proposes the use of Correlational Measurements (CMs) for improving distribution system observability by leveraging correlation between loads and between Distributed Energy Resources (DERs). Load correlation due to class-residential, commercial, industrial, etc.- is captured via correlational load measurements, the first type of CM. Injection correlation of DERs-such as wind turbines located in proximity-is captured via correlational DER measurements, the second type of CM. This thesis presents the CM formulation and derives modified node voltage and branch current based state estimators (NVSE and BCSE) accommodating CMs. Additionally, this thesis contains an on-line parameter estimation scheme for CMs that improves the accuracy of state estimates over time. The observability analysis framework with CMs is demonstrated through simulations on the IEEE 123 Node Distribution Test Feeder. It is shown that the introduction of CMs leads to improved observability of the system with only a 0.7% decrease in state estimation accuracy for NVSE and a 0.5% increase in state estimation accuracy for BCSE.

John, Yohan M↗

Structure- and Physics-Preserving Reductions of Power Grid Models

The large size of multiscale, distribution and transmission, power grids hinders fast systemwide estimation and real-time control and optimization of operations. This paper studies graph reduction methods of power grids that are favorable for fast simulations and follow-up applications. While the classical Kron reduction has been successful in reduced order modeling of power grids with traditional, hierarchical design, the selection of reference nodes for the reduced model in a multiscale, distribution and transmission, network becomes ambiguous. In this work we extend the use of the iterative Kron reduction by utilizing the electric grid's graph topology for the selection of reference nodes, consistent with the design features of multiscale networks. Additionally, we propose further reductions by aggregation of coherent subnetworks of triangular meshes, based on the graph topology and network characteristics, in order to preserve currents and build another power-flow equivalent network. Our reductions are achieved through the use of iterative aggregation of subgraphs that include general tree structures, lines, and triangles. Important features of our reduction algorithms include that (i) the reductions are either equivalent to the Kron reduction or otherwise produce a power-flow equivalent network; (ii) due to the mentioned power-flow equivalence, the reduced network can model the dynamic of the swing equations for a lossless, inductive, steady state network; (iii) the algorithms efficiently utilize hash tables to store the sequential reduction steps. The third feature allows for easy reintroduction of detailed models into the reduced, conceptual network and makes the final reduced order model backward compatible with a sequence of intermediate, partially reduced networks with varying resolution---the ordered sequence of iterative reductions corresponds to a sequence of reduced order models. The performance of our graph reduction algorithms and features of the reduced grids are discussed on a real-word transmission and distribution grid. We produce visualizations of the reduced models through open source libraries and release our reduction algorithms with example code and toy data.

24 POWER TRANSMISSION AND DISTRIBUTION↗

On the Growth and Saturation of the Gyroresonant Streaming Instabilities

The self-regulation of cosmic-ray (CR) transport in the interstellar and intracluster media has long been viewed through the lenses of linear and quasi-linear kinetic plasma physics. Such theories are believed to capture the essence of CR behavior in the presence of self-generated turbulence but cannot describe potentially critical details arising from the nonlinearities of the problem. We utilize the particle-in-cell numerical method to study the time-dependent nonlinear behavior of the gyroresonant streaming instabilities, self-consistently following the combined evolution of particle distributions and self-generated wave spectra in one-dimensional periodic simulations. We demonstrate that the early growth of instability conforms to the predictions from linear physics, but that the late-time behavior can vary depending on the properties of the initial CR distribution. We emphasize that the nonlinear stages of instability depend strongly on the initial anisotropy of CRs—highly anisotropic CR distributions do not efficiently reduce to Alfvénic drift velocities, owing to reduced production of left-handed resonant modes. We derive estimates for the wave amplitudes at saturation and the timescales for nonlinear relaxation of the CR distribution and then demonstrate the applicability of these estimates to our simulations. Bulk flows of the background plasma due to the presence of resonant waves are observed in our simulations, confirming the microphysical basis of CR-driven winds.

79 ASTRONOMY AND ASTROPHYSICS↗