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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Suppressing Screening-Current-Induced Strain in a 72-mm-Bore REBCO Insert for a 20-T Magnet: A Numerical Study

Rare-earth barium copper oxide (REBCO) high-temperature superconducting (HTS) magnets are considered a game changer for the capability of generating magnetic fields exceeding 20 T at or above liquid helium temperature, combined with the potential for substantially reduced manufacturing and operational costs. Princeton Plasma Physics Laboratory is dedicated to developing large-bore high-field superconducting magnet systems to support forefront physics research. Here, our current project focuses on a φ72-mm cold-bore REBCO insert comprising 21 dry-wound double-pancake coils, designed to generate at least 8 T when nested within a 12-T outsert magnet. A significant challenge for high-field REBCO magnets is the time-varying-magnetic-field-induced screening currents (SC) in the REBCO conductors, which can cause localized strain and stress concentrations. This paper presents a numerical study of the SC-induced strain in the REBCO insert magnet, confirming that the SC-induced strain can be substantially suppressed by energizing the REBCO insert before the outsert magnets. We discuss and reveal the mechanism behind this reduction. By applying this strategy, we expect to unleash the potential of the REBCO insert magnet to generate up to 12 T in a 12-T background field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Purpose, Design, and Analysis of the NSTX-U Centerstack Casing Lateral Support Shims

The National Spherical Torus Experiment Upgrade (NSTX-U) centerstack casing surrounds the central magnet core of NSTXU. It provides the inner vacuum boundary of the plasma chamber and supports the radially inboard plasma facing components lining the plasma chamber. The casing is supported like a cantilever, fixed at the lower flange by a connection to the TF Bundle through a G-10 ring. It is flexibly supported at the top by bellows, and insulated shims intended to react the largest part of the inventory of the disruption loads at the upper support, thereby protecting the upper and lower connections from excessive loads. A large part of the qualification of the casing and tiles requires addressing and quantifying disruption loads. These have been estimated based on extensions of the behavior of NSTX during its run period roughly 15 years ago. Structural qualification of the Upgrade used these estimates of the loading as input to calculations to demonstrate adequacy for a “most probable” disruption. The estimates are included in a requirements document to guide conservative design of all components loaded by disruptions. However, the changes made for the upgrade made the extrapolation of disruption behavior uncertain. Disruptions include axisymmetric loading from axisymmetric currents induced in the vessel shells. In addition, more complex currents called Halo currents are a result of plasma poloidal currents impinging on the casing. The halo currents add another layer of uncertain behavior and net lateral loading to the casing. As a result, NSTX-U is equipped with a number of instruments intended to assess the nature and severity of disruptions. The casing shims are required for lateral support of the casing, but these have also been converted to load cells to provide another measure of the magnitude of symmetric and asymmetric loading. Some of the logic in the development of the requirements document are presented. Electromagnetic analysis of axisymmetric and halo loading is described. Global simulations of the casing transient dynamic response, and local analysis and testing of the shim/load cell are included. The shim load cells and other instrumentation will be an important guide in planning NSTX-U experiments beyond its next run period expected in late 2026 or early 2027.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Impurity Powder Injection Experiments in the Large Helical Device

The Impurity Powder Dropper (IPD) is a device capable of injecting controlled amounts of sub-millimetre powder into the plasma under the action of gravity. In 2019 the IPD was first installed on the Large Helical Device (LHD) in Japan, with the aim of improving the plasma performances through real time boronization and assessing the compatibility of this technique with steady state operation. Extensive series of experiments have been performed using the IPD, focused on the improvement of the plasma performance via low-Z powder injection and the understanding of the underlying physical phenomena. In this article, we review the experiments that took place in the period 2019-2024. The main results include the demonstration of the improvement of the wall conditions (reduction of intrinsic impurity content, wall recycling) both on a shot-to-shot basis and in real time. Furthermore, a reduced-turbulence improved confinement regime has been observed coincident with powder injection, resulting in an increase of the plasma temperature of the order of 25%, with enhancements that can reach up to 50% for ion temperature.

