Search NASA⌕ Search

SEARCH · Search NASA

Results for “Distribution systems state 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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Integral Reactor Containment Condensation Model and Experimental Validation

This NEUP funded project, NEUP 12-3630, is for experimental, numerical and analytical studies on high-pressure steam condensation phenomena in a steel containment vessel connected to a water cooling tank, carried out at Oregon State University (OrSU) and the University of Wisconsin at Madison (UW-Madison). In the three years of investigation duration, following the original proposal, the planned tasks have been completed: (1) Performed a scaling study for the full pressure test facility applicable to the reference design for the condensation heat transfer process during design basis accidents (DBAs), modified the existing test facility to route the steady-state secondary steam flow into the high pressure containment for controllable condensation tests, and extended the operations at negative gage pressure conditions (OrSU). (2) Conducted a series of DBA and quasi-steady experiments using the full pressure test facility to provide a reliable high pressure condensation database (OrSU). (3) Analyzed experimental data and evaluated condensation model for the experimental conditions, and predicted the prototypic containment performance under accidental conditions (UW-Madison). A film flow model was developed for the scaling analysis, and the results suggest that the 1/3 scaled test facility covers large portion of laminar film flow, leading to a lower average heat transfer coefficient comparing to the prototypic value. Although it is conservative in reactor safety analysis, the significant reduction of heat transfer coefficient (50%) could under estimate the prototypic condensation heat transfer rate, resulting in inaccurate prediction of the decay heat removal capability. Further investigation is thus needed to quantify the scaling distortion for safety analysis code validation. Experimental investigations were performed in the existing MASLWR test facility at OrST with minor modifications. A total of 13 containment condensation tests were conducted for pressure ranging from 4 to 21 bar with three different static inventories of non-condensable gas. Condensation and heat transfer rates were evaluated employing several methods, notably from measured temperature gradients in the HTP as well as measured condensate formation rates. A detailed mass and energy accounting was used to assess the various measurement methods and to support simplifying assumptions required for the analysis. Condensation heat fluxes and heat transfer coefficients are calculated and presented as a function of pressure to satisfy the objectives of this investigation. The major conclusions for those tests are summarized below: (1) In the steam blow-down tests, the initial condensation heat transfer process involves the heating-up of the containment heat transfer plate. An inverse heat conduction model was developed to capture the rapid transient transfer characteristics, and the analysis method is applicable to SMR safety analysis. (2) The average condensation heat transfer coefficients for different pressure conditions and non-condensable gas mass fractions were obtained from the integral test facility, through the measurements of the heat conduction rate across the containment heat transfer plate, and from the water condensation rates measurement based on the total energy balance equation. 15 (3) The test results using the measured HTP wall temperatures are considerably lower than popular condensation models would predict mainly due to the side wall conduction effects in the existing MASLWR integral test facility. The data revealed the detailed heat transfer characteristics of the model containment, important to the SMR safety analysis and the validation of associated evaluation model. However this approach, unlike separate effect tests, cannot isolate the condensation heat transfer coefficient over the containment wall, and therefore is not suitable for the assessment of the condensation heat transfer coefficient against system pressure and noncondensable gas mass fraction. (4) The average condensation heat transfer coefficients measured from the water condensation rates through energy balance analysis are appropriate, however, with considerable uncertainties due to the heat loss and temperature distribution on the containment wall. With the consideration of the side wall conduction effects, the results indicate that the measured heat transfer coefficients in the tests is about 20% lower than the prediction of Dehbi’s correlation, mainly due to the side wall conduction effects. The investigation also indicates an increase in the condensation heat transfer coefficient at high containment pressure conditions, but the uncertainties invoked with this method appear to be substantial. (5) Non-condensable gas in the tests has little effects on the condensation heat transfer at high elevation measurement ports. It does affect the bottom measurements near the water level position. The results suggest that the heavier non-condensable gas is accumulated in the lower portion of the containment due to stratification in the narrow containment space. The overall effects of the non-condensable gas on the heat transfer process should thus be negligible for tall containments of narrow condensation spaces in most SMR designs. Therefore, the previous correlations with noncondensable gas effects are not appropriate to those small SMR containments due to the very poor mixing of steam and non-condensable gas. The MELCOR simulation results agree with the experimental data reasonably well. However, it is observed that the MELCOR overpredicts the heat flux for all analyzed tests. The MELCOR predicts that the heat fluxes for CCT’s approximately range from 30 to 45 kW/m2 whereas the experimental data (averaged) ranges from about 25 to 40 kW/m2. This may be due to the limited availability of liquid film models included in MELCOR. Also, it is believed that due to complex test geometry, measured temperature gradients across the heat transfer plate may have been underestimated and thus the heat flux had been underestimated. The MELCOR model predicts a film thickness on the order of 100 microns, which agrees very well with film flow model developed in this study for scaling analysis. However, the expected differences in film thicknesses for near vacuum and near atmospheric test conditions are not significant. Further study on the behavior of condensate film is expected to refine the simulation results. Possible refinements include but are not limited to, the followings: CFD simulation focusing on the liquid film behavior and benchmarking with experimental analyses for simpler geometries. 16 1 INTRODUCTION This NEUP funded project, NEUP 12-3630, is for experimental, numerical and analytical studies on high-pressure steam condensation phenomena in a steel containment vessel connected to a water cooling tank, carried out at Oregon State University (OrSU) and the University of Wisconsin at Madison (UW-Madison). The experimental results are employed to validate the containment condensation model in reactor containment system safety analysis code for integral SMRs. Such a containment condensation model is important to demonstrate the adequate cooling. In the three years of investigation, following the original proposal, the following planned tasks have been completed: (1) Performed a scaling study for the full pressure test facility applicable to the reference design for the condensation heat transfer process during design basis accidents (DBAs), modified the existing test facility to route the steady-state secondary steam flow into the high pressure containment for controllable condensation tests, and extended the operations at negative gage pressure conditions (OrSU). (2) Conducted a series of DBA and quasi-steady experiments using the full pressure test facility to provide a reliable high pressure condensation database (OrSU). (3) Analyzed experimental data and evaluated condensation model for the experimental conditions, and predicted the prototypic containment performance under accidental conditions (UW-Madison). The results are applicable to integral Small Modular Reactor (SMR) designs, including NuScale, mPower, Westinghouse SMR, Holtec-160 and other integral reactors with small containments of relatively high pressures under accidental conditions. Testing has been conducted at the OrSU laboratory in the existing MASLWR (Multi-Application Small Light Water Reactor) integral test facility sponsored by the US Department of Energy. Its highpressure stainless steel containment model (~2 MPa) is scaled to the NuScale SMR currently under development at NuScale Power, Inc.. Minor modifications to the model containment have been made to control the non-condensable gas fraction and to utilize the secondary loop stable steam flow for condensation testing. UW-Madison has developed a containment condensation model, which leveraged previous validated containment heat transfer work carried out at UW-Madison, and extended the range of applicability of the model to integral SMR designs that utilize containment vessels of high heat transfer efficiencies. In this final report, the research background and literature survey are presented in Chapter 2 and 3, respectively. The test facility description and modifications are summarized in Chapter 4, and the scaling analysis is introduced in Chapter 5. The tests description, procedures, and data analysis are presented in Chapter 6, while the numerical modeling is presented in Chapter 7, followed by a conclusion section in Chapter 8.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Value of Geothermal Energy Storage for Supply-Side and Demand-Side Applications

