Search NASA⌕ Search

SEARCH · Search NASA

Results for “Evolution operators”

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 1,045 records · Page 58

Systematically tracking the hourly progression of large wildfires using GOES satellite observations

Abstract. In the western United States, prolonged drought, a warming climate, and historical fuel buildup have contributed to larger and more intense wildfires as well as to longer fire seasons. As these costly wildfires become more common, new tools and methods are essential for improving our understanding of the evolution of fires and how extreme weather conditions, including heat waves, windstorms, droughts, and varying levels of active-fire suppression, influence fire spread. Here, we develop the Geostationary Operational Environmental Satellites (GOES)-Observed Fire Event Representation (GOFER) algorithm to derive the hourly fire progression of large wildfires and create a product of hourly fire perimeters, active-fire lines, and fire spread rates. Using GOES-East and GOES-West geostationary satellite detections of active fires, we test the GOFER algorithm on 28 large wildfires in California from 2019 to 2021. The GOFER algorithm includes parameter optimizations for defining the burned-to-unburned boundary and correcting for the parallax effect from elevated terrain. We evaluate GOFER perimeters using 12 h data from the Visible Infrared Imaging Radiometer Suite (VIIRS)-derived Fire Event Data Suite (FEDS) and final fire perimeters from the California's Fire and Resource Assessment Program (FRAP). Although the GOES imagery used to derive GOFER has a coarser resolution (2 km at the Equator), the final fire perimeters from GOFER correspond reasonably well to those obtained from FRAP, with a mean Intersection-over-Union (IoU) of 0.77, in comparison to 0.83 between FEDS and FRAP; the IoU indicates the area of overlap over the area of the union relative to the reference perimeters, in which 0 is no agreement and 1 is perfect agreement. GOFER fills a key temporal gap present in other fire tracking products that rely on low-Earth-orbit imagery, where perimeters are available at intervals of 12 h or longer or at ad hoc intervals from aircraft overflights. This is particularly relevant when a fire spreads rapidly, such as at maximum hourly spread rates of over 5 km h−1. Our GOFER algorithm for deriving the hourly fire progression using GOES can be applied to large wildfires across North and South America and reveals considerable variability in the rates of fire spread on diurnal timescales. The resulting GOFER product has a broad set of potential applications, including the development of predictive models for fire spread and the improvement of atmospheric transport models for surface smoke estimates. The resulting GOFER product has a broad set of potential applications, including the development of predictive models for fire spread and the improvement of atmospheric transport models for surface smoke estimates (https://doi.org/10.5281/zenodo.8327264, Liu et al., 2023).

Geology↗

Time at Temperature Experiments on Neutron Irradiated Zircaloy-2 using Conventional & Flash DSC

Dryout events in Boiling Water Reactors (BWRs) are currently treated by NRC regulations as automatic disqualification for continued fuel rod operation, even though this criterion does not account for the rate or duration of power increases, associated changes in material behavior, or the possibility of rewetting. Operational history from Anticipated Operational Occurrences (AOOs) shows that short, transient power excursions often demand only modest heat removal, and industry experience suggests that fuel can briefly enter dryout yet return to safe, stable operation. The lack of detailed understanding of the material response during such events motivates the present series of experiments. This study uses unique and innovative methods to investigate microstructural evolution in irradiated Zircaloy-2 exposed to high temperatures in inert environments. Differential Scanning Calorimetry (DSC) and FlashDSC are combined to build a comprehensive experimental framework capable of identifying the a–ß phase transformation in zirconium and examining defect annealing under steep thermal gradients. FlashDSC enables rapid heating and cooling of Focused Ion Beam (FIB)–prepared large-area lift-outs (LALOs) of irradiated Zircaloy-2 at rates of 1,000 K/s to peak temperatures of 600°C, 750°C, and 900°C. The goal is to determine whether these conditions produce measurable microstructural changes that could influence cladding performance in typical BWR environments. Microstructural characterization includes quantifying dislocation density and assessing secondary phase particle (SPP) size and distribution using Transmission Electron Microscopy (TEM). Ongoing analysis, such as diffraction pattern indexing and 4D STEM processing, will further refine these observations.

36 - MATERIALS SCIENCE↗

Geodetic imaging of thermal deformation in geothermal reservoirs - production, depletion and fault reactivation