Nespoli, F. [Princeton Plasma Physics Laboratory (↗

DOE OSTI2024

Net lithium deposition and dominant self-sputtering in lithium tokamak experiment-β with a liquid lithium wall

We observed enhanced net lithium deposition and lithium erosion, possibly dominated by physical sputtering of lithium by lithium-ion bombardment, on the outer plasma-facing surface in the Lithium Tokamak eXperiment-β (LTX-β) during liquid lithium wall operations. Silicon crystal samples with micro-trenches (30 μm × 30 μm × 2–7 μm deep) were exposed to hydrogen plasmas in LTX-β for solid and liquid lithium wall operations. Post-mortem analysis using X-ray photoelectron spectroscopy combined with argon ion sputtering measured net lithium deposition of 8.2 or 21 nm on the silicon crystal surface exposed for ~ 50 repeated shots of ~ 50-ms hydrogen plasma discharges during the liquid lithium wall operations at a vessel temperature of 475 K. Energy dispersive X-ray spectroscopy measured oxygen concentration patterns on the micro-trench floors, which were due to oxidized lithium deposition. Using the inhomogeneous oxygen concentration pattern caused by an ion-shadowing effect associated with the micro-trench’s geometric structure, we determined a polar incident ion direction of 68.4 ± 1.6° referenced to the surface normal direction. This observation was well-explained by the hypothesis that self-sputtering of Li was a dominant lithium erosion source in addition to lithium sputtering by hydrogen bombardment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Taming harmful bursts and heat flux in high-confinement tokamak plasmas

A major challenge in tokamak fusion research is first-wall erosion caused by steady heat loads and sudden energy bursts known as edge-localized modes. Divertor detachment reduces steady-state heat flux, while resonant magnetic perturbations can suppress these instabilities. However, integrating the two has been difficult because they require conflicting operating conditions. Here we demonstrate simultaneous achievement of resonant magnetic perturbations mitigated small edge-localized modes and impurity seeded partial divertor detachment in plasmas with an ITER-similar shape on the DIII-D tokamak. Experiments and simulations show that resonant magnetic perturbations facilitate detachment by redistributing particles, lowering the core density and increasing the scrape-off layer density, thereby reducing the amount of injected gas required. Cooling-gas injection eliminates the secondary heat-flux peak created by three-dimensional magnetic lobes, while edge cooling weakens the plasma response to the applied magnetic fields. These advances illustrate a viable pathway for integrating edge stability control with power exhaust in future fusion reactors.

magnetically confined plasmas↗

DOE OSTI2026

Effects of density turbulence on helicon wave propagation in the core plasmas

Radio frequency wave propagation can be significantly affected by density irregularities, such as filaments in the scrape-off layer or instabilities in the core plasma. In this study, we examine the impact of edge turbulence on helicon wave propagation using the Petra-M simulation tool. To analyze the effect of edge turbulence, we utilize a realistic background plasma derived from XGC simulations, which includes spatial density fluctuations at the edge. This focus is particularly relevant for the DIII-D configuration, characterized by edge density turbulence of the core plasmas in a wide pedestal QH-mode. We focus on helicon wave propagation in the core plasma since the slow mode cannot propagate into the core plasma when the density is higher than the lower hybrid resonance, where a wide pedestal QH-mode occurs. The simulation results indicate that edge density fluctuations have a substantial impact on helicon wave coupling due to scattering in the core. We specifically demonstrate that the toroidal mode number and the level of density fluctuations are significant factors influencing scattering. While we have minimized the excitation of slow modes from the antenna, we still observe mode-converted slow modes resulting from the incoming helicon waves. The insights derived from these simulations will inform upcoming tokamak experiments regarding helicon antenna coupling in long pulse scenarios.

Kim, Eun-Hwa [Princeton Plasma Physics Laboratory ↗

DOE OSTI2026

A MINi Powder INjector (MINPIN) for the SMART tokamak

Here, this article presents the design and performance characterization of a MINi Powder INjector (MINPIN) developed for boron wall conditioning (boronization) in the SMall Aspect Ratio Tokamak (SMART) at the University of Seville, Spain. SMART studies positive and negative triangularity plasmas within the low aspect ratio range of 1.4 < R/a < 3, with a major radius (R) of 0.40–0.60 m and a minor radius (a) of 0.20–0.33 m. MINPIN, in design, will gravitationally deliver boron (B) powder into the plasma, which ablates and transports the injected material. Mounted on the top of the vessel, the MINPIN will drop controlled amounts of powder using a piezoelectric blade driven by a signal generator that vibrates at a set frequency. The B will then be deposited in a layer on SMART’s stainless steel vacuum vessel walls, and this boronization process reduces the influx of intrinsic wall impurities, such as oxygen, carbon, nitrogen, and iron, which radiate power and cool the plasma. In the laboratory, we calibrated the drop rate for 1–10 mg/s, which will enable controlled boronization coating while achieving efficient powder ablation. We propose MINPIN’s operational scenarios based on the calibrated drop rate and estimated drop duration, assuming the distance from MINPIN to SMART’s plasmas.

Boron powder↗

DOE OSTI2026

Assessing plasma face component thermal response to rotating 3D magnetic fields for SPARC tokamak

Thermal response simulations of plasma-facing components (PFCs) in the SPARC tokamak, performed with the HEAT code, show that three-dimensional (3D) heat loads resulting from stationary n=1 perturbations require highly radiative scenarios, with up to 95% of the power crossing the separatrix (P SOL ) being radiated, to maintain PFC temperatures within acceptable operational limits, whereas the application of slowly rotating 3D fields substantially reduces the thermal loads. The HEAT module, developed to predict heat loads from non-axisymmetric plasmas, is extended to model time-dependent heat flux patterns generated by rotating 3D fields, and a comprehensive thermal analysis is performed on PFCs subjected to both the maximum and minimum power loads, as well as to rotating heat flux distributions, to evaluate the temperature evolution for varying perturbation amplitudes and rotation frequencies. The extension of this analysis to 3D fields with toroidal mode number n=2 shows that this configuration leads to weaker localized heat flux peaks relative to the n=1 case, enabling safe operation with less than 80% of the power radiated when static 3D fields of low amplitude are applied, while using slowly rotating fields at higher amplitudes. These results indicate that n=2 perturbations are generally less detrimental to divertor power exhaust, emphasizing the strong dependence of divertor power exhaust on the characteristics of the applied 3D fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

MHD Analysis of Dual-Coolant Lead-Lithium Blanket for Spherical Tokamak Advanced Reactor

The tritium breeding blanket is vital for future fusion power plants, with the Spherical Tokamak Advanced Reactor (STAR) project highlighting the dual-coolant lead-lithium (DCLL) design. The DCLL blanket performs shielding, energy exhaust, and tritium breeding using a lead-lithium alloy , with lithium as the breeder and lead as the neutron multiplier. It also serves as the primary coolant, with helium providing supplemental cooling. Reduced-activation ferritic/martensitic steel is used for the blanket structure. Magnetohydrodynamic (MHD) phenomena influence the liquid metal flow in a magnetic field, affecting heat transfer in the breeder affected by energetic neutrons. Understanding key flow parameters in such conditions is critical for efficient DCLL design. This study uses three-dimensional thermofluid MHD analysis with ANSYS CFX software, modified at Princeton Plasma Physics Laboratory, to simulate high Hartmann flows. The neutronics code MCNP, coupled with plasma equilibrium, provides heat source distribution. In conclusion, we examine electromagnetic interactions in adjacent fluid domains and analyze the magnetic field’s impact on flow distribution in the inboard and outboard blanket layout, using detailed mesh generation for accurate results.

DCLL↗

DOE OSTI2025

Investigation of radiated-power for low aspect ratio fusion plasmas

Investigating radiated power in fusion plasmas is of utmost importance to understand the effect of undesired impurities, such as metals, in present devices, or also desired impurities, such as noble gasses, to purposefully radiate a large fraction of power in future devices. These studies are especially important for high Z impurities which will play a crucial role in future generation fusion pilot plants. In this work, we have developed a power radiation analysis module, which is used to investigate 2D distributions of impurity densities and radiated power asymmetries caused due to both plasma rotation and the cooling rate dependence on temperature profile for the cases of experimental NSTX plasmas and designed scenarios for the spherical tokamak advanced reactor (STAR), both being low aspect ratio tokamaks. Two different atomic databases have been tested during this work to study their impact on the radiated power distribution, especially for high Z impurities. Also in PRAM, self-consistent calculations of two-dimensional electron, main ion, and impurity ion densities are carried out using one-dimensional input density, temperature, and rotation profiles. In the case of NSTX, discharges with high rotation of ∼ 170 km s, measured with charge exchange recombination spectroscopy, have been investigated. Rotation-induced charge separation, leading to an electrostatic potential, is calculated iteratively to a self-consistent solution while testing high Z impurities to observe any 2D asymmetry in the core radiated power due to centrifugal forces. The STAR design, being much larger (R = 4 m), is projected to have a much lower rotation, and is shown to have low rotation-induced asymmetries, on the order of ten percent or less, between the low field and high field sides. However, another effect not due to rotation but to the dependence of impurity cooling rates on temperature can lead to radiation peaking off-axis, near the plasma edge. This effect is noticeable for argon in NSTX, for example, but can also be enhanced for certain impurities at much higher temperatures projected for STAR (T e0 ~ 32 keV), for example for undesired tungsten or possibly desired xenon.

NSTX↗

DOE OSTI2026

Real-time reconstruction and control of pedestal-top electron density using RMP and gas puff at KSTAR

We report the experimental results of controlling the pedestal-top electron density by applying resonant magnetic perturbation (RMP) with in-vessel control coils and the main gas puff in the 2024-2025 KSTAR experimental campaign. The density is reconstructed using a parameterized $ψ_N$ grid and five channels of line-averaged density measured by the two-colored interferometer (TCI). The reconstruction procedure is accelerated by deploying a multi-layer perceptron to run in approximately 120 µ s and is sufficiently fast for real-time control. A proportional-integral controller was adopted, with the controller gains estimated from the system identification procedure. The experimental results demonstrate that the developed controller can follow a dynamic target while exclusively using both actuators. The absolute percentage errors between the electron density at $ψ_N$ = 0.89 and the target were approximately 1.5% median and a 2.5% average, respectively. The developed controller can even lower the density by using the pump-out mechanism under RMP, and it can follow a more dynamic range of density targets than a single actuator controller. The developed controller will enable experimental scenario exploration within a shot by dynamically setting the density target or maintaining a constant electron density within a discharge.

EFIT↗

DOE OSTI2026

Continuum radiated power density $(P^{\textrm{C}}_{\textrm{rad}})$ and effective charge $(Z_{\textrm{eff}})$ estimates from multi-energy photon-counting measurements

Multi-energy soft x-ray pinhole cameras have been designed, built, calibrated, and operated at Madison symmetric Torus, Alcator C-Mod, and more recently at Tokamak a configuration variable and Tungsten Environment in Steady-State Tokamak (WEST), to measure plasma emission across multiple energy ranges. Here we describe a new methodology to estimate the local continuum radiated power density and the plasma effective charge (Z eff ) directly from photon-counting measurements of the line-free continuum emission (Bremsstrahlung and Radiative Recombination) in several energy bands between 11 and 18 keV. This capability is particularly valuable for confinement systems using metal plasma-facing components, where x-ray losses from interactions with the sputtered wall can represent a significant fraction of the total radiated power (P rad ). The approach leverages a well-characterized detector responsivity, modeled by a complementary error function, and interprets the differential multi-energy measurements between adjacent energy levels through the probability density function of a Gaussian distribution. The implementation of this diagnostic technique is currently under development on the WEST tokamak, aiming at the goal of providing real-time P rad and Z eff measurements during long-pulse operation (up to 1000 s) in the 2026 campaign.

WEST↗

DOE OSTI2025

Development and validation of non-axisymmetric heat flux simulations with 3D fields using the HEAT code

A new comprehensive module to simulate heat fluxes from three-dimensional (3D) magnetic fields has been implemented in the HEAT code. Especially compact tokamaks like SPARC require tools to predict and manage large heat fluxes. Existing release versions of HEAT can only simulate axisymmetric heat flux on 3D plasma facing components. The new module uses an M3D-C1 perturbed equilibrium and the MAFOT code to trace field lines of the perturbed 3D magnetic field. Heat flux is then assigned to the resulting footprints via a 3D layer model. The model distinguishes between the scrape-off layer, the magnetic lobes and the private flux region, and employs only 0D parameters like the layer width, diffusive spread and the last closed flux surface position in the perturbed edge to generate a heat flux profile. The magnitude is normalized to the total input power. Resulting heat flux simulations are compared and validated against infrared measurements in the DIII-D tokamak with applied 3D fields; good agreement is found for several cases. The new module can now be applied to the SPARC tokamak; a preliminary result for applied rotating 3D fields is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

3D modeling of n = 1 RMP driven heat fluxes on the SPARC tokamak PFCs using HEAT

3D heat flux calculations at the lower outer divertor plate of SPARC using the HEAT code show that 3D fields generated from error field correction coils can lead to enhanced peak heat fluxes up to 15 times larger compared to the axisymmetric case. Previously employed to simulate axisymmetric heat flux on 3D plasma facing components, the HEAT code can now predict 3D heat flux generated by non-axisymmetric plasmas. This is achieved via a new HEAT module which leverages the 3D field line tracing capabilities of MAFOT starting from an M3D-C1 (MHD resistive code) perturbed equilibria. The resulting heat flux is assigned using the magnetic footprint and the heat flux layer model, an extension of the 2D heat flux model also known as the Eich, to 3D non-axisymmetric plasmas. For SPARC, the new capabilities of HEAT are used to calculate the 3D heat loads resulting from n = 1 perturbation fields (with n indicating the toroidal periodicity) applied through a toroidal array of six picture frame coils with different amplitude. The comparison with the unperturbed case shows significant changes in shape and intensity of the heat flux profile. The results show that the application of n = 1 3D field leads to a localized enhancement of the heat flux peak, influenced by the wetted area impacted by the magnetic footprint, and the appearance of a secondary heat flux peak, whose intensity depends on amplitude of the applied 3D field and toroidal location.

3D heat flux↗

DOE OSTI2025

Active tungsten expulsion in ELM-absent H-mode plasmas via on-demand ELM triggering with lithium granule injection

Lithium granules gravitationally injected into the upper X-point region demonstrated on-demand edge-localized modes (ELM) triggering in otherwise ELM-suppressed H-mode plasmas on the Experimental Advanced Superconducting Tokamak. Sub-mm lithium granules dropped into enhanced D-alpha H-mode plasmas achieved a high triggering efficiency, while enabling ELM frequencies from several to hundreds of hertz. Core radiation from heavy impurities, dominated by W, was reduced by up to 60%, and the normalized energy confinement increased by up to 30%. At low injection frequencies, ELMs of substantially reduced size compared to spontaneous type-I ELMs were observed. At high injection frequencies, a transition to a mixed ELM phase occurred, characterized by intermittent larger ELMs and suppression of the quasi-coherent mode, achieving the most significant W reduction and energy confinement improvement. These results highlight a promising pathway for active W control via controlled, small ELMs in long-pulse, high-performance scenarios.

EDA H-mode↗

DOE OSTI2025

HIPED: Machine learning framework for spherical tokamak pedestal prediction and optimization

We introduce a Machine Learning framework, HIPED (HeIght and width Predictor for Edge Dynamics), for predicting and optimizing pedestal and core performance in spherical tokamak plasmas. Trained on pedestal and core datasets from the third MAST-U campaign, HIPED provides accurate estimates of pedestal height and width. The results reveal notable differences compared with conventional aspect-ratio studies; for instance, a simple power-law relation between pedestal width and height has very low accuracy. Instead, additional parameters such as normalized plasma pressure, elongation, and Greenwald fraction significantly improve accuracy. HIPED can also be trained only on `control room parameters' to inform experimentalists of which controllable parameters to adjust for improving core-integrated performance. The framework further includes a multi-objective optimization scheme that helps guide experimental planning and optimization. We find Pareto-optimal discharges with respect to various features, including distance from edge-localized modes and normalized plasma pressure, track their parameter trajectories over time, and identify the control room parameters required for these Pareto-optimal discharges. This provides a framework for systematically optimizing core and edge performance according to different experimental priorities.

Parisi, Jason F. [Princeton Plasma Physics Laborat↗

DOE OSTI2026

An overview of 3D field optimization for control of transport and edge instabilities on KSTAR

An international team from several laboratories and universities has made key advances over the last few years in the control of plasma transport and edge instabilities with applied 3D fields in the KSTAR tokamak to optimize long pulse operation scenarios. This overview begins with the optimization of both core and edge resonant magnetic perturbations (RMP) to improve fast ion confinement to avoid excessive limiter heat loads due to fast ion losses and successful modeling of the experimental results. Integrated and advanced plasma control techniques with machine learning (ML) and adaptive control were then used to optimize the 3D field spectrum in real-time to control edge localized modes (ELMs) while avoiding core locked modes that could disrupt the plasma. Accelerating the offline model of 3D fields with a surrogate ML model can optimize ELM suppression in the edge while limiting the impact of the applied RMP fields deeper in the plasma core in real-time. In addition, the impact of the 3D fields on the divertor heat load has been modeled and compared with experimental measurements. An analysis of a multi-machine database including KSTAR has been performed to better understand the metrics for the observed RMP thresholds for ELM suppression and the resulting plasma performance. Predictive modeling of the operational space for ELM suppression and density pumpout due to RMP has shown the importance of magnetic islands in the plasma edge and their impact on plasma turbulence. This research has culminated in the development of successful long pulse operational scenarios on KSTAR while attempting to overcome challenges of the new tungsten divertor.

3D fields↗

DOE OSTI2026

Simultaneous ELM suppression and divertor detachment via synergistic boron powder and neon injection in EAST

A novel approach for simultaneous power exhaust and edge-localized mode (ELM) control is presented in the Experimental Advanced Superconducting Tokamak discharges, which utilize an ITER-like tungsten divertor. Real-time injection of boron (B) powder and neon (Ne) gas overcomes their limitations encountered when used separately. Pure Ne seeding leads to a narrow operational window constrained by core impurity accumulation and H-mode to L-mode back transitions, while pure solid B injection (SBI) is insufficient for effective divertor cooling. In comparison, their combined use achieves a stable, stationary, ELM-suppressed H-mode with adequate power exhaust. This synergistic scenario features partial energy detachment at the outer divertor while maintaining good plasma confinement (H 98 ∼ 1) with minimal degradation. Two key features of this scenario are: (1) the SBI triggers a persistent Edge Harmonic Mode (EHM), which provides a crucial continuous particle transport channel, preventing Ne and tungsten/molybdenum accumulation without flushing out by ELM, and (2) the B + Ne mixture allows for active optimization of the radiated power profile. Core radiation can be reduced by substituting a portion of the Ne with B, leveraging their complementary non-coronal equilibrium radiation efficiencies. This combined B + Ne injection scheme presents a promising pathway toward integrated core-edge scenarios, offering the potential to minimize total impurity throughput while leveraging an actuator (powder injection) already being considered for ITER.

Boron↗