SEARCH · Search NASA

Results for “Equilibria”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Dynamic re-equilibration controlled multi-step transformations in (Mn, Fe) 2 (P, Si) alloys

In (Mn, Fe) 2 (P, Si) alloys crystallographic first-order phase transformations enable strong coupling of magnetic and entropic properties, potentially leading to high-efficiency energy generation and refrigeration applications. Although hysteresis losses that limit these applications can be reduced through careful control of alloy composition, compositional tuning can also unfavorably influence transformation temperatures and magnetocaloric coupling strength. Hence, exploration of additional processing variables enabling independent control of transformation properties is crucial. In this work, we investigate the role of thermal history as an additional processing variable, exploiting thermally-activated mechanisms to control properties of non-diffusive transformations in (Mn, Fe) 2 (P, Si) alloys. In so doing, we report an unusual transformation-splitting phenomenon following annealing at intermediate times, where a single well-defined magneto-structural transformation evolves towards a multi-step transformation with individual steps occurring at multiple distinct temperatures. Here, on longer annealing at the same temperatures, single-step transformation behavior is recovered. Through additional magnetic and crystallographic characterization, we show that the thermal history-controlled multi-step behavior results from sluggish thermally-activated diffusion. The two-step transformation corresponds to non-equilibrium bimodal composition distributions in the transforming phase, and these develop through a dynamic re-equilibration process as the alloy passes relatively slowly between different thermal equilibria. Together, these results suggest that thermal history primarily controls the transformation properties of (Mn, Fe) 2 (P, Si) alloys indirectly through the composition of one or more transforming phases. Additional investigations are needed to develop thermal history processing for decoupling hysteresis control from other transformation properties.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

Stabilizing effect of enhanced resistivity on peeling-ballooning instabilities on EAST

Previous stability analysis of NSTX equilibrium with lithium-conditioning demonstrates that enhanced resistivity due to the increased effective charge number Z eff (i.e. increased impurity level) can provide a stabilizing effect on low-n edge localized modes. This paper extends the resistivity stabilizing effect to the intermediate-n peeling-ballooning (PB) instabilities with the linear stability analysis of EAST high-confinement mode equilibria in NIMROD two-fluid calculations. However, in this work, the resistivity stabilizing effect on PB instabilities in the EAST tokamak appears weaker than that found in NSTX. This work may provide better insight into the physical mechanism behind the beneficial effects of impurities on pedestal stability.

EAST tokamak↗

DOE OSTI2021

Computation of linear MHD instabilities with the multi-region relaxed MHD energy principle

Over the last decade, a variational principle based on a generalisation of Taylor's relaxation, referred to as multi-region relaxed magnetohydrodynamics (MRxMHDs) has been developed. The numerical solutions of the MRxMHD equilibria have been constructed using the Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502). In principle, SPEC could also be established to describe the MRxMHD stability of an equilibrium. In this work, a theoretical framework is developed to relate the second variation of the energy functional to the so-called Hessian matrix, enabling the prediction of MHD linear instabilities of cylindrical plasmas, and is implemented in SPEC. The negative and positive eigenvalues of the Hessian matrix predict the stability of an equilibrium. Verification studies of SPEC stability results with the M3D-C 1 code and the tearing mode $\Delta^{^{\prime}}$ criterion have been conducted for ideal and resistive MHD instabilities, respectively, in a pressureless cylindrical tokamak, and have shown good agreement. Overall, our stability analysis is a critical step towards understanding the MHD stability of three-dimensional MHDs where nested flux surfaces, magnetic islands and stochastic regions co-exist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

A comparative study of internal kink stability in EU DEMO designs with negative and positive triangularity

Abstract Internal kink (IK) instability is investigated for European demonstration fusion reactor (EU DEMO) plasmas in both negative triangularity (NT) and positive triangularity (PT) configurations. For NT plasmas, the IK becomes more unstable as an ideal conformal wall moves away from the plasma boundary, with the mode growth rate saturating at the wall radial location of about b / a = 1.5 , where a is the plasma minor radius and b the wall radial location. The plasma resistivity destabilizes the IK mode. The effect of sub-sonic toroidal plasma flow is sufficiently weak and can thus be ignored for these EU DEMO equilibria. These results are consistent with those for PT plasmas, albeit with larger mode growth rate in the NT configuration. Both perturbative and self-consistent magneto-hydrodynamic (MHD)-kinetic hybrid calculations predict (partial) stabilization of the IK modes in both NT and PT configurations, with inclusion of various kinetic contributions. Precessional drift motion of trapped fusion-born alphas in EU DEMO produces weak stabilization to the IK mode. Stronger stabilization occurs with the toroidal precession of trapped thermal particles (ions and electrons) and the bounce-transit motion of thermal ions. The stabilization is similar between the NT and PT configurations, due to the similarity of the mode eigenfunction (occupying a nearly circular region in the plasma core) despite the sign difference in the triangularity. The non-perturbative MHD-kinetic hybrid model predicts much less stabilization of the mode than the perturbative model, primarily due to the self-consistent determination of the mode eigenvalue in the former. Generally, no significant difference in the IK mode stability is found between the NT and PT plasmas in EU DEMO.

