Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “feasible power flow solution”

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.

71 records · Page 4

Fabrication and Testing of Solid-Solution Strengthened Corrosion Resistant Alloys For Service in Molten Fluoride Environments

The demand for higher system thermal efficiencies requires the operation of power generation cycles and heat conversion systems at progressively higher temperatures. As the system operating temperature increases, existing materials may not provide adequate mechanical properties or environmental compatibility or both. There is an increasing commercial interest in the development and deployment of liquid-fueled Molten Salt Reactors (MSRs). Hastelloy®N, the highest performing candidate MSR structural alloy, is not capable of operations at temperatures above 700°C, thus limiting the performance of these systems. Using an Integrated Computational Materials Engineering (ICME)-approach and small laboratory scale heats, ORNL developed a class of patented alloys covered by U.S. Patent 9,435,011 B2, “Creep-resistant, Cobalt-free alloys for high temperature, liquid-salt heat exchanger systems,” similar to Hastelloy®N in that they are primarily solid solution strengthened. In contrast to precipitation strengthened alloys, the microstructure of solid solution alloys and hence the high temperature mechanical properties are stable for extended periods of time allowing long reactor operating life. The new alloys have shown to possess good resistance to liquid fluorides at temperatures up to 850°C and have significantly improved creep properties when compared to Hastelloy®N. The purpose of the CRADA project was for ORNL to collaborate with Haynes International- a materials producer, MetalTek International- a foundry, and Kairos Power – an advanced reactor developer – to scale-up selected alloys, evaluate their properties, and identify one solid solution strengthened alloy that can meet the property requirements for the reactor being developed by Kairos Power and other similar liquid fluoride-salt cooled reactors. As part of the project, eight alloys were down-selected and fabricated in larger industrial scale heats by Haynes International. Resistance to molten salt was evaluated in flowing FLiNaK and FLiBe by Kairos Power using their Rotating Cage Loop (RCL) system. Accounting for iron deposition during these tests, the new alloys displayed very low net mass change showing excellent corrosion performance in molten salt. Creep properties evaluated at ORNL were found to be better than that of Hastelloy®N and 316 stainless steel. Long-term stabilities of the alloys evaluated by Haynes International showed that these alloys have excellent thermal stability in the temperature range 704.4-815.6°C, with the change in strength and ductility being less than 10-15% after a 4000-hour exposure at 815.6°C. Autogenously Gas Tungsten Arc Welding (GTAW) welded samples showed less than 10% change in yield strength / ultimate tensile strength / total elongation compared to the basemetal, indicating that the alloys have excellent weldability. Three parts were successfully investment-cast using one alloy with very little voiding showing feasibility of fabricating parts using the casting process. This project enabled extensive interaction between the material producer Haynes International, casting supplier MetalTek, and reactor developer Kairos Power. This facilitated testing of materials and components produced using the newly developed alloys by the end-user. This allowed the generation of critical dataset required for down-selection of a few promising alloys for further development. This data is also currently being shared with other reactor designers for them to evaluate the suitability of this alloy for their reactor design. The availability of this alloy will ultimately enable the design and development and deployment of MSRs with increased temperature of operation and thus, improved efficiencies.

36 MATERIALS SCIENCE↗

Forced Mixer Nozzle Optimization

Computational fluid dynamic (CFD) and computational acoustic analyses (CAA) were performed for a TFE731-40 compound nozzle, a TFE731-60 mixer nozzle and an Energy Efficient Engine (E(sup 3)) mixer nozzle for comparison with available data. The CFD analyses were performed with a three dimensional, Navier-Stokes solution of the flowfield on an unstructured grid using the RAMPANT program. The CAA analyses were performed with the NASA Glenn MGB program using a structured grid. A successful aerodynamic solution for the TFE731-40 compound nozzle operating statically was obtained, simulating an engine operating on a test stand. Analysis of the CFD results of the TFE731-40 with the MGB program produced predicted sound power levels that agree quite well with the measured data front full-scale static engine tests. Comparison of the predicted sound pressure with the data show good agreement near the jet axis, but the noise levels are overpredicted at angles closer to the inlet. The predicted sound power level for the TFE731-60 did not agree as well with measured static engine data as the TFE731-40. Although a reduction in the predicted noise level due to the mixed flow was observed, the reduction was not as significant as the measured data. The analysis of the V2 mixer from the E(sup 3) study showed that peak temperatures predicted in the mixer exit flowfield were within 5 percent of the values measured by the exit probes. The noise predictions of the V2 mixer nozzle tended to be 3-5 dB higher in peak noise level than the measurements. In addition, the maximum frequency of the noise was also overpredicted. An analysis of the 3 candidate mixer nozzle configurations demonstrated the feasibility of using centerbody lobes and porosity to improve mixing efficiency. A final configuration was designed with a predicted thermal mixing efficiency that was 5 percent higher than the 3 candidate mixers. The results of the MGB noise calculations show that the final design will exceed the design goal of a 3 dB reduction in noise as compared to the baseline TFE731-40.

Sheoran, Yogi↗

Design of an Additively Manufactured Recuperator With 800 °C Inlet Temperature for sCO2 Power Cycle Application

