Search NASA⌕ Search

Engineering topics

Finn, John M.

Publications and source records attributed to Finn, John M..

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Stability and Control of Burning Tokamak Plasmas with Resistive Walls (Final Report)

This research has focused on quantitative prediction of the stability, control, and equilibrium state evolution in toroidal burning plasmas. The stability of long pulse burning plasmas is highly sensitive to the physics of resonant layers in the plasma, sources of momentum and flow, kinetic effects of energetic particles, and boundary conditions at the wall, including feedback control and error fields. In ITER in particular, the low toroidal flow equilibrium state, sustained primarily by energetic alpha particles from fusion reactions, will require the consideration of all of these key elements to predict quantitatively the stability and evolution. The principal investigators on this proposal are leading experts in the relevant theoretical and computational areas, and aimed to perform computations guided by analytic modeling, to address this physics in realistic configurations. The overall goal is to understand the key physics mechanisms that describe resistive toroidal burning plasmas, surrounded by a resistive wall, under active feedback control. With the physics of the energetic ions, resonant layers, resistive wall, and toroidal momentum transport included, this study will extend from recent publications in theory and simulation of individual effects and move toward predictive modeling for burning plasmas.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Stability and Control of Burning Tokamak Plasmas with Resistive Walls (Final Report)

Work supported by the grant, "Stability and Control of Burning Tokamak Plasmas with Resistive Walls" in the period Aug 2017 to 2023 can be cast in three general categories: (1) linear tearing regimes and the effects on the response to error fields and resistive walls: (2) nonlinear tearing modes in the presence of error fields, with locking; and (3) energetic particle effects. In the first category, we have completed a study of tearing mode regimes with extended MHD effects, including parallel dynamics and toroidal effects, notably favorable average curvature. This work has been presented at the 2018 Sherwood Conference. This work has been published in the ArXiv and is being prepared for publication in Physics of Plasmas. A systematic study of shaping in a tokamak plasma with resistive tearing layers and a resistive wall has been presented in the PhD dissertation by Dov Rhodes, a student of Andrew Cole at Columbia University. John Finn was on Rhodes's PhD committee, and was a co-author of a paper published in Physics of Plasmas on that work. A detailed study of nonlinear tearing modes in the presence of finite beta, favorable average curvature, and an error field has been nearly completed. In this study it was found that the response to an error field with real frequency a tearing layer agreed with analytic results, and that tearing mode saturation by flattening of the current profile as well as destabilization due to flattening of the pressure profile were important factors in the locking process. This work was presented at the 2017 APS/DPP Conference and at the 2018 Sherwood Conference. The effects of energetic ions on tearing modes was presented at the 2017 APS/DPP Conference and in the PhD dissertation of Michael Halfmoon. John Finn served as a committee member on Halfmoon's PhD committee. Extensions to this work including two fluid effects were performed and presented at the APS/DPP Conference in 2017.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Improved accuracy in degenerate variational integrators for guiding centre and magnetic field line flow

First-order-accurate degenerate variational integration (DVI) was introduced in Ellison et al. ( Phys. Plasmas , vol. 25, 2018, 052502) for systems with a degenerate Lagrangian, i.e. one in which the velocity-space Hessian is singular. In this paper we introduce second-order-accurate DVI schemes, both with and without non-uniform time stepping. We show that it is not in general possible to construct a second-order scheme with a preserved two-form by composing a first-order scheme with its adjoint, and discuss the conditions under which such a composition is possible. We build two classes of second-order-accurate DVI schemes. We test these second-order schemes numerically on two systems having non-canonical variables, namely the magnetic field line and guiding centre systems. Variational integration for Hamiltonian systems with non-uniform time steps, in terms of an extended phase space Hamiltonian, is generalized to non-canonical variables. It is shown that preservation of proper degeneracy leads to single-step (one-step) methods without parasitic modes, i.e. to non-uniform time step DVIs. This extension applies to second-order-accurate as well as first-order schemes, and can be applied to adapt the time stepping to an error estimate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Machine learning methods for probabilistic locked-mode predictors in tokamak plasmas

A rotating tokamak plasma can interact resonantly with the external helical magnetic perturbations, also known as error fields. This can lead to locking and then to disruptions. We leverage machine learning (ML) methods to predict the locking events. We use a coupled third-order nonlinear ordinary differential equation model to represent the interaction of the magnetic perturbation and the plasma rotation with the error field. This model is sufficient to describe qualitatively the locking and unlocking bifurcations. Here, we explore using ML algorithms with the simulation data and experimental data, focusing on the methods that can be used with sparse datasets. These methods lead to the possibility of the avoidance of locking in real-time operations. We describe the operational space in terms of two control parameters: the magnitude of the error field and the rotation frequency associated with the momentum source that maintains the plasma rotation. The outcomes are quan- tified by order parameters that completely characterize the state, whether locked or unlocked. We use unsupervised ML methods to classify locked/unlocked states and note the usefulness of a certain normalization of the order parameters. Three supervised ML classifiers are used in suite to estimate the probability of locking in the region of control parameter space with hysteresis, i.e., the set of control parameters for which both locked and unlocked states can exist. The results show that a neural network gives the best estimate of the locking probability. An analogy of the present locking model with the van der Waals equation of state is also provided.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Measurement of the Generalized Spin Polarizabilities of the Neutron in the Low-Q{sup 2} Region