Physics↗

DOE OSTI2021

Numerical investigation of active control of tearing mode by magnetic coils and the role of Δ'

Magnetic feedback stabilization of the tearing mode (TM) is numerically investigated, utilizing the MARS-F code (Liu et al 2000 Phys. Plasmas 7 3681) for toroidal tokamak equilibria. With control coil configurations assumed in this study, magnetic feedback partially or fully stabilizes the TM, with either vanishing or finite equilibrium pressure. The best control is achieved by the combination of internal active coils and internal poloidal sensors. The internal and external tearing indices are evaluated for the close-loop system, based on the MARS-F computed mode eigenvalue and eigenfunction, respectively. In the absence of the favorable curvature effect, these two indices are real-valued and quantitatively agree well with each other. For the equilibrium with finite pressure gradient at the mode rational surface, the favorable average curvature effect becomes important and the close-loop tearing index also becomes complex-valued, partly due to interaction of the feedback system with the dissipative wall eddy current response. Isolating the inner layer and outer region response to magnetic feedback, with either proportional or proportional-derivative actions, allows to establish that feedback stabilization of the TM occurs mainly due to modification of the behavior of the external ideal solution, further confirming the analytic result reported in He et al 2021 Phys. Plasmas 28 012504.

Physics↗

DOE OSTI2023

Structure of pressure-gradient-driven current singularity in ideal magnetohydrodynamic equilibrium

Abstract Singular currents typically appear on rational surfaces in non-axisymmetric ideal magnetohydrodynamic (MHD) equilibria with a continuum of nested flux surfaces and a continuous rotational transition. These currents have two components: a surface current (Dirac δ -function in flux surface labeling) that prevents the formation of magnetic islands, and an algebraically divergent Pfirsch–Schlüter current density when a pressure gradient is present across the rational surface. On flux surfaces adjacent to the rational surface, the traditional treatment gives the Pfirsch–Schlüter current density scaling as J ∼ 1 / Δ ι , where Δ ι is the difference of the rotational transform relative to the rational surface. If the distance s between flux surfaces is proportional to Δ ι , the scaling relation J ∼ 1 / Δ ι ∼ 1 / s will lead to a paradox that the Pfirsch–Schlüter current is not integrable. In this work, we investigate this issue by considering the pressure-gradient-driven singular current in the Hahm–Kulsrud–Taylor problem, which is a prototype for singular currents arising from resonant magnetic perturbations. We show that not only the Pfirsch–Schlüter current density but also the diamagnetic current density are divergent as ∼ 1 / Δ ι . However, due to the formation of a Dirac δ -function current sheet at the rational surface, the neighboring flux surfaces are strongly packed with s ∼ ( Δ ι ) 2 . Consequently, the singular current density J ∼ 1 / s , making the total current finite, thus resolving the paradox. Furthermore, the strong packing of flux surfaces causes a steepening of the pressure gradient near the rational surface, with ∇ p ∼ d p / d s ∼ 1 / s . In general non-axisymmetric MHD equilibrium, contrary to Grad’s conjecture that the pressure profile is flat around densely distributed rational surfaces, our result suggests a pressure profile that densely steepens around them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Macroscopic trends of linear tearing stability in cylindrical current profiles

Abstract The likelihood of realising tokamak power-plants will be greatly improved by the discovery of high-gain equilibria that resist the formation of small islands and hence avoid the disruptive neoclassical tearing mode. We propose a series of studies to understand how simple tokamak design can leverage aspects of tearing onset physics to maximise passive resistance to island formation. Here we investigate the variation that current profiles can bring about in preventing tearing onset through the cylindrical linear tearing stability parameter Δ ′ . A database of 159148 realistic pilot-plant current profiles was generated with Monte Carlo sampling, and the distribution of Δ ′ values was linked with interpretable profile characteristics. In agreement with prior theoretical and experimental studies, Δ ′ was found to be strongly correlated with the existence and steepness of a local toroidal current well or hill, with the former destabilising and the latter stabilising. In the absence of these two cases, the remaining Δ ′ values were linearly bounded by the toroidal current gradient at the rational surface.

