Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional Distribution”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Assessing Photovoltaic Capacity Factor Variability Using Long-Term Satellite Derived Solar Resource Data Under Brazilian Climate

Accurate estimation of photovoltaic (PV) energy yield and its variability is essential for reducing financial risk and supporting reliable system planning for rapidly expanding PV markets. In Brazil, high solar adoption and increasing levels of distributed energy resources are beginning to introduce operational challenges such as curtailment and evolving grid requirements. Understanding how natural variability in solar resource propagates into PV system performance is therefore increasingly important for both project design and grid integration. Modern PV yield assessments commonly rely on multi-year meteorological datasets and probabilistic exceedance metrics (e.g., P50/P90) to quantify energy yield uncertainty for project financing. However, the implications of long-term solar resource variability for PV system design choices and high-adoption grid conditions remain less well characterized for rapidly expanding markets such as Brazil. In particular, understanding how weather-driven variability propagates into PV production distributions and capacity factor expectations is important for evaluating curtailment exposure, deployment strategies, and storage requirements in regions experiencing rapid growth of distributed and utility-scale PV. Seasonal and interannual variability in atmospheric conditions can produce substantial fluctuations in monthly PV energy production, which propagate into uncertainty in annual energy yield and capacity factor expectations. Characterizing this variability using long-term meteorological datasets allows probabilistic estimation of PV system performance and provides improved insight into the range of expected PV energy outcomes. This study explores the use of long-term satellite-derived meteorological data from the National Solar Radiation Database (NSRDB) to evaluate the variability of photovoltaic system performance across multiple locations in Brazil. Using a 27-year dataset (1998-2024), PV system simulations are performed to characterize the distribution of annual and seasonal capacity factors and energy yield outcomes, while propagating key sources of meteorological variability and model uncertainty through the PV modeling chain. The analysis also investigates the sensitivity of PV performance outcomes to key system design assumptions within the PV modeling chain, including tracking configuration and system sizing parameters. The resulting probabilistic performance characterization provides insight into how weather-driven variability influences PV production expectations and capacity factor distributions. These results provide a foundation for evaluating how weather-driven variability interacts with high PV adoption and potential storage or curtailment mitigation strategies.

14 SOLAR ENERGY↗
DOE OSTI2022

A Conditional Autoencoder for Galaxy Photometric Parameter Estimation

Astronomical photometric surveys routinely image billions of galaxies, and traditionally infer the parameters of a parametric model for each galaxy. This approach has served us well, but the computational expense of deriving a full posterior probability distribution function is a challenge for increasingly ambitious surveys. In this paper, we use deep learning methods to characterize galaxy images, training a conditional autoencoder on mock data. The autoencoder can reconstruct and denoise galaxy images via a latent space engineered to include semantically meaningful parameters, such as brightness, location, size, and shape. Our model recovers galaxy fluxes and shapes on mock data with a lower variance than the Hyper Suprime-Cam photometry pipeline, and returns reasonable answers even for inputs outside the range of its training data. When applied to data in the training range, the regression errors on all extracted parameters are nearly unbiased with a variance near the Cramr-Rao bound.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Impact of Time Dependent Reactor and Sensor Physics on Core Power Synthesis

