SEARCH · Search NASA

Results for “Distribution System 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Improving Trustworthiness of Data-Driven Power Grid Contingency Analysis With Bayesian Residual Graph Neural Networks

The evolving energy landscape requires novel tools to efficiently perform contingency analysis and reliability assessment of power grids, potentially in real-time. The high computational cost of traditional power flow solvers limits their applicability in practice. Machine learning (ML) surrogates such as deep neural networks (NNs) accelerate power flow solvers computations, enabling high-order contingency analysis and real-time decision-making by learning highly nonlinear functions and integrating grid topology via graph architectures. However, (graph) NNs lack predictive power away from training data and do not provide predictive confidence estimates. Here, we present a Bayesian residual graph NN that integrates knowledge from low-fidelity data via residual training and embeds granular quantification of uncertainties, improving trustworthiness critical for high-consequence decision-making. Applying Bayesian concepts to NNs is challenging due to the high-dimensionality of both the parameter space, complicating derivation of a meaningful prior, and the output space in large grid systems, requiring enhanced techniques to assess the predicted high-dimensional uncertainties. Our contributions include: (1) Deriving a prior for fully connected and graph NNs that leverages low-fidelity data to guide mean predictions and appropriately control prior predictive uncertainty. (2) Integrating this prior within an ensembling with anchoring scheme for efficient approximate posterior inference. (3) Deriving enhanced metrics to assess accuracy of both the mean and uncertainty predictions in high dimensions, appropriately accounting for correlations propagated through graph layers. The resulting Bayesian residual graph NN is tested on a contingency analysis task for 14-bus and 118-bus grids.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2015

MONTE CARLO NEUTRINO TRANSPORT THROUGH REMNANT DISKS FROM NEUTRON STAR MERGERS

We present Sedonu, a new open source, steady-state, special relativistic Monte Carlo (MC) neutrino transport code, available at bitbucket.org/srichers/sedonu. The code calculates the energy- and angle-dependent neutrino distribution function on fluid backgrounds of any number of spatial dimensions, calculates the rates of change of fluid internal energy and electron fraction, and solves for the equilibrium fluid temperature and electron fraction. We apply this method to snapshots from two-dimensional simulations of accretion disks left behind by binary neutron star mergers, varying the input physics and comparing to the results obtained with a leakage scheme for the cases of a central black hole and a central hypermassive neutron star. Neutrinos are guided away from the densest regions of the disk and escape preferentially around 45° from the equatorial plane. Neutrino heating is strengthened by MC transport a few scale heights above the disk midplane near the innermost stable circular orbit, potentially leading to a stronger neutrino-driven wind. Neutrino cooling in the dense midplane of the disk is stronger when using MC transport, leading to a globally higher cooling rate by a factor of a few and a larger leptonization rate by an order of magnitude. We calculate neutrino pair annihilation rates and estimate that an energy of 2.8 × 10{sup 46} erg is deposited within 45° of the symmetry axis over 300 ms when a central BH is present. Similarly, 1.9 × 10{sup 48} erg is deposited over 3 s when an HMNS sits at the center, but neither estimate is likely to be sufficient to drive a gamma-ray burst jet.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

A XANES and EPMA study of Fe3+ in chlorite: Importance of oxychlorite and implications for cation site distribution and thermobarometry

Chlorite is a ubiquitous product of metamorphism, alteration of magmatic rocks and hydrothermal processes owing to its large stability field and wide compositional range. Its composition is governed by several substitutions and has been used as a geothermometer, on the basis of empirical, semi-empirical, and thermodynamic models. As in some other phyllosilicates of petrological interest, the oxidation state of iron in chlorite may differ from the usually assumed divalent state. Yet, the crystal chemistry of trivalent iron in chlorite remains poorly known, and the thermodynamic properties of ferric chlorite are missing from databases used for petrological modeling. As part of an attempt to fill this gap, we present conclusions from in situ, micrometer-scale measurements of the oxidation state of iron in various chlorite-bearing samples. X-ray absorption near-edge spectroscopy (XANES) was combined with electron probe microanalysis (EPMA) on the same crystals. Results show iron oxidation states varying from ferrous to ferric; iron is in octahedral coordination in all ferromagnesian chlorites but to ~ 25% tetrahedral in the lithian chlorite cookeite (1.0 wt% Fe 2 O 3(total) ). Absolute amounts of ferric iron cover an unprecedented range (0 to ~ 30 wt% Fe 2 O 3 ). For highly magnesian, ferric chlorite, Fe concentrations are low and can be accounted for by Al = Fe 3+ substitution. In Fe-rich samples, Fe 3+ may exceed 2 atoms per formula unit (pfu, 18 oxygen basis). When structural formulas are normalized to 28 charges corresponding to the standard O 10 (OH) 8 anionic basis, these measurements define the exchange vector of a di-trioctahedral-type substitution: 3 VI (Mg, Fe 2+ ) = VI $\square$ + 2 VI Fe 3+ , as described in earlier studies. However, structural formulas calculated on the basis of the oxygen contents actually measured by EPMA show that this trend is an artifact, due to the neglect of variations in the number of protons in the structure. Our measurements indicate increasing hydrogen deficiency with increasing Fe 3+ content, up to ~ to 2 H + pfu in the Fe 3+ -rich chlorite samples, corresponding to a net exchange vector of the type R 2+ + H + = Fe 3+ . These results do not support substitutions toward di-trioctahedral ferric end-members, and highlight the need for considering substitution toward an "oxychlorite" (i.e., H-deficient) ferric component, close to tri-trioctahedral, with an O 12 (OH) 6 anionic basis, even in green, pristine-looking chlorite. The effects of iron oxidation and H deficiency on chlorite geothermometers were explored. They are deterring if H deficiency is ignored but, given the sensitivity of most thermometers to octahedral vacancy, the assumption Fe total = Fe 2+ is still safer than using high measured Fe 3+ contents and the standard 28 charge basis, which artificially increases vacancies. In such ferric chlorites, EPMA measurement of oxygen allows a fair estimate of H content if Fe 3+ /Fe 2+ is known; it should be more systematically implemented. For the same reasons, literature data reporting Fe 3+ -rich chlorite with vacancy content along the possibly artificial di-trioctahedral-type substitution should be verified. With the help of constraints from thermodynamic models, charge balance, crystal symmetry, and proton loss, a new cation site distribution is proposed for di-tri- to tri-trioctahedral chlorites in the Fe 2+ -Fe 3+ -Mg-Al-Si-O-H system, allowing a more realistic thermodynamic handling of their solid solutions.

58 GEOSCIENCES↗

DOE OSTI2019

Data Analytics for Electrical Distribution Systems with Micro PMUs (GMLC 1.4.9 Technical Report)

