A Riccati Solution for the Ideal MHD Plasma Response with Applications to Real-time Stability Control

Glasser, A.; Kolemen, E.; Glasser, A.H.
Issue date: March 2018
Cite as:
Glasser, A., Kolemen, E., & Glasser, A.H. (2017). A Riccati Solution for the Ideal MHD Plasma Response with Applications to Real-time Stability Control [Data set]. Princeton Plasma Physics Laboratory, Princeton University.
@electronic{glasser_a_2017,
  author      = {Glasser, A. and
                Kolemen, E. and
                Glasser, A.H.},
  title       = {{A Riccati Solution for the Ideal MHD Pla
                sma Response with Applications to Real-t
                ime Stability Control}},
  publisher   = {{Princeton Plasma Physics Laboratory, Pri
                nceton University}},
  year        = 2017
}
Abstract:

To effectuate near real-time feedback control of ideal MHD instabilities in a tokamak geometry, a rapid solution for stability analysis is a prerequisite. Toward this end, we reformulate the δW stability method with a Hamilton-Jacobi theory, elucidating analytical and numerical features of the generic tokamak ideal MHD stability problem. The plasma response matrix is demonstrated to be the solution of an ideal MHD matrix Riccati differential equation (MRDE). Since Riccati equations are prevalent in the control theory literature, such a shift in perspective brings to bear a range of numerical methods that are well-suited to the robust, fast solution of control problems. We discuss the usefulness of Riccati techniques in solving the stiff ODEs often encountered in ideal MHD stability analyses-—for example, in tokamak edge and stellarator physics. We then demonstrate the applicability of such methods to an existing 2D ideal MHD stability code—DCON—enabling its parallel operation in near real-time. Output is shown to match with high accuracy, and with wall-clock time ≪ 1s. Such speed may help enable active feedback ideal MHD stability control, especially in tokamak plasmas whose ideal MHD equilibria evolve with inductive timescale τ > 1s-—as in ITER.

Show More
Filename Size