Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magnetic fusion”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Tokamak ITG-KBM transition benchmarking with the mixed variables/pullback transformation electromagnetic gyrokinetic scheme

Electromagnetic gyrokinetic simulation of high temperature plasma is required to predict confinement in magnetic fusion devices and has posed challenges for existing codes. In this paper, we demonstrate successful global gyrokinetic simulation of the ion temperature gradient-driven mode-kinetic ballooning mode transition in a toroidal fusion plasma test case using the mixed variables/pullback transformation (MV/PT) scheme with the particle-in-cell codes XGC and ORB5, and compare to results from a conventional continuum code from the literature. Furthermore, the MV/PT scheme combines explicit time integration with mitigation of the well-known electromagnetic gyrokinetic “cancelation problem.” We calculate eigenmodes in the electrostatic and parallel vector potentials, and find good agreement in growth rate, real frequency, and the normalized plasma pressure of mode transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Data-driven assessment of magnetic charged particle confinement parameter scaling in magnetized liner inertial fusion experiments on Z

In magneto-inertial fusion, the ratio of the characteristic fuel length perpendicular to the applied magnetic field R to the α-particle Larmor radius $ϱ_α$ is a critical parameter setting the scale of electron thermal-conduction loss and charged burn-product confinement. Here, using a previously developed deep-learning-based Bayesian inference tool, we obtain the magnetic-field fuel-radius product BR ∝ R / $ϱ_α$ from an ensemble of 16 magnetized liner inertial fusion (MagLIF) experiments. Observations of the trends in BR are consistent with relative trade-offs between compression and flux loss as well as the impact of mix from 1D resistive radiation magneto-hydrodynamics simulations in all but two experiments, for which 3D effects are hypothesized to play a significant role. Finally, we explain the relationship between BR and the generalized Lawson parameter χ. Our results indicate the ability to improve performance in MagLIF through careful tuning of experimental inputs, while also highlighting key risks from mix and 3D effects that must be mitigated in scaling MagLIF to higher currents with a next-generation driver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Deuterium and tritium anomalous transport in the generalized Hasegawa–Wakatani resistive drift wave turbulence model with finite ion Larmor radius

Anomalous transport of mixed deuterium–tritium plasma in the edge of magnetic fusion reactors is investigated using numerical solutions of resistive drift wave turbulence model equations, including finite Larmor radius effects, that are derived within the generalized Hasegawa–Wakatani framework. The anomalous cross field diffusivities of deuterium and tritium are compared in turbulence regimes with different values of the electron adiabaticity parameter controlling the existence of zonal flow. The dependence of the tritium-to-deuterium diffusivity ratio on the deuterium and tritium densities and the logarithmic density gradients is analyzed, and a scaling relation is obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Dependence of Heat Removal Rate of Pebble-Based Rods on Inter-Pebble Matrix Fill Fraction

Carbon pebble rods are a promising candidate for use in high heat flux regions of magnetic fusion energy reactor walls. Under high (10 – 50 MW/m 2 ) heat loads, carbon pebble rods release hot pebbles from the exposed surface, carrying away heat as the pebble rod surface recedes. In this work, we show that the surface recession rate during heating can be adjusted by changing the mechanical strength of the extruded rods, modifying the heat removal rate; this is accomplished here by varying the fill fraction of the inter-pebble matrix. A three-dimensional finite element model is presented that captures many experimental observations, including the sphere temperature and the surface recession rate. The model predicts that pebble release is caused by thermally driven crack propagation through the matrix and that the matrix strength against breaking is the single most important material parameter setting the pebble release rate; this prediction is supported by experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Toward Verification of RANS Simulations of the T-Tube Modular Divertor Using Large Eddy Simulations of Impinging Turbulent Plane Jets

