Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first wall”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Density functional theory calculation of helium segregation and decohesion effect in W110/Ni111 interphase boundary

Density functional theory (DFT) is employed to study the preferential distribution and decohesion effect of He in a W-NiFe composite consisting of W particles embedded in an Ni-based solid solution matrix. A slab containing {110}<100>W//{111}<110>Ni interface is employed as a surrogate model for the W-NiFe system. Firstly, the fracture energy of the W/Ni interphase boundary (IB) (4.37 J/m 2 ) is higher than the cleavage energy of Ni{111} (3.82 J/m^2) and lower than the cleavage energy of W{110} (6.60 J/m 2 ). The comparison shows that the cohesion of the IB is stronger than the Ni{111} planes that are away from the IB. However, the cohesion between the Ni{111} planes adjacent to the IB is found to be the weakest in this system, with a cleavage energy of 3.11 J/m 2 . Subsequently, formation energy of He is calculated in the Ni slab, W slab, and various interstitial sites in the IB. The DFT calculations show that He is significantly more stable in Ni than in W, by about 1.75 eV. Interestingly, He does not prefer to segregate at the IB as compared to bulk Ni. Nevertheless, it prefers to segregate to the region between the Ni{111} planes adjacent to the IB, and decreases the cohesion of the already weakest region. Based on an estimated amount of He gas production in 5 years under first wall neutron irradiation (neutron flux of 1.04x10 15 n/cm 2 /s), the He will decrease the cleavage energy of the weakest region by 21.2%, 15.4%, and 12.2% at 800, 1000, and 1200 C, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Thermal gradient effect on helium and self-interstitial transport in tungsten

First-wall materials in a fusion reactor are expected to withstand harsh conditions, with high heat and particle fluxes that modify the materials microstructure. These fluxes will create strong gradients of temperature and concentration of diverse species. Besides the He ash and the hydrogenic species, neutron particles generated in the fusion reaction will collide with the material creating intrinsic defects, such as vacancies, self-interstitials atoms (SIAs), and clusters of such point defects. These defects and the He atoms will then migrate in the presence of the aforementioned gradients. In this study, we use nonequilibrium molecular dynamics to analyze the transport of He and SIAs in the presence of a thermal gradient in tungsten. We observe that, in all cases, the defects and impurity atoms tend to migrate toward the hot regions of the tungsten sample. The resulting species concentration profiles are exponential distributions, rising toward the hot regions of the sample, in agreement with irreversible thermodynamics analysis. For both He atoms and SIAs, we find that the resulting species flux is directed opposite to the heat flux, indicating that species transport is governed by a Soret effect (thermal-gradient-driven diffusion) characterized by a negative heat of transport that drives species diffusion uphill (from the cooler to the hot regions of the sample). Here, we demonstrate that the steady-state species profiles obtained accounting for the Soret effect vary significantly from those where temperature-gradient-driven transport is not considered and discuss the implications of such a Soret effect on the response to plasma exposure of plasma-facing tungsten.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Evidence on the effects of main-chamber neutrals on density shoulder broadening

Evidence that density shoulder broadening is dependent on high main-chamber neutral density is presented. Shoulder broadening does not occur when the sources for main-chamber neutrals are minimized using divertor baffles and wide gaps to the first wall (~3× the density decay length). Removing the baffles or reducing the gap to the inner wall both act to increase the density shoulder amplitude in otherwise identical TCV discharges. Radial turbulent transport is correlated with shoulder amplitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Development of pebble-based extruded carbon rods for extreme plasma heat flux environments

