Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium reconstruction”

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

Investigation of core transport changes in DIII-D discharges with off-axis T e profile peaks

Abstract DIII-D discharges that transition to H-mode solely with off-axis electron cyclotron heating (ECH) often exhibit strong off-axis peaking of electron temperature profiles at the heating location. Electron heat transport properties near these off-axis temperature peaks have been studied using modulated ECH. The Fourier analyzed electron temperature data have been used to infer electron thermal diffusivity. Comparisons with numerical solutions of the time-dependent electron thermal equation find that the data are consistent with a narrow region with electron diffusivity χ e an order of magnitude lower than the average value across the plasma, suggesting an electron internal transport barrier (ITB) near the ECH heating location. Detailed profile analysis and equilibrium reconstructions suggest that the formation of these ITBs are correlated with off-axis values of the safety factor q being near 1. Furthermore, the ECH driven H-mode discharges demonstrate more rapid electron heating rate near the ECH deposition location than L-mode discharges with higher auxiliary ECH heating power. Additional modeling attributes this difference to the modification of electron heat transport in the core at the L-H transition, which also sustains the off-axis electron temperature peaks.

Physics↗
DOE OSTI2025

A review of collaborative studies between the NSTX/-U and MAST/-U spherical tokamaks

The National Spherical Torus Experiment (NSTX) at the Princeton Plasma Physics Laboratory in the United States, and the mega ampere spherical tokamak (MAST) at the United Kingdom Atomic Energy Authority in the United Kingdom, and their respective upgrades (NSTX-U and MAST-U) are two MAST fusion devices that have operated roughly over the past two decades. Both devices have made significant contributions to understanding spherical tokamak (ST) plasma physics, and fusion plasmas in general, and both have contributed data to multi-machine database studies. Several diagnostics have been physically moved from one machine to the other by diagnostic teams working on both devices. Collaboration has benefited both research teams in the areas of operational expertise, scenario development, and equilibrium reconstruction techniques. More focused comparative studies between the two devices have been pursued over the years in many areas as well, including stability calculations, disruption characterization, pedestal and edge localized mode stability, confinement and transport, energetic particles, and heating and current drive modelling. Together NSTX/-U and MAST/-U set the stage for the future of STs, which is entering the phase of design of demonstration power plant devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modification of the magneto-hydrodynamic equilibrium by the lower-hybrid wave driven fast electrons on the TST-2 spherical tokamak

Presence of a large population of fast particles may qualitatively modify the tokamak equilibrium from the one that can be described by standard magneto-hydrodynamics (MHD). The kinetic modification of the MHD equilibrium was studied for a lower-hybrid (LH) wave driven plasma on the TST-2 spherical tokamak. The analysis was performed using an equilibrium reconstruction method based on an extended MHD model that considers a two-component plasma of bulk MHD fluid and kinetic fast electrons. Scrape-off-layer current appeared naturally in the extended MHD analysis because of the finite electron orbit excursion from the flux surfaces. This reduced the bulk pressure contribution to the toroidal current substantially. Furthermore, the resulting poloidal flux profile was more consistent with that inferred from the Thomson scattering measurement.

extended magneto-hydrodynamics↗
DOE OSTI2023

MHD stability of negative triangularity DIII-D plasmas

Abstract Negative triangularity (NT) experiments in DIII-D point to an emergent reactor scenario free of sawteeth, endowed with benign, nondisruptive n = 2 tearing modes, which experience q min ⩾ 1 similar to the positive triangularity hybrid scenario. Plasmas exhibiting this behavior attain β N > 3 , high enough to reconsider long held views of NT stability. Ideal MHD and tearing stability analysis of well-diagnosed equilibrium reconstructions of experimental hybrid-like plasmas predict that among shape parameters, MHD stability limits are only sensitive to average triangularity. Operation is predicted to be possible at β N relevant to upcoming tokamaks and commercial NT reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Development of high-performance long-pulse discharge in KSTAR

