2025-12-08 | | Total: 8
Magnetic shear across the polarity inversion line (PIL) plays an important role in the explosive nature of reconnection onset and in the equilibration of current sheets, acting as a source of free energy that can enhance or inhibit the onset process under certain conditions. In this study, we use a 2D PIC simulation to examine the local interaction between the reconnection guide field and thermodynamic variables during reconnection onset in a region of initially depleted thermal energy and enhanced magnetic energy in a large guide field background. We identify critical stages of the equilibration process, characterize intervals based on whether the pressure evolution is driven by changes in density or temperature, and discuss what these intervals imply about the evolution of local heat and work density. Finally, we examine power densities associated with electromagnetic field time evolution and electromagnetic energy transfer and compare to those related to thermodynamic changes.
We consider a pair of driver/witness electron bunches propagating in an ionized gas background a configuration similar to the one produced in a capillary discharge where a plasma oscillation has been excited by a driving pulse. We assume as in the plasma nonlinear regime that the plasma electrons behind the driver are completely expelled and an ellipsoidal cavity filled with ions only is formed. The fields are linear in both longitudinal and transverse directions, at least in the region of interest for particle acceleration, as the one produced by a uniform ion distribution within a uniformly charged ellipsoidal volume. The fields produced by the ions and experienced by a witness electron beam are purely electrostatic, being the ions at rest in the laboratory frame on the time scale of interest and it can be represented with the field distribution produced by a 3D charged ellipsoidal. The energy spread and emittance degradation has been studied by slicing the bunch in an array of cylinders and solving envelope equations for each bunch slice. The properties of transverse envelope and emittance oscillations and energy spread degradation have been analyzed together with the related matching conditions for optimal transport and acceleration.
There is an ongoing interest in using spectroscopy in inertial confinement fusion (ICF) experiments, where dopants such as krypton can provide vital information about the temperature and density of the imploding plasma. While the most advanced tools for calculating Stark profiles are computer simulation models (CSMs), their application to complex lineshapes under the extreme conditions of ICF experiments is computationally challenging. In this manuscript, we present results of several CSM realizations applied to the Stark shape of the krypton He-beta line and its satellites at ICF-relevant conditions (ne = 1e24 to 1e25 cm-3, Te = 3keV). We demonstrate that codes with the same underlying physics but different numerical approaches yield identical results and analyze the differences in the line profile caused by various physical effects.
Magnetohydrodynamic (MHD) instabilities and runaway electrons (REs) interact in several ways, making it important to self-consistently model these interactions for accurate predictions of RE generation and the design of mitigation strategies, such as massive gas injection (MGI). Using M3D-C1 - an extended MHD code with a RE fluid model - we investigate the effects of 3-D nonlinear MHD activity, material injection, and 2-D axisymmetric vertical displacement events (VDEs) on RE evolution during disruptions on SPARC - a high-field, high-current tokamak designed to achieve a fusion gain Q > 1. Several cases, comprising different combinations of neon (Ne) and deuterium ($\text{D}_2$) injection, are considered. Our results demonstrate key effects that arise from the self-consistent RE + MHD coupling, such as an initial increase in RE generation due to MHD instability growth, decreased saturation energies of the m/n = 1/1 mode driving sawteeth-like activity, RE losses in stochastic magnetic fields, and subsequent RE confinement and plateau formation due to re-healing of flux surfaces. Large RE plateaus (>5 MA) are obtained with Ne-only injection ($2-5 \times 10^{21}$ atoms), while combined $\text{D}_2$ + Ne injection ($2 \times 10^{21}$ Ne atoms; $1.8 \times 10^{22} \, \text{D}_2$ molecules) produces a lower RE current (<2 MA). With $\text{D}_2$ + Ne injection, a post thermal quench "cold" VDE terminates the RE beam, preventing a steady plateau. These simulations couple REs, 3-D MHD instabilities, MGI, and axisymmetric VDEs for the first time in SPARC disruption simulations and represent a crucial step in understanding RE generation and mitigation in high-current devices like SPARC.
Small-scale turbulence originating from microinstabilities limits the energy confinement time in magnetic confinement fusion. Here we develop a semi-analytical dispersion relation based on lowest-order solutions to the gyrokinetic equations in an asymptotic expansion in the ratio of transit (bounce) frequency to the mode frequency for ions (electrons), capable of describing two common instabilities: the ion temperature gradient (ITG) mode and trapped-electron mode (TEM), in the electrostatic limit. The dispersion relation, which is valid in arbitrary toroidal geometry, takes into account resonances with the magnetic ion and bounce-averaged electron drifts, incorporates non-local effects along the magnetic field line, is valid for arbitrary sign of the growth rate and magnetic curvature, and is shown to satisfy a variational property. Several common approximation models are introduced for both the magnetic drift and finite Larmor radius (FLR) damping, with the Padé approximation for FLR effect in particular resulting in remarkable agreement with the baseline dispersion relation model at significantly reduced costs. The baseline model is verified by comparing solutions of the dispersion relation model to high-fidelity linear gyrokinetic simulations, where the exact eigenfunction of the electrostatic potential from simulations is used as a trial function, showing good quantitative agreement for ITGs and TEMs in (shaped) tokamaks as well as low-magnetic-shear stellarators.
A direct quantitative comparison of Scrape-Off Layer (SOL) filament properties from fluid turbulence simulations using the GBS code and from experiments on the TCV tokamak is performed within the TCV-X21 validation case. This comparison is made possible by extending the open TCV-X21 dataset with 2D turbulence measurements obtained with Gas Puff Imaging (GPI), providing critical information on the size, velocity, and other key characteristics of turbulent filaments at the outboard midplane and in the divertor region. For the comparison, GBS simulations of TCV-X21 are analyzed using a dedicated synthetic GPI diagnostic that models the neutral helium-plasma interaction and emission processes and accounts for line-integration effects. Poloidal and radial filament velocities are found to be in good agreement between simulations and experiments, while the simulations overestimate the filament radial and poloidal sizes and underestimate the relative fluctuation levels. The simulations further indicate that filaments in the SOL are predominantly represented by density perturbations rather than temperature perturbations, consistent with previous assumptions in experimental analyses of cross-field turbulent transport from GPI data. The poloidal velocity direction of the filaments agrees with the time-averaged $\boldsymbol{E}\times\boldsymbol{B}$ direction at the outboard midplane and X-point region, but not in the divertor leg. Possible explanations are proposed and discussed, highlighting the influence of the instantaneous $\boldsymbol{E}\times\boldsymbol{B}$ velocity components in both poloidal and radial directions. This study provides new insights into turbulent filament behavior and contributes to guiding future efforts to improve first-principles simulations of the boundary plasma.
A novel, conservative discontinuous Galerkin algorithm is presented for particle kinetics on manifolds. The motion of particles on the manifold is represented using using both canonical and non-canonical Hamiltonian formulations. Our schemes apply to either formulations, but the canonical formulation results in a particularly efficient scheme that also conserves particle density and energy exactly. The collisionless update is coupled to a Bhatnagar-Gross-Krook (BGK) collision operator that provides a simplified model for relaxation to local thermodynamic equilibrium. An iterative scheme is constructed to ensure collisional invariants (density, momentum and energy) are preserved numerically. Rotation of the manifold is incorporated by modifying the Hamiltonian while ensuring a canonical formulation. Several test problems, including a kinetic version of the classical Sod-shock problem, Kelvin-Helmholtz instability on the surfaces of a sphere and a paraboloid, with and without rotations, is presented. A prospectus for further development of this approach to simulation of kinetic theory in general relativity is presented.
A number of simulations have seen the emergence of strongly-toroidally-magnetized accretion disks from interstellar medium inflows. Recently, Guo et al. 2025 (G25) studied an idealized test problem of toroidally-magnetized disks in isothermal ideal MHD with an Eulerian static-mesh method, and argued the midplane behavior changes qualitatively (with a significant loss of toroidal magnetic flux) when the the thermal scale-length is resolved ($Δx < H_{\rm thermal}$). We rerun the G25 test problem with two Lagrangian methods: meshless finite-mass, and meshless finite-volume. We show that Lagrangian methods reproduce the high-resolution ($Δx \ll H_{\rm thermal}$) Eulerian G25 results. At low resolution ($Δx \gg H_{\rm thermal}$), behaviors differ: Lagrangian methods still lose flux and evolve 'as close as possible' to the converged solution, while Eulerian methods show no evolution. We argue this difference in convergence behavior is related to the ability of Lagrangian codes to follow flows to an arbitrarily thin midplane layer, analogous to the well-studied difference in Jeans fragmentation problems. This and results from other higher-resolution simulations and different codes suggest that the sustained midplane toroidal fields seen in recent Lagrangian multi-scale, multi-physics simulations cannot be a numerical resolution effect, and some physical difference between those simulations and the G25 test problem explains their different behaviors.