Related papers: Linearized model Fokker-Planck collision operators…
A high-order piecewise field-aligned triangular finite element method is developed and implemented for global electromagnetic gyrokinetic particle-in-cell simulations of tokamak plasmas with open field lines. The approach combines locally…
Optimal control theory aims to find an optimal protocol to steer a system between assigned boundary conditions while minimizing a given cost functional in finite time. Equations arising from these types of problems are often non-linear and…
The super simple Vlasov (ssV) code was developed to study instabilities, turbulence, and reconnection in weakly magnetized plasmas, such as the solar wind in the dissipation range and the edge of fusion plasmas. The ssV code overcomes the…
As we aim to control complex systems, use of a simulator in model-based reinforcement learning is becoming more common. However, it has been challenging to overcome the Reality Gap, which comes from nonlinear model bias and susceptibility…
We present a practical numerical method for evaluating the Lagrange multipliers necessary for maintaining a constrained linear geometry of particles in dynamical simulations. The method involves no iterations, and is limited in accuracy…
The equilibrium-dispersive model of the linear GC (gas chromatography) was derived using both assumptions and the method of its derivation different from the known in literature. It was concluded that this model is a specific case of the…
In this paper the normal collision of spherical particles is investigated. The particle interaction is modelled in a macroscopic way using the Hertzian contact force with additional linear damping. The goal of the work is to develop an…
The unified gas kinetic scheme (UGKS) was initially designed to address multiscale challenges in rarefied gas dynamics and then extended to radiative transfert theory, as described by BGK like relaxation models. In this work, we extend its…
Granular systems confined in vertically vibrated shallow horizontal boxes (quasi two-dimensional geometry) present a liquid to solid phase transition when the frequency of the periodic forcing is increased. An effective model, where grains…
Under low-collisionality conditions the isotropic part of the electron velocity distribution function in a plasma becomes non-local and the electrons can be described by a single global distribution function . This is also the regime…
This paper proposes a second-order accurate numerical scheme for the Patlak-Keller-Segel system with various mobilities for the description of chemotaxis. Formulated in a variational structure, the entropy part is novelly discretized by a…
We propose fully discrete, implicit-in-time finite-volume schemes for a general family of non-linear and non-local Fokker-Planck equations with a gradient-flow structure, usually known as aggregation-diffusion equations, in any dimension.…
The application of discontinuous Galerkin (DG) schemes to hyperbolic systems of conservation laws requires a careful interplay between space discretization, carried out with local polynomials and numerical fluxes at inter-cells, and…
A combination of reaction-diffusion models with moving-boundary problems yields a system in which the diffusion (spreading and penetration) and reaction (transformation) evolve the system's state and geometry over time. These systems can be…
Linear gyrokinetic simulations covering the collisional -- collisionless transitional regime of the tearing instability are performed. It is shown that the growth rate scaling with collisionality agrees well with that predicted by a…
The gyrokinetic toroidal code (GTC) has been upgraded for global simulations by coupling the core and scrape-off layer (SOL) regions across the separatrix with field-aligned particle-grid interpolations. A fully kinetic particle pusher for…
Coulomb collisions in plasmas are typically modeled using the Boltzmann collision operator, or its variants, which apply to weakly magnetized plasmas in which the typical gyroradius of particles significantly exceeds the Debye length.…
In our previous work, numerical schemes for a simplified version of 3-wave kinetic equations, in which only the simple forward-cascade terms of the collision operators are kept, have been successfully designed, especially to capture the…
This paper proposes a novel numerical integrator for modeling multispecies Coulomb collisions in kinetic plasmas. The proposed scheme provides an energy-, momentum-, and positivity-preserving particle discretization of the nonlinear Landau…
Without access to the full quantum state, modeling dissipation in an open system requires approximations. The physical soundness of such approximations relies on using realistic microscopic models of dissipation that satisfy completely…