Search NASA⌕ Search

SEARCH · Search NASA

Results for “MAGNETIC ISLANDS”

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.

575 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Equatorial Dynamics Observed by Rocket, Radar, and Satellite During the CADRE/MALTED Campaign: Programmatics and small-scale fluctuations - 1

In August 1994, the Mesospheric and Lower Thermospheric Equatorial Dynamics (MALTED) Program was conducted from the Alcantara rocket site in northeastern Brazil as part of the International Guard Rocket Campaign to study equatorial dynamics, irregularities, and instabilities in the ionosphere. This site was selected because of its proximity to the geographic (2.3 deg S) and magnetic (approx. 0.5 deg S) equators. MALTED was concerned with planetary wave modulation of the diurnal tidal amplitude, which exhibits considerable amplitude variability at equatorial and subtropical latitudes. Our goals were to study this global modulation of the tidal motions where tidal influences on the thermal structure are maximum, to study the interaction of these tidal structures with gravity waves and turbulence at mesopause altitudes, and to gain a better understanding of dynamic influences and variability on the equatorial middle atmosphere. Four (two daytime and two nighttime) identical Nike-Orion payloads designed to investigate small-scale turbulence and irregularities were coordinated with 20 meteorological falling-sphere rockets designed to measure temperature and wind fields during a 10-day period. These in situ measurements were coordinated with observations of global-scale mesospheric motions that were provided by various ground based radars and the Upper Atmosphere Research Satellite (UARS) through the Coupling and Dynamics of Regions Equatorial (CADRE) campaign. The ground-based observatories included the Jicamarca radar observatory near Lima, Peru, and medium frequency (MF) radars in Hawaii, Christmas Island, and Adelaide. Since all four Nike-Orion flights penetrated and overflew the electrojet with apogees near 125 km, these flights provided additional information about the electrodynamics and irregularities in the equatorial ionospheric E region and may provide information on wave coupling between the mesosphere and the electrojet. Simultaneous with these flights, the CUPRI 50-MHz radar (Cornell University) provided local sounding of the electrojet region. A description of the campaign logistics and the measurements performed with the Nike-Orion instrumentation and their implications for turbulence due to gravity waves and tidal instability in the mesosphere and lower thermosphere (MLT) are presented here. From a study of electron density fluctuations measured by rocket probes, we have found evidence for equatorial mesospheric neutral-atmospheric turbulence between 85 and 90 km. Furthermore, falling-sphere data imply that gravity wave breaking was a source for this turbulence. Mean motions and the various planetary, tidal, and gravity wave structures and their coherence and variability are the subjects of a companion paper.

Goldberg, Richard A.↗
DOE OSTI2024

New Muon Campus Simulations for the Muon $g\textrm{-}2$ Experiment at Fermilab

The Muon $g\textrm{-}2$ Experiment (E989) at Fermilab conducted high-precision measurements of the muon anomalous magnetic moment $a_{\mu}$ using a storage ring from 2018 to 2023, achieving a remarkable precision of $200\:\mathrm{ppb}$ over Runs 1-3, with analyses for Runs 4-6 ongoing. A comprehensive understanding of the storage ring's beam dynamics and its accurate simulations are crucial for achieving the experiment's ambitious goals. One of the requirements for this effort is a very detailed knowledge of the phase space distribution of the beam. To address this requirement, we performed high-statistics simulations of the Muon $g\textrm{-}2$ Target Station (AP0) and the Muon Campus beamlines: M2 and M3, followed by the Delivery Ring, and then M4 and M5. The resulting muon distribution at the end of the M5 beamline from our previous $3\times10^{12}$ protons-on-target (PoT) simulation serves as an essential input for the storage ring simulations. In 2024, to facilitate the analyses of Runs 4-6, we have updated our Muon Campus models and re-optimised certain parameters to reflect the operational currents and wire chamber measurements of the beam. For these optimisations, we employed the heterogeneous island method, implemented in our evolutionary optimisation tool, glyfada. This key update addresses the need to use the best possible beam for the storage ring simulations and thus supports the experiment's overall precision. The Muon $g\textrm{-}2$ Experiment at Fermilab has successfully achieved its $70\:\mathrm{ppb}$ systematic uncertainty goal and collected 21 times more data than its predecessor at BNL. The updated and improved Muon Campus models and simulations not only facilitate the experiment's efforts to potentially resolve the current tension between experimental measurements and theoretical predictions of $a_{\mu}$, but also provide a basis for future simulations for the Mu2e Experiment (E-973), which utilises shared Muon Campus beamlines.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Fast ion transport in the quasi-single helical reversed-field pinch

