Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PPPL”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Preliminary proof-of-concept of real-time divertor heat flux control from infrared cameras with nitrogen injection in the DIII-D tokamak

In future tokamak reactors like ITER and the Fusion Pilot Plant (FPP), real-time feedback control of heat flux to the plasma-facing components (PFC) will be critical for steady-state operation. This work presents the first experimental demonstration of real-time divertor heat flux estimation with infrared thermography and feedback control with impurity seeding on the DIII-D tokamak. The flexible infrastructure of the Plasma Control System (PCS) on DIII-D makes this new capability possible. The PCS software runs on a gateway computer system, and five real-time compute nodes. An array of low latency streaming digitizers from D-TACQ Solutions connects to these real-time computers to collect and process data, and send commands to actuators during plasma discharges. This system handles the signal IO from the tokamak and allows the PCS to utilize the diagnostic data necessary to perform control in real-time. Feedback control on heat flux was accomplished by feeding infrared camera data from the “Infrared TV” (IRTV) camera to a custom-developed User Datagram Protocol (UDP) server. This server transmits infrared data to a newly developed PCS algorithm that estimates the heat flux to PFC. Here, a proportional integral derivative (PID) controller minimizes the error between a heat flux reference and the real-time estimate by injecting nitrogen gas into the divertor.

Algorithm development↗
DOE OSTI2025

A terminology for scientific workflow systems

The term “scientific workflow” has evolved over the last two decades to encompass a broad range of compositions of interdependent compute tasks and data movements. It has also become an umbrella term for processing in modern scientific applications. Today, many scientific applications can be considered as workflows made of multiple dependent steps, and hundreds of workflow systems have been developed to manage and run these scientific workflows. However, no turnkey solution has emerged from the field to address the diversity of scientific processes and the infrastructure on which they are supposed to be implemented. Instead, new research problems requiring the execution of scientific workflows with some novel feature often lead to the development of an entirely new workflow system. A direct consequence of this situation is that many existing workflow management systems (WMSs) share some salient features, offer similar functionalities, and can manage the same categories of workflows but at the same time also have some distinct capabilities that can be important for specific applications. This situation makes researchers who develop workflows face the complex question of selecting a WMS. This selection can be driven by technical considerations, to find the system that is the most appropriate for their application and for the computing and storage resources available to them, or other factors such as reputation, adoption, strong community support, or long-term sustainability. To address this problem, a group of WMS developers and practitioners joined their efforts to produce a community-based terminology of WMSs. This paper summarizes their findings and introduces this new terminology to characterize WMSs. Furthermore, this terminology is composed of fives axes: workflow structure and characteristics, composition, orchestration, data management, and metadata capture. Each axis comprises several concepts that capture the prominent features of WMSs. Based on this terminology, this paper also presents a classification of 23 existing WMSs according to the proposed axes and terms.

Community-based terminology↗
DOE OSTI2025

Thermovelocimetric characterization of liquid metal convection in a rotating slender cylinder

Rotating turbulent convection occurs ubiquitously in natural convective systems encompassing planetary cores, oceans and atmospheres, as well as in many industrial applications. While the global heat and mass transfer of water-like rotating Rayleigh–Bénard convection is well-documented, the dynamics in low-Prandtl-number liquid metals remain less understood. In this study, we experimentally investigate rotating Rayleigh-Bénard convection in liquid gallium (Prandtl number Pr ≈ 0.027) within a slender cylinder (diameter-to-height aspect ratio Γ = D/H =1/2) using novel thermovelocimetric diagnostic techniques that integrate simultaneous multi-point thermometry and ultrasonic Doppler velocity measurements. Our results reveal the formation of a stable, global-scale azimuthal wavenumber m = 2 quadrupolar vortex at low supercriticality. We propose that enhanced wall modes facilitated by the slender cylinder geometry interact with the bulk flow to create these large-scale axialized vortices. Furthermore, our findings imply a distinct scaling behavior for the wall-mode precession frequency in liquid metals, extending previous results obtained for moderate-Pr fluids. In conclusion, this provides new insights into wall-bulk coupling mechanisms of low-Pr rotating convective turbulence.

Magnetohydrodynamics↗
DOE OSTI2026

Removal of deuterium retained in boron powder by oxygen or high-temperature bakeout