This report describes how electrical utilities with residential, industrial, commercial and government customers make decisions using Data Analytics (DA) applications and discusses how OpenFMB was leveraged to share that data among disparate types of devices for DA. OpenFMB is a data standard developed by utilities for distributed communication and control of power system assets. Here we explored how OpenFMB can be used for sharing analytics and what those analytics could be for a set of representative grid devices. A review of DA applications performed by electrical utilities and customers is presented. An overview of how the testbed we used was built and configured is summarized. Several experiments which were carried out to explore the capabilities of OpenFMB alongside other data sharing mechanisms for gathering sensor measurements and sharing analytics is also discussed. And finally a comparison of OpenFMB and LBNL’s SPARK framework is summarized. Today’s electrical utilities focus on monitoring the real-time customer consumptions, offering new real-time tariff models based on customer usage habits (rewards and punishments), providing energy peak during peak load times (plan energy generation and distribution), improving the load management for transformers and power lines, and controlling the energy and grid resources efficiently. On the other hand, the DA applications for residential customers are based on modifying the customer electricity usage habit to save money, making accurate and timely decisions for electrical billing. However, the DA applications for industrial, commercial and government customers focus on creating efficient programs to improve the environmental condition (reducing CO₂ emission), developing better load forecasting for industrial and/or commercial activities. In this project, a testbed was built to explore how the OpenFMB framework can be used to distribute data among various types of representative grid edge devices. Micro PMUs (Phasor Measurement Unit) were used to collect and transmit data to an OpenPDC server alongside relays with a framework configured with a system using OpenFMB to collect and transmit data. The testbed was composed of a physical asset based on a 60 Hz AC radial power system constructed using SI-GRID components. The power system consisted of a 3-phase 14 volts line-to-ground (24 volts line-to-line) power source, three 3-phase line emulators, a programmable resistive load bank, two smart relays, two micro PMUs, and a (SCADA) and Phasor Data Concentrator (PDC) computer. This testbed configuration allowed us to perform breaker operation tests collecting data from the micro PMUs and from the relays’ sensors. During the micro PMU breaker operation tests we collected the measured and estimated fundamental voltages, currents, frequency as well as the total true, reactive and apparent power. The displacement power factor from two micro PMUs located in parallel power lines was also used. The software steps for the implementation of the Synchrophasor system are explained in detail in this report. The setting (1-2-3) and data collection (4-5) phases were performed with a free download open source software. In the setting phase, the micro PMUs were set (1), the data stream for each micro PMU was validated, and the connection and configuration files were created (2). The PDC computer was set using the connection and configuration files of the micro PMUs (3). In the data collection phase, the binary data from the micro PMUs was collected (4), and the binary data files of the micro PMUs was converted to analyze the test results (5). OpenFMB was configured to use protocol translators to read Modbus and transform the reading into the OpenFMB data model to be read for DA applications.Our experiments showed that by leveraging the OpenFMB data model and the micro PMUs together we could obtain better situational awareness of network faults, which in some cases could lead to failures in the system and would be less perceivable otherwise. For example, in one experiment the fundamental RMS current angle, and the fundamental displacement power factor (estimated by the ratio between the measured real and apparent total power) had a disturbance at the micro PMU connected to the beginning of the parallel power line because the breaker was opened without the operators or SCADA system’s knowledge, and no current was available. Based on the breaker state scenarios and micro PMU locations, the DA applications using the current angle and displacement power factor from the micro PMUs should be studied in detail to avoid a non-desired measurement that could affect the decision from the DA. Also, the sample rates for PMU, Relay and micro PMU devices were compared. The DA application improved the control and operation between the electrical utilities and their customers (residential, industrial, transportation, commercial and government). The higher the sample rate, the better the quality of the measurement for DA estimations.

42 ENGINEERING↗

NASA NTRS1994

Geometrical evidence for dark matter: X-ray constraints on the mass of the elliptical galaxy NGC 720

We describe (1) a new test for dark matter and alternate theories of gravitation based on the relative geometries of the X-ray and optical surface brightness distributions and an assumed form for the potential, of the optical light, (2) a technique to measure the shapes of the total gravitating matter and dark matter of an ellipsoidal system which is insensitive to the precise value of the temperature of the gas and to modest temperature gradients, and (3) a new method to determine the ratio of dark mass to stellar mass that is dependent on the functional forms for the visible star, gas and dark matter distributions, but independent of the distance to the galaxy or the gas temperature. We apply these techniques to X-ray data from the ROSAT Position Sensitive Proportional Counter (PSPC) of the optically flattened elliptical galaxy NGC 720; the optical isophotes have ellipticity epsilon approximately 0.40 extending out to approximately 120 sec. The X-ray isophotes are significantly elongated, epsilon = 0.20-0.30 for semimajor axis a approximately 100 sec. The major axes of the optical and X-ray isophotes are misaligned by approximately 30 deg +/- 15 deg. Spectral analysis of the X-ray data reveals no evidence of temperature gradients or anisotropies and demonstrates that a single-temperature plasma (T approximately 0.6 keV) having subsolar heavy element abundances and a two-temperature model having solar abundances describe the spectrum equally well. Considering only the relative geometries of the X-ray and optical surface brightness distributions and an assumed functional form for the potential of the optical light, we conclude that matter distributed like the optical light cannot produce the observed ellipticities of the X-ray isophotes, independent of the gas pressure, the gas temperature, and the value of the stellar mass; this comparison assumes a state of quasi-hydrostatic equilibrium so that the three-dimensional surfaces of the gas emissivity trace the three-dimensional isopotential surfaces -- we discuss the viability of this assumption for NGC 720. Milgrom's Modification of Newtonian Dynamics (MOND) cannot dispel this manifestation of dark matter. Hence, geometrical considerations require, without mention of pressure or temperature, the presence of an extended, massive dark matter halo in NGC 720. Employing essentially the technique of Buote & Canizares (1992; Buote 1992) we use the shape of the X-ray surface brightness to constrain the shape of the total gravitating matter. The total matter is modeled as either an oblate or prolate spheriod of constant shape and orientation having either a Ferrers (rho approximately r(exp -n)) or Hernquist density. Assuming the X-ray gas is in hydrostatic equilibrium, we construct a model X-ray gas distribution for various temperature profiles. We determine the ellipticity of the total gravitating matter to be epsilon approximately 0.50-0.70. Using the single-temperature model we estimate a total mass approximately (0.41-1.4) x 10(exp 12) h(sub 80) solar mass interior to the ellipsoid of semimajor axis 43.6 h(sub 80) kpc. Ferrers densities as steep as r(exp -3) do not fit the data, but the r(exp -2) and Hernquist models yield excellent fits. We estimate the mass distributions of the stars and the gas and fit the dark matter directly. For a given gas equation of state and functional forms for the visible stars, gas, and dark matter, these models yield a distance-independent and temperature-independent measurement of the ratio of dark mass to stellar mass M(sub DM)/M(sub stars). We estimate a minimum M(sub DM)/M(sub stars) greater than or equal to 4 which corresponds to a total mass slightly greater than that derived from the single-temperature models for distance D = 20h(sub 80) Mpc.

Buote, David A.↗

DOE OSTI2018

NREL Pyrheliometer Comparisons: September 24 - October 5, 2018 (NPC-2018)

Accurate measurements of direct normal (beam) solar irradiance from pyrheliometers are important for developing and deploying solar energy conversion systems, improving our understanding of the Earth's energy budget for climate change studies, and for other science and technology applications involving solar flux. Providing these measurements places many demands on the quality system used by the operator of commercially available radiometers. Maintaining accurate radiometer calibrations that are traceable to an international standard is the first step in producing research-quality solar irradiance measurements. In 1977, the World Meteorological Organization (WMO) established the World Radiometric Reference (WRR) as the international standard for the measurement of direct normal solar irradiance (Frohlich 1991). The WRR is an internationally recognized, detector-based measurement standard determined by the collective performance of six electrically self-calibrated absolute cavity radiometers comprising the World Standard Group (WSG). Various countries, including the United States, have contributed these specialized radiometers to the Physikalisch-Meteorologisches Observatorium Davos - World Radiation Center (PMOD/WRC) to establish the WSG. As with all measurement systems, Absolute Cavity Radiometers (ASR) are subject to performance changes over time. Therefore, PMOD/WRC in Davos, Switzerland, hosts an International Pyrheliometer Comparison (IPC) for transferring the WRR to participating radiometers quinquennialy by invitation. NREL has represented the U.S. Department of Energy (DOE) in each IPC since 1980. As a result, NREL has developed and maintained a select group of absolute cavity radiometers with direct calibration traceability to the WRR, and uses these reference instruments to calibrate pyrheliometers and pyranometers using the International Organization for Standardization (ISO) 17025 accredited Broadband Outdoor Radiometer Calibration (BORCAL) process (Reda et al. 2008). To fill the gap between each IPC, NREL pyrheliometer comparisons (NPCs) are held annually at the Solar Radiation Research Laboratory (SRRL) in Golden, Colorado. Open to all ACR owners and operators, each NPC provides an opportunity to determine the unique WRR transfer factor (WRR-TF) for each participating pyrheliometer. By adjusting all subsequent pyrheliometer measurements by the appropriate WRR-TF, the solar irradiance data are traceable to the WRR. NPC-2018 was September 24 through October 5, 2018. Participants operated 43 ACRs to simultaneously measure clear-sky direct normal solar irradiance during this period. The Transfer Standard Group (TSG) of reference radiometers for NPC-2018 consisted of four NREL radiometers with direct traceability to the WRR, having participated in the Twelfth International Pyrheliometer Comparisons (IPC-XII) in the fall of 2015. As a result of NPC-2018, each participating absolute cavity radiometer was assigned a new WRR-TF, computed as the reference irradiance determined by the TSG divided by the observed irradiance from the participating radiometer. The performance of the TSG during NPC-2018 was consistent with previous comparisons, including IPC-XII. The measurement performance of the TSG allowed the transfer of the WRR to each participating radiometer with an estimated uncertainty of plus/minus 0.36% with respect to the International System of Units.

