Related papers: Mixed semi-Lagrangian/finite difference methods fo…
The quasi-neutral hybrid model with kinetic ions and fluid electrons is a promising approach for bridging the inherent multi-scale nature of many problems in space and laboratory plasmas. Here, a novel, implicit, particle-in-cell based…
Stable, accurate, divergence-free simulation of magnetized supersonic turbulence is a severe test of numerical MHD schemes and has been surprisingly difficult to achieve due to the range of flow conditions present. Here we present a new,…
Magnetic confinement fusion reactors suffer severely from heat and particle losses through turbulent transport, which has inspired the construction of ever larger and more expensive reactors. Numerical simulations are vital to their design…
Understanding the behavior of interacting fermions is of fundamental interest in many fields ranging from condensed matter to high energy physics. Developing numerically efficient and accurate simulation methods is an indispensable part of…
We describe an efficient numerical method for simulating the dynamics and steady states of collective spin systems in the presence of dephasing and decay. The method is based on the Schwinger boson representation of spin operators and uses…
Solving fluid dynamics equations often requires the use of closure relations that account for missing microphysics. For example, when solving equations related to fluid dynamics for systems with a large Reynolds number, sub-grid effects…
We outline a 2D algorithm for solving incompressible flow--structure interaction problems for mixed rigid/soft body representations, within a consistent framework based on the remeshed vortex method. We adopt the one--continuum formulation…
In space and astrophysical plasmas turbulence leads to the development of coherent structures characterized by a strong current density and important magnetic shears. Using hybrid-kinetic simulations of turbulence (3D with different energy…
Computer aided engineering of multi-time-scale plasma systems which exhibit a quasi-steady state solution are challenging due to the large number of time steps required to reach convergence. Machine learning techniques combined with…
The purpose of the present paper is to compare two semi-Lagrangian methods in the context of the four-dimensional Vlasov--Poisson equation. More specifically, our goal is to compare the performance of the more recently developed…
We introduce numerical methods for simulating the diffusive motion of rigid bodies of arbitrary shape immersed in a viscous fluid. We parameterize the orientation of the bodies using normalized quaternions, which are numerically robust,…
This paper studies the nonlinear evolution of magnetic field turbulence in proximity of steady ideal MHD configurations characterized by a small electric current, a small plasma flow, and approximate flux surfaces, a physical setting that…
Perturbative guiding center theory adequately describes the slow drift motion of charged particles in the strongly-magnetized regime characteristic of thermal particle populations in various magnetic fusion devices. However, it breaks down…
Collisional and stochastic wave-particle dynamics in plasmas far from equilibrium are complex, temporally evolving, stochastic processes which are challenging to model. In this work, we extend previous methods coupling differentiable…
3D2V continuum gyrokinetic simulations of electrostatic plasma turbulence in a straight, open-field-line geometry have been performed using the full-$f$ discontinuous-Galerkin code Gkeyll. These simulations include the basic elements of a…
We propose a finite difference scheme to simulate solutions to a certain type of hyperbolic stochastic partial differential equation (HSPDE). These solutions can in turn estimate so called volatility modulated Volterra (VMV) processes and…
We present an efficient mixed finite element method to solve the fourth-order thin film flow equations using moving mesh refinement. The moving mesh strategy is based on harmonic mappings developed by Li et al. [J. Comput. Phys., 170…
This paper considers the discontinuous Galerkin (DG) methods for solving the Vlasov-Maxwell (VM) system, a fundamental model for collisionless magnetized plasma. The DG methods provide accurate numerical description with conservation and…
Understanding turbulence via simplified fluid models is crucial for optimizing magnetic confinement in tokamak devices. In this work, we propose a novel high-order turbulence law that describes the turbulent cascade at the edges of fusion…
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…