Online synthesis of the power distribution is critical in the operation and control of nuclear power reactors to ensure that the core is operating within safety margins, and to provide essential knowledge associated with the burnup of the fuel. In light water reactors (LWRs), power synthesis is achieved by using some a priori knowledge of the state of the reactor core and updating based on the signals coming from in-core sensors—namely, self-powered neutron detectors (SPNDs). This report aims to study the effects of fuel burnup and sensor degradation on the ability to accurately synthesize the power distribution in a LWR. Several modeling tools were used to simulate power synthesis based on the responses of SPNDs, with emitters made out of Rh or V. A representative pressurized water reactor low-enriched uranium (LEU) core was modeled using the Polaris/Purdue Advanced Reactor Core Simulator (PARCS) approach. The Monte Carlo N-Particle Transport 6 (MCNP6) code was used, as well, to calculate response functions between different segments of fuel to individual SPNDs; this is a crucial parameter for power synthesis. The Oak Ridge Isotope GENeration (ORIGEN) package in the Standardized Computer Analyses for Licensing Evaluation (SCALE) code was used to model the time-dependent isotopic transmutation in the SPND emitters. All these data were fed into a custom code that enacted the point-based iterative (PBI) method to simulate power synthesis. Developmental work was also performed on high-fidelity SPND models in the GEometry ANd Tracking 4 (Geant4) code, which enables higher-accuracy modeling of the current responses from SPNDs. In this work, five sets of time-dependent power synthesis test cases were conducted. In these test cases, systematic changes in the input conditions enabled an analysis of the effect of (1) slightly inaccurate a priori power distribution assumptions with respect to fuel burnup, (2) highly inaccurate a priori assumptions with respect to fuel burnup (such that burnup is not included in the a priori assumed distribution), and (3) differences between Rh and V SPNDs in terms of downstream consequences of the transmutation in the emitters. The authors discovered that one may permissibly have slightly inaccurate a priori assumptions of the fuel burnup (such that the level of burnup may be slightly under- or over-approximated by the accumulated burnup in approximately 9.3 full power days), but to not account for burnup at all in the a priori assumption leads to severe levels of error, approaching 25% at maximum. The authors also discovered that V SPNDs are extraordinarily robust in the low-enriched uranium fuel cycle considered in this modeling work, whereas Rh SPNDs undergo significant transmutation that can result in large errors in the synthesized power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Using Multiple Signatures to Improve Accuracy of Substorm Identification

Abstract We have developed a new procedure for combining lists of substorm onset times from multiple sources. We apply this procedure to observational data and to magnetohydrodynamic (MHD) model output from 1–31 January 2005. We show that this procedure is capable of rejecting false positive identifications and filling data gaps that appear in individual lists. The resulting combined onset lists produce a waiting time distribution that is comparable to previously published results, and superposed epoch analyses of the solar wind driving conditions and magnetospheric response during the resulting onset times are also comparable to previous results. Comparison of the substorm onset list from the MHD model to that obtained from observational data reveals that the MHD model reproduces many of the characteristic features of the observed substorms, in terms of solar wind driving, magnetospheric response, and waiting time distribution. Heidke skill scores show that the MHD model has statistically significant skill in predicting substorm onset times.

Haiducek, John D.↗
DOE OSTI2026

Impact of Time-Dependent Reactor and Sensor Physics on Core Power Synthesis (Rev.1)

Online synthesis of power distribution is critical in the operation and control of nuclear power reactors to ensure that the core is operating within safety margins and to provide essential knowledge associated with the burnup of the fuel. In light-water reactors, power synthesis is achieved by using some a priori knowledge of the state of the reactor core and updating based on the signals coming from in-core sensors—namely, self-powered neutron detectors (SPNDs). This report examines the effects of fuel burnup and sensor degradation on the ability to accurately synthesize the power distribution in a pressurized water reactor (PWR), considering the typical low-enriched uranium (LEU, 3%-5% enrichment) fuel cycle as well as the higher enrichment LEU+ (5%-8% enrichment) fuel cycle. Several modeling tools were used to simulate power synthesis based on the responses of SPNDs, with emitters made out of Rh or V. A representative PWR LEU core was modeled using the Polaris/Purdue Advanced Reactor Core Simulator (PARCS) approach. The Monte Carlo N-Particle Transport 6 (MCNP6) code was used as well to calculate response functions between different segments of fuel to individual SPNDs; this is a crucial parameter for power synthesis. The Oak Ridge Isotope GENeration (ORIGEN) package in the Standardized Computer Analyses for Licensing Evaluation (SCALE) code was used to model the time-dependent isotopic transmutation in the SPND emitters. All these data were fed into a custom code that enacted the point-based iterative method to simulate power synthesis. Developmental work was also performed on high-fidelity SPND models in the GEometry ANd Tracking 4 (Geant4) code, which enables higher-accuracy modeling of the current responses from SPNDs. In this work, five sets of time-dependent power synthesis test cases were conducted. In these test cases, systematic changes in the input conditions enabled an analysis of the effect of (1) slightly inaccurate a priori power distribution assumptions with respect to fuel burnup, (2) highly inaccurate a priori power distribution assumptions with respect to fuel burnup (such that burnup is not included in the a priori assumed distribution), and (3) differences between Rh and V SPNDs in terms of downstream consequences of the transmutation in the emitters and the extended nature of the LEU+ fuel cycle in comparison with LEU. The authors discovered that one may permissibly have slightly inaccurate a priori assumptions of the fuel burnup (such that the level of burnup may be slightly underapproximated or overapproximated by the accumulated burnup in approximately 9.3 full power days), but to not account for burnup at all in the a priori assumptions leads to severe levels of error, approaching 25% at maximum (for LEU). The authors also discovered that V SPNDs are extraordinarily robust in both the LEU and LEU+ fuel cycles considered in this modeling work, whereas Rh SPNDs undergo significant transmutation that can result in large errors in the synthesized power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Possibility to Identify the Contributions from Collapsars, Supernovae, and Neutron Star Mergers from the Evolution of the r-process Mass Abundance Distribution

