Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuation power flow”

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

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

At least 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High-fidelity multiphysics load following and accidental transient modeling of microreactors using NEAMS tools: Application of NEAMS codes to perform multiphysics modeling analyses of micro-reactor concepts

The feasibility of modeling microreactors using high-fidelity models with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) tools is investigated in this report. Three overarching questions guided this research: can NEAMS tools readily be applied for high-fidelity multiphysics modeling of different types of transients in microreactor designs; how accurate are the results obtained; and are improvements needed in accuracy or user experience of NEAMS tools, especially considering newly developed capabilities? This work builds upon FY-2022 work, and two microreactor concepts considering heat pipe (HP-MR) and gas-cooled (GC-MR) technologies were further analyzed using high-fidelity multiphysics simulations. The NEAMS tools considered and coupled within the MultiApp environment are Griffin for neutronics, BISON for thermo-mechanics, Sockeye for heat pipe modeling (in HP-MR), SAM for 1D Fluid – 3D solid modeling of coolant channels and system modeling of balance of plant components (in GC-MR), and the SWIFT code for hydrogen redistribution in hydride moderator. The Heat Pipe MicroReactor (HP-MR) concept was further analyzed in FY-2023 to demonstrate the stochastic TRISO failure modeling capability in BISON to check operational limits of the TRISO fuel. A new full-core Gas-Cooled MicroReactor (GC-MR) model was developed based on the initial assembly-model used in Y-2022 and used for steady-state and accidental depressurization transient simulations. Accuracy of the simulations performed was assessed through 1) verification analyses completed on the different physics with code-to-code comparison, and 2) validation of the multiphysics simulations based on modeling of the Kilopower Reactor Using Stirling Technology (KRUSTY) experiment. In FY-2023, the mesh and model of KRUSTY was updated to closely match publicly available data, and the neutronic model was verified and validated against experimental control rod worth measurements. The multiphysics model of KRUSTY was developed and used for steady-state analysis and for modeling reactivity insertion transient. The calculated power increase and stabilization agrees well with experimental data following adjustment in fuel thermal expansion coefficient. As an important component of this project, the ANL team gathered experience with a wide range of NEAMS tools: the MOOSE Mesh System, Griffin, BISON, SWIFT, Sockeye, SAM, Workbench, and the MOOSE MultiApp System, and provided assessment of new capabilities. Noteworthy are the user assessment of the “vapor-only” flow model in Sockeye and development of a multiphysics startup transient in HP-MR unit cell for use as tutorial in Sockeye. The full-core GC-MR model was used for assessment of SAM for balance of plant modeling and for demonstrating the SWIFT code capability for hydrogen redistribution modeling in multiphysics transient analyses. In this process, several bugs/issues were identified and reported to developers. Finally, the assembly GC-MR model developed in FY-2022 coupling Griffin, BISON and SAM through flow blockage and rod ejection transients was published to the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB). The Heat Pipe MicroReactor (HP-MR) concept high-fidelity multiphysics coupling of Griffin/BISON/Sockeye in load-following and heat pipe failure transients was also published on the VTB. Those submissions are enabling thorough review of these models as well as wide distribution to industry, regulator, and university users. In this analysis, several new research questions were uncovered, and follow-up analyses are recommended to further improve some models, consider additional transients, and continue development of VTB models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Studying Wind Loading on CSP Collectors to Improve Performance and Reliability

Electricity generation through Concentrating Solar-thermal Power (CSP) adds the advantage of thermal energy storage and heat production. CSP collector costs constitute nearly one-third of total plant costs. Managing wind loading, especially dynamic wind loading caused by the turbulent wind flow, is a significant design challenge. Traditional collector designs rely on wind tunnel experiments and numerical simulations, which do not fully capture dynamic effects. Here, we provide an overview of NREL's activities on wind loading on CSP collectors, aiming to improve their reliability and cost-efficiency. We present insights from field measurements at operational parabolic trough and power-tower CSP plants. Over two years, we measured turbulent wind fields and resulting structural loads in a parabolic trough plant, revealing wind, turbulence, and collector field interactions. Similar measurements are ongoing at the Crescent Dunes heliostat field. These campaigns show how atmospheric turbulence, wind direction, and collector orientation affect dynamic wind loading. Upwind collector structures can generate turbulent structures, causing fluctuating loads at downstream collectors, impacting fatigue lifetime and optical efficiency. Preliminary results indicate that interior collectors experience higher turning moments than those at the edges, suggesting higher drive wear. At Crescent Dunes, we use instrumented heliostats at the field's edge and interior to study the translation of turbulent wind to dynamic structural loads. The campaign is still ongoing, and we will present first findings. Another focus of our work is tying the experimental wind loading findings to optical performance of the collectors, using simulation tools. For example, we derive torsional errors of the parabolic troughs - the angular offset from the sun position along a collector row. This error is impacted by wind-induced mirror deformations. Future work will focus on continuing the measurement efforts and combining them with computational fluid dynamic simulations. By providing detailed field measurements and validated models, our efforts aim to enhance understanding of wind loading on CSP collectors, improving their structural integrity and optical performance.

