Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Explanability”

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

Applying simulation-based inference to spectral and spatial information from the Galactic Center gamma-ray excess

Abstract The two most favored explanations of theFermiGalactic Center gamma-ray excess (GCE) are millisecond pulsars and self annihilation of the smooth dark matter halo of the galaxy. In order to distinguish between these possibilities, we would like to optimally use all information in the available data, including photon direction and energy information.To date, analyses of the GCE have generally treated directional and energy information separately, or have ignored one or the other completely.Here, we develop a method for analyzing the GCE that relies on simulation-basedinference with neural posterior models to jointly analyze photon directional and spectral information while correctly accounting for the spatial and energy resolution of the telescope, here assumed to be theFermiLarge Area Telescope (LAT).Our results also have implications for analyses of the diffuse gamma-ray background, which we discuss.

Astronomy & Astrophysics↗
DOE OSTI2025

Confronting new NICER mass-radius measurements with phase transition in dense matter and twin compact stars

The (re)analysis of data on the X-ray emitting pulsars PSR J0030+0451 and J0740+6620, as well as new results on PSR J0437-4715 and J1231-1411, are confronted with the predictions of the equation of state (EoS) models allowing for strong first-order phase transition for the mass-radius (M-R) diagram. Here, we use models that are based on a covariant density functional (CDF) EoS for nucleonic matter at low densities and a quark matter EoS, parameterized by the speed of sound, at higher densities. To account for the variations in the ellipses for PSR J0030+0451 obtained from different analyses, we examined three scenarios to assess their consistency with our models, focusing particularly on the potential formation of twin stars. We found that in two scenarios, where the ellipses for PSR J0030+0451 and J0437-4715 with masses close to the canonical mass ∼ 1.4 M ⊙ are significantly separated, our models allow for the presence of twin stars as a natural explanation for potential differences in the radii of these stars.

X-ray binaries↗
DOE OSTI2025

High-significance detection of correlation between the unresolved gamma-ray background and the large-scale cosmic structure

Our understanding of the γ-ray sky has improved dramatically in the past decade, however, the unresolved γ-ray background (UGRB) still has a potential wealth of information about the faintest γ-ray sources pervading the Universe. Statistical cross-correlations with tracers of cosmic structure can indirectly identify the populations that most characterize the γ-ray background. In this study, we analyze the angular correlation between the γ-ray background and the matter distribution in the Universe as traced by gravitational lensing, leveraging more than a decade of observations from the Fermi-Large Area Telescope (LAT) and 3 years of data from the Dark Energy Survey (DES). We detect a correlation at signal-to-noise ratio of 8.9. Most of the statistical significance comes from large scales, demonstrating, for the first time, that a substantial portion of the UGRB aligns with the mass clustering of the Universe as traced by weak lensing. Blazars provide a plausible explanation for this signal, especially if those contributing to the correlation reside in halos of large mass (∼ 10 14 M ⊙ ) and account for approximately 30–40% of the UGRB above 10 GeV. Additionally, we observe a preference for a curved γ-ray energy spectrum, with a log-parabolic shape being favored over a power-law. We also discuss the possibility of modifications to the blazar model and the inclusion of additional γ-ray sources, such as star-forming galaxies, misalinged active galactic nuclei, or particle dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Speculations on the W-mass measurement at CDF*

Abstract The W mass determination at the Tevatron CDF experiment reported a deviation from the SM expectation at the 7 σ level. We discuss a few possible interpretations and their collider implications. We perform electroweak global fits under various frameworks and assumptions. We consider three types of electroweak global fits in the effective-field-theory framework: the S - T , S - T - , and eight-parameter flavor-universal one. We discuss the amounts of tensions between different measurements reflected in these fits and the corresponding shifts in central values of these parameters. With these electroweak fit pictures in hand, we present a few different classes of models and discuss their compatibility with these results. We find that while explaining the discrepancy, the single gauge boson extensions face strong LHC direct search constraints unless the is fermiophobic (leptophobic), which can be realized if extra vector fermions (leptons) mix with the SM fermions (leptons). Vector-like top partners can partially generate the needed shift to the electroweak observables. The compatibility with the top squark is also studied in detail. We find that the non-degenerate top squark soft masses enhance the needed operator coefficients, enabling an allowed explanation compatible with current LHC measurements. Overall, more theoretical and experimental developments are highly in demand to reveal the physics behind this discrepancy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Cause and impact of low-frequency chirping modes in DIII-D hybrid discharges

