Ideal MHD Equilibria
|
|
- Brianna Bennett
- 5 years ago
- Views:
Transcription
1 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 equivalent stationary MHD equations governed by (ρv) = 0 ρ (v ) v = p + ( B) B + ρg (v ) p = γp v (v B) = 0 and B = 0 given g and parameter γ, solve for ρ, v, p, B
2 CapSel Equil - 02 dramatic simplication when considering STATIC equilibria v = 0 leaves only p + ( B) B + ρg = 0 and B = 0 governing equations for magnetostatic equilibria without external gravity g = 0 p }{{} pressure gradient = ( B) B } {{ } Lorentz force B = 0
3 CapSel Equil - 03a magnetostatic case without gravity: choose dimensionality and geometry 1D cartesian geometry p(x), B x (x), B y (x), B z (x) constant (or zero) B x from B = 0 d p + B2 y + B 2 z dx 2 } {{ } total pressure = 0 x z B y
4 CapSel 1D cylindrical geometry p(r), B R (R), B ϕ (R), B Z (R) (constant RB R hence) zero B R from B = 0 d p + B2 ϕ + B 2 Z = B2 ϕ dr 2 Z ϕ } {{ } total pressure R }{{} radially inward tension Equil - 03b R B ODE governs 1D static equilibria: 4 functions ρ, p, B 2, B 3 3 free profiles, namely ρ and 2 from p, B 2, B 3
5 CapSel Equil D cylindrical case: dimensions from cylinder radius a, axial field strength B 0, axial density ρ 0 2 essential dimensionless parameters β 0 = 2p 0 B 2 0 and µ 0 = ab ϕ RB Z R=0 = aj Z0 2B 0 axial plasma beta β 0 inverse pitch µ = ab ϕ /RB Z Latter uses axial current J Z0 = J Z (R = 0) = 1 R d dr (RB ϕ) R=0
6 CapSel Equil - 05 introduce dimensionless quantities f and unit profiles f R R a R p p B 2 0 ρ ρ ρ 0 ρ = β 0 2 p p 0 β 0 2 p B Z B Z B 0 B Z B ϕ B ϕ B 0 µ 0 µ R B Z
7 CapSel Equil D cylindrical equilibrium fully specified by β 0, µ 0 and two (three with ρ) unit profiles p( R) and µ = µ/µ 0 determine B Z from equilibrium relation d d R β 0 B Z 2 2 p + 2 ( R 2 µ 2 0 µ ) find B ϕ R µ B Z and B ϕ = µ 0 Bϕ BZ 2 = Rµ 2 0 µ 2 core problem, freedom in parameter values and unit profile variations for cylinder of finite length L 2πR 0 : alternatively introduce ɛ = a/r 0 and work with q 0 = ɛ/µ 0 instead q-profile from q = RB Z /R 0 B ϕ
8 CapSel Grad-Shafranov equation Equil - 07a consider axisymmetric ( ϕ = 0), 2D static equilibria assume right-handed coordinate system (R, Z, ϕ) to solve for 5 two-dimensional functions ρ(r, Z), p(r, Z), B R (R, Z), B Z (R, Z), B ϕ (R, Z) density profile arbitrary, not in equilibrium p = ( B) B from B = 0 or 1 R R (RB R) + B Z Z = 0 flux function ψ(r, Z) with B = 1 R eˆ ϕ ψ + B ϕ eˆ ϕ such that B R = 1 R ψ Z and B Z = 1 ψ R R
9 CapSel Equil - 07b note that B ψ = 0 constant ψ surfaces (or contours in (R, Z)-plane) are flux surfaces field lines lie on these flux surfaces 3 orthogonal directions from p = ( B) B Poloidal R-Z plane ϕ J J p B
10 CapSel Equil - 08 since J = B we have J = 0 analogously define current stream function I(R, Z) J = 1 R eˆ ϕ I + J ϕ eˆ ϕ this yields J R = 1 I R Z and J Z = 1 I R R from J = B we have as well ϕ component of p = ( B) B J R = B ϕ Z I = RB ϕ 0 = J R B Z J Z B R = 1 I R 2 Z since commutator vanishes, find I(ψ) ψ R I ψ R Z
11 CapSel Equil - 09 R and Z component of force balance yield p R = J ZB ϕ J ϕ B Z = p Z = J ϕb R J R B ϕ = again find p(ψ) and conclude II R J ϕ 2 R II p = II R J ϕ 2 R R J ϕ 2 R ψ R ψ Z
12 CapSel Equil - 10a use J = B to find J ϕ = B Z R B R Z in terms of ψ we get RJ ϕ = R R Grad-Shafranov equation thus reads R R ( 1 R ψ R ) + 2 ψ Z 2 1 ψ + 2 ψ R R Z = 2 II p R 2 to solve for ψ(r, Z) under given profiles p(ψ) and I(ψ)
13 CapSel Equil - 10b GS: 2nd order PDE character from aψ RR + 2bψ RZ + cψ ZZ +... find b 2 ac = 1 < 0: elliptic Boundary value problem well-posed for elliptic Laplace equation: solution reaches extremum on boundary here we have (different from Laplace ψ = 0): R 2 1 R 2 ψ = II p R 2 normally have ψ [0, ψ 1 ]: zero at magnetic axis, fixed BV ψ 1
14 CapSel Equil - 11 Grad-Shafranov for tokamak equilibria tokamak = toroidal vessel for achieving controlled nuclear fusion task to achieve steady state magnetically confined plasma within poloidal cross-section of vessel: ideal MHD equilibrium static plasma and axisymmetry: solve GS within given boundary overall configuration: donut with nested flux surfaces Z δ ϕ R R 0 2a
15 CapSel Equil - 12 Scaling for Grad-Shafranov natural to set units by a (length) and B 0 (magnetic and pressure) from 2a horizontal diameter of vessel geometric center of vessel at R 0 ɛ a/r 0 inverse aspect ratio use strength of vacuum field at R 0 for B 0 value of flux on outer boundary ψ 1 use to define a unit flux label ψ = ψ/ψ1 define dimensionless inverse flux α a 2 B 0 /ψ 1 change to dimensionless poloidal coordinates (x, y) R Z x = R R 0 a = R 1 ɛ y = Z = Z a
16 CapSel Equil - 13 LHS of GS changes to R R 1 ψ + 2 ψ ψ 1 R R Z 2 a 2 2 ψ x 2 ɛ ψ 1 + ɛx x + 2 ψ y 2 RHS of GS becomes I di dψ R2 dp dψ a2 B 2 0 ψ 1 Ĩ dĩ 2 d p + R d ψ d ψ using Ĩ RB ϕ/ab 0 and p p/b 2 0
17 CapSel Equil - 14 introduce dimensionless profiles scaled pressure P = P ( ψ) α2 ɛ p Q = Q( ψ) ɛα2 G = G( ψ) 1 ɛ 2a 2 B 2 0 [ I 2 (ψ) R 2 0B 2 0 ] [ Q( ψ) P ( ψ) ] under these definitions, GS becomes 2 ψ x 2 ɛ ψ 1 + ɛx x + 2 ψ y = dg dp x(2 + ɛx) 2 d ψ d ψ two arbitrary profiles G dg/d ψ and P dp/d ψ
18 CapSel Equil - 15 core problem now identified, separate freedom in G and P magnitude and shape: amplitude and unit profile G AΓ( ψ) with Γ(0) = 1 P 1 2 ABΠ( ψ) with Π(0) = 1 roughly related to current profile and pressure gradient profile Final form of GS is then 2 ψ x 2 boundary conditions ψ = 1 on boundary ɛ ψ 1 + ɛx x + 2 ψ y = A 2 ψ = ψ x = ψ y = 0 at magnetic axis (δ, 0) Γ( ψ) + Bx(1 + ɛx 2 2nd order PDE with 4 BCs: overdetermined!!! )Π( ψ)
19 CapSel Equil - 16 A (and B if input δ) to be determined with solution ( eigenvalues ) once solution ψ, A, and B identified reconstruct pressure from p = ɛ α 2 ψ magnetic field components from 1 2 ABΠ d ψ + P (0) B R = 1 ɛ ψ α 1 + ɛx y B Z = 1 ɛ ψ α1 + ɛx x B ϕ = ɛ ψ A(ɛΓ 1 1/2 BΠ)d ψ 1 + ɛx α 2 2 different realization of same core problem solution for different α!
20 CapSel Equil - 17 The Soloviev solution analytic solution to Grad-Shafranov R R 1 ψ + 2 ψ R R Z = 2 II p R 2 for linear profiles I 2 /2 = I 2 0/2 Eψ and p = p 0 F ψ E, F free profile parameters polynomial solution derived by Soloviev ψ = ( C + DR 2) two additional parameters C, D verify by insertion [ E + ( F 8D 2 ) R 2] Z 2
21 CapSel Equil - 18 use Soloviev for creating tokamak equilibria instructive to change to ψ and (x, y) coordinates Soloviev solution rewritten as ɛ ψ = x 2 (1 2 x2 ) + 1 ɛ2 [1 + τɛx(2 + ɛx)] y2 4 σ 2 unit contour through 4 points (x = ±1, y = 0) (0, y = ±σ) ɛ = a/r 0 inverse aspect ratio wrt geometric center τ measures triangularity (ellips-like for τ = 0) ψ = 0 at magnetic axis (δ, 0) with δ = 1 ɛ [ 1 + ɛ2 1 ] outward Shafranov shift of magnetic axis geometric effect required for equilibrium in torus
22 CapSel Equil - 19 did not consider boundary as given here determined as unit contour from solution cases ɛ = 0.333, τ = 0 and σ = (0.5, 1, 1.5) elliptic flux contours
23 CapSel Equil - 20 varying the triangularity parameter τ case τ = 0.5 versus τ = 0 and τ = 1.5 triangular flux contours, sign of τ for direction
24 CapSel Equil - 21 Force-Free equilibria So far: 1D and 2D static ideal MHD equilibria (no gravity, no flow) pressure gradient balancing Lorentz force in low beta plasma β 1 only solve ( B) B = 0 plus solenoidal constraint current J = B must be B force-free magnetic configuration, can be 1D, 2D, or 3D
25 CapSel Equil - 22 simplest case: vanishing current J = 0 or potential fields vacuum magnetic field solution since B = 0 define φ from B = φ since B = 0 solve 2 φ = 0 elliptic Laplace equation, BVP well-posed frequently used for reconstructing coronal magnetic field normal B component from magnetograms of solar photosphere B n photosphere = φ n photosphere given
26 CapSel more general: field-aligned currents J = αb from ( B) = 0 find B α = 0 Equil - 23a α is constant along field line 1D examples are slab: B y (x) = sin(αx) and B z (x) = cos(αx) cylinder: B Z (R) = J 0 (αr) and B ϕ (R) = J 1 (αr)
27 CapSel Equil - 23b example taken from Régnier & Amari three reconstructions: potential, constant α, varying α best agreement with observations for nonlinear force-free case
28 CapSel Equil - 24 Stationary equilibria static (v = 0) case only left momentum balance stationary solutions v 0 involve all equations 1D case still fairly trivial slab: when v x = 0 same as before 1D slab with arbitrary background flow v y (x) and v z (x) obeys d dx p + B2 y + Bz 2 2 } {{ } total pressure = ρg with gravity g = gê x density enters example solution: force-free (B y, B z ) and ρ = 1 Dx on x [0, 1] density contrast D and pressure p = p 0 g(x Dx2 2 )
29 CapSel Equil D cylindrical equilibria take v R = 0 then again only force balance remains assume gravitational line source g = GM d dr p + B2 ϕ + BZ 2 2 } {{ } total pressure = B2 ϕ R }{{} radially inward tension + R 2 ê R ρvϕ 2 } R{{} centrifugal force ρ GM suitable for accretion disks R [R inner, R outer ]: midplane Z ϕ R 2 }{{} gravity M * R for jets: neglect gravity and take R [0, R jet ]
30 CapSel Equil D (axisymmetric) stationary ideal MHD equilibria allow for both toroidal and poloidal flow in tokamak steady-state core problem is again 2nd order PDE for flux ψ(r, Z) BUT additional algebraic equation for M 2 (R, Z) ρvp/b 2 p 2 latter expresses energy conservation along fieldline much more complex problem: character of governing PDE can change from location to location elliptic or hyperbolic nature depends on M 2 turns out: 3 elliptic regimes: subslow, slow, fast elliptic flow regions ε ss ε s ε f 0 1 M 2 slow Alfven fast
31 CapSel Equil - 27 References J.P. Goedbloed, Transactions of Fusion Technology 37, 2T, march 2000 B.C. Low, Lecture series HAO 1994 J.P. Goedbloed & S. Poedts, Advanced Magnetohydrodynamics, Lecture series K.U. Leuven (CPA, Doctoraatsopleiding Wiskunde) J.P. Goedbloed & S. Poedts, Principles of Magnetohydrodynamics, Cambridge University Press, to appear 2003 J.P. Goedbloed, Magnetohydrodynamics of astrophysical plasmas, Lectures at Utrecht University 2001 J.P. Goedbloed & A.E. Lifschitz, Phys. Plasmas 4, 1997, 3544 J.P. Goedbloed, Physica Scripta T98, 2002, 43 S. Régnier web site solar.physics.montana.edu/regnier T. Amari web site
The 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 informationSimple examples of MHD equilibria
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
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 informationGeneralized Solovev equilibrium with sheared flow of arbitrary direction and stability consideration
Generalized Solovev equilibrium with sheared flow of arbitrary direction and stability consideration D.A. Kaltsas and G.N. Throumoulopoulos Department of Physics, University of Ioannina, GR 451 10 Ioannina,
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 informationUn problème inverse: l identification du profil de courant dans un Tokamak
Un problème inverse: l identification du profil de courant dans un Tokamak Blaise Faugeras en collaboration avec J. Blum et C. Boulbe Université de Nice Sophia Antipolis Laboratoire J.-A. Dieudonné Nice,
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 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 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 informationConfiguration Optimization of a Planar-Axis Stellarator with a Reduced Shafranov Shift )
Configuration Optimization of a Planar-Axis Stellarator with a Reduced Shafranov Shift ) Shoichi OKAMURA 1,2) 1) National Institute for Fusion Science, Toki 509-5292, Japan 2) Department of Fusion Science,
More informationMODELLING TWISTED FLUX TUBES PHILIP BRADSHAW (ASTROPHYSICS)
MODELLING TWISTED FLUX TUBES PHILIP BRADSHAW (ASTROPHYSICS) Abstract: Twisted flux tubes are important features in the Universe and are involved in the storage and release of magnetic energy. Therefore
More informationSOLAR MHD Lecture 2 Plan
SOLAR MHD Lecture Plan Magnetostatic Equilibrium ü Structure of Magnetic Flux Tubes ü Force-free fields Waves in a homogenous magnetized medium ü Linearized wave equation ü Alfvén wave ü Magnetoacoustic
More informationEMAFF: MAgnetic Equations with FreeFem++
EMAFF: MAgnetic Equations with FreeFem++ The Grad-Shafranov equation & The Current Hole Erwan DERIAZ Bruno DESPRÉS Gloria FACCANONI š Kirill Pichon GOSTAF Lise-Marie IMBERT-GÉRARD Georges SADAKA Remy
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 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 informationMagnetohydrodynamics (MHD) Philippa Browning Jodrell Bank Centre for Astrophysics University of Manchester
Magnetohydrodynamics (MHD) Philippa Browning Jodrell Bank Centre for Astrophysics University of Manchester MagnetoHydroDynamics (MHD) 1. The MHD equations. Magnetic Reynolds number and ideal MHD 3. Some
More informationMagnetohydrodynamics (MHD) Philippa Browning Jodrell Bank Centre for Astrophysics University of Manchester
Magnetohydrodynamics (MHD) Philippa Browning Jodrell Bank Centre for Astrophysics University of Manchester MagnetoHydroDynamics (MHD) 1. The MHD equations. Magnetic Reynolds number and ideal MHD 3. Some
More informationMagnetohydrodynamic Waves
Magnetohydrodynamic Waves Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 17, 2016 These slides are largely based off of 4.5 and 4.8 of The Physics of
More informationForce-Free Magnetosphere of an Accreting Kerr Black Hole
Force-Free Magnetosphere of an Accreting Kerr Black Hole Dmitri A. Uzdensky Princeton University, Princeton, NJ 08540, USA I consider a stationary axisymmetric force-free degenerate magnetosphere of a
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 informationPLASMA EQUILIBRIUM IN TOKAMAKS
PLAMA EQUILIBRIUM IN TOKAMAK H.J. de Blank FOM Institute DIFFER Dutch Institute for Fundamental Energy Research, Association EURATOM-FOM, P.O. Box 27, 343 BE Nieuwegein, The Netherlands, www.differ.nl.
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 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 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 informationExact solutions for magnetic annihilation in curvilinear geometry
Exact solutions for magnetic annihilation in curvilinear geometry E. Tassi b,, V.S. Titov and G. Hornig Theoretische Physik IV, Ruhr-Universität Bochum, 44780 Bochum, Germany b Theoretische Physik IV,
More informationarxiv:astro-ph/ v1 27 May 2005
2D stationary resistive MHD flows: borderline to magnetic reconnection solutions D.H. Nickeler a,, H.-J. Fahr b arxiv:astro-ph/0505554v1 27 May 2005 a Astronomical Institute, Utrecht University, Princetonplein
More informationStellar structure Conservation of mass Conservation of energy Equation of hydrostatic equilibrium Equation of energy transport Equation of state
Stellar structure For an isolated, static, spherically symmetric star, four basic laws / equations needed to describe structure: Conservation of mass Conservation of energy (at each radius, the change
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 informationVertical angular momentum transfer from accretion disks and the formation of large-scale collimated jets
Author manuscript, published in "Plasma Physics and Controlled Fusion 50 (2008) 4020" DOI : 10.1088/0741-3335/50/12/124020 35th EPS Conference on Plasma Physics Vertical angular momentum transfer from
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 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 informationLinear stability of MHD configurations
Linear stability of MHD configurations Rony Keppens Centre for mathematical Plasma Astrophysics KU Leuven Rony Keppens (KU Leuven) Linear MHD stability CHARM@ROB 2017 1 / 18 Ideal MHD configurations Interested
More informationProf. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit
Prof. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit Rough breakdown of MHD shocks Jump conditions: flux in = flux out mass flux: ρv n magnetic flux: B n Normal momentum flux: ρv n
More informationDT Fusion Ignition of LHD-Type Helical Reactor by Joule Heating Associated with Magnetic Axis Shift )
DT Fusion Ignition of LHD-Type Helical Reactor by Joule Heating Associated with Magnetic Axis Shift ) Tsuguhiro WATANABE National Institute for Fusion Science, 322-6 Oroshi-cho, Toki 509-5292, Japan (Received
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 informationMagnetohydrodynamic stability of negative central magnetic shear, high pressure ( pol 1) toroidal equilibria
Magnetohydrodynamic stability of negative central magnetic shear, high pressure ( pol 1) toroidal equilibria Robert G. Kleva Institute for Plasma Research, University of Maryland, College Park, Maryland
More informationS3: Description of Applied Focusing Fields S3A: Overview
S3: Description of Applied Focusing Fields S3A: Overview The fields of such classes of magnets obey the vacuum Maxwell Equations within the aperture: Applied fields for focusing, bending, and acceleration
More informationConservation Laws in Ideal MHD
Conservation Laws in Ideal MHD Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 3, 2016 These lecture notes are largely based on Plasma Physics for Astrophysics
More informationBefore you begin read these instructions carefully:
NATURAL SCIENCES TRIPOS Part IB & II (General Friday, 30 May, 2014 9:00 am to 12:00 pm MATHEMATICS (2 Before you begin read these instructions carefully: You may submit answers to no more than six questions.
More informationReconstruction of the Pressure Profile of LDX High Beta Plasma
Reconstruction of the Pressure Profile of LDX High Beta Plasma M. Mauel, I. Karim, A. Boxer, J. Ellsworth, D. Garnier, A. Hansen, J. Kesner, and E. Ortiz Dept. of Applied Physics and Applied Mathematics,
More informationIssues in Neoclassical Tearing Mode Theory
Issues in Neoclassical Tearing Mode Theory Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Austin, TX Tearing Mode Stability in Tokamaks According to standard (single-fluid)
More informationMHD Linear Stability Analysis Using a Full Wave Code
US-Japan JIFT Workshop on Progress of Extended MHD Models NIFS, Toki,Japan 2007/03/27 MHD Linear Stability Analysis Using a Full Wave Code T. Akutsu and A. Fukuyama Department of Nuclear Engineering, Kyoto
More informationProblem set: solar irradiance and solar wind
Problem set: solar irradiance and solar wind Karel Schrijver July 3, 203 Stratification of a static atmosphere within a force-free magnetic field Problem: Write down the general MHD force-balance equation
More informationImproved Analytical Flux Surface Representation and Calculation Models for Poloidal Asymmetries. T. G. Collart, W. M. Stacey
Improved Analytical Flux Surface Representation and Calculation Models for Poloidal Asymmetries T. G. Collart, W. M. Stacey Georgia Institute of Technology Atlanta, GA 3332 USA December, 215 Abstract An
More informationNIMEQ: MHD Equilibrium Solver for NIMROD
NIMEQ: MHD Equilibrium Solver for NIMOD E.C.Howell, C..Sovinec University of Wisconsin-Madison 5 th Annual Meeting of Division of Plasma Physics Dallas, Texas, Nov. 17-Nov. 1,8 1 Abstract A Grad-Shafranov
More informationPlasma formation in MAST by using the double null merging technique
1 Plasma formation in MAST by using the double null merging technique P. Micozzi 1, F. Alladio 1, P. Costa 1, A. Mancuso 1, A. Sykes 2, G. Cunningham 2, M. Gryaznevich 2, J. Hicks 2, M. Hood 2, G. McArdle
More informationNumerical calculation of the Hamada basis vectors for three-dimensional toroidal magnetic configurations
PHYSICS OF PLASMAS 12, 072513 2005 Numerical calculation of the Hamada basis vectors for three-dimensional toroidal magnetic configurations J. N. Talmadge and S. P. Gerhardt a HSX Plasma Laboratory, University
More informationMagnetohydrodynamics (MHD)
Magnetohydrodynamics (MHD) Robertus v F-S Robertus@sheffield.ac.uk SP RC, School of Mathematics & Statistics, The (UK) The Outline Introduction Magnetic Sun MHD equations Potential and force-free fields
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 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 informationF q. Gas at radius R (cylindrical) and height z above the disk midplane. F z. central mass M
Accretion Disks Luminosity of AGN derives from gravitational potential energy of gas spiraling inward through an accretion disk. Derive structure of the disk, and characteristic temperatures of the gas.
More informationApproximate Solutions of the Grad-Schlüter-Shafranov Equation
Approximate Solutions of the Grad-Schlüter-Shafranov Equation Gerson Otto Ludwig Associated Plasma Laboratory, National Space Research Institute 17-010, São José dos Campos, SP, Brazil ludwig@plasma.inpe.br
More informationOptimization approach for the computation of magnetohydrostatic coronal equilibria in spherical geometry ABSTRACT. 2.
A&A 475, 701 706 (2007) DOI: 10.1051/0004-6361:20078244 c ESO 2007 Astronomy & Astrophysics Optimization approach for the computation of magnetohydrostatic coronal equilibria in spherical geometry T. Wiegelmann
More informationSimulation Study of Interaction between Energetic Ions and Alfvén Eigenmodes in LHD
1 Simulation Study of Interaction between Energetic Ions and Alfvén Eigenmodes in LHD Y. Todo 1), N. Nakajima 1), M. Osakabe 1), S. Yamamoto 2), D. A. Spong 3) 1) National Institute for Fusion Science,
More informationWaves and characteristics: Overview 5-1
Waves and characteristics: Overview 5-1 Chapter 5: Waves and characteristics Overview Physics and accounting: use example of sound waves to illustrate method of linearization and counting of variables
More informationSolution Set Eight. 1 Problem #1: Toroidal Electromagnet with Gap Problem #4: Self-Inductance of a Long Solenoid. 9
: Solution Set Eight Northwestern University, Electrodynamics I Wednesday, March 9, 6 Contents Problem #: Toroidal Electromagnet with Gap. Problem #: Electromagnetic Momentum. 3 3 Problem #3: Type I Superconductor.
More informationMHD Simulation of High Wavenumber Ballooning-like Modes in LHD
1 TH/P9-16 MHD Simulation of High Wavenumber Ballooning-like Modes in LHD H. Miura and N. Nakajima National Institute for Fusion Science, 322-6 Oroshi, Toki, Gifu 509-5292, JAPAN e-mail contact of main
More informationThe correlation between magnetic flux and jet power
The correlation between magnetic flux and jet power Frontiers in Astronomy and Space Sciences 2017 Elena Nokhrina Moscow Institute of Physics and Technology IAU Symposium 342 - Perseus in Sicily: from
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 informationEdge Momentum Transport by Neutrals
1 TH/P3-18 Edge Momentum Transport by Neutrals J.T. Omotani 1, S.L. Newton 1,2, I. Pusztai 1 and T. Fülöp 1 1 Department of Physics, Chalmers University of Technology, 41296 Gothenburg, Sweden 2 CCFE,
More informationTwo Fluid Dynamo and Edge-Resonant m=0 Tearing Instability in Reversed Field Pinch
1 Two Fluid Dynamo and Edge-Resonant m= Tearing Instability in Reversed Field Pinch V.V. Mirnov 1), C.C.Hegna 1), S.C. Prager 1), C.R.Sovinec 1), and H.Tian 1) 1) The University of Wisconsin-Madison, Madison,
More informationStellar Winds: Mechanisms and Dynamics
Astrofysikalisk dynamik, VT 010 Stellar Winds: Mechanisms and Dynamics Lecture Notes Susanne Höfner Department of Physics and Astronomy Uppsala University 1 Most stars have a stellar wind, i.e. and outflow
More informationarxiv: v1 [astro-ph] 29 Jan 2008
Contrib. Astron. Obs. Skalnaté Pleso?, 1 6, (2018) Non-dipolar magnetic fields in Ap stars arxiv:0801.4562v1 [astro-ph] 29 Jan 2008 J.Braithwaite 1 Canadian Institute for Theoretical Astrophysics 60 St.
More informationLow Beta MHD Equilibrium Including a Static Magnetic Island for Reduced MHD Equations in a Straight Heliotron Configuration
Low Beta MHD Equilibrium Including a Static Magnetic Island for Reduced MHD Equations in a Straight Heliotron Configuration Kinya SAITO 1,a), Katsuji ICHIGUCHI 1,2) and Ryuichi ISHIZAKI 1,2) 1) The Graduate
More informationGravitational steady states of solar coronal loops. Abstract
Gravitational steady states of solar coronal loops Linda E. Sugiyama Laboratory for Nuclear Science, Massachusetts Institute of Technology, Cambridge MA 02139 M. Asgari-Targhi Harvard-Smithsonian Center
More informationPAPER 68 ACCRETION DISCS
MATHEMATICAL TRIPOS Part III Tuesday, 2 June, 2009 9:00 am to 11:00 am PAPER 68 ACCRETION DISCS There are THREE questions in total Full marks can be obtained by completing TWO questions The questions carry
More informationL Aquila, Maggio 2002
Nonlinear saturation of Shear Alfvén Modes and energetic ion transports in Tokamak equilibria with hollow-q profiles G. Vlad, S. Briguglio, F. Zonca, G. Fogaccia Associazione Euratom-ENEA sulla Fusione,
More informationD.J. Schlossberg, D.J. Battaglia, M.W. Bongard, R.J. Fonck, A.J. Redd. University of Wisconsin - Madison 1500 Engineering Drive Madison, WI 53706
D.J. Schlossberg, D.J. Battaglia, M.W. Bongard, R.J. Fonck, A.J. Redd University of Wisconsin - Madison 1500 Engineering Drive Madison, WI 53706 Concept Overview Implementation on PEGASUS Results Current
More informationInfluence of Beta, Shape and Rotation on the H-mode Pedestal Height
Influence of Beta, Shape and Rotation on the H-mode Pedestal Height by A.W. Leonard with R.J. Groebner, T.H. Osborne, and P.B. Snyder Presented at Forty-Ninth APS Meeting of the Division of Plasma Physics
More informationExperimental Study of Hall Effect on a Formation Process of an FRC by Counter-Helicity Spheromak Merging in TS-4 )
Experimental Study of Hall Effect on a Formation Process of an FRC by Counter-Helicity Spheromak Merging in TS-4 ) Yasuhiro KAMINOU, Michiaki INOMOTO and Yasushi ONO Graduate School of Engineering, The
More informationIs the Troyon limit a beta limit?
Is the Troyon limit a beta limit? Pierre-Alexandre Gourdain 1 1 Extreme State Physics Laboratory, Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627, USA The plasma beta,
More informationMAGNETOHYDRODYNAMIC 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 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 informationTheoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions
Theoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions Tom Elsden 1 Andrew Wright 1 1 Dept Maths & Stats, University of St Andrews DAMTP Seminar - 8th May 2017 Outline Introduction Coordinates
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 informationGeneral Solution of the Incompressible, Potential Flow Equations
CHAPTER 3 General Solution of the Incompressible, Potential Flow Equations Developing the basic methodology for obtaining the elementary solutions to potential flow problem. Linear nature of the potential
More informationD.J. Schlossberg, D.J. Battaglia, M.W. Bongard, R.J. Fonck, A.J. Redd. University of Wisconsin - Madison 1500 Engineering Drive Madison, WI 53706
D.J. Schlossberg, D.J. Battaglia, M.W. Bongard, R.J. Fonck, A.J. Redd University of Wisconsin - Madison 1500 Engineering Drive Madison, WI 53706 Non-solenoidal startup using point-source DC helicity injectors
More informationª 10 KeV. In 2XIIB and the tandem mirrors built to date, in which the plug radius R p. ª r Li
Axisymmetric Tandem Mirrors: Stabilization and Confinement Studies R. F. Post, T. K. Fowler*, R. Bulmer, J. Byers, D. Hua, L. Tung Lawrence Livermore National Laboratory *Consultant, Presenter This talk
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 informationMHD equilibrium calculations for stellarators. Antoine Cerfon, MIT PSFC with F. Parra, J. Freidberg (MIT NSE)
MHD equilibrium calculations for stellarators Antoine Cerfon, MIT PSFC with F. Parra, J. Freidberg (MIT NSE) March 20, 2012 MAGNETIC FIELD LINE HAMILTONIAN Consider a general toroidal coordinate system
More informationarxiv: v2 [physics.plasm-ph] 8 Mar 2018
Quasi-axisymmetric magnetic fields: weakly non-axisymmetric case in a vacuum G.G. Plunk 1 and P. Helander 1 1 Max-Planck-Institut für Plasmaphysik, EURATOM Association, 17491 Greifswald, Germany arxiv:1801.02990v2
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 informationChapter 1 Hamiltonian Representation of Magnetic Field
Chapter 1 Hamiltonian Representation of Magnetic Field In this chapter we present the Hamiltonian formulation of the equations for magnetic field lines. We specifically consider the magnetic field corresponding
More informationHeat Flux Management via Advanced Magnetic Divertor Configurations and Divertor Detachment.
Heat Flux Management via Advanced Magnetic Divertor Configurations and Divertor Detachment E. Kolemen a, S.L. Allen b, B.D. Bray c, M.E. Fenstermacher b, D.A. Humphreys c, A.W. Hyatt c, C.J. Lasnier b,
More informationMHD modeling of the kink double-gradient branch of the ballooning instability in the magnetotail
MHD modeling of the kink double-gradient branch of the ballooning instability in the magnetotail Korovinskiy 1 D., Divin A., Ivanova 3 V., Erkaev 4,5 N., Semenov 6 V., Ivanov 7 I., Biernat 1,8 H., Lapenta
More information02. Multipole Fields *
02. Multipole Fields * Prof. Steven M. Lund Physics and Astronomy Department Facility for Rare Isotope Beams (FRIB) Michigan State University (MSU) US Particle Accelerator School Accelerator Physics Steven
More informationGA A22863 PLASMA PRESSURE AND FLOWS DURING DIVERTOR DETACHMENT
GA A22863 PLASMA PRESSURE AND FLOWS DURING DIVERTOR DETACHMENT by M.J. SCHAFFER, J.A. BOEDO, N.H. BROOKS, R.C. ISLER, and R.A. MOYER AUGUST 1998 DISCLAIMER This report was prepared as an account of work
More informationToroidal plasmas: Overview T-1
Toroidal plasmas: Overview T-1 Toroidal plasmas Overview Introduction: Context is Alfvén wave dynamics, statement of equilibrium problem, cylindrical equilibrium, rotational transform, importance for stability;
More informationMHD Simulation of Solar Chromospheric Evaporation Jets in the Oblique Coronal Magnetic Field
MHD Simulation of Solar Chromospheric Evaporation Jets in the Oblique Coronal Magnetic Field Y. Matsui, T. Yokoyama, H. Hotta and T. Saito Department of Earth and Planetary Science, University of Tokyo,
More informationBogoyavlenskij Symmetries of Isotropic and Anisotropic MHD Equilibria as Lie Point Transformations
Proceedings of Institute of Mathematics of NAS of Ukraine 2004, Vol. 50, Part 1, 69 76 Bogoyavlenskij Symmetries of Isotropic and Anisotropic MHD Equilibria as Lie Point Transformations Alexei F. CHEVIAKOV
More informationRole of the mean curvature in the geometry of magnetic confinement configurations
arxiv:1008.1567v2 [physics.plasm-ph] 11 Aug 2010 Role of the mean curvature in the geometry of magnetic confinement configurations A.A. Skovoroda I.A. Taimanov Abstract Examples are presented of how the
More informationADVANCES IN NOVEL FRC PLASMAS
GA A23824 ADVANCES IN NOVEL FRC PLASMAS by M.J. SCHAFFER, M.R. BROWN, J.A. LEUER, P.B. PARKS, and C.D. COTHRAN APRIL 22 DISCLAIMER This report was prepared as an account of work sponsored by an agency
More informationCoronal Magnetic Field Extrapolations
3 rd SOLAIRE School Solar Observational Data Analysis (SODAS) Coronal Magnetic Field Extrapolations Stéphane RÉGNIER University of St Andrews What I will focus on Magnetic field extrapolation of active
More informationMulti-D MHD and B = 0
CapSel DivB - 01 Multi-D MHD and B = 0 keppens@rijnh.nl multi-d MHD and MHD wave anisotropies dimensionality > 1 non-trivial B = 0 constraint even if satisfied exactly t = 0: can numerically generate B
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 information6.1 Momentum Equation for Frictionless Flow: Euler s Equation The equations of motion for frictionless flow, called Euler s
Chapter 6 INCOMPRESSIBLE INVISCID FLOW All real fluids possess viscosity. However in many flow cases it is reasonable to neglect the effects of viscosity. It is useful to investigate the dynamics of an
More informationSheared Flow Stabilization in the Z-Pinch
1 IF/P7-28 Sheared Flow Stabilization in the Z-Pinch U. Shumlak, C.S. Adams, J.M. Blakely, B.J. Chan, R.P. Golingo, S.D. Knecht, B.A. Nelson, R.J. Oberto, M.R. Sybouts, and G.V. Vogman Aerospace & Energetics
More informationNoninductive Formation of Spherical Tokamak at 7 Times the Plasma Cutoff Density by Electron Bernstein Wave Heating and Current Drive on LATE
1 EX/P6-18 Noninductive Formation of Spherical Tokamak at 7 Times the Plasma Cutoff Density by Electron Bernstein Wave Heating and Current Drive on LATE M. Uchida, T. Maekawa, H. Tanaka, F. Watanabe, Y.
More informationDIAGNOSTICS FOR ADVANCED TOKAMAK RESEARCH
DIAGNOSTICS FOR ADVANCED TOKAMAK RESEARCH by K.H. Burrell Presented at High Temperature Plasma Diagnostics 2 Conference Tucson, Arizona June 19 22, 2 134 /KHB/wj ROLE OF DIAGNOSTICS IN ADVANCED TOKAMAK
More information