Abstract Recuperators with design temperatures at and above 800 °C can further increase the thermal efficiency of supercritical CO2 power cycles by enabling higher turbine exhaust temperatures. Mar-M247 is a well-suited nickel-based superalloy for high temperature service due to its high creep strength that prevents excessive material thickness being required for pressure containment. Additive manufacturing using a high-speed laser-directed energy deposition (L-DED) process presents a promising solution, with build trials demonstrating the ability to produce nonconventional flow channels for enhanced heat transfer. A design process is presented that includes aerothermal and mechanical evaluation to maximize performance within the constraints of the manufacturing process. A 2-D heat transfer network and pressure drop code allows prediction of flow distribution and its effect on overall thermal performance. Established literature correlations, along with CFD simulation, inform the prediction of heat transfer coefficients and friction factors for the flowpaths and enhancement features in the heat exchanger core. Mechanical evaluation using finite element analysis (FEA) modeling with the intent of the ASME Boiler & Pressure Vessel Code (BPVC) Section VIII, Div. 2 assesses the operational safety of the design. The detailed design features annular finned passages that take advantage of helical flow paths to distribute the flow from separated headers to shared heat transfer surfaces. Performance predictions for the recuperator at a 50 kW scale provide insights into the feasibility of the additively manufactured (AM) process to produce recuperators on a commercial scale that extend existing operating envelopes.

Engineering↗

Automated Euler and Navier-Stokes Database Generation for a Glide-Back Booster

The past two decades have seen a sustained increase in the use of high fidelity Computational Fluid Dynamics (CFD) in basic research, aircraft design, and the analysis of post-design issues. As the fidelity of a CFD method increases, the number of cases that can be readily and affordably computed greatly diminishes. However, computer speeds now exceed 2 GHz, hundreds of processors are currently available and more affordable, and advances in parallel CFD algorithms scale more readily with large numbers of processors. All of these factors make it feasible to compute thousands of high fidelity cases. However, there still remains the overwhelming task of monitoring the solution process. This paper presents an approach to automate the CFD solution process. A new software tool, AeroDB, is used to compute thousands of Euler and Navier-Stokes solutions for a 2nd generation glide-back booster in one week. The solution process exploits a common job-submission grid environment, the NASA Information Power Grid (IPG), using 13 computers located at 4 different geographical sites. Process automation and web-based access to a MySql database greatly reduces the user workload, removing much of the tedium and tendency for user input errors. The AeroDB framework is shown. The user submits/deletes jobs, monitors AeroDB's progress, and retrieves data and plots via a web portal. Once a job is in the database, a job launcher uses an IPG resource broker to decide which computers are best suited to run the job. Job/code requirements, the number of CPUs free on a remote system, and queue lengths are some of the parameters the broker takes into account. The Globus software provides secure services for user authentication, remote shell execution, and secure file transfers over an open network. AeroDB automatically decides when a job is completed. Currently, the Cart3D unstructured flow solver is used for the Euler equations, and the Overflow structured overset flow solver is used for the Navier-Stokes equations. Other codes can be readily included into the AeroDB framework.

Chaderjian, Neal M.↗

Demonstrate FARM supervisory capabilities for a thermal energy storage problem for the DETAIL facility: IES Simulation Ecosystem Control System Development

The goal of the power dispatch problem for an Integrated Energy System (IES) is to adjust the power output and the heat flow of each component to maximize the profitability of the whole unit. Facilities that can integrate real-time digital signals, mock nuclear power, thermal energy storage and industrial heat use via high-temperature electrolysis were constructed at INL to support the research activities. The Dynamic Energy Technology and Integration Laboratory (DETAIL) houses the Microreactor Agile Non-nuclear Experimental Test Bed (MAGNET) and the Thermal Energy Distribution System (TEDS). In this report, the hierarchical control system architecture proposed in June 2023 milestone for the flexible operation of DETAIL facility is finalized and demonstrated. A brief description of the components and the corresponding Dymola models from the HYRBID repository is first provided. Then, the current control strategy is presented. In particular, the approach for generating the set-point trajectories to be fed to the PI controllers is analyzed, and its limits were identified. To preserve safe operation over both long-time and real-time horizons, the integration of a Supervisory Control layer embedding a modified version of FARM (Feasible Actuator Range Modifier) module is proposed. FARM is a component of the RAVEN-based FORCE framework designed to support HERON module at optimizing the operation of IES units. The proposed control system for DETAIL foresees FARM to be applied twice, i.e., the original version (“FARM-Validator”) aiding the solution of the power dispatch problem, and a modified version (“FARM-Supervisory”) coordinating the PID controllers. Despite the kernel of the two modules is the same, their tasks are quite different. The former intervenes at the beginning of each hour to prevent constraint violations over long time periods, the latter addresses real-time control tasks and monitors the response of constrained variables at a much finer time resolution. A tentative procedure for training the embedded Digital Twins with the experimental data is also proposed. Finally, the capabilities of the designed architecture and the impact of the added Supervisory Control layer are demonstrated by simulating a representative power dispatch scenario.

25 ENERGY STORAGE↗

Throughput Optimization of Molybdenum Carbide Nanoparticle Catalysts in a Continuous Flow Reactor Using Design of Experiments