Turbulent impinging jets have been proposed to cool high heat flux plasma-facing components such as the solid tungsten target plates of the divertor in long-pulse magnetic fusion energy reactors. In particular, the T-tube modular divertor, originally developed by the ARIES Team, consists of two concentric cylindrical tubes where helium flows through a slot in the inner tube, forming an approximately planar jet that impinges upon and cools the inner surface of the pressure boundary (namely, the outer tube) and the ~15-cm 2 plasma-facing W target. The objective of this work is to demonstrate that large eddy simulations (LESs) accurately simulate the thermal transport in canonical flows that comprise the cooling flow in the T-tube, as well as validate temperatures from LES with experimental measurements in a simplified T-tube geometry. Wall‑resolved LESs, validated by experimental data and verified by direct numerical simulations (DNSs), provide benchmark data for two canonical flows in the T‑tube, namely, planar impinging and wall jets, for Reynolds numbers Re B = 4 × 10 3 to 2 × 10 4 . Our LES results are within 4% to 12% root-mean-square error (RMSE) of surface Nusselt number distributions (Nu) from experiments and DNSs. The validated LES results are then used as the ground truth to evaluate four Reynolds‑averaged Navier-Stokes (RANS) turbulence closures, namely, the k‑ω SST, realizable k‑ε, GEKO, and γ‑SST models. The k‑ω SST model has the best overall performance in terms of heat transfer, giving surface Nu within 12% RMSE of the LES results for high‑ReB impinging jets and reduced overprediction in the wall‑jet region. The GEKO model with default constants has the next best performance, providing slightly better Nu predictions for low ReB impinging jets (versus k-ω SST) but worse overall performance over the full range of ReB studied here. The realizable k‑ε turbulence model significantly overestimates turbulence near the stagnation point, while the γ‑SST model suppresses near‑wall production, biasing the simulations toward simulating laminar surface heat transfer. Simulations of the simplified T‑tube show that LES and RANS simulations with the k‑ω SST model give nearly identical average heat transfer coefficients (HTCs) over the impingement surface. The realizable k‑ε model predicts significantly lower wall temperatures due to overestimation of HTC in the outlet flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

One-dimensional simulation and validation of divertor detachment induced through nitrogen seeding on HL-2A

