SEARCH · Search NASA

Results for “Distribution systems state estimation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 991 records · Page 55

All sources NASA NTRS DOE OSTI

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 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 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↗

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 OSTI2015

Port Graham Community Building Biomass Heating Design Project

Native Village of Port Graham completed preconstruction activities to prepare for construction and operations of a cord wood biomass heating system to five or more community buildings in Port Graham, Alaska. Project Description Native Village of Port Graham (NVPG) completed preconstruction activities that pave the way towards reduced local energy costs through the construction and operations of a cord wood biomass heating system. NVPG plans include installation of a GARN WHS 3200 Boiler that uses cord wood as fuel source. Implementation of the 700,000 Btu per hour output biomass community building heat utility would heat 5-community buildings in Port Graham, Alaska. Heating system is estimated to displace 85% of the heating fuel oil or 5365 gallons of fuel on an annual basis with an estimated peak output of 600,000 Btu per hour. Estimated savings is $15,112.00 per year. The construction cost estimate made to install the new biomass boiler system is estimated $251,693.47 with an additional Boiler Building expansion cost estimated at $97,828.40. Total installed cost is estimated $349,521.87. The WHS 3200 Boiler would be placed inside a new structure at the old community Water Plant Building site that is controlled by NVPG. Design of the new biomass heat plant and hot water loop system was completed by Richmond Engineering, NVPG contractor for the project. A hot water heat loop system running off the boiler is designed to be placed underground on lands controlled by NVPG and stubbed to feed hot water to existing base board heating system in the following community buildings: 1. Anesia Anahonak Moonin Health and Dental Clinic 2. Native Village of Port Graham offices 3. Port Graham Public Safety Building/Fire Department 4. Port Graham Corporation Office Building which also houses the Port Graham Museum and Head Start Center 5. North Pacific Rim Housing Authority Workshop/Old Fire Hall Existing community buildings fuel oil heating systems are to be retro-fitted to accommodate hot water from the proposed wood-burning GARN Boiler, once installed, and rely on the existing fuel oil-fired hot water heating equipment for backup. The boiler would use an estimated 125 bone dry tons, equivalent to 100 cords, woody biomass feedstock obtained from local lands per year. Project would use local labor as described in the Port Graham Biomass Project, report completed by Chena Power, Inc. and Winters and Associates as part of the in-kind support to the U. S. Department of Energy (DOE) project for work on a project for State of Alaska’s Alaska Energy Authority (AEA). NVPG will likely initiate operations of the biomass boiler system even though several operational variations were studied. Obtaining the fuel source could be done by contractors, PGVC employees, or NVPG employees. Feeding the system would likely be done by NVPG employees. A majority of the buildings heated would be owned by NVPG. The PGVC office would be heated as well as the Old Fire Hall used as a workshop and storage area for North Pacific Rim Housing Authority. One methodology studied to charge for cost of utilizing the community building biomass system would use a percentage of use of hot water generated by the biomass hot water system based on past heating oil usage in relation to all buildings heated by biomass hot water. The method is better described in the Port Graham Biomass Project report. Fuel source agreements have been drafted to enter into agreements with area landowners. One Native allotment owner has asked Chugachmiut Forestry to begin a timber sale process to sell timber off her lands, specifically wind thrown timber that was determined to be of sufficient quantity to supply to the proposed biomass heating system for approximately 5-years. On NVPG’s behalf, Chugachmiut has presented to PGVC three different documents, attached, that could lead to a sale of woody biomass fuel for the project for up to 25-years, the expected life of the project. PGVC has signed a letter of intent to negotiate a sale of woody biomass material April 30, 2015. Chugachmiut Forestry has conducted two different field forest measurements of Native allotment lands and PGVC forest and timber lands. Lands deemed road accessible for biomass harvest were analyzed for this project. Forestry then conducted three different analyses and developed two reports to determine forest biomass on a tons per acre basis in addition to timber volume measurements taken for forest management purposes. Permits required were limited. For the biomass building, the Kenai Peninsula Borough did not require a permit. State of Alaska, Department of Public Safety, Division of Fire and Life Safety requires a plan review for fire and life safety requirements called an application for Fire and Life Safety Plan Review that would require a registered design professional to sign the document. State of Alaska State Forest Practices Act is required to be followed for any timber sale or harvest. This Act also requires consultation with Alaska Department of Fish and Game when operations are in close proximity or cross anadromous waters. Native allotment lands require following U. S. Bureau of Indian Affairs timber sale contracting process and approval.

09 BIOMASS 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 OSTI2018

Field Demonstration of The Krauklis Seismic Wave In A Novel Mva Method For Geologic CO 2 Storage