Understanding the nucleon spin structure in the regime where the strong interaction becomes truly strong poses a challenge to both experiment and theory. At energy scales below the nucleon mass of about 1 GeV, the intense interaction among the quarks and gluons inside the nucleon makes them highly correlated. Their coherent behaviour causes the emergence of effective degrees of freedom, requiring the application of non-perturbative techniques such as chiral effective field theory1. Here we present measurements of the neutron’s generalized spin polarizabilities that quantify the neutron’s spin precession under electromagnetic fields at very low energy-momentum transfer squared down to 0.035 GeV2. In this regime, chiral effective field theory calculations2–4 are expected to be applicable. Our data, however, show a strong discrepancy with these predictions, presenting a challenge to the current description of the neutron’s spin properties.

Sulkosky, Vincent↗
DOE OSTI2020

Stability and Control of Burning Tokamak Plasmas with Resistive Walls (Second Year Report/ Final Report)

This research will focus on quantitative prediction of the stability, control, and equilibrium state evolution in toroidal burning plasmas. The stability of long pulse burning plasmas is highly sensitive to the physics of resonant layers in the plasma, sources of momentum and ow, kinetic effects of energetic particles, and boundary conditions at the wall, including feedback control and error fields. In ITER in particular, the low toroidal ow equilibrium state, sustained primarily by energetic alpha particles from fusion reactions, will require the consideration of all of these key elements to predict quantitatively the stability and evolution. The principal investigators on this proposal are leading experts in the relevant theoretical and computational areas, and aim to perform computations guided by analytic modeling, to address this physics in realistic configurations. The overall goal is to understand the key physics mechanisms that describe resistive toroidal burning plasmas, surrounded by a resistive wall, under active feedback control. With the physics of the energetic ions, resonant layers, resistive wall, and toroidal momentum transport included, this study will extend from recent publications in theory and simulation of individual effects and move toward predictive modeling for burning plasmas.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Nonlinear error field response in the presence of plasma rotation and real frequencies due to favorable curvature

This works describes the interaction of an phenomenon known to arise in tokamak fusion experiments with the magnetic perturbations that can arise at the edge of such experiments. We present nonlinear NIMROD resistive MHD simulations of the response of a rotating plasma to an error field when the plasma has weakly damped linear tearing modes (TMs), stabilized by a pressure gradient and favorable curvature. The favorable curvature leads to the Glasser effect: the occurrence of real frequencies and stabilization with positive stability index Δ'.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Final Technical Report Under SC0018278

Work supported by the grant "Stability and Control of Burning Tokamak Plasmas with Resistive Walls" in the period Aug 2017 to Aug 2018 can be cast in three general categories: (1) linear tearing regimes and the effects on the response to error fields and resistive walls: (2) nonlinear tearing modes in the presence of error fields, with locking; and (3) energetic particle effects. In the first category, we have completed a study of tearing mode regimes with extended MHD effects, including parallel dynamics and toroidal effects, notably favorable average curvature. This work has been presented at the 2018 Sherwood Conference in Ref.[1]. This work has been published in the ArXiv in Ref.[2] and is being prepared for publication in Physics of Plasmas. A systematic study of shaping in a tokamak plasma with resistive tearing layers and a resistive wall has been presented in the Ph D dissertation[3] by Dov Rhodes, a student of Andrew Cole at Columbia University. John Finn was on Rhodes's Ph D committee, and was a co-author of a paper published in Physics of Plasmas on that work. A detailed study of nonlinear tearing modes in the presence of finite beta, favorable average curvature, and an error field has been nearly completed. In this study it was found that the response to an error field with real frequency a tearing layer agreed with analytic results, and that tearing mode saturation by flattening of the current profile as well as destabilization due to flattening of the pressure profile were important factors in the locking process. This work was presented at the 2017 APS/DPP Conference and at the 2018 Sherwood Conference. The effects of energetic ions on tearing modes was presented at the 2017 APS/DPP Conference and in the Ph D dissertation of Michael Halfmoon. John Finn served as a committee member on Halfmoon's Ph D committee. Extensions to this work including two-fluid effects were performed and presented at the APS/DPP Conference in 2017.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Properties of the electrostatically driven helical plasma state

