Simple examples of MHD equilibria
|
|
- Allen Skinner
- 5 years ago
- Views:
Transcription
1 Department of Physics Seminar. grade: Nuclear engineering Simple examples of MHD equilibria Author: Ingrid Vavtar Mentor: prof. ddr. Tomaž Gyergyek Ljubljana, 017 Summary: In this seminar paper I will represent the magnetohydrodynamics description (MHD) of the plasma. In the first part I will focus on the MHD equations and their properties. In the second part of the seminar paper the basic properties of the MHD equilibrium model will be presented; I will focus on radial pressure balance in cylindrical system, namely θ-pinch, Z-pinch and screw pinch. 1
2 TABLE OF CONTENTS 1 INTRODUCTION... MHD equations Basic scaling relations for MHD Summary of the MHD equations THE MHD EQUILIBRIUM MODEL General properties Flux surfaces Current surfaces Magnetic pressure and tension Radial pressure balance The θ-pinch The Z-pinch The screw pinch CONCLUSION SOURCES AND LITERATURE INTRODUCTION Fusion is a form of nuclear energy which involves the merging (i.e., the fusing) of light elements, mainly hydrogen (H) and its isotopes deuterium (D) and tritium (T). The main application of fusion energy is the production of electricity, because the laws of physics have shown that a large amount of kinetic energy is released every time a nuclear fusion reaction occurs. To achieve fusion reaction it is necessary that the material (fuel) is in the form of plasma (plasma is hot gas made of electrons and bare nuclei) and that fusion plasma achieves a temperature of about 15 kev and a pressure of approximately 7 atm. Fusion plasma is a fully ionized gas whose behaviour is dominated by long-range electric and magnetic fields. A major consequence of this behaviour is that plasma is an exceptionally good conductor of electricity. The reason is associated with the fact that at a high temperature and because of the low density of plasma, there are rarely any Coulomb collisions between electrons and ions and thus there is very little resistance to the current flow. In terms of fusion reactor it is the high conductivity that opens up the possibility of confinement by magnetic fields. Carrying out a workable fusion reactor requires good knowledge of plasma properties, therefore my focus in this seminar paper will be on the description of the plasma. I will introduce a magnetic hydrodynamic description of the plasma, the adequacy thereof and equilibrium magnetic configuration on fusion reactors. Plasma can be described with a self-consistent plasma model, where it is considered that the charges in the plasma which move under the influence of electric and magnetic fields affect
3 those fields back. In developing self-consistent models we should be aware that various levels of description are possible. The most accurate models are kinetic models which strive to determine the distribution functions of ions and electrons and then the macroscopic quantities are determined as the moments of the distribution functions. They are the most accurate and describe almost all physical phenomena in plasma but are mathematically very demanding. The next level of description corresponds to macroscopic fluid models, where the great simplification in comparison to kinetic models is that all the unknowns are functions of only space and time, but not of the velocity. In comparison with the kinetic models they are less accurate and they cannot describe some phenomena which can be described in kinetic models, but they are mathematically less demanding. In a simplified two-fluid model it is assumed that plasma contains two fluids ; electrons and ions. The two-fluid model consists of the conservation of mass, momentum, energy for electrons and ions plus Maxwell s equations and describes the self-consistent behaviour of a fusion plasma including all of the main phenomena such as macroscopic equilibrium and stability, transport, heating and current drive. The third level of description is magnetohydrodynamics (MHD) on which I will focus in the following. The MHD model is a reduction of the two-fluid model to one-fluid model derived by focusing attention on the length and time scales characteristic of macroscopic behavior. MHD equations The major issue in which self-consistency plays a crucial role is the macroscopic equilibrium. It needs to be known how a magnetic field can produce forces to hold plasma in a stable, macroscopic equilibrium thereby allowing fusion reactions to take place in a continuous, steady state mode of operation. The analysis of macroscopic equilibrium and stability is based on a single-fluid model known as MHD. The derivation of the MHD model from the two-fluid model is relatively simple, but requires a few steps. In the following I will present MHD equations and their characteristics..1 Basic scaling relations for MHD It is necessary to define the characteristic spatial and temporal criteria to allow comparison orders of magnitude of individual parts of equations in the two-fluid model and to allow neglecting small parts in equations. The appropriate length scale L is the plasma radius (L~a) while the appropriate time scale τ is the ion thermal transit time across the plasma ( τ~a/v Ti ), which leads to a characteristic velocity u (u~l/τ~v Ti ). v Ti is the fastest macroscopic speed that the plasma can achieve, a.k.a. the ion sound speed. The ion sound speed v Ti assimilated to the average ion velocity in two-fluid model u i. In a fusion plasma there is E B drift (because of a uniform perpendicular electric field superimposes a constant drift velocity on the gyro motion known as the E B drift and is perpendicular to both field, electric and magnetic; it is also independent of the mass and charge, so electrons and ions drift with same velocity.). Even E B drift speed should not exceed the ion sound speed (u e ~u i ~E /B ~v Ti ). Since the E B drift velocity is identical for both electrons and ions the current densities J = e(n i u i n e u e ) cancels to leading order. It follows J /env Ti 1. The next request is that the plasma pressure is assumed to be finite compared to the magnetic pressure, which requests that β p/( B μ 0 )~1. 3
4 Inequalities which must be satisfied are also small ion gyro radius r Li a and low collision frequency ν ei v Ti /a.. Summary of the MHD equations The MHD equations which are collected together in the following are the basis for studies of macroscopic balance in fusion plasma studies. First equation of MHD model is mass equation: dρ dt + ρ v = 0. (1) The momentum equation is the most important one in MHD as it describes the basic force balance of the plasma: ρ dv dt = J B p. () The left side of the equation is a product of density and acceleration of plasma which represents the inertial force which is important in determining the dynamical behaviour of the plasma. On the right side of the equation are the magnetic field force used to confine the plasma and pressure gradient force that causes a hot core of plasma to expand outwards. In a steady state situation without flow the inertial force is zero and equilibrium is achieved when the magnetic force balances the pressure gradient force. In terms of scaling it can be seen that the current density J must satisfy J~ p J for the two forces to balance. This implies that ~ r Li 1 as ab e 0 nv Ti a previously asserted. [1] If the MHD equilibrium is achieved, the resistive Ohm s law can be derived from equations of conservation of momentum in two-fluid model: E + v B = ηj. (3) A more detailed kinetic analysis shows that the resistance of the plasma varies considerably in the direction which is parallel to B, as in the direction perpendicular to it. The resistance η in fact corresponds to the rectangular part of the resistance; parallel part of the resistance is about half the size. Taking into account the resistance of plasma in Ohm s law brings two important effects. First, use of the resistive Ohm s law allows a treatment of a large number of plasma instability. The time scale of these instabilities, however, is much larger than the time scale for MHD, which was defined at the beginning, so that instabilities do not lead to loss of plasma in its entirety, but only to increased transport losses. Secondly, even if there were no instabilities as a result of the final resistance of the plasma it is a fact that the resistance of the plasma is only a dissipative mechanism in the equation of motion. This causes the diffusion of the particles and the magnetic field which are the two main mechanisms of transport losses. However, these appear as slower than the characteristic MHD time constant. In the MHD regime the energy equation is simplified considerably, because all sources and sinks of energy can be neglected, since all the processes of heating and cooling are much slower than the MHD time constant τ, which leads to the single-fluid energy equation: 4
5 d dt ( p ρ γ) = 0, (4) where γ gives the ratio of specific heats and his value is γ = 5 3. Some simplifications of Maxwell s equations can be made. Taking into account that the speed of macroscopic plasma flows is much smaller than the speed of lights (v Ti c), the displacement current 1 E can be neglected. Because characteristic length scale of MHD easily c t satisfies the inequality λ De a, the E term in Poisson s equation can be neglected leading to the quasi-neutrality relation: n e n i. Quasi-neutrality implies that ε 0 E en e (not that E = 0). These two simplifications reduce Maxwell s equations from the exact relativistic Lorentz-invariant form to a self-consistent, low-frequency, Galilean-invariant form: E = B t, (5) B = μ 0 J, B = 0. In the equations above single fluid variables, mass density ρ, macroscopic velocity v, and the pressure p appear. J denotes the current densities, B the magnetic field and E the electric field. 3 THE MHD EQUILIBRIUM MODEL An important application of the MHD model is concerned with the calculation of equilibrium. In MHD equilibrium model we are looking for an answer to the question of how does a combination of externally applied and internally induced magnetic fields need to act to provide an equilibrium force balance that holds the plasma together at the desired location in the vacuum chamber. The equilibria of interest must correspond to the confined equilibria, where the purpose of the magnetic field is to isolate the plasma from the wall of the vacuum chamber which keeps the plasma hot and the wall cool. The MHD equilibrium in a toroidal geometry can be separated into two pieces, namely the radial pressure balance (along the minor radius r) and toroidal force balance (along the major radius R) (see Figure 1). The plasma is a hot core of gas that tends to expand uniformly along the minor radius r. Two basic magnetic configurations that can produce radial pressure balance are the θ-pinch and Z-pinch plus combinations thereof. These configurations will be discussed below. A problem of toroidal force balance arises solely because of the toroidal geometry. Because of toroidicity, unavoidable forces are generated by both the toroidal and poloidal magnetic fields that tend to push the plasma horizontally outwards along the direction of the major radius R. The toroidal force balance is outside the frameworks of this seminar paper. 5
6 Figure 1: On the left radial MHD balance and on the right torodal MHD balance. Reproduced from [1]. To focus only on the equilibrium the MHD model needs to be simplified. The simplifications result from two basic assumptions. First assumption is that for plasma in equilibrium all quantities are independent of time, which results in = 0. The second assumption is that the t plasma is treated as a stationary liquid, v = 0. With these assumptions the MHD equations are converted in the following form: J B = p, B = μ 0 J, (6) B = 0. These equations describe the equilibrium properties of all magnetic configurations of fusion interest. 3.1 General properties Flux surfaces The first property concerns the concept of flux surfaces. Vector p is, by definition, perpendicular to the surface where p is a constant. So if the first equation from simplified MHD model (6) is multiplied with the magnetic field B it immediately follows that: B p = 0. (7) In equation (7) the rule for mixed product of three vectors is taken into account. The implication is that the magnetic field lines must lie on isobaric surfaces, so that there is no component of B perpendicular to the surface. These surfaces are called flux surfaces and are forming a set of closed, nested, toroidal surfaces Current surfaces Current surfaces need to be proceeded similarly as flux surfaces above. Momentum equation (6) forming the dot product with current densities J results in: J p = 0. (8) The current lines, like magnetic field, lie in the constant pressure surfaces (see Figure ) and there is also no component of J perpendicular to the pressure contours, which implies that 6
7 current flows between flux surfaces and not across them. Vectors B and J thus both lie on the isobaric surfaces, but we cannot conclude anything about what is the angle between them. Figure : Vectors B and J both lie on the isobaric surfaces. On right reproduced from [1], on left reproduced from [] Magnetic pressure and tension Magnetic pressure and tension are two ways in which the magnetic field can act to hold the plasma in equilibrium force balance. Magnetic pressure uniformly presses on the plasma from all sides and tension extends plasma in the direction of magnetic field lines. Combining first two equations (6) and the vector identity ( B ) = B ( B) + (B )B yields to: (p + B ) B κ = 0, (9) μ 0 μ 0 where is the perpendicular component of the gradient operator and κ is the curvature vector. Equation (9) describes pressure balance perpendicular to the magnetic field. The quantity p represents the plasma pressure, quantity B B μ 0 represents the magnetic pressure and the last term μ 0 κ represents the tension force created by the curvature of the field lines. So the magnetic field distresses plasma in the direction perpendicular to the magnetic field and stretches it in the parallel direction (see Figure 3). [1] Figure 3: Magnetic field distresses plasma in the direction perpendicular to the magnetic field and stretches it in the parallel direction. Reproduced from []. 7
8 3. Radial pressure balance For further analysis simpler cylindrical geometry will be used, where the torus will be stretched in the cylinder. This way the problem is much more simplified because it becomes onedimensional (it is necessary to observe only the radial compression of the plasma) The θ-pinch Consider a cylindrical column of plasma in a conductive, thin, metal panel, where the axis of the plasma column and metal panel coincide. Electric current is running through this panel in the poloidal direction. In the plasma a current is induced in the opposite direction and it creates a magnetic field in the axial direction, i.e., toroidal direction. Tags poloidal and toroidal are understood in the terms of a toroid which is extended into the cylinder. The geometry and field components of a θ-pinch are illustrated in Figure 4. Figure 4: Schematic diagram of a θ-pinch on left and typical θ-pinch profiles on right. Reproduced from [1]. The non-trivial unkowns in the problem are p = p(r), B = B z (r)e z and J = J θ (r)e θ, where current density in the θ direction is the origin of the name θ-pinch. First, by symmetry the condition. B = 0 is automatically satisfied. Second, Ampère s law reduces to μ 0 J θ = db z dr. Substituting J θ into the momentum equation leads to a simple relation between p and B z, which can be written as d (p + B z ) = 0 and can be easily integrated: dr μ 0 p(r) + B z (r) = B 0. (10) μ 0 μ 0 In equation (10) the part B 0 = const. is the externally applied magnetic pressure. Equation μ 0 (10) represents the basic radial pressure balance relation for a θ-pinch. It states that at any radial position r the sum of the local plasma pressure plus internal magnetic pressure is equal to the applied magnetic pressure. Equation (10) can be solved with different radial dependence of the pressure in a magnetic field. We are interested in a balance where the plasma is confined which means that p(r) has the maximum at r = 0 and then monotonically decreases towards the edge of the plasma column at r = a. The radius of the metal panel should be b and b > a. A typical example is illustrated in Figure 4. [1] 3.. The Z-pinch The Z-pinch is the complementary configuration of the θ-pinch. Like in the case of θ-pinch, it is considered a cylindrical column of plasma, but now the electric current in the column flows 8
9 in a toroidal direction, i.e., in z direction, which produces a purely poloidal magnetic field in the θ direction. It is given the name Z-pinch because of the current density in the z direction. The geometry and field components of a Z-pinch are illustrated in Figure 5 on the left side. Figure 5: On the left is shematic diagram o a Z-pinch, in the middle typical Z-pinch profiles, on the right the radial forces on a Z-pinch. Reproduced from [1]. In this configuration the non-trivial unknowns are given by p = p(r), B = B θ (r)e θ and J = J z (r)e z. Even in this case it is so because the symmetry. B = 0. Since the magnetic field B θ on the closed circular lines is constant, Ampère s law yields a relation for J z given by μ 0 J z = 1 d r dr (rb θ). Z-pinch pressure balance relation can be written as: d dr (p + B θ ) + B θ = 0, (11) μ 0 rμ 0 if we expand rb θ derivative from pressure balance relation which is obtained by substituting the expression p = J z B θ into momentum equation. In order, from left to right, the terms correspond to the plasma pressure, the magnetic pressure and the magnetic tension. In this case the equation (11) cannot be directly integrated, therefore we need to invent a physically acceptable, radial dependence density of the electric current flow J z (r) through the plasma column from which poloidal magnetic field B θ (r) and pressure profile p(r) can be then calculated. Typical Z-pinch profiles are illustrated in Figure 5 in the middle. In equation (11) the positive force dp wants to expand the plasma column in the radial direction and the negative dr forces on the right side are pressing on it (see Figure 5 on the right). It can be seen that only the tension force is confining the outer portion of the plasma. The conclusion is that the Z-pinch configuration is capable of providing radial pressure balance by means of the magnetic tension force. [1] 3..3 The screw pinch The screw pinch is a combination of a θ-pinch and Z-pinch. In its configuration the magnetic lines twisted around the surface giving the appearance of a screw thread and hence the name screw pinch. Using a very similar method to the one in θ-pinch and Z-pinch, the desired expression for radial pressure balance in a generalized screw pinch yields: d dr (p + B θ + B z ) + B θ = 0. (1) μ 0 μ 0 rμ 0 9
10 It is a differential equation with three variables, two of which we can choose on our own. Both are usually a magnetic field. This can be done in an experiment, for example, if we create magnetic fields in axial and poloidal direction. Both can be independently modified. [1] 4 CONCLUSION In the seminar paper, magnetic hydrodynamic description (MHD) of the plasma is discussed. The first part of the seminar paper is focused on the MHD equations and its properties. The MHD model is reduction of the two-fluid model to the one-fluid model derived by focusing attention on the length and time scales characteristic of macroscopic bahavior. MHD model contains the mass and momentum equations, the resistive Ohm s law, energy equation and three Maxwell s equations. An important application of the MHD model is the analysis of equilibrium. The second part of this seminar paper was focused on it. To focus only on the equilbrium, the MHD model was simplified. The MHD equilibrium in a toroidal geometry can be separated into two pieces, namely radial pressure balance and toroidal force balance. In the seminar paper, only radial pressure balance was presented. The analysis of general properties like flux surfaces, current surfaces and magnetic pressure and tension points to the notion that pressure contours from a set of closed-nested toroidal surfaces. Both the magnetic field lines and current density lines lie on these surfaces. The magnetic forces that hold the plasma together arise from a combination of magnetic pressure and magnetic tension. Finally, several basic configurations like θ-pinch, Z-pinch and screw pinch were researched using the simplified MHD model. The conclusion is that the θ-pinch and Z-pinch configurations are capable of providing radial presure balance and that the screw pinch is a configuration that consists of an arbitrary combination of θ-pinch and Z-pinch fields. 5 SOURCES AND LITERATURE [1] J. P. Freidberg. Plasma physics and fusion energy. Cambridge University Press, New York 007. [] J. A. Bittencourt. Fundamentas of plasma physics 3 rd edition. Springer-Verlag, New York 004. [3] Glossary of terms in nuclear technology and radiation protection. Nuclear society of Slovenia, Ljubljana 01. [4] T. Gyergyek. Lecture notes: Physics and engineering of fusion reactors. Faculty of mathematics and physics, University of Ljubljana, Ljubljana 016. [5] J. P. Freidberg. Ideal MHD. Cambridge University Press, New York 014. [6] H. P. Goedbloed, S. Poedts. Princiles of magnetohydrodynamics. Cambridge University Press, New York
MAGNETOHYDRODYNAMIC EQUILIBRIUM AND STABILITY OF PLASMA
University of Ljubljana Faculty of Mathematics and Physics Seminar 1 b -1st year, II. cycle MAGNETOHYDRODYNAMIC EQUILIBRIUM AND STABILITY OF PLASMA Author: Lino alamon Advisor: prof. dr. Tomaº Gyergyek
More information0 Magnetically Confined Plasma
0 Magnetically Confined Plasma 0.1 Particle Motion in Prescribed Fields The equation of motion for species s (= e, i) is written as d v ( s m s dt = q s E + vs B). The motion in a constant magnetic field
More informationMAGNETOHYDRODYNAMICS
Chapter 6 MAGNETOHYDRODYNAMICS 6.1 Introduction Magnetohydrodynamics is a branch of plasma physics dealing with dc or low frequency effects in fully ionized magnetized plasma. In this chapter we will study
More informationThe Virial Theorem, MHD Equilibria, and Force-Free Fields
The Virial Theorem, MHD Equilibria, and Force-Free Fields Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 10 12, 2014 These lecture notes are largely
More informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More informationIdeal Magnetohydrodynamics (MHD)
Ideal Magnetohydrodynamics (MHD) Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 1, 2016 These lecture notes are largely based on Lectures in Magnetohydrodynamics
More informationIdeal MHD Equilibria
CapSel Equil - 01 Ideal MHD Equilibria keppens@rijnh.nl steady state ( t = 0) smoothly varying solutions to MHD equations solutions without discontinuities conservative or non-conservative formulation
More informationFormation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )
Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University
More informationMagnetohydrodynamic waves in a plasma
Department of Physics Seminar 1b Magnetohydrodynamic waves in a plasma Author: Janez Kokalj Advisor: prof. dr. Tomaž Gyergyek Petelinje, April 2016 Abstract Plasma can sustain different wave phenomena.
More informationMacroscopic plasma description
Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion
More informationFundamentals of Magnetohydrodynamics (MHD)
Fundamentals of Magnetohydrodynamics (MHD) Thomas Neukirch School of Mathematics and Statistics University of St. Andrews STFC Advanced School U Dundee 2014 p.1/46 Motivation Solar Corona in EUV Want to
More informationCharged particle motion in external fields
Chapter 2 Charged particle motion in external fields A (fully ionized) plasma contains a very large number of particles. In general, their motion can only be studied statistically, taking appropriate averages.
More informationTurbulence and Transport The Secrets of Magnetic Confinement
Turbulence and Transport The Secrets of Magnetic Confinement Presented by Martin Greenwald MIT Plasma Science & Fusion Center IAP January 2005 FUSION REACTIONS POWER THE STARS AND PRODUCE THE ELEMENTS
More informationEffect of magnetic reconnection on CT penetration into magnetized plasmas
Earth Planets Space, 53, 547 55, 200 Effect of magnetic reconnection on CT penetration into magnetized plasmas Yoshio Suzuki, Takaya Hayashi 2, and Yasuaki Kishimoto Naka Fusion esearch Establishment,
More informationReduced MHD. Nick Murphy. Harvard-Smithsonian Center for Astrophysics. Astronomy 253: Plasma Astrophysics. February 19, 2014
Reduced MHD Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 19, 2014 These lecture notes are largely based on Lectures in Magnetohydrodynamics by Dalton
More informationIntroduction to Fusion Physics
Introduction to Fusion Physics Hartmut Zohm Max-Planck-Institut für Plasmaphysik 85748 Garching DPG Advanced Physics School The Physics of ITER Bad Honnef, 22.09.2014 Energy from nuclear fusion Reduction
More informationINTRODUCTION TO MAGNETIC NUCLEAR FUSION
INTRODUCTION TO MAGNETIC NUCLEAR FUSION S.E. Sharapov Euratom/CCFE Fusion Association, Culham Science Centre, Abingdon, Oxfordshire OX14 3DB, UK With acknowledgments to B.Alper for use of his transparencies
More informationTURBULENT TRANSPORT THEORY
ASDEX Upgrade Max-Planck-Institut für Plasmaphysik TURBULENT TRANSPORT THEORY C. Angioni GYRO, J. Candy and R.E. Waltz, GA The problem of Transport Transport is the physics subject which studies the physical
More informationDirect drive by cyclotron heating can explain spontaneous rotation in tokamaks
Direct drive by cyclotron heating can explain spontaneous rotation in tokamaks J. W. Van Dam and L.-J. Zheng Institute for Fusion Studies University of Texas at Austin 12th US-EU Transport Task Force Annual
More informationToroidal flow stablization of disruptive high tokamaks
PHYSICS OF PLASMAS VOLUME 9, NUMBER 6 JUNE 2002 Robert G. Kleva and Parvez N. Guzdar Institute for Plasma Research, University of Maryland, College Park, Maryland 20742-3511 Received 4 February 2002; accepted
More informationIntroduction to MagnetoHydroDynamics (MHD) Antoine Cerfon, Courant Institute, New York University
Introduction to MagnetoHydroDynamics (MHD) Antoine Cerfon, Courant Institute, New York University Email: cerfon@cims.nyu.edu SULI Introductory Course in Plasma Physics, June 6, 2016 PART I: DESCRIBING
More informationThe Linear Theory of Tearing Modes in periodic, cyindrical plasmas. Cary Forest University of Wisconsin
The Linear Theory of Tearing Modes in periodic, cyindrical plasmas Cary Forest University of Wisconsin 1 Resistive MHD E + v B = ηj (no energy principle) Role of resistivity No frozen flux, B can tear
More informationTokamak Fusion Basics and the MHD Equations
MHD Simulations for Fusion Applications Lecture 1 Tokamak Fusion Basics and the MHD Equations Stephen C. Jardin Princeton Plasma Physics Laboratory CEMRACS 1 Marseille, France July 19, 21 1 Fusion Powers
More informationPlasmoid Motion in Helical Plasmas
Plasmoid Motion in Helical Plasmas Ryuichi ISHIZAKI and Noriyoshi NAKAJIMA National Institute for Fusion Science, Toki 509-5292, Japan (Received 12 December 2009 / Accepted 18 May 2010) In order to explain
More informationLecture # 3. Introduction to Kink Modes the Kruskal- Shafranov Limit.
Lecture # 3. Introduction to Kink Modes the Kruskal- Shafranov Limit. Steve Cowley UCLA. This lecture is meant to introduce the simplest ideas about kink modes. It would take many lectures to develop the
More informationStabilization of sawteeth in tokamaks with toroidal flows
PHYSICS OF PLASMAS VOLUME 9, NUMBER 7 JULY 2002 Stabilization of sawteeth in tokamaks with toroidal flows Robert G. Kleva and Parvez N. Guzdar Institute for Plasma Research, University of Maryland, College
More informationWaves in plasma. Denis Gialis
Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.
More informationIntroduction to Fusion Physics
Introduction to Fusion Physics J. W. Haverkort October 15, 2009 Abstract This is a summary of the first eight chapters from the book Plasma Physics and Fusion Energy by Jeffrey P. Freidberg. Contents 2
More informationChapter IX: Nuclear fusion
Chapter IX: Nuclear fusion 1 Summary 1. General remarks 2. Basic processes 3. Characteristics of fusion 4. Solar fusion 5. Controlled fusion 2 General remarks (1) Maximum of binding energy per nucleon
More informationASTR-3760: Solar & Space Physics...Spring 2017
ASTR-3760: Solar & Space Physics...Spring 2017 Review material for midterm exam (March 22, 2017) Although I m not recommending full-on memorization of everything in this document, I do think it s important
More informationSMR/ Summer College on Plasma Physics. 30 July - 24 August, Introduction to Magnetic Island Theory.
SMR/1856-1 2007 Summer College on Plasma Physics 30 July - 24 August, 2007 Introduction to Magnetic Island Theory. R. Fitzpatrick Inst. for Fusion Studies University of Texas at Austin USA Introduction
More informationTokamak Divertor System Concept and the Design for ITER. Chris Stoafer April 14, 2011
Tokamak Divertor System Concept and the Design for ITER Chris Stoafer April 14, 2011 Presentation Overview Divertor concept and purpose Divertor physics General design considerations Overview of ITER divertor
More informationChapter 5 MAGNETIZED PLASMAS. 5.1 Introduction. 5.2 Diamagnetic current
Chapter 5 MAGNETIZED PLASMAS 5.1 Introduction We are now in a position to study the behaviour of plasma in a magnetic field. In the first instance we will re-examine particle diffusion and mobility with
More informationThe physics of fusion power. A.G. Peeters,
The physics of fusion power A.G. Peeters, May 8, 2008 2 PREFACE These lecture notes give a first introduction into the physics processes of importance to fusion research. You are looking at the first version
More informationSensors Plasma Diagnostics
Sensors Plasma Diagnostics Ken Gentle Physics Department Kenneth Gentle RLM 12.330 k.gentle@mail.utexas.edu NRL Formulary MIT Formulary www.psfc.mit.edu/library1/catalog/ reports/2010/11rr/11rr013/11rr013_full.pdf
More informationNeutral beam plasma heating
Seminar I b 1 st year, 2 nd cycle program Neutral beam plasma heating Author: Gabrijela Ikovic Advisor: prof.dr. Tomaž Gyergyek Ljubljana, May 2014 Abstract For plasma to be ignited, external heating is
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationParticle Pinch Model of Passing/Trapped High-Z Impurity with Centrifugal Force Effect )
Particle Pinch Model of Passing/Trapped High-Z Impurity with Centrifugal Force Effect ) Yusuke SHIMIZU, Takaaki FUJITA, Atsushi OKAMOTO, Hideki ARIMOTO, Nobuhiko HAYASHI 1), Kazuo HOSHINO 2), Tomohide
More informationReal Plasma with n, T ~ p Equilibrium: p = j B
Real Plasma with n, T ~ p Equilibrium: p = j B B lines must lie in isobaric surfaces. Since B = 0, only possible if isobaric surfaces are topological tori. Magnetic field lines must form nested tori. Equilibrium
More informationSingle particle motion
Single particle motion Plasma is a collection of a very large number of charged particles moving in, and giving rise to, electromagnetic fields. Before going to the statistical descriptions, let us learn
More informationCLASSICAL ELECTRICITY
CLASSICAL ELECTRICITY AND MAGNETISM by WOLFGANG K. H. PANOFSKY Stanford University and MELBA PHILLIPS Washington University SECOND EDITION ADDISON-WESLEY PUBLISHING COMPANY Reading, Massachusetts Menlo
More informationPLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT. Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford
1 PLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford 2 Overview Why do we need plasmas? For fusion, among other things
More informationMHD RELATED TO 2-FLUID THEORY, KINETIC THEORY AND MAGANETIC RECONNECTION
MHD RELATED TO 2-FLUID THEORY, KINETIC THEORY AND MAGANETIC RECONNECTION Marty Goldman University of Colorado Spring 2017 Physics 5150 Issues 2 How is MHD related to 2-fluid theory Level of MHD depends
More informationNeoclassical transport
Neoclassical transport Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 28 th January 2013 Dr Ben Dudson Magnetic Confinement Fusion (1 of 19) Last time Toroidal devices
More informationIntroduction to Plasma Physics
Introduction to Plasma Physics Hartmut Zohm Max-Planck-Institut für Plasmaphysik 85748 Garching DPG Advanced Physics School The Physics of ITER Bad Honnef, 22.09.2014 A simplistic view on a Fusion Power
More informationELECTROMAGNETIC LIQUID METAL WALL PHENOMENA
ELECTROMAGNETIC LIQUID METAL WALL PHENOMENA BY BOB WOOLLEY 15-19 FEBRUARY 1999 APEX-6 MEETING LIQUID WALLS A sufficiently thick, flowing, liquid first wall and tritium breeding blanket which almost completely
More informationNIMROD simulations of dynamo experiments in cylindrical and spherical geometries. Dalton Schnack, Ivan Khalzov, Fatima Ebrahimi, Cary Forest,
NIMROD simulations of dynamo experiments in cylindrical and spherical geometries Dalton Schnack, Ivan Khalzov, Fatima Ebrahimi, Cary Forest, 1 Introduction Two experiments, Madison Plasma Couette Experiment
More informationMAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT
MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT ABSTRACT A. G. Tarditi and J. V. Shebalin Advanced Space Propulsion Laboratory NASA Johnson Space Center Houston, TX
More informationToroidal confinement devices
Toroidal confinement devices Dr Ben Dudson Department of Physics, University of York, Heslington, York YO10 5DD, UK 24 th January 2014 Dr Ben Dudson Magnetic Confinement Fusion (1 of 20) Last time... Power
More informationwhere the index 2 indicates the particle species (electrons and different ion species if present). With
hapter 3 Magnetohydrodynamics luid equations are probably the most widely used equations for the description of inhomogeneous plasmas. While the phase fluid which is governed by the Boltzmann equation
More informationBurn Stabilization of a Tokamak Power Plant with On-Line Estimations of Energy and Particle Confinement Times
Burn Stabilization of a Tokamak Power Plant with On-Line Estimations of Energy and Particle Confinement Times Javier E. Vitela vitela@nucleares.unam.mx Instituto de Ciencias Nucleares Universidad Nacional
More informationIntroduction. Chapter Plasma: definitions
Chapter 1 Introduction 1.1 Plasma: definitions A plasma is a quasi-neutral gas of charged and neutral particles which exhibits collective behaviour. An equivalent, alternative definition: A plasma is a
More informationDispersive Media, Lecture 7 - Thomas Johnson 1. Waves in plasmas. T. Johnson
2017-02-14 Dispersive Media, Lecture 7 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasmas as a coupled system Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas
More informationMotion of a charged particle in an EM field
Department of physics Seminar - 4 th year Motion of a charged particle in an EM field Author: Lojze Gačnik Mentor: Assoc. Prof. Tomaž Gyergyek Ljubljana, November 2011 Abstract In this paper I will present
More informationDerivation of dynamo current drive in a closed current volume and stable current sustainment in the HIT SI experiment
Derivation of dynamo current drive and stable current sustainment in the HIT SI experiment 1 Derivation of dynamo current drive in a closed current volume and stable current sustainment in the HIT SI experiment
More information12. MHD Approximation.
Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal
More information26. Non-linear effects in plasma
Phys780: Plasma Physics Lecture 26. Non-linear effects. Collisionless shocks.. 1 26. Non-linear effects in plasma Collisionless shocks ([1], p.405-421, [6], p.237-245, 249-254; [4], p.429-440) Collisionless
More informationSources of Magnetic Field I
Sources of Magnetic Field I Physics 2415 Lecture 17 Michael Fowler, UVa Today s Topics Forces between currents Ampère s Law Fields inside wire and solenoid Magnetic Field from a Current in a Long Straight
More information22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 14: Formulation of the Stability Problem
.65, MHD Theory of Fusion Systems Prof. Freidberg Lecture 4: Formulation of the Stability Problem Hierarchy of Formulations of the MHD Stability Problem for Arbitrary 3-D Systems. Linearized equations
More informationHeating and current drive: Radio Frequency
Heating and current drive: Radio Frequency Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 13 th February 2012 Dr Ben Dudson Magnetic Confinement Fusion (1 of 26)
More informationSingle particle motion and trapped particles
Single particle motion and trapped particles Gyromotion of ions and electrons Drifts in electric fields Inhomogeneous magnetic fields Magnetic and general drift motions Trapped magnetospheric particles
More informationRecapitulation: Questions on Chaps. 1 and 2 #A
Recapitulation: Questions on Chaps. 1 and 2 #A Chapter 1. Introduction What is the importance of plasma physics? How are plasmas confined in the laboratory and in nature? Why are plasmas important in astrophysics?
More informationThe Field-Reversed Configuration (FRC) is a high-beta compact toroidal in which the external field is reversed on axis by azimuthal plasma The FRC is
and Stability of Field-Reversed Equilibrium with Toroidal Field Configurations Atomics General Box 85608, San Diego, California 92186-5608 P.O. APS Annual APS Meeting of the Division of Plasma Physics
More informationThe Levitated Dipole Experiment: Towards Fusion Without Tritium
The Levitated Dipole Experiment: Towards Fusion Without Tritium Jay Kesner MIT M.S. Davis, J.E. Ellsworth, D.T. Garnier, M.E. Mauel, P.C. Michael, P.P. Woskov MCP I3.110 Presented at the EPS Meeting, Dublin,
More informationModeling of ELM Dynamics for ITER
Modeling of ELM Dynamics for ITER A.Y. PANKIN 1, G. BATEMAN 1, D.P. BRENNAN 2, A.H. KRITZ 1, S. KRUGER 3, P.B. SNYDER 4 and the NIMROD team 1 Lehigh University, 16 Memorial Drive East, Bethlehem, PA 18015
More informationVII. Hydrodynamic theory of stellar winds
VII. Hydrodynamic theory of stellar winds observations winds exist everywhere in the HRD hydrodynamic theory needed to describe stellar atmospheres with winds Unified Model Atmospheres: - based on the
More informationFast Ion Confinement in the MST Reversed Field Pinch
Fast Ion Connement in the MST Reversed Field Pinch Gennady Fiksel B. Hudson, D.J. Den Hartog, R.M. Magee, R. O'Connell, S.C. Prager MST Team - University of Wisconsin - Madison Center for Magnetic Self-Organization
More informationSingle Particle Motion
Single Particle Motion Overview Electromagnetic fields, Lorentz-force, gyration and guiding center, drifts, adiabatic invariants. Pre-requisites: Energy density of the particle population smaller than
More informationCreation and destruction of magnetic fields
HAO/NCAR July 30 2007 Magnetic fields in the Universe Earth Magnetic field present for 3.5 10 9 years, much longer than Ohmic decay time ( 10 4 years) Strong variability on shorter time scales (10 3 years)
More informationSelf-consistent modeling of ITER with BALDUR integrated predictive modeling code
Self-consistent modeling of ITER with BALDUR integrated predictive modeling code Thawatchai Onjun Sirindhorn International Institute of Technology, Thammasat University, Klong Luang, Pathumthani, 12121,
More informationImpact of neutral atoms on plasma turbulence in the tokamak edge region
Impact of neutral atoms on plasma turbulence in the tokamak edge region C. Wersal P. Ricci, F.D. Halpern, R. Jorge, J. Morales, P. Paruta, F. Riva Theory of Fusion Plasmas Joint Varenna-Lausanne International
More informationDenis Bytschkow. Investigation of MHD mode rotation in the tokamak ASDEX Upgrade
Denis Bytschkow Investigation of MHD mode rotation in the tokamak ASDEX Upgrade IPP 1/341 September, 2010 P H Y S I K - D E P A R T M E N T TECHNISCHE UNIVERSITÄT MÜNCHEN Investigation of MHD mode rotation
More informationThe RFP: Plasma Confinement with a Reversed Twist
The RFP: Plasma Confinement with a Reversed Twist JOHN SARFF Department of Physics University of Wisconsin-Madison Invited Tutorial 1997 Meeting APS DPP Pittsburgh Nov. 19, 1997 A tutorial on the Reversed
More informationSimulation of alpha particle current drive and heating in spherical tokamaks
Simulation of alpha particle current drive and heating in spherical tokamaks R. Farengo 1, M. Zarco 1, H. E. Ferrari 1, 1 Centro Atómico Bariloche and Instituto Balseiro, Argentina. Consejo Nacional de
More informationSW103: Lecture 2. Magnetohydrodynamics and MHD models
SW103: Lecture 2 Magnetohydrodynamics and MHD models Scale sizes in the Solar Terrestrial System: or why we use MagnetoHydroDynamics Sun-Earth distance = 1 Astronomical Unit (AU) 200 R Sun 20,000 R E 1
More information35. RESISTIVE INSTABILITIES: CLOSING REMARKS
35. RESISTIVE INSTABILITIES: CLOSING REMARKS In Section 34 we presented a heuristic discussion of the tearing mode. The tearing mode is centered about x = 0, where F x that B 0 ( ) = k! B = 0. Note that
More informationMHD. Jeff Freidberg MIT
MHD Jeff Freidberg MIT 1 What is MHD MHD stands for magnetohydrodynamics MHD is a simple, self-consistent fluid description of a fusion plasma Its main application involves the macroscopic equilibrium
More information2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1. Waves in plasmas. T. Johnson
2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasma physics Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas Transverse waves
More informationLow Temperature Plasma Technology Laboratory
Low Temperature Plasma Technology Laboratory CENTRAL PEAKING OF MAGNETIZED GAS DISCHARGES Francis F. Chen and Davide Curreli LTP-1210 Oct. 2012 Electrical Engineering Department Los Angeles, California
More informationCalculation of alpha particle redistribution in sawteeth using experimentally reconstructed displacement eigenfunctions
Calculation of alpha particle redistribution in sawteeth using experimentally reconstructed displacement eigenfunctions R. Farengo, H. E. Ferrari,2, M.-C. Firpo 3, P. L. Garcia-Martinez 2,3, A. F. Lifschitz
More informationAP Physics C Mechanics Objectives
AP Physics C Mechanics Objectives I. KINEMATICS A. Motion in One Dimension 1. The relationships among position, velocity and acceleration a. Given a graph of position vs. time, identify or sketch a graph
More informationPHYSICS BASIS FOR THE GASDYNAMIC MIRROR (GDM) FUSION ROCKET. Abstract
PHYSICS BASIS FOR THE GASDYNAMIC MIRROR (GDM) FUSION ROCKET T. Kammash Department of Nuclear Engineering and Radiological Sciences University of Michigan Ann Arbor, MI 48109 (313) 764-0205 W. Emrich Jr.
More informationSingle Particle Motion
Single Particle Motion C ontents Uniform E and B E = - guiding centers Definition of guiding center E gravitation Non Uniform B 'grad B' drift, B B Curvature drift Grad -B drift, B B invariance of µ. Magnetic
More informationCurrent-driven instabilities
Current-driven instabilities Ben Dudson Department of Physics, University of York, Heslington, York YO10 5DD, UK 21 st February 2014 Ben Dudson Magnetic Confinement Fusion (1 of 23) Previously In the last
More information(a) (b) (c) (d) (e) (f) r (minor radius) time. time. Soft X-ray. T_e contours (ECE) r (minor radius) time time
Studies of Spherical Tori, Stellarators and Anisotropic Pressure with M3D 1 L.E. Sugiyama 1), W. Park 2), H.R. Strauss 3), S.R. Hudson 2), D. Stutman 4), X-Z. Tang 2) 1) Massachusetts Institute of Technology,
More informationModule 3: Electromagnetism
Module 3: Electromagnetism Lecture - Magnetic Field Objectives In this lecture you will learn the following Electric current is the source of magnetic field. When a charged particle is placed in an electromagnetic
More informationFundamentals of Magnetic Island Theory in Tokamaks
Fundamentals of Magnetic Island Theory in Tokamaks Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Austin, TX, USA Talk available at http://farside.ph.utexas.edu/talks/talks.html
More informationTHE SHEARED-FLOW STABILIZED Z-PINCH
THE SHEARED-FLOW STABILIZED Z-PINCH U. Shumlak* J. Chadney R.P. Golingo D.J. Den Hartog M.C. Hughes S.D. Knecht W. Lowrie V.S. Lukin B.A. Nelson R.J. Oberto J.L. Rohrbach M.P. Ross and G.V. Vogman *University
More informationFIG. 1. "Flower-like" configuration of filaments used for modelling. Magnetic field values for this configuration can be described analytically. Induc
Ion Motion Modelling within Dynamic Filamentary PF-Pinch Column A. Gaψlkowski 1), A. Pasternak 2), M. Sadowski 2) 1) Institute of Plasma Physics and Laser Microfusion, Warsaw, Poland 2) The Andrzej Soltan
More informationMagnetically Confined Fusion: Transport in the core and in the Scrape- off Layer Bogdan Hnat
Magnetically Confined Fusion: Transport in the core and in the Scrape- off Layer ogdan Hnat Joe Dewhurst, David Higgins, Steve Gallagher, James Robinson and Paula Copil Fusion Reaction H + 3 H 4 He + n
More informationIntroduction to Magnetohydrodynamics (MHD)
Introduction to Magnetohydrodynamics (MHD) Tony Arber University of Warwick 4th SOLARNET Summer School on Solar MHD and Reconnection Aim Derivation of MHD equations from conservation laws Quasi-neutrality
More informationRelativistic magnetohydrodynamics. Abstract
Relativistic magnetohydrodynamics R. D. Hazeltine and S. M. Mahajan Institute for Fusion Studies, The University of Texas, Austin, Texas 78712 (October 19, 2000) Abstract The lowest-order description of
More informationJet Stability: A computational survey
Jet Stability Galway 2008-1 Jet Stability: A computational survey Rony Keppens Centre for Plasma-Astrophysics, K.U.Leuven (Belgium) & FOM-Institute for Plasma Physics Rijnhuizen & Astronomical Institute,
More informationCreation and destruction of magnetic fields
HAO/NCAR July 20 2011 Magnetic fields in the Universe Earth Magnetic field present for 3.5 10 9 years, much longer than Ohmic decay time ( 10 4 years) Strong variability on shorter time scales (10 3 years)
More informationW.A. HOULBERG Oak Ridge National Lab., Oak Ridge, TN USA. M.C. ZARNSTORFF Princeton Plasma Plasma Physics Lab., Princeton, NJ USA
INTRINSICALLY STEADY STATE TOKAMAKS K.C. SHAING, A.Y. AYDEMIR, R.D. HAZELTINE Institute for Fusion Studies, The University of Texas at Austin, Austin TX 78712 USA W.A. HOULBERG Oak Ridge National Lab.,
More informationChapter 4. Electrostatic Fields in Matter
Chapter 4. Electrostatic Fields in Matter 4.1. Polarization 4.2. The Field of a Polarized Object 4.3. The Electric Displacement 4.4. Linear Dielectrics 4.5. Energy in dielectric systems 4.6. Forces on
More informationChapter 8. Conservation Laws. 8.3 Magnetic Forces Do No Work
Chapter 8. Conservation Laws 8.3 Magnetic Forces Do No Work 8.2 Momentum of EM fields 8.2.1 Newton's Third Law in Electrodynamics Consider two charges, q 1 and q 2, moving with speeds v 1 and v 2 magnetic
More informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More informationEquilibrium and transport in Tokamaks
Equilibrium and transport in Tokamaks Jacques Blum Laboratoire J.-A. Dieudonné, Université de Nice Sophia-Antipolis Parc Valrose 06108 Nice Cedex 02, France jblum@unice.fr 08 septembre 2008 Jacques Blum
More information