This work presents first experiments toward the development of continuously renewable (extrudable) pebble-based carbon rods for use as plasma-facing components in extreme steady-state plasma flux environments. The primary envisioned application of this work is a first wall that can survive long-term in future magnetic fusion power reactors while also improving recovery of the reactor fuel (tritium and deuterium atoms). Bench tests applying extreme steady-state front-surface heat loads of up to 50 MW/m 2 are presented. Continuous pebble rod front-surface recession and intact pebble recovery are successfully demonstrated, at a rate of order 0.2 cm/s. Numerical simulations of the pebble rod front-surface recession are able to match observations reasonably well, indicating that the recession mechanism can be understood as occurring due to pebble thermal expansion and resulting shock and cracking of the inter-pebble binder. Tests of the pebble rod extrusion demonstrate that friction between the rods and the stainless-steel extrusion channel is tolerably low (<50 N for the expected channel length) over a wide range of temperatures. Front-surface outgassing rates below 1000 Torr L/s/m 2 are achieved, believed to be sufficiently low for use in magnetic fusion reactors. In conclusion, initial parametric scans over pebble rod size and binder fraction to vary front-surface recession rates are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Filamentary velocity scaling validation and spin dynamics in the DIII-D tokamak

Measured filament velocities in the DIII-D tokamak are compared against theoretical scalings, finding that the latter often represents an upper limit on experimental velocity distributions with most filaments possessing lower velocity. Filament spin from internal E × B drift is experimentally demonstrated to alter filament radial velocity. A critical spin velocity, where filament radial velocity peaks, is observed and corresponds to approximately 5 km/s. This transition is corroborated using a less direct measure of filament spin in the form of a temperature ratio. These techniques are combined to find that the critical spin velocity closely aligns with transport times along and across filaments becoming comparable. The normalized filament size distribution is consistent with the most stable size as dictated by Kelvin–Helmholtz and curvature-driven instabilities. Overall, the findings suggest filament stability and spin alter filamentary transport that may threaten the integrity of first walls in fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Design Overview of a High-Pressure Helium Flow Visualization Apparatus for Blanket Cooling Studies

Cooling of the fusion blanket first wall remains a significant challenge given the adverse conditions of heat and particle flux encountered near the plasma. Helium emerges as an attractive cooling candidate because of its chemical and neutronic inertness and separability from hydrogenic species (e.g. tritium). Because of the low thermal mass of helium, optimization of these coolant channels is warranted to provide high heat transfer performance at low pumping costs. Increasingly, computational fluid dynamics (CFD) simulations are employed to model and optimize these flow channels, and accompanying experimental data are needed to validate the predictions of these models. To provide the aforementioned experimental data, a high-pressure helium flow visualization upgrade has been designed for the Helium Flow Loop Experiment facility. This apparatus was built to American Society of Mechanical Engineers boiler and pressure vessel standards to withstand operating pressure of 4 MPa and mated to high-pressure glass windows. Seedless flow visualization is performed via high-speed background oriented schlieren (BOS), with image correlation used for time-resolved two-dimensional velocimetry at frequencies in excess of 60 kHz. Rectangular flow channel test articles are additively manufactured via laser powder bed fusion and installed into this visualization apparatus, with one-sided heating supplied by resistive heaters. In conclusion, the chosen test geometries were informed by prior CFD simulations, and the helium flow structures observed via BOS (detachment, recirculation, etc.) will be used for the validation of these accompanying models, in support of the design and optimization of blanket cooling channel configurations.

Helium flow↗
DOE OSTI2023

Scenario optimization for the tokamak ramp-down phase in RAPTOR: Part B. safe termination of DEMO plasmas

