Search NASA⌕ Search

SEARCH · Search NASA

Results for “RF simulation”

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

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗
DOE OSTI2024

Analysis of fusion alphas interaction with RF waves in D-T plasma at JET

This work studies the influence of RF waves in ICRH range of frequency on fusion alphas during the recent JET D-T campaign. Fusion alphas from D-T reactions are born with energies of about 3.5MeV and therefore have significant Doppler shift enabling synergistic interaction between them and RF waves at broad range of frequencies including the ones foreseen for future fusion machines ITER and SPARC. Resonant interaction between RF waves and alphas, also called synergistic effects, will modify the alpha distribution and ultimately will have an impact on alpha orbit losses and heating. Data from JET 3.43T/2.3MA pulses based on the hybrid scenario during the DTE2 campaign were used for the analysis in this study. The impact of synergistic effects on alpha orbit losses and alpha heating is assessed. Conclusions are based on analysis of experimental data for fast alphas losses, i.e. measurements from neutral particle analyser, fast ion losses scintillator detector, Faraday cups, and TRANSP simulations. Experimental data and TRANSP analysis indicate that there are indeed changes in the alphas' distribution function due to interaction with RF waves. Data from the scintillator detector and the Faraday cups were compared for pulses with and without ICRH power and versus cases with enhanced alpha losses due to MHD activities. The trends from these diagnostics consistently show no additional alpha losses due to interaction with RF waves. TRANSP predictions for the impact of the synergistic effects on alpha heating show up to 42% increase in alpha electron heating and up to 25% increase in alpha ion heating. These effects however become negligibly small, less than 1%, when alpha heating is compared to the total auxiliary hearting power in the investigated JET pulses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Radio Frequency Algorithms

Final Technical Report for DOE SC0019287 award titled Radio Frequency Algorithms. This project supported the design, simulation, analysis, and optimization of robust algorithms to configure, optimize, and control Radio Frequency and Low-Level RF systems for several modern accelerator systems.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Multiple Beam Triode Driven RF Sources for Accelerator Applications Phase I Final Report

Calabazas Creek Research, Inc, (CCR), in collaboration with Microwave Power Products, Inc. (MPP), formerly Communications & Power Industries, LLC (CPI,) and JP Accelerator Works, Inc. (JPAW), embarked on a program to develop multiple beam triodes to produce RF power from 350 – 800 MHz with an average power exceeding 200 kW. The effort was motivated by the performance of a triode-based RF source which produces 25 kW of UHF power at 90% efficiency. The CCR effort focused on implementing this technology into a multiple beam device to increase the output power while retaining the low cost, compact size, and high efficiency. The program performed extensive simulations indicating that the goals could be achieved, and a prototype multiple beam triode was built, baked, and tested. Unfortunately, a grid to cathode short terminated the testing before the tube could generate RF power. Nevertheless, the effort demonstrated that a multiple beam triode could be designed, built, baked, and energized to high voltage. The multiple beam triode used oxide cathodes, which are only capable of pulsed operation. The multiple beam triode will be rebuilt using dispenser cathodes, which will allow high duty or continuous operation. The grids were also modified to be more robust to avoid previous issues. The MB triode will provide the beam power for RF generation. The RF is generated by surrounding the triode with input and output cavities to convert beam power to RF power. RF cavities to generate 200 kW CW at 350-450 MHz using the MB triode with dispenser cathodes was assembled during the program. The next Phase of this effort is to assemble the multiple beam triode using the subassemblies built in the Phase I program and test with the RF cavities. The Phase I program also initiated design of a higher frequency, higher power multiple beam triode. That design is forecast to produce approximately 500 kW CW from 350 - 500 MHz.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Machine Learning for Multipactor Susceptibility Prediction in Planar RF Gaps