CSP↗
NASA NTRS2004

NASA Multidimensional Stirling Convertor Code Developed

A high-efficiency Stirling Radioisotope Generator (SRG) for use on potential NASA Space Science missions is being developed by the Department of Energy, Lockheed Martin, Stirling Technology Company, and the NASA Glenn Research Center. These missions may include providing spacecraft onboard electric power for deep space missions or power for unmanned Mars rovers. Glenn is also developing advanced technology for Stirling convertors, aimed at substantially improving the specific power and efficiency of the convertor and the overall power system. Performance and mass improvement goals have been established for second- and third-generation Stirling radioisotope power systems. Multiple efforts are underway to achieve these goals, both in house at Glenn and under various grants and contracts. These efforts include the development of a multidimensional Stirling computational fluid dynamics code, high-temperature materials, advanced controllers, an end-to-end system dynamics model, low-vibration techniques, advanced regenerators, and a lightweight convertor. Under a NASA grant, Cleveland State University (CSU) and its subcontractors, the University of Minnesota (UMN) and Gedeon Associates, have developed a twodimensional computer simulation of a CSUmod Stirling convertor. The CFD-ACE commercial software developed by CFD Research Corp. of Huntsville, Alabama, is being used. The CSUmod is a scaled version of the Stirling Technology Demonstrator Convertor (TDC), which was designed and fabricated by the Stirling Technology Company and is being tested by NASA. The schematic illustrates the structure of this model. Modeled are the fluid-flow and heat-transfer phenomena that occur in the expansion space, the heater, the regenerator, the cooler, the compression space, the surrounding walls, and the moving piston and displacer. In addition, the overall heat transfer, the indicated power, and the efficiency can be calculated. The CSUmod model is being converted to a two-dimensional model of the TDC at NASA Glenn. Validation of the multidimensional Stirling code is an important part of the grant effort. UMN has been generating data in an oscillating-flow test facility using two different test sections: a 90 turn and a cooler/regenerator/heater test section. CSU has created computational fluid dynamics models of both these test sections and has been making comparisons with the data, then improving their models to improve the agreement with the test data. CSU has also been using data available in the literature for code validation. UMN is now preparing to begin fabrication of a new 180 turn test section that will be more representative of certain portions of the Stirling engine geometry. Simulations to almost periodic steady state with the two-dimensional CSUmod model indicate that, to reach periodic steady state on a single 2-GHz desktop computer, 75 to 100 complete simulation cycles would be required and between 1 and 2 months of computer time. Therefore, Glenn has purchased the first 8 computers, of a 64-computer cluster, to be run in parallel to accelerate the simulation. On the basis of CFD Research Corp.'s experience with running the parallelized version of CFD-ACE on their clusters, we estimate that the complete 64-computer cluster will reduce simulation computing time by a factor of about 40. Plans are to continue development of these multidimensional Stirling codes and to use them to study the fluid-flow and heat-transfer phenomena that occur inside Stirling convertors. This is expected to lead to improved thermodynamic loss understanding, onedimensional design and performance codes, and engine performance.

Tew, Roy C.↗
NASA NTRS2011

An Examination of the Effect of Boundary Layer Ingestion on Turboelectric Distributed Propulsion Systems

A Turboelectric Distributed Propulsion (TeDP) system differs from other propulsion systems by the use of electrical power to transmit power from the turbine to the fan. Electrical power can be efficiently transmitted over longer distances and with complex topologies. Also the use of power inverters allows the generator and motors speeds to be independent of one another. This decoupling allows the aircraft designer to place the core engines and the fans in locations most advantageous for each. The result can be very different installation environments for the different devices. Thus the installation effects on this system can be quite different than conventional turbofans where the fan and core both see the same installed environments. This paper examines a propulsion system consisting of two superconducting generators, each driven by a turboshaft engine located so that their inlets ingest freestream air, superconducting electrical transmission lines, and an array of superconducting motor driven fan positioned across the upper/rear fuselage area of a hybrid wing body aircraft in a continuous nacelle that ingests all of the upper fuselage boundary layer. The effect of ingesting the boundary layer on the design of the system with a range of design pressure ratios is examined. Also the impact of ingesting the boundary layer on off-design performance is examined. The results show that when examining different design fan pressure ratios it is important to recalculate of the boundary layer mass-average Pt and MN up the height for each inlet height during convergence of the design point for each fan design pressure ratio examined. Correct estimation of off-design performance is dependent on the height of the column of air measured from the aircraft surface immediately prior to any external diffusion that will flow through the fan propulsors. The mass-averaged Pt and MN calculated for this column of air determine the Pt and MN seen by the propulsor inlet. Since the height of this column will change as the amount of air passing through the fans change as the propulsion system is throttled, and since the mass-average Pt and MN varies by height, this capture height must be recalculated as the airflow through the propulsor is varied as the off-design performance point is converged.

Felder, James L.↗
DOE OSTI2022

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

ExaSGD: 2022 Kernel Thrust Activities

