Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “general physics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Boosted Dark Matter Search with the ICARUS Detector at the Gran Sasso Underground National Laboratory

The quest to discover dark matter is one of the most pressing topics in modern physics. Evidence of gravitational anomalies is present at different cosmological scales, spanning the stability and rotation of galaxies to the overall density profile of the entire observable Universe. Due to the discrepancy between the mass measured and the gravitational influence on the mass, if General Relativity is to correctly describe gravity at different scales, extra electromagnetically invisible matter termed dark matter must exist. In recent decades, many experiments have been performed to detect dark matter with masses that span many orders of magnitude, from Mass Compact Halo Objects (MACHOs) with masses several orders of magnitude that of our sun to very light particles on the order of $\sim\mu$eV. However, experiments have only been able to constrain the possible parameter space of dark matter models and, in general, probe dark matter models with single constituent dark matter like Weakly Interacting Massive Particles (WIMPs). In the Boosted Dark Matter (BDM) paradigm, there are four particles in the dark sector: $\chi_0$, $\chi_1$, $\chi_2$, and the dark photon $X$. WIMP-like properties are ascribed to $\chi_0$ to describe the physics of the gravitational anomalies, while relativistic components are assigned to $\chi_1$ and $\chi_2$. The inelastic boosted dark matter (iBDM) channel $\chi_1 e^-\rightarrow \chi_2 e^-\rightarrow \chi_1 Xe^-\rightarrow \chi_1 (e^+ e^- ) e^-$, where $\chi_2$ is an excited state of $\chi_1$, has a unique signature that is distinguishable from neutrino interactions, making neutrino experiments a viable search environment for BDM. Using data collected by the ICARUS collaboration with ICARUS T-600 Liquid Argon Time Projection Chamber (LArTPC) during its operation at the Gran Sasso Underground Laboratory (LNGS), the first-ever LArTPC dark matter search using the ICARUS detector is presented. With a total detector exposure of 0.13 kton$\cdot$year, a total of 4 iBDM candidate events were identified by several stages of data scanning. All were rejected due to defined selection criteria and a wire-by-wire $dE/dx$ evaluation, giving a null result. Exclusion plots at 90\%C.L. for the dark photon ($m_X$, $\epsilon$) parameter space for fixed ($m_0,m_1,m_2$) DM mass sets are presented, setting new limits on the dark photon visible decay parameter space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Integrated modeling of boron powder injection for real-time plasma-facing component conditioning

An integrated modeling framework for investigating the application of solid boron (B) powder injection for real-time surface conditioning of plasma-facing components (PFCs) in tokamak environments is presented. Utilizing the DIII-D impurity powder dropper (IPD) setup, this study simulates B powder injection scenarios ranging from milligrams to tens of milligrams per second, corresponding to boron flux rates of 10 20 –10 21 B/s in standard L-mode conditions. The comprehensive modeling approach combines EMC3-EIRENE for simulating the deuterium plasma background and the Dust Injection Simulator (DIS) for the ablation and transport of the boron powder particles. EMC3 trace impurity fluid modeling results show substantial boron transport to the inboard lower divertor, predominantly influenced by the main ion plasma flow. The dependency on powder particle size (5-250 µm) was found to be insignificant for the scenario considered. The effects of erosion and redeposition were considered to reconcile the discrepancies with experimental observations, which saw substantial deposition on the outer divertor plasma-facing components. For this purpose, the WallDYN3D code was updated to include boron sources within the plasma domain and integrated into the modeling framework. The mixed-material migration modeling shows evolving boron deposition patterns, suggesting the formation of mixed B-C layers or predominantly B coverage depending on the powder mass flow rate. While the modeling outcomes at lower B injection rates tend to align with DIII-D experimental observations, the prediction of near-pure boron layers at higher rates has yet to be experimentally verified in the carbon environment of the DIII-D tokamak. The extensive reach of boron layers found in the modeling suggests the need for modeling that encompasses the entire wall geometry for more accurate experimental correlations. This integrated approach sets a precedent for analyzing and applying real-time in-situ boron coating techniques in advanced tokamak scenarios, potentially extendable to ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Nonlinear causality of Israel-Stewart theory with diffusion

