Search NASA⌕ Search

SEARCH · Search NASA

Results for “Large Plasma Device”

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

Open‐Air Combustion Synthesis with Rapid Plasma Processing of Large‐Area Transparent Conducting Oxides

A vacuum-free, high-throughput synthesis of indium tin oxide (ITO) via Combustion Oxidation with Rapid Plasma Processing (CORP) utilizes a solution-based exothermic combustion reaction to generate the oxide with tunable control of either amorphous or crystalline phases. A subsequent open-air, forming gas plasma treatment is used to introduce oxygen vacancies and promote crystallization. Here, the evolution of the oxide structure is elucidated by extended X-ray absorption spectroscopy fine structure analysis. Using CORP, fabrication of 300 cm 2 of ITO possessing a champion sheet resistance of 38 Ω sq. −1 , visible transmission of 89%, conductivity stability for over 250 days, roughness < 2nm, and Haacke figure of merit (%T 550nm 10 /R s ) of 0.012 Ω −1 is achieved. Cost modeling of CORP demonstrates up to a 67% reduction in price for TCOs using fully continuous, in-line unit operations compared with vacuum sputtering. The work shows a path toward a low-cost, vacuum-free manufacturing method for TCOs at commercial scales.

42 ENGINEERING↗
DOE OSTI2021

Fiber-optic silicon Fabry–Perot interferometric bolometer with improved detection limit for magnetic confinement fusion

Fiber-optic bolometers (FOBs) intended for plasma radiation measurement in magnetically confined fusion have been previously developed using a silicon pillar that functions as both a Fabry–Perot interferometer (FPI) for temperature measurement and an absorber for the radiation. We report an FOB design that can significantly improve the detection sensitivity over earlier designs by engineering the absorber of the FOB. Our design uses the fact that, compared with a silicon pillar, a gold film with the same x-ray absorption thickness will show a much higher temperature rise from a given power density of the radiation. Therefore, the responsivity of an FOB can be improved by attaching a large gold disk to the silicon FPI as the absorber. We have developed a fabrication method for FOBs of such design and obtained an FOB with a 4-µm-thick, 0.6-mm-diameter gold disk attached to a 200-µm-diameter, 100-µm-thick silicon FPI. We have characterized the noise, responsivity, response time, and noise-equivalent power density (NEPD) and compared these with the earlier design where the absorber is mainly the silicon FPI itself. The experimental result suggests that the FOB with the gold disk achieves a responsivity of ∼2.8 mK/(W/m2) and a noise-equivalent-power-density of 0.13 W/m2, which are, respectively, more than nine times larger and six times smaller compared to the FOB using a previous design. Improved NEPD and good absorption over a broad frequency range will make the FOB more attractive for applications in magnetic-confinement fusion devices.

Uddin, Nezam↗
DOE OSTI2024

The CXSFIT spectral fitting code: Past, present and future

Magnetically confined plasma experiments generate a wealth of spectroscopic data. The first step toward extracting physical parameters is to fit a spectral model to the often complex spectra. The CXSFIT (Charge eXchange Spectroscopy FITting) spectral fitting code was originally developed for fitting charge exchange spectra on JET from the late 1980s onward and has been further developed over decades to keep up with the needs of the users. The primary use is to efficiently fit a large number of spectra with many constrained Gaussian spectral lines of which the physical parameters can be coupled in a user-friendly manner. More recent additions to the code include time-dependent couplings between parameters, flexible background subtraction, and a non-linear coupling scheme between fit parameters. The latter was a pre-requisite for implementing Zeeman and motional Stark effect multiplets in the library of spectral features. The ability to save and replay “fit recipes,” even when multiple iterations are required, has ensured the traceability of the results and is one of the keys to the longevity and success of the code. The code is also in use on other tokamaks (AUG, ST-40) and to fit data from other spectroscopic diagnostics on JET. In this paper, we document the current capabilities and philosophy behind the structure of the code, including some of the algorithms used to calculate spectral features numerically efficiently. We also provide an outline of how CXSFIT could be transferred into a framework that would be able to meet the spectral fitting requirements of future devices, such as ITER.

Delabie, Ephrem G.↗
DOE OSTI2022

Hybrid simulation of Alfvén wave parametric decay instability in a laboratory relevant plasma

