Search NASA⌕ Search

SEARCH · Search NASA

Results for “radio frequency heating”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Role of microstructure on flux expulsion of superconducting radio frequency cavities

The trapped residual magnetic flux during the cool-down due to the incomplete Meissner state is a significant source of radio frequency losses in superconducting radio frequency cavities. Here, in this study, we clearly correlate the niobium microstructure in elliptical cavity geometry and flux expulsion behavior. In particular, a traditionally fabricated Nb cavity half-cell from an annealed poly-crystalline Nb sheet after an 800 °C heat treatment leads to a bi-modal microstructure that ties in with flux trapping and inefficient flux expulsion. This non-uniform microstructure is related to varying strain profiles along the cavity shape. A novel approach to prevent this non-uniform microstructure is presented by fabricating a 1.3 GHz single cell Nb cavity with a cold-worked sheet and subsequent heat treatment leading to better flux expulsion after 800 °C/3 h. Microstructural evolution by electron backscattered diffraction-orientation imaging microscopy on cavity cutouts, and flux pinning behavior by dc-magnetization on coupon samples confirms a reduction in flux pinning centers with increased heat treatment temperature. The heat treatment temperature-dependent mechanical properties and thermal conductivity are reported. The significant impact of cold work in this study demonstrates clear evidence for the importance of the microstructure required for high-performance superconducting cavities with reduced losses caused by magnetic flux trapping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

The 2022 High-energy Outburst and Radio Disappearing Act of the Magnetar 1E 1547.0–5408

We report the radio and high-energy properties of a new outburst from the radio-loud magnetar 1E 1547.0–5408. Following the detection of a short burst from the source with Swift-BAT on 2022 April 7, observations by NICER detected an increased flux peaking at (6.0 ± 0.4) × 10 –11 erg s –1 cm –2 in the soft X-ray band, falling to a baseline level of 1.7 × 10 –11 erg s –1 cm –2 over a 17 day period. Joint spectroscopic measurements by NICER and NuSTAR indicated no change in the hard nonthermal tail despite the prominent increase in soft X-rays. Observations at radio wavelengths with Murriyang, the 64 m Parkes radio telescope, revealed that the persistent radio emission from the magnetar disappeared at least 22 days prior to the initial Swift-BAT detection and was redetected two weeks later. Such behavior is unprecedented in a radio-loud magnetar, and may point to an unnoticed slow rise in the high-energy activity prior to the detected short bursts. Finally, our combined radio and X-ray timing revealed the outburst coincided with a spin-up glitch, where the spin frequency and spin-down rate increased by 0.2 ± 0.1 μHz and (–2.4 ± 0.1) × 10 –12 s –2 , respectively. A linear increase in the spin-down rate of (–2.0 ± 0.1) × 10 –19 s –3 was also observed over 147 days of postoutburst timing. Our results suggest that the outburst may have been associated with a reconfiguration of the quasi-polar field lines, likely signaling a changing twist, accompanied by spatially broader heating of the surface and a brief quenching of the radio signal, yet without any measurable impact on the hard X-ray properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Investigation of Materials for Radio Frequency Antenna Plasma Facing Components

Here, the interaction of radio frequency (RF) sheaths with fusion reactor relevant materials (e.g., tungsten and titanium diboride) is being studied on the RF Plasma Interaction Experiment (RF PIE). The RF PIE consists of an electron cyclotron resonance (ECR) plasma source (2.45 GHz, 5 kW) with a biased and heated RF electrode that is used to simulate antenna surfaces in contact with the edge plasma. Helium plasmas (density of ~1e18/m 3 , electron temperature of 4–5 eV) are being used to explore sheath formation on material surfaces with biases up to 500 V. The erosion of a tungsten surface is being studied spectroscopically using a mirror-linked 1 m Czerny-Turner UV imaging spectrometer with a spectral resolution of 0.012 nm for measuring plasma emission in and near the sheath. Tungsten line emission intensity is higher for RF versus dc biasing for similar plasma conditions and average ion energy. RF biasing causes a broadening of the ion energy distribution function (IEDF) due to the RF sheath, as determined from the hPIC2 code, and results in enhanced sputtering. Calculations of the expected sputtering yield for dc and RF biasing are consistent with experimental observations of changes in the 400.9 nm tungsten line emission intensity as a function of ion energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The multiwavelength view of shocks in the fastest nova V1674 Her