14 SOLAR ENERGY↗

NASA NTRS2017

NASA's Space Launch System: Deep-Space Delivery for Smallsats

Designed for human exploration missions into deep space, NASA's Space Launch System (SLS) represents a new spaceflight infrastructure asset, enabling a wide variety of unique utilization opportunities. While primarily focused on launching the large systems needed for crewed spaceflight beyond Earth orbit, SLS also offers a game-changing capability for the deployment of small satellites to deep-space destinations, beginning with its first flight. Currently, SLS is making rapid progress toward readiness for its first launch in two years, using the initial configuration of the vehicle, which is capable of delivering 70 metric tons (t) to Low Earth Orbit (LEO). On its first flight test of the Orion spacecraft around the moon, accompanying Orion on SLS will be small-satellite secondary payloads, which will deploy in cislunar space. The deployment berths are sized for "6U" CubeSats, and on EM-1 the spacecraft will be deployed into cislunar space following Orion separate from the SLS Interim Cryogenic Propulsion Stage. Payloads in 6U class will be limited to 14 kg maximum mass. Secondary payloads on EM-1 will be launched in the Orion Stage Adapter (OSA). Payload dispensers will be mounted on specially designed brackets, each attached to the interior wall of the OSA. For the EM-1 mission, a total of fourteen brackets will be installed, allowing for thirteen payload locations. The final location will be used for mounting an avionics unit, which will include a battery and sequencer for executing the mission deployment sequence. Following the launch of EM-1, deployments of the secondary payloads will commence after sufficient separation of the Orion spacecraft to the upper stage vehicle to minimize any possible contact of the deployed CubeSats to Orion. Currently this is estimated to require approximately 4 hours. The allowed deployment window for the CubeSats will be from the time the upper stage disposal maneuvers are complete to up to 10 days after launch. The upper stage will fly past the moon at a perigee of approximately 100km, and this closest approach will occur about 5 days after launch. The limiting factor for the latest deployment time is the available power in the sequencer system. Several NASA Mission Directorates were involved in the development of programs for the competition, selection, and development of EM-1 payloads that support directorate priorities. CubeSat payloads on EM-1 will include both NASA research experiments and spacecraft developed by industry, international and potentially academia partners. The Human Exploration and Operations Mission Directorate (HEOMD) Advanced Exploration Systems (AES) Division was allocated five payload opportunities on the EM-1 mission. Near Earth Asteroid (NEA) Scout is designed to rendezvous with and characterize a candidate NEA. A solar sail, an innovation the spacecraft will demonstrated for the CubeSat class, will provide propulsion. Lunar Flashlight will use a green propellant system and will search for potential ice deposits in the moon's permanently shadowed craters. BioSentinel is a yeast radiation biosensor, planned to measure the effects of space radiation on deoxyribonucleic acid (DNA). Lunar Icecube, a collaboration with Morehead State University, will prospect for water in ice, liquid, and vapor forms as well as other lunar volatiles from a low-perigee, highly inclined lunar orbit using a compact Infrared spectrometer. Skyfire, a partnership with Lockheed Martin, is a technology demonstration mission that will perform a lunar flyby, collecting spectroscopy, and thermography data to address questions related to surface characterization, remote sensing, and site selection. NASA's Space Technology Mission Directorate (STMD) was allocated three payload opportunities on the EM-1 mission. These slots will be filled via the Centennial Challenges Program, NASA's flagship program for technology prize competitions, which directly engages the public, academia, and industry in open prize competitions to stimulate innovation. The NASA Science Mission Directorate (SMD) was allocated two payload opportunities on the EM-1 mission. The CubeSat Mission to Study Solar Particles (CuSP) payload will study the sources and acceleration mechanisms of solar and interplanetary particles in near-Earth orbit, support space weather research by determining proton radiation levels during Solar Energetic Particle (SEP) events and identifying suprathermal properties that could help predict geomagnetic storms. The LunaH-Map payload will help scientists understand the quantity of H-bearing materials in lunar cold traps (~10 km), determine the concentration of H-bearing materials with 1m depth, and constrain the vertical distribution of H-bearing materials. The final three payload opportunities for the EM-1 mission were allocated for NASA's international space agency counterparts. The flight opportunities are intended to benefit the international space agency and NASA as well as further the collective space exploration goals. ArgoMoon is sponsored by ESA/ASI and will fly along with the ICPS on its disposal trajectory to perform proximity operations with the ICPS post-disposal, take external imagery of engineering and historical significance, and perform an optical communications demonstration. EQUULEUS, sponsored by JAXA, will fly to a libration orbit around the Earth-Moon L2 point and demonstrate trajectory control techniques within the Sun-Earth- Moon region for the first time by a nano spacecraft. The mission will also contribute to the future human exploration scenario by understanding the radiation environment in geospace and deep space, characterizing the flux of impacting meteors on the far side of the moon, and demonstrating the future deep space exploration scenario using the "deep space port" at Lagrange points. OMOTENASHI, also sponsored by JAXA, will land the smallest lunar lander to date on the lunar surface to demonstrate the feasibility of the hardware for distributed cooperative exploration system. Small landers will enable multi-point exploration, which is complimentary with large-scale human exploration. Once on the lunar surface, the OMOTENASHI spacecraft will observe the radiation and soil environments of the lunar surface by active radiation measurements and soil shear measurements. Following EM-1, Space Launch System will evolve to the more-powerful Block 1B configuration, which uses a new Exploration Upper Stage to increase the vehicle's LEO payload capability from 70 t to 105 t. With that transition, the Orion Stage Adapter, which will carry the secondary payloads on EM-1, will be phased out, and a new Universal Stage Adapter will be introduced, creating opportunities for flying larger secondary payloads. This paper will provide a brief status of SLS progress toward first launch; an overview of smallsat accommodations, integration, and operations on EM-1; information about the specific payloads flying on that launch; and a discussion of future accommodations and opportunities for secondary payloads on SLS for Exploration Mission-2 and beyond.

Robinson, Kimberly F.↗

NASA NTRS2017

NASA's Space Launch System: Deep-Space Opportunities for SmallSats