Transition metal carbides (TMCs) have attracted significant attention because of their applications toward a wide range of catalytic transformations. However, the practicality of their synthesis is still limited because of the harsh conditions in which most TMCs are prepared. Recently, a solution-phase synthesis of phase-pure a-MoC1-x nanoparticles was presented. While this synthetic route yielded nanoparticles with exceptional catalytic performance, the reaction parameter space was not explored, and catalyst throughput was not optimized for scale-up. Continuous flow platforms coupled with statistical design of experiments (DoE) can provide a powerful method for understanding the reaction parameter space for optimizations. Here, we demonstrate the use of statistical DoE in tandem with response surface methodology for a parametric screening analysis to optimize the throughput of a MoC1-x nanoparticle synthesis utilizing a millifluidic flow reactor. A full factorial design was implemented to evaluate four input variables (reaction temperature, flow rate, solvent fraction of oleylamine, and precursor concentration) that carry statistically significant effects on three responses (throughput, residence time, and isolated yield). A Doehlert matrix was implemented to investigate each significant variable at a higher number of levels to optimize throughput. Our results give a nonintuitive set of experimental conditions that resulted in an optimized throughput of 2.2 g h-1. This translates to a 50-fold increase in throughput compared to the previously reported batch method. The catalytic performance of the MoC1-x nanoparticles produced under optimized throughput was demonstrated in the CO2 hydrogenation reaction. This DoE screening analysis and throughput optimization of MoC1-x synthesis open the door to an increased feasibility for scale-up.

design of experiments↗

Numerical investigation of internal high-speed viscous flows using a parabolic technique

A feasibility study has been conducted to assess the applicability of an existing parabolic analysis (ADD-Axisymmetric Diffuser Duct), developed previously for subsonic viscous internal flows, to mixed supersonic/subsonic flows with heat addition simulating a SCRAMJET combustor. A study was conducted with the ADD code modified to include additional convection effects in the normal momentum equation when supersonic expansion and compression waves are present. A set of test problems with weak shock and expansion waves have been analyzed with this modified ADD method and stable and accurate solutions were demonstrated provided the streamwise step size was maintained at levels larger than the boundary layer displacement thickness. Calculations made with further reductions in step size encountered departure solutions consistent with strong interaction theory. Calculations were also performed for a flow field with a flame front in which a specific heat release was imposed to simulate a SCRAMJET combustor. In this case the flame front generated relatively thick shear layers which aggravated the departure solution problem. Qualitatively correct results were obtained for these cases using a marching technique with the convective terms in the normal momentum equation suppressed. It is concluded from the present study that for the class of problems where strong viscous/inviscid interactions are present a global iteration procedure is required.

Anderson, O. L.↗

Solving the Pulsed-Power Driven Inertial Fusion Energy Standoff Problem (Abbreviated Final Report)

The standoff problem in pulsed-power driven IFE (Inertial Fusion Energy) refers to the problem of electrically connecting a fusion target with the driver in a way that preserves the driver/target interface in the presence of a large fusion energy release (100’s of megajoules to 1000 megajoules). High yield fusion drivers for stockpile stewardship applications have similar concerns, but without the complication of the high repetition rate needed for energy production. All proposed IFE or high yield systems have the challenge of protecting the reactor vessel and driver from the energy release, but pulsed power drivers have the additional challenge of establishing an electrical connection in a way that does not create excessive debris or require costly/massive structures that must be expendable. Fortunately, there are conceptual solutions in the form of very low mass transmission lines, usually in the form of wires or foils, since very little mass is needed to overcome the magnetic forces driving the transmission lines apart or provide a low resistance pathway given the short pulse duration of the driver. Several important aspects of the standoff problem were analyzed in the course of this LDRD project. Conceptual means of introducing the low mass transmission lines into the reactor chamber were proposed, including analysis of methods to retain the power flow gap between the electrodes in the presence of the chamber environment. The energy losses due to ohmic dissipation and magnetic acceleration of the low mass transmission lines were evaluated for candidate driver pulses and different transmission line masses, and the fluence of x-ray and neutron energy on the persistent power flow structures that must survive the pulse of energy were evaluated. Finally, means of rapidly pumping down the chamber were explored as a way to return the system to low pressure following the multi-atmosphere post-pulse pressure rise. The overall conclusion from the analysis is that, while the reactor environment creates stressful conditions for a low mass transmission line standoff system, there are concepts expending 10 kg or less of electrode mass each pulse with the potential to mitigate the stresses and successfully drive fusion targets. Whether such a system could lead to a practical and economical fusion reactor would require extensive research and development beyond the scope of the feasibility study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DEM and MELT In-Can Vitrification Process for Fukushima Daiichi Water Treatment Secondary Waste - 20034