The reversed-field pinch (RFP) can spontaneously transition from an axisymmetric magnetic topology to a 3D-helical geometry. Investigations on fast ion transport associated with energetic particle driven Alfvén instabilities, tearing mode induced stochasticity, and neoclassical effects have been performed on the Madison Symmetric Torus. STELLGAP produced shear-Alfvén continua seeded with V3FIT 3D-equilibrium reconstructions describe the response of Alfvénic bursting activity as a direct consequence of the equilibrium change on the fast ion resonance. Furthermore, far infrared interferometry resolved electron density perturbations associated with the bursts provide a spatial measurement of the mode structure and support the reconstructions. The bursts produce no global resonant fast ion transport; however, their disappearance at a high core-resonant amplitude implies other transport mechanisms at play. Additionally, neutral particle analysis and neutron signals suggest fast ion losses at sufficient core tearing mode strength, supporting the lack of Alfvénic activity. The guiding-center code ORBIT corroborates rapid fast ion loss times in the helical state largely as a consequence of remnant tearing modes. Additionally, ORBIT simulations demonstrate little neoclassical enhancement of particle transport. While superbanana orbits may exist, the growth in the core-resonant fast ion island and the associated secondary mode overlap govern the largest transport process, leading to robust fast ion losses in the 3D-RFP.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Fast Ion Transport in the Three-Dimensional Reversed-Field Pinch

Here, we report on the first comprehensive experimental and numerical study of fast ion transport in the helical reversed-field pinch (RFP). Classical orbit effects dominate the macroscopic confinement properties. The strongest effect arises from growth in the dominant fast ion guiding-center island, but substantial influence from remnant subdominant tearing modes also plays a critical role. At the formation of the helical RFP, neutron flux measurements indicate a drastic loss of fast ions at sufficient subdominant mode amplitudes. Simulations corroborate these measurements and suggest that subdominant tearing modes strongly limit fast ion behavior. Previous work details a sharp thermal transport barrier and suggests the helical RFP as an Ohmically heated fusion reactor candidate; the enhanced transport of fast ions reported here identifies a key challenge for this scheme, but a workable scenario is conceivable with low subdominant tearing mode amplitudes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Numerical investigation of plasma edge transport and limiter heat fluxes in Wendelstein 7-X startup plasmas with EMC3-EIRENE

Here the results of a first systematic assessment of plasma edge transport processes for the limiter startup configuration at Wendelstein 7-X are presented. This includes an investigation of transport from intrinsic and externally injected impurities and their impact on the power balance and limiter heat fluxes. The fully 3D coupled plasma fluid and kinetic neutral transport Monte Carlo code EMC3-EIRENE is used. The analysis of the magnetic topology shows that the poloidally and toroidally localized limiters cause a 3D helical scrape-off layer (SOL) consisting of magnetic flux tubes of three different connection lengths L C . The transport in the helical SOL is governed by L C as topological scale length for the parallel plasma loss channel to the limiters. A clear modulation of the plasma pressure with L C is seen. The helical flux tube topology results in counter streaming sonic plasma flows. The heterogeneous SOL plasma structure yields an uneven limiter heat load distribution with localized peaking. Assuming spatially constant anomalous transport coefficients, increasing plasma density yields a reduction of the maximum peak heat loads from 12 MWm -2 to 7.5 MWm -2 and a broadening of the deposited heat fluxes. The impact of impurities on the limiter heat loads is studied by assuming intrinsic carbon impurities eroded from the limiter surfaces with a gross chemical sputtering yield of 2 . The resulting radiative losses account for less than 10% of the input power in the power balance with marginal impact on the limiter heat loads. It is shown that a significant mitigation of peak heat loads, 40–50%, can be achieved with controlled impurity seeding with nitrogen and neon, which is a method of particular interest for the later island divertor phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