Divertor detachment is a promising method to solve the power exhaust problem in tokamak devices or even in future magnetic fusion reactors. In this work, a detachment experiment (HL-2A shot #38008) with mixed gas seeding (60% nitrogen and 40% deuterium) is simulated using the SD1D module in BOUT++. In the process from attachment to detachment, the target electron temperature and the target ion saturation current in simulations are found to be consistent with the experimental results measured by Langmuir probes on the target plate. In order to understand the underlying detachment mechanism on HL-2A, this work analyses the role of different particle species in the cases with different seeding rates. It shows that the plasma density varies little and the density of neutrals (D and D 2 ) slightly decreases with the increase of nitrogen seeding rate, such that plasma–neutral interactions cannot effectively reduce plasma energy and plasma momentum in the divertor. The case with a high seeding rate predicts that increasing seeding rate cannot reach a target temperature lower than 2.5 eV, which is the required temperature for strong plasma–molecule interactions. Thus, the plasma–molecule interactions may not be important in the divertor during nitrogen seeding. Here this work also studies the parallel forces including (1) force due to the parallel electric field, (2) friction forces, (3) ion- and electron-thermal forces, and (4) collisional reactions (e.g. charge exchange recombination and ionisation). It is found that the friction force between nitrogen ions and other particle species (primarily D + ) is the dominant force towards the target, while ion- and electron-thermal forces are the dominant force pushing nitrogen ions back to upstream. Parallel forces determine the parallel distribution of nitrogen impurities, and therefore decide the region of nitrogen radiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Radial drift of plasma blobs in a toroidal magnetic field with fully kinetic and reduced fluid models

Abstract In curved magnetic geometries, field-aligned regions of enhanced plasma pressure and density, termed ‘blobs,’ move as coherent filaments across the magnetic field lines. Coherent blobs account for a significant fraction of transport at the edges of magnetic fusion experiments and arise in naturally-occurring space plasmas. This work examines the dynamics of blobs with a fully kinetic electromagnetic particle-in-cell code and with a drift-reduced fluid code. In low-beta regimes with moderate blob speeds, good agreement is found in the maximum blob velocity between the two simulation schemes and simple analytical estimates. The fully kinetic code demonstrates that blob speeds saturate near the initial sound speed, which is a regime outside the validity of the reduced fluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Interpretative transport modeling of the WEST boundary plasma: main plasma and light impurities

Understanding impurity transport in tokamak plasmas is crucial to control radiative losses and material migration in future magnetic fusion reactors. In this work we deploy the SolEdge2D-EIRENE code to model the boundary plasma in a WEST discharge, satisfactorily reproducing measurements of both upstream and divertor plasma conditions. Additionally, the spatial distribution of oxygen, studied here as a representative light impurity, is compared to vacuum ultraviolet spectroscopy measurements acquired with an oscillating line of sight. The simulation captures a key feature of the experiment, namely a factor of sime2 higher oxygen brightness in the inner divertor region compared to the outer one. This spatial asymmetry in oxygen concentration is interpreted by analyzing the balance of friction forces and thermal gradient forces that the light impurity exchanges with the main plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Compressional Alfvén eigenmodes excited by runaway electrons

Compressional Alfvén eigenmodes (CAEs) driven by energetic ions have been observed in magnetic fusion experiments. In this paper, we show that the modes can also be driven by runaway electrons formed in post-disruption plasma, which may explain kinetic instabilities observed in DIII-D disruption experiments with massive gas injection. The spatial structure is calculated, as are the frequencies which are in agreement with experimental observations. Using a runaway electron distribution function obtained from a kinetic simulation, the mode growth rates are calculated and found to exceed the collisional damping rate when the runaway electron density exceeds a threshold value. The excitation of CAEs poses a new possible approach to mitigate seed runaway electrons during the current quench and surpassing the avalanche.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The dependence of tungsten fuzz layer thickness and porosity on tungsten deposition rate and helium ion fluence

Abstract Fuzz formation on a heated tungsten surface in the presence of a helium-containing plasma and tungsten deposition source was investigated. Tungsten samples were exposed at 1123 K to pure helium plasma with ion incident energy of 76 eV, W/He ion flux ratio of ∼ 0.4 × 10 − 4 , and varied helium ion fluence from 0.18 to 3.4 × 10 26 m −2 . Fuzz thickness was measured by cross-sectional scanning electron microscopy to increase from 0.22 to 15 µ m with increasing helium ion fluence. No indication of saturation in fuzz thickness at high fluence was observed, in contrast to fuzz produced on a tungsten surface without tungsten deposition. Additional tungsten samples were exposed at 1123 K to pure helium plasma with ion incident energy of 76 eV, helium ion fluence of ∼ 3.4 × 10 26 m −2 , and varied W/He ion flux ratio from 0.26 to 3.0 × 10 − 4 . Fuzz thickness increased from 7.5 to 120 µ m with increasing W/He ion ratio. A final sample exposed at 1123 K to a mixed helium-deuterium plasma with ion incident energy of 76 eV, helium ion fluence of 0.18 × 10 26 m −2 , and W/He ion flux ratio of 2.2 × 10 − 4 developed nearly identical fuzz structures to that developed in a pure He plasma. As a function of deposited tungsten fluence, all results were found to trace out a single layer-growth curve given by a power law relation, indicating that fuzz thickness is independent of the W/He ion flux ratio in the range investigated and independent of any deuterium present in the plasma. As a result, for tungsten plasma facing walls in magnetic fusion devices at 1000–2000 K with 10 −4 W/He ion flux ratio, fuzz with thicknesses greater than hundreds of microns may form in as little as 10 4 s (in the absence of ELM-induced erosion or annealing), and may more significantly affect its thermophysical properties than fuzz generated without a tungsten deposition source.

Physics↗
DOE OSTI2025

Efficient calculation of self magnetic field, self-force, and self-inductance for electromagnetic coils with rectangular cross-section

Abstract For designing high-field electromagnets, the Lorentz force on coils needs to be computed in order to design suitable support structures, and the inductance should be computed to evaluate the stored energy and dynamics. Also, the magnetic field and its variation inside the conductor is of interest for computing stress and strain, and due to superconducting quench limits. For these force, inductance, energy, and internal field calculations, the coils cannot be naively approximated as infinitesimally thin filaments due to divergences when the source and evaluation points coincide, so more computationally demanding calculations are usually required, resolving the finite cross-section of the conductors. Here, we present a new alternative method that enables the internal magnetic field vector, self-force, and self-inductance to be computed rapidly and accurately within a 1D filament model. The method is applicable to coils for which the curve center-line can have general noncircular shape, as long as the conductor width is small compared to the radius of curvature. This paper extends a previous calculation for circular-cross-section conductors (Hurwitz et al 2024 IEEE Trans. Magn. ) to consider the case of rectangular cross-section. The reduced model is derived by rigorous analysis of the singularity, regularizing the filament integrals such that they match the true high-dimensional integrals at high coil aspect ratio. The new filament model exactly recovers analytic results for a circular coil, and is shown to accurately reproduce full finite-cross-section calculations for a non-planar coil of a stellarator magnetic fusion device. Due to the efficiency of the model here, it is well suited for use inside design optimization.

Landreman, Matt (ORCID:000000027233577X)↗
DOE OSTI2025

MHD, disruptions and control physics: Chapter 4 of the special issue: on the path to tokamak burning plasma operation

In this chapter, we review the progress in MHD stability, disruptions and control in magnetic fusion research that has occurred over the past (more than) one and a half decades since the publication by Hender et al in 2007 on the same topic as part of the update of ITER Physics Basis. During this period, remarkable progress has been achieved in the understanding of the basic physics and overall control of MHD instabilities through a wide spectrum of dedicated experiments, theory and modeling. The sawtooth activities are probably today one of the best understood of MHD events and very robust control schemes have been developed for reliable operation of tokamaks through core heating. Similarly, significant improvements have been achieved in understanding and control of neoclassical tearing modes, resistive wall modes or locked modes and their control through ECCD or error field control. The field of disruption prediction through application of artificial intelligence, machine learning or deep learning methods, which had already started at the time of the 2007 review, has progressed significantly due to general progress in these fields and application of newer, more sophisticated algorithms. However, although remarkable progress has been achieved in the field of Disruptions, their understanding, prediction, possible avoidance and mitigation still remain probably the most active fields of R&D globally in this field. This is especially because reactor grade machines like ITER and DEMO will be much less tolerant in respect of disruptions and runaway currents, and their occurrences must be either avoided altogether or minimized to an acceptable value without causing any significant hindrance to robust machine operations. This review is intended to present a broad spectrum of the R&D that has occurred in this field in support of ITER, which will also be of immense significance for all future machines, especially reactors like DEMO.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Bayesian parameter estimation and evaluation of the K -ω shear stress transport model for plane impinging jets

Numerical simulations with semi-empirical turbulence models are commonly used to model impinging jets, often used for cooling solid surfaces. In this work, the constants in the k-ω shear stress transport model in ANSYS FLUENT are calibrated to experimental velocity and heat transfer data for a plane turbulent impinging air jet to determine if Kennedy-O'Hagan calibration (Kennedy and O'Hagan 2001 J. R. Stat. Soc. B 63 425–64) can improve predictions of near-surface velocities and surface Nusselt numbers for similar flows. Impinging jets have been proposed to cool the target plates of the divertor in future magnetic fusion energy reactors, where simulations are used to estimate divertor performance. The flat-plate divertor (Wang et al 2009 Fusion Sci. Technol .56 1023–7) uses a plane jet of helium issuing from a B = 0.5 mm slot to cool a surface with radius of curvature of 44 B at a distance 4 B from the slot. Predictions from the calibrated numerical model are compared with independent experimental data at different flow conditions, as well as surface temperature data for a flat plate divertor test section. The contribution of this work is evaluation of the accuracy of a calibrated turbulence model for modest extrapolations in flow geometry and flow conditions for a plane impinging jet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Explicit structure-preserving geometric particle-in-cell algorithm in curvilinear orthogonal coordinate systems and its applications to whole-device 6D kinetic simulations of tokamak physics

