Hinton 1976
Hinton 1976
Hinton 1976
The dissipation induced by coulomb-collisional scattering provides an irreducible minimum, and thus a
useful standard for comparison, for transport processes in a hot, magnetically confined plasma. The kinetic
description of this dissipation is provided by an equation of the Fokker — Planck form. As in the standard
transport theory for a neutral gas, approximate solution of the Fokker — Planck equation permits the
calculation of transport coefficients, which linearly relate the fluxes of particles, energy, and electric charge,
to the density and temperature gradients, and to the electric field. The transport relations are useful in
studying the confinement properties of present and future experimental devices for research in controlled
thermonuclear fusion. The transport theory for a magnetized plasma (in which the Larmor radius is much
smaller than gradient scale lengths describing the plasma fluid) departs from the theory for a neutral gas in
several fundamental ways. Thus, transport coefficients for a magnetized plasma can be calculated even
when the collisional mean free path is much longer than the gradient scale length (as would pertain in
thermonuclear regimes). Such transport coefficients are generally nonlocal, being defined in terms of
averages over surfaces with macroscopic dimensions. Furthermore, when the mean free path is long, the
magnetized-plasma transport coefficients depend crucially upon the magnetic field geometry, the effects of
which must be treated at the kinetic level of the Fokker — Planck equation. The results display several novel
couplings between collisional dissipation and the electromagnetic field. The present review of magnetized-
plasma transport theory is intended to be as widely accessible as possible. Thus the relevant features of
magnetic confinement in closed (toroidal) systems, and of charged particles in spatially varying fields, are
derived, at least in outline, from first principles. Although consideration is given to "classical" transport in
which most field geometric effects are omitted, major emphasis is placed on the "neoclassical" theory
which has been developed over the last decade, Neoclassical transport coefficients are specifically relevant
to a magnetically confined plasma, rather than to just a magnetized plasma; their unusual features, such as
nonlocality and geometry dependence, become particularly important in the high temperature regime of
proposed thermonuclear reactors, The area of neoclassical theory which seems most complete —
its
application to axisymmetric tokamak-type confinement systems —
is correspondingly stressed.
Reviews of Modern Physics, Vol. 48, No. 2, Part I, April 1976 Copyright 1976 American Physical Society 239
240 F. L. Hinton and R. D. Hazeltine: Theory of plasma transport
priately defined surface, even as, in response to trans- roidal plasma confinement has been presented by Artsi-
port, the surface configuration changes in time. movich (19'l2).
It is useful to distinguish two phases in the develop-
ment of magnetized-plasma transport theory. In the A. History
earlier, classical phase, the appropriate form of the In less than ten years, the theory of neoclassical
Fokker-Planck collision operator was established, and transport has acquired a vast literature. We survey the
then used to analyze transport properties both parallel major contributions here, assuming (conveniently, but
and perpendicular to the magnetic field. It is a char- only temporarily) that the rea. der is familiar with neo-
acteristic of classical transport theory that the spatial classical terminology.
variation of the magnetic field has no effect, at least at The kinetic theory of plasma transport in tokamaks
the level of the kinetic equation; classical transport co- was initiated by Galeev and Sagdeev(1968). They showed
efficients can be, and usually were, calculated for the that trapped particles, with "banana" orbits, are re-
case of a uniform field. The second phase of develop- sponsible for a significant enhancement of the calcu-
ment, which is now generally referred to as rseoelzz- lated diffusion coefficient and thermal conductivity, at
sieal, began in the 1960's, especially after the pioneer- very small values of collision frequency (i. e. , in the
ing work of Galeev and Sagdeev (1968). Neoclassical banana regime). ' The enhancement was explained qual-
transport theory is also based on the use of the Fokker- itatively on the basis of a random walk of the banana
Planck collision operator; it departs from classical orbits, mith a step size significantly larger than the
theory in recognizing that the magnetic field of a toroi- Larmor radius. Qaleev and Sagdeev also demonstrated
dal confinement system is necessarily nonuniform, and the existence of an intermediate collision frequency re-
that in the long mean-free-path regime appropriate to gime, in mhich the diffusion coefficient is independent
thermonuclear temperatures, the spatial variation of of collision frequency (the "plateau regime").
the field has crucial effects, which must be treated at Plasma diffusion in a stellerator was investigated by
the kinetic level. Galeev, Sagdeev, Furth, and Rosenbluth (1969). They
The crucial effects of field variation are particle found a maximum in the diffusion coefficient at low col-
trapping and the gradient-& and curvature drifts of gy- lision frequencies, due to the existence of "localized"
rating particles across field lines, resulting from grad- particles (Gibson and Taylor, 196'l; Gibson and Mason,
ual distortion of the Larmor orbit. Thus the magnetic 1968), with "superbanana" orbits.
inhibition of free particle motion, mentioned previously. , The neoclassical theory of plasma transport generated
has complementary weakness. Because it provides lo- considerable interest, because of its novel character
calization of charged particles in only two dimensions, and because of its possible importance in determining
it requires a toroidal system for containment; but the toroidal confinement properties. Kovrizhnykh (1969)
field curvature resulting from toroidicity leads to per- and Rutherford (1970) independently pointed out that the
pendicular drifts, i. e. , to a weakening of the original ion and electron diffusion rates should be equal and in-
inhibition. Neoclassical transport is a prominent result dependent of the radial electric field, in axisymmetric
of these weaknesses. It should be noted that, while clas- systems, such as tokamaks. Kovrizhnykh calculated
sical and neoclassical transport processes are additive, the transport coefficients for a variety of confinement
neoclassical effects are typically much larger; recall devices (tokamak, levitron, stellerator, and bumpy
the previous remark (ii). torus), and considered both weakly ionized and fully
The present review considers classical theory only ionized plasmas. The details of the calculations, for
briefly, for purposes of comparison and for complete- axisymmetric systems, are particularly clear in Ruth-
ness. A relatively comprehensive treatment of neoclas- erford's paper. Frieman (1970) considered diffusion in
sical theory is attempted. We have tried to make the general nonaxisymmetric systems. His paper contains
argument self-contained, and accessible to as wide an a derivation of the drift kinetic equation, which mas the
audience as possible, including researchers unfamiliar starting point for the Qaleev-Sagdeev calculation.
with plasma transport theory. Especially for the bene- Other neoclassical transport properties mere also
fit of the latter, we draw attention to the following pre- soon discovered. The effect of trapped electrons on the
vious reviews. Ohmic conductivity of a tokamak plasma, was consid-
The famous monograph of Spitzer (1967) discusses, ered by Hinton and Oberman (1969) who showed that, at
among other relevant topics, the nature of Coulomb low collision frequencies, the conductivity is reduced by
collisions in a plasma, and a number of classical trans- an amount proportional to the fraction of trapped elec-
port processes. More detailed treatments of the clas-' trons. The possibility of a new kind of pinch effect, due
sical theory, with extensive references, may be found to trapped particles, was discovered independently by
in Braginskii (1965) and Kaufman (1966). The micro- Ware (1970) and Galeev (1971). They pointed out that
scopic processes underlying neoclassical transport, as the toroidal electric field, which is necessarily present
related to field geometry, are considered by Kadomtsev in tokamaks, causes a radially inward motion of trapped
and Pogutse (1971), who also discuss a related class of particles, at a velocity which is much larger than the
plasma instabilities. Neoclassical transport is also
considered in a review by Galeev and Sagdeev (1975), The orders of magnitude of these coefficients in the high
and in a relatively elementary reviem by Hazeltine collision frequency regime had been worked out several years
—
(1975). Finally, a coherent if unfortunately somewhat earlier by Pfirsch and ScMutter (f962) and Shafranov (f965)
—
dated account of the e&Pemmental observations on to- from a fluid model.
Ex B drift in the classical pinch effect. The existence gimes by Tsang and Frieman (1975).
of a "diffusion driven" current ("bootstrap current" ) The banana regime analysis of Rosenbluth, Hazeltine
was predicted independently by Galeev (1971) [see also and Hinton (1972) was generalized to include noncircular
Sagdeev and Galeev (1970), and Galeev and Sagdeev cross section tokamaks, such as doublet (Ohkawa. , 1968)
(1971)] and Bickerton, Connor and Taylor (1971). The by Glasser and Thompson (1973). A further generaliza-
latter authors considered the possibility of a steady- tion, to include all collisionality regimes, was consid-
state "bootstrap" tokamak, in which the poloidal con- ered by Bernstein (1974). He generalized the variational
fining magnetic field would be provided by this "boot- principles used by Hinton and Rosenbluth (1973) to in-
strap" current, and in which the toroidal electric field clude arbitrary axisymmetric systems, and introduced
would be zero; this possibility was also considered in- a convenient inner-product formalism, which has been
dependently by Kadomtsev and Shafranov (1972). The adopted in this review.
physical origin of the bootstrap current is the banana In trying to present a coherent account of the subject
orbits of the trapped electrons, in the presence of a under review, we frequently use arguments quite differ-
density gradient. These result in a current along the ent from those of the original contributors. Yet we have
magnetic field lines, analogous to the classical diamag- attempted proper assignment of credit, not only to those
netic current due to particle gyration about the field authors mentioned in the brief history above, but to
lines. Through collisions with ions, these particles many others whose contributions have been significant.
become untrapped; through collisions with untrapped To those authors whose work is not included, we offer
electrons, the latter are "entrained. "
The result of our apologies; some omissions were necessary, in
these collisions is a current carried by untrapped elec- order that we might explore the chosen topics in some
trons. depth.
All of the calculations referred to used either simplified
collision operators to describe like-particle collisions, or B. Synopsis
used nonrigorous approximations to deal with the exact This review is concerned with the calculation of trans-
Fokker-Planck operators. The resulting numerical port coeff icients, f rom the Fokker- Planck equation, in
transport coefficients differed by as much as an order the small gyroradius limit. The transport coefficients
of magnitude from one paper to the next. Rosenbluth, linearly relate the fluxes (of particles, energy, and
Hazeltine and Hinton (1972) carried out an accurate and electric charge) to the thermodynamic forces (the pres-
systematic evaluation of the tokamak neoclassical coef- sure and temperature gradients for each charged par-
ficients, starting with the full Fokker- Planck collision ticle species in the plasma, and the electric field).
operator. By using a variational principle to treat the When combined with the exact conservation laws for
low collision frequency banana regime, the e'ffects of particles and energy, and with Maxwell's equations, the
like-particle collisions were included rigorously. Their linear transport relations provide a closed set of equa-
results are correct to lowest order in the inverse aspect tions, which predict the temporal evolution of the plasma
ratio, in the large aspect ratio, circular cross section fluid, from a properly chosen set of initial and boundary
case. conditions.
Although the correct transport coefficients were known The Fokker-Planck collision operator is introduced at
in the banana regime (Rosenbluth, Hazeltine and Hinton, the end of this section. In Sec. II, we postpone further
1972) and the plateau regime (Galeev, 1971), it was the kinetic analysis to consider the velocity moments of the
transition between these regimes which was thought to Fokker-Planck equation. The moment equations permit
be more relevant to the tokamak experiments of the late a relatively straightforward introduction of the small
1960's and early 1970's. Hinton and Rosenbluth (1973) gyroradius parameter (Sec. II.B). Then the toroidal na-
calculated the transport coefficients for a large aspect- ture of a magnetically confined plasma equilibrium—
ratio tokamak in the banana-plateau transitional regime. which is critical to neoclassical theory— can be estab-
They used numerical methods to make use of general- lished (Sec. II.C). We find that this equilibrium is nec-
izations of a variational principle due to Rutherford essarily dynamic, i. e. , it can exist only in the presence
(1970). The results showed that diffusion and thermal of significant, divergence-free, flows of particles and
conductivity coefficients increase monotonically with energy. A crucial averaging operator — the so-called
collision frequency, even in the plateau regime. —
flux surface average which is local with respect to mi-
None of the previous papers answered the question of nor toroidal radius, is introduced. We also derive con-
how the radial electric field is determined in the quasi- venient, perspicuous expressions for the cross-field
steady state of ambipolar diffusion in axisymmetric sys- particle and energy fluxes, and compare the classical
terns. Rosenbluth, Rutherford, Taylor, Frieman, and and neoclassical contributions (Sec. II.D). The changing
Kovrizhnykh (1971) studied this question for the banana character of diffusion, for different regimes of collision
regime. By going to fourth order in the gyroradius, frequency, is considered. These results are derived for
they obtained an equation which determines the time the case of a completely general confinement geometry;
rate of change of the toroidal plasma angular momen- simplifications permitted in the axisymmetric case, in
tum, in terms of radial diffusion of angular momentum. which the major toroidal axis is a symmetry axis, are
The radial electric field is thus determined, through its considered separately (Sec. II.E). Finally, the relevant
relation to the plasma toroidal angular velocity. This dissipation mechanisms are summarized, by considera-
question was also studied for the collisional regime by tion of the heat production rate (Sec. II.F).
Hazeltine (1974), and for the plateau and transitional re- As remarked previously, the slow drift motion of gy-
rating particles —the so-called guiding center drift— species. The linearized kinetic equations so obtained
across the nonuniform confining field, plays a funda- are novel also in their more general and rigorous treat-
mental role in neoclassical transport. Hence in Sec. III ment of certain ion-electron coupling effects. Next, we
we adopt a microscopic viewpoint to review the relevant show that the kinetic equations can be obtained from a
conclusions of the guiding center theory. As in Sec. II, variational principle, which is valid independently of the
an initially general formalism (Secs. III.A, III. B) is later size of the collision frequency, i. e. , in both short and
specialized to the imports. nt case of axisymmetry (Sec. long mean-free-path regimes (Sec. V. C). As a result,
III. C). Once the character of the collisionless guiding variational expressions for the transport coefficients
center orbits is established, conventional random walk are obtained, in terms of integrals of the distribution
arguments are used to estimate diffusion coefficients function. A simpler variational principle is derived for
under various circumstances (Sec. III. D). (These argu- the long mean-free-path limit (Sec. V. D). In.the oppo-
ments provide useful understanding and their conclu- site, collision-dominated limit, the kinetic equations
sions are approximately correct; but it should be noted can be solved directly (Sec. V. E); one thus obtains, by
that the details of the simplified random walk picture a somewhat different argument, the axisymmetric ver-
are not always borne out by kinetic theory. ) Next we sion of the results of Sec. IV. D.
consider the predictions of guiding center- theory with In long mean-free-path regimes, the linearized kinet-
regard to the single-particle distribution function, ic equations involve a quite complicated average of the
f(x, v, f), describing the expected density of particles at already formidable Fokker-Planck collision operator,
the point (x, v) in phase space (Sec. III. E). This is shown so that rigorous solutions have been obtained only in
to consist of two parts, f
=f +f, with the following prop- certain limits. For example, if the ionic charge Z is
erties: fis determined by Larmor gyration, and yields very large, the electron collision operator reduces to a
f
classical transport, while is determined by guiding Lorentz gas operator (to lowest order in Z '), and the
center motion, as modified by collisions, and yields kinetic theory is analytically tractable; this case is con-
neoclassical transport. A small gyroradius approxima- sidered in Sec. V. D. A more generally important limit
tion to the Fokker-Planck equation —
the so-called drift is that of a "thin" (large aspect-ratio) torus, whose mi-
kinetic equation —
which f
must satisfy, is derived here. nor radius is much smaller than its major radius. This
Having established, in Secs. II and III, precisely which case is considered in Sec. VI; we show that, for suffi-
moments of the distribution function are required, we ciently small collision frequency, the distribution func-
begin in Sec. IV the actual calculation of transport coef- tion is localized in velocity space, in such a way that
ficients. The classical perpendicular transport problem the collision operator may be expanded in powers of the
is considered first, and solved by straightforwardly square root of the inverse aspect ratio. Since the low-
combining some of our previous results (Sec. IV. B). We est-order version resembles a Lorentz gas operator,
then turn to the problem of collisional transport parallel the kinetic equations are readily solved. The velocity
to the magnetic field (Sec. IV. C). The variational meth-
od used to treat this problem has some resemblance to
space region in which f is localized changes as the col-
lision frequency increases (while remaining smaller than
techniques used in later sections, to which it therefore it is in the collision-dominated regime, where f is not
provides a useful introduction. Furthermore, we show localized). Hence somewhat different analytical treat-
that the parallel transport equations are readily com- ments are required, for the nearly collisionless limit
bined with certain conclusions of Sec. II, to yield the (Secs. VI. B, VI. C), and for the regime of intermediate
neoclassical (cross-field) transport coefficients for an collisionality (Sec. VI. D). Results for the intermediate
arbitrary toroidal geometry, in the large collision fre- f
regime are affected by the behaviour of near a certain
quency limit (Sec. IV. D). boundary layer in velocity space. Explicit neoclassical
Unfortunately, under the experimental conditions typ- transport coefficients, which are rigorous in the con-
ical of most present and planned toroidal confinement text of the aspect-ratio expansion, are obtained for the
sy stems, these short mean-f ree-path results pertain to low-to-intermediate transitional regime by substituting
only a small fraction of the plasma volume. To treat a numerical solutions to the large aspect-ratio kinetic
less collision-dominated plasma (including, in particu- equations into the general variational formalism of Sec.
lar, a thermonuclear plasma), more sophisticated kine- V.
tic theory, based on the drift kinetic equation, is re- For the regime of intermediate to high collisionality
quired. This theory, specialized to the case of an axi- (Sec. VI. E) transport coefficients are obtained by numer-
symmetric system, is presented in Sec. V. Fir st we ical solution of the large aspect-ratio kinetic equations.
present a refined, and slightly modified, version of the Approximate, but very convenient, interpolation formu-
small gyroradius ordering of Sec. II. In the refined ver- las are then used to connect smoothly the numerical val-
sion, the electron to ion mass ratio is ordered with re- ues of the coefficients in the various collision frequency
spect to the gyroradius, in such a way as to allow the regimes. The final expressions for the transport coef-
maximum number of transport effects to be considered ficients may be found at the end of Sec. VI (Sec. VI. F).
simultaneously (Sec. V.A). The ordering yields an ex- In Sec. VII we return our attention to the moment equa-
pansion of the drift-kinetic equation in powers of the tions, which express the conservation of particles and
(poloidal) gyroradius. The expansion is treated slightly energy for each plasma species. The flux-surface aver-
unconventionally here (Sec. V. B), in order to unify as ages of these equations, describing transport in minor
much as possible the analysis of two transport prob- radius, are seen to involve the particle and energy flux-
lems: that of a pure plasma, with electrons and one es, and the electric current, in precisely the same
species of ion, and that of a plasma with several ion forms as were calculated from transport theory in Secs.
EV-VI. Then, to obtain a closed set of equations for the U„B(x) = x '(x' 5„~ —x„x~) .
plasma fluid, the conservation laws and transport rela-
tions are combined with Maxwell's equations. The elec- The Coulomb logarithm is formally given by
tromagnetic coupling is especially interesting, because lnA = In(9N), (1 4)
the explicit form of the flux-surface average depends
where N» 1 is the number of particles in a sphere of
upon the magnetic field configuration, which can change,
in general, on the same time scale as other processes
radius &» the Debye length.
Equation (1.2) was first derived by Landau (1936).
of interest. Thus the global nature of neoclassical
Chandrasekhar (1943) used a different argument to de-
transport profoundly affects the problem of closing the
rive a special case of Eq. (1.2), which was later ex-
set of moment equations. The closed set of equations
tended, somewhat, by Cohen, Spitzer and Routly (1950).
derived here differs from previously published ver-
Rosenbluth, MacDonald and Judd (1957) derived the ex-
sions, and is intended to be more generally convenient.
pression (for e,' = e,')
However, the closure problem is treated only in axisym-
metric geometry.
C, = —1", f, ' ——
82
f, BVcf 8VB
An exhaustive treatment of neoclassical transport the- ~ Vn Vn 2 Vo Vg
ory is not attempted in this review. In particular, we
consider only briefly the large literature concerning
transport in asymmetrical. confinement systems, such where I; =4ve'InA/m, ', and
as the stellerator. [Early work on transport in steller-
ators is reviewed by Galeev and Sagdeev (19'75).J /z, =Q (1+m./nz, ) d'v'f', (v )/~v —v ~, (1.6)
Transport coefficients for the asymmetrical case are
presented, in terms of averages involving the field ge-
ometry, only in the short mean-free-path regime. Our g=g 1 d'v'f (v')Iv —v'I
discussion of rotational relaxation is similarly incom-
plete, and several topics relevant to a contaminated
It is not hard to show that Eqs. (1.2) and (1.5) are equiv-
plasma, such as charge exchange effects, are barely
mentioned. Most of the omitted topics relate, in our alent; for many applications, the latter is more conven-
ient. The functions h, and g are called "Rosenbluth po-
tentials, " because they satisfy the differential equations
opinion, to elements of the theory which, although po-
tentially very important, have yet to be thoroughly and
rigorously developed.
Although the axisymmetric theory is developed, as
much as possible, in the context of arbitrary magnetic
geometry, this is mainly to emphasize the fact that the
transport theory does not depend upon any specific as-
sumptions concerning the magnetic field configuration.
Explicit numerical results are thus given only for the where (V, )' = &'/Bv &v, etc.
large aspect ratio, circular cross section case, which We review neither the derivations of Eq. (1.2) nor the
is most generally useful. attempts to derive improved operators [see the review
by Fried (1966)]. However, some comments on the val-
C. Collision operator idity of the Fokker-Planck formalism are appropriate.
A basic element in all derivations of the collision op-
All of the transport literature under review is based erator is that, because of the long-range nature of the
on the assumption that the distribution function f, (x, v, t) Coulomb interaction, the great majority of scattering
satisfies the Fokker-Planck equation, events result in quite small deflections: Av/v« l. In
the absence of coherent wave propagation, the effect of
&f /&t+v V f, +(e, /m, ) (E+c 'vxB) &f /&v =C, (f), many such deflections is to cause each particle to per-
form a random walk in velocity space; the correspond-
where e, is the charge and m, the mass of particles of ing evolution of the velocity distribution function is then
species a, E and B are, respectively, the (macroscopic) determined by an operator having the Fokker-Planck
electric and magnetic fields, and C, is the Fokker- form of Eq. (1.5).
Planck collision operator: Because the range of the "bare" Coulomb force is in-
finite, the accumulated effect of very small deflections
8,
is divergent, unless Debye shielding is taken into ac-
count. In its simplest form, this shielding —
the result
2ve.'e', , f. (v) sf, (v') of space-charge polarization —
appears in the Coulomb
ab logarithm: the Coulomb interaction is cut off at a dis-
tance corresponding to the Debye length. Thus a more
general expression for the Coulomb logarithm is given
by (Spitzer, 196'7; Braginskii, 1965)
(1.8)
Here, a sum over repeated Cartesian indices (n, P) is where & and ~ .
are, respectively, the largest and
lmpl leds and smallest values of the impact parameter for the class of
some confinement devices of interest, this assumption with v~ being the thermal speed and
is not strictly justified, The modification of the colli-
sion operator which is required to tr~at the case 0» ~~
has been considered [see, for example, Montgomery, If the temperatures of the two species are roughly com-
et al. (1974), and references cited therein]; it appears parable, and if u; is not much larger than v, „, , then the
+mdivis to replace the Debye length cut-
that a dominant effect ratio
off by a gyroradius cut-off, in the Coulomb logarithm. v, h; /ott„-(m, /m;)'~' «1
D. Conservation laws provides a natural small parameter, with respect to
which Eq. (1.13) represents a first order expansion.
A realistic collision operator must conserve particles, The linearized electron-ion collision operator is ob-
momentum, and energy. These conservation laws are tained by assuming uq/&~; -f, and retaining only first
expressed by the following relations between moments order, O(f), terms in Eq, (1.14). We find (Braginskii,
of C,~: 1965)
(1.10) C g
—I ll;
1
— —'U(v)'
& &f,
+
2u; v
fgI (1.16)
Here the l superscript refers to the linearization.
v [m~v Cgy Cyg] = 0 The first term of Eq. (1.16), which involves only the
(1.12)
electron distribution function, describes angular scat-
tering at fixed I vI, and vanishes if f, is isotropic
is because, in collisions with the much more massive
This.
It is easily verified, by straightforward
substitution and ions, the energy change of an electron is very small (of
Integrat&on by pax'ts that the colll, sion opex'atox' of Eq. order m, /mq). Thus neither Eq. (1.14) nor Eq. (1.16) is
(1.2) satisfies Eqs. (1.10)-(1.12). Other basic proper- sufficiently accurate for calculating the rate of energy
ties of the collision operator (e. g. , that it satisfies exchange between the two species.
Boltzmann's II theorem) are verified in later sections. Momentum exchange between ions and electrons is
C;, =
9= f S v(m/2) (v'—uPC(f); (2. 9)
and the collisional change in energy flux
(1.18)
6= d'v (m v'/2) v C( f) . (2. 10)
where T, =m, v,'„, /2. In Eq. (1.18), the first term, in
square brackets, determines the rate of energy ex- Making explicit the species subscript a, we recall that
change. The last term, involving F„ is responsible for the conservation of momentum and energy in collisions
momentum conservation, as in Eq. (1.11). require s
Finally, it is convenient to present here the form of
the like-species, lineaxized collision operator. Equa- (2. 11)
tions (1.2) and (1.15) readily yield
c'. = —
'
8 /. (vlfd'u /. lv )& 8(v — v')'' Q (Q, +F. u, ) =0. (2. 12)
-sf.(vg af, (v ) Now the even (~ v~' and v') moments of Eq. (1.1), which
(1.19)
O'U8 O'Ug express the conservation of particles and energy, re-
spectively, may be written as
II. MOMENT EQUATIONS
Bn/St+ V. (nu) = 0, (2. 13)
A. Oefinitions
(8/&t) 3p/2+ V'. Q = Q+u. (F +enE) . (2. 14)
To provide a framework for later discussion, we con-
sider here the velocity moments of Eq. (1.1). First we The odd (v and v'jv) moments, which express the con-
~
define, omitting species subscripts for convenience, the servation of momentum and of energy flux, respectively,
particle flux take the form
nu, =
(&/&t) mnu+ 0 P —en(E + c 'u xB) = F,
(2. 15)
(~ /& t ) Q+ V ~
R —(3/2) (e/m) Ep —(e/m) E P ~
Transport theory of a magnetically confined plasma One further refinement is of general importance. In
has the principal object of providing a closed set of low-P toroidal confinement systems, the magnetic field
equations for the time evolution of the densities and can be decomposed into poloidal (B~) and toroidal (Br)
pressures of each species in the plasma fluid. Thus '
components such that B =B~+B» with, typically
only the even moment equations appear in the final
closed set. The odd moment equations are useful main- IB,/B. I«1, (2. 23)
ly because they alone involve the magnetic field. Hence as shown in Sec. II.C, below. Hence we introduce the
an ordering argument based on the strong field limit can poloidal gyroradius
be understood at the fluid, rather than kinetic, level = v~/Qp, = eBp/mc,
pp Qp (2. 24)
only through consideration of Eqs. (2. 15) and (2. 16).
This strong field, or small gyroradius, ordering is and let
made precise in the following subsection. It is then
~~ =(B/Bp) 6 (2. 25)
used, with Eqs. (2. 13)-(2.16), to study the confined plas-
ma equilibrium, to provide physical understanding of Thus two measures of the gyroradius may be distin-
several transport processes, and to provide expressions guished for each plasma species. We adopt the strongest
for the cross-field fluxes which are significantly more version of Eq. (2. 19)
convenient than the definitions of Eqs. (2. 1) and (2. 3).
(2. 26)
B. Small gyroiadius ordering since this applies to most of the literature under review.
The subscripts are frequently omitted for the sake of
We denote the scale length for changes in macroscopic simplicity. However, we point out that certain steller-
parameters, such as the pressure, by ator experiments do not satisfy Eq. (2. 26), and that it is
t=-IvlnPI '; only marginally satisfied on some present tokamak de-
vices.
the thermal speed by
„—
v, = (2T/m)'
'. 2. Lowest-order consequences
and the transit frequency by As shown by kinetic argument in Sec. V, the order-
ings (2. 19)-(2.21) force the distribution functions to be
~= va /I (2. 17) approximately Maxwellian
measure of a particle gyroradius (I armor
A convenient
radius) is then
f -f~+ 0(~), (2. 27)
f„=-n,(m "v~) ' exp[- (v/v, h)'] . (2. 28)
=
p —v~/Q =mcv, h/(eB) . (2. 18)
Hence in lowest order each plasma species is described
The basic ordering assumption is that this length is by a density n, and a temperature T =mvt'h/2. [This cir-
much smaller than l: cumstance is so crucial to the linear transport process-
6—= p/l = &u/Q«1 . (2. 19) es under consideration that, if it could not be deduced
from the orderings, Eq. (2. 27) would be assumed to hold
We further assume anyway. ] Some elementary consequences are that
& 1nP/& t = 0(6'(u) (2. 20) (nu, Q, F,(P —Ip), [R —(5pT/2m) l]j=o(6), (2.29)
and where is the unit dyadic. From these results, and
I
cE/(Bv, „)= 0(6), (2. 21) from the definitions of Eqs. (2. 3) and (2. 6), it follows
that
where E is the electric field. Equation (2. 20), in which
P could be replaced by n, T, ete. , serves to identify the q =[Q —(5/2) pu][1+0(&)] . (2. 30)
(diffusion) time scale under consideration. It is assumed Note also that, in view of the slow time variation de-
that all faster processes have evolved to equilibrium. scribed by Eq. (2. 20), the electric field is predominantly
Equation (2. 21), called the dhift o~dexinI„serves to rule electrostatic
out certain rapid fluid motions associated with the case
cE/(Bv, „)=0(1) (magnetohydrodynamical ordering). (c/Bv, „) (E+ VC) =0(&'), (2. 31)
where 4 is the electrostatic potential.
Ref inernents We now turn to the odd moments, Eqs. (2. 15) and
Because of the small electron-to-ion mass ratio (2. 16). Using Eq. (2.29) to identify the zeroth-order
(m. /m;) &10 ', terms, we f ind
B [VP+enV4] =0(&),
it is occasionally useful to include the species sub-
scripts; we assume B [v(PT)+ePve ] =0(6) .
Hence
so that B vT =0(6), (2. 32)
6./6, - v,„;jv ~, —(m, /m;) ' ~' «1 . (2. 22) B V[nexp( —eC/T)] =O(6) . (2. 33)
In combination with the quasi-neutrality condition 271 upon one traversal of the minor (major) toroidal cir-
cumf erence.
g e. n. =o, (2. 34) (iii) 0 and g are chosen such that the ratio B VO/
B V'f depends only on 4.
(iv) 4, 6 and f are chosen such that the metric de-
Eq. (2. 33) yields the constraints terminant
B Vn=O(&) =B VP, ~
(2. 35) 0 =-
ve (ve x vg) '
g(@ g) I I (2. 43)
B VC =0(5) . (2. 36)
is independent of &
Kul s rud1958, ) .
Gra, d (1967) has emphasized that the equilibrium rela-.
(2. 45) and (2.46), one can express B in terms of its flux
functions
tions need not possess solutions in general. The exis-
tence of solutions can be guaranteed only for the axisym- B =(2v) ' VQ XVO+V/XVg (2. 48)
metric case and, in an approximate sense, for near
axisymmetry. We therefore restrict our attention to Here, the first term is called the toroidal field B2, and
the second term the poloidal field B~. These fields are
configurations in which any asymmetries are sufficiently
tangent to the (not necessarily orthogonal) & and 0 coor-
weak to permit the existence of nested flux surfaces.
dinate lines, respectively. An equivalent expression is
a. Flux coordinates —
B& —Vf X V|tp, Br —
—qVJX Ve, (2.49)
Figure I shows a typical flux surface, the field line where the quantity
from which it is constructed, and the magnetic axis.
The directions of variation of so-called "flux coordi- q(0) =(2~) 'dy /de (2. 50)
nates" (4, 0, f) are also shown in the figure. The func- is the reciprocal of the rotational transform, and de-
tions 4 (x), e(x), and g(x) have the following properties. scribes the average pitch of the field lines: g is the
(i) 4', the radial coordinate, is constant on each flux number (not rational, in general) of toroidal circuits
surface. for one complete poloidal circuit, of a field line around
(ii) e(p), the poloidal (toxoidal) angIe, increases by the magnetic axis.
A similar argument yields, for the parallel heat flow, is called the classical flux, and (Wimmel, 1970)
9')), = —(5/2) (mQ) 'IPdT/dy+I. (g) B . (2. 68) nuNc=- (mn)-'nx [v' (P —ip) —e(nE +zz c)] (2. 78)
While the functions K(P) and I(g) remain to be deter- is called the neoclassicaL flux. (This terminology, which
mined, Eqs. (2. 66) and (2. 68) are useful because all the is roughly conventional, assumes that both fluxes are
(0, f) dependence is given explicitly in terms of the field finally evaluated under the assumption of ordinary Cou-
geometr y. lomb collisions, as discussed in Sec. 1.)
Equation (2. 66) also provides the first-order current
J =K(HALI)B —c(dP/dg) g' Vrjr x VO 1. Classical diffusion
or To interpret the fldx nu„we may observe from Eq.
j~, = jC(g)B —c(I/B) dP/dP, (2. 70) (2, 61) that the gyration induced diamagnetic drifts for
oppositely charged species are in opposite directions.
where The resulting counter-streaming between ions and elec-
trons is accompanied by collisional friction,
If =+ e, IC. (2. 71) F~, = v„m, n, (u,„~ —u, , ) where v„. is the Coulomb colli-
sion frequency for electron-ion momentum exchange.
is a function of g only. Note that Eq. (2. 69) can also be (An additional contribution to F~, coming from the tem-
derived from Eq. (2. 40) and the quasi-neutrality condi- perature gradient, is irrelevant to the present, qualita-
tion tive, argument. ) Equation (2. 77) then yields a radial
flux nu, = —(v/Q)(mQ) '&P, and correspondingly a, clas-
(2. 72) sical diffusion coefficient (Rosenbluth and Kaufman,
1958)
3. F irst-order stress
Equation (2. 14) can be considered as the trace of the
more general moment
electron
', (2. 79)
gyroradius defined by Eq. (2.18).
Equation (2. 79) (which happens to be exact when T, =T,
has been derived for the case of electron diffusion, but
in fact in the case of a single ion species, it holds for
).
Without writing Eq. (2.73) any more explicitly, we note the ions as well. This is because collisional momentum
that it can be reduced to the form [see, for example, conservation, Eq. (2. 11), guarantees ambipolarity of
Chew, Goldberger and Low (1956)] the classical flux
11 x P + P x B = 0(5') .
Hence, in a coordinate system with one axis along B. the
g e, n, u„=0. (2. 80)
first-order stress must be diagonal. Notice that Eq. (2.61), and therefore Eq. (2. 79), re-
quire only
P =R~Pii + (~ —nn) P + 0(&'), (2. 74)
v&&Q,
where P, [ and P are, respectively, the parallel and per-
pendicular stresses. Equation (2.29) requires the stress i.e. , that the l. owest-order gyromotion be unperturbed
anisotropy to be at most first order in ~. by collisions. As shown in Sec. III, such nonrandom
This stress tensor describes a plasma in which the particle motion drives a non-Maxwellian perturbation
various directions of velocity in the plane perpendicular on the distribution function, f
—f~ = 0(&), and this per-
to B appear uniformly populated, when viewed on time turbation in turn drives F~.
scales long compared to the gyroperiod.
A similar argument yields, for the tensor R 2. Neoclassical transport
(2.75) In our definition the neoclassical flux has two sources,
P —ip and re +n «'.
The latter pertains mainly to mo-
D. Second-order floxes tion of the magnetic surfaces, and is discussed in Sec,
As suggested by Eqs. (2. 20) and (2.63), motion of par- II.E below. Here we consider the typically more impor-
ticles and energy across flux surfaces appears as a sec-
tant source, P —lp. It's form and origin depend upon the
ond-order process. While closed form expressions for size of the collision frequency: ahvays assuming v «Q,
the radial particle and energy fluxes can be obtained only we must also specify the relative magnitudes of v and ~.
In the collisional regime
by solving the kinetic equation, useful information is
nonetheless available at the fluid level.
Consider first the particle flux. Using Eqs. (2.37) and
(2.61), we may write the stress anisotropy is kept small by collisional ran-
domization,
nu~ =nu» +nu, +nuNc+ 0(&'), (2. 76)
(2.81)
where
but the mean free path is short enough to allow pressure
— '~xp
~u, = —(m~) (2. 77) variation
F. L. Hinton and R. D. Hazeltine: Theory of plasma transport
dy
(P') (2.91)
are the classical and neoclassical contributions, re-
where
spectively.
F ~ VV = O(52), (2. 92)
E. Axisymmetiic systems and we have used Eqs. (2.57) and (2. 74). The remaining
The transport theory for a rigorously axisymmetric terms in Eq. (2.15) yield simply
torus is relatively complete, not only because of its
(nu Vg) = —(c/e) (R2Vg (F + enE)) + O(5') . (2. 93)
comparative theoretical simplicity, but also because this
is the more relevant case: asymmetry is purposefully This expression contains both classical and neoclassi-
minimized in the design of most present confinement cal contributions. These may be distinguished by noting
systems. Axisymmetry provides essentially two sources from Eq. (2.77) that
of simplification in. our previous formulae. The first is
&nu. . Vy& =-(c/e)&R'Vg F ), (2. 94)
related to the axisymmetric version of Ampere's law,
and pertains mainly to the plasma equilibrium. The sec- and therefore that
After substituting the conservation laws, Eqs. (2.13) and Note that in Eq. (3.1) we measure p, in a fixed reference
(2. 14), into Eq. (2. 111), and using Eq. (2. 30), we obtain frame, rather than in one moving with the E x B drift.
This reflects our choice of the drift ordering, Eq. (2. 21),
e, =-u, VP, /T, —q, VT, /T', rather than the magnetohydrodynamical ordering, cE/II
+ Q, /T. +u. (F, + e.n, E)/T.
~
. (2.112) '
vg„. Note al s o that the dr if t order ing, with Far aday s
law, requires only d lnB/dt = O(«u); but the stronger
Therefore, using Eqs. (2. 11) and (2. 12), we have
assumption
~ =-gfu, VP, +q, VT/T I+5 E. (2.113) d lnB/dt = 0(&'cu)
a
which is implied by Eq. (2.20), is more relevant to
In this form, the irreversible nature of ~ is apparent:
transport theory, and is adopted here.
Eq. (2. 113) has the characteristic form of a sum of The equation of motion
products of forces (VP„VT,/T„E) and fluxes (u„q„T).
Next, we specialize to an axisymmetric system, and
consider the flux-surface average of &. After making [E (x, I ) + c ' v B (x, I)] && (3.6)
use of Eqs. (2. 69), (2.86), (2. 89), and (2. 93), we obtain
can be written in terms of p, , e, and &, with the result
{K) = —g((q, Vg) d lnT. /dP+ (c/e, ){RIr. )n. 'dp. /dy} (3.7)
a
—(e/4~){E a}dI/dy. (2. 114)
dc e d4
Here, the first term requires no interpretation. The (3.8)
second term involves the radial particle flux, relative
to the moving poloidal flux surface, as shown by Eq.
(2. 108). To understand the last term of Eq. (2. 114), note
=Q+e, . +, +XV
Vj
v~,
dt
—
m
E, (3 9)
that the Poynting flux relative to the poloidal flux sur- where, on the right-hand sides, d/dt=— s/st+v V. ~
For small gyroradius, the particle velocity, v =v~~+v, After substituting Eq. (3.11}into Eq. {3.14), and using
is most conveniently parametrized in terms of the mag- Eqs. {3.7) —(3.9), we find
netic moment —
—n(v()+ +vg, (3.16)
vgc v())
g = v'/2a, where
the energy v I) = (v~/2Q) n ' V x n, (3.16)
e = v'/2+e4/I, (3.2)
and the gyrophase angle, v~ =cE &n/II+n&& (/J VR+ v~~&z V~Ã+ V~~ Bn/at)/Q (3.17)
=-tan '(v e, /v e, ),
g
~ ~
(3.3) 2A lucid treatment of the averaging is presented by Banos
where the unit vectors e, and e, are chosen to make (1967) who however, uses the magnetohydrodynamical order-
(n, , e„e,) a right-handed orthogonal triplet at each point. ing. More extensive treatments of the g. c. theory Inay be
Thus found in, for example, Northrop (1963) Hastie et al. (1967),
and Morosov and Solov'ev (1966); the last of these is most
v = (2pB)'~'[e, cosf —e, sing]. (3.4) relevant to transport theory.
are, respectively, the parallel and perpendicular g. c. formally similar to Eq. (3.18). To pursue this analogy,
drifts. The first term of Eq. (3.17) is the familiar E x B one presumes the existence of a function 4~(x, p, , e)
drift, the second term is called the gradient-B drift, s atisf y ing
and the third term is called the curvature drift. The
third and fourth terms together comprise the "accelera-
Bg. VC'~ =0, (3.28)
tion drift, " in the terminology of Northrop (1963). The analogously to 8 V4 =0. The surfaces corresponding
various perpendicular drifts may be understood in terms to fixed values of p, , e, and 4 „are called drift surfaces;
of distortions of the gyro-orbit; see, for example, Eq. (3.28) states that the first order g. c. motion is con-
Rosenbluth and I ongmire (1957). fined to a single drift surface, just as the zeroth-order
Equations (3.15)—(3.17) are exact. The small & as- motion is confined to a single flux surface.
sumption enters in the statement that when viewed on While the family of trajectories 8+(p. , e) are generally
time scales long compared to the gyroperiod, the exact well-defined, the existence of drift surfaces ean be
particle trajectory may be approximated by v, . guaranteed only for confinement geometries possessing
Thus the lowe st- order g. c. mot ion consists of a appropriate symmetry. Then 4+ appears as the dynami-
streaming along the magnetic field lines. cal invariant, in a g. c. Lagrangian sense, corresponding
to the ignorable coordinate (Morosov and Solev'ev,
vg, =nv I, [1+ O(&)], (3.18)
1966). In an axisymmetric toroidal system, for exam-
where ple, 4+ is realized as the gyrophase-averaged, canoni-
v II = ~[2 (e —p. B —e 4'/m) p~' (3.1 9)
.
cal angular momentum. We consider this case in further
detail.
In a nonuniform field, even this motion displays inter-
esting features, because Eqs. (3.7) and (3.8) imply C. Axisymrnetric systems
dy/dt =-(Iu, /B) aB/at, (3.20) The predominant magnetic field variation in an axisym-
dc/dt =(e/m)[a4/at — c 'vIIn aA/at], (3.21) metric (or nearly axisymmetric) toroidal system re-
sults from
i.e. , that the averaged changes in and e are relatively jU,
vn ~
Vy = I v ~ ~ V(v~~/Q)+ 0(&') .
~~
(3.32) '+ell —untrapped orbit
~ VllO ~O~
The agreement between this expression, which holds for
any P, and the axisymmetric version of Eq. (3.24), which Magnetic surface
assumes P = 0(&), may seem surprising. But Eq. (3.24)
omits only a term proportional to n&& VP, as shown by
Eq. (3.22), and axisymmetry implies Vg nx VP
In VP = 0(6), for any P, as noted in Eq. (2. 35).
-~
ell- trapped orbits:
A more interesting result follows from the exact in-
variant of Eq. (3.30): { Vila & o~
~ Vila»~
dg/dt = (cm/e) d(Rv r)/dt, (3.33)
where the derivative is taken along the particle trajec-
tory. Considering only the trapped particles, we inte- FIG. 2. Trapped and untrapped guiding center orbits, in an
grate Eq. (3.33) between two successive times t, and t, axisymmetric system, projected onto a constant-& plane. The
at which ~ ~ vanishes, and obtain dashed line is the projection of the magnetic surface to which
each orbit shown is tangent; the parallel velocity at the posi-
tion of that magnetic surface is labeled v()p, The signs indi-
g(t, ) —y(t, ) = (sq/~t+v Vq)(t, —t, ) =0, (3.34) . cated for v correspond to a positively charged particle,
~~0
when the toroidal magnetic field and current are in the direc-
where the overbar denotes an orbital average. Equation tion out of the paper. The designations "well-untrapped" and
(2.106) then provides the result "well-trapped" refer to the cases (v~~o && ~v(DoB/B)' and )
)v~~&
~
«v ~o(DB/B)t'~, respectively, while "marginally
v cREr ~—
Vg =u, Vg = (3.35) trapped" refers to the case )v~~0 ] R v~o(AB/B)'~ .
effect"), which bears significantly on transport in a where Q~ =eB~/mc is the poloidal gyrofrequency, bv, ~
tokamak, was discovered by Ware (1970) and indepen- is to be determined from Eq. (3.19), and we have ap-
dently by Galeev (1971). proximated, for large aspect ratio
In Sec. II we pointed out that, without collisions, the
induced, parallel electric field must be presumed to ~t II/ai=t R/R =a/R«1, (3.37)
vanish [cf. the discussion concerning Eq. (2. 109)]. In where a is the minor radius of the torus. Equation
this case the velocity u, coincides with an E &&B drift in (3.36) further simplifies in the following two regions of
the induced electric field [recall Eq. (2. 98}], and Eq. phase space. (i) For the far untrapped particles with
(3.35) appears trivial. It is not trivial because it re- small p. , Eq. (3.19) yields ~Av~~~ vthQ/R, whence
mains approximately correct in the presence of colli-
sions, when. u, is typically much faster than the E &B zr = p~(a/R), v~ «v~~, (3.38)
drift, provided the collision frequency is sufficiently where p~ is the gyroradius in the poloidal field. (ii} For
small (Rutherford, et at. , 1970). This surprising cir- the trapped and nearly trapped particles we find 4m~~
cumstance can be thoroughly explained only from the —vth(a/R)'I' and therefore
kinetic considerations of Sec. V. However, a partial
explanation is possible here. Note first that, when act- ar = p~(a/R)'/' v —v . (3.39)
ing on the untrapped particles, the (parallel components
of the) accelerating toroidal electric field and deceler- Hence these particles make the largest excursions.
That the g.c. excursion is proportional to p~, rather than
ating toroidal friction force tend to balance. Untrapped
particles therefore experience, roughly, only the rela- p, was first noted by Tamm (1 959).
The trapped particle g. c. orbits have been named
tively slow E&&B drift; see Eq. (2.109). For the trapped
banana orbits (Furth and Rosenbluth, 1969), for reasons
particles, on the other hand, the parallel electric field,
evident from Fig. 2, and the trapped (or nearly trapped)
which always acts in the same direction, cannot be ba-
lanced by the parallel friction force, which changes
particle excursion width
direction each time the trapped particle is reflected. b. , = p~(a/R)'/' (3.40)
The kinetic analysis of Sec. V shows that for these par-
ticles the effect of friction, when averaged over the is called the banana width.
collisionless orbit, is so small that Eq. (3.35) remains
D. Dif fusion
roughly valid.
Equation (3.31) may be used to estimate the radial ex- The preceding discussion neglects the microscopic
cursion of a g. c. from a given flux surface. I, et 6 de- f Coulomb interaction between particles. Such "colli-
note the change in f
which occurs on the lowest-order sions" change each particle's energy and magnetic mo-
orbit. Then A4~ =0 implies that the poloidal flux change ment, causing diffusion in both velocity- and coordinate-
associated with the excursion is A/ =Id, (v~~/Q); the ex- space. As noted in Sec. II, the character of this diffu-
cursion in minor radius, r, is Ar =b. g/~VQ~ or sion depends upon the size of the collision frequency.
length" qR (where q is the safety factor and R the major of trapped banana orbits (superbanana diffusion). How-
toroidal radius). Therefore q'R' = (v,'„/v, )i; with ever, with the neglect of such extremely collisionless
u~ = pu, „/R, Eq. (3.42) provides the Pfirsch —Schliiter effects, the classification outlined above is quite gen. —
coefficient erally adequate [see, for example, Connor and Hastie
D q2p2 p q2D (3.43) (1974)] .
We now estimate (neoclassical) diffusion coefficients
for these two regimes, specializing first to the axisym-
2. Col l I sl on I ess reg I m Bs metric, large aspect ratio case.
In the banana regime, the dominant contribution to
Let 8 be the pitch angle of the gyro-orbit
diffusion comes from the trapped and nearly trapped
sin'0 = 2pB/u'. (3.44) particles, whose g. e. 's make the largest radial exeur-
In a random walk description, the change in 6 due to
sions. In Fig. 4 we show two electron orbits, both of
which pass through the point E'. (Although banana orbits
Coulomb scattering is proportional to the square-root of
the time. Since the classical collision frequency, v, ,
are depicted, the following argument applies equally to
nearly trapped particles with small v~~. ) Because the
pertains to roughly 90 changes in 0, the collision fre-
orbits are being traversed in opposite directions, one
quency for some arbitrary change AO is
lies inside the indicated flux surface and the other lies
(3.45) outside. It is clear that a reversal of the g. c. parallel
For example, the collision frequency for scattering a velocity at I' induces a radial step of the orbit center by
a distance -~~. From Fig. 3, we see that pitch angle
diffusion by as little as Ae = (a/R)' ' is sufficient for
magnetically trapped particle onto an untrapped orbit is
roughly
such reversals, which therefore occur with frequency
v, = (B/aB) v„ (3.46) Hence the trapped particle diffusion coefficient is
D = (a/R)''v a'
(~/R )1/2 p p2 (3.53)
where the second form follows from Eqs. (3.40) and
(3.46). The banana coefficient exceeds D, by as much as
two orders of magnitude in typical tokamak devices, be-
cause p~» p. FIG. 3. Trapped and untrapped regions of velocity space, in
In the plateau regime, collisions are too frequent for terms of the parallel (v!!) and perpendicular (vj} velocities.
undisturbed banana motion, but the collisionless g. c. The trapped region is shaded; it occupies the region outside
of a cone whose angle depends upon the field variation,
motion of particles with smaLL v, ! remains significant. as indica, ted.
Consider the class of "resonant" particles with v!!& v„
where the speed uo is to be determined (Stringer, 1970).
From Eq. (3.45), the effective collision frequency for v~/v„, and they comprise a fraction (b/R)'~' of al]. the
particles. Hence our usual argument yields the diffu-
scattering particles out of this class is roughly sion coefficient
v(v, ) = (v, „/v, )'v, . (3. 54)
(3.58)
v(v, ) & v~~/qR, (3.55)
appropriate to the banana regime
then resonant particles will suffer effective collisions
before they are able to sample upper and lower halves
of the torus. The uncompensated radial drift of such in an asymmetric system (Gibson and Mason, 1969;
particles, between effective collisions, leads to a ran- Galeev, et al. , 1969; Frieman, 1970; Hazeltine and
dom walk with step-size v~/u(U, ). The corresponding Rosenbluth, 1972; Connor and Hastie, 1974). As re-
diffusion coefficient is marked previously, diffusion across the asymmetrical
(3.56) confining field of a stellerator is more complicated, in
that several low collision frequency regimes can be dis-
since only a fraction v, /u, „of all the particles are reso- tinguished. A detailed discussion of g. c. orbits and dif-
nant in this sense. Finally, Eqs. (3.54) and (3.55) imply fusion in stellerator geometry may be found in the re-
u, /qR = v, (v, „/v )2, so that Eq. (3.56) reduces to the view of Kadomtsev and Pogutse (1971), and in the re-
plateau diffusion coefficient view of Galeev and Sagdeev (1975).
A serious disadvantage of asymmetry, with regard to
Dp = U~/(dr (3.57)
transport in a high temperature, nearly collisionless
so-called because it it independent of v, . plasma, is evident from the &, ' behavior of D~.
In experimentally realistic geometries true plateau
behavior does not occur, because the necessary inequal- E. DISCI'lbUtlQA fUAGtlOA
ities, (3.52), cannot simultaneously be satisfied very
strongly unless the aspect ratio is extremely large. 1. Velocity coordinates
Thus, for typical aspect ratios, the diffusion coefficient It is generally convenient to express the distribution
increases monotonically with &, and only a mild flat- function, f, in terms of the velocity variables (p. , e, &),
tening is observed in the regime of Eq. (3.52) (Hinton defined by Eqs. (3.1)-(3.4). Some technical properties
and Rosenbluth, 1973). The plateau argument nonethe-
of the transformation v- (p. , e, g) are considered here.
less provides a useful estimate for this regime. First note that the variables (g, e, g) do not determine
3. Asymmetry
magnetic field perturbation 5(&) «R
1
An asymmetrizing
l
has no effect on the form of the diffusion coefficients, I flux
as estimated above, for the collisional or plateau re- I surface
gimes. However the banana formula, Eq. (3.53), depends l
upon the periodic nature of the collisionless orbits, !
which pertains only in the axisymmetric case. In gener-
symmetry
al, we must consider particles trapped by local maxima
of b. Because the g. c. orbits of such "localized" parti-
cles (Gibson and Mason, 1969) may be confined to the
upper or lower halves of the torus, they drift away from FIG. 4. Two banana orbits, both of which pass through the
the magnetic surface, between effective collisions, in. a same point, &, on a flux surface, but which differ in the di-
manner similar to resonant particles in the plateau re- rection in which they are being traversed, at H. A small-
gime. The effective collision frequency for localized angle collision at 2, .which causes a change in the sign of v!!,
particles is v~ = (R/b)v„ their step-size is roughly leads to a radial step of the orbit center.
v uniquely, because the sign of v~~ is not specified [cf. surface, i.e. , ~, depends only on 4. Explicitly,
Eq. (3, 19)]. One must therefore introduce the additional, (3.67)
discrete variable
where R is the smallest distance between the sym-
(3.69)
metry axis and the flux surface under consideration. The
Evidently v is a constant of the motion for untrapped par- size of the trapped region still depends upon position,
ticles, but not for trapped particles, whose parallel through A. . The bounce surface, S~(A), consists of points
velocities change sign at each "bounce. " The signifi- . such that
cance of this fact can be appreciated after the trapped
R =XI . (s.66)
and untrapped regions of velocity space have been more
explic itly distinguished. Finally, we take note of the Jacobians for the two
In general, the distinction between trapped and un- velocity- coordinate systems prese nted above:
trapped regions must be made relatively to a specific
field line, about which a particle with given (p, , e) is d v=B djLdcdg Ull (s.69)
gyrating. (The gyrophase f is clearly irrelevant to ques- G
tions of trapping. ) If q„or c —p, B(x) —eC'(x)/m, can vanish where the integration domain is given by eC'/m & e &~,
for points x on a chosen field line, then p. and e are in 0& p. & (e —eC/m)/B; and
tQe trapped region relative to that field line. The points
at which oil vanishes will be called bounce points, and d'v =Bg dA. dude (3.70)
denoted by x, . The set of all bounce points, correspond-
ing to different field lines but to the same values of p, where the integration domain is given by 0 & u & ~,
and c, is called the bounce surface and denoted by 0&A. &A .
S, (p. , e) —
= (xl e —p, B(x) = e4 (x)/m}. (s. 60)
f
Since the value of at points x, on the bounce surface 2. Gyrophase dependence
measures the number of particles being reflected at x„
we must have
We now consider the dependence of upon the gyro- f
phase angle g. Recalling the gyrophase average of Eq.
f(x„p., e, v=+1)= f(x~, p, , e, (x= —1), (s.61) (3.13), we define a function by f
for all (p, , c) in trapped regions. Equation (3.61) may be f(x, p, e, g, t) = f(x, p, , e, t)+ f(x, p. , e, g, t) . (3.71)
considered as a boundary condition, in the sense that,
for fixed p. and e, the bounce surface constitutes a An approximate expression for is easily derived. We f
boundary of the domain of definition of f.
Of course the
note that, in the presence of rapid gyration, can satis- f
fy the assumed orderings,
v dependence of f
is not constrained in untrapped re-
gions. S lnf/St «0, lV infl«p-', (3. 72)
We have already noted [Eq. (2. 36)] that the electro-
only if it is nearly constant on the gyrotrajectory
static potential is nearly constant along B. For particles
with large ionic charge, g»1, the variation of potential f(x, p, &, K, t) =a(x —p(K), p. , &, t),
erterI, y along field lines remains important (Hazeltine,
where x —p =R, the position of the g. c., and p is given by
et al. , 1974), but more typically,
Eq. (3.11). A first-order Taylor expansion of the func-
v„V(e@/m) «vg, V(p, B) . ~
(s.62) tion g provides f
=g(x, p, , e, t) and =-p V But f ~
f. f
must be nearly Maxwellian, so we have
In this case the kinetic energy is a lowest-order con-
stant of the motion, and the variables (&, w, r), where f, = p'Vfu+O(6— '). (3.73)
A, = p. /w, l8 =E' —e4/m = v /2 (s.63) Later in this section, we outline the derivation, from
are most convenient. Since kinetic theory, of the more accurate expression
v = 0'[2 to (1
~[
—XB)]' (s.64) f = —p Vf+p Bf
8ji B
—pVB+v
1
ll
n'Vn+v ll
BB.
at
—e E—
nZ
the trapped and untrapped velocity space regions are
given by ~ »»,
and 0 & ~ & ~, , respectively. Here +
Bf e sA
96 Rzc Bt
(s.66)
is the largest A, for which f(x, A, , w) is defined, and the (s. 74)
critical & for trapping is
where A is the vector potential, p =p/p, and v =v /v
~c=Bmax ~
(s.66) It is instructive to compute the perpendicular particle
where B,
„ is a local maximum of B occurring on a given flux from Eqs. (2.1) and (3.74)
field line.
As an example, we consider the axisymmetric geom- d'vv~ f= d'vvi f. (3.75)
etry of Eq. (2.87), with B =Br. The ergodic property of
field lines then implies that the distinction between After considerable manipulation, using Eqs. (3.4),
trapped and untrapped regions is uniform on each flux (3.61), and (3.69), one obtains
d'v f (vD+v~, ) =cnE xB/B'+ (c/eB)(P (nx VB)/B netization can enter the energy conservation law, as
shown below.
+ P))(V x n)
Finally, it is instructive to consider the lowest-order
+ P~(Vxn)~~]. ,
flow. Substituting Eq. (3.73) into Eq. (3.75) we quickly
find
(s. 78)
where (s. as)
(Vxn)~ =nxKVxn)xn] =nxn Vn. ~
(3.79) where the right-hand side is the diamagnetic (and elec-
We interpret the perpendicular fluxes of Eqs. (3.76) trostatic) drift given by Eq. (2.61), and we have recalled
and (3.78) by comparison with the rigorous flux, nu Eq. (2. 36). Thus, as remarked in Sec. II, this first-
obtained from the moment equation, Eq. (2. 37). Through order flow is a result of I armor gyration. Similarly,
O(&'), the latter may be written as
Jl d'v(vD +v~, ) fu =nu~, + V x M, , (s.a4)
nu~ =cnE x 8/B'+(c/eB)7P~(nxVB)/B+ P~Vxn
+ (P~~ P )nx (n V)n]' where M, =cjn/eB is the lowest-order magnetization.
We next turn our attention to the energy flux, Q~. By
+nu, —VxM, (3.80) manipulations similar to those used in deriving Eq.
where nu, is the classical flux given by Eq. (2. 77), and (3.81), we find
we have used Eq. (2. 74) for the stress. From Eqs.
(3.78)-(3.80) we conclude d'v f (m v'/2)(vv + v ~~) = @~ —Q, + V x (mcBR ~/eB')
d'vv~ f =nu~ —nu, . (s.82) terms on the right-hand side of Eq. (3.85) are closely
analogous to those on the right-hand side of Eq. (3.81),
Hence the g. c. results, Eqs. (3.15) and (3.74), are not and only the last term, which apparently describes an
sufficiently accurate to describe classical collisional electrostatic transport of energy, requires comment. It
transport. This is not surprising, because the collision enters the energy conservation law through its diver-
frequency does not appear explicitly in either Eq. (3.15) gence; notice that
or Eq. (3.74). The O(v&') terms in which produce nu,f cP, nx V4/B = ed@ V xM,
V ~ ~
(3.86)
have been evaluated by Rosenbluth and Kaufman (1958),
for the case of a uniform magnetic field. Alternatively, using Eq. (3.77). Hence the last term of Eq. (3.85) takes
the classical transport coefficients may be evaluated into account electrostatic energy changes, enu V4, ~
from moments of the collision operator, Eqs. (2.77) and associated with the plasma magnetization. From Eq.
(2.84). In Sec. IV, we show that this method yields the (2. 59) we observe that the lowest-order potential C'(4)
same results, and requires knowledge of only the low- cannot contribute to the flux surface average of Eq.
est-order f, given by Eq. (3.'73). (3.86). In fact, in most theoretical treatments the
On. the other hand, neoclassical transport is fully in- poloidal variation of @ turns out to be so small that the
cluded in Eqs. (3.81) and (3.82); the relevant collisional last term of Eq. (3.85) is omitted.
effects enter implicitly, through non-Maxwellian correc- The lowest-order energy flux is
tions to f.
In other words, while classical transport,
which results from collisional interaction with the gyro- d'v va+v~j ~= ~&+V& 2 3 87
motion, is neglected in Eq. (3.74) for neoclassical f,
transport, which results from collisional interaction where Q» is given by Eq. (2.62), and we have recalled
with the g. c. motion, can be included, in principle, by Eq. (2.29).
sufficiently accurate calculation of f. In conclusion, we have shown that the solution to the
Eq. (3.4), and dp/dt, de/d, t and dg/dt are given by Eqs. f/St+ (v, n+v„) Vf +
i
~
n ~
V
B
(3.7)—(3.9). Note that all terms in Eq. (3.88) are at most
of order =&0, except the 0 term in dt/dt
cu To em. pha- + (de/dt)g, Bf/Be =Q(f, f), (3.95)
size this, and to simplify notation, we introduce the where
operator ~*, defined by
(de/dt)g, =de/dt+ p. ea/—
at
Bf s f
djl Bf de
( ) = (e/m)aC /st+ BR/at —(e/mc)v~~ n BA/Bt,
Bt 9p, dt
~
p
(3.96)
using Eq. (3.21).
Equation (3.95) is easily interpreted. The coefficient
Then we have
of Vf is the g. c. velocity; the coefficient of af/eg dis-
'dfls& =6 * f (s. 89) plays the well-known (Hastie, et a/. , 1968) departure of
p, from the exact, adiabatically nvariant, magnetic mo-
~,
which implies
ment; and the coefficient of Bf/Be, given by Eq. (3.96),
5gf-0 (s.90) is the rate of change of energy of a g. c. with magnetic
moment p. and velocity v„=v~~[1+ 0(&)]. Thus appears f
in Eq. (3.95) as the phase-space distribution of g. c. 's.
9
Qgf Qgf To further simplify the equation, we begin by observ-
(3.91)
Bg ing that the coefficient of Bf/Bp, is second order in &,
since n Vxn =n Vxl3/R- p-6. From consideration of
~
Eq. (2. 31), which implies and by requiring Eq. (2.27) to hold,
we find that
(e/mc) v([n ~
BA/Bt = 0(&(o u,'„) . (s. 98) f —
=f~(4, &) +g(0, p, &) (»[[/&)Bfd(/Bp, (3.104)
With regard to the BB/Bt term in Eq. (3.96), it must where the function g, which must be constant on flux
be noted that our previous estimate, Eq. (3.5), does not surfaces, represents an 0(&) correction to the Maxwel-
take 'nto account the low-P assumption of Eq. (3.94). lian.
Thus, restricting our attention to situations in which It turns out that Eq. (3.104) accurately describes the
changes in the external current occur much more slowly first-order distribution function even in the presence of
than changes in the plasma current, and noting, from collisions (Rutherford, 1970). Of course, solution of
Eq. (2.40), that the plasma current contr'ibutes to B only the first-order drift-kinetic equation, including colli-
in order P, we have sions, is required to determine the function g.
- pB Inp/Bt The main content of Eq. (3.104) is that the first order
B lnB/Bt = 0(p&'cu) . (3.99) distribution function. can depen. d upon the poloidal angle
Hence for second-order accuracy, the BB/Bt term in 0 only through the v[[/0 factor appearing in the last
Eq. (3.96) may be neglected. term. In this sense, Eq. (3.104) is similar to Eqs.
In axisymmetric systems, a more explicit and instruc- (2.66), (2.68), and (2.88) for the first-order flows,
tive derivation of Eq. (3.99) is possible. We observe which were derived from consideration of the moment
from the equilibrium relation, Eq. (2. 90), that equations. Indeed, after substituting Eq. (3.104) into the
definitions, Eqs. (2.1), (2. 3), and (2.6), of the particle
B', = ~Vy/R('-&, and heat fluxes, and using Eq. (3.69), we find that Eqs.
whence (2.66) and (2.68) are reproduced; the flux functions K(g)
and 1.(g) appear in this way as integrals of g (Hazeltine,
(3.100) et al. , 1973)
Since B~, which is driven mainly by external currents,
is nearly constan. t, we have K (d ) = 2 n gaJ d d d E d, (s.lo5)
= Br —
(1 + B~~/B~~) —B (Bq/B2r ), — /(d) = 2m
got ddd~ m[s —ed/m —(5/2)T/m]d
«
assumes v, „/v I, and, in the sense described at the end mentum of that species. Thus only the unlike-species
of Sec. III, is nearly local. Thus the flux-surface geom- collision operator enters Eq. (4. 1). For simplicity, we
etry must be considered (to obtain the average particle consider the case of a single ion species, and calculate
or energy flux across a surface) only at the fluid, rather the friction force exerted on the electrons
than at the kinetic, level. This partial decoupling of
collisional and geometrical considerations is in sharp F„=Jl d'vm, v, C„(f„f,). (4. 3)
contrast to the kinetic theory for collisionless regimes,
in. which a relatively complicated, flux-surface-aver- As noted in Sec. II, to compute the flux through 0(&'),
aged, version of the drift-kinetic equation must be only the 0(&) form of F is required. We therefore ex-
solved. pand
Hence the perpendicular transport problem for the
collisional regime is comparatively simple, and only in fa = fa~ + fai + ' ' ' &
this regime have exact transport coefficients, for arbi- where f, ~ is a Maxwellian, as in Eq. (2. 28), and f, = 0(5)
trary toroidal geometry, been calculated. can be decomposed into its gyrophase-averaged and
For convenience, we begin this section by calculating gyr ophase- dependent parts
the classical perpendicular transport coefficients. How-
ever, it should be recalled that classical transport al- f$ f] fj
ways occurs (so long as Q~& v), and that the same clas- The corresponding expansion of the collision operator
sical results pertain in all collision frequency regimes. is written as
We next consider, in Sec. IV. C, the classical Parallel
transport problem. While this topic has definite intrin- c., (f., f,) = c„(f,„,f, )+c'„(f,)+ 0(&'), (4.4)
sic interest, it would be beyond the subject of this re- where C' denotes the linearized operator
view if it did not bear on the problem of neoclassical,
perpendicular transport. Indeed, from the discussion c.', (f, ) = c.,(f, ~, f~, ) + c., (f.„f, „). (4. 5)
of collisional transport in Sec. III, it is evident that the
neoclassical (Pfirsch —Schluter) enhancement of colli- The first term of Eq. (4.4) vanishes in general only if
sion-dominated transport depends upon diffusion and heat all plasma species have the same temperature
conduction in the direction of the magnetic field.
Finally, in Sec. IV.D, we combine the results for the
C.~(f.~, fute) = (&. —7'&) &(I vl) . (4 6)
parallel transport problem with some conclusions of The specific form of the function P, is not required
Sec. II, to obtain the neoclassical transport coefficients here. We need only observe that, since the collision
for a general toroidal geometry. operator is rotationally symmetric —
it does not involve
a preferred direction in velocity space and since the —
B. Classical perpendicular transport Maxwellians are similarly independent of direction, 6
We have already noted that, in a toroidal confinement can depend upon only the magnitude of the velocity.
system, classical transport typically presents only a Hence the first term of Eq. (4.4) cannot contribute to the
small correction to neoclassical transport. Hence the integral of Eq. (4. 3).
following discussion, which is mainly intended to illus- The rotationally symmetric property of the linea+ized
trate some remarks of previous sections, is brief, and collision operator implies that, for example,
only the major results are presented. More detailed
discussions of classical transport may be found in, for C'(gsin&) ~ sing,
example, the reviews of Kaufman (1966) and Braginskii for any function g which is independent of the gyrophase
(1965). Hence, after recalling Eq. (3.69), we see that the
perpendicular friction force comes entirely from f„and
Particle fluxes „a,
f, nd we have.
By our definition, Eq. (2. 77), the classical perpendicu-
lar particle flux is proportional to the perpendicular F~, = d'vm, v C,', (f, )+0(&'). (4. 7)
friction force
Note, by a similar argument, that the Paralle/ friction
F~~ = d v ~~v~Cg . force involves only : f,
Here E))e c/ v ~g UI[ C~~ y
+ 0 (4.8)
c. = g c., (f., f,) Equations (4. 7) and (4.8), when considered together with
Eqs. (2. 94) and (2.95), make explicit the previously
is the bilinear Fokker-Planck collision operator, and emphasized relation between classical transport, re-
v~ is given by Eq. (3.4). Any momentum-conserving sulting from gyro-oscillation (f ), and neoclassical
collision operator must satisfy transport, which results from g. c. motion ( ). f
We next calculate F~, explicitly, using the lowest-
d~v m, vc„(f„f, ) = 0, (4.2) order f, given by Eq. (3.73) and the small mass-ratio
collision operator given by Eq. (1.16). After substituting
since like-species collisions cannot change the bulk mo- the latter into Eq. (4. 3), we obtain
m In,. 8 Sf„2m, u; v Here we have noted that the collision operator must be
Galilean invariant. Since, as noted previously, the
right-hand side of Eq. (4.17) cannot contribute to the
(4. 9) friction force, we have shown that the electrostatic po-
tential cannot appear in the classical particle flux. This
The first term here, coming from C«(f„,f, ~), ay be m. conclusion is evidently independent of the specific form
simplified by integrating by parts, twice, and using the of the collision operator, provided C„. is Galilean in-
identity variant and rotationally symmetric. In particular, the
=-2v/v'. conclusion does not depend upon any small mass-ratio
s/sv u (4. 10)
expansion.
The integral of the second term, which comes from After substituting Eq. (4.14) into Eq. (2. VV), we obtain
C,';(f, ~, f, , ), is straightforward to evaluate, because of the classical perpendicular particle flux (Rosenbluth
the rotational symmetry of f, ~. Thus Eq. (4. 9) reduces and Kaufman, 1958)
to
nu„=-(m, Q,'v, ) '(V P —'n, V T, )—
, (4. 18)
v. ' v,'„, '+ (m. m. /7, )u, ,
F„=—(Sm'~'/4)m. J d'u/. , v, v =nu, ,
The first term of Eq. (4. 18) has been interpreted in Sec.
(4.19)
(4. 11)
II; a physical interpretation of the second, temperature
Here gradient term may be found in the review of Braginskii
(1965). It will be recalled that the ambipolarity condi-
Te vthe ~e (4. 12) tion of Eq. (4.19) follows from collisional momentum
]6&&/2 ~2&4+ lnA
i
conservation.
is the electron-ion momentum exchange time (or "elec-
tron collision time") and u, ~ is the first-order ion per-
pendicular velocity, given by Eq. (2. 61) (the 1-subscript 2. Energy floxes
is suppressed). From Eqs. (3.11) and (3.73) we have According to Eqs. (2.10) and (2.84), the classical
energy flux is proportional to the collision operator mo-
nx Vf. „, ment
where IQ, ~= eel/m, c (e&0), and the gradient to be taken d'v(mv'/2) v, C'(f, ) .
at constant energy
As in the case of diffusion, only the linearized operator
V f I = f, „[Vln P, + (v'/v, '„, —5/2)V lnT, —eV C/T] . a.ppears. But unlike the diffusion case, both like and
unlike species collisions contribute. In particular, ion
(4. 13)
thermal conduction is dominated by the effects of ion-
The integral in the first term of Eq. (4.11) is now read- ion collisions; to lowest order in (m, /m, . )'/', we have
ily evaluated, and after substitution from Eq. (2.61) we
f lnd G~= d v~ ~ v 2 vjC ~ ~
(4.20)
F =,
Q,
nx [VP —', n, VT, ],-
7-, I
,
I
(4.14) Since the like-particle collision operator. must vanish
when acting on a displaced Maxwellian [c.f. Eq. (4.17)],
since Q, = —Zm, Q, /m, , for ions with charge Ze. Recall it is clear that only the temperature gradient term in
that P = Q, P, is the total pressure. f„contributes to G~,. : to this order in the mass ratio,
An important feature of Eq. (4.14) is that the electro- the classical ion energy flux is driven solely by the ion
static potential does not appear. Since this fact has not temperature gr adient.
always been appreciated, we provide an alternate deri- The corresponding electron moment, on the other
vation. Suppose VP, = VP,. = V'T, =0. Then the electron hand,
distribution function may be written as
f, =2f, //(v~ ui /v, („, (4. 15) G~, = d'v nz, v' 2 v~ C,', „+C,',- (4.21)
Q„= —(I, 027, ) '[0.91p,V~T, —T, V~P) . (4. 23) (4. 28)
in Eq. (4.22), T, is the ion-ion collision time, defined
by
From Eq. (3.97), we see that f~ must be constant
along field lines, so that parallel gradients of n and T
T =3.m' 'T~'/47r' 'z e n ln A (4.24) enter the kinetic theory through f, . To estimate the
The numerical coefficient in the first term of Eq. (4.23) steepness of these gradients, we consider the parallel
pertains only to the case of unit ionic charge; the value force balance law, Eq. (2. 96). In the collision-domi-
of this coefficient for other values of Q may be found in nated case, the stress is nearly isotropic, so we have
Braginskii (1965). &
A
' &P' =-&IIP=&iI+ «&ii (4. 29)
Equations (4.18), (4. 19), (4.22), and (4.23) complete
the solution to the classical perpendicular transport = O(6~-'), (4. 30)
pr oblem. Notice that: because ~ v [recall Eqs. (2.81) and (2.82)]. It follows
E~~
Q, ,
—(m,. /m, )"' Q. . . (4. 25) f =f '+f '+' (4. 31)
while For a maximal ordering, the parallel electric field
term in Eq. (3.102) is also treated as O(66 '). More
Qec Pe +ec Ps +ic (4.26)
specifically, we assume
The ordering of Eq. (4.25} is not surprising, because we
would expect Q„/Q„- (p2/p2)(v, /v, . ), and Eqs. (2.18), (e/mc)v, ~6 'dA/st = O(5a '~U~}, (4. 32)
(4. 12), and (4.24) yield (p',. /p', )(v, /~, . ) —(m,. /m, )'~'. How- as a refinement of Eq. (3.98). [In small collision fre-
ever, the disparity between the ion energy flux and the quency regimes, Eq. (4. 32} is restrictive, and excludes
other cross-field fluxes yields two distinct time scales electron runaway effects. See Sec. V. ]
for collisional dissipation; if T, 4T, , some delicacy is Turning our attention now to Eq. (3.102), we change
required to obtain equations which simultaneously de- the energy variable, for convenience, from e to zo =e
scribe processes occurring on both time-scales. This —eC'/m. This transformation yields
matter is considered in detail in Sec. V.
(ii) to lowest order in the mass ratio, Q, , involves
~f/&t+ (&ii +~v~)'&f + ~/mEii&ii fs/sw =C(f f }
only &T, The reason for this has been discussed pre- (4. 33)
viously. Of course there is a correction term, Q, , ~
where the time and space derivatives are taken at con-
-Q„, which involves the other gradients; to our knowl- stant p. and m. The first-order version is evidently
edge, this correction has not been calculated.
It shall be seen that remarks similar to (i) and (ii) v~~n &f, -C'(f, ) =-v„'Vf„+eE~ f„/T, (4. 34)
also apply to the neoclassical fluxes. where C' is the linearized operator defined by Eq. (4.5).
Before expanding Eq. (4. 34) i.n 6, it is convenient to
C. Parallel transport recall, from Sec. I, some properties of the linearized,
unlike species, collision operators.
1. Kinetic equation
The basic kinetic equations describing classical paral- 2. Electron version
lel transport can be derived from the Fokker-Planck
The operator C,', is given, to low'est order in the mass
equation (Braginskii, 1965), even in the absence of a
ratio, by Eq. (1.16). When the distribution functions are
magnetic field (Cohen, Spitzer and Routly, 1950). For
independent of gyrophase, as in Eq. (4. 33), this expres-
reasons of consistency, our derivation proceeds from
sion reduces to
the drift-kinetic equation, Eq. (3.102).
All derivations are based upon an ordering in the (~)I&f ~ +»ii(~alii/~h )f u] (4. 33)
Chapman-Enskog parameter —
the mean free path di-
vided by the scale length —
which is denoted by
where 7, is given by Eq. (4.12)
where f~ is the usual Maxwellian, and f, is further & = 2 (&/& &)(I —8)(&/& 5) .
expanded (4. 39)
Equations (4. 35)-(4.39) are written for the case of a confined, quiescent plasma, we show in Sec. V that
single-ion species with charge ge; however, if all ion -u e
u.j e (4.49)
species have (to sufficient accuracy) the same parallel il ll
flow velocity u,. ii, then the equations remain correct for Therefore, for the purposes of this review, the parallel
any number of ion species, provided we make the re- flow terms may be omitted; we thus obtain from Eq.
placement (4.47)
n g ~ n ~
g ~ C~(a.') = f.uv Il(8+II/T. +II. Va. '),
which becomes, after substitution from Eq. (4.48),
jn Eq. (4.12). The quantity
g ff n~ g] n (4.40)
v'
C'(g, )o=f, v A, v — 5 A, +, (4.50)
the
—
Here the "forces" A. „, which drive the electron distribu-
is called the "effective g." tion away from equilibrium, are given by
Note that the electron-ion collision operator contains
a term, proportional to u, ll, which is independent of f„. A, =VII lnp, + ezII/T. , (4. 51)
The convenient way to treat this inhomogeneous term is
to introduce a function g'„defined by
2 =&ii »&. ~ (4.52)
2u,.2'ii v
fel = feN + 8'e (4.41) we have suppressed the subscripts of Eq. (4.48), with
Vthe
the understanding that A„= O(&A '). The solution to Eq.
Equations (4.35) and (4. 39) then imply (4. 50) is considered later.
' ; .
C. (f, ) = ~.; &g. . (4.42)
3. l on version
Evidently g, represents the electron distribution function
as measured in the rest frame of the ions [recall Eq. The collision operator describing the effect on the ion
(4. 16)]. Since distribution function of scattering by electrons (or by
relatively light ions) is given by Eq. (1.18). For small
.
C. (f.„(v u), f.„(v —u)) =o, — (4.43) mass ratio and f;, l- f„l, the dominant. term is
l l
to all orders in u/v, h„we have C,', (f) = (E„~V II/p, )f, . (4. 53)
Cee (2&II V II Vthe feII) = 0 . (4 44) where E, is the parallel component of the ion friction
li
Vfe therefore introduce the notation force. Evidently a transformation of the form of Eq.
(4.41) is not useful for the ion kinetic equation, which is
(4.45) therefore written in terms of f, , In other respects the
i.e. , C,', is the homogeneous electron collision operator, reduction of the ion kinetic equation for parallel trans-
corresponding to the case u,. =0. Equation (4.42) implies port proceeds analogously to the electron case. In par-
that ticular, the dominant perturbation to the ion distribution
,', has the form given by Eq. (4.48), so the
f,
c,'(f, .
) = c,'. (g;), (4.46)
function,
O(&b. ') version of Eq. (4. 34) becomes, in this case
so that Eq. (4. 34) becomes
VII n Vg, —C~(g) = eVII EII f, II/T,
A
— —ve Vf, II
c,', (f;, ) = c,' , (f, , ). -.
+ f IIVII VII lnpI +
v2
2
—
— Vll lnT. —eFII/TI
—f &V II
& V(2RIII V II/Vshe) (4.47) thi
With regard to the O(&A ') terms in Eq. (4.47), observe 9, =— Vii lnT; . (4. 56)
first that from Eq. (3.24) This relation describes ion parallel heat conduction.
v, . Vf„= O(&~0) + O(&'~-') . 4. Variational principle
Hence the g. c. drift is irrelevant in this order. The
%e have reduced the kinetic equations for both species
parallel flow terms would contribute if we assumed to the general form
u, il-u, , l-~D; in this way classical viscosity coeffi-
cients could be computed. However, for a toroidally c'(f) =q(v, (), (4. 57)
where C' is a linearized, homogeneous, Fokker —Planck (4. 55), omit the collisional energy exchange term, com-
collision operator, f
is the function to be determined, ing from lnf~. The result is (8S/Bt), =S.
and Q is a. given function of u and f, = v()/1). Because spe- That the entropy production be positive (Boltzmann's
cial cases of Eq. (4. 57) were first solved by Spitzer and "H-theorem" ) is guaranteed by the self-adjointness
collaborators (Cohen, Spitzer and Routly, 1950; Spitzer property
and Harm, 1953), it is often referred to as the "Spitzer
problem. " In the special case of Eq. (4. 50) with V))P,
(4.63)
=V)(T, =0, the solution is called the SPitzew function, and Equation (4.63) is readily verified from the definitions
denoted by v() E((f2 of the operators C,'. ; and C,'o.
Now consider the first order variation of the function
(4. 56)
Oe( I)f2 ) )lfeM/ e ~
f
[6l]. Letting g= +&f in Eq. (4.61), and using Eq. (4.62),
Spitzer's method of solution was based on numerical in- we find
6[61] =26[6 ] - 5[6]
tegration; the function f2, tabulated by Spitzer and
Harm (1953), is nearly exact.
Braginskii (1965) has solved Eq. (4.57) analytically, by where 6[6'] =6'[6f] and, according to Eq. (4.63)
expanding f
in a (truncated) series of Laguerre poly-
nomials in v (we shall see that the (-dependence of is f
readily determined, because of the rotational symmetry It immediately follows from Eq. (4. 57) that [6t] is varia-
of C'). This method is found to give quite accurate re- tion al,
sults, at least with regard to the relevant velocity mo-
(4.64)
ments of f.
A third method for solving the Spitzer problem, pre- Furthermore, by considering the second-order varia-
sented by Robinson and Bernstein (1962), is based on a, tion of [6l], one can show that [(R] ~(R, i.e. , the varia-
variational principle. We outline the variational pro- tional principle is rnaxsmal.
cedure here, because it resembles, in some respects, It is instructive to evaluate (R for the specific cases of
techniques used in neoclassical theo'ry, and thus pro- Eqs. (4. 50) and (4. 55). Using Eqs. (4.41) and (5.62), we
vides a useful introduction to the analysis of Sec. V. have, for the electron case,
[See also Kaufman (1960).] v2 5
Since Robinson and Bernstein allowed for arbitrary
relative magnitudes of the collision frequency and gyro- the
evaluation of the transport coefficients, we consider first TABLE I. A comparison of the parallel transport coefficients,
the ion Spitzer problem given by Eq. (4. 55). The solu- K~ as evaluated by Spitzer and Harm (1 953), by Braginskii
tion, f„
is represented by a trial function f~; for exam- (1965), and by the variational procedure outlined in the text.
Only the case of Z, ff =1 is included.
ple
f~ =a, v'(1+a, v')g, (4. 7s) Spit zer- Har m Br agin s k1 1 V ar iational
(4.as)
Finally, we note that the g dependence of f„, as given
in Eq. (4. 73), is exact. This can be seen by expanding
f and Q in series of Legendre polynomials, P„(g) (P~ =1, TABLE II. The Spitzer —Harm transport coefficients K~, for
several values of the effective ionic charge &,ff.
P, = g, etc. )
eff jeff jeff
1.98 2 32 2. 67 3.39
Rotational symmetry of the collision operator implies K(2 1.38 2.05 2.91 5.16
4. 17 6.75 10.01 20. 32
(4.79)
These equations permit a straightforward generalization radial fluxes of interest. The two ion species case, in
of Eqs. (2.93) and (2. 100) to the non-axisymmetric case. which one species is much more massive than the other,
Considering first the radial particle flux, we multiply the can be treated by a closely analogous argument (Ruther-
momentum conservation law, Eq. (2. 15), by ford, 1974).
—(c/e)g''&Qx&0 = (c/e)g q Br 4. 84 2. Neoclassical transport coefficients
where g is the metric determinant of Eq. (2. 43) and q is
In order to obtain transport equations from Eqs.
the safety factor of Eq. (2. 50). In view of Eq. (4.81), the
(4.89) —(4. 91), we must express the parallel gradients
stress tensor term is given to sufficient accuracy by (A.„A, and A,. ) in terms of radial gradients (dP/dg and
g'~2%(x VO ~
VP = g~ 2V) x %8 V( &P/9(
~
dT/dP) and the electric field. But our previous results
= &p/&r.
make this task very simple. Consider first the neoclas-
sical ion heat flux of Eq. (4.91). According to Eq. (4. 69),
The flux-surface average of this term vanishes because the Pm&/lel ion heat flux is
of Eq. (2.44), and, since,
g ~ VgxVO'uxB =u'Vg, to lowest order in ~ and &. Qn the other hand, it was
we obtain the lowest-order radial particle flux shown as Sec. II that this flux must have the general
fol m
(nu Vg) = —(c/e)q '(g'~'Br (F+enE))
~
(4.85)
in close analogy to Eq. (2. 93). The classical particle q~~,
= c/Zea)Ip; dT;/dq+I,
(5/2)(— (y)a.
flux, Therefore
(nu ~
ig)
(c/e)—
=q '(g'~'Br F, ), ~
(4. 86) A, =(5/2)(c/Zea)K, 'In, dT, /dq+. L, .(q) B/K, T, 4. 93)
(.
could be computed, for a given field geometry, from The function L, (g) can now be determined by observing
Eq. (4. 14). We restrict our attention here to the neo- that Eqs. (2. 53) and (4. 56) require
classical contribution, which is typically much larger.
Recalling Eq. (2.67), which may be written as (A,. a)=o
I($, O, g}=q g' B'B (4.87)
'
thus
we find that Eq. (4.85) yields A, = (5/2)(c/gea)K, 'n, dT, /dg[I —B'(I)/(B2)] . (4. 94)
(nu„c Vg) = —(c/e)(I(F~~/B+ enBr E/Br B)) . (4.88) Equations (4. 76) and (4.91) then provide the neoclassical
ion heat flux
Finally, we specialize to the electron version of Eq.
(4.88). Parallel momentum balance [Eq. (4.29)] then re- (q, Nc Vg) = —1.6(c/Z e)'(nz, . /7, . )P,. dT,. /dg
~
tron heat flux and dT/dg, can be traced back to our diffusion-time
scale ordering, Eq. (2. 20). Thus, for example, an ex-
(q, „c Vg) = (5/2)(c/e)P, T, (I A, /B),
~
(4. 90) cessive variation of T,. along B, would induce an ion
where A, =n VT, /T„as in Eq. (4.52); and, for the radial parallel heat flux which could not be consistent with Eq.
ion heat flux (2.68), i.e. , with the assumed equilibrium condition
~ q,. =O(5').
(q, „, &tt') =-(5/2)(c/e)P; T;(I A;/B), (4. 91) Electron transport equations are obtained from Eqs.
where A. , =n V T,. /T,
~ (4. 89) and (4. 90) in a. similar way. Thus we combine
One consequence of Eq. (4.85), of course, is that the Eqs. (2. 70) and (2.68) of Sec. II
lowest-order particle flux is automatically ambipolar Jii = c(I/B) dP/dg +K(—
g)B,
(i.e. , ambipolar for any electrostatic field), even in the
presence of asymmetry, for ~«1. This can also be ~~
= (5/2)(c/e) (P, /B) dT, /dP + L, (g)a,
seen from Eqs. (4.29) and (4.89), since the former im-
plies with the electron parallel transport equations Eqs.
(4.67) and (4.68). After the geometrical identities,
6 vp=R. . v p a =0 (4. 92) (aA, ) =0, (aA, ) =e(Z a)/T. ,
through O(&A '); thus the total pressure is more nearly are used to eliminate K(P) and L, (P), we obtain
uniform on each flux surface than are the pressures of K„A., +K,~A2 = —(c/e) dP/dg[I/B —B(I )/(B') ]
the individual species. It follows that, in the case of a
single ion species, Eqs. (4.89)—(4. 91) provide all the + (eK„/T. )a(Z a)/(a'),
Rev. Mod. Phys. , Vol. 48, No. 2, Part I, April t976
F. L. Hinton and R. D. Hazeltine: Theory of plasma transport
K»A, +K»A2=. (5— /2)(c/e)n, dT, /dg[I/B —B(I &/(B'&] as previously estimated. The original Pfirseh-Schliiter
result ("Pfirsch —Schliiter factor") was slightly different,
+ (eSC„/T, )B(Z B)/(B'& . (4. 97)
Dp2 /D, = 2q2, because the early studies did not properly
These equations can be solved for the A. „, in terms of allow for' Vi) T,
the radial gradients and (E„B&. When the results are The large aspect-ratio, axisymmetric version of Eq.
substituted into Eqs. (4.89) and (4.90), we obtain the (4.95) can be obtained in a similar way; this result was
r adial transport equations first derived by Shafranov (1965).
c &0& =-(c/e)'( ./&. )[&I'/B'& — &I &'/(B'&)
Finally, we note that the relation between the parallel
and radial g2adienls, expressed by Eqs. (4.96) and
~
[(~»/~) dP/dg —', (~12/-~)n, d T, /d g] (4. 97), can (in the axisymmetric case) be used to obtain
similar relations between the parallel and radial fluxes.
+cn, [(I &(E B&/(B') —(IB ' E/B B&], In particular, the parallel current can be written in
(4. 98) terms of the radial particle and heat fluxes (Hazeltine
and Hinton, 1973). That such a relation also obtains,
approximately, in other collision frequency regimes„
(q, c ig& = —(5/2)(c/e)'(m, T, /7, )
~
was first pointed out by Ware (1973).
[&I'/B') —&I &'/(B') ] V. KINETIC THEORY OF TRANSPORT IN
~
[(5/2)(z„/z)n, dT, /dg —(e12/z) dP/dg]. AXISYIVIMETR IC SYSTEMS
(4. 99) A. Small gyroradius expansion
Here The kinetic theory of transport in toroidal systems is
2
(4.100) based on a generalization of the well-known Chapman-
K K11K22- K12 ~
Enskog. method in. the kinetic theory of gases (Chapman
and the &'s are given, for several values of ionic charge, and Cowling, 1952), in which the motion along field lines
in Table II. Again, the axisymmetrie case is obtained is not necessarily localized by collisions. The basis of
by taking J to be constant on flux surfaces. this kinetic theory is an expansion in powers of the
Note also that Eqs. (4.67) and (4. 96) provide the fol- Poloidal gyxoxadius p~, which measures the degree of
lowing expression for the parallel current localization in the direction perpendicular to a flux sur-
0'2 B(E2 B &/(B & c dP/dg[1 /B —(I &B/(B ] &
face.
We begin by rewriting the drift-kinetic equation. , Eq.
(4.101) (3.102), in terms of the independent variable e =—e
The first term here gives the Spitzer current, as modi- + (e/m)((C & —4) (instead of e) and i1
fied by toroidicity and V~~P effects. Note than an elec- sf — e B(4& Bf
trostatic field cannot contribute, because of Eq. (2. 53). Btl
+(vii+vu) If + . Vl Btl
+v2 E2 — =C(f
Bc f) ~
ence of the diamagnetic current. The return current cannot order v, h(B~/B)/l, while the third term v2 Vf, is of
vanish (for nonuniform B), and will generally change order v, „p/l', and smaller than the second by a factor
sign across the plasma cross section. p~/l, where p~ =—pB/B~ is the poloidal gyroradius, men-
Equation (4.98) can be seen to agree with the estimate tioned in Sec. II. We choose the basic small parameter
of Sec. III, to be 5~,. —
= p~,. /l, the ion poloidal gyroradius divided by a
characteristic length (which may be the major or minor
use radius of the plasma). In this section, we use the as-
for large aspect-ratio, axisymmetric geometry. This sumption that B~/B«1, so that p&&~ p, and we may ex-
case is considered in Sec. VI; it is convenient to antici- pand the drift kinetic equation to any or der in p~.
pate some elements of the Sec. VI discussion here. %'e make a noncommittal choice of ordering for the
With B =Br =I/R, the flux-surface averages in Eq. like-particle collision terms, that these are of the same
(4.98) reduce to order as the term vt~ 7f. That is, the collision fre-
quency is of the same order as the bounce or transit fre-
(I 2/B2 (I )2/(B2 (Il2 (B-2 &-1 2~2
& & &
quency. Subsidiary expansions will be carried out later
where & is the minor radius of the flux surface, which for the limiting cases vl/v, „(B2/B) small (Sec. V.D) or
is assumed to be nearly circular. Using ~VJ~=B~Il, we large (Sec. V.E) compared with unity. Before expanding
find that the flux u~c. & is related to the density gradient f in powers of 6~„we must also choose the ordering of
dn, /dr by the Pfirsch-Schliiter diffusion coefficient E~], and specify the collision terms in more detail. We
make use of the fact that (m, /m, . )'2 is small, in order to
Dps =2(m, 02', ) '(~„/a)q2(T, +T, ), (4.102) simplify the unlike particle collision terms.
where the safety factor q is roughly given by
q = (r/Il)(Br/B2), for this geometry. Referring to Eq. 1. E lectron-ion col lision term
(4.18), we see that %'e recall that the linear ized electron. -ion. collision
DP2/D, = 2(@22/e)q2, operator is given by Eq. (1.16), to lowest order in
B'V d v Olog
of the drift perpendicular to flux surfaces, the deriva-
tive of v~~ may be taken at constant m = v'/2 =e
—(q/m)(4'), i.e., we may use v, =[2(m —pB)P~' or
d'v log f, C„(f„,f, ) + d'v log f„C„(f,o, f, ,) . - v(1 —&&) ', where ~=g/w. We also use A, and v as
v~~—
~
U(w)
Bw
-.
we have already written down the (somewhat more ac-
~,
stant along field lines, at constant p. , e. Hence, curate) moments of the exact Fokker-Planck equation,
Eqs. (2.13) and (2.14), which need only to be averaged
f,o =N, (m, /2wT, )~' exp( m, e/T, ), (5.17) - over a flux surface to yield the desired equations. This
f, O=N, (m, /2mT, )+2 e.xp. (. m, e/T, ), (5.18) topic will be pursued in more detail in Sec. VII.
We now give a physical interpretation of some of the
where N, =n, oexp[ e(& )/T—, ], N~ ——n,. expO[Z, . e(4')/T, ], above equations. If the time derivatives were replaced
with n,-~, n;o, T„and T, functions of g and t. Charge in Eqs. (5.11) and (5.12), their solutions would approach
neutrality requires that m, o =g,. n,.o.
The fir st- or der equations are
v~~ Vf., —(C.', +C.',. )f., = v, . Vf., —(e/T. -)v~~&, A„'
(5.19)
the Maxwellian steady states, given by Eqs. (5.17) and
(5.18), in a time of the order of the 90'-scattering colli-
sion time. The decay of small departures from these
steady states, is described by linear equations. In the
short mean-free-path limit, this decay is related to the
classical parallel transport processes, considered in
v~~ Vf —Cia fgi =— '
v~~ Vf~o ~ (5.20)
See. IV. Thus heat conduction makes the temperatures
where the superscript l denotes the linearized collision uniform on a magnetic surface, resistivity reduces the
operators. Qn the right-hand sides of these equations, relative parallel flow velocity to zero, and parallel vis-
the following expression for the drift velocity is to be cosity eliminates the parallel mass motion through re-
used [recall Eq. (3.24) or Eq. (3.32)] gions of spatially varying magnetic field strength.
vq=-v~~nxV(v~~/n) vq, A similar interpretation of the time-dependent version
of Eq. (5.20) can be given if the ion temperature gradient
= Iv n V(v)(/0),
I) is zero. In this sPecial case, the steady-state solution
where the spatial derivative acting on v~] is taken at con- ls
stant e, p. , neglecting the 0(&) term [contained in Eqs. = 2v zt
the
fg. , (5.25)
flux surface. Then, since we only need the component where
).
vii Vf; —C'; » f; = -vq; Vf. i + v (n; Z, eEii + I'ii,. ) o/P,
'
where 0 is independent of 8 (and f); This means that 6 v~ Vp = (BzJB~,) v „n V(v), /Q), (5.36)
may be obtained from a knowledge of the poloidal field
where we have used RB~ RpB~ Further, we write
'
'&
dE, (5.31) (5.39)
0
so that
where the integration path is tangent to the poloidal pro-
jection of a field line (on a plane of constant g). The re- Vp = (v), /Q~) n V(v, )k), (5.40 )
quirement that 0 increase by 2m on one poloidal circuit
where
around the magnetic axis implies that
Q, =— (eBo/mc)(Bp, /BT, ) = ego/mc (5.41)
is a typical value for the gyrofrequency in the poloidal
magnetic field. We shall usually use the approximation
Note that 0 is also related to the derivative with respect
B =Br, so that it =R/Ro.
to poloidal flux of the volume enclosed by a flux surface
dV d9 dip 2n
(5. 32) 2. Transformation of the kinetic equations
2w dg B V0 ~
and to the metric determinant of Eq. (2.43): k =g Rather than work with Eqs. (5.21) and (5.29), it is use-
The natural choice of & as the symmetry angle, in ful to transform them by introducing the functions II„II„
axisymmetric systems, requires that B VP =Br/R defined as follows (now using —e for the electron charge
which is not constant. Therefore, the local pitch of the and Z,. e for the ion charge}
field lines, B Vg/B VO, is not constant, i.e. , the field
~
(5. 33) v, 0
)
s lnp; Z,. e s(4&
Q~p Bp Tg sp
With the above choice of poloidal angle 8, this is simply
'~ dl'
(5. 34) +(Z, e/T, ) f o
0
.
P
[BE~ —B2(EgB&/(B~&]+Hi 1
8 9
v = (vv;/20 v)IV;//I— lnP, + (T;/Z;T, ) ln p, , (5.56)
„~(
Bp Bp
A, = Bp
lnT, , A„=B,(E B&/(B'&, (5.5V)
which follows from Eq. (5.43) and the definition of u, )),
with A„= —(T, /Z, T, )p, , (5.58)
d'v v)) H, =-
,
(n, u('„,./2n, .~)u,. /h. (5.47) A. „=—(n, Z, e. / m.
Q,. )(E B&/B, A„=
9
Bp
lnT, .
The fact that this parallel ion flow is proportional to B (5.59)
(=Bo/l() follows from the fact that it must be divergence- The driving term A. „. has been expressed as an inte-
less, since the other terms in Eq. (5.46) provide the gral of H, in Eq. (5.50); this equation can also be written
parallel return flows discussed in Sec. II. Note that the as
—
radial electric field B(4 &/sp has cancelled and does not
appear in Eq. (5.45).
The equation satisfied by the function II; is
A„= ((/H, )
f d'v H, q, (5.60)
(5.61)
The term (Ed, B&, which contains the friction force, can From the linearity of Eqs. (5.51) and (5.52), it follows
be expressed as an integral of Iz„as follows. By multi- that the right-hand sides of Eqs. (5.60) and (5.61) may be
plying Eq. (5.21) by m, v)) and integrating over all veloc- written as linear combinations of the six forces A. „„A.„,
ities, we obtain the parallel electron momentum equa- Then A„. and A. 4, may be eliminated, by expressing them
tion: in terms of the four independent forces A. A. A. 3„
and A. „.. Before doing this, we introduce the transport
„, „,
d'v m, v))v)) Vf„Fj)H —n, () eE)) . — coefficients.
= o B (Zii B)/(B'),
Jii, — These inner products are the transport coefficients, re-
ii
4 gives
(1/Te}(( Ii +lie)/~) =g(net
n=J,
gne)Ane ~ (5.v2)
~;/T; = (e., g„)A.. . (5.82)
while Eq. (5.76) yields
e;/T; =Q(P. , g.i)A. i . (5.73)
n=l = (T;/Z;T, ) p. ; = —(T-;/Z, T, )(P„g„)A„,
A„— (5.s3)
Here, following Bernstein (1974), we have introduced where A„. =s(lnT, . )/sp. These are the lowest-order re-
inner products, defined, for example, by
sults, in which the ion fluxes are unaffected by the elec-
(~„g. ) = ( . f ' a, d. .Ug (5.74)
tron forces, and the ion force A„. appears in the expres-
sions for the electron fluxes in a simple additive way,
through A. 4, . 8 8
A first correction to Eq. (5.82) may be obtained by A, = lnT, , A2~ = 1nT,
substituting Eq. (5.83) into Eq. (5.75); we obtain
Neglecting the effect of collisions with the electrons, the
(5.84) species "b" particle flux is given by the ambipolar con-
dition
with
3 Z. r. +~I =0-
A„= (1/Z;) Q(n„g„,)A„, The relations between the fluxes and forces are very
n=l.
similar to Eq. (5. 70)-(5.73)
—(»/&;»)(~. , z. .)(p„g,;)&.;I (5.s5)
l. =Q(n„g„.)A„. ,
Ily comparison with Eq. (5.77), we see that the weak-
coupling approximation is equivalent to assuming I' = I,
or, more precisely, V./&. =Q(n.
n=&
, g. )A. .. .
I' —1 =- (T, /Z2 T, )(n„g„)(p„g„.) «1. (5.88)
eb/7'b =Q
n=g
(P„g.b) A.b,
More specific conditions will be given in Sec. VI, for the
various collisionality regimes, after the transport coef- where the auxiliary forces A„and A», defined by
ficients have been obtained for the large aspect ratio Ass (+a b/Zb~ )I b Al b (+II +)/ b bb 0
case. We will find that the weak-coupling approximation
is satisfied, in the electron-ion problem, under most appear because of the frictional couP/ing of the two spe-
conditions. cies. [Note that, neglecting the effect of E)), there is no
force corresponding to A„, in Eq. (5.57); the force cor-
responding to A4, has therefore been labeled A„.] These
5. Two ion-species problem auxiliary forces are eliminated just as in the electron-
As a first step in the solution of the problem of trans- ion problem, to give
r. ~) —
port in the presence of impurities, we consider the z'
special case of one-imPuxity sPecies. We consider a re-
stricted problem, in which the two ion species have
A, .=-(Z.r, /Z, (p„g„)g(n„g„.)A„.
2
S-=—
e
where H, = H, /.f, — .
,
Steady-state ion heat conduction is governed by the
law that [S,. ] is a minimum (where the square brackets be stationary with respect to variation. s in II, , where
indicate the functional of H;), subject to the constraints ~, and ~ are Lagrange multipliers. These variational
imposed by the guiding center drifts and the existence of principles give Euler equations for II,' which are con-
a nonzero temperature gradient. The subsystem con- sistent with the constraints, provided that ~, =2, & =-2.
sisting of the plasma between two nearby magnetic sur- By noticing that only [&a) contains both H, and H, , we
faces is maintained away from thermal equilibrium, may combine the two variational principles. The func-
since the drifts and the temperature gradient prevent tional
the distribution function from becoming a Maxwellian. (5.93)
If the subsystem were a closed system, its entropy
would increase to the maximum value consistent with the is thus station. ary with respect to independent variations
constraints of constant total number of particles and of II, and II; . In fact, by using the negative-definite
energy, and the distribution function would then be Max- property of the collision term
wellian. Aetua1. ly, the entropy increases indefinitely J «Q (5.94)
because a steady state difference is maintain, ed, between
the distribution function and the equilibrium Maxwellian. . (with equality if and only if H,'=—0), i. t
is easily shown
This difference must adjust so that the rate of increase that [S,] is maximal with respect to H, and minimal with
of entropy is as small as possible, however. respect to H, Its extremal value is given by Eq. (5.92),
The mathematical form of the const~ainfs, imposed by or
the temperature gradient and the guiding center motion,
may be expressed in the form of integrals obtained from S, = &D = -(q, /T, } (5.95)
Eq. (5.52). We first define the even and odd [in o
=sign(u~, )] parts of the distribution function using Eq. (5.91) and the definition of the heat flux, Eq.
'
H, =--.'[H,. (~ =+) +H,. (v = -)),
(5.84). From Eq. (5.94) we have S, &0, so that a positive
coefficient z, (the ion thermal conductivity}, defined by
(5.88)
H, = ,'[II, (o =+) —H, ((x=-—)]—. q, =-~,. (aT, /ap),
Then Eq. (5.52) yields the following equations for the odd may be obtained from the extremal value of the func-
and even parts of H tional [S,.]:
I" gl~ 7H; -Cg;H, =-vg; ' tp(~~/~~g —5/2)
+.
'f;„ S,. =~,. (s lnT, /ap)'.
This variational principle was first introduced by
(5.89)
re-.
Rutherford (1970) in a form which is somewhat differ-
Iv~, lie &H, —C„H,
+
=0. (5.90) ent, but equivalent to Eq. (5.93}.
+ A somewhat more convenient form of the variational
By multiplying Eqs. (5.89) and (5.90) by H, and H, principle is obtained when g„ is used in. place of II,. ; the
spectively, and averaging over a magnetic surface, we extremal value is then equal to the inner product
find the constraints to be
(P. , u;;) = ~;. -
J~ —J+ —JD =Q, er~+ J =Q, (5.91) All of the inner products required to obtain all of the
fluxes of both particle species may be calculated from
where the following variational expressions:
where c stands for either e or i, y„stands for either collision term in Eq. (5.80) or (5.81) may be treated
o, „or P„, and the superscripts a indicate the parity in o. formally as a small perturbation. For systems with ar-
[Note that the n„'s are all even, and P, is even, but P, bitrary aspect ratio (or an arbitrary fraction of trapped
is odd, so that o.„=ti2 =P, =0.] The errors in these particles) the effective collision frequency for trapped
variational expressions are quadratic in the errors in electrons is of the order of the momentum transfer fre-
the functions g„',. The diagonal (m =n) terms in the ma- quency v„. (v), defined by Eq. (5.4). The bounce fre-
trix [U '„] provide mini-max variational principles which quency is given roughly by
are useful in solving for the functions g„„g„,. The ex-
trema1. values of all of the variational expressions give ,
~, —(2m/I. ) ~ v,
the required inner products where I., (=2nqB) is
the connection length. The "bana-
"
na, or "collisionless" regime for the electrons is then
V&'„'=- d'v g„'. y -g„-, y-. (5. 97) defined by
(P„a„)= (P„a„). - (5, 99) The lowest-order equation, obtained by neglecting the
It should be evident that the diagonal matrix element collision term is then
corresponding to the even force P, is negative, as in the
electron case
"
v[, n Vg„(, = -o. „f„, (5.102)
where we have treated all of the n„'s uniformly, as being
of order unity in the collision frequency parameter. The
It is easily shown that the other diagonal matrix element special property of the ~„'s, for the axisymmetric case
is also negative considered here, is that they have the form v„n . V(. ).
We shall consider here only the n =1 equation in detail;
(p„g„.) = v t[' = (f
«'« [ g,',. c,
g„+ g, ; c z, ;() ~ « . g
the general solution of Eq. (5.102) is [recall Eq. (3.104)]
The fact, that all of the diagonal matrix elements are (5.103)
negative, corresponds to the definition of the forces as where
being positive in the direction opposite to the direction
of the fluxes which they drive [e.g. , diffusion due to a (5.104)
density gradient goes in the direction opposite to the The function 6„ is determined by the first-order equa-
density gradient]. This is due to the positive-definite tion
nature of the entropy production.
(5.105)
D. Banana regime
"
The unknown function gz', can be eliminated by the oper-
ation g d8%/v„, where the integration is a round trip be-
Definitions tween turning points (at which v„changes sign) for the
Vixen the effective collision frequencies of the trapped trapped particles (h;„~ XBO~ h,
„), or a complete cir-
particles are smaller than their bounce frequencies, the cuit around the minor circumference, for the untrapped
Untropped
" d68 v„h I
(5.106) I
VII ep
o, () } o,'() }
The perturbation analysis thus leads to a reduction of —Bounce Points
Eq. (5.80) (for n= 1) to this constraint equation, for a
function G„which is independent of 0. This analysis FIG. 5. Phase space, showing trapped and untrapped regions,
the boundary layer, and the trapped-particle bounce points.
was first carried out by Galeev and Sagdeev (1968).
The values 1+ e and l —e for hnl» and Am~ apply only to the
Very similar constraint equations may be obtained for large aspect-ratio, circular cross section case, where
the functions g„, g3„and g4„by following this same e = r/R, .
procedure. These constraint equations have not been
solved, in general, because of the complexity of the
Fokker-Planck collision terms. %'e will subsequently
82(x) dg~
give an explicit solution for a special case which is
tractable: the Lorentz model. v "~(x)
Io"
Aep
i~"' ")=' (5.111)
We now show that Eq. (5.106) is equivalent to a, vari- h;„~ LB~~ „(trapped). Since these equations are
ational principle (Rosenbluth, Ha. zeltine and Hinton,
for h,
equivalent to Eq. (5.106), the variational property of
1972) which is useful in obtaining approximate solutions.
We first note that Eq. (5.104) and the boundary condi-
"]
[Q,', has been demonstrated. That [Q,', has a. mini- "]
mum follows from the negative-definite property of the
tions imply that &„ is even in o, in the trapped region collision operator.
of phase space, since the odd part of g,'o' (and hence of A general property of G„ follows from Eqs. (5.111)
G„) must vanish at the trapped pa, rticle bounce points, and (5.110) (when summed over a) and the fact that the
,
& = II, (X), and also be independent of II. Thus, collision term preserves the parity in v [e.g. , C, (v„f,o)
G„(o' =+ ) —G„((r = —) = 0, f or I1,„&XB, & I1 „. 1S odd 111 (r]:
(5.10V) G„(o =+ ) + G„((r= —) = 0 for all X,
3. Variational principle i.e.,the function t"„ is odd in 0. It follows from Eq.
(5.107) that G„ is identically zero in the trapped parti-
The constraint equation, Eq. (5.106), is equivalent to cle region.
the requirement that the functional The variational expression for [Q,", '] may also be ob-
tained directly from Eq. (5.96). We note that g,', is first
(e) G~
2 77
) 3 v))k
IQ ep I
~
" ' v))k
IQ ep I " ") order in v, ~qR/v, and that the trial functions must satis-
fy the constraints of Eqs. (5.103) and (5.104), which im-
(5.108) ply the cancellation of g,'Jv„ln Vg„", and g,', a."„ to first
be a minimum, subject to the constraints of Eqs. (5.104) order. Therefore, sj.nce o.', =0, we have
and (5.107). To prove this, we note that variation of G„ [Q(e)] [y(e)]
and use of the sen-adjointness property of C„Eq.
neglecting terms quadratic in v„qR/v. Furthermore,
(5.67), gives from Eq. (5.9V), it follows that the extremal value is the
"d88 ~ VPg required transport coefficient, in this case the diffusion
6[Q,", '] = —
-
v'dv
fy & v II IA ep I ")' coefficient
(5.109) Q„- -(o(, g(.). (5.112)
where we have used The other transport coefficients may be obtained from
variational principles also:
~h
., 2lv ~'
~dy
d v=2+ I
vdv
vp
[Q.".'] =— J (5.113)
l
tionals are the transport coefficients: term. The first term on the right-hand side, of Eq.
(5.115) is the negative of the Spitzer distribution func-
tion [except for factors of E„and h, which are accounted
4. Lorentz model for in the definition of g3„Eq. (5.69).j Comparing with
We now consider the special case Z; 1, in which » Eq. (5.42), we may thus identify G~, with f„/E„, when
electron-electron all other driving terms except E„are absent. [For pur-
collisions may be neglected. A sim-
ple solution of the constraint equation, Eq. (5.106),
poses of interpretation, we ignore the distinction be-
tween E„and B(E„B)/(&').j The time-averaged effect of
may be obtained, in this case. The results are, pre-
sumably, qualitatively similar to those of the more gen-
collisions, on the difference between the distribution
function f„and the Spitzer function (proportional to E„),
eral case Z,. -1, although the latter have not been cal-
must be zero. For the trapped elect ons, the time-av-
culated, except for large aspect ratio. An analysis of
the I.orentz model for general geometry was first given eraged effect of &, is zero, since the work done by
~
where X' =—XBp The second expression for l, y was ob- 3 c 3 nepe~
tained by integrating by parts in A. ', and is a convenient ( 3 tg4e) 4 0~PO I33t (c'4, g4. ) = -8
Te
form for calculations with specific geometries; for the
case of large aspect ratio and circular cross section, (5.127)
discussed in Sec. VI, we have
I, = 1.38(2E) + O(C ) ( ., g..)=-—n,' e'&,Z'I. .
8
'ti' Bl e e
(5.128)
where e=—3'/R0 is the inverse aspect ratio. Although the
assumption of a single minimum in B is implicit in Eq. where I» is given by Eq. (5.123), and
(5.123), it is easily generalized to include multiple x'dh' ~ &min
integral expression given by Eq. (5. 120), the factor v„. " (v'/v'„„;
cancels out of the integrand, leaving an energy-depen- -(P., g. ;) = d'v —5/2) C, , g,'P . (5.135)
ip
dent term proportional to v (—,m, v' ——,T,). This inte-
grates to zero: the average electron energy carried in- The banana-regime variational principle is
ward in the trapped particle pinch effect is exactly —,'T, .
When the electron-electron collisions are included, (5.13 6)
there is no longer an exact cancellation.
Because of the Onsager relations, we have a term the minimum value of which is the thermal conductivity
proportional to the pressure gradients, in the parallel (5.137)
current density, whose coefficient is (n„g„)= (n„g„).
Such a. term was first pointed out by Galeev (1971), and In the case e = 1, the perturbation analysis is some-
independently by Bickerton, Connor and Taylor (1971). what different. The leading term is of order (v;;qR/v) ',
The corresponding current is called the "bootstrap cur- and satisfies
rent, " since it makes possible the idea of a "bootstrap
tokamak, "
operating in a steady state with E~= 0. This
possibility was first investigated by Blckerton et Ql, It is determined by the constraint equation
(1971), and independently by Kadomtsev and Shafranov
(1972). A limit, on plasma pressure in a tokamak, was
dOB
C iigii ', &
= 2QieBofio
2 ~&
&& B2dg
II i tbi {B2) ~t
tion. The constraint equation is found to be Ceg(ie =v(( ' &(gne +v(('4feo)& (5.140)
tl
(o)
Czi&2i =0r (5.13S) "
where Z„ is defined by Eq. (5.114). Equation (5.138) im-
plies that g„', is the perturbation of a Maxwellian due
to pressure and temperature perturbations
where -"
g„', = [a„+ (v'/v'„„, —5/2)f „]f . (5.141)
g2'P = —(v„h/A~e)(v2/v2„„; —5/2) f,o+ Go, , (5.134) Equation (5.139) is then the Spitzer problem [c.f., Sec.
with BG2;/88=0. The ion thermal conductivity is given, IV.C] and the parallel gradients of pressure and tem-
in terms of the solution of this equation, by perature drive parallel current and heat flows deter-
VS gf7 g
/ 1 +g g V j Off + Kg 2 V }5 ff
[ f
(5.143)
where the K „'s are the parallel transport coefficients, FIG. 6. Toroidal coordinate system, showing its relation to
defined by Eq. (4.77) and Table II. Equation (5.140) has the cylindrical coordinates R, Z, g.
solubility conditions which determine the function g„', ""
uniquely
After solving Eqs. (5.142) and (5.143) for V„a„and V„b„,
(5.144) with the left-hand sides of those equations given by Eqs.
(5.146) —(5.149), and substituting into Eq. (5.150), we ob-
tain
(5.145)
The solution of these differential equations may be writ-
.
(c'., g. ) = —— "~.(p'./, /7. ) [&k'& —&k '& '], (5.151)
ten as
r
(5.152)
d'vv„g„", ' = — d2vv', jy„ f„+k, (p)/k,
., g, .) —" . .'/,
II
VI. TRANSPORT GQE F F IGI ENTS FQR LARGE so that the coordinate 8, introduced in Sec. V. B, is rela-
ASP EGT- RATIO TQKAlVIAKS ted to ~ by
A. Geomet~ 8 = (d —A(Y/R, ) Sin% + O(6'). (6. 10)
The problem of specifying the coordinate system, in a We will continue to work with the coordinate 8; in fact,
concrete way, is greatly simplified in the large aspect- the -approximation
ratio limit
(6. 11)
q = x/Ro « i. (6. 1)
will be adequate.
Here y is the minor radius, and R, the major radius; The toroidal flux is given by
these will be defined subsequently. In this limit, many
m'+ O(e')
solutions of the flux surface equilibrium equations are (6. 12)
known, the most generally useful being given as follows so that the minor radius coordinate p, introduced in Sec.
(Shafranov, 1966). V. B, may be identified with x,
We define the toroidal coordinates y, (d, f in terms of
the usual cylindrical coordinates R, Z, f (see Fig. 6) by B. Electron transport in the banana regime
cosa, Z =+ sin~,
R= R, (x)+x 1. Large aspect-ratio expansion
where the function R, (x) gives the distance from the sym- In the large aspect-ratio limit, the solution of the ki-
metry axis to the center of the (circular) flux surface netic equations is simplified by making use of the small-
cross section whose minor radius is x. The distance to ness of the fraction of trapped particles, which is ap-
the magnetic axis is Ro(0) =Ro. The poloidal flux function proximately (2e)'~'. Assuming that
is then (2&)'i' «1, (6. 13)
r
dyBp. (~)+ O(e'), the trapped particle region of phase space is very nar-
y = R, (6.3)
0 row, consisting of the small range of pitch angles such
where Bpo(x) is the effective poloidal field magnitude, (the that 1 —&» XBo» 1+ &, where A. is the pitch-angle vari-
time dependence is not written explicitly here) which able, defined by Eq. (3.63). A variational method (Rosen-
must be related to the pressure gradient by bluth, Hazeltine and Hinton, 1972) can then be used to
calculate the transport coefficients exactly, in the limit
1 d / z 1
(6.4) of small e (and small collision frequency).
t- dt Bo d'v d'v We begin the analysis with the banana regime varia-
Here we have written the toroidal field as tional principle of Sec. V. D. We wish to minimize the
functional
B, = 1,(~)/R. (6. 5)
d
Equation (6.4) is the equation of pressure balance in a, [q„]-
(e) =—
) d 3vg„C, g„, (6. 14)
cylindrical plasma of circular cross section, and is a
consequence of the equilibrium equation, Eq. (2. 90), to subject to the constraint that g„has the form given by
zeroth order in &. A term of order & is not present in Eq. (5.103), where G, is independent of 8 and is an even
function of o (= —signv„) in the trapped particle region of
Eq. (6. 3) provided that Ro(x) satisfies
phase space. [Note that we have dropped the superscript
~Bp'0 [R', (y)/2]+ ~BE= 8mB (6. 6) (0); the functions g„and G, are here understood to be
zeroth order in collision frequency. ] The minimum val-
That is, the relative shift, of the flux surface centers is ue of [Q,", ] then gives the desired transport coefficient,
determined by the profiles of pressure and poloidal field. according to Eq. (5. 112).
In general, dRO/Ch= O(c), provided that P&— = 8nP/B, tbe '„ The function G„defined by
ratio of plasma pressure to poloidal magnetic pressure, ie = (vn@/'I +eg I )feo+ Gi (6. i5)
is not much larger than unity (it is typically somewhat
less than one, in present tokamak experiments). is assumed to have an expansion of the form
Neglecting terms of order &', the poloidal field is 1
(0)+ g(&)+. .. (6.i6)
therefore given by
where the superscript gives the order, in q',
of the
RPv~(&) derivative, XBG/81, in tbe trapped particle region. Be-
" R[1+ (dR, /d~) cosa)] ' cause of our assumption concerning the smallness of
this region in the A. direction (see Fig. 5) the predomi-
nant effect of the collision operator C, is contained in
where 0„, P„are unit vectors in the x and ~ directions, the X derivatives, i. e. , pitch-angle scattering is more
respectively. The magnitude of the poloidal field is important than the other collisional processes contained
~=Bp.(Br)[1 — (A~)(~ /R) cos(u] '+ O(e'), (6.8) in C, . To make use of this assumption, we write
G"' = -(v„m/I fl„l )f., + I "'. (6. 17)
where
The function 0' ' is assumed to be "localized" in X, in
A~/R, = -(~/ +RdR, /dr ), — (6.9) tbe sense that X&h'o'/&X is of order unity in the trapped
region [much larger than vBIr(0)/Bv, which is O(q't')], I = 1.38(2~)'t ' (6. 25)
but is much smaller, of order &, in the untrapped re-
which is the large aspect-ratio limit of I», defined by
gion. The consistency of this assumption can be verified
Eq. (5. 123).
after the solution is obtained. By making use of Eq.
(1.19) for C, and expressing the velocity integral in Eq. By comparison with Eq. (5. 121), the effect of elec-
(6. 14) as a symmetrical quadratic form, we find that tron —electron collisions is seen to be an inc~ease, in
the value of the diffusion coefficient, by a numerical fac-
[Q "] — "'
=
dg
2n
d'vh' 'C")Ir(')+
e O(&) (6. 18) tor of 1.53 for S;=1, 1.26 for Z,. =2, etc. This may be
surprising, in view of the fact that neoclassical diffu-
where the first term on the right-hand side is of order sion, like classical diffusion, is due to unlike-species
&' ', and is the dominant term, containing the pitch- collisions. In the classical case, the form of the dis-
angle scattering part of C, : tribution function, which enters the integral for the (per-
pendicular) friction force, is not affected by like-species
C,"'=—(v„(v)+ v„(v))S. (6.19) collisions. In the neoclassical case, on the other hand,
The electron —ion term v„.Z was defined by Eqs. (5. 3) the distribution function which goes into the calculation
and (5. 5); the elect+on electr-on collision frequency for of the (parallel) friction force is quantitatively affected
pitch-angle scattering is defined by by electron —electron collisions, even though it is of ze-
roth order in the collision frequency.
(6. 20) In order to demonstrate this, we return to Eq. (5. 119),
and try to discover the distribution function g„which
where must be used in this expression. [Note that Eq. (6.23)
contains the pitch-angle scattering operator C,' ', rather
((x)-=(1 —, 1
2x' ex((x)(.
X
(6. 21) than the exact C„which appears in Eq. (5. 119).] A bet-
ter approximation to the function G is needed, so we re-
and v~, =(2T,/m, )'t'. The O(e) correction in Eq. (6. 18) turn to the exact constraint equation, Eq. (5. 106), and
contains
sidered here.
6"'+ ~ ~ ~
C;
as well as C, — (0) it will not be con- subtract from it the zeroth-order constraint equation,
Eq. (6. 22). To first order in the expansion of Eq. (6. 16),
By using Eq. (6.17) and the properties of G"' (that it we obtain
is independent ()f 9 and even in o in the trapped region), dOB
we may verify that minimizing the zeroth order part of (C G(1) C(1)I (0)) 0
[Q,' ], Eq. (6. 18) yields the zexoth olde~ const~ai-nt
equation, where C,"'—
ution of
= C, —C, "'. We define a function f, as the sol-
dgB C(0)
In„l
II
" + g(P) O (6. 22) f C (1 )Ir(0 ) (6.26)
By using this equation, Eq. (6.17) and the above-men- Then the first-order constraint equation becomes
tioned properties of G"', the minimum value of [Q'„"] B (y)
may be written as II
0.69(26} z
lg
O(E, }
XBo
I+K,
FIG 7. Banana regime distribution functions: (a} the "localized" function O'Oi [divided by vf o/~g„~] vs p„~o, for s=o. i,
f~~, defined by Eq. (6.26) [divided by the same factor]. Both are evaluated at 0=0, with the O(c) corrections
(b) the function
ignored. Also shown is (c) the smoothing of the function g&~ (=4 ) in the boundary layer, as described by the solution of Eqs.
(6.9O) and (6.91).
effects. In this connection, Kovrizhnykh (1970) has poin- C. IOA tl'BASPOI't IA th8 baAIA8 f89IAl8
ted out that the pinch effect coefficient for a weakly ion- 't. Ion thermal conductivity
ized plasma is smaller by roughly a factor & than in a
fully ionized plasma. This is due to the much stronger When the ions are in the banana regime, the calcula-
effect of velocity diffusion, in the latter. tion of the ion thermal conductivity —(p„g„) proceeds
The effective reduction in conductivity, given in Eq. in a way very similar to the calculations in the previous
(6. 33), is also enhanced by electron —electron collisions. subsection. The main difference is that the operator C;;
This is because the current-carrying untrapped electrons conserves momentum, while C, does not. The conse-
see the trapped electrons (in a time-average sense) as quences of this fact will become clear subsequently.
fixed scatterers, which are present in addition to the ion The banana regime var iational expression is given by
scatterers. The latter point has been discussed in a
physical way by Coppi and Sigmar (1973). [Q )= — fd vR' C (6.36)
One might be tempted to calculate the current density
directly from the definition, Eq. (5.65), since the odd (in where
o) part of H, is presumably given by Eqs. (5. 118) and A'„= —(v)(h/0;i) (v'/v, '„,. —5/2) f;() + G„ (6. 37)
(5.69) (Daybelge, 1971). This does not give the effect of
electron —electron collisions correctly, however. It is with G„ independent of 8, and' even (in v) in the trapped
necessary to include functions similar to introduced f„„ particle region. We use an expansion G„. =G +Q '
above in connection with the calculation of the diffusion +. , similar to Eq. (6. 16), and further assume
coefficient. This was first pointed out by Sigmar (1972),
and was demonstrated in a way similar to the above ana-
G'" =(v))Iz/& p) (v'/v(h; —y)f;. +Iz,
(6. 38) '",
lysis by Hazeltine, Hinton and Rosenbluth (1973). The where h; is "localized" in the same sense as was k„,
point is that, since the functions analogous to lz'0) and
(in the above analysis) are both O(e'~') they must both be
f, and y is a parameter to be determined. The lowest-
order variational expression is similar to Eq. (6. 18)
included in the "direct" calculation, Eq. (5.65). On. the
other hand, in the "variational" calculation [analogous (6.39)
to Eq. (6.23)] the X derivatives appear, and since &f„,/BX
is smaller than BIz' )/BA. by one order in E' ', f„, is not The zeroth-order collision operator is
needed. The inclusion of the Spitzer function in the CI;~ = v;;(v)Z, (6.40)
transformation Eq. (5.42), is necessary for obtaining this
"'
result, since the function h. would not otherwise have with the pitch-angle scattering operator given by Eq.
the property of "localization. " (5. 6), and
The electrical conductivity reduction coefficient,
(n„g„) has a numerical magnitude such that (for Z,. = 1) (6.41)
the effective conductivity
where v~, =(2T;/nz;)' ', Q(x) is given by Eq. (6.21), and
(6.34)
unphysically changes sign for aspect ratios less than 3 1/2T3/2 (6. 42)
about 4. In order to remedy this, the O(e) correction to i
y by minimizing Q», which is quadratic in y, as given where 0;, =Z; eB,/m;c. This poloidal rotation is a result
by Eq. (6.44) with Eq. (6.45). The value of y is then of collisions between trapped and untrapped ions, in the
found to be determined by presence of a temperature gradient.
The ion parallel flow velocity represents an additional
d3vv C(o) I, (0) =0 (6.46) source of friction on the electrons: the fourth electron
thermodynamic force is given by
This is equivalent to choosing y consistent with momen-
tum conservation (not automatically satisfied by the
pitch-angle scattering approximation CI; ). The numer-
ical value of y so obtained, and the minimum value of (6. 51)
') are
Q2[,
The O(e' ') correction may be neglected, since A„ is
y=1. 33, —(P„g„)=0.66n; . e'~'p', ()/~; (6.47) multiplied, in Eqs. (5.'lO)-(5. 72), by transport coeffi-
These same results can be obtained in a more straight-
cients which are themselves of order e' This force '.
forward way by replacing C;; by a model collision oper-
can be combined with A. „,
in the large aspect-ratio limit
because of the fact that n~ differs from a, only by terms
ator, suggested by Kovrizhnykh (1970), which consists
of order c. Hence,
of the pitch-angle scattering term C;; plus an ad hoc
momentum-conserving term. The details may be found (6.52)
in Appendix A of the paper by Rosenbluth, Hazeltine and
Hinton (1972). The Kovrizhnykh model operator can be and the Onsager relations, Eq. (5.98) imply further that
used to obtain the correct result, while other similar (6. 53)
operators can not [see, e. g. , Furth, el al. , (1970)] be-
cause it is self-adjoint, and therefore leads to a varia-
tional principle, which is equivalent to Eq. (6. 39).
In Eqs. (5.'lO)-(5. 72), therefore,
the same factor as A„,
„
A. is multiplied by
and we may use the combined
The ion thermal conductivity, as given by Eq. (6.47), force
is the result of pitch-angle scattering of trapped ions by
collisions with untrapped ones. [Recall, however, that A. „+%4,= —inn,
8
(1+ T;/Z; T,) + —lnT,
9
tion of the numbers of ions at a given value of r, with s;=-Qe'vg P,. v~) . ~
1
d(Reg, =
Oo
dPReg, = z, P, =( —v'/v', +5/2) 0;A —
2
(1 y(') . (6.'79)
we obtain
Recalling that the 8 derivative is at constant (, we note
(n„g„)= -(w ' ~'/ )4'enp', ~(B, v, „ /rB), (6. 71) that a term (v(u ~~/v,'h;) f;, may be added to the function g„,
which is independent of collision frequency (and hence without contributing to the left-hand side of Eq. (6. '78),
independent of Z;). This result was first obtained by provided that Su„/&8=0. In order that this equation have
Galeev and Sagdeev (1968), by using a Krook model col- a solution at all, the right-hand side must be orthogonal
lision operator. [The numerical value is correctly given to the above null solution. That is, multiplication of Eq.
in the paper by Galeev (1971).] The above derivation, (6. 78) by v'( and integration over 8 and v yields the con-
using the pitch-angle scattering model, was first given straint
by Rutherford (19"70).
The other electron transport coefficients may be ob- 2 e d8sin8 d'v v' $(1 —P) " =0, (6.80)
tained in a very similar way; the results are
~/2
which serves to determine the arbitrary parameter u]].
(o„g„) = —-~ e'n, p,'~ (v, „, B,/rB), (6. 72) This parameter can be introduced, in terms of the pa-
rameter y just as in the banana regime [c.f. Eqs. (6. 3'7)
(6.73) and (6.38)]:
' e'(n, c/B~, ) (v, 7, B„,/rB), (6.81)
(u„g„) = —(0.26+0. 24/Z;) „,
(6.'74) where II; may be assumed to be localized. Substitution
' into Eq. (6. 78) and use of Eq. (6.40) for the collision
( og )3= —(0. 13+0.27/Z;) e'(n, c/B„,) (v,„,T, B,/rB), term gives
which is equivalent to the condition of zero ion flux, lute maximum of V„at the maximum possible rate.
when I+(' (in the expression for the drift velocity) is The value of V', ",' so obtained (using the arbitrary con-
replaced by unity. Comparing this with Eq. (6.84), we vergence criterion that the fractional change in I/'yy
see that the 2 in the integrand can be omitted; alterna- after one iteration, be less than 10 ) gave the numer-
tively, the & can be replaced by y, making the ion ther- ical values of the transport coefficients. These mere
mal conductivity quadratic in y, as in the banana regime calculated for three inverse aspect ratios, a=0. 01, 0.04,
variational calculation. The condition which determines and 0. 09, and for a range of values of collisionality,
y, Eq. (6. 85), is then equivalent to minimizing the ther- 10 '& v, & 10, where the electron collisionality param-
mal conductivity with respect to y. The results of carry- eter is defined as
ing out the integrals in Eq. (6.84) and determining y in —v 2 xBo/(B, v, z, q' ~'),
v~, = „, (6. 87)
this way, are
y=3, (P„g„)= ——,' v'~' e'n; p&'e(v, „B„/~.B) . (6.86) (the banana-plateau transition occurs nominally for
v„, -1). These results were then'fitted to analytical
The thermal conductivity given by this expression is formulas (which will be presented in Sec. VI. F) using
the result of scattering of "resonant" ions (in the sense the method of least squares.
described in Sec. III.D) by nonresonant ones. The fric- The starting point for the ion thermal conductivity cal-
tional force of reaction, back on the nonresonant ions, culation, is a variational principle equivalent to Eq.
produces a parallel f lorn related to y in the same may (6.82), which may be obtained by substituting Eq. (6.81)
as described in Sec. VI. C; Equations (6.49)-(6.51) of into Eq. (5. 96) for m=n=2, c=i, and using the pitch-
that section remain true in the plateau regime as well angle scattering operator, Eq. (6.40), on the localized
as the banana regime, with only the numerical values of function II;. The numerical method of calculation, using
y being different in the two regimes. Likewise, Eqs. the maximal variational principle, is the same as for
(6. 54) and (6.57) hold also in the plateau regime. the electron transport coefficients. The value of the
parameter y is obtained, as in the banana regime cal-
2. Variational calculation culation, by minimizing the ion thermal conductivity,
which is quadratic in y. Analytical formulas, obtained
The collision-frequency dependence of the transport from least squares fits to y and (P„g„), as functions of
coefficients, in the entire transitional range of collision ion collisionality, mill be presented in Sec. VI. F.
frequencies between the banana regime and plateau re- Analytical checks on the accuracy of the numerical
gime, is quite relevant to present experiments. This calculations, were made in the following way. We de-
dependence was calculated by Hinton and Rosenbluth fine the functions Q „(v) in terms of the variational ex-
(1973). They carried out a numerical solution of Eqs. pressions, Eq. (5.96), before integrating over particle
(6.63) and (6. 82) by a finite difference relaxation method, energy; for example,
which makes use of a maximal variational principle to
gain rapid convergence. The transport coefficients were
calculated from variational expressions, so that only the
localized parts of the distribution functions were needed,
and the pitch-angle scattering approximation mas ade-
quate. Good agreement was thus obtained with the ba- (6.88)
nana regime limiting values which had been calculated where v is defined as a function of & by Eq. (6.65).
analytically, as described in Secs. VI. B and C. (Since v appears in Eqs. (6. 63) and (6. 82) only as a pa-
The starting point for the calculation of the electron rameter, the matrix elements Q „ themselves provide
transport coefficients, is Eq. (5.96). The functions g„, variational principles for these equations. ) In the ba-
are assumed to be localized in the variable X, so that nana regime limit, v«e' ', Q»(v) must approach the
the collision operator can be replaced by the pitch-angle constant value (w/2) Ie ' ' =3.07 [where I is given by
scattering operator, Eq. (6. 19). [This variational prin- Eq. (6.25)]. In the plateau regime limit, e' '«v«1,
ciple is then equivalent to Eq. (6.63), when the trans- the function Q„(v) is given approximately by
formation to the variables (8, g) is made. ] Using the
pitch-angle scattering approximation, it is possible to vQ„(v) ~ (v' e'~'/8) [1+0.32(v/e' ') ~~'] (6.89)
eliminate the two functions g„, g, , (for n=1, for ex- where the factor in square brackets includes a correc-
ample) in favor of a single function C„. This function
is like a "potential, " in that both g,', may be derived
tion to the plateau limit. This was obtained by solving
Eq. (6.68) to one higher order in e/v' ', and substituting
from it by differentiation. the result into the variational principle. Good agree-
The variational principle, mhen expressed in terms of ment with these analytical results was obtained, as
this single function, is maximal. That is V,', has an shown in Fig. 9. Good agreement was also obtained mith
absolute maximum when the function C„satisfies the the result of a boundary layer calculation, mhich we
appropriate partial differential equation and boundary now describe.
conditions. This property was found to be quite useful,
in combination with a finite-difference relaxation meth-
od of solving the partial differential equations on a com- 3. Boundary layer calculation
puter. By choosing the relaxation parameter to maxi- Corrections to the banana regime result mere obtained,
mize V,', after each iteration, an increasing sequence by a careful treatment of the boundary layer between the
of values was generated, which converged to the abso- trapped and untrapped regions of phase space. As we
dx sg'„/s y, dx xsg,', /s y,
.50
to the exact value s, whic h can be obtained f rom Eq.
(6. 92). The result is
(6. 94)
nana regime limit, Q~~ =3.07, and the plateau limit v'Q&& an amount proportional to the shift. (The actual broad-
= 71 ~/8. The correction to the plateau limit given by Eq. (6.89), ening of the function g'„contributes negligibly. ) We thus
and the result of the boundary-layer calculation, Eq. (6.97), find the boundary-layer contribution to the diffusion co-
are also shown. efficient to be given by Eq. (6. 88), with
at the outset. Then Eq. (6. 63), with v„=O, can be re- v only. The next-order equation is
not forced to have the form of a Maxwellian with per- collisions conserve momentum. By expanding in powers
turbed (8-dependent) values of density and temperature. of
Thus, Eq. (6. 101) is not a perturbed Maxwellian, and
this is the reason for the numerical discrepancy in the
diffusion coefficient: in a sense, the result is a "non- and writing the 0 dependence of the first-order term ex-
equilibrium" diffusion coefficient. plicitly
We conclude that the numerical, results shown in Fig.
10 are unphysical in the large v limit for two reasons.
Firstly, the pitch-angle scattering representation of elec- we obtain
tron —electron collisions [contained in Q»(v) if v is given
by Eq. (6. 65)] is not adequate in the collision-dominated igh„—(xB/B~, v) C„h„=(v'/v, '„,. —y) ev(1+ (') f;, .
regime. We recognized this by specializing to the case 2Q;p
Z; »1. Secondly, even when Z; »1, the absence of any (6. 105)
process, which would tend to Maxwellianize the distri-
An approximate method of solution will be described
bution function, invalidates the results for the trans-
port coefficients. Thus, the numerical results of Hinton subsequently.
and Rosenbluth (1973) are useful only for the banana- The ion thermal conductivity is given by
plateau regime.
(P, r Z2i) = 2g f,
1
o 0
d' v —(v '/v,
2
'z; —5/2) (1 + P) Heh„
E. The plateau-collisional transition
(6. 106)
We now consider the regime of intermediate to high
The constraint, Eq. (6. 80), which determines y, may
collisionality, v~»1, where v„, is given by Eq. (6.87)
be written as follows, after integrating by parts
and v~; by a similar expression. In the large aspect-
ratio, circular cross section case which we are con- d'v v'(1 —3P) Reh„= 0 . (6. 107)
sidering, the 0 dependence of the distribution functions
is simply sinusoidal, so that only the dependence on the We note that the mirror force term (proportional to
velocity variables v, $ remains to be determined.
esin8) in Eq. (6. 82) has been neglected, in obtaining Eq.
(6. 105); it has been used in obtaining the constraint, Eq.
1. E ectron transport
I
(6. 107). In the intermediate-to-high collisionality re-
We begin with Eq. (6. 63), with the exact collision term gime, the majority of particles (the well-untrapped
C„ in place of the pitch-angle scattering term. By ex- ones) which contribute to diffusion are only slightly af-
panding in the inverse aspect ratio fected by the mirror force, so it can be treated as a
+(x) ++(2) +. .. perturbation. In the absence of any parallel magnetic
field gradient, however, the ion parallel flow velocity
and writing the 6 dependence of the first-order term would be arbitrary, corresponding to the choice of
explic itly reference frame. The small mirror force removes this
g~" = Im(h„, e'~), degeneracy, in a toroidal confinement system, and en-
ables the parallel velocity [and hence y, in Eq. (6. 81)]
we have (for n = 1, for example) to be uniquely determined.
ev(1+ P)
igg„, —(rB/B„,v) C,e 0„=
je — 2~g f«,
eo ~ (6. 103) 3. Model like-particle collision operator
~
where C, is given by Eq. (1.19). The expressions for the The above kinetic equations [Eq. (6. 103) and its n =2
inner products now involve integration only over velocity; and n =3 counterparts, and Eq. (6. 105)] were solved by
for example, the diffusion coefficient is Rawls, Chu, and Hinton (1975) by using a model like-
particle collision operator. Since this operator may be
-(o.„g„)=
Q„ Bo
2
d'v —(1+ $') Heh„,
2
(6. 104) generally useful, we give some of the details
C;;f;, =-~;, (~) [&f~, -f;, +~&(~ii/~, ~) Pf~.
I l
Eq. (6, 103) would still be difficult, because of the com- + pf;, +s(v'/v', . —5/2) f;, ], (6.108)
plexity of the Fokker-Planck collision operator for elec- where Z is the pitch-angle scattering operator given by
tron-electron collisions. An approximate method of Eq. (5. 5). The collision frequency is a function of vel-
solution will therefore be discussed subsequently. ocity o, def ined here as
lisions with other ions (Trubnikov, 1965). The values of 1. E lectron fluxes
P and & were determined by using the corresponding
model electron-electron collision operator to solve the
Rather than fit the inner products g„,), it was (a,
found more convenient to fit a dimensionless set of elec-
classical parallel transport problem of Spitzer and tron transport coefficients K „, defined as follows
Harm (1953). By finding a best fit to the values of x„,
II.'„, and ~» given in Table II, for Z; =1, the values —(a„g, ) =n e (peel&e)A&» (6. 114)
P =0.96, A. =0. 55 were obtained. The parameters P, p,
2(au gee) = e
(aI& g3e) (pe e/ e) +(2 (6. 115)
and s are to be determined by the requirements of con-
servation of particle number, momentum, and energy
in like-particle collisions
. . .
—(a„g. ) —5(a„g. ) —~(a„g, ) = n. ~"(p'. e/~. ) ~„,
(6. 116)
(6. 110) —(a„g„)= e'/'(n, c/B„,) K», (6. 117)
—(a„g„)——,' (a „g„)= e' /'(n, c/B~, ) K», (6. 118)
It can be shown that this model collision operator is
self-adjoint and satisfies the II theoreme
A convenient method of solution of Eqs. (6. 103) and
—(a„g. ) = . ~"(~((/T. ) A;. . (6. 119)
where e =r/Bo, p', e =2m, T, c'/e'B~'„v'e is given by Eq.
(6. 105), with the model operator C„or C;;, involves
expansion of the unknown functions, h„and k„, in a
(5.4), and o(( is given by Eq. (5.66).
In terms of these coefficients, the electron transport
series of Legendre polynomials P,((). The expansion relations, Eqs. (5.70)-(5.72) may be written as follows,
coefficients satisfy recursion relations which can be
using the weak-coupling approximation and Eq. (5. 83)
solved by a computer, assuming the series to be trun-
I.
cated at a large value of l. The self-consistent solution 8
l2
for the parameters P, p, and s involves only the coef-
ficients of P, (() and P, (g), while the transport coeffi-
cients are obtained from the coefficients of P, (g) and —K, n, e'/'c(E„/h)/B, , (6. 120)
P, (().
The numerical results obtained by this method agree
fairly well in the plateau limit with those given in Sec. e. + T. ('. = —n. T. E'~'(p'. ~/v, ) K„A, +K„—
—, '.
22 ar (nT
e' ~~'(
2
c~'„/5
~ e'~ ~~'~ K ~~(K )tl+c
n
«~
~ m no
~I eI
2 ~
~~ ~e
Ie
3 2 ~
~/~ „v~ „v~ E ] (6. 129)
~e
The results of Rawls, Chu, and Hinton (1975) for e«1
~
mn m~
~~
(6. 125) were fitted to such expressions; their values for & „
and c „are listed in Table III (except for b» and c», as
For n=3, mentioned previously). Again, note that these values
were obtained from least-squares fits, so they do not
exactly agree with the collisional regime values given
(6. 126) by Eqs. (5. 151)—(5. 153).
The numerical coefficients, K „, & „, & „, and c „, were Since these two limiting cases represent regimes of
obtained by considering the following two limiting cases. collisionality which overlap, when e «1, we may regard
(i) v„, «e ' ' |banana —
p/a/eau regime) the expressions given in Eqs. (6. 125) and (6. 126) to be
valid for all values of the collisionality parameter v~, .
(6. 127) It then seems reasonable to extrapolate the results to
The results of Hinton and Rosenbluth (1973) for e« 1, finite values of &. In Fig. 11, we show the result of
were fitted to such functions. The term involving v~, ' using the interpolation formula, Eq. (6. 125), for a value
reflects the role of the boundary layer, as mentioned a=0. 2, which makes it possible to see all three colli-
previously. The values of K o and a „so obta. ined are sionality regimes on the same linear scale. The diffu-
listed in Table III. Note that the values of K '„were ob- sion coefficient —(o, , g„), which is proportional to
tained from least-squares fits, in the range 10 '& v~, v~, K»(v„,), is plotted; the ordinate gives the diffusion
& 10; they do not agree exactly with the banana regime coefficient normalized to its plateau value (called D in
values, from Eqs. (6.28) —(6.33), which apply in the limit the figure). Also shown for comparison are the banana,
v~, -0. Also listed in Table DI are the values of b33 and plateau, and collisional regime asymptotes, the ba-
c» obtained by Hinton and Rosenbluth (1973), since these nana-plateau transition formula, Eq. (6. 127), and the
are the only known values of these coefficients. plateau —collisional transition formula, Eq. (6. 128).
(ii) v~, »1 (p/a/eau —eo//isiona/ regime): For m or
n=1 or 2, 2. ion fluxes
A dimensionless ion thermal conductivity coefficient
TABLE III. Numerical coefficients in Eqs. (6.125) and (6.126):
the results of least-squares fits to the transport coefficients. K, may be defined as
(P2 g21) =n e' (P e/&~)~2 (6. 130)
A. Zi= 1
'
where pPe =2m; T; c'/Z, e'R„'„adnT; is given by Eq.
&mn
(6.42). In terms of this coefficient, the ion heat flux,
11 1.04 2.01 1.53 0.89 Eq. (5.84), may be written as follows (in the weak-cou-
12 1.20 0.76 0.67 0.56 pling approximation)
22
13
2. 55
2. 30
0.45
1.02
0.43
1.07
0.43
i.o7 2 i K~iei » g+ (P„g,;) &; I./Z;
23 4. 19 0.57 0.61 0.61 1 + +2 +3
33 1.83 0.68 0.32 0.66 (6. 131)
B. Zi =2 Here we have used the Onsager relation (P„g„)
= —(P„g„.). We have used the approximation of Eq.
(5.75) for A„, which gives, on comparison with Eq.
11 0.86 2. 18 1.17 0.79 (5.70),
12 0.95 0.78 0.50 o. 51
22 1.99 0.46 0.26 0.34
1
I./Z,
2 (6. 132)
13 1.87 0.89 0.62 0.69 1 + p ge Q3
low s
Z«~ =0.66, a, =1.03, &, =0.31, c, =0.74 .
The mean ion parallel velocity, averaged over a mag-
netic surface, is obtained from Eqs. (5.46) and Eq.
(5.83), in the weak-coupling approximation:
cT, Z;e S(4)
Zea p"" BlnT;
sr
B)nP;
sx r, ar I
(6. 134)
The coefficient (P„g„), is equal to 1.17 in the banana
regime, —0. 5 in the plateau limit, and —2. 1 in the col- 0 I
. .
&./&=[(&. g2. ) —(~„g. ) (~, g. ) (P„g,~) (Z'. /Zl &)1, '- Z'
+[ (~, g„) —(%, g„)' (P „g») (Z', /Z» I ) + ( a.„g„)(P „g„)(Z, /Zn E) ] (6. 139)
lnP, Z, lnP,
qb/T =(o'„g„)(P „g,~) (Z, /Zn I') V
+ [ (P2, g, n) + (W, g„) (P „g») (Z, /Z, 5') + (o'„g„)(P „g'»)' (Z ', /Z'n P) ] (6. 140)
where 1 = 1+(Z2/Z~)(o3, g~, )(P~, g»), and we have used by transposing the results given earlier in this section,
T, = T, =—T, and the Onsager relations, Eqs. (5.98) and i.e. , by replacing m, by m„m; by m„etc.
(We note,
(5.99). however, that the required inner products are not yet
The inner products may be obtained in special cases available in the plateau-collisional transition regime. )
TABLE IV. Estimates of H —f, whose smallness is required for the weak-coupling approxi-
mation to be valid.
Electrons in Electrons in
banana r egime collisional r egime
.
(~, «~) )) ~ -3/2)
+b Zb e IQA
Tbb 3 +
JVg b
i/2 m 3/2 (6. 142)
)
(o.'„g„) (E- 1) slnT
(6. 151)
The charge-neutrality condition n, =Z; n; has no counter-
Bg
part in the two ion species problem; the ion charge num- &b=(P„g.«) 8
(6. 152)
ber Z; in certain expressions must be replaced by the
parameter (Connor, 1973) where only the larger, O(mb/m, )' ' term has been re-
n = n, Z,'/n. Z', . tained in qb. These expressions agree with those given
(6. 143)
by Hinton and Moore (1974), when the proper identifica-
The banana regime coefficients, for example, Eqs. tion of coefficients is made. For the required inner
(6.28)-(6. 30), are expressed as continuous functions of products, it is convenient to use the approximate fits
n by replacing 1+0.53/Z; by 1+0.53/o. , etc. obtained by Hinton and Moore (1974), in which the re-
Eqllatlolls (6. 138)—(6. 140) Illay be sllllpllf led collslder- suIts are expressed as continuous functions of e. We
ably when both ion species are in the collisional regime. first define coefficients A'„and A „, as follows:
Using Eqs. (5. 155), (5. 156), and (5. 164), we find
p2 ~i/2
(6. 153)
)I Blnp4 Z4 &lnpb ~lnT
) (6 144)
( (
.
(~„g, ) =- Oa Pap+abE
i/2
(~i2 —~ I~»)
a 5 a
(6. 154)
g,
T' =(~„g, . ), 1np, ' —Z, ~lnpb'
Z', +(~, g. ) . 1nT (6. 145) ( Z~ g24) =
2
Tab
~1/2
(~22 5~12 + 4 ~IX) ~ (6. 155)
and (6. 151), was found to be coefficient much smaller than the standard one, by the
& =1+(2m. /m„)" Z,', /nSC, ', .
factor I/E= (n, /n. ) (m„/m. )'~'.
(6. 158) The continuity equation, in which I; appears, has the
The coefficients are given approximately as follows: form
0. 73 (1 + 0. 53/o. )
1+1.26 v~, (1 +0. 53/o) (6. 159)
0. 73 (I + 0. 71/o. } where 8, is the difference between the ionization and re-
' (6, 160) combination rates which affect species &. Neglecting
1+0.37 v„(1+0.71/o. )
these effects, and in the absence of a temperature grad-
1.46 (1 + 0. 8/o. ) ient (& T/&r = 0), diffusion leads to a final state (1, = 0)
(6. 161)
1 + 0. 2 v„, (1 + 0. 8/n) in which ~Inn, /sr = (Z, /Z, ) sinn~/&r or
0. 55 g 0.73 n~ const (n ) b~
1+094v ' 1+037 v
(Spitzer, 1952; Taylor, 1961; Braginskii, 1965). If
1.65 species & is hydrogen, and Z&» 1, then the impurity
K (6. 162)
species "&" is much more highly concentrated in the
where center of the plasma (+=0) than the hydrogen ions, in
this final state. The effect of a temperature gradient is
~*a = ~2~+0/(+po Uu 'ra~ & (6. 163) to increase the inward diffusion rate for impurities, if
both ion species are in the collisional regime. As
pointed out by Rutherford (1974), however, the coeffi-
u„~ =- v 2 xylo/(B„, v, » v'» e'~') = Z', v„,/v 2 Z', , (6. 164) cient &T/9r, in the banana regime expression for F„ is
with w, ~ and 7'» given by Eqs. (6. 141) and (6. 142), and positive, so that a normal temperature gradient
=(2&/~, )' ', vug, =(2T/m~)'~'. The numerical coef- (&T/&+&0} causes inward diffusion of the light species,
ficients in Eqs. (6. 159)-(6.162) were obtained by: (i) re- and hence ouAeaxd diffusion of the impuxity species,
quiring the limits v~, -0, v~~-0 to give the exact banana since F~ = —(Z, /Z, ) I",. As the collisionality is increased
regime values, and (ii) fitting the remaining parameter in the banana-plateau transition, this temperature grad-
to numerical calculations similar to those of Hinton and ient "screening" effect, changes sign (Hinton and Moore,
Rosenbluth (19'l3), with the constraint that the plateau 1974), leading to' results qualitatively like those in the
limiting values be independent of n. The terms pro- collisional regime.
portional to 1/o. in Eqs. (6. 159)-(6.161) are the contri-
butions from like-species collisions. In Eq. (6. 162}, the Vll. Closure of the IVloment Equations
assumption o. »(m, /m~)'~' was used; the effect, on the
The moment equations, discussed in Sec. II, become a
heavy species, of unlike-species collisions is thus ne-
closed set of equations only when the appropriate kinetic
glected here, although it is retained in Eq. (6. 158}.
equations have been solved and the required moments
It is now easy to check on the weak-coupling approxi-
have been calculated from the distribution functions. Xn
mation, using Eqs. (6. 158), (6. 159), (6. 162), and (6. 164).
this section, we discuss the closure problem for tokamak
If both ion species are in the banana regime, then
Ii —1-(m, /m„)' '/o. , which is typically small [except in devices, in which there is an additional difficulty: the
the uninteresting case o. & (m, /m~)'~'j so that the weak- magnetic field is not completely determined by currents
external to the plasma, but must be determined self-con-
coupling approximation is valid. However, if the light
sistently with the transport processes in the plasma
species is in the banana regime, but the heavy species
is in the plateau regime (note that v~, vq„assuming » (Grad and Hogan, 1970; Grad, 1970). The goal is to ob-
tain a set of equations which, with appropriate initial
Zg && Z~) then
data and boundary conditions, completely determine the
S —1=1.25(n. /n, ) (~/m, )'~' v, (1+0.53/u) . . spatial and temporal dependence of the electron density,
the electron and ion temperatures„and the magnetic fieM
Hence, if
v„. (1+0. 53/o. ) ~ (n, /n. ) (m, /m. )" (6. 165)
in the plasma.
As we have seen in Secs. IV-VI, the magnetic surfaces
the weak-coupling approximation is not valid. The provide the basis for the coordinate system in which the
right-hand side of the inequality is typically small; for transport relations —
linear relations between "fluxes"
example, with as much as 5% oxygen in a hydrogen plas- —
and "forces" are given. In fact, the fluxes are aver-
ma, it has the value 0.2. If both ion species are in the ages of certain odd moments over these magnetic sur-
plateau regime, then I'-
1= (n, /n~) (m, /m~)' ', which is faces, while the forces are derivatives of even moments
typically large: the weak-coupling approximation is with respect to a distance coordinate perpendicular to
again not valid. As shown by Hinton and Moore (19'l4), the surfaces. (Although the transport relations were not
the right-hand side of Eq. (6.165) gives the value of col- derived in precisely this form, the quantity (E ~~B) having
lisionality at which the banana-plateau transition actu- been treated as a "force", we will show, shortly, how
ally takes place (i. e. , where the diffusion coefficient be- they may be put in this form. ) Since the magnetic field
comes approximately independent of v~, ). This value configuration in tokamaks depends on the plasma pres-
may be much smaller than the "standard" value for this sure, which changes on the diffusion time scale, we need
transition (v~, -1), which makes the "plateau" diffusion an equation to determine the time dependence of the mag-
„„Bt
magnetic fields, in the presence of a constant external times which differ by &t, divided by &t. By using dS
toroidal magnetic field. It is then more natural to con- = (dS/IV%'I)V4, we may write the above result as
sider the poloidal magnetic field, or g, to be determined
by Faraday's law, Eq. (2.103), with an electric field con-
8
at ... d'x5'(x, t) = d'x —
&F
+
dS
vei
I u~ V4.
~
v(t
—
E B
~
(7.11)
tive minor radius coordinate p, we must turn to the de- I I
rivation of an equation for the poloidal flux, ()). First, Combining Eqs. (7.9) —(V. ll) gives, therefore,
however, we shall derive a theorem which will be
needed in the following. 2m —
8$
at, = —c
dS
vyi
E ~
B,
This theorem gives the time derivative of an integral I
J, defined by = —c BV (Z — a) (7.12)
(7.4) where the second equality follows from Eq. (2. 56). From
where the integral is carried out over the region inter-
Eq. (V.8) we finally obtain
ior to a moving surface 4'(x, t) = const. The result is
(7.13)
= lj d'x —+ &p dS u F'
~
(7.5) where we have also used the relation
where the surface integral is taken over the surface 4
=const, and where u~ is the velocity of the constant +
which follows from the definition, Eq. (2. 50). The rea- c &I 2 C 8 t
BV
son for using p (which is equivalent to using P) as the 4~ eV By
spatial coordinate should now be clear: Eq. (7.13) is a
simple and exact equation, which contains the electric which may be written as
field only through (E))B&, which is obtainable from the c 2B Xeg
transport relations. "»B&=4.'BV I BV
It is actually more convenient, in applications, to use
an expression for the time derivative of the effective po- where X, as defined by Eq. (7.3), may also be written
loidal field B„(p, t), defined by Eq. (5.36). To obtain as
this expression, we simply differentiate Eq. (7.13) with
respect to p
(V. 14) When the effective minor radius p is used as the spatial
coordinate, we have
where cI' B ' X Bg)
4)TV'(p) Bp IV'(p) Bpg '
Eo = (E))B&/IRo(R (7.15)
This equation has the form of Faraday's law for a cylin- where V'(p) =—
BV/Bp, or
drical system, with Eo as the axial component of the el- cI B XRO
ectric field. 4m V'(p) Bp IV' (p) (7.20)
Next, we derive a useful expression for the velocity,
u, , of a surface of constant p (or constant Q). We note where B„is the effective poloidal field, defined by Eq.
that (5.36).
.V C. Transport relations
u P Vp =u P ~ V
~ B@=u
Bp
We now wish to make two points, with regard to the
where relationships expressed by Eqs. (5.70)-(5.73). The first
8 is quite simple: although the electric field appears in
u ~ V@=— (7.17) the role of a force, through A„— = (E B)&B /(oB ), and the
current density appears as a flux, these roles may
and the time derivative is taken at constant x. Also, easily be interchanged. That is, by solving Eq. (5.72)
since „,
for A. in terms of (J„/h) and A„, A, „A~,, and sub-
stituting the result into Eqs. (5.70) and (5.71), we obtain
transport relations in which (E„B& appears as a flux,
and (J„/11& (=(J„B&/Bo) appears as a force. Note that the
we have Spitzer current term in Eq. (5.72) is, by definition
(J„,/I &
= ~„(E„B&/B,.
(7.18)
It is natural to regard (O„B& as a force', since it is
Finally, combining Eqs. (7.16), (7.18), (7.10), and given in terms of spatial derivatives by Eq. (7.20). It is
(7.13), we have also natural that (E„B& should be regarded as a flux:
recall Eq. (2.115), which gives the Poynting flux rela-
u, Vp = [c(E„B&/I(R &
—cRE z, ]/R~B, . (V. 19) tive to a constant-P surface, averaged over the surface.
The second point which we wish to make is that the
B. Ampere's law particle flux, relative to a, constant-@ surface, is simp-
ly related to the flux I', , defined by Eq. (5.62). [Recall
The transport theory described in earlier sections al- that Eq. (2.108) gives the particle flux relative to a
lows us to relate (E„B& to (J, B& (as well as to spatial constant-g surface. ] To find this relation, we begin
derivatives of the density and temperatures). We now with Eq. (3.81). Noting that the essential difference be-
show how (J B) can be expressed in terms of spatial tween (the radial components of) vD, Eq. (3.17), and v„,
derivatives of |t). We have Eq. (3.24), is the inductive E xB drift, we obtain, from
(J„B&=(J, B,&+P, .5, &
Eq. (3.81),
where the subscripts p and T denote the poloidal and (&,u, ' Vp&— o Vp + 'fIac
toroidal components, respectively. Using the poloidal
and toroidal components of Ampere's law, Eq. (2.42), (7.22)
+(n, u, Vp&,
where E'~'— = -c 'BA/Bt is the inductive electric field.
From now on, we shall neglect the classical flux
(nu, . Vp), although it is easily included.
J ~ vg= —v x B
4m ~
~ vr„= —v (vg/R
4m
) The first term on the right-hand side of Eq. (7.22) can
be transformed, by using Eq. (5.42), in which the indef-
we have inite integral can be eliminated by using Eq. (3.32) and
The Ohmic heating term in Eq. (V.34) can be trans- ' I',
formed to a more useful form by using Eqs. (2.69) and {u, (F, +n,. Z, eE)) =
~ '+(y ——,')n, /n, . (7;43)
(2.89) This expression was first obtained by Hinton and Rosen-
J E = —(c/4z)(BI/Bt))E„B —c(BP/B)t))E R . (7.38) bluth (1973); the banana regime limit, in which y =1.33,
was given by Rosenbluth, et al. , (1972).
The derivative BI/Bp may be eliminated by using Eq.
(ii) In the collisional regime, A„, as given by Eq.
(2.90), in which the right-hand side is (4m/c)R2J. &7f; (5.85), is negligibly small [cf. Eq. {5.156) and the sub-
multiplying by R ' and averaging over a magnetic sur-
sequent discussion].
face yields We now collect results, and write the general surface-
BI (J Vg)+ cBP/BP
~ averaged electron energy equation as
') (7.39)
B&))
By manipulations
I(R
similar to those which lead to Eq.
3 Q 8 pl
p)p, ]+~, Bf + —V'(p)[q, + 'T, (,(, — —
)~
(7.44)
where X is defined by Eq. (7.3). By averaging Eq.
(7.38) and using Eq. (7.19) to eliminate (ErR), we find where Q,. is given by Eq. (7.37), and (J Vg) is given by
(J E) =(J. Vf)(E„B)/I(R-')+(u, . Vp)BP/Bp. (7.41) Eq. (V.40). The ion equation is very similar, the differ-
ences in the right-hand sides being an overall sign
The Ohmic heating term is thus related to the relative change and the absence of the (E, B) term.
Poynting flux [recall Eq. (2.115)] and to the work done In the large aspect-ratio, circular cross section case,
by the moving surface against the plasma pressure with magnetic surface motion neglected, Eq. (7.44)
gradient. takes a more familiar form
The remaining terms on the right-hand side of Eq.
(7.34), involving u, , are evaluated as follows. Noting
that only the O(6') terms are required, and making ex-
1 8
3 Bpe'+ — —~[q + ,'T
e 2 r]-
e e
for
and
the
U y ln Sec. II. C.2; the result may be writ-
—
Q.i + ~ '
+ ,'T
I'e P'- 8
r.i ]—i i.T.i A li. (7.47)
equilibrium rotation. It follows in particular that the e gp
radial electric field never appears in the closed set of
second- order moment equations. The meaning of the factor (V')'~' in the time derivatives
The right-hand side of Eq. (V.42) appears in both elec- can be made clear by considering an example in which
tron and ion energy equations, with opposite signs in transport processes are negligible.
each equation, and represents a transfer of energy be-
tween the species. It can be simplified in either of two F. Adiabatic compression
limiting cases: Gne successful method of heating tokamak plasmas
(i) In the large aspect-ratio circular cross-section has been compression in major radius (Bol, et al. ,
case, with electrons in the banana-plateau regime, A„ 1972). By changing the external poloidal magnetic field,
is given by Eq. (6.55), in terms of the particle flux. with a time-independent toroidal field, the equilibrium
If the ions are also in the banana-plateau regime, then position of the plasma column can be shifted to smaller-
p,. is given by Eq. (6.50), in terms of the ion tempera- values of major radius. U this is done fast enough that
ture gradient and the parameter y. Neglecting motion dissipation and transport effects are negligible, but
of the magnetic surfaces, (u, 7'p) is given by Eq. (7.27) slowly enough that magnetohydrodynamic equilibrium is
divided by n„so that Eq. (V.42) reduces to approximately maintained, the plasma is adiaba3tical/y
compressed. In this case, we assume that the time de- obtained by first solving the kinetic equations, Eqs.
rivative terms dominate in Eqs. {7.13), (7.30), (7.46), (5.80) and (5.81). Thus, the closure problem is solved
and (7.47), which become in principle, although explicit results for the transport
coefficients are known only in special cases.
(7.48) The most useful special case is the large aspect-ra-
tio, circular cross- section case, for which the geometry
—
8
8t
(n e
V') =0, (7.49)
was discussed in Sec. VI.A. Only in this case are the
transport coefficients known explicitly, for all regimes
of collisionality. In order to facilitate the use of these
—[P, (V')' '] = &
Lt',.( V')' '] = 0 . (7.50) results in computer programs, we now summarize the
complete set of equations and transport coefficients for
If the temperature equilibration occurs more quickly the large aspect ratio, circular cross-section case.
than the compression, Eq. (7.50) should be replaced by The magnetic surface configuration is assumed here to
be time- independent.
—
Bt P
(V')'"] =0, (7.51) The equations which (with appropriate boundary and in-
i—
8
n„T,
, T, , and functions of
& and t are as follows (with p, = n, T, and—P, =n, T,
B„as
):.
constant toroidal flux. By differentiating with respect '+
e — xl" =0,
8t
to Q and using the definition of the safety factor q, Eq,
(2.50), we have q =const. , following the flux surface
motion. Equation (7.49) gives n, V' =const. which ex-
3 8pe'+ 1 8
—.
8&
—
~(q e +-,~'T e I"e )
presses the conservation of the number of particles be-
tween two nearby magnetic surfaces, while Eq. (7.51)
yields P(V')'~' = const. which is the law for the change
gpss
—
3~e' Qe'(T —T.)+ cEp
e
&e i
8—
47t~ 8y
(~B&0) —
—Ie'
+e
8p&
8&
+ p. i T i A. .
at each instant of time by Eq. (2.90), using Eqs. (7.39) and (6.136), respectively.
and (7.40) to express 8I/8$ in terms of spatial deriva. —
ACKNOWLEDGMENTS
tives of the pressure' P and the poloidal magnetic field
B„. These equations describe the radial profiles and This work was supported in part by U. S. Energy Re-
search and Development Administration through Con-
magnetic surfaces on time scales long compared with
the time needed to establish magnetohydrodynamic equi- tract AT-(40-1)-4478 and also in part by the Alfred P.
librium. We have shown that, when the transport is Sloan Foundation, through a Research Fellowship held
negligible, the equations are consistent with the scaling by one of the authors (F.L.H. ). Figure 11 was kindly
laws for adiabatic compression. supplied by D. W. Ross.
We have shown that these equations are a closed set REFERENCES
.of equations when the following quantities are related to
Alfven, H. , 1950, CosmicaE Electrodynamics (Oxford U. P. ,
the density, temperatures and poloidal field and their Fairlawn, N. J.).
spatial derivatives: I', , q, , q, , (E B), p, „A„.. But these Artsirnovich, L. A. , 1972, Nucl. Fusion 12, 215.
relations are given by Eqs. (5.70)-(5.73), (5.75), and Banos, A. , Jr. , 1967, J. Plasma Phys. 1, 305.
(5.76), in terms of transport coefficients '.ch must be Berk, H. L. , and &. A. Galeev, 1967, Phys. Fluids 10, 441.
Rutherford, P. H. , 1970, Phys. Fluids i3, 482. Spitzer, L. , Jr. , 1967, Physics of 5'ully Ionized Gases (»«r-
Rutherford, P. H. , 1974, Phys. Fluids 17, 1782. science, New York).
Rutherford, P. H. , and E. A. Frieman, 1968, Phys, Fluids 11, Spitzer, L. , Jr. , and R. Harm, 1953, Phys. Rev. 89, 977.
569. Stix, T. H. , 1973, Phys. Fluids 16, 1260.
Rutherford, P. H. , L. M, Kovrizhnikh, M. N. Rosenbluth, and Stringer, T. E., 1970, Phys. Fluids 13, 810.
F. L. Hinton, 1970, Phys. Rev. Lett. 25, 1090. Su, C. H. , and C. Oberman, 1968, Phys. Rev. Lett. 20, 427.
Sagdeev, R. Z. , and A. A. Galeev, 1970, Sov. Phys. —Dokl. 14, Tamm, I. E., 1959, in Plasma Physics and the Pxoblem of
1198. Contxolled Tkexmonucleax Pusion (Pergamon, New York),
Shafranov, V. D. , 1958, Sov.. Phys. —JETP 8, 545. Vol. 1, p. 35.
Shafranov, V. D. , 1965, At. Energ. 19, 120. Taylor, J. B., 1961, Phys. Fluids 4, 1142.
Shafranov, V. D. , 1966, in Reviewers of Plasma Physics, edited Trubnikov, B. A. , 1965, in geviezes of Plasma Physics,
by M. A. Leontovich (Consultants Bureau, New York), Vol. edited by M. A. Leontovich (Consultants Bureau, New York),
2, p. 103. Vol. 1, p. 105.
Sigmar, D. J., 1972, talk given at the Sherwood Theory Meet- Tsang, K. T. , and E. A. Frieman, 1975, to be published in
ing, Los Alamos, March 23 (unpublished). Phys. Fluids.
Sigmar, D. J., J. F. Clarke, R. V. Neidigh, and K. L. Vander Ware, A. A. , 1970, Phys. Rev. Lett. 25, 916.
Sluis, 1974, Phys. Rev. Lett. 33, 1376. Ware, A. A. , 1973, Nucl. Fusion t3, 793.
Spitzer, L. , Jr. , 1952, Astrophys. J. 116, 299. Wimmel, H. K. , 1970, Nucl. Fusion 10, 117.