Classical novae are shock-powered multiwavelength transients triggered by a thermonuclear runaway on an accreting white dwarf. V1674 Her is the fastest nova ever recorded (time to declined by two magnitudes is t 2 = 1.1 d) that challenges our understanding of shock formation in novae. We investigate the physical mechanisms behind nova emission from GeV γ-rays to cm-band radio using coordinated Fermi-LAT, NuSTAR, Swift, and VLA observations supported by optical photometry. Fermi-LAT detected short-lived (18 h) 0.1–100 GeV emission from V1674 Her that appeared 6 h after the eruption began; this was at a level of (1.6 ± 0.4) × 10 −6 photons cm −2 s −1 . Eleven days later, simultaneous NuSTAR and Swift X-ray observations revealed optically thin thermal plasma shock-heated to k $T$ shock = 4 keV. The lack of a detectable 6.7 keV Fe Kα emission suggests super-solar CNO abundances. The radio emission from V1674 Her was consistent with thermal emission at early times and synchrotron at late times. The radio spectrum steeply rising with frequency may be a result of either free-free absorption of synchrotron and thermal emission by unshocked outer regions of the nova shell or the Razin–Tsytovich effect attenuating synchrotron emission in dense plasma. The development of the shock inside the ejecta is unaffected by the extraordinarily rapid evolution and the intermediate polar host of this nova.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Evaluation of predictive correlation between flux expulsion and grain growth for superconducting radio frequency cavities

Abstract A series of experiments were carried out in an effort to develop a simple method for predicting the magnetic flux expulsion behavior of high purity niobium used to fabricate superconducting radio frequency (SRF) cavities. Using conventional metallographic characterizations in conjunction with high spatial resolution electron backscattered diffraction-orientation imaging microscopy, we find that the flux expulsion behavior of 1.3 GHz single cell SRF Nb cavities is significantly associated with the grain growth of the Nb material during heat treatment. Most of the Nb grains rapidly grew during 900 °C heat treatment, and were likely full-recrystallized with 1000 °C HT. When compared to the magnetic flux expulsion ratio ( B sc / B nc ) at Δ T = 5 K, the flux expulsion efficiency of the cavities increases with the increase of grain size. Most interestingly, 900 °C HT shows a roughly linear trend that suggests that this criterion could be used to predict the appropriate heat treatment temperature for sufficient flux expulsion behavior in SRF-grade Nb. This result can be used to see if flux expulsion can be predicted by examining the materials coming from Nb vendors, prior to cavity fabrication.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Investigating the effects of electron bounce-cyclotron resonance on plasma dynamics in capacitive discharges operated in the presence of a weak transverse magnetic field

Recently, Patil et al. [Phys. Rev. Res. 4, 013059 (2022)] have reported the existence of an enhanced operating regime when a low-pressure (5 mTorr) capacitively coupled discharge (CCP) is driven by a very high radio frequency (60 MHz) source in the presence of a weak external magnetic field applied parallel to its electrodes. Their particle-in-cell simulations show that a significantly higher bulk plasma density and ion flux can be achieved at the electrode when the electron cyclotron frequency equals half of the applied radio frequency for a given fixed voltage. In the present work, we take a detailed look at this phenomenon and further delineate the effect of this “electron bounce-cyclotron resonance (EBCR)” on the electron and ion dynamics of the system. We find that the ionization collision rate and stochastic heating are maximum under resonance condition. The electron energy distribution function also indicates that the population of tail-end electrons is highest for the case where EBCR is maximum. Formation of electric field transients in the bulk plasma region is also seen at lower values of applied magnetic field. Finally, we demonstrate that the EBCR-induced effect is a low-pressure phenomenon and weakens as the neutral gas pressure increases. The potential utility of this effect to advance the operational performance of CCP devices for industrial purposes is discussed in this report.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Jefferson Laboratory C100 Superconducting Radio-Frequency Cavity Fault Data, 2020