Large-amplitude Alfvén waves are subject to parametric decays, which can have important consequences in space, astrophysical, and fusion plasmas. Although the Alfvén wave parametric decay instability (PDI) was predicted decades ago, observational evidence is limited, stimulating considerable interest in laboratory demonstration of the instability and associated numerical modeling. Here, we report an important step toward direct hybrid simulation of the Alfvén wave PDI in a laboratory plasma, using antenna-like wave injection of a circularly polarized wave and realistic wave-plasma parameters. Considering collisionless damping, we identify the threshold Alfvén wave amplitudes and frequencies required for triggering the instability. These threshold behaviors are corroborated by simple theoretical analysis, considering the balance between PDI growth and ion Landau damping and PDI development in a bounded plasma. Other effects not included in the present model such as finite transverse wave scale and ion–neutral collisions are briefly discussed. These hybrid simulations promise to be a useful tool for investigating laboratory Alfvén wave dynamics and may provide guidance for future laboratory demonstration of the PDI process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Exploring the transition from continuous turbulence fluctuations to bursting ELMs in high SOL density regimes

BOUT++ turbulence simulations of the DIII-D reveal that the density profile between the separatrix and pedestal plays a crucial role in the dynamics of edge localized modes (ELMs) and edge plasma turbulent transport. Nonlinear simulations demonstrate that small ELMs in the DIII-D hybrid scenario under high SOL density conditions are predominantly driven by local ballooning modes near the separatrix, stabilizing global instabilities while enhancing localized pressure fluctuations. A key control parameters for ELM dynamics is the separatrix-to-pedestal density ratio, n e,sep /n e,ped . A high ratio indicates a shallow gradient, favoring small ELMs, while a lower ratio signals a steep gradient, which increases the likelihood of large ELMs. Comprehensive parameter scans, including n e,sep /n e,ped , density gradient profiles near the separatrix, and resistivity, reveal the critical role of these parameters in shaping transitions between turbulence-driven transport and ELM bursting. The scans demonstrate that in high SOL density regimes, small ELMs can result from either global resistive MHD instabilities or local ballooning modes near the separatrix, depending on the steepness of the separatrix density gradient. These findings also highlight the transition from continuous turbulence to small ELMs. The post-crash peak in pressure fluctuations, δP rms serves as a critical metric for identifying transition from continuous turbulence fluctuations to ELM bursting. Larger δP rms values correlate with ELM bursts driven by local or global instabilities, whereas smaller values indicate turbulence-dominated transport. Drift-Alfvén and resistive ballooning turbulence enhance the entrainment of fluctuations from the pedestal to the SOL, contributing to the complex interplay of dynamics in this regime. These findings emphasize the importance of separatrix density shaping and pedestal gradient control for optimizing ELM behavior in ITER and future fusion devices.

