Search NASA⌕ Search

SEARCH · Search NASA

Results for “ARISE”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

High-Current Density Durability of Pt/C and PtCo/C Catalysts at Similar Particle Sizes in PEMFCs

The durability of carbon supported PtCo-alloy based nanoparticle catalysts play a key role in the longevity of proton-exchange membrane fuel cells (PEMFC) in electric vehicle applications. To improve its durability, it is important to understand and mitigate the various factors that cause PtCo-based cathode catalyst layers (CCL) to lose performance over time. These factors include i) electrochemical surface area (ECSA) loss, ii) specific activity loss, iii) H + /O 2 -transport changes and iv) Co 2+ contamination effects. We use a catalyst-specific accelerated stress test (AST) voltage cycling protocol to compare the durability of Pt and PtCo catalysts at similar average nanoparticle size and distribution. Our studies indicate that while Pt and PtCo nanoparticle catalysts suffer from similar magnitudes of electrochemical surface area (ECSA) losses, PtCo catalyst shows a significantly larger cell voltage loss at high current densities upon durability testing. The distinctive factor causing the large cell voltage loss of PtCo catalyst appears to be the secondary effects of the leached Co 2+ cations that contaminate the electrode ionomer. A 1D performance model has been used to quantify the cell voltage losses arising from various factors causing degradation of the membrane electrode assembly (MEA).

08 HYDROGEN↗
DOE OSTI2023

Visualization Quality Assessment

Understanding how inaccuracies in visualizations affect users’ perception and understanding of scientific data is hard. Inaccuracies in visualizations are quite common and could arise from a range of sources such as errors in the original dataset arising from compression artifacts, errors in the capturing device, noise during transmission of the data, effects due to the algorithm being used to convert data to visualization images, images generated from neural networks, and sources we have yet to discover. Many image quality assessment metrics have been developed to quantify image errors. However, these are usually focused on “natural images” rather than visualizations of scientific data. Common image quality assessment metrics (IQAs) include MSE, PSNR, perceptual metrics such SSIM, FSIM as well as perceptual metrics using deep learning approaches. However, a critical part of understanding how errors are perceived by humans, and subsequently developing more accurate quality assessment metrics, is through user evaluation studies. The goal of this software is to develop a visualization quality assessment (VQA) process that will enable the generation of VQAs that can be used to quantify errors in scientific data visualizations. The VQA development process will include software to support user evaluation experimental design, analysis of visualization differences against standard quality metrics, and the ability to develop additional VQA metrics specific to scientific visualization images.

Grosset, Andre↗
DOE OSTI2024

Earth System Reanalysis in Support of Climate Model Improvements

Recent climate model developments, established through increased model resolution, have led to substantial improvements in model simulations of the time-evolving, coupled Earth system and its subcomponents. However, regardless of resolution, climate models will always produce climate features and variability that differ from the real world and will be prone to biases. This is due to many remaining uncertainties, such as in parametric and structural model uncertainty, in the initial conditions prescribed, and in the prescribed (scenario) forcing which varies on decadal to centennial timescales. Further model improvements are expected to arise specifically from improved representation of physical processes realized through model-data fusion. This will create an unprecedented opportunity to better exploit a large array of Earth observations, from in situ measurements to weather radars and satellite observations, as the resolved scales of the models approach those of the observations. For this, climate DA will be the central tool to bring models and observations into consistency, by improving initial conditions, inferring uncertain model parameters and structure, and quantifying uncertainty. Generally, there will be advantages and complementarities of adjoint-based smoother approaches, ensemble-based filter approaches, or new ML-inspired approaches. Yet, the ever-increasing model resolution will present growing challenges arising from computational cost, calling for new ways of performing data assimilation and model optimization. Using the complementarity in a hybrid approach, blending tools and concepts from variational, ensemble and ML methods might be what is required in the future. In this context ML could be important to handle non-linear responses, and to better approximate non-Gaussian distributions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Theory for the Balance between Warm Rain and Ice Crystal Processes of Precipitation in Mixed-Phase Clouds