Fukushima Daiichi Nuclear Power Station was severely damaged by a tsunami induced by the earthquake on March 11, 2011. A large amount of water treatment secondary waste has been generated from the purification work of contaminated water that consists of the water flowing into the reactor building. The Mid and Long-Term Road map specifies that solid waste at Fukushima Daiichi Nuclear Power Station should be stored safely and reasonably according to its characteristics. In order to meet this goal, a versatile treatment solution is needed, and the In-Can vitrification seems to be a suitable solution. The DEM and MELT In-Can vitrification process is being developed and qualified in France to treat waste arising from decommissioning and dismantling operations. This qualification is carried out through the consortium of the DEM and MELT project coordinated by CEA which gathers Orano, ECM technologies and Andra. The DEM and MELT In-Can vitrification process is a robust, simple and versatile in situ vitrification process. It is designed for high or intermediate level waste and is compact enough to be implemented in an existing facility or close to the waste to be treated. It is developed to treat liquid and solid waste, to produce a small amount of secondary waste and to minimize investment and operating costs. The DEM and MELT In-Can vitrification process is also developed with a modular design in order to be adapted to nuclear operators needs and requirements. The In-Can melting process has already been developed and fully qualified by the CEA for Alpha liquid waste and can be rated TRL 7 for this scope. The In-Can vitrification technology is considered as a qualified solution that could be applied for the treatment of waste resulting from the D and D work in Fukushima Daiichi Nuclear Power Station. The feasibility evaluation of the DEM and MELT In-Can vitrification process for Fukushima Daiichi water treatment secondary waste is a project led by Orano, CEA and ANADEC. This R and D project is performed by the funds from the Ministry of Economy, Trade and Industry as The Subsidy Program 'Project of Decommissioning and Contaminated Water Management'. This paper presents the first results of this feasibility evaluation, from the scenarios and glass formulation definition to the industrial implementation feasibility. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

HI-STORM Overpack and MPC-32 Thermal-Hydraulic Model with MOOSE Framework

Nuclear power is a significant source of electricity in the United States, but the average age of nuclear power plants is around 40 years old. The safe management of the spent nuclear fuel (SNF) is a key aspect of the back-end of the nuclear fuel cycle. Spent fuel dry storage systems are becoming a popular and effective solution in this area, given the absence of a final disposal system. The spent fuel cask system (dry cask method) provides a feasible solution to maintain spent fuel for 60 years before final disposal. Dry cask storage has many characteristics that make it attractive. It fulfills the safety requirements of the Nuclear Regulatory Commission (NRC) while providing modularity and flexibility to contractors. The HI-STORM overpack and MPC-32 canister are the main parts of the HI-STORM 100 dry cask storage system. These components remove heat from the system using natural circulation, requiring no human intervention. This is the characteristic that provides passive heat removal and low maintenance features in dry cask storage systems. To develop a thermal model for a dry cask storage system, the physics behind the system should be defined clearly. There are two natural circulation loops in the system; circulation of helium cools down the nuclear assemblies in the MPC, while circulation of air cools down the walls of the MPC. This work aims to develop a thermal model of the MPC-32 canister and HI-STORM overpack using the Multiphysics Object-Oriented Simulation Environment (MOOSE). MOOSE is an open-source framework developed by Idaho National Laboratory (INL) for multiscale, multiphysics simulations. In this study, we will investigate and demonstrate the thermal-hydraulics modeling capabilities of the MOOSE framework, including natural circulation, heat transfer, and porous flow.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of a MOOSE thermal model of the MPC-32 canister and HI-STORM overpack

Nuclear power is a significant source of electricity in the United States, but the average age of nuclear power plants is around 40 years old. The safe management of the spent nuclear fuel (SNF) is a key aspect of the back-end of the nuclear fuel cycle. Spent fuel dry storage systems are becoming a popular and effective solution in this area, given the absence of a final disposal system. The spent fuel cask system (dry cask method) provides a feasible solution to maintain spent fuel for 60 years before final disposal. Dry cask storage has many characteristics that make it attractive. It fulfills the safety requirements of the Nuclear Regulatory Commission (NRC) while providing modularity and flexibility to contractors. The HI-STORM overpack and MPC-32 canister are the main parts of the HI-STORM 100 dry cask storage system. These components remove heat from the system using natural circulation, requiring no human intervention. This is the characteristic that provides passive heat removal and low maintenance features in dry cask storage systems. To develop a thermal model for a dry cask storage system, the physics behind the system should be defined clearly. There are two natural circulation loops in the system; circulation of helium cools down the nuclear assemblies in the MPC, while circulation of air cools down the walls of the MPC. This work aims to develop a thermal model of the MPC-32 canister and HI-STORM overpack using the Multiphysics Object-Oriented Simulation Environment (MOOSE). MOOSE is an open-source framework developed by Idaho National Laboratory (INL) for multiscale, multiphysics simulations. In this study, we will investigate and demonstrate the thermal-hydraulics modeling capabilities of the MOOSE framework, including natural circulation, heat transfer, and porous flow.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Parallel Simulation of Unsteady Turbulent Flames