We present the first fully nonlinear causality constraints in D = 3 + 1 dimensions for Israel-Stewart theory in the presence of energy and number diffusion in the Eckart and Landau hydrodynamic frames, respectively. These constraints are algebraic inequalities that make no assumption on the underlying geometry of the spacetime or the equation of state. In order to highlight the distinct physical and structural behavior of the two hydrodynamic frames, we discuss the special ultrarelativistic ideal gas equation of state considered in earlier literature in D = 1 + 1 dimensions, and show that our general D = 3 + 1 constraints reduce to their results upon an appropriate choice of angles. For this equation of state in both D = 1 + 1 and D = 3 + 1 dimensions one can show that: (i) there exists a region allowed by nonlinear causality in which the baryon current transitions into a spacelike vector in the Landau frame, and (ii) an analogous argument shows that the solutions of the Eckart frame equations of motion never violate the dominant energy condition, assuming nonlinear causality holds. Furthermore, we then compare our results with those from linearized Israel-Stewart theory and show that the linear causality bounds fail to capture the new physical constraints on energy and number diffusion that are successfully obtained through our nonlinear causality approach.

Quark-gluon plasma↗
DOE OSTI2025

Towards a general subtraction formula for NNLO QCD corrections to processes at hadron colliders: final states with quarks and gluons

We describe the calculation of integrated subtraction terms in the nested soft-collinear subtraction scheme for hadron collider processes with quarks and gluons, thereby extending the results presented in ref. [1]. Although this extension eventually proves to be straightforward, it requires a more careful treatment of certain collinear limits to achieve a compact and physically-transparent final result. We also show that the cancellation of infrared divergences can be organized in such a way that, once soft contributions are removed, it occurs independently for each of the external partons. We consider these results to be important stepping stones on the way to deriving finite remainders of the integrated subtraction terms for fully-general hadron collider processes in the context of the nested soft-collinear subtraction scheme.

NNLO calculations↗
DOE OSTI2026

Overview of the MAST Upgrade physics programme: testing novel concepts at low aspect ratio to inform future devices

The research programme performed on the Mega Amp Spherical Tokamak (MAST) Upgrade experiment has made significant advances in developing the physics understanding of low aspect ratio tokamaks in support of the operation of ITER and design of fusion powerplants. High performance plasma scenarios have been developed to facilitate a broad programme of experiments, in which confinement is constrained by the presence of m/n = 2/1 modes that cause substantial losses of fast ions. The onset of these modes coincides with the q = 2 surface residing in a local minimum in the toroidal current density profile. The maximum electron temperature at the pedestal top, T e,ped is limited with gas fuelling to ∼350 eV to maintain regular ELMs; higher T e,ped results in a transition to a non-stationary ELM-free regime. The operational space of spherical tokamaks has been expanded into small and ELM-free regimes. Strong shaping of the last closed flux surface can induce a transition from large to small ELMs, and ELM suppression with resonant magnetic perturbations has been observed for the first time in a low aspect ratio tokamak. Negative triangularity shaping has induced a transition from ELMy H-mode to a high-performance L-mode regime for the first time in a low aspect ratio tokamak. In studies of fast ion confinement, losses of fast particles due to Global Alfvén Eigenmodes have been identified. Interactions between fast ions generated by off-axis neutral beam injection and thermal neutrals can result in significant losses of fast ions. Experiments with on- and off-axis neutral beam injection exhibit a flux pumping mechanism, where the central safety factor is held to ∼1 in the absence of sawteeth. In studies of pedestal physics, it has been found that elevated main chamber neutral pressures result in an increase in the electron density and reduction in the temperature at the pedestal top. Advances in understanding plasma exhaust include the integration of a high-performance plasma core with detached outer divertors in the X-point target configuration. A newly commissioned lower divertor cryopump reduces the lower divertor neutral pressure by up to 50%, with minimal effect on the main chamber or upper divertor. New measurements and SOLPS-ITER simulations emphasise the importance of plasma–neutral interactions on divertor detachment in the conditions accessible in experiments. Real-time control of the ionisation front location in both divertor chambers independently has been demonstrated in double null experiments, enabled by the tightly baffled divertor chambers.

MAST Upgrade↗
DOE OSTI2025

Fractonic coset construction for spontaneously broken translations