This report presents the results of a study examining the value potential for geothermal energy storage (GES), a long-duration energy storage resource that stores thermal and/or geomechanical energy in the subsurface. GES could benefit the overall U.S. power system by temporally shifting electricity generation (supply-side) or meeting building heating and cooling load (demand-side). This report analyzes supply-side and demand-side opportunities independently because of differences in applications and models. Currently there is significant uncertainty about the development costs for GES, with only a limited number of demonstration plants for electric energy storage and building heating and cooling storage developments. In this report, we estimate the value of supply-side and demand-side GES to the bulk power system in the contiguous United States. Because of the significant uncertainty about GES development costs, this analysis does not consider GES deployment costs but instead focuses on the value of GES to the U.S. electricity system. The estimated values of GES provide reference points for economically competitive commercial cost targets. Supply-side GES is modeled as part of an enhanced geothermal system (EGS) generation plant in NREL's Regional Energy Deployment System (ReEDS) capacity expansion model (Ho et al. 2021). In contrast to conventional geothermal plants, which generate constant power, EGS plants have unique features that may allow for in-reservoir energy storage for flexible generation. Demand-side GES for heating and cooling, including seasonal hot and cold storage and short-duration heat pump storage, is incorporated into a price-taker model using Cambium electricity marginal cost projections. To establish an upper bound for the value of GES, analysis focused on favorable scenarios for storage with high generation from zero marginal cost, variable renewable energy resources. High penetrations of variable renewable energy generation can increase hourly electricity price variability, which increases the value of temporal energy arbitrage for storage technologies like GES.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Optimization-Based Dynamic Voltage Support of Microgrids Using Energy Storage Systems

A microgrid network is characterized by a high R/X ratio, making the voltage more sensitive to active power changes compared to bulk power systems, where the voltage is regulated primarily by reactive power. Due to its sensitivity, voltage control approaches for microgrids should also consider the active power input coupling, making it very different from conventional power systems. Additionally, as the energy costs associated with active and reactive powers are different and the operational conditions of microgrids connected to active distribution systems vary over time, the ideal controller to provide voltage support must be flexible enough to handle these technical and operational constraints. This paper proposes a model predictive control approach to provide dynamic voltage support using energy storage systems. This approach uses a simplified predictive model of the system to solve the model predictive control problem. By proper selection of model predictive control weighting parameters, the quality of service provided can be adjusted to achieve the desired performance. A simulation study in MATLAB/Simulink validates the proposed approach for the Cordova, Alaska microgrid. Results show that the performance of the voltage support can be adjusted depending on the choice of weight and constraints of the controller.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Electronic Spin Relaxation and Clustering in High-Pressure High-Temperature Synthesized Microcrystalline Diamond Particles with Reduced Nitrogen Content

