Engineering topics

Catto, Peter J.

Publications and source records attributed to Catto, Peter J..

12 records

DOE OSTI2023

Merging of the superbanana plateau and $\sqrt{v}$ transport regimes in nearly quasisymmetric stellarators

We report alpha particle confinement is one of the most demanding issues for stellarators. It now seems clear that it is possible to design optimized stellarators that confine the background plasma at near tokamak radial transport levels. Moreover, adequate collisionless alpha particle confinement is possible in the core of a highly optimized stellarator. Here, the collisional confinement of barely trapped alphas in an optimized stellarator is considered by accounting for the resonance due to the reversal in direction of the drift within a flux surface and investigating the sensitive role of magnetic shear in keeping this resonance close to the passing boundary in some nearly quasisymmetric stellarator configurations. The treatment relies on a narrow collisional boundary layer formulation that combines the responses of both these resonant pitch angle alphas and the remaining barely trapped alphas. A novel merged regime treatment leads to explicit expressions for the energy diffusivity for both superbanana plateau (or resonant plateau) and $\sqrt{v}$ transport in the large aspect ratio limit for a slowing down tail alpha distribution function, where v is the effective pitch angle scattering collision frequency of the trapped alphas off the background ions. Depending on the details of the optimization scheme and the sign of the magnetic shear, modest magnetic shear can be used to reduce superbanana (or resonant) plateau transport to below the $\sqrt{v}$ transport level. In addition, a quasilinear equation retaining spatial diffusion is derived for a general alpha distribution function that allows the radial alpha transport to modify the distribution so it is no longer isotropic in velocity space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Lower hybrid current drive in a tokamak for correlated passes through resonance

Standard quasilinear descriptions are based on the constant magnetic field form of the quasilinear operator so improperly treat the trapped electron modifications associated with tokamak geometry. Moreover, successive poloidal transits of the Landau resonance during lower hybrid current drive in a tokamak are well correlated, and these geometrical details must be properly retained to account for the presence of trapped electrons that do not contribute to the driven current. The recently derived quasilinear operator in tokamak geometry accounts for these features and finds that the quasilinear diffusivity is proportional to a delta function with a transit or bounce averaged argument (rather than a local Landau resonance condition). The new quasilinear operator is combined with the Cordey (Nucl. Fusion, vol. 16, 1976, pp. 499–507) eigenfunctions to properly derive a rather simple and compact analytic expression for the trapped electron modifications to the driven lower hybrid current and the efficiency of the current drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Collisional broadening of nonlinear resonant wave–particle interactions

A general procedure for understanding plasma behaviour when resonant wave–particle interactions are the sole destabilizing and transport mechanism or only heating and/or current drive source is highlighted without recourse to involved numerical or analytical treatments. These phenomena are characterized by transport that appears to be collisionless even though collisions play a central role in narrow collisional boundary layers. The order of magnitude estimates, which include nonlinear effects, are shown to provide expressions in agreement with the principal results of recent toroidal Alfvén eigenmode (TAE), toroidal magnetic field ripple, and heating and current drive treatments. More importantly, the retention of nonlinearities leads to new estimates of the alpha particle energy diffusivity at saturation for TAE modes, and the ripple threshold at which superbanana plateau evaluations of alpha particle transport are modified by nonlinear radial drift effects. In addition, the estimates indicate when quasilinear descriptions for heating and current drive will begin to fail. The phenomenological procedure demonstrates that in magnetic fusion relevant plasmas, narrow collisional boundary layers must be retained for resonant wave–particle interactions as they enhance the role of collisions, and make stochastic particle motion unlikely to be more important than other nonlinear processes.

Physics↗

DOE OSTI2021

Collisional broadening of nonlinear resonant wave-particle interactions

