Engineering topics
Berk, H. L.
Publications and source records attributed to Berk, H. L..
Resonances between high energy particles and ideal magnetohydrodynamic modes in tokamaks
Particle trajectory surfaces in an ideal magnetohydrodynamic high energy particle resonance are studied using kinetic Poincaré plots and through a calculation by perturbing near the resonance and finding canonical variables in the resonance, allowing the study of the distortion of the structure from that of a simple pendulum and to assist in the construction of models for the modification of particle distributions due to the modes. As a result, it is found that the narrow structure of an ideal mode eigenfunction can lead to a significant decrease in the resonance width compared to a case in which the eigenfunction does not vary within the resonant island.
Resonance broadened quasi-linear (RBQ) model for fast ion distribution relaxation due to Alfvénic eigenmodes
The burning plasma performance is limited by the confinement of the super-alfvénic fusion products such as alpha particles and the auxiliary heating ions capable of exciting the Alfvénic eigenmodes (AEs) (Gorelenkov et al 2014 Nucl. Fusion 54 125001). Here in this work the effect of AEs on fast ions is formulated within the quasi-linear (QL) theory generalized for this problem recently (Duarte 2017 PhD Thesis University of São Paulo, Brazil). The generalization involves the resonance line broadened interaction of energetic particles (EP) with AEs supplemented by the diffusion coefficients depending on EP position in the velocity space. A new resonance broadened QL code (or RBQ1D) based on this formulation allowing for EP diffusion in radial direction is built and presented in details. In RBQ1D applications we reduce the wave particle interaction (WPI) dynamics to 1D case when the particle kinetic energy is nearly constant. The diffusion equation for EP distribution evolution is then solved simultaneously for all particles along the angular momentum direction. We make initial applications of the RBQ1D to a DIII-D plasma with elevated q-profile where the beam ions show stiff transport properties (Collins et al (The DIII-D Team) 2016 Phys. Rev. Lett. 116 095001). In conclusion, AE driven fast ion profile relaxation is studied for validations of the QL approach in realistic conditions of beam ion driven instabilities in DIII-D.
Study of the likelihood of Alfvénic mode bifurcation in NSTX and predictions for ITER baseline scenarios
Rare Alfvénic wave transitions between fixed-frequency and chirping phases are identified in NSTX, where Alfvénic waves are normally observed to exhibit either chirping or avalanching responses. For those transitions, we apply in this paper a criterion (Duarte et al 2017 Nucl. Fusion 57 054001) to predict the nature of fast ion redistribution in tokamaks to be in the convective or diffusive nonlinear regimes. For NSTX discharges in which the transition is not accompanied by changes in the beam deposited power or modifications in the injected radiofrequency power, it has been found that the anomalous fast ion transport is a likely mediator of the bifurcation between the fixed-frequency mode behavior and rapid chirping. For a quantitative assessment, global gyrokinetic simulations of the effects of electrostatic ion temperature gradient turbulence and trapped electron mode turbulence on chirping were pursued using the GTS code. The investigation is extended by means of predictive studies of the probable spectral behavior of Alfvénic eigenmodes for baseline ITER cases consisting of elmy, advanced and hybrid scenarios. Finally, it has been observed that most modes are found to be borderline between the steady and the chirping phases.
Theory and observation of the onset of nonlinear structures due to eigenmode destabilization by fast ions in tokamaks
Alfvén waves can induce the ejection of fast ions in different forms in tokamaks. In order to develop predictive capabilities to anticipate the nature of fast ion transport, a methodology is proposed to differentiate the likelihood of energetic-particle-driven instabilities to produce frequency chirping or fixed-frequency oscillations. Furthermore, the proposed method employs numerically calculated eigenstructures and multiple resonance surfaces of a given mode in the presence of energetic ion drag and stochasticity (due to collisions and micro-turbulence). Toroidicity-induced, reversed-shear and beta-induced Alfvén-acoustic eigenmodes are used as examples. Waves measured in experiments are characterized, and compatibility is found between the proposed criterion predictions and the experimental observation or lack of observation of chirping behavior of Alfvénic modes in different tokamaks. It is found that the stochastic diffusion due to micro-turbulence can be the dominant energetic particle detuning mechanism near the resonances in many plasma experiments, and its strength is the key as to whether chirping solutions are likely to arise. We proposed a criterion that constitutes a useful predictive tool in assessing whether the nature of the transport for fast ion losses in fusion devices will be dominated by convective or diffusive processes.
Explosive attractor solutions to a universal cubic delay equation
New explosive attractor solutions have been found in a universal cubic delay equation that has been studied in both the plasma and the fluid mechanics literature. Through computational simulations and analytic approximations, it is found that the oscillatory component of the explosive mode amplitude solutions are described through multi-frequency Fourier expansions with respect to a pseudo-time variable. The spectral dependence of these solutions as a function of a system parameter, ϕ, is studied. The mode amplitude that is described in the explosive regime has two main features: a well-known envelope (t0−t)−5/2, with t0 the blow-up time of the amplitude, and a spectrum of discrete oscillations with ever-increasing frequencies, which may give experimental information about the properties of a system's equilibrium.
Calculation of the fast ion tail distribution for a spherically symmetric hot spot
The fast ion tail for a spherically symmetric hot spot is computed via the solution of a simplified Fokker-Planck collision operator. Emphasis is placed on describing the energy scaling of the fast ion distribution function in the hot spot as well as the surrounding cold plasma throughout a broad range of collisionalities and temperatures. It is found that while the fast ion tail inside the hot spot is significantly depleted, leading to a reduction of the fusion yield in this region, a surplus of fast ions is observed in the neighboring cold plasma region. The presence of this surplus of fast ions in the neighboring cold region is shown to result in a partial recovery of the fusion yield lost in the hot spot.
Erratum: “1.5D quasilinear model and its application on beams interacting with Alfvén eigenmodes in DIII-D” [Phys. Plasmas 19, 092511 (2012)]
No abstract prepared.
Comparing the line broadened quasilinear model to Vlasov code
The Line Broadened Quasilinear (LBQ) model is revisited to study its predicted saturation level as compared with predictions of a Vlasov solver BOT [Lilley et al., Phys. Rev. Lett. 102, 195003 (2009) and M. Lilley, BOT Manual. The parametric dependencies of the model are modified to achieve more accuracy compared to the results of the Vlasov solver both in regards to a mode amplitude's time evolution to a saturated state and its final steady state amplitude in the parameter space of the model's applicability. However, the regions of stability as predicted by LBQ model and BOT are found to significantly differ from each other. The solutions of the BOT simulations are found to have a larger region of instability than the LBQ simulations.
Nonlinearly driven harmonics of Alfvén modes
In order to study the leading order nonlinear magneto-hydrodynamic (MHD) harmonic response of a plasma in realistic geometry, the AEGIS code has been generalized to account for inhomogeneous source terms. These source terms are expressed in terms of the quadratic corrections that depend on the functional form of a linear MHD eigenmode, such as the Toroidal Alfvén Eigenmode. The solution of the resultant equation gives the second order harmonic response. Preliminary results are presented here.
Magnetically driven jets and winds
Four equations for the origin and propagation of nonrelativistic jets and winds are derived from the basic conservation laws of ideal MHD. The axial current density is negative in the vicinity of the axis and positive at larger radii; there is no net current because this is energetically favored. The magnetic field is essential for the jet solutions in that the zz-component of the magnetic stress acts, in opposition to gravity, to drive matter through the slow magnetosonic critical point. For a representative self-consistent disk/jet solution relevant to a protostellar system, the reaction of the accreted mass expelled in the jets is 0.1, the ratio of the power carried by the jets to the disk luminosity is 0.66, and the ratio of the boundary layer to disk luminosities is less than about 0.13. The star's rotation rate decreases with time even for rotation rates much less than the breakup rate.
Linearized theory of inhomogeneous multiple 'water-bag' plasmas
Equations are derived for describing the inhomogeneous equilibrium and small deviations from the equilibrium, giving particular attention to systems with trapped particles. An investigation is conducted of periodic systems with a single trapped-particle water bag, taking into account the behavior of the perturbation equations at the turning points. An outline is provided concerning a procedure for obtaining the eigenvalues. The results of stability calculations connected with the sideband effects are considered along with questions regarding the general applicability of the multiple water-bag approach in stability calculations.
Stability of a trapped particle equilibrium
Stability of trapped particle equilibrium, considering contribution to plasma wave dispersion
Nonlinear study of Vlasov's equation for a special class of distribution functions
Nonlinear study of Vlasov equation for distribution function for one-dimensional rarefied plasma