## Brazilian Journal of Physics

##
*Print version* ISSN 0103-9733

*On-line version* ISSN 1678-4448

### Braz. J. Phys. vol.31 no.1 São Paulo Mar. 2001

#### http://dx.doi.org/10.1590/S0103-97332001000100006

**Effect of plasma subsonic toroidal flows induced by Alfvén waves on transport processes in the edge of elongated tokamaks **

V. S. Tsypin, J. H. F. Severo, I. C. Nascimento,

R. M. O. Galvão, and Yu. K. Kuznetsov

*Institute of Physics, University of São Paulo,*

*Rua do Matão, Travessa R, 187, 05508-900, São Paulo, Brasil*

**Received on 4 May, 2000 **

There is a renewed interest in using Alfvén waves (AW) in tokamak plasmas. Previously, AW were actively explored mostly for current drive and plasma heating in tokamaks. Presently, the possibility of the anomalous and neoclassical transport suppression by AW in tokamak plasmas is being vividly discussed. AW can also induce poloidal and toroidal plasma rotation. Toroidal plasma rotation can reach the subsonic level. These flows can substantially affect neoclassical transport both in collisional and weakly collisional plasmas. In this paper, the effect of plasma subsonic toroidal flows induced by Alfvén waves on transport processes in the edge of elongated tokamak is investigated. The dependence of poloidal plasma rotation and ion heat conductivity on the elongation parameter and the ratio of induced toroidal velocity to the sonic speed are analytically obtained.

**I Introduction**

It is well-known that Alfvén waves are one of the basic methods of current drive (AWCD) and plasma heating (AWPH) in tokamaks.^{1-5} These methods were clearly demonstrated in numerous experiments^{6-8} and will be anew applied in forthcoming experiments in TCABR (Tokamak Chauffage Alfvén Brasil) tokamak.^{9,10} The comparatively recent applications of AW in tokamaks are creation of transport barriers and suppression of anomalous and neoclassical transport by these waves.^{11-14} The main idea is to induce a sheared radial electric field that can take apart the turbulent plasma eddies and thus suppress anomalous transport.^{15-17} Suppression of neoclassical transport in weakly collisional tokamak plasmas is connected with ion banana orbit squeezing by the sheared radial electric field.^{18,19} Alfvén waves can be used to efficiently induce the strongly sheared radial electric field^{11-14} thus achieving suppression of anomalous and neoclassical transport. The first experimental evidence that AW can be useful to create plasma rotation, and, consequently, the radial electric field, was obtained in the Phaedrus-T tokamak.^{20}

Besides these very attractive topics, there is also the additional physical mechanism where AW can play a substantial role, i.e., the ion neoclassical transport in the tokamak edge (collisional) plasma. The importance of this mechanism was emphasized in the papers by Rogister,^{21-23} where the experimental results in modern tokamaks were carefully analyzed. Rogister considers the transport processes in tokamak edge plasmas to be in the H-mode (the so-called high regime), when the characteristic length of plasma macroscopic parameters *L _{n}* is of the order of the ion Larmor radius, r

_{iq }=

*cv*/

_{Ti}M_{i}*e*

_{i}B_{qp}, calculated on the poloidal magnetic field

*B*

_{qp}, i.e.,

*L*~ r

_{n}_{iq}. Here

*e*and

_{i}*M*are the ion charge and mass, respectively,

_{i}*B*

_{qp}is the ''physical'' component of the magnetic field, to differ it from covariant or contravariant components that appear below, and

*v*= is the ion thermal velocity. However, there are some relevant problems that need to be investigated in collisional tokamak plasmas, in the L-mode (the low regime), when the condition

_{Ti }*L*> r

_{n}_{iq}is fulfilled.

