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

Sensitivity of the magnetic axis to variations in the magnetic field

Lagrangian variational principles are used to locate the magnetic axis and to determine its sensitivity to variations in the magnetic field. Only the change in the magnetic field on axis is required to determine, to lowest order, the location of the new magnetic axis. Using a simple example, namely, the magnetic field produced by a set of current-carrying filaments (a.k.a. coils), the variation in the coil geometry that induces the largest change in the magnetic axis is determined using singular value decomposition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Enhancement to the conductivity of surface transfer-doped (111) diamond through thermochemical surface etching

The use of a transition metal catalyzed thermochemical etching method for improving the carrier transport properties of the near-surface two-dimensional (2D) hole gas in surface transfer-doped hydrogen-terminated (111) diamond is demonstrated. Using Ni 0.8 Cr 0.2 films deposited and annealed to a temperature of 900 °C, with up to three etch cycles, preferential (111) surface etching produces large terraces exceeding 10 μm in size with a surface microroughness, σ 2 RMS λ, that is two orders of magnitude lower than for the pre-etched (111) surface. Magnetotransport measurements on hydrogen-terminated Hall bars engineered on the pre- and post-etched surfaces and rendered conductive by the adsorbed water layer formed on exposure to ambient conditions demonstrate that this etching causes an improvement in the hole mobility by an order of magnitude, resulting in a measured sheet resistivity of 1.04 kΩ/sq at a temperature of 4.2 K without gating.

42 ENGINEERING↗
DOE OSTI2024

Plasma sources for advanced semiconductor applications

Semiconductors are the foundation of modern technology, used in our personal, industrial, and military-grade devices. Further, every aspect of U.S. society is closely tied to semiconductors, and our economy cannot progress at the current pace with existing chip manufacturing methods as chip features approach an atomistic scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An experimental study of the existence regions and non-linear interactions of drift wave and Kelvin–Helmholtz instabilities in a linear magnetized plasma

Experimental observations of the intrinsic excitation and non-linear interactions of drift wave (DW) and Kelvin–Helmholtz (KH) instabilities in a linear magnetized plasma column are presented. The experiments are carried out in the inverse mirror plasma experimental device (IMPED)—a cylindrical, magnetized, linear plasma machine designed to study low-frequency waves and instabilities in plasma. A novel feature of IMPED is the ability to control plasma profiles, such as the density n(r)⁠, electron temperature T e (r)⁠, and plasma potential V p (r) by varying the ratio Rm of the magnetic field in the main chamber to that in the source chamber. At high values of Rm, higher-density gradient scale length promotes the drift wave (DW) instability while lower Rm value results in a higher radial electric field, inducing a sheared poloidal flow that enhances the dominance of the Kelvin–Helmholtz (KH) mode. The background and fluctuating plasma parameters are characterized using various configurations of multiple in situ electric probes at different spatial locations to quantify the local gradients that excite the low-frequency primary instabilities. Statistical, spectral, and bispectral analysis of the density and potential signals help identify these modes in terms of wave number, frequency, phase, and amplitude and also delineate the nature of their non-linear interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A generalized effective potential for differentially rotating plasmas

Global stability of differentially rotating plasma is investigated using a generalized effective potential. We first, for a current-free system, obtain a general form of an effective potential in terms of the free energies of global curvature and gradients of rotation for non-axisymmetric disturbances. We then examine the stability of differentially rotating disks for several rotation profiles and present the associated effective potential for the onset of these instabilities in the MHD regime. In particular, results for global axisymmetric magnetorotational instability as well as local and global non-axisymmetric modes are presented. The latter constitute two distinct non-axisymmetric modes, a high frequency local MRI and a global low-frequency non-axisymmetric mode (the magneto-curvature mode, introduced in Ebrahimi and Pharr [Astrophys. J. 936, 145 (2022)]), confined either between two Alfvénic resonances or an Alfvénic resonance and a boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Near-ideal relaxed MHD in slab geometry

