Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPERIMENTAL REACTORS”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Model-free stabilization via Extremum Seeking using a cost neural estimator

In this paper, a fully model-free architecture for vertical stabilization of thermonuclear plasmas in tokamak experimental reactors is presented. For the first time, an Extremum Seeking control algorithm is combined with neural networks to estimate the Lyapunov function to be minimized, resulting in a fully data-driven control architecture. The performance of different neural networks are compared. Specifically, Multilayer Perceptrons and Extreme Learning Machines are considered. The proposed architecture is tested in simulation to show that it can counteract relevant plasma disturbances, resulting in a significant improvement in terms of the achievable operative space compared to the Extremum Seeking algorithm, which still relies on model-based cost estimator.

42 ENGINEERING↗
DOE OSTI2022

Demonstration of RELAP5-3D for transient analysis of a dual coolant lead lithium fusion blanket concept

Our research aims to build a multiphysics framework for transient analysis of the Dual- Coolant Lead-Lithium (DCLL) blanket design of the proposed Fusion Nuclear Science Facility (FNSF) for all Onticipated Operational Occurrences (AOOs) using Reduced Order Models (ROMs). This paper presents an assessment of using RELAP5-3D for transient thermal-hydraulic blanket analysis of the DCLL design. Preliminary validation studies for the RELAP5-3D properties were conducted using a representative vertical flow loop model and a comparison study using experimental data from the thermal convection corrosion loops at ORNL. These studies inspected the thermal hydraulic response of Lead-Lithium Eutectic (PbLi) systems within the RELAP5-3D code; mainly pressure drop and heat transfer, as compared to analytical calculations and experimental data. The Magnetohydrodynamic (MHD) pressure drop effect for liquid metal under the effects of a uniform magnetic field was implemented into our models using an equivalent forms loss coefficient based on the phasic momentum equations. We verified our implementation of MHD pressure drop within RELAP5-3D using a systematic approach consisting of two studies. Models were developed to perform comparison studies based on experimental data from the Argonne Liquid metal EXperiment (ALEX) facility and the Magnetohydrodynamic PbLi Experiment (MaPLE) facility. Using the developed validation basis, we built a simplified DCLL channel model within RELAP5-3D, inspected MHD pressure drop and performed transient temperature analysis. The DCLL channel model used radial heating data from the International Thermonuclear Experimental Reactor (ITER) which was normalized to the FNSF design. Here, a representative startup transient was performed using a relative power curve from the initial ramp up of a long-term power pulse of the Tore Supra reactor. This study demonstrated that RELAP5-3D ROMs are capable of accurate transient analysis of the DCLL blanket. We have since utilized MCNP to develop heating profiles based on the FNSF design and are working toward implementing them into our DCLL model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Estimation of the Thermal-Fluids and Thermal-Structural Performance of Helium-Cooled Modular Finger-Type Divertors

Over the past decade, our group has investigated the thermal-fluid performance of the helium-cooled modular divertor with multiple jets (HEMJ) and a simplified “flat” design of the HEMJ for long-pulse magnetic fusion energy (MFE) reactors. Experimental studies were performed in a helium (He) loop at the prototypical pressure of 10 MPa, nearly prototypical He temperatures and incident heat fluxes using test sections made from stainless steel and tungsten alloys. Correlations for average Nusselt numbers and pressure loss coefficients were developed from the data and are used to validate computational fluid dynamics (CFD) models. Here, this work presents updated thermal-fluids performance curves based on these correlations that estimate the maximum heat flux that can be accommodated by the plasma-facing surface and coolant pumping power requirements at prototypical operating conditions. Thermal-structural performance curves developed from ITER structural design criteria are introduced, which include protection against ductile and non-ductile failure, ratcheting fatigue, and creep fatigue. The performance design curves for these finger-type divertors demonstrate that the “flat” design, with a significantly less complicated geometry than the HEMJ, has thermal-fluid and thermal-structural performance comparable to the original HEMJ concept, and may be superior in terms of non-ductile failure criteria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Highest fusion performance without harmful edge energy bursts in tokamak

