SEARCH · Engineering Papers
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.
When a fusion plasma suddenly intercepts a solid surface (wall or pellet), thermal collapse is distinctly kinetic & has novel physics. VPIC simulations & theory revealed the fascinating dynamics of four propagating fronts controlling parallel electron temperature cooling.
Sheared flow increases the reactivity of fusion plasma. In unmagnetized plasma with flow gradients comparable to the mean free path of reacting ions, fusion reactivity can be more than doubled. The effect is of particular relevance to inertial confinement fusion (ICF), where it allows implosion kinetic energy to contribute to the fusion burn even before thermalizing. Finally, because colder fuel stops alpha particles more quickly, ignition is possible in a smaller volume, substantially reducing energy requirements in fast-ignition designs.
This grant facilitated the APS-DPP Community Planning Process (DPP-CPP) which was a strategic planning process with the goals of identifying scientific and technological opportunities in the fields of Plasma Physics and Fusion Energy Science, and making consensus recommendations for a strategy to address these opportunities to the DOE Office of Fusion Energy Science Advisory Committee (FESAC). The DPP-CPP was initiated by the APS-DPP executive committee in order to help FESAC fulfill a charge from DOE to develop a strategic plan for the DOE Office of Fusion Energy Science (FES). In this charge, the DOE Deputy Director for Science, Dr. Stephen Binkley, requests “that the Fusion Energy Sciences Advisory Committee (FESAC) undertake a new long-range strategic planning activity for the Fusion Energy Sciences (FES) program. The strategic planning activity - to encompass the entire FES research portfolio (namely, burning plasma science and discovery plasma science) - should identify and prioritize the research required to advance both the scientific foundation needed to develop a fusion energy source, as well as the broader FES mission to steward plasma science.” The process as set out in the proposal (and below) was successful in meeting all of its stated goals culminating in the final report to the FESAC subcommittee given the charge.
We report that in steady state, the fuel cycle of a fusion plasma requires inward particle fluxes of fuel ions. These particle flows are also accompanied by heating. In the case of classical transport in a rotating cylindrical plasma, this heating can proceed through several distinct channels depending on the physical mechanisms involved. Some channels directly heat the fuel ions themselves, whereas others heat electrons. Which channel dominates depends, in general, on the details of the temperature, density and rotation profiles of the plasma constituents. However, remarkably, under relatively few assumptions concerning these profiles, if the α particles, the by-products of the fusion reaction, can be removed directly by other means, then a hot-ion mode tends to emerge naturally.
Toroidal current profile measurements in the tokamak plasma edge are critical for fusion plasma physics research and model validation. A three-wave Faraday-effect polarimeter-interferometer with a sub-centimeter spatial resolution is proposed on the DIII-D tokamak to determine the edge current profile via Abel inversion. By using probe beams with 316 µm wavelength, a low-field-side, vertical-view, single-pass optical layout covering the plasma edge region (R = 2.15–2.27 m) is assessed. Measurements with no greater than 0.1° polarimetric systematic uncertainty, no greater than 0.01° polarimetric root-mean-square noise (1 kHz bandwidth), and a 0.8 cm radial chord spacing are considered feasible based on the achieved performance of existing systems using similar wavelengths on fusion devices. Furthermore, synthetic diagnostic calculations taking various factors into account, such as diagnostic uncertainty and quality of magnetic flux surfaces, find that the edge current profile can be determined with up to 0.12 MA/m 2 uncertainty, or about 10% of the peak current density in the pedestal of an investigated high-confinement plasma.
Advanced operation scenarios such as high poloidal beta ( β P ) or high q min are promising concepts to achieve the steady-state high-performance fusion plasmas. However, those scenarios are prone to substantial Alfvénic activity, causing fast-ion transport and losses. Recent experiments with the advanced operation scenario on KSTAR tokamak have shown that the electron cyclotron current drive (ECCD) is able to mitigate and suppress the beam-ion driven toroidal Alfvén eigenmodes (TAEs) for over several tens of global energy confinement time. Co-current directional intermediate off-axis ECCD lowers the central safety factor slightly and tilts the central q -profile shape so that the continuum damping in the core region increases. Besides, the rise of central plasma pressure and increased thermal-ion Landau damping contribute to TAE stabilization. While the TAEs are suppressed, neutron emission rate and total stored energy increase by approximately 45% and 25%, respectively. Fast-ion transport estimated by TRANSP calculations approaches the classical level during the TAE suppression period. Substantial reduction in fast-ion loss and neutron deficit is also observed. Subsequently, enhancement of fast-ion confinement by suppressing the TAEs leads to an increase of non-inductive current fraction and will benefit the sustainment of the long-pulse high-performance discharges.
In the dynamic environment of burning, thermonuclear deuterium–tritium plasmas, diagnosing the time-resolved neutron energy spectrum is of critical importance. Strategies exist for this diagnosis in magnetic confinement fusion plasmas, which presently have a lifetime of ~10 12 longer than inertial confinement fusion (ICF) plasmas. In this work, we present a novel concept for a simple, precise, and scale-able diagnostic to measure time-resolved neutron spectra in ICF plasmas. The concept leverages general tomographic reconstruction techniques adapted to time-of-flight parameter space, and then employs an updated Monte Carlo algorithm and National Ignition Facility-relevant constraints to reconstruct the time-evolving neutron energy spectrum. Reconstructed spectra of the primary 14.028 MeV n DT peak are in good agreement with the exact synthetic spectra. The technique is also used to reconstruct the time-evolving downscattered spectrum, although the present implementation shows significantly more error.
Magnetized target fusion (MTF) attempts to combine the favorable attributes of magnetic confinement fusion (MCF) for energy confinement with the attributes of inertial confinement fusion (ICF) for efficient compression heating and wall-free containment of the fusing plasma. It uses a material liner to compress and contain a magnetized plasma. For practical applications, standoff drivers to deliver the imploding momentum flux to the target plasma remotely are required. Spherically converging plasma jets have been proposed as standoff drivers for this purpose. The concept involves the dynamic formation of a spherical plasma liner by the merging of plasma jets, and the use of the liner so formed to compress a spheromak or a field reversed configuration (FRC).
Magnetized target fusion (MTF) attempts to combine the favorable attributes of magnetic confinement fusion (MCF) for energy confinement with the attributes of inertial confinement fusion (ICF) for efficient compression heating and wall-free containment of the fusing plasma. It uses a material liner to compress and contain a magnetized plasma. For practical applications, standoff drivers to deliver the imploding momentum flux to the target plasma remotely are required. Spherically converging plasma jets have been proposed as standoff drivers for this purpose. The concept involves the dynamic formation of a spherical plasma liner by the merging of plasma jets, and the use of the liner so formed to compress a spheromak or a field reversed configuration (FRC). For the successful implementation of the scheme, plasma jets of the requisite momentum flux density need to be produced. Their transport over sufficiently large distances (a few meters) needs to be assured. When they collide and merge into a liner, relative differences in velocity, density and temperature of the jets could give rise to instabilities in the development of the liner. Variation in the jet properties must be controlled to ensure that the growth rate of the instabilities are not significant over the time scale of the liner formation before engaging with the target plasma. On impact with the target plasma, some plasma interpenetration might occur between the liner and the target. The operating parameter space needs to be identified to ensure that a reasonably robust and conducting contact surface is formed between the liner and the target. A mismatch in the "impedance" between the liner and the target plasma could give rise to undesirable shock heating of the liner leading to increased entropy (thermal losses) in the liner. Any irregularities in the liner will accentuate the Rayleigh-Taylor instabilities during the compression of the target plasma by the liner.
Charge transfer at electron-volt energies between multiply charged atomic ions and neutral atoms and molecules is of considerable importance in astrophysics, plasma physics, and in particular, fusion plasmas. In the year covered by this report, several major tasks were completed. These include: (1) the re-calibration of the ion gauge to measure the absolute particle densities of H2, He, N2, and CO for our current measurements; (2) the analysis of data for charge transfer reactions of N(exp 2 plus) ion and He, H2, N2, and CO; (3) measurement and data analysis of the charge transfer reaction of (Fe(exp 2 plus) ion and H2; (4) charge transfer measurement of Fe(exp 2 plus) ion and H2; and (5) redesign and modification of the ion detection and data acquisition system for the low energy beam facility (reflection time of flight mass spectrometer) dedicated to the study of state select charge transfer.
One goal of the Numerical Turbulent Transport Project is to model a tokamak (fusion) plasma with 10(sup 8 - 10)(sup 9) particles, to explain anomalous transport of particles and energy.
Extracting reliable information from diagnostic data in tokamaks is critical for understanding, analyzing, and controlling the behavior of fusion plasmas and validating models describing that behavior. Recent interest within the fusion community has focused on the use of principled statistical methods, such as Gaussian process regression (GPR), to attempt to develop sharper, more reliable, and more rigorous tools for examining the complex observed behavior in these systems. While GPR is an enormously powerful tool, there is also the danger of drawing fragile, or inconsistent conclusions from naive GPR fits that are not driven by principled treatments. Here we review the fundamental concepts underlying GPR in a way that may be useful for broad-ranging applications in fusion science. We also revisit how GPR is developed for profile fitting in tokamaks. We examine various extensions and targeted modifications applicable to experimental observations in the edge of the DIII-D tokamak. Finally, we discuss best practices for applying GPR to fusion data.
Small-amplitude, symmetry-breaking magnetic field perturbations, including resonant magnetic perturbations (RMPs) and error fields, can profoundly impact plasma properties in both tokamaks and stellarators. In this work, we perform the first comparison between the Stepped Pressure Equilibrium Code (SPEC) (a comparatively fast and efficient equilibrium code based on energy-minimisation principles) and M3D-C 1 (a high-fidelity albeit computationally expensive initial-value extended-magnetohydro-dynamic (MHD) code) to assess the conditions under which SPEC can be used to model the nonlinear, non-ideal plasma response to an externally applied (m = 2,n = 1) RMP field in an experimentally relevant geometry. We find that SPEC is able to capture the plasma response in the weakly nonlinear regime – meaning perturbation amplitudes below the threshold for break up of the separatrix and onset of secondary magnetic island formation – when around half of the total toroidal flux is enclosed in the volume containing the q = 2 resonant surface. Further, the observed dependence of SPEC solutions on input parameters, including toroidal flux and the number of volumes into which the plasma is partitioned, indicates that additional exploration of the underlying Multi-Region Relaxed MHD physics model is needed to constrain the choice of parameters. Nonetheless, this work suggests promising applications of SPEC to optimisation and fusion plasma design.
Fusion is an attractive clean-energy solution, thanks to its various advantages, such as reduced radioactivity, little high-level nuclear waste, ample fuel supplies, and increased safety. However, the harsh operating environment introduced by a complex fusion plasma system makes design and integration of fusion blankets incredibly challenging and time-consuming. This work focuses on developing a fully integrated multiphysics simulation framework based on an advanced open-source platform—the Multiphysics Object-Oriented Simulation Environment (MOOSE)—to alleviate the difficulties in fusion blanket design and integration. MOOSE is a massively parallel finite element/volume multiphysics simulation platform that has been widely adopted within the nuclear fission community. Even though fission and fusion are fundamentally different, they involve similar multiphysics phenomena. A fully integrated open-source multiphysics simulation framework tailored for the fusion blanket design will be implemented by leveraging the well-established multiphysics capabilities in MOOSE. Once successfully developed, this fully integrated framework will rapidly evaluate a blanket design concept and offer insights for subsequent iterations. As the first step, we will mainly aim to integrate neutronics analysis, system thermal hydraulics simulation, and full 3-D heat transfer calculations. The efficacy of the integrated framework will be verified using an innovative solid ceramic blanket design. While the project’s final goal is to enable a fully integrated multiphysics simulation platform for various fusion blanket concepts, here this work, as a preliminary step, will mainly focus on a solid ceramic breeder helium-cooled blanket.
We report the ITER Diagnostic Residual Gas Analyzer (DRGA) will measure the distribution of gas species, i.e., deuterium (D), tritium (T), and impurities, in the divertor exhaust stream and in the plasma periphery, with time resolution relevant to fusion plasma–wall particle dynamics. The uniqueness of the DRGA, over previous implementations of plasma dynamics residual gas analysis, is an integrated approach, combining mass and low-temperature plasma-activated optical spectroscopy, in a differentially pumped analysis station. A further unique feature of the ITER divertor-specific DRGA is an ~8-m separation of the analysis station from the sampled pumping duct, while still maintaining a ~1-s response time for hydrogen isotopic concentrations. ITER DRGA final design activities are strongly benefiting from testing of prototypical DRGA components and methods on present fusion devices, most currently on JET and W7-X. DRGA systems are implemented on both these devices and include sensors (and pumping methods) that are directly relevant to the ITER DRGA design. The recent JET-DTE2 campaign has provided the first experience on operating the combined ITER DRGA sensors with D-T plasmas. While enhancing system design for ITER, this experience on operating devices has also revealed additional engineering challenges, which further guide the continuing final design project. Meanwhile, the recent determination that the ITER DRGA, with slight optimization, will resolve the helium isotopes well enough to support an ITER pre-DT, He-3-based heating scheme, has greatly increased ITER Research Program interest in the DRGA and its implementation well ahead of the DT phase.
Tokamak fusion plasmas benefit from high pressures but are then susceptible to modes of instability. These magnetohydrodynamic (MHD) modes are macroscopic distortions of the plasma, but certain collective motions of individual particles can provide stabilizing effects opposing them. The presence of a resistive wall slows the mode growth, converting a kink to a resistive wall mode (RWM). A kinetic MHD model includes Maxwell's equations, ideal MHD constraints, and kinetic effects included through the pressure tensor, calculated with the perturbed drift-kinetic distribution function of the particles. The kinetic stabilizing effects on the RWM arise through resonances between the plasma rotation and particle drift motions: precession, bounce, and transit. A match between particle motions and the mode allows efficient transfer of energy that would otherwise drive the growth of the mode, thus damping the growth. The first approach to calculating RWM stability is to write a set of equations for the complex mode frequency in terms of known quantities and then to solve the system. The “energy principle” approach, which has the advantage of clarity in distinguishing the various stabilizing and destabilizing effects, is to change the force balance equation into an equation in terms of changes of kinetic and potential energies, and then to write a dispersion relation for the mode frequency in terms of those quantities. These methods have been used in various benchmarked codes to calculate kinetic effects on RWM stability. Importantly, the theory has illuminated the important roles of plasma rotation, energetic particles, and collisions in RWM stability.
This is the Final Progress Report of the Helicity Injected Torus (HIT) program covering May 2016 through December 2020. During this final period the HIT program continued to make advancements in computational modeling, theory, and experimental studies of global self-organization and inductively sustained spheromaks. This included many valuable contributions to the fusion, plasma and physics community, through peer-reviewed publications, PhD theses and conference presentations and posters.
The parasitic interaction between ion cyclotron range of frequencies (ICRF) waves and the edge plasma has long limited high-power ICRF utilization in magnetically confined fusion plasmas. In this paper, we report on the successful mitigation of this interaction by power tapering a four-strap field-aligned (FA) antenna in the Alcator C-Mod tokamak. With the antenna operated in dipole phasing ($0\pi0\pi$), the ratio of the power coupled by the central two straps ($P_\mathrm {cent}$) to the power coupled by the outer two straps was varied at fixed total coupled ICRF power ($P_\mathrm {tot}$). With $P_\mathrm {cent}$/$P_\mathrm {tot} \approx 0.8-0.9$, no enhancement of the plasma potential was measured despite 1 MW of coupled ICRF power. When the antenna is powered, energy deposition increases on its four corners and on neighboring limiter tiles magnetically connected to some of those corners. The energy deposited on the antenna and limiter is minimized together with near-field plasma potentials. Additionally, we show that deviation from $P_\mathrm {cent}$/$P_\mathrm {tot} \approx 0.8$ enhances wave-plasma interactions in places not magnetically connected to the antenna. These far-field effects are correlated with changes in the antenna spectrum. Hence, on C-Mod, near-field and far-field interactions were minimized for about the same range of $P_\mathrm {cent}$/$P_\mathrm {tot}$. Optimal power tapering was also found to improve ICRF heating efficiency and facilitate H-mode access. Finally, we show that the heat flux pattern on the toroidally aligned antennas reverses with the direction of the RF-induced E×B flow in front of the antenna, a correlation not observed for the FA antenna. We argue that the RF-induced E×B flow, which carries high-density plasma from the near-SOL, is aligned with the FA antenna and does not intercept its corners, where the plasma potential is large. This could explain the previously reported smaller impurity generation by the FA antenna (Wukitch et al 2013 Phys. Plasmas 20 056117).