Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid protection”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Evaluation of Modeling NO2 Concentrations Driven by Satellite-Derived and Bottom-Up Emission Inventories Using In-Situ Measurements Over China

Chemical transport models together with emission inventories are widely used to simulate NO2 concentrations over China, but validation of the simulations with in situ measurements has been extremely limited. Here we use ground measurements obtained from the air quality monitoring network recently developed by the Ministry of Environmental Protection of China to validate modeling surface NO2 concentrations from the CHIMERE regional chemical transport model driven by the satellite-derived DECSO and the bottom-up MIX emission inventories. We applied a correction factor to the observations to account for the interferences of other oxidized nitrogen compounds (NOz), based on the modeled ratio of NO2 to NOz. The model accurately reproduces the spatial variability in NO2 from in situ measurements, with a spatial correlation coefficient of over 0.7 for simulations based on both inventories. A negative and positive bias is found for the simulation with the DECSO (slopeD0.74 and 0.64 for the daily mean and daytime only) and the MIX (slopeD1.3 and 1.1) inventories, respectively, suggesting an underestimation and overestimation of NOx emissions from corresponding inventories. The bias between observed and modeled concentrations is reduced, with the slope dropping from 1.3 to 1.0 when the spatial distribution of NOx emissions in the DECSO inventory is applied as the spatial proxy for the MIX inventory, which suggests an improvement of the distribution of emissions between urban and suburban or rural areas in the DECSO inventory compared to that used in the bottom-up inventory. A rough estimate indicates that the observed concentrations, from sites predominantly placed in the populated urban areas, may be 10-40% higher than the corresponding model grid cell mean. This reduces the estimate of the negative bias of the DECSO-based simulation to the range of -30 to 0% on average and more firmly establishes that the MIX inventory is biased high over major cities. The performance of the model is comparable over seasons, with a slightly worse spatial correlation in summer due to the difficulties in resolving the more active NOx photochemistry and larger concentration gradients in summer by the model. In addition, the model well captures the daytime diurnal cycle but shows more significant disagreement between simulations and measurements during nighttime, which likely produces a positive model bias of about 15% in the daily mean concentrations. This is most likely related to the uncertainty in vertical mixing in the model at night.

NO2 modelling↗
DOE OSTI2025

Solid State Transformer Architecture and Control Compensation for Common Mode Currents

A high-altitude electromagnetic pulse (HEMP) or similar geomagnetic disturbance (GMD) has the potential to impact the operation of large-scale electric power grids. By introducing low-frequency common-mode (CM) currents, these events can degrade the performance of critical system components, such as large power transformers by introducing CM currents which can lead to magnetic saturation of the transformer core. In this work, a solid-state transformer (SST) is developed to replace susceptible equipment and improve grid resiliency by safely absorbing these CM disturbances. This device will be referred to as a common-mode solid-state transformer (CM-SST). An SST architecture based on a four-legged AC/DC converter is developed. This architecture enables active control of CM signals without disturbing the AC voltages or the real and reactive power delivery capabilities. A system-level model of this architecture is created, and time-domain simulations are performed to evaluate the SST’s performance in response to simulated CM disturbances. A control strategy for mitigating CM current is also investigated. Hamiltonian surface shaping and power flow control (HSSPFC) is used to design a nonlinear controller for the SST’s output inverter. The objectives of the controller are to suppress CM-induced AC current offsets and regulate AC currents to desired setpoints. Nonlinear system analysis is applied to design and validate the controller. Two cases are tested: (a) the proposed four-leg inverter and (b) a standard three-leg inverter. The results show that the proposed controller rapidly mitigates CM disturbances while maintaining high-quality AC current waveforms in the four-leg configuration. Finally, the hardware performance of an SST prototype is evaluated. In particular, the ability of the SST to safely redirect and absorb CM currents is demonstrated, showing how it can protect neighboring conventional transformers in the system. The study confirms that appropriate control laws allow the SST to protect both itself and adjacent transformers during a HEMP or GMD event.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Securing Grid-interactive Efficient Buildings (GEB) through Cyber Defense and Resilient System (CYDRES)