An optimized plasma current ramp-down strategy is critical for safe and fast termination of plasma discharges in a tokamak demonstration fusion reactor (DEMO), both in planned and emergency scenarios, avoiding plasma disruptions and excessive heat loads to the first wall. Plasma stability limits and machine-specific technical requirements constrain the stable envelope through which the plasma must be navigated. Large amounts of auxiliary heating are required throughout the ramp-down phase, to avoid a radiative collapse in the presence of intrinsic tungsten and seeded xenon impurities, as quantitatively estimated in this work. As the plasma current is reduced, the current density becomes increasingly peaked, reflected by a growing value of the internal inductance $\ell$ $i3$ , resulting in reduced controllability of the vertical position of the plasma. The feasibility of different plasma current ramp-down rates is tested by applying an automated optimization framework embedding the RAPTOR core transport solver. Optimal time traces for plasma current $I$ $p$ $(t)$ and plasma elongation $κ(t)$ are proposed, to satisfy an $I$ $p$ -dependent upper limit on the plasma internal inductance, as obtained from vertical stability studies using the CREATE-NL code, as well as a constraint on the time evolution of $q$ 95 , to avoid an ideal MHD mode. A negative current density near the plasma edge is observed in our simulations, even for the most conservative $I$ $p$ ramp-down rate, indicating significant transient dynamics due to a large resistive time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Measurement of post-disruption runaway electron kinetic energy and pitch angle during final loss instability in DIII-D

Post-disruption runaway electron (RE) kinetic energy K and pitch angle sin$\vartheta$ are critical parameters for determining resulting first wall material damage during wall strikes, but are very challenging to measure experimentally. During the final loss instability, confined RE K and sin$\vartheta$ are reconstructed during center-post wall strikes for both high impurity (high-Z) and low impurity (low-Z) plasmas by combining soft x-ray, hard x-ray, synchrotron emission, and total radiated power measurements. Deconfined (wall impacting) RE sin$\vartheta$ is then reconstructed for these shots by using time-decay analysis of infra-red imaging. Additionally, deconfined RE K and sin$\vartheta$ are reconstructed for a low-Z downward loss shot by analyzing resulting damage to a sacrificial graphite dome limiter. The damage analysis uses multi-step modeling simulating plasma instability, RE loss orbits, energy deposition, and finally material expansion (MARS-F, KORC, GEANT-4, and finally COMSOL). Overall, mean kinetic energies are found to be in the range $\langle$$K$$\rangle$ $≈ 3-4$ MeV for confined REs. KORC simulations indicate that the final loss instability process does not change individual RE kinetic energy K. In conclusion, confined RE pitch angles are found to be fairly low initially pre-instability, $\langle$sin$\vartheta$$\rangle$ $≈ 0.1 – 0.2$, but appear to increase roughly 2x, to $\langle$sin$\vartheta$$\rangle$ $≈ 0.3 – 0.4$ for both confined and deconfined REs during instability onset in the low-Z case; this increase is not observed in the high-Z case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Routes to high-performance operation in Wendelstein 7-X: turbulence suppression with shaping of the density profile

Steep density gradients generally lead to improved plasma performance in the neoclassically optimized stellarator W7-X. This is evident in the global energy confinement time as well as in the ion temperature and can be explained by a strong reduction of the ion temperature gradient turbulence. Such conditions can experimentally be realized by several methods: injection of cryogenic hydrogen pellets, appropriate combination of neutral beam and electron cyclotron resonance heating (ECRH) and, in some cases, with low power of ECRH after preconditioning of the first wall. The duration of the improved phases is determined by the ability to sustain the steep density gradient, by technical limitations of the involved systems and, eventually, by the plasma stability. This paper gives an overview of relevant experimental results and presents example discharges where the improved confinement conditions could be extended to multiple seconds: up to 4 s using neutral beam injection and from 14 to 40 s with steady state pellet injection. In these plasmas the turbulent thermal diffusivity is reduced by a factor of 3 to 4 in a broad radial range, which allows high ion temperatures of up to 3 keV at the densities of about 1.5 • 10 20 m −3 .

high performance plasma↗
DOE OSTI2026

Post-puff SOL broadening on MAST-U under high-recycling conditions: evidence consistent with cross-field transport changes

