Related papers: A fully implicit, conservative, non-linear, electr…
This paper describes a numerical scheme for multi-fluid hydrodynamics in the limit of small mass densities of the charged particles. The inertia of the charged particles can then be neglected, which makes it possible to write an evolution…
The heliosphere is predominantly a partially ionized plasma that consists of electrons, ions and significant neutral atoms. Nonlinear interactions amongst these species take place through direct collision or charge exchange processes. These…
An algorithm for a family of self-starting high-order implicit time integration schemes with controllable numerical dissipation is proposed for both linear and nonlinear transient problems. This work builds on the previous works of the…
A macroscopic model to describe the dynamics of ion transport in ion channels is the Poisson-Nernst-Planck(PNP) equations. In this paper, we develop a finite-difference method for solving PNP equations, which is second-order accurate in…
We present an efficient and accurate energy-conserving implicit particle-in-cell~(PIC) algorithm for the electrostatic Vlasov system, with particular emphasis on its high robustness for simulating complex plasma systems with multiple…
Certain features of the method of characteristics are of considerable interest in relation with Vlasov simulation [H. Abbasi {\it et al}, Phys. Rev. E \textbf{84}, 036702 (2011)]. A Vlasov simulation scheme of this kind can be recurrence…
The fundamental "two-fluid" model for describing plasma dynamics is given by the Euler-Maxwell system, in which compressible ion and electron fluids interact with their own self-consistent electromagnetic field. We prove global stability of…
In this work, we suggest an easy-to-code higher-order finite volume semi-discrete scheme to analyze the nonlinear behavior of the electron-plasma oscillations by solving electron fluid equations numerically. The present method employs a…
A second-order-accurate finite volume method, hybridized by blending an extended double-flux algorithm and a traditionally conservative scheme, is developed. In this scheme, hybrid convective fluxes as well as hybrid interpolation…
The $\delta f$ particle-in-cell algorithm has been a useful tool in studying the physics of plasmas, particularly turbulent magnetized plasmas in the context of gyrokinetics. The reduction in noise due to not having to resolve the full…
Many interesting physical problems described by systems of hyperbolic conservation laws are stiff, and thus impose a very small time-step because of the restrictive CFL stability condition. In this case, one can exploit the superior…
The Energy Conserving semi-implicit method (ECsim), presented by Lapenta in 2017, is a Particle in Cell (PIC) algorithm for the simulation of plasmas. Energy conservation is achieved within a semi-implicit formulation that does not require…
In this paper, we consider the numerical approximation for a diffuse interface model of the two-phase incompressible inductionless magnetohydrodynamics problem. This model consists of Cahn-Hilliard equations, Navier-Stokes equations and…
Phenomena taking place in capacitively coupled plasmas with large electrodes and driven at very high frequencies are studied numerically utilizing a novel energy- and charge-conserving implicit fully electromagnetic particle-in-cell / Monte…
An explicit high-order noncanonical symplectic algorithm for ideal two-fluid systems is developed. The fluid is discretized as particles in the Lagrangian description, while the electromagnetic fields and internal energy are treated as…
We consider the solution of the fully kinetic (including electrons) Vlasov-Amp\`ere system in a one-dimensional physical space and two-dimensional velocity space (1D-2V) for an arbitrary number of species with a time-implicit Eulerian…
Most plasmas are only partially ionized. To better understand the dynamics of these plasmas, the behaviors of a mixture of neutral species and plasma in ideal magnetohydrodynamic states are investigated. The current approach is about the…
The numerical solution of the many-body problem of interacting electrons and ions is a key challenge in condensed matter physics, chemistry, and materials science. Traditional methods to solve the multi-component quantum Hamiltonian are…
A basic model for describing plasma dynamics is given by the Euler-Maxwell system, in which compressible ion and electron fluids interact with their own self-consistent electromagnetic field. In this paper we consider the "one-fluid"…
Conventional explicit electromagnetic particle-in-cell (PIC) algorithms do not conserve discrete energy exactly. Time-centered fully implicit PIC algorithms can conserve discrete energy exactly, but may introduce large dispersion errors in…