The path of tokamak fusion and International thermonuclear experimental reactor (ITER) is maintaining high-performance plasma to produce sufficient fusion power. This effort is hindered by the transient energy burst arising from the instabilities at the boundary of plasmas. Conventional 3D magnetic perturbations used to suppress these instabilities often degrade fusion performance and increase the risk of other instabilities. This study presents an innovative 3D field optimization approach that leverages machine learning and real-time adaptability to overcome these challenges. Implemented in the DIII-D and KSTAR tokamaks, this method has consistently achieved reactor-relevant core confinement and the highest fusion performance without triggering damaging bursts. This is enabled by advances in the physics understanding of self-organized transport in the plasma edge and machine learning techniques to optimize the 3D field spectrum. The success of automated, real-time adaptive control of such complex systems paves the way for maximizing fusion efficiency in ITER and beyond while minimizing damage to device components.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Novel dual-reflection design applied for ITER core x-ray spectrometer

Here. a novel dual-reflection configuration is introduced for the International Tokamak Experimental Reactor (ITER) core x-ray spectrometer to fit the allocated space where it will be placed accompanied by moving the detectors backward to reduce the incident radiation dose. The highly oriented pyrolytic graphite, which has a mosaic structure of microscopic crystallites, is chosen for the front reflector motivated by higher x-ray throughput and stronger misalignment tolerance compared to the perfect crystal reflector. In the ITER core x-ray spectrometer, a combination of several reflector-deflected Lines of Sight (LOSs) and a direct LOS is proposed for the first time named X-Ray Crystal Spectroscopy Core (XRCS-Core). The system is optimized to observe lines from externally seeded xenon and the intrinsic tungsten impurity, meeting both port integration needs and measurement requirements. Its spectral performance is simulated using an analytical-raytracing mixed code——XRSA, showing good imaging quality with a spectral resolution higher than 8000. The XRCS-Core system is thought to be applicable in various ITER scenarios through the assessment taking into account the spectrometers’ specifications and the chosen lines’ emissivity in different plasma parameters.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Models of resistive wall tearing mode disruptions

Disruptions are a serious issue in tokamaks. In a disruption, the thermal energy is lost by means of an instability which could be a resistive wall tearing mode (RWTM). During precursors to a disruption, the plasma edge region cools, causing the current to contract. Model sequences of contracted current equilibria are given, and their stability is calculated. A linear stability study shows that there is a maximum value of edge q a ≈ 3 for RWTMs to occur. Further, this also implies a minimum rational surface radius normalized to plasma radius from RWTMs to be unstable. Nonlinear simulations are performed using a similar model sequence derived from an equilibrium reconstruction. There is a striking difference in the results, depending on whether the wall is ideal or resistive. With an ideal wall, the perturbations saturate at moderate amplitude, causing a minor disruption without a thermal quench. With a resistive wall, there is a major disruption with a thermal quench, if the edge q a ≤ 3. There is a sharp transition in nonlinear behavior at q a = 3. This is consistent with the linear model and with experiments. If disruptions are caused by RWTMs, then devices with highly conducting walls, such as the International Tokamak Experimental Reactor will experience much milder, tolerable, disruptions than presently predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The motional Stark effect diagnostic for ITER

An overview of the plans for the motional Stark effect (MSE) diagnostic installation on the International Thermonuclear Experimental Reactor (ITER) is presented. The MSE diagnostic uniquely provides spatially localized magnetic field measurements inside the plasma. These are used to constrain equilibrium reconstructions to determine q(r), the safety factor as a function of minor radius. Meeting the system requirements to deliver q-profiles and related quantities with the specified radial resolution of 20 points over the minor radius, 10 ms time resolution, and better than 10% accuracy is challenging. MSE systems observe the D/H-α emission near 656.3 nm from neutral beams. As the beam atoms traverse the magnetic field, B⃗, at high velocity, v⃗, they experience a Lorentz electric field, v⃗×B⃗, which causes the spectral emission to be split and polarized due to the Stark effect. Traditional MSE-LP (line polarization) measurements determine the direction of the magnetic field in the observation volume using polarimetric analysis of the detected light. The harsh conditions of ITER are expected to deposit thin films of contaminants on the first mirror, which would alter the polarization state of reflected light significantly. On ITER, the combination of high magnetic field strength and high energy beams makes the Stark spectrum resolution suitable for the determination of the magnetic field magnitude from the line shift, so this approach has been selected. Every aspect of the measurement system must be planned for the burning plasma environment and carefully analyzed ahead of time. Current status and plans for the system are presented.

Instruments & Instrumentation↗
DOE OSTI2024

The slow wave resonance cone in the collisional regime

