2026-10-16 | | Total: 5
Quantitative comparison between SMBI fueling simulations and experiments has long been hampered by the poor constraint on the effective source. Here, the particle increment per pulse is extracted experimentally and used to calibrate the simulation source. For EAST shot 71438, the density evolution from microwave reflectometry is analyzed for nine SMBI pulse groups. The SMBI-induced density increase is isolated from the background by comparing the rising and falling phases around each pulse group. Its radial profile peaks in the mid-pedestal region for every pulse group. Volume integration of this increment provides the experimental reference. The simulations use a multi-fluid model in the BOUT++ framework in which two transport operators are recast into conservative flux form, with the particle balance residual reduced to about 2% of the injected amount, and the perpendicular transport coefficients are constructed from the discrete radial flux balance, with the unphysical changes of the initial density and temperature profiles in source-free evolution reduced by an order of magnitude. For three cases sharing the same nominal molecular flux, the simulated deposition peak agrees with the experimental one in the relative pedestal coordinate. The simulated conversion of the injected particles at the end of the injection is slightly lower than the value inferred from the experiment, while the fraction retained at the end of the pulse is higher. The perturbation thus decays more slowly in the simulation, which we attribute to fast loss channels not represented in the model, such as parallel losses to the scrape-off layer. The edge electron temperature stays on the high-temperature plateau of the dissociation and ionization rate coefficients throughout the pulse, and the inward shift of the deposition with edge cooling reported earlier does not occur under the present conditions.
In a Z-pinch dynamic hohlraum, an aluminum liner collides with a central foam and converts kinetic energy into radiation, and the liner-foam response depends on the aluminum equation of state (EOS). We construct a wide-range aluminum equation-of-state table from density functional theory (DFT) simulations over solid-liquid-plasma regions, denoted DFT-eos. Mixed deterministic-stochastic finite-temperature DFT (mDFT) gives the electron thermal contribution; quantum molecular dynamics (QMD) with classical calibration gives the ion part. We perform one-dimensional cylindrical radiation magnetohydrodynamics (MHD) simulations at a peak current of 50 MA and rise times of 220-300 ns on a single-layer light foam load. Opacity and electrical conductivity are fixed; only the liner Al EOS is changed between DFT-eos and the quotidian EOS (QEOS). Impact time and peak liner kinetic energy change little between the two tables. Larger differences appear in flight-segment density, the contact state at impact, and the time history of outgoing radiation power at the outer interface.
We examine the consistency of the equilibrium magnetic field $\mathbf{B}_0$ used in the analytical calculation of the sideways wall force by A. A. Martynov and S. Yu. Medvedev [Phys. Plasmas 27, 012508 (2020)] and the subsequent reformulation presented in [Phys. Plasmas 27, 054702 (2020)]. In these attempts, the current density $\mathbf{j}_0 = \nabla \times \mathbf{B}_0$ must vanish in the vacuum regions adjacent to the wall, making $\mathbf{B}_0$ curl-free as a complement to the natural $\nabla \cdot \mathbf{B}_0 = 0$. Since the cited works incorporate selected toroidal corrections without explicitly specifying the corresponding differential operators, we test $\mathbf{B}_0$ adopted in these models against Maxwell's equations using both toroidal and quasi-cylindrical formulations. We show that the divergence-free requirement is violated in the toroidal interpretation, while the curl-free condition is not satisfied in the quasi-cylindrical limit. The reformulation likewise fails to define a self-consistent equilibrium field with toroidal effects because the accompanying field corrections required for consistency with Maxwell's equations are not specified. Since the sideways force calculations are built upon these inconsistent field prescriptions, their analytical derivations are consequently compromised. The resulting arguments therefore cannot serve as a valid basis for revising either the earlier results established by D. V. Mironov and V. D. Pustovitov [Phys. Plasmas 24, 092508 (2017)] or those obtained within the single-harmonic analysis [Phys. Plasmas 22, 052502 (2015)].
During relativistic magnetic reconnection, X-points ($E>B$ regions at vanishing guide field) have been proposed to play an essential role in injecting particles into the nonthermal population. Using fully kinetic particle-in-cell simulations at a range of magnetizations $σ\equiv B_0^2/(4πn_0 m_e c^2) = 12$--$800$, we study in detail how the particle energy gains up to the injection threshold $γ_{\rm inj} = σ/4$ are partitioned into work done inside $E>B$ regions (unmagnetized), within one gyration after exit (re-magnetizing), and at all other times (magnetized). The magnetized channel supplies the largest part of the energy gain for $\sim 95\%$ of injected tracers at $σ= 12$ and $\sim 65\%$ at $σ= 800$. Most of the tracers in this magnetized majority never enter an $E>B$ region before injection. We also compare these injection shares with the earlier studies that proposed an essential role for $E>B$ regions. We conclude that $E>B$ regions are neither necessary for injection nor the dominant direct contributors to pre-injection work by injection share in two-dimensional, nonradiative pair-plasma reconnection at vanishing guide field.
Cosmic rays (CRs) streaming faster than the Alfvén speed generate Alfvén waves that scatter them. The saturated wave amplitude is regulated by a damping mechanism. This damping controls the CR scattering rate and drift velocity, and thus the rate at which CRs transfer momentum and energy to the ambient plasma, which sets the strength of CR feedback. In the ionized interstellar medium, nonlinear Landau damping (NLLD) is commonly invoked as the dominant damping mechanism for Alfvén waves. We test this assumption with fluid-PIC simulations in which thermal ions and electrons are modeled as fluids with optional Landau closures, allowing us to compare otherwise identical runs with and without NLLD. When waves are self-generated by monoenergetic CRs, runs with and without Landau closures are indistinguishable: NLLD is inefficient and the CRs isotropize in the wave frame. Once high-amplitude Alfvén waves mimicking external turbulence are seeded, NLLD curbs wave growth and equalizes forward- and backward-propagating waves. We show that the saturated drift velocity, which falls below the Alfvén speed, can be accurately predicted from the wave spectra alone. If the CR distribution is undersampled, a numerical instability corrupts the largest resonant scales: unphysical, unresolved pitch-angle gradients resonantly generate backward-propagating waves, driving an inverse cascade that should be absent from purely streaming-driven waves. While power-law CR distributions may permit stronger NLLD, our results indicate that NLLD is unlikely to be an intrinsic saturation mechanism of the streaming instability; instead, its effectiveness hinges on the interplay of streaming-generated modes with those of external origin.