We study the homogeneous breaking of spatial translation symmetry concomitantly with the spontaneous breaking of other internal and spacetime symmetries, including dilations. We use the symmetry-breaking pattern as the only input to derive, via the coset construction, general effective field theories for the symmetry-originated modes associated with Goldstone’s theorem, namely the Nambu-Goldstone candidates. Through explicit computations, we show that integrating out the explicit massive Nambu-Goldstone candidates or imposing symmetric constraints, namely the inverse Higgs constraints, to express massive modes in terms of the massless ones leads to physically distinct effective field theories. This sensitivity to the chosen method can be traced back to the homogeneous breaking of translations, the homogeneous aspect of the breaking induces a mixing between internal and spacetime symmetries at the level of the Lie algebra. This, in turn, leads to subtle discussions about the inverse Higgs constraints, in particular that they lead to a loss of generality in our specific examples. The derived general effective field theories also give rise to a broad class of theories exhibiting emergent enhanced shift symmetries, which constrain the mobility of the modes. The latter are referred to as fractonic modes. Published by the American Physical Society 2025

Chavda, Ameya (ORCID:0000000211731605)↗
DOE OSTI2025

An energetic link between order and strength in metals: A nanocrystalline strength limit in high-entropy alloys and intermetallic compounds

The metallurgy and materials communities have long understood and exploited fundamental links between chemical and structural ordering in metallic solids to tailor their mechanical properties. We extend these ideas to include prediction of the nanocrystalline strength limit in high-entropy alloys and intermetallic compounds, where a breakdown occurs in the classical Hall-Petch strengthening behavior. The highest reported strength achievable through alloying has rapidly climbed and given rise to new classifications of materials with extraordinary properties, with a notable case being nanocrystalline metals. High-entropy alloys (chemically disordered, concentrated solid solutions) and intermetallic compounds are two boundary cases of how tailored order can be used to manipulate mechanical behavior. Here, we show that the complex electronic-structure mechanisms governing the peak strength of alloys and pure metals can be reduced to a few physically meaningful parameters based on their atomic arrangements and used – with no fitting parameters – to predict the maximum strength of these materials. This includes a generalized energy-based accounting for the degree of structural and chemical ordering that allows for rapid and reasonably accurate prediction of peak strength (validated in the nanocrystalline limit) as a function of temperature. Predictions of maximum strength based on the activation energy (with all materials properties derived from DFT calculations or experiments) for a stress-driven transition to an amorphous state is shown to accurately describe the breakdown in Hall-Petch behavior at the smallest crystallite sizes for pure metals, intermetallic compounds, high-entropy alloys, and metallic glasses. Further, this activation energy is also shown to be directly proportional to interstitial electronic charge density, which is a good predictor of ductility, stiffness (moduli), and phase stability in high-entropy alloys and solid metals generally. The proposed framework suggests the possibility of coupling ordering and intrinsic strength to mechanisms like dislocation nucleation, hydrogen embrittlement, and transport properties, such as through correlations between the activation energies for amorphization with stacking-fault and grain boundary energies. It additionally opens the prospect for greatly accelerated structural materials design and development to address materials challenges limiting more sustainable and efficient use of energy.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Emulating 2D Materials with Magnons

Spin waves (magnons) in two-dimensional (2D) materials have received increasing interest due to their unique states and potential for tunability. However, many interesting features of these systems, including Dirac points and topological states, occur at high frequencies, where experimental probes are limited. Here, we study a crystal formed by patterning a hexagonal array of holes in a perpendicularly magnetized thin film. Through simulation, we find that the magnonic band structure imitates that of graphene, but additionally has some kagomelike character and includes a few flat bands. Surprisingly, its nature can be understood using a nine-band tight-binding Hamiltonian. This clear analogy to 2D materials enables band-gap engineering in 2D, topological magnons along 1D phase boundaries, and spectrally isolated modes at 0D point defects. Interestingly, the 1D phase boundaries allow access to the valley degree of freedom through a magnonic analog of the quantum valley Hall insulator. These approaches can be extended to other magnonic systems, but are potentially more general due to the simplicity of the model, which resembles existing results from electron, phonon, photon, and cold-atom systems. This finding brings the physics of spin waves in 2D materials to more experimentally accessible scales, augments it, and outlines a few principles for controlling magnonic states.

Ferrimagnets↗
DOE OSTI2026

Numerical Modeling of a Two-Stage Ocean Current Turbine