Designed for human exploration missions into deep space, NASA's Space Launch System (SLS) represents a new spaceflight infrastructure asset, enabling a wide variety of unique utilization opportunities. While primarily focused on launching the large systems needed for crewed spaceflight beyond Earth orbit, SLS also offers a game-changing capability for the deployment of small satellites to deep-space destinations, beginning with its first flight. Currently, SLS is making rapid progress toward readiness for its first launch in two years, using the initial configuration of the vehicle, which is capable of delivering 70 metric tons (t) to Low Earth Orbit (LEO). On its first flight test of the Orion spacecraft around the moon, accompanying Orion on SLS will be small-satellite secondary payloads, which will deploy in cislunar space. The deployment berths are sized for "6U" CubeSats, and on EM-1 the spacecraft will be deployed into cislunar space following Orion separate from the SLS Interim Cryogenic Propulsion Stage. Payloads in 6U class will be limited to 14 kg maximum mass. Secondary payloads on EM-1 will be launched in the Orion Stage Adapter (OSA). Payload dispensers will be mounted on specially designed brackets, each attached to the interior wall of the OSA. For the EM-1 mission, a total of fourteen brackets will be installed, allowing for thirteen payload locations. The final location will be used for mounting an avionics unit, which will include a battery and sequencer for executing the mission deployment sequence. Following the launch of EM-1, deployments of the secondary payloads will commence after sufficient separation of the Orion spacecraft to the upper stage vehicle to minimize any possible contact of the deployed cubesats to Orion. Currently this is estimated to require approximately 4 hours. The allowed deployment window for the cubesats will be from the time the upper stage disposal maneuvers are complete to up to 10 days after launch. The upper stage will fly past the moon at a perigee of approximately 100km, and this closest approach will occur about 5 days after launch. The limiting factor for the latest deployment time is the available power in the sequencer system. Several NASA Mission Directorates were involved in the development of programs for the competition, selection, and development of EM-1 payloads that support directorate priorities. CubeSat payloads on EM-1 will include both NASA research experiments and spacecraft developed by industry, international and potentially academia partners. The Human Exploration and Operations Mission Directorate (HEOMD) Advanced Exploration Systems (AES) Division was allocated five payload opportunities on the EM-1 mission. Near Earth Asteroid (NEA) Scout is designed to rendezvous with and characterize a candidate NEA. A solar sail, an innovation the spacecraft will demonstrated for the CubeSat class, will provide propulsion. Lunar Flashlight will use a green propellant system and will search for potential ice deposits in the moon's permanently shadowed craters. BioSentinel is a yeast radiation biosensor, planned to measure the effects of space radiation on deoxyribonucleic acid (DNA). Lunar Icecube, a collaboration with Morehead State University, will prospect for water in ice, liquid, and vapor forms as well as other lunar volatiles from a low-perigee, highly inclined lunar orbit using a compact Infrared spectrometer. Skyfire, a partnership with Lockheed Martin, is a technology demonstration mission that will perform a lunar flyby, collecting spectroscopy, and thermography data to address questions related to surface characterization, remote sensing, and site selection. NASA's Space Technology Mission Directorate (STMD) was allocated three payload opportunities on the EM-1 mission. These slots will be filled via the 2 Centennial Challenges Program, NASA's flagship program for technology prize competitions, which directly engages the public, academia, and industry in open prize competitions to stimulate innovation. The NASA Science Mission Directorate (SMD) was allocated two payload opportunities on the EM-1 mission. The CubeSat Mission to Study Solar Particles (CuSP) payload will study the sources and acceleration mechanisms of solar and interplanetary particles in near-Earth orbit, support space weather research by determining proton radiation levels during Solar Energetic Particle (SEP) events and identifying suprathermal properties that could help predict geomagnetic storms. The LunaH-Map payload will help scientists understand the quantity of H-bearing materials in lunar cold traps (~10 km), determine the concentration of H-bearing materials with 1m depth, and constrain the vertical distribution of H-bearing materials. The final three payload opportunities for the EM-1 mission were allocated for NASA's international space agency counterparts. The flight opportunities are intended to benefit the international space agency and NASA as well as further the collective space exploration goals. ArgoMoon is sponsored by ESA/ASI and will fly along with the ICPS on its disposal trajectory to perform proximity operations with the ICPS post-disposal, take external imagery of engineering and historical significance, and perform an optical communications demonstration. EQUULEUS, sponsored by JAXA, will fly to a libration orbit around the Earth-Moon L2 point and demonstrate trajectory control techniques within the Sun-Earth-Moon region for the first time by a nano spacecraft. The mission will also contribute to the future human exploration scenario by understanding the radiation environment in geospace and deep space, characterizing the flux of impacting meteors on the far side of the moon, and demonstrating the future deep space exploration scenario using the "deep space port" at Lagrange points. OMOTENASHI, also sponsored by JAXA, will land the smallest lunar lander to date on the lunar surface to demonstrate the feasibility of the hardware for distributed cooperative exploration system. Small landers will enable multi-point exploration, which is complimentary with large-scale human exploration. Once on the lunar surface, the OMOTENASHI spacecraft will observe the radiation and soil environments of the lunar surface by active radiation measurements and soil shear measurements. Following EM-1, Space Launch System will evolve to the more-powerful Block 1B configuration, which uses a new Exploration Upper Stage to increase the vehicle's LEO payload capability from 70 t to 105 t. With that transition, the Orion Stage Adapter, which will carry the secondary payloads on EM-1, will be phased out, and a new Universal Stage Adapter will be introduced, creating opportunities for flying larger secondary payloads. This paper will provide a brief status of SLS progress toward first launch; an overview of smallsat accommodations, integration, and operations on EM-1; information about the specific payloads flying on that launch; and a discussion of future accommodations and opportunities for secondary payloads on SLS for Exploration Mission-2 and beyond.

Robinson, Kimberly F.↗

DOE OSTI2013

An Analytical Impact Assessment Framework for Wildlife to Inform the Siting and Permitting of Wind Energy Facilities

In the United States overall electrical generation capacity is expected to increase by 10-25 gigawatts (GW) per year to meet increases in demand. Wind energy is a key component of state and federal renewable energy standards, and central to the Department of Energy’s 20% by 2030 wind production goals. Increased wind energy development may present increased resource conflict with avian wildlife, and environmental permitting has been identified as a potential obstacle to expansion in the sector. ICF developed an analytical framework to help applicants and agencies examine potential impacts in support of facility siting and permitting. A key objective of our work was to develop a framework that is scalable from the local to the national level, and one that is generalizable across the different scales at which biological communities operate – from local influences to meta-populations. The intent was to allow natural resource managers to estimate the cumulative impacts of turbine strikes and habitat changes on long-term population performance in the context of a species demography, genetic potential, and life history. We developed three types of models based on our literature review and participation in the scientific review processes. First, the conceptual model was developed as a general description of the analytical framework. Second, we developed the analytical framework based on the relationships between concepts, and the functions presented in the scientific literature. Third, we constructed an application of the model by parameterizing the framework using data from and relevant to the Altamont Pass Wind Resource Area (APWRA), and an existing golden eagle population model. We developed managed source code, database create statements, and written documentation to allow for the reproduction of each phase of the analysis. ICF identified a potential template adaptive management system in the form of the US Fish & Wildlife Service (USFWS) Adaptive Harvest Management (AHM) program, and developed recommendations for the structure and function of a similar wind-facility related program. We provided a straw-man implementation of the analytical framework based on assumptions for APWRA-wide golden eagle fatalities, and presented a statistical examination of the model performance. APWRA-wide fatality rates appear substantial at all scales examined from the local APWRA population to the Bird Conservation Region. Documented fatality rates significantly influenced population performance in terms of non-territorial non-breeding birds. Breeder, Juvenile, Subadult, and Adult abundance were mostly unaffected by Baseline APWRA-wide fatality rates. However, increased variability in fatality rates would likely have impacts on long-term population performance, and would result in a substantially larger loss of resources. We developed four recommendations for future study. First, we recommend establishment of concept experts through the existing system of non-profits, regulatory agencies, academia, and industry in the wind energy sector. Second, we recommend the development of a central or distributed shared data repository, and establish guidelines for data sharing and transparency. Third, we recommend development a forum and process for model selection at the local and national level. Last, we recommend experimental implementation of the prescribed system at broader scales, and refinement the expectations for modeling and adaptive management.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Application of NEAMS Codes to Capture MSR Phenomena