A novel plasma state has been found [Akçay et al ., Phys. Plasmas 24 , 052503 (2017)] in the presence of a uniform applied axial magnetic field in periodic cylindrical geometry. This state is driven by external electrostatic fields provided by helical electrodes with a ($m=1,n=1$) (helical) symmetry where $m$ and $n$ represent the poloidal and axial harmonics. The resulting plasma is a function of the cylinder radius $r < r_w$ and helical angle $u=mθ-nζ$, where θ is the poloidal angle, and $ζ=z/R$ is the normalized axial coordinate in the context of a periodic cylinder. In this reference, the strongly driven form of the state was found to have a strong axial mean current density, with a mean-field line safety factor $q_0(r)$ just above the pitch of the electrodes m/n=1 in the interior, where the plasma is nearly force-free. However, at the edge the current density has a component perpendicular to the magnetic field B. This perpendicular current density drives nearly Alfvénic helical plasma flows, a notable feature of these states. This state is being studied for its possible application in $DC$ electrical transformers. We present results on several issues of importance for these applications: the transient leading to the steady state; the twist and writhe of the field lines and their relation with the current density; the properties of the current density streamlines and length of the current density lines connected to the electrodes; the sensitivity to changes in the velocity boundary conditions; the effect of varying the radial resistivity profile; and the effects of a concentrated electrode potential

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Geometric scalings for the electrostatically driven helical plasma state

A new plasma state has been investigated [Akcay et al., Phys. Plasmas 24, 052503 (2017)], with a uniform applied axial magnetic field in a periodic cylinder of length $L=2πR$, driven by helical electrodes. The drive is single helicity, depending on $mθ+kz=mθ-nζ$, where $ζ=z/R$ and $k=-n/R$. For strong $(m,n)=(1,1)$ drive, the state was found to have a strong axial mean current density, with a mean-field safety factor $q_0(r)$ just above the pitch of the electrodes $m/n=1$ in the interior. This state has possible applications to DC electrical transformers and tailoring of the current profile in tokamaks. Here, we study two geometric issues of interest for these applications: (i) scaling of properties with the plasma length or aspect ratio and (ii) behavior for different helicities, specifically $(m,n)=(1,n)$ for n > and $(m,n)=(2,1)$.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Electrostatically driven helical plasma state

A novel plasma state has been found in the presence of a uniformly applied axial magnetic field in periodic cylindrical geometry. This state is driven electrostatically by helical electrodes, providing a driving field that depends on the radius and mθ−nζ, where θ is the poloidal angle and ζ=z/R is the toroidal angle. We focus on m=n=1. The radial magnetic field at the wall is taken to be zero. With weak driving, the resulting distortion is very small, but for stronger driving, the mean field of the state has field line safety factor q0(r) just above the pitch of the electrodes m∕n = 1 except near the edge, where q0 increases monotonically. This state is characterized as a single helicity Ohmic equilibrium with the helical symmetry of the applied field. The plasma appears to be close to force-free in the interior, but current density crosses the magnetic flux surfaces near the edge, where current must enter and exit through the helical electrodes. This perpendicular current density drives large helical plasma flows. The sensitivity of this state to flow boundary conditions, plasma resistivity profile, the strength of electrostatic driving, and parameters such as the loop voltage and the Lundquist number is explored. The magnetic helicity is calculated for both the transient period and time-asymptotic state. Possible applications to current drive in toroidal confinement devices and to electrical transformers are discussed.

Akçay, Cihan↗
DOE OSTI2015

Error field penetration and locking to the backward propagating wave

In this letter we investigate error field penetration, or locking, behavior in plasmas having stable tearing modes with finite real frequencies w r in the plasma frame. In particular, we address the fact that locking can drive a significant equilibrium flow. We show that this occurs at a velocity slightly above v = w r /k, corresponding to the interaction with a backward propagating tearing mode in the plasma frame. Results are discussed for a few typical tearing mode regimes, including a new derivation showing that the existence of real frequencies occurs for viscoresistive tearing modes, in an analysis including the effects of pressure gradient, curvature and parallel dynamics. The general result of locking to a finite velocity flow is applicable to a wide range of tearing mode regimes, indeed any regime where real frequencies occur.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Issues in measure-preserving three dimensional flow integrators: Self-adjointness, reversibility, and non-uniform time stepping