Multipactor discharge is a nonlinear electron avalanche that limits the performance of high-power radio-frequency (RF) and vacuum electronic devices. Predicting multipactor susceptibility traditionally relies on Monte Carlo or particle-in-cell (PIC) simulations, which become computationally expensive for large parametric studies. In this work, we present a supervised machine-learning (ML) framework for prediction of multipactor susceptibility in a two-surface planar geometry. The models are trained using high-fidelity PIC simulation generated susceptibility data and learn the relationship between operational parameters, geometry, and material-dependent secondary electron emission properties. The proposed approach enables rapid reconstruction of susceptibility charts while preserving the physical structure of multipactor growth regions.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Study of ponderomotive instability in the FRIB beta=0.53 half-wave resonator

Superconducting radio-frequency niobium cavities are susceptible to deformations caused by external or internal forces, leading to shifts in the cavity resonant frequency. One source of deformation comes from the radiation pressure of the cavity fields, producing the so-called Lorentz force detuning effect. This effect can couple to the cavity mechanical modes in generator driven mode, leading to ponderomotive instability. In the FRIB 322 MHz, β=0.53 Half-wave Resonators (HWR), the instability appeared when the cavity was detuned, with thresholds depending on low-level RF control parameters, such as closed loop gain, as well as accelerating gradient. Using a measured Lorentz Transfer Function, Simulink simulations were conducted to predict the instability thresholds, which were then compared with the experimental results. We will discuss how these thresholds can be broadened to enable stable operation at higher gradients.

Accelerator Physics↗
DOE OSTI2025

Fermilab's controls development with virtual accelerator

Control Systems development is often the last thing considered when designing and building new equipment, e.g. a new detector or superconducting RF LINAC; however when the new equipment is installed, it is the first thing desired to be operational for testing. Due to frequent delays in building new equipment and project deadlines, control system development and testing is often curtailed. A way to alleviate this problem is to simulate the control system, though this will be challenging for complex systems.The Fermilab PIP-II (proton improvement plan - II) project is being constructed at Fermilab to deliver $800\,MeV$ protons of $>1\,MW$ beam power to replace the present LINAC for the remainder of the existing accelerator complex. The new LINAC consists of a warm front end (WFE), 23 superconducting RF cryomodules (of 5 types), and a beam transfer line (BTL) to the existing complex.The accelerator physics group has a parallel project to create a digital twin (DT) of the PIP-II accelerator. We have coupled the EPICS controls to this DT and are developing both the DT and EPICS software in parallel. This will allow us to develop the EPICS software framework, the HMIs, sequences, high level physics applications, and other services for use in a fully functional control system.This presentation will detail the work that we have performed to date and show demonstrations of controlling and monitoring the status of the accelerator, as well as future plans for this work.

Hanlet, Pierrick [Fermilab]↗
DOE OSTI2025

A fully ro-vibrationally resolved corona model for the molecular hydrogen Fulcher-𝛼 system

A ro-vibrationally resolved corona model for the molecular hydrogen Fulcher-𝛼 transition based on the Yacora solver is introduced. The model couples 1365 ro-vibrational levels of the 𝑋 1 ⁢𝛴$^+_𝑔$, 𝑑 3⁢ 𝛱 𝑢 and 𝑎 3 ⁢𝛴$^+_𝑔$ states via electron impact excitation from the ground state into the 𝑑 3⁢ 𝛱 𝑢 state and subsequent spontaneous emission into the 𝑎 3 ⁢𝛴$+_𝑔$ state. For the process of electron impact excitation a set of 45260 fully ro-vibrationally resolved cross sections calculated with the molecular convergent close-coupling (MCCC) method in the adiabatic-nuclei formulation is applied. The MCCC cross sections are compared with results of a scaling method used in other works, demonstrating the need for dedicated ro-vibrationally resolved cross sections. By post-processing the model output, entire Fulcher spectra can be simulated. These spectra are benchmarked with measurements from an inductively coupled plasma discharge at a pressure between 1.1–10 Pa and a RF power of 700 W. The model results agree very well with the experiment both in relative shape and in absolute value. A non-ro-vibrationally resolved, purely electronically resolved collisional radiative model for molecular hydrogen (including further states and reaction channels) is applied to investigate the validity of the corona approximation for the benchmark plasmas and the relevance of possible process extensions of the corona model. Furthermore, the influence of collisional quenching, autoionization and predissociation is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