Li, Nami [Lawrence Livermore National Laboratory (↗
DOE OSTI2024

Relative sensitivity of plastic scintillator: A comparative analysis with 60 Co gamma rays, deuterium–deuterium, and deuterium–tritium neutrons

A plastic scintillator has found extensive application in the realm of high-energy physics and national security science. Many applications in those fields often involve the simultaneous production of photons, neutrons, and charged particles, which makes the relative sensitivity information for these different radiation types important. In this study, we have adopted a multi-head detector comprised of a plastic scintillator and high gain phototubes, which provides a large dynamic range and linearity. A comparative study on the relative sensitivities of plastic scintillators was facilitated by adopting three distinct radiation calibration sources (i.e., 60 Co γ rays, DD neutrons, and DT neutrons). Neutrons from a DD source generate a comparable level of scintillation to gamma rays emitted by 60 Co (i.e., 60 Co-γ/DD-n = 0.92 ± 16%). DT neutrons induce ~3.5 times the scintillation observed with DD neutrons (i.e., DT-n/DD-n = 3.5 ± 28%). In addition, the Geant4 simulation granted us valuable insights into the relative sensitivity of the scintillator. This comparative study will provide a useful database for users in diverse applications.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Enhanced Pedestal H-mode at low edge ion collisionality on NSTX

The Enhanced Pedestal (EP) H-mode regime is an attractive wide-pedestal ELM-free high-betap scenario for NSTX-U and next-step devices as it achieves enhanced energy confinement (H98y,2 > 1.5), large normalized pressure (betaN > 5) and significant bootstrap fraction (f_BS > 0.6) at I_p/B_T = 2 MA/T. This regime is realized when the edge ion collisionality becomes sufficiently small that a positive feedback interaction occurs between a reduction in the ion neoclassical energy transport and an increase in the particle transport from pressure-driven edge instabilities. EP H-mode was most often observed as a transition following a large ELM in conditions with low edge neutral recycling. It is hypothesized that the onset of pressure-driven instabilities prior to the full recovery of the neutral density leads to a temporary period with elevated ion temperature gradient that triggers the transition to EP H-mode. Linear CGYRO and M3D-C1 calculations are compared to beam emission spectroscopy (BES) and magnetic spectroscopy in order to describe the evolution of the edge particle transport mechanisms during the ELM recovery and the saturated EP H-mode state. The observations are consistent with the hypothesis that the onset of pressure-driven edge instabilities, such as the KBM and kink-peeling, can be responsible for the increased particle transport in EP H-mode.

edge and boundary physics↗
DOE OSTI2023

Microstability of β ~ 1 tokamak equilibria

High-power-density tokamaks offer a potential solution to design cost-effective fusion devices. One way to achieve high power density is to operate at a high β value (the ratio of thermal to magnetic pressure), i.e. β~1 . However, a β~1 state may be unstable to various pressure- and current-driven instabilities or have unfavourable microstability properties. To explore these possibilities, we generate β~1 equilibria and investigate their stability. First, we demonstrate the generation of high-β equilibria with the computer code VMEC. We then analyse these equilibria to determine their stability against the infinite-n ideal-ballooning mode. We follow that by engaging in a detailed microstability study using the GS2 code, beginning with assessments of electrostatic ion-temperature-gradient and trapped election mode instabilities. We observe interesting behaviour for the high-β equilibria – stabilization of these modes through two distinct mechanisms – large negative local shear and reversal of electron precession drift. Finally, we perform electromagnetic gyrokinetic simulations and observe enhanced stability in the outer core of high-β equilibria and absence of kinetic ballooning modes in the negative-triangularity, high-β equilibria. The enhanced outer-core stability of high-β equilibria is different from their lower-β counterparts and offers an alternative, potentially favourable regime of tokamak operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Exploration of the parameter space of quasisymmetric stellarator vacuum fields through adjoint optimisation

Optimising stellarators for quasisymmetry leads to strongly reduced collisional transport and energetic particle losses compared with unoptimised configurations. Although stellarators with precise quasisymmetry have been obtained in the past, it remains unclear how broad the parameter space is where good quasisymmetry may be achieved. We study the range of aspect ratios and rotational transform values for which stellarators with excellent quasisymmetry on the boundary can be obtained. A large number of Fourier harmonics is included in the boundary representation, which is made computationally tractable by the use of adjoint methods to enable fast gradient-based optimisation and by the direct optimisation of vacuum magnetic fields, which converge more robustly compared with solutions from magnetohydrostatics. Several novel configurations are presented, including stellarators with record levels of quasisymmetry on a surface, three field period quasiaxisymmetric stellarators with substantial magnetic shear, and compact quasisymmetric stellarators at low aspect ratios similar to tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Global gyrokinetic simulations of electrostatic microturbulent transport in LHD stellarator with boron impurity

Abstract Global gyrokinetic simulations of electrostatic microturbulent transport for discharge # 166256 of the Large Helical Device stellarator in the presence of boron impurity show the co-existence of the ion temperature gradient (ITG) turbulence and trapped electron mode (TEM) turbulence before and during boron powder injection. ITG turbulence dominates in the core, whereas TEM dominates near the edge, consistent with the experimental observations. Linear TEM frequency increases from ∼ 80 kHz to ∼ 100 kHz during boron injection, and ITG frequency decreases from ∼ 20 kHz to ∼ 13 kHz, consistent with the experiments. The poloidal wave number spectrum is broad for both ITG (0–0.5 mm −1 ) and TEM (0–0.25 mm −1 ). The nonlinear simulations with boron impurity show a reduction in the heat conductivity compared to the case without boron. The comparison of the nonlinear transport before and during boron injection shows that the ion heat transport is substantially reduced in the region where the TEM is dominant. However, the average electron heat transport throughout the radial domain and the average ion heat transport in the region where the ITG is dominant are similar. The simulations with boron show the effective heat conductivity values qualitatively agree with the estimate obtained from the experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Selective CW Laser Synthesis of MoS2 and Mixture of MoS2 and MoO2 from (NH4)2MoS4 Film

Very recently, the synthesis of 2D MoS2 and WS2 through pulsed laser-directed thermolysis can achieve wafer-scale and large-area structures, in ambient conditions. In this paper, we report the synthesis of MoS2 and MoS2 oxides from (NH4)2MoS4 film using a visible continuous-wave (CW) laser at 532 nm, instead of the infrared pulsed laser for the laser-directed thermolysis. The (NH4)2MoS4 film is prepared by dissolving its crystal powder in DI water, sonicating the solution, and dip-coating onto a glass slide. We observed a laser intensity threshold for the laser synthesis of MoS2, however, it occurred in a narrow laser intensity range. Above that range, a mixture of MoS2 and MoO2 is formed, which can be used for a memristor device, as demonstrated by other research groups. We did not observe a mixture of MoS2 and MoO3 in the laser thermolysis of (NH4)2MoS4. The laser synthesis of MoS2 in a line pattern is also achieved through laser scanning. Due to of the ease of CW beam steering and the fine control of laser intensities, this study can lead toward the CW laser-directed thermolysis of (NH4)2MoS4 film for the fast, non-vacuum, patternable, and wafer-scale synthesis of 2D MoS2.

(NH4)2MoS4↗
DOE OSTI2023

The DESC stellarator code suite Part 3: Quasi-symmetry optimization

The DESC stellarator optimization code takes advantage of advanced numerical methods to search the full parameter space much faster than conventional tools. Only a single equilibrium solution is needed at each optimization step thanks to automatic differentiation, which efficiently provides exact derivative information. A Gauss–Newton trust-region optimization method uses second-order derivative information to take large steps in parameter space and converges rapidly. With just-in-time compilation and GPU portability, high-dimensional stellarator optimization runs take orders of magnitude less computation time with DESC compared to other approaches. This paper presents the theory of the DESC fixed-boundary local optimization algorithm along with demonstrations of how to easily implement it in the code. Example quasi-symmetry optimizations are shown and compared to results from conventional tools. Three different forms of quasi-symmetry objectives are available in DESC, and their relative advantages are discussed in detail. In the examples presented, the triple product formulation yields the best optimization results in terms of minimized computation time and particle transport. This paper concludes with an explanation of how the modular code suite can be extended to accommodate other types of optimization problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2022

Linear simulation of magnetohydrodynamic plasma response to three-dimensional magnetic perturbations in high-β P plasmas

In this work, we report the numerical analyses of linear magnetohydrodynamics (MHD) plasma response to applied three-dimensional magnetic perturbations (MPs) in a joint DIII-D/EAST collaboration on high-β P (poloidal beta) plasmas, utilizing the extended-MHD code M3D-C1, with the purpose of realizing a better understanding of the existing experiment in which the n=3 MPs were applied to such high-β P plasmas attempting to control large amplitude type-I ELMs. Such high-β P plasmas obtained at the DIII-D tokamak feature an upper-biased double null configuration, a high edge safety factor q 95 ~6.4, and a stable internal transport barrier (ITB) leading to relatively high core pressures. Single-fluid simulations show that the plasma response to n=3 MPs, including both non-resonant/kinking and resonant components, is significantly weaker than that to n=1 or 2 MPs. To survey the impact of q 95 on plasma response to applied MPs, the SEGWAY (Self-consistent Equilibrium Generating Workflow for AnalYsis) module, developed in the OMFIT integrated modelling framework, is employed to generate a series of equilibria with a wide range of q 95 while other key parameters including the normalized beta, electron density at pedestal top, and plasma shape are kept fixed. Compared to the vacuum response, single-fluid M3D-C1 simulations predict a much more significant decrease of resonant plasma response to the applied n=3 MPs at the maximum penetration radii as q 95 increases. In contrast to single-fluid simulation results showing resonant penetration occurs only near the pedestal top where the E×B toroidal rotation frequency is zero, two-fluid simulations show two comparable resonant penetrations locating near the pedestal top and the ITB foot, where the perpendicular electron rotation frequency is zero. Such resonant field penetration near the ITB foot may be responsible for the observed formation of a staircase structure in both electron density and temperature profiles and thereby a considerable deterioration of global plasma performance when MPs are applied in high-β P plasmas. Motivated by this numerical work, we provide some ideas for the future research, with the purpose of realizing effective ELM control in such high-β P plasmas on the DIII-D and EAST devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Roadmap on electronic structure codes in the exascale era

Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources devoted to these calculations. As we enter the exascale era, exciting new opportunities to increase simulation numbers, sizes, and accuracies present themselves. In order to realize these promises, the community of electronic structure software developers will however first have to tackle a number of challenges pertaining to the efficient use of new architectures that will rely heavily on massive parallelism and hardware accelerators. This roadmap provides a broad overview of the state-of-the-art in electronic structure calculations and of the various new directions being pursued by the community. It covers 14 electronic structure codes, presenting their current status, their development priorities over the next five years, and their plans towards tackling the challenges and leveraging the opportunities presented by the advent of exascale computing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Grad–Shafranov equilibria via data-free physics informed neural networks

A large number of magnetohydrodynamic (MHD) equilibrium calculations are often required for uncertainty quantification, optimization, and real-time diagnostic information, making MHD equilibrium codes vital to the field of plasma physics. In this paper, we explore a method for solving the Grad–Shafranov equation by using physics-informed neural networks (PINNs). For PINNs, we optimize neural networks by directly minimizing the residual of the partial differential equation as a loss function. We show that PINNs can accurately and effectively solve the Grad–Shafranov equation with several different boundary conditions, making it more flexible than traditional solvers. This method is flexible as it does not require any mesh and basis choice, thereby streamlining the computational process. We also explore the parameter space by varying the size of the model, the learning rate, and boundary conditions to map various tradeoffs such as between reconstruction error and computational speed. Additionally, we introduce a parameterized PINN framework, expanding the input space to include variables such as pressure, aspect ratio, elongation, and triangularity in order to handle a broader range of plasma scenarios within a single network. Parameterized PINNs could be used in future work to solve inverse problems such as shape optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Characterization of self-magnetic pinch radiographic diode performance on RITS-6 at Sandia National Laboratories. II. Coupling between the inductive voltage adder and the SMP load

The self-magnetic pinch (SMP) diode is a type of radiographic diode used to generate an intense electron beam for radiographic applications. At Sandia National Laboratories, SMP was the diode load for the six-cavity radiographic integrated test stand inductive voltage adder (IVA) driver operated in a magnetically insulated transmission line (MITL). The MITL contributes a flow current in addition to the current generated within the diode itself. Extensive experiments with a MITL of 40Ω load impedance [T. J. Renk et al., Phys. Plasmas 29, 023105 (2022)] indicate that the additional flow current leads to results similar to what might be expected from a conventional high-voltage interface driver, where flow current is not present. However, when the MITL flow impedance was increased to 80 Ω, qualitatively different diode behavior was observed. This includes large retrapping waves suggestive of an initial coupling to low impedance as well as diode current decreasing with time even as the total current does not. A key observation is that the driver generates total current (flow + diode) consistent with the flow impedance of the MITL used. The case is made in this paper that the 80Ω MITL experiments detailed here can only be understood when the IVA-MITL-SMP diode is considered as a total system. Additionally, the constraint of fixed total current plus the relatively high flow impedance limits the ability of the diode (whether SMP or other type) to act as an independent load. An unexpected new result is that in tracking the behavior of the electron strike angle on the converter as a function of time, we observed that the conventional cIV x “Radiographic” radiation scaling (where x~2.2) begins to break down for voltages above 8 MV, and cubic scaling is required to recover accurate angle tracking.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A new tubular ironless linear induction motor using non-interrupted winding method

The ironless linear induction motor device can be used for any application which requires a large motion span or rotation restrictedly by the ironless feature due to the strong magnetic field environment. One of the possible applications is ITER diagnostics shutters. The tokamak reactor systems have a severe radiation issue during the operation. The traditional Kapton-insulated wires have significantly lower radiation resistance than an expected large radiation dose at ITER Motion Stark Effect (MSE) shutter location. Furthermore, the traditional single-sided linear induction motor winds the wire to make the N-turn pitch coils independently and then connect them in series to make a one-phase winding. Therefore, it is a big challenge to make so many radiation-protection connectors in such a limited space as inside the vacuum vessel of fusion reactors. The approved mineral-insulated kinds of cables need a large bend radius for the end turns, so the traditional single-sided linear motor with two-layer three-phase windings need to be upgraded to the tubular linear motors. The new proposed tubular ironless linear induction motor can not only overcome the large bend radius issue with the mineral-insulated cables but also allow the motor coils wound continuously around the tubular-shaped stator, forward and backward N times to make a N-turn winding for each phase, in order to resolve the radiation protection issue with the winding connection. The 3D transient Maxwell electromagnetic model with a tubular three-phase linear induction motor has been analyzed. Furthermore, this paper will present the detailed design and winding method of ironless linear induction motor actuators, and results of the 3D transient numerical simulation of the actuator performance.

Ironless↗