This report documents work completed in FY21 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program’s Molten Salt Reactor (MSR) Application Drivers activity at Argonne. The common focus was on identifying the modeling and simulation functional requirements for designing, licensing, and operating MSRs and applying those capabilities already developed in NEAMS tools to example problems of interest. The four main parts of this report each focuses on a specific area of simulation physics as it relates to MSR phenomena: fuel evolution, chemistry, computational fluid dynamics (CFD), and systems analysis. In terms of fuel evolution, which includes depletion, decay, on-line separations, and transmutation, the current state of computational capabilities for modeling this behavior in liquid-fueled molten salt reactor is discussed. Some of the functional requirements to accomplish the various applications of MSR fuel depletion modeling are highlighted, followed by a summary of recent approaches and code development activities. The chemistry functional requirements were discussed in the context of several applications of high importance for MSRs, such as corrosion, salt chemistry, and species transport. Each of these types of chemistry modeling have considerable impact on various aspects of reactor applications, including informing on reactor designs, improving operational efficiencies, analyzing safety and reactivity concerns, and estimating the mechanistic source term of the reactor. A brief overview is also provided on code development activities ongoing under NEAMS relevant to chemistry modeling of MSRs. In terms of CFD applications, the state-of-the-art spectral element code, Nek5000 was used to model the fluid dynamics within a full core of the Molten Salt Fast Reactor (MSFR) concept designed as part of the Euratom EVOL project. This concept was selected as the challenge problem because of its similar features to several U.S. industry concepts. The goal was to model some of the fast MSR design challenges, including potential large internal re-circulations, the need of accurate tracking of delayed neutron precursors (DNP), and the lack of relevant thermal-hydraulics models/correlations, etc. Therefore, a series of CFD models were created for the MSFR core cavity using a k – τ model two-equation model for the turbulence. These first full core results demonstrated that any potential recirculation zones could be properly identified with the current NEAMS CFD capabilities. The development of these models will also set the stage for future testing of Nek5000’s functionalities to model other MSR thermal-fluid phenomena. Lastly, the validation of SAM against experimental data from the Molten Salt Reactor Experiment, which started in FY19, continues with the inclusion of modeling the reactivity insertion tests. This involved using SAM and its point kinetics model for flowing fuel salt to recreate the time dependent power changes and response after positive reactivity insertions at the 1, 5, and 8 MWt power levels. Through these exercises, several code modifications were suggested to the SAM development team and accommodated to enable closer agreement with solid technical and physical justifications. These include adding a moderator reactivity feedback coefficient as an available input and modifying the solution approach for the point kinetics model. Additional SAM development suggestions for flowing fuel MSRs include adding the capability to allow the moderator power change proportionally with the reactor power and enabling specification of the power and DNP distributions separately.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2013

Virtual Design Studio (VDS) - Development of an Integrated Computer Simulation Environment for Performance Based Design of Very-Low Energy and High IEQ Buildings

Executive Summary The objective of this study was to develop a “Virtual Design Studio (VDS)”: a software platform for integrated, coordinated and optimized design of green building systems with low energy consumption, high indoor environmental quality (IEQ), and high level of sustainability. This VDS is intended to assist collaborating architects, engineers and project management team members throughout from the early phases to the detailed building design stages. It can be used to plan design tasks and workflow, and evaluate the potential impacts of various green building strategies on the building performance by using the state of the art simulation tools as well as industrial/professional standards and guidelines for green building system design. Engaged in the development of VDS was a multi-disciplinary research team that included architects, engineers, and software developers. Based on the review and analysis of how existing professional practices in building systems design operate, particularly those used in the U.S., Germany and UK, a generic process for performance-based building design, construction and operation was proposed. It distinguishes the whole process into five distinct stages: Assess, Define, Design, Apply, and Monitoring (ADDAM). The current VDS is focused on the first three stages. The VDS considers building design as a multi-dimensional process, involving multiple design teams, design factors, and design stages. The intersection among these three dimensions defines a specific design task in terms of “who”, “what” and “when”. It also considers building design as a multi-objective process that aims to enhance the five aspects of performance for green building systems: site sustainability, materials and resource efficiency, water utilization efficiency, energy efficiency and impacts to the atmospheric environment, and IEQ. The current VDS development has been limited to energy efficiency and IEQ performance, with particular focus on evaluating thermal performance, air quality and lighting environmental quality because of their strong interaction with the energy performance of buildings. The VDS software framework contains four major functions: 1) Design coordination: It enables users to define tasks using the Input-Process-Output flow approach, which specifies the anticipated activities (i.e., the process), required input and output information, and anticipated interactions with other tasks. It also allows task scheduling to define the work flow, and sharing of the design data and information via the internet. 2) Modeling and simulation: It enables users to perform building simulations to predict the energy consumption and IEQ conditions at any of the design stages by using EnergyPlus and a combined heat, air, moisture and pollutant simulation (CHAMPS) model. A method for co-simulation was developed to allow the use of both models at the same time step for the combined energy and indoor air quality analysis. 3) Results visualization: It enables users to display a 3-D geometric design of the building by reading BIM (building information model) file generated by design software such as SketchUp, and the predicted results of heat, air, moisture, pollutant and light distributions in the building. 4) Performance evaluation: It enables the users to compare the performance of a proposed building design against a reference building that is defined for the same type of buildings under the same climate condition, and predicts the percent of improvements over the minimum requirements specified in ASHRAE Standard 55-2010, 62.1-2010 and 90.1-2010. An approach was developed to estimate the potential impact of a design factor on the whole building performance, and hence can assist the user to identify areas that have most pay back for investment. The VDS software was developed by using C++ with the conventional Model, View and Control (MVC) software architecture. The software has been verified by using a simple 3-zone case building. The application of the VDS concepts and framework for building design and performance analysis has been illustrated by using a medium-sized, five story office building that received LEED Platinum Certification from USGBC.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2017

Geomechanics-Based Stochastic Analysis of Injection- Induced Seismicity

The production of geothermal energy from dry and low permeability reservoirs is achieved by water circulation in natural and/or man-made fractures, and is referred to as enhanced or engineered geothermal systems (EGS). Often, the permeable zones have to be created by stimulation, a process which involves fracture initiation and/or activation of discontinuities such as faults and joints due to pore pressure and the in-situ stress perturbations. The stimulation of a rock mass is often accompanied by multiple microseismic events. Micro-seismic events associated with rock failure in shear, and shear slip on new or pre-existing fracture planes and possibly their propagations. The microseismic signals contain information about the sources of energy that can be used for understanding the hydraulic fracturing process and the created reservoir properties. Detection and interpretation of microseismic events is useful for estimating the stimulated zone, created reservoir permeability and fracture growth, and geometry of the geological structures and the in-situ stress state. The process commonly is referred to as seismicity-based reservoir characterization (SBRC). Although, progress has been made by scientific & geothermal communities for quantitative and qualitative analysis of reservoir stimulation using SBRC several key questions remain unresolved in the analysis of micro-seismicity namely, variation of seismic activity with injection rate, delayed micro-seismicity, and the relation of stimulated zone to the injected volume and its rate, and the resulting reservoir permeability. In addition, the current approach to SBRC does not consider the full range of relevant poroelastic and thermoelastic phenomena and neglects the uncertainty in rock properties and in-situ stress in the data inversion process. The objective of this research and technology developments was to develop a 3D SBRC model that addresses these shortcomings by taking into account hydro-thermo-poro-mechanical mechanisms associated with injection and utilizing a state-of-the-art stochastic inversion procedure. The approach proposed herein is innovative and significantly improves the existing SBCR technology (e.g., Shapiro et al. 2003) for geothermal reservoirs in several ways. First, the current scope of the SBRC is limited with respect to the physical processes considered and the rock properties used. Usually, the geomechanics analyses within SBRC is limited to the pore pressure diffusion in the rock mass, which is modeled using a time-dependent parabolic equation and solved using a finite element algorithm with either a line or a point source. However, water injection induces both poroelastic and thermoelastic stresses in the rock mass which affect the stress state. In fact, it has been suggested that thermoelastic stresses can play a dominant role in reservoir seismicity (Ghassemi et al., 2007). We include these important effects by using a fully-coupled poro-thermoelastic constitutive equations for the rock mass which will be solved using a 3D finite element model with more realistic injection geometries such as multiple injection/extraction sources (and in fractures), uncertainty in the material parameters and the in-situ stress distribution to better reflect the pore pressure and stress distributions. In addition, we developed a 3D stochastic fracture network model to study MEQ generation in fracture rocks. The model was verified using laboratory experiments, and calibrated and applied to Newberry EGS stimulation. In previous SBRC approaches, the triggering of micro-seismicity is modeled base on the assumption that the prior stochastic criticality model of the rock mass is a valid and adequate description. However, this assumption often does not hold in the field. Thus, we improved upon the current SBRC approach by using the micro-seismic responses to estimate the hydraulic diffusivity as well as the criticality distribution itself within the field. In this way, instead of relying on our a priori knowledge of criticality distribution, we combine an initial probabilistic description of criticality with the information contained in microseismic measurements to arrive at criticality solutions that are conditioned on both field data and our prior knowledge. Previous SBRC have relied upon a deterministic inversion approach to estimate the permeability, and the extent of the stimulated zone, whereas a stochastic inversion algorithm that recognizes and quantifies the uncertainties in the prior model, the time evolution of pore pressure distributions (modeling errors), and the observed seismic events is developed and used herein to realistically assess the quality of the solution. Finally, we developed a technique for processing discrete MEQ data to estimate fracture network properties such as dip and dip directions. The approach was successfully applied to the Fenton Hill HRD experiment and the Newberry EGS with results in good agreement with field observations.