Here, we investigate thermally induced surface deformation in geothermal systems. To define source mechanisms at depth, we assess the mechanical process of subsurface deformation by assuming a spherically cooled fractured reservoir in an infinite medium and derive relations that define magnitudes of thermal contraction, stress change and permeability evolution. The magnitude of thermal deformation in typical geothermal system is larger than anticipated and suggests two different modalities of surface subsidence – thermal contraction and fault reactivation. Here, surface deformation (vertical displacement, surface tilt and horizontal strain) induced by the two different modalities are assessed with Mogi (contraction) and Okada (slip) models and compared with instrumental sensitivity of high precision surface geodetic tools. Here, we show that 1 year of geothermal operation at 10 MW with a power plant conversion efficiency of 12% can yield ~ 3.0 × 10 4 m 3 of subsurface volume change. For a reservoir at 2000 m depth, this induces ~ 1.7 mm of vertical surface displacement, ~ 800 nano-radians of surface tilt and ~ 900 nano-strains of surface strain. This result implies that typically observed magnitudes of surface subsidence (order of cm/year) are naturally expected in massive (100 MW scale) geothermal operations and observed surface subsidence may largely be the result of thermal contraction. Conversely, thermal unloading can trigger fault reactivation. Analysis with an Okada slip model shows these shear offsets on pre-existing faults can also result in surface deformations of considerable magnitude. Our analysis of field operational data from various geothermal projects suggests that both thermal contraction and slow fault reactivation may contribute to the observed large surface deformation. Comparison of predicted deformation with instrumental sensitivity of high precision surface tools confirms that geodetic signals, especially tilt and strain, are indeed sufficiently large to describe reservoir evolution and to potentially deconvolve reservoir parameters of interest, such as permeability.

15 GEOTHERMAL ENERGY↗

Evolution of the Rate-Limiting step: From Thin Film to Thick Ni-rich Cathodes

Ni-rich NMC (LiNixMnyCo1-x-yO2, x = 0.6) has the attributes of high specific capacity (greater than 180 mAh g-1), high operating voltage (ca. 3.8 V) and low cost, which is therefore deemed as one of the most promising cathode candidates for next-generation high-energy Li-ion batteries (LIBs). However, challenges exist before the large-scale commercialization of Ni-rich cathode such as gas generation, capacity degradation especially at high rate or elevated temperatures. From thin film to densely pack high loading cathode, this work investigates the rate-limiting steps in Ni-rich cathodes e.g., phase structure, electronic and ionic conductivity to understand the reaction kinetics at different scales. This work provides new insights in understanding the electrochemical properties of Ni-rich cathode at different levels to inspire revolutionary ideas to address the fundamental issues of Ni-rich cathode from materials to relevant electrode levels.

Hu, Jiangtao↗

Fault-tolerant resource comparison of qudit and qubit encodings for diagonal quadratic operators

Finite local Hilbert-space truncations arise naturally in quantum simulations of lattice field theories and motivate qudit encodings, but their fault-tolerant advantage over qubit encodings remains unclear. We compare the non-Clifford cost of implementing quadratic diagonal evolutions, exemplified by 𝑈 = 𝑒$^{−𝑖⁢𝑡⁢𝜙^2_𝑥}$ in a uniform field-amplitude discretization of a real scalar field, using either one logical 𝑑-level qudit or 𝑛 𝑏 = ⌈log 2⁡ 𝑑⌉ logical qubits. We analyze two standard settings: product-formula simulation and linear combination of unitaries (LCU) per block encoding, taking the resource metric to be the number of non-Clifford gates after synthesis into a discrete logical gate set. Because tight synthesis bounds for general single-qudit rotations are not known, we express the qudit constructions in terms of embedded two-level SU⁡(2) rotations and derive explicit finite-𝑑 break-even conditions for their synthesis cost; these serve as compiler targets for when qudit encodings can outperform the qubit baseline. Within the constructive models studied here, product-formula implementations would require an exponentially stronger per-primitive synthesis advantage for qudits to win asymptotically, while in the LCU setting the qubit encoding is asymptotically cheaper in 𝑑. Nevertheless, the finite-𝑑 threshold analysis identifies low-dimensional regions in which qudits can yield meaningful constant-factor savings, particularly for LCU-based implementations. As a secondary analysis of the LCU construction, we use an idealized negligible-overhead qubit-qudit code-switching model to give an absolute 𝑇-count comparison and reinterpret the savings as an allowable per-switch overhead budget.

