Search NASA⌕ Search

SEARCH · Search NASA

Results for “multipliers”

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 OSTI2022

Multi-site assessment of the potential of fine resolution red-edge vegetation indices for estimating gross primary production

Gross primary production (GPP) models driven by fine resolution remote sensing data characterize the spatial and temporal heterogeneities in plant photosynthesis, which is largely dependent on biome-specific maximum photosynthetic capacity. The red-edge reflectance, sensitive to leaf chlorophyll content, is a good proxy of maximum photosynthetic capacity. More importantly, studies show that the red-edge reflectance-related chlorophyll content index (CIr) multiplied by the incident photosynthetic active radiation (PAR in ) strongly correlates to GPP estimated at carbon flux towers (GPP flux ). Yet, to the best of our knowledge, there is no systematic study investigating the general relationship between fine spatial resolution CIr and GPP among biomes and the relationship between CIr and maximum photosynthetic capacity in GPP models. To provide an overview on incorporating space-borne CIr into a GPP model, we applied fine resolution Sentinel-2-derived CIr and GPP flux over 57 flux sites representative of 10 biomes. We investigated the relationship between CIr and GPP flux , and the spatio-temporal relationship between CIr and ecosystem maximum photosynthetic capacity indicated by the potential ecosystem light use efficiency (LUE pot ). We also evaluated the relationship between other five vegetation indices (VIs) and GPP flux . Results showed that the CIr multiplied by PARin has a higher agreement (R 2 > 0.5) with GPP than other VIs. A universal relationship exists between the CIr multiplied by PARin and GPP, except for forest biomes. The CIr also strongly (R 2 > 0.5) relates to the LUE pot during the peak of the growing season. The CIr has a low spatial variance (CV = 0.25) among biomes, highlighting that CIr can be a proxy of maximum photosynthetic capacity in GPP models that do not require biome-dependent coefficients. This study provides insight for incorporating CIr into GPP models and better quantifying global terrestrial photosynthesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Stress-Induced Structural Transformations in Au Nanocrystals

Nanocrystals can exist in multiply twinned structures like icosahedron or single crystalline structures like cuboctahedron. Transformations between these structures can proceed through diffusion or displacive motion. Experimental studies on nanocrystal structural transformations have focused on high-temperature diffusion-mediated processes. Limited experimental evidence of displacive motion exists. We report structural transformation of 6 nm Au nanocrystals under nonhydrostatic pressure of 7.7 GPa in a diamond anvil cell that is driven by displacive motion. X-ray diffraction and transmission electron microscopy were used to detect the structural transformation from multiply twinned to single crystalline. Single crystalline nanocrystals were recovered after unloading, then quickly reverted to the multiply twinned state after dispersion in toluene. Furthermore, the dynamics of recovery was captured using TEM which showed surface recrystallization and rapid twin boundary motion. Molecular dynamics simulations showed that twin boundaries are unstable due to defects nucleated from the interior of the nanocrystal.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

W-band system-on-chip electron cyclotron emission imaging system on DIII-D

Monolithic, millimeter wave “system-on-chip” (SoC) technology has been employed in heterodyne receiver integrated circuit radiometers in a newly developed Electron Cyclotron Emission Imaging (ECEI) system on the DIII-D tokamak for 2D electron temperature profile and fluctuation evolution diagnostics. Here, a prototype module operating in E-band (72-80 GHz) was first employed in a 2 x 10 element array that demonstrated significant improvements over the previous quasi-optical Schottky diode mixer arrays during the 2018 operational campaign of the DIII-D tokamak. For compatibility with International Thermonuclear Experimental Reactor (ITER) relevant scenarios on DIII-D, the SoC ECEI system was upgraded with 20 horn-waveguide receiver modules. Each individual module contains a University of California at Davis designed W-band (75 -110 GHz) receiver die that integrates a broadband low noise amplifier (LNA), a double balanced down-converting mixer, and a x4 multiplier on the local oscillator (LO) chain. A x2 multiplier and two IF amplifiers are packaged and selected to further boost the signal strength as well as downconvert the signal frequency. The upgraded W-band array exhibits > 30 dB additional gain and 20x improvement in noise temperature comparing with the previous Schottky diode radio frequency (RF) mixer input systems; an internal 8 times multiplier chain is used to bring down the LO frequency below 12 GHz, thereby obviating the need for a large aperture for quasi-optical LO coupling and replacing it with coaxial connectors. The horn-waveguide shielding housing avoids out-of-band noise interference on each individual module. The upgraded ECEI system plays an important role for absolute electron temperature evolution and fluctuation measurements for edge and core region transport physics studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Sensitivity Analysis of Transient Critical Heat Flux by RIA Under High-Pressure Flow Boiling Conditions in TRTL