15 GEOTHERMAL ENERGY↗

DOE OSTI2026

Spatially Resolved Domicile Charging Demands for Light-, Medium-, and Heavy-Duty Electric Vehicles in Virginia

The use of plug-in electric vehicles (PEVs) and resulting grid impacts are likely to grow rapidly, and evaluation of optimal smart charge management and grid integration strategies is warranted now. Evaluating distribution grid impacts requires fine-grained models of PEV operations to estimate charging loads across diverse vehicles at high spatial resolution. We propose such a model and consider a high-electrification scenario in Richmond and Newport News, Virginia. Our framework considers four categories of vehicle that are amenable to early aggressive electrification: light-duty passenger vehicles (LDV), trucks and vans with a focus on delivery or other local operations, school buses, and transit buses. These vehicles have a relatively consistent domicile, reducing the need for public charging infrastructure rollout to electrify. We apply a recent LDV model and propose new models for each vocation of medium- and heavy-duty vehicle, leveraging telematics data. We demonstrate our framework in Virginia and find energy demands in the region may total 15 GWh day, with most consumed by LDV. However, considering power demand at high spatial resolution reveals a different trend: LDVs have relatively small peak loads at specific sites (peak site demand around 800 kW) compared to average and high demand medium- and heavy-duty vehicle charging sites (peak site demand around 6,000 kW at a transit bus depot, 1,500 kW at a local freight hub, and 1,000 kW at a school). Our framework yields insights on the relative impacts of each vocation and enables future work to tailor grid integration strategies to each vehicle category.

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS1999

A Study of Direct and Cloud-Mediated Radiative Forcing of Climate Due to Aerosols

The Intergovernmental Panel on Climate Change (IPCC) has reported that in the southeastern US and eastern China, the general greenhouse warming due to anthropogenic gaseous emissions is dominated by the cooling effect of anthropogenic aerosols. To verify this model prediction in eastern China and southeastern US, we analyzed regional patterns of climate changes at 72 stations in eastern China during 1951- 94 (44 years), and at 52 stations in the southeastern US during 1949-94 (46 years) to detect the fingerprint of aerosol radiative forcing. It was found that the mean rates of change of annual mean daily, maximum, minimum temperatures and diurnal temperature range (DTR) in eastern China were 0.8, -0.2, 1.8, and -2.0 C/100 years respectively, while the mean rates of change of annual mean daily, maximum, minimum temperatures and DTR in the southeastern US were -0.2, -0.6, 0.2, and -0.8 C/100 years, respectively. This indicates that the high rate of increase in annual mean minimum temperature in eastern China results in a slightly warming trend of daily temperature, while the high rate of decrease in annual mean maximum temperature and low rate of increase in annual mean minimum temperature lead to the cooling trend of daily temperature in the southeastern US. We found that the warming from the longwave forcing due to both greenhouse gases and aerosols was completely counteracted by the shortwave aerosol forcing in the southeastern US in the past 46 years. A slightly overall warming trend in eastern China is evident; winters have become milder. This finding is explained by hypothesizing that increasing energy usage during the past 44 years has resulted in more coal and biomass burning, thus increasing the emission of absorbing soot and organic aerosols in eastern China. Such emissions, in addition to well-known Asia dust and greenhouse gases, may be responsible for the winter warming trend in eastern China that we have reported here. The sensitivity of aerosol radiative properties to aerosol composition, size distribution, relative humidity (RH) is examined for the following aerosol systems: inorganic and organic ions (Cl-, Br-, NO3 -, SO4 2-, Na+, NH4 +, K+, Ca2+, Mg2+, HCOO-, CH3COO-, CH3CH2COO-, CH3COCOO-, OOCCOO2-, MSA-1); water-insoluble inorganic and organic compounds (elemental carbon, n-alkanes, SiO2, Al2O3, Fe2O3 and other organic compounds). The partial molar refraction method was used to calculate the real part of the refractive index. It was found that the asymmetry factor increased by approximately 48% with the real part varying from 1.40 to 1.65, and the single scattering albedo decreased by 24% with the imaginary part varying from -0.005 to -0.1. The asymmetry factor increased by 5.4 times with the geometric standard deviation varying from 1.2 to 3.0. The radiation transmission is very sensitive to the change in size distribution; other factors are not as significant. To determine the aerosol direct radiative forcing (ADRF), the aerosol optical depth (AOD) values at the three operational wavelengths (415, 500 and 673 nm) were determined at a regionally representative site, namely, Mt. Gibbs (35.78 deg N, 82.29 deg W, elevation 2006 m) in Mt. Mitchell State Park, NC, and a site located in an adjacent valley (Black Mountain, 35.66 deg N, 82.38 deg W, elevation 951 m) in the southeastern US. The two sites are separated horizontally by 10 km and vertically by 1 km. It was found that the representative total AOD values at 500 nm at the valley site for highly polluted (HP), marine (M) and continental (C) air masses were 0.68 +/- 0.33, 0.29 +/- 0.19 and 0.10 +/- 0.04, respectively. A search-graph method was used to retrieve the columnar size distribution (number concentration N, effective radius reff and geometric standard deviation=?g) from the optical depth observations at three operational wavelengths. The ground albedo, single scattering albedo and imaginary part of the refractive index were calculated using a mathematically unique procedure involving a Mie code and a radiative transfer code in conjunction with the retrieved aerosol size distribution, AOD, and diffuse-direct irradiance ratio. It was found that N, r(eff) and sigma(g) were in the ranges of 10 to 1.7 x 10(exp 4)/cubic cm, 0.09 to 0.68 micrometers and 1.12 to 2.95, respectively. The asymmetry factor and single scattering albedo were in the ranges of 0.63 to 0.75 and 0.74 to 0.97 respectively. The ground albedo for the forested terrain and imaginary part of refractive index were found to be in the ranges of 0.06 to 0.29 and 0.005 to 0.051 respectively. On the basis of these aerosol radiative properties obtained at the research sites and computations using the Column Radiation Model (CRM) of National Center of Atmospheric Research (NCAR) Community Climate Model (CCM3), it was found that the average cloud-free 24-hour ADRF values were -13 +/- 8, -8 +/- 3, -33 +/- 16 W/square m for marine, continental, and polluted air masses, respectively. On the assumption that the fractional coverage of clouds is 0.61, it was estimated that the annual mean ADRF was 7 +/- 2 W/square m in the southeastern US. The review with respect to the current knowledge of organic acids shows that aerosol formate and acetate concentrations range from 0.02 to 5.3 nmol/cubic m and from 0.03 to 12.4 nmol/cubic m respectively, and that between 34% to 77% of formate and between 21% to 66% of acetate are present in the fine fraction of aerosols. It was found that although most (98-99%) of these volatile organic acids were present in the gas phase, their concentrations in the aerosol particles were sufficient to make them a good candidate for cloud condensation nuclei (CCN). It is hypothesized that organic acids are at least one of the primary sources of CCN in the atmosphere due to their ubiquitous presence in the troposphere, especially over the continental forested areas. The results of our measurements at Palmer Station, Antarctica show that the daily average CCN concentrations at 0.3% and 1% supersaturations ranged from 0.3 to 160/cubic cm and from 4 to 168/cubic cm, respectively, during the period from 17 January to 26 February, 1994. New evidence for substantial and definitive CCN enhancement near and within cloud has been observed at Mt. Mitchell, North Carolina. The results show that the average monthly CCN concentrations were 460 +/- 217, 386 +/- 286, 429 +/- 228 and 238 +/- 134/cubic cm for in-cloud, overcast, clear and rainy conditions, respectively. The typical CCN spectra show that there were a lot of small CCN produced and the ion concentrations (especially H+ and SO4 2-) were very high during the CCN enhancement period. The significantly positive correlation between black carbon (BC) and CCN at 1% supersaturation indicates that a percentage of the BC measured at the site may be in the form of an internal mixture and participated in the formation of CCN.