The Energy & Environmental Research Center (EERC) conducted a feasibility study of a new, low-impact geophysics monitoring method designed to incrementally track the saturation front of CO 2 injected into a reservoir. The method leverages a new way of transmitting energy from the surface to the reservoir and employs a new subsurface signal called the Krauklis wave (K-wave) and other guided waves for seismic monitoring. The EERC partnered with Seismos, Inc., the company that holds the exclusive technology license for this method, to advance development of K-wave technology, and Denbury Onshore, LLC, to field-test the system. The system was field-tested at Denbury’s Bell Creek oil field, currently undergoing CO 2 -based enhanced oil recovery (EOR) in southeastern Montana. The K-wave system uses a wellhead-mounted source to propagate energy down the wellbore in the form of a tube wave and through the well perforations into the reservoir. Part of the energy transmitted through the perforations into the reservoir converts to slow moving K-waves and/or other guided waves that are expected to propagate through the formation to neighboring wells. The propagated signal would induce tube waves in the neighboring wells that return to the surface to be recorded with a wellhead-mounted receiver. The propagation time of the signal is directly affected by the presence of CO 2 between the wells, so the difference in transit time of the returning signal between pre- and postinjection data sets could be used to help define the lateral extent of a CO 2 saturation front. Objectives of the study were to demonstrate, validate, and evaluate the K-wave method’s ability to monitor the morphology of injected CO 2 and estimate saturation distribution. Results from the demonstration were to be compared to conventional 4-D surface seismic methods, dynamic reservoir simulations, and CO 2 breakthrough data to validate the technology and provide useful information that may result in improvements. Prior to mobilizing for the full-scale K-wave baseline survey, a field test was conducted to better understand wellhead hardware configurations and test the response from the Bell Creek reservoir using three wells. Useful operational information resulted but signal from the first-generation gas release/gas injection source was not recorded in neighboring wells. After the initial field test, numerical modeling of the K-wave method was initiated using parameters specific to the Bell Creek Field. Modeling results indicated that the fundamental concepts behind the K-wave method exhibited problems and revealed a significantly greater attenuation of the signal at the wellbore perforations than expected, totaling on the order of 10-8 from source well to receiver well over a 400-m well spacing. Modeling also showed that greater source frequency bandwidth was needed to produce guided waves within a 30-foot-thick reservoir layer such as that at Bell Creek and that reservoir layer seismic velocity would need to be slower than the velocities of the bounding layers in order to act as a wave guide opposite to what exists at Bell Creek. Modeling results prompted a second field test with a stronger, second-generation displacement source developed in the interim. Concurrent with the second field test, a 3-D seismic survey was acquired to serve as a baseline for project validation and to serve as a monitor data set for the portion of the project area where CO 2 injection had already commenced. The new source was tested using a four-well geometry and failed to produce a detectable signal, which appeared to confirm modeling results and indicated that the Bell Creek Field is not an ideal site for the K-wave test. It was concluded that a source with sufficient strength and bandwidth could not be built to overcome the factors discovered by the modeling within a reasonable amount of time. Several other sites undergoing CO 2 -based EOR were also evaluated; however, none of these sites was amenable because of geology, limitations of source and receiver hardware, or the advanced state of the CO 2 flood. In addition, two other configurations of the source and receiver layout were considered; however, these design modifications either changed the scope of the project or did not meet the goals of the Carbon Storage Program. In consultation with the U.S. Department of Energy (DOE) Project Manager and Technology Manager, a no-go decision was invoked to forego additional expenditures as the technology would not meet current project objectives. Project activities to date have moved the technology forward by developing modeling equations that can be applied elsewhere, advanced the state of wellhead sensors and sources, and spurred improvements to field operations. The 3-D seismic survey in the project area and the dynamic reservoir simulations done as part of this study have provided value for ongoing and future synergistic activities at Bell Creek related to DOE carbon storage goals.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2017

Investigation of Microbunching Instabilities in Modern Recirculating Accelerators