A reactivity-initiated accident (RIA) is a design-basis accident under which critical heat flux (CHF) is likely to be exceeded. The operational margin for RIAs is currently determined using steady-state CHF lookup tables, which provide conservative estimates relative to transient CHF phenomena. The Transient Reactor Test Loop (TRTL) facility at Oregon State University is capable of performing out-of-pile rapid heating experiments representative of a RIA at conditions representative of a pressurized water reactor (PWR). Here, to further our understanding of and ability to predict transient CHF under PWR conditions, we performed a sensitivity analysis on a RELAP5-3D model of the TRTL facility coupled to the RAVEN code framework to define a proposed experimental test matrix to be performed at the TRTL facility. We then implemented a flow boiling CHF correlation into RELAP5-3D and performed a secondary sensitivity analysis inspecting the impact of the built-in RELAP5-3D CHF and heat transfer multipliers on both the prediction of CHF and key safety parameters, such as peak cladding temperature and heat flux. The results show that the multiplier with the highest influence toward the prediction of CHF occurrence and the safety parameters is the transient CHF multiplier. Operational performance envelopes have been developed for each of the test matrix cases and will be used for validation once the experiments are performed. The TRTL facility is currently performing shakedown testing to verify system performance prior to proceeding with the experimental campaign. Restart testing results include pump curve restart testing, pressure tests, and heater rod thermocouple transients.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Reference Cell Performance and Modeling on a One-Axis Tracking Surface: Preprint

Performance of five silicon-based reference cells is examined on a single-axis tracking surface. The reference cells' output are modeled using one-minute sampled spectral irradiance and reference cell temperature measurements. The model also incorporates spectral responsivity data for the reference cell. The transmission of light through the glazing multiplies the sum over all appropriate wavelengths of temperature adjusted reference cell responsivity times the measured spectral irradiance. Modeled reference cell output is compared with measured reference cell output under clear sky and totally cloudy sky conditions. For each reference cell the ratio of modeled to measured output varies by less than 2% over the year.

modeling↗
DOE OSTI2021

Ice Multiplication by Fragmentation during Quasi-Spherical Freezing of Raindrops: A Theoretical Investigation

Ice multiplication by fragmentation during collision–freezing of supercooled rain or drizzle is investigated. A zero-dimensional dynamical system describes the time evolution of number densities of supercooled drops and ice crystals in a mixed-phase cloud. The characteristic time scale for this collision–freezing ice fragmentation is controlled by the collision efficiency, the number of ice fragments per freezing event, and the available number concentration of supercooled drops. The rate of the process is proportional to the number of supercooled drops available. Thus, ice may multiply extensively, even when the fragmentation number per freezing event is relatively small. The ratio of total numbers of ice particles to those from the first ice, namely, the “ice-enhancement factor,” is controlled both by the number of fragments per freezing event and by the available number concentration of supercooled drops in a similar manner. Especially, when ice fragmentation by freezing of supercooled drops is considered in isolation, the number of originally existing supercooled drops multiplied by the fragmentation number per freezing event yields the eventual number of ice crystals. When supercooled drops are continuously generated by coalescence, ice crystals from freezing fragmentation also continuously increase asymptotically at a rate equal to the generation rate of supercooled drops multiplied by the fragmentation number per freezing event. All these results are expressed by simple analytical forms, thanks to the simplicity of the theoretical model. These parameters can practically be used as a means for characterizing observed mixed-phase clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Failure Analysis–Informed Risk Assessment Framework for Geological Carbon Storage Using Numerical Simulation and Machine Learning

