Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cyclotrons”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Full-wave modeling of RF waves in fusion plasmas with finite element method: Progress in past decades and its future role

This paper reviews the progress in computing radio frequency (RF) wave fields in fusion plasmas, specifically focusing on simulations utilizing the finite element method (FEM) over the past few decades. Computing RF wave fields in fusion plasmas presents unique challenges due to large simulation domains, complex antenna structures, non-local dielectric properties, and wide ranges of spatial scales. It highlights key developments and outlines future directions, primarily addressing waves in the ion cyclotron (IC) to lower hybrid (LH) frequency range. We begin with briefly revisiting earlier developments before the widespread availability of modern computer-aided engineering (CAE) software based on FEM. This historical perspective illuminates early progress and the physics difficulties that motivated ongoing work within the community. Modern wave simulations for RF antennas based on FEM are characterized by the use of detailed 3D antenna model geometry generated from engineering CAD software and localized wave dielectric model. Significant advancements have also been made in improving physics models to include phenomena such as RF sheath rectification and wave scattering. FEM-based RF simulations have also been applied to compute wave propagation in the core region, where the inclusion of non-local dielectric response is crucial. This is a challenging goal, and several promising approaches have been proposed in this area. Additionally, RF simulation development initiatives based on open-source libraries have gained popularity, demonstrating scalability and flexibility in extending physics models. This paper will discuss the advantages and disadvantages of using such a publicly available FEM library.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Preliminary Result of Neoclassical Impurity Transport Affected by ICRH with a 4-D Fokker-Planck Code

To understand the role of ion cyclotron resonance heating (ICRH) in mitigating tungsten accumulation, we evaluate the impact of ICRH-induced poloidal asymmetries on neoclassical impurity transport using a newly developed 4-D Fokker–Planck code (FP4D). Unlike conventional approaches that separate RF heating and transport modelling, FP4D solves both fast-ion distribution formation and impurity transport self-consistently without bounce-averaging. Benchmarking against the NEO code confirms the accuracy of the neoclassical physics in FP4D. Using FP4D, we analyze anisotropic and non-Maxwellian minority distributions under ICRH and their effects on the equilibrium potential and tungsten transport. Results show that anisotropy induced poloidal potential variations can reduce inward tungsten flux, depending on the modelling approach. Our findings highlight the importance of resolving the full velocity and poloidal structure of minority species in predicting impurity transport under strong RF heating conditions.