The Kernel Thrust milestone ADSE22-407 covers the development of device-capable optimization algorithms and solvers technologies required by the ExaSGD project’s software stack in order to solve security-constrained alternating current optimal power flow (SC-ACOPF) problems on emerging exascale architectures. To this extent, in FY22 the main objective of the Kernel Thrust was (i) provide sparse optimization solver that runs efficiently on hardware accelerator devices (i.e., NVIDIA and AMD GPUs) to perform intra-node computations, (ii) strengthen the reliability and increase the performance of the mixed-dense sparse (MDS) solver of HiOp for deployment on the FY22 target architectures, Summit and Crusher, and (iii) increase performance by improving the mathematical algorithm and refining the parallel MPI-based implementation of the coarse-grain parallel solver HiOp-PriDec for capabilities deployment on the FY22 target architectures, Summit and Crusher. This document presents the developments and contributions done by the Kernels Thrust Team in FY22 toward completion of the above-mentioned objectives. These contributions progressed along four main development (sub)thrusts: (1) Design and implementation of a sparse optimization solver for use on hardware accelerators; (2) Improvement of the mathematical algorithm and of the parallel implementation of HiOp-PriDec to ensure readiness and efficient coarse-grain parallelism for FY23 target exascale machine; and (3) Support Software and Application Development Thrusts of the exaSGD project in their deployment of the project’s software stack on AMD- and NVIDIA-based architectures. The development of the sparse optimization solver (thrust 1 above) was new in FY22 and resulted in a new sparse solver in HiOp (available as of version 0.6). The second development thrust was a continuation of the efforts from FY21 and improved the mathematical algorithm and the communication strategy of the HiOp-PriDec solver. The last developement thrust is a large collaborative effort. Namely, the project’s teams from multiple labs (LLNL, PNNL, ORNL, and NREL) performed large-scale demonstration of the ExaSGD software stack, namely the optimization solvers of HiOp interfaced with the modeling front-end ExaGO and the stochastic sampler PowerScenarios. These demonstration efforts solved large-scale instances of the SC-ACOPF challenge problem of medium network sizes (10, 000-bus system) and large number of contingencies on Summit (NVIDIA accelerators) and Crusher (AMD accelerators) systems at ORNL.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

Towards Sheaf Theoretic Analyses for Delay Tolerant Networking

The goal of Delay Tolerant Networking (DTN) is to take a collection of heterogeneous, disparate connections between satellites, space assets, ground stations, and ground infrastructure and bring it together into a cohesive, functioning overlay network. Depending on the systems being considered, one can find links with a one-way light time exceeding minutes (and hours),periodic links which can sometimes be predicted by orbital mechanics, and restrictions based on the variety of capabilities built into these systems. These characteristics preclude traditional network models and routing techniques and have classically led to either rigid routing tables or purely probabilistic models. As the deeper underlying structures remain unknown, development of more DTN-optimized algorithms has lacked the necessary foundation. In a continuation of previous work, the goal of this paper is to identify and study these fundamental structures that exist in delay tolerant networks (DTN), with a focus on space networks. The current routing methodology has been to use contact graph routing (CGR) algorithms. CGR models a series of known contacts as a static graph. For CGR to work, this graph must be globally consistent and must have an accurate picture of the network. Because this is a globally controlled structure, there is little room for flexibility in the event of changes to the network which would naturally occur as the network grows. As a response to the desire for flexibility as the network changes, we introduced the mathematical structure known as sheaves to DTNs last year. The tag-line for sheaves is that they are a mathematically precise way of gluing local data together into unique global data. Thus, sheaves lend extra power to traditional models(and routing algorithms) by taking additional information and merging it, in as consistent a manner as possible, with the representation itself. The clearest example of how Earth-bound networks exhibit behavior that is “sheafy” is link state routers, which build a local-to-global picture of their network by gluing local information together into a global network, exactly as a sheaf would do. For routing within delay tolerant networks to truly exploit this structure, a deeper structure than a graph is required. In this paper, we develop sheaves that can work over directed graphs such as temporal flow networks, we construct a sheaf representation for Dijkstra’s algorithm, and we outline a construction for routing sheaves capable of modeling multicast scenarios. Finally, there is a section of future work suggesting follow-on research.

Robert Short↗
DOE OSTI2021

Analysis of Lobe Power Calculator and Indication System with Physics and Cycle Based Models