We investigate the solutions of the relaxed magnetohydrodynamic (MHD) model (RxMHD) of R. Dewar and Z. Qu. This model generalizes Taylor relaxation by including the ideal Ohm's law constraint using an augmented Lagrangian method, providing a pathway to extend the multi-region relaxed MHD (MRxMHD) model. We present the first numerical solution of the RxMHD model by Dewar and Qu, demonstrating that it is mathematically well-defined and computationally feasible for constructing MHD equilibria in slab geometry. We also show that a cross-field flow can exist without enforcing an arbitrary constraint on the angular momentum, as is done in the case of MRxMHD with flow. Our results also demonstrate the self-organization of fully relaxed regions during the optimization, which was an important motivation behind developing this model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Verification of the Shafranov shift in free-boundary VMEC, DESC, and SPEC calculations with elliptical geometry

Here, given a set of current-carrying filaments that create a magnetic field with rotating elliptical flux surfaces, we present verification calculations using the VMEC, DESC, and SPEC free-boundary magnetohydrostatic equilibrium codes for both the “vacuum” equilibrium, for which there is no plasma pressure and no plasma currents, and for a non-zero pressure, zero toroidal current equilibrium. For the vacuum case, these codes are quantitatively compared to the magnetic field produced by the coils. For the non-zero pressure case, they are compared to each other. As the stepped pressure profile used in SPEC approaches the continuous profile used in VMEC and DESC, the SPEC magnetic axis approaches the VMEC and DESC axes. To within an error given by the difference between the pressure profiles, the VMEC, DESC, and SPEC calculations give the same pressure-induced geometric shift of the magnetic axis, known as the Shafranov shift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Effects of oxidation and impurities in lithium surfaces on the emitting wall plasma sheath

The use of lithium as a surface coating in fusion devices improves the plasma performance, but the change in the wall properties affects the secondary electron emission properties of the material. Lithium oxidizes easily, which drives the emission yield well above unity. We present here simulations demonstrating the change in the sheath structure from monotonic to the nonmonotonic space-charge limited sheath using an energy-dependent data-driven emission model, which self-consistently captures both secondary emission and backscattering populations. Increased secondary electron emission from the material has ramifications for the degradation and erosion of the wall. The results show that the oxidation leads to an increased electron energy flux into the wall and a reduced ion energy flux. The net transfer of energy to the surface is significantly greater for the oxidized case than for the pure lithium case. High backscattering rates of low-energy particles lead to a high re-emission rate at the wall.

Fusion reactors↗
DOE OSTI2025

Experimental demonstration of transient CHI start-up using a floating biased electrode configuration

Transient coaxial helicity injection (CHI) current start-up in a spherical tokamak (ST) has been demonstrated for the first time using a reactor-relevant floating single biased (FSB) electrode configuration on QUEST. Unlike previous applications of CHI on HIT-II and on NSTX in which the entire ST vessel was split in two using two toroidal ceramic rings as vacuum breaks, in the FSB electrode configuration on QUEST, only the lower outer divertor plate is electrically insulated from the rest of the vessel. This is the CHI injector region. A CHI discharge is generated by driving current, using an external power supply, along magnetic field lines that connect the inner and outer divertor plates in the CHI injector region. The CHI plasma evolution and closed flux formation on QUEST are very similar to those on NSTX and have generated toroidal currents of 150 kA with 100 kA of current flowing on closed flux surfaces. An alternate, force-free model approach has been used to derive the CHI scaling relations. The resulting CHI scaling relations from the model are consistent with the conventional CHI scaling relations, which state that the CHI generated toroidal current is proportional to the CHI injector flux and the CHI injector current is inversely proportional to the toroidal field. These results from QUEST bode well for the application of CHI in a new generation of compact high field STs and low aspect ratio tokamaks in which the space for the central solenoid is very restricted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Global gyrokinetic simulations of kinetic-ballooning-mode turbulence in Wendelstein 7-X

We apply the global gyrokinetic code GENE-3D to perform kinetic-ballooning-mode (KBM) turbulence simulations in the Wendelstein 7-X (W7-X) stellarator using experimentally relevant parameters. In contrast to the flux-tube approach, we observe saturation of the turbulent fluctuations at finite amplitude. Analyzing the turbulent characteristics, we identify the presence of weakly driven KBM turbulence hybridizing with ion-temperature-gradient and trapped-electron-mode turbulence in the scenario under consideration, showing that it is possible that KBM activity can be driven in W7-X, opening up the path for studying this type of turbulence, which is highly relevant for future fusion power plants.

Electrostatics↗
DOE OSTI2025