The ion neoclassical transport in tokamak collisional plasmas^{24} was originally investigated by Shafranov. He used the well-known approach by Pfirsch and Schlüter to calculate the ion heat conductivity in tokamak plasmas with circular magnetic surfaces. These results were later generalized for arbitrary geometry of the magnetic field.^{25} The next step in this direction was taken by Mikhailovskii and Tsypin,^{26} who considered the effect of sonic toroidal flows on transport processes in collisional plasmas. The very interesting effect of changing the sign of the ion poloidal velocity at a specific value of the ratio of toroidal velocity *U _{i}*z to sound velocity = (

*T*+

_{e }*T*) /

_{i}*M*, and the importance of the collisional parameter

_{i}*b*=

*q*

^{2}

*R*

^{2}/l

_{i}

^{2}were brought up in this paper. Here,

*q*is the inverse rotational transform,

*R*is the tokamak major radius, l

*=*

_{i }*v*/n

_{Ti}_{i}is the ion mean free path,

*T*and

_{e}*T*are the electron and ion temperature, respectively, and n

_{i}_{i}is the ion-ion collision frequency.

^{27}This study was later continued to determine also the ion heat conductivity taking into account of large values of the collisional parameter

*b*, 1 <

*b*<

*M*/

_{i}*M*,

_{e}^{28,29}external forces,

^{30}and elongation.

^{31}In this paper we carry on this investigation studying the effect of plasma subsonic toroidal flows induced by Alfvén waves on poloidal plasma rotation and ion heat conductivity in collisional plasmas of elongated tokamaks, with smooth profiles of macroscopic plasma parameters,

*L*> r

_{n}_{iq}.

**II Magnetic field and metric **

To solve our problem, we need to choose the coordinate system to carry out the calculations for a toroidal plasma column with elliptic cross-section. As the calculations of the metrics have been repeatedly presented in many previous publications, we here show only main steps to obtain it. Initially, we take the coordinates r¢, q^{¢¢} and z to be attached to the geometrical center of the magnetic surfaces, where r¢ is the radial coordinate in the tokamak cross-section and q^{¢¢} and z are the poloidal and toroidal angles, respectively. The length element in these coordinates is

The next step is the ''encircling'' of the magnetic surfaces by means of the transformation