EIC hadron ring longitudinal stability simulations at injection

At the EIC injection energy (23.8 GeV for protons), the hadron beam experiences a relatively strong space charge field. The corresponding wakefield, in conjunction with the wakes from the RF cavities (i.e. the geometrical wakes), may become the cause of single-bunch instabilities, both longitudinal and transverse. In this report, we focus on the longitudinal stability aspect of this picture, using a simplified model for the ring.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Integration Development and Testing of Rear Transition Monitor for Beam Current Monitoring System

Addressing baseline effects in accelerator environments is crucial for accurate data acquisition and analysis, since baseline effects can obscure signal clarity and impact the reliability of beam current monitoring systems. There are many potential contributors to baseline noise, such as variations in beam dynamics, electromagnetic interference from nearby equipment, or RF interference. Previous applications of noise reduction systems don t sufficiently filter sources of asynchronous noise, so a new algorithm was implemented. A simulation dataset was created to replicate beam conditions and a Red Pitaya FPGA was used to collect data through the streaming application. A Python script was developed to implement noise reduction algorithms and efforts were made to integrate real-time data streaming with the Redis platform and Acnet Front End infrastructure.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Technical Report for Bayesian Optimization and Reinforcement Learning for Beam Polarization Increase in the BNL Hadron Injectors

This project developed and evaluated physics-informed Bayesian learning and machine learning (ML)-based optimization methods for improving beam polarization preservation in the BNL hadron injector chain. The work focused on uncertainty-aware digital twin modeling, Bayesian calibration of accelerator simulations using beam measurements, and data-efficient optimization strategies including Bayesian optimization and reinforcement learning. These methods were applied to injector tuning and RF control problems in realistic accelerator settings to support improved operational robustness and readiness for RHIC operations and future Electron–Ion Collider facilities. No subject inventions were disclosed under this award.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Convergence study of wakefield simulations with GdfidL and ECHO3D

The interaction of charged particle beams with vacuum chamber components gives rise to electromagnetic wakefields, whose frequency-domain representation is known as beam coupling impedance. Geometric impedance arising from discontinuities and transitions in the vacuum chamber is the focus of this study. Minimizing this impedance is essential to mitigate adverse collective effects in modern storage rings operating with high-intensity particle beams. Accurate and reliable impedance simulations is a key factor of the vacuum chamber design. This paper presents the results of a convergence study of two widely used electromagnetic solvers, GdfidL and ECHO3D, applied to key vacuum-chamber components of the National Synchrotron Light Source II (NSLS-II) storage ring. Detailed comparisons are performed for several geometries, including flange absorbers, RF bellows, button-type beam position monitors, and an in-vacuum undulator (IVU). The results show notable differences in convergence and computational efficiency between the two codes. While GdfidL provides highly resolved results and serves as a common benchmark tool, ECHO3D yields consistent results with coarser meshes, significantly reducing simulation time and memory demands. Simulations with a full-geometry IVU model demonstrate that simplified taper-transition models can miss important impedance contributions. In conclusion, these findings provide practical guidelines for efficient and accurate impedance modeling to optimize design of vacuum chamber components for accelerators.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Broadband cross polarization for ultra-wideline magic-angle spinning NMR

