Related papers: Mixed semi-Lagrangian/finite difference methods fo…
This work explores the potential of using differentiable simulation for learning quadruped locomotion. Differentiable simulation promises fast convergence and stable training by computing low-variance first-order gradients using robot…
We develop new numerical schemes for Vlasov--Poisson equations with high-order accuracy. Our methods are based on a spatially monotonicity-preserving (MP) scheme and are modified suitably so that positivity of the distribution function is…
Discrete variational methods show excellent performance in numerical simulations of different mechanical systems. In this paper, we introduce an iterative procedure for the solution of discrete variational equations for boundary value…
A Hamiltonian two-field gyrofluid model is used to investigate the dynamics of an electron-ion collisionless plasma subject to a strong ambient magnetic field, within a spectral range extending from the magnetohydrodynamic (MHD) scales to…
The inherent complexity of boundary plasma, characterized by multi-scale and multi-physics challenges, has historically restricted high-fidelity simulations to scientific research due to their intensive computational demands. Consequently,…
This paper presents a novel approach for simulating plasma instabilities in tokamak plasmas using the piecewise field-aligned finite element method in combination with the particle-in-cell method. Our method traditionally aligns the…
A general framework for performing event-driven simulations of systems with semi-flexible or rigid bodies interacting under impulsive torques and forces is outlined. Two different approaches are presented. In the first, the dynamics and…
Variational integrators are a special kind of geometric discretisation methods applicable to any system of differential equations that obeys a Lagrangian formulation. In this thesis, variational integrators are developed for several…
An algorithm to obtain data-driven models of oscillatory phenomena in plasma space propulsion systems is presented, based on sparse regression (SINDy) and Pareto front analysis. The algorithm can incorporate physical constraints, use data…
The numerical analysis for the small amplitude motion of an elastic beam with internal damping is investigated in domain with moving ends. An efficient numerical method is constructed to solve this moving boundary problem. The stability and…
Self-consistent modeling of turbulence-driven transport is critical for optimizing confinement in magnetically confined fusion plasmas, such as in tokamaks and stellarators. In particular, capturing the long-term co-evolution of turbulence,…
We study particle acceleration in strongly turbulent pair plasmas using novel 3D Particle-in-Cell simulations, featuring particle injection from an external heat bath and diffusive escape. We demonstrate the formation of steady-state,…
We present a Vlasov-DArwin numerical code (ViDA) specifically designed to address plasma physics problems, where small-scale high accuracy is requested even during the non linear regime to guarantee a clean description of the plasma…
A conservative scheme has been formulated and verified for gyrokinetic particle simulations of electromagnetic waves and instabilities in magnetized plasmas. An electron continuity equation derived from drift kinetic equation is used to…
A surrogate model for particle-in-cell plasma simulations based on a graph neural network is presented. The graph is constructed in such a way as to enable the representation of electromagnetic fields on a fixed spatial grid. The model is…
Controlling instability in plasma is one of the central challenges in fusion energy research. Among the various sources of instability, kinetic effects play a significant role. In this work, we aim to suppress the instability induced by…
A method of solution of the collisionless Vlasov equation, by following collisionless phase point trajectories in phase space, is presented. It is shown that by increasing the number of phase points, without enhancing the resolution of…
We investigate the drift wave -- zonal flow dynamics in a shearless slab geometry with the new flux-balanced Hasegawa-Wakatani model. As in previous Hasegawa-Wakatani models, we observe a sharp transition from a turbulence dominated regime…
We have developed a deterministic conservative solver for the inhomogeneous Fokker-Planck-Landau equation coupled with the Poisson equation, which is a {classical mean-field} primary model for collisional plasmas. Two subproblems, i.e. the…
The majority of studies on multi-scale vortex motions employ a two-dimensional geometry by using a variety of observational and numerical data. This approach limits the understanding the nature of physical processes responsible for vortex…