A general procedure for understanding plasma behavior when resonant wave-particle interactions are the sole destabilizing and transport mechanism or only heating and/or current drive source is highlighted without recourse to involved numerical or analytical treatments. These phenomena are characterized by transport that appears to be collisionless even though collisions play a central role in narrow collisional boundary layers. The order of magnitude estimates, which include nonlinear effects, are shown to provide expressions in agreement with the principle results of recent toroidal Alfvén eigenmode (TAE), toroidal magnetic field ripple, and heating and current drive treatments. More importantly, the retention of nonlinearities leads to new estimates of the alpha particle energy diffusivity at saturation for TAE modes, and the ripple threshold at which superbanana plateau evaluations of alpha particle transport are modified by nonlinear radial drift effects. In addition, the estimates indicate when quasilinear descriptions for heating and current drive will begin to fail. The phenomenological procedure demonstrates that in magnetic fusion relevant plasmas narrow collisional boundary layers must be retained for resonant wave-particle interactions as they enhance the role of collisions, and make stochastic particle motion unlikely to be more important than other nonlinear process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Reimagining full wave rf quasilinear theory in a tokamak

The velocity dependent resonant interaction of particles with applied radio frequency (rf) waves during heating and current drive in the presence of pitch angle scattering collisions gives rise to narrow collisional velocity space boundary layers that dramatically enhance the role of collisions as recently shown by Catto (J. Plasma Phys., vol. 86, 815860302, 2020). The behavior is a generalization of the narrow collisional boundary layer that forms during Landau damping as found by Johnston (Phys. Fluids, vol. 14, 1971, pp. 2719-2726) and Auerbach (Phys. Fluids, vol. 20, 1977, pp. 1836-1844). For a wave of parallel wave number k|| interacting with weakly collisional plasma species of collision frequency ν and thermal speed vth , the effective collision frequency becomes of order ν(k_||v_th /ν)^2/3>> ν . The narrow boundary layers that arise because of the diffusive nature of the collisions allows a physically meaningful wave-particle interaction time to be defined that is the inverse of this effective collision frequency. The collisionality implied by the narrow boundary layer results in changes in the standard quasilinear treatment of applied rf fields in tokamaks while remaining consistent with causality. These changes occur because successive poloidal interactions with the rf are correlated in tokamak geometry and because the resonant velocity space dependent interactions are controlled by the spatial and temporal behavior of the perturbed full wave fields rather than just the spatially local Landau and Doppler shifted cyclotron wave-particle resonance condition associated with unperturbed motion of the particles. The correlation of successive poloidal circuits of the tokamak leads to the appearance in the quasilinear operator of transit averaged resonance conditions localized in velocity space boundary layers that maintain negative definite entropy production.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Lower hybrid current drive in a tokamak for correlated passes through resonance

Standard quasilinear descriptions are based on the constant magnetic field form of the quasilinear operator so improperly treat the trapped electron modifications associated with tokamak geometry. Moreover, successive poloidal transits of the Landau resonance during lower hybrid current drive in a tokamak are well correlated, and these geometrical details must be properly retained to account for the presence of trapped electrons that do not contribute to the driven current. The recently derived quasilinear operator in tokamak geometry accounts for these features and finds that the quasilinear diffusivity is proportional to a delta function with a transit or bounce averaged argument (rather than a local Landau resonance condition). The new quasilinear operator is combined with the Cordey ( Nucl. Fusion , vol. 16, 1976, pp. 499–507) eigenfunctions to properly derive a rather simple and compact analytic expression for the trapped electron modifications to the driven lower hybrid current and the efficiency of the current drive.

Physics↗

DOE OSTI2021

Collisional effects on resonant particles in quasilinear theory

