Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

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

At least 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
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.↗
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 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 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↗
NASA NTRS2013

Dynamics and Control of a Disordered System in Space

In this paper, we present some ideas regarding the modeling, dynamics and control aspects of granular spacecraft. Granular spacecraft are complex multibody systems composed of a spatially disordered distribution of a large number of elements, for instance a cloud of N grains in orbit, with N greater than 10(exp 3). These grains can be large (Cubesat-size) or small (mm-size), and can be active, i.e., a fully equipped vehicle capable sensing their own position and attitude, and enabled with propulsion means, or entirely passive. The ultimate objective would be to study the behavior of the single grains and of large ensembles of grains in orbit and to identify ways to guide and control the shape of a cloud composed of these grains so that it can perform a useful function in space, for instance, as an element of an optical imaging system for astrophysical applications. This concept, in which the aperture does not need to be continuous and monolithic, would increase the aperture size several times compared to large NASA observatories such as ATLAST, allowing for a true Terrestrial Planet Imager that would be able to resolve exo-planet details and do meaningful spectroscopy on distant world. In the paper, we address the modeling and autonomous operation of a distributed assembly (the cloud) of large numbers of highly miniaturized space-borne elements (the grains). A multi-scale, multi-physics model is proposed of the dynamics of the cloud in orbit, as well as a control law for cloud shape maintenance, and preliminary simulation studies yield an estimate of the computational effort, indicating a scale factor of approximately N(exp 1.4) as a function of the number of grains. A granular spacecraft can be defined as a collection of a large number of space-borne elements (in the 1000s) designed and controlled such that a desirable collective behavior emerges, either from the interactions among neighboring grains, and/or between the grains and the environment. In this paper, each grain is considered to be a highly miniaturized spacecraft which has limited size and mass, hence it has limited actuation, limited propulsive capability, limited power, limited sensing, limited communication, limited computational resources, limited range of motion, limited lifetime, and may be expendable. The modeling and dynamics of clouds of vehicles is more challenging than with conventional vehicles because we are faced with a probabilistic vehicle composed of a large number of physically disconnected vehicles. First, different scales of motion occur simultaneously in a cloud: translations and rotations of the cloud as a whole (macro-dynamics), relative rotation and translation of one cloud member with respect to another (meso-dynamics), and individual cloud member dynamics (micro-dynamics). Second, the control design needs to be tolerant of the system complexity, of the system architecture (centralized vs. decentralized large scale system control) as well as robust to un-modeled dynamics and noise sources. Figure 1, top left, shows the kinematic parameters of a 1000 element cloud in orbit. The motion of the system is described with respect to a local vertical-local horizontal (LV-LH) orbiting reference frame (x,y,z)=F(sub ORF) of origin O(sub ORF) which rotates with mean motion omega and orbital semi-major axis R(sub 0). The orbital geometry at the initial time is defined in terms of its six orbital elements, and the orbital dynamics equation for point O(sub ORF) is propagated forward in time under the influence of the gravitational field of the primary and other external perturbations, described below. The origin of this frame coincides with the initial position of the center of mass of the system, and the coordinate axes are z along the local vertical, x toward the flight direction, and y in the orbit normal direction. The assumptions we used to model the dynamics are as follows: 1) The inertial frame is fixed at Earth's center. 2) The orbiting Frame ORF follows Keplerian orbit. 3) the cloud system dynamics is referred to ORF. 4) the attitude of each grain uses the principal body frame as body fixed frame. 5) the atmosphere is assumed to be rigidly rotating with the Earth. Regarding the grains forming the cloud: 1) each grain is modeled as a rigid body; 2) a simple attitude estimator provides attitude estimates, 3) a simple guidance logic commands the position and attitude of each grain, 4) a simple local feedback controller based on PD control of local states is used to stabilize the attitude of the vehicle. Regarding the cloud: 1) the cloud as a whole is modeled as an equivalent rigid body in orbit, and 2) an associated graph establishes agent connectivity and enables coupling between modes of motion at the micro and macro scales; 3) a simple guidance and estimation logic is modeled to estimate and command the attitude of this equivalent rigid body; 4) a cloud shape maintenance controller is based on the dynamics of a stable virtual truss in the orbiting frame. Regarding the environmental perturbations acting on the cloud: 1) a non-spherical gravity field including JO (Earth's spherical field) zonal component, J2 (Earth's oblateness) and J3 zonal components is implemented; 2) atmospheric drag is modeled with an exponential model; 3) solar pressure is modeled assuming the Sun is inertially fixed; and 4) the Earth's magnetic field is model using an equivalent dipole model. The equations of motion are written in a referential system with respect to the origin of the orbiting frame and the state is propagated forward in time using an incremental predictor-corrector scheme. A representative cloud with varying number of grains is simulated to identify the limitations in computation time as the number of grains grows. We derive a control law to track a desired surface in the ORF (equivalently to maintain a reference cloud shape) by defining an error from a desired surface shape, and designing a control law that is exponentially stable and reduces the tracking error to zero. Figure 1 (top right) shows a comparison of various requirements for simulation of single spacecraft vs. granular spacecraft, indicating the high degree of complexity that needs to be taken into consideration. The ORF components of control force required by one of the grains is, for this particular case, in the micro-Newton range. However, no attempt has been made yet to reconfigure (or re-orient) the cloud configuration internally, for which forces in the milli-Newton level are expected, depending on the time required to do the reconfiguration. Figure 1, bottom, shows the computation time as a function of the number of grains, indicating an order N(exp 1.43) scaling on a 8 Gb, 1067 MHz RAM MacOSX computer with a 3.06 GHz Intel Core 2 Duo processor. With this metric, the same simulation for a system of N=1000 grains would take 5.4 hours, and 146 hours (i.e., 6 days) for a system with N=10,000 grains. Therefore, efficient ways to simulate this complex system, where not only the time scales of natural system dynamics, but also the sampling times of the Guidance, Navigation, and Control are included, remain to be explored. Additional details on the cloud modeling, dynamics, and control will be described in the paper.

simulation↗
DOE OSTI2022

Production of Fischer-Tropsch Synfuels at Nuclear Plants

A case study analysis was performed to evaluate nuclear-powered synthetic fuel production in the midwestern United States (U.S.). A Fischer-Tropsch (FT) fuel synthesis plant design was used as the basis for the analysis. The FT plant design was configured to produce a product slate consisting of diesel fuel, jet fuel, and motor gasoline blend stocks from carbon dioxide (CO 2 ) and hydrogen (H 2 ) feedstocks. The CO 2 feedstock for the FT plant was assumed to be sourced from biorefineries in the region around a Midwest light water reactor (LWR) nuclear power plant (NPP). The analysis specifies that power from the LWR is used to produce H 2 via high-temperature steam electrolysis and to operate the FT synfuel production plant. Capital costs were estimated for the FT plant while capital costs for the electrolysis plant were based on previous Idaho National Laboratory (INL) studies. In addition to labor and maintenance costs for the FT and electrolysis plants, operating costs also include the costs for CO 2 feedstock transport. An analysis was performed to determine the cost of transporting CO 2 from the distributed biorefinery sources to the centralized fuel synthesis plant as a function of the synfuel plant capacity and corresponding CO 2 demand. The primary revenue streams are associated with sales of the synthetic fuel products. The synthetic fuel products will likely follow the same market trends as the conventional fuel products. The synfuel price data was thus based on projections made by the U.S. Energy Information Administration (EIA) 2021 Annual Energy Outlook (AEO) for conventional fuel products minus federal and state taxes, as well as marketing and distribution costs. The economic analysis also considered cases that included and excluded revenues from the 2022 Inflation Reduction Act (IRA) clean hydrogen production tax credit (PTC) of $\$ $3.00/kg for the first ten years of operation. The economic analysis calculated the net present value (NPV) for cases involving steady-state synfuel production for comparison with the NPV for a business-as-usual case in which NPP continues to sell only electric power to the grid. A synfuel production “Reference Case” was considered in addition to sensitivity cases in which the plant capacity, electricity price, and synthetic fuel product prices were perturbed. The synfuel production Reference Case considered a scenario in which the electrolysis and synfuel plants utilized a combined electrical load of 1000 megawatt electrical (MWe) from the LWR with the balance of the LWR power output being sold to the electric grid. The economic analysis suggests that the synfuel production Reference Case evaluated in this analysis would lead to considerable economic potential for near-term deployment of a nuclear-based synfuel production plant. Specifically, the economic analysis suggests that the deployment of a 1000 megawatt (MW) nuclear-powered synfuel plant could result in a NPV increase of approximately $\$ $1.7 billion for a case with no clean synfuel price premium relative to conventional petroleum fuels when accounting for the additional revenues from the 2022 IRA clean hydrogen PTCs of $\$ $3/kg. Sensitivity analysis was performed to evaluate the effect of perturbation of selected model input parameters on the NPV for the synfuel production Reference Case. The sensitivity analysis indicates that the plant capacity has the largest impact on the differential NPV, with a smaller synfuel production capacity resulting in a decrease in revenue when a larger fraction of the power from the NPP is sold to the grid and a smaller fraction of the power is used to produce synthetic fuel products. The synfuel product pricing has the next largest impact on the differential NPV, with lower synfuel prices resulting in decreased NPV from decreased synfuel sales revenue while higher synfuel prices result in increased NPV from increased synfuel sales revenue. Electricity pricing has a smaller effect on the NPV than the fuel sales price since, in the Reference Case, most of the energy from the NPP is used for synfuel production and a smaller amount of the system revenues are associated with electrical power sales. However, the electricity price sensitivity does indicate that the Synfuel Integrated Energy System (IES) would have a greater NPV than the business-as-usual case (e.g., grid power sales only) when electricity market prices are low, suggesting that synfuel production could provide a strategy for decreasing the economic risks to NPPs posed by a loss of revenues attributed to falling electricity market prices.

10 SYNTHETIC FUELS↗
DOE OSTI2020

Catalytic Bio-crude Production in a Novel, Short-Contact Time Reactor

RTI has developed a catalytic fast pyrolysis (CFP) process to produce biomass-based hydrocarbon replacement fuels. From a technology perspective, this advanced biofuels process produces liquid transportation fuels that can leverage capital expenditures in the existing petroleum distribution infrastructure. From a business perspective, this technology does not face potential market limitations caused by oxygenate blend limitations for light duty vehicle applications. A promising catalyst has been developed with excellent potential for deoxygenation to produce a hydrocarbon-rich biocrude that is more suitable for upgrading with traditional hydroprocessing technology. The role of the catalyst is to control the chemistry during biomass pyrolysis to produce a biocrude in a single step that has lower oxygen content and is more thermally stable than conventional biomass fast pyrolysis oil. Catalyst properties are optimized to minimize gas and coke production and improve catalytic deoxygenation and biocrude yields. Results from this project have demonstrated that biocrudes from catalytic fast pyrolysis contain <20 wt% oxygen, with the organic fraction as low as 12% oxygen. These oxygen numbers are > 50% lower compared to traditional fast pyrolysis oils. The biocrude produced from our CFP process can be upgraded to gasoline range hydrocarbons. A very conservative approach for upgrading included a preliminary stabilization step followed by a hydroprocessing step. The oxygen content in the upgraded product was decreased by 97% (to less than 3wt %) from the oxygen content of the biocrude sample. Additionally, the acid content was reduced by 99% from the acid content of the biocrude sample. In an optimized integrated process, the CFP step and the hydroprocessing (upgrading) step need to be carefully balanced to maximize the input biomass energy that is recovered in the finished biofuels. Hydrogen utilization and carbon recovery during both catalytic biomass pyrolysis and biocrude upgrading need to be optimized to maximize energy recovery in the finished biofuel. A nominal 1 ton per day (TPD) catalytic biomass pyrolysis system was designed based on the laboratory experimental data. The system was designed to achieve the short residence times (0.5–2 sec) and high heat transfer rates required for maximum liquid biocrude yields while optimizing process integration to maintain catalyst activity by continuous regeneration. The catalytic biomass pyrolysis reactor is a continuously circulating single-loop transport reactor design that is flexible enough to allow sensitivity studies around temperature, residence time, biomass feed rate, and catalyst-to-biomass ratio (i.e., catalyst circulation rate) for process optimization. For a nominal 100 lb/hr feed rate of dry biomass, the products include ~25 lb/hr char and coke, ~20 lb/hr permanent gases, and ~55 lb/hr condensable vapors. Pyrolysis temperatures are 350–600°C, and the system is operated close to ambient pressure. Production of the catalyst developed within the project was scaled up to provide enough material (200kg) for bench-scale testing. The 1TPD bench-scale unit was fabricated, installed and commissioned with this project. Biocrude samples from wood were produced in this unit with the CFP catalyst with material balance > 90%. The intermediate biocrude yields and chemical composition compare well with previous laboratory results and a hydroprocessing strategy for upgrading the resulting biocrude was developed. The experimental results obtained within this project were used to develop a process model for a 2000 TPD commercial facility. Detailed flow sheets were developed for all major unit operations in the process and provided a basis for an economic analysis to develop capital and operating costs. Estimates of the biofuel production costs were used to develop internal rates of return for the yields and compositions measured for the current state of technology at the end of this project and for the yields anticipated for an optimized (nth plant) scenario as additional process improvements are made.

09 BIOMASS FUELS↗
DOE OSTI2022

WISP: Watching grid Infrastructure Stealthily through Proxies (Final Technical Report)

The complex interdependencies of cyber systems (sensors and communications), physical grids and associated electricity market operations make protecting electric power grids a significant challenge. The energy sector is constantly under new, targeted, advanced and dangerous cyber-attacks that have the potential to result in the loss of human life. These threats are further exacerbated by our need to modernize the grid. One focus of cyber security research in smart grids is the securing of the SCADA system through advanced intrusion detection systems (IDS) and bad data detection algorithms in state estimation. These methods either require full knowledge of the system topology and parameters or fail to understand the physical behaviors under attack. WISP (Watching grid Infrastructure Stealthily through Proxies) is designed to provide additional protection to the power grid using only publicly available data. In particular, WISP exploits the spatio-temporal nature of the real time locational marginal prices (LMPs), in conjunction with other information such as bids, weather, outages and load data to analyze anomalous power pricing behaviors and then correlate those observations to localize regions of interest and identify potential cyber events. WISP is non-intrusive as the tool is deployed as a service in the Cloud or on premise and provides reliable information to system operators for enhanced situational awareness, without impeding energy delivery functions. The WISP technology comprises three modules: the data-driven anomaly detection core, the vulnerability and risk analysis and the root cause analysis. The data-driven anomaly detection core performs the tasks of feature selection, anomaly detection and attack region localization. The vulnerability and risk analysis module provides system level information of the vulnerable variables and times, assisting the operators in selecting monitoring and protection nodes. The root cause analysis module takes the detection results and identifies potential operational conditions that contribute to the detected anomalies. In Phase I, we have demonstrated the feasibility and effectiveness of WISP. We developed a realistic electricity market simulator capable of generating normal and attack market data under various operational conditions. We developed a series of cyber-attack detection and analysis algorithms and evaluated them under multiple data sources. Finally, we integrated all modules into an end-to-end software, providing functions for data management, data analytics and visualization. Specifically, we have achieved: (i) real-time data acceptance from external utility interfaces with >99% acceptance rate; (ii) high performance anomaly detection algorithms with >98% detection accuracy and <0.1% false alarm rate; and (iii) ultra-low computing delay <50 milliseconds. Additionally, our team developed algorithms to identify the vulnerable variables in electricity market operations and root cause analysis functions to identify major contributors to the price spikes. These ancillary modules are necessary when deploying WISP in real world industry environment. In Phase II, we have demonstrated the effectiveness of WISP software on realistic largescale power systems. We performed red team testing for the Phase I WISP software and identified software vulnerabilities and implemented corresponding mitigation solutions. We adapted the electricity market simulator for the Texas synthetic 2000-bus system and generated datasets for the false data injection attacks. We created database and visualization interfaces for the Texas system and the ISO New England system. We performed software optimization in terms of operation efficiency, computing speed and detection accuracy. Finally, we tested the software on the Texas system and the ISO New England system and evaluated the detection performance. Overall, we achieved above 89% detection rate, below 3% false alarm rate and below 37 seconds of end-to-end detection delay.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRSDate not supplied

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand [1] while simultaneously maintaining one of the safest forms of transportation [2], [3]. One of the reasons for this success is the ability of the system and the operators to adapt and accommodate to situations that routinely disrupt optimal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators ability to control. These factors can lead to states where automation is unable to properly handle these issues and therefore air traffic controllers and pilots have to intervene, ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions this can be viewed as an increase in complexity. The reasoning for this is because humans are now required to make tactical decisions in response to external factors, resulting in a departure from the strategic plan where operations would be more efficiently managed. Human operators control airspace complexity under rigid regulations that are constantly changing. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. There has been past work that devised airspace complexity metrics in commercial aviation and related these metrics to controller workload (e.g., [4],[5]). The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic, including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly [6]. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to our proposed effort that identifies such contributing factors or precursor patterns. To define the scope we are proposing to measure complexity from the viewpoint of the Terminal Radar Approach Control Facilities (TRACON) controller’s perspective. In particular we are analyzing arrivals into KSFO. With safety as the top concern for airspace operators, it is important to recognize that as density and heterogeneity grow, the focus of the system will change. Times of the day when the airspace has low density and heterogeneity, the flights will follow more efficient paths where the aircraft move on established routes that are more or less directly to the destination. However, when density and heterogeneity increases, the system will begin changing focus to avoiding conflicts and collisions and route the flights in a more flexible way. Higher flexibility requires more communication and coordination between controllers and pilots which the current automation is unable to handle. This paper proposes a novel approach that monitors airspace complexity at multiple scales, uses a Machine Learning-based tool that predicts when operations will transition to a regime of greater complexity, and identifies actions that can reduce the complexity while still maintaining efficient and safe operations. We demonstrate our proposed approach using data from multiple complementary sources. This includes, but is not limited to: historical aircraft surveillance data from NASA’s Sherlock Data Warehouse [7], METAR weather data, and airport configuration data from Aviation System Performance Metrics (ASPM). The surveillance data flight paths are sampled at a variable sample rate — increasing as the aircraft approaches the airport. This is due to how Sherlock manages flight track stitching between different radar facilities which have different sampling rates. The weather and performance data are logged at defined intervals throughout the day at a courser refresh rate. In addition to the logged data and metrics, we leverage pre-defined Standard Terminal Arrival Routes (STARs) procedures to characterize the path of each flight. Each flight files for one of these routes in the flight plan well before entering the terminal airspace, and approximately follows the route until it leaves the STAR, typically on the final fix of a runway transition. However, most flights do not always fly the full STAR procedure to completion [8], but the majority do adhere to the fixes within the common route of the procedure. Our approach leverages fixes in the common route of each of the STARs to build a reference path to the airport. This allows us to characterize the flight paths in what we are defining as the “maneuvering area” (the airspace between the STAR and before the flight is lined up on the runway’s final approach) to determine how off nominal the flights are to calculate its complexity score. Determining airspace complexity is a concept that does not have a concrete answer. In designing this metric, we consider what increases the workload for the air traffic controllers. Consequently more specialized vectoring maneuvers results in higher workload. Accordingly, we start with a theory: each flight has a direct path it takes from the STAR’s common route to the final approach’s outer marker fix for the flight’s landing runway. It is important to note that the direct path is only used as a reference. If the majority of the flights have a large consistent offset as compared to other routes it does not necessarily mean that those flights have higher complexity. We are merely building a distribution based on this direct path for that particular STAR and runway pair to determine the normal mode of operations for that route. Flights that are in the upper tail of these distributions will result in higher complexity scores and flights that fly in the median will represent the normal mode of operations and therefore will have lower complexity scores. Since flights following each STAR route take different paths to the airport, we have a different distribution for each STAR route and therefore can model these distributions to compute a complexity score from their respective normalized distributions. To evaluate the effectiveness of our proposed airspace complexity metric we will compare against an established approach based on trajectory clustering [9]. This unsupervised learning technique consists of the following steps: (1) identify the general maneuvering areas (waypoints) by performing $\kappa$-means or DBSCAN clustering on locations where aircraft frequently turn based on the surveillance radar track data, (2) map flight trajectories onto sequences of waypoints, and (3) cluster the sequences based on their common subsequences. From a high-level perspective, this baseline model learns nominal operations in the airspace through the sequence of waypoints that are representative of where aircraft change direction and defines deviations from the nominal operations as “complex.” Therefore, more deviations from the nominal operations correspond to higher complexity values. For our validation, we re-implemented this technique and tune model hyper-parameters to correctly detect waypoints for the arrival traffic into the San Francisco bay area. We will compute the complexity measure over a one-year period using our proposed technique as well as the baseline. Our validation will be based on each technique’s ability to detect a set of undesirable outcomes (e.g., go-arounds, holding patterns, average time in the airspace, etc.). Since our current complexity metric is derived from the offset from the direct reference path, it’s important to understand what causes these offsets. In many of the flights with high offset distance, flights performing holding patterns and S turns can be observed. These maneuvering tactics are utilized to add distance between the aircraft and the destination runway to prevent multiple flights from having conflicting arrival times. In order to predict a rise in complexity (or the precursor to complexity), it’s necessary to be able to identify these potential conflicts (which in turn, result in higher offsets). To do this, we define a “representative flight” for each STAR route and runway pair. This flight is approximately the path the flight would take if there was a clear path with no other flights in the airspace — including the time remaining to the airport. We first identify the flights for a given STAR runway pair using the offset to the reference path distributions that fall between the 44-55 percentiles. This yields the flights that conform to the most normal mode of operation. Each of these flights is partitioned based on the percent complete from the entry point into the maneuvering areas from 0\% – 100\% complete. Then for each percent “bin”, we take the median value of the flight’s latitude/longitude coordinates, airspeed, and (non causal) time remaining to the airport to construct a lookup table for each percent complete bin on a given route. As a flight enters the maneuvering area, we can find the estimated arrival time of a flight to the airport by finding the closest point to the representative path’s percent complete bin (relative to the flight’s current position at any snapshot in the airspace) and therefore retrieve the corresponding remaining time left on the “representative path”. We assume that the flight will follow the representative path to completion when deriving these estimates. We can then compare these estimated arrival times against other flights for the same snapshot in time to identify potential conflicts. If more flights are estimated to arrive within a tolerance window than there are runways available, then we have a potential conflict. We can use this derived measure along with other factors expected to add disruption to the operation such as weather and runway configuration changes as an input to machine learning tools to detect precursors that increases in our complexity measure. This novel method will assist in uncovering insights into the contributing factors that lead to increased complexity that may allow for in-time responses to avoid reaching a high complexity state in the airspace.

complexity↗
NASA NTRS2017

Validation of a Manually Oscillating Chair for In-The-Field Assessment of Dynamic Visual Acuity on Crewmembers Within Hours of Returning From Long-Duration Spaceflight

Long-duration spaceflight results in sensorimotor adaptations, which cause functional deficits during gravitational transitions, such as landing on a planetary surface after long-duration microgravity exposure. Both the vestibular system and the central nervous system are affected by gravitational transitions. These systems are responsible for coordinating head and eye movements via the vestibulo-ocular reflex (VOR) and go through an adaptation period upon exposure to microgravity. Consequently, they must also re-adapt to Earth's gravitational environment upon landing. This re-adaptation causes decrements in gaze control and dynamic visual acuity, with crewmembers reporting oscillopsia and blurred vision caused by retinal slip, or the inability to keep an image focused on their retina. This is thought to drive motion sickness symptoms experienced by most crewmembers following landing. Retinal slip can be estimated by dynamic visual acuity (DVA); visual acuity while in motion. Previously, DVA has been assessed in the laboratory where subjects walked at 6.4 km/hr on a motorized treadmill. Using this method, Peters et al. (2011) found that DVA is worsened in astronauts by an average of 0.75 eye-chart lines one day after landing. However, it is believed that re-adaptation occurs quickly and that DVA might be worse immediately upon re-exposure to a gravitational environment. Since many crewmembers are unable to walk safely upon landing, it was necessary to develop a method for replicating the vertical head movements associated with walking. In addition, the use of a chair to imitate the head displacement caused by walking isolates eye-head interactions without allowing for trunk and lower-body compensation, as seen with treadmill walking (Mulavara & Bloomberg 2003). Therefore, a modality for assessing DVA in the field within a few hours of landing was developed. In this study, we validated the ability of a manually operated oscillating chair to reproduce the oscillatory frequency of walking on a treadmill. Healthy non-astronaut subjects (n=14) participated in one test session and completed three static (seated) and three dynamic (walking/oscillated) visual acuity tests. DVA was assessed using a motorized treadmill, an automated oscillating chair, and a manually operated chair, both developed in the Neuroscience Laboratory at JSC. The automated chair was motor-driven and set to oscillate vertically at 2 Hz with a vertical displacement of +/- 5 cm to simulate vertical translation while walking. The manually operated chair was oscillated vertically by a test operator to the beat of a metronome at 120 beats/min (2 Hz) and a vertical displacement of approximately +/- 5 cm. As the subject was oscillated, they were asked to discern the direction gap of Landolt-C optotypes of varying sizes and verbally reported the direction while an operator recorded their response using a gamepad. Subjects were outfitted with accelerometers (sampling rate = 128 Hz) on their head, trunk and lumbar spine. A fast Fourier transform was performed on the vertical trunk acceleration to compare the peak and spread of the distribution of oscillation frequencies for each oscillating condition. The spread of the frequency distribution for the manual chair was not significantly different from either the treadmill or the automated chair. However, all three conditions had similar non-zero standard error values, suggesting a variance in head movement frequency which may affect DVA. The average oscillation frequency of the manual chair (1.85 Hz) was significantly different (α=0.05) from that of treadmill walking (2.24 Hz), but not significantly different from that of the automated chair (1.85 Hz) and all three conditions had small standard errors (SEM = 0.04, 0.06, and 0.08 Hz for manual, treadmill, and automated respectively). This implies that both chairs oscillate at a frequency below that of treadmill walking, but are comparable to each other and reproducible across sessions. Additionally, DVA scores did not vary significantly across conditions. The smaller spread values of the oscillating chairs' frequencies indicated mitigation of variation induced by locomotor strategies, which enables better examination of the issue of VOR adaptation. Furthermore, due to the deconditioned state of crewmembers in the initial hours after landing, it is easier to transport a manual bouncing chair into the field and safer to perform a vision test while seated in a chair versus walking on a treadmill. Therefore, the manually oscillating chair has been deemed to meet and exceed the DVA testing capabilities previously obtained by treadmill walking.

Kreutzberg, G. A.↗
DOE OSTI2020

Distribution Development for Residual Inventory at the New York West Valley Site - 20400

The New York State Energy Research and Development Authority (NYSERDA) is the owner of the Western New York Nuclear Service Center (WNYNSC), a 1,351 ha site located approximately 48 km south of Buffalo, New York. In 1962, Nuclear Fuel Services, Inc. (NFS) entered into Agreements with the Atomic Energy Commission and New York State to construct the first commercial reprocessing plant of nuclear fuel in the United States. NFS, a private company, built and operated the spent fuel reprocessing plant and waste disposal facilities, processing 640 Mg of spent nuclear fuel from 1966 to 1972 under an Atomic Energy Commission license. Nuclear fuel reprocessing operations ended in 1972 and never reopened, leaving behind radioactive and chemical wastes. Operations led to contamination in a number of facilities and locations. Some of that contamination has migrated from waste disposal zones to other layers, formations, and features on and off the WNYNSC. Phase I decommissioning activities are ongoing and involve the removal of a number of areas and structures that have been associated with contamination. The purpose of this work is to outline the approach for characterizing contamination not associated with disposed wastes, contaminated structures, or specific releases. In this work, the term, residual radiological activity, is used to describe environmental contamination that exists subsequent to the completion of Phase I decommissioning activities, that is not associated with disposed wastes, contaminated structures, or specific releases. Contamination from the Site was quantified relative to data that characterize the concentrations of radionuclides that exist in background. Background concentrations are those present in the area but having no influence from Site related activities. The existence of residual radiological activity that is elevated relative to background has the potential to contribute to future risks to human health and the environment. As a consequence, the residual inventory information is used to inform the West Valley Probabilistic Performance Assessment (PPA) model to characterize potential future risks to human health and the environment. The centralized West Valley Data Management System (DMS) was the source of information for the data assembled in this analysis. The DMS is a fairly large compilation consisting of thousands of records from investigation studies, with sample dates ranging from 1990 to present. Samples from monitoring wells, boreholes, geoprobe studies, surface water, surface soils, storm water outfalls, ventilation stack filters, plant and animal tissues, and more are included in the DMS. Results are typically reported in units of activity per unit volume. For the purpose of the analyses presented here, all results were converted into consistent units of pCi per unit volume. Since 1990, data have been collected from various locations across the WNYNSC at different times with varying frequency over the course of several decades. As a consequence, a number of potential issues can arise with respect to the assembly of a dataset that is deemed adequate for the characterization of residual radiological activity. These issues were assessed and resolved to the extent possible through careful consideration of the properties of the distributions. The intent was to use data which characterize the current state of the Site. Radionuclides can be designated to one of several groups depending on their origin. In this work the groups considered were 1) Naturally Occurring Radioactive Material (NORM), 2) fallout, and 3) Other (including power plant, medical research, etc). This grouping is a useful construct with respect to the interpretation of fixed laboratory results. For example, NORM radionuclides that exist within a decay chain should have approximately equivalent distributions of concentrations if they are representative of background conditions. Insights such as these can be used as a check to identify sample results that need to be further investigated or omitted due to issues associated with reported results from fixed laboratory analyses. This type of analysis provided a foundation for the assessment of the adequacy of sample results for use in subsequent components of an assessment. The general process for the assessment of residual radiological contamination at the Site consists of a sequence of several steps. First, for each analyte, several statistical tests were performed to assess the weight of evidence against the null hypothesis that the mean of the distribution of concentrations was equal to zero. If the mean of the distribution of concentrations for a given radionuclide was not found to be greater than zero, then it was removed from consideration as a component of the residual radiological contamination. If there was significant evidence to reject the hypothesis of the mean being equal to zero, the second step was to compare the distribution of the data from the Site to that of the corresponding background. A suite of tests was used to compare the distributions of the site and background data. The results of these tests were collectively used to determine if site data are elevated relative to background. The third step was to develop distributions using a Bayesian framework to characterize the distribution of mean of the increment present above background for each of the radionuclides. The Bayesian model implemented allowed for the comparison of site-specific records to background concentrations to better approximate contamination attributed to the Site. A final screening step was employed for radionuclides that exceed background. This screening step compared 95% upper confidence limits (UCLs) from the increment distribution developed in the previous step to the risk screening levels. This approach yields a list of analytes that were determined to be elevated relative to background.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Quantifying the Solar Energy Resource for Puerto Rico