Time-accurate simulation of turbulent flames in high Reynolds number flows is a challenging task since both fluid dynamics and combustion must be modeled accurately. To numerically simulate this phenomenon, very large computer resources (both time and memory) are required. Although current vector supercomputers are capable of providing adequate resources for simulations of this nature, the high cost and their limited availability, makes practical use of such machines less than satisfactory. At the same time, the explicit time integration algorithms used in unsteady flow simulations often possess a very high degree of parallelism, making them very amenable to efficient implementation on large-scale parallel computers. Under these circumstances, distributed memory parallel computers offer an excellent near-term solution for greatly increased computational speed and memory, at a cost that may render the unsteady simulations of the type discussed above more feasible and affordable.This paper discusses the study of unsteady turbulent flames using a simulation algorithm that is capable of retaining high parallel efficiency on distributed memory parallel architectures. Numerical studies are carried out using large-eddy simulation (LES). In LES, the scales larger than the grid are computed using a time- and space-accurate scheme, while the unresolved small scales are modeled using eddy viscosity based subgrid models. This is acceptable for the moment/energy closure since the small scales primarily provide a dissipative mechanism for the energy transferred from the large scales. However, for combustion to occur, the species must first undergo mixing at the small scales and then come into molecular contact. Therefore, global models cannot be used. Recently, a new model for turbulent combustion was developed, in which the combustion is modeled, within the subgrid (small-scales) using a methodology that simulates the mixing and the molecular transport and the chemical kinetics within each LES grid cell. Finite-rate kinetics can be included without any closure and this approach actually provides a means to predict the turbulent rates and the turbulent flame speed. The subgrid combustion model requires resolution of the local time scales associated with small-scale mixing, molecular diffusion and chemical kinetics and, therefore, within each grid cell, a significant amount of computations must be carried out before the large-scale (LES resolved) effects are incorporated. Therefore, this approach is uniquely suited for parallel processing and has been implemented on various systems such as: Intel Paragon, IBM SP-2, Cray T3D and SGI Power Challenge (PC) using the system independent Message Passing Interface (MPI) compiler. In this paper, timing data on these machines is reported along with some characteristic results.

Menon, Suresh↗

A pushered capsule implosion as an alternate approach to the ignition regime for inertial confinement fusion

We report in inertial confinement fusion, the threshold for ignition is a highly dynamic quantity as the sources and sinks of power in the hot spot can vary rapidly. In this article, we consider the ignition condition as a race between heating and disassembly rates and make use of a prior solution to the fusion hot-spot thermodynamics to develop a Lawson-like ignition criteria for pressure × confinement time (p-τ) vs temperature. Low-Z capsule designs reach the temperature for this threshold using as much of the shell as feasible as ablator but then are limited in τ by low stagnated mass. An alternate approach, the pushered single shell (PSS) design [D. D.-M. Ho, S. MacLaren, and Y. Wang, “High-yield implosions via radiation trapping and high rho-R,” paper presented at the 60th Annual Meeting of the APS Division of Plasma Physics, 2018], introduces a dense inner layer of Mo-Be alloy that is smoothly graded outward to pure Be, increasing the confinement time at stagnation and lowering the temperature requirement at the ignition threshold. Here, we describe a PSS ignition design for the National Ignition Facility and use the theory as well as simulations to compare it with the low-Z capsule approach. Additionally, we show how an adjustment to the design is used to anticipate the effects of mixing at the fuel–ablator interface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Particle-based high-temperature thermochemical energy storage reactors

Solar and other renewable energy driven gas-solid thermochemical energy storage (TCES) technology is a promising solution for the next generation energy storage systems due to its high operating temperature, efficient energy conversion, ultra-long storage duration, and potential high energy density. Experimental and theoretical studies suggest that the respective gravimetric and volumetric TCES energy storage densities vary from 200 to 3000 kJ kg –1 and 1–3 GJ m –3 . Solar radiation or heat generated from electric furnaces powered by renewable electricity can be stored in the form of chemical energy through endothermic reactions, while the stored chemical energy can be converted to thermal energy via an exothermic reaction when needed. The design of highly effective reactors requires a deep understanding of materials, thermodynamics, chemical kinetics, and transport phenomena. At time of writing, TCES reactors are yet to be deployed at commercially relevant scales, leaving a substantial gap between development efforts and commercial feasibility. Therefore, this review aims to examine the state-of-the-art design and performance of particle-based TCES reactors with different reactive materials. Fundamentals related to TCES reactive materials, reaction conditions, thermodynamics and kinetics, and transport phenomena are reviewed in detail to provide a comprehensive understanding of the reactor design and operation. Five major types of TCES reactors have been comprehensively reviewed and compared, including fixed, moving, rotary, fluidized, and entrained bed reactors. Most reported prototype reactors in the literature operate at lab scale with thermal inputs below 40 kW, and scaled TCES reactors (e.g., at megawatt level) are yet to be demonstrated. The nominal reactor operating temperatures range from 300 to 1500 °C, depending on the selected chemistry, reactive material, and heat sources. To evaluate their designs, the reactors are assessed in aspects of performance, cost, and durability. Discrepancies in performance indicators of energy storage density, extent of reaction, and various energy efficiencies are highlighted. The scale-up of reactors and power block integration, which hold the key to the successful commercialization of TCES systems, are critically analyzed. Furthermore, advanced materials (both reactive materials and ceramic reactor housing materials), effective particle flow control, advanced modeling tools, and novel system design may bring significant improvement to the energy efficiency, storage density and cost competitiveness of particle-based TCES reactors.

25 ENERGY STORAGE↗

AGGREGATE: dAta-driven modelinG preservinG contRollable dEr for outaGe mAnagemenT and rEsiliency (Final Report)