The negatively charged nitrogen-vacancy (NV – ) color center in diamonds is widely studied because of numerous applications of this unique quantum system in sensing and quantum information sciences. While substitutional nitrogen is required to form the NV – centers in diamond, it also yields other paramagnetic defects─primarily the neutrally charged substitutional nitrogen centers (P1)─that decrease NV – spin coherence, which in turn degrades performance in applications. Herein, we investigate high-pressure high-temperature synthesized diamond microparticles (ca. 140–185 μm) having lower─ranging from 3 to 38 ppm─than the typical nitrogen content of type 1b diamond ( ca. 100 ppm and higher) typically used for the production of fluorescent diamond particles with NV – centers. A suite of electron paramagnetic resonance, optically detected magnetic resonance, and nuclear magnetic resonance methods are used to characterize spin properties of P1 and NV – centers in the particles. Upon decreasing the nitrogen content from 29 to 3 ppm, the ensemble NV – T 2 relaxation time increased by about 3-fold as measured directly in the Hahn Echo experiment at magnetic field of 1.2 T. Analysis of electronic relaxation of P1 centers revealed the existence of at least two distinct populations of P1 centers, consisting of fast and slower relaxing spins and allowed for an estimation of local concentrations. Even with <10 ppm nitrogen contents, the analysis indicated a highly heterogeneous distribution of P1 centers, suggesting the possibility of P1 spin clustering even at low nitrogen concentrations. The combined data demonstrate that the particles prepared from HPHT diamond with a low nitrogen content offer improved spin properties that are beneficial for NV – sensing applications.

Carbon↗
DOE OSTI2014

Hadron collider potential for excited bosons search

The e {sup +} e {sup −} and μ{sup +}μ{sup −} dilepton final states are the most clear channels for a new heavy neutral resonance search. Their advantage is that usually in the region of expected heavy-mass resonance peak the main irreducible background, from the Standard Model Drell-Yan process, contributes two orders of magnitude smaller than the signal. In this paper we consider the future prospects for search for the excited neutral Z*-bosons. The bosons can be observed as a Breit-Wigner resonance peak in the dilepton invariant mass distributions in the same way as the well-known extra gauge Z′ bosons. However, the Z* bosons have unique signatures in transverse momentum, angular and pseudorapidity distributions of the final leptons, which allow to distinguish them from the other heavy neutral resonances. At present only the ATLAS Collaboration has looked for such new excitations at the Large Hadron Collider and has published its results for 7 TeV collision energy. After successful comparison of our evaluation with these official results we present our estimations for the discovery potential and the exclusion limits on the Z*-boson search in pp collisions at higher centre-of-mass energies and different luminosities. In particular, LHC Run 2 can discover Z*-boson with its mass up to 5.3 TeV, while the High Luminosity LHC can extend that reach to 6.2 TeV. The High Energy LHC (with collision energy of 33 TeV) will be able to probe two times heavier resonance masses at the same integrated luminosities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Aerosol Vertical Profiling at Oliktok Point (AVPOP) Field Campaign Report