The dataset was created to train machine learning models for the task of identifying the (1) cavity and (2) fault type from C100-type cryomodules at the Thomas Jefferson National Accelerator Facility (Jefferson Lab), thereby replacing the time-consuming efforts of a subject matter expert. Superconducting radio-frequency (SRF) cavity trips represent a significant source of accelerator downtime. Real-time – rather than post-mortem – identification of the offending cavity and classification of the fault type would give control room operators valuable feedback for corrective action planning. The anticipated benefit is increased beam-on-target time for users and provides performance metrics that can be used to improve future cavity designs. A series of 17 RF signals are recorded for each of the 8 cavities in a C100 cryomodule every time a cavity trips. These time-series signals are written to file using a specially designed data acquisition system. The dataset represents fault events recorded during Continuous Electron Beam Accelerator Facility (CEBAF) beam operations between January 18, 2019 and March 9, 2020. The following filtering steps were applied to collected data; (1) only 4 of the 17 signals per cavity are retained (GMES, GASK, CRFP, DETA2) (2) only events with data from each of the eight cavities in the cryomodule are kept, (3) only events that were sampled at 5 kHz were kept, (4) events from cryomodule 0L04 were neglected, (5) events occurring between February 4, 3PM and February 5, 12PM were neglected. As a result of preprocessing, the dataset is comprised of 2,375 unique events. The full dataset is comprised of three files: features.csv, cavity_labels.csv, fault_labels.csv. Each instance in faults.csv includes a timestamp (“date_time”), a label for the cryomodule which experienced the trip (“zone_label”), and 192 features (“feature_1”, “feature_2”... “feature_192”). The features correspond to 6 autoregressive features for each of 4 signals per cavity for each of the 8 cavities (6 × 4 signals/cavity × 8 cavities/cryomodule = 192). To deal with the large variation of signal amplitudes, time-series standardization via the z-score (standard score) function was applied prior to computing the features. For each instance, there is an associated label for the (1) cavity which faulted first (cavity_labels.csv) and (2) the type of fault that caused the trip (fault_labels.csv). The cavity identification can take values of [0, 1, 2, 3, 4, 5, 6, 7, 8] and the fault type can take values of [‘Microphonics’, ‘Quench_100ms’, ‘Controls_Fault’, ‘E_Quench’, ‘Quench_3ms’, ‘Single_Cav_Turn_Off’ , ‘Heat_Riser_Choke’, ‘Multi_Cav_Turn_Off’].

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Design and Prototyping of the 162.5 MHz PIP-II RF Reference-Line Station

The Proton Improvement Plan II (PIP-II) Reference Line distributes phase-stable radio-frequency signals so accelerating cavities transfer energy to the proton beam at the correct point in each RF cycle. The full system contains 162.5, 325, and 650 MHz sections; this project focuses on the 162.5 MHz station. Commercial filters, amplifiers, splitters, mixers, and hardware were assembled into a compact, serviceable heat-plate layout and tested to verify performance metrics at 162.5 MHz.

Subedi, Harsheet [Unlisted, US]↗
DOE OSTI2022

Dehydration of gas-phase benzyl amine alcohols studied at atmospheric pressure

Gas-phase ion chemistry is a useful approach for the investigation of physical organic chemistry and the study of reactivity, structure, and thermochemical properties of ionic and neutral organic substrates. Herein, one-step dehydration of gas-phase benzyl amine alcohols to give benzazetidines was discovered without the use of catalysts. Mechanistic investigations of the gas-phase dehydration reactions were explored to study the possible influence of charged microdroplet acceleration during electrospray ionization (ESI) and collisions occurring during ion transfer through the atmospheric pressure interface (API) of the mass spectrometer. The product ion distribution was observed to be less sensitive to droplet effects (spray distance and voltage), and temperature of transfer capillary in the API. However, the product ion distribution exhibit high sensitivity to different S-lens radio frequency (RF) voltage of the API, providing evidence that the dehydration reaction is driven by collisions during ion transfer in the API. With this insight, we developed atmospheric pressure thermal dissociation platform that allowed thermal-induced collisions outside of the mass spectrometer when the benzyl amine alcohols are carried through a heated coiled tube. Up to 98% dissociation efficiency was achieved with the coiled tube heated to a temperature of 200 °C, yielding only the desired benzazetidine product. Finally, by performing the atmospheric pressure thermal dissociation experiment outside of the mass spectrometer, we believe it will be straightforward to collect the dehydration product.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

First-principles study and experimental characterization of metal incorporation in germanium telluride