The AGGREGATE project team successfully developed and validated various modules for outage management. Brief summaries of each module are provided to showcase their strength for outage management and restoration for a distribution system with a high penetration of connected distribution energy resources (DERs). In recent years, inverter-based DERs have been widely deployed in distribution system. A most of behind-the-meter (BTM) solar power generation is not visible to the utility. The data-driven DER and load estimation modules are using machine learning (ML) and artificial intelligence (AI) to manage this issue, which provides an opportunity for distribution system operators (DSOs) to operate systems and make decisions in real-time for a distribution system with a high penetration of DERs deployed. Also, the estimated DER and true load can be further leveraged in network aggregation and cold-load pick up estimation for reducing the computing complexity and providing for fast restoration. After load demand and DER power generations have been estimated, the information will support topology and state estimation (SE). The topology estimation module demonstrated the viability of mixed integer linear programming (MILP) formulation to estimate the most likely operational radial topology and outage sections using power flow measurements, historical/estimated load and DERs data and smart meter ping measurements. Formulation includes continuous (power flow, load and DERs data) and binary measurements (smart meter ping measurements) in a single formulation. Errors in continuous data and binary data are modeled as normal distribution and Bernoulli distribution, respectively. In the future distribution grid, the power injection from controllable DERs will be essential for efficient and resilient grid operation. However, determining the optimal DER injections and restoration actions is dependent on knowledge of the system states. State estimation (SE), already the cornerstone of transmission energy management systems, will become commonplace in distribution management systems as more measurements become available from deployment of automated metering infrastructure (AMI). Observability analysis is the first step in SE, as it determines the sufficiency of the available measurements for accurately estimating the current system states. A new type of pseudo-measurement called a Correlational Measurement (CM) is introduced in this module, to enhance the observability of the system to enable more accurate SE. CMs encapsulate knowledge of correlation between demand patterns for similar classes of loads as well as injection patterns for same-technology renewable DERs. During grid contingency scenarios, DERs have been traditionally disconnected, without any fault ride-through capabilities. However, with new regulations and better technology, it is feasible for these resources to contribute to the grid’s restoration after an adverse event and hence enhance resilience. The controllability module proposes a two-step restoration scheme for the power system restoration process by leveraging additional degrees of freedom in power electronics interfaced DERs for mitigating voltage problems. In a resilience mode without the utility system, the distribution grid relies on DERs to serve critical load. In such a severe event with multiple faults on the distribution feeders, actuation of various protective devices (PDs) divides the distribution system into electrical islands. The undetected actuated PDs due to fault current contributions from DERs can delay the restoration process, thereby reducing the system resilience. The Advanced Outage Management (AOM) and the Advanced Feeder Restoration (AFR) modules developed in this project provide improved system resilience with multiple DERs. AOM identifies the faulted sections and actuated PDs in a distribution system with DERs by incorporating smart meter data. The most credible outage scenario including fault locations, PD actuations, and fault indicator (FI) failures is identified by a set of binary integer linear programming incorporating hypotheses. The AFR module serves to restore a distribution system with available energy resources taking into consideration the availability of utility sources and DERs. By partitioning the system into islands, critical load will be served with the available generation resources within islands based on the solution of a MILP. When the utility systems become available, the optimal path will be determined by a spanning tree search algorithm that reconnects these islands back to substations and restores the remaining load. The transmission and distribution (T&D) co-simulation module was used to validate the effect of a control action performed on the distribution side assets as it propagates to the transmission side. This ensures that the control action performed results in a feasible operating point on both the transmission and the distribution system. In addition to validation, the team used the T&D co-simulation module to demonstrate how distribution system assets can be used to mitigate issues on the transmission system. Specifically, the team demonstrated that appropriate switching operations on the distribution side can alleviate the line overload condition on the transmission side without causing new operational constraint violations.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Community-Centric Microgrid Feasibility Analysis Framework to Accelerate Resilience of Island Communities

Marine energy offers a reliable energy solution for island and coastal communities, which often lack traditional local generation, to support their transition to energy independence and reduce reliance on externally imported fuels. Successful deployment of new technologies in these isolated locations requires community acceptance and approval from the outset, as these communities typically lack the financial and technical resources to operate and maintain new systems. This report presents a community-centric microgrid planning framework for remote coastal and island communities. Community engagement is integrated as the first step in the planning process, incorporating community profiles and visions into energy development scenarios. A case study was conducted in St. George, Pribilof Islands, Alaska, which relies entirely on diesel yet has significant wind and wave energy potential. Community engagement revealed a unique history and current economic status, with an interest in adopting advanced energy technologies despite past failures. Various microgrid configurations were optimized, considering different technologies to meet current and future energy needs while balancing cost and energy resilience. Wave energy converters (WECs) were a key component, integrated with other energy sources using the Xendee optimization tool. The Marine Energy Microgrid Toolkit, developed as part of this work, uses commercial power system analysis tools to optimize and analyze microgrid scenarios. The developed framework and toolkit can be applied to island and coastal communities to enhance resilience and support microgrid deployments. Future enhancements will include incorporating new marine resources, developing dynamic models, and automating the integration of Xendee and PowerFactory simulations.