Transient broadening of the scrape-off layer (SOL) density profile can modify main-chamber first-wall particle fluxes and divertor loading, yet its control parameters remain debated between divertor-regime transitions, neutral dynamics and changes in cross-field transport. We investigate fueling-driven SOL density-profile evolution and post-fueling relaxation on MAST-Upgrade (MAST-U) in ohmic L-mode, using two otherwise similar double-null Conventional Divertor discharges (I p = 450 kA, B T = 0.33 T) with identical 50 ms low-field side gas puffs; the only intentional difference is the puff start time. Upstream Thomson scattering shows that both discharges develop a transient far-SOL density profile modification, expressed as an increased SOL-width metric λ n e and a far-SOL enhancement consistent with a shoulder-like signature. In the earlier-puff case, the SOL-width metric remains elevated after puff termination and the post-puff decay time scales are systematically longer across the analysed radii. Outer-target Langmuir probes indicate high-recycling conditions during the analysed window (few-eV T e with radially peaked j sat,∥ and q ∥ profiles, without signatures of deep detachment). The outer-divertor collisionality proxy Λ div is elevated in both cases and does not discriminate between the different post-puff persistence. The target-integrated ion flux proxy ∫ J sat dA evolves nearly identically in both discharges when aligned to the puff start. Taken together, these observations suggest that the late/post-puff upstream SOL evolution is not set by parallel exhaust to the outer targetalone and is more consistent with upstream cross-field redistribution, with a possible role for plasma–neutral coupling and fueling geometry.

MAST-U↗
DOE OSTI2023

Initial study on thermal stability of cold spray tantalum coating irradiated with deuterium for fusion applications

Abstract Removal of neutral hydrogen atoms in the plasma edge reduces the number of charge exchange events and thus, the net energy losses in the plasma, significantly improving performance of fusion devices. Effective control of the residual pressure of hydrogen isotopes (HIs) in the plasma edge may be achieved by utilizing a hydrogen absorbing first wall interface capable of withstanding the harsh fusion environment. In this study, we have investigated tantalum (Ta) coating deposited by cold spray technology on 316L stainless steel substrate as a potential plasma-facing material surface. High fluence low energy deuterium plasma irradiation experiments and subsequent thermal annealing cycles associated with thermal desorption spectrometry (TDS) demonstrated superior structural stability of the Ta coating. TDS experiments revealed the outgassing of deuterium (as measure of its retention) for cold spray Ta coatings to be three times higher than bulk Ta and two orders of magnitude greater than bulk polycrystalline W. X-ray photoelectron spectroscopy revealed evolution of oxidation states upon deuterium irradiation and a partial recovery of the metallic signature of Ta after the thermal treatment at 1100 K.

Physics↗
DOE OSTI2020

Simulation of MHD instabilities with fluid runaway electron model in M3D- C 1

