Princeton Plasma Physics Laboratory

We present nonlinear magnetohydrodynamic (MHD) simulations of high-β Wendelstein 7-X plasmas using the stellarator extension of the M3D-C1 code, building on the recent work that shows benign saturation of ideal ballooning modes above the designed β limit in the standard configuration [Zhou et al., Phys. Rev. Lett. 133, 135102 (2024)]. First, we examine the results' sensitivity to the parallel thermal conductivity. It is found that while an increased parallel conductivity reduces the linear growth rate, the saturated pressure profile is barely affected. Second, we consider the dependence on the profile shape. It is shown that an equilibrium with a peaked pressure profile and lower β is subject to more significant change than a broad profile with higher β and a larger growth rate, suggesting that benign saturation, or nonlinear stability, is not guaranteed and not dictated by linear growth. Third, we study the influence of the magnetic configuration, with the equilibrium rotational transform varied by adjusting the planar coil current. With similar growth rates, similar magnitudes of profile change are found regardless of the presence of a low-order resonance, which implies that the saturation mechanism is not specific to a resonant or non-resonant mode. These results indicate that MHD stability should still be treated seriously in stellarator operation and design, for which nonlinear modeling using tools like M3D-C1 can play an instrumental role.

Abstract This paper reports on the generic prediction capability of full electromagnetic plasma initiation modelling with DYON, which was carried out for the first time in fusion research by the joint modelling of the International Tokamak Physics Activity—Integrating Operation Scenario group. The following devices were included in the experiment database: VEST (spherical torus, copper coils, Stainless steel wall, R / a = 0.3 m/0.2 m, V v = 3.7 m 3 ), MAST-U (spherical torus, copper coils, C wall, R / a = 0.7 m/0.5 m, V v = 55 m 3 ), EAST (conventional tokamak, superconducting coils, metallic wall, R / a = 1.85 m/0.5 m, V v = 38 m 3 ), DIII-D (conventional tokamak, copper coils, C wall, R / a = 1.67 m/0.65 m, V v = 35 m 3 ), and KSTAR (conventional tokamak, superconducting coils, C wall, R / a = 1.8 m/0.5 m, V v = 55 m 3 ). Despite the different hardware features of the devices, the required operating spaces of the loop voltage induction and prefill gas pressure for inductive plasma initiation in each device were successfully reproduced by the predictive simulations with DYON using only the individual hardware design and the control room input data for each discharge. This successful validation across multiple machines demonstrates that the full electromagnetic DYON modelling can capture the essential physics of inductive plasma initiation. The simulation settings commonly employed for all modelling and the modifications necessary to account for the discrepancies between individual devices are reported. Predictions for ITER based on the multi-machine validation indicate that a wide range of prefill gas pressures exists for the Townsend breakdown and the plasma burn-through (0.01–1.5 mPa).

Abstract Linear and quasilinear magnetohydrodynamic (MHD) modeling is carried out for two DIII-D discharges that both featured a transition from quiescent H-mode (QH) to wide-pedestal QH (WPQH). The MHD perturbations, associated with the edge harmonic oscillations (EHOs) observed during the QH-phase in both discharges, are identified as low- n ( n is the toroidal mode number) resistive kink-peeling instabilities, with the ideal MHD counterpart remaining stable. The quasilinear model successfully simulates EHO-like perturbations during the QH phase in both discharges, confirming the experimental observations. A less intuitive finding is the EHO-like behavior involving the n = 2 perturbation, simulated for one of the discharges during the WPQH-phase. This result, while consistent with experimental observations, is obtained despite the fact that the initial perturbation is linearly stable. The eventual growth of the perturbation is solely due to nonlinear interaction between the MHD perturbation and the plasma toroidal flow. The occurrence of the EHO-like perturbation during the WPQH-phase is found to be sensitive to the initial profile of the plasma edge rotation. For DIII-D plasmas considered, the neoclassical toroidal viscosity, generated by three-dimensional low- n perturbations, is found to play a dominant role in modifying the edge flow during EHOs.