The Advanced Test Reactor (ATR) at INL measures reactor power through two methods, thermal and Nitrogen-16 (N-16) activity. Water power calculator (WPC) is a thermal power system that measures flow and temperature to determine the thermal quadrant powers. The N-16 system utilizes a beta chamber detector that outputs reactivity levels to calculate lobe power through an algorithm called lobe power calculation and indication system (LPCIS). The LPCIS utilizes the N-16 system and the WPC system to determine reactor core power levels. The WPC provides accurate calculations of quadrant and total reactor thermal power. With the use of WPC measurements, thermal-to-N-16 (T2N) power ratios are produced to determine if the two indication systems agree on core power. Relative magnitude equations are used to utilize N-16 coefficients and multipliers to improve the indications of the LPCIS. This is crucial for maintaining safety limits. Currently, the LPCIS system uses linear equations and matrices to calculate lobe power through multipliers and coefficients. Advancements in technology and system upgrades have assisted system engineers at INL with the objective to reach a more dynamic system. In return the system demonstrates an increase in the accuracy of power reading while maintaining safety margins. The new proposed coefficient and multiple method implements a cycle specific, and a physics-based model intended for changing coefficients during operation. This calibration experiment focuses on power splitting, outer shim and neck shim, as well as fuel burning into the RDAS weighting factors. Results show that the Physics learning method yields a smaller error margin inside of the desired power range for the data set 166-A. This proved true for both constrained and unconstrained testing. This is most likely due to the physics data fitting approach that resulted in favorable coefficients and multipliers for 166-A. Continuing to improve the physics-based model will help improve the power accuracy of the LPCIS system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Polymer Macrocycles: A novel topology to control dynamics of rubbery materials (Final Report)

One of the current challenges in polymer physics and engineering is to understand the behavior of circular macromolecules. Though such molecules have been made by dilute solution ring closure methods in the past, such procedures have been limited to relatively small, unentangled rings as well as to small sample quantities. Initial studies were also hampered by issues surrounding purity of the rings, though this problem has been somewhat overcome by the use of liquid chromatography at critical conditions (LCCC) methods that separate rings and linear chains. Furthermore, in spite of the continued progress a full understanding of ring dynamics remains elusive. The present work has addressed multiple aspects of the behavior of ring molecules by using a novel synthetic route: reversible radical recombination redox polymerization (R3P) to produce poly(3,6-dioxa-1,8-octanedithiol) (PolyDODT) molecules. The synthesis conditions were chosen to produce linear (LDODT), cyclic (RDODT) and linear-cyclic mixtures (LRDODT) in order to investigate dilute solution and rheological responses of the materials. Importantly, the R3P synthesis makes larger quantities and higher molecular weights than do the dilute solution ring closure methods and also has the potential to be scaled to industrial quantities, thus becoming commercializable. We emphasize two important outcomes of the work. The first is that the rheological investigation shows that very large molecular weight the PolyDODT rings permitted the investigation of ring dynamics from glassy towards terminal flow for molecular sizes corresponding to approximately 300 entanglements of the linear counterpart. At the same time, these very large sized rings could be diluted in a solvent in order to reduce the molecular entanglement density to less than unity. The results demonstrate that ring viscosity follows a Rouse-like linear dependence on molecular weight or entanglement number to approximately 15 entanglements. This range of "Rouse-like" response is greater than what was achieved in all prior works where synthesis conditions have limited the molecular sizes achievable. After the entanglement number $Z_w$ increases beyond approximately 15, we find that the viscosity-molecular weight scaling undergoes a rapid change from the Rouse-like $η\sim{Z}^1_w$ to a very strong power-law scaling that is approximately $η\sim{Z}^{5.8}_w$. This aspect of the work strongly suggests that prior investigations on highly pure rings have invariably been performed on unentangled systems simply because ring molecules have very high entanglement thresholds. The results are further supported by plateau modulus measurements that show that in the "Rouse-like" viscosity regime, for the RDODT samples there is a corresponding lack of a rubbery entanglement plateau and once entanglement coupling is seen, the rubbery plateau also appears in the dynamic modulus data. Furthermore, upon appropriate scaling, we show that the unentangled polyDODT rings (i.e., the diluted high molecular weight rings) show virtually identical dynamics (upon appropriate scaling to account for the chain dilutions) to the dynamics of similarly unentangled LCCC fractionated polystyrene rings.

36 MATERIALS SCIENCE↗
NASA NTRS1997

Scale Model Icing Research Tunnel

NASA Lewis Research Center's Icing Research Tunnel (IRT) is the world's largest refrigerated wind tunnel and one of only three icing wind tunnel facilities in the United States. The IRT was constructed in the 1940's and has been operated continually since it was built. In this facility, natural icing conditions are duplicated to test the effects of inflight icing on actual aircraft components as well as on models of airplanes and helicopters. IRT tests have been used successfully to reduce flight test hours for the certification of ice-detection instrumentation and ice protection systems. To ensure that the IRT will remain the world's premier icing facility well into the next century, Lewis is making some renovations and is planning others. These improvements include modernizing the control room, replacing the fan blades with new ones to increase the test section maximum velocity to 430 mph, installing new spray bars to increase the size and uniformity of the artificial icing cloud, and replacing the facility heat exchanger. Most of the improvements will have a first-order effect on the IRT's airflow quality. To help us understand these effects and evaluate potential improvements to the flow characteristics of the IRT, we built a modular 1/10th-scale aerodynamic model of the facility. This closed-loop scale-model pilot tunnel was fabricated onsite in the various shops of Lewis' Fabrication Support Division. The tunnel's rectangular sections are composed of acrylic walls supported by an aluminum angle framework. Its turning vanes are made of tubing machined to the contour of the IRT turning vanes. The fan leg of the tunnel, which transitions from rectangular to circular and back to rectangular cross sections, is fabricated of fiberglass sections. The contraction section of the tunnel is constructed from sheet aluminum. A 12-bladed aluminum fan is coupled to a turbine powered by high-pressure air capable of driving the maximum test section velocity to 550 ft/sec (Mach 0.45). The air turbine and instrumentation are housed inside a fiberglass nacelle. Total and static pressure measurements can be taken around the loop, and velocity and flow angularity measurements can be taken with hot-wire and five-hole probes at specific locations. The Scale Model Icing Research Tunnel (SMIRT) is undergoing checkout tests to determine how its airflow characteristics compare with the IRT. Near-term uses for this scale-model tunnel include determining the aerodynamic effects of replacing the 52-yearold W-shaped heat exchanger with a flat-faced heat exchanger. SMIRT is an integral part of the improvements planned for the IRT because testing the proposed IRT improvements in a scale-model tunnel will lower costs and improve productivity.