In the low-density edge plasma of tokamaks, ion cyclotron range of frequencies actuators may parasitically emit slow waves. If the density is sufficiently low, which may be common in large future devices such as international thermonuclear experimental reactor (ITER), these slow waves take the form of so-called resonance cones. The traditional theoretical description of this wave mode relies on formally relating an electrostatic approximation of the frequency-domain wave equation to a time-domain wave equation and relating the cone angle to the wave speed in the time-domain wave equation. In the cold plasma collisional regime, that wave speed is complex. We investigate that scenario in this work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Power balance analysis at the L-H transition in JET-ILW NBI-heated deuterium plasmas

The understanding of the physics underlying the L-H transition has strong implications for ITER experimental reactor and demonstration power plant (DEMO). In many tokamaks, including JET, it has been observed that, at a particular plasma density, n e,min , the power necessary to access H-mode P L-H is minimum. In the present work, L-H transitions of JET deuterium plasmas heated by neutral beam injection (NBI) are studied for the first time by means of a power balance analysis to characterize the main contributions in the transition, through integrated transport modelling. In the pulses analysed, we do observe a minimum of the L-H power threshold in density, indicating the presence of density branches and of n e,min . Electron and ion heat fluxes at the transition are estimated separately. The electron/ion equipartition power results in favour of the ions, as shown by QuaLiKiz quasilinear gyrokinetic simulations, which predict a larger ion transport that causes T e > T i . The resulting edge ion heat flux also shows a clear change of slope below n e,min , similarly to ASDEX-Upgrade (AUG) NBI pulses (Ryter et al 2014 Nucl. Fusion 54 083003). JET NBI data are compared to radio-frequency heated AUG and Alcator C-mod pulses (Schmidtmayr et al 2018 Nucl. Fusion 58 056003), showing a different trend of the power, coupled to ions at the L-H transition with respect to the linearity observed in the radio-frequency heated plasmas. Finally, the presence of n e,min and the role of the ion heat flux is discussed in the paper, although it seems it is not possible to explain the presence of a P L-H minimum in density by a critical ion heat flux and by the equipartition power for the JET NBI-heated plasmas analysed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Investigation of the compatibility of pellet fueling with ELM-free H-mode plasmas in EAST tokamak

Abstract Experiments on pellet fueling have been carried out in edge localized mode (ELM)-free high-confinement mode (H-mode) plasmas with q 95 ∼ 6 in the EAST tokamak. Cryogenic deuterium pellets were injected into the ELM-free plasmas at a frequency of 10 or 5 Hz from ∼45 cm above the mid-plane on the low-field side. It is found that the ELM-free H-mode plasmas are still sustained even if both the edge and core plasma are impacted by the pellet injections (PIs). Several small ELMs would appear and the edge coherent mode accompanying the ELM-free phase fades or even disappears just after the pellet events, but the plasma would rapidly recover to the ELM-free state. Although the ELMing phase is very short, it may be an issue that still needs to be resolved in the future. Meanwhile, the 3/2 tearing mode often appearing in ELM-free discharges would be stabilized by the PIs, and the high-Z impurity concentration would be reduced during the PIs. All these results will be meaningful to International Thermonuclear Experimental Reactor and future fusion reactors.

Physics↗
DOE OSTI2023

ORBIT simulations of fast ion power loads on the wall of the Divertor Tokamak Test

Abstract Neutral beam injection in tokamaks produces a population of fast ions, which interact with 3D magnetic fields in a variety of ways, often resulting in energetic particle losses in very short times to the wall. Careful design of neutral beams and active control of error fields helps to keep these losses to a minimum. Nevertheless, past experience in tokamaks in the 1980s and detailed simulations for future machines, such as International Thermonuclear Experimental Reactor, suggest that very localized fast ion losses (‘hot spots’) can be present, even if the overall losses are low. In this paper, we discuss this issue in the Divertor Tokamak Test (DTT) project, and we show that in the standard single-null full-power scenario of the DTT, fast ions produce two hot spots, corresponding approximately to the beam injection and exit toroidal angles: the former being mainly due to prompt losses/passing particles, while the latter is due to trapped ions. However, the maximum power load in these spots is of the order ∼ 63 kW m −2 , below the tolerance of plasma-facing components of the machine.

Physics↗
DOE OSTI2025

Ultra-low magnetization and hysteresis loss in APC Nb 3 Sn superconductors

