2026-03-02 | | Total: 7
An automated algorithm to construct island divertors for stellarators is presented and is used to find divertors that meet heat flux requirements determined by engineering and material limits. The algorithm uses just two initial conditions: two starting coordinates on the island separatrix. We leverage the simplicity of the algorithm to explore the divertor parameter space in a fixed magnetic equilibrium. Heat loads are approximated using the field line diffusion model implemented in the FLARE code. Optimal divertor solutions that satisfy engineering requirements are found using a parameter scan and a Bayesian optimization routine. The optimization achieves a 95% reduction in computational cost compared to the parameter scan. The resulting divertors are proven to be robust to varying plasma parameters through simulations with different cross-field heat diffusivities. This work represents a first step towards island divertor optimization for stellarators.
We propose a method for predicting the nonlinear saturation level of convective instabilities in neutral and magnetized fluids. The method combines Gardner's restacking algorithm, which computes the available energy and ground states of collisionless plasmas in phase space, and Lagrangian relaxation, where fluid elements find lower-energy equilibria while preserving local invariants. For the incompressible Rayleigh-Taylor instability, the problem is formally equivalent to Gardner's and the restacking algorithm directly applies in configuration space. To treat compressibility, we follow restacking with Lagrangian relaxation to obtain the ground state, and the results show excellent agreement with direct numerical simulations. Successful extension to the $m=0$ interchange instability in a Z-pinch demonstrates the method's potential as a general framework for estimating the nonlinear extent of convective instabilities, which can facilitate the design and operation of fusion reactors.
We propose a coupled system for the nonlinear interaction between high-frequency, circularly polarized, intense electromagnetic (EM) waves and low-frequency electron-density perturbations, driven by the EM-wave ponderomotive force, in an unmagnetized plasma composed of fully degenerate relativistic electrons and stationary positive ions, including a higher-order correction to the nonlocal nonlinearity. We show that the modulational instability (MI) growth rate associated with the generation of EM envelope solitons gets significantly reduced with a slight increase in either the nonlocal nonlinear correction or the degeneracy parameter. Furthermore, a three-wave temporal model predicts the existence of quasiperiodic and chaotic states of EM solitons while interacting with longitudinal electron density perturbations. We show that the greater the degeneracy (or the higher the contribution from the nonlocal correction), the smaller the instability domain of modulation wave numbers; thus, degeneracy favors the stability of EM soliton evolution. The existence of temporal chaos in a low-dimensional model could be a signature of the development of spatiotemporal chaos in the complete nonlinear model, in which many electromagnetic solitons can be excited and saturated as they interact with electron plasma waves.
Magnetic reconnection in partially ionized plasmas plays a crucial role in a wide range of solar, astrophysical, and laboratory environments. While reconnection in such plasmas is commonly characterized by the ion-neutral coupling strength and the ionization fraction $χ=n_{i}/(n_{i}+n_{n})$, most previous studies have focused primarily on the former. A systematic exploration of the ionization fraction, particularly in combination with ion-neutral coupling, is still lacking. This study presents the first systematic scan of the two-dimensional parameter space defined by ion-neutral collisionality and ionization fraction, enabling investigation of the transition from strongly coupled reconnection to faster, decoupled reconnection. To achieve this, we employ a new three-fluid, five-moment numerical model that treats electrons, ions, and neutrals as separate species on an equal footing. We find that in the strongly coupled regime, the reconnection rate is consistent with a $χ^{1/4}$ scaling. As collisionality decreases, the system transitions to a fast, ionization-independent regime. On the other hand, in all simulations, the current sheet thins down to the ion inertial length di, rather than the expanded hybrid scale $d_{i}χ^{-1/2}$ predicted by analytic fluid theories. The identified critical thickness and the resulting onset of fast reconnection agree reasonably well with recent fully kinetic simulations and laboratory experiments. In addition, we show that, over a wide range of coupling strengths, the ion outflow velocities remain Alfvénic, scaling with the appropriate ion or hybrid Alfvén speed, while the hybrid outflow velocity scales as $χ^{1/2}$ when normalized by ion Alfvén speed.
The impact of magnetic geometry on zonal-flow generation in ion temperature gradient driven turbulence is investigated by means of linear and nonlinear gyrokinetic simulations. The modulation of geodesic curvature on various configurations has revealed amplification of the zonal-flow intensity in relatively smaller geodesic curvature. Based on these findings, a nonlinear proxy model for explorations of a novel magnetic geometry to activate the zonal-flow dynamics is proposed.
We report on the experimental detection of shockwave dynamics produced in a spark discharge, using a nonresonant four-wave mixing optical technique. In particular, we observe the spark-induced local density perturbation across a millimeter-range probe volume, centered on the discharge, via single-shot coherent Rayleigh-Brillouin scattering. We detect the emergence of shock-induced flow velocities, which appear as distinct features in the spectrum, and monitor their dynamic evolution from a few hundred nanoseconds to microseconds after the spark. Finally, we benchmark our measurements against simulations based on a one-dimensional compressible flow model. Our results pave the way for quantitative measurements of highly non-uniform transient flows in challenging environments featuring non-equilibrium gas kinetics.
Plasma-based accelerators are compact and provide high gradients, yet their practical use has been limited by energy gain, stability, beam quality, and energy transfer efficiency. Here, we address several of these challenges simultaneously using a hybrid scheme in which an electron bunch from a laser wakefield accelerator (LWFA) drives a subsequent plasma wakefield accelerator (PWFA) stage with internal witness injection. Close to driver depletion in the PWFA stage, we obtain witness bunches with higher electron energy, reduced energy spread and divergence, and higher angular-spectral charge density compared to LWFA alone. We report energy transformer ratios approaching~2, and about 20\% of the initial energy in the drive beam was transferred to the witness bunch, thereby achieving a driver-to-witness energy transfer efficiency that largely surpasses that of all previous PWFA experiments.