Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PLASMA CONFINEMENT”

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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Progress towards a practical multicell positron trap

Described here is progress in an experimental program to develop a 21 cell multicell trap for the accumulation and storage of {approx} 10{sup 12} positrons. The basic architecture is an arrangement of multiple Penning-Malmberg (PM) trapped plasmas (i.e., cells) arranged in parallel in a common vacuum system and magnetic field. Experiments are described that are intended to address several key issues, including the effects of large space charge potentials and high plasma densities on: plasma heating, deterioration of confinement, and decreased efficiency of rotating electric fields in producing plasma compression. The confinement of PM plasmas displaced both radially and toward the ends of the uniform magnetic field region will also be investigated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Understanding and controlling capsule symmetry in near vacuum hohlraums at the National Ignition Facility

The near vacuum hohlraum platform is an inertial confinement fusion design at the National Ignition Facility (NIF) that uses the lowest practical density of helium gas of 30 μg/cc to fill the hohlraum, which is ten times lower than now used routinely. This has several advantages, such as high laser coupling; however, the inability to understand and simulate the symmetry of the imploded capsule has limited the use of this platform. This work presents the first simulations that are able to accurately capture the highly prolate implosion seen experimentally without unphysical, ad hoc model changes. While previous investigations attributed this asymmetry to multi-species interpenetration in the hohlraum, we find that this alone has little effect on symmetry. Instead, it is the presence of crossed-beam energy transfer (CBET), occurring with no applied wavelength shift between the laser beams, that increases the laser power to the inner cones and causes a more prolate implosion. The effect of CBET is increased in the simulation model when the hohlraum laser entrance hole hardware is included. Using this understanding, CBET is exploited by shifting the inner-beam wavelength by –0.75 Å (at 1ω) with respect to the outer-beams. This transfers laser power to the outer-beams in contrast to positive wavelength shifts as done routinely on NIF and produces a round capsule implosion in our simulations. Furthermore, this work shows the possibility of the near vacuum hohlraum as a viable experimental platform.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An experimental and computational study of thin-layer Rayleigh–Taylor instability development during deceleration with and without an externally applied magnetic field

The importance of mitigating the Rayleigh–Taylor instability (RTI) in inertial confinement fusion (ICF) is critical to successfully achieve high gain fusion yield. Consequently, understanding the seed mechanisms of RTI and the potential evolution of RTI in ICF relevant conditions is crucial. Single feature perturbations consistently demonstrate non-linear RTI evolution, for which an experimental platform on OMEGA-EP is developed. Manufacturing defects introduced into the target design require exploration of unanticipated changes to RTI development and an identification of targets that will still render quantifiable physics results. Consequently, it is presented that the inherent 3D nature of experimental targets necessitates 3D modeling for accurate design work and predictive modeling of experimental targets, especially when high resolution imaging diagnostics, like Fresnel Zone Plates, are utilized. A study of the morphology of the RTI evolution due to changing initial conditions and the presence of an externally applied magnetic field are also explored. Experimental data show thin-layer RTI morphology comparable to resistive magneto hydrodynamic 3D results. A discussion on the impacts of an externally applied magnetic field makes the case for continued efforts to observe a magnetic field's impact on RTI morphology.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Local measurements of the pedestal magnetic field profile throughout the ELM cycle on DIII-D

We report new high speed localized measurements of the pedestal magnetic field during the edge localized mode (ELM) cycle of a DIII-D High confinement mode (H-mode) discharge indicate a temporally and spatial complex redistribution of the edge current density profile, $j_{edge}$. The measurement technique extracts the magnetic field magnitude, $\textit{B}$, via the spectral separation of Stark-split neutral beam radiation in the pedestal. Single spatial channel measurements from a novel spatial heterodyne spectrometer are validated in discharges with core current profile changes. The technique measures Stark-splitting changes that imply $\textit{B}$ changes as small as 1 mT with high time resolution (50 μs). At normalized poloidal flux $ψ_n$ = 1.0, $\textit{B}$ appears saturated in the inter-ELM period and then rapidly decreases in <200 μs by ~1%, before edge recycling emission begins to increase. Radially inboard of $j_{edge}$, $\textit{B}$ increases at the ELM crash. The behavior is consistent with a rapid collapse of $j_{edge}$ at the ELM crash and subsequent pedestal recovery. In some discharges, at $ψ_n$ < 0.96, changes in $\textit{B}$ are observed throughout the ELM cycle. In others, $\textit{B}$ recovers and is relatively stable until a few ms leading up to the next crash. Measurements of $\textit{B}$ during the H-mode transition show a large increase at $ψ_n$ = 1 with little change at $ψ_n$ = 0.9, consistent with the formation of the edge bootstrap current density peak. The $ψ_n$ = 0.9 spectrum is complicated by predicted changes to the Stark component intensities with density at the L–H transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