High-performance long-pulse plasma operation is essential for producing economically viable fusion energy in tokamak devices. To achieve such discharges in KSTAR, firstly, the rapid increase in the temperature of plasma-facing components was mitigated. The temperature increase of the poloidal limiter, especially, was associated with beam-driven fast ion orbit loss and the discrepancy of the equilibrium reconstructed with heated magnetic probes of signal drift. The fast ions lost to the poloidal limiter were reduced by optimizing the plasma shape and the composition of neutral beam injection (NBI). This nonlinear signal drift was successfully reduced by a new thermal shielding protector on the magnetic probes. Secondly, a lower loop voltage approach was implemented to reduce a poloidal flux consumption rate. A plasma current of 400 kA and a line-averaged electron density of ~2.0 × 10 19 m –3 were chosen by considering the L–H power threshold, fast ion orbit loss, and beam shine-through power loss for low loop voltage in KSTAR. In addition, the application of electron cyclotron heating also helped maintain the plasma with low loop voltage (~25 mV) by enhancing the NBI-driven current and achieving a high poloidal beta (β P ) state. KSTAR has achieved a long pulse (~90 s) operation with the high performance of β P $\leqslant$ 2.7, thermal energy confinement enhancement factor (H 98y2 ) ~ 1.1, and fraction of non-inductive current (f NI ) $\leqslant$ 0.96. Still, gradual degradation of the plasma performance has been observed over time in the discharges. In one of the long-pulse discharges, β P reduced by ~18% over the time of ~8τ R (current relaxation time, τ R ~5 s) and ~1067τ E,th (thermal energy confinement time, τ E,th ~45 ms). The degradation may be closely associated with weak, yet growing, and persistent toroidal Alfvén eigenmodes and their effect on fast ion confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

NSTX-U research advancing the physics of spherical tokamaks

Abstract The objectives of NSTX-U research are to reinforce the advantages of STs while addressing the challenges. To extend confinement physics of low- A , high beta plasmas to lower collisionality levels, understanding of the transport mechanisms that set confinement performance and pedestal profiles is being advanced through gyrokinetic simulations, reduced model development, and comparison to NSTX experiment, as well as improved simulation of RF heating. To develop stable non-inductive scenarios needed for steady-state operation, various performance-limiting modes of instability were studied, including MHD, tearing modes, and energetic particle instabilities. Predictive tools were developed, covering disruptions, runaway electrons, equilibrium reconstruction, and control tools. To develop power and particle handling techniques to optimize plasma exhaust in high performance scenarios, innovative lithium-based solutions are being developed to handle the very high heat flux levels that the increased heating power and compact geometry of NSTX-U will produce, and will be seen in future STs. Predictive capabilities accounting for plasma phenomena, like edge harmonic oscillations, ELMs, and blobs, are being tested and improved. In these ways, NSTX-U researchers are advancing the physics understanding of ST plasmas to maximize the benefit that will be gained from further NSTX-U experiments and to increase confidence in projections to future devices.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Investigation of pedestal parameters and divertor heat fluxes in small ELM regimes in DIII-D

Abstract Divertor heat flux and its correlation with pedestal parameters within various small edge localized mode (ELM) regimes, including high beta poloidal, type-II and ELMs with negative triangularity H-modes were investigated in DIII-D. The parallel energy fluences of type-II and high beta poloidal small ELM regimes fall below the linear scaling with pedestal electron pressure for type-I ELMs put forward in Eich et al 2017 ( Nucl. Mater. Energy 12 84–90). The negative triangularity of H-mode ELMs follow the Eich scaling for type-I ELMs. The parallel heat flux and total heat loads to the divertor were determined using high-time resolution infrared thermography, while pedestal parameters were obtained through self-consistent kinetic equilibrium reconstructions. Linear regressions for the type-II and high beta poloidal regimes demonstrate that an equivalent 7.5 MA small ELM scenario in ITER would fall below the ~5 MJ m − 2 leading edge melting limit for tungsten (Gunn et al 2017 Nucl. Fusion 57 046025). Utilizing fast thermography, the scrape-off layer power fall-off length for both inter-ELM and intra-ELM was determined and compared to the Eich scaling with poloidal magnetic field in Eich et al (ASDEX Upgrade Team and JET EFDA Contributors 2013 Nucl. Fusion 53 093031). Except for the high beta poloidal scenario, all the small ELM regimes during both inter- and intra-ELM periods had power fall-off lengths ( λ q ) larger then would be expected from the B pol , MP − 1 scaling associated with type-I ELMs, signifying their potential in managing heat loads and offering a solution for core–edge integration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Metrics and extrapolation of resonant magnetic perturbation thresholds for ELM suppression