Poor confinement in stellarators at high energy II

This work discusses particle resonances in stellarators. There are two important effects. First, resonances provide locations for Alfvén modes, which can cause significant particle loss. Second, if a resonance matches the period of the modulation of the equilibrium field magnitude it produces islands in particle orbits increasing in size with particle energy, also possibly causing loss. A previous publication suffered from a restriction to the consideration of deeply passing particle orbits, clearly showing the resonance islands. Large particle magnetic moment increases the effect of particle cross field drift, essential for this phenomenon, so passing particles with larger values of magnetic moment give larger islands. It is possible to avoid designs exhibiting such resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

General kinetic ion-induced electron emission model for metallic walls applied to biased Z-pinch electrodes

A kinetic ion-induced electron emission (IIEE) model for general applications is developed to obtain the emitted electron energy spectrum for a distribution of ion impacts on a metallic surface. We assume an ionization cascade mechanism and use empirical models for the ion and electron stopping powers. The emission spectrum and the secondary electron yield (SEY) are validated for a variety of materials. The IIEE model is used to study the effect of IIEE on the plasma-material interactions of Z-pinch electrodes. Un-magnetized Boltzmann-Poisson simulations are performed for a Z-pinch plasma doubly bounded by two biased copper electrodes with and without IIEE at bias potentials from 0 to 9 kV. At the anode, the SEY decreases from 0 to 1 kV, but then increases at higher bias potentials. At the cathode, the SEY is much larger due to higher energy ion bombardment and grows with bias potential. As the bias potential increases, the emitted cathode electrons are accelerated to higher energies into the domain, collisionally heating the plasma. Above 1 kV, the heating is strong enough to increase the plasma potential. Despite SEY greater than 1, only a classical sheath forms as opposed to a space-charge limited or inverse sheath due to the emitted electron flux not reaching the space charge current saturation limits. Furthermore, the current in the emissionless cases saturates to a value lower than experiment. With IIEE, the current does not saturate and continues to increase with the 4 kV case, matching most closely with the experiment.

Carbon based materials↗
DOE OSTI2025

Validation of a comprehensive first-principles-based framework for predicting the performance of future stellarators

This paper presents the validation of the GENE-KNOSOS-Tango framework for recovering both the steady-state plasma profiles in the considered radial domain and selected turbulence trends in a stellarator. This framework couples the gyrokinetic turbulence code GENE, the neoclassical transport code KNOSOS, and the transport solver Tango in a multi-timescale simulation feedback loop. Ion-scale kinetic-electron and electron-scale adiabatic-ion flux-tube simulations were performed to evolve the density and temperature profiles for four OP1.2b W7-X scenarios. The simulated density and temperature profiles showed good agreement with the experimental data using a reasonable set of boundary conditions. Equally important was the reproduction of observed trends for several turbulence properties, such as density fluctuations and turbulent heat diffusivities. Key effects were also touched upon, such as electron-scale turbulence and the neoclassical radial electric field shear. The validation of the GENE-KNOSOS-Tango framework enables credible predictions of physical phenomena in stellarators and reactor performance based on a given set of edge parameters.

Gyrokinetic simulations↗
DOE OSTI2025

Velocity-space signatures of shock-drift acceleration at quasi-perpendicular collisionless shocks

The shock-drift acceleration of ions at quasi-perpendicular shocks is a well-known kinetic mechanism for the acceleration of a small fraction of incoming ions to high energy. Here, we use a suite of sixteen hybrid simulations of quasi-perpendicular collisionless shocks over the range of Alfvén Mach number 4.3 ≤ M A ≤ 15.8 (corresponding to a range of fast magnetosonic Mach numbers 2.6 ≤ M f ≤ 9.4) and shock-normal angle 45°≤ θ Bn ≤ 90° to identify the velocity-space signature of shock-drift acceleration using the field-particle correlation technique. We show that the features of the ion velocity distribution in the shock foot and ramp regions can be clearly interpreted by analysis of the single-particle trajectory of a reflected ion through the full 3D-3V phase space. The characteristic features of the velocity-space signature of shock-drift acceleration remain qualitatively robust over the full parameter range of our simulations, providing a potential means for its identification using single-point spacecraft measurements. At higher Alfvén Mach numbers M A ≳8 (M f ≳5), kinetic instabilities generate fluctuations of the electromagnetic fields within the shock transition region, leading to fluctuations in and smearing out of the resulting velocity-space signatures, but the signature remains generally robust and identifiable. The results on the shock-drift acceleration of ions presented here represent a novel means to determine more completely the partitioning of upstream bulk flow kinetic energy into plasma heating, particle acceleration, and electromagnetic fields in collisionless shocks.