Abstract Mixed-phase clouds contain both supercooled cloud liquid and ice crystals. In principle, precipitation may be initiated either by the liquid phase or by the ice phase. Ice crystals may grow by vapor diffusion to become snow (“ice crystal process”), forming “cold” precipitation. Equally, cloud droplets, when large enough, coalesce to form “warm” precipitation by the “warm rain process.” Warm rain could be supercooled and freeze as “warm” graupel. In the present paper, a new simplified theoretical analysis is provided to examine the microphysical system consisting of three species of hydrometeor, namely, cloud liquid, “cold ice” (crystals, snow), and “warm rain” (frozen or supercooled). This is obtained by nondimensionalizing and simplifying the evolution equations for the mass of each species. Analytical formulas are given for equilibria. Feedback analysis shows that the sign of the feedback is linked to the abundance of precipitation, with a neutral surface in the 3D phase space. The system’s precipitation amount explodes while in the initial unstable regime, crossing the neutral surface and approaching the equilibrium point that is a stable attractor. Positive and negative feedbacks are elucidated. In a standard case, the cold ice mass is about 1000 times larger than the warm rain mass. To illustrate the physical behavior of the theory, sensitivity tests are performed with respect to environmental conditions (e.g., aerosol, updraft speed) and microphysical parameters (e.g., riming and sedimentation rates for cold ice). Cold ice prevails, especially in fast ascent, due to its low bulk density, favoring slow sedimentation and a wide cross-sectional area for riming. Significance Statement The theory elucidates how the ice phase can prevail in the precipitation from any mixed-phase clouds with supercooled cloud liquid and crystals. The ice phase radically suppresses cloud liquid by riming when active and “wins” the competition against coalescence. This prevalence of ice is shown to arise from the low bulk density of snow. The cloud is viewed as a system of negative and positive feedbacks that prevail in realms of stability and instability in a 3D phase space.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Seasonally and Regionally Dependent Shifts of the Atmospheric Westerly Jets under Global Warming

A distinct feature of the atmospheric circulation response to increasing greenhouse gas forcing is the poleward shift of the zonal-mean westerly jet. The dynamical mechanisms of the zonal-mean poleward jet shift have been extensively studied in literature. At seasonal/regional scales, however, the westerly jets can shift equatorward, such as in the early-summer Asia–Pacific region, the late-winter America–Atlantic region, and the winter/spring east Pacific. These equatorward jet shifts imply climate impacts distinct from those of the poleward shifts, yet their causes are not well understood. In this work, based on a hierarchy of coupled, prescribed-SST, and aquaplanet simulations, we attribute the seasonal/regional equatorward jet shifts to the enhanced tropical upper-level warming (ETUW), which arises from both the tropical moist adiabat and the enhanced equatorial surface warming. By steepening the meridional temperature gradient in the subtropical upper-to-middle level and assisted by positive eddy feedback, the ETUW increases the zonal wind equatorward of the climatological jet. When the regional/seasonal meridional temperature gradients (or equivalently the westerly jets) are weak and peak close to the tropics, the ETUW effect overcomes the poleward jet-shift mechanisms and leads to the equatorward jet shifts. This climatological-state dependency is consistently seen in the decomposed jet responses to uniform warming and surface warming pattern, and further demonstrated through idealized aquaplanet experiments with designed climatological states. For uniform warming, the ETUW arising from moist adiabat makes the general poleward jet shifts insignificant in the aforementioned favorable regions/seasons. For warming pattern, the ETUW from enhanced equatorial warming drives substantial equatorward jet shifts in these favorable seasons/regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Observational Constraints on the Cloud Feedback Pattern Effect