where r = , h = ln(*l*_{2}/*l*_{1}), and *l*_{2} and *l*_{1} are semimajor and semiminor axes, respectively, of the magnetic surface. Then, we choose new coordinates to be attached to the geometric center of the tokamak cross-section, and straighten the magnetic lines of force:

where D(*r*) is the shift of the geometric center of the magnetic surface, determined from equilibrium conditions.

The straightening parameter d can be found from the condition that the radial component of the current density vanishes, *j ^{r}*=0, by taking into account the expression for the magnetic field with the straight lines of force,

(here c and f are the poloidal and toroidal magnetic fluxes, respectively, and *g* is the metric tensor determinant), and after finding the metric tensor components *g _{ik}*, after the substitution of Eqs. (2) and (3) into Eq. (1),

As a result, the parameter d is equal to

Then, we find the metric covariant tensor components and their determinant

where the parameter e^{*} is e^{*}=eexp(-h/2), and e = *r*/*R* is the tokamak inverse aspect ratio. The angle dependence of the magnetic field can be obtained from the expression *B*_{0 }= ; we have

where

and *q *= f¢/c¢ is the inverse rotational transform defined in general geometry.

We also need a component of the Cristoffel's symbol

namely,

The non-orthogonal coordinates are *x ^{i}*=(

*r*, q, z).

**III Starting equations**

As we mentioned in the Introduction, our main goal is to calculate the surface averaged ion heat flux G_{Ti} and the ion poloidal rotation velocity *U _{i}*q as functions of the ratio a = / (the squared Mach number) and the elongation parameter

*A*. The ion heat flux is defined by the expression

where

and the heat flux radial contravariant component *q _{i}^{r}* is equal to

Here, *p _{i}* is the ion isotropic pressure, w

*=*

_{ci }*e*/

_{i}B*cM*is the ion cyclotron frequency,

_{i}*g*

^{11}is the contravariant component of the metric tensor, and

*h*

^{z }=

*B*

^{z}/

*B*is the z-contravariant component of the unit vector

**h = B**/

*B*. The conventional relation between vector covariant and contravariant components

*B*

_{a }=

*g*

_{ab}

*B*

^{b}is also used below. The physical poloidal component of velocity

*U*q is defined by the expression

_{i}*U*q =

_{i}*rV*

_{i}^{q}, where

*V*

_{i}^{q}is the q-contravariant component of the ion macroscopic velocity

**V**

_{i}.

From Eqs. (15)-(17) we find the surface averaged ion heat flux,

Here, *n*_{0} is the particle density.

Further, to find the poloidal rotation velocity *U _{i}*q we employ the one fluid momentum equation

where *p *= *p _{i }*+

*p*is the plasma pressure,

_{e}**p**is the ion viscosity,

^{27}in which we take into account only parallel viscosity p

_{||}, therefore,

and **F**^{h} is the radio frequency force affecting electrons and ions. Using the ambipolarity condition,

which can be obtained from the current continuity equation · **j** » **0** averaged over the magnetic surface, one obtains from Eq. (19)

Here, m_{iz} is the coefficient describing the ion toroidal or anomalous viscosity or friction with neutrals, and and are ''physical'' poloidal and toroidal components of the radio frequency force **F**^{h}, respectively.

The plasma density poloidal perturbations _{0} in Eq. (18) can be found from the momentum equation (19), taking into account the ion inertial terms connected with the subsonic toroidal flow, *U _{i}*z <

*c*, and the ion and electron temperature perturbations,

_{s}Using Eq. (24), the expression for the Cristoffel's symbol Eq. (14), and the parallel component of Eq. (19), one obtains from Eq. (22),

Thus, our basic equations for the problem under study are Eqs. (18) and (25).

**IV Poloidally perturbed quantities **

We see from Eqs.(18) and (25) that we should find the perturbed ion and electron temperatures and ion parallel viscosity. To find the particle temperatures, we proceed from the heat transport equations,^{27}

where a = *i*, *e*, the heat exchange *Q*_{a} between ions and electrons is equal to

and particle heat fluxes are

Then, we obtain equations for perturbed quantities from Eqs. (26)-(28)

The terms with the poloidal velocities *U _{i}*q and

*U*q in Eqs. (29) and (30) can be obtained using the quasistationary continuity equations

_{e}and the frozen-in condition

As a result, we have

Using equation (33) and solving the system of Eqs.(29) and (30), we find the perturbed particle temperatures

where *d*(*b*) =1 + 2.2*b*, *d*_{1}(*b*) = 1 + 0.54*b*, and *b *= *q*^{2}*R*^{2}/l_{i}^{2}. We have also used

with *T _{e}*0 »

*T*0 =

_{i}*T*

_{0}. The expression for

*U*q in Eq. (36) can be obtained from momentum equations for electrons and ions.

_{e}^{27}

We also need to calculate the parallel viscosity _{||} in Eqs. (24) and (25). Note, that the Braginskii's approach^{27} for viscosity (the Navier-Stokes type of the viscosity) is not relevant for problems to be under study in this paper. Therefore, we take the parallel viscosity p_{||} in the form derived in Refs. 26, 32, and 33. This type of viscosity is the so-called Burnett type of the viscosity.^{34} Besides the spatial derivatives of the ion velocity, it also includes the derivatives of the ion thermal fluxes **q**_{i}. These additional terms are similar to the thermal force terms in the ion-electron friction in the momentum equation.^{27} Thus, we have,^{26,32,33}

where

Using Eqs.(28) and (33), we simplify Eqs. (38) and (39)

Finally, we obtain the expression for the parallel viscosity Eqs. (37), (40), and (41), employing Eqs. (8), (11), and (24),

This expression gives us the possibility to find the ion ambipolar poloidal velocity *U _{i}*q.

**V Ion fluxes and their analyses **

Now, we can derive the final expressions for plasma rotation velocity and the radial ion heat flux and analyze them. The poloidal velocity *U _{i}*q can be obtained from Eqs. (25), (35), and (42),

where

To plot functions *G _{u}*1(a,

*A*),

*G*2(a,

_{u}*A*),

*G*3(a,

_{u}*A*) in Fig. 1, we choose the parameter

*b*= 50. One can observe that function

*G*1(a,

_{u}*A*) (figure 1), as in the case of an axially-symmetric tokamak,

^{30}changes sign at a

_{0}» 2

*d*(

*b*)/

*b*. This quantity coincides with our previously obtained results.

^{26,28-31}It is a very interesting and important result that inertial forces, affecting ions, enforce ion poloidal rotation velocity

*U*q to change its sign in the absence of radio frequency waves in a plasma. This effect, which was originally obtained in Ref. 26, is also relevant for weakly collisional plasmas but sign changing takes place at the Mach number of the order unity.

_{i}^{35}The maximum of function

*G*1(a,

_{u}*A*) occurs at

*A*» 0 and is

*G*1(a, 0) » 2.4. Increasing the parameter

_{u}*A*, the maximum of the quantity

*G*1(a,

_{u}*A*) decreases and shifts to supersonic values, a > 1. Nevertheless, only the case a < 1 is of the practical interest.

The function *G _{u}*2(a,

*A*) characterizes the role of the parallel current and diamagnetic drift in the ion poloidal rotation velocity. This function increases with the parameter a, but stays at a small level,

*G*2(a,

_{u}*A*) << 1 for a < 1.

The function *G _{u}*3(a,

*A*) describes the role of external forces in inducing the ion poloidal rotation. This function decreases with the growth of the plasma toroidal rotation, and, correspondingly, so does the role of external forces on the ion poloidal rotation. The induced radial electric field, which is currently supposed to be the dominant quantity for transport barrier creation in tokamaks,

^{12-17}can be found using the expression which is defined via the ion velocity from the ion motion equation,

^{27}

Using Eqs. (11), (24), (34), (42), and (43), we derive the magnetic surface averaged radial ion heat flux from Eq. (18). It is convenient to write the final equation in the Shafranov form,^{24,25}

where

In a supposition that the collisional parameter *b* satisfies the inequality 1 < *b* < , one obtains the well-known result for the radial ion heat flux.^{24,25}

The coefficients *G _{T}*1(a,

*A*) and

*G*3(a,

_{T}*A*) are plotted in Fig. 2. The quantity

*G*2(a,

_{T}*A*) is small under the condition

*b*£

*M*/

_{i}*M*, therefore, it is not presented in Fig. 2. The functions

_{e}*G*1(a,

_{T}*A*) and

*G*3(a,

_{T}*A*) characterize the neoclassical contribution and the influence of external forces, respectively. The function

*G*1(a,

_{T}*A*) is a growing function of the elongation

*A*and the squared Mach number a. The role of external forces increases with the Mach number, nevertheless, less intensively for larger elongations [see the function

*G*3(a,

_{T}*A*) in Fig. 2].

**VI Conclusion**

We have studied the effect of plasma subsonic toroidal flows induced by Alfvén waves on poloidal rotation and ion heat conductivity of collisional plasmas of elongated tokamaks with smooth profiles of macroscopic plasma parameters. This study is important as ion neoclassical fluxes are sometimes in a satisfactory coincidence with the experimental results,^{36} as opposed to the electron

Figure. 1 The dependence of the functions G1(a, _{u}A), G2(a, _{u}A), and G3(a, _{u}A) on a, for different magnitudes of the parameter A : A_{1 }= 0, A_{2 }= 1, A_{3 }= 2.. |

Figure. 2 The dependence of the functions G1(a, _{T}A) and G3(a, _{T}A) on a, for different magnitudes of the parameter A : A_{1 }= 0, A_{2 }= 1, A_{3 }= 2.. |

anomalous transport. This paper is a generalization of our previous results for the case of elongated tokamaks.^{26,28-31} It is shown that a tokamak elongation can substantially affect ion neoclassical fluxes. These fluxes are presented in the paper as functions of three quantities to be proportional to 1) the conventional ion temperature gradient, 2) the parallel current and diamagnetic drift, and 3) radio frequency forces. The proper coefficients are plotted on the Fig. 1 and Fig. 2. Once more, the importance of the collisional parameter *b* is emphasized in the paper (see also Refs. 21-23, 26, 28-31). It is also demonstrated that the very interesting effect of changing the sign of the ion poloidal rotation velocity in collisional tokamak plasmas, which depends on the square of the Mach number a, obtained previously,^{26} also takes place in elongated tokamaks. The poloidal velocity has a maximum as a function of the Much number. This sign changing and maximum of the ion poloidal rotation is the consequence of taking into account inertial forces affecting ions. It is also shown that the radio frequency forces can in some cases determine the ion neoclassical fluxes in elongated tokamaks.