Fuel retention in dust accumulated in fusion devices is a potential operational hazard. Boronization is one of the leading wall-conditioning candidates for future fusion devices. Fuel retention characteristics of boron dust must be studied to accurately understand the potential hazards of boron dust in fusion devices. We evaluated the retention of deuterium (D) in commercially available boron powder as a proxy for tritium retention in boron dust. Diffuse reflectance infrared Fourier transform spectroscopy revealed that boron powder between 300 and 500 K exposed to neutral hydrogen (H 2 ) gas retains H via B-H bonding. We utilized D 2 exposures of 17,400 L on the boron powder, which led to retention on the order of 0.5 g D per 100 kg B or 1.2 g D per m 2 B, as estimated from temperature programmed desorption measurements. We evaluated how powder bakeouts under ultra-high vacuum (UHV), or in an oxygen (O 2 ) or H 2 gas backfilled environment can remove the retained fuel. Bakeouts below 523 K under UHV exhibited poor D removal, promoting recapture of D 2 present as background gas in the chamber. Higher temperature bakeouts led to the removal of over 80% of the retained D within 3 days (623 K) or 1 h (723 K). Bakeouts at 423 K under an O 2 gas environment removed ∼150% more retained D from the boron powder, when compared to bakeouts conducted under UHV. These results suggest that efficient fuel removal from boron dust can be achieved by high-temperature bakeout (>623 K) or by bakeouts at lower temperature (423 K) under an O 2 environment.

Bakeout↗
DOE OSTI2025

Tungsten erosion during L-mode discharges in the DIII-D SAS-VW divertor

A series of L-mode plasma discharges was performed in the DIII-D tokamak to assess the impact of outer strike point (OSP) position and toroidal magnetic field direction on erosion and core contamination potential of the recently-installed, tungsten-coated Small Angle Slot (SAS-VW) divertor. In one discharge, in-slot emission spectroscopy measured an<48 % increase in the W gross erosion rate when the OSP was moved 3 cm outwards, away from the V-shaped vertex of the slot divertor. However, the effective W yield (erosion rate divided by the incident D flux) was, overall, insensitive to changes in OSP location. Consistently low estimates of the effective W yield based on measurements taken a few cm outwards from the vertex suggest potentially significant C surface contamination. No W emission signal was detected when orienting the toroidal magnetic field such that the ion B×∇B drift direction is pointed away from the X-point. However, measurements of W content in the plasma core for both toroidal magnetic field directions suggest the presence of additional, unmeasured sources of erosion. The difference in the measured core W density with OSP position is much greater than the difference in the measured erosion rates, which may suggest that the leakage of eroded impurities out of the divertor is governed primarily through the parallel ion temperature gradient and friction forces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Experiment-modeling studies comparing energy dissipation in the DIII-D SAS and SAS-VW divertors

Recent DIII-D experiments on Small Angle Slot (SAS) divertors have confirmed that a combination of divertor closure and target shaping can enhance cooling across the divertor target and increase energy dissipation, but with significant dependence on B T (toroidal magnetic field) direction. In these novel divertors, the roles of closure, target shaping, drifts, and scale lengths are all interconnected in optimizing dissipation, with the separatrix electron density n eSEP being the key parameter associated with the level of dissipation/detachment. After modifying the original flat-targeted graphite SAS to include a V shape with a tungsten coating on the outer side of the divertor (SAS-VW), matched series of discharges were run to compare to detailed SOLPS-ITER modeling. Experimentally, when run as designed with the outer strike point at the slot vertex, SAS-VW requires nearly identical n eSEP for detachment as the original SAS, with little difference in dissipation for the new geometry. This is in contrast to (1) earlier modeling predictions that a small change of the SAS geometry to a V shape should enhance dissipation at the same n eSEP for magnetic configurations having better H-mode access (ion B × ∇B drift directed into the divertor), and (2) despite the achievement of significantly higher (2-7x) neutral pressures and compression in the SAS-VW slot. Comparisons of experimental density scans to the most recent SOLPS-ITER modeling with ExB drifts show reasonable agreement for dissipation/detachment onset when using separatrix density as the independent parameter. In order to help understand the discrepancy in modeled vs actual performance for the new configuration, additional measurements varying gas injection location and impurity injection were undertaken. In-slot D 2 gas fueling is more effective (5–22 %) in promoting detachment, in accord with modeling. In-slot impurity injection (N 2 or Ne) can yield 30 % lower core Z eff and 15 % less confinement degradation after detachment compared to main chamber puffing, as well as relatively lower tungsten leakage from the divertor. Modeling can also reproduce the improved detachment seen as the strike point moves inboard of the slot vertex. While we can explain the effects of the most important parameters causing energy dissipation in these slot divertors, it remains that many aspects of their behavior cannot be accurately modeled using state-of-art codes such as SOLPS-ITER. This is of concern for future model-driven designs utilizing similar V-shaped geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Recollections for the 50th anniversary of the plasma surface interactions (PSI) in controlled fusion devices conference