Godwood, Samuel [Univ. of Liverpool (United Kingdo↗

Co-dependent microstructural evolution pathways in metastable d-ferrite in cast austenitic stainless steels during thermal aging

Cast austenitic stainless steels (CASS) are excellent alloys because they combine high corrosion resistance with high strength and toughness. However, they are susceptible to embrittlement upon long-term thermal aging at elevated temperatures. Thus, the microstructural evolution pathways during thermal aging need to be understood to predict and potentially prevent embrittlement. Atom probe tomography was used to identify and quantify the microstructural evolution pathways of the d-ferrite in different CASS alloys aged for up to 10,000 hours at temperatures between 290 °C and 400 °C. The four steels – CF8, CF8M, CF3, and CF3M – which vary by Mo and C concentration, each experienced spinodal decomposition of the d-ferrite, and precipitation of G-phase clusters and Cu clusters attached to the G-phase. There were large differences in the extent of these features due to their Mo and C concentration. Using radial distribution function analysis, the interactions of constituent elements was found to determine the evolution of these features, with Mo and C specifically influencing the movement of Cr, Ni, Si, Mn, and Cu atoms due to their relative miscibility with these elements. The results will help inform predictive models for the use of duplex stainless steels for extended operation at high temperatures.

Duplex stainless steel, Aging, Spinodal decomposit↗

Synthesis and hydrogen-evolution activity of tungsten selenide thin films deposited on tungsten foils

Thin films of WSe 2 have been deposited onto a conductive substrate (tungsten foil) using a relatively simple chemical-vapor-transport technique. X-ray photoelectron spectroscopy, energy-dispersive X-ray spectroscopy, X-ray powder diffraction, scanning electron microscopy, and high-resolution transmission electron microscopy indicated that the films consisted of micron-sized single crystals of WSe 2 that were oriented perpendicular to the surface of the tungsten foil substrate. Linear sweep voltammetry was used to assess the ability of the WSe 2 films to catalyze the hydrogen-evolution reaction and chronopotentiometry was used to gauge the temporal stability of the catalytic performance of the films under cathodic conditions. A 350 mV overpotential (η) was required to drive the hydrogen-evolution reaction at a current density of –10 mA cm –2 in aqueous 0.5 M H 2 SO 4 , representing a significant improvement in catalytic performance relative to the behavior of macroscopic WSe 2 single crystals. The WSe 2 thin films were relatively stable under catalytic conditions, with the overpotential changing by only ~10 mV after one hour and exhibiting an additional change of ~5 mV after another hour of operation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding the Angular Momentum Evolution of T Tauri and Herbig Ae/Be Stars

We investigate a sample of six Herbig Ae/Be stars belonging to the Orion OB1 association, as well as 73 low-mass objects, members of the σ Orionis cluster, in order to explore the angular momentum evolution at early stages of evolution, and its possible connection with main-sequence Ap/Bp magnetic stars. Using FIES and HECTOCHELLE spectra, we obtain projected rotational velocities through two independent methods. Individual masses, radii, and ages are computed using evolutionary models, distance, and cluster extinction. Under the assumption that similar physical processes operate in both T Tauri and Herbig Ae/Be stars, we construct snapshots of the protostar’s rotation against mass during the first 10 Myr with the aid of a rotational model that includes a variable disk lifetime, changes in the stellar moment of inertia, a dipolar magnetic field with variable strength, and angular momentum loss through stellar winds powered by accretion. We use these snapshots, as well as the rotational data, to infer a plausible scenario for the angular momentum evolution. We find that magnetic field strengths of a few kilo-Gauss at 3 Myr are required to match the rotational velocities of both groups of stars. Models with masses between 2 M {sub ⊙} and 3 M {sub ⊙} display larger angular momentum values by a factor of ~3, in comparison to stars of similar spectral types on the main sequence. Even though some quantitative estimates on this dramatic decrease with age for Ap/Bp magnetic main-sequence stars are presented, the results obtained for the angular momentum evolution do not explain their low rotation rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Forecasting for ESCAPE: A Multi-Institution Hybrid Forecasting and Nowcasting Operation for Sea-Breeze Convection Supporting a Ground-Based and Airborne Field Campaign

The Experiment of Sea-Breeze Convection, Aerosols, Precipitation and Environment (ESCAPE) field project deployed two aircraft and ground-based assets in the vicinity of Houston, Texas, between 27 May and 2 July 2022, examining how meteorological conditions, dynamics, and aerosols control the initiation, early growth stage, and evolution of coastal convective clouds. To ensure that airborne- and ground-based assets were deployed appropriately, a forecasting and nowcasting team was formed. Daily forecasts guided real-time decision-making by assessing synoptic weather conditions, environmental aerosol, and a variety of atmospheric modeling data to assign a probability for meeting specific ESCAPE campaign objectives. During the research flights, a small team of forecasters provided “nowcasting” support by analyzing radar, satellite, and new model data in real time. The nowcasting team proved invaluable to the campaign operation, as sometimes changing environmental conditions affected, for example, the timing of convective initiation. In addition to the success of the forecasting and nowcasting teams, the ESCAPE campaign offered a unique “testbed” opportunity where in-person and virtual support both contributed to campaign objectives. The forecasting and nowcasting teams were each composed of new and experienced forecasters alike, where new forecasters were given invaluable experience that would otherwise be difficult to attain. Both teams received training on forecast models, map analysis, Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT), and thermodynamic sounding analysis before the beginning of the campaign. In this article, the ESCAPE forecasting and nowcasting teams reflect on these experiences, providing potentially useful advice for future field campaigns requiring forecasting and nowcasting support in a hybrid virtual/in-person framework.

54 ENVIRONMENTAL SCIENCES↗

Predicting nonequilibrium Green’s function dynamics and photoemission spectra via nonlinear integral operator learning

Understanding the dynamics of nonequilibrium quantum many-body systems is an important research topic in a wide range of fields across condensed matter physics, quantum optics, and high-energy physics. However, numerical studies of large-scale nonequilibrium phenomena in realistic materials face serious challenges due to intrinsic high-dimensionality of quantum many-body problems and the absence of time-invariance. The nonequilibrium properties of many-body systems can be described by the dynamics of the correlator, or the Green's function of the system, whose time evolution is given by a high-dimensional system of integro-differential equations, known as the Kadanoff–Baym equations (KBEs). The time-convolution term in KBEs, which needs to be recalculated at each time step, makes it difficult to perform long-time numerical simulation. In this paper, we develop an operator-learning framework based on recurrent neural networks (RNNs) to address this challenge. We utilize RNNs to learn the nonlinear mapping between Green's functions and convolution integrals in KBEs. By using the learned operators as a surrogate model in the KBE solver, we obtain a general machine-learning scheme for predicting the dynamics of nonequilibrium Green's functions. Besides significant savings per each time step, the new methodology reduces the temporal computational complexity from $O(N_t^3)$ to $O(N_t)$ where N t is the number of steps taken in a simulation, thereby making it possible to study large many-body problems which are currently infeasible with conventional KBE solvers. Through various numerical examples, we demonstrate the effectiveness of the operator-learning based approach in providing accurate predictions of physical observables such as the reduced density matrix and time-resolved photoemission spectra. Moreover, our framework exhibits clear numerical convergence and can be easily parallelized, thereby facilitating many possible further developments and applications.

97 MATHEMATICS AND COMPUTING↗

Shapeshifting Nanocatalyst for CO2 Conversion

The conversion of CO2 into high-value chemicals through a photoreduction reaction in water is a promising route to reduce the dependence on fossil fuels. Enhancing selectivity toward hydrocarbons or alcohols can be achieved by Ag-Cu alloys. However, the stabilized surface state created by Ag-Cu interactions is still poorly understood. In this work, multi-modal in situ X-ray experiments reveals underlying mechanisms and the evolution of Ag-Cu nanoparticles under CO2 reduction reaction (CO2RR) conditions. Both morphological and chemical changes of Ag and Cu species induced by diffusion mechanics are tracked during nanocatalyst operation. The initial spheroid Ag-Cu nanoparticles are composed of a Cu-rich shell and Ag-rich core. The reduction treatment promotes Ag migration toward the surface. During photocatalytic CO2 reduction reaction, Cu atoms migrate back to the surface, forming Ag-Cu-O species. The study observes the surface oxidation of Cu(0) to Cu+ and the presence of Ag at the sub-surface region. Furthermore, nanoparticles change their shape, decreasing their specific surface area, driven by Cu diffusion during the CO2 photoreduction reaction. The results provide invaluable insights into the dynamic restructuring of the catalyst under reaction conditions and into the active species responsible for CO2 conversion.

CO2 reduction reaction↗

Search, engineering, and applications of new oxidative biocatalysts

Abstract Most industrial enzymes are hydrolases, such as glycosidases and esterases. However, oxidoreductases have an unexploited potential for substituting harsh (and scarcely selective) chemical processes. A group of basidiomycetes are the only organisms degrading the aromatic lignin polymer, enabling the subsequent use of plant polysaccharides. Therefore, these fungi and their ligninolytic peroxidases are the biocatalysts of choice for industrial delignification and oxidative biotransformations of aromatic and other organic compounds. The latter also include oxygenation reactions, which are catalyzed with high regio/stereo selectivity by fungal peroxygenases. In search for novel and more robust peroxidases/peroxygenases, basidiomycetes from unexplored habitats were screened, and hundreds of genes identified in basidiomycete genomes (in collaboration with the DOE JGI). The most interesting genes were heterologously expressed, and the corresponding enzymes structurally‐functionally characterized. The information obtained enabled us to improve the enzyme operational and catalytic properties by directed mutagenesis. However, the structural‐functional relationships explaining some desirable properties are not established yet and, therefore, their introduction was addressed by ‘non‐rational’ directed evolution. Then, over 100 oxidative biotransformations were analyzed. Among them, it is noteworthy to mention the regio/stereo selective hydroxylation of long/short‐chain alkanes (a chemically challenging reaction), epoxidation of alkenes, and production of hydroxy‐fatty acids. Concerning aromatic oxygenations, the regioselective hydroxylation of flavonoids, and stereoselective hydroxylation/epoxidation of alkyl/alkenyl‐benzenes were among the most remarkable reactions, together with enzymatic hydroxylation of benzene (as an alternative for harsh chemical process). Finally, peroxidases and peroxygenases also showed a potential as delignification biocatalysts and in the decolorization of contaminant dyes from textile industries. © 2014 The Authors. Biofuels, Bioproducts and Biorefi ning published by Society of Chemical Industry and John Wiley & Sons, Ltd.

Martínez, Angel T.↗

High purity Mn 5 O 8 nanoparticles with a high overpotential to gas evolution reactions for high voltage aqueous sodium-ion electrochemical storage

Here, developing electrodes with high specific energy by using inexpensive manganese oxides is of great importance for aqueous electrochemical energy storage (EES) using non-Li charge carriers such as Na-or K-ions. However, the energy density of aqueous EES devices is generally limited by their narrow thermodynamic potential window (~1.23 V). In this paper, the synthesis of high purity layered Mn 5 O 8 nanoparticles through solid state thermal treatment of Mn3O4 spinel nanoparticles, resulting in a chemical formula of [Mn 2 2+ ] [Mn 3 4+ O 8 2– ], evidenced by Rietveld refinement of synchrotron- based X-ray diffraction, has been reported. The electro-kinetic analyses obtained from cyclic voltammetry measurements in half-cells have demonstrated that Mn5O8 electrode has a large overpotential (~ 0.6 V) towards gas evolution reactions, resulting in a stable potential window of 2.5 V in an aqueous electrolyte in half-cell measurements. Symmetric full-cells fabricated using Mn5O8 electrodes can be operated within a stable 3.0 V potential window for 5000 galvanostatic cycles, exhibiting a stable electrode capacity of about 103 mAh/g at a C-rate of 95 with nearly 100% coulombic efficiency and 96% energy efficiency.

25 ENERGY STORAGE↗

In-situ exsolution of bimetallic CoFe nanoparticles on (La,Sr)FeO 3 perovskite: Its effect on electrocatalytic oxidative coupling of methane

This study presents in-situ reduction of lanthanum strontium cobalt ferrite (LSCF) perovskite as an effective method for modifying its surface properties and enhancing its electrocatalytic reactivity for oxidative coupling of methane (OCM). The evolution of hetero-phases during the reduction of LSCF resulted in CoFe nanoparticles being formed at the surface. The in-situ reduced LSCF cell for OCM could be operated in either an ion pump or a fuel cell mode. High selectivity of 63% and 10.2% were reported for C 2+ hydrocarbons and C 3 H 6 , respectively. DFT calculations on LSCF and CoFe revealed that the high selectivity of C 2+ hydrocarbons on the LSCF primarily stems from the presence of CoFe nanoparticles. In-situ DRIFTS conducted under CH 4 proved that complete oxidation of CH 4 can be effectively inhibited by reducing LSCF, and control of oxygen supply is an important parameter for selective conversion of CH 4 to higher order hydrocarbon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Surface Charge Predicts the Presence of Cation Effects in Electrocatalysis

The identity of electrolyte cations has an important influence on the rates of many electrocatalytic reactions, but these effects are not always observed. Recently, we found that the surface charge of the catalyst, quantified by the potential of zero total charge (PZTC), is a useful heuristic for predicting when cation effects will be observed for the oxygen reduction reaction (ORR). Here, we demonstrate that this descriptor allows us to rationalize the observation or absence of cation effects across a range of conditions, reactions (the hydrogen evolution reaction, the ORR, methanol oxidation, ethylene glycol oxidation, glycerol oxidation, and glucose reduction) and metal surfaces (Pt, Pd, Ag, and Au). These results suggest that when the reaction’s operating potential is negative of the metal’s PZTC, electrolyte cations accumulate at the catalyst surface and influence reaction rates. As a result, when reactions occur positive of the PZTC, cation effects are not observed.

Cations↗

Quantum logic gates from time-dependent global magnetic field in a system with constant exchange

We propose a method that implements a universal set of one- and two-quantum-bit gates for quantum computation in a system of coupled electron pairs with constant non-diagonal exchange interaction. In our proposal, suppression of the exchange interaction is performed by the continual repetition of single-spin rotations. A small g-factor difference between the electrons allows for addressing qubits and avoiding strong magnetic field pulses. Numerical experiments were performed to show that, to implement the one- and two-qubit operations, it is sufficient to change the strength of the magnetic field by a few Gauss. This introduces one and then the other electron in a resonance. To determine the evolution of the two-qubit system, we use the algorithms of optimal control theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Origin of enhanced water oxidation activity in an iridium single atom anchored on NiFe oxyhydroxide catalyst

The efficiency of the synthesis of renewable fuels and feedstocks from electrical sources is limited, at present, by the sluggish water oxidation reaction. Single-atom catalysts (SACs) with a controllable coordination environment and exceptional atom utilization efficiency open new paradigms toward designing high-performance water oxidation catalysts. Here, using operando X-ray absorption spectroscopy measurements with calculations of spectra and electrochemical activity, we demonstrate that the origin of water oxidation activity of IrNiFe SACs is the presence of highly oxidized Ir single atom (Ir 5.3+ ) in the NiFe oxyhydroxide under operating conditions. We show that the optimal water oxidation catalyst could be achieved by systematically increasing the oxidation state and modulating the coordination environment of the Ir active sites anchored atop the NiFe oxyhydroxide layers. Based on the proposed mechanism, we have successfully anchored Ir single-atom sites on NiFe oxyhydroxides (Ir 0.1 /Ni 9 Fe SAC) via a unique in situ cryogenic–photochemical reduction method that delivers an overpotential of 183 mV at 10 mA · cm – 2 and retains its performance following 100 h of operation in 1 M KOH electrolyte, outperforming the reported catalysts and the commercial IrO 2 catalysts. These findings open the avenue toward an atomic-level understanding of the oxygen evolution of catalytic centers under in operando conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermoeconomic cost optimization of superconducting magnets for proton therapy gantries

A compact gantry delivering 70-220 MeV protons with fixed field in the superconducting magnets could reduce the cost and improve the adoption of proton therapy. While a number of magnet and cryogenics designs have been proposed, the combined capital and operating costs of state-of-the-art superconducting materials have not been analyzed. In response, we develop a thermoeconomic model of a multi-stage, conduction cooled gantry lattice and analyze the cryocooler operating cost, cryocooler capital cost and conductor capital cost for Nb-Ti, Nb 3 Sn, REBCO and Bi-2223 over a continuous range of magnet temperatures, and a differential evolution algorithm is used to identify the optimal combination of thermal intercept temperatures. Although Nb3Sn yields the lowest Net Present Value (NPV) of $111.7k at a magnet temperature of 9.4 K, the optimized Bi-2223 design at 12.8 K approaches the realm of commercial feasibility by offering improved thermal stability and forgoing the need for costly conductor heat treatment and magnet quench training. Furthermore, it was found that Nb3Sn was more cost effective than Nb-Ti and that REBCO was not economically viable for the parameters of this investigation. Overall, the thermoeconomic model developed herein can optimize conductor choices, magnet temperatures and thermal staging which has value for any conduction-cooled superconducting magnet.

BSCCO↗