Canacci, Victor A.↗
DOE OSTI2020

Progression to Compatibility Evaluations in Flowing Molten Salts

Molten salt compatibility with structural alloys has been identified as a key issue for the development of Generation 3 concentrating solar power (CSP) systems with thermal storage. To accelerate this evaluation to pumped systems, the goal of this project was to conduct thermal convection loop (TCL) experiments with a peak temperature of ≥700°C and a typical temperature gradient of ~100°C. The experiments indicated that conventional ~16wt.%Cr Ni-based alloys are compatible up to 700°C with purified (i.e. low O) or dried (low H₂O) industrial-sourced Mg-K-Na chloride salt with <10 μm/yr loss. This two-year project was conducted based on the experimental and mechanistic understanding developed more than 60 years ago at Oak Ridge National Laboratory (ORNL). The first TCL experiment met the <15 μm/yr corrosion metric for this project with specimens of Ni-based alloy 600 exposed at 580°-700°C for 1000 h and post-exposure room temperature tensile tests showed minimal degradation. The second TCL experiment successfully deployed an electrochemical sensor from Argonne National Laboratory and had a peak temperature of 750°C but only ran for ~110 h due to a furnace failure. Both experiments used highly purified industrial-sourced salt with an O content of ~3 μg O/g salt and a Mg addition of 0.04 wt.%. In the second year, the Chloride Collective developed a more economical drying procedure such that the O content was much higher (>20,000 μg O/g salt). A third TCL experiment was conducted with a sensor and a peak temperature of 700°C using dried salt from the same industrial source and increased NaCl content (~20 wt.%). In addition to a 0.05% Mg addition to the salt, a Mg coupon was added in the coldest part of the loop which dissolved during the experiment. Again, small mass changes were noted for specimens of alloys 600 and C276 but the values were slightly higher than those measured in purified salt with Mg. A thin non-continuous and non-adherent oxide layer was deposited on most specimens containing Mg, Si and Al but Cr depletion also was observed. Both years included facilities qualification crucible experiments and then capsule experiments to confirm a baseline isothermal reaction rate. The first year capsule experiment led to the conclusion that the two-stage ORNL purification process using NH₄Cl and CCl₄ left the salt with a high Cl potential and a Mg addition (~0.05wt.%) was needed to lower the potential. The second year capsule experiments at 600° and 700°C found little difference in depth of attack for 0-0.25%Mg additions. In general, the complex reactions where salt can be trapped in the porous surface layer of metal indicated that mass change is an unreliable metric and average depth of Cr depletion is a better metric for assessing the extent of attack. These results have created a new baseline that is contrary to the recent published literature for chloride salts where mass losses have been reported that can be extrapolated to very significant annual metal loss rates. Chloride salt corrosion can be controlled and the results also indicate that salt purification to low O levels may not be necessary. However, additional TCL experiments are needed at different times and temperatures to generate important engineering information such as temperature dependent corrosion rates and reaction rate laws for extrapolation to long-term behavior and isolate the effect of salt additives and impurities.

36 MATERIALS SCIENCE↗
NASA NTRS2004

Chemical Reaction and Flow Modeling in Fullerene and Nanotube Production

