Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dissipative systems”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Plasma heating by electron cyclotron wave and the temperature effects on lower hybrid current drive on EAST

This paper presents the progress in the long—pulse operation of the electron cyclotron (EC) system and the achievements in high—electron temperature plasmas by the combined EC and lower hybrid (LH) waves heating since the EC system was built in 2015. An electron temperature of up to 12 keV with a duration over 100 s was realized by the simultaneous heating of EC and LH waves at the line-averaged density ${\bar n_{\text{e}}}$ ~1.8 × 10 19 m –3 . The plasma heating effect strongly depends on the location of EC power deposition. H-mode plasmas with solely EC wave auxiliary heating have been obtained on EAST for the first time. These H-mode discharges show an enhanced confinement factor H 98(y,2) around 1.0, which is higher than the previous H-mode using LH power alone (Xu et al 2011 Nucl. Fusion 51 072001). The total heating power is very close to the threshold value for L–H transition according to the international tokamak scaling. In addition to the temperature effects inside the separatrix, higher electron temperature produced by EC wave is found to reduce the LH power loss in the scrape-off layer due to collisional absorption, which is beneficial to further increase the LH current drive efficiency. Ray-tracing/Fokker–Planck modeling results indicate that higher electron temperature can shorten the LH wave propagation on EAST in a multiple-pass regime, thus decreasing the collisional dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Polarization-type potential-induced degradation in bifacial PERC modules in the field

This study examines the susceptibility of bifacial glass/glass passivated emitter and rear cell (PERC) modules to potential-induced degradation-polarization (PID-p) in the field. While there are several studies showing PID-p occurring on both front and back faces of bifacial PERC in accelerated tests, we address the yet unclarified behavior in fielded modules. We examine the effects of mounting configuration; specifically, comparing modules mounted near ground and in elevated ground rack configurations. Modules with the cell circuit in -1500 V system voltage configuration, whether mounted on racks about 30 cm above the ground or elevated 2 m high showed mean degradation of 4.5% to 6% in power under standard test conditions over about 2.5 weeks as measured from the front side of the module. This extent of degradation remained sustained for a duration of about 6 months analyzed. Average daytime temperatures of modules in the various mounting configurations were similar and therefore judged to be insufficient to be a primary influence for the modest PID-p rate differences that we observed among mounting configurations. Increased leakage current in the morning suggests morning dew was sustained longer on modules near the ground measured over six months which would be expected to increase the PID-p rate over the long term. However, the main difference seen between the modules on the various mountings during the initial period with up to 6% mean degradation by PID-p was the approximately two times the irradiance from albedo on the rear of modules mounted in elevated ground rack compared to those on the near ground rack. This difference in incident albedo led to a modestly reduced rate of the development of PID-p of the modules on the elevated ground rack. The difference is attributed to the dissipation of PID-p-causing electrical charge by the albedo incident on the module rear. The behavior could be modeled by a sigmoidal equation with consideration of the differences in the insolation on the module rear.

14 SOLAR ENERGY↗
DOE OSTI2020

Enhancement of Thermal Conductivity of Bentonite Buffer Materials with Copper Wires/Meshes for High-Level Radioactive Waste Disposal

In high-level radioactive waste disposal, a heat-generating waste canister is generally encased with a layer of bentonite-based buffer material acting as an engineered barrier to limit water percolation and radionuclide release. The low thermal conductivity of bentonite (~0.5 W/m∙K) combined with a high thermal loading waste package may result in a high surface temperature on the package that can potentially impact the structural integrity of the package itself as well as the surrounding buffer material. We show in this paper that the thermal conductivity of bentonite can be effectively enhanced by embedding copper wires/meshes across the buffer layer to form fully connected high heat conduction pathways. A simple calculation based on Rayleigh’s model shows that a required thermal conductivity of 5 W/m∙K for effective heat dissipation can be achieved simply by adding ~1 vol % of copper wires/meshes into bentonite. As a result, the peak surface temperature on a large waste package such as a dual-purpose canister can be reduced by up to 300°C, leading to a significant reduction in the surface storage time for waste cooling and therefore the overall cost for direct disposal of such waste packages. Because of the ensured full thermal percolation across the buffer layer, copper wires/meshes turn out to be much more effective than any other materials currently suggested (such as graphene or graphite) in enhancing the thermal conductivity of buffer material. Furthermore, the embedded copper wires/meshes can help reinforce the mechanical strength of the buffer material, thus preventing the material from a potential erosion by a possible intrusion of dilute groundwater.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Field-Dependent Ionic Conductivities from Generalized Fluctuation-Dissipation Relations