Geological carbon storage (GCS) is recognized as a critical technology for achieving large-scale reductions in anthropogenic carbon dioxide (CO 2 ) emissions. Ensuring long-term containment and safety requires robust risk assessment frameworks that account for geological uncertainty and identify potential failure scenarios. Among various indicators, the area of review (AoR) serves as a key metric for evaluating storage performance, regulatory compliance, and monitoring design, as it delineates the spatial extent impacted by pressure buildup and plume migration. However, conventional AoR-based risk assessments typically perturb parameters within narrow uncertainty bounds, potentially overlooking rare but high-impact events arising from extreme geological conditions. In this study, we present a failure analysis–informed risk assessment framework for large-scale GCS projects to improve site prescreening and monitoring design. A suite of 300 numerical simulations was generated using stochastic geological models that vary five key parameters: net-to-gross ratio, anisotropy azimuth, porosity multiplier, permeability multiplier, and vertical-to-horizontal permeability ratio. Among these, 200 realizations represent normal geological uncertainty, while 100 additional cases explore extreme yet plausible conditions for failure-case analysis. The AoR was simulated and computed from pressure and CO 2 saturation fields, where the baseline AoR boundary, representing the extent predicted under typical geological uncertainty, was defined as the union of 200 normal-range simulations, and failure was identified when extreme-range cases exceeded this baseline. Results show that incorporating broader parameter uncertainty produces significantly larger AoR extents, underscoring the potential underestimation of risk under conventional uncertainty ranges. Furthermore, spatial probability maps derived from failure-induced AoR exceedance identify regions requiring enhanced monitoring attention. Various machine learning (ML)–based classifiers were developed to predict failure occurrence from geological parameters, with the random forest model achieving the highest performance (F1-score of 0.986). Consistent findings from correlation coefficient, feature importance, and Sobol sensitivity analyses reveal that low net-to-gross ratios and permeability multipliers are the dominant risk drivers, reflecting reduced reservoir connectivity and limited pressure dissipation. Altogether, these results provide a novel framework for risk-informed site prescreening and monitoring design that explicitly considers rare but high-impact geological scenarios in GCS projects.

25 ENERGY STORAGE↗
DOE OSTI2023

A New Plasma Radar Concept for Simultaneous Magnetic and Density Measurements

An innovative, compact 288GHz interferometer has been fabricated, tested, installed and successfully demonstrated on the LAPD-U magnetized plasma at UCLA. The system takes advantage of frequency modulated (FM) radar techniques to deliver a compact heterodyne system. In addition, the reflected power from the source is taken advantage of to eliminate the need for additional quasi-optical components. Electron density in LAPD-U plasma has recently been increased substantially thereby requiring a higher frequency/shorter wavelength interferometer to avoid deleterious refractive effects. This system satisfies those needs. The system uses a 96GHz varactor tuned Gunn oscillator which passes to a passive tripler. This tripler has ~3% conversion efficiency. The 288GHz radiation is then coupled to free space using a so-called dual-mode or Pickett horn. The output 288GHz beam is then coupled to an aspheric lens manufactured from low-loss, high-density polyethylene. This lens is employed to collimate the emerging beam. Small axial adjustment of the lens position can also be used to create a slowly focusing beam so as to optimize the measured signal. In addition, up-down or side-to-side adjustment of the lens can be utilized to steer the beam vertically or horizontally – again to optimize alignment. The propagating beam passes through a beam splitter and then through a water-free, bubble-free fused quartz window into the LAPD-U vacuum vessel. The beam-splitter is a thin sheet of G10 which reflects a small fraction of the incident power (~5 %) towards a zero-bias detector optimized for the frequency range from 220 to 300GHz. Note that waveguides at this frequency have dimensions of ~0.9mm x 0.45mm and so have very large conductive losses. This drives the use of quasi-optical propagation. The detector requires no DC bias and is very responsive (> 1V/mW into 1MΩ). Radiation is coupled to the detector via a similar lens-horn arrangement used for the launch. This reflected beam acts as the local oscillator or reference millimeter-wave beam for the detector. The remainder of the launched source beam then enters the LAPD-U vacuum vessel and passes through the plasma at the mid-plane until reaching the opposing port which is closed off with an aluminum flange. This flange is used as a mirror to retroreflect the incident 288GHz beam back along its path. The retroreflected beam exits the input port but does NOT couple directly into the zero-bias detector. Instead, the majority of the return power continues towards the 288GHz source. As mentioned above the transmitted beam enters the source a second time. This would appear undesirable. However, at these frequencies multipliers are highly non-linear elements which results in a significant portion of the return beam (~20%) re-emerging from the multiplier and horn and then coupling via the G10 beam-splitter to the zero-bias detector. This approach eliminated the need for a second quasi-optical beam-splitter. The system is extremely compact measuring approximately 28 inches x 20 inches. The above did not explain how heterodyne operation was achieved. As mentioned above the Gunn oscillator is able to be varactor tuned. This allows a low voltage to be applied to control the operating frequency of the Gunn oscillator. During heterodyne operation a sawtooth shaped voltage is applied to the varactor at 750kHz using an 80MHz Arbitrary Waveform Generator (AWG). This voltage changes the Gunn frequency linearly during the up-sweep which is then reset abruptly at the sawtooth crash to be immediately followed by another linear sweep. Passage through the 288GHz multiplier triples the frequency change experienced by the electromagnetic wave. These frequency changes are small – tens of megahertz. This FM radar approach results in the launched electromagnetic wave frequencies at the detector for the reference and plasma wave to be different. The approximately 10 ns delay propagation delay for the plasma beam results in the local oscillator and plasma beams NOT having an identical frequency – there is in fact a fixed difference frequency. The frequency tuning level of the Gunn oscillator is then adjusted so that there is ONE cycle of this difference frequency during each linear ramp. During the sawtooth crash or downward re-sweep this one cycle replays in reverse but on a very fast timescale. The process then repeats. Low-pass filtering eliminates the fast re-sweep to leave a pure sine wave heterodyne signal. When the plasma is present it introduces a phase delay in the sine wave (caused by the extremely small Doppler shift resulting from the optical path length change). Of course, to measure this phase change we need a reference. This is simply obtained from an arbitrary waveform generator which provides a synchronized output pulse train which again is low pass filtered to obtain a 750kHz sinusoidal voltage reference for the interferometer. The interferometer was installed on LAPD-U where it has worked reliably and has established that electron densities exceeding 1x10 13 cm -3 are routinely achieved. In addition, the system sensitivity was able to easily observe density fluctuation at frequencies up to 50kHz. FM Radar techniques have enabled a full demonstration of a compact, sensitive, high frequency (288GHz/1mm) heterodyne interferometer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Analysis of Lobe Power Calculator and Indication System with Physics and Cycle Based Models

