Learning from Each Other: Cross-Cutting Diagnostic Development Activities Between Magnetic and Inertial Confinement Fusion
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Accurately incorporating non-local transport into radiation-hydrodynamics codes, and indeed any fluid system, has long been elusive. To date, a simplified and accurate theory that can be easily integrated has not been available. This limitation affects modeling in inertial confinement fusion (ICF) and magnetic confinement fusion systems, among others, where non-local transport is well-known to be present. Here, we present a coupling methodology between a full Vlasov–Fokker–Planck (VFP) electron kinetic code and radiation-hydrodynamics (rad-hydro) codes. The VFP code is used to adjust native electron transport in the rad-hydro code, thus enabling improved transport without the need to integrate a full electron VFP solver into the rad-hydro code. This approach necessitates only occasional invocation of the VFP code, reducing computational intensity compared to following the dynamic evolution entirely with the VFP code on fluid time scales. We illustrate that the methodology is more accurate than other simplified methods in thermal decay systems relevant to ICF and can replicate standard theoretical results with high accuracy.
Plasma radiation plays an important role in the operation and control of magnetic confinement fusion and is typically measured by bolometers. Here, we study the influence of mechanical vibration and magnetic field, which are typically present in magnetic confinement fusion devices, on the performance of a fiber-optic bolometer (FOB) based on silicon Fabry-Perot interferometer (FPI) demodulated by a scanning diode laser. It is found that both the vibration and magnetic field could cause large noise due to the birefringence of the silicon FPI. To mitigate the birefringence effects, we demonstrated two effective methods: polarization maintaining FOB and polarization scrambling. We show that two methods can be effective in mitigating the effect caused by the sensor birefringence: (1) replacing the regular single-mode fiber in the FOB system with polarization-maintaining fiber and (2) using a polarization scrambler after the laser source. For both methods, no significant increase in the noise was observed in the vibration and magnetic field tests. Finally, our research highlights the importance of the polarization management for the FOB systems toward their practical applications in magnetic-confinement fusion systems.
Not Available
Magnetizing a cryogenic deuterium–tritium (DT)-layered inertial confinement fusion (ICF) implosion can improve performance by reducing thermal conduction and improving DT-alpha confinement in the hot spot. A room-temperature, magnetized indirect-drive ICF platform at the National Ignition Facility has been developed, using a high-Z, high-resistivity AuTa 4 alloy as the hohlraum wall material. Experiments show a 2.5× increase in deuterium–deuterium (DD) neutron yield and a 0.8-keV increase in hot-spot temperature with the application of a 12-T B-field. Further, for an initial 26-T B-field, we observed a 2.9× yield increase and a 1.1-keV temperature increase, with the inferred burn-averaged B-field in the compressed hot spot estimated to be 7.1 ± 1.8 kT using measured primary DD-n and secondary DT-n neutron yields.
Accurately predicting redeposition is vital for high-Z plasma-facing component (PFC) survivability in magnetic confinement fusion. In this study, we categorize high-Z redeposition into three mechanisms: geometric-driven (prompt), sheath-driven (local), and scrape-off-layer-driven (far) redeposition. To investigate these mechanisms, we employ Monte Carlo transport codes to simulate azimuthally symmetric tungsten source erosion and redeposition in a tokamak. By iteratively analyzing critical parameters, we evaluate redeposition scaling for each mechanism. Specifically, we investigate the impact of magnetic-field-to-PFC pitch angle assumptions on PFC losses into the scrape-off layer. Our findings reveal significant pitch angle sensitivity due to an asymmetric prompt vs local redeposition trade-off. These results enhance our understanding of redeposition phenomena in fusion plasma environments.
The dynamics of edge-core coupling is critically important to the optimization of magnetically confined fusion plasmas. Since early proposals, there has been persistent speculation that inward propagation of turbulence from the boundary is a possible means to energize the edge-core coupling region. However, the detailed mechanism of this process has remained a mystery until recent experiments observed that regular, intense gradient relaxation events generated blob-void pairs very close to the last closed flux surface. Blobs ($\tilde{n}$ >0) propagate outward and detach from the bulk plasma, while voids ($\tilde{n}$ <0) propagate inward, and so stir the core plasma. Here, in this work, we demonstrate that this heretofore ignored process of void emission can drive a broad turbulent layer of width ∼100 𝜌 𝑠 , for typical parameters. The mechanism is the Cherenkov emission of drift waves from inward-propagating voids. The model shows promise to resolve several questions surrounding the shortfall problem and the strong turbulence in the edge-core coupling region.
Experiments on the DIII-D tokamak have identi fied a novel regime in which applied resonant magnetic perturbations (RMPs) increase the particle con finement and overall performance. This work details a robust range of counter-current rotation over which RMPs cause this density pump-in effect for high con finement (H-mode) plasmas. The pump-in is shown to be caused by a reduction of the turbulent transport and to be correlated with a change in the sign of the induced neoclassical transport. Furthermore, this novel reversal of the RMP induced transport has the potential to signifi cantly improve reactor relevant, three-dimensional magnetic con finement scenarios.
Turbulent transport is regarded as one of the key issues in magnetic confinement nuclear fusion, both for tokamaks and stellarators. Here, in this work, we show that a significant decrease in a microstability-based proxy, as opposed to a geometric one, for the turbulent heat flux, namely the quasilinear heat flux, can be obtained in an efficient manner by coupling stellarator optimization with linear gyrokinetic simulations. This is accomplished by computing the quasilinear heat flux at each step of the optimization process, as well as the deviation from quasisymmetry, and minimizing their sum, leading to a balance between neoclassical and the turbulent transport proxy.
Diagnosing plasma magnetization in inertial confinement fusion implosions is important for understanding how magnetic fields affect implosion dynamics and to assess plasma conditions in magnetized implosion experiments. Secondary deuterium–tritium (DT) reactions provide two diagnostic signatures to infer neutron-averaged magnetization. Magnetically confining fusion tritons from deuterium–deuterium (DD) reactions in the hot spot increases their path lengths and energy loss, leading to an increase in the secondary DT reaction yield. In addition, the distribution of magnetically confined DD-triton is anisotropic, and this drives anisotropy in the secondary DT neutron spectra along different lines of sight. Implosion parameter space as well as sensitivity to the applied B-field, fuel ρR, temperature, and hot-spot shape will be examined using Monte Carlo and 2D radiation-magnetohydrodynamic simulations.
The viability of magnetic confinement fusion as an energy source depends on achieving the high ion temperatures required for D-T fusion. Among the available techniques, plasma heating with waves in the ion cyclotron range of frequencies (ICRF) is a prominent method for bulk ion heating in fusion plasmas. Furthermore, a detailed understanding of the non-linear physics of alpha heating and the complex impact of MeV-range fast ions on plasma dynamics becomes progressively more important. This paper provides a comprehensive overview of recent developments with the three-ion ICRF scenarios on Alcator C-Mod, ASDEX Upgrade and JET tokamaks. The results demonstrate the flexibility of these novel scenarios for heating bulk ions in D-T ≈ 50%-50% plasmas and efficient generation of MeV-range fast ions in multi-ion species plasmas. In conclusion, several key results relevant for ITER and future fusion reactors are highlighted.
The formation of the plasma current spike at the end of the thermal quench (TQ) phase is studied systematically, which is found to strongly correlate with the onset of the 1/1 kink-tearing reconnection only in the simulation results presented here. Here, the magnetohydrodynamic activity on the q = 1 surface plays a critical role in the spike formation and the disruption process. In particular, when the safety factor in the magnetic axis q 0 exceeds 1, the plasma major disruption transits into successive minor disruptions, and the start of the TQ phase is delayed.
Transient magneto-hydrodynamic (MHD) events like edge localized modes (ELMs) or disruptions are a concern for magnetic confinement fusion power plants. Research with the MHD code JOREK towards understanding control of such instabilities is reviewed here in a concise way to provide a complete overview, while we refer to the original publications for details. Experimental validation for unmitigated vertical displacement events progressed. The mechanism of vertical force mitigation by impurity injection was identified. Two-way eddy current coupling to CARIDDI was completed. Shattered pellet injection was simulated in JET, KSTAR, ASDEX Upgrade (AUG) and ITER. Benign runaway electron beam termination in JET and ITER was studied. Coupling of kinetic REs to the MHD is ongoing and a virtual RE synchrotron radiation diagnostic was developed. Regarding pedestal physics, regimes devoid of large ELMs in AUG were simulated and predictive JT60-SA simulations are ongoing. For ELM suppression by resonant magnetic perturbations (RMPs), AUG, ITER and EAST simulations were performed. A free boundary RMP model was validated against experiments. Evidence for penetrated magnetic islands at the pedestal top based on AUG experiments and simulations was found. Simulations of the naturally ELM-free quiescent H-mode in AUG and HL-3 show external kink mode formation prevents pedestal build-up towards an ELM within windows of the edge safety factor. With kinetic neutral particles, high field side high density formation in ITER was simulated and with kinetic impurities, tungsten transport in AUG RMP plasmas was studied. To capture turbulent transport, electro-static full-f particle in cell models for ion temperature gradient and trapped electron modes were established and benchmarked. Application to RMP plasmas shows enhanced turbulence in comparison to unperturbed states. Energetic particle interactions with MHD were studied. Flux pumping that prevents the safety factor on axis from dropping below unity was simulated. First non-linear stellarator applications include current relaxation in $l$ = 2 stellarators, while verification for advanced stellarators progresses.
It is well documented that the central electron temperature in the national spherical torus experiment (NSTX) remains largely unchanged as the external heating power, and hence the normalized volume averaged plasma pressure β increases [Stutman, Phys. Rev. Lett. 102, 115002 (2009)]. Herein we present a hypothesis that low n, pressure driven ideal magnetohydrodynamic (MHD) instabilities that are nondisruptive, can break magnetic surfaces in the central region and thereby flatten the electron temperature profiles. We demonstrate this mechanism in a 3D resistive MHD simulation of a NSTX discharge. By varying the toroidal magnetic field strength, and/or the heating power, we show that there is a critical value of β, above which the central temperature profile no longer peaks on axis.
Applying an external magnetic field to laser-driven inertial confinement fusion implosions is a promising approach for enhancing fusion yield. The field is compressed with the plasma, producing a magnetized hotspot that anisotropically suppresses thermal losses and traps alpha particles, making performance sensitive to the compressed field orientation. Here, we derive a simple, readily applicable analytic model that enables rapid evaluation of the compressed field topology and show that ablation into the hotspot amplifies the central field, while the ablated ice near the hotspot edge develops a decaying, radially bent field, with a discontinuity in the field direction. The radially bent field renders thermal insulation at the hotspot edge negligible and largely independent of the applied field strength, whereas insulation in the hotspot core still depends strongly on the applied field. Applying the model to non-axial initial field configurations, we find that an initially applied mirror field provides the greatest suppression, followed by the standard axial field.
By producing localized wave regions at the ends of an open-field-line magnetic confinement system, ponderomotive walls can be used to differentially confine different species in the plasma. Furthermore, if the plasma is rotating, this wall can be magnetostatic in the laboratory frame, resulting in simpler engineering and better power flow. However, recent work on such magnetostatic walls has shown qualitatively different potentials than those found in the earlier, nonrotating theory. Here, using a simple slab model of a ponderomotive wall, we resolve this discrepancy. We show that the form of the ponderomotive potential in the comoving plasma frame depends on the assumption made about the electrostatic potential in the laboratory frame. If the laboratory-frame potential is unperturbed by the magnetic oscillation, one finds a parallel-polarized wave in the comoving frame, while if each field line remains equipotential throughout the perturbation region, one finds a perpendicularly polarized wave. As a result, this in turn dramatically changes the averaged ponderomotive force experienced by a charged particle along the field line, not only its scaling, but also its direction.
Highly accurate simulation of plasma transport is needed to drive the successful design and operation of magnetically confined fusion reactors. Unfortunately, the extreme anisotropy present in magnetized plasmas results in thin boundary layers that are expensive to resolve. Here, this work investigates how various mesh refinement strategies might reduce that expense to allow for more efficient simulation by comparing standard variable refinement approaches that use a field quantity to an adaptive approach that uses an error estimator. It is first verified that higher order discretization only realizes the proper rate of convergence once the mesh resolves the thin boundary layer, therefore motivating the focusing of refinement on the boundary layer. For three two-dimensional test cases that contain characteristic features of tokamak magnetic fields, an exponential refinement strategy based on the magnetic flux function, which is the standard refinement approach in the field, is compared to an adaptive strategy utilizing the established Zienwiekicz and Zhu error estimator. The adaptive mesh refinement strategy consistently achieves the same accuracy using orders of magnitude less degrees of freedom than either exponential or uniform refinement. This result makes the adaptive refinement strategy more efficient than the exponential refinement strategy while also being more generalizable to problems with complex magnetic geometries. Scaling laws are derived that quantify the improvement in cost of the adaptive refinement strategy over other refinement approaches in 2D and 3D.
Recent simulation and experimental results suggest that the magnetic island and flow on a resonant surface often do not satisfy the “no-slip” condition in the steady state. A new theoretical model on a nonlinear plasma response to external magnetic perturbation in the absence of a no-slip condition is proposed. The model is composed of the equations for the evolution of both the width and phase of the magnetic island due to forced reconnection driven by the external magnetic perturbation, and the force-balance equation for the plasma flow. When the island width is much less than the resistive layer width, the island growth is governed by the linear Hahm–Kulsrud–Taylor solution in the presence of time-dependent plasma flow. In the other regime when the island width is much larger than the resistive layer width, the evolution of both the island width and phase can be described using the Rutherford theory. The island solution is used to construct the quasilinear electromagnetic force, which, together with the viscous one, contributes to the nonlinear variation in plasma flow. The no-slip condition assumed in the conventional error field theory is not imposed here, where the island oscillation frequency depends on but does not necessarily equal to the plasma flow frequency at the rational surface.