We study the evolution of rapid neutron-capture process (r-process) isotopes in the galaxy. We analyze relative contributions from core-collapse supernovae (CCSNe), neutron star mergers, and collapsars under a range of astrophysical conditions and nuclear input data. Here we show that, although the r-process in each of these sites can lead to a similar (universal) elemental distribution, the detailed isotopic abundances can differ from one site to another. These differences may allow for the identification of which sources contributed to the early evolution of r-process material in the galaxy. Our simulations suggest that the early evolution was dominated by CCSNe and collapsar r-process nucleosynthesis. This conclusion may be testable if the next generation of observatories can deduce isotopic r-process abundances.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The onset distribution of rotating m , n = 2 , 1 tearing modes and its consequences on the stability of high-confinement-mode plasmas in DIII-D

Abstract Analysis of a multi-scenario database of over 13 000 DIII-D H-mode discharges shows that the m , n = 2 , 1 magnetic islands are dominantly pressure gradient driven, stochastically triggered non-linear instabilities at all edge safety factor ( q 95 ) values. The instability onset time closely follows the exponential distribution in intermediate and high q 95 scenarios and is characterized by near constant onset rate ( λ ), in accordance with Poisson-point processes. This implies that the plasmas are operated in marginally stable conditions, characterized by a small threshold for instability growth and variations in the trigger amplitude and/or the stabilizing mechanisms with temporally uniform random distribution in this database. While the majority of the tearing modes occur in the first current-profile relaxation time of the β N flattop, constant λ throughout the β N flattop shows that the tearing onset is insensitive to the evolution of the equilibrium current profile. In low q 95 scenarios, where a large fraction of the plasmas are operated at low torque, λ increases over the course of the β N flattop, showing that these plasmas evolve toward more unstable conditions. The onset rate rapidly increases with β N , while it does not show a clear dependence on the current gradient at the mode rational surface. Overall, these observations support that the majority of the analyzed 2,1 tearing modes are non-linear, neoclassically driven instabilities and classical stability does not play a dominant role in their onset.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Boundary condition effects on runaway electron mitigation coil modeling for the SPARC and DIII-D tokamaks