Explicit structure-preserving geometric particle-in-cell (PIC) algorithm in curvilinear orthogonal coordinate systems is developed. The work reported represents a further development of the structure-preserving geometric PIC algorithm achieving the goal of practical applications in magnetic fusion research. The algorithm is constructed by discretizing the field theory for the system of charged particles and electromagnetic field using Whitney forms, discrete exterior calculus, and explicit non-canonical symplectic integration. In addition to the truncated infinitely dimensional symplectic structure, the algorithm preserves exactly many important physical symmetries and conservation laws, such as local energy conservation, gauge symmetry and the corresponding local charge conservation. As a result, the algorithm possesses the long-term accuracy and fidelity required for first-principles-based simulations of the multiscale tokamak physics. The algorithm has been implemented in the SymPIC code, which is designed for high-efficiency massively-parallel PIC simulations in modern clusters. The code has been applied to carry out whole-device 6D kinetic simulation studies of tokamak physics. A self-consistent kinetic steady state for fusion plasma in the tokamak geometry is numerically found with a predominately diagonal and anisotropic pressure tensor. The state also admits a steady-state sub-sonic ion flow in the range of 10 km s -1 , agreeing with experimental observations and analytical calculations Kinetic ballooning instability in the self-consistent kinetic steady state is simulated. It is shown that high-n ballooning modes have larger growth rates than low-n global modes, and in the nonlinear phase the modes saturate approximately in 5 ion transit times at the 2% level by the E × B flow generated by the instability. These results are consistent with early and recent electromagnetic gyrokinetic simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Staircase resiliency in a fluctuating cellular array