02 PETROLEUM↗

Investigating Low-Altitude Constellations of Ad-Hoc Lunar PNT System for Distributed Spacecraft Autonomy

In this study, we examine a low-altitude Lunar Position, Navigation, and Timing (LPNT) constellations and the localization performance of Centralized Extended Kalman Filter (CEKF) and Decentralized Extended Kalman Filter (DEKF) algorithms. The primary investigation involves a 100-node swarm operating at a 100 km altitude, in contrast to previous studies that examined a 21-node asset in a frozen-orbit at 5,500 km. The autonomous operation of large-scale swarm is based on two-way Inter-Satellite Link (ISL) measurements, which involve pseudoranges and relative velocities among swarm nodes. We perform a numerical assessment of the two filtering approaches, utilizing ‘fully sampled’ measurements from all available assets as well as ‘two ISL’ measurements where each spacecraft is restricted to only two antennas. This research includes an analysis of CEKF under 2-ISL constraints and evaluates the performance of DEKF in a 100-node swarm, which has not been explored in previous studies. In addition, we examine the impact of increasing the sampling frequency for DEKF, showing that the update cycle can be shortened from a 10-minute interval. A novel approach for ‘2-ISL limited’ DEKF will also be introduced, using a matching formulation that exhaustively enumerates all potential matches. This study provides valuable insights into large-scale distributed swarm operations, considering various filter configurations, sampling frequencies, matching strategies, and scalability of CEKF and DEKF for low-altitude LPNT applications. The Lunar PNT technology plays a key role in providing reliable and robust navigation services on the Moon's surface and the South pole, where the primary Lunar missions are planned. To support upcoming Lunar missions, including small satellites from NASA's Commercial Lunar Payload Services program, the Lunar PNT system must be adaptable to smaller platforms like CubeSats. Driven by the growing involvement of public and private exploration partnerships, the traditional low Earth orbit missions are shifting to beyond geosynchronous orbit [1]. These upcoming missions aim to foster a sustainable and innovative exploration program, in collaboration with commercial and international partners, to facilitate human expansion throughout the solar system and return new knowledge and opportunities to Earth [2]. As part of this trend, there are increasing efforts to utilize science missions in Lunar orbit to develop a non-dedicated and ad-hoc PNT network system. Two traditional approaches, the Deep Space Network (DSN) and the weak signal Global Positioning System (GPS), are established deep-space navigation technologies for missions beyond the geosynchronous orbit. Beginning in 1958, the DSN was developed to communicate with the Explorer 1 spacecraft based on the use of radiometric tracking in spacecraft navigation [3]. The DSN is capable of providing nearly unfettered coverage to spacecraft beyond low-Earth orbit (LEO), however, increased space mission volume has created concerns about future expectations of DSN usage for spacecraft navigation [4]. For cislunar mission applications, the position accuracy using DSN achieves 100 m (3σ) with at least three geometrically diverse ground stations when using radiometric tracking alone [5]. The DSN's dependence on Earth-based ground stations restricts its operational capabilities to periods of Earth visibility. This limitation, coupled with its poor localization performance, renders the DSN unsuitable for future lunar missions that demand continuous tracking and precise positioning. To satisfy the increasing requirements of DSN in Lunar applications, spacecrafts are also required to improve their onboard antenna power and efficiency of the transmission. However, there is an important aggregate cost trade between adding capabilities to every spacecraft and adding to a capacity on the ground that serves multiple spacecraft [6]. A weak GPS system can provide PNT service while the user spacecraft is bound to the Moon, leveraging a single, steerable high gain antenna with the relatively narrow beam which includes all the sources in its field of view [7]. However, the higher the altitude the receiver is above the GPS constellations, the poorer and the weaker are the relative geometry and the received signal powers, respectively, leading to a significant navigation accuracy reduction [8]. The transmitted power becomes weaker with increasing distance from the Earth as well as signals tracked from one of the side lobes of the GPS antenna pattern. As a results, the number of visible satellites and relative geometric condition of the GPS satellites at very high altitude drops dramatically and reduces the navigation solution accuracy. Therefore, the weak GPS system is also not an ideal way to provide PNT service to upcoming Lunar missions when considering its limited geometric condition and the recued navigation accuracy. Another navigation approach on the Moon is being developed, similar to the Global Navigation Satellite System (GNSS) on Earth, aiming to offer navigation service with continuous 24/7 coverage across the entire Lunar surface. For example, lunar communications relay and navigation systems (LCRNS) by NASA and Lunar navigation satellite systems (LNSS) by JAXA are designed to serve as dedicated Position, Navigation, and Timing (PNT) systems for the Moon. However, designing a dedicated LNSS and PNT service involves additional challenges, which are unique to the lunar environment, including limited payload capacity for the CubeSat platform, i.e., the size, weight, and power (SWaP) of the onboard clock, limited lunar ground monitoring stations, and limited financial investment as compared to the legacy Earth-GPS [9]. NASA’s focus on utilizing CubeSat platforms on the Moon leads to an alternative Lunar navigation platform that leverages the existing Lunar science and exploration assets. The small satellites used in Lunar missions can be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions [10]. As upcoming Lunar missions will often operate at low-altitude about 30 km to 100 km for scientific observations and mapping purposes, the low-altitude orbital constellations could be employed to create an ad-hoc Lunar PNT system. However, several issues must be addressed, such as the instability of these orbits, which often require maintenance or are only suitable for short-duration missions, operating for fewer than 90 days. Additionally, at an altitude of 100 km, the satellites have a limited period during which they are above the horizon and capable of providing PNT service to users. The implementation of a non-dedicated, ad-hoc Lunar navigation constellation facilitates on-demand PNT services. A preliminary study of ad-hoc Lunar PNT system was conducted using 21 spacecraft in 5,5000 km altitude frozen orbits to test its feasibility and a basic performance of orbital asset localization among ad-hoc Lunar constellations in small satellites format [10]. These swarm assets are designed for autonomous localization with minimal Earth interaction, reducing dependency on bandwidth and ground resources. The design in [10] demonstrated the feasibility of a decentralized PNT approach, specifically employing a DEKF approach for state estimation, which helps minimize onboard operating costs. The DEKF method distributes computation across individual satellites, which lightens the computational load while maintaining accuracy in orbit ephemeris and clock offsets, similar to centralized systems [11]. In a follow-on study [12], each spacecraft was limited to 2 communications antennae, forcing the selection of measurements and scheduling spacecraft activities to perform the measurements. A matching algorithm is implemented to select the best measurements and schedule position estimation updates. The decentralized localization performance is also investigated with increasing levels of network degradation for swarm assets considering the impact of intermittent and permanent communication failure, to demonstrate the robustness and fidelity of the decentralized Lunar PNT service [13]. This study confirmed that the ad-hoc PNT constellations in frozen orbit are highly robust and resilient to communication failures. However, unlike frozen orbit swarm assets, the low-altitude satellites have a limited ground view at an altitude of 100 km, where the ad-hoc Lunar constellation consists of 98 low-altitude satellites, evenly distributed across seven circular polar orbital planes, alongside two satellites in a frozen orbit at an altitude of 5,500 km (Figure 1). Therefore, the number of satellites visible to ground users is significantly limited in low-altitude orbit constellations. As each visibility of a spacecraft remains intact for only a few ticks before it moves out of the field of view, the ground user encounters challenges in maintaining continuous navigation service, resulting in sparse availability and provision of Lunar PNT system. Consequently, service availability is primarily restricted to the Lunar South Pole region (Figure 2). Given these limitations and concerns, the localization performance of low-altitude swarm assets will be assessed in this study. We focus on the investigation of the localization performance of low-altitude swarm assets and ground users near the Lunar South Pole. The overall flow of the Lunar PNT simulation incorporates the DEKF approach of asset localization and the weighted least-squares approach in user localization (Figure 3). The autonomous Lunar PNT simulation is primarily implemented in MATLAB, where the DEKF based on the matching scheduler is implemented with Google’s OR-tools as a model builder and Gurobi optimization tool as a backend solver. The General Mission Analysis Tool (GMAT) is utilized to generate ephemeris data for swarm assets, and accounts for satellite orbital details, mass, and perturbations like solar radiation pressure and drag coefficients. Each ephemeris dataset is produced in the Moon International Celestial Reference Frame (ICRF) inertial coordinate system. For state estimation, the distributed swarm assets rely on two-way Inter-Satellite Link (ISL) measurements, which involve tracking pseudoranges and relative velocities between visible satellites and anchor nodes during each observation. Numerical evaluations of the decentralized localization process are conducted to demonstrate the feasibility of the low-altitude PNT system in providing reliable navigation services. The main approach involves using DEKF and CEKF to localize 100 satellites in low-altitude constellations, where the CEKF is implemented to serve as a baseline for comparing the performance of distributed algorithms. In both cases, we evaluate ‘fully sampled’ measurements from all available assets, and ‘two ISL’ measurements when spacecraft are constrained to have only two antennas. We test four estimation techniques: CEKF fully sampled, CEKF two ISL, DEKF fully sampled, and DEKF two ISL filters. As the DEKF update cycle is comprised of network setup, communication, and computations, a global broadcast network and 2-way ISL network setup will take from 4 to 6 minutes as maximum [12]. In this simulation, the DEKF update cycle is set to 10 minutes, including a 4-minute latency for obtaining and computing the actual measurement updates. We experiment an increased update cycle to demonstrate the feasibility and evaluate the impact on localization performance using various tuning values for measurement noise covariances (Figures 4 and 5). By comparing centralized and decentralized approaches using a matching algorithm, we analyze the influence of cross-correlation factors in the covariance matrix, assuming 100% reliability of all assets and measurements. The increased frequency and the adjustments of tuning parameters reveal distinct error patterns between the two scenarios. The localization accuracy of the swarm assets and ground users is assessed by taking the median error across 100 assets and one ground user (84.9°S, 137.5°E) over 7-day simulation period (Table 1). Since the user localization accuracy is significantly affected by the performance of the swarm assets, it is crucial to maintain high localization accuracy within the swarm. This study will continue to explore decentralized filtering for autonomous LPNT operations, with further investigation of an 'iterative' matching approach which enumerates every valid matching pair, planned for the following month.

Yeji Kim↗