The Advanced Test Reactor (ATR) at INL measures reactor power through various methods, two of them being thermal and Nitrogen-16 (N-16) activity. Water power calculator (WPC) is a thermal power system that measures flow and temperature to determine thermal quadrant powers. The N-16 system utilizes a beta detector that outputs Nitrogen activation levels to calculate lobe power through an algorithm called Lobe Power Calculation and Indication system (LPCIS). The LPCIS utilizes the N-16 system and the WPC system to determine reactor core power levels. The WPC provides accurate calculations of quadrant and total reactor thermal power. With the use of WPC measurements, thermal-to-N-16 (T2N) power ratios are produced to determine if the two indication systems agree on core power. Relative magnitude equations are used to utilize N-16 coefficients and multipliers to improve the indications of the LPCIS. These correct indications are crucial for maintaining safety limits because operators rely on this information for decision making. Currently, the LPCIS system uses linear equations and matrices to calculate lobe power through multipliers and coefficients. Advancements in technology and system upgrades have increased the accuracy of power readings by making the system more dynamic. The new proposed coefficient and multiplier method implements a cycle specific and physics-based model intended for changing coefficients during operation. This calibration experiment focused on power splitting, outer shim and neck shim, as well as fuel burning into the reactor digital acquisition system (RDAS) weighting factors. Results demonstrated that the physics learning method yields a smaller error margin inside of the desired power range for the data set 166-A. Continuing to improve the physics-based model will help improve the power accuracy of the LPCIS system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

HarDWR - Cumulative Water Rights Curves