Significant variations in MHD activity and fast-ion transport are observed in the DIIID high-beta, steady-state hybrid discharges with a mixture of Electron Cyclotron (EC) waves and Neutral Beam Injection (NBI). When Electron Cyclotron Heating (ECH) or Current Drive (ECCD) is applied, Alfven Eigenmodes (AEs) are usually suppressed and replaced by low-frequency bursting modes. The analysis of a recently compiled database of hybrid discharges suggests that the change of the fast-ion pressure especially the perpendicular pressure is the main factor responsible for the instability transition. The lower ratio of fast-ion injection speed v inj to Alfven speed v alfven and slight drop of safety factor q min during ECCD may also facilitate the transition. The database shows that AEs mainly occur when the fast-ion fraction P f /P total is less than 0.53 and v inj /v alfven is greater than 0.50, while low-frequency bursting modes appear in the opposite regime. In this research, P f and P total are the central fast-ion pressure from classical prediction and total plasma pressure, respectively. The reason why the instability transition correlates with P f /P total and v inj /v alfven is that they can significantly modify the drive of low-frequency bursting modes and AEs. The explanation is supported by the observation that low-frequency bursting modes are rarely seen in the hybrids with NBI only, with EC waves and counter-NBI, or with high plasma density. A careful check of the low-frequency bursting modes suggests that they are mainly chirping (neoclassical) tearing modes [referred to as chirping (N)TM], i.e. the mode frequency firstly jumps up from the steady (N)TM frequency, then chirps down, and finally returns to the steady (N)TM frequency. Sometimes, the (N)TMs are fully stabilized and replaced with pure fishbones. Lastly, the resonance condition calculation and ’Kick’ model simulations suggest that (N)TMs and fishbones can interact through modification of the fast ion distribution in phase space, which influences the drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Microtearing modes as the source of magnetic fluctuations in the JET pedestal

Abstract We report on a detailed study of magnetic fluctuations in the JET pedestal, employing basic theoretical considerations, gyrokinetic simulations, and experimental fluctuation data to establish the physical basis for their origin, role, and distinctive characteristics. We demonstrate quantitative agreement between gyrokinetic simulations of microtearing modes (MTMs) and two magnetic frequency bands with corresponding toroidal mode numbers n = 4 and 8. Such disparate fluctuation scales, with substantial gaps between toroidal mode numbers, are commonly observed in pedestal fluctuations. Here we provide a clear explanation, namely the alignment of the relevant rational surfaces (and not others) with the peak in the ω * profile, which is localized in the steep gradient region of the pedestal. We demonstrate that a global treatment is required to capture this effect. Nonlinear simulations suggest that the MTM fluctuations produce experimentally-relevant transport levels and saturate by relaxing the background electron temperature gradient, slightly downshifting the fluctuation frequencies from the linear predictions. Scans in collisionality are compared with a simple MTM dispersion relation. At the experimental points considered, MTM growth rates can either increase or decrease with collision frequency depending on the parameters thus defying any simple characterization of collisionality dependence.

Physics↗
DOE OSTI2021

Heat pulse propagation and anomalous electron heat transport measurements on the optimized stellarator W7-X

The optimized stellarator Wendelstein 7-X (W7-X) is designed to have an approximately quasi-isodynamic magnetic configuration with reduced neoclassical transport in comparison to a classical stellarator, and turbulent transport is expected to be a significant source of anomalous heat transport across the plasma minor radius. The ion temperature gradient driven mode and the trapped electron mode (TEM) are thought to be responsible for the ion-scale turbulence in W7-X plasmas with volume averaged pressure below 1%. In this work, the electron temperature gradient driven turbulence is shown to be a good candidate for the explanation of the observed electron heat flux, in the inner plasma region where the density gradient is weak (in the outer region, a relatively stronger density gradient would drive additional TEM turbulence). Further, the experimental electron heat transport measured during electron cyclotron resonant heating power and plasma density scans is compared to neoclassical predictions, and the stiffness in the electron heat transport measured during transient transport experiments is presented in three common magnetic configurations of W7-X. In low-(β) plasma discharges, the stiffness in the electron heat flux, quantified by the ratio of the heat pulse to power balance diffusivity, ${\chi }_{\mathrm{e}}^{\mathrm{H}\mathrm{P}}/{\chi }_{\mathrm{e}}^{\mathrm{P}\mathrm{B}}$, is measured to be less than 2, and trend downwards with increasing collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Understanding detachment of the W7-X island divertor

