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

Analysis of neoclassical tearing mode stabilization experiment by electron cyclotron injection in KSTAR

We report the neoclassical tearing mode (NTM) stabilization experimental results by injection of the electron cyclotron (EC) beam in KSTAR, and its analysis with integrated modelling. In the KSTAR experiment, NTM was intentionally induced and the EC beam was injected by the plasma control system (PCS) to stabilize it. Here, the EC angle was controlled based on the minimum seeking algorithm for the island growth rate. Although the EC angle was not changed significantly, the island width significantly decreased towards the end of the experiment followed by the recovery of plasma performance. To assess the NTM stabilization experiment from a stability perspective, we performed integrated modelling with TRIASSIC utilizing the reconstructed equilibrium with magnetic diagnostics and motional stark effect diagnostics. As a result, we observed that the NTM is stabilized by EC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Evaluation of DIII-D plasmas for Future Measurements of the Driven Current Density in the DIII-D High Field Side Lower Hybrid Experiment

High field side (HFS) lower hybrid current drive (LHCD) is a promising method for efficiently driving off-axis current in steady state tokamak power plants. The first test of this technology is underway at the DIII-D tokamak. The initial physics goal is to measure the LH-driven current density profile in order to validate the ray-tracing/Fokker-Plank codes GENRAY/CQL3D and gain confidence in their ability to predict LH deposition in future DIII-D experiments and fusion power plants. DIII-D’s motional Stark effect (MSE) diagnostic will be used to constrain the magnetic equilibrium reconstructions. From these reconstructions, the ohmic, bootstrap, and driven non-inductive components of the current density profile can be extracted. Several recent DIII-D plasmas have been identified in which GENRAY/CQL3D predict measurable amounts of current with only 100 kW of coupled power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The magnetic structure of long-wavelength magnetohydrodynamic modes in current-carrying stellarator plasmas

Magnetic field fluctuations are observed in current-carrying stellarator plasmas when the rotational transform is close to a rational value at the edge of the plasma. At low plasma pressure, these fluctuations are associated with perturbed currents parallel to the equilibrium magnetic field lines. A model for these magnetohydrodynamic modes in a low-β, three-dimensional stellarator equilibria has been developed. A set of helical current filaments are constrained to mimic the structure of magnetic field lines on rational surfaces derived from three-dimensional (3D) equilibrium reconstructions. Transformation to straight field line coordinates then allows fitting of the poloidal magnetic sensor data to a single harmonic function, which fixes the modeled toroidal mode structure via the field line flow geometry. The developed procedure accurately captures phase and amplitude variation for m/n = 3/2, 3/1, and 4/1 modes in the 3D equilibria of the compact toroidal hybrid experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

On the frequency bifurcations of the MHD startup modes in NSTX

The observed bifurcations of the low frequency (<50 kHz) and low toroidal periodicity (n < 5) magnetohydrodynamic (MHD) activity often present in the initial part of the National Spherical Tokamak Experiment (NSTX) discharges can be explained by the evolution of the radial profile of the safety factor (q=rB φ /RB θ ) crossing multiple rational surfaces in the core. Important performance limiting instability mechanisms in the NSTX spherical tokamak are often linked to low frequency and low-n MHD activity. They are quite common in long-pulse NSTX plasmas. They can be present at the beginning of the plasma current flat-top, at the end of the discharge or during the whole duration, and they have been observed to deleteriously impact performance over a wide range of q 95 . An interesting feature observed in some NSTX discharges is the presence of a bifurcation in the frequency of the low n modes, as low as n = 1, that have frequencies comparable to the plasma core rotation divided by n. Equilibrium reconstructions constrained by magnetic diagnostics data and motional stark effect pitch angle radial profiles suggest that the observed bifurcations are linked to a fast evolving minimum value of q. Finally, 3D non-linear resistive MHD simulations show that these modes are ideal and exist as non-resonant before the correspondent rational surface enters the plasma.

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

Improved accuracy and robustness of electron density profiles from JET’s X-mode frequency-modulated continuous-wave reflectometers