**Acknowledgments**

This work was supported by the Research Support Foundation of the State of São Paulo (FAPESP), National Council of Scientific and Technological Development (CNPq), and Excellence Research Programs (PRONEX) RMOG 50/70 grant from the Ministry of Science and Technology, Brazil.

**References**

[1] Ya. I. Kolesnichenko, V. V. Parail, and G. V. Pereverzev, in *Reviews of Plasma Physics,* edited by B. B. Kadomtsev, (Consultants Bureau, New York, 1990) Vol. 17, p. 3. [ Links ]

[2] A. G. Elfimov, A. G. Kirov and V. P. Sidorov, in *High-Frequency Plasma Heating*, edited by A. G. Litvak (American Institute of Physics Translation series, New-York, 1992), p.239. [ Links ]

[3] R. Klima, Czech. J. Phys. B **18**, 1280 (1968). [ Links ]

[4] N. J. Fish and F. F. Karney, Phys. Fluids **24**, 27 (1981). [ Links ]

[5] D. W. Faulconer, Fusion Technology **33** (2T) 219 (1998). [ Links ]

[6] R.A.Demirkhanov, A.G.Kirov, L.F.Ruchko, and A.V.Sukachev, Soviet Phys. JETP Lett. **33**, 28 (1981). [ Links ]

