The Virial Theorem, MHD Equilibria, and Force-Free Fields
|
|
- Sydney Williams
- 6 years ago
- Views:
Transcription
1 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 based on Plasma Physics for Astrophysics by Russell Kulsrud, Lectures in Magnetohydrodynamics by the late Dalton Schnack, Ideal Magnetohydrodynamics by Jeffrey Freidberg, Hydrodynamic and Hydromagnetic Stability by S. Chandrasekhar, Classical Electrodynamics by J. Jackson, and examples by A. Savcheva and A. Spitkovsky.
2 Outline We will look at the properties and key characteristics of MHD equilibria. Our discussion will focus on: The virial theorem 0 = 2E V + 3 (γ 1) E p + E B + E g MHD equilibria J B Force-free fields c = p J B = 0 B = αb
3 The Virial Theorem for MHD (following Kulsrud 4.6) The Virial Theorem allows us to understand broadly the equilibrium properties of a system in terms of energies Suppose there exists a magnetized plasma within a finite volume. The scalar moment of inertia is I = 1 ρr 2 dv (1) 2 where r is the position vector about some arbitrary origin Our strategy: Calculate di/dt and d 2 I/dt 2 Ignore surface integrals by assuming the volume is large Put the result in terms of energies Set d 2 I/dt 2 = 0 for an equilibrium Determine the conditions under which the resulting equation can be satisfied V
4 Take the first time derivative of I Use the continuity equation, the radial form r 2 = 2r, the identity (f A) f A + A f, and Gauss theorem. di = 1 ρ dt 2 V t r 2 dv = 1 (ρv) r 2 dv 2 V = 1 (ρvr 2) dv + 1 ρv r 2 dv 2 V 2 = 1 r 2 ρv ds + ρv r dv (2) 2 S The moment of inertia changes by mass entering the system, or mass moving toward or away from the origin Consider a volume large enough so that no mass enters or leaves. The surface integral then vanishes: di dt = ρv r dv (3) V V
5 Take the second time derivative of I Use the momentum equation, r = I, the tensor identity A T = (A T) + T : A, and Gauss theorem d 2 I dt 2 = r (ρv) dv = ( T) r dv V t V = (T r) dv + T : r dv V V = ds T r + trace (T) dv (4) S The first term represents surface stresses. If we assume that surface stresses are negligible, then we are left with V d 2 I dt 2 = trace (T) dv (5) V
6 Let s look again at the stress tensor T (with gravity) We need to take T = ρvv + pi + B2 8π I BB 4π + ( φ)2 I 8πG ( φ)( φ) 4πG (6) The Reynolds stress is ρvv = V x V x V x V y V x V z V x V y V y V y V y V z V x V z V y V z V z V z (7) Then take its trace by adding up the diagonal elements: trace (ρvv) = ρ ( Vx 2 + Vy 2 + Vz 2 ) = ρv 2 (8) It s just twice the kinetic energy density!
7 Now evaluate the traces of the other terms in T The traces of the other stress tensors yield energy densities times constants! trace ( B 2 8π I BB trace ( ( φ) 2 I 8πG ( φ)( φ) 4πG trace (pi) = 3p (9) ) = B2 (10) 4π 8π ) = ( φ)2 (11) 8πG Recall that the internal energy density is given by p/(γ 1)
8 Now let s put these back into the volume integral By replacing trace (T) in Eq. 5 we arrive at d 2 ( ) I dt 2 = ρv 2 + 3p + B2 8π ( φ)2 dv 8πG where V = 2E V + 3 (γ 1) E p + E B + E g (12) Kinetic energy: EV 0 Internal energy: Ep 0 Magnetic energy: EB 0 Gravitational energy: E g 0 (only possible negative term!) In an equilibrium, this expression must equal zero: 0 = 2E V + 3 (γ 1) E p + E B + E g (13) We often consider averages of confined systems over a long time.
9 What happens when we neglect magnetic and internal energy? If E B = E g = 0, then we recover E V = 1 2 E g (14) The kinetic energy must equal half the magnitude of the gravitational energy. This is a well-known result in self-gravitating systems such as star clusters, galaxies, and galaxy clusters This result has been used to infer the presence of dark matter
10 What happens when we drop gravity in a static system? In the absence of gravity and bulk motions, we are left with 0 = 3 (γ 1) E p + E B (15) But E p 0 and E B 0! We have a contradiction! A magnetized plasma cannot be in MHD equilibrium under forces generated solely by its own internal currents. Equilibria are possible if there are external currents as in laboratory plasmas Accounted for from the surface integrals we dropped In astrophysics, this might not be satisfied
11 What limits on the magnetic energy does the Virial Theorem imply? If E V = E p = 0, then 2E V + 3 (γ 1) E p + E B + E g = 0 (16) For a stable equilibrium, the magnetic energy must not exceed the magnitude of the gravitational energy. The virial theorem provides broad insight into the equilibrium properties of a system without having to worry about the details
12 MHD Equilibria We often care about systems that are in equilibrium Let s look at the momentum equation (neglecting gravity) ( ) ρ t + V V = J B p (17) c For a static equilibrium, the configuration must have J B c in the absence of other forces = p (18)
13 Properties of MHD equilibria Take B of the equilibrium equation: B ( p) = ( ) J B B c (19) B p = 0 (20) where we use that B is orthogonal to J B. B p is the directional derivative of p in the direction of B Plasma pressure is constant along magnetic field lines Similiarly, if we take J the equilibrium equation then J p = 0 (21) since J is orthogonal to J B also.
14 Effects of fast thermal conduction on equilibria Ideal MHD does not include thermal conduction However, thermal conduction is very fast along field lines! If temperature is constant along field lines, then For p = nkt, then we also have B T = 0 (22) B n = 0 (23) Eqs. 22 and 23 are not exact results, but rather arise from the approximation of infinitely fast parallel thermal conduction.
15 Example: equilibria with a unidirectional magnetic field Consider a configuration where the magnetic field is purely in the ẑ direction The equilibrium condition is then p + B z 8π = p tot (24) where the total pressure, p tot, is a constant The tension forces disappear because the field lines are straight
16 Example: consider 1D cylindrical equilibria A Z-pinch (above) has current flowing in the ẑ direction so that B is purely azimuthal A θ-pinch has current flowing in the ˆθ direction so B is purely axial A screw pinch has components of J and B in both the axial and azimuthal directions For these configurations, we look for solutions of the form p = p(r) ; J = J θ (r) ˆθ + J z (r)ẑ ; B = B θ (r) ˆθ + B z (r)ẑ (25) for which B = 0 is trivially satisfied
17 Finding a Z-pinch 1D equilibrium Set J θ = 0 and B z = 0 since current is purely axial Ampere s law becomes J z (r) = c 1 d 4π r dr (rb θ) (26) The ˆr component of the momentum equation is We then apply Eq. 26 J z B θ = dp dr (27) dp dr + c B θ d 4π r dr (rb θ) = 0 (28)
18 Finding a Z-pinch 1D equilibrium This can be rearranged to ( ) d p + B2 θ + B2 θ = 0 (29) dr 8π 4πr }{{}}{{} total pressure tension Or, putting this in terms of the curvature vector, ( ) p + B2 B2 8π 4π κ = 0 (30)
19 A Z-pinch equilibrium can be found by specifying B θ (r) and then solving for p(r) Shown above is the Bennett pinch with B θ r r 2 + r 2 0 ; p, J z r 2 0 (r 2 + r 2 0 )2 (31) If the domain is r [0, ] then the magnetic energy diverges! Need an outer wall, which is not present in astrophysics.
20 Axisymmetric equilibria are found by solving the Grad-Shafranov equation Astrophysical applications include flux ropes in the corona/solar wind, compact object magnetospheres, etc. Fundamentally important for fusion devices like tokamaks
21 To derive the Grad-Shafranov equation, we first define a flux function We introduce a flux function ψ such that B r = 1 ψ r z 1 ψ B z = r t which satisfies the divergence constraint for any B θ (r, z). Contours of constant ψ represent the projection of magnetic field lines into the poloidal (r-z) plane (32) (33)
22 We apply Ampere s law with this magnetic field configuration The current density for this configuration is where the operator is J r = c 1 4π r z (rb θ) (34) c 1 J z = 4π r r (rb θ) (35) c ψ J θ = (36) 4π r r r ( 1 r r ) + 2 z 2 (37)
23 Now look again at the equilibrium relation Equilibrium requires J B c By symmetry, p/ θ = 0 so that This then becomes 1 r 2 = p (38) J z B r J r B z = 0 (39) [ r (rb θ) ψ z z (rb θ) ψ r ] = 0 (40) To satisfy this, rb θ must be a function of ψ alone, so we define rb θ F (ψ) (41)
24 Deriving the Grad-Shafranov equation Likewise, p is constant along field lines and is also a function of ψ: p = p (ψ) (42) The r component of the equilibrium relation is j θ B z j z B θ = p r Combining these results we arrive at ψ 4πr 1 ψ r r 1 df ψ F πr dψ r r = dp ψ dψ r (43) (44) which simplifies to the Grad-Shafranov equation ψ + F df dψ = 4πr 2 dp dψ (45)
25 How do we find a solution of the Grad-Shafranov equation? The Grad-Shafranov equation is given by where ψ + F df dψ = 4πr 2 dp dψ, (46) r ( ) r r r z 2 (47) F (ψ) rb θ (48) To solve the Grad-Shafranov equation, we need to Specify p(ψ) Specify F (ψ) Solve for ψ The Grad-Shafranov equation is usually solved numerically
26 Equilibria and the Virial Theorem The natural state of a flux rope is to try to expand to infinity In laboratory plasmas, conducting wall outer boundaries and externally applied magnetic fields prevent this These show up as surface integrals in the Virial Theorem In astrophysical plasmas, a flux rope can be held in place from J B and p forces from the surrounding medium The surrounding medium, in turn, can be held in place by gravitational forces Example: the solar corona
27 How does gravity change things? The equilibrium condition becomes 0 = J B c p + ρg (49) B and J are no longer necessarily orthogonal to p B p = ρb g (50) J p = ρj g (51) A radially outward flux tube reduces to the case of hydrostatic equilibrium: p = ρg The Grad-Shafranov equation can be generalized to include gravitational forces I decided against assigning this as a homework problem.
28 Force-free fields When pressure is constant or in the limit of β 0, the pressure gradient force vanishes. Equilibria then have J B = 0 (52) Such configurations are called force-free because there is no Lorentz force and no plasma pressure gradient force Using Ampere s law, this reduces to where α is constant along field lines B = αb (53) Vector fields parallel to their own curl are called Beltrami fields
29 Linear force-free fields have constant α Start with the condition for a linear force-free field: Take its curl: B = αb (54) ( B) = α B (55) Use Eq. 54 for the RHS and vector identities ( B) 2 B = α 2 B (56) Linear force-free fields obey the Helmholtz equation 2 B + α 2 B = 0 (57) which can be solved using separation of variables, Green s functions, Fourier series, or numerically.
30 Nonlinear force-free fields (NLFFF) have non-constant α Finding nonlinear force-free fields requires solving B = α (x) B (58) where α must still remain constant along field lines Some analytical examples exist in 2D 3D solutions must generally be found numerically
31 Example: pulsar magnetospheres I Force-free field models are able to capture features that do not exist in a vacuum solution I Figure of oblique rotator by Anatoly Spitkovsky. Color marks direction of toroidal field.
32 Example: NLFFF modeling of sigmoidal solar active regions Hinode/XRT NLFFF model/photospheric magnetic flux contours Significant free energy is stored in the sheared/twisted magnetic field of sigmoids (S-shaped structures in the corona) The modeling strategy is to introduce a flux rope into a potential field model, let it relax, and see if it matches the observed X-ray emission (Savcheva et al. 2012)
33 Potential magnetic fields represent the minimum energy condition for a configuration s boundary conditions These are a special case of force-free fields with α = 0 Potential fields have J = 0, so Ampere s law becomes 0 = B (59) The curl of a gradient is zero, so this is satisfied if where ζ is a scalar magnetic potential B = ζ (60) Taking the divergence of both sides shows that ζ can be found by solving Laplace s equation B = ζ (61) 0 = 2 ζ (62)
34 Potential magnetic fields can be calcuated using Laplace s equation for the vector potential From Ampere s law with no currents we know that From these two expressions we arrive at B = A (63) 0 = c 4π B (64) ( A) = 0 (65) Using vector identities and choosing the Coulomb gauge ( A = 0) we arrive at 2 A = 0 (66)
35 Determining the magnetic free energy of a system The magnetic free energy of a system is given by ( ) B 2 E M 8π B2 p dv (67) 8π V where B p is the potential field solution for the system s boundary conditions However, the quantity B 2 8π B2 p 8π (68) should not be considered a free energy density since free energy is an integral quantity. At times this could be negative! Current density is a good proxy for how much stress exists in a magnetic field
36 Summary The Virial Theorem allows us to understand key requirements for equilibria even before we know the details A plasma cannot be in an MHD equilibrium generated solely by its internal currents 2D MHD equilibria are found using the Grad-Shafranov equation However, many equilibria are unstable 3D MHD equilibria are difficult to find and/or generalize Force-free fields have J B = 0 so that B = αb Linear force-free fields have constant α and can be solved for using the Helmholtz equation Nonlinear force-free fields have nonconstant α and usually are solved numerically Potential fields are the minimum energy state for given boundary conditions
Conservation 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 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 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 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 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 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 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 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 informationBeyond Ideal MHD. Nick Murphy. Harvard-Smithsonian Center for Astrophysics. Astronomy 253: Plasma Astrophysics. February 8, 2016
Beyond Ideal MHD Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 8, 2016 These lecture notes are largely based on Plasma Physics for Astrophysics by
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 informationProblems in Magnetostatics
Problems in Magnetostatics 8th February 27 Some of the later problems are quite challenging. This is characteristic of problems in magnetism. There are trivial problems and there are tough problems. Very
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 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 informationCHAPTER 8 CONSERVATION LAWS
CHAPTER 8 CONSERVATION LAWS Outlines 1. Charge and Energy 2. The Poynting s Theorem 3. Momentum 4. Angular Momentum 2 Conservation of charge and energy The net amount of charges in a volume V is given
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 informationShell Balances in Fluid Mechanics
Shell Balances in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University When fluid flow occurs in a single direction everywhere in a system, shell
More information1 Fundamentals. 1.1 Overview. 1.2 Units: Physics 704 Spring 2018
Physics 704 Spring 2018 1 Fundamentals 1.1 Overview The objective of this course is: to determine and fields in various physical systems and the forces and/or torques resulting from them. The domain of
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 informationElectromagnetic energy and momentum
Electromagnetic energy and momentum Conservation of energy: the Poynting vector In previous chapters of Jackson we have seen that the energy density of the electric eq. 4.89 in Jackson and magnetic eq.
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 informationMomentum transport from magnetic reconnection in laboratory an. plasmas. Fatima Ebrahimi
Momentum transport from magnetic reconnection in laboratory and astrophysical plasmas Space Science Center - University of New Hampshire collaborators : V. Mirnov, S. Prager, D. Schnack, C. Sovinec Center
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 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 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 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 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 informationFundamental Stellar Parameters. Radiative Transfer. Stellar Atmospheres
Fundamental Stellar Parameters Radiative Transfer Stellar Atmospheres Equations of Stellar Structure Basic Principles Equations of Hydrostatic Equilibrium and Mass Conservation Central Pressure, Virial
More informationLecture 8. Engineering applications
Lecture 8 Engineering applications In Lecture 6 we saw one classic example of the application of vector calculus to Maxwell s equation. In this lecture we explore a few more examples from fluid mechanics
More informationAY202a Galaxies & Dynamics Lecture 7: Jeans Law, Virial Theorem Structure of E Galaxies
AY202a Galaxies & Dynamics Lecture 7: Jeans Law, Virial Theorem Structure of E Galaxies Jean s Law Star/Galaxy Formation is most simply defined as the process of going from hydrostatic equilibrium to gravitational
More informationThe Magnetorotational Instability
The Magnetorotational Instability Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics March 10, 2014 These slides are based off of Balbus & Hawley (1991), Hawley
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 informationThe Physics of Fluids and Plasmas
The Physics of Fluids and Plasmas An Introduction for Astrophysicists ARNAB RAI CHOUDHURI CAMBRIDGE UNIVERSITY PRESS Preface Acknowledgements xiii xvii Introduction 1 1. 3 1.1 Fluids and plasmas in the
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 informationTutorial 3 - Solutions Electromagnetic Waves
Tutorial 3 - Solutions Electromagnetic Waves You can find formulas you require for vector calculus at the end of this tutorial. 1. Find the Divergence and Curl of the following functions - (a) k r ˆr f
More informationComputations of Vector Potential and Toroidal Flux and Applications to Stellarator Simulations
Computations of Vector Potential and Toroidal Flux and Applications to Stellarator Simulations NIMROD Team Meeting Torrin Bechtel April 30, 2017 Outline 1 Project Goals and Progress 2 Vector Potential
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 informationA SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations
A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations Hervé Lemonnier DM2S/STMF/LIEFT, CEA/Grenoble, 38054 Grenoble Cedex 9 Ph. +33(0)4 38 78 45 40 herve.lemonnier@cea.fr, herve.lemonnier.sci.free.fr/tpf/tpf.htm
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 informationPotential/density pairs and Gauss s law
Potential/density pairs and Gauss s law We showed last time that the motion of a particle in a cluster will evolve gradually, on the relaxation time scale. This time, however, is much longer than the typical
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 informationPlasma spectroscopy when there is magnetic reconnection associated with Rayleigh-Taylor instability in the Caltech spheromak jet experiment
Plasma spectroscopy when there is magnetic reconnection associated with Rayleigh-Taylor instability in the Caltech spheromak jet experiment KB Chai Korea Atomic Energy Research Institute/Caltech Paul M.
More informationRadiation Integrals and Auxiliary Potential Functions
Radiation Integrals and Auxiliary Potential Functions Ranga Rodrigo June 23, 2010 Lecture notes are fully based on Balanis [?]. Some diagrams and text are directly from the books. Contents 1 The Vector
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 informationOn existence of resistive magnetohydrodynamic equilibria
arxiv:physics/0503077v1 [physics.plasm-ph] 9 Mar 2005 On existence of resistive magnetohydrodynamic equilibria H. Tasso, G. N. Throumoulopoulos Max-Planck-Institut für Plasmaphysik Euratom Association
More information2. Equations of Stellar Structure
2. Equations of Stellar Structure We already discussed that the structure of stars is basically governed by three simple laws, namely hyostatic equilibrium, energy transport and energy generation. In this
More informationChapter 1. Maxwell s Equations
Chapter 1 Maxwell s Equations 11 Review of familiar basics Throughout this course on electrodynamics, we shall use cgs units, in which the electric field E and the magnetic field B, the two aspects of
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 informationIndiana University Physics P331: Theory of Electromagnetism Review Problems #3
Indiana University Physics P331: Theory of Electromagnetism Review Problems #3 Note: The final exam (Friday 1/14 8:00-10:00 AM will be comprehensive, covering lecture and homework material pertaining to
More informationSpace Plasma Physics Thomas Wiegelmann, 2012
Space Plasma Physics Thomas Wiegelmann, 2012 1. Basic Plasma Physics concepts 2. Overview about solar system plasmas Plasma Models 3. Single particle motion, Test particle model 4. Statistic description
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 informationThe Distribution Function
The Distribution Function As we have seen before the distribution function (or phase-space density) f( x, v, t) d 3 x d 3 v gives a full description of the state of any collisionless system. Here f( x,
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 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 informationStellar Winds. Star. v w
Stellar Winds Star v w Stellar Winds Geoffrey V. Bicknell 1 Characteristics of stellar winds Solar wind Velocity at earth s orbit: Density: Temperature: Speed of sound: v 400 km/s n 10 7 m 3 c s T 10 5
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 informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationPhysics 505 Fall 2005 Homework Assignment #7 Solutions
Physics 505 Fall 005 Homework Assignment #7 Solutions Textbook problems: Ch. 4: 4.10 Ch. 5: 5.3, 5.6, 5.7 4.10 Two concentric conducting spheres of inner and outer radii a and b, respectively, carry charges
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
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 informationSummary of the Equations of Fluid Dynamics
Reference: Summary of the Equations of Fluid Dynamics Fluid Mechanics, L.D. Landau & E.M. Lifshitz 1 Introduction Emission processes give us diagnostics with which to estimate important parameters, such
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 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 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 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 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 informationPhysical Processes in Astrophysics
Physical Processes in Astrophysics Huirong Yan Uni Potsdam & Desy Email: hyan@mail.desy.de 1 Reference Books: Plasma Physics for Astrophysics, Russell M. Kulsrud (2005) The Physics of Astrophysics, Frank
More informationElectrodynamics of Magnetized Rotators Anatoly Spitkovsky,, UC Berkeley
Electrodynamics of Magnetized Rotators Anatoly Spitkovsky,, UC Berkeley Magnetized rotators are ubiquitous: pulsars, AGN, GRBs (?) Rotation very efficient at long-term energy storage Extraction of rotational
More informationPEAT SEISMOLOGY Lecture 3: The elastic wave equation
PEAT8002 - SEISMOLOGY Lecture 3: The elastic wave equation Nick Rawlinson Research School of Earth Sciences Australian National University Equation of motion The equation of motion can be derived by considering
More information1 Shapes of elliptical galaxies
Lecture 3 : Elliptical galaxies II 1 Shapes of elliptical galaxies Figure 1: Observed distribution of apparent axial ratios (Ryden, 1992, ApJ, 396, 445) The observed ellipticity of a galaxy is a combination
More informationDIVERGENCE AND CURL THEOREMS
This document is stored in Documents/4C/Gausstokes.tex. with LaTex. Compile it November 29, 2014 Hans P. Paar DIVERGENCE AND CURL THEOREM 1 Introduction We discuss the theorems of Gauss and tokes also
More informationLecturer: Bengt E W Nilsson
9 3 19 Lecturer: Bengt E W Nilsson Last time: Relativistic physics in any dimension. Light-cone coordinates, light-cone stuff. Extra dimensions compact extra dimensions (here we talked about fundamental
More informationSummary: Applications of Gauss Law
Physics 2460 Electricity and Magnetism I, Fall 2006, Lecture 15 1 Summary: Applications of Gauss Law 1. Field outside of a uniformly charged sphere of radius a: 2. An infinite, uniformly charged plane
More informationThe Madison Dynamo Experiment: magnetic instabilities driven by sheared flow in a sphere. Cary Forest Department of Physics University of Wisconsin
The Madison Dynamo Experiment: magnetic instabilities driven by sheared flow in a sphere Cary Forest Department of Physics University of Wisconsin February 28, 2001 Planets, stars and perhaps the galaxy
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 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 informationMeridional Flow, Differential Rotation, and the Solar Dynamo
Meridional Flow, Differential Rotation, and the Solar Dynamo Manfred Küker 1 1 Leibniz Institut für Astrophysik Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany Abstract. Mean field models of rotating
More informationκ = f (r 0 ) k µ µ k ν = κk ν (5)
1. Horizon regularity and surface gravity Consider a static, spherically symmetric metric of the form where f(r) vanishes at r = r 0 linearly, and g(r 0 ) 0. Show that near r = r 0 the metric is approximately
More informationAstrophysical Dynamos
Astrophysical Dynamos Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics April 19, 2016 These lecture notes are based off of Kulsrud, Cowling (1981), Beck et al.
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 informationChapter 4. MHD Equilibrium and Stability. 4.1 Basic Two-Dimensional Equilibrium Equations and Properties. Resistive Diffusion
Chapter 4 MHD Equilibrium and Stability Resistive Diffusion Before discussing equilibrium properties let us first consider effects of electric resistivity. Using the resistive form of Ohm s law with constant
More informationComment about Didactical formulation of the
Comment about Didactical formulation of the Ampère law Hendrik van Hees Institute for Theoretical Physics, Goethe University Frankfurt, Max-von-Laue-Str. 1, D-60438 Frankfurt, Germany Frankfurt Institute
More information2nd Year Electromagnetism 2012:.Exam Practice
2nd Year Electromagnetism 2012:.Exam Practice These are sample questions of the type of question that will be set in the exam. They haven t been checked the way exam questions are checked so there may
More informationAim: Understand equilibrium of galaxies
8. Galactic Dynamics Aim: Understand equilibrium of galaxies 1. What are the dominant forces? 2. Can we define some kind of equilibrium? 3. What are the relevant timescales? 4. Do galaxies evolve along
More informationA cylinder in a magnetic field (Jackson)
Problem 1. A cylinder in a magnetic field (Jackson) A very long hollow cylinder of inner radius a and outer radius b of permeability µ is placed in an initially uniform magnetic field B o at right angles
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 informationwhere G is Newton s gravitational constant, M is the mass internal to radius r, and Ω 0 is the
Homework Exercise Solar Convection and the Solar Dynamo Mark Miesch (HAO/NCAR) NASA Heliophysics Summer School Boulder, Colorado, July 27 - August 3, 2011 PROBLEM 1: THERMAL WIND BALANCE We begin with
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 informationPart III. Flow and dissipation
Part III Flow and dissipation in this web service 1 Waves and instabilities of stationary plasmas 1.1 Laboratory and astrophysical plasmas 1.1.1 Grand vision: magnetized plasma on all scales In Chapter
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 informationStellar Magnetospheres part deux: Magnetic Hot Stars. Stan Owocki
Stellar Magnetospheres part deux: Magnetic Hot Stars Stan Owocki Key concepts from lec. 1 MagRe# --> inf => ideal => frozen flux breaks down at small scales: reconnection Lorentz force ~ mag. pressure
More information1 Introduction to Governing Equations 2 1a Methodology... 2
Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................
More information4Divergenceandcurl. D ds = ρdv. S
4Divergenceandcurl Epressing the total charge Q V contained in a volume V as a 3D volume integral of charge density ρ(r), wecanwritegauss s law eamined during the last few lectures in the general form
More informationElectromagnetic Theory Prof. D. K. Ghosh Department of Physics Indian Institute of Technology, Bombay
Electromagnetic Theory Prof. D. K. Ghosh Department of Physics Indian Institute of Technology, Bombay Module - 4 Time Varying Field Lecture - 30 Maxwell s Equations In the last lecture we had introduced
More informationdt 2 = 0, we find that: K = 1 2 Ω (2)
1 1. irial Theorem Last semester, we derived the irial theorem from essentially considering a series of particles which attract each other through gravitation. The result was that d = K + Ω (1) dt where
More informationIntermission Page 343, Griffith
Intermission Page 343, Griffith Chapter 8. Conservation Laws (Page 346, Griffith) Lecture : Electromagnetic Power Flow Flow of Electromagnetic Power Electromagnetic waves transport throughout space the
More informationPLASMA ASTROPHYSICS. ElisaBete M. de Gouveia Dal Pino IAG-USP. NOTES: (references therein)
PLASMA ASTROPHYSICS ElisaBete M. de Gouveia Dal Pino IAG-USP NOTES:http://www.astro.iag.usp.br/~dalpino (references therein) ICTP-SAIFR, October 7-18, 2013 Contents What is plasma? Why plasmas in astrophysics?
More informationThe structure and evolution of stars. Introduction and recap
The structure and evolution of stars Lecture 3: The equations of stellar structure 1 Introduction and recap For our stars which are isolated, static, and spherically symmetric there are four basic equations
More informationCHAPTER 7 ELECTRODYNAMICS
CHAPTER 7 ELECTRODYNAMICS Outlines 1. Electromotive Force 2. Electromagnetic Induction 3. Maxwell s Equations Michael Faraday James C. Maxwell 2 Summary of Electrostatics and Magnetostatics ρ/ε This semester,
More information