henry2004y / TestParticle.jl

Test particle tracing in electromagnetic field
https://henry2004y.github.io/TestParticle.jl/dev/
MIT License
15 stars 3 forks source link

Curvature and grad-B drifts #61

Open henry2004y opened 1 year ago

henry2004y commented 1 year ago

From the fluid theory, the curvature drift exists while the grad-B drift does not. Are we able to demonstrate this with our test particle model?

TCLiuu commented 1 year ago

It is very interesting, I am not familiar with the fluid theory. But according to the kinetic theory, the effect of grad-B drift has an influence on the distribution function and we call it kinetic effect (include curvature drift, grad-B drift and FLR).

I have found a thread about this question.

TCLiuu commented 1 year ago

I found two papers relating to this question. Paper 1: calculating the drift flow by gyrokinetic moments. Paper 2: also gyrokinetic model.

henry2004y commented 1 year ago

It is very interesting, I am not familiar with the fluid theory. But according to the kinetic theory, the effect of grad-B drift has an influence on the distribution function and we call it kinetic effect (include curvature drift, grad-B drift and FLR).

I have found a thread about this question.

The linked thread is slightly biased, but I get his point. In the fluid theory, which I quoted from Chen's book, there is a strong assumption

If f(v) remains Maxwellian in a nonuniform B field, and there is no density gradient, then the net momentum carried into any fixed fluid element is zero.

Obviously even if only $\nabla B$ exists, an initially Maxwellian distributed plasma will be distorted since

$$ \mathbf{v}{\nabla B} = \pm \frac{1}{2}v\perp rL\frac{\nabla B\times\mathbf{B}}{B^2} \propto v\perp^2 $$

However, I think this will NOT generate anisotropies in the distribution since it is symmetric with $|v\perp|$? If there are still the same number of particles with the same $|v\perp|$, then Chen's argument is still valid?


I haven't read those two papers. Do you have a quick scan of what their conclusion is about our question @TCLiuu ?

TCLiuu commented 1 year ago

It is very interesting, I am not familiar with the fluid theory. But according to the kinetic theory, the effect of grad-B drift has an influence on the distribution function and we call it kinetic effect (include curvature drift, grad-B drift and FLR). I have found a thread about this question.

The linked thread is slightly biased, but I get his point. In the fluid theory, which I quoted from Chen's book, there is a strong assumption

If f(v) remains Maxwellian in a nonuniform B field, and there is no density gradient, then the net momentum carried into any fixed fluid element is zero.

Obviously even if only ∇B exists, an initially Maxwellian distributed plasma will be distorted since

v∇B=±12v⊥rL∇B×BB2∝v⊥2

However, I think this will NOT generate anisotropies in the distribution since it is symmetric with |v⊥|? If there are still the same number of particles with the same |v⊥|, then Chen's argument is still valid?

I haven't read those two papers. Do you have a quick scan of what their conclusion is about our question @TCLiuu ?

To summarize, from the viewpoint of gyrokinetics, the fluid perpendicular current in isotropic plasmas, $c\mathbf{b}\times\nabla p{\perp}/B$, contains both the diamagnetic current and part of the current caused by the gradient drift. The other part of the current caused by the gradient drift cancels the current generated by the curvature drift for isotropic distribution functions. From the viewpoint of gyrokinetics, the diamagnetic current is $c\nabla\times (p{\perp}\mathbf{b}/B)$ actually.

TCLiuu commented 1 year ago

We can calculate the perpendicular current by the method of moment integral. It can be write as: $$\mathbf{J}_\perp =\mathbf{J}_d+\mathbf{J}_M= \int e(\dot{\mathbf{X}}+\dot{\mathbf{\rho}})\delta(\mathbf{X}+\mathbf{\rho}-\mathbf{r})F(\mathbf{Z})\mathrm{d}^6\mathbf{Z}$$ $$\mathbf{J}_M= \int e\dot{\mathbf{\rho}}\delta(\mathbf{X}+\mathbf{\rho}-\mathbf{r})F(\mathbf{Z})\mathrm{d}^6\mathbf{Z}$$ $$\mathbf{J}_d= \int e\dot{\mathbf{X}}\delta(\mathbf{X}+\mathbf{\rho}-\mathbf{r})F(\mathbf{Z})\mathrm{d}^6\mathbf{Z}$$

We can separate naturally the perpendicular current to two terms (the so-called diamagnetic drift term and the guiding center drift term). The part of diamagnetic drift term will cancel the effect of the grad-B drift and reduce to the result of the fluid theory.

TCLiuu commented 1 year ago

Recently, I have been preparing for this part of the content and therefore studied it again. I think ffchen's statement here is not very rigorous.

If considering anisotropic plasma, then the fluid velocity in the direction perpendicular to the magnetic field should be written as

image

where

image

And the actual requirements that need to be met is ignoring the parallel flow, that is $u_\parallel$≈0. If this item is not ignored, there will be a corresponding curvature drift.

After ignoring the parallel flow, according to my previous post, magnetic field gradient drift and curvature drift still exist, but all the gradient drift parts and some of the curvature drift parts and diamagnetic drift parts are cancelled out. In isotropic plasma, only the remaining diamagnetic drift part can be seen, while in anisotropic plasma, we can still see curvature drift. This does not mean that there is no gradient drift. From the perspective of kinetic theory, it is clear that diamagnetic drift is a result of FLR effect. As for FFChen's discussion on distribution, I don't think it's rigorous enough but it's acceptable in introductory courses.

I strongly recommend that you take a look at paper2 mentioned above, which clearly demonstrates the entire process of this issue from gyrokinetic theory to MHD.

I am curious about your current understanding of this issue.

henry2004y commented 1 year ago

My current understanding of this issue is based on Section 3.6 in Paul Bellan's Fundamentals of Plasma Physics: Relations of Drift Equations to the Double Adiabatic MHD equations. He showed that by taking the diamagnetic current, grad B current, curvature current, and polarization current into account, the drift equations for phenomena with characteristic frequencies $\omega << \omega_{ci}$ and the double adiabatic MHD equations are equivalent descriptions of plasma dynamics.