Abstract Extended-MHD modeling of planned Runaway Electron Mitigation Coils (REMC) for SPARC and DIII-D is performed with the NIMROD code. A coil has been designed for each machine, with the two differing in shape and location, but both having n = 1 symmetry (with n the toroidal mode number). Compared to previous modeling efforts, three improvements are made to the simulations boundary conditions. First a resistive wall model is used in place of an ideal wall. Second, the ThinCurr code is used to compute the time-dependent 3D fields used as magnetic boundary conditions for the simulations. Third, the simulation boundary is moved from the first-wall location to the Vacuum Vessel (VV), which extends the boundary past the location of the internal REMC. To remove the 3D coil from the simulation domain, an equivalent set of 3D fields is calculated at the VV boundary that produce approximately the same field distribution at the last closed flux surface assuming vacuum between the two. Each of these three boundary condition improvements leads to an improvement in the predicted performance of the REMC for both machines. The resistive wall alone primarily effects the resonance of the coil with the plasma after the TQ, affecting the q-profile evolution in the SPARC modeling, and allowing the applied spectrum to be modified in response to the plasma in the DIII-D modeling. The movement of the simulation boundary has the most significant effect on the RE confinement overall, including in the early stages, particularly for a DIII-D inner wall limited equilibrium, where the RE loss fraction increases from 90% to > 99%, with SPARC RE losses also occurring much earlier when the boundary is placed at the VV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Ultrasonic Measurements of Temperature Profile and Heat Fluxes in Coal-Fired Power Plants (Final Report)

Many industrial processes are inaccessible or inhospitable to characterization by traditional temperature measurement methods, such as thermocouples, especially over prolonged exposure to harsh environments. Ultrasound is an established characterization technology with diverse applications ranging from medical imaging to therapies to flaw detection to nondestructive evaluation. Ultrasound may characterize solid materials and components noninvasively as a nondestructive evaluation modality and obtain internal measurements of material properties. For example, the speed of ultrasound propagation changes with Young’s modulus and Poisson’s ratio, which can be found from its measurements. Traditional ultrasonic characterization assumes all material properties remain constant with the position. When this assumption holds, a property of interest may be measured by relating it to the speed of ultrasound propagation (or a speed of sound, SOS) and measuring the SOS by timing the ultrasound propagation through a known distance. However, when a property of interest is spatially distributed, the propagation time depends on the SOS changing with the position along the ultrasound propagation path. The multiple temperature distributions may lead to an identical time of flight (TOF). Temperature is one property that impacts the speed of ultrasound and often cannot be assumed to remain constant with the position. Previously, in the context of temperature, we addressed the challenge of ultrasonic characterization of spatially distributed properties by developing a method for measuring segmental temperature distributions (MSTD). This method divides the ultrasonic propagation into segments bound by echogenic features. These features provide ultrasonic interfaces where some energy is reflected toward the receiving transducer, and the rest continues through the medium. The time-of-flight between the echoes reflected from echogenic features characterizes the spatial distribution in the properties of interest in the corresponding segment of the ultrasonic propagation path. This project demonstrated the application of the MSTD method in industrial conditions of the coal-fired power plant. We implemented the MSTD using metals and alloys waveguides, which may be the existing structure for which the temperature distribution is characterized or purposefully designed waveguides added to the structure by welding or other means specifically to quantify thermal properties using the MSTD method. Previous iterations of the MSTD method used ceramic and cementitious waveguides, which significantly attenuate ultrasound. On the other hand, low attenuation in metallic waveguides creates interactions between echogenic features which compilates the signal analysis in the segmental TOF measurements. We have established the WG design principles that minimize the interferences between trailing and primary echoes and, in some cases, eliminate them. The waveguides in which echoes do not interfere improve the timing accuracy and the robustness of ultrasonic measurements of the spatial distributions in material properties. Our emphasis remained on the estimation of the temperature distributions. We have developed general recommendations for designing ultrasonically segmented waveguides with the reduced influence of trailing echoes. Two of our waveguide designs were tested in the industry. The first waveguide was designed for insertion into a combustion zone of the utility-scale coal-fired power plant boiler. The second design allows the characterization of temperature distribution in the direction normal to the boiler’s water wall, a large heat exchanger converting the chemical energy released during combustion to the steam driving the electrical power generation turbines. These waveguides were designed to operate within a restrictive space of thermally insulated water wall and incorporate densely located echogenic features while combatting the influence of trailing echoes. The project has successfully demonstrated the feasibility of using the developed method for accurate, continuous, and robust temperature measurements in extreme environments of power generation and other industrial processes. It, therefore, has achieved its overarching goal of advancing the technology readiness level of the novel Ultrasound Measurements of Segmental Temperature Distribution (US-MSTD) method for real-time measurements of the temperature distribution and heat fluxes closer to commercial availability, developing a prototype multipoint measurement system, and validating its performance on coal-fired utility boilers. The success of this project was achieved in collaboration with the power generator, Rocky Mountain Power, and set the stage for the transfer of this technology from the laboratory to the industry.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Ribosome-inactivation by a class of widely distributed C-tail anchored membrane proteins