Abstract Model evidence for the “pattern effect” assumes that global climate models (GCMs) faithfully simulate how clouds respond to varying sea surface temperature (SST) patterns and associated meteorological perturbations. We exploit time-invariant satellite-based estimates of the sensitivity of marine low clouds to meteorological perturbations to estimate how these clouds responded to time-varying SST patterns and meteorology between 1870 and 2014. GCMs and reanalyses provide estimates of the historical meteorological changes. Observations suggest that increasing estimated inversion strength (EIS) between 1980 and 2014 produced a negative low cloud feedback, opposite to the positive feedback expected from increasing CO 2 . This indicates that the processes responsible for marine cloud changes from 1980 to the near present are distinct from those associated with an increase in CO 2 . We also observationally constrain the difference between the historical near-global marine low cloud feedback, , and that arising from increasing CO 2 , . We find that this cloud feedback pattern effect depends strongly on time period and reanalysis dataset, and that varying changes in EIS and SST with warming explain much of its variability. Between 1980 and 2014, we estimate that (90% confidence) assuming meteorological changes from the Multiple Reanalysis Ensemble, implying a total pattern effect (that arising from all climate feedbacks) of 1.86 ± 0.45 W m −2 K −1 . This observational evidence corroborates previous quantitative estimates of the pattern effect, which heretofore relied largely upon GCM-based cloud changes. However, disparate historical meteorological changes across individual reanalyses contribute to considerable uncertainty in its magnitude.

58 GEOSCIENCES↗
DOE OSTI2024

Decadal Thermal Variability of the Upper Southern Ocean: Zonal Asymmetry

Abstract As the major sink of anthropogenic heat, the Southern Ocean has shown quasi-symmetric, deep-reaching warming since the mid-twentieth century. In comparison, the shorter-term heat storage pattern of the Southern Ocean is more complex and has notable impacts on regional climate and marine ecosystems. By analyzing observational datasets and climate model simulations, this study reveals that the Southern Ocean exhibits prominent decadal (>8 years) variability extending to ∼700-m depth and is characterized by out-of-phase changes in the Pacific and Atlantic–Indian Ocean sectors. Changes in the Pacific sector are larger in magnitude than those in the Atlantic–Indian Ocean sectors and dominate the total heat storage of the Southern Ocean on decadal time scales. Instead of heat uptake through surface heat fluxes, these asymmetric variations arise primarily from wind-driven heat redistribution. Pacemaker and preindustrial simulations of the Community Earth System Model version 1 (CESM1) suggest that these variations in Southern Ocean winds arise primarily from natural variability of the tropical Pacific, as represented by the interdecadal Pacific oscillation (IPO). Through atmospheric teleconnection, the positive phase of the IPO gives rise to higher-than-normal sea level pressure and anticyclonic wind anomalies in the 50°–70°S band of the Pacific sector. These winds lead to warming of 0–700 m by driving the convergence of warm water. The opposite processes, involving cyclonic winds and upper-layer divergence, occur in the Atlantic–Indian Ocean sector. These findings aid our understanding of the time-varying heat storage of the Southern Ocean and provide useful implications on initialized decadal climate prediction.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Opposing Hysteresis Patterns in Flow and Outflow Macroscopic Fundamental Diagrams and Their Implications

Two key aggregated traffic models are the relationship between average network flow and density (known as the network or flow macroscopic fundamental diagram [flow-MFD]) and the relationship between trip completion and density (known as network exit function or the outflow-MFD [o-FMD]). The flow- and o-MFDs have been shown to be related by average network length and average trip distance under steady-state conditions. However, recent studies have demonstrated that these two relationships might have different patterns when traffic conditions are allowed to vary: the flow-MFD exhibits a clockwise hysteresis loop, while the o-MFD exhibits a counter-clockwise loop. One recent study attributes this behavior to the presence of bottlenecks within the network. The present paper demonstrates that this phenomenon may arise even without bottlenecks present and offers an alternative, but more general, explanation for these findings: a vehicle’s entire trip contributes to a network’s average flow, while only its end contributes to the trip completion rate. This lag can also be exaggerated by trips with different lengths, and it can lead to other patterns in the o-MFD such as figure-eight patterns. A simple arterial example is used to demonstrate this explanation and reveal the expected patterns, and they are also identified in real networks using empirical data. Then, simulations of a congestible ring network are used to unveil features that might increase or diminish the differences between the flow- and o-MFDs. Finally, more realistic simulations are used to confirm that these behaviors arise in real networks.

42 ENGINEERING↗
DOE OSTI2023

Fast truncated SVD of sparse and dense matrices on graphics processors