Aerosols are an important component of the atmospheric system through their various impacts on climate. Notably, the largest uncertainty of the energy budget (i.e., radiative forcing estimate) is the aerosol indirect effect. Aerosol-induced microphysical modifications influence cloud lifetime and albedo. However, constraining aerosol-cloud impacts in climate models, specifically in regions such as the Arctic, remains a significant challenge due to limited observations. The efficacy of an aerosol to impact the arctic energy budget largely depends on its concentration, composition, and source. While ground-based measurements will provide a consistent, high-resolution time series of aerosol observations, the arctic atmosphere can be stratified, and thus what is observed on the ground may not be representative of the aerosol aloft. A key gap in arctic aerosol research is routine vertical profiling of abundance and source characterization. Previous work outside of the Arctic has focused on using launched or tethered balloon systems (LBS and TBS, respectively) to measure vertical profiles of aerosol size distributions and has demonstrated the utility of such observations (e.g., Hofmann et al. 1989; Hofmann 1993; Iwasaka et al. 2003; Kim et al. 2003; Tobo et al. 2007). More specifically, particle spectrometers have also been deployed via TBS systems (de Boer et al. 2018; Greenberg et al. 2009; Maletto et al. 2003; Renard et al. 2016; Wehner et al. 2007; Siebert et al. 2004) affording information on aerosol layer locations and evolution by means of multiple profiles. Filter samples for post-processing of aerosol composition and ice nucleation properties have also been conducted in mid-latitudes (Creamean et al. 2018). However, while these studies present novel measurements of aerosol vertical profiling via balloon systems, few exist to date in the Arctic. The U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility TBS is composed of a variety of winches, balloons, and instrumentation that flies within ARM’s Restricted Airspace (R-2204) at the third ARM Mobile Facility (AMF3) at Oliktok Point, Alaska (Figure 1). The system payload comprises instruments that can simultaneously collect in situ data on horizontal winds, ice microphysics, turbulence, thermodynamic state, aerosols, and the cloud-top environment (https://www.arm.gov/capabilities/instruments/tbs). Recent campaigns at AMF3 have involved the use of the ARM helikite, which is a balloon-kite hybrid that uses lighter-than-air principles to obtain its initial lift and a kite to achieve stability and dynamic lift. Additional information on the specifications of the helikite TBS can be found in the ARM TBS handbook (Dexheimer 2018).

Oliktok Point, North Slope of Alaska,, tethered ba↗
NASA NTRS1998

Acquisition and Early Losses of Rare Gases from the Deep Earth

Direct observations show that the deep Earth contains rare gases of solar composition distinct from those in the atmosphere. We examine the implications of mantle rare gas characteristics on acquisition of rare gases from the solar nebula and subsequent losses due to a large impact. Deep mantle rare gas concentrations and isotopic compositions can be obtained from a model of transport and distribution of mantle rare gases. This model assumes the lower mantle closed early, while the upper mantle is open to subduction from the atmosphere and mass transfer from the lower mantle. Constraints are derived that can be incorporated into models for terrestrial volatile acquisition: (1) Calculated lower-mantle Xe-isotopic ratios indicate that the fraction of radiogenic Xe produced by I-129 and Pu-244 during the first about 10(exp 8) yr was lost, a conclusion also drawn for atmospheric Xe. Thus, either the Earth was made from materials that had lost >99% of rare gases about (0.7-2) x 10(exp 8) yr after the solar system formed, or gases were then lost from the fully formed Earth. (2) Concentrations of 3He and 20Ne in the lower mantle were established after these losses. (3) Neon-isotopic data indicates that mantle Ne has solar composition. The model allows for solar Ar/Ne and Xe/Ne in the lower mantle if a dominant fraction of upper mantle Ar and Xe are subduction-derived. If Earth formed in the presence of the solar nebula, it could have been melted by accretional energy and the blanketing effect of a massive, nebula-derived atmosphere. Gases from this atmosphere would have been sequestered within the molten Earth by dissolution at the surface and downward mixing. It was found that too much Ne would be dissolved in the Earth unless the atmosphere began to escape when the Earth was only partially assembled. Here we consider conditions required to initially dissolve sufficient rare gases to account for the present lower mantle concentrations after subsequent losses at 10(exp 8) yr. It is assumed that equilibration of the atmosphere with a thoroughly molten mantle was rapid, so that initial abundances of gases retained in any mantle layer reflected surface conditions when the layer solidified. For subsequent gas loss of 99.5% and typical solubility coefficients, a total pressure of 100 atm was required for an atmosphere of solar composition. Calculations of the pressure at the base of a primordial atmosphere indicate that this value might be exceeded by an order of magnitude or more for an atmosphere supported by accretional energy. Surface temperatures of about 4000 K would have been produced, probably high enough to melt the deep mantle. Initial distributions of retained rare gases would then be determined by the history of surface pressure and temperature during mantle cooling and solidification, i.e., the coupled cooling of Earth and atmosphere. The Earth's thermal state was determined by its surface temperature and the efficiency of convection in the molten mantle, estimated to be sufficient to maintain an adiabatic gradient. Because the melting curve is steeper than the adiabat, solidification of the mantle proceeded outward from the interior. Incorporation of atmospheric gases in the mantle therefore occurred over a range in surface temperature of a few thousand degrees Kelvin. The thermal state of the atmosphere was controlled by total luminosity of the Earth (energy) released by accreting planetesimals and the cooling Earth), nebular temperature and pressure, and atmospheric opacity. The energy released by accretion declined with time as did nebular pressure. Analytical solutions for an idealized (constant opacity radiative atmosphere show that declining energy sources under constant nebular conditions result in slowly diminishing surface temperature but dramatically increasing surface pressure. For such an atmosphere with declining nebular pressure but constant total luminosity, surface pressure decreases gradually with decreasing temperaure. A decline in accretion luminosity might be compensated by energy released as the mantle cools for about 10(exp 5) year, after which luminosity must decline. The total complement of dissolved rare gases will depend on the particular evolutionary path determined by the declining accretional luminosity, the Earth thermal history, removal of the nebula, and opacity variations of the atmosphere. Models for these coupled evolutionary histories for Earth's acquisition of nebular-derived noble gases are in progress. The later losses required at about 10(exp 8) yr (depleting the interior concentrations of the sequestered solar gases by a factor of > 100) were presumably related to the major impact in which the Moon formed.

Porcelli, D.↗
DOE OSTI2026

Real-Time Multiregional Market-to-Market Congestion Management Through Exchange of Relief Cost Curve

This paper introduces a novel method for multiregional market-to-market (M2M) coordinated congestion management. It identifies shortcomings in existing M2M approaches, where Regional Transmission Organizations (RTOs) exchange shadow prices and relief requests to optimize congestion relief allocations across interconnected regions. Two methods are proposed to enhance flow and price convergence. The first method proposes that both Regional Transmission Organizations (RTOs) use state-estimator flows directly to determine relief requirements, eliminating delays and potential oscillations caused by using market flows calculated from the prior period under existing M2M approach. The second method involves exchanging transmission relief cost curves, enabling each RTOs to integrate other RTOs' relief costs curve into its real-time security-constrained economic dispatch (SCED). This method can effectively extend the coordination to multiple transmission lines and across more than two RTOs. The alternating direction method of multipliers (ADMM) is also applied to the M2M coordination problem and compared with the proposed methods. Case studies on small and large-scale systems demonstrate the effectiveness of these approaches.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Heavy-Duty Vehicle Activity Updates for MOVES Using NREL Fleet DNA and CE-CERT Data

The U.S. Environmental Protection Agency's (EPA's) Motor Vehicle Emission Simulator (MOVES) is a publicly available tool used by researchers and policymakers to help understand motor vehicle emission sources at a national, county, and project level. Estimates of heavy-duty activity in the most recent version of the model at the time this work was conducted, MOVES2014, was identified as an area in need of improvement. The start activity in MOVES2014 is based on a limited and dated data set. In addition, MOVES2014 relies on drive cycles that represent on-network activity but do not account for idling activity that occurs on off-network roads, such as at a distribution center, while the truck is queuing or during loading and unloading. As a result, MOVES2014 may currently underestimate the number of starts and idle and soak time for heavy-duty trucks in real-world operation. The National Renewable Energy Laboratory (NREL) has previously leveraged its expansive Fleet DNA database of heavy-duty vehicles to idle and start activity for six of the nine heavy-duty vehicle source types of classes in the MOVES model. The data available in Fleet DNA from 416 conventional, diesel-powered vehicles provided activity estimates from more than 120,000 hours of operation throughout 14,682 vehicle days between October 2006 and January 2016. NREL calculated start fraction, starts per day, soak fraction, and idle fraction by hour of the day for each vehicle type, state, and vocation, and provided results in .CSV files that can be translated to MOVES table inputs. The idle and start activity from this initial analysis of Fleet DNA data was used to develop default idle and start data for heavy-duty vehicles in MOVES3. Satisfied with the results from the Fleet DNA data used for MOVES3, the EPA asked NREL to extend this start/soak/idle analysis using additional data from a larger number of vehicles for a potential future update to the MOVES model. Such a data set was achieved from a project led by the University of California at Riverside, College of Engineering, Center for Environmental Research & Technology (CE-CERT) and funded by California Air Resources Board. Specifically, this data set consists of 90 heavy-duty vehicles operated mainly in California, which can be separated into five of the nine heavy-duty vehicle classes in the MOVES model. In addition, the heavy-duty activity database collected by CE-CERT provided activity estimates from more than 44,000 hours of operation throughout 4,724 vehicle days between November 2014 and September 2016. This report details the analysis of the heavy-duty activity database collected from the University of California at Riverside by providing graphical analysis and context for the start, soak, and idle distributions. The comparison of the related results from both the Fleet DNA and CE-CERT data sets are documented as well.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Herbaceous Feedstock 2018 State of Technology Report