The DOE CYDRES project is driven by the urgent need to address critical research gaps in the domain of cyber-physical security of smart buildings, including Grid-interactive Efficient Buildings (GEBs). CYDRES, a real-time advanced building resilient platform, aims to enhance the cyber-attack-immune capabilities of buildings through multi-layered prevention, detection, and adaptation mechanisms. CYDRES consists of five key modules: a multi-layer network analyzer, an Automatic Fault Detection, Diagnosis, and Prognosis (AFDDP) framework, an intelligent mode selector, a cyber-resilient control framework, and a situation awareness platform. The Network Analyzer employs a data-driven framework that includes a protocol state learning tool and a CRF (Conditional Random Field) command validator. In Hardware-In-the-Loop (HIL) testbeds, it achieved 100% detection accuracy with a false alarm rate of 3%, validating its efficacy in identifying selected cyber-attacks. The AFDDP framework leverages pattern matching, PCA (Principal Component Analysis)-based strategies, and a DBN (Dynamic Bayesian Network)-based fault diagnosis approach to pinpoint the causes of physical system abnormalities using Building Automation System (BAS) data. In HIL experiments, the AFDDP module attained a detection accuracy of over 95% with a false alarm rate below 7%. Additionally, the fault detector utilized machine learning (Random Forest) and deep learning (Multi-Layer Perceptron) methods with acoustic sensor data to achieve a 100% fault detection accuracy in Heating, Ventilation, and Air-Conditioning (HVAC) equipment. The Mode Selector offered real-time impact analysis, allowing immediate actions to protect BASs in the face of emerging threats. The cyber-resilient control framework included an adaptive Model Predictive Control (MPC) and a measurement compensator, reducing temperature violations by up to 94% and improving the total demand flexibility by up to 70% in HIL experiments. Such HIL experiments covered a cyber-attack case and a physical fault case, showcasing CYDRES’ efficiency in maintaining operational continuity during threats. The situation awareness platform in Grafana enhanced real-time threat detection and response visualization, augmenting the operational awareness for building operators. CYDRES demonstrated high technical effectiveness in various test scenarios, particularly in HIL environments. The project's phased development approach ensured efficient use of resources, highlighting its practical feasibility and readiness for commercialization. By enhancing the security and resilience of building operations, CYDRES represents a significant advance in mitigating risks associated with cyber-physical systems, thereby enhancing public confidence in the safety of modern building infrastructure. Future directions for the project include expanding testing protocols, refining AFDDP methodologies, exploring more comprehensive resilient control strategies, and testing in real commercial buildings.

42 ENGINEERING↗
DOE OSTI2025

A census of fish passage facilities at US hydropower developments across the conterminous United States

Hydropower provides reliable and secure electricity and contributes significantly to the flexibility and stability of the US electricity grid. Hydropower dams generally are aquatic barriers and hazards to migratory fish in rivers. Fish passage facilities mitigate risks from hydropower to migratory fish, but information on these facilities is incomplete for the conterminous United States (CONUS). Here, we present the first CONUS-scale dataset of fish passage facilities at US hydropower developments in over 30 years. The existence of fish passage facilities (presence or absence) was specified for 1909 hydropower features, 390 of which had at least one facility. Most features had a single passage facility that provided passage in only one direction, with downstream being more common than upstream passage. Bypasses and ladders accounted for 60 % of the fish passage facility types. While we documented 659 fish passage facilities across CONUS, facilities were most common in the New England, Pacific Northwest, Mid-Atlantic, and Great Lakes regions. In general, fish passage facilities were more common at features that were closer to the ocean, at lower elevations, and at shorter dams, but not related to installed electrical generation capacity. Hydrologic sub-basins containing salmonids also contained the largest number of hydropower features, but the proportion of features with passage was generally higher in sub-basins containing multiple migratory taxa. This census provides valuable information on existing fish passage mitigation and is a benchmark to gauge progress toward a modernized hydropower fleet that provides affordable, reliable energy while protecting fishery resources and river ecosystems.

13 HYDRO ENERGY↗
DOE OSTI2025

Package Data for CERF-Data Centers