A dataset within the Harmonized Database of Western U.S. Water Rights (HarDWR). For a detailed description of the database, please see the meta-record v2.0. Changelog v2.0 - Recalculated based on Harmonized Water Rights Records v2.0 sourced from WestDAAT - Added "Unspecified" was a water source category, and files associated with this category v1.01 - Updated the names of each file with an ID number less than 3 digits to include leading 0s v1.0 - Initial public release Description This product an updated version of the database used as input to the WBM model (Grogan et al. in review; Grogan et al. 2022), and is the result of the step creating cumulative water rights curves described in Lisk et al. (2024). This database contains 2,667 individual .csv files, three for each Water Management Area (WMA) in the 11-state region. File Naming WMA_[###]_[X]W.csv, where [###] is the unique identifier for each WMA, and [X] is either "S" for surface water rights, "G" for groundwater rights, or "U" for unspecified. Column headers in each file: Year: the priority date year CUML: the total cumulative water rights allocated up to this priority date year (ft3s-1) Irrigation: The percent of cumulative water rights allocated to the Irrigation category up to this priority date year (%) Domestic: The percent of cumulative water rights allocated to the Domestic category up to this priority date year (%) Livestock: The percent of cumulative water rights allocated to the Livestock category up to this priority date year (%) Fish: The percent of cumulative water rights allocated to the Fish category up to this priority date year (%) Industrial: The percent of cumulative water rights allocated to the Industrial category up to this priority date year (%) Environmental: The percent of cumulative water rights allocated to the Environmental category up to this priority date year (%) Other: The percent of cumulative water rights allocated to the Other category up to this priority date year (%) In addition to the database files, there is a companion .csv file, called stateWMAs_ID.csv. The main purpose of this file is to provide the means of translating between the [###] unique identifier and the various other id of the WMA the file is attached to. This translation file has five columns, which are: basinNum: The official state given alpha-numeric identifier of the WMA basinName: the state provided English name of the area, where applicable state: the state name uniID: a unique identifier we created by concatenating the state name, and underscore, and the state numerical ID ID: a unique numeric identifier we created as a requirement for the files to be used within WBM (Grogan et al. in review) For those whom may be interested in exploring our code more in depth, we are also making available an internal data file for convenience. The file is in .RData format and contains everything described above as well as some minor additional objects used within the code calculating the cumulative curves. For completeness, here is a detailed description of the various objects which can be found within the .RData file: states: A character vector containing the state names for those states in the study region. More importantly, the index of the state name is also the index in which that state's data can be found in the various following list objects. For example, if California is the second index in this object, the data for California will also be in the second index for each accompanying list. wmasRightsPersGrd: A list of nested lists of data frames. Each state has one entry in the top level list and consist of a list with one data frame per WMA. The data frames in this object are the same as the above described .csv files, specifically for groundwater. wmasRightsPersSur: A list of nested lists of data frames. Each state has one entry in the top level list and consist of a list with one data frame per WMA. The data frames in this object are the same as the above described .csv files, specifically for surface water. wmasRightsPersUnsp: A list of nested lists of data frames. Each state has one entry in the top level list and consist of a list with one data frame per WMA. The data frames in this object are the same as the above described .csv files, specifically for water from an unspecified source. wmasRightsTotsGrd: A list of nested lists of data frames. Each state has one entry in the top level list and consist of a list with one data frame per WMA. The values in the data frames of this object would be the equivalent to multiplying the CUML column by the sector percentage columns to create the sector total ft3s-1, specifically for groundwater. This object is provided for convenience to check out original total values, if desired. wmasRightsTotsSur: A list of nested lists of data frames. Each state has one entry in the top level list and consist of a list with one data frame per WMA. The values in the data frames of this object would be the equivalent to multiplying the CUML column by the sector percentage columns to create the sector total ft3s-1, specifically for surface water. This object is provided for convenience to check out original total values, if desired. wmasRightsTotsUnsp: A list of nested lists of data frames. Each state has one entry in the top level list and consist of a list with one data frame per WMA. The values in the data frames of this object would be the equivalent to multiplying the CUML column by the sector percentage columns to create the sector total ft3s-1, specifically for water from an unspecified source. This object is provided for convenience to check out original total values, if desired. wmaIDByState: A list of data frames which contain the by state allocation by WMA matrix. There is a matrix for each state in a different data frame within the object.

Economics↗
DOE OSTI2025

Integrated GW Farm ABM

This Data Repository includes data used for the integrated groundwater- farm ABM model, raw model output from scenario ensemble, and processed outputs that isolate the groundwater storage depletion outcomes for the 35,000 farm cells. Model Inputs: Farm ABM Inputs: This folder contains the input data used by the integrated groundwater - farm ABM modelling script (Python file) used for the high performance computing (HPC) experiments. The sub-folder "data inputs" contains all of the farm attribute data, while the three files in the folder have the hydrogeological data lookup table (NLDAS Cost Curve Attributes.csv), a lookup table (Theis well function table.csv) for the groundwater cost curve function, and the farm indexes and corresponding NLDAS ids for all of the cells run in this experiment (nldas farms subset final.csv). NLDAS Cost curve hydrogeological data: Hydrogeological data aggregated to 1/8 degree resolution and aligned with the NLDAS grid. Parameters include: water depth below ground surface [meters], subsurface porosity [unitless], aquifer depth from ground surface to aquifer bottom [meters], annual average recharge (USGS: mm, Doll: meters), and three different hydraulic conductivity (K) values (meters/day). The three K values represent the mean value from Gleeson et al. (2018), one standard deviation above the mean from Gleeson et al. (2018), and the de Graaf et al. 2020 modifications to certain lithologies. Additional information about these datasets and their processing are documented in the supplement to Yoon et al. 2025 (in review). Output: Raw outputs: This folder contains a .zip file that has model outputs for the entire scenario ensemble. There is one csv for each farm id, using the format "farm farmid cases.csv". The relationship between the farm id and NLDAS id is defined by the "nldas farms subset final.csv" located in the Farm ABM Inputs folder. Each csv has 625 rows, corresponding to 625 combinations of different scenario parameter values. Each row (scenario) represents the outcome of a 100 year simulation. Columns define scenario settings and summary statistics for each scenario. The first four columns define the scenario settings: "hydro ratio," "econ ratio," "K scenario," and "gamma scenario." The hydro and econ ratios are values passed to the modeling script that influence multipliers for other model parameters, as documented in the supplement to Yoon et al. 2025 (in review). The gamma multiplier is a coefficient multiplier applied to the baseline gamma values (values below 1 represent lower unobserved costs compared to baseline, values above 1 represent higher costs). The K scenario names represent K values of: "low": 0.5 m/d, "int 1": 2.5 m/d, "int 2": 10 m/d, "high": 50 m/d, and "gleeson": mean Gleeson K value. "Perc vol depleted" is the fraction of groundwater depleted at the end of the 100 simulation. Processed Output: Derived depletion outcomes from raw outputs: All of the individual csv files from the Raw outputs were aggregated into a single file that has the scenario settings and fraction depletion "Perc vol depleted" for every farm cell, for every scenario. The other two files define relationships between the farm id, NLDAS id, and local and major aquifer units, used for aquifer-level depletion analysis.

Agent based modeling↗
DOE OSTI2026

Sensitivity Analysis of Drivers Water Shortage in the Los Angeles Region During Drought

The code and detailed step-by-step instructions for generating the model output data, processing results, and analysis and plotting are provided at https://github.com/IMMM-SFA/Ferencz_et_al_2026_ER_Water. The PyArtes model is a python adaptation of the Artes model. PyArtes uses many of the same input data and optimization model architecture as Artes. Documentation for the PyArtes model is provided in the Supplement to the paper. The primary data product are simulated monthly water shortages for indoor and outdoor demand under a large ensemble of drought scenarios (>13,000). The droughts are hypothetical and are not based on historical time series data of supply sources - though historical data did help inform ranges explored for supply parameters. Demands are informed by recent 2017-2021 water supply data. Demands used for the model can be accessed at https://github.com/IMMM-SFA/Ferencz_et_al_2026_ER_Water. Simulations resolve demand for over 90 water providers in the study region. The results report 36 months of water shortage data for each indoor and outdoor demand node. The study also developed a multilayer perceptron (MLP) neural network trained on a subset of the simulated shortage ensemble to emulate worst annual water shortage for a given set of parameter multipliers -- provided the parameter values fall within the ranges sampled in the ensemble. Emulated water shortages for synthetic ensembles are in the MLP-generated shortages folder. The MLP model was used to generate larger ensembles to support Sobol analysis that would have been extremely computationally expensive to simulate. Datasets provided in this repository*: Simulated shortages. These results are used for the analysis for Figures 5, 8, and 9 in the paper, and also to train the MLP emulator. .zip file containing outputs for the 13,312 scenario ensemble. Separate .csv files for indoor and outdoor shortage for each scenario. Rows = demand ids (~100), Columns = months (36) Units = acre-feet/month of shortage (shortage = monthly demand - supply). 1 acft = 1233.48 m^3 .csv files of aggregated shortages derived from the 13,312 ensemble Rows = scenarios (13,312), Columns = demand ids (~100) Units = acre-feet/year (either worst annual shortage or total shortage over the 3-year drought) .csv file of the parameter multipliers scenarios for the ensemble .csv file of the parameter ranges and baseline values the multipliers were applied to MLP-generated shortages. These results are used for Figures 4, 6, and 7 in the paper. mwd higher folder: scenario ensembles, emulated worst year total shortages (acft), and Sobol results Emulated shortages. Rows = scenarios, columns = demand ids, units acft Sobol results. Rows = demand ids, columns Sobol (S1, ST, or 95% confidence interval) value for each parameter mwd lower folder: scenario ensembles, emulated worst year total shortages (acft), and Sobol results same organization as mwd higher MLP performance: performance metrics (R^2, RMSE, BIAS, MAPE) for the testing subset (20% or 2,662 scenarios) and simulated vs emulated worst year shortage (acre-feet/year) for every demand node, MWD wholesale regions, and the entire study region (LAC). Supporting data for figures. Figure plotting scripts in the associated GitHub repo. These files support analysis and visualization. Geospatial Data used for plotting simulated water shortages and Sobol results. Dictionary of full names for demand nodes in the model and estimates of water supply by source type informed by Artes input files and California Urban Water Management Planning data: https://water.ca.gov/Programs/Water-Use-And-Efficiency/Urban-Water-Use-Efficiency/Urban-Water-Management-Plans *Readme files provided for each folder.

drought↗
DOE OSTI2021

Distributed control of electric power grids

A method for determining power output levels of a plurality of nodes in an electric power system includes receiving, at a first node of the plurality of nodes, voltage information and multipliers of all neighboring nodes of the first node within the electric power system, determining, by the first node, a local power generation and a local voltage using the voltage information and the multipliers of the neighboring nodes and distributing the local power generation and the local voltage to the neighboring nodes, determining, by the first node, an estimated voltage of each of the neighboring nodes and distributing the estimated voltage to each of the neighboring nodes, and updating, by the first node, a local multiplier using the voltage information received from the neighboring nodes and the estimated voltage of each of the neighboring nodes determined by the node.

Ghosh, Soumyadip↗
DOE OSTI2021

Analysis of Lobe Power Calculator and Indication System with Physics and Cycle Based Models

The Advanced Test Reactor (ATR) at INL measures reactor power through two methods, thermal and Nitrogen-16 (N-16) activity. Water power calculator (WPC) is a thermal power system that measures flow and temperature to determine the thermal quadrant powers. The N-16 system utilizes a beta chamber detector that outputs reactivity levels to calculate lobe power through an algorithm called lobe power calculation and indication system (LPCIS). The LPCIS utilizes the N-16 system and the WPC system to determine reactor core power levels. The WPC provides accurate calculations of quadrant and total reactor thermal power. With the use of WPC measurements, thermal-to-N-16 (T2N) power ratios are produced to determine if the two indication systems agree on core power. Relative magnitude equations are used to utilize N-16 coefficients and multipliers to improve the indications of the LPCIS. This is crucial for maintaining safety limits. Currently, the LPCIS system uses linear equations and matrices to calculate lobe power through multipliers and coefficients. Advancements in technology and system upgrades have assisted system engineers at INL with the objective to reach a more dynamic system. In return the system demonstrates an increase in the accuracy of power reading while maintaining safety margins. The new proposed coefficient and multiple method implements a cycle specific, and a physics-based model intended for changing coefficients during operation. This calibration experiment focuses on power splitting, outer shim and neck shim, as well as fuel burning into the RDAS weighting factors. Results show that the Physics learning method yields a smaller error margin inside of the desired power range for the data set 166-A. This proved true for both constrained and unconstrained testing. This is most likely due to the physics data fitting approach that resulted in favorable coefficients and multipliers for 166-A. Continuing to improve the physics-based model will help improve the power accuracy of the LPCIS system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Optimal Control of the Energy-Saving Hybrid Hydraulic-Electric Architecture (HHEA) for Off-Highway Mobile Machines

Most off-highway constructions and agriculture equipment use hydraulics, which has unmatched power density, for power transmission and throttling as a means for control. A novel hybrid hydraulic-electric architecture (HHEA) has recently been proposed to improve efficiency for high-power machines that would have been cost-prohibitive to electrify directly. HHEA uses a set of common pressure rails (CPRs) to transmit the majority of power hydraulically and small electric motor drives to modulate that power and to achieve precise control. This article proposes a computationally efficient Lagrange multiplier method (LMM) for computing the optimal sequence of pressure rail selections to minimize energy use. This is needed to evaluate HHEA's energy-saving potential and for iterative architecture design and sizing. An interesting complication is that the cost function is not fully defined until the candidate control sequence is fully specified. This issue is dealt with by decomposing the original problem into a set of sub-problems with additional constraints that can be solved efficiently. Computational effort can be further reduced if actuators are optimized individually instead of together. However, additional steps are required to prevent the constraint functions from becoming discontinuous with respect to the Lagrange multipliers, which is necessary for meeting the constraints. Lastly, a case study of a construction machine demonstrates the efficacy of the method and shows that the HHEA reduces energy consumption by 68%-73% compared to the baseline load-sensing architecture.

Lagrange multiplier↗
DOE OSTI2021

Sensitivity analysis of in-pile critical heat flux experiments in TREAT for characterization of RIA power-transient effects

A reactivity-initiated accident (RIA) is one type of postulated design basis accident (DBA) that can cause a departure from nucleate boiling (DNB) event in pressurized-water reactors (PWRs). A DNB occurrence and its consequences depend on the thermophysical properties of the fuel components and coolant, characteristics of the transient energy insertion into the fuel rod, and the onset of the critical heat flux (CHF) phenomenon. To leverage the restart of the Transient Reactor Test (TREAT) Facility, an effort is currently underway to better understand the cladding-to-coolant heat transfer mechanisms and the CHF phenomenon under fast-transient irradiation conditions. This paper characterizes the impact of power transients on the thermal-hydraulic behavior of a TREAT Facility reactor heater rodlet CHF experiment to provide the priority of parameters that need to be investigated for an improved CHF model. Sobol sensitivity analysis methods and the Reactor Excursion and Leak Analysis Program (RELAP5-3D) code were used to identify key input parameters on the uncertainty in the prediction of peak outer- and inner-surface temperatures of the heater tube, as well as the time of the DNB event. A series of sensitivity analyses revealed the total energy deposition on the tube and the transient effects of power pulse had large impacts on the maximum temperatures. The CHF multiplier had the largest impact on the time occurrence of CHF. The overall results show the energy deposition rate in the tube is the most influencing factor to the manifestation of CHF and the resulting thermal-hydraulic behaviors of the tube. The multiplier for the CHF, which is interpreted as the predicted CHF value, has the largest Sobol indices for the time of the CHF in all cases, since it directly determines the occurrence of CHF. It is inferred that the uncertainties in the thermal-hydraulic behaviors of fuels increase with respect to the key parameters as the power pulse becomes broader, and an accurate estimation of the energy deposition rate is required to reduce the uncertainty in the evaluation of the integrity of fuel if the CHF is expected to occur near the peak power. Therefore, the outputs are expected to provide rigorous interpretation of ongoing in-pile CHF experiments in the TREAT Facility reactor regarding thermal-hydraulic behavior of the fuel system aiming for a new transient in-pile CHF model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

On restricting first order form of gauge theories to one-loop order

Highlights: • We prove the equivalence of the first and second-order formalism. • We use Lagrange multiplier fields to limit radiative effects to one-loop order. • We prove that the quantized theory is both renormalizable and unitary. The first order form of the Yang–Mills and Einstein–Hilbert actions are quantized, and it is shown how Green’s functions computed using the first and the second order form of these theories are related. Next we show how by use of Lagrange multiplier fields (LM), radiative effects beyond one-loop order can be eliminated. This allows one to compute Green’s functions exactly without loss of unitarity. The consequences of this restriction on radiative effects are examined for the Yang–Mills and Einstein–Hilbert actions. In these two gauge theories, we find that the quantized theory is both renormalizable and unitary once the LM field is used to eliminate effects beyond one-loop order.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Restricting loop expansions in gauge theories coupled to matter

Highlights: • Lagrange multiplier fields are used to restrict radiative effects to one-loop order. • This approach is exemplified in the Yang-Mills theory coupled to a scalar matter field. • We showed that it can be used with the Einstein–Hilbert action coupled to matter. • The resulting quantum gravity theory is both renormalizable and unitary. Quantizing any model in which a Lagrange multiplier (LM) field is used to restrict field configurations to those that satisfy the classical equations of motion, leads to at most one-loop radiative corrections. This approach can be used with both the Yang–Mills (YM) and Einstein–Hilbert (EH) action; the resulting theory is both renormalizable and unitary, has a positive energy spectrum and has no negative norm states contributing to physical processes. Although this approach cannot be consistently used with scalar fields alone, scalar fields can be coupled to gauge fields so that loop effects in the gauge sector are restricted to one-loop order in a way that satisfies the usual criterion for a consistent quantum field theory. The tree-level diagrams are those of the classical theory in which the metric couples to the energy–momentum tensor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