After Hurricane Maria, multiple U.S. Department of Energy laboratories studied the state of the electric grid in Puerto Rico and analyzed grid resilience and grid integration of renewable energy. As part of the work done at the National Renewable Energy Laboratory, researchers created new solar resource data, conducted a technical potential and supply curve analysis, and studied the interannual variability of the solar resource. A new methodology was developed to downscale solar resource data from the National Solar Radiation Data Base (NSRDB) from a 4-km x 4-km spatial and 30-minute temporal resolution to a 2 km x 2 km and 5-minute resolution. This methodology primarily used simple physical principles to develop high-resolution cloud properties which were then used to compute solar radiation. The high-resolution datasets were validated against ground measurements and the error metrics were found to be similar to the original lower resolution dataset. Using 20 years of downscaled data from the NSRDB multi-year capacity factors for photovoltaics (PV) were developed for both single-axis tracking and fixed latitude-tilt configurations. Use of the multi-year data provides the ability to understand variability in capacity factors due to variability in weather over a long period of time. For Puerto Rico the coastal regions were found to have significant higher capacity factors than inland. Using land-use and terrain information a technical potential analysis was conducted for Puerto Rico. This analysis restricted single PV plant development to a maximum of 100 MW nameplate capacity. The nameplate capacity for each municipality were then determined. Based on our assumptions, 56 of the 78 total municipalities of Puerto Rico contain some level of solar capacity. Most of the interior municipalities did not have any capacity because of the geographic exclusions used in this study. The lowest capacity for a PV plant observed in a municipality was 10 MW. The maximum capacity within a county was 2,000 MW. Further a supply curve analysis was conducted by taking the results of the technical potential and quantifying system and transmission costs. The levelized cost of energy (LCOE) was calculated for each theoretical PV plant site, and the levelized cost of transmission was added to the LCOE to produce a total cost estimate for each site. The results of the supply curve analysis allow for a relative comparison of the cost for integrating new PV capacity into the grid. This analysis indicates that cheaper total LCOE sites tend to be larger in capacity. The total capacity in this study was found to be far beyond the maximum peak load for the island. However, this study does not consider the economic and market potential for development. The cumulative capacity presented in this study assumes that the best locations are developed first and ignores the complex decision paths for new power plant development. Therefore, this analysis can only be treated as illustrative. Finally this study investigates the impact of inter-annual variability of resource using a variety of metrices including probability of exceedance and variation in capacity factor and LCOE. This study demonstrates that the capacity factor or LCOE could vary by over 10% year to year. This clearly indicates the risks involved in using any particular year of data and clearly points to the use of multi-year data to reduce some of the risks related to variability in weather.