The U.S. Department of Energy (DOE) promotes the production of advanced liquid transportation fuels from lignocellulosic biomass by funding fundamental and applied research that advances the state of technology (SOT). As part of its involvement with this mission, Idaho National Laboratory (INL) completes an annual SOT report for biomass feedstock logistics. This report summarizes supply system impacts of Bioenergy Technologies Office (BETO)-funded research and development efforts at INL and elsewhere (such as the High-Tonnage Feedstock Logistics projects (Webb et al. 2013a, Webb et al. 2013b, Webb et al. 2013c, Webb and Sokhansanj 2014, Sokhansanj et al. 2014) that lead to improvements in feedstock supply systems. These include improvements to and observed performance of innovative harvest and collection methods, storage technologies, transportation and handling approaches, and advanced preprocessing technologies. Biomass quality and variability, and the interface between feedstock quality and conversion performance are key drivers in addition to delivered feedstock cost. In this report, we estimate the benefits of R&D improvements to individual supply system unit operations, and present the status of feedstock logistics technology development for converting biomass into biofuels. These analyses are supported by experimental data where possible, and help to align the SOT relative to the cost goals defined in the Multi-Year Program Plan. The 2018 Herbaceous SOT aligned feedstock logistic design with current biorefinery’s design capacity utilized by biochemical conversion platform. Currently biochemical conversion platform utilizes a 725,000 dry ton/year biorefiney design for the techno economic analysis. Hence, feedstock delivered cost in the 2018 Herbaceous SOT is calculated based on biorefinery’s 725,000 dry ton design capacity instead of 800, 000 dry ton capacity utilized in the 2017 Herbaceous SOT. Biomass availabilities in this SOT were updated to year 2018 data from the 2016 Billion-Ton Report (BT16) (DOE 2016a), with the exception of switchgrass, for which the 2018 Herbaceous SOT utilized the 2019 switchgrass availability data from BT16. The BT16 report (DOE 2016a) does not project switchgrass availability in 2018; the soonest switchgrass is available in the BT16 report is 2019. Therefore, availability of switchgrass for this analysis was that projected for 2019. The 2018 Herbaceous SOT incorporates same technologies utilized in the 2017 Herbaceous SOT. However, a sensitivity analysis is performed to understand the impact of variation of process parameters on those technologies on feedstock logistic cost. New R&D data that shows the variations of process parameters affecting process performance is incorporated in the 2018 SOT to measure the variations in delivered feedstock cost. The 2018 Herbaceous SOT has also provided projected delivered feedstock of 2022 design case based on near term technical target under BETO funded R&D project. Finally, updated biorefinery size of 725,000 dry ton/year was incorporated within least-cost formulation model to select optimal siting and depot scales during optimization of the least cost blend. This modification to the optimization algorithm allows the trade-off between the cost of increased supply radius and the savings from selecting biomass from higher producing counties to be assessed. Such optimization has also showed the economic benefit of decentralized depots in comparison to centralized preprocessing co-located with the biorefinery by decoupling the biorefinery and feedstock locations. The 2018 Herbaceous SOT report documents the current modeled cost of a herbaceous feedstock supply system (from harvest to the pretreatment reactor throat, including grower payment) for hydrocarbon fuel production via biochemical conversion, based on equipment and processes now available or potentially available in the near term. The modeled cost also considers both the required quality and the availability of the biomass resources. The 2018 Herbaceous SOT predicts a modeled delivered feedstock cost of $83.67/dry ton (2016$); this is a $0.23/dry ton (2016$) decrease from the 2017 Herbaceous SOT. The modification of biorefinery’s designed capacity and increased projected biomass availability in the same supply shed contributed to this modeled cost reduction. The least-cost formulation model to optimally site and scale local distributed preprocessing depots also contributed to the cost reduction by considering county-level grower payment and distance from the biorefinery as variables in the optimization algorithm. Sensitivity analysis on various process parameters that affect delivered feedstock cost in the 2018 Herbaceous SOT shows that the delivered cost could varies from $80.45-$88.83/dry ton. The top factors that causes such variations are: effective baling rate, bale density, hammer mill throughput, interest rate and storage dry matter loss.

09 BIOMASS FUELS↗
DOE OSTI2021

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs Results

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells - increasing or decreasing the fluid flow rates across the wells - and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. Data and supporting literature from a study describing a new approach combining reservoir modeling and machine learning to produce models that enable strategies for the mitigation of decreased heat and power production rates over time for geothermal power plants. The computational approach used enables translation of sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy and discovery of optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an "open-source" reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 hours, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 seconds. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs. Includes a synthetic, yet realistic, model of a geothermal reservoir, referred to as open-source reservoir (OSR). OSR is a 10-well (4 injection wells and 6 production wells) system that resembles Brady Hot Springs (a commercially operational geothermal field in Nevada, USA) at a high level but has a number of sufficiently modified characteristics (which renders any possible similarity between specific characteristics like temperatures and pressures as purely random). We study OSR through CMG simulations with a wide range of flow allocation scenarios. Includes a dataset with 101 simulated scenarios that cover the period of time between 2020 and 2040 and a link to the published paper about this project, where we focus on the Machine Learning work for predicting OSR's energy production based on the simulation data, as well as a link to the GitHub repository where we have published the code we have developed (please refer to the repository's readme file to see instructions on how to run the code). Additional links are included to associated work led by the USGS to identify geologic factors associated with well productivity in geothermal fields. Below are the high-level steps for applying the same modeling + ML process to other geothermal reservoirs: 1. Develop a geologic model of the geothermal field. The location of faults, upflow zones, aquifers, etc. need to be accounted for as accurately as possible 2. The geologic model needs to be converted to a reservoir model that can be used in a reservoir simulator, such as, for instance, CMG STARS, TETRAD, or FALCON 3. Using native state modeling, the initial temperature and pressure distributions are evaluated, and they become the initial conditions for dynamic reservoir simulations 4....