Runaway electrons are generated in a tokamak during the start up, during normal operation and during a plasma disruption. During a disruption, runaway electrons can be accelerated to high energies, potentially damaging the first wall. To predict the consequences of runaway generation during a disruption, it is necessary to consider resonant interactions of runaways with the bulk plasma. Here we consider the interactions of runaways on low mode number tearing modes. We have developed a fluid runaway electron model for the 3D MHD code M3D-C1[Jardin,et al. J Comput. Sci Discovery 6 014002 (2012)]. To benchmark, we have reproduced the MHD linear tearing mode results (with runaway electrons) in a circular cylinder presented in previous analytic studies[[Helander, P., et al, Phys. Plasmas 14 144102 (2007)] and extended here with a numerical eigenvalue calculation. Furthermore, we find that the low mode number tearing mode has a rotation caused by the MHD - runaways interaction, and the toroidal current scale length is much smaller with runaways than that for without and decreases as the runaway speed increases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An atom probe tomography and inventory calculation examination of second phase precipitates in neutron irradiated single crystal tungsten

Tungsten is the prime candidate material for the plasma facing divertor system and first wall armor of future nuclear fusion reactors. However, the understanding of the microstructural and chemical evolution of pure tungsten under neutron irradiation is relatively unknown, in part due to a lack of experimental data on this topic. Here, single crystal tungsten has been irradiated in the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory at a nominal temperature of 700–800 °C to damage levels of 0.1 and 1.8 displacements per atom (dpa). Inventory calculations of the neutron irradiation experiments have been used to track isotope generation and decay to inform atom probe tomography (APT) results in the determination of the transmutation-induced precipitate compositions. Furthermore, APT crystallography has been used to show the relationship between precipitates and matrix. The composition of the precipitates is shown to progress towards that of the σ-phase at the highest dose studied, with those precipitates lying along crystallographic planes similar to those of displacement-induced dislocations. This work also sets the framework for APT studies of materials that contain isotopic ratios far from those observed in the natural state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Beryllium-driven structural evolution at the divertor surface

Erosion of the beryllium first wall material in tokamak reactors has been shown to result in transport and deposition on the tungsten divertor. Experimental studies of beryllium implantation in tungsten indicate that mixed W–Be intermetallic deposits can form, which have lower melting temperatures than tungsten and can trap tritium at higher rates. To better understand the formation and growth rate of these intermetallics, we performed cumulative molecular dynamics (MD) simulations of both high and low energy beryllium deposition in tungsten. In both cases, a W–Be mixed material layer (MML) emerged at the surface within several nanoseconds, either through energetic implantation or a thermally-activated exchange mechanism, respectively. While some ordering of the material into intermetallics occurred, fully ordered structures did not emerge from the deposition simulations. Targeted MD simulations of the MML to further study the rate of Be diffusion and intermetallic growth rates indicate that for both cases, the gradual re-structuring of the material into an ordered intermetallic layer is beyond accessible MD time scales(≤1 μs). However, the rapid formation of the MML within nanoseconds indicates that beryllium deposition can influence other plasma species interactions at the surface and begin to alter the tungsten material properties. Therefore, beryllium deposition on the divertor surface, even in small amounts, is likely to cause significant changes in plasma-surface interactions and will need to be considered in future studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

3D radiation, density, and MHD structures following neon shattered pellet injection into stable DIII-D Super H-mode discharges

Six nominally repeat neon shattered pellet injection (SPI) shutdowns of stable DIII-D Super H-modes are studied to understand the 3D properties of the radiation and impurity transport. The radiation efficiency and radiation peaking determine whether first wall melting is expected following disruption mitigation in ITER. Previous studies make use of axisymmetric approximations to infer radiation efficiencies, but validating the high efficiency required by ITER necessitates improved accuracy, and this work contributes by exploring the 3D radiation and density structures that will inform forward modeling. When the neon shatter plume produced by the SPI reaches the plasma edge, m/n = 3/1 and 2/1 island O-points are observed to align with the injection trajectory in five out of six cases, suggesting that the injected material seeds the island O-points. Field aligned neon structures emitting Ne-I line radiation drift at 1 km/s in the ion diamagnetic drift direction during the pre-TQ, tracking the motion of the m/n = 2/1 island O-point. Neon fragments penetrate to the q = 2 surface by the time of the TQ. Techniques to constrain the 3D emissivity are explored, and one method constrains a 3D flux tube that is consistent with the radiation data, and when mapped to the interferometers, intersects the lasers that measure the highest density. The resulting structure derived from the radiation measurements exists near the 2/1 island X-point. In five repeatable discharges, the peak of the radiation in the toroidal direction exists in a 120° toroidal sector where the injection occurs, in contrast with the outlier discharge where the toroidal peak exists in the complementary 240° toroidal sector far from the injector, and where a 50% lower density rise is observed. In conclusion, the n = 1 phase behavior is markedly different in the outlier discharge, suggesting a possible dependence of the radiation structure and the assimilation efficiency on MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Collisionless losses of fast ions in the Divertor Tokamak Test due to toroidal field ripple

In this paper we analyze fast ion motion in the Divertor Tokamak Test device. DTT is planned to be heated through a mix of 45 MW heating power, including 15 MW negative-ion-based Neutral Beam heating (NNBI) which is currently being developed by Consorzio RFX in Padova, Italy. An issue for DTT is that a Toroidal Field (TF) ripple with a maximum value of about 0.42% (with respect to the on-axis magnetic field $B_0$) is expected on the low-field side (LFS), and this ripple interacts with fast ions through rather well-known phenomena of ripple-precession resonances, in addition to prompt losses of ions which do not complete a full orbit in the poloidal plane. We will show that, with the planned geometry of NNBI, prompt losses are negligible, and ripple-precession losses amount to a maximum 0.15%. The calculations are performed with the guiding center code ORBIT using two different equilibria, and a beam with energy 400 keV and injection angle $\alpha_{inj} = 40^{\circ}$ (measured w.r.t. the first wall), which corresponds to a pitch of injected particles $\lambda = v_{\parallel}/v \approx \sin{\alpha_{inj}} = 0.65$. Main resonances are of the form $\omega_b - n N \omega_d = 0$, $\omega_b$ and $\omega_d$ being the bounce and precession frequency, respectively, $N=18$ the ripple periodicity and $3 \le n \le 6$ are the toroidal wavenumbers of the resonances. Although collisionless interaction with the TF ripple does not pose a serious threat to the NNBI project, an open question remains whether the presence of these resonances will interact with fast particles accelerated by \alfven eigenmodes, and if stochastization of the resonances is possible in DTT, as it was observed in the past in TORE SUPRA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Self-consistent core-pedestal ITER scenario modeling

The purpose of these integrated ITER simulations is to identify dependencies that can impact the performance of ITER. The dependence of fusion power production, temperature and density pedestals, and core profiles on varying magnetic-q, edge density fueling strength, neoclassical transport, magnetic field strength, and alpha heating, are examined. It is noted that ion modes usually dominate the interior, while electron modes dominate the barrier. The electron mode without flow shear is substantially stronger. The q increase tends to cause more kinetic ballooning modes at the edge and the q decrease means that the peeling is getting worse, but the ballooning is getting better. The large edge q tends to provide hollow density while small edge q gives normal density profile. This is consistent with the observation that the density of particles is inversely proportional to the density of the plasma current. The rise in the source of the edge particle increases the density of the edge and the reaction rate close to the edge increases temperature fluxes, resulting in weak pedestal barriers. When neoclassical transport is turned off or the strength of B-field is increased, the pronounced edge barriers are identified. In addition, a distinct effect on the magnetohydrodynamic (MHD) modes is observed due to the increased B-field. The interesting aspect is that the barrier collapses about q = 2 as the stability of MHD has decreased to replace the barrier with a continuous slope but with roughly the same fusion power. This would actually be good for a reactor. The slope of the H-mode pedestal is found to be reduced due to the alpha heating. It is expected that the problem of edge-localized modes causing damage to the first wall will be eased. The effects of circular, elongated, and general geometry on plasma profiles are also compared and contrasted. The circular geometry has a hollow density profile and no edge or internal thermal transport barrier (ETB or ITB). However, both ETBs and ITBs can be found in the elongated geometry. In general geometry, there is a weaker ETB, comparable ITBs, but a higher edge particle transport barrier than in elongated geometry. Higher density near the edge, on the other hand, causes more wall erosion, so elongated geometry might be the best option.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Improved scaling law for the prediction of deuterium retention in beryllium co-deposits

Hydrogen isotope co-deposition with Be eroded from the first wall is expected to be the main fusion fuel retention mechanism in ITER. Since good fuel accounting is crucial for economic and safety reasons, reliable predictions of hydrogen isotope retention are needed. Furthermore, this study builds upon the well-established empirical De Temmerman scaling law that predicts D/Be ratios in co-deposited layers based on deposition temperature, deposition rate, and deuterium particle energy. Expanding the data used in the original development of the scaling law with an additional dataset obtained with more recent measurements using a different technique to the original De Temmerman approach, allows us to obtain new values for free parameters and improve the prediction capabilities of the new scaling law. In an effort to improve the model even further, scaling with D 2 background pressure was included and a new two-term model derived, describing D retention in low- and high-energy traps separately.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