This large database study of resonant magnetic perturbation (RMP) edge localized mode (ELM) suppression thresholds in the AUG, DIII-D, EAST, and KSTAR tokamaks details the key strengths and weaknesses of RMP metrics. The RMP ELM suppression database used for this work contains plasma information at the time of transition from ELMing to ELM suppressed states where a clear experimental threshold is identified. The experimental threshold distributions are compared for five metrics: (1) the island overlap width, (2) pedestal top Chirikov overlap, (3) peeling edge displacement, (4) pedestal top resonant drive, and (5) edge dominant mode overlap. The distributions, the regularity of the dependence on RMP coil currents, and the sensitivities of a given metric to equilibrium reconstruction details are compared. The overlap metric proves to be a good compromise between including the appropriate plasma response physics and maintaining a numerical robustness. This quantity does not exhibit clear power-law scalings for projection, but machine learning can assist in predicting thresholds within the existing parameter ranges and providing uncertainty quantification of those predictions. Two new first-principles models, one utilizing a threshold from the non-linear Modified Rutherford equation evaluated at the pedestal top and one utilizing the SLAYER code to calculate the linear tearing threshold from torque balance, offer possible paths to extrapolation beyond the existing database parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Conceptual study on using Doppler backscattering to measure magnetic pitch angle in tokamak plasmas

We introduce a new approach to measure the magnetic pitch angle profile in tokamak plasmas with Doppler backscattering (DBS), a technique traditionally used for measuring flows and density fluctuations. The DBS signal is maximised when its probe beam's wavevector is perpendicular to the magnetic field at the cutoff location, independent of the density fluctuations [Hillesheim \emph{et al} 2022 \emph{Nucl. Fusion} \textbf{55} 073024]. Hence, if one could isolate this effect, DBS would then yield information about the magnetic pitch angle. By varying the toroidal launch angle, the DBS beam reaches cutoff with different angles with respect to the magnetic field, but with other properties remaining similar. Hence, the toroidal launch angle which gives maximum backscattered power is thus that which is matched to the pitch angle at the cutoff location, enabling inference of the magnetic pitch angle. We performed systematic scans of the DBS toroidal launch angle for repeated DIII-D tokamak discharges. Experimental DBS data from this scan were analysed and combined with Gaussian beam-tracing simulations using the Scotty code [Hall-Chen \emph{et al} 2022 \emph{Plasma Phys. Control. Fusion} \textbf{64} 095002]. The pitch-angle inferred from DBS is consistent with that from magnetics-only and motional-Stark-effect-constrained (MSE) equilibrium reconstruction in the edge. In the core, the pitch angles from DBS and magnetics-only reconstructions differ by one to two degrees, while simultaneous MSE measurements were not available. The uncertainty in these measurements was under a degree; we show that this uncertainty is primarily due to the error in toroidal steering, the number of toroidally separated measurements, and shot-to-shot repeatability. We find that the error of pitch-angle measurements can be reduced by optimising the poloidal launch angle and initial beam properties. Since DBS has high spatial and temporal resolutions, is non-perturbative, does not require neutral beams, and is likely robust to neutron damage of and debris on the first mirrors, using DBS to measure the pitch angle in future fusion energy systems is especially appealing.

beam tracing↗
DOE OSTI2026

Variability of MHD instabilities in benign termination of high-current runaway electron beams in the JET and DIII-D tokamaks

Benign termination, in which magnetohydrodynamic (MHD) instabilities deconfine runaway electrons (REs) following hydrogenic injections, is a promising strategy for mitigating dangerous RE loads after disruptions. Recent experiments on the Joint European Torus (JET) have explored this scenario at higher pre-disruptive plasma currents than are achievable on other devices, revealing challenges in obtaining benign terminations at I p ≥2.5 MA. This work analyzes the evolution of these high-current RE beams and their terminating MHD events using fast magnetic sensor measurements and EFIT equilibrium reconstructions for approximately 40 JET and 20 DIII-D tokamak discharges. On JET, unsuccessful non-benign terminations occur at low edge safety factor (q edge ≈ 2), and are preceded by intermittent, non-terminating MHD events at higher rational qedge. Trends in the internal inductance I i indicate more peaked RE current profiles in the high-I p non-benign population, which may hinder successful recombination through re-ionization of the companion plasma. In contrast, benign terminations on JET typically occur at higher q edge ≥3 and exhibit less peaked RE current profiles. DIII-D displays a broader range of terminating edge safety factors, again correlated with the measured values. Across both tokamaks, the RE current peaking is therefore found to determine which MHD instability boundary is encountered, a result confirmed by linear resistive MHD modeling with the CASTOR3D code. Measured growth rates are similar for benign and non-benign cases, indicating that ideal MHD timescales at low density after hydrogenic injection do not alone explain efficient RE deconfinement. Instead, non-benign cases are most readily characterized by their comparably lower overall MHD perturbation amplitudes δB. These observations suggest that the interplay between ideal and resistive dynamics governs the termination process, with implications for extrapolating benign RE termination to high-I p reactor scenarios.

