2026-07-07 | | Total: 15
SoPlasmaFoam is an open-source, multi-region plasma-dielectric solver built on OpenFOAM, integrated with the PETSc linear-algebra suite (CPU and GPU back-ends), the blastAMR adaptive-mesh-refinement library (hexahedral and polyhedral meshes), and the ROUND family of high-resolution convective schemes. It solves drift-diffusion-reaction transport for charged species, coupled self-consistently to Poisson's equation explicitly or semi-implicitly, with plasma and dielectric regions joined by a monolithic multi-domain coupling for arbitrary curved interfaces. This work makes three contributions. First, a systematic assessment of convective schemes on a stiff scalar-advection problem and the positive-streamer benchmark shows that Scharfetter-Gummel is stable but excessively diffusive on coarse meshes, while ROUNDF outperforms all tested TVD limiters and is recommended for streamer transport. Second, an analysis of Poisson-transport coupling shows that fixed-point correction loops critically control accuracy, that a semi-implicit Poisson formulation does not remove this requirement, and that coupling must be tightened even when Courant and dielectric-relaxation numbers are well below unity. Third, a drift-robust wall boundary condition acting on discretized matrix coefficients is introduced, remaining accurate in the drift-dominated limit where conventional mixed-boundary mappings fail. The solver is validated against a low-pressure DC glow discharge and the positive-streamer benchmark, and its multi-region capability is demonstrated on a nanosecond surface dielectric barrier discharge. Performance analysis confirms memory-bound finite-volume scaling and shows that with AMR the solver is competitive with the fastest reported plasma codes. The framework provides a modular foundation for multiphysics simulations in plasma-assisted combustion, plasma processing, and plasma-based flow control.
The interaction of ultraintense lasers with magnetized overdense plasmas reveals a fundamentally new regime of laser-driven particle acceleration. Particle-in-cell simulations demonstrate the generation of directional, quasi-monoenergetic electrons in the MeV energy range superimposed on a broad thermal electron background with the estimated acceleration gradient of 3.6 MeV/μm, which is the highest till date. In contrast to conventional laser-plasma accelerators, which rely on underdense plasmas and are therefore constrained to relatively low plasma densities and limited beam charge, the present scheme operates in plasmas with densities orders of magnitude higher, opening new possibilities for the generation of high-flux energetic electron beams. A central result of this work is the demonstration of the excitation of electron Bernstein waves during relativistic laser interaction with magnetized overdense plasmas. The subsequent Landau damping of these electrostatic warm-plasma modes selectively transfers energy to resonant electrons, leading to the emergence of quasi-monoenergetic spectral peaks at energies that can be tuned through the applied magnetic field. To support the simulation results, we experimentally demonstrate the directional emission of energetic electrons from a simple permanent-magnet target irradiated by an ultraintense laser pulse, highlighting the practical feasibility of controlled electron-beam generation in dense plasma environments. These findings establish electron Bernstein waves as an efficient mediator of laser energy coupling in overdense plasmas and introduce a new paradigm for controlled particle acceleration and energy deposition in high-energy-density plasma systems.
This work presents deep learning models to predict magnetohydrodynamic equilibrium parameters and profiles for the ADITYA-U tokamak. A synthetic free-boundary equilibrium dataset consisting of 100,760 cases was generated using the pyIPREQ Grad-Shafranov solver, with inputs derived from 766 ADITYA-U plasma discharges and constrained to experimentally relevant circular limiter plasmas near the flat-top phase. Several deep learning approaches were investigated for predicting scalar equilibrium quantities, one-dimensional safety factor profiles and two-dimensional poloidal flux profiles. These approaches included Dense neural networks, principal component analysis based reduced-order models, one-dimensional and two-dimensional convolutional neural networks, and physics-informed neural networks incorporating Grad-Shafranov residual constraints. In addition, an inverse model was developed to estimate poloidal field coil currents from desired plasma equilibrium conditions. The results demonstrate that key equilibrium parameters and profiles can be accurately estimated within the operational domain represented by the dataset. The developed models provide a computationally efficient alternative to conventional equilibrium estimation and can be useful for real-time plasma control, rapid equilibrium analysis, and experimental planning in ADITYA-U operations.
We present a linear gyrokinetic theory for the electromagnetic collisionless trapped-electron mode (CTEM). It is found that the weak electromagnetic effects of CTEMs originate from the particle dynamics. Theoretical analysis reveals that the kinetic and fluid-like components of the trapped-electron parallel current cancel at leading order. The ion parallel current is also negligible due to the weak ion transit resonance. Consequently, the perturbed parallel current in the electromagnetic CTEM is dominated by passing electrons. We demonstrate that these characteristics of particle dynamics decouple the CTEM from the shear Alfvén wave branch, rendering the electromagnetic effects subdominant. Both eigenmode analyses and gyrokinetic simulations validate these findings.
Using fully kinetic simulations that capture unprecedentedly large (from electron to ion) scales, we study magnetogenesis driven by continuous large-scale forcing until nonlinear dynamo saturation. We uncover a two-stage mechanism in collisionless ion-electron plasmas whose dynamics diverge dramatically from the pair-plasma case. In the first phase, electron pressure anisotropy triggers electron-Weibel modes, seeding small-scale magnetic fields. Then, a second growth phase emerges when the more massive ions develop their own strong anisotropy and drive ion-Weibel-type modes; concurrently, a Biermann-battery mechanism contributes to amplifying the magnetic field. This combined dynamics provides a tenfold amplification of the magnetic field in comparison to the pair-plasma case. Over long times, dynamo action continues until the system reaches a statistical steady state. This self-consistent kinetic mechanism provides a plausible explanation for robust magnetogenesis wherever an external forcing continuously stirs the plasma.
Plasma in nuclear fusion reactors is only partially observable: diagnostics are constrained by limited access, cost, and the harsh plasma environment, while high-fidelity simulations remain prohibitively costly at reactor-relevant scales to address the observability gap. This paper presents an ML model for reconstructing full-domain plasma states from a small number of accessible measurements. The model combines temporal encoding of sparse sensor histories with spatial decoding into complete plasma-field maps. Demonstrations using high-fidelity kinetic simulation data show that multiple coupled plasma quantities can be reconstructed from only a few density sensors, with robustness to sparse and randomly located probes. The approach provides a route toward diagnostic augmentation, real-time state estimation, and data-driven digital twins for fusion-relevant plasmas.
We investigate the dynamics of charged particles interacting with ultra-intense electromagnetic X-modes in strongly magnetized plasmas. We demonstrate that particle motion becomes chaotic for relative wave intensities $δ= B_w/B_0 \gtrsim 0.25$ (not above the field reversal threshold $δ\geq 1$). The transition to chaos occurs via the Chirikov resonance overlap mechanism and the related destruction of Kolmogorov-Arnold-Moser (KAM) tori. The maximum Lyapunov exponent increases logarithmically with $δ$, even though the unmagnetized $δ\to \infty$ limit is strictly integrable. In the $δ\gg 1$ regime, incomplete re-laminarization of the phase space flow leads to two distinct populations: (i) the majority of particles undergoing stochastic diffusion, and (ii) a fraction of particles that become phase-locked with the wave, experiencing macroscopic intermittent surfing (Lévy flights). The 1D Particle-In-Cell simulations using the EPOCH code in the highly magnetized ($σ\gg 1$) and under-dense regime are generally consistent with the Hamiltonian single-particle theory. The dissipation fraction of the initial EM energy remains mild.
The moving multi-mirror (MMM) concept mitigates axial losses in magnetic mirrors by using inward-propagating multi-mirror sections to transport escaping particles back toward the central cell. Single-particle simulations are used to establish the underlying dynamics, while a modified rate-equation model provides quantitative estimates of the resulting confinement enhancement. It is found that the steady-state outgoing flux can be robustly suppressed by several orders of magnitude. The analysis uncovered confinement challenges in the MMM concept, indicating that additional scattering processes are required in both the central cell and the MMM sections to achieve the desired confinement.
Robust control of plasma current ramp-up is an absolute necessity, as an efficient and uncontaminated plasma current ramp-up is essential for achieving prolonged, high-pressure tokamak plasma discharges. In conventional tokamaks with Ohmic breakdown, the plasma current ramp-up is achieved primarily with pre-fixed temporal profiles of the applied toroidal electric field and the equilibrium magnetic field (Bv). The pre-fixed temporal profiles of these fields are often insufficient to maintain a successful plasma current ramp-up, as several unquantified dynamical variables, such as the condition of the vessel wall and plasma-facing components, influence the plasma current rise. Fuel gas injection in an appropriate quantity at a suitable time during the current ramp-up is therefore used to control the plasma current rise rate, ensuring successful plasma current start-up in Aditya-U. The gas injection time and gas quantity are controlled based on real-time measurement of plasma current rise rate using a digital signal processor (DSP) controller. This special control scheme is capable of achieving the plasma current to rise nearly at the desired rate, resulting in a successful start-up and a stable plasma discharge.
The state-of-the-art concerning the modeling of the thermomechanical response of graphite to runaway electron (RE) impact is based on one-way coupled linear thermoelasticity combined with Rankine's strength-based failure criterion and limited to the onset of material failure. Here, the predictive capabilities are extended to the nonlinear damage phase characterized by material fragmentation and debris expulsion. This is achieved by introducing plasticity via the Johnson-Holmquist constitutive model, adopting an effective plastic strain-based failure criterion and coupling finite element analysis with smoothed-particle hydrodynamics. The extended thermomechanical model is successfully benchmarked against the results of two controlled RE-induced damage experiments recently carried out in DIII-D. This constitutes an important step towards the final objective of quantitatively describing the thermomechanical response of tungsten to REs.
This dissertation focuses on Atmosphere-Breathing Electric Propulsion (ABEP) for propellant-less continuous orbiting in very low Earth orbits (VLEO) and very low Mars orbits (VLMO). ABEP utilizes an intake to collect atmospheric particles as propellant for an electric thruster, offsetting aerodynamic drag. An analytical system-level analysis assesses collectible mass flow, drag, and power requirements, using GOCE as a case study. Three intake designs are developed based on gas-surface interactions. Two diffuse-reflection designs yield collection efficiencies $η_c < 0.5$, while a specular-reflection design achieves $η_c < 0.95$. The specular intake demonstrates superior robustness against flow misalignment, maintaining high efficiency at large angles. The proposed ABEP thruster features a contactless, neutralizer-free design using a helicon plasma discharge, preventing erosion from aggressive species like atomic oxygen. Numerical simulations via the HELIC code determined key design parameters (frequency, magnetic field, plasma density), leading to the implementation of a resonant birdcage antenna with high electrical efficiency. The thruster was manufactured and experimentally validated using Ar, N_2, and O_2. Tests demonstrated successful ignition and operation across relevant mass flows with low power consumption ($P_f < 60$ W). Additionally, a B-dot probe was developed to detect helicon waves in the plasma plume. Ultimately, this work delivers verified intake designs and a novel contactless plasma thruster with high electrical efficiency (>99% in vacuum). Both components represent promising technologies for near-future ABEP applications. Future development outlooks, diagnostics, and spacecraft configurations are also discussed.
We present a quantum-kinetic linear theory of the gravitational instability in the context of fuzzy dark matter universe. Starting from the Wigner transport equation, we apply Landau's approach to the linearized Wigner--Poisson system and derive a kinetic dispersion relation that incorporates quantum effects exactly by introducing the plasma dispersion function. The growth rate as a function of wavenumber is characterized by a dimensionless quantum-to-thermal ratio $α= k_{\mathrm{qJ}}/k_{\mathrm J}$, where $k_{\mathrm{qJ}}$ and $k_{\mathrm J}$ represent the quantum and thermal Jeans wavenumbers, respectively. We derive an analytic expression for the spectral slope at the cutoff wavenumber, revealing that the spectral shape undergoes a sharp transition across $α\sim 0.5$. This implies a crossover from a thermally dominated kinetic regime, in which collisionless damping occurs via phase mixing and Landau resonance, to a regime dominated by quantum pressure. By applying these results to fuzzy dark matter, we show that the cutoff scale and its spectral shape depend sensitively on both the particle mass and the initial velocity dispersion, suggesting a method for simultaneously constraining these parameters through observations of the matter power spectrum. This framework provides a theoretical basis for future studies on the transition from early-phase thermal states to the formation of Bose-Einstein condensates in galactic structures.
Particle-in-cell (PIC) simulation is one of the most important research tools in theoretical plasma physics. To solve the motion of charged particles, the Boris method (a.k.a. the Boris integrator/pusher/solver) has been used for more than a half century. Although the Boris solver has good accuracy, the demand for high-accuracy numerical solvers has been increasing, and new integrators have been actively developed. In this contribution, we present novel high-accuracy particle integrators, the hyper Boris integrators, for nonrelativistic kinetic simulations. We further discuss their connection to 3D rotation representations.
Small-scale dynamos (SSDs) are believed to be the primary source of magnetic fields in all turbulent astrophysical systems, especially those with weak rotation such as elliptical galaxies and galaxy clusters. The initial kinematic phase of these dynamos is relatively well understood. Here we demonstrate analytically and numerically that, in an appropriate limit, the eddy-damped quasi-normal Markovian (EDQNM) closure for incompressible magnetohydrodynamic turbulence is strictly equivalent to the earlier models of kinematic dynamos. Moreover, it allows the extension of the kinematic dynamo framework to multi-scale turbulent flows and into the nonlinear regime. The EDQNM closure also enables us to explore a wide parameter range which is inaccessible to direct numerical simulations of the SSD. Using nonhelical EDQNM simulations, we identify several asymptotic regimes of nonlinear dynamo action when the system is highly turbulent with fluid Reynolds number $Re \gtrsim 10^6$ for magnetic Prandtl number $Pm > 1$ and magnetic Reynolds number $Rm \gtrsim 10^6$ for $Pm < 1$: 1) the kinematic growth rate approaches a value independent of $Pm$, 2) the saturated magnetic to kinetic energy ratio similarly converges to $\simeq 0.55$ across $Pm$, while the ratio of magnetic to kinetic integral wavenumbers asymptotes to $\simeq 3$. For all $Pm$, we further find strong feedback between magnetic field and velocity field largely via Alfvénisation leading to a saturated kinetic and magnetic spectra with almost the same inertial range with a slope of $-3/2$. These findings could provide guidance for future global simulations and for modeling the nonlinear regime of astrophysical systems living in these extreme limits.
Magnetic switchbacks are large-amplitude fluctuations in the interplanetary magnetic field, and appear frequently in the near-Sun solar wind explored recently by Parker Solar Probe: these new observations have prompted many new studies into their properties and origins. Here, we first review what is known about how switchbacks evolve as they travel away from the Sun: both in terms of their expansion-driven growth and their decay due to various processes like turbulence, reconnection, dispersion, parametric instability, and interaction with interplanetary shocks. We then review the current state of knowledge on how switchbacks impact the physics of the solar wind as a whole: in terms of the turbulent cascade, acceleration and heating of the wind, modification of the open solar flux and scattering of energetic particles. Finally, we suggest future studies to further our understanding of switchback evolution and impacts on the heliosphere.