Related papers: A fully implicit, conservative, non-linear, electr…
We discuss the interpolation of the electric and magnetic fields within a charge-conserving Particle-In-Cell scheme. The choice of the interpolation procedure for the fields acting on a particle can be constrained by analyzing conservation…
Linear kinetic Monte Carlo particle transport models are frequently employed in fusion plasma simulations to quantify atomic and surface effects on the main plasma flow dynamics. Separate codes are used for transport of neutral particles…
Several methods in nonadiabatic molecular dynamics are based on Madelung's hydrodynamic description of nuclear motion, while the electronic component is treated as a finite-dimensional quantum system. In this context, the quantum potential…
This paper studies mechanism of preconcentration of charged particles in a straight micro-channel embedded with permselective membranes, by numerically solving coupled transport equations of ions, charged particles and solvent fluid without…
Particle-in-cell merging algorithms aim to resample dynamically the six-dimensional phase space occupied by particles without distorting substantially the physical description of the system. Whereas various approaches have been proposed in…
A scheme which allows to compute the dynamics of strongly correlated classical ions embedded into a partiallzy ionized quantum plasma by first principles molecular dynamics is presented. The dynamically screened dust approach of Joyce and…
In this article we develop an algorithm for the efficient simulation of electrolytes in the presence of physical boundaries. In previous work the Discrete Ion Stochastic Continuum Overdamped Solvent (DISCOS) algorithm was derived for triply…
We consider the periodic problem for two-fluid non-isentropic Euler-Maxwell systems in plasmas. By means of suitable choices of symmetrizers and an induction argument on the order of the time-space derivatives of solutions in energy…
We describe kinetic simulations of transient problems in partially ionized weakly-collisional plasma around spherical bodies absorbing or emitting charged particles. Numerical solutions of kinetic equations for electrons and ions in 1D2V…
This article studies the vortical nature and structure of phase-space holes -- nonlinear B.G.K. trapping modes found in the phase-space collision-free plasmas. A fluid-like outlook of the particles' phase-space is explored, which makes it…
The trapezoidal rule, which is a special case of the Newmark family of algorithms, is one of the most widely used methods for transient hyperbolic problems. In this work, we show that this rule conserves linear and angular momenta and…
In this work, an implicit scheme for particle-in-cell/Fourier electromagnetic simulations is developed and applied to studies of Alfv\'en waves in one dimension and three-dimensional tokamak plasmas. An analytical treatment is introduced to…
We present a six-moment multi-fluid model, which solves the governing equations for both ions and electrons, with pressure anisotropy along and perpendicular to the magnetic field direction, as well as the complete set of Maxwell equations.…
Fredholm integral equations of the first kind are the prototypical example of ill-posed linear inverse problems. They model, among other things, reconstruction of distorted noisy observations and indirect density estimation and also appear…
While magnetic confinement fusion (MCF) and inertial confinement fusion (ICF) remain the primary routes toward controlled fusion, progress is still constrained by energy loss, plasma instabilities, and the cost and complexity of large-scale…
In this paper, we present and analyze an energy-conserving and linearly implicit scheme for solving the nonlinear wave equations. Optimal error estimates in time and superconvergent error estimates in space are established without time-step…
A conservative finite-volume framework, based on a collocated variable arrangement, for the simulation of flows at all speeds, applicable to incompressible, ideal-gas and real-gas fluids is proposed in conjunction with a fully-coupled…
Common time-explicit numerical methods for kinetic simulations of plasmas in the low-collisions limit fall into two classes of algorithms: momentum conserving and energy conserving. Each has certain drawbacks. The PIC algorithm does not…
Radiation hydrodynamical transport of ionization fronts in the next generation of cosmological reionization simulations holds the promise of predicting UV escape fractions from first principles as well as investigating the role of…
In this paper, we propose a meshfree approximation method for the implicit filter developed in [2], which is a novel numerical algorithm for nonlinear filtering problems. The implicit filter approximates conditional distributions in the…