The Equinox Ocean Turbines (EQOT) current energy converter has a unique design with power generation in two small-diameter turbines attached to the tips of a large-diameter passive rotor. This configuration offers some key advantages for capturing ocean currents. With no centrally placed generator, almost no reaction torque is required at the nacelle of the main large-diameter rotor, and the small-diameter tip turbine generators operate at a higher speed and lower torque. The physics that determine the performance and loads on the turbine are also unique. The interactions of the flow field between the two stages and the general architecture of the system cannot be captured with traditional mid-fidelity modeling tools. For design iterations and large sets of load cases, it is important to have mid-fidelity models that can capture the important phenomenon with enough accuracy to identify global trends. This work uses a limited set of high-fidelity computational fluid dynamics (CFD) simulations to help inform the selection of and construction of a custom mid-fidelity model. Mid-fidelity modeling approaches were verified by comparing key turbine performance quantities to those found with the CFD model. Hydrodynamic interactions of the two-stage rotor were identified through high-fidelity CFD modeling. This highlighted the impact of the main rotor tip vortex and wake on the secondary rotor apparent inflow. This results in a relative flow rotation and sharp deficit, that change the optimal secondary rotor rotation speed and adds unsteadiness to the blade loading respectively. Multiple mid-fidelity approaches were evaluated for their ability to capture these effects. A simple approximation of the combined-stage performance based on single-stage BEM provides a reasonable rough prediction, especially near the peak TSR values, with some larger discrepancy at higher TSRs. Predicting the combined-stage performance based on single-stage CFD data improves this prediction across the TSR range. Although the combined-stage modeling in OLAF was not successful in this stage of the project, it showed promise as a mid-fidelity method, assuming the parameters can be tuned to account for the significant differences in time and length scales between the main and secondary rotors. This may be addressed through code changes in future work. A significant finding from the OLAF work was the agreement between the vortex core radius values found independently via a parameter space search and via CFD. The technique of using single-stage secondary rotor BEM, with a custom inflow taken from single-stage main rotor CFD or OLAF, provides an efficient method to capture one-way coupled flow interactions. This method provided generally good predictions of the impact of the flow rotation on the secondary rotor but struggled to accurately predict the peaks of the unsteady load progression. Future work could include some superposition of a tuned main rotor trailing edge viscous wake into the custom inflow to better predict this interaction.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Accidental suppression of Wilson coefficients in Higgs coupling

Higgs couplings are essential probes for physics beyond the Standard Model (BSM) since they can be modified by new physics, such as through the Higgs portal interaction |H| 2 Ο. These modifications influence Higgs interactions via dimension-6 operators of the form (∂|H| 2 ) 2 and |H| 6 , which are generally expected to be of comparable size. This paper discusses a phenomenon of accidental suppression, where the |H| 6 coupling is significantly smaller than (∂|H| 2 ) 2 . This suppression, arising from the truncation of the tree-level effective potential, lacks a clear symmetry explanation but persists in portal models. This paper aims to inspire further studies on additional instances of accidental suppression without symmetry explanations or a general framework to characterize such suppression. We also discuss constraints, at the HL-LHC and future colliders, on the Wilson coefficients of the two dimension-6 operators for various benchmark scenarios of the concrete model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Sensor Reduction for Diversion Detection in a Realistic Heat Pipe Microreactor Using Supervised Machine Learning

Microreactors are designed as a smaller, cheaper, and safer alternative to traditional nuclear power plants. Their non-traditional characteristics and prospect of mass production and deployment will likely require new approaches to nuclear safeguards. The primary proliferation concern with microreactors is the diversion of fuel material. Such diversion may produce measurable defects in key physical attributes like neutron flux, which may in turn be detectable using machine learning models. Preliminary work has demonstrated this ability for modeled nominal and diversion scenarios using large quantities of energy integrated neutron flux data. In practice, the number of available sensors for such measurements will be limited and energy integrated flux information will not be available. This work explores the ability of tree-based gradient boosted ensemble models to classify a given microreactor core is nominal or diversion, and determine the number of fuel pins diverted in the case of diversion with reduced numbers of sensors and more realistic detector responses. Classification accuracy of greater than 98% and regression errors as low as 5% of the total number of fuel pins were achieved with as few as 15 sensors, compared to 99% and 4.1% with a maximum of 240 sensors.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Pedestal formation via different trajectories in the stability space in response to the timing scan of neutral beam heating in DIII-D