Germanium telluride is a well-known phase change material (PCM) used in non-volatile memory cells and radio frequency switches. Controlling the properties of GeTe for improved PCM device performance has sometimes been achieved by doping and/or alloying with metals, often at concentrations greater than 10 at. % and using non-equilibrium methods. Here, since switching PCMs between the low-resistance crystalline and high-resistance amorphous states requires a heating cycle, the stability of metal-incorporated GeTe (Ge 0.5–x M x Te 0.5 ) films is also critical to practical implementation of these materials in electronic and optoelectronic devices. In this work, we use both density-functional theory and experimental characterization methods to probe the solubility and critical properties of Ge 0.5–x M x Te 0.5 films. Using first-principles calculations, we determine the enthalpy of formation for GeTe with 2.08, 4.17, and 6.25 at. % of Cu, Fe, Mn, Mo, and Ti and show trends between the stability of the Ge 0.5–x M x Te 0.5 systems and the atomic position, composition, and distribution of the metal atoms in the GeTe matrix.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A ruthenium oxide thermometer for dilution refrigerators operating down to 5 mK

At the lowest temperatures achieved in dilution refrigerators, ruthenium oxide resistance thermometers often saturate and therefore lose their sensitivity. In an effort to extend the range of such temperature sensors, we built a thermometer which maintains sensitivity to 5 mK. A key component of this thermometer is an in situ radio frequency filter which is based on a modern rf absorption material. Here, we show that the use of such a filter is only effective when it is encased in the same rf-tight enclosure as the ruthenium oxide sensor. Our design delivers an attenuation level that is necessary to mitigate the effects of parasitic heating of a fraction of pW present in our circuit. Furthermore, we show that the likely origin of this parasitic heating is the black body radiation present within the experimental space of the refrigerator. We found that the equilibration time of the thermometer increases rapidly as the temperature is lowered; below 5 mK this thermometer becomes impractical because of the prohibitively long equilibration times.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Fluid-kinetic modeling of a high power density radio frequency inductively coupled positive hydrogen ion source

High power density radio-frequency (RF) inductively coupled positive ion sources are attractive candidates for next-generation neutral beam injection (NBI) systems, where higher injected power and longer pulse lengths are desired without sacrificing source reliability. Operating at absorbed power densities of order $\gt 1~\mathrm{W\,cm}^{-3}$ places these sources in a regime with stronger gas heating, higher dissociation, and non-Maxwellian electron energy distributions. The Large Uniform Plasma for Ionizing Neutrals (LUPIN) is an RF inductively coupled plasma source designed to explore this high power density regime and to provide guidance for a positive ion source upgrade for the DIII-D NBI system. LUPIN is designed to operate at up to 20 kW of RF power at 2 MHz, coupling energy through a cylindrical quartz vessel to achieve target ion current densities of $2100\,\mathrm{A\,m}^{-2}$ . This paper presents fluid-kinetic modeling of LUPIN using the hybrid plasma equipment model where electrons are treated kinetically, and the simulations reveal that electron energy distribution function transitions from nearly Maxwellian in the core to bi-Maxwellian towards the edge. Parametric simulations investigate the effects of RF power, gas pressure, and frequency on plasma density, ion flux, and uniformity. Parametric sweeps reveal that increasing power shifts the primary ionization channel from molecular to atomic with diminishing flux gains due to skin-depth contraction and gas rarefaction. Higher frequency localizes heating and increases $\mathrm{H}_2^+$ and $\mathrm{H}_3^+$ delivery to the grid, while elevated pressure boosts ionization yet hinders ion transport due to increase in collisionality.

inductively coupled plasma↗
DOE OSTI2023

Rotating spokes, potential hump and modulated ionization in radio frequency magnetron discharges