The fundamental behavior of the W7-X island divertor under detached conditions, which has been theoretically predicted with the EMC3-Eirene code, is re-examined here under the experimental conditions achieved so far and compared with the first experimental results. Both simulations and experiments cover a range of divertor configurations and plasma parameters, and show the following common trends: (1) with rising impurity radiation, the target heat load decreases 'uniformly' over the entire target surface in the sense that both the peak and average heat loads can drop by an order of magnitude. Impurity radiation (mainly from intrinsic carbon) occurs primarily at the plasma edge and the resulting negative impact on the stored energy is less than 10%. (2) When the total radiation exceeds a critical level, the target particle flux (the recycling flux Γrecy) begins to fall and can drop by a factor of 3–5 at high radiation levels without an obvious indication of significant volume recombination. (3) While Γ recy decreases, the divertor neutral pressure continues to build up and reaches a maximum, at which point Γ recy has declined significantly. (4) During detachment, the electron temperature at the last closed flux surface falls in a way that is not quantitatively understandable from parallel classical heat conduction processes. This paper presents a physical explanation of the numerical/experimental results described above. Furthermore, using the EMC3-Eirene code as a diagnostic tool, we are able, apparently for the first time, to provide a full quantitative analysis of each transport channel in the island divertor, aiming to clarify how the island divertor plasma self-regulates to maintain particle, energy, and momentum balance under detached conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Interstitial hydrogen enhances the mobility of some grain boundaries in tungsten

Segregation of interstitials at a grain boundary (GB) is known to generally lower its mobility. This phenomenon, called 'solute-drag', has important ramifications on the process of recrystallization and microstructural evolution. Here, we present predictions from molecular dynamics (MD) simulations which demonstrate that interstitial hydrogen in tungsten can in fact increase the mobility of some GBs which exhibit shear coupling. Assuming a disconnection-based mechanism, activation energies and pre-factors for disconnection nucleation are predicted from simulations of shear-coupled motion. In GBs where enhanced mobility is predicted, interstitial H reduces both the activation energy and the pre-factor for disconnection nucleation, thus effectively increasing the mobility. For GBs with diminished mobility, MD predicts that presence of interstitial H reduces the pre-factor and, in some cases, increases the activation energy. The reduction in the activation energy inferred from MD simulations are confirmed by nudged elastic band calculations. Temperature-dependent structural transitions are observed for some GBs, and the effect of interstitial H is found to change with the changes in structure. The effect of interstitial H is predicted to be complex and highly variable, providing some plausible explanations for experimental observations on the recrystallization of tungsten in presence of H-loaded plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Simulation study of particle transport by weakly coherent mode in the Alcator C-Mod tokamak

A simulation study has been conducted of the physical mechanisms behind the weakly coherent mode (WCM) and its produced particle transport in the I-mode edge plasmas by using the BOUT++ code. The WCM is identified in our simulations by its poloidal and radial distributions as well as its frequency and wavenumber spectra. Its produced radial particle flux is calculated and compared with the experimental value. The good agreement indicates that the WCM is an important particle transport channel in the I-mode pedestal. It is found that the WCM can transport particles across the strong outer shear layer of the Er well established in the formation of I-mode, based on which a possible explanation is provided why I-mode does not feature a density pedestal. The key point lies in the change of the cross-phase between the electric potential and density fluctuations induced by the E × B Doppler shift. In the strong shear layer, although the electric potential fluctuation is significantly suppressed, the cross-phase is close to π/2, resulting in a strong drive of the density fluctuation and particle transport. To identify the physical nature of the WCM, a linear dispersion relation for drift Alfvén modes is derived in the slab geometry. A drift Alfvén wave instability is found to have similar dependence to the simulated linear instability behind the WCM on the resistivity and the parallel electron pressure gradient and thermal force terms in the parallel Ohm's law.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A reduced model for the ITER divertor based on SOLPS solutions for ITER Q = 10 baseline conditions: A. identifying options for the control parameters*