The frequency of type-I ELMs decreases as the initiation of the neutral beam injection (NBI) heating is delayed with respect to the time when plasma current (I p ) reaches flat-top in the ITER Baseline Scenario discharges in DIII-D. Henceforth, the time gap between the NBI initiation and I p flat-top will be referred to as “heating delay.” As the heating delay is modified, pedestal formation follows different trajectories in the edge current density–pedestal pressure gradient (j edge -∇p e ped ) space from the L-H transition toward the first ELM event. During the stationary phase after the first ELM, the ELM frequency (f ELM ) decreases by a factor of ~2 as the heating delay is increased. A longer pedestal recovery time in the inter-ELM period is observed for the low f ELM discharges as compared to the high f ELM discharges. Both low and high f ELM discharges show nearly identical profiles of electron density and temperature and have a similar MHD stability just before an ELM crash. However, a marked difference is observed in the magnetic spectrogram of the high and low f ELM discharges in response to the variation in the heating delay. The main difference is in the 200–400 kHz range of the magnetic spectra. A quasi-coherent mode (QCM) at 220 kHz and weaker broadband fluctuations are observed in the high f ELM discharges, while only strong broadband fluctuations are prevalent in the low f ELM discharges. ELM-synchronized analysis shows that the time evolution of these modes is different for the high and low f ELM discharges. The localization of both these modes is confirmed at the maximum gradient region of the pedestal. We hypothesize that these modes cause important pedestal transport and that the difference in the pedestal recovery of the high and low f ELM discharges is a result of the difference in transport driven by these modes, as they change with changes in the heating delay. It is demonstrated experimentally for the first time that discharges with similar pedestal parameters can carry the history of the heating delay into the stationary phase and that changes in turbulent-driven transport are a likely cause of changes in f ELM observed with variations of heating delay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Reduced-Order Modeling of Hydrogen Releases from Vent Stacks and with Wind Effects

Here, the physical release behavior of hydrogen is important to understand from a safety and design perspective. The consequences of unignited pressurized gaseous hydrogen plumes exiting vent stacks were considered by extending and modifying existing general hydrogen plume models. Entrainment, vent stack backpressure, and the flow regime of hydrogen exiting the vent were found to be significant factors affecting plume shape and size, but further investigation and validation with unchoked, low-Froude-number flows is recommended to improve the model’s robustness. Additionally, models for the effects of wind on unignited plume momentum and entrainment were added to explore this behavior. Wind was assumed to increase mixing of hydrogen with the ambient air, and to affect the momentum of the released jet. Introducing wind into the plume model led to a shorter plume for all wind and jet directions. A high counter-flowing wind led to non-physical results and challenges in interpreting the visualization. The proposed jet plume wind sub-models (specifically entrainment coefficients) were fit and compared to experimental data of different releases of hydrogen into a wind tunnel, but the quantity of data available and experimental conditions were limited. Thus, collection of more empirical data and for a wider range of conditions is recommended for improvement of the proposed computational models. Developing reduced-order models for these physical phenomena can improve accessibility to predicted physical behavior and the rate at which hydrogen systems can be safely designed and deployed.

entrainment↗
DOE OSTI2025

Dynamic realization of emergent high-dimensional optical vortices

The dimensionality of vortical structures has recently been extended beyond two dimensions, providing additional topological complexity and robustness for high-capacity information processing and turbulence control. The generation of high-dimensional vortical structures has mostly been demonstrated in classical systems through the complex interference of fluidic, acoustic, or electromagnetic waves. However, natural materials rarely support three- or higher-dimensional vortical structures and their physical interactions. Here, we experimentally demonstrate a high-dimensional gradient thickness optical cavity (GTOC) in which the optical coupling of planar metal-dielectric multilayers implements topological interactions across multiple dimensions. At non-trivial topological phases, high-dimensional GTOC induces high-dimensional vortical structures in generalized parameter space in three, four dimensions, and beyond. These emergent high-dimensional vortical structures are observed under electro-optic tomography as optical vortex dynamics in two-dimensional real-space, employing the optical thicknesses of dielectric layers as synthetic dimensions. Our findings hold significant promise for emulating high-dimensional physics and developing active topological photonic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Directional finite difference method for directly solving 3D gyrokinetic field equations with enhanced accuracy