Physics↗

DOE OSTI2024

Core inductive electric field during sawtooth crashes on DIII-D

Sawtooth crashes on tokamak plasmas exhibit relaxation much faster than resistive time scales via a mechanism not fully understood. Using core magnetic measurements from the Radial Interferometer-Polarimeter (RIP) diagnostic on the DIII-D tokamak, Grad–Shafranov equilibria constrained by internal magnetic measurements that have high time resolution (<1μs) can be computed, allowing analysis of how equilibrium parameters such as safety factor q, current density J, and parallel electric field E ∥ , particularly on-axis, evolve. At the sawtooth crash, on-axis safety factor q 0 is observed to rise by 5% but remain below 1 throughout the cycle, and on-axis current density J 0 is observed to drop by 5%. On-axis parallel electric field E ∥ (0) is found to be balanced by ηJ 0 (resistivity times on-axis current density) except during the 200 µs crash period, where E ∥ (0) reaches 22 V m –1 , exceeding ηJ 0 by a factor of more than 2000. These first measurements in tokamak plasmas verify that generalized Ohm's law is not balanced during the crash by resistive effects alone; this is a finding expected due to the relaxation being much faster than resistive timescales. As a result, measurement of the electric field during the tokamak sawtooth serves to illuminate the physical mechanisms at work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Accuracy of kinetic equilibrium reconstruction of NSTX and NSTX-U plasmas and its impact on the transport and stability analysis

An accurate magnetohydrodynamic (MHD) equilibrium reconstruction is an essential starting point for stability and transport plasma analysis. Herein this work describes an approach for obtaining kinetic equilibrium reconstructions using the OMFIT framework, which has been applied for the first time to spherical tokamak data from NSTX and NSTX-U. The EFIT equilibrium solver is integrated with experimental data analysis procedures and subsequent TRANSP transport simulations to enhance the accuracy of the reconstruction, in particular, at the edge region, by adding constraints on the total pressure and current density profiles, based on the transport code solution. The accuracy of the equilibrium reconstruction depends on the uncertainty and number of constraints, as well as the choice of basis functions to represent the pressure and current density profiles. Improved fidelity of the equilibrium reconstruction is demonstrated by reducing the variability of the magnetic axis and boundary locations from several centimeters, for reconstructions based on magnetic and experimental pressure constraints, to only several millimeters, for kinetic reconstructions based on transport code constraints, when different representations of basis functions were tested. The variability of the safety factor on axis was reduced ten times in the same sensitivity study. The accuracy of the equilibrium reconstruction and subsequent mapping of the experimental kinetic profile data have a significant impact on the trapped gyro Landau fluid and linear CGYRO turbulence simulations, which predict different spectra of unstable modes and turbulent fluxes for cases with different numbers of constraints in the equilibrium reconstruction. Conversely, the stability analysis performed using the GATO code shows plasmas that are stable to n = 1 MHD modes in both equilibria using magnetic and experimental pressure constraints as well as the transport code constrained equilibrium. However, a scan of parameters away from these conditions shows considerable deviation in the threshold of unstable modes between these reconstructions. Therefore, for reliable plasma analysis and use in turbulence and stability calculations, a high-fidelity equilibrium reconstruction with accurate kinetic constraints based on transport code solutions is necessary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Resilient stellarator divertor characteristics in the Helically Symmetric eXperiment

Resilient divertor features connected to open chaotic edge structures in the Helically Symmetric eXperiment are investigated. For the first time, an expanded vessel wall was considered that would give space for implementation of a physical divertor target structure. The analysis was done for four different magnetic configurations with very different chaotic plasma edges. A resilient plasma wall interaction pattern was identified across all configurations. This manifests as qualitatively very similar footprint behavior across the different plasma equilibria. Overall, the resilient field lines of interest with high connection length L C lie within a helical band along the wall for all configurations. This resiliency can be used to identify the best location of a divertor. The details of the magnetic footprint’s resilient helical band is subject to specific field line structures which are linked to the penetration depth of field lines into the plasma and directly influence the heat and particle flux patterns. The differences arising from these details are characterized by introducing a new metric, the minimum radial connection min $(\delta_N)$ of a field line from the last closed flux surface. The relationship, namely the deviation from a scaling law, between min $(\delta_N)$ and L C of the field lines in the plasma edge field line behavior suggests that the field lines are associated with structures such as resonant islands, cantori, and turnstiles. This helps determine the relevant magnetic flux channels based on the radial location of these chaotic edge structures and the divertor target footprint. These details will need to be taken into account for resilient divertor design.