We investigate the solution of low-rank matrix approximation problems using the truncated singular value decomposition (SVD). For this purpose, we develop and optimize graphics processing unit (GPU) implementations for the randomized SVD and a blocked variant of the Lanczos approach. Our work takes advantage of the fact that the two methods are composed of very similar linear algebra building blocks, which can be assembled using numerical kernels from existing high-performance linear algebra libraries. Furthermore, the experiments with several sparse matrices arising in representative real-world applications and synthetic dense test matrices reveal a performance advantage of the block Lanczos algorithm when targeting the same approximation accuracy.

Computer Science↗
DOE OSTI2021

Orbital Conflict: Cutting Planes for Symmetric Integer Programs

Cutting planes have been an important factor in the impressive progress made by integer programming (IP) solvers in the past two decades. However, cutting planes have had little impact on improving performance for symmetric IPs. Rather, the main breakthroughs for solving symmetric IPs have been achieved by cleverly exploiting symmetry in the enumeration phase of branch and bound. In this work, we introduce a hierarchy of cutting planes that arise from a reinterpretation of symmetry-exploiting branching methods. There are too many inequalities in the hierarchy to be used efficiently in a direct manner. However, the lowest levels of this cutting-plane hierarchy can be implicitly exploited by enhancing the conflict graph of the integer programming instance and by generating inequalities such as clique cuts valid for the stable set relaxation of the instance. We provide computational evidence that the resulting symmetry-powered clique cuts can improve state-of-the-art symmetry-exploiting methods. Furthermore, the inequalities are then employed in a two-phase approach with high-throughput computations to solve heretofore unsolved symmetric integer programs arising from covering designs, establishing for the first time the covering radii of two binary-ternary codes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Data for An End-to-End Pipeline for Succinic Acid Production at an Industrially Relevant Scale Using Issatchenkia orientalis

Microbial production of succinic acid (SA) at an industrially relevant scale has been hindered by high downstream processing costs arising from neutral pH fermentation for over three decades. Here, we metabolically engineer the acid-tolerant yeast Issatchenkia orientalis for SA production, attaining the highest titers in sugar-based media at low pH (pH 3) in fed-batch fermentations, i.e. 109.5 g/L in minimal medium and 104.6 g/L in sugarcane juice medium. We further perform batch fermentation using sugarcane juice medium in a pilot-scale fermenter (300×) and achieve 63.1 g/L of SA, which can be directly crystallized with a yield of 64.0%. Finally, we simulate an end-to-end low-pH SA production pipeline, and techno-economic analysis and life cycle assessment indicate our process is financially viable and can reduce greenhouse gas emissions by 34–90% relative to fossil-based production processes. We expect I. orientalis can serve as a general industrial platform for production of organic acids.

Metabolomics↗
DOE OSTI2020

Data from Improved Annotation of Untargeted Metabolomics Data through Buffer Modifications That Shift Adduct Mass and Intensity

Annotation of untargeted high-resolution full-scan LC-MS metabolomics data remains challenging due to individual metabolites generating multiple LC-MS peaks arising from isotopes, adducts, and fragments. Adduct annotation is a particular challenge, as the same mass difference between peaks can arise from adduct formation, fragmentation, or different biological species. To address this, here we describe a buffer modification workflow (BMW) in which the same sample is run by LC-MS in both liquid chromatography solvent with 14NH3–acetate buffer and in solvent with the buffer modified with 15NH3–formate. Buffer switching results in characteristic mass and signal intensity changes for adduct peaks, facilitating their annotation. This relatively simple and convenient chromatography modification annotated yeast metabolomics data with similar effectiveness to growing the yeast in isotope-labeled media. Application to mouse liver data annotated both known metabolite and known adduct peaks with 95% accuracy. Overall, it identified 26% of ∼27 000 liver LC-MS features as putative metabolites, of which ∼2600 showed HMDB or KEGG database formula match. This workflow is well suited to biological samples that cannot be readily isotope labeled, including plants, mammalian tissues, and tumors.

Conversion↗
DOE OSTI2024

Statistical mechanics and pressure of composite multimoded weakly nonlinear optical systems