The Plasma Surface Interactions in Controlled Fusion Devices (PSI) conference reached an important milestone in 2024 with its 50th anniversary. It was celebrated at its venue in Marseille by a special round table discussion gathering 6 former chairmen of its Programme Committees, who gave some highlights presented at the conference during the five decades. The article provides a summary of this overview.

Plasma Surface Interactions Conference↗
DOE OSTI2026

Deuterium retention characteristics during lithium powder and granule injection in EAST

Fuel retention remains a critical challenge for magnetic confinement fusion devices. As a low-Z material, lithium plays a key role in tokamak wall conditioning and impurity control. In this study, the fuel retention behaviors associated with two lithium forms—powder and granule—are systematically examined using the gas balance method. Both forms of lithium injection significantly enhance fuel retention, shifting the wall behavior from net outgassing to net absorption. For lithium powder injection, compared with the reference discharge, a suppression efficiency exceeding 40% is achieved relative to the number of injected lithium atoms. Moreover, discharges with different injection rates show a monotonic increase in retained fuel with the injected lithium quantity. In contrast, repeated granule injections exhibit a pronounced cumulative effect, with the retention ratio varying from 0.16 to 0.65 over successive discharges, indicating progressive enhancement of wall absorption and deuterium retention. These results advance the understanding of wall behavior modification induced by solid material injection and provide insights for developing effective wall-conditioning strategies using low-Z materials in ITER and future fusion reactors, while the expected high retention may ultimately exclude the use of Li wall conditioning in future fusion devices.

Fuel retention↗
DOE OSTI2026

Influence of preheating on tungsten damage under pulsed hydrogen plasma exposure

Preheated tungsten plate behavior under multiple pulsed hydrogen plasma irradiations in the PW-7 plasma accelerator was investigated. The preheating temperature was 973 K, which is above the ductile to brittle transition temperature (DBTT ≈ 673 K). The plate was exposed to three consecutive series of 13 plasma pulses each. Preheating sufficiently suppressed surface cracking at the early irradiation stage: after 13 pulses, cracks were rare, localized, and did not form extensive fractures. During subsequent irradiation series, the area fraction of fine cracks increased by ∼ 2% due to accumulation of dislocations around the pre-existing defects. However, their width remained nearly constant, averaging 0.2–0.3 μm, with rare cracks up to 1 μm in size. Repeated plasma exposure resulted in pronounced surface roughening associated with multiple melting cycles, crater formation and blister rupture, with roughness reaching ∼ 204 nm and a recorded sample mass loss of ∼ 2.8 mg. These results indicate that controlling tungsten preheating well above DBTT effectively reduces damage. However, the surface condition during initial irradiation significantly influences its subsequent morphological evolution. Preheating to 773 K, selected regarding the minimum DBTT value, has previously been demonstrated to be an effective approach for limiting crack growth. A further increase in temperature of approximately ∼ 1000 K did not qualitatively improve results: crack still formed and widths did not decrease. At the same time, the absence of further crack coarsening indicates that the limiting effect of preheating has been reached. Considering the increased energy costs, further increasing the preheating temperature seems unjustified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Geostrophic to magnetically-damped turbulence in liquid metal rotating magnetoconvection

Understanding planetary core convection dynamics requires the study of convective flows in which the Coriolis and Lorentz forces attain a leading-order, so-called magnetostrophic balance. Experimental investigations of rotating magnetoconvection (RMC) in the magnetostrophic regime are therefore essential to broadly characterize the properties of local-scale planetary core flow. Toward this end, we present here the first thermovelocimetric measurements of magnetostrophic, liquid metal convection, which are made using liquid gallium as the working fluid, at moderate rotation rates (Ekman numbers 10 -4 ≤ Ek ≤ 10 -5 ) and in the presence of dynamically strong magnetic fields (Elsasser number Λ = 1). Complementary rotating convection (RC) experiments are performed at the same rotation rates to serve as reference cases. Our RMC velocity measurements adequately follow a geostrophic turbulent scaling for cases in which local-scale convective inertial forces exceed the Lorentz forces in the fluid bulk. In cases where Lorentz forces exceed local-scale inertia (local or convective interaction parameter N ℓ or N C ≳ 3), the root-mean-square RMC velocities are magnetically damped, yielding values below the geostrophic turbulent RC scaling prediction. An enhancement in heat transfer is observed, which we attribute to the increased coherence of vertically aligned magnetostrophic convective flow. Finally, our results and comparisons to dynamo simulations suggest that the system is in the strong-field regime, where convection-scale flows are magnetically damped with N C ≳ 3).