Edge codes such as SOLPS coupled to neutral codes such as EIRENE have become so comprehensive and sophisticated that they now constitute, in effect, 'code-experiments' that, as for actual experiments, can benefit from interpretation using simple models and conceptual frameworks, i.e. reduced models. The first task is the identification of options for the reduced model control parameters that are best suited for control of the action of the divertor, i.e. for control of target power loading and sputter-erosion, primarily. A strong correlation between the electron temperature at the divertor target, T e,t , and the neutral deuterium D 2 density at the target, n D2,t , flux-tube resolved, has recently been reported for a number of code studies including SOLPS-4.3 modeling of a set of ~50 ITER baseline cases: Q DT = 10, q 95 = 3, P SOL = 100 MW, metallic walls, and Ne seeding. This part A of the present study reports new results for largely the same ITER cases, confirming the strong correlations reported earlier between local values of T e,t , and (i) n D2,t , and (ii) normalized volumetric losses of power and pressure in the divertor. Strong correlations have now also been found, and are reported here for the first time, between T e,t and all of the divertor target quantities of practical interest. A physical explanation for this surprising result has not been fully identified; nevertheless it has encouraging implications for reduced modeling of the ITER divertor. For such ITER conditions, (i) the global Ne injection rate, Inj Ne [Ne/s], and (ii) the electron temperature at the location on the target where the peak power deposition occurs, T e,t @q ⊥,pk (eV), are found to be promising reduced model control parameters. In the companion report, part B, a reduced model for the ITER divertor is developed and described in detail, based on reversed-direction 2 point modelling, Rev2PM. The input to the reduced model is a value of the variable pair (T et @q ⟂,pk ,Inj Ne ) and the output are values of the various target as well as divertor-entrance quantities of practical interest, e.g. q ⊥,pk , n e,Xpt (the electron density at the poloidal location of the X-point), etc. In part B the reduced model is quantitatively characterized using one half of the code cases; it is then used to successfully predict (replicate) the code values of e.g. n e,Xpt for the other half of the cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A reduced model for the ITER divertor based on SOLPS solutions for ITER Q = 10 baseline conditions: B. A reduced model based on reversed-direction two point modeling

Edge codes such as SOLPS coupled to neutral codes such as EIRENE have become so comprehensive and sophisticated that they now constitute, in effect, 'code-experiments' that, as for actual experiments, can benefit from interpretation using simple models and conceptual frameworks, i.e. reduced models. The first task is the identification of options for the reduced model control parameters that are best suited for control of the action of the divertor, i.e. for control of target power loading and sputter–erosion, primarily. A strong correlation between the electron temperature at the divertor target, T e,t , and the neutral deuterium D 2 density at the target, n D2,t , flux-tube resolved, has recently been reported for a number of code studies including SOLPS-4.3 modeling of a set of ~50 ITER baseline cases: Q DT = 10, q 95 = 3, P SOL = 100 MW, metallic walls, and Ne seeding (Pitts et al 2019). Part A of the present study reports new results for largely the same ITER cases, confirming the strong correlation reported earlier between local values of T e,t , and (i) n D2,t , and (ii) normalized volumetric losses of power and pressure in the divertor. Strong correlations have now also been found, and are reported here for the first time, between T e,t and all of the divertor target quantities of practical interest. A physical explanation for this surprising result has not as yet been fully identified; nevertheless it has encouraging implications for reduced modeling of the ITER divertor. For such ITER conditions, (i) the global Ne injection rate, Inj Ne (Ne s –1 ), and (ii) the electron temperature at the location on the target where the peak power deposition occurs, T e,t @q ⊥,pk (eV), are found to be promising reduced model control parameters. In this part B, a reduced model for the ITER divertor is developed and described in detail, based on reversed-direction two point modeling, Rev2PM. The input to the reduced model is a value of the variable pair (T e,t @q ⊥,pk , Inj Ne )for a chosen case and the output are values of the various target as well as divertor-entrance quantities of practical interest, e.g. q ⊥,pk , the electron density at the X-point, n e,Xpt , etc. The reduced model was quantitatively characterized using one half of the code cases; it was then used to successfully predict (replicate) the code values of e.g. n e,Xpt for the other half.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The physics basis to integrate an MHD stable, high-power hybrid scenario to a cool divertor for steady-state reactor operation

