2026-08-14 | | Total: 9
The macroscopic ideal-MHD stability of an axisymmetric mirror device with sonic levels of plasma rotation is analyzed by approximating the plasma equilibrium as a rotating theta pinch possessing an artificial gravity. An eigenmode equation is derived that governs the stability of the equilibrium to small perturbations in the case of an arbitrary plasma angular velocity profile. The stability of the m=1 and m=2 modes is investigated. The plasma is found to be stable to these two modes provided that it is sufficiently short in the axial direction. The critical axial length of the device below which the modes are stabilized first decreases with increasing plasma rotation, attains a minimum value when the rotation is roughly sonic, and then increases with increasing plasma rotation. The value of the plasma rotation off the magnetic axis is found to have a significantly stronger effect on the stability of the modes than the value on the magnetic axis.
In this work, a three-dimensional computational fluid dynamics model is developed for a vortex-stabilized microwave CO$_2$ plasma reactor operating over the pressure range of 100-400 mbar. The model combines experimentally constrained, emission-based plasma sizes and volumetric heat-source distributions with thermally dominated finite-rate heavy-particle chemistry for a multi-component mixture. Turbulent flow and transport are described using the SST k-omega model. The model reproduces the measured radial gas-temperature profiles in the plasma core and the non-monotonic pressure dependence of CO$_2$ conversion, including a maximum at 150 mbar and a pronounced decrease at 400 mbar. A vortex-driven recirculation region redistributes gas upstream. Turbulent mixing and cooling are strongest near the upper reactor boundary, but their contribution decreases as pressure increases. The pressure dependence of conversion is determined by the competition between CO$_2$ dissociation and CO recombination. CO$_2$ direct dissociation reaction dominates in the high-temperature plasma core, whereas O-assisted conversion reaction contributes near the plasma edges and in the surrounding hot region. At 150 mbar, enhanced CO$_2$ dissociation is accompanied by limited CO loss, resulting in the highest conversion. With pressure increasing to 400 mbar, slower cooling and more frequent three-body collisions promote CO recombination in the afterglow, causing more than 60% of the CO formed near the plasma to be lost downstream. Moreover, additional CO loss occurs in the upper region of the reactor at higher pressures because of the reduced cooling rate.
In previous work, we introduced a semi-analytical method for computing time-domain solutions of linearized Vlasov problems. Rather than representing the plasma response as a sum of residues associated with Landau poles, the method constructs a regularized frequency-domain response spectrum that is subsequently inverted numerically. Explicit applications have so far been limited to unmagnetized plasmas. In this work, we extend the construction to plasmas with a uniform background magnetic field and a Maxwellian equilibrium distribution. We apply the resulting formulation to the electrostatic ion-Bernstein density response of a plasma with kinetic ions and adiabatic electrons, and provide error estimates for truncating both the spectral integration domain and the cyclotron-harmonic expansion. The solution serves as a time-domain reference for damped finite-$k_\parallel$ regimes, in which dispersion-relation roots alone are insufficient for pointwise-in-time verification of simulation codes. Finally, we indicate how the framework can be extended to fully electromagnetic problems.
Nuclear batteries powered by alpha decay have been deployed successfully for over 60 years, on a worldwide $^{238}$Pu supply of kilograms per year. We show that the 14 MeV neutrons of a single deuterium-tritium fusion plant can produce alpha emitter battery fuels up to tons per year, in three classes: fuels with completely new production pathways ($^{236}$Pu, $^{227}$Ac, $^{210}$Pb), fuels previously proposed whose scarce feedstock the same pathways now breed at scale ($^{232}$U, $^{228}$Th), and the established $^{238}$Pu. OpenMC simulations of actinide channels in a tokamak blanket give, per GW yr of fusion: 11 to 57 kg of $^{236}$Pu, whose chain releases 18 GJ per gram over a century, ending at stable $^{208}$Pb, plus up to 5.2 t of co-product $^{238}$Pu; up to 1.4 t of $^{231}$Pa from thorium, and, from channel with $^{231}$Pa feedstock, up to $\sim$15 t of $^{232}$U or $\sim$122 kg of $^{210}$Pb, with $^{227}$Ac produced at 21 g/yr per tonne of $^{231}$Pa. Neutron capture also upgrades $^{241}$Am to a $^{242}$Cm/$^{242m}$Am/$^{241}$Am/$^{238}$Pu blend with up to 10 times higher power density. The same $^{236}$Pu and $^{232}$U also serve as proliferation safeguards: the $^{237}$Np, $^{232}$Th, and $^{231}$Pa channel products are self-protecting, the plutonium by $^{236}$Pu and $^{238}$Pu decay heat and the 2.6 MeV gammas from $^{208}$Tl content, and similarly the uranium from its $^{232}$U. Many of these fuels ($^{236}$Pu, $^{232}$U, $^{228}$Th, $^{227}$Ac) have an order of magnitude higher power and energy density than current alpha emitters, and at human spaceflight-relevant doses the $^{227}$Ac and $^{210}$Pb chains need less shield mass than $^{238}$Pu or $^{241}$Am above a few hundred watts. Fusion neutrons could therefore enable nuclear batteries at the kilowatt to megawatt scale and unlock new possibilities for power sources requiring exceptionally high energy density.
Compact heavy ion accelerators have numerous applications, ranging from heavy ion fusion to carbon ion radiotherapy, and testing radiation-hardened electronics. The demand could be met by developing high-gradient traveling wave plasma accelerators of high-charge ($\simμ\mathrm{C}$) relativistic ion beams. We will discuss a novel ion acceleration regime -- Counter-propagating ionization Front Acceleration (CFA) -- utilizing counter-propagating Ionization Front (IF) and high-current Relativistic Electron Beam (REB). Theoretical modeling and 3D PIC simulations demonstrate the possibility of using typical REBs produced by induction voltage adders propagating through a gas-filled tube undergoing laser ionization to achieve acceleration gradients in excess of $\sim 250 {\rm MeV/m}$ while accelerating micro-Coulombs of ions over meters distance. A unique energy conversion mechanism -- from the REB to electromagnetic fields to the ions is discussed, as well as the limits on the accelerated ions charge and the degree of its neutralization, acceleration gradient, and ion energy spread.
Accurate numerical simulation of rarefied hypersonic plasmas is increasingly important for optimization of re-entry spacecraft design and the development of advanced aerospace technologies. For kinetic simulation methods, it is convention to enforce ions and electrons to diffuse at the same rate, known as the ambipolar diffusion approximation. This approach circumvents costly resolution of fast electron motion, but neglects the complex plasma dynamics of ions and electrons. Almost all studies that investigated the efficacy of the ambipolar diffusion approximation in hypersonics report noticeable differences in flowfield properties when electrostatic modeling is used, including increases in vehicle surface heat flux and decreases in electron temperature. However, it is unknown whether these reported differences originate directly from acceleration and deceleration of charged species through the electric fields and momentum-exchange collisions between charged and neutral species, defined as first-order effects, or from subsequent interactions with particles experiencing first-order effects, defined as second-order effects. Kinetic hypersonic flow simulations with electrostatic modeling are performed with argon to quantify the validity of the ambipolar diffusion approximation in terms of capturing first-order plasma effects along a one-dimensional stagnation streamline. Three different plasma diffusion regimes are studied under two sets of rarefied freestream flow conditions. The approximation is evaluated in terms of predicting plasma density distributions, electron temperature, and stagnation point heat flux. New criteria are proposed for identification of plasma diffusion regimes in hypersonic flows and use of the ambipolar diffusion approximation.
Type Ia supernovae (SNe Ia) explosively destroy carbon-oxygen white dwarfs (WDs) in multiple stellar systems. They produce approximately 50% of the iron-group elements in the Universe, synthesize electron-capture (EC) elements, drive nuclear physics experiments, and underpin high-precision cosmology. To first order, the outcome is governed by nuclear physics, a property often described as stellar amnesia. Recently, this stellar amnesia has begun to be broken by the nearly universal detection of EC elements with JWST. These elements trace high-density burning, largely ruling out the currently popular helium-triggered, sub-Mch detonation models as the dominant channel. Instead, the ubiquitous presence of EC is shifting back the focus to dynamical and secular mergers, and near-Mch explosions similar to the deflagration model W7, but in which the nuclear flame undergoes a deflagration-to-detonation transition. The early deflagration phase is especially important because spherical simulations identify the central WD density, and thus the WD mass, as a key parameter governing the explosion. Here, we present detailed magneto-hydrodynamical simulations. We find that small-scale, pre-existing turbulence expected from the pre-explosion smoldering phase is essential for overcoming the fundamental challenges imposed by the intrinsic 3D physics. This turbulence systematically reduces the production of EC elements by about a factor of two, implying the need for WD central densities closer to those associated with accretion-induced collapse to a neutron star. We also demonstrate the effect of magnetic fields near the saturation field strength and highlight the need for higher-precision EC rates at low Ye.
Diffusive shock acceleration, at shocks from coronal mass ejections to supernova-remnant blast waves, presupposes a scattering wave field that the accelerated particles themselves maintain. This self-regulation has not been resolved in situ. We report Parker Solar Probe observations of a fast (~2800 km/s), near-parallel interplanetary shock at 0.24 AU on 2023 March 13 and separate its upstream wave field into four families, a classification not made before at a fast shock near the Sun. Right-hand and left-hand circularly polarized families over a common wavenumber band, with a field-aligned linearly polarized family, are cyclotron-resonant with the suprathermal-to-MeV protons streaming from the shock: the beam drives the field that scatters it, and the measured mean free path, half the precursor scale, leaves the beam anisotropic enough to sustain the drive. Outside this loop lies a weak, oblique, linearly polarized component, a few per cent of the wave power, resolved here for the first time at an in situ foreshock. Its in-phase density and field-magnitude fluctuations identify the compressive part as fast magnetosonic and shift the cyclotron-resonance energies of the resonant families by up to 13 % along the precursor. Acceleration at shocks inside 0.3 AU is governed upstream, in a foreshock the shock builds for itself.
At every scale they occupy, magnetic fields affect various phenomena, including star formation, cosmic ray transport, charged particle acceleration, space weather, transport in planetary atmospheres, and laboratory plasmas. These fields are often generated and sustained by turbulent flows in a process called the dynamo. In 1955, E. N. Parker parameterized the effects of small-scale turbulence to propose a mean-field dynamo theory. The widely used theory reproduces observed large-scale fields but suffers from difficulty in tuning parameters as they are not justified from first principles: Studies of turbulent flows show tangled magnetic fields, which are folded and fragmented into small-scale structures due to shear-flow straining. Here, considering a shear flow that is unstable and driven, we develop analytic theory and perform three-dimensional (3D), advanced computer simulations of turbulence with up to 4096 x 4096 x 8192 grid points, showing ab initio generation of quasi-periodic, large-scale magnetic fields. The generation occurs via the mean-vorticity effect---an additional mean-field dynamo process postulated in 1990. Crucial to this dynamo is the prior generation of large-scale 3D jets, robustly produced as topologically protected and exact nonlinear solutions of the magnetohydrodynamic equations. The jet-driven dynamo applies to shear-driven laboratory and astrophysical systems. These include binary neutron star mergers, where the reported dynamo likely operates on microsecond timescales to produce in milliseconds some of the strongest magnetic fields in the Universe, providing signals for multimessenger astronomy.