Over the past decade, there has been a sustained interest in using frequency-swept (FS) pulses for the efficient acquisition of wideline and ultra-wideline (UW) NMR powder patterns. Such experiments are typically conducted under static conditions, employing both direct- and indirect-excitation methods (i.e., WCPMG and BRAIN-CP/WCPMG, respectively). Recently, Koppe et al. demonstrated that the WCPMG pulse sequence can be used to efficiently acquire wideline and UW NMR spectra with spinning sideband (SSB) manifolds under magic-angle spinning (MAS) conditions, capitalizing on the increased signal-to-noise ratios (SNR) afforded by MAS. To date, there have been only a few instances of broadband cross-polarization (CP) experiments using FS pulses under MAS conditions and no applications to systems exhibiting wideline and/or ultra-wideline powder patterns, despite the clear advantages these experiments could offer. Herein, we demonstrate that FS pulses selectively applied to a single sideband of the S spin can be used for efficient 1 H-S polarization transfer to S = 1/2 nuclides with large anisotropic chemical shift interactions at slow to moderate MAS rates. The Hartmann–Hahn matching conditions in BRAIN-CP/WCPMG-MAS experiments bear similarity to those of standard CP sequences, yet operate over UW frequency ranges and only require low-amplitude RF pulses on the S channel. Crucial to the success of the BRAIN-CP/WCPMG-MAS experiment is careful calibration of the RF amplitude, transmitter offset, and effective frequency sweep of the FS pulse applied to the S spins at a given MAS rate. Thus, by means of numerical simulations and experimental testing, we provide recommendations for the parameterization and setup of BRAIN-CP/WCPMG-MAS experiments for their most efficient use. Results showcasing the capability of the BRAIN-CP/WCPMG-MAS pulse sequence are presented, including applications to 119 Sn, 195 Pt, and 103 Rh NMR.