The development of processes to produce fullerenes and carbon nanotubes has largely been empirical. Fullerenes were first discovered in the soot produced by laser ablation of graphite [1]and then in the soot of electric arc evaporated carbon. Techniques and conditions for producing larger and larger quantities of fullerenes depended mainly on trial and error empirical variations of these processes, with attempts to scale them up by using larger electrodes and targets and higher power. Various concepts of how fullerenes and carbon nanotubes were formed were put forth, but very little was done based on chemical kinetics of the reactions. This was mainly due to the complex mixture of species and complex nature of conditions in the reactors. Temperatures in the reactors varied from several thousand degrees Kelvin down to near room temperature. There are hundreds of species possible, ranging from atomic carbon to large clusters of carbonaceous soot, and metallic catalyst atoms to metal clusters, to complexes of metals and carbon. Most of the chemical kinetics of the reactions and the thermodynamic properties of clusters and complexes have only been approximated. In addition, flow conditions in the reactors are transient or unsteady, and three dimensional, with steep spatial gradients of temperature and species concentrations. All these factors make computational simulations of reactors very complex and challenging. This article addresses the development of the chemical reaction involved in fullerene production and extends this to production of carbon nanotubes by the laser ablation/oven process and by the electric arc evaporation process. In addition, the high-pressure carbon monoxide (HiPco) process is discussed. The article is in several parts. The first one addresses the thermochemical aspects of modeling; and considers the development of chemical rate equations, estimates of reaction rates, and thermodynamic properties where they are available. The second part addresses modeling of the arc process for fullerene and carbon nanotube production using O-D, 1-D and 2-D fluid flow models. The third part addresses simulations of the pulsed laser ablation process using time-dependent techniques in 2-D, and a steady state 2-D simulation of a continuous laser ablation process. The fourth part addresses steady state modeling in O-D and 2-D of the HiPco process. In each of the simulations, there is a variety of simplifications that are made that enable one to concentrate on one aspect or another of the process. There are simplifications that can be made to the chemical reaction models , e.g. reduction in number of species by lumping some of them together in a representative species. Other simulations are carried out by eliminating the chemistry altogether in order to concentrate on the fluid dynamics. When solving problems with a large number of species in more than one spatial dimension, it is almost imperative that the problem be decoupled by solving for the fluid dynamics to find the fluid motion and temperature history of "particles" of fluid moving through a reactor. Then one can solve the chemical rate equations with complex chemistry following the temperature and pressure history. One difficulty is that often mixing with an ambient gas is involved. Therefore, one needs to take dilution and mixing into account. This changes the ratio of carbon species to background gas. Commercially available codes may have no provision for including dilution as part of the input. One must the write special solvers for including dilution in decoupled problems. The article addresses both ful1erene production and single-walled carbon nanotube (SWNT) production. There are at least two schemes or concepts of SWNT growth. This article will only address growth in the gas phase by carbon and catalyst cluster growth and SW T formation by the addition of carbon. There are other models that conceive of SWNT growth as a phase separation process from clusters me up carbon and metal catalyst, with the carbon precipitating from the cluster as it cools. We will not deal with that concept in this article. Further research is needed to determine the rates at which these composite clusters form, evaporate, and segregate.

Scott, Carl D.↗
DOE OSTI2022

Development of an Integrated Multidiagnostic to Assess the High-Z Impurity Fluxes in the Metallic Environment of WEST Using IMAS

WEST is an actively cooled, long-pulse tokamak with nearly all plasma-facing components (PFC) made of tungsten. One of the aims of WEST is to study plasma operations with tungsten PFCs in preparation for long-pulse operations on high-Z divertor devices, such as ITER. For long-pulse operation, the high-Z impurity content and transport to the core plasma are critical concerns that require further measurement and interpretation to improve plasma performance and PFC durability. This work details the impurity influxes in WEST during a series of discharges in which the lower hybrid (LH) injected power was incrementally increased. An analysis has been performed of measurements collected from an array of edge diagnostics. Visible spectroscopy was utilized to measure the spectral radiances generated by fuel particles (D) and impurities (W, O, and C) at the divertor and at the antennas with a newly developed spectral peak-fitting tool used to analyze the data in WEST. The scrape-off layer (SOL) plasma conditions (n e and T e ) measured at the divertor target with flush Langmuir probes and near the outer mid-plane (OMP) using reciprocating Langmuir probes (RCPs) are used to evaluate the number of ionizations per photons (S/XB) coefficients required to estimate the impurity fluxes obtained with the collisional-radiative model ColRadPy. Here, the array of edge diagnostics discussed in this work, coupled with SOL plasma modeling tools, represents a multidiagnostic interpretative modeling workflow that will continue to be applied to upcoming experimental campaigns on the WEST experiment to assess high-Z impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2009

High-Efficiency Artificial Photosynthesis Using a Novel Alkaline Membrane Cell

A new cell designed to mimic the photosynthetic processes of plants to convert carbon dioxide into carbonaceous products and oxygen at high efficiency, has an improved configuration using a polymer membrane electrolyte and an alkaline medium. This increases efficiency of the artificial photosynthetic process, achieves high conversion rates, permits the use of inexpensive catalysts, and widens the range of products generated by this type of process. The alkaline membrane electrolyte allows for the continuous generation of sodium formate without the need for any additional separation system. The electrolyte type, pH, electrocatalyst type, and cell voltage were found to have a strong effect on the efficiency of conversion of carbon dioxide to formate. Indium electrodes were found to have higher conversion efficiency compared to lead. Bicarbonate electrolyte offers higher conversion efficiency and higher rates than water solutions saturated with carbon dioxide. pH values between 8 and 9 lead to the maximum values of efficiency. The operating cell voltage of 2.5 V, or higher, ensures conversion of the carbon dioxide to formate, although the hydrogen evolution reaction begins to compete strongly with the formate production reaction at higher cell voltages. Formate is produced at indium and lead electrodes at a conversion efficiency of 48 mg of CO2/kilojoule of energy input. This efficiency is about eight times that of natural photosynthesis in green plants. The electrochemical method of artificial photosynthesis is a promising approach for the conversion, separation and sequestration of carbon dioxide for confined environments as in space habitats, and also for carbon dioxide management in the terrestrial context. The heart of the reactor is a membrane cell fabricated from an alkaline polymer electrolyte membrane and catalyst- coated electrodes. This cell is assembled and held in compression in gold-plated hardware. The cathode side of the cell is supplied with carbon dioxide-saturated water or bicarbonate solution. The anode side of the cell is supplied with sodium hydroxide solution. The solutions are circulated past the electrodes in the electrochemical cell using pumps. A regulated power supply provides the electrical energy required for the reactions. Photovoltaic cells can be used to better mimic the photosynthetic reaction. The current flowing through the electrochemical cell, and the cell voltage, are monitored during experimentation. The products of the electrochemical reduction of carbon dioxide are allowed to accumulate in the cathode reservoir. Samples of the cathode solution are withdrawn for product analysis. Oxygen is generated on the anode side and is allowed to vent out of the reservoir.