[7] G. A. Collins, F. Hofmann, B. Joye, R. Keller, A. Lietti, J. B. Lister, and A. Pochelon, Phys. Fluids, **29**, 2260 (1986). [ Links ]

[8] N. Hershkovitz, P. Moroz, P. Probert, E. Y. Wang, T. Intrator, D. Diebold, R. Breun, D. Brouchous, M. Doczy, D. Edgel, A. Elfimov, R. Fonk, M. Kishinevsky, C. Litwin, P. Nonn, S. Regan, J. Sorensen, J. Tataronis, M. Vucovich, X. Wang, and S. Wukitch, in *Proceedings of the 15th International Atomic Energy Agency Conference on Plasma Physics and Controlled Fusion*, Seville, 1994 (International Atomic Energy Agency, Vienna, 1995), Vol. 2, p. 245. [ Links ]

[9] I. C. Nascimento, R. M. O.Galvão and TCABR team, in Proc. of X Intern. Conf. on Plasma Physics. Contr. Papers. Brazil, Foz do Iguaçu, 31.10.94-04.11.94, (AIP, Woodbury, New York, 1996), Vol. 1, p. 69. [ Links ]

[10] R. M. O. Galvão, A. G. Elfimov, G. Amarante-Segundo, V. S. Tsypin, L. F. Ruchko, I. C. Nascimento, M. Tendler. Plasma Phys. Contr. Fusion **41**, A487 (1999). [ Links ]