Statistical mechanics can provide a versatile theoretical framework for investigating the collective dynamics of weakly nonlinear-wave settings that can be utterly complex to describe otherwise. In optics, composite systems arise due to interactions between different frequencies and polarizations. The purpose of this work is to develop a thermodynamic theory that takes into account the synergistic action of multiple components. We find that the type of the nonlinearity involved can have important implications in the thermalization process and, hence, can lead to different thermal equilibrium conditions. Importantly, we derive closed-form expressions for the actual optomechanical pressure that is exerted on the system. In particular, the total optomechanical pressure is the sum of the partial pressures due to each component. Our results can be applied to a variety of weakly nonlinear optical settings such as multimode fibers, bulk waveguides, photonic lattices, and coupled microresonators. We present two specific examples, where two colors interact in a one-waveguide array with either a cubic or quadratic nonlinearity.

Efremidis, Nikolaos K. (ORCID:0000000298300268)↗
DOE OSTI2022

Analysis, identification and confirmation of synthetic opioids using chloroformate chemistry: Retrospective detection of fentanyl and acetylfentanyl in urine and plasma samples by EI-GC-MS and HR-LC-MS

Electron Impact Gas Chromatography-Mass Spectrometry (EI-GC-MS) and High Resolution Liquid Chromatography-Mass Spectrometry (HR-LC-MS) have been used in the analysis of products arising from the trichloroethoxycarbonylation of fentanyl and acetylfentanyl in urine and plasma matrices. The method involves the initial extraction of both synthetic opioids separately from the matrices followed by detection of the unique products that arise from their reaction with 2,2,2-trichloroethoxycarbonyl chloride (Troc-Cl), namely Troc-norfentanyl and Troc-noracetylfentanyl. The optimized protocol was successfully evaluated for its efficacy at detecting these species formed from fentanyl and acetylfentanyl when present at low and high levels in urine (fentanyl: 5 and 10 ng/mL and acetylfentanyl: 20 and 100 ng/mL) and plasma (fentanyl: 10 and 20 ng/mL and acetylfentanyl: 50 and 200 ng/mL), values that reflect levels reported in overdose victims. The HR-LC-MS method’s LOQ (limit of quantitation) for the Troc-norfentanyl and Troc-noracetylfentanyl products was determined to be ~10 ng/mL for both species. Even though the superiority in the detection of these species by HR-LC-MS over EI-GC-MS, the latter method proved to be important in the detection of the second product from the reaction, namely 2-phenylethyl chloride that is crucial in the determination of the original opioid. This observation highlights the importance of using complimentary analytical techniques in the analysis of a sample, whether biological or environmental in nature. The method herein serves as a complementary, qualitative confirmation for the presence of a fentanyl in collected urine, plasma and by extension other biological samples amenable to the common extraction procedures described for opioid analysis. More importantly, the method’s main strength comes from its ability to react with unknown fentanyls to yield products that can be not only detected by EI-GC-MS and HR-LC-MS but can then be used to retrospectively identify an unknown fentanyl.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

MFNets: Multifidelity data-driven networks for Bayesian learning and prediction

This paper presents a multifidelity uncertainty quantification framework called MFNets. We seek to address three existing challenges that arise when experimental and simulation data from different sources are used to enhance statistical estimation and prediction with quantified uncertainty. Specifically, we demonstrate that MFNets can (1) fuse heterogeneous data sources arising from simulations with different parameterizations, e.g simulation models with different uncertain parameters or data sets collected under different environmental conditions; (2) encode known relationships among data sources to reduce data requirements; and (3) improve the robustness of existing multi-fidelity approaches to corrupted data. MFNets construct a network of latent variables (LVs) to facilitate the fusion of data from an ensemble of sources of varying credibility and cost. These LVs are posited as explanatory variables that provide the source of correlation in the observed data. Furthermore, MFNets provide a way to encode prior physical knowledge to enable efficient estimation of statistics and/or construction of surrogates via conditional independence relations on the LVs. We highlight the utility of our framework with a number of theoretical results which assess the quality of the posterior mean as a frequentist estimator and compare it to standard sampling approaches that use single fidelity, multilevel, and control variate Monte Carlo estimators. We also use the proposed framework to derive the Monte Carlo-based control variate estimator entirely from the use of Bayes rule and linear-Gaussian models -- to our knowledge the first such derivation. Finally, we demonstrate the ability to work with different uncertain parameters across different models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