Narayan, Sri↗
NASA NTRS1994

Alkali Metal Rankine Cycle Boiler Technology Challenges and Some Potential Solutions for Space Nuclear Power and Propulsion Applications

Alkali metal boilers are of interest for application to future space Rankine cycle power conversion systems. Significant progress on such boilers was accomplished in the 1960's and early 1970's, but development was not continued to operational systems since NASA's plans for future space missions were drastically curtailed in the early 1970's. In particular, piloted Mars missions were indefinitely deferred. With the announcement of the Space Exploration Initiative (SEI) in July 1989 by President Bush, interest was rekindled in challenging space missions and, consequently in space nuclear power and propulsion. Nuclear electric propulsion (NEP) and nuclear thermal propulsion (NTP) were proposed for interplanetary space vehicles, particularly for Mars missions. The potassium Rankine power conversion cycle became of interest to provide electric power for NEP vehicles and for 'dual-mode' NTP vehicles, where the same reactor could be used directly for propulsion and (with an additional coolant loop) for power. Although the boiler is not a major contributor to system mass, it is of critical importance because of its interaction with the rest of the power conversion system; it can cause problems for other components such as excess liquid droplets entering the turbine, thereby reducing its life, or more critically, it can drive instabilities-some severe enough to cause system failure. Funding for the SEI and its associated technology program from 1990 to 1993 was not sufficient to support significant new work on Rankine cycle boilers for space applications. In Fiscal Year 1994, funding for these challenging missions and technologies has again been curtailed, and planning for the future is very uncertain. The purpose of this paper is to review the technologies developed in the 1960's and 1970's in the light of the recent SEI applications. In this way, future Rankine cycle boiler programs may be conducted most efficiently. This report is aimed at evaluating alkali metal boiler technology for space Rankine cycle systems. Research is summarized on the problems of flow stability, liquid carryover, pressure drop and heat transfer, and on potential solutions developed, primarily those developed by the NASA Lewis Research Center in the 1960's and early 1970's.

Stone, James R.↗
NASA NTRS1993

A long-term space astrophysics research program: The evolution of the quasar continuum

The research program supported by this grant now has great momentum. Numerous papers are in progress, and a strong multi-wavelength observing program is rapidly accumulating data on samples of high redshift quasars across the spectrum. ROSAT spectra of quasars continue to yield surprises. Of four z = 3 quasars with X-ray spectra, three show strong absorption. This contrasts strongly with the situation for luminous AGN at low redshifts where fewer than 1 in 20 show X-ray absorption. A new site for this absorption is probably needed, either around the quasar (e.g. in a cluster cooling flow) or along the line of sight (e.g. in a Damped Lyman-alpha system). The unabsorbed quasar allows limits on the physical conditions in a damped Lyman-alpha cloud to be calculated, and will allow a X-ray Gunn-Peterson test to be applied that will limit the fraction of the closure mass in an intergalactic medium. The X-ray spectral indices of these z = 3 quasars show no change from those of similar objects at low z, suggesting that 'short-lifetime' models apply. Eight other z = 3-4 quasars have been detected and their energy distributions from X-rays to Infrared (using new infrared spectrographs) have been compiled. These are now being compared with the low z continua from the 'Atlas of Quasar Energy Distributions' to search for evolutionary changes. The discovery of a likely warm absorber in 3C351 made recognition of another example simple. Also, modeling of the conditions in the absorber in 3C351 using the OVI absorption line from HST and the high ionization emission lines, suggests that the broad line region is indeed the origin of the warm absorber in this quasar, and by extension, others. Warm absorbers can now be used as a new diagnostic of this region. The X-ray spectrum of a 'Red Quasar', 3C212, has a cut-off spectrum, which could be fitted by an absorbed power-law, or more remarkably, by an unabsorbed black body. Using our quasi-simultaneous optical data and photoionization modeling we are able to rule out this latter possibility. The long-sought key to understanding red quasars has now been found. (Reddening was earlier ruled out because no 2175A absorption feature was seen. This is now seen to have been misleading.)

Elvis, M.↗
NASA NTRS2013

Passivation of Flexible YBCO Superconducting Current Lead With Amorphous SiO2 Layer