divertor↗

DOE OSTI2025

Large-scale tearing-mode hazard function analysis with standard matched equilibrium reconstructions

The association between features from standard tokamak equilibrium reconstructions and the onset of n = 1 tearing modes (TMs) is analyzed at scale. The TM onset rate is directly modeled with a ‘hazard’ function which gives the expected number of onsets (per unit time spent) in a given equilibrium parameter region. In particular the different statistical modeling performance achieved for magnetics-only reconstructions and motional Stark effect (MSE) enhanced reconstructions is studied. It is observed that a better hazard model for the TM onset rate can be built with the MSE-enhanced equilibria compared to the matched magnetics-only situation. This advantage disappears if internal profile details are withheld from the matched analysis. Plausibility of the hazard function is further demonstrated with visualizations of global trends in the operational space, and time-traces from specific tokamak discharges. As a result, TMs typically degrade tokamak plasma performance and may lead to plasma termination, motivating this statistical study.

equilibrium↗

DOE OSTI2025

Omnigenous stellarators with improved ideal and kinetic ballooning stability

Omnigenity is a property of a magnetic field which ensures confinement of trapped particles. It is a necessary requirement for any high-performance stellarator. After creating an omnigenous equilibrium, one must also ensure reduced transport resulting from kinetic and magnetohydrodynamic (MHD) instabilities. To this end, we leverage the GPU-accelerated DESC optimization suite, which is used to design stable, finite-β omnigenous equilibria with poloidal, toroidal, and helical symmetry, achieving Mercier, ideal ballooning, and as a consequence, improved kinetic ballooning stability. We discover stellarators with second stability, a regime of large pressure gradient where an equilibrium becomes ideal ballooning stable, and demonstrate and explore both using theory and gyrokinetic simulations the connection between ideal and kinetic ballooning stability.

optimization↗

DOE OSTI2020

Physics-guided machine learning approaches to predict the ideal stability properties of fusion plasmas

One of the biggest challenges to achieve the goal of producing fusion energy in tokamak devices is the necessity of avoiding disruptions of the plasma current due to instabilities. The disruption event characterization and forecasting (DECAF) framework has been developed in this purpose, integrating physics models of many causal events that can lead to a disruption. Two different machine learning approaches are proposed to improve the ideal magnetohydrodynamic (MHD) no-wall limit component of the kinetic stability model included in DECAF. First, a random forest regressor (RFR), was adopted to reproduce the DCON computed change in plasma potential energy without wall effects for a large database of equilibria from the national spherical torus experiment (NSTX). This tree-based method provides an analysis of the importance of each input feature, giving an insight into the underlying physics phenomena. Secondly, a fully-connected neural network has been trained on sets of calculations with the DCON code, to get an improved closed form equation of the no-wall β limit as a function of the relevant plasma parameters indicated by the RFR. The neural network has been guided by physics theory of ideal MHD in its extension outside the domain of the NSTX experimental data. The estimated value has been incorporated into the DECAF kinetic stability model and tested against a set of experimentally stable and unstable discharges. Moreover, the neural network results were used to simulate a real-time stability assessment using only quantities available in real-time. Finally, the portability of the model was investigated, showing encouraging results by testing the NSTX-trained algorithm on the mega ampere spherical tokamak (MAST).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

The role of edge resonant magnetic perturbations in edge-localized-mode suppression and density pump-out in low-collisionality DIII-D plasmas