An adjoint method for neoclassical stellarator optimization

Stellarators remain an enticing route to steady-state fusion power. However, to achieve the required confinement, the magnetic geometry must be highly optimized. This optimization requires navigating high-dimensional spaces, often necessitating the use of gradient-based methods. The gradient of the neoclassical fluxes is expensive to compute with classical methods, requiring $O(N)$ flux computations, where $N$ is the number of parameters. To reduce the cost of the gradient computation, we present an adjoint method for computing the derivatives of moments of the neoclassical distribution function for stellarator optimization. The linear adjoint method allows derivatives of quantities which depend on solutions of a linear system, such as moments of the distribution function, to be computed with respect to many parameters from the solution of only two linear systems. This reduces the cost of computing the gradient to the point that the finite-collisionality neoclassical fluxes can be used within an optimization loop. With the neoclassical adjoint method, we compute solutions of the drift kinetic equation and an adjoint drift kinetic equation to obtain derivatives of neoclassical quantities with respect to geometric parameters. When the number of parameters in the derivative is large ( $O(10^{2})$ ), this adjoint method provides up to a factor of 200 reduction in cost. We demonstrate adjoint-based optimization of the field strength to obtain minimal bootstrap current on a surface. With adjoint-based derivatives, we also compute the local sensitivity to magnetic perturbations on a flux surface and identify regions where tight tolerances on error fields are required for control of the bootstrap current or radial transport. Moreover, the solve for the ambipolar electric field is accelerated using a Newton method with derivatives obtained from the adjoint method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Application of cross-beam energy transfer to control drive symmetry in ICF implosions in low gas fill Hohlraums at the National Ignition Facility

Cross beam energy transfer (CBET), invoked by setting a wavelength difference, Δλ, between inner and outer beam cones, can be used to increase the drive on the waist in indirectly driven inertial confinement fusion experiments at the National Ignition Facility (NIF). Historically, hot spot symmetry control in capsule implosions in high (≥0.9 mg/cm 3 4 He) gas fill Hohlraums was enabled by substantial CBET. However, these implosion designs suffered from inflight symmetry swings, high SRS backscatter on the inner cones, and significant hot electron generation posing a threat to DT fuel preheat. Subsequent experiments in larger, low (≤0.6 mg/cm 3 4 He) gas fill Hohlraums demonstrated round implosions by varying the inner cone fraction throughout the laser drive at Δλ = 0 Å while keeping backscatter and hot electron generation very low. To enable driving larger capsules at a given Hohlraum size, additional tools for implosion symmetry control are required. With this goal in mind, here we present a detailed experimental study of using CBET in low gas fill Hohlraums near NIF's current peak power capability. We find a ~2.5× higher sensitivity of the P 2 Legendre mode with respect to Δλ changes compared to that of high gas fill designs. We attribute this observation to the fact that backscatter remains very low and that CBET remains in a linear regime, as suggested by simulations. As a result, a much smaller Δλ of order 1 Å is sufficient for sustaining implosion symmetry while keeping laser-to- Hohlraum coupling high and hot electron generation very low. While this study used plastic ablator capsules, our findings can be generalized to other ablator materials and, hence, show great promise for using wavelength detuning as a strong lever for implosion symmetry control in future low gas fill designs that require smaller case to capsule ratios in order to increase the energy coupled to the capsule.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An experimentally informed design process for future inertial confinement fusion facilities