Convective turbulence↗
DOE OSTI2025

Methanol adsorption and dissociation on GaP(110) studied by ambient pressure X-ray photoelectron spectroscopy

Ambient pressure X-ray photoelectron spectroscopy (AP-XPS) was used to investigate methanol (CH 3 OH) adsorption and reaction on the GaP(110) surface. Exposure of CH 3 OH to GaP(110) at room temperature led to the formation of at least four different surface species as indicated by analysis of C 1s and O 1s XPS features. By combining AP-XPS data with density functional theory calculations, the surface species were identified as methoxy (CH 3 O*), formaldehyde (CH 2 O*), and paired methanol (p-CH 3 O*H) and methoxy (p-CH 3 O*) species, where “paired” means that they belong to a hydrogen-bonded methoxy-methanol complex. Asterisk * here indicates an adsite. The formation of CH 2 O* via the dehydrogenation of CH 3 O* was shown to be limited by the availability of vacant phosphorus (P) sites on GaP(110). With an increase in CH 3 OH pressure, the fractional coverage of CH 3 O* species reached 0.55, and the surface P sites were completely saturated with hydrogen. Under a constant CH 3 OH pressure of 0.5 Torr, the surface concentration of the paired species and of CH 2 O* remained constant until 400 K. At higher temperatures, thermally driven reactions led to a significant increase in the concentration of surface CH x * species, which suggests that C-O bond cleavage of the CH 3 O group is the dominant decomposition mechanism on GaP(110). In conclusion, based on the reactivity of GaP(110) toward CH 3 OH dehydrogenation, elevated temperatures and CH 3 OH pressures may be used to functionalize this surface.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effect of insulator end cap thickness on time-dependent Hartmann flow in a rotating mirror

We present a framework for analysing plasma flow in a rotating mirror. By making a series of physical assumptions, we reduce the magnetohydrodynamic (MHD) equations in a three-dimensional cylindrical system to a one-dimensional system in a shallow, cuboidal channel within a transverse magnetic field, similar to the Hartmann flow in ducts. We then solve the system both numerically and analytically for a range of values of the Hartmann number and calculate the dependence of the plasma flow speed on the thickness of the insulating end cap. We observe that the mean flow overshoots and decelerates before achieving a steady-state value, a phenomenon that the analytical model cannot capture. This overshoot is directly proportional to the thickness of the insulating end cap and the external electric field, with a weak dependence on the external magnetic field. Our simplified model can act as a benchmark for future simulations of the supersonic mirror device CMFX (centrifugal magnetic fusion experiment), which will employ more sophisticated physics and realistic magnetic field geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A quasi-linear model of electromagnetic turbulent transport and its application to flux-driven transport predictions for STEP

A quasi-linear reduced transport model is developed from a database of high-β electromagnetic nonlinear gyrokinetic simulations performed with spherical tokamak for energy production (STEP) relevant parameters. The quasi-linear model is fully electromagnetic and accounts for the effect of equilibrium flow shear using a novel approach. Its flux predictions are shown to agree quantitatively with predictions from local nonlinear gyrokinetic simulations across a broad range of STEP-relevant local equilibria. This reduced transport model is implemented in the T3D transport solver that is used to perform the first flux-driven simulations for STEP to account for transport from hybrid kinetic ballooning mode turbulence, which dominates over a wide region of the core plasma. Nonlinear gyrokinetic simulations of the final transport steady state from T3D return turbulent fluxes that are consistent with the reduced model, indicating that the quasi-linear model may also be appropriate for describing the transport steady state. Within the assumption considered here, our simulations support the existence of a transport steady state in STEP with a fusion power comparable to that in the burning flat top of the conceptual design, but do not demonstrate how this state can be accessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Constructing nested coordinates inside strongly shaped toroids using an action principle