We derive a relationship for the electric field dependent ionic conductivity in terms of fluctuations of time integrated microscopic variables. We demonstrate this formalism with molecular dynamics simulations of solutions of differing ionic strength with implicit solvent conditions and molten salts. These calculations are aided by a novel nonequilibrium statistical reweighting scheme that allows for the conductivity to be computed as a continuous function of the applied field. In strong electrolytes, we find the fluctuations of the ionic current are Gaussian, and subsequently, the conductivity is constant with applied field. In weaker electrolytes and molten salts, we find the fluctuations of the ionic current are strongly non-Gaussian, and the conductivity increases with applied field. Furthermore, this nonlinear behavior, known phenomenologically for dilute electrolytes as the Onsager-Wien effect, is general and results from the suppression of ionic correlations at large applied fields, as we elucidate through both dynamic and static correlations within nonequilibrium steady states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Hourly dynamic line ratings for existing transmission across the contiguous United States (preliminary results) [Slides]

The thermal rating of a transmission line depends on the rate at which heat generated by current flow through the line can be dissipated into the environment by convective and radiative cooling, keeping the conductor below a maximum operating temperature. Static or seasonal line ratings assume nearly-worst-case weather conditions (high ambient air temperature, high solar irradiance, and low wind speed), resulting in a low risk of "overrating" a line (i.e., allowing more current to flow than the line is designed for) but a high fraction of time when the line is "underrated" (i.e., able to flow more current while staying below the maximum temperature). Dynamic transmission line ratings (DLR) use real-time measurements of weather conditions, line sag, or line temperature to determine the thermal rating of transmission lines at hourly or higher time resolution, allowing for increased ratings during some (cold, windy, and/or dark) periods and correcting for overrating during other (hot, calm, and/or sunny) periods. Many studies have explored the benefits of DLR, but typically for small collections of transmission lines under a limited sample of weather conditions, leaving open the question of how widely applicable the benefits of DLR might be over diverse geographies and weather conditions. Here, we combine geospatial data on tens of thousands of transmission line routes across the contiguous United States with high-resolution hourly historical weather data for 2007-2013 to estimate the impact of DLR and other weather-adjusted rating systems. We compare static and dynamic line ratings to intermediate ratings adjusted by ambient temperature, day/night irradiance, clear-sky irradiance, and/or measured irradiance, illustrating the benefits of incorporating progressively more detailed weather data into the calculated line ratings. Moving from nighttime-adjusted ratings to clear-sky-adjusted ratings provides more benefit than moving from clear-sky-adjusted ratings to measured-irradiance-adjusted ratings. Depending on the assumed static windspeed for irradiance-adjusted-ratings, DLR either corrects for overrated hours or provides substantially increased ratings in the majority of hours on the majority of lines.

14 SOLAR ENERGY↗
DOE OSTI2020

Magnetic reconnection and kinetic waves generated in the Earth's quasi-parallel bow shock

Magnetic reconnection in quasi-parallel shocks, relevant to the Earth’s bow shock, is studied by means of two-dimensional full particle-in-cell simulations. As the Alfvénic Mach number increases, the propagation direction of the waves excited in the transition region changes, and the shock becomes more turbulent with more reconnection sites. In the higher Mach number shock, abundant electron-only reconnection sites are generated with scales on the order of the ion skin depth or less. Non-reconnecting current sheets can also generate electron jets and energy dissipation can occur there as well. However, non-reconnecting current sheets with the magnetic field reversal typically show a smaller energy dissipation rate than reconnecting current sheets. In the shock transition region, two types of waves are responsible for driving reconnection: one has a wavelength around three ion skin depths (d i ), and the other has a wavelength less than 1 d i . Here, electron and ion distribution functions show that in regions where the former type of waves is excited, there are two ion beams and a single-peaked electron distribution. In contrast, in regions where the latter type of waves is excited, there are multiple electron and ion beams. The waves propagating obliquely to the magnetic field bend the magnetic field lines, and magnetic reconnection occurs where oppositely directed field lines come into contact.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A review of thermal physics and management inside lithium-ion batteries for high energy density and fast charging