The achievement of ignition in the laboratory has renewed interest in defining the requirements for a future high-gain inertial confinement fusion (ICF) facility. Our best chance of predicting future ICF performance is with 3-D radiation hydrodynamic simulations that have been benchmarked against experimental data, but their high computational cost is prohibitive for use in practical design studies. We introduce a hierarchical approach where 3-D simulations are tuned to match experimental measurements and used to train 3-D degradation models in 1-D simulations allowing for accurate predictions over the entire OMEGA direct-drive database. A genetic algorithm was used in combination with the trained 1-D simulations to search for optimal direct-drive implosion designs at driver energies ranging from 20 kJ to 10 MJ. As the fidelity of 3-D codes improves, this approach will provide a viable experimentally informed tool for defining the next ICF facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Current density distributions and sputter marks in electron cyclotron resonance ion sources

Most electron cyclotron resonance ion sources use hexapolar magnetic fields for the radial confinement of the plasma. The geometry of this magnetic structure is then-induced by charged particles-mapped onto the inner side of the plasma electrode via sputtering and deposition. The resulting structures usually show two different patterns: a sharp triangular one in the central region which in some cases is even sputtered deep into the material (referred to as thin groove or sharp structure), and a blurred but still triangular-like one in the surroundings (referred to as broad halo). Therefore, both patterns seem to have different sources. To investigate their origins we replaced the standard plasma electrode by a custom-built plasma electrode acting as a planar, multi-segment current-detector. For different biased disc voltages, detector positions, and source biases (referred to the detector) we measured the electrical current density distributions in the plane of the plasma electrode. The results show a strong and sharply confined electron population with triangular shape surrounded by less intense and spatially less confined ions. Observed sputter- and deposition marks are related to the analysis of the results. Our measurements suggest that the two different patterns (thin and broad) indeed originate from different particle populations. The thin structures seem to be caused by the hot electron population while the broad marks seem to stem from the medium to highly charged ions. In this paper we present our measurements together with theoretical considerations and substantiate the conclusions drawn above. The validity of these results is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2011

Power loss of a single electron charge distribution confined in a quantum plasma

The dielectric tensor for a quantum plasma is derived by using a linearized quantum hydrodynamic theory. The wave functions for a nanostructure bound system have been investigated. Finally, the power loss for an oscillating charge distribution of a mixed state will be calculated, using the dielectric function formalism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Numerical investigation of edge plasma phenomena in an enhanced D-alpha discharge at Alcator C-Mod: Parallel heat flux and quasi-coherent edge oscillations

Reduced-model scrape-off layer turbulence (SOLT) simulations of an enhanced D-alpha (EDA) H-mode shot observed in the Alcator C-Mod tokamak were conducted to compare with observed variations in the scrape-off-layer (SOL) width of the parallel heat flux profile. In particular, the role of the competition between sheath- and conduction-limited parallel heat fluxes in determining that width was studied for the turbulent SOL plasma that emerged from the simulations. The SOL width decreases with increasing input power and with increasing separatrix temperature in both the experiment and the simulation, consistent with the strong temperature dependence of the parallel heat flux in balance with the perpendicular transport by turbulence and blobs. The particularly strong temperature dependence observed in the case analyzed is attributed to the fact that these simulations produce SOL plasmas which are in the conduction-limited regime for the parallel heat flux. A persistent quasi-coherent (QC) mode dominates the SOLT simulations and bears considerable resemblance to the QC mode observed in C-Mod EDA operation. The SOLT QC mode consists of nonlinearly saturated wave-fronts located just inside the separatrix that are convected poloidally by the mean flow, continuously transporting particles and energy and intermittently emitting blobs into the SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Simulations of drift resistive ballooning L-mode turbulence in the edge plasma of the DIII-D tokamak

Results from simulations of electromagnetic drift-resistive ballooning turbulence for tokamak edge turbulence in realistic single-null geometry are reported. The calculations are undertaken with the BOUT three-dimensional fluid code that solves Braginskii-based fluid equations [X. Q. Xu and R. H. Cohen, Contrib. Plasma Phys. 36, 158 (1998)]. The simulation setup models L-mode edge plasma parameters in the actual magnetic geometry of the DIII-D tokamak [J. L. Luxon et al., Fusion Sci. Technol. 48, 807 (2002)]. The computations track the development of drift-resistive ballooning turbulence in the edge region to saturation. Fluctuation amplitudes, fluctuation spectra, and particle and thermal fluxes are compared to experimental data near the outer midplane from Langmuir probe and beam-emission-spectroscopy for a few well-characterized L-mode discharges in DIII-D. The simulations are comprised of a suite of runs in which the physics model is varied to include more fluid fields and physics terms. The simulations yield results for fluctuation amplitudes, correlation lengths, particle and energy fluxes, and diffusivities that agree with measurements within an order of magnitude and within factors of 2 or better for some of the data. The agreement of the simulations with the experimental measurements varies with respect to including more physics in the model equations within the suite of models investigated. The simulations show stabilizing effects of sheared E × B poloidal rotation (imposed zonal flow) and of lower edge electron temperature and density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Validation of transport models using additive flux minimization technique