On the Effect of Beating during Nonlinear Frequency Chirping

Spectroscopic analyses of energetic particle (EP) driven bursts of MHD fluctuations in magnetically confined plasmas often exhibit chirps that occur simultaneously in groups of two or more. While the superposition of oscillations at multiple frequencies necessarily causes beating in the signal acquired by a localized external probe, self-consistent hybrid simulations of chirping EP modes in a JT-60U tokamak plasma have demonstrated the possibility of global beating, where the mode's electromagnetic field vanishes globally between beats and reappears with opposite phase [Bierwage et al., Nucl. Fusion 57, 016036 (2017)]. This implies that there can be a single coherent field mode that oscillates at multiple frequencies simultaneously when it is resonantly driven by multiple density waves in EP phase space. Conversely, this means that the EP density waves are mutually coupled and interfere with each other via the jointly driven field, a mechanism ignored in some theories of chirping. In this thesis-style treatise, we study the role of field pulsations in general and beating in particular using the Hamiltonian guiding center orbit-following code ORBIT with a reduced wave-particle interaction model in realistic geometry. Beating is found to drive the evolution of EP phase space structures. Additionally, a key mechanism is the pulsation of effective phase space islands combined with the alternation of their effective O- and X-points due to phase jumps between each beat. Observations: (1) Beating causes density wave fronts to advance radially in a pulsed manner and the resulting chirps become staircase-like. (2) The pulsations facilitate convective transfer of material between neighboring layers of phase space density waves. On the one hand, this may inhibit the early detachment of solitary phase space vortices. On the other hand, it facilitates the accumulation of hole and clump fragments into larger structures. (3) Long-range chirping is observed when massive holes or clumps detach and drift away from the turbulent belt around the seed resonance. It is remarkable that the detached vortices remain robust and,on average, maintain their concentric nested layers while being visibly perturbed by the field's continued beating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Gradient-based surface nuclear magnetic resonance for groundwater investigation

In medical magnetic resonance imaging, spatial localization (imaging) is based upon the application of controlled magnetic field gradients on top of the main magnetic field to spatially modulate the frequency and/or phase of the nuclear magnetic resonance (NMR) signal across the volume of investigation. In this work, we have applied similar physical principles to produce controlled magnetic field gradients during surface NMR-based groundwater investigations. In this approach, a gradient pulse of variable amplitude or duration is applied immediately after the excitation pulse to cause predictable phase encoding of the NMR signal as a function of depth. This approach is also applicable to emerging surface NMR detection methods that use a prepolarization field with fast nonadiabatic turn-off to generate detectable NMR signals from the shallow subsurface. In this case, the gradient pulse is applied after terminating the prepolarization field and provides a heretofore unavailable means of localizing the NMR response as a function of depth. The application of gradients can also be combined with tip-angle-based modulation to yield higher imaging resolution than can be achieved through either gradient- or tip-angle-based imaging alone. We implemented this new gradient-based capability into a surface NMR gradient generation accessory that is compatible with the GMR-Flex instrument and developed surface NMR-specific forward modeling and linear inverse models. We validated the accuracy of this novel gradient-based sNMR technology using computer simulations, experiments using a small pool filled with a discrete layer of bulk water, and field experiments at well-characterized groundwater test sites along Ebey Island, WA, and Larned, KS. The gradient-based sNMR imaging observations were compared with high-resolution direct push NMR results observed at these sites. The results of computer simulations and field experiments indicate improvements in both the detection (signal-to-noise ratio) and spatial resolution of shallow subsurface water content using gradient-based surface NMR, compared with traditional surface NMR imaging methods.

