DOCSLIB.ORG

Fundamentals of Plasma Physics and Controlled Fusion Kenro Miyamoto

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text

Load more

Fundamentals of Plasma Physics and Controlled Fusion Kenro Miyamoto (Received Sep. 18, 2000) NIFS-PROC-48 Oct. 2000 Fundamentals of Plasma Physics and Controlled Fusion by Kenro Miyamoto i Fundamentals of Plasma Physics and Controlled Fusion Kenro Miyamoto Contents Preface 1 Nature of Plasma ....................................................................1 1.1 Introduction .......................................................................... 1 1.2 Charge Neutrality and Landau Damping ..............................................1 1.3 Fusion Core Plasma ...................................................................3 2 Plasma Characteristics ..............................................................7 2.1 Velocity Space Distribution Function, Electron and Ion Temperatures ..................7 2.2 Plasma Frequency, Debye Length ......................................................8 2.3 Cyclotron Frequency, Larmor Radius ..................................................9 2.4 Drift Velocity of Guiding Center ..................................................... 10 2.5 Magnetic Moment, Mirror Confinement, Longitudinal Adiabatic Constant ............12 2.6 Coulomb Collision Time, Fast Neutral Beam Injection ................................14 2.7 Runaway Electron, Dreicer Field .....................................................18 2.8 Electric Resistivity, Ohmic Heating ...................................................19 2.9 Variety of Time and Space Scales in Plasmas .........................................19 3 Magnetic Configuration and Particle Orbit ......................................21 3.1 Maxwell Equations ...................................................................21 3.2 Magnetic Surface .....................................................................23 3.3 Equation of Motion of a Charged Particle ............................................24 3.4 Particle Orbit in Axially Symmetric System ..........................................26 3.5 Drift of Guiding Center in Toroidal Field .............................................28 a Guiding Center of Circulating Particles b Guiding Center of Banana Particles 3.6 Orbit of Guiding Center and Magnetic Surface .......................................31 3.7 Effect of Longitudinal Electric Field on Banana Orbit ................................32 4 Velocity Space Distribution Function and Boltzmann’s Equation .............................................................35 4.1 Phase Space and Distribution Function ...............................................35 4.2 Boltzmann’s Equation and Vlasov’s Equation ........................................ 36 5 Plasma as MHD Fluid ..............................................................39 5.1 Magnetohydrodynamic Equations for Two Fluids .....................................39 5.2 Magnetogydrodynamic Equations for One Fluid ......................................41 5.3 Simplified Magnetohydrodynamic Equations ..........................................42 5.4 Magnetoacoustic Wave ...............................................................44 6 Equilibrium ..........................................................................47 6.1 Pressure Equilibrium .................................................................47 6.2 Equilibrium Equation for Axially Symmetric and Translationally Symmetric Systems . 48 ii Contents 6.3 Tokamak Equilibrium ................................................................50 6.4 Poloidal Field for Tokamak Equilibrium ..............................................55 6.5 Upper Limit of Beta Ratio ...........................................................58 6.6 Pfirsch-Schl¨uter Current ..............................................................59 6.7 Virial Theorem .......................................................................61 7 Diffusion of Plasma, Confinement Time ..........................................65 7.1 Collisional Diffusion (Classical Diffusion) .............................................66 a Magnetohydrodynamic Treatment b A Particle Model 7.2 Neoclassical Diffusion of Electrons in Tokamak .......................................69 7.3 Fluctuation Loss, Bohm Diffusion, and Stationary Convective Loss ...................71 7.4 Loss by Magnetic Fluctuation ........................................................75 8 Magnetohydrodynamic Instabilities ...............................................77 8.1 Interchange, Sausage and Kink Instabilities ...........................................77 a Interchange Instability b Stability Criterion for Interchange Instability, Magnetic Well c Sausage Instability d Kink Instability 8.2 Formulation of Magnetohydrodynamic Instabilities ...................................85 a Linearization of Magnetohydrodynamic Equations b Energy Principle 8.3 Instabilities of a Cylindrical Plasma ..................................................90 a Instabilities of Sharp-Boundary Configuration: Kruskal-Shafranov Condition b Instabilities of Diffuse-Boundary Configurations c Suydam’s Criterion d Tokamak Configuration e Reversed Field Pinch 8.4 Hain-L¨ust Magnetohydrodynamic Equation .........................................102 8.5 Ballooning Instability ...............................................................103 8.6 ηi Mode due to Density and Temperature Gradient .................................. 107 9 Resistive Instability ............................................................... 111 9.1 Tearing Instability ..................................................................111 9.2 Resistive Drift Instability ...........................................................115 10 Plasma as Medium of Waves ....................................................121 10.1 Dispersion Equation of Waves in a Cold Plasma ....................................121 10.2 Properties of Waves ................................................................124 a Polarization and Particle Motion b Cutoff and Resonance 10.3 Waves in a Two-Components Plasma .............................................. 126 10.4 Various Waves .....................................................................129 a Alfv´en Wave b Ion Cyclotron Wave and Fast Wave c Lower Hybrid Resonance d Upper Hybrid Resonance e Electron Cyclotron Wave 10.5 Conditions for Electrostatic Waves .................................................134 11 Landau Damping and Cyclotron Damping .....................................137 11.1 Landau Damping (Amplification) ..................................................137 11.2 Transit-Time Damping .............................................................140 11.3 Cyclotron Damping ................................................................140 11.4 Quasi-Linear Theory of Evolution in the Distribution Function .....................142 12 Wave Propagation and Wave Heating ..........................................145 12.1 Energy Flow .......................................................................145 Contents iii 12.2 Ray Tracing .......................................................................148 12.3 Dielectric Tensor of Hot Plasma, Wave Absorption and Heating ....................149 12.4 Wave Heating in Ion Cyclotron Range of Frequency ................................154 12.5 Lower Hybrid Wave Heating .......................................................157 12.6 Electron Cyclotron Heating ........................................................159 13 Velocity Space Instabilities (Electrostatic Waves) .............................163 13.1 Dispersion Equation of Electrostatic Wave .........................................163 13.2 Two Streams Instability ........................................................... 164 13.3 Electron Beam Instability ..........................................................165 13.4 Harris Instability ..................................................................165 14 Instabilities driven by Energetic Particles ......................................169 14.1 Fishbone Instability ................................................................169 a Formulation b MHD Potential Energy c Kinetic Integral of Hot Component d Growth Rate of Fishbone Instability 14.2 Toroidal Alfv´en Eigenmodes .......................................................176 a Toroidicity Inuced Alfv´en Eigenmode b Instabilities of TAE Driven by Energetic Particles c Various Alfv´en Modes 15 Development of Fusion Researches .............................................189 16 Tokamak ...........................................................................197 16.1 Tokamak Devices ..................................................................197 16.2 Equilibrium ........................................................................201 a Case with Conducting Shell b Case without Conducting Shell c Equilibrium Beta Limit of Tokamaks with Elongated Plasma Cross Section 16.3 MHD Stability and Density Limit ..................................................203 16.4 Beta Limit of Elongated Plasma ...................................................205 16.5 Impurity Control, Scrape-Off Layer and Divertor ...................................206 16.6 Confinement Scaling of L Mode ....................................................210 16.7 H Mode and Improved Confinement Modes ........................................212 16.8 Noninductive Current Drive ........................................................217 a Lower Hybrid Current Drive b Electron Cyclotron Current Drive c Neutral Beam Current Drive d Bootstrap Current 16.9
Recommended publications
  • Analysis of Plasma Instabilities and Verification of the BOUT Code for The

    Analysis of Plasma Instabilities and Verification of the BOUT Code for The

    Analysis of plasma instabilities and verification of the BOUT code for the Large Plasma Device P. Popovich,1 M.V. Umansky,2 T.A. Carter,1, a) and B. Friedman1 1)Department of Physics and Astronomy and Center for Multiscale Plasma Dynamics, University of California, Los Angeles, CA 90095-1547 2)Lawrence Livermore National Laboratory, Livermore, CA 94550, USA (Dated: 8 August 2018) The properties of linear instabilities in the Large Plasma Device [W. Gekelman et al., Rev. Sci. Inst., 62, 2875 (1991)] are studied both through analytic calculations and solving numerically a system of linearized collisional plasma fluid equations using the 3D fluid code BOUT [M. Umansky et al., Contrib. Plasma Phys. 180, 887 (2009)], which has been successfully modified to treat cylindrical geometry. Instability drive from plasma pressure gradients and flows is considered, focusing on resistive drift waves, the Kelvin-Helmholtz and rotational interchange instabilities. A general linear dispersion relation for partially ionized collisional plasmas including these modes is derived and analyzed. For LAPD relevant profiles including strongly driven flows it is found that all three modes can have comparable growth rates and frequencies. Detailed comparison with solutions of the analytic dispersion relation demonstrates that BOUT accurately reproduces all characteristics of linear modes in this system. PACS numbers: 52.30.Ex, 52.35.Fp, 52.35.Kt, 52.35.Lv, 52.65.Kj arXiv:1004.4674v4 [physics.plasm-ph] 26 Jan 2011 a)Electronic mail: [email protected] 1 I. INTRODUCTION Understanding complex nonlinear phenomena in magnetized plasmas increasingly relies on the use of numerical simulation as an enabling tool.
  • Stellarators. Present Status and Future Planning

    Stellarators. Present Status and Future Planning

    STELLARATORS Present Status and Future Planning H. Wobig, Garching IMax-Planck-Institut für Plasmaphysikl magnetic axis Fig. 1 — Period of toroidally dosed magnetic surfaces. The stellarator is a closed toroidal device open-ended configurations like mirror ma­ designed to confine a hot plasma In a ma­ chines (see EN, 12 8/9) there are always gnetic field. It is one of the oldest concepts plasma particles which escape the confine­ magnetic surfaces is investigated by field to have been investigated in the search for ment volume along the magnetic field and line integration. The topology of the field controlled thermonuclear fusion. an isotropic distribution function cannot be lines is determined by the rotational trans­ The basic idea of the stellarator was maintained. If this anisotropy is too strong form or twist number 1/27t (or+), which is developed by Lyman Spitzer, Professor of it gives rise to instabilities and enhanced the number of revolutions of a field line Astronomy at Princeton in 19511). Experi­ plasma losses. It is mainly to avoid these around the magnetic axis during one toroi­ mental studies began in 1952 and after the disadvantages that toroidal confinement dal revolution. declassification of fusion research in 1958 has been preferred. The need for a helical field can be under­ the idea was picked up by other research In toroidal configurations, the currents stood from a look at the particle orbits. groups. In Europe, the first stellarators which generate the confining magnetic Charged particles tend to follow the field were built in the Max-Planck Institute for fields can be classified under three lines and unless these are helical, because Physics and Astrophysics in Munich and categories : of inhomogeneities, the particles will drift later at Culham, Moscow and Karkhov.
  • 10. Collisions • Use Conservation of Momentum and Energy and The

    10. Collisions • Use Conservation of Momentum and Energy and The

    10. Collisions • Use conservation of momentum and energy and the center of mass to understand collisions between two objects. • During a collision, two or more objects exert a force on one another for a short time: -F(t) F(t) Before During After • It is not necessary for the objects to touch during a collision, e.g. an asteroid flied by the earth is considered a collision because its path is changed due to the gravitational attraction of the earth. One can still use conservation of momentum and energy to analyze the collision. Impulse: During a collision, the objects exert a force on one another. This force may be complicated and change with time. However, from Newton's 3rd Law, the two objects must exert an equal and opposite force on one another. F(t) t ti tf Dt From Newton'sr 2nd Law: dp r = F (t) dt r r dp = F (t)dt r r r r tf p f - pi = Dp = ò F (t)dt ti The change in the momentum is defined as the impulse of the collision. • Impulse is a vector quantity. Impulse-Linear Momentum Theorem: In a collision, the impulse on an object is equal to the change in momentum: r r J = Dp Conservation of Linear Momentum: In a system of two or more particles that are colliding, the forces that these objects exert on one another are internal forces. These internal forces cannot change the momentum of the system. Only an external force can change the momentum. The linear momentum of a closed isolated system is conserved during a collision of objects within the system.
  • Controlled Nuclear Fusion

    Controlled Nuclear Fusion

    Controlled Nuclear Fusion HANNAH SILVER, SPENCER LUKE, PETER TING, ADAM BARRETT, TORY TILTON, GABE KARP, TIMOTHY BERWIND Nuclear Fusion Thermonuclear fusion is the process by which nuclei of low atomic weight such as hydrogen combine to form nuclei of higher atomic weight such as helium. two isotopes of hydrogen, deuterium (composed of a hydrogen nucleus containing one neutrons and one proton) and tritium (a hydrogen nucleus containing two neutrons and one proton), provide the most energetically favorable fusion reactants. in the fusion process, some of the mass of the original nuclei is lost and transformed to energy in the form of high-energy particles. energy from fusion reactions is the most basic form of energy in the universe; our sun and all other stars produce energy through thermonuclear fusion reactions. Nuclear Fusion Overview Two nuclei fuse together to form one larger nucleus Fusion occurs in the sun, supernovae explosion, and right after the big bang Occurs in the stars Initially, research failed Nuclear weapon research renewed interest The Science of Nuclear Fusion Fusion in stars is mostly of hydrogen (H1 & H2) Electrically charged hydrogen atoms repel each other. The heat from stars speeds up hydrogen atoms Nuclei move so fast, they push through the repulsive electric force Reaction creates radiant & thermal energy Controlled Fusion uses two main elements Deuterium is found in sea water and can be extracted using sea water Tritium can be made from lithium When the thermal energy output exceeds input, the equation is self-sustaining and called a thermonuclear reaction 1929 1939 1954 1976 1988 1993 2003 Prediction Quantitativ ZETA JET Project Japanese Princeton ITER using e=mc2, e theory Tokomak Generates that energy explaining 10 from fusion is fusion.
  • Inertial Electrostatic Confinement Fusor Cody Boyd Virginia Commonwealth University

    Inertial Electrostatic Confinement Fusor Cody Boyd Virginia Commonwealth University

    Virginia Commonwealth University VCU Scholars Compass Capstone Design Expo Posters College of Engineering 2015 Inertial Electrostatic Confinement Fusor Cody Boyd Virginia Commonwealth University Brian Hortelano Virginia Commonwealth University Yonathan Kassaye Virginia Commonwealth University See next page for additional authors Follow this and additional works at: https://scholarscompass.vcu.edu/capstone Part of the Engineering Commons © The Author(s) Downloaded from https://scholarscompass.vcu.edu/capstone/40 This Poster is brought to you for free and open access by the College of Engineering at VCU Scholars Compass. It has been accepted for inclusion in Capstone Design Expo Posters by an authorized administrator of VCU Scholars Compass. For more information, please contact [email protected]. Authors Cody Boyd, Brian Hortelano, Yonathan Kassaye, Dimitris Killinger, Adam Stanfield, Jordan Stark, Thomas Veilleux, and Nick Reuter This poster is available at VCU Scholars Compass: https://scholarscompass.vcu.edu/capstone/40 Team Members: Cody Boyd, Brian Hortelano, Yonathan Kassaye, Dimitris Killinger, Adam Stanfield, Jordan Stark, Thomas Veilleux Inertial Electrostatic Faculty Advisor: Dr. Sama Bilbao Y Leon, Mr. James G. Miller Sponsor: Confinement Fusor Dominion Virginia Power What is Fusion? Shielding Computational Modeling Because the D-D fusion reaction One of the potential uses of the fusor will be to results in the production of neutrons irradiate materials and see how they behave after and X-rays, shielding is necessary to certain levels of both fast and thermal neutron protect users from the radiation exposure. To reduce the amount of time and produced by the fusor. A Monte Carlo resources spent testing, a computational model n-Particle (MCNP) model was using XOOPIC, a particle interaction software, developed to calculate the necessary was developed to model the fusor.
  • Law of Conversation of Energy

    Law of Conversation of Energy

    Law of Conservation of Mass: "In any kind of physical or chemical process, mass is neither created nor destroyed - the mass before the process equals the mass after the process." - the total mass of the system does not change, the total mass of the products of a chemical reaction is always the same as the total mass of the original materials. "Physics for scientists and engineers," 4th edition, Vol.1, Raymond A. Serway, Saunders College Publishing, 1996. Ex. 1) When wood burns, mass seems to disappear because some of the products of reaction are gases; if the mass of the original wood is added to the mass of the oxygen that combined with it and if the mass of the resulting ash is added to the mass o the gaseous products, the two sums will turn out exactly equal. 2) Iron increases in weight on rusting because it combines with gases from the air, and the increase in weight is exactly equal to the weight of gas consumed. Out of thousands of reactions that have been tested with accurate chemical balances, no deviation from the law has ever been found. Law of Conversation of Energy: The total energy of a closed system is constant. Matter is neither created nor destroyed – total mass of reactants equals total mass of products You can calculate the change of temp by simply understanding that energy and the mass is conserved - it means that we added the two heat quantities together we can calculate the change of temperature by using the law or measure change of temp and show the conservation of energy E1 + E2 = E3 -> E(universe) = E(System) + E(Surroundings) M1 + M2 = M3 Is T1 + T2 = unknown (No, no law of conservation of temperature, so we have to use the concept of conservation of energy) Total amount of thermal energy in beaker of water in absolute terms as opposed to differential terms (reference point is 0 degrees Kelvin) Knowns: M1, M2, T1, T2 (Kelvin) When add the two together, want to know what T3 and M3 are going to be.
  • NEWSLETTER Issue No. 7 September 2017

    NEWSLETTER Issue No. 7 September 2017

    NEWSLETTER September 2017 Issue no. 7 Nuclear Industry Group Newsletter September 2017 Contents Notes from the Chair ................................................................................... 3 IOP Group Officers Forum .......................................................................... 4 NIG Committee Elections ............................................................................ 6 Nuclear Industry Group Career Contribution Prize 2017 .......................... 7 Event – Gen IV Reactors by Richard Stainsby (NNL) ................................ 8 Event – Nuclear Security by Robert Rodger (NNL) and Graham Urwin (RWM) ......................................................................................................... 12 Event – The UK’s Nuclear Future by Dame Sue Ion ................................ 13 Event – Regulatory Challenges for Nuclear New Build by Mike Finnerty. .................................................................................................................... 15 Event – European Nuclear Young Generation Forum ............................. 18 Event – Nuclear Fusion, 60 Years on from ZETA by Chris Warrick (UKAEA), Kate Lancaster (York Plasma Institute), David Kingham (Tokamak Energy) and Ian Chapman (UKAEA) ....................................... 19 IOP Materials and Characterisation Group Meetings .............................. 25 “Brexatom” – the implications of the withdrawal for the UK from the Euratom Treaty. ........................................................................................
  • Shear Flow-Interchange Instability in Nightside Magnetotail Causes Auroral Beads As a Signature of Substorm Onset

    Shear Flow-Interchange Instability in Nightside Magnetotail Causes Auroral Beads As a Signature of Substorm Onset

    manuscript submitted to JGR: Space Physics Shear flow-interchange instability in nightside magnetotail causes auroral beads as a signature of substorm onset Jason Derr1, Wendell Horton1, Richard Wolf2 1Institute for Fusion Studies and Applied Research Laboratories, University of Texas at Austin 2Physics and Astronomy Department, Rice University Key Points: • A wedge model is extended to obtain a MHD magnetospheric wave equation for the near-Earth magnetotail, including velocity shear effects. • Shear flow-interchange instability to replace shear flow-ballooning and interchange in- stabilities as substorm onset cause. • WKB applicability conditions suggest nonlinear analysis is necessary and yields spatial scale of most unstable mode’s variations. arXiv:1904.11056v1 [physics.space-ph] 24 Apr 2019 Corresponding author: Jason Derr, [email protected] –1– manuscript submitted to JGR: Space Physics Abstract A geometric wedge model of the near-earth nightside plasma sheet is used to derive a wave equation for low frequency shear flow-interchange waves which transmit E~ × B~ sheared zonal flows along magnetic flux tubes towards the ionosphere. Discrepancies with the wave equation result used in Kalmoni et al. (2015) for shear flow-ballooning instability are dis- cussed. The shear flow-interchange instability appears to be responsible for substorm onset. The wedge wave equation is used to compute rough expressions for dispersion relations and local growth rates in the midnight region of the nightside magnetotail where the instability develops, forming the auroral beads characteristic of geomagnetic substorm onset. Stability analysis for the shear flow-interchange modes demonstrates that nonlinear analysis is neces- sary for quantitatively accurate results and determines the spatial scale on which the instability varies.
  • Plasma Stability in a Dipole Magnetic Field Andrei N. Simakov MS Physics, Moscow Institute of Physics and Technology, Russia

    Plasma Stability in a Dipole Magnetic Field Andrei N. Simakov MS Physics, Moscow Institute of Physics and Technology, Russia

    Plasma Stability in a Dipole Magnetic Field by Andrei N. Simakov M.S. Physics, Moscow Institute of Physics and Technology, Russia, 1997 Submitted to the Department of Physics in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Physics at the MASSACHUSETTS INSTITUTE OF TECHNOLOGY June 2001 @ Massachusetts Institute of Technology 2001. All ri0h FTECNOOG 8cience JUN 2 2 2001 LIBRARIES Author................................................. "-partment of Physics April 26, 2001 C ertified by ..... ................... Peter J. Catto Senior Research Scientist, Plasma Science and Fusion Center Thesis Supervisor Certified by Miklos Porkolab Professor, Department of Physics Thesis Supervisor Accepted by ........ .......... Thomas J. Greytak Professor, Associate Depart nt Head for Education Ii 2 Plasma Stability in a Dipole Magnetic Field by Andrei N. Simakov M.S. Physics, Moscow Institute of Physics and Technology, Russia, 1997 Submitted to the Department of Physics on April 26, 2001, in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Physics Abstract The MHD and kinetic stability of an axially symmetric plasma, confined by a poloidal magnetic field with closed lines, is considered. In such a system the stabilizing effects of plasma compression and magnetic field compression counteract the unfavorable field line curvature and can stabilize pressure gradient driven magnetohydrodynamic modes provided the pressure gradient is not too steep. Isotropic pressure, ideal MHD stability is studied first and a general interchange stability condition and an integro-differential eigenmode equation for ballooning modes are derived, using the MHD energy principle. The existence of plasma equilibria which are both interchange and ballooning stable for arbitrarily large beta = plasma pressure / magnetic pressure, is demonstrated.
  • [Insert Your Title Here]

    [Insert Your Title Here]

    Laboratory Investigation of the Dynamics of Shear Flows in a Plasma Boundary Layer by Ami Marie DuBois A dissertation submitted to the Graduate Faculty of Auburn University in partial fulfillment of the requirements for the Degree of Doctor of Philosophy Auburn, Alabama December 14, 2013 Copyright 2013 by Ami Marie DuBois Approved by Edward Thomas, Jr., Chair, Professor of Physics William Amatucci, Staff Scientist at Naval Research Laboratory Uwe Konopka, Professor of Physics David Maurer, Professor of Physics Abstract The laboratory experiments presented in this dissertation investigate a regime of instabilities that occur when a highly localized, radial electric field oriented perpendicular to a uniform background magnetic field gives rise to an azimuthal velocity shear profile at the boundary between two interpenetrating plasmas. This investigation is motivated by theoretical predictions which state that plasmas are unstable to transverse and parallel inhomogeneous sheared flows over a very broad frequency range. Shear driven instabilities are commonly observed in the near-Earth space environment when boundary layers, such as the magnetopause and the plasma sheet boundary layer, are compressed by intense solar storms. When the shear scale length is much less than the ion gyro- radius, but greater than the electron gyro-radius, the electrons are magnetized in the shear layer, but the ions are effectively un-magnetized. The resulting shear driven instability, the electron-ion hybrid instability, is investigated in a new interpenetrating plasma configuration in the Auburn Linear EXperiment for Instability Studied (ALEXIS) in the absence of a magnetic field aligned current. In order to truly understand the dynamics at magnetospheric boundary layers, the EIH instability is studied in the presence of a density gradient located at the boundary layer between two plasmas.
  • Highlights in Early Stellarator Research at Princeton

    Highlights in Early Stellarator Research at Princeton

    J. Plasma Fusion Res. SERIES, Vol.1 (1998) 3-8 Highlights in Early Stellarator Research at Princeton STIX Thomas H. Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08540, USA (Received: 30 September 1997/Accepted: 22 October 1997) Abstract This paper presents an overview of the work on Stellarators in Princeton during the first fifteen years. Particular emphasis is given to the pioneering contributions of the late Lyman Spitzer, Jr. The concepts discussed will include equilibrium, stability, ohmic and radiofrequency plasma heating, plasma purity, and the problems associated with creating a full-scale fusion power plant. Brief descriptions are given of the early Princeton Stellarators: Model A, Model B, Model B-2, Model B-3, Models 8-64 and 8-65, and Model C, and also of the postulated fusion power plant, Model D. Keywords: Spitzer, Kruskal, stellarator, rotational transform, Bohm diffusion, ohmic heating, magnetic pumping, ion cyclotron resonance heating (ICRH), magnetic island, tokamak On March 31 of this year, at the age of 82, Lyman stellarator was brought to the headquarters of the U.S. Spitzer, Jr., a true pioneer in the fields of astrophysics Atomic Energy Commission in Washington where it re- and plasma physics, died. I wish to dedicate this presen- ceived a favorable reception. Spitzer chose the name tation to his memory. "Project Matterhorn" for the project which was to be Forty-six years ago, in early 1951, Spitzer, then sited in the Princeton area, on the newly acquired For- chair of the Department of Astronomy at Princeton restal tract, and funding began on July 1 of that year University, together with Princeton physicist John 121- Wheeler, had been thinking about the physics of ther- Spitzer's earliest stellarator papers comprise a truly monoclear processes.
  • Heat and Energy Conservation

    Heat and Energy Conservation

    1 Lecture notes in Fluid Dynamics (1.63J/2.01J) by Chiang C. Mei, MIT, Spring, 2007 CHAPTER 4. THERMAL EFFECTS IN FLUIDS 4-1-2energy.tex 4.1 Heat and energy conservation Recall the basic equations for a compressible fluid. Mass conservation requires that : ρt + ∇ · ρ~q = 0 (4.1.1) Momentum conservation requires that : = ρ (~qt + ~q∇ · ~q)= −∇p + ∇· τ +ρf~ (4.1.2) = where the viscous stress tensor τ has the components = ∂qi ∂qi ∂qk τ = τij = µ + + λ δij ij ∂xj ∂xi ! ∂xk There are 5 unknowns ρ, p, qi but only 4 equations. One more equation is needed. 4.1.1 Conservation of total energy Consider both mechanical ad thermal energy. Let e be the internal (thermal) energy per unit mass due to microscopic motion, and q2/2 be the kinetic energy per unit mass due to macroscopic motion. Conservation of energy requires D q2 ρ e + dV rate of incr. of energy in V (t) Dt ZZZV 2 ! = − Q~ · ~ndS rate of heat flux into V ZZS + ρf~ · ~qdV rate of work by body force ZZZV + Σ~ · ~qdS rate of work by surface force ZZX Use the kinematic transport theorm, the left hand side becomes D q2 ρ e + dV ZZZV Dt 2 ! 2 Using Gauss theorem the heat flux term becomes ∂Qi − QinidS = − dV ZZS ZZZV ∂xi The work done by surface stress becomes Σjqj dS = (σjini)qj dS ZZS ZZS ∂(σijqj) = (σijqj)ni dS = dV ZZS ZZZV ∂xi Now all terms are expressed as volume integrals over an arbitrary material volume, the following must be true at every point in space, 2 D q ∂Qi ∂(σijqi) ρ e + = − + ρfiqi + (4.1.3) Dt 2 ! ∂xi ∂xj As an alternative form, we differentiate the kinetic energy and get De