Yu, Shao-Cai↗

NASA NTRSDate not supplied

Advancement of Entry System Modeling to Support Exploration of Giant Planets

This paper describes NASA’s efforts to advance entry system modeling and simulation capabilities to support future exploration of Giant planets. The Giant planets are key destinations of interest to the planetary science community for their potential to provide insight into the formation and evolution of our Solar System, as well as extrasolar planetary systems. To date, the Galileo atmospheric probe is the only purpose-built entry probe to a Giant planet. Post-flight analysis of Galileo’s performance showed that there was significant recession of the thermal protection system (TPS), well beyond what was anticipated on the flank, and this was due in part to insufficiently accurate capability for estimating the flight environment and TPS response. While Galileo ultimately survived its flight, the example serves to highlight the great challenge of designing successful missions for environments that are poorly understood or where models have not yet been validated. An important means to reduce mission risks is the incorporation of physics-based modeling with well-quantified uncertainties. The emphasis on physics-based modeling – in contrast to empirically-driven models – is motivated by the fact that it is impossible to completely replicate entry environments through ground tests and, therefore, extrapolation to the flight environment is required. Basing analysis in fundamental physics removes the bias of ground test limitations, though one must then be careful to properly characterize model inputs, simplifying assumptions, and the limits wherein the model is valid. NASA’s Entry Systems Modeling (ESM) Project is tasked with investigating such considerations for planetary science missions across the Solar System, and in recent years has begun to do so for Giant planets. The most distinctive features of the Giant planets, from an entry system perspective, are the atmospheres composed primarily of hydrogen and helium. The entry velocities of proposed missions are generally very large and can therefore be expected to result in significant convective and radiative heating generated by the vehicle’s shock layer. Yet thermochemical behavior of the hydrogen-helium system is not well understood under such conditions. The ESM project is leading efforts to develop accurate thermochemical databases based on state-of-the-art measurements in the Electric Arc Shock Tube and detailed computational chemistry. The large heat fluxes anticipated by missions has driven interest in new TPS materials, in particular woven materials, which may be enabling but have never been flown before. Consequently, multiscale models are in development to describe properties and performance of the materials from micro- to system-scale. The goal is to not only provide accurate thermal response but also to inform thermostructural reliability predictions for extreme entries. Additionally, new computational models have been developed to evaluate performance of non-destructive evaluation techniques which are vital to establishing acceptance of systems to be free of manufacturing faults like material cracking, voids, and debonding. Finally, in the area of guidance and control, aerocapture has been shown conceptually to provide a number of mission benefits, including reducing transit time and increasing payload fraction. The ESM project is building a launch-to-landing trajectory simulation capability to enable detailed studies of aerocapture maneuvers in the context of Giant planets missions. The final presentation and paper will describe each of these topics in detail, including discussion of specific gaps and the technical approach to solving them. In addition, the final paper will briefly discuss ongoing coordination between ESM project work and an ESA-funded technology development activity comprised of validation testing in the Oxford T6, IRS PWK and IST ESTHER tunnels, as well as state-to-state modeling of the shock layer to better represent non-Boltzmann energy distributions leading to non-equilibrium radiation.

Entry systems↗

DOE OSTI2025

The value of integrating a geothermal district heating system into a microgrid

As electrical grids increasingly rely on variable renewable energy, maintaining reliability and cost efficiency becomes more complex. To address these challenges, this study analyzed the integration of geothermal district heating as a grid-responsive thermal resource within a microgrid in Tuttle, Oklahoma. Building energy modeling using EnergyPlus estimated annual district heating demand at 2.9 GWh, with a peak load of 2.8 MW th . Techno-economic analyses were conducted to meet the heating demand under three geothermal scenarios, varying by production depth, flow rate, and thermal output, each supplemented by natural gas peaking boilers. In parallel, equivalent electrical load profiles were developed using typical coefficients of performance (COPs) for air-source heat pumps and electric boilers to establish an electrified baseline scenario. A complete end-use electrical load profile was also developed for the microgrid using Cambium dataset. The modeling results demonstrated reliable and economic operation of the geothermal systems over 30 years, with COPs ranging from 2.6 to 8.9 and the lowest levelized heating cost at $\$$54.6/MWh. Geothermal integration reduced electricity consumption by up to 94.7 % compared to the non-geothermal base case, yielding annual energy savings of up to $\$$803 k. Avoided grid costs ranged from $\$$65 k–$\$$147 k per year, with individual events avoiding up to $\$$4,863 per hour. Grid-responsive operation further reduced wholesale energy costs by 53–56 %. These findings demonstrate geothermal heating, traditionally treated as a non-grid-responsive thermal resource, can be reconfigured to support dynamic grid services, offering a scalable pathway to enhance reliability and reduce costs in renewable-rich microgrids and district heating networks.

15 GEOTHERMAL ENERGY↗

DOE OSTI2025

Evaluating the potential of short-term instrument deployment to improve distributed wind resource assessment

Distributed wind projects, which are connected at the distribution level of an electricity system or in off-grid applications to serve specific or local energy needs, often rely solely on wind resource models to establish wind speed and energy generation expectations. Historically, anemometer loan programs have provided an affordable avenue for more accurate onsite wind resource assessment, and the lowering cost of lidar systems has shown similar advantages for more recent assessments. While a full 12 months of onsite wind measurement is the standard for correcting model-based long-term wind speed estimates for utility-scale wind farms, the time and capital investment involved in gathering onsite measurements must be reconciled with the energy needs and funding opportunities that drive expedient deployment of distributed wind projects. Much literature exists to quantify the performance of correcting long-term wind speed estimates with 1 or more years of observational data, but few studies explore the impacts of correcting with months-long observational periods. This study aims to answer the question of how short you can go in terms of the observational time period needed to make impactful improvements to model-based long-term wind speed estimates. Three algorithms, multivariable linear regression, adaptive regression splines, and regression trees, are evaluated for their skill at correcting long-term wind resource estimates from the European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5) using months-long periods of observational data from 66 locations across the US. On average, correction with even 1 month of observations provides significant improvement over the baseline ERA5 wind speed estimates and produces median bias magnitudes and relative errors within 0.22 m s −1 and 4 percentage points of the median bias magnitudes and relative errors achieved using the standard 12 months of data for correction. However, in cases when the shortest observational periods (1 to 2 months) used for correction are not well correlated with the overlapping ERA5 reference, the resultant long-term wind speed errors are worse than those produced using ERA5 without correction. Summer months, which are characterized by weaker relative wind speeds and standard deviations for most of the evaluation sites, tend to produce the worst results for long-term correction using months-long observations. The three tested algorithms perform similarly for long-term wind speed bias; however, regression trees perform notably worse than multivariable linear regression and adaptive regression splines in terms of correlation when using 6 months or less of observational data for correction. Translating the analysis to wind energy, median relative errors in the capacity factor are on average within 10 % using 1 month of training. If the observation period used for correction is not well correlated with the reference data, however, misrepresentation of the observed capacity factor can be substantial. The risk associated with poor correlation between the observed and reference datasets decreases with increasing training period length. In the worst-correlation scenarios, the median capacity factor relative errors from using 1, 3, and 6 months are within 47 %, 26 %, and 16 %, respectively.

17 WIND ENERGY↗

DOE OSTI2013