Traditionally it has been assumed that battery thermal management systems should be designed to maintain the battery temperature around room temperature. That is not always true as lithium-ion battery (LIB) R&D is pivoting towards the development of high energy density and fast charging batteries. Therefore, it is necessary to have a comprehensive review of thermal considerations for LIBs targeted for high energy density and fast charging, i.e., the optimal thermal condition, thermal physics (heat transport and generation) inside the battery, and thermal management strategies. As the energy density and charge rate increases, the optimal battery temperature can shift to be higher than room temperature. In the first part of the review various sources of heat generation inside LIBs and various approaches to minimizing battery heat generation are summarized. The importance of heat of mixing due to ion diffusion during fast charging is also highlighted. To improve the temperature uniformity and avoid excessive internal temperature rise, heat transfer inside the battery needs to be enhanced, and reducing the thermal contact resistance between the electrodes and separator can significantly increase the effective thermal conductivity of batteries. In the second part of the review various challenges and latest developments related to thermal transport and properties of LIBs are discussed. Finally, a summary of latest advancement on smart control of internal temperature of LIBs is discussed as depending on the ambient temperature and the optimal temperature; the battery heat needs to be retained or dissipated to elevate or avoid temperature rise.

25 ENERGY STORAGE↗
DOE OSTI2025

Dissipative ground state preparation in ab initio electronic structure theory

Dissipative engineering is a powerful tool for quantum state preparation, and has drawn significant attention in quantum algorithms and quantum many-body physics in recent years. In this work, we introduce a novel approach using the Lindblad dynamics to efficiently prepare the ground state for general ab initio electronic structure problems on quantum computers, without variational parameters. These problems often involve Hamiltonians that lack geometric locality or sparsity structures, which we address by proposing two generic types of jump operators for the Lindblad dynamics. Type-I jump operators break the particle number symmetry and should be simulated in the Fock space. Type-II jump operators preserves the particle number symmetry and can be simulated more efficiently in the full configuration interaction space. For both types of jump operators, we prove that in a simplified Hartree-Fock framework, the spectral gap of our Lindbladian is lower bounded by a universal constant. For physical observables such as energy and reduced density matrices, the convergence rate of our Lindblad dynamics with Type-I jump operators remains universal, while the convergence rate with Type-II jump operators only depends on coarse grained information such as the number of orbitals and the number of electrons. To validate our approach, we employ a Monte Carlo trajectory-based algorithm for simulating the Lindblad dynamics for full ab initio Hamiltonians, demonstrating its effectiveness on molecular systems amenable to exact wavefunction treatment.

Quantum chemistry↗
DOE OSTI2026

A New Theoretical Framework for Parameterizing Nonequilibrium Fractionation During Evaporation From the Ocean