Here, in this work, the transition from the gradient drift instability (GDI) into an m = 1 rotating spoke in the radio frequency magnetron discharge was studied by means of the two-dimensional axial-azimuthal (z–y) particle-in-cell/Monte Carlo collision method. The kinetic model combined with the linear analysis of the perturbation revealed that the cathode sheath (axial) electric field E z triggers the GDI, deforming the local potential until the instability condition is not fulfilled and the fluctuation growth stops in which moment the instability becomes saturated. The potential deformation consequently leads to the formation of the potential hump, surrounding which the azimuthal electric field E y is present. The saturation level of E y was found to be synchronized with and proportional to the time-changing voltage applied at the cathode, resulting in the RF-modulation of the electron heating in the E y due to $\nabla B$ drift. In the saturated stage of the instability, it was shown that the rotation velocity and direction of the spoke present in the simulations agree well with the experimental observation (Panjan 2019 J. Appl. Phys.125 203303). In the instability linear stage, the instability mode wavelength and the growth rate were also found to be in good agreement with the prediction of the GDI linear fluid theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

S-band acoustoelectric amplifier in an InGaAs-AlScN-SiC architecture

Here, we report on an acoustoelectric slab waveguide heterostructure for phonon amplification using a thin Al 0.58 Sc 0.42 N film grown directly on a 4H-SiC substrate with an ultra-thin In 0.53 Ga 0.47 As epitaxial film heterogeneously integrated onto the surface of the Al 0.58 Sc 0.42 N. The aluminum scandium nitride film grown directly on silicon carbide enables a thin (∼850 nm thick) piezoelectric film to be deposited on a thermally conductive bulk substrate (370 W/m K for 4H-SiC); the high thermal conductivity of the substrate, large mobility of the semiconductor (∼7000 cm 2 /V s), and low carrier concentration (∼5 × 10 15 cm −3 ) yield low self-heating. A Sezawa mode with optimal overlap between the peak of its evanescent electric field and the semiconductor charge carriers is supported. The high velocity of the heterostructure materials allows us to operate the Sezawa mode amplifier at 3.05 GHz, demonstrating a gain of 500 dB/cm (40 dB in 800 μm). Additionally, a terminal end-to-end radio frequency gain of 7.7 dB and a nonreciprocal transmission of 52.6 dB are achieved with a dissipated DC power of 2.3 mW. In conclusion, the power added efficiency and acoustic noise figure are also characterized.

Hackett, L.↗
DOE OSTI2022

Star formation and AGN feedback in the local Universe: Combining LOFAR and MaNGA

The effect of active galactic nuclei (AGN) on their host galaxies – in particular their levels of star formation – remains one of the key outstanding questions of galaxy evolution. Successful cosmological models of galaxy evolution require a fraction of energy released by an AGN to be redistributed into the interstellar medium to reproduce the observed stellar mass and luminosity function and to prevent the formation of over-massive galaxies. Observations have confirmed that the radio-AGN population is energetically capable of heating and redistributing gas at all phases, however, direct evidence of AGN enhancing or quenching star formation remains rare. With modern, deep radio surveys and large integral field spectroscopy (IFS) surveys, we can detect fainter synchrotron emission from AGN jets and accurately probe the star-forming properties of galaxies, respectively. Here, in this paper, we combine data from the LOw Frequency ARray Two-meter Sky Survey (LoTSS) with data from one of the largest optical IFS surveys, Mapping Nearby Galaxies at Apache Point Observatory (MaNGA), to probe the star-forming properties of 307 local (z < 0.15) galaxies that host radio-detected AGN (RDAGN). We compare our results to a robust control sample of non-active galaxies that each match the stellar mass, redshift, visual morphology, and inclination of a RDAGN host. We find that RDAGN and control galaxies have broad star-formation rate (SFR) distributions, typically lie below the SFMS, and have negative stellar light-weighted age gradients. These results indicate that AGN selected based on their current activity are not responsible for suppressing their host galaxies’ star formation. Rather, our results support the maintenance mode role that radio AGN are expected to have in the local Universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Integrated modeling of RF-induced tungsten erosion at ICRH antenna structures in the WEST tokamak *