Geochemistry & Geophysics↗
DOE OSTI2016

Hurricane Resilient Wind Plant Concept Study Final Report

Hurricanes occur over much of the U.S. Atlantic and Gulf coasts, from Long Island to the U.S.-Mexico border, encompassing much of the nation's primary offshore wind resource. Category 5 hurricanes have made landfall as far north as North Carolina, with Category 3 hurricanes reaching New York with some frequency. Along the US West coast, typhoons strike with similar frequency and severity. At present, offshore wind turbine design practices do not fully consider the severe operating conditions imposed by hurricanes. Although universally applied to most turbine designs, International Electrotechnical Commission (IEC) standards do not sufficiently address the duration, directionality, magnitude, or character of hurricanes. To assess advanced design features that could mitigate hurricane loading in various ways, this Hurricane-Resilient Wind Plant Concept Study considered a concept design study of a 500-megawatt (MW) wind power plant consisting of 10-MW wind turbines deployed in 25-meter (m) water depths in the Western Gulf of Mexico. This location was selected because hurricane frequency and severity provided a unique set of design challenges that would enable assessment of hurricane risk and projection of cost of energy (COE) changes, all in response to specific U.S. Department of Energy (DOE) objectives. Notably, the concept study pursued a holistic approach that incorporated multiple advanced system elements at the wind turbine and wind power plant levels to meet objectives for system performance and reduced COE. Principal turbine system elements included a 10-MW rotor with structurally efficient, low-solidity blades; a lightweight, permanent-magnet, direct-drive generator, and an innovative fixed substructure. At the wind power plant level, turbines were arrayed in a large-scale wind power plant in a manner aimed at balancing energy production against capital, installation, and operation and maintenance (O&M) costs to achieve significant overall reductions in COE.

17 WIND ENERGY↗
DOE OSTI2021

MHD simulations of cold bubble formation from 2/1 tearing mode during massive gas injection in a tokamak

Massive gas injection (MGI) experiments have been carried out in many tokamaks to study disruption dynamics and mitigation schemes. Two events often observed in those experiments are the excitation of the m = 2, n = 1 magnetohydrodynamic mode, and the formation of cold bubble structure in the temperature distribution before the thermal quench (TQ). Here m is the poloidal mode number, n the toroidal mode number. The physics mechanisms underlying those phenomena, however, have not been entirely clear. In this work, our recent NIMROD simulations of the MGI process in a tokamak have reproduced the main features of both events, which has allowed us to examine and establish the causal relation between them. In these simulations, the 3/1 and 2/1 islands are found to form successively after the arrival of impurity ion cold front at the corresponding q = 3 and q = 2 rational surfaces. At the interface between impurity and plasma, a local thin current sheet forms due to an enhanced local pressure gradient and moves inward following the gas cold front, this may contribute to the formation of a dominant 2/1 mode. Following the growth of the 2/1 tearing mode, the impurity penetration into the core region inside the q = 2 surface gives rise to the formation of the cold bubble temperature structure and initiates the final TQ. Here, a subdominant 1/1 mode developed earlier near the q = 1 surface alone does not cause such a cold bubble formation, however, the exact manner of the preceding impurity penetration depends on the nature of the 1/1 mode: kink-tearing or quasi-interchange.

2/1 tearing mode↗
DOE OSTI2024

Numerical simulations of laser-driven experiments of ion acceleration in stochastic magnetic fields