The evaporation model for water isotopes proposed by Craig and Gordon (1965, https://books.google.co.in/books?id=6wIKAQAAIAAJ) is used in most isotope-enabled atmospheric models for the parameterization of nonequilibrium fractionation during evaporation from the ocean. In this model, one of the most uncertain parameters is the nonequilibrium fractionation factor 𝛼 𝑘 . Many isotope models use the formulation of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029), which parameterizes 𝛼 𝑘 as a function of wind speed and distinguishes between a smooth and a rough regime to account for the effect of ocean waves. The resulting discontinuity in 𝛼 𝑘 between smooth and rough regimes has been disputed by several empirical studies. Here, we present a new approach to parameterizing 𝛼𝑘 by explicitly accounting for the influence of wave drag on the momentum flux near the surface. Following Cifuentes-Lorenzen et al. (2018, https://doi.org/10.1007/s10546-018-0376-0), we add a third wave-induced component to the total momentum flux, in addition to the viscous and turbulent components, and extend the definition of the eddy viscosity to account for the momentum flux due to waves and turbulent dissipation near the surface. The new scheme predicts a slight decrease of 𝛼 𝑘 with wind speed, similar to the smooth-regime parameterization of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029). This new parameterization is incorporated into the isotope-enabled Community Atmosphere Model, where it improves the correlation of simulated and measured vapor deuterium excess relative to the default version and a version with constant 𝛼 𝑘 , suggesting that it may be used as a valid representation of fractionation during evaporation from the ocean in future isotope models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Low-dissipation centred schemes for hyperbolic equations in conservative and non-conservative form

We propose a one-parameter family of low-dissipation centred numerical schemes for solving hyperbolic equations in conservative or non-conservative form, using finite volume or discontinuous Galerkin finite element methods. The new schemes spring out from the multi-dimensional FORCE method and are determined by a single parameter α ≥ 1. Given an increasing sequence of real numbers 1 ≤ α 1 < α 2 < . . . < α K , there corresponds a sequence of numerical schemes with stability restriction associated to a decreasing sequence of Courant numbers 1 > c 1 > c 2 > . . . > c K > 0 and a decreasing sequence of corresponding numerical viscosity functions d 1 > d 2 > . . . > d K . For a given Courant number c k ≤ 1 there is a real number α k > 0 and a corresponding stable scheme with minimal numerical viscosity d k . The proposed schemes suit very well the family of high-order discontinuous Galerkin finite element methods and the recently proposed class of ADER-TR schemes, whose orders of accuracy define decreasing sequences of Courant numbers, as the order of accuracy increases. The centred methods of this paper are stable in 2D and 3D in the frame of simultaneous updating formulae, unlike other centred methods, such as 1D FORCE, which are not. Furthermore, the schemes are highly accurate for slowly-moving waves, which is precisely the kind of waves that traditional centred methods smear disastrously. Additional features of the proposed schemes include ease of implementation and applicability to any hyperbolic system either in conservative or non-conservative form. Here, the proposed schemes are analysed and computationally assessed through a suite of test problems for a linear model system, for the Euler equations in one and two space dimensions, and for the Baer-Nunziato equations for compressible two-phase flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Fluctuation dynamo in a weakly collisional plasma

The turbulent amplification of cosmic magnetic fields depends upon the material properties of the host plasma. In many hot, dilute astrophysical systems, such as the intracluster medium (ICM) of galaxy clusters, the rarity of particle–particle collisions allows departures from local thermodynamic equilibrium. These departures – pressure anisotropies – exert anisotropic viscous stresses on the plasma motions that inhibit their ability to stretch magnetic-field lines. We present an extensive numerical study of the fluctuation dynamo in a weakly collisional plasma using magnetohydrodynamic (MHD) equations endowed with a field-parallel viscous (Braginskii) stress. When the stress is limited to values consistent with a pressure anisotropy regulated by firehose and mirror instabilities, the Braginskii-MHD dynamo largely resembles its MHD counterpart, particularly when the magnetic field is dynamically weak. If instead the parallel viscous stress is left unabated – a situation relevant to recent kinetic simulations of the fluctuation dynamo and, we argue, to the early stages of the dynamo in a magnetized ICM – the dynamo changes its character, amplifying the magnetic field while exhibiting many characteristics reminiscent of the saturated state of the large-Prandtl-number ( ${Pm}\gtrsim {1}$ ) MHD dynamo. We construct an analytic model for the Braginskii-MHD dynamo in this regime, which successfully matches simulated dynamo growth rates and magnetic-energy spectra. A prediction of this model, confirmed by our numerical simulations, is that a Braginskii-MHD plasma without pressure-anisotropy limiters will not support a dynamo if the ratio of perpendicular and parallel viscosities is too small. This ratio reflects the relative allowed rates of field-line stretching and mixing, the latter of which promotes resistive dissipation of the magnetic field. Finally, in all cases that do exhibit a viable dynamo, the generated magnetic field is organized into folds that persist into the saturated state and bias the chaotic flow to acquire a scale-dependent spectral anisotropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

In this work, we employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant–matrix pair correlation function and the penetrant–matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Equipping Neural Network Surrogates with Uncertainty for Propagation in Physical Systems

Coarse-grained or filtered models typically rely on closure models to account for unresolved scales. For instance, large eddy simulation for modeling turbulent fluid flows explicitly resolves the largest scales, but requires modeling closure terms to account for the sub-filter scales. With the vast amount of data available from high-fidelity simulations, there are unique opportunities to leverage data-driven modeling techniques to formulate expressive and flexible closure models. Despite their flexibility, data-driven models struggle in domain shift settings, i.e. when deployed in configurations not captured in the training dataset. In particular, the efficacy of neural network surrogates is difficult to assess a priori due to the deterministic, point-estimate nature of predictions. In high-consequence applications, such models require reliable uncertainty estimates in the data-informed and out-of-distribution regimes. To quantify uncertainties in both regimes, we employ Bayesian neural networks which are able to capture both epistemic and aleatoric uncertainties. We will discuss challenges associated with the training and evaluation of these networks. Furthermore, we will discuss uncertainty embedding strategies to enable efficient sampling and propagation of uncertainty through high-fidelity simulations.

Bayesian neural networks↗
DOE OSTI2025

The Remarkable 2024 North Atlantic Mid‐Season Hurricane Lull

Abstract The 2024 North Atlantic (hereafter Atlantic) hurricane season started quickly, with the earliest Category 5 on record (Beryl) and three hurricanes forming through 14 August. Following Ernesto's dissipation on 20 August, the Atlantic hurricane season became extremely quiet during the climatological peak of hurricane season, with only one Category 2 hurricane (Francine) and one tropical storm through 23 September. Several environmental factors likely contributed to this unexpected, prolonged lull. During mid‐to‐late August, subseasonal conditions were broadly favorable for Atlantic hurricanes, but a northward shift in African easterly wave emergence latitude yielded fewer tropical cyclone seed disturbances that also traversed unfavorably cool ocean water. During early‐to‐mid September, subseasonal variability driven by the Madden‐Julian oscillation was less conducive to hurricane activity, with several bouts of increased vertical wind shear across the central Atlantic. Throughout most of the lull, the tropical Atlantic was anomalously dry and subsident, suppressing hurricane formation chances. Plain Language Summary The 2024 North Atlantic (hereafter Atlantic) hurricane season was anticipated to be extremely active, with seasonal forecasts calling for a potentially historic season. Following a busy start including the earliest Category 5 Atlantic hurricane on record (Beryl), the season became very quiet between 20 August–23 September, with only one Category 2 hurricane (Francine) and one tropical storm forming. Following the lull, the Atlantic became extremely busy again, with seven hurricanes forming from 25 September through the end of the season–the most on record for this period. This paper focuses on the surprising lull and attributes it to several different factors, including a late‐August northward shift in African easterly wave tracks (organized thunderstorm complexes spanning thousands of kilometers that often serve as Atlantic hurricane seeds). This northward shift brought these systems out over colder water and into a drier airmass, suppressing their hurricane formation chances. During the first 3 weeks of September, a large‐scale tropical phenomenon known as the Madden‐Julian oscillation, resulted in increased sinking motion, drier mid‐levels and increased vertical wind shear–all of which typically reduce Atlantic hurricane activity. Overall, the lull demonstrates how seed disturbances and other environmental factors can strongly modulate aggregate seasonal tropical cyclone activity. Key Points Following an active start, the 2024 Atlantic hurricane was quiet during the climatological peak with only 1 tropical storm and 1 hurricane The lull was not anticipated by seasonal forecasting agencies, which called for an extremely active 2024 Atlantic hurricane season Factors reducing storm count included a northward shift in African easterly waves, broad‐scale subsidence, and reduced mid‐level moisture

Klotzbach, P J↗
DOE OSTI2025

DC-Link Capacitor Design for a Neutral-Point-Less Three-Level Dual-Phase Inverter for Traction Application

Conventional multilevel inverters, such as neutral￾point-clamped and T-type inverters, have gained popularity in electric vehicle applications due to their advantages, including high voltage range, high power capability, low switching losses, low total harmonic distortion, and low electromagnetic interference. However, these traditional multilevel inverters require a neutral point connection to generate a zero-voltage vector. The neutral point current oscillates at three times the fundamental frequency, leading to voltage imbalance and overvoltage stress on the power modules. Additionally, the use of two stacked DC-link capacitors increases the volume required for the same overall capacitance and complicates packaging due to ripple current and heat dissipation from separate components. This is a significant concern for traction drive units, where space is limited. In this paper, a neutral-point-less multilevel dual three￾phase inverter topology is investigated for traction inverter applications. Simulation results demonstrate that the proposed topology effectively retains the benefits of multilevel operation while utilizing a single DC-link capacitor. The inverter model was simulated in conjunction with an industry-standard battery model to evaluate the potential for capacitor size reduction compared to a conventional three-level inverter.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Tungsten–carbon surface evolution and erosion modeling for a small angle slot divertor in DIII-D

Abstract We modeled tungsten–carbon mixed surface evolution, sputtering erosion, and transport for the tungsten coated region of a small angle slot (SAS) divertor design for the DIII-D tokamak. This divertor concept aims to achieve a closed slot dissipative plasma to minimize heat load and surface erosion, and to study high- Z material performance. Our advanced simulations use coupled ITMC-DYN material mixing/response and 3D full kinetic REDEP/WBC erosion/redeposition code packages, with divertor plasma solution from the SOLPS-ITER package with 4 MW power input. The SAS design geometry and resulting in-slot plasma parameters cause significant differences in predicted sputter/transport from a conventional divertor. For 2% C/D incident plasma ratio, an equilibrium mixed C/W surface is attained at ∼30 s of discharge, from wall sputtered carbon transported to the 10 cm long tungsten divertor region. Tungsten remains exposed to the plasma, but the evolved surface composition varies with different C/D ratios. Tungsten is primarily sputtered from the mixed surface by impinging carbon ions in the +1 to +4 charge states, with some self-sputtering. Redeposition of sputtered tungsten to the divertor is significant, ∼80% along the higher plasma temperature attached plasma SAS entrance region, but this is less than the typically near-unity values for a conventional divertor. Plasma-incident carbon is highly backscattered (∼50%) from the mixed surface, with little redeposition (<10%); this helps maintain tungsten in the surface sputter zone. Carbon is mainly sputtered from the mixed surface by D + ions, also with low redeposition (∼10%–30%). Finally, the modeling shows non-zero but low sputtered tungsten current from the divertor to the core plasma direction. These results appear favorable for effective testing of a tungsten-containing SAS divertor in DIII-D, and extrapolation of mixed-material evolution/response findings to the analogous low- Z /high- Z , Be/W, ITER plasma facing system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Langevin dynamics of generalized spins as SU( N ) coherent states

Classical models of spin systems traditionally retain only the dipole moments, but a quantum spin state will frequently have additional structure. Spins of magnitude S have N = 2⁢S + 1 levels. Alternatively, the spin state is fully characterized by a set of N 2 - 1 local physical observables, which we interpret as generalized spin components. For example, a spin with S = 1 has three dipole components and five quadrupole components. Further, these components evolve under a generalization of the classical Landau-Lifshitz dynamics, which can be extended with noise and damping terms. In this paper, we reformulate the dynamical equations of motion as a Langevin dynamics of SU(N) coherent states in the Schrödinger picture. This viewpoint is especially useful as the basis for an efficient numerical method to sample spin configurations in thermal equilibrium and to simulate the relaxation and driven motion of topological solitons. To illustrate the approach, we simulate a nonequilibrium relaxation process that creates CP 2 skyrmions, which are topological defects with both dipole and quadrupole character.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Imaging a terahertz superfluid plasmon in a two-dimensional superconductor

The superconducting gap defines the fundamental energy scale for the emergence of dissipationless transport and collective phenomena in a superconductor. In layered high-temperature cuprate superconductors, in which the Cooper pairs are confined to weakly coupled two-dimensional (2D) copper–oxygen (CuO 2 ) planes, terahertz (THz) spectroscopy at subgap millielectronvolt (meV) energies has provided crucial insights into the collective superfluid response perpendicular to the superconducting layers. However, within the CuO 2 planes, the collective superfluid response manifests as plasmonic charge oscillations at energies far exceeding the superconducting gap, obscured by strong dissipation. Here, in this study, we present spectroscopic evidence of a below-gap, 2D superfluid plasmon in few-layer Bi 2 Sr 2 CaCu 2 O 8+x and spatially resolve its deeply subdiffractive THz electrodynamics. By placing the superconductor in the near field of a spintronic THz emitter, we reveal this distinct resonance—absent in bulk samples and observed only in the superconducting phase—and determine its plasmonic nature by mapping the geometric anisotropy and dispersion. Crucially, these measurements offer a direct view of the momentum-dependent and frequency-dependent superconducting transition in two dimensions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