Toy model

ToyDrift

Electrostatic drift toy model for a single ion species in a given background magnetic field.

Overview

ToyDrift is a stripped-down guiding-center model used to isolate the

electrostatic drift part of the dynamics. It is intended for algorithm prototyping and reduced verification problems rather than production drift-kinetic studies.

Use cases

This model is appropriate for:

reduced electrostatic guiding-center benchmarks
testing the field solve plus $E \times B$ pusher
algorithm prototyping before moving to the full drift-kinetic model

Governing equations

PDEs solved by model:

Drift equation:

\frac{\partial f}{\partial t} + \frac{E \times b _{0}}{B _{∥}^{*}} \cdot \frac{\partial f}{\partial X} = 0

Poisson equation:

- \nabla \cdot \nabla ϕ = \int f B_{∥}^{*} d v_{∥} d μ

where $f (X, v_{∥}, μ, t)$ is the guiding center distribution and

E = - \nabla ϕ, B^{*} = B_{0} + ε v_{∥} \nabla \times b_{0}, B_{∥}^{*} = B^{*} \cdot b_{0}

The control variate method can be activated in the Poisson equation; if enabled, the following Poisson equation is solved:

- \nabla \cdot \nabla ϕ = \int (f - f_{0}) B_{∥}^{*} d v_{∥} d μ

Normalization

The reference speed is the ion thermal velocity:

v^= v^_{i}, E^= v^_{i} B^, ϕ^= E^x^.

Discretization

Time integration is performed by the following propagators (in sequence):

struphy.propagators.poisson_solve.PoissonSolve
struphy.propagators.push_guiding_center_bx_estar.PushGuidingCenterBxEstar

Diagnostics

The following scalars are tracked during simulation:

Electrostatic field energy: en_phi
Particle kinetic energy: en_particles
Total energy: en_tot

Example

Create and initialize the toy drift model:

from struphy.models import ToyDrift

model = ToyDrift()
model.em_fields.phi

model.kinetic_ions.var