Properties of integration schemes for solenoidal fields in three dimensions are studied, with a focus on integrating magnetic field lines in a plasma using adaptive time stepping. It is shown that implicit midpoint (IM) and a scheme we call three-dimensional leapfrog (LF) can do a good job (in the sense of preserving KAM tori) of integrating fields that are reversible, or (for LF) have a “special divergence-free” (SDF) property. We review the notion of a self-adjoint scheme, showing that such schemes are at least second order accurate and can always be formed by composing an arbitrary scheme with its adjoint. We also review the concept of reversibility, showing that a reversible but not exactly volume-preserving scheme can lead to a fractal invariant measure in a chaotic region, although this property may not often be observable. We also show numerical results indicating that the IM and LF schemes can fail to preserve KAM tori when the reversibility property (and the SDF property for LF) of the field is broken. We discuss extensions to measure preserving flows, the integration of magnetic field lines in a plasma and the integration of rays for several plasma waves. The main new result of this paper relates to non-uniform time stepping for volume-preserving flows. We investigate two potential schemes, both based on the general method of Feng and Shang [Numer. Math. 71, 451 (1995)], in which the flow is integrated in split time steps, each Hamiltonian in two dimensions. The first scheme is an extension of the method of extended phase space, a well-proven method of symplectic integration with non-uniform time steps. This method is found not to work, and an explanation is given. The second method investigated is a method based on transformation to canonical variables for the two split-step Hamiltonian systems. This method, which is related to the method of non-canonical generating functions of Richardson and Finn [Plasma Phys. Controlled Fusion 54, 014004 (2012)], appears to work very well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Variational principles with Padé approximants for tearing mode analysis

Tearing modes occur in several distinct physical regimes, and it is often important to compute the inner layer response for these modes with various effects. There is a need for an approximate and efficient method of solving the inner layer equations in all these regimes. In this paper, we introduce a method of solving the inner layer equations based on using a variational principle with Padé approximants. For all the regimes considered, the main layer equations to be solved are inhomogeneous, and Padé approximants give a convenient and efficient method of satisfying the correct asymptotic behavior at the edge of the layer. Results using this variational principle—Padé approximant method in three of these regimes is presented. These regimes are the constant-ψ resistive-inertial (RI) regime, the constant-ψ viscoresistive regime, and the non-constant-ψ inviscid tearing regime. The last regime includes the constant-ψ RI regime and the inertial regime. The results show that reasonable accuracy can be obtained very efficiently with Padé approximants having a small number of parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Three-dimensional force-free looplike magnetohydrodynamic equilibria

Computations of three-dimensional force-free magnetohydrodynamic (MHD) equilibria, del x B = lambdaB with lambda = lambda(sub 0), a constant are presented. These equilibria are determined by boundary conditions on a surface corresponding to the solar photosphere. The specific boundary conditions used correspond to looplike magnetic fields in the corona. It is found that as lambda(sub 0) is increased, the loops of flux become kinked, and for sufficiently large lambda(sub 0), develop knots. The relationship between the kinking and knotting properties of these equilibria and the presence of a kink instability and related loss of equilibrium is explored. Clearly, magnetic reconnection must be involved for an unknotted loop equilibrium to become knotted, and speculations are made about the creation of a closed hyperbolic field line (X-line) about which this reconnection creating knotted field lines is centered.

Finn, John M.↗
NASA NTRS1993

Long-term containment of energetic particles in coronal loops

Recent observation from the Compton Gamma-Ray Observatory shows that gamma-ray emission after a solar flare can last for as long as 8 hours. There is also evidence that electrons and protons are accelerated only during the impulsive phase of the flare and are subsequently mirror trapped in coronal magnetic loops. This poses the following dilemma: if the magnetic field lines in the loop are simple plane arches, the protons will drift across the cross section of the loop in seconds to minutes, rather than hours. To solve the dilemma, we use guiding center theory to show that long-term containment of energetic protons in a coronal loop is possible if magnetic field lines have enough twist. We also find that in the trapped region of the loop, the twist angle of field lines between the mirror points of a bounce orbit is approximately 2 pi.

Lau, Yun-Tung↗
NASA NTRS1993

Loss of equilibrium and reconnection in tearing of two-dimensional equilibria

Two-dimensional tearinglike behavior is studied in reduced resistive magnetohydrodynamics (MHD) with flux conserving boundary conditions on a rectangular box. The tearinglike perturbations do not destroy the symmetries of the initial state, either discrete or continuous. In such cases, in which the perturbation does not break a symmetry of the equilibrium, linear instability is typically not directly observed. However, there can be a loss of equilibrium associated with the existence of a tearing unstable state. These ideas are illustrated with three examples: a very elongated tokamak, a tokamak with pinching coils to elongate its flux surfaces, and a model for the magnetotail or for solar arcades. The loss of equilibrium is demonstrated by means of a nonlinear energy functional. The importance of the fact that the dynamics shows a loss of equilibrium is that a large amount of free energy can be released, in the form of reconnection, and that there is a possibility of hysteresis.

Finn, John M.↗