[11] G. G. Craddock and P. H. Diamond, Phys. Rev. Lett. **67**, 1535 (1991). [ Links ]

[12] V. S. Tsypin, A. G. Elfimov, M. Tendler, C. A. de Azevedo, and A. S. de Assis, Phys. Plasmas **5**, 7 (1998). [ Links ]

[13] V. S. Tsypin, R. M. O. Galvão, I. C. Nascimento, A. G. Elfimov, M. Tendler, C. A. de Azevedo, and A. S. de Assis, Phys. Rev. Lett. **81**, 3403 (1998). [ Links ]

[14] V. S. Tsypin, I. C. Nascimento, R. M. O. Galvão, A. G. Elfimov, G. S. Amarante Segundo, and M. Tendler Phys. Plasmas **6**, 3548 (1999). [ Links ]

[15] H. Biglary, P. H. Diamond, and P. W. Terry, Phys. Fluids B **2, **1 (1990). [ Links ]

[16] T. S. Hahm and K. H. Burrel. Phys. Plasmas **2**, 1648 (1995) [ Links ]

[17] K. H. Burrell, Phys. Plasmas **4**, 1499 (1997). [ Links ]

[18] H. L. Berk, A. A. Galeev. Phys. Fluids **10**, 44 (1967). [ Links ]

[19] K. C. Shaing and R. D. Hazeltine, Phys. Fluids B **4**, 2547 (1992). [ Links ]

[20] S. Wukitch, C. Litwin, M. Harper, R. Parker, and N. Hershkovitz, Phys. Rev. Lett. **77**, 294 (1996). [ Links ]

[21] A. Rogister, Phys. Plasmas **1, **619 (1994). [ Links ]

[22] A. Rogister, Phys. Rev. Lett. **81, **3663 (1998). [ Links ]

[23] A. Rogister, Phys. Plasmas **6, **200 (1999). [ Links ]

[24] V. D. Shafranov, Sov. Atom. Energy **19**, 1008 (1965). [ Links ]

[25] V. D. Shafranov, J. W. Connor, and C. J. Watson, Sov. J. Plasma Phys. **2**, 99 (1976). [ Links ]

[26] A. B. Mikhailovskii and V. S. Tsypin, Sov. Phys.-JETP **56**, 75 (1982). [ Links ]

[27] S. I. Braginskii, in *Reviews of Plasma Physics*, edited by M. A. Leontovich (Consultants Bureau, New York, 1965), Vol. 1, p. 205. [ Links ]

[28] V. S. Tsypin, Sov. J. Plasma Phys. **11**, 661 (1985). [ Links ]

[29] V. S. Tsypin, Sov. J. Plasma Phys. **14**, 466 (1988). [ Links ]

[30] V. S.Tsypin, D. Kh.Morozov, J. J. E.Herrera, J. J. Martinell, M. Tendler, I. F. Potapenko, A. S. de Assis, and C. A. de Azevedo. Plasma Phys. Control. Fusion **39**, 1681 (1997). [ Links ]

[31] V. S. Tsypin, C. A. de Azevedo, and A. S. de Assis, Physics Letters A. **219**, 282 (1996). [ Links ]

[32] A. B. Mikhailovskii and V. S. Tsypin, Plasma Physics **13**, 785 (1971). [ Links ]

[33] A. B. Mikhailovskii and V. S. Tsypin, Beiträge Plasma Physik 24, 335 (1984). [ Links ]

[34] D. Burnett, Proc. Lond. Math. Soc. **39** 385; **40** 382 (1935). [ Links ]

[35] S. K. Wong and K. H. Burrel, Phys. Fluids **25 **1863 (1982). [ Links ]

[36] F. Wagner and U. Strouth, Plasma Phys. Contr. Fusion **34** 1803 (1992). [ Links ]