734 research outputs found
Mixed semi-Lagrangian/finite difference methods for plasma simulations
In this paper, we present an efficient algorithm for the long time behavior
of plasma simulations. We will focus on 4D drift-kinetic model, where the
plasma's motion occurs in the plane perpendicular to the magnetic field and can
be governed by the 2D guiding-center model.
Hermite WENO reconstructions, already proposed in \cite{YF15}, are applied
for solving the Vlasov equation. Here we consider an arbitrary computational
domain with an appropriate numerical method for the treatment of boundary
conditions.
Then we apply this algorithm for plasma turbulence simulations. We first
solve the 2D guiding-center model in a D-shape domain and investigate the
numerical stability of the steady state. Then, the 4D drift-kinetic model is
studied with a mixed method, i.e. the semi-Lagrangian method in linear phase
and finite difference method during the nonlinear phase. Numerical results show
that the mixed method is efficient and accurate in linear phase and it is much
stable during the nonlinear phase. Moreover, in practice it has better
conservation properties.Comment: arXiv admin note: text overlap with arXiv:1312.448
On the scaling of entropy viscosity in high order methods
In this work, we outline the entropy viscosity method and discuss how the
choice of scaling influences the size of viscosity for a simple shock problem.
We present examples to illustrate the performance of the entropy viscosity
method under two distinct scalings
A conservative implicit multirate method for hyperbolic problems
This work focuses on the development of a self adjusting multirate strategy
based on an implicit time discretization for the numerical solution of
hyperbolic equations, that could benefit from different time steps in different
areas of the spatial domain. We propose a novel mass conservative multirate
approach, that can be generalized to various implicit time discretization
methods. It is based on flux partitioning, so that flux exchanges between a
cell and its neighbors are balanced. A number of numerical experiments on both
non-linear scalar problems and systems of hyperbolic equations have been
carried out to test the efficiency and accuracy of the proposed approach
- …