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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Measurement of Zero-Frequency Fluctuations Generated by Coupling between Alfvén Modes in the JET Tokamak

We report the first experimental detection of a zero-frequency fluctuation that is pumped by an Alfvén mode in a magnetically confined plasma. Core-localized Alfvén modes of frequency inside the toroidicity-induced gap (and its harmonics) exhibit three-wave coupling interactions with a zero-frequency fluctuation. The observation of the zero-frequency fluctuation is consistent with theoretical and numerical predictions of zonal modes pumped by Alfvén modes, and is correlated with an increase in the deep core ion temperature, temperature gradient, confinement factor H 89,P , and a reduction in the main ion heat diffusivity. Despite the energetic particle transport induced by the Alfvén eigenmodes, the generation of a zero-frequency fluctuation that can suppress the turbulence leads to an overall improvement of confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Electron-Only Magnetic Reconnection and Inverse Magnetic-Energy Transfer at Subion Scales

We derive, and validate numerically, an analytical model for electron-only magnetic reconnection applicable to strongly magnetized plasmas. Our model predicts subion-scale reconnection rates significantly higher than those pertaining to large-scale reconnection, aligning with recent observations and simulations. Here, we apply this reconnection model to the problem of inverse magnetic energy transfer at subion scales. We derive time-dependent scaling laws for the magnetic energy decay and the typical magnetic structure dimensions that differ from those previously found in the magnetohydrodynamics regime. These scaling laws are validated via two- and three-dimensional simulations, demonstrating that subion-scale magnetic fields can reach large, system-size scales via successive coalescence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Weyl semimetal phases and intrinsic spin-Hall conductivity in SbAs ordered alloys

Here, using density functional theory calculations, we investigated possible Weyl semimetal (WSM) phases in antimony arsenide ordered alloys Sb 1-x ⁢As x (x=0, 1/6, 1/3, 1/2, 2/3, 5/6, 1). We find WSM phases for all As compositions of Sb 1-x ⁢As x with broken inversion symmetry, in contrast to Bi 1-x⁢ Sb x where only compositions x=1/2 and 5/6 were predicted to exhibit WSM phases. The WSM phases in Sb 1-x ⁢As x are characterized by the presence of 12 Weyl points, located within 55 meV from the Fermi level in the case of x = 1/2. The robust spin-orbit coupling strength and Berry curvature in these alloys produce large spin-Hall conductivity in the range of 176–602 ($\hslash$/e)(S/cm), comparable to that in the BiSb alloys. Finally, Sb 0.5 ⁢As 0.5 is predicted to be almost lattice matched to GaAs(111), with the Fermi level within the gap of the semiconductor, facilitating growth and characterization, and thus, offering promising integration with conventional semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structural inversion asymmetry in epitaxial ultrathin films of Bi(111)/InSb(111)B

Bismuth (Bi) films hold potential for spintronic devices due to strong spin-orbit coupling. Understanding the growth, surface states, and interactions with the substrate is key to their functionalization. Large-area high-quality (111) Bi ultrathin films were grown on InSb (111)B substrates by molecular beam epitaxy (MBE). Strong film-substrate interactions epitaxially stabilize the (111) orientation and lead to nonequivalent interface potentials. Analysis of angle-resolved photoemission spectroscopy (ARPES) measurements, employed to characterize the evolution of the surface states with film thickness, indicate a crossing at the $\overline{M}$ point, suggesting a topologically trivial phase in the thin film. In conclusion, the results show the presence of interfacial bonds to the substrate breaks inversion symmetry, preventing the semimetal-to-semiconductor transition predicted for freestanding bismuth layers, highlighting the importance of controlled functionalization and surface passivation of two-dimensional materials.