Kimball, James J. [Florida State University, Talla↗
DOE OSTI2026

Cryogenic-Refined MOSFET Modeling for Oscillator, Frequency Divider, and Amplifier Designs Below 4 K

Capturing device characteristic changes at cryogenic temperatures is crucial for cryo-CMOS circuit designs. In this work, we present an isothermal cryogenic-refined modeling approach for CMOS transistors that is simple, low overhead, and easy to implement while offering the required accuracy for predicting circuit performance at the designated temperatures. Guided by die-level measurement data and circuit design principles, the model introduces corrections to only five critical parameters: threshold voltage, carrier mobility, elevated low-frequency flicker noise, dominant high-frequency shot noise, and subthreshold swing (SS). These refinements are implemented around the foundry-provided SPICE model, which is typically validated only down to about 200 K. With these adjustments, the proposed cryogenic-refined model achieves less than 5% error in both large-signal metrics (I–V characteristics) and small-signal parameters (e.g., transconductance) when compared with device measurements at deep-cryogenic temperatures. The methodology is validated in two advanced technologies: TSMC 40-nm CMOS and GlobalFoundries (GF) 45-nm RF-SOI. We further demonstrate its applicability in three representative RF circuits: a 30-GHz LC oscillator, a high-speed current-mode-logic (CML) frequency divider (FD), and a subthreshold Gb/s amplifier, all showing close agreement between simulated predictions and measurements performed at 4 and 2.5 K. Finally, we believe that the proposed approach is implementation-friendly and can significantly accelerate the development of cryo-CMOS integrated circuits.

circuit modeling↗
DOE OSTI2025

Results and lessons learned from accelerating radio frequency modeling using machine learning [slides]

The “advanced tokamak” reactor concept is a leading candidate for a steady state fusion pilot plant. An advanced tokamak (AT) sustains a majority of the required plasma current with effects resulting from maintenance of the peaked pressure at the device center. This current is augmented by auxiliary current drive sources. These auxiliary actuators may consist of neutral particle beams and/or radio frequency (RF) systems such as lower hybrid current drive (LHCD) and high harmonic fast wave (HHFW) current drive using radio and microwaves from antennas. The primary focus of this work is to develop models of RF current profile control suitable for use in integrated modeling frameworks and for real-time control in experiments. Direct physics models of RF current drive can be computationally intensive. In order to achieve predictive times appropriate for the thousands of calls needed in real-time control of experiments and for use in integrated models, we will apply modern machine learning (ML) techniques to accelerate these models and interpolate their results. To generate the fast and accurate models for use in control level algorithms and integrated modeling we need to replace present models with high dimensional interpolation of their results. We will perform additional simulations across a broader parameter range for EAST and other tokamaks in different physics regimes (Alcator C-Mod, DIII-D, WEST, CFETR, ARC, ITER) and combine them into a larger database for training and testing of the ML models. Further testing of the control level models with experimental current profile data from EAST and C-Mod tokamaks will provide additional confirmation of the control level model before integration in a tokamak control system or integrated modeling suite. ML will be used to optimize the selection of training data consisting of RF current driven at different values of density profile, temperature profile, plasma current, and wavenumber. ML will also be used to facilitate classification of current drive from these input data. The output of this effort will be a validated classifier capable of determining the current drive profiles for HHFW CD and LHCD on a mille-second timescale. This will provide a breakthrough capability enabling real-time control of RF driven current profiles in experiments including ITER ICRF and use integrated modeling frameworks requiring thousands of current profile calculations in discharge simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Diagnosing Ghost Bunches with the Upstream Extinction Monitor in the Mu2e Experiment

The Mu2e experiment has a stringent requirement for extinction of the pulsed proton beam, referring to the elimination of particles between proton bunches to a relative level of 10$^-10$, which means a single out-of-time particle in the inter-pulse gaps for every 250 complete proton pulses. As the construction of the Mu2e experiment nears completion, it is crucially important to make an early measurement of the beam extinction in its current condition. Hence the upstream extinction monitor was constructed and operated to probe for problems in the proton pulse structure or a higher than expected incidence rate of out-of-time particles. The analysis in this work comes from data taken in March 2026. The long data run showed a significant presence of out-of-time particles from ghost bunches in the Delivery Ring approximately 388 ns after the centers of the main proton pulses. These are hypothesized to be the result of a combination of a RF frequency mismatch, particle space charge, and machine impedance during the rebunching sequence in the Recycler Ring, which can lead to particles leaking into adjacent buckets, but further studies and simulations are needed to confirm this

Hensley, R. [UC, Davis] (ORCID:0000000216064485)↗
DOE OSTI2024

A high-frequency, low-power resonant radio-frequency neutron spin flipper for high-resolution spectroscopy

Here, we present a resonant-mode, transverse-field, radio-frequency (rf) neutron spin flipper design that uses high-temperature superconducting films to ensure sharp transitions between uniform magnetic field regions. Resonant mode allows for low-power, high-frequency operation but requires strict homogeneity of the magnetic fields inside the device. This design was found to efficiently flip neutrons at 96.6 ± 0.6% at an effective frequency of 4 MHz in bootstrap configuration with a beam size of 2.4 × 2.5 cm 2 and a wavelength of 0.4 nm. The high frequency and efficiency enable this device to perform high-resolution neutron spectroscopy with comparable performance with currently implemented rf flipper designs. The limitation of the maximum frequency was found due to the field homogeneity of the device. We numerically analyze the maximum possible efficiency of this design using a Bloch solver simulation with magnetic fields generated from finite-element simulations. We also discuss future improvements of the efficiency and frequency to the design based on the experimental and simulation results.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Two-Frequency RF Fields Induced Multipactor in Coaxial Transmission Lines

Multipactor is a nonlinear discharge phenomenon that occurs in vacuum RF systems, potentially leading to signal distortion, power loss, and even permanent damage to high-power components. This study presents a detailed investigation of two-surface multipactor in coaxial transmission lines under two-frequency excitation, using one-dimensional (1D) Monte Carlo simulations validated by three-dimensional (3D) Particle-in-Cell (PIC) results and experimental data. Introducing a second carrier mode is shown to suppress multipactor by reshaping and shrinking the susceptibility region, with the extent and location of suppression strongly dependent on the device aspect ratio and the relative phase of the second mode. Distinct suppression patterns emerge across different frequency–gap distance (fd) regimes, and in certain cases, susceptibility expansion is also observed. The study identifies and distinguishes pure and mixed multipactor modes in coaxial geometry, where analytical mode boundaries are not readily defined. Unlike planar systems, pure-mode regions in coaxial structures overlap with mixed-mode domains, complicating classification. Image charge forces are found to have minimal effect on susceptibility thresholds but do influence electron growth rates. These findings provide new insights into waveform-driven control of multipactor in high-power RF systems.

43 PARTICLE ACCELERATORS↗