1887
banner image
No data available.
Please log in to see this content.
You have no subscription access to this content.
No metrics data to plot.
The attempt to load metrics for this article has failed.
The attempt to plot a graph for these metrics has failed.
Linearized model Fokker–Planck collision operators for gyrokinetic simulations. II. Numerical implementation and tests
Rent:
Rent this article for
USD
10.1063/1.3155085
/content/aip/journal/pop/16/7/10.1063/1.3155085
http://aip.metastore.ingenta.com/content/aip/journal/pop/16/7/10.1063/1.3155085

Figures

Image of FIG. 1.
FIG. 1.

(Left) Solid line indicates the scaling of the leading order error, averaged over all grid points, of the conservative finite difference scheme for a Gauss–Legendre grid (the grid used in GS2). The slope of the dotted line corresponds to a first order scheme. (Right) Factor by which the conservative finite difference scheme of Eq. (43) amplifies the true collision operator amplitude at the boundaries of the Gauss–Legendre grid.

Image of FIG. 2.
FIG. 2.

Plots showing evolution of the perturbed local density, parallel momentum, and energy over fifty collision times. Without the conserving terms (9)–(11), both parallel momentum and energy decay significantly over a few collision times (long dashed lines). Inclusion of conserving terms with the conservative scheme detailed in Sec. III leads to exact moment conservation (solid lines). Use of a nonconservative scheme leads to inexact conservation that depends on grid spacing (short dashed lines).

Image of FIG. 3.
FIG. 3.

Plot of the evolution of entropy generation for the homogeneous plasma slab over 20 collision times. Our initial distribution in velocity space is random noise, and we use a grid with 16 pitch angles and 8 energies. The entropy generation rate is always non-negative and approaches zero in the long-time limit.

Image of FIG. 4.
FIG. 4.

Evolution of for the electromagnetic plasma slab with , , and . Inclusion of the ion drag term in the electron-ion collision operator leads to the theoretically predicted damping rate for the parallel current given in Eq. (50). Without the ion drag term, the parallel current decays past zero (at ) and converges to a negative value as the electron flow damps to zero.

Image of FIG. 5.
FIG. 5.

Evolution of perturbed parallel flow for the electromagnetic plasma slab with , , and . Without inclusion of the ion drag term in Eq. (12), the electron flow is erroneously damped to zero (instead of to the ion flow).

Image of FIG. 6.
FIG. 6.

Damping rate of the slow mode for a range of collisionalities spanning the collisionless to strongly collisional regimes. Dashed lines correspond to the theoretical prediction for the damping rate in the collisional and collisionless limits. The solid line is the result obtained numerically with GS2. Vertical dotted lines denote approximate regions (collisional and collisionless) for which the analytic theory is valid.

Image of FIG. 7.
FIG. 7.

Evolution of ion particle and heat fluxes for an electrostatic, two-species Z-pinch simulation. We are considering and . The particle flux is indicated by the solid line and is given in units of . The heat flux is indicated by the dashed line and is given in units of . We see that a steady state is achieved for both fluxes without artificial dissipation.

Image of FIG. 8.
FIG. 8.

Linear growth rate spectrum of the entropy mode in a Z pinch for , where is major radius and is density gradient scale length. The solid line is the collisionless result, and the other lines represent the result of including collisions. The short dashed line corresponds to using only the Lorentz operator; the dotted line corresponds to the full model collision operator without conserving terms; and the long dashed line corresponds to the full model collision operator with conserving terms. All collisional cases were carried out with .

Tables

Generic image for table
Table I.

Sherman–Morrison variable definitions for Lorentz and energy diffusion operator equations

Loading

Article metrics loading...

/content/aip/journal/pop/16/7/10.1063/1.3155085
2009-07-14
2014-04-19
Loading

Full text loading...

This is a required field
Please enter a valid email address
752b84549af89a08dbdd7fdb8b9568b5 journal.articlezxybnytfddd
Scitation: Linearized model Fokker–Planck collision operators for gyrokinetic simulations. II. Numerical implementation and tests
http://aip.metastore.ingenta.com/content/aip/journal/pop/16/7/10.1063/1.3155085
10.1063/1.3155085
SEARCH_EXPAND_ITEM