Inbar, Hadass S. [Univ. of California, Santa Barba↗
DOE OSTI2026

Tokamak level performance in the optimized stellarator Wendelstein 7-X with stable peaked density profiles

Progress towards achieving net gain from nuclear fusion is typically characterized by the triple product 𝑛⁢𝑇⁢𝜏. Although stellarators have several benefits for a reactor, such as the inherent steady-state capability and the lack of disruptions, the maximum achieved triple product has historically lagged behind that of tokamaks. In the large optimized stellarator Wendelstein 7-X, a triple product of (1.10 ± 0.15) × 10 20 m −3 ⁢ keVs has recently been achieved and held stable for 1.9 s. This marks an important step in stellarator core plasma performance and puts it in line with 𝐻-mode performance in similar sized tokamaks. The progress in Wendelstein 7-X performance is due to the suppression of ion-temperature-gradient turbulence by strong density gradients. In the record case, these are maintained by a fine balance of the core particle transport and central fueling using neutral beam injection. The duration of the scenario was limited in time only by technical restrictions of the neutral beam injectors and control of the divertor heat loads. This paper presents the result in the context of the leading magnetic confinement devices along with details of the scenario.

Bannmann, Sebastian [Max Planck Institute for Plas↗
DOE OSTI2026

Understanding Large-Scale Dynamos in Unstratified Rotating Shear Flows

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

accretion disk & black-hole plasma↗
DOE OSTI2025

Forced 3D Reconnection in an Exponentially Separating Magnetic Field

We present a solvable scenario for 3D reconnection in a sheared magnetic field. We consider a localized external force that is applied slowly to a flux tube and then maintained, generating an Alfvénic perturbation that spreads along the field lines. Separation of the sheared field lines reduces the scale of the perturbation across the field, enhancing magnetic diffusion. For a fusion-motivated equilibrium with exponential field-line separation, we find a reconnection timescale proportional to 𝒮/ ln⁡ 𝒮 under magnetohydrodynamics (MHD) and to 𝒮 1/3 for semicollisional electron-only reconnection, where 𝒮 is the Lundquist number of the perturbed flux tube. We generalize these results to arbitrary magnetic geometries, showing that the semicollisional case is geometry independent. Interestingly, we find that slower field-line separation yields an increased reconnection rate in MHD.

Magnetic confinement fusion↗
DOE OSTI2025

Nonideal stability analysis of differentially rotating plasmas with global curvature effects

The linear stability of global nonaxisymmetric modes in differentially rotating, magnetized, nonideal plasma is critical to classifying turbulence and transport phenomena. We investigate the competition between the local magneto-rotational instability (MRI) and the magneto-curvature instability (MCI)—a distinct nonaxisymmetric low-frequency curvature-driven global branch that appears alongside MRI. Here, to accomplish this, we developed a nonideal global spectral method, which is validated against NIMROD code simulations. This spectral approach allows for the direct derivation of an extended effective potential formalism and a resistive Alfvénic resonance condition, providing a framework for direct analysis of energy contributions and confinement mechanisms. Our study reveals that the global, low-frequency MCI persists at low magnetic Reynolds numbers (Rm), whereas the localized, high-frequency MRI is stabilized by diffusive broadening of its structure around its Alfvénic resonances. Consequently, we identify the global MCI branch as the primary onset mechanism for nonaxisymmetric magnetohydrodynamic instability in systems with finite curvature, e.g., astrophysical rotators. We establish distinct parameter regimes for mode dominance: MCI prevails in geometrically moderate-thickness disks with intermediate curvature and radial gaps, while MRI dominates in thin, low-curvature disks with large radial gaps. Mode competition is also highly sensitive to the flow profile, particularly vorticity and its gradient, with nonuniform shear profiles exhibiting more robust instability due to flow curvature (i.e., the second derivative of the flow profile) and shear contributions. A key outcome is the development of spectral diagrams derived from the global spectral method. These diagrams comprehensively map dominant instabilities and their characteristics, offering a predictive tool for critical onset parameters (i.e., flow curvature, magnetic field, and Rm) and facilitating the interpretation of experimental and simulation results. Notably, these diagrams demonstrate that the global MCI is generally the sole unstable mode at the initial onset of nonaxisymmetric instability.

Haywood, Alexander [Princeton Univ., NJ (United St↗
DOE OSTI2025

Piecewise omnigenous stellarators with zero bootstrap current

Until now, quasi-isodynamic magnetic fields have been the only known stellarator configurations that, at low collisionality, give small radial neoclassical transport and zero bootstrap current for arbitrary plasma profiles, the latter facilitating control of the magnetic configuration. The recently introduced notion of piecewise omnigenous fields has enormously broadened the space of stellarator configurations with small radial neoclassical transport. Here, in this Letter, the existence of piecewise omnigenous fields that give zero bootstrap current is proven analytically and confirmed numerically. These results establish piecewise omnigenity as an alternative approach to stellarator reactor design.

Calvo, Iván [Research Centre for Energy, Environme↗
DOE OSTI2026

Band-gap reduction and band alignments of dilute bismide III–V alloys

Adding a few atomic percent of Bi to III–V semiconductors leads to significant changes in their electronic structure and optical properties. Bismuth substitution on the pnictogen site leads to a large increase in spin-orbit splitting Δ SO at the top of the valence band (Γ 8⁢𝑣 −Γ 7⁢𝑣 ) and a large reduction in the band gap, creating unique opportunities in semiconductor device applications. Quantifying these changes is key to the design and simulation of electronic and optoelectronic devices. Using hybrid functional calculations, we predict the band gap of III–Vs (III = Al, Ga, In and V = As, Sb) with low concentrations of Bi (3.125% and 6.25%), the effects of adding Bi on the valence- and conduction-band edges, and the band offset between these dilute alloys and their III–V parent compounds. As expected, adding Bi raises the valence-band maximum (VBM). However, contrary to previous assumptions, the conduction-band minimum (CBM) is also significantly lowered, and both effects contribute to the sizable band-gap reduction. Changes in band gap and Δ SO are notably larger in the arsenides than in the antimonides. In conclusion, we also predict cases of band-gap inversion (Γ 6⁢𝑐 below Γ 8⁢𝑣 ), and Δ SO larger than the band gap, which are key parameters for designing topological materials and for minimizing losses due to Auger recombination in infrared lasers.

Saboor, Abdul [Univ. of Delaware, Newark, DE (Unit↗
DOE OSTI2026

Commissioning of the Power Supplies and Coils of the SMART Tokamak

The small aspect ratio tokamak (SMART) is a spherical tokamak (ST) that offers unique capabilities for studying the potential of negative triangularity. It has been designed, constructed, and is currently being operated by the Plasma Science Fusion Technology (PSFT) Laboratory at the University of Seville. SMART has a total of 21 coils, organized into seven independent circuits and driven by five modular power supplies (PS). The PS operation relies on switching converter technology based on IGBT and supercapacitors (SCs). The PS delivers predefined current waveforms to the copper coil system consisting of the central solenoid (CS), 12 toroidal field (TF) coils, three series pairs of poloidal field (PF) coils, and two independent PF coils. This study details the commissioning and validation of the PS toroidal and solenoid coils, as well as the assembly of coils in SMART. The maximum rated current and slope were measured, along with the series impedance of the coils, the output current ripple, and the noise levels. The internal parameters of the SCs were measured, and optimized current profiles were proposed to enhance overall performance. A comparison has been made between the theoretical values and the experimental results, providing insight into the performance of the PS and areas for improvement.

Power electronics↗
DOE OSTI2026

Spectroscopic Characterization of Plasmoid Properties During Pellet Fueling in W7-X

This study utilizes a spectroscopic approach to investigate the properties of plasmoids that are formed during the process of cryogenic hydrogen pellet fueling in the Wendelstein 7-X (W7-X) stellarator. An analysis of the Balmer series emissions was conducted using a diagnostic that was installed during the 2024 operational campaign. Electron temperature, density, and plasma beta ( β ) values can be inferred from the emissions of radiation from the ablation plasmoid. These values are essential for validating pellet ablation models and, in the future, optimizing fueling strategies in steady-state fusion devices.

Cryogenic hydrogen pellets↗
DOE OSTI2026

Multiphysics Analysis of Li Cooled Divertor Substrate During Loss of Coolant Accident (LOCA)

In the ongoing study of potential designs for liquid metal (LM) plasma-facing components (PFCs), so-called “slow” and “fast” Li flow divertor concepts are under investigation. In the previous studies on design and analysis of the slow Li flow divertor and comparison with the fast Li flow divertor, the magnetohydrodynamics (MHD)/heat transfer effects of the Li flowing inside the substrate as a second coolant were comprehensively investigated under the normal steady-state operation conditions. Here, in the present study, the multiphysics analysis is extended to the unsteady abnormal divertor scenario where the Li layer on top of the substrate does not provide full coverage or even totally disappears for a certain period of time, so that the substrate becomes directly exposed to the incident high plasma heat flux. Such an unwanted event may happen regardless of the concept of the divertor and is worth detailed investigations, typically referred to as a loss of coolant accident (LOCA). To address this situation, a simplified scoping analysis is conducted first in 2-D, and then an integrated 3-D modeling is performed using a time-dependent multiphysics model in COMSOL Multiphysics that integrates LM MHD, heat transfer, and solid mechanics. The main goal is to evaluate conditions under which the major material limits, such as the maximum allowable temperature, stress, and displacement of the substrate, can still be met. It was shown that the maximum time over which the substrate of RAFM steel can retain structural integrity during the LOCA is around 0.2 ∼ 0.3 s. Any divertor concept that utilizes RAFM steel as a substrate material and liquid Li as a second coolant should take such a permitted time into consideration.

divertor↗
DOE OSTI2025

Development and transferability of neural-network models for plasma-surface interactions

Plasma-surface interactions are increasingly critical to modern technologies; yet, accurate molecular dynamics simulations remain limited by the capabilities of interatomic potentials. Deep Potentials (DPs) promise to revolutionize the field by providing a systematic method for producing accurate interatomic potentials. The primary challenge of DP development is selecting a dataset, which efficiently spans the set of atomic environments one expects to encounter in the subsequent molecular dynamics simulations. The computational cost of density functional theory calculations, which are the typical basis for DP development, makes it impossible to directly verify the quality of a given DP. To address this challenge, we explore the development of a deep-learned interatomic potential, “DeepREBO,” trained to reproduce the behavior of the REBO2 empirical potential, enabling direct validation of training methodology and transferability. Using an active learning framework, we begin with a minimal dataset and iteratively expand it to train a Deep Potential-Smooth Edition model that faithfully reproduces REBO2 results for 25 eV hydrogen bombardment of diamond (001), a particularly challenging case. We show that small, carefully curated datasets can outperform large, unguided ones, with effective models requiring fewer than 15 000 snapshots. Subsequent transferability tests demonstrate that while DeepREBO generalizes well to diamond (111) surfaces, performance degrades for amorphous carbon or higher-energy impacts, highlighting the need for use-case-specific training data. We also evaluate methods to improve short-range repulsion. This study outlines best practices for training robust deep potentials and underscores the importance of dataset design for predictive plasma simulations.

Ab-initio molecular dynamics↗
DOE OSTI2025

What to Support When You’re Compressing

Over the last nearly 20 years, lossy compression has become an essential aspect of HPC applications’ data pipelines, allowing them to overcome limitations in storage capacity and bandwidth and, in some cases, increase computational throughput and capacity. However, with the adoption of lossy compression comes the requirement to assess and control the impact lossy compression has on scientific outcomes. In this work, we take a major step forward in describing the state of practice and by characterizing workloads. We examine applications’ needs and compressors’ capabilities across 9 different supercomputing application domains. We present 24 takeaways that provide best practices for applications, operational impacts for facilities achieving compressed data, and gaps in application needs not addressed by production compressors that point towards opportunities for future compression research.

Error-Bounded Lossy Compression↗
DOE OSTI2025

Integrating and Characterizing HPC Task Runtime Systems for hybrid AI-HPC workloads

Scientific workflows increasingly involve both HPC and machine-learning tasks, combining MPI-based simulations, training, and inference in a single execution. Launchers such as Slurm’s srun constrain concurrency and throughput, making them unsuitable for dynamic and heterogeneous workloads. We present a performance study of RADICAL-Pilot (RP) integrated with Flux and Dragon, two complementary runtime systems that enable hierarchical resource management and high-throughput function execution. Using synthetic and production-scale workloads on Frontier, we characterize the task execution properties of RP across runtime configurations. RP+Flux sustains up to 930 tasks/s, and RP+Flux+Dragon exceeds 1,500 tasks/s with over 99.6% utilization. In contrast, srun peaks at 152 tasks/s and degrades with scale, with utilization below 50%. For IMPECCABLE.v2 drug discovery campaign, RP+Flux reduces makespan by 30–60% relative to srun/Slurm and increases throughput more than four times on up to 1,024. These results demonstrate hybrid runtime integration in RP as a scalable approach for hybrid AI-HPC workloads.

HPC-AI↗
DOE OSTI2023

XGC-6D

XGC-6D is a geometric full-f Particle-In-Cell code on unstructured meshes including XGC-like unstructured mesh, tetrahedral mesh, etc. The code is particularly suitable to study boundary plasma with steep gradient and frequency higher than ion gyro-frequency. References: Wang Z, Qin H, Sturdevant B, Chang CS. Geometric electrostatic particle-in-cell algorithm on unstructured meshes. Journal of Plasma Physics. 2021;87(4):905870406.

Wang, Zhenyu [Princeton Plasma Physics Laboratory ↗
DOE OSTI2025

Exascale workflow applications and middleware: An ExaWorks retrospective

Exascale computers offer transformative capabilities to combine data-driven and learning-based approaches with traditional simulation applications to accelerate scientific discovery and insight. However, these software combinations and integrations are difficult to achieve due to the challenges of coordinating and deploying heterogeneous software components on diverse and massive platforms. Here, we present the ExaWorks project, which addresses many of these challenges. We developed a workflow Software Development Toolkit (SDK), a curated collection of workflow technologies that can be composed and interoperated through a common interface, engineered following current best practices, and specifically designed to work on HPC platforms. ExaWorks also developed PSI/J, a job management abstraction API, to simplify the construction of portable software components and applications that can be used over various HPC schedulers. The PSI/J API is a minimal interface for submitting and monitoring jobs and their execution state across multiple and commonly used HPC schedulers. We also describe several leading and innovative workflow examples of ExaWorks tools used on DOE leadership platforms. Furthermore, we discuss how our project is working with the workflow community, large computing facilities, and HPC platform vendors to address the requirements of workflows sustainably at the exascale.

97 MATHEMATICS AND COMPUTING↗