MHD instabilities↗
DOE OSTI2020

STELLOPT

The STELLOPT code is designed to optimize 3D MHD equilibria to a set of target physics parameters encompassing stellarator design and 3D equilibrium reconstructions.

Lazerson, Samuel↗
DOE OSTI2025

Real-time Electron Temperature and Density Profile Measurements for NSTX-U

This project set out to design, build, and deploy a real-time multi-point Thomson scattering (rt MPTS) system for NSTX-U, capable of providing electron temperature and density profiles—and their uncertainties—in real time for use in equilibrium reconstruction and plasma control. The long-term goal was to move Thomson scattering from a purely post-shot diagnostic to a central element of real-time situation awareness and feedback control on NSTX-U.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Accelerating discoveries at DIII-D with the Integrated Research Infrastructure

DIII-D research is being accelerated by leveraging high performance computing (HPC) and data resources available through the National Energy Research Scientific Computing Center (NERSC) Superfacility initiative. As part of this initiative, a high-resolution, fully automated, whole discharge kinetic equilibrium reconstruction workflow was developed that runs at the NERSC for most DIII-D shots in under 20 min. This has eliminated a long-standing research barrier and opened the door to more sophisticated analyses, including plasma transport and stability. These capabilities would benefit from being automated and executed within the larger Department of Energy Advanced Scientific Computing Research program’s Integrated Research Infrastructure (IRI) framework. The goal of IRI is to empower researchers to meld DOE’s world-class research tools, infrastructure, and user facilities seamlessly and securely in novel ways to radically accelerate discovery and innovation. For transport, we are looking at producing flux matched profiles and also using particle tracing to predict fast ion heat deposition from neutral beam injection before a shot takes place. Our starting point for evaluating plasma stability focuses on the pedestal limits that must be navigated to achieve better confinement. This information is meant to help operators run more effective experiments, so it needs to be available rapidly inside the DIII-D control room. So far this has been achieved by ensuring the data is available with existing tools, but as more novel results are produced new visualization tools must be developed. In addition, all of the high-quality data we have generated has been collected into databases that can unlock even deeper insights. This has already been leveraged for model and code validation studies as well as for developing AI/ML surrogates. The workflows developed for this project are intended to serve as prototypes that can be replicated on other experiments and can be run to provide timely and essential information for ITER, as well as next stage fusion power plants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Reconstruction of plasma equilibrium and eddy current induced error fields in the Lithium Tokamak eXperiment-β (Final Report)

This is the final report for the DOE award DE-SC0019239, titled “Reconstruction of plasma equilibrium and eddy current induced error fields in the Lithium Tokamak eXperiment-b” for the period July 1, 2018 – July 2, 2023. This award supported research into computational models for plasma equilibrium and 3D conducting structures in the Lithium Tokamak eXperiment-b (LTX-b) using the PSI-Tri/TokaMaker and PSI-Tet codes. Additional detail on progress in support of these tasks is provided below.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Plasma flux expansion control on the DIII-D tokamak

A new controller has been developed with help of the flexible divertor poloidal-field coil set of the DIII-D tokamak, to aid in the precise control of the flux expansion in the scrape-off layer (SOL). Here, the single-input multiple-output (SIMO) architecture ensures flexibility through a complementary set of orthogonal actuator direction to guarantee minimum effect on existing controlled variables, e.g. radial and vertical position control of the X-point. A non-linear free-boundary simulation code (GSevolve) is used for simulating the closed-loop response and for verifying the implementation of the control algorithm on the DIII-D plasma control system (PCS). First results of the experimental commissioning of the new controller during 2020 DIII-D campaign are also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Single Gaussian process method for arbitrary tokamak regimes with a statistical analysis