This paper introduces STRIPE (Simulated Transport of RF Impurity Production and Emission), an advanced modeling framework developed to analyze material erosion and the global transport of eroded impurities originating from radio-frequency (RF) antenna structures in magnetic confinement fusion devices. STRIPE integrates multiple physics modules: SolEdge3x for scrape-off-layer plasma profiles, COMSOL for 3D RF rectified sheath potentials, RustBCA for erosion yields and surface interactions, and global impurity transport for 3D ion energy-angle distributions and impurity transport. The framework is applied to an ion cyclotron RF-heated L-mode discharge (#57877) in the WEST tokamak, where it predicts a thirty-fold increase in gross tungsten erosion at antenna limiters during the transition from ohmic to ICRH operation. Additionally, under ICRH conditions, a tenfold enhancement in erosion is observed when comparing RF sheath effects to purely thermal sheath conditions. High-charge-state oxygen ions ($\mathrm{O}$ 6+ and above) are identified as the dominant contributors to tungsten sputtering. To validate the model, a synthetic diagnostic tool based on inverse photon efficiency (S/XB coefficients) from the ColRadPy collisional-radiative model enables direct comparison with spectroscopic measurements. Model predictions using a plasma composition of 1% oxygen and 99% deuterium show good agreement with observed W − I (400.9 nm) emission for discharge #57877, supporting the accuracy of the STRIPE framework. This study focuses specifically on gross erosion calculations to demonstrate STRIPE’s capabilities. Future extensions of this work will incorporate net erosion, re-deposition, self-sputtering effects, and whole-device modeling of sputtered tungsten impurity transport. STRIPE is also being applied to other RF-heated linear and toroidal devices, offering valuable insights for antenna design, impurity control, and performance optimization in next-generation fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Fast-ion physics in SPARC

Potential loss of energetic ions including alphas and radio-frequency tail ions due to classical orbit effects and magnetohydrodynamic instabilities (MHD) are central physics issues in the design and experimental physics programme of the SPARC tokamak. The expected loss of fusion alpha power due to ripple-induced transport is computed for the SPARC tokamak design by the ASCOT and SPIRAL orbit-simulation codes, to assess the expected surface heating of plasma-facing components. We find good agreement between the ASCOT and SPIRAL simulation results not only in integrated quantities (fraction of alpha power loss) but also in the spatial, temporal and pitch-angle dependence of the losses. If the toroidal field (TF) coils are well-aligned, the SPARC edge ripple is small (0.15–0.30 %), the computed ripple-induced alpha power loss is small ( ~0.25% ) and the corresponding peak surface power density is acceptable ( 244 kW m -2 ). However, the ripple and ripple-induced losses increase strongly if the TF coils are assumed to suffer increasing magnitudes of misalignment. Surface heat loads may become problematic if the TF coil misalignment approaches the centimetre level. Ripple-induced losses of the energetic ion tail driven by ion cyclotron range of frequency (ICRF) heating are not expected to generate significant wall or limiter heating in the nominal SPARC plasma scenario. Because the expected classical fast-ion losses are small, SPARC will be able to observe and study fast-ion redistribution due to MHD including sawteeth and Alfvén eigenmodes (AEs). SPARC's parameter space for AE physics even at moderate Q is shown to reasonably overlap that of the demonstration power plant ARC (Sorbom et al., Fusion Engng Des., vol. 100, 2015, p. 378), and thus measurements of AE mode amplitude, spectrum and associated fast-ion transport in SPARC would provide relevant guidance about AE behaviour expected in ARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Metaplectic geometrical optics for ray-based modeling of caustics: Theory and algorithms

The optimization of radio frequency-wave (RF) systems for fusion experiments is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at caustics such as cutoffs and focal points, erroneously predicting the wave intensity to be infinite. This is a critical shortcoming of GO, since the caustic wave intensity is often the quantity of interest, e.g., RF heating. Full-wave modeling can be used instead, but the computational cost limits the speed at which such optimizations can be performed. Here, we have developed a less expensive alternative called metaplectic geometrical optics (MGO). Instead of evolving waves in the usual x (coordinate) or k (spectral) representation, MGO uses a mixed X$\equiv$Ax+Bk representation. By continuously adjusting the matrix coefficients A and B along the rays, one can ensure that GO remains valid in the X coordinates without caustic singularities. The caustic-free result is then mapped back onto the original x space using metaplectic transforms. Here, we overview the MGO theory and review algorithms that will aid the development of an MGO-based ray-tracing code. We show how using orthosymplectic transformations leads to considerable simplifications compared to previously published MGO formulas. We also prove explicitly that MGO exactly reproduces standard GO when evaluated far from caustics (an important property that until now has only been inferred from numerical simulations), and we relate MGO to other semiclassical caustic-removal schemes published in the literature. Finally this discussion is then augmented by an explicit comparison of the computed spectrum for a wave bounded between two cutoffs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