A new additive flux minimization technique is proposed for carrying out the verification and validation (V and V) of anomalous transport models. In this approach, the plasma profiles are computed in time dependent predictive simulations in which an additional effective diffusivity is varied. The goal is to obtain an optimal match between the computed and experimental profile. This new technique has several advantages over traditional V and V methods for transport models in tokamaks and takes advantage of uncertainty quantification methods developed by the applied math community. As a demonstration of its efficiency, the technique is applied to the hypothesis that the paleoclassical density transport dominates in the plasma edge region in DIII-D tokamak discharges. A simplified version of the paleoclassical model that utilizes the Spitzer resistivity for the parallel neoclassical resistivity and neglects the trapped particle effects is tested in this paper. It is shown that a contribution to density transport, in addition to the paleoclassical density transport, is needed in order to describe the experimental profiles. It is found that more additional diffusivity is needed at the top of the H-mode pedestal, and almost no additional diffusivity is needed at the pedestal bottom. The implementation of this V and V technique uses the FACETS::Core transport solver and the DAKOTA toolkit for design optimization and uncertainty quantification. The FACETS::Core solver is used for advancing the plasma density profiles. The DAKOTA toolkit is used for the optimization of plasma profiles and the computation of the additional diffusivity that is required for the predicted density profile to match the experimental profile.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Resistive wall mode active control physics design for KSTAR

As KSTAR H-mode operation approaches the region where the resistive wall mode (RWM) can be unstable, an important issue for future long pulse, high beta plasma operation is to evaluate RWM active feedback control performance using a planned active/passive RWM stabilization system on the device. In particular, an optimal design of feedback sensors allows mode stabilization up to the highest achievable β{sub N} close to the ideal with-wall limit, β{sub N}{sup wall}, with reduced control power requirements. The computed ideal n = 1 mode structure from the DCON code has been input to the VALEN-3D code to calculate the projected performance of an active RWM control system in the KSTAR three-dimensional conducting structure device geometry. Control performance with the midplane locked mode detection sensors, off-midplane saddle loops, and magnetic pickup coils is examined. The midplane sensors measuring the radial component of the mode perturbation is found to be strongly affected by the wall eddy current. The off-axis saddle loops with proper compensation of the prompt applied field are computed to provide stabilization at β{sub N} up to 86% of β{sub N}{sup wall} but the low RWM amplitude computed in the off-axis regions near the sensors can produce a low signal-to-noise ratio. The required control power and bandwidth are also estimated with varied noise levels in the feedback sensors. Further improvements have been explored by examining a new RWM sensor design motivated by the off-midplane poloidal magnetic field sensors in NSTX. The new sensors mounted off of the copper passive stabilizer plates near the device midplane show a clear advantage in control performance corresponding to achieving 99% of β{sub N}{sup wall} without the need of compensation of the prompt field. The result shows a significant improvement of RWM feedback stabilization using the new sensor set which motivates a future feedback sensor upgrade.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

The importance of matched poloidal spectra to error field correction in DIII-D

Optimal error field correction (EFC) is thought to be achieved when coupling to the least-stable “dominant” mode of the plasma is nulled at each toroidal mode number (n). The limit of this picture is tested in the DIII-D tokamak by applying superpositions of in- and ex-vessel coil set n = 1 fields calculated to be fully orthogonal to the n = 1 dominant mode. In co-rotating H-mode and low-density Ohmic scenarios, the plasma is found to be, respectively, 7× and 20× less sensitive to the orthogonal field as compared to the in-vessel coil set field. For the scenarios investigated, any geometry of EFC coil can thus recover a strong majority of the detrimental effect introduced by the n = 1 error field. Despite low sensitivity to the orthogonal field, its optimization in H-mode is shown to be consistent with minimizing the neoclassical toroidal viscosity torque and not the higher-order n = 1 mode coupling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Integrated modeling of temperature profiles in L-mode tokamak discharges

