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.
In this project, we have further developed GTC gyrokinetic particle simulations to study plasma confinement properties in the stellarators and tokamaks with 3D magnetic fields, where both neoclassical and turbulent transport are self-consistently treated on an equal footing.
Multimachine empirical scaling predicts an extremely narrow heat exhaust layer in future high magnetic field tokamaks, producing high power densities that require mitigation. Here, in the experiments presented, the width of this exhaust layer is nearly doubled using actuators to increase turbulent transport in the plasma edge. This is achieved in low collisionality, high confinement edge pedestals with their gradients limited by turbulent transport instead of large-scale, coherent instabilities. The exhaust heat flux profile width and divertor leg diffusive spreading both double as a high frequency band of turbulent fluctuations propagating in the electron diamagnetic direction doubles in amplitude. The results are quantitatively reproduced in electromagnetic XGC particle-in-cell simulations which show the heat flux carried by electrons emerges to broaden the heat flux profile, directly supported by Langmuir probe measurements.
A steady state solution for the self-consistent electrostatic potential due to a plasma confined in a magnetic flux tube is considered. A steady state distribution function is constructed for the trapped particles from the constants of the motion, in the absence of waves and collisions. Using Liouville's theorem, the particle density along the geomagnetic field is determined and found to depend on the local magnetic field, self-consistent electric potential, and the equatorial plasma distribution function. A hot anisotropic magnetospheric plasma in steady state is modeled by a bi-Maxwellian at the equator. The self-consistent electric potential along the magnetic field is calculated assuming quasineutrality, and the potential drop is found to be approximately equal to the average kinetic energy of the equatorially trapped plasma. The potential is compared with that obtained by Alfven and Faelthammar (1963).
The separatrix electron density is an important parameter for core-edge scenario integration in tokamak devices, as it influences plasma confinement, divertor detachment and disruption avoidance. This quantity has been measured in H-mode discharges on JET, ASDEX Upgrade and Alcator C-Mod by applying the same fitting function to Thomson scattering measurements, and by employing the same analysis technique based on scrape-off layer power balance. To estimate the power crossing the separatrix, the inter-ELM time derivative of the plasma energy d W /d t has been experimentally evaluated and found to be approximately a constant fraction of the absorbed heating power. Correlations between n e,sep and engineering parameters have been investigated, revealing that n e,sep scales with the divertor neutral pressure p 0,div in a similar manner across all devices. Additionally, when n e,sep is normalized to the obtained n e,sep dependency, no clear correlation with the plasma current is found. These observations are in agreement with the 2-point model, which suggests that the upstream separatrix density is mainly set by the recycling at the divertor target.
This document details the work accomplished during the period 04/15/2012-01/31/2020 for the collaborative project “Interactions Between Microturbulence and Magnetohydrodynamic Instabilities in Magnetized Plasmas,” and serves as the final report on the University of California, San Diego part of this collaborative project for Department of Energy Office of Science - Fusion Energy Sciences Grant No. DE- SC0007783. This project focused on the development of analytic and computational tools, seeking to better understand the mutual interactions between small-scale temperature-gradient driven microturbulence and macroscopic magnetohydrodynamic (MHD) instabilities, particularly tearing modes, in magnetically confined plasmas relevant for the generation of fusion energy. This project advanced understanding through numerical simulations of different sets of model equations, corresponding to different microturbulence instabilities. Questions addressed included how the magnitudes, spatial structures, and scaling's of turbulent forcing's and drags impact magnetic island dynamics, as well as quantifying the “region of influence” of an island on the microturbulence dynamics. A key finding of the project is that the microturbulence fluctuations and associated cross-field particle and energy fluxes exhibit spatial asymmetries for sufficiently large islands. The results have improved our fundamental understanding of plasma dynamics, motivated new experimental investigations, and will help inform future fusion predictive modeling efforts.
We propose to continue a comprehensive and aggressive collaborative study of the critical issues of RFP and fusion physics at MST by enhancing the inferred ion distribution function and the high (spatial and temporal) resolution of magnetic fluctuation measurements. For this, we propose to make novel use of a retarding potential before a neutral particle analyzer as well as a novel millimeter magnetic probe. The measurement of relationships between ion and electron temperatures and magnetic field fluctuations can provide otherwise unobtainable insights on the role of turbulence in fusion plasmas with special focus on the impact of magnetic field turbulence on plasma confinement and scaling. These insights can be particularly valuable when they are coupled with an increased understanding of the issues which drive the turbulence onset, enhance steady-state sustainment, and enable turbulent dissipation. With these new insights, MST will be particularly well positioned to provide important insights in a manner portable to a broad class of fusion devices. We propose to exploit this opportunity and, at the same time, take full advantage of the associated educational enhancements for plasma physics students as FAMU continues to expand the diversity of the fusion scientific workforce.
The field of nuclear diagnostics for Inertial Confinement Fusion (ICF) is broadly reviewed from its beginning in the seventies to present day. During this time, the sophistication of the ICF facilities and the suite of nuclear diagnostics have substantially evolved, generally a consequence of the efforts and experience gained on previous facilities. As the fusion yields have increased several orders of magnitude during these years, the quality of the nuclear-fusion-product measurements has improved significantly, facilitating an increased level of understanding about the physics governing the nuclear phase of an ICF implosion. The field of ICF has now entered an era where the fusion yields are high enough for nuclear measurements to provide spatial, temporal and spectral information, which have proven indispensable to understanding the performance of an ICF implosion. At the same time, the requirements on the nuclear diagnostics have also become more stringent. To put these measurements into context, this review starts by providing some historical remarks about the field of ICF and the role of nuclear diagnostics, followed by a brief overview of the basic physics that characterize the nuclear phase and performance of an ICF implosion. A technical discussion is subsequently presented of the neutron, gamma-ray, charged-particle and radiochemistry diagnostics that are, or have been, routinely used in the field of ICF. This discussion is followed by an elaboration of the current view of the next-generation nuclear diagnostics. Since the seventies, the overall progress made in the areas of nuclear diagnostics and scientific understanding of an ICF implosion has been enormous, and with the implementation of new high-fusion-yield facilities world-wide, the next-generation nuclear diagnostics will play an even more important role for decades to come.
Designing magnets for three-dimensional plasma confinement is a key task for advancing the stellarator as a fusion reactor concept. Stellarator magnets must produce an accurate field while leaving adequate room for other components and being reasonably simple to construct and assemble. In this paper, a framework for coil design and optimization is introduced that enables the attainment of sparse magnet solutions with arbitrary restrictions on where coils may be located. The solution space is formulated as a 'wireframe' consisting of a mesh of interconnected wire segments enclosing the plasma. Two methods are developed for optimizing the current distribution on a wireframe: Regularized Constrained Least Squares, which uses a linear least-squares approach to optimize the currents in each segment, and Greedy Stellarator Coil Optimization, a fully discrete procedure in which loops of current are added to the mesh one by one to achieve the desired magnetic field on the plasma boundary. Examples are presented of solutions obtainable with each method, some of which achieve high field accuracy while obeying spatial constraints that permit easy assembly.
FIDASIM is a code that models signals produced by charge-exchange reactions between neutrals and ions (both fast and thermal) in magnetically confined plasmas. With the ion distribution function as input, the code predicts the efflux to a neutral particle analyzer diagnostic and the photon radiance of Balmer-alpha light to a fast-ion D α diagnostic, in addition to many other related quantities. A new, parallelized version of the Monte Carlo code FIDASIM has been developed in Fortran90 that is substantially faster than the original interactive data language version. Modified algorithms include more accurate treatments of the time dependent collisional-radiative equations that describe neutral energy levels, of the cloud of ‘halo’ neutrals that surround the injected neutral beam, and of finite Larmor radius effects. Enhanced physics capabilities include modelling ‘passive’ signals from cold edge neutrals, the ability to treat general three-dimensional magnetic confinement configurations, and calculations of diagnostic-specific weight functions that enable tomographic reconstructions of the fast-ion distribution function. Neutral beam attenuation, beam emission, and fast-ion birth profiles are also modelled. Finally, the new algorithms have been successfully validated against experimental data and new features have been tested through benchmarks between two independently developed versions of the code.
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.
Chaotic transport is a subject of paramount importance in a variety of problems in plasma physics, specially those related to anomalous transport and turbulence. On the other hand, a great deal of information on chaotic transport can be obtained from simple dynamical systems like two-dimensional area-preserving (symplectic) maps, where powerful mathematical results like KAM theory are available. In this work we review recent works on transport barriers in area-preserving maps, focusing on systems which do not obey the so-called twist property. For such systems KAM theory no longer holds everywhere and novel dynamical features show up as nonresistive reconnection, shearless curves and shearless bifurcations. After presenting some general features using a standard nontwist mapping, we consider magnetic field line maps for magnetically confined plasmas in tokamaks.
We present results of a quantitative analysis of structured plasma outflows above a polar coronal hole observed by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) spacecraft. In a 6 hr interval of continuous high-cadence SDO/AIA images, we identified more than 2300 episodes of small-scale plasma flows in the polar corona. The mean upward flow speed measured by the surfing transform technique is estimated to be 122 ± 34 km s−1, which is comparable to the local sound speed. The typical recurrence period of the flow episodes is 10–30 minutes, and the mean duration and transverse size of each episode are about 3–5 minutes and 3–4 Mm, respectively. The largest identifiable episodes last for tens of minutes and reach widths up to 40 Mm. For the first time, we demonstrate that the polar coronal-hole outflows obey a family of power-law probability distributions characteristic of impulsive interchange magnetic reconnection. Turbulent photospheric driving may play a crucial role in releasing magnetically confined plasma onto open field. The estimated occurrence rate of the detected self-similar coronal outflows is sufficient for them to make a dominant contribution to the fast-wind mass and energy fluxes and to account for the wind's small-scale structure.
Abstract It is estimated that pilot plants and reactors may experience rates of net erosion and deposition of solid plasma facing component (PFC) material of 10 3 –10 5 kg yr −1 . Even if the net erosion (wear) problem can be solved, the redeposition of so much material has the potential for major interference with operation, including disruptions due to so-called ‘unidentified flying objects (UFOs)’ and unsafe dust levels. The potential implications appear to be no less serious than for plasma contact with the divertor target: a dust explosion or a major UFO-disruption could be as damaging for an actively-cooled deuterium-tritium (DT) tokamak as target failure. It will therefore be necessary to manage material deposits to prevent their fouling operation. This situation appears to require a fundamental paradigm shift with regard to meeting the challenge of taming the plasma–material interface: it appears that any acceptable solid PFC material will in effect be flow-through , like liquid–metal PFCs, although at far lower mass flow rates. Solid PFC material will have to be treated as a consumable , like brake pads in cars. ITER will use high-Z (tungsten) armor on the divertor targets and low-Z (beryllium) on the main walls. The ARIES-AT reactor design calls for a similar arrangement, but with SiC cladding on the main walls. Non-metallic low-Z refractory materials such as ceramics (graphite, SiC, etc) used as in situ replenishable, relatively thin—of order mm—claddings on a substrate which is resistant to neutron damage could provide a potential solution for the main walls, while reducing the risk of degrading the confined plasma. Separately, wall conditioning has proven essential for achieving high performance. For DT devices, however, standard methods appear to be unworkable, but recently powder droppers injecting low-Z material ∼continuously into discharges have been quite effective and may be usable in DT devices as well. The resulting massive generation of low-Z debris, however, has the same potential to seriously disrupt operation as noted above. Powder droppers provide a unique opportunity to carry out controlled studies on the management of low-Z slag in all current tokamaks, independent of whether their protection tiles use low-Z or high-Z material.
Abstract The companion part A paper (Stangeby et al 2022) reports a number of independent estimates indicating that high-duty-cycle DT tokamaks starting with pilot plants will likely experience rates of net erosion and deposition of solid PFC, plasma facing component, material in the range of 10 3 to 10 4 kg yr −1 , regardless of the material used. The subsequent redeposition of such large quantities of material has the potential for major interference with tokamak operation. Similar levels and issues will be involved if ∼continuous low-Z powder dropping is used for surface conditioning of DT tokamaks, independent of the material used for the PFC armor. In Stangeby et al (2022) (part A) it is proposed that for high-duty-cycle DT tokamaks, non-metallic low-Z refractory materials such as ceramics (graphite, SiC, etc) used as in situ replenishable, relatively thin—of order mm—claddings on a substrate which is resistant to neutron damage could provide a potential solution for protecting the main walls, while reducing the risk of degrading the confined plasma. Assessment of whether such an approach is viable will require information, much of which is not available today. Section 6 of part A identifies a partial list of major physics questions that will need to be answered in order to make an informed assessment. This part B report describes R&D needed to be done in present tokamaks in order to answer many of these questions. Most of the required R&D is to establish better understanding of low-Z slag generation and to identify means to safely manage it. Powder droppers provide a unique opportunity to carry out controlled studies on the management of low-Z slag in current tokamaks, independent of whether their protection tiles use low-Z or high-Z material.
An initial high fluence campaign was performed in WEST, in 2023, on the newly installed actively cooled tungsten divertor composed of ITER-grade monoblocks. The campaign consisted in the repetition of a 60 s long Deuterium L-mode pulse in attached divertor conditions, cumulating over 10000s of plasma exposure. A maximum deuterium fluence of approximately 5⋅1026 m−2 was reached in the outer strike point region, representative of a few high performance ITER pulses. Gross tungsten erosion inferred from visible spectroscopy shows that the most eroded plasma facing component is the inner divertor target with rates ten times larger than on the outer divertor target. The outer midplane tungsten bumpers, located a few centimeters from the plasma, show gross erosion rates two times lower than at the outer divertor. We conclude that the outer midplane bumpers have a negligible contribution to the long range tungsten migration and deposition onto the lower divertor. The cumulated gross erosion rate on the inner divertor translates in an effective gross erosion thickness of about 20μm, while it is about 2μm for the outer divertor. Strikingly, these orderings coincide with the thickness of deposits found locally on the divertor: the exposed surfaces of high field side monoblocks are covered with several tens of μm tungsten deposits, while on the lower field side, few μm thin tungsten deposits are only found on the magnetically shadowed parts of monoblocks. The strong impact of those deposits on WEST operation, namely perturbation of surface temperature measurement with infra-red thermography, and the emission of flakes causing radiative perturbation of the confined plasma, calls for anticipating similar issues in ITER. In particular, the start of research operation shall consider the definition of a divertor erosion budget in order to anticipate the formation of deleterious deposits.
This work addresses discrepancies in the alignment of the H-mode pedestal profiles of the electron and ion properties in the DIII-D tokamak as measured by Thomson Scattering (TS) and Charge Exchange Recombination Spectroscopy (CER) diagnostics. While the alignment of these profiles is key for accurate studies of tokamak physics and plasma confinement, misalignments can occur due to inaccuracies, such as in magnetic equilibrium reconstructions required to map measurements in different poloidal and toroidal locations. Both FIDASIM, an established simulation package, and a simplified collisional radiative model are used to simulate neutral beam state densities and neutral beam emission. Simulated neutral beam emissions are calculated based on shifted TS profiles and compared to beam emission measurements from the Main Ion CER system to determine the best shift for aligning TS with CER. This analysis is performed on various DIII-D discharges.
Experiments with extrinsic impurity seeding in strongly negative triangularity shapes in DIII-D achieved radiated power fractions (relative to input power) of up to ≈85% total radiation and ≈55% core radiation in steady-operating conditions. The relationship between core and total radiation was sensitive to impurity species and input power. Attempts to reach higher radiation levels via higher impurity flows resulted in radiative collapse disruptions. Nitrogen, neon, argon, and krypton were tested. Injection was by gas puffing, usually controlled by feeding back real-time estimates for core or total radiating fractions ($P_\textrm{rad} / P_\textrm{input}$). Argon and krypton were controlled by feeding back total radiating fraction, whereas neon was controlled by feeding back the core radiating fraction, due to the lower efficiency of neon as a divertor radiator. Nitrogen flows were pre-programmed. Reasonable total $P_\textrm{rad}$ control target following was achieved with argon or krypton. Control was more challenging with the neon/core radiation configuration, which was more prone to slow response and overshooting of the control target. Poor particle removal contributed to the control challenge: neon particle inventory within the last closed flux surface was roughly constant for up to 1 s (the longest duration tested) after neon injection was halted. However, separate experiments with laser-blow off of non-recycling impurities measured a short impurity confinement time, on the scale of the energy confinement time of ∼100 ms. Modeling with the Aurora impurity transport simulation matched experimental neon density profiles with full recycling (R = 1.0) and weak pumping but predicted rapid decreases in neon inventory if pumping were increased or recycling decreased. This indicates that changes outside the confined plasma (adding a well-placed pump) would improve controllability for all highly recycling species.
In this paper, we show the expansion of the high-temperature regime and related physical phenomena in the Large Helical Device (LHD). At the LHD, we have developed a high-temperature operational regime for both ion (T i ) and electron (T e ) temperatures through understanding the physical characteristics of high-temperature plasma, the enhancement of plasma heating, and the development of technical methods such as wall recycling control and impurity control. In particular, the formation of the ion internal transport barrier (i-ITB) and electron internal transport barrier (e-ITB) in the plasma core region is the key to improving plasma performance. In addition, the deuterium experiment that began in 2017 showed that, due to the isotope effect, the thermal transport of both ions and electrons is suppressed in deuterium plasma compared to light hydrogen plasma, and that a high central temperature can be achieved efficiently. Furthermore, the combination of neutral beam injection (NBI) and electron cyclotron heating (ECH) has extended the operating conditions that simultaneously maintain high T i and T e , and comprehensive research on plasma confinement characteristics with an eye to future fusion reactor conditions has progressed.