Ribosome hibernation is a commonly used strategy that protects ribosomes under unfavorable conditions and regulates developmental processes. Multiple ribosome-hibernation factors have been identified in all domains of life, but due to their structural diversity and the lack of a common inactivation mechanism, it is currently unknown how many different hibernation factors exist. Here, we show that the YqjD/ElaB/YgaM paralogs, initially discovered as membrane-bound ribosome binding proteins in E. coli, constitute an abundant class of ribosome-hibernating proteins, which are conserved across all proteobacteria and some other bacterial phyla. Our data demonstrate that they inhibit in vitro protein synthesis by interacting with the 50S ribosomal subunit. In vivo cross-linking combined with mass spectrometry revealed their specific interactions with proteins surrounding the ribosomal tunnel exit and even their penetration into the ribosomal tunnel. Thus, YqjD/ElaB/YgaM inhibit translation by blocking the ribosomal tunnel and thus mimic the activity of antimicrobial peptides and macrolide antibiotics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Environmental control on the distribution of metabolic strategies of benthic microbial mats in Lake Fryxell, Antarctica

Ecological theories posit that heterogeneity in environmental conditions greatly affects community structure and function. However, the degree to which ecological theory developed using plant- and animal-dominated systems applies to microbiomes is unclear. Investigating the metabolic strategies found in microbiomes are particularly informative for testing the universality of ecological theories because microorganisms have far wider metabolic capacity than plants and animals. We used metagenomic analyses to explore the relationships between the energy and physicochemical gradients in Lake Fryxell and the metabolic capacity of its benthic microbiome. Statistical analysis of the relative abundance of metabolic marker genes and gene family diversity shows that oxygenic photosynthesis, carbon fixation, and flavin-based electron bifurcation differentiate mats growing in different environmental conditions. The pattern of gene family diversity points to the likely importance of temporal environmental heterogeneity in addition to resource gradients. Overall, we found that the environmental heterogeneity of photosynthetically active radiation (PAR) and oxygen concentration ([O2]) in Lake Fryxell provide the framework by which metabolic diversity and composition of the community is structured, in accordance with its phylogenetic structure. The organization of the resulting microbial ecosystems are consistent with the maximum power principle and the species sorting model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Nonequilibrium spherulitic magnetite in the Ryugu samples

Here, we have investigated several particles collected during each of two touchdowns of the Hayabusa2 spacecraft at the surface of the C-type asteroid 162173 Ryugu using various electron microscope techniques. Our detailed transmission electron microscopy study shows the presence of magnetite with various morphologies coexisting in close proximity. This is characteristic of CI chondrite-like materials and consistent with the mineral assemblages and compositions in the Ryugu parent body. We describe the microstructural characteristics of magnetite with different morphologies, which could have resulted from the chemical conditions (growth vs. diffusion rate) during their formation. Furthermore, we describe the presence of magnetites with a spherulitic structure composed of individual radiating fibers that are characterized by pervasive, homogeneously distributed euhedral to subhedral pores that have not been described in previous chondrite studies. This particular spherulitic structure is consistent with crystallization under nonequilibrium conditions. Additionally, the presence of a high density of defects within the magnetite fibers, the high surface/volume ratio of this morphology, and the presence of amorphous materials in several pores and at the edges of the acicular fibers further support their formation under nonequilibrium conditions. We suggest that the growth processes that lead to this structure result from the solution reaching a supersaturated state, resulting in an adjustment to a lower free energy condition via nucleation and rapid growth.

58 GEOSCIENCES↗
DOE OSTI2024

Loss cone effects and monotonic sheath conditions of a partially magnetized plasma sheath