15 GEOTHERMAL ENERGY↗
DOE OSTI2019

High Energy Storage Capacity Low Cost Iron Flow Battery

A new flow battery was proposed that utilizes low cost materials: iron as the only active element, cheap aqueous electrolytes, and inexpensive separators. During charging, ferrous iron (Fe 2+ ) is oxidized to ferric iron (Fe 3+ ) at the positive electrode while it (Fe 2+ ) is reduced to form iron metal (Fe 0 ) at the negative electrode. Because iron is plated at the negative electrode during charging, conventional electrode structures couple the energy storage capacity and the power rating of the battery. In order to decouple the energy and power ratings and regain the economic advantages of a flow battery, a slurry electrode design was proposed for use as the negative electrode. The slurry electrode is made by flowing electrically conductive particles in an electrolyte containing the dissolved iron species. On charging, iron is plated onto the particles. The particles can then carry the iron metal out of the cell to be stored in external reservoirs. This project looked to develop a slurry for the negative electrode of the all-iron flow battery that can operate at 200 mA/cm 2 with <100 mV overpotential, be pumped with parasitic energy costs < 10%, and promote plating on the particles. A scalable and stackable cell design of 1100 cm 2 was developed. The project was initially under contract in mid-January 2013. The project was converted from a one-year proof of concept seedling to an additional two-year program, then extended with an 18-month plus-up, and finally with an additional 6-month approved cost-overrun. During the seeding project, the concept of plating on carbon particles in a slurry was demonstrated. However, this was accomplished with high-cost multi-walled carbon nanotubes as the slurry material. The following two-year program discovered several low-cost carbon blacks that could also be used in this application. These findings were supported by an analysis that demonstrated the effectiveness of the carbon that takes into account not just the effective electronic conductivity of the slurry, but also the surface area per volume of the slurry. It was demonstrated that stable, flowable slurries could be created with these carbons, even when the carbons where fully loaded with iron during charging. During the two-year program, laboratory-scale 50 cm 2 single-cells demonstrated charge-discharge cycling at current densities of 150 mA/cm 2 and energy efficiencies of >50%. Meanwhile, cost modeling continued to show promise for this energy storage approach. Following the conversion period, the overall objective of the nominally 18-month plus up period was to close the technical gap of a 25X scale-up and reduce risk to enable commercial investment required to develop a slurry iron flow battery product. A specific objective was to demonstrate 150 cycles of performance exceeding 70% energy efficiency with a 10-cell stack of 1156 cm 2 cells operating at 100 mA/cm 2 . The 10-cell short stack was designed to have a power delivery rating of approximately 1 kW, which represents a size commensurate with the requirements of our commercialization partner at that time. A complete battery support system for the full-size cell was designed, fabricated, and tested. This includes pumps (multiple types were considered, and constant volume hose pumps were found to be acceptable for performance and cost), electrolyte rebalance reactor, and sensors and control systems. The support system proved to be well-behaved and worked well. A new low-cost membrane was developed for this system based on a microporous polymer support with a layer of hydrogel polymer. The membrane had low resistivity, and significantly reduced pressure driven fluid flow across the membrane (as compared to the un-coated microporous support). This membrane approach was scaled-up in conjunction with a commercial toll coating company, and a large amount of membrane was manufactured to support large cell and stack testing. A new concept for an electrolyte rebalance reactor to react hydrogen gas from the negative electrolyte (from plating inefficiency) with excess ferric ion in the positive electrolyte was developed. This design was very simple, low cost, and effective. The full scale 1156 cm 2 single-cell was designed, modeled, fabricated and tested. Modeling showed that uniform flow-distribution should be obtained with the designed flow distributor. The design took into consideration the need to stack the cells, and components were fabricated for a ten-cell stack. Testing of the cell for charge and discharge demonstrated current densities of 50 mA/cm 2 with efficiencies of 50%. However, optimization and long cycling could not be accomplished because the slurry electrode would plug within the cell and block the flow. This would shut down the cell after several charge-discharge cycles. Various approaches were pursued to address this problem with some progress. However, there was not adequate time to find solutions to the problem although several promising approaches were identified for future work. The final deliverable for this proposed 18-month effort was to be a 1 kW, 6 kWh slurry flow battery, with complete balance of plant, ready for testing at an ARPA-E ‘Charges’ site with an estimated capital cost less than $\$ $45/kWh for a six-hour system (excluding power electronics). The deliverable milestone was not achieved for reasons mentioned above related to the slurry plugging. At the completion of this project, it was determined that the cost targets are achievable. The cost estimate for a 1 MW system with a one-hour capacity was estimated at $\$ $140/kWh (including power electronics) and for a six-hour system was estimated at $\$ $40/kWh (including power electronics). These costs include the battery, the balance of plant, the power conditioning and a concrete pad. This program resulted in three patents including the slurry-based iron flow battery, the rebalance reactor, and the membrane separator. At the time of this report, slurry-based iron flow battery patents have been issued in the US, China, Japan, Korea, and Europe. The other patent applications were under various stages of examination. A licensee was identified and the company initiated a scale-up and development program to take this technology to the marketplace, and CWRU staff have transferred their know-how and knowledge base to the company. This technology will enhance the economic and energy security of the United States by enabling the use of intermittent renewable energy technology, such as wind and solar power, as dispatchable resources. This report summarizes the results of this research and development program.