We combine simulations with new analyses that overcome previous pitfalls to explicate how nonhelical mean-field dynamos grow and saturate in unstratified, magnetorotationally driven turbulence. Shear of the mean radial magnetic field amplifies the azimuthal component. Radial fields are regenerated by velocity fluctuations that induce shear of radial magnetic fluctuations, followed by Lorentz and Coriolis forces that source a negative off-diagonal component in the turbulent diffusivity tensor. We present a simple schematic to illustrate this dynamo growth. A different part of the Lorentz force forms a third-order correlator in the mean electromotive force that saturates the dynamo.

Diamond defects are among the most promising qubits. Modeling their properties through accurate quantum mechanical simulations can further their development into robust units of information. We use the recently developed capped density functional embedding theory (capped-DFET) with the multiconfigurational n-electron valence second-order perturbation theory to characterize the electronic excitation energies for different spin manifolds of the well-characterized negatively charged substitutional N defect adjacent to a vacancy (VC) in diamond (NCVC−). We successfully reproduce vertical excitation energies for both triplet and singlet states of NCVC− with errors < 0.1 eV. Unlike other embedding methods, capped-DFET exhibits robust predictions that are approximately independent of the embedded cluster size: it only requires a cluster to contain the defect atoms and their nearest neighbors (as small as a 40-atom capped cluster). Furthermore, our method is free from slowly converging Coulomb interactions between charged defects, and thus also only weakly dependent on supercell size.

In preparation for high field side lower hybrid current drive (HFS LHCD) experiments in DIII-D, the HFS LHCD was to be commissioned and physics experiments commence once the system operated up to 300 kW for 0.5 s. The initial physics experiments sought to characterize coupling, wave propagation, and driven current measurements. In HFS LHCD first campaign, the maximum power was limited to &lt;200 kW due to waveguide pressure leaks limiting the number of available modules and power per module. Here we summarize commissioning progress and initial physics observations. The HFS LHCD coupler was optimized for high qmin, DIII-D discharges where efficient off-axis current at r/a~0.6-0.8 is desired. The coupler n|| spectrum is peaked at 2.7 and is predicted to generate ~0.14 MA/MW coupled for 1.6 MW injected. In preparation, the HFS scrape-off layer density profile was characterized and found to have steeper profiles and lower fluctuation levels than the low field side. Furthermore, the HFS SOL density profile can be accurately predicted using global plasma quantities using machine learning. Thus far, one module has injected ~100 kW for 0.5 s with &lt;5% reflected power. Nonthermal electrons have been observed on lower frequency channels of the electron cyclotron emission radiometer correlated with the LH power indicating core wave absorption. To avoid 30R neutral beam heat flux, a split launcher is proposed to avoid the high heat flux region while maintaining power spectrum and directivity.

Abstract Extended periods of radio pulsations have been observed for six magnetars, displaying characteristics different from those of ordinary pulsars. In this Letter, we argue that radio emission is generated in a closed, twisted magnetic flux bundle originating near the magnetic pole and extending beyond 100 km from the magnetar. The electron–positron flow in the twisted bundle has to carry electric current and, at the same time, experiences a strong drag from the radiation field of the magnetar. This combination forces the plasma into a “radiatively locked” state with a sustained two-stream instability, generating radio emission. We demonstrate this mechanism using novel first-principles simulations that follow the plasma behavior by solving the relativistic Vlasov equation with the discontinuous Galerkin method. First, using one-dimensional simulations, we demonstrate how radiative drag induces the two-stream instability, sustaining turbulent electric fields. When extended to two dimensions, the system produces electromagnetic waves, including superluminal modes capable of escaping the magnetosphere. We measure their frequency and emitted power and incorporate the local simulation results into a global magnetospheric model. The model explains key features of the observed radio emission from magnetars: its appearance after an X-ray outburst, wide pulse profiles, luminosities ∼10 30 erg s −1 , and a broad range of frequencies extending up to ∼100 GHz.

Abstract None