Abstract Coupling a high-performance core to a low heat flux divertor is a crucial step for ITER and a Fusion Pilot Plant or DEMO. Experiments in DIII-D recently expanded the steady-state hybrid scenario to high density and divertor impurity injection to study the feasibility of a radiating mantle solution. This work presents the physics basis for trade-offs between density, current drive and stability to tearing modes (TMs) at high β. EC power is crucial to tailor the plasma profiles into a passively stable state, and to eject impurities from the core. Off-axis EC depositions decrease the heating efficiency, but calculated electron heat transport coefficients show that this effect is partially mitigated by improved confinement inside the EC deposition. Additionally, the reduction in pressure is recovered by increasing the density. This favourable scaling of confinement with density was observed in high power plasmas for years, and this work provides a comprehensive explanation. ELITE predictions indicate that a path in peeling-ballooning stability opens up for certain conditions of density, power, q 95 and shaping, allowing the edge pressure to continue increasing without encountering a limit. In the core, calculated anomalous fast-ion diffusion coefficients are consistent with density fluctuation measurements in the toroidicity-induced Alfvén eigenmode range, showing that smaller fast-ion losses contribute to the enhanced confinement at high density. The edge integration study shows that divertor heat loads can be reduced with Ne and Ar injection, but this eventually triggers a cascade of n = 1, 2, 3 core TMs. We can now show that impurity radiation in the core is small and it is not the cause for the drop in confinement at high Ar and Ne injection rates. The overlap between the core TMs is consistent with the loss of pressure as estimated by the Belt model for the coupled rational surfaces. Optimization of these trade-offs has achieved plasmas with sustained H 98y2 = 1.7, f GW = 0.7 and ∼85% mantle radiation. The scenario and its variations at higher density and on- vs off-axis EC heating has been studied as a candidate for an integrated solution for several reactor designs, such as ITER, ARC, and the ARIES-ACT1 case, showing promising results in terms of fusion power and gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Transport barrier and spinning blob dynamics in the tokamak edge

Abstract In this work, we investigate the dynamics of plasma blobs in the edge of magnetic confinement devices using a full-f gyrokinetic particle-in-cell code with X-point geometry. In simulations, the evolution of a seeded blob is followed as it approaches a naturally-forming zonal shear layer near the separatrix, where the blob is stabilized by a large spin induced by the self-consistent adiabatic electron response, and blob bifurcation and trapping are observed during the cross-field propagation of blobs. A new theoretical explanation in both the zonal free and zonal shear layer is constructed, where the dominant E × B spin motion is included. A theoretical condition for a transport barrier induced by the interaction between spinning blobs and the zonal shear layer is obtained, and its scaling is verified with simulations. The new theoretical framework, especially the transport barrier, is applicable to explain and predict various experimental phenomena. In particular, the transport barrier condition calculated with experimental parameters demonstrates that the blob radial transport for H mode is smaller than L mode in experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Observation of long-radial-range-correlation in turbulence in high-collisionality high-confinement plasmas on DIII-D

Abstract We report on the observation of spatially asymmetric turbulent structures with a long radial correlation length in the core of high-collisionality H -mode plasmas on DIII-D tokamak. These turbulent structures develop from shorter wavelength turbulence and have a radially elongated structure. The envelope of turbulence spans a broad radial range in the mid-radius region, leading to streamer-like transport events. The underlying turbulence is featured by intermittency, long-term memory effect, and the characteristic spectrum of self-organized criticality. The amplitude and the radial scale increase substantially when the shearing rate of the mean flow is reduced below the turbulent scattering rate. The enhanced long-radial-range-correlated (LRRC) transport events are accompanied by apparent degradation of normalized energy confinement time. The emergence of such LRRC transport events may serve as a candidate explanation for the degrading nature of H-mode core plasma confinement at high collisionality on DIII-D tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Observation of a new pedestal stability regime in MAST Upgrade H-mode plasmas