We present numerical simulations used to interpret laser-driven plasma experiments at the GSI Helmholtz Centre for Heavy Ion Research. The mechanisms by which non-thermal particles are accelerated in astrophysical environments, e.g., the solar wind, supernova remnants, and gamma ray bursts, is a topic of intense study. When shocks are present, the primary acceleration mechanism is believed to be first-order Fermi, which accelerates particles as they cross a shock. Second-order Fermi acceleration can also contribute, utilizing magnetic mirrors for particle energization. Despite this mechanism being less efficient, the ubiquity of magnetized turbulence in the universe necessitates its consideration. Another acceleration mechanism is the lower-hybrid drift instability, arising from gradients of both density and magnetic field, which produce lower-hybrid waves with an electric field that energizes particles as they cross these waves. With the combination of high-powered laser systems and particle accelerators, it is possible to study the mechanisms behind cosmic-ray acceleration in the laboratory. In this work, we combine experimental results and high-fidelity three-dimensional simulations to estimate the efficiency of ion acceleration in a weakly magnetized interaction region. We validate the FLASH magneto-hydrodynamic code with experimental results and use OSIRIS particle-in-cell code to verify the initial formation of the interaction region, showing good agreement between codes and experimental results. We find that the plasma conditions in the experiment are conducive to the lower-hybrid drift instability, yielding an increase in energy ΔE of ~ 264 keV for 242 MeV calcium ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Nonlinear modeling of the scaling law for the m/n = 3/2 error field penetration threshold

The scaling law for the m/n=3/2 error field (EF) penetration threshold is predicted numerically based on nonlinear single-fluid and two-fluid modeling using the TM1 code. The simulated penetration threshold of radial magnetic field b r at the plasma edge is scaled to the electron density n e , temperature T e , viscous time τ μ , toroidal field B t and the natural frequency ω in the form of b r /B t ∝n e αn T e αT τ μ αμ B t αB ω αω by scanning these parameters separately. Here, α n , α T , α μ , α B and α ω are the scaling coefficients on n e , T e , τ μ , B t and ω, respectively. Single-fluid modeling shows that the 3/2 EF threshold scales as b r /B t ∝n e 0.56 T e 0.6 τ μ -0.59 B t -1.15 ω, which is similar with the analytical scaling law in both the Rutherford and visco-resistive regimes. Yet, two-fluid modeling shows that the scaling law differs significantly in particular regarding the dependence on plasma rotation. In detail, the scaling coefficient α n on density decreases from 0.67 to 0.56 and α T on temperature decreases from 0.67 to 0.32, while α μ on viscous time is around -0.45 and α B on toroidal field decreases slightly from -1.15 to -1, when the ratio |ω E /ω *e | between plasma rotation frequency ω E and diamagnetic drift frequency ω *e varies from 0 to 10. Scans of the plasma rotation reveals that the penetration threshold linearly depends on the perpendicular rotation frequency (or natural frequency) ω ⊥e =ω E +ω *e , and there is a minimum in the required field amplitude when ω ⊥e 0. In addition, the enduring mystery of non-zero penetration threshold at zero plasma natural frequency in EF experiments is resolved by two-fluid simulations. We report that the very small island and smooth bifurcation in EF penetration near zero frequency is hard to detect in the experiment, leading to a finite penetration threshold within the capability of the experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Particle Acceleration and Fractional Transport in Turbulent Reconnection

We consider a large-scale environment of turbulent reconnection that is fragmented into a number of randomly distributed unstable current sheets (UCSs), and we statistically analyze the acceleration of particles within this environment. We address two important cases of acceleration mechanisms when particles interact with the UCS: (a) electric field acceleration and (b) acceleration by reflection at contracting islands. Electrons and ions are accelerated very efficiently, attaining an energy distribution of power-law shape with an index 1–2, depending on the acceleration mechanism. The transport coefficients in energy space are estimated from test-particle simulation data, and we show that the classical Fokker–Planck (FP) equation fails to reproduce the simulation results when the transport coefficients are inserted into it and it is solved numerically. The cause for this failure is that the particles perform Levy flights in energy space, while the distributions of the energy increments exhibit power-law tails. We then use the fractional transport equation (FTE) derived by Isliker et al., whose parameters and the order of the fractional derivatives are inferred from the simulation data, and solving the FTE numerically, we show that the FTE successfully reproduces the kinetic energy distribution of the test particles. We discuss in detail the analysis of the simulation data and the criteria that allow one to judge the appropriateness of either an FTE or a classical FP equation as a transport model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Relevance of anisotropy in the Kondo effect: Lessons from the symplectic case