Jeong, Yunho [Hanyang Univ., Seoul (Korea, Republi↗
DOE OSTI2026

Overview of the ICRF heating system in SPARC

Ion Cyclotron Range of Frequencies (ICRF) heating is a critical system for the SPARC mission of demonstrating net positive fusion gain. In the first campaign, 20 MW of installed power at 120 MHz will be supplied to 10 four-strap antennas, arranged in poloidal pairs. The matching network is designed to accommodate for a range of loading conditions and a fast controller adjusts the matching during a pulse with frequency modulation of ±1 MHz. For the main L-mode scenario of the first experimental campaign, the power coupling and absorption is analysed in this work. Several feeding schemes are compared and limits on the amount of the total coupled power are evaluated. Coupling of the full 20 MW power is possible well within the limit of the peak electric field in the transmission line of 15 kV/cm. The 3 He minority and 2 nd harmonic T heating scheme is validated in full-wave modelling to be efficient for heating plasma in the SPARC L-mode Q>1 scenario, with dominant ion heating and efficient single-pass absorption. A range of ICRF wave parameters and minority concentration would be suitable for operation and allow tailoring the heating properties for plasma conditions and experimental objectives.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

ICRF near-field effects sensitivity to magnetic tilt-angle

This paper complements the study reported in G. Urbanczyk et al 2025 Nucl. Fusion 65 046018 by investigating how near-field effects from Ion Cyclotron Range of Frequency (ICRF) antennas depend on the magnetic geometry, particularly focusing on the sensitivity of rectified sheath potentials and associated impurity production to variations in magnetic field tilt-angle. Simulations were conducted using Petra-M, employing a flat model of the ASDEX Upgrade (AUG) three-strap antenna, and COMSOL, utilizing a curved model of the WEST antenna. Results demonstrate indications of sensitivity of the rectified sheath potentials and impurity sputtering to the alignment of Faraday screen bars with the magnetic field lines, emphasizing the necessity for careful consideration of geometric alignment in antenna design for fusion devices. This trend is consistent with experimental data measured with various diagnostics.

Urbanczyk, Guillaume [Universite de Lorraine, Nanc↗
DOE OSTI2026

Simulations of ICRF Heating for SPARC during First Campaign and Primary Reference-Like Discharge using the Stix Code

High magnetic field tokamaks, like SPARC, rely on ion cyclotron radio frequency heating (ICRF) to reach fusion relevant temperatures. The SPARC tokamak will have 14 ICRF antennas in 7 toroidal locations delivering > 20 MW of power to the plasma. New capabilities with the full wave cold plasma solver, Stix, now allow for resolving the wave-particle resonances using lower order thermal corrections to capture core absorption of Landau damping and ion resonances in devices like SPARC. Favorable comparisons to the TORIC codes give confidence in the single pass absorption of this model to accurately capture the strength of edge interactions of the RF. Using this new dielectric formulation in the Stix code, simulations of the 2D poloidal cross section of SPARC are completed for the first campaign and primary reference-like discharges (PRD-like). A scan of the minority ion concentrations of helium-3 is performed and shows the expected behavior that as the helium-3 decreases the amount of single pass absorption also decreases which is seen in both scenarios. Additionally, both scenarios show only slight differences in single-pass absorption for the range of 3% to 5% helium-3 allowing for more flexibility in experiments. This study also highlights the differences between the first campaign and PRD-like with the first campaign discharges showing much more multi-pass absorption and an effect of confining the wave to a smaller portion of the cross-section due to the fast wave cut-off. This latter result suggests that far-field sheath rectification at the high-field side would be minimal for the first campaign scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

DIII-D High Field Side Lower Hybrid Current Drive Experiment Overview

In preparation for high field side lower hybrid current drive (HFS LHCD) experiments in DIII-D, the HFS LHCD was to be commissioned and physics experiments commence once the system operated up to 300 kW for 0.5 s. The initial physics experiments sought to characterize coupling, wave propagation, and driven current measurements. In HFS LHCD first campaign, the maximum power was limited to <200 kW due to waveguide pressure leaks limiting the number of available modules and power per module. Here we summarize commissioning progress and initial physics observations. The HFS LHCD coupler was optimized for high qmin, DIII-D discharges where efficient off-axis current at r/a~0.6-0.8 is desired. The coupler n|| spectrum is peaked at 2.7 and is predicted to generate ~0.14 MA/MW coupled for 1.6 MW injected. In preparation, the HFS scrape-off layer density profile was characterized and found to have steeper profiles and lower fluctuation levels than the low field side. Furthermore, the HFS SOL density profile can be accurately predicted using global plasma quantities using machine learning. Thus far, one module has injected ~100 kW for 0.5 s with <5% reflected power. Nonthermal electrons have been observed on lower frequency channels of the electron cyclotron emission radiometer correlated with the LH power indicating core wave absorption. To avoid 30R neutral beam heat flux, a split launcher is proposed to avoid the high heat flux region while maintaining power spectrum and directivity.

Wukitch, Stephen J. [Massachusetts Inst. of Techno↗
DOE OSTI2020

Analytic stability boundaries for compressional and global Alfvén eigenmodes driven by fast ions. II. Interaction via Landau resonance

We report that conditions for net fast ion drive are derived for beam-driven, co-propagating, sub-cyclotron compressional (CAE) and global (GAE) Alfvén eigenmodes driven by the Landau resonance with super-Alfvénic fast ions. Approximations applicable to realistic neutral beam distributions and mode characteristics observed in spherical tokamaks enable the derivation of marginal stability conditions for these modes. Such conditions successfully reproduce the stability boundaries found from numerical integration of the exact expression for local fast ion drive/damping. Coupling between the CAE and GAE branches of the dispersion due to finite ω/ω ci and |k ∥ /k ⊥ | is retained and found to be responsible for the existence of the GAE instability via this resonance. Finally, encouraging agreement is demonstrated between the approximate stability criterion, simulation results, and a database of NSTX observations of co-CAEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Deep convolutional neural networks for multi-scale time-series classification and application to tokamak disruption prediction using raw, high temporal resolution diagnostic data

In this paper we discuss recent advances in deep convolutional neural networks (CNN) for sequence learning, which allow identifying long-range, multi-scale phenomena in long sequences, such as those found in fusion plasmas. We point out several benefits of these deep CNN architectures, such as not requiring experts such as physicists to hand-craft input data features, the ability to capture longer range dependencies compared to the more common sequence neural networks (recurrent neural networks like long short-term memory (LSTM) networks), and the comparative computational efficiency. We apply this neural network architecture to the popular problem of disruption prediction in fusion energy tokamaks, utilizing raw data from a single diagnostic, the Electron Cyclotron Emission imaging (ECEi) diagnostic from the DIII-D tokamak. Initial results trained on a large ECEi dataset show promise, achieving an F 1 -score of ~91% on individual time-slices using only the ECEi data. This indicates the ECEi diagnostic by itself can be sensitive to a number of pre-disruption markers useful for predicting disruptions on timescales not only for mitigation but also avoidance. Future opportunities for utilizing these deep CNN architectures with fusion data are outlined, including impact of recent upgrades to the ECEi diagnostic.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Measurement and modeling of the radio frequency sheath impedance in a large magnetized plasma

The DC and RF properties of radio frequency (RF) driven sheaths were studied in the Large Plasma Device (LAPD) at the University of California, Los Angeles. The experiments diagnosed RF sheaths on field lines connected to a grounded plate at one end and an ion cyclotron range of frequencies (ICRF) antenna at the other end. The experimental setup permitted measurement of the RF sheath impedance at the plate as a function of DC sheath voltage, with the latter controlled by varying the RF current applied to the antenna. Here, the DC current-voltage characteristics of these sheaths and the RF sheath impedance measurements were compared with modeling. Hot electrons, present in the LAPD plasma, were inferred to contribute significantly to both the DC and RF currents and hence the RF impedance. It was postulated that at very low power hot electrons could not access the region of the plasma subject to RF waves, resulting in an increased RF impedance. Within some experimental limitations and significant assumptions, an RF sheath impedance model was verified by the experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Normal modes, rotational inertia, and thermal fluctuations of trapped ion crystals

The normal modes of a trapped ion crystal are derived using an approach based on the Hermitian properties of the system’s dynamical matrix. This method is equivalent to the standard Bogoliubov method, but for classical systems it is arguably simpler and more general in that canonical coordinates are not necessary. The theory is developed for stable, unstable, and neutrally-stable systems. The method is then applied to ion crystals in a Penning trap. Reduced eigenvalue problems for the case of large applied magnetic field are developed, for which the spectrum breaks into E x B drift modes, axial modes, and cyclotron modes. Thermal fluctuation levels in these modes are analyzed and shown to be consistent with the Bohr-van-Leeuwen theorem, provided that neutrally-stable modes associated with crystal rotations are included in the analysis. An expression for the rotational inertia of the crystal is derived, and a magnetic contribution to this inertia, which dominates in large magnetic fields, is described. An unusual limit is discovered for the special case of spherically-symmetric confinement, in which the rotational inertia does not exist and changes in angular momentum leave the rotation frequency unaffected.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Evidence of E × B staircase in HL-2A L-mode tokamak discharges

Theory and simulation suggest the existence of an E ×B staircase in the plasma core consisting of a series of nested m/n = 0/0 E×B shear layers that regulate turbulent transport across the layers together with mesoscale transport events occurring between them. Here, we show evidence for these phenomena in HL-2A L-mode discharges. Both high resolution electron cyclotron emission Te profiles and frequency modulated continuous wave reflectometer ne profiles show multiple gradient corrugations. The analysis of simultaneous poloidally and radially separated beam emission spectroscopy turbulent density fluctuation measurements over most of the minor radius of the plasma provides the evidence of eddy tilting and propagation direction reversals largely coincident with the profile corrugations, with long-range transport events occurring between these locations. The results provide significant additional evidence for the existence of an E × B staircase that influences particle and heat transport in L-mode discharges.

Liu, Wenbin (ORCID:0000000310424117)↗
DOE OSTI2021

Transport coefficients for magnetic-field evolution in inviscid magnetohydrodynamics

The magnetized resistivity and electrothermal tensors when substituted into the induction equation lead to electrothermal magnetic field generation, resistive magnetic diffusion, and magnetic field advection due to resistivity gradients, temperature gradients and currents. Here, the advection terms driven by temperature gradient and current have cross-field components (perpendicular to both the magnetic field and the driving term) that depend on significantly modified versions of Braginskii’s transport coefficients. The improved fits to Braginskii’s coefficients given by Epperlein and Haines and Ji and Held give physically incorrect results for cross-field advection at small Hall parameters (product of cyclotron frequency and collision time). The errors in Epperlein and Haines’ fits are particularly severe, giving increasing advection velocities below a Hall parameter of one when they should decrease linearly to zero. Epperlein and Haines’ fits can also give erroneous advection terms due to variations in effective atomic number. The only serious error in Braginskii’s fits is an overestimate in advection due to perpendicular resistivity. New fits for the cross-field advection terms are obtained from a direct numerical solution of the Fokker-Planck equation and Ji and Held’s higher order expansion approach that are continuous functions of the effective atomic number.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Performance demonstration of vacuum microwave components critical for the operation of the ITER low-field side reflectometer

Final design studies in preparation for manufacturing have been performed for functional components of the vacuum portion of the ITER Low-Field Side Reflectometer (LFSR). These components consist of an antenna array, electron cyclotron heating (ECH) protection mirrors, phase calibration mirrors, and vacuum windows. Evaluation of these components was conducted at the LFSR test facility and DIII-D. The antenna array consists of six corrugated-waveguide antennas for simultaneous profile, fluctuation, and Doppler measurements. A diffraction grating, incorporated into the plasma-facing miter bend, provides protection of sensitive components from stray ECH at 170 GHz. For in situ phase calibration of the LFSR profile reflectometer, an embossed mirror is incorporated into the adjacent miter bend. Measurements of the radiated beam profile indicate that these components have a small, acceptable effect on mode conversion and beam quality. Baseline transmission characteristics of the dual-disk vacuum window are obtained and are used to guide ongoing developments. Preliminary simulations indicate that a surface-relief structure on the window surfaces can greatly improve transmission. The workability of real-time phase measurements was demonstrated on the DIII-D profile reflectometer. The new automated real-time analysis agrees well with the standard post-processing routine.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

System-on-chip upgrade of millimeter-wave imaging diagnostics for fusion plasma

Monolithic, millimeter wave “system-on-chip” technology has been employed in chip heterodyne radiometers in a newly developed Electron Cyclotron Emission Imaging (ECEI) system on the DIII-D tokamak for 2D electron temperature and fluctuation diagnostics. The system employs 20 horn-waveguide receiver modules each with customized W-band (75–110 GHz) monolithic microwave integrated circuit chips comprising a W-band low noise amplifier, a balanced mixer, a ×2 local oscillator (LO) frequency doubler, and two intermediate frequency amplifier stages in each module. Compared to previous quasi-optical ECEI arrays with Schottky mixer diodes mounted on planar antennas, the upgraded W-band array exhibits >30 dB additional gain and 20× improvement in noise temperature; an internal eight times multiplier chain is used to provide LO coupling, thereby eliminating the need for quasi-optical coupling. The horn-waveguide shielding housing avoids out-of-band noise interference on each module. The upgraded ECEI system plays an important role for absolute electron temperature and fluctuation measurements for edge and core region transport physics studies. An F-band receiver chip (up to 140 GHz) is under development for additional fusion facilities with a higher toroidal magnetic field. Visualization diagnostics provide multi-scale and multi-dimensional data in plasma profile evolution. Finally, a significant aspect of imaging measurement is focusing on artificial intelligence for science applications.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Design of microwave broadband CMOS transmitter and receiver circuits for MIR and ECEI plasma diagnostics

To efficiently determine the plasma electron density fluctuations using the MIR diagnostic technique, a 55–75 GHz 65 nm-CMOS transmitter has been developed where four separate intermediate frequency (IF) signals are up-converted, amplified, and then combined to generate an 8-tone RF output; a broadband 90 nm-CMOS receiver has also been constructed, which consists of an RF-low noise amplifier (LNA), mixer, and IF amplifier. The circuits and their corresponding modules will soon be deployed on the DIII-D and NSTX-U fusion devices. A 110–140 GHz 65 nm-CMOS receiver has also been designed, which is suitable for measuring the deep-core temperature fluctuations in the DIII-D tokamak using the electron cyclotron emission imaging diagnostic system. In addition to the RF-LNA/balun, mixer, and IF amplifier, an LO balun/tripler and driving amplifier are now included in this highly integrated circuit chip. By adopting the microwave and millimeterwave system-on-chip concept in the front-end system design, this paper demonstrates that compact transmitter and receiver modules can be easily built, which, in turn, facilitates array implementation and maintenance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Alfvénic modes excited by the kink instability in PHASMA

Magnetic flux ropes have been successfully created with plasma guns in the newly commissioned PHAse Space MApping (PHASMA) experiment. The flux ropes exhibit the expected m = 1 kink instability. The observed threshold current for the onset of this kink instability is half of the Kruskal–Shafranov current limit, consistent with predictions for the non-line tied boundary condition of PHASMA. The helicity, paramagnetism, and growth rate of the observed magnetic fluctuations are also consistent with kink instability predictions. The observed fluctuation frequency appears to be a superposition of a real frequency due to a Doppler shift of the kink mode arising from plasma flow (~2kHz) and a contribution from a wave mode (~5 kHz). The dispersion of the wave mode is consistent with an Alfvén wave. Distinct from most previ-ous laboratory studies of flux ropes, the working gas in PHASMA is argon. Thus, the ion cyclotron frequency in PHASMA is quite low and the frequency of the Alfvénic mode plateaus at ~0.5 of the ion gyro frequency with increasing background magnetic field strength.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

New, improved analysis of correlation ECE data to accurately determine turbulent electron temperature spectra and magnitudes (invited)

Turbulent electron temperature fluctuation measurement using a correlation electron cyclotron emission (CECE) radiometer has become an important diagnostic for studying energy transport in fusion plasmas, and its use is widespread in tokamaks (DIII-D, ASDEX Upgrade, Alcator C-Mod, Tore Supra, EAST, TCV, HL-2A, etc.). The CECE diagnostic typically performs correlation analysis between two closely spaced (within the turbulent correlation length) ECE channels that are dominated by uncorrelated thermal noise emission. This allows electron temperature fluctuations embedded in the thermal noise to be revealed and fluctuation level and spectra determined. We have demonstrated a new, improved CECE coherency-based analysis for calculating the temperature fluctuation frequency spectrum and level, which has been verified both numerically through the simulation of synthetic ECE radiometer data and through analysis of experimental data from the CECE system on DIII-D. The new formulation places coherency-based analysis on a firm foundational footing and corrects some currently published methodologies. This new method accurately accounts for bias error in the coherence function and correctly calculates noise levels for a fixed data record length. It provides excellent accuracy in determining temperature fluctuation level (e.g., <10% error) even for a small realization number in the ensemble average. The method also has a smaller uncertainty (i.e., error bar) in the power spectrum when compared to the more standard cross-power method when evaluated at low coherency. Direct calculation of system noise level using correlation between randomized intermediate frequency signals is recommended.

Wang, G. (ORCID:0000000225739827)↗
DOE OSTI2021

Energy-angle distribution of the ions in the RF sheath of ICRH antennas

Radio-frequency sheaths forming at the surface of ion cyclotron resonance heating (ICRH) actuators in fusion experiments are associated with enhanced impurity sputtering from ICRH plasma-facing components (PFCs). The minimization of impurity fluxes from the ICRH PFC is a vital task for the usage of ICRH systems. Capturing the ion kinetics is important for an accurate description of plasma–material interactions, because the ion dynamics plays a crucial role in RF sheaths. Here, we present a hybrid particle-in-cell (hPIC) model able to capture the kinetic behavior of the ions. We analyze the kinetic ion energy-angle distributions (IEADs) impacting the RF antenna and its dependence on different plasma and RF sheath parameters. In particular, the IEAD dependence on RF frequency and magnetic field alignment is investigated. Using hPIC, we simulated a case emulating the latest experimental campaign from JET. The simulation showed that under specific plasma and RF parameters, the kinetic motion of ions results in a cusp formation in the IEAD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