745 research outputs found
GEMPIC: Geometric ElectroMagnetic Particle-In-Cell Methods
We present a novel framework for Finite Element Particle-in-Cell methods
based on the discretization of the underlying Hamiltonian structure of the
Vlasov-Maxwell system. We derive a semi-discrete Poisson bracket, which retains
the defining properties of a bracket, anti-symmetry and the Jacobi identity, as
well as conservation of its Casimir invariants, implying that the semi-discrete
system is still a Hamiltonian system. In order to obtain a fully discrete
Poisson integrator, the semi-discrete bracket is used in conjunction with
Hamiltonian splitting methods for integration in time. Techniques from Finite
Element Exterior Calculus ensure conservation of the divergence of the magnetic
field and Gauss' law as well as stability of the field solver. The resulting
methods are gauge invariant, feature exact charge conservation and show
excellent long-time energy and momentum behaviour. Due to the generality of our
framework, these conservation properties are guaranteed independently of a
particular choice of the Finite Element basis, as long as the corresponding
Finite Element spaces satisfy certain compatibility conditions.Comment: 57 Page
A Lagrangian kinetic model for collisionless magnetic reconnection
A new fully kinetic system is proposed for modeling collisionless magnetic
reconnection. The formulation relies on fundamental principles in Lagrangian
dynamics, in which the inertia of the electron mean flow is neglected in the
expression of the Lagrangian, rather then enforcing a zero electron mass in the
equations of motion. This is done upon splitting the electron velocity into its
mean and fluctuating parts, so that the latter naturally produce the
corresponding pressure tensor. The model exhibits a new Coriolis force term,
which emerges from a change of frame in the electron dynamics. Then, if the
electron heat flux is neglected, the strong electron magnetization limit yields
a hybrid model, in which the electron pressure tensor is frozen into the
electron mean velocity.Comment: 15 pages, no figures. To Appear in Plasma Phys. Control. Fusio
A discontinuous Galerkin method for the Vlasov-Poisson system
A discontinuous Galerkin method for approximating the Vlasov-Poisson system
of equations describing the time evolution of a collisionless plasma is
proposed. The method is mass conservative and, in the case that piecewise
constant functions are used as a basis, the method preserves the positivity of
the electron distribution function and weakly enforces continuity of the
electric field through mesh interfaces and boundary conditions. The performance
of the method is investigated by computing several examples and error estimates
associated system's approximation are stated. In particular, computed results
are benchmarked against established theoretical results for linear advection
and the phenomenon of linear Landau damping for both the Maxwell and Lorentz
distributions. Moreover, two nonlinear problems are considered: nonlinear
Landau damping and a version of the two-stream instability are computed. For
the latter, fine scale details of the resulting long-time BGK-like state are
presented. Conservation laws are examined and various comparisons to theory are
made. The results obtained demonstrate that the discontinuous Galerkin method
is a viable option for integrating the Vlasov-Poisson system.Comment: To appear in Journal for Computational Physics, 2011. 63 pages, 86
figure
Hamiltonian approach to hybrid plasma models
The Hamiltonian structures of several hybrid kinetic-fluid models are
identified explicitly, upon considering collisionless Vlasov dynamics for the
hot particles interacting with a bulk fluid. After presenting different
pressure-coupling schemes for an ordinary fluid interacting with a hot gas, the
paper extends the treatment to account for a fluid plasma interacting with an
energetic ion species. Both current-coupling and pressure-coupling MHD schemes
are treated extensively. In particular, pressure-coupling schemes are shown to
require a transport-like term in the Vlasov kinetic equation, in order for the
Hamiltonian structure to be preserved. The last part of the paper is devoted to
studying the more general case of an energetic ion species interacting with a
neutralizing electron background (hybrid Hall-MHD). Circulation laws and
Casimir functionals are presented explicitly in each case.Comment: 27 pages, no figures. To appear in J. Phys.
- …