In this Letter, we propose the conditions for monotonic plasma sheaths adjacent to a floating wall in the presence of an applied, oblique magnetic field. The electron velocity distribution function (VDF) at the sheath edge obtained from a kinetic model exhibits a loss cone shaped truncation. Using an approximation of the truncated VDF, we derive an analytical framework of the sheath edge condition (i.e., Bohm condition), namely, the relation of ion injection velocity and sheath potential drop as a function of the magnetic field angle. The results show that the sheath edge velocity and total potential drop decrease for a steady-state sheath, which eventually collapses when the magnetic field lines become parallel to the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A Poisson equation method for prescribing fully developed non-Newtonian inlet conditions for computational fluid dynamics simulations in models of arbitrary cross-section

Prescribing inlet boundary conditions for computational fluid dynamics (CFD) simulations of internal flow in complex geometries such as anatomical vascular models is challenging. In the absence of patient-specific inlet velocity data, a common approach for long blood vessels is to assume that the inlet flow is fully developed. In vessels of irregular cross section, however, prescribing fully developed conditions is complicated due to the lack of a general closed-form analytical solution. In this study, we develop a simple Poisson equation method for prescribing fully developed inlet conditions for the flow of either Newtonian or non-Newtonian fluids in CFD models of arbitrary cross-section. We first derive the generalized Poisson equation for fully developed flow of a non-Newtonian fluid and we then develop and verify a methodology for numerically computing the solution on any planar boundary domain. In addition, we develop a simple extension of the method for prescribing a non-orthogonal inlet velocity that represents fully developed flow from an upstream tube that is connected to the CFD inlet at a non-orthogonal angle. This may be used to investigate a common source of uncertainty in CFD simulations of internal flow that is due to a lack of information concerning the exact streamwise flow direction at the inlets. Comparison to several Newtonian and non-Newtonian benchmark verification solutions shows the method to be extremely accurate. As a practical demonstration case, we use the method to prescribe fully developed conditions on multiple non-circular inlets for the non-Newtonian flow of blood in a patient-specific model of the inferior vena cava (IVC). Finally, we further demonstrate the utility of the method by performing a sensitivity study using the patient-specific IVC model, wherein we investigate the influence of inlet velocity flow direction on the non-Newtonian IVC hemodynamics. Given its simplicity and computational efficiency, the method is shown to be far superior to alternative approaches for prescribing fully developed inlet conditions in such complicated geometries. In conclusion, to facilitate the adoption of our Poisson equation method, we have distributed our OpenFOAM source code and the associated test cases from this study as open-source software.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Drop clustering and drop size correlations from holographic imagery suggest cloud droplet spectral broadening via entrainment-mixing

The question of how droplets rapidly grow large enough to initiate collision-coalescence has persisted for decades. Many theories explaining the production of sufficiently large drops (i.e., those in the “bottleneck” size range; ∼ 25–50 µm diameters) involve drop clustering on millimeter scales. A novel method is introduced to evaluate drop clustering trends particle-by-particle (i.e., the number/proximity of neighboring drops for given droplets; defined as drop clustering fields) which are diagnosed relative to drops within their shared drop environments – in contrast to previous studies which diagnose drop clustering of defined sample volumes, or in terms of absolute length scales. Specifically, this study evaluates the statistical likelihood that drops of a given size are associated with either a significant number of neighboring drops, or are significantly isolated from neighboring drops. Observations are acquired from the HOLODEC during the Cloud System Evolution in the Trades campaign, which sampled subtropical marine clouds. The HOLODEC measures drop size distributions and the 3D spatial coordinates of droplets. Results show drops within the bottleneck size range (diameters of ∼ 25–50 µm) are most likely to be significantly isolated from neighboring drops. This “isolated large drop trend” is primarily observed at subsaturated conditions, suggesting entrainment is the contributing factor. Holograms associated with this trend are more likely to have broader drop size distributions, larger maximum drop sizes and overly regions where precipitation reaches the lowest altitudes from the sampled cloud, suggesting entrainment-mixing drop size distribution broadening is a relevant precipitation-initiation mechanism.

D'Alessandro, John J. [Univ. of Washington, Seattl↗
DOE OSTI2020