Adiabatic demagnetization refrigerators (ADR) are operated in space to cool detectors of cosmic radiation to a few 10s of mK. A key element of the ADR is a superconducting magnet operating at about 0.3 K that is continually energized and de-energized in synchronism with a thermal switch, such that a piece of paramagnetic salt is alternately warm in a high magnetic field and cold in zero magnetic field. This causes the salt pill or refrigerant to cool, and it is able to suck heat from an object, e.g., the sensor, to be cooled. Current has to be fed into and out of the magnets from a dissipative power supply at the ambient temperature of the spacecraft. The current leads that link the magnets to the power supply inevitably conduct a significant amount of heat into the colder regions of the supporting cryostat, resulting in the need for larger, heavier, and more powerful supporting refrigerators. The aim of this project was to design and construct high-temperature superconductor (HTS) leads from YBCO (yttrium barium copper oxide) composite conductors to reduce the heat load significantly in the temperature regime below the critical temperature of YBCO. The magnet lead does not have to support current in the event that the YBCO ceases to be superconducting. Cus - tomarily, a normal metal conductor in parallel with the YBCO is a necessary part of the lead structure to allow for this upset condition; however, for this application, the normal metal can be dispensed with. Amorphous silicon dioxide is deposited directly onto the surface of YBCO, which resides on a flexible substrate. The silicon dioxide protects the YBCO from chemically reacting with atmospheric water and carbon dioxide, thus preserving the superconducting properties of the YBCO. The customary protective coating for flexible YBCO conductors is silver or a silver/gold alloy, which conducts heat many orders of magnitude better than SiO2 and so limits the use of such a composite conductor for passing current across a thermal gradient with as little flow of heat as possible to make an efficient current lead. By protecting YBCO on a flexible substrate of low thermal conductivity with SiO2, a thermally efficient and flexible current lead can be fabricated. The technology is also applicable to current leads for 4 K superconducting electronics current biasing. A commercially available thin-film YBCO composite tape conductor is first stripped of its protective silver coating. It is then mounted on a jig that holds the sample flat and acts as a heat sink. Silicon dioxide is then deposited onto the YBCO to a thickness of about 1 micron using PECVD (plasma-enhanced chemical vapor deposition), without heating the YBCO to the point where degradation occurs. Since SiO2 can have good high-frequency electrical properties, it can be used to coat YBCO cable structures used to feed RF signals across temperature gradients. The prime embodiment concerns the conduction of DC current across the cryogenic temperature gradient. The coating is hard and electrically insulating, but flexible.

Johannes, Daniel↗
DOE OSTI2023

Hybrid Tandem Photovoltaics

Tandem solar cell structures are the only strategy demonstrated to surpass the detailed balance efficiency limit of high-quality single-junction solar cells. To continue to improve the efficiencies of cost-effective terrestrial solar power, hybrid tandems of dissimilar subcells are being considered by many around the world, especially designs that incorporate silicon solar cells as a bottom subcell. In this project, we studied a wide variety of tandem design possibilities including those with three-terminal (3T) and four-terminal (4T) configurations. The use of 3T and 4T designs could be useful for efficient and economical hybrid tandem designs that utilize the best available subcell materials such as emerging perovskite materials. Three-terminal configurations, in particular, have not been sufficiently studied previously. We have laid the foundational groundwork in this project for understanding the operation of 3T tandems: developing a taxonomy for naming, a methodology for measuring and interconnecting, and models for simply characterizing 3T tandems. Electrical and optical subcell coupling between the subcells was also measured and modeled. An important part of this work was the fabrication of novel example tandem structures, including 4T GaAs/Si, 3T GaInP/Si, 3T GaAs/Si, and 3T GaInP/GaAs devices. Using these high-quality tandem cells, we have been able to clearly demonstrate the achievability of high-efficiencies, and subtle physical effects such as photon recycling and luminescent coupling. We have developed and demonstrated essential building-block tools such as transparent conductive adhesives (TCA) and 3T silicon bottom cells with interdigitated back contacts (IBC) that can also be used in many other tandem designs. We have tested the reliability of these tools and devices under standardized testing and outdoor measurements. We have found 4T GaAs/Si tandems to be relatively straightforward to fabricate and robust in real-world outdoor conditions. While we have demonstrated working hybrid 3T III-V/TCA/Si IBC tandems, we experienced low yields even with our best process flows yet. Further work is still needed to improve the processing yield of these devices. We therefore also created tandem cells using an all-III-V 3T tandem process which was very robust with high yields, allowing for the creation of voltage-matched strings in many different configurations using 8 nearly identical 3T tandems. Using these robust 3T tandem examples, we were able measure and precisely characterize 3T tandem behaviors to predict their operation under changing spectrum and temperature. The optoelectronic equivalent-circuit model was shown to be very general and applicable to hybrid tandems, and encompassed the operation 3T Si IBC cells. This general model has been distributed to the public in as open-source Python-based software called PVcircuit. We have calculated the implications of these new tandem device designs on the real-world energy production and shown how the relative performance of different tandem configurations is situational and can be engineered using the tools developed here.

14 SOLAR ENERGY↗