A Kondo model with symplectic symmetry was recently put forward as the effective low-energy theory of a superconducting-island device coupled to multiple leads. This model, which possesses non-Fermi liquid physics and effective anyons, was argued to belong to the class of topological Kondo effects. Here, we clarify the extent of stability of its exotic fixed point using perturbative and numerical renormalization group in conjunction with bosonization and conformal field theory. In contrast to previous claims, we show that asymmetry in the coupling to the leads destabilizes the non-Fermi liquid. Other destabilizing perturbations include asymmetry in the superconducting pairing or internal energy of the individual quantum dots in the island. Nevertheless, these perturbations all generate the same relevant operators. Thus, only a small number of couplings need to be tuned individually, and these can be selected according to experimental convenience. Our results highlight a common misconception that anisotropy in single-channel Kondo couplings is always irrelevant. As demonstrated, relevant terms will emerge whenever the group generators do not span the full space of impurity operators. This calls for a more detailed inspection of models that exhibit this property, such as large-spin impurities and SO⁡( M ) Kondo models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Loss of energetic particles due to resistive wall mode instability in ITER

Effects of an unstable n = 1 (n is the toroidal mode number) resistive wall mode (RWM) on the energetic particle (EP) confinement and loss are numerically investigated, for an ITER steady state scenario with 10 MA plasma current and 5.3 T toroidal field. The eigenfunction of the RWM is computed, with the associated three-dimensional magnetic field perturbation superposed with the 2D equilibrium field for tracing the EP drift orbits. Considered are mono-energetic EPs at 0.5 MeV and 1 MeV for deuterium ions, and 3.5 MeV for fusion-born alphas, with a range of distribution in the particle pitch angle. Modeling finds that less than 20% of EPs can be lost to the limiting surface in ITER assuming a source distribution uniform in minor radius, due to an unstable RWM that produces 100 Gauss poloidal field perturbation at the outboard mid-plane just inside the (effective) resistive wall surface. On top of the initial prompt drift orbit loss for counter-current EPs, the RWM induced particle loss occurs on a one second time scale, which is comparable to the RWM growth time in ITER. The 'wetted' area, due to the lost EPs striking the limiting surface, is generally found to be large due to the RWM. This is a favorable prediction for ITER. Here, the loss distribution in the poloidal angle is more uniform for co-current EPs. Counter-current EPs experience outward orbit drift when launched from the low-field side, and tend to more often hit the bottom region of the limiting surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Tuning of the rotational transform in Wendelstein 7-X

The control of rotational transform in Wendelstein 7-X (W7-X) is key to the island divertor operation and safety of plasma facing components. The island divertor concept in W7-X relies on an edge flux surface with rotational transform of $t$ = 1 resonating with an intrinsic n/m = 5/5 resonance to form a five lobed island chain. This island chain intersects with divertor plates to give rise to the island divertor. Changes in the relative position of the rational surface and the divertor plates can result in changes in divertor performance, thus the control of the rotational transform is essential to operation of the W7-X device. During the first divertor campaign electromagnetic loads resulted in elastic deformations of the shaped modular stellarator coils. Such deformations made these coils more planar, reducing the vacuum rotational transform, subsequently shifting the $t$ = 1 resonance outward. Unintended plasma wall interactions provided the first clear evidence of this effect during plasma operation. Flux surface measurements were utilized to estimate the correct level of current in the planar coils for correction of $t$, and found to be around ~-100 A . Scans of the planar coil current for iota correction were performed during plasma operation. These measurements suggest planar coil currents between -250 and ~ 0 A would place the strike lines at the designed distance from the pumping gaps. Divertor Langmuir and upstream probe measurements confirm these estimates along with measurements of divertor neutral gas pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Sanibel Symposium in the Petascale-Exascale Computational Era