Particle accelerators are machines to accelerate and store charged particles, such as electrons or protons, to the energy levels for various scientific applications. A collection of charged particles usually forms a particle beam. There are three basic types of particle accelerators: linear accelerators (linac), storage-ring (or circular) accelerators, and recirculating accelerators. In a linac, particles are accelerated and pass through once along a linear or straight beamline. Storage-ring accelerators propel particles around a circular track and repetitively append the energy to the stored beam. The third type, also the most recent one in chronology, the recirculating accelerator, is designed to accelerate the particle beam in a short section of linac, circulate the beam, and then either continue to accelerate for energy boost or decelerate it for energy recovery. The beam properties of a linac machine are set at best by the initial particle sources. For storage rings, the beam equilibria are instead determined by the overall machine design. The modern recirculating machines share with linacs the advantages to both accelerate and preserve the beam with high beam quality, as well as efficiently reuse the accelerating components. The beamline design in such a machine configuration can however be much more complicated than that of linacs. As modern accelerators push toward the high-brightness or high-intensity frontier by demanding particles in a highly charged bunch (about nano-Coulomb per bunch) to concentrate in an ever-decreasing beam phase space (transverse normalized emittance about 1 μm and relative energy spread of the order of 10^-5 in GeV beam energy), the interaction amongst particles via their self-generated electromagnetic fields can potentially lead to coherent instabilities of the beam and thus pose significant challenges to the machine design and operation. In the past decade and a half, microbunching instability (MBI) has been one of the most challenging issues for such high-brightness or high-intensity beam transport, as it would degrade lasing performance in the fourth-generation light sources, reduce cooling efficiency in electron cooling facilities, and compromise the luminosity of colliding beams in lepton or lepton-hadron colliders. The dissertation work will focus on the MBI in modern recirculating electron accelerators. It has been known that the collective interactions, the coherent synchrotron radiation (CSR) and the longitudinal space charge (LSC) forces, can drive MBI. The CSR effect is a collective phenomenon in which the electrons in a curved motion, e.g. a bending dipole, emit radiation at a scale comparable to the micro-bunched structure of the bunch distribution. The LSC effect stems from non-uniformity of the charge distribution, acts as plasma oscillation, and can eventually accumulate an amount of energy modulation when the beam traverses a long section of a beamline. MBI can be seeded by non-uniformity or shot noise of the beam, which originates from granularity of the elementary charge. Through the aforementioned collective effects, the modulation of the bunch sub-structure can be amplified and, once the beam-wave interaction formed a positive feedback, can result in MBI. The problem of MBI has been intensively studied for linac-based facilities and for storage-ring accelerators. However, systematic studies for recirculation machines are still very limited and form a knowledge gap. Because of the much more complicated machine configuration of the recirculating accelerators than that of linacs, the existing MBI analysis needs to be extended to accommodate the high-brightness particle beam transport in modern recirculating accelerators. This dissertation is focused on theoretical investigation of MBI in such machine configuration in the following seven themes: (1) Development and generalization of MBI theory The theoretical formulation has been extended so as to be applicable to a general linear beamline lattice including horizontal and vertical transport bending elements, and beam acceleration or deceleration. These featured generalizations are required for MBI analysis in recirculation accelerators. (2) Construction of CSR impedance models In addition to the steady-state CSR interaction, it has been found that the exit transient effect (or CSR drift) can even result in more serious MBI in high-brightness recirculation arcs. The onedimensional free-space CSR impedances, especially the exit transients, are derived. The steady-state CSR impedance is also extended to non-ultrarelativistic beam energy for MBI analysis of low-energy merger sections in recirculating accelerators. (3) Numerical implementation of the derived semi-analytical formulation This includes the development of a semi-analytical Vlasov solver for MBI analysis, and also benchmarking of the solver against massive particle tracking simulations. (4) Exploration of multistage amplification behavior of CSR microbunching development The CSR-induced MBI acts as an amplifier, which amplifies the sub-bunch modulation of a beam. The amplification is commonly quantified by the amplification gain. A beam transport system can be considered as a cascaded amplifier. Unlike the two-stage amplification of four-dipole bunch compressor chicanes employed in linacs, the recirculation arcs, which are usually constituted by several tens of bending magnets, show a distinguishing feature of up to six-stage microbunching amplification for our example arc lattices. That is, the maximal CSR amplification gain can be proportional to the peak bunch current up to sixth power. A method to compare lattice performance has been developed in terms of gain coefficients, which nearly depend on the lattice properties only. This method has also proven to be an effective way to quantify the current dependence of the maximal (5) Control of CSR MBI in multibend transport or recirculation arcs The existing mitigation schemes of MBI mostly aim to linac-based accelerators and may not be practical to the recirculating accelerator facilities. Thus a set of conditions for suppression of CSR MBI was proposed and examined for example lattices from low (~100 MeV) to high (~1 GeV) energies. (6) Study of more aspects of microbunched structures in beam phase spaces For a cascaded amplifier in circuit electronics, the total amplification gain can be estimated as the product of individual gains. In a beam transport line of an accelerator, the (scalar) gain multiplication was examined and found to under-estimate the overall microbunching amplification. The concept of gain matrix was developed, which includes the density, energy and transverse-longitudinal modulations in a beam phase space, and used to analyze MBI for a proposed recirculating machine. Throughout the gain matrix approach, it reasonably gives the upper limit of spectral MBI gain curves. This extended analysis can be employed to study multi-pass recirculation. (7) Study of MBI for magnetized beams Driven by a recent energy-recovery-linac based cooler design for electron cooling at Jefferson Lab Electron-Ion Collider Project, the generalized theoretical formulation for MBI to a transversely coupled beam has been developed and applied to this study. A magnetized beam in general features non-zero canonical angular momentum, thus considered to be a transversely coupled beam. A novel idea of utilizing magnetized beam transport was proposed for improvement of cooling efficiency and possible mitigation of collective effects. A concern of MBI regarding this design was studied and excluded. The large transverse beam size associated with the beam magnetization is found to help suppress MBI via the transverse-longitudinal correlation.

43 PARTICLE ACCELERATORS↗