This dataset contains sample input 100m resolution raster files for running the CERF-DC python package (see https://github.com/IMMM-SFA/cerf_data_centers) at the state level across the CONUS. Due to data availability constraints, some of the items included in this dataset are proxies or assumptions for siting factors used in the model. These are individually noted in the item descriptions and can be exchanged with more detailed information upon availability. Data Descriptions The following raster files are included in the data download: state_siting_region.tif — State areas identified by state FIPS code composite_siting_suitability.tif — Value of 1 indicates suitable siting location, 0 otherwise. The following areas are excluded from siting: Areas within 300m of a federal airport runway Waterbodies Areas with slope >16% Areas susceptible to sinkholes High coastal or inland flood risk areas Local, state, and federal parks, leisure areas, and cemeteries Areas >2 km away from electric substations Areas >5 km away from a municipal water supplier service area Areas >2 km away from high-speed fiber provider service territory Protected Areas Database of the United States (PAD-US) areas Railroads, major roadways, and minor roadways Military areas and training grounds Developed lands Areas >0.8 km (0.5 miles) from developed lands land_value_dollar_per_sqft.tif — USD per square foot (sqft) derived from USDA $/acre land cost personal_property_tax_rate.tif — Personal property tax rate by state. Uses an assumed 0.0125 personal property tax rate for states with personal property tax, 0 for states without personal property tax. real_property_tax_rate.tif — Real property tax rate. Based on county level residential real estate property tax rates. sales_tax_rate.tif — Sales tax rate by state. mechanical_cooling_fraction.tif — Fraction of year (values between 0 and 1, inclusive) that the data center would be cooled through mechanical processes based on local water stress and humidity levels. water_cooling_fraction.tif — Fraction of year (values between 0 and 1, inclusive) that the data center would be cooled through evaporative (water cooled) processes based on local water stress and humidity levels. distance_to_substation.tif — Distance to nearest substation in hundreds of meters (i.e., value of 1 equals a distance of 100m). Offshore areas have a value of 0. industrial_electricity_rates_dollar_per_kwh.tif — USD/kWh industrial electricity rates. Represents the average industrial rate across all utilities that operate within a given county. Values are derived from the US Utility Rate Database. commercial_electricity_rates_dollar_per_kwh.tif — USD/kWh commercial electricity rates. Represents the average commercial rate across all utilities that operate within a given county. Values are derived from the US Utility Rate Database. data_center_market_locations.tif — Grid cells with positive values represent the centroid of existing data center market clusters. The value of non-zero grid cells represents the number of data centers in the market cluster. All other grid cells have a value of 0. Geospatial Metadata CRS: Albers Equal Area Conic (ESRI:102003) Extent: -2415585.0000000023283064,-1441981.2605773280374706 : 2384414.9999999976716936,1708018.7394226719625294 Dimensions: X: 48000 Y: 31500 Bands: 1 Origin: -2415585.0000000023283064,1708018.7394226719625294 Pixel Size: 100,-100 Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program. License This data is made available under a CCBY4.0 License Disclaimer This material was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor the United States Department of Energy, nor the Contractor, nor any or their employees, nor any jurisdiction or organization that has cooperated in the development of these materials, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness or any information, apparatus, product, software, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or Battelle Memorial Institute. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. PACIFIC NORTHWEST NATIONAL LABORATORYoperated byBATTELLEfor theUNITED STATES DEPARTMENT OF ENERGYunder Contract DE-AC05-76RL01830

Mongird, Kendall↗
DOE OSTI2021

Specific absorbed fractions for a revised series of the UF/NCI pediatric reference phantoms: internal photon sources

Assessment of radiation absorbed dose to internal organs of the body from the intake of radionuclides, or in the medical setting through the injection of radiopharmaceuticals, is generally performed based upon reference biokinetic models or patient imaging data, respectively. Biokinetic models estimate the time course of activity localized to source organs. The time-integration of these organ activity profiles are then scaled by the radionuclide S-value, which defines the absorbed dose to a target tissue per nuclear transformation in various source tissues. S-values are computed using established nuclear decay information (particle energies and yields), and a parameter termed the specific absorbed fraction (SAF). The SAF is the ratio of the absorbed fraction—fraction of particle energy emitted in the source tissue that is deposited in the target tissue—and the target organ mass. While values of the SAF may be computed using patient-specific or individual-specific anatomic models, they have been more widely available through the use of computational reference phantoms. In this study, we report on an extensive series of photon SAFs computed in a revised series of the University of Florida and the National Cancer Institute pediatric reference phantoms which have been modified to conform to the specifications embodied in the ICRP reference adult phantoms of Publication 110 (e.g. organs modeled, organ ID numbers, blood contribution to elemental compositions). Following phantom anatomical revisions, photon radiation transport simulations were performed using MCNPX v2.7 in each of the ten phantoms of the series—male and female newborn, 1 year old, 5 year old, 10 year old, and 15 year old—for 60 different tissues serving as source and/or target regions. A total of 25 photon energies were considered from 10 keV to 10 MeV along a logarithm energy grid. Detailed analyses were conducted of the relative statistical errors in the Monte Carlo target tissue energy deposition tallies at low photon energies and over all energies for source–target combinations at large intra-organ separation distances. Based on these analyses, various data smoothing algorithms were employed, including multi-point weighted data smoothing, and log–log interpolation at low energies (1 keV and 5 keV) using limiting SAF values based upon target organ mass to bound the interpolation interval. The final dataset is provided in a series of ten electronic supplemental files in MS Excel format. Finally, the results of this study were further used as the basis for assessing the radiative component of internal electron source SAFs as described in our companion paper (Schwarz et al 2021) for this same pediatric phantom series.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS2019

Lightning Strike Distance Distribution Beyond a Preexisting Lightning Area

The 45th Weather Squadron (45 WS) asked the Applied Meteorology Unit (AMU) to review the 30-year-old, lightning stand-off distances of 5 nautical miles (nmi) for applicability to today's operations. This was based on the realization that previous lightning strike distance studies did not match how 45 WS issues lightning warnings (Roeder, 2008). The previous lightning distance studies were from the point of origin of the lightning or the average starting location that would tend to be in the core of the thunderstorm. However, the 45 WS issues lightning warnings based on the edge of a preexisting lightning area. Before beginning the AMU project, it took several years to develop a method to calculate a distance distribution beyond a preexisting area (Roeder, 2015). The AMU pulled Lightning Detection and Ranging (LDAR) sensor data from 1/1/2013 to 12/31/2013. This dataset consisted of 37 million individual source data points from the LDAR sensors. Only sources within 50 km north, south, east or west of the LDAR grid center were included in the dataset. This limited the use of LDAR data to that with the greatest accuracy of source detection and increased data processing. Points were grouped into flashes based on spatial and temporal criteria. Based on the sensitivity analysis the AMU performed on the flash clustering algorithm, a time value of 0.3 seconds was found to model flashes adequately. Distance parameters were tested from 1,500 to 7,500 meters (m) in 500 m increments. Distance parameters of both 3,000 m and 4,000 m produced results in the plotting tool that were most representative of the physical behavior of lightning. Thus statistics were gathered for the most representative of these spatial and temporal criteria on the flash size and the polygon expansion distance in order to find the correct data distributions. The best fit curves for the LDAR polygon expansion frequency vs. distances for both the 3 kilometer (km) and 4 km distance threshold values were exponential decay functions and had R2 values of > 0.998, indicating good model fits. The equations of the best fit curves were then used to calculate a desired safety radius of 4 nmi for either 3 km or 4 km distance threshold criteria. The AMU analysis concludes the safe reduction of the 5 nmi lightning warning circles to 4 nmi should improve the operational impact by 36% if based on distance from the center of the property area being protected. If based on the edge of the property being protected, then the reduction is 4.5 nmi to 4 nmi and the operational impact is decreased by 21%. For the 6 nmi lightning warning circles, the recommended 4 nmi stand-off distance will result in a safe reduction of operational impact of 31% if based on the center of the area being protected, or 16% if based on the edge of the property.

Flash clustering algorithm↗
DOE OSTI2026

Decentralised Reinforcement Learning for Dynamic Cyberattack Response in Microgrid Networks

Microgrids rely on communication networks for reliable operation, which makes them inherently vulnerable to cyberattacks. Such attacks can destabilise system dynamics and drive states away from their nominal operating trajectories. Although several physics-informed and machine learning-based strategies have been developed to counter these threats, the rapidly evolving cyber landscape enables adversaries to bypass static defences or rules-based mitigation approaches. This paper proposes a dynamic, online-trained and fully decentralised reinforcement learning (RL)-based cyberattack response framework to protect microgrids from evolving cyberattacks. The proposed framework deploys multiple deep Q-networks (DQNs), each associated with a distributed energy resource (DER), to enable localised and adaptive attack mitigation. In this framework, each DQN processes local voltage and frequency measurements—combined with intrusion detection system (IDS) alerts—as observations and rewards to guide decision-making. Extensive simulation studies demonstrate the robustness of the proposed framework under diverse attack scenarios and varying IDS-induced detection delays. Comparative analysis highlights its superiority over existing static or preexisting rules-based mitigation approaches. Finally, we present an analysis that shows the framework's scalability to real-life microgrids with more interacting agents.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2014

Simulation of Fluid Flow and Collection Efficiency for an SEA Multi-element Probe

Numerical simulations of fluid flow and collection efficiency for a Science Engineering Associates (SEA) multi-element probe are presented. Simulation of the flow field was produced using the Glenn-HT Navier-Stokes solver. Three-dimensional unsteady results were produced and then time averaged for the heat transfer and collection efficiency results. Three grid densities were investigated to enable an assessment of grid dependence. Simulations were completed for free stream velocities ranging from 85-135 meters per second, and free stream total pressure of 44.8 and 93.1 kilopascals (6.5 and 13.5 pounds per square inch absolute). In addition, the effect of angle of attack and yaw were investigated by including 5 degree deviations from straight for one of the flow conditions. All but one of the cases simulated a probe in isolation (i.e. in a very large domain without any support strut). One case is included which represents a probe mounted on a support strut within a finite sized wind tunnel. Collection efficiencies were generated, using the LEWICE3D code, for four spherical particle sizes, 100, 50, 20, and 5 micron in diameter. It was observed that a reduction in velocity of about 20% occurred, for all cases, as the flow entered the shroud of the probe. The reduction in velocity within the shroud is not indicative of any error in the probe measurement accuracy. Heat transfer results are presented which agree quite well with a correlation for the circular cross section heated elements. Collection efficiency results indicate a reduction in collection efficiency as particle size is reduced. The reduction with particle size is expected, however, the results tended to be lower than the previous results generated for isolated two-dimensional elements. The deviation from the two-dimensional results is more pronounced for the smaller particles and is likely due to the reduced flow within the protective shroud. As particle size increases differences between the two-dimensional and three dimensional results become negligible. Taken as a group, the total collection efficiency of the elements including the effects of the shroud has been shown to be in the range of 0.93 to 0.99 for particles above 20 microns. The 3D model has improved the estimated collection efficiency for smaller particles where errors in previous estimates were more significant.

Water Content↗
NASA NTRS2017

Tortuosity Computations of Porous Materials using the Direct Simulation Monte Carlo

Low-density carbon fiber preforms, used as thermal protection systems (TPS) materials for planetary entry systems, have permeable, highly porous microstructures consisting of interlaced fibers. Internal gas transport in TPS is important in modeling the penetration of hot boundary-layer gases and the in-depth transport of pyrolysis and ablation products. The gas effective diffusion coefficient of a porous material must be known before the gas transport can be modeled in material response solvers; however, there are very little available data for rigid fibrous insulators used in heritage TPS.The tortuosity factor, which reflects the efficiency of the percolation paths, can be computed from the effective diffusion coefficient of a gas inside a porous material and is based on the micro-structure of the material. It is well known, that the tortuosity factor is a strong function of the Knudsen number. Due to the small characteristic scales of porous media used in TPS applications (typical pore size of the order of 50 micron), the transport of gases can occur in the rarefied and transitional regimes, at Knudsen numbers above 1. A proper way to model the gas dynamics at these conditions consists in solving the Boltzmann equation using particle-based methods that account for movement and collisions of atoms and molecules.In this work we adopt, for the first time, the Direct Simulation Monte Carlo (DSMC) method to compute the tortuosity factor of fibrous media in the rarefied regime. To enable realistic simulations of the actual transport of gases in the porous medium, digitized computational grids are obtained from X-ray micro-tomography imaging of real TPS materials. The SPARTA DSMC solver is used for simulations. Effective diffusion coefficients and tortuosity factors are obtained by computing the mean-square displacement of diffusing particles.We first apply the method to compute the tortuosity factors as a function of the Knudsen number for computationally designed materials such as random cylindrical fibers and packed bed of spheres with prescribed porosity. Results are compared to literature values obtained using random walk methods in the rarefied and transitional regime and a finite-volume method for the continuum regime. We then compute tortuosity factors for a real carbon fiber material with a transverse isotropic structure (FiberForm), quantifying differences between through-thickness and in-plain tortuosities at various Knudsen regimes.

Tortuosity↗
DOE OSTI2026

Package Data for CERF-Data Centers

This dataset contains sample input 100m resolution raster files for running the CERF-DC python package (see https://github.com/IMMM-SFA/cerf_data_centers) at the state level across the CONUS. Due to data availability constraints, some of the items included in this dataset are proxies or assumptions for siting factors used in the model. These are individually noted in the item descriptions and can be exchanged with more detailed information upon availability. Data Descriptions The following raster files are included in the data download: state_siting_region.tif — State areas identified by state FIPS code composite_siting_suitability.tif — Value of 1 indicates suitable siting location, 0 otherwise. The following areas are excluded from siting: Areas within 300 m of a federal airport runway or within an airport area boundary Waterbodies Areas with slope >16% Areas susceptible to sinkholes High coastal or inland flood risk areas Local, state, and federal parks, leisure areas, and cemeteries Areas >2 km away from electric substations Areas >5 km away from a municipal water supplier service area Areas >2 km away from high-speed fiber provider service territory USGS Protected Areas Database of the United States (PAD-US) GAP status 1, 2, or 3 areas US National Parks Wetlands USFWS critical habitats BIA land areas Railroads, major roadways, and minor roadways Military areas and training grounds NLCD developed lands Areas >0.8 km (0.5 miles) from NLCD developed lands land_value_dollar_per_sqft.tif — USD per square foot (sqft) derived from USDA $/acre land cost personal_property_tax_rate.tif — Personal property tax rate by state. Uses an assumed 0.0125 personal property tax rate for states with personal property tax, 0 for states without personal property tax. real_property_tax_rate.tif — Real property tax rate. Based on county level residential real estate property tax rates. sales_tax_rate.tif — Sales tax rate by state. mechanical_cooling_fraction.tif — Fraction of year (values between 0 and 1, inclusive) that the data center would be cooled through mechanical processes based on local water stress and humidity levels. water_cooling_fraction.tif — Fraction of year (values between 0 and 1, inclusive) that the data center would be cooled through evaporative (water cooled) processes based on local water stress and humidity levels. distance_to_substation.tif — Distance to nearest substation in hundreds of meters (i.e., value of 1 equals a distance of 100m). Offshore areas have a value of 0. industrial_electricity_rates_dollar_per_kwh.tif — USD/kWh industrial electricity rates. Represents the average industrial rate across all utilities that operate within a given county. Values are derived from the US Utility Rate Database. commercial_electricity_rates_dollar_per_kwh.tif — USD/kWh commercial electricity rates. Represents the average commercial rate across all utilities that operate within a given county. Values are derived from the US Utility Rate Database. data_center_market_locations.tif — Grid cells with positive values represent the centroid of existing data center market clusters. The value of non-zero grid cells represents the number of data centers in the market cluster. All other grid cells have a value of 0. Geospatial Metadata CRS: Albers Equal Area Conic (ESRI:102003) Extent: -2415585.0000000023283064,-1441981.2605773280374706 : 2384414.9999999976716936,1708018.7394226719625294 Dimensions: X: 48000 Y: 31500 Bands: 1 Origin: -2415585.0000000023283064,1708018.7394226719625294 Pixel Size: 100,-100 Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program. License This data is made available under a CCBY4.0 License Disclaimer This material was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor the United States Department of Energy, nor the Contractor, nor any or their employees, nor any jurisdiction or organization that has cooperated in the development of these materials, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness or any information, apparatus, product, software, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or Battelle Memorial Institute. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. PACIFIC NORTHWEST NATIONAL LABORATORYoperated byBATTELLEfor theUNITED STATES DEPARTMENT OF ENERGYunder Contract DE-AC05-76RL01830

Mongird, Kendall↗
DOE OSTI2023

Simulation of compound flooding in Japan using a nationwide model

For this study, a high-resolution, summit-to-sea unstructured-grid model was used to simulate two compound flooding events in different geomorphic settings in Japan: the July 2012 flood in Kumamoto city and July 2018 flood in Okayama City, both caused by a torrential rainfall during seasonal Meiyu front but otherwise exhibiting quite different compound flood characteristics. The model shows good performance in simulating flooding extent; e.g., the Hit Rate for the inundation event in 2018 exceeds 0.9. Sensitivity tests were conducted to determine the roles/significance of levee structures along major rivers in coastal regions. Our results indicate that riverine and oceanic factors as well as levees significantly contribute to the observed flooding extents, and we are able to quantify the contributions from each factor at different locations. The model has been implemented operationally as the first nationwide 3-day compound flooding forecast system for Japan.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1992

A High Angle of Attack Inviscid Shuttle Orbiter Computation

An accurate description of the aerothermal environment is required to minimize the weight of the Thermal Protection System required on the leeside of winged reentry vehicles. The inability of ground-test facilities to reproduce the high enthalpy, separated flow present during reentry flight conditions, coupled with the prohibitive expense of flight tests, leads to the use of an analytical method - namely Computational Fluid Dynamics (CFD) - to describe the flow. While the ultimate goal of this work is to accurately predict the leeside flow and its associated thermal environment, an essential and reasonable first step towards that goal is a comparison of pressure predictions by a code with wind-tunnel data. Until such CFD pressure predictions agree with wind-tunnel test cases, there is little hope of accurately predicting the thermal environment at flight conditions. Thus, the objective of this study is to compare the pressure distributions predicted by inviscid, perfect gas CFD to Shuttle Orbiter wind-tunnel data and to address any significant issues encountered during the computation. While flight data is available for the Shuttle Orbiter, a wind-tunnel case is chosen for this study to allow a tractable problem for preliminary investigation. A wind-tunnel case allows the perfect gas assumption for the flow chemistry. This provides a significant computational savings over a several species finite-rate chemistry model which would be necessary if high-temperature effects present at flight conditions were to be included. In addition, by concentrating on the surface pressures, the analysis need only consider inviscid flow for general evaluation of the code capability. This further reduces the computational expense due to the absence of viscous terms and the associated decrease in the number of points required for the computational grid. Previous computational efforts (such as STEIN and HALIS) have been directed toward the windward surface quantities, primarily due to restrictions in treating either the winged geometry or its associated subsonic regions at high angle of attack. The code used for this study, the LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) code of Gnoffo, represents a state-of-the-art code for computing the flow over complex configurations at hypersonic speeds. In the study, the LAURA code is applied to a wind-tunnel condition to initiate the assessment of the code's ability to predict the flow over a relatively complex hypersonic vehicle at high angles of attack. This presentation is used to highlight the pertinent results of a more detailed investigation of the inviscid calculation over the Shuttle Orbiter with the LAURA code.

William L. Kleb↗
NASA NTRS2008

Classification of Small Negative Lightning Reports at the KSC-ER

The NASA Kennedy Space Center (KSC) and Air Force Eastern Range (ER) operate an extensive suite of lightning sensors because Florida experiences the highest area density of ground strikes in the United States, with area densities approaching 16 fl/sq km/yr when accumulated in 10x10 km (100 sq km) grids. The KSC-ER use data derived from two cloud-to-ground (CG) lightning detection networks, the "Cloud-to-Ground Lightning Surveillance System" (CGLSS) and the U.S. National Lightning Detection Network (TradeMark) (NLDN) plus a 3-dimensional lightning mapping system, the Lightning Detection and Ranging (LDAR) system, to provide warnings for ground operations and to insure mission safety during space launches. For operational applications at the KSC-ER it is important to understand the performance of each lightning detection system in considerable detail. In this work we examine a specific subset of the CGLSS stroke reports that have low values of the negative inferred peak current, Ip, i.e. values between 0 and -7 kA, and were thought to produce a new ground contact (NGC). When possible, the NLDN and LDAR systems were used to validate the CGLSS classification and to determine how many of these reported strokes were first strokes, subsequent strokes in a pre-existing channel (PEC), or cloud pulses that the CGLSS misclassified as CG strokes. It is scientifically important to determine the smallest current that can reach the ground either in the form of a first stroke or by way of a subsequent stroke that creates a new ground contact. In Biagi et al (2007), 52 low amplitude, negative return strokes ([Ip] < or = 10 kA) were evaluated in southern Arizona, northern Texas, and southern Oklahoma. The authors found that 50-87% of the small NLDN reports could be classified as CG (either first or subsequent strokes) on the basis of video and waveform recordings. Low amplitude return strokes are interesting because they are usually difficult to detect, and they are thought to bypass conventional lightning protection that relies on a sufficient attractive radius to prevent "shielding failure" (Golde, 1977). They also have larger location errors compared to the larger current events. In this study, we use the estimated peak current provided by the CGLSS and the results of our classification to determine the minimum Ip for each category of CG stroke and its probability of occurrence. Where possible, these results are compared to the findings in the literature.

Ward, Jennifer G.↗
DOE OSTI2022

Large-scale Hydrogen Storage – Risk Assessment Seattle City Light and Port of Seattle [Abstract]

This CRADA presents the strategy that Pacific Northwest National Laboratory (PNNL) and Sandia National Laboratories (SNL) will take to support Seattle City Light (SCL), and the Port of Seattle (Port) in performing a risk assessment of large-scale hydrogen storage. Risk assessment is often used to ensure that adequate measures are taken to protect workers and the public, the environment, infrastructure, and assets. A detailed risk assessment can also be used to direct funding and upgrades, to specific components and sub-systems in order to mitigate risks to the larger system. In this way risk assessments are often employed as a part of a larger risk management strategy, with the goal of minimizing the occurrence of hazards and to identify means to limit their consequences. Risk assessment is often used to engage and inform regulators, and to communicate how specific regulations are being met. However, it is important to note that the proposed work is not intended for SCL and the Port to use in order to gain regulatory acceptance for their proposed activities. The work performed as a part of this effort will be a preliminary risk assessment for early-stage component and system designs and should be considered research and development (R&D). As such, the proposed work will be performed to a quality level and design maturity consistent with R&D and is not considered appropriate for final safety analysis and regulatory compliance purposes. Previous and on-going work at SCL and the Port demonstrated the utility of deploying hydrogen systems at the Port. The deployment of hydrogen at the Port is a part of a larger vision of using hydrogen to address a range of issues for SCL and the Port. These include large-scale fueling of MD/HD vehicles, cargo-handling equipment (CHE), and harbor vessels to reduce emissions; support of adjacent LD vehicles; support of critical port operations during extreme events (i.e., resiliency); deferral of more capital- and time-intensive electrical distribution system upgrades while still supporting evolving port operations and decarbonization efforts; facilitating electrification by establishing energy storage as a grid resource, starting at strategic port locations; creation of a flexible market resource that can be used by SCL to generate revenue via arbitrage; support of planned future maritime operations that involve heavy use of hydrogen for ocean-going vessels; and future end-use applications involving natural gas pipeline hydrogen injection. Ultimately, the success of these activities is underpinned by the deployed storage capacity. Large-scale deployment of hydrogen systems will require hydrogen storage at a scale that has not been demonstrated. In addition, the ideal location for such multi-use systems is near the end user which will often necessitate deploying into urban and/or industrial areas. A detailed risk assessment using the Port as a test case is necessary to ensure the deployment of large-scale hydrogen is successful. Many technologies have been proposed for hydrogen storage; however, these technologies need to be analyzed as they apply to an actual site. The physical infrastructure and hydrogen use cases for the Port will be analyzed, and a risk assessment for compressed hydrogen, liquified hydrogen, and Liquid Organic Hydrogen Carrier (LOHC) storage will be performed. These risk assessments will be useful for understanding how each of these technologies would perform in terms of facility and public safety. The operating states of the proposed hydrogen systems at the Port will be analyzed and incorporated into the storage risk assessment. Scalability will also be analyzed to understand how future port uses would affect the overall risk assessment. Finally, using the risk assessment as a tool to inform engagement and to gain stakeholder acceptance will be explored.

08 HYDROGEN↗
DOE OSTI2023

A review of privacy in energy applications

As the distribution system continues to experience an increase in distributed energy resource (DER) and electric vehicle (EV) penetration, so does the need for new solutions that can help grid operators manage and leverage their capabilities. This will undoubtedly lead to new operational schemes and business opportunities that will transform the traditional consumer into a prosumer who will be more actively engaged in grid operations. Although the field is still under active development, many of the potential use cases presented in literature or industry are built upon edge computing, two-way communications, and other innovative computational constructs to attain their goals. However, at their core, many use cases assume a great level of data access to aid with the decision-making process, an assumption that may need to be revised to ensure fair and equitable operational processes are maintained. This may be particularly true as edge resources are predicted to participate in retail-side, many-to-many, or peer-to-peer markets and thus may lead to financial impacts if data access considerations are ignored. The need to revise data access mechanisms can be further justified by the introduction of new participants into the operational process, who do not have the same level of trust, nor the incentives to focus on energy delivery as their primary objective. At the same time, more consumers are becoming aware of their own data, and the potential impacts of its abuse. To help solution developers better understand these risks, this report has been developed to offer an initial introduction to the topic of privacy. This is achieved by 1) Highlighting the need for privacy-aware solutions; 2) Encouraging system designers to be inquisitive about the status quo; 3) Documenting the existing threat space; 4) Presenting and evaluating tools that may be helpful towards enabling better privacy postures; and 5) Making recommendations to encourage the adoption of better practices. From a technical perspective, the report focuses on evaluating two potential techniques by applying them to the Transactive Energy Space. Based on the obtained results, it can be established that differential privacy (DP) methods may have limited applicability when highly correlated, time-series data records need to be protected. However, DP may be a powerful tool when it is used to aggregate and analyze mid-size and large-size data sets in a more traditional statistical environment. The second tool under evaluation is threshold cryptography, which can guarantee complete secrecy (and thus privacy) but requires the establishment of key management procedures and dedicated communication channels for key coordination. Therefore, due to its increased computational overhead, the use of threshold cryptography must be weighted using a cost/benefit analysis on a per-application basis.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Characterization of MC&A for the Molten Salt Fuel Cycle