A careful examination of the effects of collisions on resonant wave-particle interactions leads to an alternate interpretation and deeper understanding of the quasilinear operator originally formulated by Kennel and Engelmann (Phys. Fluids vol. 9, 1966, pp. 2377- 2388) for collisionless, magnetized plasmas, and widely used to model radio frequency heating and current drive. The resonant and nearly resonant particles are particularly sensitive to collisions that pitch angle scatter them out of and into resonance. As a result, the resonant particle-wave interactions occur in the center of a narrow collisional boundary when the collision frequency nu is very small compared to the wave frequency omega. The diffusive nature of the pitch angle scattering combined with the wave-particle resonance condition enhances the collision frequency by (omega/nu)2/3 >>1, resulting in an effective resonant particle collision time of tau_int ~ (nu /omega)2/3 nu <<1/ nu . A rigorous collisional boundary layer analysis generalizes the standard quasilinear operator to a form that is fully consistent with Kennel-Englemann, but allows replacing the delta function appearing in the diffusivity with a simple integral (having the appropriate delta function limit) retaining the new physics associated with the narrow boundary layer, while preserving the entropy production principle. The limitations of the collisional boundary layer treatment are also estimated, and indicate that substantial departures from Maxwellian are not permitted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Reimagining full wave rf quasilinear theory in a tokamak

The velocity dependent resonant interaction of particles with applied radiofrequency (rf) waves during heating and current drive in the presence of pitch angle scattering collisions gives rise to narrow collisional velocity space boundary layers that dramatically enhance the role of collisions as recently shown by Catto ( J. Plasma Phys. , vol. 86, 2020, 815860302). The behaviour is a generalization of the narrow collisional boundary layer that forms during Landau damping as found by Johnston ( Phys. Fluids , vol. 14, 1971, pp. 2719–2726) and Auerbach ( Phys. Fluids , vol. 20, 1977, pp. 1836–1844). For a wave of parallel wave number ${k_{||}}$ interacting with weakly collisional plasma species of collision frequency $\nu$ and thermal speed ${v_{\textrm{th}}}$ , the effective collision frequency becomes of order $\nu {({k_{||}}{v_{th}}/\nu )^{2/3}} \gg \nu $ . The narrow boundary layers that arise because of the diffusive nature of the collisions allow a physically meaningful wave–particle interaction time to be defined that is the inverse of this effective collision frequency. The collisionality implied by the narrow boundary layer results in changes in the standard quasilinear treatment of applied rf fields in tokamaks while remaining consistent with causality. These changes occur because successive poloidal interactions with the rf are correlated in tokamak geometry and because the resonant velocity space dependent interactions are controlled by the spatial and temporal behaviour of the perturbed full wave fields rather than just the spatially local Landau and Doppler shifted cyclotron wave–particle resonance condition associated with unperturbed motion of the particles. The correlation of successive poloidal circuits of the tokamak leads to the appearance in the quasilinear operator of transit averaged resonance conditions localized in velocity space boundary layers that maintain negative definite entropy production.

Physics↗

DOE OSTI2021

Drift kinetic theory of alpha transport by tokamak perturbations

We report that upcoming tokamak experiments fuelled with deuterium and tritium are expected to have large alpha particle populations. Such experiments motivate new attention to the theory of alpha particle confinement and transport. A key topic is the interaction of alpha particles with perturbations to the tokamak fields, including those from ripple and magnetohydrodynamic modes like Alfvén eigenmodes. These perturbations can transport alphas, leading to changed localization of alpha heating, loss of alpha power and damage to device walls. Alpha interaction with these perturbations is often studied with single-particle theory. In contrast, we derive a drift kinetic theory to calculate the alpha heat flux resulting from arbitrary perturbation frequency and periodicity (provided these can be studied drift kinetically). Novel features of the theory include the retention of a large effective collision frequency resulting from the resonant alpha collisional boundary layer, correlated interactions over many poloidal transits and finite orbit effects. Heat fluxes are considered for the example cases of ripple and the toroidal Alfvén eigenmode (TAE). The ripple heat flux is small. The TAE heat flux is significant and scales with the square of the perturbation amplitude, allowing the derivation of constraints on mode amplitude for avoidance of significant alpha depletion. A simple saturation condition suggests that TAEs in one upcoming experiment will not cause significant alpha transport via the mechanisms in this theory. However, saturation above the level suggested by the simple condition, but within numerical and experimental experience, which could be accompanied by the onset of stochasticity, could cause significant transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Collisional effects on resonant particles in quasilinear theory