25 ENERGY STORAGE↗
DOE OSTI2026

NLR's FleetDNA Speed Distributions for EPA MOVES

The U.S. Environmental Protection Agency (EPA) MOtor Vehicle Emission Simulator (MOVES) models mobile sources of air toxics at small and large scales across the Unites States. Vehicle speed distributions based on the type of vehicle, type of road traveled, and hour of the day are inputs into MOVES. In conjunction with vehicle miles traveled (VMT) inputs, they are used to estimate total operation time and select the relevant driving cycles from which running operating modes and emissions are estimated. In an ongoing research effort, the National Laboratory of the Rockies (NLR) partnered with the EPA to provide default average speed data for heavy-duty vehicles based on NLR's FleetDNA database and the Bourns College of Engineering - Center for Environmental Research and Technology (CE-CERT) data set from the University of California, Riverside.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2014

Sub-Second Parallel State Estimation

This report describes the performance of Pacific Northwest National Laboratory (PNNL) sub-second parallel state estimation (PSE) tool using the utility data from the Bonneville Power Administrative (BPA) and discusses the benefits of the fast computational speed for power system applications. The test data were provided by BPA. They are two-days’ worth of hourly snapshots that include power system data and measurement sets in a commercial tool format. These data are extracted out from the commercial tool box and fed into the PSE tool. With the help of advanced solvers, the PSE tool is able to solve each BPA hourly state estimation problem within one second, which is more than 10 times faster than today’s commercial tool. This improved computational performance can help increase the reliability value of state estimation in many aspects: (1) the shorter the time required for execution of state estimation, the more time remains for operators to take appropriate actions, and/or to apply automatic or manual corrective control actions. This increases the chances of arresting or mitigating the impact of cascading failures; (2) the SE can be executed multiple times within time allowance. Therefore, the robustness of SE can be enhanced by repeating the execution of the SE with adaptive adjustments, including removing bad data and/or adjusting different initial conditions to compute a better estimate within the same time as a traditional state estimator’s single estimate. There are other benefits with the sub-second SE, such as that the PSE results can potentially be used in local and/or wide-area automatic corrective control actions that are currently dependent on raw measurements to minimize the impact of bad measurements, and provides opportunities to enhance the power grid reliability and efficiency. PSE also can enable other advanced tools that rely on SE outputs and could be used to further improve operators’ actions and automated controls to mitigate effects of severe events on the grid. The power grid continues to grow and the number of measurements is increasing at an accelerated rate due to the variety of smart grid devices being introduced. A parallel state estimation implementation will have better performance than traditional, sequential state estimation by utilizing the power of high performance computing (HPC). This increased performance positions parallel state estimators as valuable tools for operating the increasingly more complex power grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Grid Resiliency using Hydro-based Microgrids During Wildfires