Grassland Bird Distribution and Raptor Flight Patterns in the Competitive Renewable Energy Zones of the Texas Panhandle [Thesis]

The consistent wind resource in the Great Plains of North America has encouraged the development of wind energy facilities across this region. In the Texas Panhandle, a high quality wind resource is only one factor that has led to the expansion of wind energy development. Other factors include federal tax incentives and the availability of subsidies. Moreover, the State Renewable Portfolio Standards (RPS), mandating production of 10,000 mega-watts of renewable energy in the state by 2025, has contributed to an amicable regulatory and permitting environment (State Energy Conservation Office 2010). Considering the current rate of development, the RPS will be met in coming years (American Wind Energy Association 2011) and the rate of development is likely to continue. To meet increased energy demands in the face of a chronically constrained transmission grid, Texas has developed a comprehensive plan that organizes and prioritizes new transmission systems in high quality wind resource areas called Competitive Renewable Energy Zones (CREZ). The CREZ plan provides developers a solution to transmission constraints and unlocks large areas of undeveloped wind resource areas. In the northern Texas panhandle, there are two CREZs that are classified Class 3 wind (Class 5 is the highest) and range from 862,725 to 1,772,328 ha in size (Public Utility Commission of Texas 2008). Grassland bird populations have declined more than any other bird group in North America (Peterjohn and Sauer 1999, Sauer et al. 2004). Loss of grassland habitat from agricultural development has been the greatest contributor to the decline of grassland bird populations, but development of non-renewable (i.e., oil, coal, and gas) and renewable energy (i.e., wind, solar, biomass, and geothermal) sources have contributed to the decline as well (Pimentel et al. 2002, Maybe and Paul 2007). The effects of wind energy development on declining grassland bird populations has become an area of extensive research, as we attempt to understand and minimize potential impacts of a growing energy sector on declining bird populations. Based on data from post-construction fatality surveys, two grassland bird groups have been the specific focus of research, passerines (songbird guild) and raptors (birds of prey). The effects of wind energy development on these two groups of birds, both of conservation concern, have been examined over the last decade. The primary focus of this research has been on mortality resulting from collision with wind turbines (Kuvlesky et al. 2007). Most studies just quantify post-construction fatality levels (e.g., Erickson et al. 2002) while very few studies provide a comparison of bird populations prior to development through a Before-After-Control-Impact (BACI) study design. Before-After-Control-Impact studies provide powerful evidence of avian/wind energy relationships (Anderson et al. 1999). Despite repeated urgency on conducting these types of studies (Anderson et al. 1999, Madders and Whitfield 2006, Kuvelsky et al. 2007), few have been conducted in North America. Although several European researchers (Larsson 1994, de Lucas et al. 2007) have used BACI designs to examine whether wind facilities modified raptor behavior, there is a scarcity of BACI data relating to North America grassland ecosystems that examine avian-wind energy relationships. There are less than a handful of studies in the entire United States, let alone the southern short grass prairie ecosystem, that incorporate preconstruction data to form the baseline for post-construction impact estimates (Johnson et al. 2000, Erickson et al. 2002). Although declines in grassland bird populations are well-documented (Peterjohn and Sauer 1999, Sauer et al. 2004), the causal mechanisms affecting the decline of grassland birds with increasing wind energy development in the southern short grass prairie are not well-understood (Kuvlesky et al. 2007, Maybe and Paul 2007). Several factors may potentially affect the bird population when wind turbines are constructed in areas with high bird densities (de Lucas et al. 2007). Habitat fragmentation, noise from turbines, physical movement of turbine blades, and increased vehicle traffic have been suggested as causes of decreased density of nesting grassland birds in Minnesota (Leddy et al. 1999), Oklahoma (O’Connell and Piorkowski 2006), and South Dakota (Shaffer and Johnson 2008). Similarly, constructing turbines in areas where bird flight patterns place them at similar heights of turbine blades increases the potential for bird collisions (Johnson et al. 2000, Hoover 2002). Raptor fatalities have been associated with topographic features such as ridges, saddles and rims where birds use updrafts from prevailing winds (Erickson et al. 2000, Johnson et al. 2000, Barrios and Rodriquez 2004, Hoover and Morrison 2005). Thus, wind energy development can result in indirect (e.g., habitat avoidance, decreased nest success) and direct (e.g., collision fatalities) impacts to bird populations (Anderson et al. 1999). Directly quantifying the level of potential impacts (e.g., estimated fatalities/mega-watthour) from wind energy development is beyond the scope of this study. Instead, I aim to quantify density, occupancy and flight behavior for the two bird groups mentioned earlier: obligate grassland songbirds and raptors, respectively, predict where impacts may occur, and provide management recommendations to minimize potential impacts. The United States Department of Energy (DOE), through the Office of Energy Efficiency and Renewable Energy Allocation, contracted Texas Tech University to investigate grassland bird patterns of occurrence in the anticipated CREZ in support of DOE’s 20% Wind Energy by 2030 initiative. In cooperation with Iberdrola Renewables, Inc., studies initiated by Wulff (2010) at Texas Tech University were continued at an area proposed for wind energy development and a separate reference site unassociated with wind energy development. I focused on four primary objectives and this thesis is accordingly organized in four separate chapters that address grassland bird density, grassland bird occupancy, raptor flight patterns, and finally I summarize species diversity and composition. The following chapters use formatting from the Journal of Wildlife Management guidelines (Block et al. 2011) with modifications as required by the Texas Tech University Graduate School. 1) I estimate pre-construction bird density patterns using methods that adjust for imperfect detection. I used a distance sampling protocol that effectively accounts for incomplete detection in the field where birds are present but not detected (Buckland et al. 2001). I improved density estimates with hierarchical distance sampling models, a modeling technique that effectively incorporates the detection process with environmental covariates that further influence bird density (Royle et al. 2004, Royle and Dorazio 2008). Covariates included road density and current oil and gas infrastructure to determine the relationship between existing energy development and bird density patterns. Further, I used remote sensing techniques and vegetation field data to investigate how landcover characteristics influenced bird density patterns. I focused species-specific analyses on obligate grassland birds with >70 detections per season namely grasshopper sparrow (Ammodramus savannarum) and horned lark (Eremophila alpestris). Chapter II focuses on hierarchical models that model and describe relationships between grassland bird density and anthropogenic and landscape features. 2) A large number of bird detections (>70) are needed to estimate density using distance sampling and collection of such quantity are often not feasible, particularly for cryptic species or species that naturally occur at low densities (Buckland et al. 2001). Occupancy models operate with far fewer data and are often used as a surrogate for bird abundance when there are fewer detections (MacKenzie and Nichols 2004). I used occupancy models that allow for the possibility of imperfect detection and species abundance to improve estimates of occurrence probability (Royal 2004). I focused species-specific analyses on grassland birds with few detections: Cassin’s sparrow (Peucaea cassinii), eastern meadowlark (Sturnella magna), and upland sandpiper (Bartramia longicauda). Chapter III uses a multi-season dynamic site occupancy model that incorporates bird abundance to better estimate occurrence probability. 3) When I considered the topographic relief of the study sites, the proposed design of the wind facility and its location within the central U.S. migratory corridor, I expanded the study to investigate raptor abundance and flight behavior (Hoover 2002, Miller 2008). I developed a new survey technique that improved the accuracy of raptor flight height estimates and compared seasonal counts and flight heights at the plateau rim and areas further inland. I used counts and flight behaviors to calculate species-specific collision risk indices for raptors based on topographic features. I focused species-specific analyses on raptors with the highest counts: American kestrel (Falco sparverius), northern harrier (Circus cyaneus), red-tailed hawk (Buteo jamaicensis), Swainson’s hawk (Buteo swainsoni), and turkey vulture (Cathartes aura). Chapter IV describes patterns of seasonal raptor abundance and flight behavior and how topography modulates collision risk with proposed wind energy turbines. 4) Finally, for completeness, in Chapter V I summarize morning point count data for all species and provide estimates of relative composition and species diversity with the Shannon-Wiener Diversity Index (Shannon and Weaver 1949).

17 WIND ENERGY↗