The 56 th Sanibel Symposium was held February 14-19 2016 at the King and Prince Hotel, St. Simons Island, GA. It successfully brought quantum chemists and chemical and condensed matter physicists together in presentations, posters, and informal discussions bridging those two communities. The Symposium has had a significant role in preparing generations of quantum theorists. As computational potency and algorithmic sophistication have grown, the Symposium has evolved to emphasize more heavily computationally oriented method development in chemistry and materials physics, including nanoscience, complex molecular phenomena, and even bio-molecular methods and problems. Given this context, the 56 th Sanibel meeting systematically and deliberately had sessions focused on exascale computation. A selection of outstanding theoretical problems that need serious attention was included. Five invited sessions, two contributed sessions (hot topics), and a poster session were organized with the exascale theme. This was a historic milestone in the evolution of the Symposia. Just as years ago linear algebra, perturbation theory, density matrices, and band-structure methods dominated early Sanibel Symposia, the exascale sessions of the 56 th meeting contributed a transformative influence to add structure and strength to the computational physical science community in an unprecedented way. A copy of the full program of the 56 th Symposium is attached. The exascale sessions were Linear Scaling, Non-Adabatic Dynamics, Interpretive Theory and Models, Computation, Software, and Algorithms, and Quantum Monte Carlo. The Symposium Proceedings will be published in Molecular Physics (2017). Note that the Sanibel proceedings from 2015 and 2014 were published as Molecular Physics vol. 114, issue 3-4 (2016) and vol. 113, issue 3-4 (2015) respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Disruption event characterization and forecasting in tokamaks

Disruption prediction and avoidance is a critical need for next-step tokamaks, such as ITER. Disruption Event Characterization and Forecasting (DECAF) research fully automates analysis of tokamak data to determine chains of events that lead to disruptions and to forecast their evolution allowing sufficient time for mitigation or complete avoidance of the disruption. Disruption event chains related to local rotating or global magnetohydrodynamic (MHD) modes and vertical instability are examined with warnings issued for many off-normal physics events, including density limits, plasma dynamics, confinement transitions, and profile variations. Along with Greenwald density limit evaluation, a local radiative island power balance theory is evaluated and compared to the observation of island growth. Automated decomposition and analysis of rotating tearing modes produce physical event chains leading to disruptions. A total MHD state warning model comprised of 15 separate criteria produces a disruption forecast about 180 ms before a standard locked mode detector warning. Single DECAF event analyses have begun on KSTAR, MAST, and NSTX/-U databases with thousands of shot seconds of device operation using from 0.5 to 1 × 10 6 tested sample times per device. An initial multi-device database comparison illustrates a highly important result that plasma disruptivity does not need to increase as βN increases. Global MHD instabilities, such as resistive wall modes (RWMs), can give the briefest time period of warning before disruption compared to other physics events. In an NSTX database with unstable RWMs, the mode onset, loss of boundary and current control, and disruption event warnings are found in all cases and vertical displacement events are found in 91% of cases. An initial time-dependent reduced physics model of kinetic RWM stabilization created to forecast the disruption chain predicts instability 84% of the time for experimentally unstable cases with a relatively low false positive rate. Instances of the disruption event chain analysis illustrate dynamics including H–L back transitions for rotating MHD and global RWM triggering events. Furthermore, disruption warnings are issued with sufficient time before the disruption (on transport timescales) to potentially allow active profile control for disruption avoidance, active mode control, or mitigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