For the accelerator magnets of the next hadron collider, reducing superconductor persistent-current magnetization is not only important for achieving the desired field quality, but also crucial for its sustainability because the magnetization loss is the major heat load to the magnet cold mass. For conventional Nb 3 Sn conductors this requires reduction of effective subelement size (D eff ). For the restacked-rod-process (RRP ® ) conductors a physical subelement size (D sub ) as small as 35 µm (corresponding to a D eff close to 45 µm) can be reached, but at a significant price in J c . Another way to reduce the magnetization is by introducing artificial pinning centers (APC) using the internal oxidation approach. APC conductors outperform conventional Nb 3 Sn wires in two aspects: 1) higher J c at high fields, and 2) much lower J c and magnetization at low fields (e.g. below 5 T). Here, in this work we explored the fabricability of APC wires with small D sub . A 180-stack APC wire was produced and drawn to 0.7- and 0.5 mm diameters with good quality, with D sub s of 34 and 24 µm (D eff s of 36 and 25 µm), respectively. For the 34 µm-D sub wire, its non-Cu J c is higher than that of an RRP ® wire used for the High-Luminosity Large Hadron Collider (HL-LHC) project above 13 T (e.g. 36% higher at 4.2 K, 18 T), while its non-Cu magnetization at 1 T, ΔM(1 T), is only 29% of the RRP ® wire. Its non-Cu hysteresis loss for a cycle between 1 and 14 T, Q h (1–14 T), is 37% of the RRP ® wire. For the 24 µm-D sub wire, its non-Cu J c surpasses the HL-LHC RRP ® wire above 17.5 T, while its ΔM(1 T) and Q h (1–14 T) are only 17% and 23% of the RRP® wire, respectively. Its non-Cu Q h (±3 T) even meets the specification of the International Thermonuclear Experimental Reactor project.