Abstract Gaussian process regression is a Bayesian method for inferring profiles based on input data. The technique is increasing in popularity in the fusion community due to its many advantages over traditional fitting techniques including intrinsic uncertainty quantification and robustness to over-fitting. This work investigates the use of a new method, the change-point method, for handling the varying length scales found in different tokamak regimes. The use of the Student’s t-distribution for the Bayesian likelihood probability is also investigated and shown to be advantageous in providing good fits in profiles with many outliers. To compare different methods, synthetic data generated from analytic profiles is used to create a database enabling a quantitative statistical comparison of which methods perform the best. Using a full Bayesian approach with the change-point method, Matérn kernel for the prior probability, and Student’s t-distribution for the likelihood is shown to give the best results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Late Holocene glacier and climate fluctuations in the Mackenzie and Selwyn mountain ranges, northwestern Canada

Over the last century, northwestern Canada experienced some of the highest rates of tropospheric warming globally, which caused glaciers in the region to rapidly retreat. Our study seeks to extend the record of glacier fluctuations and assess climate drivers prior to the instrumental record in the Mackenzie and Selwyn mountains of northwestern Canada. We collected 27 10 Be surface exposure ages across nine cirque and valley glacier moraines to constrain the timing of their emplacement. Cirque and valley glaciers in this region reached their greatest Holocene extents in the latter half of the Little Ice Age (1600–1850 CE). Four erratic boulders, 10–250 m distal from late Holocene moraines, yielded 10 Be exposure ages of 10.9–11.6 ka, demonstrating that by ca. 11 ka, alpine glaciers were no more extensive than during the last several hundred years. Estimated temperature change obtained through reconstruction of equilibrium line altitudes shows that since ca. 1850 CE, mean annual temperatures have risen 0.2–2.3 °C. We use our glacier chronology and the Open Global Glacier Model (OGGM) to estimate that from 1000 CE, glaciers in this region reached a maximum total volume of 34–38 km 3 between 1765 and 1855 CE and had lost nearly half their ice volume by 2019 CE. OGGM was unable to produce modeled glacier lengths that match the timing or magnitude of the maximum glacier extent indicated by the 10 Be chronology. However, when applied to the entire Mackenzie and Selwyn mountain region, past millennium OGGM simulations using the Max Planck Institute Earth System Model (MPI-ESM) and the Community Climate System Model 4 (CCSM4) yield late Holocene glacier volume change temporally consistent with our moraine and remote sensing record, while the Meteorological Research Institute Earth System Model 2 (MRI-ESM2) and the Model for Interdisciplinary Research on Climate (MIROC) fail to produce modeled glacier change consistent with our glacier chronology. Finally, OGGM forced by future climate projections under varying greenhouse gas emission scenarios predicts 85 % to over 97 % glacier volume loss by the end of the 21st century. The loss of glaciers from this region will have profound impacts on local ecosystems and communities that rely on meltwater from glacierized catchments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Ringing Artifact Reduction Method for Ultrasound Reconstruction Using Multi-Agent Consensus Equilibrium

Non-destructive characterization of multi-layered structures that can be accessed from only a single side is important for applications such as well-bore integrity inspection. Existing methods related to Synthetic Aperture Focusing Technique (SAFT) rapidly produce acceptable results but with significant artifacts. Recently, ultrasound model-based iterative reconstruction (UMBIR) approaches have shown significant improvements over SAFT. However, even these methods produce ringing artifacts due to the high fractional-bandwidth of the excitation signal.In this paper, we propose a ringing artifact reduction method for ultrasound image reconstruction that uses a multi-agent consensus equilibrium (RARE-MACE) framework. Our approach integrates a physics-based forward model that accounts for the propagation of a collimated ultrasonic beam in multi-layered media, a spatially varying image prior, and a denoiser designed to suppress the ringing artifacts that are characteristic of reconstructions from high-fractional bandwidth ultrasound sensor data. We test our method on simulated and experimental measurements and show substantial improvements in image quality compared to SAFT and UMBIR.

Alanazi, Abdulrahman M↗