14 SOLAR ENERGY↗
DOE OSTI2026

Generator Interconnection Costs to the Transmission System in non-ISO Balancing Authorities [Slides]

Electric transmission system operators—including Independent System Operators (ISOs), Regional Transmission Organizations (RTOs), and utilities—require proposed power plants to undergo a series of interconnection studies before connecting to the grid. These studies assess what transmission upgrades or new infrastructure may be necessary and assign the associated costs to the project. Lawrence Berkeley National Laboratory has compiled, aggregated, and cleaned interconnection cost data, originally for ISOs/RTOs, and now for five non-ISO Balancing Authorities: PacifiCorp, Bonneville Power Authority, Duke Energy Progress, Duke Energy Carolinas and Duke Energy Florida. Insufficient transparency in interconnection cost data may contribute to rapidly expanding interconnection queues, with active queue capacities tripling between 2020 and 2024 in the studied BAs. Most projects withdraw after receiving high interconnection cost estimates. Interconnection costs have increased since the early 2000s, with average costs for "complete" projects reaching $194/kW between 2018 and 2024. Active queue projects and withdrawn projects incur substantially higher costs, primarily due to rising network upgrade costs. Recent interconnection costs in non-ISO balancing authorities are higher than in ISO regions, potentially due to a greater willingness to pay among developers. Utility-scale solar, wind, and storage projects have interconnection costs that exceed those for natural gas. However, when focusing on projects that do not withdraw from the queue, the interconnection costs for these technologies are more similar to natural gas projects. Other key findings include: (1) Larger generation projects benefit from lower proportional interconnection costs, (2) capacity transmission service (NRIS) often requires additional network investments, and (3) projects with high network upgrade costs are often clustered geographically. The dataset includes results from 2,104 interconnection studies conducted between 2000 and 2024, covering projects that are operational, withdrawn, or still progressing through the study process. The Excel file contains (a) the complete project-level interconnection cost dataset, and (b) seven additional tabs summarizing cost metrics across dimensions such as time, market structure, cost category (point of interconnection vs. broader network upgrades), fuel type, service type (ERIS vs. NRIS), generator size, and geography.

24 POWER TRANSMISSION AND DISTRIBUTION↗