How nitrogen and oxygen shape SRF cavity performance

Nitrogen and oxygen-based surface treatments have revolutionized the performance of superconducting radiofrequency (SRF) cavities, enabling them to reach higher gradients and lower losses. However, the exact mechanisms by which these treatments improve cavity performance remain largely unknown. This work provides new insights into the role of nitrogen and oxygen in SRF cavity performance by using time-of-flight secondary ion mass spectrometry (TOF-SIMS) to precisely quantify the concentrations and depth profiles of these impurities within niobium cutouts. We correlate these impurity profiles with detailed cavity performance measurements, including surface resistance and quality factor, and compare our findings with predictions from BCS theory. The results demonstrate that while both nitrogen and oxygen enhance performance, ten times more oxygen is required to achieve the same reduction in BCS resistance as interstitial nitrogen. We present a potential model in which the observed variation arises from nitrogen's greater effectiveness in trapping hydrogen, thus reducing the formation of niobium hydrides and enhancing superconducting gap.

Hu, Hannah [Chicago U.]↗
DOE OSTI2019

Enhancement of EM Pump Performance through Modeling and Testing

Liquid sodium has seen development as a nuclear reactor coolant in the United States for over 50 years. While sodium has many desirable characteristics as a heat transfer fluid, its reactivity poses numerous engineering challenges for even the most basic plant components like pumps. Despite the challenges, mechanical pumps have seen success in sodium service albeit at the expense of simplicity and flexibility. In fact, significant failures of these pumps have been recorded which can be attributed to their complexity. Notwithstanding the complexity, mechanical pumps are still widely proposed for the main circuit pump of sodium cooled reactors. However, an alternative to mechanical pumps are Electromagnetic Induction Pumps (EMIPs) which use the principle of the Lorentz Force. EMIPs avoid mechanically moving parts in contact with sodium which simplifies the design by avoiding mechanical wear and seals. Therefore, these pumps are generally safer and more reliable than their mechanical counterparts. That being said, EMIPs suffer from relatively low efficiencies that range from 5% to 45%. Some of the inefficiencies are inherent to the EMIPs pumps; namely, inefficiencies arise in the form of resistive heating losses in the fluid itself. However, other inefficiencies arise from Edge Effects which include Finite Length and Finite Width effects. Due to the significant safety advantages EMIPs have compared to mechanical pumps, there is a strong motivation to integrate these technologies into the design of Generation IV liquid-metal-cooled nuclear power plants. However, EMIPs have low efficiencies partly attributed to Edge Effects. Much of the literature focusing on Edge Effects is dated or simply does not exist in the public domain. More recently, numerical models of EMIPs have shown promising results. However, these models are often custom solvers or use complex software suits which limit their application. This report describes the work completed to understand Edge Effects in two EMIPs designs. These designs are the Annular Linear Induction Pump (ALIP) and Permanent Magnet Induction Pump (PMIP) also known as a Moving Magnet Pump (MMP). The methodology developed during this project aim at measuring key performance parameters (pressure-flowrate and efficiency-slip performance) using high temperature liquid-sodium. Baseline experimental measurements were then used to verify the accuracy of numerical models developed in the software suit ANSYS/FLUENT. On the experimental side, instrumentation was developed to characterize the pressure-flowrate and efficiency-slip performance of both the ALIP and PMIP in high-temperature liquid-sodium. Additionally, instrumentation was developed to measure potential output pressure pulsations and the magnetic field distributions. Moreover, a specialized ALIP was constructed to study the impacts magnetic field tapering, also known as coil grading, on the Finite Length Effect. Lastly, experimental performance curves of the PMIP were collected at 300 oC and 600 oC and compared to analytical theory which accounted for Finite Width Effects. These results showed that the Finite Width Effect and non-uniform magnetic field in a PMIP impact the performance. On the numerical side, models were developed using the software suit FLUENT/ANSYS. These methods were initially compared against a simple MHD flow in a Permanent Magnet Electromagnet Flowmeter (PMEMFM). The results were significant in that they showed agreeance between the numerical model and demonstrated that PMEMFM’s could be calibrated numerically. These modeling methods were then applied to a prototypic PMIP. Results showed that this simple model provided reasonably accurate results of the PMIP p-Q performance.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Characteristics of Ultralow Aspect Ratio Toroidal Plasmas and Non-Solenoidal Startup and Edge Stability Studies at Near-Unity Aspect Ratio (Final Report)