Xu, X [Fermi National Accelerator Laboratory (FNAL↗
DOE OSTI2023

Turbulence and E × B flow correlations across the L-H transition in DIII-D deuterium and hydrogen plasmas

The isotope dependence of the low- to high confinement-mode (L- to H-mode) transition power threshold P LH presents significant challenges for the initial (non-nuclear) hydrogen operations phase of the international thermonuclear experimental reactor (ITER). Here, we examine the isotope dependence of turbulence and E × B flow correlation properties in the L-mode edge plasma, leading up to the L-H transition. At marginal auxiliary power (near P LH ), turbulence is initially suppressed periodically during limit cycle oscillations (LCO) that precede the transition to sustained H-mode confinement. Here, we present evidence that the long-range (toroidal) correlation of the E × B edge plasma flow across the LCO phase is much weaker in hydrogen than in deuterium in the DIII-D tokamak, congruent with the higher threshold power P LH required to access LCO and H-mode in hydrogen. Concomitantly, the time required to initially quench edge turbulence via localized edge E × B flow shear is significantly longer in hydrogen (1–1.5 ms) than in deuterium (approx. 100 µs). No toroidal long-range correlation of the turbulence amplitude is observed, in agreement with expectations based on the relatively short poloidal turbulence correlation length. Radial edge turbulence and flow correlation lengths are longer in deuterium than in hydrogen plasmas as one would expect from ‘naïve’ gyro-Bohm isotope transport scaling, despite the substantially higher thermal flux across the last closed flux surface in hydrogen before the L-H transition. This article is part of a discussion meeting issue 'H-mode transition and pedestal studies in fusion plasmas'.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Turbulence-Driven Edge-Localized-Mode-Free High-Confinement Mode with Divertor Detachment in a Metal-Wall Tokamak

We report the first demonstration of a minute-scale, edge-localized-mode-free high-confinement plasma regime compatible with divertor partial detachment and enhanced pedestal performance in a metal-wall tokamak, the Experimental Advanced Superconducting Tokamak. This regime is enabled by a newly identified mechanism: during divertor partial detachment, reduced ionization and enhanced pumping in a closed divertor lead to less cooling of the pedestal by recycling neutrals and seeding impurities, resulting in an increased pedestal temperature gradient, which excites high-frequency broadband turbulence. Gyrokinetic simulations identify the high-frequency broadband turbulence as a temperature-gradient-driven trapped electron mode (𝜂 𝑒 -TEM), which drives outward transport of particles and heat, thereby maintaining the edge-localized-mode-free state. This pedestal regime is particularly promising for the International Thermonuclear Experimental Reactor, where the anticipated lower density gradient, reduced E×B shear, and lower collisionality in the pedestal will further facilitate 𝜂 𝑒 -TEM excitation. The achieved integrated scenario with a detached divertor and turbulence-dominated pedestal thus offers a compelling solution for managing heat loads and metal impurity sources for long-pulse high-performance operation in future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Simulation of Liquid Lithium Divertor Geometry Using SOLPS-ITER

Here, plasma-facing component (PFC) geometries are evaluated for a Fusion Nuclear Science Facility (FNSF) design to find solutions that are compatible with flowing liquid lithium components while satisfying requirements from the perspective of plasma and neutral particle transport. Flowing liquid metal (LM) divertor systems offer important advantages due to continual regeneration of the surface material, but may require modifications from standard optimized tokamak divertor designs. Scrape Off Layer Plasma Simulator for the International Thermonuclear Experimental Reactor (SOLPS-ITER) simulations are used to compare standard vertical target, open, and balanced baffled geometries. It is found that the open geometry offers little control of neutral particles and has narrow operating windows where the divertor fluxes can be lowered to acceptable levels, while maintaining the required upstream density. The addition of neutral baffles allows greater control over the divertor conditions with less impact on the upstream density. The divertor neutral pressure in the baffled geometry is slightly higher than the open divertor, but not as large with a vertical target. A first-pass coupling of the SOLPS solution with sheath and implantation/erosion models to compute the lithium emission shows that the lithium is largely confined to the near-surface region due to the small ionization mean free path and large main ion flow toward the surface. The lithium source was varied by three orders of magnitude. At the highest rate, when the sourced lithium approaches the level of the main ion recycling flux, a high core lithium concentration (~10%) is observed with only moderate power dissipation to affect the divertor plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Fast and Invertible Simplicial Approximation of Magnetic‐Following Interpolation for Visualizing Fusion Plasma Simulation Data

We introduce a fast and invertible approximation for fusion plasma simulation data represented as 2D planar meshes with connectivities approximating magnetic field lines along the toroidal dimension in deformed 3D toroidal spaces. Scientific variables (e.g., density and temperature) in these fusion data are interpolated following a complex magnetic-field-line-following scheme in the toroidal space represented by a cylindrical coordinate system. This deformation in the 3D space poses challenges for root-finding and interpolation. To this end, we propose a novel paradigm for visualizing and analyzing such data based on a newly developed algorithm for constructing a 3D simplicial mesh within the deformed 3D space. Our algorithm generates a tetrahedral mesh that connects the 2D meshes using tetrahedra while adhering to the constraints on node connectivities imposed by the magnetic field-line scheme. Specifically, we first divide the space into smaller partitions to reduce complexity based on the input geometries and constraints on connectivities. Then, we independently search for a feasible tetrahedralization of each partition, considering nonconvexity. We demonstrate our method with two X-Point Gyrokinetic Code (XGC) simulation datasets on the International Thermonuclear Experimental Reactor (ITER) and Wendelstein 7-X (W7-X), and use an ocean simulation dataset to substantiate broader applicability of our method. An open source implementation of our algorithm is available at https://github.com/rcrcarissa/DeformedSpaceTet.

Ren, Congrong [The Ohio State Univ., Columbus, OH ↗
DOE OSTI2024

Advanced Closures for Nonlinear Fluid Simulations of Plasmas (Final Report)

Advanced Closures for Nonlinear Fluid Simulations of Plasmas E. Held, Utah State University (Principal Investigator) J.-Y. Ji, Utah State University (Co-Investigator) Utah State University (USU) continued to develop and apply closures for hybrid fluid/kinetic simulations of magnetized fusion plasmas. The closures developed by USU's fusion theory and computation group use the drift ordering and are referred to as "continuum" and "general moment" approaches. They differ in their representation for the distribution functions which probabilistically describe the density of particles in 5 dimensional phase space. The two chosen velocity variables, needed in studies of kinetic physics for magnetized plasmas, are the speed and the component of a particle's velocity along the magnetic field. In magnetized plasmas, the charged particle trajectories are tied closely to the magnetic field lines. Both approaches account for binary interactions between particles using the Coulomb collision operator. Theoretical and computational development along both lines provided avenues for verification between the two approaches. This funding allowed USU to further develop physically accurate and numerically efficient models that advance our understanding of kinetic physics in experiments like the International Thermonuclear Experimental Reactor (ITER). The ITER device will likely prove the physical practicality of thermonuclear fusion by producing 500 megawatts of clean, abundant fusion power. This work helped to realize this goal by advancing our simulation capability using hybrid fluid/kinetic models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