Howes, Gregory G. [Univ. of Iowa, Iowa City, IA (U↗
DOE OSTI2025

Guiding center equations for the magnetic dipole

Since the discovery of Van Allen radiation belts in the 1960s, observations of energetic ions trapped in the Earth's dipole magnetic field have illustrated the remarkable confinement properties of this configuration. As such, it has been used for confining a hot plasma for nuclear fusion studies, starting from the pioneering work of Bo Lehnert and Akira Hasegawa, in the Levitated Dipole Experiment (LDX) at MIT until 2011 and in the RT-1 experiment at the University of Tokyo. More recently, the dipole has been subject to a renewed interest for fusion studies by a couple of startups and for smaller applications as a cold plasma source. While the equilibrium and magneto-hydrodynamic stability of the dipole have been investigated quite in detail, neoclassical properties of the dipole are comparatively much less known: the dipole is more known in geophysics than in fusion science. For this reason, in this paper, we propose a set of Hamiltonian, guiding-center equations to describe the motion of electrons and ions in a magnetic dipole configuration. We also developed a code, and we show the main features of particle motion, benchmarking our results with the analytical solutions for the bounce and precession motion, which are well documented in the literature. We also draw some general conclusions for the neoclassical transport in usual toroidal confinement schemes, such as the tokamak and the stellarator, pointing out the unique advantages of the dipole in confining energetic particles.

Hamiltonian mechanics↗
DOE OSTI2025

Embedded random phase approximation for magnetic systems: H 2 dissociative adsorption on Fe(110)

The random phase approximation (RPA), a method for treating electron correlation, has been shown to be superior to standard density functional theory (DFT) approximations in numerous cases. However, the RPA’s computational cost is substantially higher than that of DFT, particularly restricting its application to extended surfaces. The recently introduced embedded RPA (emb-RPA) approach [Wei et al., J. Chem. Phys. 159(19), 194108 (2023)] reduces this computational cost by approximately two orders of magnitude. While previous applications of emb-RPA focused on non-spin-polarized systems, here we extend the approach to ferromagnetic ones. Unlike other embedded correlated wavefunction methods, such as embedded complete active space self-consistent field theory, emb-RPA is advantageous for spin-polarized systems because the RPA is compatible with unrestricted DFT solutions, which are eigenfunctions of the spin angular momentum operator S z but not the total spin-squared operator S 2 . By applying emb-RPA with specific magnetization constraints, we achieved a speedup of two to three orders of magnitude (one order when accounting for the one-time embedding potential optimization cost) with only small errors (∼50 meV) compared to full periodic RPA. Moreover, emb-RPA significantly reduces the over-binding errors of DFT approximations. In conclusion, we anticipate that the acceleration enabled by the spin-polarized emb-RPA approach will broaden the applicability of RPA to magnetic materials.

Density functional theory↗
DOE OSTI2025

Fast ion stabilization of tilt in large radius FRCs

The field reversed configuration (FRC) has been a curious case in plasma physics research in that early MHD analysis suggested FRCs should be grossly unstable, while experimental results contradicted that prediction. Later, this theory was able to resolve this contradiction by understanding that finite Larmor radius effects largely negated the MHD predictions. Similarly, previous theoretical studies of beam driven FRCs predicted that such system would be unstable to beam driven modes while, again, experimental results indicated the contradiction. In this paper, we reconcile the theoretical understanding of beam driven modes with experimental observations of stability in these systems. By self-consistently capturing fast ion generation from neutral beam injection and its impact on the plasma equilibrium, we show that low amplitude perturbations in the magnetic field, driven by betatron particles, modify the precession frequencies of the betatron particles such that the drive for compressional Alfvén waves in the thermal plasma is reduced. Finally, we are able to demonstrate, for the first time, stable beam driven FRC evolution at high S*/E in 3D kinetic simulations.

Betatrons↗