A careful examination of the effects of collisions on resonant wave–particle interactions leads to an alternate interpretation and deeper understanding of the quasilinear operator originally formulated by Kennel & Engelmann ( Phys. Fluids , vol. 9, 1966, pp. 2377–2388) for collisionless, magnetized plasmas, and widely used to model radio frequency heating and current drive. The resonant and nearly resonant particles are particularly sensitive to collisions that scatter them out of and into resonance, as for Landau damping as shown by Johnston ( Phys. Fluids , vol. 14, 1971, pp. 2719–2726) and Auerbach ( Phys. Fluids , vol. 20, 1977, pp. 1836–1844). As a result, the resonant particle–wave interactions occur in the centre of a narrow collisional boundary when the collision frequency $\unicode[STIX]{x1D708}$ is very small compared to the wave frequency $\unicode[STIX]{x1D714}$ . The diffusive nature of the pitch angle scattering combined with the wave–particle resonance condition enhances the collision frequency by $(\unicode[STIX]{x1D714}/\unicode[STIX]{x1D708})^{2/3}\gg 1$ , resulting in an effective resonant particle collisional interaction time of $\unicode[STIX]{x1D70F}_{\text{int}}\sim (\unicode[STIX]{x1D708}/\unicode[STIX]{x1D714})^{2/3}/\unicode[STIX]{x1D708}\ll 1/\unicode[STIX]{x1D708}$ . A collisional boundary layer analysis generalizes the standard quasilinear operator to a form that is fully consistent with Kennel–Englemann, but allows replacing the delta function appearing in the diffusivity with a simple integral (having the appropriate delta function limit) retaining the new physics associated with the narrow boundary layer, while preserving the entropy production principle. The limitations of the collisional boundary layer treatment are also estimated, and indicate that substantial departures from Maxwellian are not permitted.

Physics↗

DOE OSTI2019

Bootstrap current and parallel ion velocity in imperfectly optimized stellarators

A novel derivation of the parallel ion velocity, and the bootstrap and Pfirsch-Schlüter currents in an imperfectly optimized (that is, almost omnigenous) stellarator magnetic field, \vec B, is presented. It is shown that, when the conventional radially local form of the drift kinetic equation is employed, the flow velocity and the bootstrap current acquire a spurious contribution proportional to ω /ν , where ω denotes the \vec E × \vec B rotation frequency (due to the radial electric field E ) and ν the collision frequency. This contribution is particularly large in the squareroot ν regime and at smaller collisionalities, where ω /ν > 1 , and is presumably present in most numerical calculations, but it disappears if a more accurate drift kinetic equation is used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Collisional alpha transport in a weakly non-quasisymmetric stellarator magnetic field

Alpha particle confinement is a serious concern in stellarators and provides strong motivation for optimizing magnetic field configurations. In addition to the collisionless confinement of trapped alphas in stellarators, excessive collisional transport of the trapped alpha particles must be avoided while they tangentially drift due to the magnetic gradient (the ∇B drift). The combination of pitch angle scatter off the background ions and the ∇B drift gives rise to two narrow boundary layers in the trapped region. The first is at the trapped-passing boundary and enables the finite trapped response to be matched to the vanishing passing response of the alphas. The second layer is a region that encompasses the somewhat more deeply trapped alphas with vanishing tangential ∇B drift. Away from (and between) these boundary layers, collisions are ineffective and the alpha ∇B drift simply balances the small radial drift of the trapped alphas. As this balance does not vanish as the trapped-passing boundary is approached, the first collisional boundary layer is necessary and gives rise to √ν transport, with ν the collision frequency. The vanishing of the tangential drift results in a separate, somewhat wider boundary layer, and significantly stronger superbanana plateau transport that is independent of collisionality. The constraint imposed by the need to avoid significant energy depletion loss in the slowing down tail distribution function sets the allowed departure of a stellarator from an optimal quasisymmetric configuration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