Inhomogeneous mixing by stationary convective cells set in a fixed array is a particularly simple route to layering. Layered profile structures, or staircases, have been observed in many systems, including drift-wave turbulence in magnetic confinement devices. The simplest type of staircase occurs in passive-scalar advection, due to the existence and interplay of two disparate timescales, the cell turn-over (τ H ), and the cell diffusion (τ D ) time. In this simple system, we study the resiliency of the staircase structure in the presence of global transverse shear and weak vortex scattering. The fixed cellular array is then generalized to a fluctuating vortex array in a series of numerical experiments. The focus is on regimes of low-modest effective Reynolds numbers, as found in magnetic fusion devices. By systematically perturbing the elements of the vortex array, we learn that staircases form and are resilient (although steps become less regular, due to cell mergers) over a broad range of Reynolds numbers. The criteria for resiliency are (a) τ D >>τ H and (b) a sufficiently high profile curvature (κ ≥ 1.5). We learn that scalar concentration travels along regions of shear, thus staircase barriers form first, and scalar concentration "homogenizes" in vortices later. The scattering of vortices induces a lower effective speed of scalar concentration front propagation. The paths are those of the least time. We observe that if background diffusion is kept fixed, the cell geometric properties can be used to derive an approximation for the effective diffusivity of the scalar. Furthermore, the effective diffusivity of the fluctuating vortex array does not deviate significantly from that of the fixed cellular array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

NSTXU Diagnostic Disruption Dynamic Loading Represented by Response Spectra

This article presents the results of transient dynamic simulations of loads due to disruption eddy currents on the NSTXU vacuum vessel. Dynamic loading at diagnostic mounting locations is expressed as response spectra derived from the time history results of the dynamic structural simulations of a variety of disruption scenarios. The disruption simulations draw on a history of the project assessments of worst case disruptions for specific components. Major efforts to assess disruption loading have included the vacuum vessel which is the major structural support for the machine, as well as the passive plates (PPs), high harmonic fast wave (HHFW) antenna, and centerstack casing. Each one of these efforts included transient electromagnetic simulations producing time-dependent eddy current Lorentz loads (and in some cases halo loads) which then were applied to time-dependent structural dynamic analyses intended to obtain the proper dynamic amplification factors. In some instances, the EM model and structural model were identical allowing direct transfer of EM forces to the structural model. In other cases, the EM and structural model were not identical and the vector potential (VP) transfer method was used. The results files from these analyses were available (or re-run) to post process in ANSYS Classic time history postprocessor. In conclusion, the ANSYS command is used to create response spectra from time history data at desired points on the vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC

The global total-f gyrokinetic particle-in-cell code XGC, used to study transport in magnetic fusion plasmas or to couple with a core gyrokinetic code while functioning as an edge gyrokinetic code, implements a 5-dimensional (5D) continuum grid to perform the dissipative operations, such as plasma collisions, or to exchange the particle distribution function information with a core code. To transfer the distribution function between marker particles and a rectangular 2D velocity-space grid, XGC employs a bilinear mapping. The conservation of particle density and momentum is accurate enough in this bilinear operation, but the error in the particle energy conservation can become undesirably large and cause non-negligible numerical heating in a steep edge pedestal. In the present work we update XGC to use a novel mapping technique, based on the calculation of a pseudo-inverse, to exactly preserve moments up to the order of the discretization space. We describe the details of the implementation and we demonstrate the reduced interpolation error for a tokamak test plasma by using 1st- and 2nd-order elements with the pseudo-inverse method and comparing to the bilinear mapping.

Fusion Plasma↗