The gyrokinetic (GK) field equation is a three-dimensional (3D) elliptic equation, but it is often simplified to a set of two-dimensional (2D) equations by assuming that the field does not vary along a specific direction. However, this simplification can introduce inevitable 0th-order numerical errors, as nonlinear mode coupling in toroidal geometry can produce undesirable harmonic modes that violate the assumption. In this work, we propose a novel directional finite difference method (FDM) with a local coordinate transformation to better resolve the target field of interest. The directional FDM can accurately solve 3D GK field equations without simplifications, which can overcome the limitations of conventional methods. The accuracy and efficiency of different FDMs are analyzed in great detail for a variety of geometries, from simple 2D Cartesian coordinates to realistic 3D curvilinear coordinates. The 0th-order numerical errors of simplified 2D GK equations were found to be more problematic for low-harmonic modes and low aspect ratio geometries such as spherical tokamaks. On the other hand, the directional 3D FDM can accurately resolve a much wider range of harmonic modes aligned to the direction of interest, including the low-harmonic modes. In conclusion, we demonstrate that the directional 3D FDM is a highly effective algorithm for solving the 3D GK field equations, achieving accuracy improvements of 10 to 100 times or more, particularly for low-harmonic modes in spherical tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Error field identification through torque balance on a saturated island in DIII-D

Measurement of the electromagnetic torque on a magnetic island could be an attractive method for error field identification in the early phase of ITER operation. Previous DIII-D experiments (Strait 2014 Nucl. Fusion 54 073004; Shiraki et al 2015 Plasma Phys. Control. Fusion 57 025016) have demonstrated the principle of this approach using a stationary or slowly rotating island, while recent developments in magnetic data analysis (Sweeney and Strait 2019 Phys. Plasmas 26 012509) allow the field of a rapidly rotating island to be readily distinguished from that of the wall currents induced by its rotation. In a recent experiment, a rotating n = 1 magnetic perturbation forced a saturated magnetic island to rotate, thus sampling all toroidal phases periodically in a single discharge. The phase and amplitude of the error field are inferred from analysis of the time-dependent torque balance on the island, including torques from the error field, the applied magnetic perturbation, and the wall currents induced by rotation of the applied perturbation and the island. Furthermore, results agree well with those from more conventional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Access and limits of RMP ELM suppression with n = 1 fields in DIII-D

This work reports on DIII-D experiments aimed at extending resonant magnetic perturbation (RMP) suppression of edge localized modes (ELMs) to n = 1 fields, where n is the toroidal mode number. Modeling of the 3D ideal MHD plasma response to the RMPs using the GPEC code is used to quantify edge and core resonant fluxes, guiding experimental strategies to increase plasma resilience against core error field penetration, optimize multicoil phasing, and explore higher q 95 operation. In DIII-D, ELM mitigation is regularly observed across a wide range of n = 1 RMP scenarios. A ∼100 ms phase of complete ELM suppression was achieved at q 95 ∼ 3.9 using an odd-parity coil configuration. The suppressed phase exhibited clear signatures of RMP ELM suppression, including the elimination of Dα spikes, increased pedestal rotation, enhanced magnetic response, and elevated broadband density turbulence. An optimized coil configuration for edge-to-core resonant flux did show increased edge resonance indicated by increased density pumpout, but did not yield RMP ELM suppression. At q 95 ∼ 5.1, a bifurcation to a grassy-like ELM regime occurred, while large type-I ELMs persisted. These results demonstrate progress in experimental access to n = 1 RMP ELM suppression in DIII-D, motivating further study for robust access. This work also highlights the potential role of 3D edge stability as well as rational surface alignment in RMP ELM suppression access, which has important implications for the use of low-n RMPs in future reactor-scale devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Linear simulations of a wide pedestal quiescent H-mode plasma with the extended-MHD code NIMROD

We present linear simulations of a Wide Pedestal QH (WPQH)-mode DIII-D plasma with the NIMROD code to address the role of ion multispecies collisionality on WPQH-modes. We use a full extended-MHD model that captures Alfvénic and magnetoacoustic waves, including two-fluid, ion gyroviscosity and cross heat flux. We show that two-fluid and ion gyroviscous effects destabilize electron-directed peeling–ballooning modes in the pedestal. The linear growth rates of the modes decrease when carbon (C) is added to a deuterium plasma because the multispecies collisional effects increase the plasma resistivity. When replacing C with tungsten (W), the multispecies collisional effects increase the resistivity further, as well as the stabilizing effect on the instabilities. The resistivity impact on mode stability is reversed when both two-fluid and ion gyroviscous effects are ignored, which is consistent with previous works, proving that multispecies collisional effects, together with two-fluid and ion gyroviscosity, should be considered to address the stability of plasmas with significant impurity content.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