Advanced reactor developers are exploring diverse reactor designs, including molten salt reactors (MSRs). These advanced reactors are considered for wider applications and a range of deployment locations, including supporting the integration of renewable energy sources in the grid. There are three main types of MSRs: (1) reactors in which the fuel salt freely circulates within the core; (2) reactors with the fuel salt contained within vented fuel tubes; and (3) reactors that use molten salt solely as a coolant, with the fuel in a separate, solid form. In this document, the term MSR refers specifically to the first two types, which use fuel salt—special nuclear material (enriched uranium, plutonium, and 233 U) in chloride or fluoride form mixed with chloride- or fluoride-based carrier salt in a peritectic mixture—as the primary medium for fission. The composition of fuel salt, both at startup and for makeup or refueling, varies depending on the MSR design and the chosen fuel cycle approach, which can be either once-through or closed. For MSRs, a variety of fuel cycle approaches (e.g., U, U–Pu, U–Pu–TRU, U–Th, U–Pu–Th) are being considered. Fuel in MSRs is much different than traditional solid fuel, including its preparation. The uniqueness warrants investigation into characterizing fuel preparation processes, known as fuel salt synthesis . This effort characterized major fuel preparation and synthesis processes, identifying temperature, equipment, and environmental requirements for uranium-, plutonium-, and thorium-based fuel preparation and synthesis. Because MSR fuel salt synthesis facilities handle special nuclear material in loose, bulk form, a material control and accounting plan will be required for licensing from the US Nuclear Regulatory Commission or under the US Department of Energy authorization. This effort serves as a foundation to investigate material control and accounting approaches for synthesis facilities, including determining measurement points and techniques. Because several MSR developers are planning demonstration facilities in the coming years, this effort will support stakeholders with preparing or reviewing material control and accounting plans for providing assurance that all special nuclear material is accounted for at fuel salt synthesis facilities. This report was produced for Materials Protection, Accounting, and Control Technologies (MPACT) program under the US Department of Energy (DOE), Office of Nuclear Energy, Nuclear Fuel Cycle and Supply Chain.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRSDate not supplied

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd [1] and associated criteria (KnGLL > 0.05) can be accurately used to predict zones of breakdown in the flow, and associated errors for surface quantities. Furthermore, when studying reacting flow, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown in minimal (see Fig. 1). Finally, we studied the impact of the grid refinement on the surface quantities in DSMC.

DSMC↗