Distribution System Segmented Model Simplification with Independent Dynamically Changing End-Use Loads

With more emerging technologies in power distribution systems, to support quasi-static time-series (QSTS) simulation over short time frames, accurately modeling the dynamic behavior of end-use load under various conditions is im- portant. However, time and cost can be a burden for the QSTS simulation of thousands of houses representing dynamic load behaviors on large distribution feeders. There are two main contributions of this paper in order to reduce the QSTS simulation complexity. First, the paper presents a method to stochasti- cally model the dynamic behavior of end-use load for QSTS simulations. Second, using the stochastically modeled load behavior as an input to a distribution sys- tem, and building on our previous QSTS simulation simplification/acceleration method that relies on uniform load behavior, this paper presents a segment substitution approach for the simplification of distribution system models with independent dynamically changing end-use loads. Our numerical results show that compared to GridLAB-D, a software tool that simulates end-use load be- havior through physical models, the proposed stochastic modeling approach achieved an 89.3% simulation time reduction for end-use loads, while accurately capturing the temporally changing total load behavior. Moreover, the QSTS simulation of a large test distribution system through OpenDSS showed that the proposed simplification method decreased the feeder simulation time by a factor of three at less than 0.3% voltage error. Compared to our previous method with independent dynamic loads, the voltage error is reduced by a factor of four.

Power distribution system modeling, stochastic pro↗
DOE OSTI2026

Development of physics-consistent conditional diffusion model to overcome data scarcity in critical heat flux

Deep generative modeling provides a powerful pathway to overcome data scarcity in energy-related applications where experimental data are often limited. By learning the underlying probability distribution of the training dataset, deep generative models, such as the diffusion model, can generate high-fidelity synthetic samples that statistically resemble the training data. Such synthetic data generation can significantly enrich the size and diversity of the available training data, and more importantly, improve the robustness of downstream machine learning models in predictive tasks. The objective of this paper is to investigate the effectiveness of diffusion models for overcoming data scarcity in nuclear energy applications. By leveraging a public dataset on critical heat flux which covers a wide range of commercial nuclear reactor operational conditions, we developed a diffusion model that can generate an arbitrary amount of synthetic samples. Since a vanilla diffusion model can only generate samples randomly, we also developed a conditional diffusion model capable of generating targeted critical heat flux data under user-specified thermal-hydraulic conditions. The performance of the diffusion model was evaluated based on its ability to capture empirical feature distributions and pair-wise correlations, as well as to maintain physical consistency. The results showed that both the diffusion model and conditional diffusion model can successfully generate realistic and physics-consistent critical heat flux data. Furthermore, uncertainty quantification results demonstrate that the conditional diffusion model is highly effective in augmenting critical heat flux data while maintaining acceptable levels of uncertainty.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

BI-LEVEL OPTIMIZATION FOR ELECTRICITY TRANSACTION IN SMART COMMUNITY WITH MODULAR PUMP HYDRO STORAGE

Grid integration of the increasing distributed energy resources could be challenging in terms of new infrastructure investment, power grid stability, etc. To resolve more renewables locally and reduce the need for extensive electricity transmission, a community energy transaction market is assumed with market operator as the leader whose responsibility is to generate local energy prices and clear the energy transaction payment among the prosumers (followers). The leader and multi-followers have competitive objectives of revenue maximization and operational cost minimization. This non-cooperative leader-follower (Stackelberg) game is formulated using a bi-level optimization framework, where a novel modular pump hydro storage technology (GLIDES system) is set as an upper level market operator, and the lower level prosumers are nearby commercial buildings. The best responses of the lower level model could be derived by necessary optimality conditions, and thus the bi-level model could be transformed into single level optimization model via replacing the lower level model by its Karush-Kuhn-Tucker (KKT) necessary conditions. Several experiments have been designed to compare the local energy transaction behavior and profit distribution with the different demand response levels and different local price structures. The experimental results indicate that the lower level prosumers could benefit the most when local buying and selling prices are equal, while maximum revenue potential for the upper level agent could be reached with non-equal trading prices.

Chen, Yang↗