A Mathematical Introduction to Magnetohydrodynamics

Total Page:16

File Type:pdf, Size:1020Kb

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. 119 5.2 Nonlinear shear Alfv´enwaves. 128 5.3 Magnetic field diffusion. 129 A Energy conservation in extended MHD models 137 B Magnetic vector potential in MHD 141 References 145 2 Preamble Magnetohydrodynamics is the theory of electrically conducting, neutral fluids in the low-frequency regime. The basic equations of magnetohydrodynamics (MHD) have been proposed by Hannes Alfv´en[1, 2], who realized the importance of the electric currents carried by a plasma and the magnetic field they generate. Alfv´encombined the equations of fluid dynamics with Faraday's and Amp`ere's laws of electrodynamics, thus obtaining a novel mathematical theory, which helped understanding space plasmas in Earth and planetary magnetospheres, as well as the physics of the Sun, solar wind, and stellar atmospheres. In fusion research, MHD is crucial to the understanding of plasma equilibria and their stability. Liquid metals and electrolytes, like salt water, can also be modeled by MHD equations. Besides the important physical applications, MHD equations exhibit a re- markably beautiful mathematical structure, with connections to geometry and topology that allows us to understand some of the dynamics of magnetic fields in plasmas in terms of topological ideas [3, 4, 5, 6, 7, 8, 9]. As a dynamical system MHD is an example of infinite-dimensional Hamiltonian system [10, 11]. The scope of this lecture is, in this regard, extremely limited. The goal is to introduce MHD equations in a reasonably self-contained way and to discuss some of their most important features. The style of this lectures is quite similar to the mathematical introduction to fluid dynamics by Chorin and Marsden [12] as the title of this note suggests. Particularly, we shall attempt to introduce the physics modeling in a mathematically precise albeit not always rigorous way. The physics literature on the subject is vast. As a reference for further read- ing, the book by Biskamp [13] provides a clear and comprehensive exposition, while the lectures by Schnack [14] offer a more gradual learning curve. For MHD of the solar atmosphere one can refer to Priest [15] as well as to Aschwanden's book on the solar corona [16]. For an introduction to magnetohydrodynamics with emphasis on equilibria and stability of fusion plasmas one can refer to the books by Freidberg [17] and Zohm [18], while Goedbloed, Poedts and Keppens address applications to both astrophysical and fusion plasmas [19, 20]. A nice introduction to MHD with a broader perspective which includes applications to metals can be found in Davidson's book [21]. On the mathematical side, MHD has received considerable attention from ap- plied mathematicians. Its rich mathematical structure has become a paradigm for the application of geometry and topology [9, 22] as well as for structure preserving discretization [23, 24, 25, 26]. As a system of partial differential equations, well-posedness of the Cauchy problem for MHD equations subject to appropriate boundary conditions have been studied first by Duvaut and Lions [27, in French], Sermange and Temam [28], Secchi [29] and more recently, by Chen and coworkers [30] and by Fefferman and coworkers [31, and references therein]. 3 1 Basic elements of fluid dynamics The basic understanding of fluid dynamics is an essential prerequisite to the study of MHD. We shall start by recalling the basic elements thereof, following Chorin and Marsden [12], cf. also Marsden and Hughes [32]. First, we define the physical quantities that describe the dynamical state of a fluid (kinematics) and continue with the equations of motion (dynamics). 1.1 Kinematics of fluids. Fluid dynamics is built on the basis of the con- tinuum hypothesis: A fluid is a distribution of matter occupying a certain region of the continuous three-dimensional space. The considered region of space is a domain (i.e., an open and connected subset) Ω ⊆ R3. We neglect the fact that any fluid is ultimately made of atoms and molecules as we are interested in studying its collective motion on a much larger spatial scale. With the continuum hypothesis, one needs to quantify how matter is dis- tributed in Ω at any time t in a certain time interval I ⊆ R. Thus, the first physical quantity of interest is the mass density, which is a non-negative time- dependent scalar field ρ : I × Ω ! R, such that ρ(t; x) ≥ 0 (1.1) gives the mass per unit of volume at time t 2 I at the spatial location x 2 Ω. By definition, the amount of mass contained in an arbitrary volume W ⊆ Ω (to be referred to as a control volume) is given by Z (mass in W at time t) = ρ(t; x)dx; W which implies that the mass density must be at least locally integrable. Physically we think of an infinitesimal volume of fluid centered around a point x 2 Ω. The volume of this infinitesimal region is mathematically repre- sented by the Lebesgue measure dx in Ω and the mass is represented by the measure ρ(t; x)dx. Such infinitesimal portions of fluid are referred to as fluid elements. In plasma physics the mass density is often replaced by an equivalent non- negative scalar field referred to as the particle number density or simply number density which is defined in terms of the mass density by n(t; x) = ρ(t; x)=m: Here the fluid is regarded as a collection of particles that all have the same mass m; thus, the number of particles contained in a control volume W is given by 1 Z Z (number of particles in W at time t) = ρ(t; x)dx = n(t; x)dx: m W W According to this definition, the number of particles does not need to be an integer, due to the continuum hypothesis. Next we need to describe the motion of the fluid. We introduce a velocity field defined as a time-dependent vector field, u : I × Ω ! R3 such that, u(t; x) 2 R3 (1.2) 4 gives the velocity of the fluid element at the point x 2 Ω and time t 2 I. The vector u(t; x) is referred to as the fluid velocity. The mass times the velocity of the fluid element, namely, ρ(t; x)u(t; x)dx, gives the linear momentum of the fluid element, hence Z (momentum in W at time t) = ρ(t; x)u(t; x)dx: W In addition to the quantities ρ(t; x) and u(t; x), that can be regarded as the counterpart in fluid dynamics of mass and velocity of a particle in mechan- ics, we need to specify another scalar field for the internal energy of the fluid element. Differently from a point-mass particle, a fluid element is a thermody- namical system that can undergo expansions and compressions, thus absorbing and releasing energy. The thermodynamic status of a fluid element is specified by the internal energy density U : I × Ω ! R≥0, where R≥0 denotes the set of non-negative real numbers. If each fluid element is regarded as an ideal gas composed by n(t; x)dx particles, we can equivalently express the internal energy density in terms of a new variable. Specifically the laws of thermodynamics for a perfect gas allow us to write the internal energy of a fluid element as 3 U(t; x)dx = n(t; x)k T (t; x)dx; 2 B where kB is the Boltzmann constant and T : I × Ω ! R is a strictly positive scalar field, such that T (t; x) > 0 (1.3) represents the local temperature of the fluid element at time t 2 I and position x 2 Ω. Therefore, the total energy carried by a fluid element is the sum of the kinetic energy associated to its motion plus the internal energy associated to its thermodynamics, namely, Z 1 2 3 (energy in W at time t) = ρ(t; x)u(t; x) + n(t; x)kBT (t; x) dx: W 2 2 When the law of a perfect gas does not apply we can still define T as above, but now it has the meaning of an effective temperature, which is just a measure of the internal energy and does not imply thermodynamical equilibrium.
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]
  • 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]
  • 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]
  • Plasma Physics and Fusion Energy
    This page intentionally left blank PLASMA PHYSICS AND FUSION ENERGY There has been an increase in worldwide interest in fusion research over the last decade due to the recognition that a large number of new, environmentally attractive, sustainable energy sources will be needed during the next century to meet the ever increasing demand for electrical energy. This has led to an international agreement to build a large, $4 billion, reactor-scale device known as the “International Thermonuclear Experimental Reactor” (ITER). Plasma Physics and Fusion Energy is based on a series of lecture notes from graduate courses in plasma physics and fusion energy at MIT. It begins with an overview of world energy needs, current methods of energy generation, and the potential role that fusion may play in the future. It covers energy issues such as fusion power production, power balance, and the design of a simple fusion reactor before discussing the basic plasma physics issues facing the development of fusion power – macroscopic equilibrium and stability, transport, and heating. This book will be of interest to graduate students and researchers in the field of applied physics and nuclear engineering. A large number of problems accumulated over two decades of teaching are included to aid understanding. Jeffrey P. Freidberg is a Professor and previous Head of the Nuclear Science and Engineering Department at MIT. He is also an Associate Director of the Plasma Science and Fusion Center, which is the main fusion research laboratory at MIT. PLASMA PHYSICS AND FUSION ENERGY Jeffrey P. Freidberg Massachusetts Institute of Technology CAMBRIDGE UNIVERSITY PRESS Cambridge, New York, Melbourne, Madrid, Cape Town, Singapore, São Paulo Cambridge University Press The Edinburgh Building, Cambridge CB2 8RU, UK Published in the United States of America by Cambridge University Press, New York www.cambridge.org Information on this title: www.cambridge.org/9780521851077 © J.
    [Show full text]
  • LECTURE NOTES on IDEAL MAGNETOHYDRODYNAMICS By
    ASSOCIATIE EURATOM-FOM FOM-INSTITUUT VOOR PLASMAFYSICA RIJNHUIZEN - NIEUWEGEIN - NEDERLAND LECTURE NOTES ON IDEAL MAGNETOHYDRODYNAMICS by J.P. GoedbJoed Rijnhuizen Report 83-145 ASSCOCIATIE EURATOM-FOM Maart 1983 FOM-INSTITUUT VOOR PLASMAFYSICA RIJNHUIZEN - NIEUWEGEIN - NEDERLAND LECTURE NOTES ON IDEAL MAGNETOHYDRODYNAMICS by J.P. Goedbloed Rijnhuizen Report 83-145 Corrected version of the notes of March 1979, originally printed as internal report at Instituto de Ffsica, Universidade Estadual de Campinas, Campinas, Brazil This work was supported by the "Stichting voor Fundamentaal Onderzoek der Materie" (FOM), the "Nederlandse Organisatie voor Zuiw-WetenschappelijK Onderzoek" (ZWO), EURATOM, the "Fundacio de Amparo i Pesquiw do Ettado de Sao Paulo" (FAPÉSP), and the "Conaelho Nacione) de PesquisM" (CNPQ, Brazil). "Then I saw that all toil and skill in work come from a man's envy of his neighbour. This also is vanity and a striving after wind." Ecclessiastes 4:4 "Ever since the creation of the world his invisible nature, namely, his eternal power and deity, has been clearly perceived in the things that have been made." Romans I:20 "Remember then to sing the praises of his work, as men have always sung them." Job 36:24 PREFACE These notes were prepared for a course of lectures for staff and students of the Instituto de Flsica, Universidade Estadual de Campinas, Brazil. The course consisted of two-hour lectures twice a week during a period of 9 weeks in the months June-August 1978. It has been my intention to make the subject- matter as much as possible self-contained, so that all needed physical and mathematical techniques and derivations were pre­ sented in detail.
    [Show full text]
  • Two-Fluid Physics and Field-Reversed Configurations
    PHYSICS OF PLASMAS 14, 055911 ͑2007͒ Two-fluid physics and field-reversed configurationsa… ͒ A. Hakimb Tech-X Corporation, 5621 Arapahoe Avenue – Suite A, Boulder, Colorado 80303 U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle, Washington 98195-2600 ͑Received 6 November 2006; accepted 30 April 2007; published online 24 May 2007͒ In this paper, algorithms for the solution of two-fluid plasma equations are presented and applied to the study of field-reversed configurations ͑FRCs͒. The two-fluid model is more general than the often used magnetohydrodynamic ͑MHD͒ model. The model takes into account electron inertia, charge separation, and the full electromagnetic field equations, and it allows for separate electron and ion motion. The algorithm presented is the high-resolution wave propagation scheme. The wave propagation method is based on solutions to the Riemann problem at cell interfaces. Operator splitting is used to incorporate the Lorentz and electromagnetic source terms. The algorithms are benchmarked against the Geospace Environmental Modeling Reconnection Challenge problem. Equilibrium of FRC is studied. It is shown that starting from a MHD equilibrium produces a relaxed two-fluid equilibrium with strong flows at the FRC edges due to diamagnetic drift. The azimuthal electron flow causes lower-hybrid drift instabilities ͑LHDI͒, which can be captured if the ion gyroradius is well resolved. The LHDI is known to be a possible source of anomalous resistivity in many plasma configurations. LHDI simulations are performed in slab geometries and are compared to recent experimental results. © 2007 American Institute of Physics. ͓DOI: 10.1063/1.2742570͔ I. INTRODUCTION electromagnetism.
    [Show full text]