This grant supported a long-term research and development program to explore the physics and technology of spherical tokamak (ST) plasma behavior at near-unity aspect ratio ( A ). The PEGASUS Toroidal Experiment is a university-scale ( R ) ~ 0.40 m, α ~ .35 m, I p ≤ 0.3 MA, B T ≤ 0.15 T) experiment that accesses A → 1, to provide large values of the normalized plasma current, I N = I p /α B T , or equivalently the toroidal field utilization factor I p / I TF . The research activities focused on: 1) developing the technique of Local Helicity Injection (LHI) as a viable means of initiating and growing an ST plasma without the use of a central solenoid; 2) exploring the equilibrium and stability properties of spherical tokamak plasmas at extremely low aspect ratio; and 3) obtaining new insights in the properties of the enhanced confinement H-mode regime at A ~ 1. The initial (Phase I) facility employed a unique high-stress central solenoid Ohmic for current drive and heating with simple magnet current waveforms. A subsequent major upgrade (Phase II) offered more flexible control of the plasma evolution along with a powerful LHI system for plasma initiation, heating and sustainment. The LHI technique consists of injection of high-density current streams into vacuum toroidal and vertical magnetic fields. These streams become unstable and support magnetic reconnection to allow the system to relax into a tokamak-like magnetic configuration. This program focused on: 1) development of electron current sources that provide sufficiently high injection currents in the challenging edge plasma region; 2) validation of theoretical limits to achievable I p arising from magnetic helicity and energy conservation plus absolute limits arising from relaxation to a minimal energy states; 3) tests of a power balance model to describe the evolution of the plasma current I p (t); and 4) study of magnetic fluctuations and their role in the I p evolution. This culminated in the production of toroidal currents of I p ≥ 0.2 MA with only ~8 kA of injected current. Operation with Ohmic induction alone in Phase I created plasmas with soft limits on the achievable I p and field utilization ( I p / I TF ~ 1) due to the onset of large low-order (2/1 or 3/2) internal resonant tearing modes. With the expanded capabilities in Phase II operation, methods employed to mitigate those modes and achieve I p / I TF > 1 were: 1) I TF ramp-downs; 2) non-solenoidal startup via LHI; and 3) extremely low I TF Ohmic operations assisted by LHI current injector pre-ionization. Extensive use of LHI-driven startup and drive for non-solenoidal operation produced favorable current profiles to stably operate tokamak plasmas with unprecedentedly high I N = I p / aB T ~ 15 and I p / I TF > 2. With strong reconnection-driven anomalous ion heating, this in turn provided access to tokamak plasmas with B t up to 100%, a world record 2 times higher than that previously achieved in advanced or spherical tokamaks. This high B t regime had a large fraction (~50%) of the plasma volume in a region of an absolute minimum- B configuration. The MHD stability limit in this regime was the ideal external kink mode without a stabilizing wall. These results provided the first demonstration of access to tokamak plasmas with very high B t , high i N , high elongation κ, and low internal inductance ℓ i , and as such served as a proof-ofprinciple for access to this long-sought ST regime. The naturally low B T needed for stability at A ≤ 1.2 reduces the L-H mode power threshold, P LH , so that H-mode plasmas were achieved with Ohmic heating alone. The power threshold for the transition to the H-mode regime was found to be 10× the accepted values from the ITER international scaling, confirming initial trends indicated by experiments on MAST and NSTX, and suggesting that the physics behind this transition is not yet fully understood. The density dependence of the transition and the observed insensitivity to whether the plasma is bounded by a limiter or magnetic divertor are consistent with the recent FM 3 model. Multiple unstable toroidal magnetic modes were measured during ELM (Edge Localized Mode) crashes, with generally lower toroidal mode numbers than seen at higher A , which was attributed to the larger peeling mode instability drive at low A . First-ever edge current profile measurements through an ELM crash showed current-carrying filaments were generated, as seen in ELM simulations and in electromagnetic blob transport theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