Two-fluid nonlinear MHD simulations using the TM1 code demonstrate that the formation of magnetic islands at the top and bottom of the H-mode pedestal, together with the strong screening of resonant fields in the gradient region of the pedestal, can account for ELM suppression and density pump-out by n = 2 resonant magnetic perturbations (RMPs) in low-collisionality DIII-D ITER similar shape (ISS) plasmas. Using experimentally relevant transport coefficients, neoclassical resistivity, electron collisionality, and RMP amplitudes, nonlinear MHD simulations reproduce the observed level of density reduction (density pump-out) in DIII-D due to the formation of narrow magnetic islands in the resistive foot of pedestal. For large amplitude RMPs ( B r / B t > 1 × 10 -4 ) simulations predict resonant field penetration and resulting in pressure reduction at the top of the pedestal, consistent with experimental observations at the onset of ELM suppression. The predicted reduction in the height and width of the pedestal by magnetic island enhanced collisional transport provides a quantitative mechanism for the stabilization of the peeling-ballooning mode (PBM). Importantly, these simulations predict strong screening of resonant fields in the steep gradient region of the pedestal due to strong E × B rotation and diamagnetic flows. Yet, if the plasma resistivity is made artificially larger (~ 10X) than neoclassical, the simulations predict magnetic stochasticity throughout the plasma edge and the collapse of the pedestal due to the reduction in the penetration threshold with increasing resistivity. A scaling relation for resonant field penetration at the pedestal top, using several hundred nonlinear simulations, reproduces the density and E × B rotation dependence of the ELM suppression threshold observed in DIII-D. Simulations using ITER model equilibria indicate that the penetration threshold at the top of the ITER pedestal will be lower than in DIII-D due to the anticipated lower perpendicular flow velocities in ITER, resulting in the weaker screening of resonant fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Pedestal collapse by resonant magnetic perturbations

Pedestal collapse (i.e., the complete loss of the edge transport barrier (ETB)) in DIII-D H-mode plasmas occurs when resonant magnetic perturbations (RMPs) penetrate the steep gradient region at the plasma edge. Normally, RMP driven magnetic islands can occur at the top and bottom of the H-mode pedestal and these islands generate conditions consistent with edge-localized-mode (ELM) suppression and density pump-out, respectively, based on nonlinear two-fluid MHD simulations. In contrast, MHD simulations show that the steep pressure gradient region between the top and bottom of the DIII-D pedestal is generally immune to resonant field penetration due to large local E × B and diamagnetic flows. By this fortuitous circumstance, the edge-transport-barrier and H-mode confinement can be maintained while achieving ELM suppression. However, pedestal collapse can occur in DIII-D when the screening flows are inadequate to prevent field penetration in the steep gradient region of the pedestal. Non-linear two-fluid MHD simulations support the role of resonant field penetration in pedestal collapse for DIII-D H-mode plasmas with weak edge E × B and diamagnetic screening flows. ITER will likely have weaker edge screening flows than present experiments due to its much larger size, making it more susceptible to resonant field penetration in the steep gradient region of the pedestal. Analysis of model ITER equilibria demonstrates that resonant field penetration in the steep pressure gradient region is possible for RMP levels of the order required for ELM suppression. The effect of such penetration on the ITER pedestal will depend sensitively on the resulting degree of island overlap.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Pedestal stability analysis on MAST in preparation for MAST-U

In preparation for the upcoming MAST-U campaign, pedestal stability of spherical tokamaks is revisited by investigating standard H-mode discharges on MAST. As a step beyond previous studies, both ion and electron profiles are used for obtaining equilibria and a diverse set of pedestals is evaluated. Stability analysis with the ELITE and CGYRO codes shows that MAST pedestals are constrained by kinetic ballooning modes and medium toroidal mode number peeling-ballooning modes, with most unstable modes ranging from n = 25 to n = 45. In discharges with a steep q profile at the edge a larger number of poloidal harmonics is excited for each toroidal mode. In conclusion, a comparison with discharges on DIII-D with matched shape and similar non-dimensional parameters indicates that the increased shear at lower aspect ratio stabilizes low n peeling modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Approach to nonlinear magnetohydrodynamic simulations in stellarator geometry

The capability to model the nonlinear magnetohydrodynamic (MHD) evolution of stellarator plasmas is developed by extending the M3D-C 1 code to allow non-axisymmetric domain geometry. We introduce a set of logical coordinates, in which the computational domain is axisymmetric, to utilize the existing finite-element framework of M3D-C 1 . A C 1 coordinate mapping connects the logical domain to the non-axisymmetric physical domain, where we use the M3D-C 1 extended MHD models essentially without modifications. We present several numerical verifications on the implementation of this approach, including simulations of the heating, destabilization, and equilibration of a stellarator plasma with strongly anisotropic thermal conductivity, and of the relaxation of stellarator equilibria to integrable and non-integrable magnetic field configurations in realistic geometries.

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↗