JET’s frequency-modulated continuous wave (FMCW) reflectometers have been operating well with the current design since 2005, and density profiles have been automatically calculated intershot since then. However, the calculated profiles had long suffered from several shortcomings: poor agreement with other diagnostics, sometimes inappropriately moving radially by several centimeters, elevated levels of radial jitter, and persistent wriggles (strong unphysical oscillations). Here, in this research, several techniques are applied to the reflectometry data analysis, and the shortcomings are significantly improved. Starting with improving the equilibrium reconstruction that estimates the background magnetic field, adding a ripple correction in the reconstructed magnetic field profile, and adding new inner-wall reflection positions estimated through ray-tracing, these changes not only improve the agreement of reconstructed profiles to other diagnostics but also solve density profile wriggles that were present during band transitions. Other smaller but also persistent wriggles were also suppressed by applying a localized correction to the measured beat frequency where persistent oscillations are present. Finally, the burst analysis method, as introduced by Varela et al., has been implemented to extract the beat frequency from stacked spectrograms. Due to the strong suppression of spurious reflections, the radial jitter that sometimes would span several centimeters has been strongly reduced. The stacking of spectrograms has also been shown to be very useful for stacking recurring events, like small gas puff modulations, and extracting transport coefficients that would otherwise be below the noise level.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Aligning the Thomson scattering and charge exchange recombination diagnostics using neutral beam emission at DIII-D

This work addresses discrepancies in the alignment of the H-mode pedestal profiles of the electron and ion properties in the DIII-D tokamak as measured by Thomson Scattering (TS) and Charge Exchange Recombination Spectroscopy (CER) diagnostics. While the alignment of these profiles is key for accurate studies of tokamak physics and plasma confinement, misalignments can occur due to inaccuracies, such as in magnetic equilibrium reconstructions required to map measurements in different poloidal and toroidal locations. Both FIDASIM, an established simulation package, and a simplified collisional radiative model are used to simulate neutral beam state densities and neutral beam emission. Simulated neutral beam emissions are calculated based on shifted TS profiles and compared to beam emission measurements from the Main Ion CER system to determine the best shift for aligning TS with CER. This analysis is performed on various DIII-D discharges.

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

A magnetic diagnostic suite for the Pegasus-III experiment

Pegasus-III is an ultralow aspect ratio spherical tokamak providing a dedicated US experiment for comparative solenoid-free startup studies. A new magnetic diagnostic suite for equilibrium and low frequency (<200 kHz) magnetohydrodynamic mode analysis has been installed. These new diagnostics address the significant challenges of measuring magnetic field in a high noise environment with the majority constrained to fit in an 8 mm diagnostic gap on the high field side. Electrostatic switching noise generated by the 16 independent current feedback-controlled power supplies produces dV cm /dt ~ 1 kV/μs and volt level common mode noise on the magnetics. Immunity to this switching noise is accomplished through differential signal runs and signal processing, along with end-to-end electromagnetic interference shielding. The magnetic measurements are simultaneously digitized at 1 MHz and conditioned by precision 8 pole Butterworth filters with a corner frequency of 200 kHz to prevent aliasing down to the 16-bit level over the full passband. Ex-vessel calibrations of the B p coils were completed with a typical uncertainty of <0.5%. Stray toroidal field pickup from coil misalignment or positioning errors is corrected using a physics-based model. Comparisons of the corrected measurements to modeling agree to within 1.3% on average. This is within the 1.5% measurement uncertainty that a sensitivity analysis determined is needed for accurate fast boundary and equilibrium reconstruction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Forming a database to study reversed magnetic shear from the National Spherical Torus eXperiment using machine learning