A new approach for constructing polar-like boundary-conforming coordinates inside a toroid with strongly shaped cross-sections is presented. A coordinate mapping is obtained through a variational approach, which involves identifying extremal points of a proposed action in the mapping space from [0,2π] 2 ×[0,1] to a toroidal domain in $\mathbb{R}$ 3 . This approach employs an action built on the squared Jacobian and radial length. Extensive testing is conducted on general toroidal boundaries using a global Fourier–Zernike basis via action minimisation. The results demonstrate successful coordinate construction capable of accurately describing strongly shaped toroidal domains. The coordinate construction is successfully applied to the computation of three-dimensional magnetohydrodynamic equilibria in the GVEC code where the use of traditional coordinate construction by interpolation from the boundary failed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Metastability of stratified magnetohydrostatic equilibria and their relaxation

Motivated by explosive releases of energy in fusion, space and astrophysical plasmas, we consider the nonlinear stability of stratified magnetohydrodynamic equilibria against two-dimensional interchanges of straight magnetic-flux tubes. We demonstrate that, even within this restricted class of dynamics, the linear stability of an equilibrium does not guarantee its nonlinear stability: equilibria can be metastable. We show that the minimum-energy state accessible to a metastable equilibrium under non-diffusive two-dimensional dynamics can be found by solving a combinatorial optimisation problem. These minimum-energy states are, to good approximation, the final states reached by our simulations of destabilised metastable equilibria for which turbulent mixing is suppressed by viscosity. To predict the result of fully turbulent relaxation, we construct a statistical mechanical theory based on the maximisation of Boltzmann's mixing entropy. This theory is analogous to the Lynden-Bell statistical mechanics of collisionless stellar systems and plasma, and to the Robert–Sommeria–Miller theory of two-dimensional vortex turbulence. Our theory reproduces well the results of our numerical simulations for sufficiently large perturbations to the metastable equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

MagNetUS: a magnetized plasma research ecosystem

MagNetUS is a network of scientists and research groups that coordinates and advocates for fundamental magnetized plasma research in the USA. Its primary goal is to bring together a broad community of researchers and the experimental and numerical tools they use in order to facilitate the sharing of ideas, resources and common tasks. Discussed here are the motivation and goals for this network and details of its formation, history and structure. An overview of associated experimental facilities and numerical projects is provided, along with examples of scientific topics investigated therein. Finally, a vision for the future of the organization is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Beam focusing and consequences for Doppler backscattering measurements

The phenomenon of focusing of microwave beams in a plasma near a turning-point caustic is discussed by exploiting the analytical solution to the Gaussian beam-tracing equations in the two-dimensional (2-D) linear-layer problem. The location of maximum beam focusing and the beam width at that location are studied in terms of the beam initial conditions. This focusing must be taken into account to interpret Doppler backscattering (DBS) measurements. We find that the filter function that characterises the scattering intensity contribution along the beam path through the plasma is inversely proportional to the beam width, predicting enhanced scattering from the beam focusing region. We show that the DBS signal enhancement for decreasing incident angles between the beam path and the density gradient is due to beam focusing and not due to forward scattering, as was originally proposed by (Gusakov et al., (Plasma Phys. Contr. Fusion, vol. 56, 2014, p. 0250092014, 2017); Plasma Phys. Rep. vol. 43(6), 2017, pp. 605–613). The analytic beam model is used to predict the measurement of the k y density-fluctuation wavenumber power spectrum via DBS, showing that, in an NSTX-inspired example, the spectral exponent of the turbulent, intermediate-to-high k y density-fluctuation spectrum might be quantitatively measurable via DBS, but not the spectral peak corresponding to the driving scale of the turbulent cascade.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Suppression of the collisionless tearing mode by flow shear: implications for reconnection onset in the Alfvénic solar wind

We analyse the collisionless tearing mode instability of a current sheet with a strong shear flow across the layer. The growth rate decreases with increasing shear flow, and is completely stabilised as the shear flow becomes Alfvénic. We also show that, in the presence of strong flow shear, the tearing mode growth rate decreases with increasing background ion-to-electron temperature ratio, the opposite behaviour to the tearing mode without flow shear. We find that even a relatively small flow shear is enough to dramatically alter the scaling behaviour of the mode, because the growth rate is small compared with the shear flow across the ion scales (but large compared with shear flow across the electron scales). Our results may explain the relative absence of reconnection events in the near-Sun Alfvénic solar wind observed recently by NASA’s Parker Solar Probe.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