The primary objective of this project is to assess the techno-economic viability of hydropower resources with supplementary investments, including microgrid controllers, storage solutions, load banks, additional generation like diesel generators, communication links, and transmission infrastructure. The aim is to enable these systems to operate autonomously (in islanded mode) and ensure grid resiliency during and after wildfire events. The project seeks to establish a comprehensive framework for setting up hydropower-based microgrids, from risk estimation to final investment decisions. These microgrids can provide critical support during emergencies and offer additional grid stability and revenue opportunities during regular operations and other system contingencies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

A neural-network-enhanced parameter-varying framework for multi-objective model predictive control applied to buildings

Management of the electrical grid is becoming more complex due to the increased penetration of alternative energy generation technologies and a broadening diversity of electric loads. This complexity creates challenges in balancing demand and generation that can increase the potential for grid instabilities. One effective way to address this issue is to leverage previously unexploited demand flexibility through advanced control strategies. In this work, we propose an advanced control method, called adaptive neural parameter-varying model predictive control (ANPV-MPC), to control the temperature and energy consumption of a building via its Heating, Ventilation, and Air Conditioning system. ANPV-MPC combines key ideas in parameter-varying control, adaptive control, and online learning strategies to bridge the gap between computationally efficient linear model predictive control and more accurate nonlinear model predictive control. The novelty in ANPV-MPC is the use of a physics-inspired Bayesian neural network to estimate the coefficients of the parameter-varying linear control model. The Bayesian neural network additionally provides uncertainty estimates, triggering online training to capture evolving building system conditions. We show that ANPV-MPC can approximate the building system dynamics with a 28.39% higher accuracy than traditional linear model predictive control, resulting in 36.23% better control performance without increasing complexity of the optimal control problem. ANPV-MPC also adapts in real time to previously unseen conditions using online learning, further improving its performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

Boil-off losses along LH2 pathway

Losses along the LH 2 pathway are intrinsic to the utilization of a cryogenic fluid. They occur when the molecule is transferred between 2 vessels (liquefaction plant to trailer, trailer to station storage, station storage to pump or compressor, then fuel cell electric vehicles …) and when the fluid is warmed up due to heat transfer with the environment. Those losses can be estimated with good accuracy using thermodynamic models based on conservation of mass and energy, providing the thermodynamic states are correctly described. Indeed, the fluid undergoes various changes as it moves along the entire pathway (2 phase transition, super-heated warming, non-uniform temperature distributions across the saturation film) and accurate equations of state and 2 phase behavior implementations are essential. The balances of mass and energy during the various dynamics processes then enable to quantify the boil-off losses. In this work, a MATLAB code previously developed by NASA to simulate rocket loading is used as the basis for the LH 2 transfer model. This code implements complex physical phenomena such as the competition between condensation and evaporation and the convection vs. conduction heat transfer as a function of the relative temperatures on both sides of the saturated film. The original code was modified to consider real gas equations of state, and some semi-empirical relationships, such as between the heat of vaporization and the critical temperature, were also replaced by a REFPROP equivalent expression, assumed to be more accurate. Non-constant liquid temperature equations were added to simulate sub-cooled conditions. It is shown that although trailer depressurization and transfer losses from the receiving vessel during a LH 2 delivery may constitute the largest sources of losses, those can be greatly reduced and even eliminated if CFR/DOT regulations are properly applied and if no-vent fill methods are implemented, particularly for a stationary LH2 storage operating at reasonably low pressure (45 psia or so). As such, it is expected that the only remaining boil-off losses for a refueling station would come from the LH 2 pump (utilization and idling), pump vessel cool down/warm up, and environment heat transfer. Based on experimental data measured at LLNL on the Linde 875 bar LH 2 cryo-pump and results from the model, and extrapolating for refueling stations of various sizes, it can be shown that boil-off losses can vary from 15% of delivered LH 2 for a 100 kg/day, 5% at 400 kg/day, and down to less than 2% for stations above 1,800 kg/day (for a system similar to the one at LLNL). Less boil-off is to be expected for LH 2 pumps dispensing at 350 bar, so that less than 0.7% can be expected above 1,800 kg/day station capacities. At last, boil-off recovery solutions are briefly analyzed, including compressors, cryo-coolers and stationary fuel cell. Their economic value is strongly related to the costs of H2 and of the electricity needed to perform the recovery, on top of the initial capital expenditures. Effective boil-off could be reduced to 0.6-1% for a 2,000 kg/day station if adequate recovery solutions are used.

08 HYDROGEN↗