Achieving a long-lived reversed magnetic shear (RMS) target plasma in the National Spherical Torus eXperiment Upgrade will require developing various sustainment scenarios. To help with the ongoing plasma control efforts, the development of a new analysis for the motional Stark effect (MSE) diagnostic using a machine learning algorithm, namely, MSE-ML, is described. MSE-ML will be used to identify patterns during RMS discharges, some of which suffer magnetohydrodynamic (MHD) events resulting in current redistribution and monotonic q-profiles. A database consisting of q and magnetic shear profiles is being constructed primarily based on the existing National Spherical Torus eXperiment data with equilibrium reconstructions constrained by the magnetic field pitch angle profile measured using the multi-channel MSE diagnostic. An unsupervised k-means clustering of the data is developed to study the RMS formation as a function of time. The initial clustering from the q-profiles shows significant differences in both amplitude and the duration of the RMS period. As a goal, the clustering results that detect and distinguish shots with substantial and sustained RMS are to be used as a preprocessing step in a supervised algorithm to identify the underlying conditions that lead to long-lasting improved confinement with RMS. Another aim of the MSE-ML study is to identify precursors of RMS-destroying MHD events in either derived data such as the q-profile or directly measured data such as the magnetic field pitch angle profile.

Uzun-Kaymak, I. U. (ORCID:0000000276251493)↗
DOE OSTI2025

Integrated modeling of advanced tokamaks for tearing mode avoidance on DIII-D

Understanding the physics of the onset and evolution of tearing modes (TMs) in tokamak plasmas is crucial for high-performance steady-state operations. The onset condition of (m,n) = (2,1) tearing stability, $Δ′$ > $Δ_c^′$> 0, has been studied in DIII-D steady-state hybrid discharges with accurate equilibrium reconstruction and well-measured plasma profiles. Here, m is the poloidal mode number and n is the toroidal mode number. The onset of n = 1 TMs is carefully identified using the Mirnov array to determine in-phase signals for the onset time and using electron cyclotron emission measurements to trace back from the saturated island width in the experiments for both initial island width and its growth rate. The tearing stability index $Δ′$ at TM onset is calculated in multiple ways, using the modified Rutherford equation and using the PEST-III or resistive DCON magnetohydrodynamic codes, and the results are in reasonable agreement with each other. The calculated $Δ′$ is compared with an analytical formula of the tearing stability threshold $Δ_c^′$ to determine the condition for the mode onset, $Δ′$ > $Δ_c^′$. Finally, the local gradient effect of the plasma current and the pressure profiles on the tearing stability is investigated in terms of the difference between $Δ′$ and $Δ_c^′$ through a series of modeled equilibria with the parametric variations. A database of experimental mode instability onsets is compatible with the simulated regions of instability obtained from the parametric variations.

Linear stability analysis↗
DOE OSTI2025

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗
DOE OSTI2022

A New Integrated Analysis Suite for Fast-Ion Study in KSTAR

Here, an integrated workflow for fast-ion analysis was developed by adapting the One Modeling Framework for Integrated Task (OMFIT) workflow manager to support a standard and unified analysis platform for KSTAR users. The newly established analysis suite offers a graphical user interface–based workflow to enable users to readily access and handle experimental data archived in various data formats and servers. Further, users can analyze the data by importing modules designed for conducting certain tasks, such as profile fitting, equilibrium reconstruction, and postprocessing of tokamak data. The procedures for preparing the inputs for fast-ion simulations are streamlined by a common workflow manager, which enables the parallel processing of various tasks to efficiently analyze large fast-ion datasets. The OMFIT platform comprises a flexible Python-based application that enables users to freely manipulate the Python scripts for applications that are unavailable in the standard workflow. The framework also offers mapping tools to translate the output data into the Integrated Modeling and Analysis Suite format to maintain application compatibility for future ITER burning plasma experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Doppler Backscattering Data Analysis and Integrated Modeling with OMFIT

One Modeling Framework for Integrated Tasks (OMFIT) is a widely used software tool in the magnetic fusion research community. OMFIT provides magnetic fusion energy researchers with a framework for the development of special-purpose physics modules. This paper describes an OMFIT physics module pertaining to the Doppler Backscattering (DBS) fusion plasma diagnostic. DBS measures density fluctuations and flow velocity through plasma scattering of electromagnetic waves. The OMFIT DBS module was developed to analyze experimental DBS data and facilitate modeling of DBS systems installed on multiple tokamak devices. The OMFIT DBS module is designed to support several analysis workflows: detailed analysis of experimental data, experimental planning, and theory-based synthetic diagnostic modeling. The DBS module uses integrated modeling by leveraging other OMFIT physics modules to perform tasks related to DBS, e.g. ray/beam–tracing simulations, edge-localized mode–synchronized data analysis, magnetic equilibrium reconstruction, and fitting kinetic profile data. Furthermore, this paper describes several supported workflows and serves a reference for the OMFIT DBS module.

Doppler backscattering↗
DOE OSTI2022

Fast calculation of the tokamak vertical instability

There has been recent interest in fast calculations of the tokamak axisymmetric vertical instability for real time feedback control purposes. It is shown that the maximum eigenvalue for the basic rigid version of this stability problem can be obtained by finding the positive root to a simple scalar function. This function can be generalized to include plasma mass and has complexity linear in the number of conductive elements. The formulation is based on standard matrix decompositions of the fixed-geometry part of the eigenproblem. The calculation bottleneck is the summary of mutual inductances from the reconstructed equilibrium current density. The with-mass spectrum can be made fully real-valued by the addition of a critical amount of damping with negligible effect on the vertical growth rate. Furthermore, the calculation has been implemented in the plasma control system at the DIII-D tokamak and used in experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

ECEI characterization of pedestal fluctuations in quiescent H-mode plasmas in DIII-D

Here, electron cyclotron emission imaging (ECEI) is employed to characterize the magneto hydraulics dynamics (MHD) fluctuations at the quiescent H-mode pedestals in DIII-D. Pedestal MHD fluctuations cause ECE radiation temperature fluctuations $\delta {T_{e{\text{,rad}}}}$ in both the pedestal and scrape-off-layer (SOL). A synthetic ECE platform is utilized for detailed interpretation of the ECEI signals in the SOL and pedestal regions. It is observed that the ECE radiation $\delta {T_{e{\text{,rad}}}}$, which is located in the SOL region according to the cold and optically thick plasma resonance assumption, is extremely sensitive to MHD radial displacements near the separatrix, exhibiting radiation inversion to $\delta {T_{e{\text{,rad}}}}$ at the pedestal. Here, the radiation inversion refers to the opposite phase between the radiation fluctuation at the pedestal and the radiation fluctuation at the SOL. Consequently, the quasi-coherent MHD (QCM), which displays a radiation inversion, is found to be consistent with an MHD radial structure that has a strong displacement near the separatrix. In contrast, the edge harmonic oscillation (EHO), which displays weak or no inversion, is found to be consistent with an MHD radial displacement structure peaking at the pedestal top. The ECEI data, interpreted with synthetic ECE, are in qualitative agreement with beam emission spectroscopy measurements on DIII-D for the relative radial extent and localization of the EHO and QCM. The high sensitivity of ECE radiation to separatrix displacements can be used to detect turbulence or MHD fluctuations near the separatrix, which may affect the transport across the separatrix and the wetted area in the divertor. The $\delta {T_{e{\text{,rad}}}}$ inversion measured with an ECE or ECEI system potentially provides important information on the magnetic field $\left| {{B_{{\text{sep}}}}} \right|$ at the separatrix, which helps constrain the pedestal equilibrium reconstruction and achieve an unambiguous mapping of the ECE/ECEI system with respect to the separatrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Operational space and performance limiting events in the first physics campaign of MAST-U

Abstract The MAST-U fusion plasma research device, an upgrade to the Mega Amp Spherical Tokamak, has recently completed its first campaign of physics operation. MAST-U operated with Ohmic, or one or two neutral beams for heating, at 400–800 kA plasma current, in conventional or ‘SuperX’ divertor configurations. Equilibrium reconstructions provide key plasma physics parameters vs. time for each discharge, and diagrams are produced which show where the prevalence of operation occurs as well as the limits in various operational spaces. When compared to stability limits, the operation of MAST-U so far has generally stayed out of the low q , low density instability region, and below the high density Greenwald limit, high beta global stability limit, and high elongation vertical stability limit. MAST-U still has the potential to reach higher elongation, which could benefit the plasma performance. Despite the majority of operations happening below established stability limits, disruptions do occur in the flat-top phase of MAST-U plasmas. The reasons for these disruptions are highlighted, and possible strategies to avoid them and to extend the operational space of MAST-U in future campaigns are discussed.

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