Drift Equations

The drift equations take the form:

$$\frac{\partial n}{\partial t} + \nabla\cdot(n\,{\bf V}_e) =0,$$ (4.219)

$$m_e\, n\,\frac{\partial {\bf V}_e}{\partial t} + m_e\, n\,({\bf V... ...cdot\nabla){\bf V}_e+ [\delta^{\,-2}]\,\nabla p_e+ [\zeta^{\,-1}]\, \nabla\cdot \mbox{\boldmath$\pi$} _e$$ (4.220)

$$+ [\delta^{\,-2}]\,e \,n\, ({\bf E} + {\bf V}_e\times {\bf B}) =[\delta^{\,-2}\,\zeta]\, {\bf F}_U + [\delta^{\,-2}]\, {\bf F}_T,$$

$$\frac{3}{2}\frac{\partial p_e}{\partial t} + \frac{3}{2} \,({\bf V}_e\cdot\nabla)\, p_e + \frac{5}{2}\,p_e\,\nabla\cdot{\bf V}_e$$

$$+ [\zeta^{\,-1}]\,\nabla\cdot{\bf q}_{Te} +\nabla\cdot{\bf q}_{Ue} =-[\delta^{\,-2}\,\zeta\,\mu^2]\, W_i$$

$$\phantom{=}+[\zeta]\,W_U + W_T,$$ (4.221)

and

$$\frac{\partial n}{\partial t} + \nabla\cdot(n\,{\bf V}_i) =0,$$ (4.222)

$$m_i \,n\,\frac{\partial {\bf V}_i}{\partial t} + m_i \,n\,({\bf V... ...ot\nabla) {\bf V}_i + [\delta^{\,-2}]\,\nabla p_i + [\zeta^{\,-1}]\,\nabla\cdot \mbox{\boldmath$\pi$} _i$$ (4.223)

$$- [\delta^{\,-2}]\,e\, n\, ({\bf E} + {\bf V}_i\times {\bf B}) = -[\delta^{\,-2}\,\zeta]\, {\bf F}_U - [\delta^{\,-2}]\,{\bf F}_T,$$

$$\frac{3}{2}\frac{\partial p_i}{\partial t} + \frac{3}{2} \,({\bf V}_i\cdot\nabla) \,p_i + \frac{5}{2}\,p_i\,\nabla\cdot{\bf V}_i$$

$$+ [\zeta^{\,-1}]\,\nabla\cdot{\bf q}_i = [\delta^{\,-2}\,\zeta\,\mu^2]\,W_i.$$ (4.224)

In the drift limit, the motions of the electron and ion fluids are sufficiently different that there is little to be gained in rewriting the drift equations in terms of the center-of-mass velocity and the plasma current. Instead, let us consider the ${\cal O}(\delta^{\,-2})$ components of Equations (4.220) and (4.223):

$${\bf E} + {\bf V}_e\times{\bf B} \simeq - \frac{\nabla p_e}{e\,n} -\frac{ 0.71\,\nabla_\parallel T_e}{e},$$ (4.225)

$${\bf E} + {\bf V}_i\times {\bf B} \simeq +\frac{\nabla p_i}{e\,n} - \frac{ 0.71\,\nabla_\parallel T_e}{e}.$$ (4.226)

In the previous equations, we have neglected all ${\cal O}(\zeta)$ terms for the sake of simplicity. Equations (4.225)-(4.226) can be inverted to give

$${\bf V}_{\perp\,e} \simeq {\bf V}_E + {\bf V}_{\ast\,e},$$ (4.227)

$${\bf V}_{\perp\,i} \simeq {\bf V}_E + {\bf V}_{\ast\,i}.$$ (4.228)

Here, ${\bf V}_E\equiv {\bf E}\times{\bf B}/B^{\,2}$ is the ${\bf E}\times{\bf B}$ velocity, whereas

$${\bf V}_{\ast\,e} \equiv \frac{\nabla p_e \times{\bf B}} {e\,n\,B^{\,2}},$$ (4.229)

and

$${\bf V}_{\ast\,i} \equiv -\frac{\nabla p_i \times{\bf B}} {e\,n\,B^{\,2}},$$ (4.230)

are termed the electron diamagnetic velocity and the ion diamagnetic velocity, respectively.

According to Equations (4.227)-(4.228), in the drift approximation the velocity of the electron fluid perpendicular to the magnetic field is the sum of the ${\bf E}\times{\bf B}$ velocity and the electron diamagnetic velocity. A similar statement can be made for the ion fluid. By contrast, in the MHD approximation the perpendicular velocities of the two fluids consist of the ${\bf E}\times{\bf B}$ velocity alone, and are, therefore, identical to lowest order. The main difference between the two orderings lies in the assumed magnitude of the electric field. In the MHD limit

$$\frac{E}{B} \sim v_t,$$ (4.231)

whereas in the drift limit

$$\frac{E}{B} \sim \delta\,v_t\sim v_d.$$ (4.232)

Thus, the MHD ordering can be regarded as a strong electric field ordering, whereas the drift ordering corresponds to a weak electric field ordering.

The diamagnetic velocities are so named because the diamagnetic current,

$${\bf j}_\ast \equiv - e\,n\,({\bf V}_{\ast\,e}-{\bf V}_{\ast\,i}) =-\frac{\nabla p\times{\bf B}}{B^{\,2}},$$ (4.233)

generally acts to reduce the magnitude of the magnetic field inside the plasma.

The electron diamagnetic velocity can be written

$${\bf V}_{\ast\,e} = \frac{T_e\,\nabla n\times{\bf b}}{e\,n\,B} + \frac{ \nabla T_e\times{\bf b}}{e\,B}.$$ (4.234)

In order to account for this velocity, let us consider a simplified case in which the electron temperature is uniform, there is a uniform density gradient running along the $x$ -direction, and the magnetic field is parallel to the $z$ -axis. (See Figure 4.3.) The electrons gyrate in the $x$ -$y$ plane in circles of radius $\rho_e\sim v_e/\vert{\mit\Omega}_e\vert$ . At a given point, coordinate $x_0$ , say, on the $x$ -axis, the electrons that come from the right and the left have traversed distances of approximate magnitude $\rho_e$ . Thus, the electrons from the right originate from regions where the particle density is approximately $\rho_e\,\partial n/\partial x$ greater than the regions from which the electrons from the left originate. It follows that the $y$ -directed particle flux is unbalanced, with slightly more particles moving in the $-y$ -direction than in the $+y$ -direction. Thus, there is a net particle flux in the $-y$ -direction: that is, in the direction of $\nabla n \times{\bf b}$ . The magnitude of this flux is

$$n\,V_{\ast\,e} \sim \rho_e\,\frac{\partial n}{\partial x}\,v_e \sim \frac{T_e}{e\,B} \frac{\partial n}{\partial x}.$$ (4.235)

There is no unbalanced particle flux in the $x$ -direction, because the $x$ -directed fluxes are associated with electrons that originate from regions where $x=x_0$ . We have now accounted for the first term on the right-hand side of Equation (4.234). We can account for the second term using similar arguments. The ion diamagnetic velocity is similar in magnitude to the electron diamagnetic velocity, but is oppositely directed, because ions gyrate in the opposite direction to electrons.

Figure 4.3: Origin of the diamagnetic velocity in a magnetized plasma.

The most curious aspect of diamagnetic flows is that they represent fluid flows for which there is no corresponding motion of the particle guiding centers. Nevertheless, the diamagnetic velocities are real fluid velocities, and the associated diamagnetic current is a real current. For instance, the diamagnetic current contributes to force balance inside the plasma, and also gives rise to ohmic heating.