Three-Dimensional Stellarator Codes

Total Page:16

File Type:pdf, Size:1020Kb

Three-Dimensional Stellarator Codes Three-dimensional stellarator codes P. R. Garabedian* Courant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New York, NY 10012 Contributed by P. R. Garabedian, June 3, 2002 Three-dimensional computer codes have been used to develop anomalous transport can be modeled by iterating on small quasisymmetric stellarators with modular coils that are promising variations of the electric potential within the magnetic surfaces candidates for a magnetic fusion reactor. The mathematics of to achieve quasineutrality between the distributions of ions and plasma confinement raises serious questions about the numerical electrons. The method simulates complicated transport in stel- calculations. Convergence studies have been performed to assess larators remarkably well, and numerical results have been ob- the best configurations. Comparisons with recent data from large tained for the large helical device (LHD) experiment at the stellarator experiments serve to validate the theory. National Institute for Fusion Studies in Toki, Japan that are in excellent agreement with recent observations of the energy odular stellarators can be viewed as an advanced tokamak confinement time at high temperatures (6, 7). A vectorized Mhybrid appropriate for implementation as a fusion reactor version of the TRAN code runs efficiently on standard work (1). Quasisymmetry of the magnetic spectrum is predicted to stations. give good confinement at high temperatures, and adequate Convergence Studies rotational transform from the external magnetic field is expected to stabilize the plasma. New configurations have been designed Calculation of toroidal equilibrium of a plasma without two- by making imaginative use of three-dimensional computer codes. dimensional symmetry is a problem in mathematics that is not Because the mathematics of these stellarators is complicated, we well posed. In terms of the toroidal flux s and a pair of angular have performed convergence studies applicable to proof of flux coordinates ␪ and ␾, the parallel current has a Fourier principle experiments that are being planned. expansion of the form We are primarily concerned with the NSTAB equilibrium and J⅐B m Ϫ In ͑m Ϫ In͒␪ ϩ ͑n Ϫ ␫m͒␾ stability code and the TRAN Monte Carlo transport code ϭ pЈ ͸ B cosͫ ͬ , developed at New York University by Octavio Betancourt and B2 n Ϫ ␫m mn 1 Ϫ ␫I Mark Taylor (2–5). The codes are applied to a compact stel- larator called the MHH2 that has two field periods and excellent where Bmn is the magnetic spectrum of Fourier coefficients of ͞ 2 ␫ quasiaxial symmetry. For the calculations we selected a config- 1 B , is the rotational transform, and I is the net current (8). uration with realistic physical parameters that provide good A configuration is called quasisymmetric if a single row, column, convergence, enabling us to perform long runs and make or diagonal of the double array Bmn dominates that spectrum. Ϫ␫ estimates of numerical errors. Theoretical conclusions can be The small denominators n m vanish at resonant surfaces ␫ drawn that are relevant to a wider range of examples, such as the where is rational, so smooth solutions of the partial differential optimized stellarator specified in Table 1. equations describing MHD equilibrium do not in general exist in The NSTAB code is a computer implementation of the three dimensions. In numerical work this leads us to construct variational principle of ideal magnetohydrodynamics (MHD). If weak, discontinuous solutions of discrete equations that are B is the magnetic field and p is the pressure, solutions of the expressed in conservation form, but even the best methods only magnetostatics equations are found by minimizing the potential converge in an asymptotic sense. Enough spectral terms must be energy included to eliminate significant truncation error, but not so many that the results become meaningless. Our convergence studies for the NSTAB code clarify how this can be accom- E ϭ ͵͵͵͓B2͞2 Ϫ p͔dV plished in practice. For convergence studies we have selected an example of the MHH2 compact stellarator, shown in Fig. 1, that has two field in a coordinate system compatible with toroidal geometry in periods, plasma aspect ratio three, excellent quasiaxial symme- three dimensions. An accurate finite difference scheme is used try, and a limit near 5% for the average value of ␤ ϭ 2p͞B2. The in the radial direction, and dependence on the poloidal and separatrix has a smooth shape that facilitates making long runs SCIENCES toroidal angles is handled by the spectral method. It is assumed of the NSTAB code. The coordinate system for the computa- that there are nested toroidal flux surfaces, and the differential tions is chosen to rotate once in the poloidal direction over a full APPLIED PHYSICAL equations are written in a conservation form that captures circuit of the device in the toroidal direction. Zoning of the islands and current sheets. The resolution is so good that poloidal angle has been adjusted to provide good resolution on questions of stability can be settled by a mountain pass theorem crude grids and reasonable spacing of the mesh both at the asserting that when more than one solution of the problem can separatrix and the magnetic axis, which is determined by a robust be found then an unstable equilibrium must exist corresponding algorithm. However, numerical errors appear at the magnetic to a saddle point in the energy landscape. Bifurcated equilibria axis if the radial mesh is too fine. are calculated whose magnetic surfaces have Poincare´sections Applying an accelerated method of steepest descent to the displaying the structure of the most unstable modes. variational principle, we test for stability by examining runs of The TRAN code uses a split time algorithm to calculate the the NSTAB code in which some dangerous mode has been confinement time of test particles by alternately tracking guiding triggered by introducing temporarily an appropriate forcing center orbits and applying a random walk that represents colli- term. A run predicts stability if the mode decays during further sions. The magnetic field B and the flow field U of the plasma in a background obtained by using NSTAB are held fixed during iterations that impose quasineutrality. Conservation of momen- Abbreviations: LHD, large helical device; MHD, magnetohydrodynamics; W7-AS, Wendel- tum might be enforced by the selection of U, but because it is not stein 7-AS. thermal U has little effect on the collision operator. However, *E-mail: [email protected]. www.pnas.org͞cgi͞doi͞10.1073͞pnas.162330399 PNAS ͉ August 6, 2002 ͉ vol. 99 ͉ no. 16 ͉ 10257–10259 Downloaded by guest on September 27, 2021 Table 1. Fourier coefficients ⌬mn of an MHH2 stellarator defined Table 2. Relationship between the highest degree N of the spectral terms used in a run of NSTAB and the resulting critical ؍ in cylindrical and toroidal coordinates by the formula r ؉ iz iu ؊imu؉2inv e ¥ ⌬mn e value of ␤ that is calculated n Ϫ10123N 16 20 24 32 48گm ␤ 0.060 0.050 0.045 0.040 0.039 Ϫ1 0.190 0.130 Ϫ0.015 0.000 0.000 0 0.000 1.000 0.000 0.000 0.000 1 0.150 3.000 0.250 0.050 0.000 2 0.000 Ϫ0.090 Ϫ0.420 Ϫ0.070 0.000 residuals so that the test of stability would become more reliable. 3 0.000 0.000 Ϫ0.040 0.080 0.000 Runs were continued as long as possible to achieve accuracy 4 0.000 0.015 0.000 Ϫ0.015 Ϫ0.015 sufficient for a convergence study. Recent advances in computer technology have enabled us to do this economically. In Table 2 we compare the degree N of the spectral calcula- tions with the corresponding estimate of an average ␤ limit based iterations, but if it grows then the equilibrium is unstable. A more on the mountain pass theorem. If the degree is too low the convincing conclusion can be drawn from the mountain pass method does not provide meaningful results about stability. theorem if the iterations converge to a bifurcated solution whose However, for n ϭ 24 one obtains efficiently answers that are of stellarator symmetry is visibly broken by magnetic surfaces that sufficient accuracy to optimize the design of a quasisymmetric exhibit the structure of the dangerous mode. The numerical stellarator like the MHH2. The results do not change signifi- results depend on the maximum degree N of factors in each of cantly for N Ͼ 32, which is a good value to choose from the point two angular coordinates that specify the spectral terms included of view of asymptotic convergence. At N ϭ 48 a stage is reached in the computation. The purpose of our convergence study is to where the numerical method may fail in long runs because of the decide whether the prediction about stability approaches a singular behavior of the solution. meaningful limit as N increases. Table 3 shows how the convergence of the iterative scheme For the example of the MHH2 we have performed equilibrium used in the NSTAB code depends on the degree N of the spectral and stability runs of the NSTAB code with between 14 and 28 terms that are used. For crude grids it is easy to reduce the mesh intervals in the radial flux coordinate s and with a largest residuals to the level of round-off error in the computer. degree N of the spectral terms ranging as high as 48. Because of However, as N increases the effectiveness of the scheme dete- the nested surface hypothesis, islands of small width are captured riorates, and the method may only converge in the asymptotic better on crude grids, and the calculations are relatively insen- sense that at first the errors become smaller, but later they sitive to the radial mesh size because the finite difference scheme increase without limit.
Recommended publications
  • Magnetohydrodynamics - Overview
    Center for Turbulence Research 361 Proceedings of the Summer Program 2006 Magnetohydrodynamics - overview Magnetohydrodynamics (MHD) is the study of the interaction of electrically conduct- ing fluids and electromagnetic forces. MHD problems arise in a wide variety of situations ranging from the explanation of the origin of Earth's magnetic field and the prediction of space weather to the damping of turbulent fluctuations in semiconductor melts dur- ing crystal growth and even the measurement of the flow rates of beverages in the food industry. The description of MHD flows involves both the equations of fluid dynamics - the Navier-Stokes equations - and the equations of electrodynamics - Maxwell's equa- tions - which are mutually coupled through the Lorentz force and Ohm's law for moving electrical conductors. Given the vast range of Reynolds numbers and magnetic Reynolds numbers characterising astrophysical, geophysical, and industrial MHD problems it is not feasible to devise an all-encompassing numerical code for the simulation of all MHD flows. The challenge in MHD is rather to combine a judicious application of analytical tools for the derivation of simplified equations with accurate numerical algorithms for the solution of the resulting mathematical models. The 2006 summer program provided a unique opportunity for scientists and engineers working in different sub-disciplines of MHD to interact with each other and benefit from the numerical capabilities available at the Center for Turbulence Research. In materials processing such as the growth of semiconductor crystals or the continuous casting of steel turbulent flows are often undesirable. The application of steady magnetic fields provides a convenient means for the contactless damping of turbulent fluctuations and for the laminarisation of flows.
    [Show full text]
  • STATUS of FUSION ENERGY Impact & Opportunity for Alberta Volume II
    STATUS OF FUSION ENERGY Impact & Opportunity for Alberta Volume II Appendices Prepared by Alberta/Canada Fusion Energy Program March 2014 ALBERTA COUNCIL OF TECHNOLOGIES Gratefully acknowledges the support of: Alberta Energy Stantec Corporation University of Alberta Alberta/Canada Fusion Energy Advisory Committee Gary Albach Nathan Armstrong Brian Baudais Will Bridge Robert Fedosejevs Peter Hackett Chris Holly Jerry Keller Brian Kryska Axel Meisen Rob Pitcairn Klaas Rodenburg John Rose Glenn Stowkowy Martin Truksa Gary Woloshyniuk Perry Kinkaide Allan Offenberger A special thank you is extended to the institutions (identified in this report) that were visited and to the many persons who so graciously hosted our site visits, provided the briefing material presented in this status report and thereby assisted our fusion assessment. Report Authors Allan Offenberger Robert Fedosejevs Klaas Rodenburg Perry Kinkaide Contact: Dr. Perry Kinkaide [email protected] 780-990-5874 Dr. Allan Offenberger [email protected] 780-483-1740 i TABLE OF CONTENTS Page List of Acronyms ………………………………………………………………………….. iii List of Figures……………………………………………………………………………… iv Appendix A: Assessment of Major Global Fusion Technologies 1.0 Context - Global Energy Demand……………………………………………………… 1 1.0.1 Foreward ……………………………………………………………………… 1 1.0.2 Energy Trends………………………………………………………………… 2 1.0.3 Energy From Fusion Reactions……………………………………………… 4 1.1 Major Approaches to Fusion Energy………………………………………………….. 7 1.1.1 Introduction……………………………………………………………………. 7 1.1.2 Fusion Reactions & the Fuel Cycle………………………………………….. 8 1.1.3 IFE Approaches to Fusion…………………………………………………… 11 1.1.3.1 Introduction………………………………………………………….. 11 1.1.3.2 Indirect Drive…………………………………………………………14 1.1.3.3 Direct Drive…………………………………………………………. 16 1.1.3.4 Fast Ignition………………………………………………………… 17 1.1.3.5 Shock Ignition………………………………………………………..19 1.1.3.6 IFE Power Reactor Systems……………………………………….20 1.1.3.7 Modeling Codes…………………………………………………….
    [Show full text]
  • A Mathematical Introduction to Magnetohydrodynamics
    Typeset in LATEX 2" July 26, 2017 A Mathematical Introduction to Magnetohydrodynamics Omar Maj Max Planck Institute for Plasma Physics, D-85748 Garching, Germany. e-mail: [email protected] 1 Contents Preamble 3 1 Basic elements of fluid dynamics 4 1.1 Kinematics of fluids. .4 1.2 Lagrangian trajectories and flow of a vector field. .5 1.3 Deformation tensor and vorticity. 14 1.4 Advective derivative and Reynolds transport theorem. 17 1.5 Dynamics of fluids. 20 1.6 Relation to kinetic theory and closure. 24 1.7 Incompressible flows . 32 1.8 Equations of state, isentropic flows and vorticity. 34 1.9 Effects of Euler-type nonlinearities. 35 2 Basic elements of classical electrodynamics 39 2.1 Maxwell's equations. 39 2.2 Lorentz force and motion of an electrically charged particle. 52 2.3 Basic mathematical results for electrodynamics. 57 3 From multi-fluid models to magnetohydrodynamics 69 3.1 A model for multiple electrically charged fluids. 69 3.2 Quasi-neutral limit. 74 3.3 From multi-fluid to a single-fluid model. 82 3.4 The Ohm's law for an electron-ion plasma. 86 3.5 The equations of magnetohydrodynamics. 90 4 Conservation laws in magnetohydrodynamics 95 4.1 Global conservation laws in resistive MHD. 95 4.2 Global conservation laws in ideal MHD. 98 4.3 Frozen-in law. 100 4.4 Flux conservation. 105 4.5 Topology of the magnetic field. 109 4.6 Analogy with the vorticity of isentropic flows. 117 5 Basic processes in magnetohydrodynamics 119 5.1 Linear MHD waves.
    [Show full text]
  • Stability Study of the Cylindrical Tokamak--Thomas Scaffidi(2011)
    Ecole´ normale sup erieure´ Princeton Plasma Physics Laboratory Stage long de recherche, FIP M1 Second semestre 2010-2011 Stability study of the cylindrical tokamak Etude´ de stabilit´edu tokamak cylindrique Author: Supervisor: Thomas Scaffidi Prof. Stephen C. Jardin Abstract Une des instabilit´es les plus probl´ematiques dans les plasmas de tokamak est appel´ee tearing mode . Elle est g´en´er´ee par les gradients de courant et de pression et implique une reconfiguration du champ magn´etique et du champ de vitesse localis´ee dans une fine r´egion autour d’une surface magn´etique r´esonante. Alors que les lignes de champ magn´etique sont `al’´equilibre situ´ees sur des surfaces toriques concentriques, l’instabilit´econduit `ala formation d’ˆıles magn´etiques dans lesquelles les lignes de champ passent d’un tube de flux `al’autre, rendant possible un trans- port thermique radial important et donc cr´eant une perte de confinement. Pour qu’il puisse y avoir une reconfiguration du champ magn´etique, il faut inclure la r´esistivit´edu plasma dans le mod`ele, et nous r´esolvons donc les ´equations de la magn´etohydrodynamique (MHD) r´esistive. On s’int´eresse `ala stabilit´ede configurations d’´equilibre vis-`a-vis de ces instabilit´es dans un syst`eme `ala g´eom´etrie simplifi´ee appel´ele tokamak cylindrique. L’´etude est `ala fois analytique et num´erique. La solution analytique est r´ealis´ee par une m´ethode de type “couche limite” qui tire profit de l’´etroitesse de la zone o`ula reconfiguration a lieu.
    [Show full text]
  • Magnetohydrodynamics 1 19.1Overview
    Contents 19 Magnetohydrodynamics 1 19.1Overview...................................... 1 19.2 BasicEquationsofMHD . 2 19.2.1 Maxwell’s Equations in the MHD Approximation . ..... 4 19.2.2 Momentum and Energy Conservation . .. 8 19.2.3 BoundaryConditions. 10 19.2.4 Magneticfieldandvorticity . .. 12 19.3 MagnetostaticEquilibria . ..... 13 19.3.1 Controlled thermonuclear fusion . ..... 13 19.3.2 Z-Pinch .................................. 15 19.3.3 Θ-Pinch .................................. 17 19.3.4 Tokamak.................................. 17 19.4 HydromagneticFlows. .. 18 19.5 Stability of Hydromagnetic Equilibria . ......... 22 19.5.1 LinearPerturbationTheory . .. 22 19.5.2 Z-Pinch: Sausage and Kink Instabilities . ...... 25 19.5.3 EnergyPrinciple ............................. 28 19.6 Dynamos and Reconnection of Magnetic Field Lines . ......... 29 19.6.1 Cowling’stheorem ............................ 30 19.6.2 Kinematicdynamos............................ 30 19.6.3 MagneticReconnection. 31 19.7 Magnetosonic Waves and the Scattering of Cosmic Rays . ......... 33 19.7.1 CosmicRays ............................... 33 19.7.2 Magnetosonic Dispersion Relation . ..... 34 19.7.3 ScatteringofCosmicRays . 36 0 Chapter 19 Magnetohydrodynamics Version 1219.1.K.pdf, 7 September 2012 Please send comments, suggestions, and errata via email to [email protected] or on paper to Kip Thorne, 350-17 Caltech, Pasadena CA 91125 Box 19.1 Reader’s Guide This chapter relies heavily on Chap. 13 and somewhat on the treatment of vorticity • transport in Sec. 14.2 Part VI, Plasma Physics (Chaps. 20-23) relies heavily on this chapter. • 19.1 Overview In preceding chapters, we have described the consequences of incorporating viscosity and thermal conductivity into the description of a fluid. We now turn to our final embellishment of fluid mechanics, in which the fluid is electrically conducting and moves in a magnetic field.
    [Show full text]
  • Application of High Performance Computing for Fusion Research Shimpei Futatani
    Application of high performance computing for fusion research Shimpei Futatani Department of Physics, Universitat Politècnica de Catalunya (UPC), Barcelona, Spain The study of the alternative energy source is one of the absolutely imperative research topics in the world as our present energy sources such as fossil fuels are limited. This work is dedicated to the nuclear fusion physics research in close collaboration with existing experimental fusion devices and the ITER organization (www.iter.org) which is an huge international nuclear fusion R&D project. The goal of the ITER project is to demonstrate a clean and safe energy production by nuclear fusion which is the reaction that powers the sun. One of the ideas of the nuclear fusion on the earth is that the very high temperature ionized particles, forming a plasma can be controlled by a magnetic field, called magnetically confined plasma. This is essential, because no material can be sustained against such high temperature reached in a fusion reactor. Tokamak is a device which uses a powerful magnetic field to confine a hot plasma in the shape of a torus. It is a demanding task to achieve a sufficiently good confinement in a tokamak for a ‘burning plasma’ due to various kinds of plasma instabilities. One of the critical unsolved problems is MHD (MagnetoHydroDynamics) instabilities in the fusion plasmas. The MHD instabilities at the plasma boundary damages the plasma facing component of the fusion reactor. One of techniques to control the MHD instabilities is injection of pellets (small deuterium ice bodies). The physics of the interaction between pellet ablation and MHD dynamics is very complex, and uncertainties still remain regarding the theoretical physics as well as the numerical modelling point of view.
    [Show full text]
  • Identification of the Sequence of Steps Intrinsic to Spheromak Formation
    Identification of the Sequence of Steps Intrinsic to Spheromak Formation P.M. Bellan, S. You, and G. S. Yun California Institute of Technology Pasadena, CA 91125, USA Abstract. A planar coaxial electrostatic helicity source is used for studying the relaxation process intrinsic to spheromak formation. Experimental observations reveal that spheromak formation involves: (1) breakdown and creation of a number of distinct, arched, filamentary, plasma-filled flux loops that span from cathode to anode gas nozzles, (2) merging of these loops to form a central column, (3) jet-like expansion of the central column, (4) kink instability of the central column, (5) conversion of toroidal flux to poloidal flux by the kink instability. Steps 1 and 3 indicate that spheromak formation involves an MHD pumping of plasma from the gas nozzles into the magnetic flux tube linking the nozzles. In order to measure this pumping, the gas puffing system has been modified to permit simultaneous injection of different gas species into the two ends of the flux tube linking the wall. Gated CCD cameras with narrow-band optical filters are used to track the pumped flows. Keywords: spheromak, magnetohydrodynamics, kink, jet PACS: 52.55.Ip, 52.35.Py, 52.55.Wq, 52.59.Dk INTRODUCTION which is rich in toroidal field energy (typical poloidal currents are 100 kA) but poor in poloidal field energy. Spheromaks are toroidal MHD equilibria incorporating Spheromak formation involves conversion of toroidal nested poloidal magnetic flux surfaces suitable for magnetic flux into poloidal magnetic flux while confining a fusion plasma. Unlike tokamaks, conserving magnetic helicity (but not conserving spheromaks are a self-organizing, minimum-energy magnetic energy).
    [Show full text]
  • Significance of MHD Effects in Stellarator Confinement
    Ms-313705 Final Paper Significance of MHD Effects in Stellarator Confinement A. Weller1, S. Sakakibara2, K.Y. Watanabe2, K. Toi2, J. Geiger1, M.C. Zarnstorff3, S.R. Hudson3, A. Reiman3, A. Werner1, C. Nührenberg1, S. Ohdachi2, Y. Suzuki2, H. Yamada2, W7-AS Team1, LHD Team2 1Max-Planck-Institut für Plasmaphysik, EURATOM-IPP Association, D-17491 Greifswald, Germany 2National Institute for Fusion Science, Toki 509-5292, Japan 3Princeton Plasma Physics Laboratory, Princeton, NJ 08543, USA Corresponding Author: Dr. Arthur Weller ( [email protected] ) Max-Planck-Institut für Plasmaphysik, EURATOM-IPP Association, Wendelsteinstr. 1, D-17491 Greifswald, Germany Fax: +49 (0)3834 88 2509 Content of Paper: Main Text: 27 pages (including references and figure captions) 12 Figures, no Table - 1 - Ms-313705 Final Paper Abstract Substantial progress has been achieved to raise the plasma beta in stellarators and helical systems by high power neutral beam heating, approaching reactor relevant values [1-3]. The achievement of high-β operation is closely linked with configuration effects on the confinement and with magnetohydrodynamic (MHD) stability. The magnetic configurations of the Wendelstein W7-AS stellarator and of the Large Helical Device (LHD) and their optimization for high-β operation within the flexibility of the devices are characterized. A comparative description of the accessible operational regimes in W7-AS and LHD is given. The finite-β effects on the flux surfaces depend on the degree of configuration optimization. In particular, a large Shafranov shift is accompanied by formation of islands and stochastic field regions as found by numerical equilibrium studies [2,4]. However, the observed pressure gradients indicate some mitigation of the effects on the plasma confinement, presumably because of the high collisionality of high-β plasmas and island healing effects (LHD [5-7]).
    [Show full text]
  • Peaceful Uses of Atomic Energy, Held in Point of Discussions Back in 1958
    INTERNATIONAL ATOMIC ENERGY AGENCY Wagramer Strasse 5, P.O. Box 100 A-1400 Vienna, Austria www.iaea.org CONTENT Guenter Mank, Head, Physics Section Alan Aqrawi, Physics Section Division of Physical and Chemical Sciences Division of Physical and Chemical Sciences Department of Nuclear Sciences and Applications Department of Nuclear Sciences and Applications CREDITS Photography: UNOG Archives Geneva (page 4,5,6,16,17,20,24,25,26), Nuclear Fusion Institute RRC "Kurchatov Institute" Russia (page 9,18,19), KFA Kernforschungsanlage Juelich Germany (page 10,11,15), ORNL Oak Ridge National Laboratory USA (page 6), IAEA Archives Vienna (page 8). DESIGN Alan Aqrawi © IAEA, 2008 INDEX FOREWORD Foreword “Celebrating fifty years of fusion… Foreword 2 …entering into the burning plasma era.” Energy in all its forms has always driven been achieved every 1.8 years. New human development. New technologies in developments in science and engineering energy production, starting from the use of fire have led to an optimized magnetic prototype The 2nd Geneva Conference - Introduction 4 itself, have driven economic and social reactor, with corresponding cost savings. development. In the mid-1950s, nuclear Inertial confinement experiments have energy created new hope for an abundant achieved similar progress. The culmination of source of that energy for the world. international collaborative efforts in fusion is the start of construction of the International Declassification 6 To promote this groundbreaking technology, Thermonuclear Experimental Reactor (ITER), and to host a neutral ground for substantive the biggest scientific endeavour the world has scientific debate, the United Nations in 1955 ever seen involving States with more than half organized the first of a series of conferences of the world’s population.
    [Show full text]
  • Plasma Physics and Controlled Fusion Research During Half a Century Bo Lehnert
    SE0100262 TRITA-A Report ISSN 1102-2051 VETENSKAP OCH ISRN KTH/ALF/--01/4--SE IONST KTH Plasma Physics and Controlled Fusion Research During Half a Century Bo Lehnert Research and Training programme on CONTROLLED THERMONUCLEAR FUSION AND PLASMA PHYSICS (Association EURATOM/NFR) FUSION PLASMA PHYSICS ALFV N LABORATORY ROYAL INSTITUTE OF TECHNOLOGY SE-100 44 STOCKHOLM SWEDEN PLEASE BE AWARE THAT ALL OF THE MISSING PAGES IN THIS DOCUMENT WERE ORIGINALLY BLANK TRITA-ALF-2001-04 ISRN KTH/ALF/--01/4--SE Plasma Physics and Controlled Fusion Research During Half a Century Bo Lehnert VETENSKAP OCH KONST Stockholm, June 2001 The Alfven Laboratory Division of Fusion Plasma Physics Royal Institute of Technology SE-100 44 Stockholm, Sweden (Association EURATOM/NFR) Printed by Alfven Laboratory Fusion Plasma Physics Division Royal Institute of Technology SE-100 44 Stockholm PLASMA PHYSICS AND CONTROLLED FUSION RESEARCH DURING HALF A CENTURY Bo Lehnert Alfven Laboratory, Royal Institute of Technology S-100 44 Stockholm, Sweden ABSTRACT A review is given on the historical development of research on plasma physics and controlled fusion. The potentialities are outlined for fusion of light atomic nuclei, with respect to the available energy resources and the environmental properties. Various approaches in the research on controlled fusion are further described, as well as the present state of investigation and future perspectives, being based on the use of a hot plasma in a fusion reactor. Special reference is given to the part of this work which has been conducted in Sweden, merely to identify its place within the general historical development. Considerable progress has been made in fusion research during the last decades.
    [Show full text]
  • Simulation Magnéto-Hydro-Dynamiques Des Edge-Localised-Modes Dans Un Tokamak
    Laboratoire de Physique des Interactions Institut de Recherche sur la Fusion par Ioniques et Moléculaires (PIIM) confinement Magnétique (IRFM) Université de Provence CEA Cadarache Thèse de doctorat de l’Université de Provence Simulation Magnéto-Hydro-Dynamiques des Edge-Localised-Modes dans un tokamak Présentée par Stanislas PAMELA Soutenue le 27 Septembre 2010 devant le jury composé de : Pr. Howard WILSON Professeur au Department of Physics Rapporteur University of York, United-Kingdom Dr. Hinrich LUTJENS Docteur d’état, Chargé de Recherche au CNRS Rapporteur CPHT, Ecole Polytechnique, Palaiseau Dr. Xavier LITAUDON Chef-Adjoint du Service de Chauffage et Président du Jury Confinement du Plasma, IRFM, CEA Cadarache Pr. Sadruddin BENKADDA Professeur au Laboratoire PIIM, Directeur de Recherche Directeur de thèse CNRS, Université de Provence, Aix-Marseille I Dr. Guido HUYSMANS Chercheur au Service de Chauffage et Responsable CEA Confinement du Plasma, IRFM, CEA Cadarache Pr. Peter BEYER Professeur au Laboratoire PIIM Invité Université de Provence, Aix-Marseille I Dr. Pierre RAMET Maître de Conférences, LABRI, INRIA Invité Université de Bordeaux I TABLE OF CONTENTS 1. Introduction ----------------------------------------------------------------------------------- 1 1.1 Nuclear Fusion ----------------------------------------------------------------------- 1 1.1.1 Fusion Energy ------------------------------------------------------------- 1 1.1.2 Energy Confinement ----------------------------------------------------- 2 1.1.3 Inertial and Magnetic
    [Show full text]
  • Seven Ways to Model Fusion Plasmas
    Seven Ways To Model Fusion Plasmas Introduction A few years ago a man named John Hedditch quit his coding job at Google to follow his dream. He wanted to work in fusion research. This was not a sound financial choice. Fusion is under resourced and it certainly does not pay as well as Google. However, Dr. Hedditch had a vision that most people in this space share: a green energy future for the human race. John is an Austrian, so he went to work at the University of Sydney under Dr. Joe Khachan. The team worked hard - and on August 2nd 2018, Dr. Hedditch submitted a paper laying out some new ideas on how to model plasma [2]. John Hedditch (Left) and his boss Joe Khachan (right). Hedditch wants to use field theory to model plasma equilibrium. This is stealing an idea that is commonly used in the world of quantum physics. Field theory asks how charged particles behave in a field. This is like modeling a spider web. The charges on the particles create a web of electromagnetic forces – that push and pull on one another. The math used, is based on this concept. It is a complex idea that would be tough to code in a computer. However, Hedditch is just using the math to uncover some new insights about plasmas in places where there is no magnetic field. Current Problems: You might think we have already mastered plasma modelling – but that is wrong. There is a body of evidence that our models consistently fall short of reality.
    [Show full text]