Abstract The first pedestal stability and structure analysis on the new MAST Upgrade (MAST-U) spherical tokamak H-mode plasmas is presented. Our results indicate that MAST-U pedestals are close to the low toroidal mode number ( n ) peeling branch of the peeling-ballooning instability, in contrast with MAST H-mode pedestals which were deeply in the high- n ballooning branch. This offers the possibility of reaching the ELM-free quiescent H-mode (Burrell et al 2005 Plasma Phys. Control. Fusion 47 B37–B52) or high-performance super H-mode (Snyder et al 2015 Nucl. Fusion 55 083026; Snyder et al 2019 Nucl. Fusion 59 086017) regimes. In addition, the coupling between the peeling and ballooning branches is weak in MAST-U, suggesting that a path to very high pedestal pressure gradient at high density may exist with sufficient heating power. A possible explanation for the differences between MAST and MAST-U pedestal stability is given in terms of plasma shaping parameters, in particular squareness and elongation, as well as the pedestal top temperature and collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Nonlinear modeling of ELM mitigation with RMP on HL-2A

Abstract Nonlinear modeling of mitigation of the edge localized mode (ELM) with resonant magnetic perturbation (RMP) is performed for the HL-2A tokamak, utilizing the three-dimensional (3D) magnetohydrodynamic code JOREK. Based on the 3D equilibrium established after application of the n = 1 ( n is the toroidal mode number) RMP at 4.9 kAt coil current with odd parity, ELM mitigation is successfully simulated consistent with the experimental result. Nonlinear simulations show strong mode coupling among toroidal Fourier harmonics, allowing redistribution of the magnetic energy such that the most unstable toroidal mode saturates at a lower level. This magnetic energy cascade offers an explanation of the RMP-induced ELM mitigation achieved in HL-2A. Detailed examination of the simulation results shows persistent resonant field screening even during the ELM mitigation phase. Finite plasma resistivity however does enable partial penetration of the resonant field thus modifying the edge magnetic topology and characteristics of the edge transport. Plasma radial profiles undergo pronounced changes around the pedestal region, when the magnetic energy of the most unstable toroidal mode reaches the maximum value. Systematic scans of the applied RMP coil current with the JOREK simulations find a threshold value of around 4.5 kAt required for achieving the ELM mitigation on HL-2A.

Physics↗
DOE OSTI2025

Designing a validation experiment for radio frequency condensation

Abstract Theoretical studies have suggested that nonlinear effects can lead to ‘radio frequency (RF) condensation’, where an initially broad current profile can coalesce in islands when they reach sufficient width. In suitable conditions, RF condensation can ‘self-focus’ the driven current to the center of an island, improving stabilization efficiency and reducing control complexity. In unsuitable conditions, the effect can prematurely deplete the RF energy before it reaches the island center, impairing stabilization. It is predicted that the RF condensation effect can significantly impact reactor-scale tokamaks. This paper presents a set of simulations investigating the conditions under which RF condensation might be encountered in present-day tokamaks. For concreteness, the calculations use equilibrium reconstructions for two shots from DIII-D and AUG. The Current Condensation Amid Magnetic Islands (OCCAMI) simulation code has been used for this investigation. The code takes as its input a numerically specified axisymmetric EFIT equilibrium solution, and it perturbatively constructs a 3D field with an island embedded at the appropriate rational surface. In the OCCAMI code, the GENRAY code is used for ray tracing and for calculating the power deposition along a ray trajectory, and GENRAY is coupled self-consistently to a solution of the thermal diffusion equation in the island. The simulation results described in the paper illuminate the conditions required for experimental validation of the theory of RF condensation. The simulations also provide an explanation of why the effect was not noticed in experiments prior to the publication of theoretical papers on the subject.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