Simulations of doublet III-D, the joint European tokamak, and the tokamak fusion test reactor L-mode tokamak plasmas are carried out using the PTRANSP predictive integrated modeling code. The simulation and experimental temperature profiles are compared. The time evolved temperature profiles are computed utilizing the Multi-Mode anomalous transport model version 7.1 (MMM7.1) which includes transport associated with drift-resistive-inertial ballooning modes (the DRIBM model [T. Rafiq et al., Phys. Plasmas 17, 082511 (2010)]). The tokamak discharges considered involved a broad range of conditions including scans over gyroradius, ITER like current ramp-up, with and without neon impurity injection, collisionality, and low and high plasma current. The comparison of simulation and experimental temperature profiles for the discharges considered is shown for the radial range from the magnetic axis to the last closed flux surface. The regions where various modes in the Multi-Mode model contribute to transport are illustrated. In the simulations carried out using the MMM7.1 model it is found that: The drift-resistive-inertial ballooning modes contribute to the anomalous transport primarily near the edge of the plasma; transport associated with the ion temperature gradient and trapped electron modes contribute in the core region but decrease in the region of the plasma boundary; and neoclassical ion thermal transport contributes mainly near the center of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

The targeted heating and current drive applications for the ITER electron cyclotron system

We present a 24 MW Electron Cyclotron (EC) system operating at 170 GHz and 3600 s pulse length is to be installed on ITER. The EC plant shall deliver 20 MW of this power to the plasma for Heating and Current Drive (H&CD) applications. The EC system is designed for plasma initiation, central heating, current drive, current profile tailoring, and Magneto-hydrodynamic control (in particular, sawteeth and Neo-classical Tearing Mode) in the flat-top phase of the plasma. A preliminary design review was performed in 2012, which identified a need for extended application of the EC system to the plasma ramp-up, flattop, and ramp down phases of ITER plasma pulse. The various functionalities are prioritized based on those applications, which can be uniquely addressed with the EC system in contrast to other H&CD systems. An initial attempt has been developed at prioritizing the allocated H&CD applications for the three scenarios envisioned: ELMy H-mode (15 MA), Hybrid (~12 MA), and Advanced (~9 MA) scenarios. In conclusion, this leads to the finalization of the design requirements for the EC sub-systems.

42 ENGINEERING↗
DOE OSTI2015

Fusion power production in International Thermonuclear Experimental Reactor baseline H-mode scenarios

Self-consistent simulations of 15 MA ITER H-mode DT scenarios, from ramp-up through flat-top, are carried out. Electron and ion temperatures, toroidal angular frequency, and currents are evolved, in simulations carried out using the predictive TRANSPort and integrated modeling code starting with initial profiles and equilibria obtained from tokamak simulation code studies. Studies are carried out examining the dependence and sensitivity of fusion power production on electron density, argon impurity concentration, choice of radio frequency heating, pedestal temperature without and with E × B flow shear effects included, and the degree of plasma rotation. The goal of these whole-device ITER simulations is to identify dependencies that might impact ITER fusion performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

X-point position dependence of edge intrinsic toroidal rotation on the Tokamak à Configuration Variable

Recent theoretical work predicts intrinsic toroidal rotation in the tokamak edge to depend strongly on the normalized major radial position of the X-point. With this motivation, we conducted a series of Ohmic L-mode shots on the Tokamak à Configuration Variable, moving the X-point from the inboard to the outboard edge of the last closed flux surface in both lower and upper single null configurations. The edge toroidal rotation evolved from strongly co-current for an inboard X-point to either vanishing or counter-current for an outboard X-point, in agreement with the theoretical expectations. The whole rotation profile shifted roughly rigidly with the edge rotation, resulting in variation of the peak core rotation by more than a factor of two. Core rotation reversals had little effect on the edge rotation. Edge rotation was slightly more counter-current for unfavorable than favorable ∇B drift discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