SIESTA: A scalable iterative equilibrium solver for toroidal applications

Hirshman, S. P.; Sanchez, R.; Cook, C. R.
June 2011
Physics of Plasmas;Jun2011, Vol. 18 Issue 6, p062504
Academic Journal
A new solver for rapidly obtaining magnetohydrodynamic (MHD) equilibria in toroidal systems in the presence of islands and stochastic regions is described. It is based on the Kulsrud-Kruskal MHD energy minimization principle. To carry out this minimization, small displacements are made around a convenient set of curvilinear coordinates obtained from a nearby three-dimensional equilibrium that assumes nested surfaces. Because the changes of the magnetic fields and pressure are small, corresponding to small changes in the initial magnetic and kinetic energies, solutions for the linearized perturbations can be used to rapidly and iteratively find lower energy states with magnetic islands. A physics-based preconditioner is developed to accelerate the convergence of the iterative procedure to obtain an ideal MHD equilibrium with broken magnetic surfaces (islands).


Related Articles

  • Peturbation: An Efficient Technique for the Solution of Very Large Instances of the Euclidean TSP. Codenotti, Bruno; Manzini, Giovanni; Margara, Luciano // INFORMS Journal on Computing;Spring96, Vol. 8 Issue 2, p125 

    Presents a study that introduced a technique for developing efficient iterated local search procedures and applied it to solve instances of the Euclidean Traveling Salesman Problem (TSP). Description of the Euclidean TSP; Iterating local search and the method of perturbation; Discussion of...

  • Constructing fast-growing perturbations for the nonlinear regime. Barkmeijer, Jan // Journal of the Atmospheric Sciences;10/1/96, Vol. 53 Issue 19, p2838 

    Presents an iterative method to determine perturbations that maintain their growth rate in the medium-range period. Model and time evolution of singular vectors; Modification of perturbations; Numerical experiments.

  • Plastic Collapse Load of Crown-Hinged Steel Circular Arches: A Theoretical Method. Spoorenberg, R. C.; Snijder, H. H.; Hoenderkamp, J. C. D. // Advances in Structural Engineering;Apr2013, Vol. 16 Issue 4, p721 

    For construction purposes and to avoid detrimental influences of foundation settlements arches are not always made from a single arch-rib but are built by connecting two curvilinear elements at the crown with a hinge. These arches are also known as crown-hinged arches. This paper presents an...

  • On the continua in two-dimensional nonadiabatic magnetohydrodynamic spectra. De Ploey, A.; Van der Linden, R. A. M.; Belie¨n, A. J. C. // Physics of Plasmas;Jun2000, Vol. 7 Issue 6 

    The equations for the continuous subspectra of the linear magnetohydrodynamic (MHD) normal modes spectrum of two-dimensional (2D) plasmas are derived in general curvilinear coordinates, taking nonadiabatic effects in the energy equation into account. Previously published derivations of...

  • On the perturbation of multiconfiguration wave functions. Rolik, Zoltán; Szabados, Ágnes; Surján, Péter R. // Journal of Chemical Physics;7/22/2003, Vol. 119 Issue 4, p1922 

    A simple variant of perturbation theory is used to correct reference states of a general multiconfigurational character. The full solution of an active space is not required, and no iterative procedure is applied to construct the resolvent operator. The perturbed wave function is expanded in a...

  • RANDOM FRACTALS AND PROBABILITY METRICS. Hutchinson, John E.; Ruschendorf, Ludger // Advances in Applied Probability;Dec2000, Vol. 32 Issue 4, p925 

    Introduces metrics in the space of random measures to prove the existence and uniqueness results for self-similar random fractal measures. Random measure; Self-similar; Probability metric; Random fractal; Monge-Kantorovich metric; Scaling operator; Iterated function system; Minimal metric.

  • Iterates of expanding maps. Barbour, A. D.; Gerrard, R. M.; Reinert, G. // Probability Theory & Related Fields;2000, Vol. 116 Issue 2, p151 

    Abstract. The iterates of expanding maps of the unit interval into itself have many of the properties of a more conventional stochastic process, when the expanding map satisfies some regularity conditions and when the starting point is suitably chosen at random. In this paper, we show that the...

  • Simulation of Laminar Natural Convection in a Cavity with Cylindrical Obstacles. Mushatet, Khudheyer S. // Australian Journal of Basic & Applied Sciences;2011, Vol. 5 Issue 6, p636 

    The laminar natural convection inside a rectangular cavity containing two cylindrical obstacles has been numerically investigated. A curvilinear coordinates system was used to transfer the physical space into a computational domain. The governing partial differential equations are solved using...

  • Return to equilibrium. Bach, Volker; Fro¨hlich, Ju¨rg; Sigal, Israel Michael // Journal of Mathematical Physics;Jun2000, Vol. 41 Issue 6 

    We study an atom with finitely many energy levels in contact with a heat bath consisting of photons (blackbody radiation) at a temperature T>0. The dynamics of this system is described by a Liouville operator, or thermal Hamiltonian, which is the sum of an atomic Liouville operator, of a...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics