Conservation Laws in Ideal MHD
|
|
- Morris Payne
- 5 years ago
- Views:
Transcription
1 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 by Russell Kulsrud, Lectures in Magnetohydrodynamics by the late Dalton Schnack, Ideal Magnetohydrodynamics by Jeffrey Freidberg, Magnetic Reconnection by Eric Priest and Terry Forbes, and various papers by M. Berger & J. Perez. The description of Stokes Theorem and the Divergence Theorem come from Stewart s Multivariable Calculus with figures adapted from Wikimedia Commons and R. Fitzpatrick s online notes.
2 Outline Mathematical interlude: Stokes Theorem Divergence Theorem Dyadic tensors Conservation of mass, momentum, and energy The frozen-in condition Helicity and cross-helicity
3 Stokes Theorem Let S be an oriented piecewise-smooth surface that is bounded by a simple, close, piecewise-smooth boundary curve C with positive orientation. Let F be a vector field whose components have continuous partial derivatives on an open region in R 3 that contains S. Then F dr = F ds (1) C where ds = ˆn ds and ˆn is a unit normal vector to S S
4 The Divergence Theorem (also called Gauss s Theorem) Let V be a simple solid region and let S be the boundary surface of V, given with positive (outward) orientation. Let F be a vector field whose component functions have continuous partial derivatives on an open region that contains V. Then F ds = F dv (2) S The interior contributions to the divergence integral cancel, so only the exterior contributions remain V
5 While somewhat horrible, dyadic tensors 1 allow the equations of MHD to be written in a compact/useful form A dyadic tensor is a second rank tensor that can be written as the tensor product of two vectors: A ij = B i C j A = BC = (B 1 ê 1 + B 2 ê 2 + B 3 ê 3 ) (C 1 ê 1 + C 2 ê 2 + C 3 ê 3 ) = B 1 C 1 ê 1 ê 1 + B 1 C 2 ê 1 ê 2 + B 1 C 3 ê 1 ê 3 + B 2 C 1 ê 2 ê 1 + B 2 C 2 ê 2 ê 2 + B 2 C 3 ê 2 ê 3 + B 3 C 1 ê 3 ê 1 + B 3 C 2 ê 3 ê 2 + B 3 C 3 ê 3 ê 3 (3) where ê 1... are unit vectors and ê 1 ê 3... are unit dyads, e.g.: ê 1 ê 3 = (4) Good resources: Lectures in Magnetohydrodynamics by Dalton Schnack and rect.html
6 Dot products and double dot products with dyadic tensors In general, BC CB Dot products from the left and right are different: ê 1 A = B 1 C 1 ê 1 + B 1 C 2 ê 2 + B 1 C 3 ê 3 (5) A ê 1 = B 1 C 1 ê 1 + B 2 C 1 ê 2 + B 3 C 1 ê 3 (6) Double-dot notation is ambiguous but we use: A : T A ij T ij where we sum over repeated indices
7 The divergence of a dyadic tensor is a vector The divergence of a dyadic tensor is This is a vector whose jth component is T = i T ij (7) ( T) j = T 1j x 1 + T 2j x 2 + T 3j x 3 (8) We can even generalize Gauss theorem! T dv = ds T (9) S
8 What about the dyad T V? The gradient operator is given by ê i ê 1 + ê 2 + ê 3 (10) x i x 1 x 2 x 3 Then T V is given in Cartesian coordinates by ( V) ij = V j x i (11) V = V j ê i ê j (12) x i x V x x V y x V z = y V x y V y y V z (13) z V x z V y z V z where i is the row index and j is the column index
9 Conservative form As mentioned previously, conservative form is usually given by (stuff) + (flux of stuff) = 0 (14) t where source and sink terms go on the RHS If stuff is a scalar, then flux of stuff is a vector If stuff is a vector, then flux of stuff is a dyadic tensor If Eq. 14 is satisfied, then stuff is locally conserved
10 The continuity equation is already in conservative form The continuity equation 2 is ρ t where ρv is the mass flux + (ρv) = 0 (15) Mass is locally conserved the universe is safe once again! 2 This expression is similar to the Law of Conservation of Space Wombats
11 Is mass conserved globally? Integrate the continuity equation over a closed volume V bounded by a surface S [ ] ρ V t + (ρv) dv = 0 (16) ρ dv + (ρv) dv = 0 (17) t V V Now define M as the mass within V and use Gauss theorem dm + ds (ρv) = 0 (18) }{{} dt S }{{} change of mass in V mass flow through S Yep, we re still good. Mass is conserved globally. Recall that this was how we derived the continuity equation in the first place!
12 The momentum equation in conservative form The momentum equation can be written as ρv t + T = 0 (19) where the ideal MHD stress tensor T is ) T = ρvv + (p + B2 I BB 8π 4π and I ê 1 ê 1 + ê 2 ê 2 + ê 3 ê 3 is the identity dyadic tensor The quantity p + B 2 /8π is the total pressure (20) One can also include the viscous and gravitational stress tensors, or use a pressure tensor rather than a scalar pressure
13 Physical meaning of the stress tensor Imagine you have an infinitesimal box Forces are exerted on each face by the outside volume The forces on each side of the box could have components in three directions The stress tensor includes the nine quantities needed to describe the total force exerted on this box
14 The Reynolds stress tensor The Reynolds stress is given by T V ρvv (21) T V represents the flux of momentum. Think: T V ρv }{{} Momentum density V }{{} times a velocity (22) (ρv x )V y is the rate at which the x component of momentum is carried in the y direction (and vice versa)
15 The Maxwell stress tensor The Maxwell stress is T B B2 8π I BB 4π Momentum is transported by the magnetic field The quantity BB 4π motions is the hoop stress which resists shearing (23)
16 We can t always drop gravity! (Kulsrud, p. 77) Gravity is often treated as a body force where g is constant More generally, g is found using the gradient of a scalar potential, φ, g = φ (24) where Poisson s equation is used to solve for the potential, The gravitational stress tensor is 2 φ = 4πGρ (25) T g = ( φ)2 I 8πG where φ φ is a dyadic tensor φ φ 4πG The usual form of the force may be recovering using (26) ρg = T g (27)
17 CGL theory allows for p p For a scalar pressure, the corresponding part of the stress tensor is symmetric p 0 0 T p = pi = 0 p 0 (28) 0 0 p Chew, Goldberger, & Low (1956) allow for different perpendicular and parallel pressures in the double adiabatic approximation T CGL = p p p (29) Used more frequently in astrophysics than in laboratory or space plasma physics
18 Conservation of energy The energy equation can be written as The energy density is w t + s = 0 (30) The energy flux is w = ρv 2 }{{} 2 kinetic + B2 8π }{{} magnetic + p γ 1 }{{} internal (31) s = ρv 2 2 V + p }{{} γ 1 V }{{} kinetic internal + pv }{{} work + ce B 4π }{{} Poynting (32) pv represents the work done on the plasma from p forces
19 We now have a set of three equations describing conservation of mass, momentum, and energy Conservation of mass, momentum, & energy are given by ρ + (ρv) t = 0 (33) ρv + T t = 0 (34) w t + s = 0 (35) where the stress tensor, energy density, and energy flux are ) T = ρvv + (p + B2 I BB (36) 8π 4π w = ρv B2 8π + p (37) γ 1 ( ρv 2 s = 2 + γ ) γ 1 p V + E B (38) 4π
20 The most important property of ideal MHD is the frozen-in condition Define S as a surface bounded by a closed curve C that is co-moving with the plasma The magnetic flux through S is Ψ = B ds (39) There are contributions to dψ dt from: Changes in B with S held fixed, dψ1 dt The flux swept out by C as it moves with the plasma, dψ 2 dt S
21 How does the flux change due to B changing? Take the time derivative of the magnetic flux with S fixed: dψ 1 B = dt S t ds = c ( E) ds S = c E dl (40) where we used Faraday s law and Stokes theorem. C
22 How does the area of S change as it is swept with the plasma? C 2 V ds dl C 1 The lateral area swept out by C as it flows with the plasma is ds = Vdt dl (41) where the line element dl is tangent to C The flux through this area is B ds = B Vdt dl, so that dψ 2 = B V dl dt C = V B dl (42) C
23 Putting it all together The total change in magnetic flux through S is then dψ dt But E + V B c = dψ 1 + dψ 2 dt dt = c E dl C }{{} changes in B = c C ( E + V B c = 0 in ideal MHD, so that V B dl C }{{} changes in S ) dl (43) dψ dt = 0. (44) The magnetic flux through any co-moving fluid element remains constant in ideal MHD.
24 In ideal MHD, the magnetic field and plasma are frozen-in to each other If a parcel of plasma moves, the magnetic field attached to the parcel moves along with it More rigorously: if two plasma elements are initially connected by a magnetic field line, they will remain connected by a magnetic field line at future times. Magnetic topology (e.g., connectivity) is preserved in ideal MHD The plasma cannot move across magnetic field lines (though it remains free to move along field lines)
25 A useful concept in ideal MHD is a field line velocity More generally, flux will be frozen-in if there exists some velocity U such that E + U B c = 0 (45) The velocity U is a field line velocity. In ideal MHD, U is the bulk fluid velocity V A field line velocity does not need to be a fluid velocity In Hall MHD, the field line velocity is the electron velocity
26 What is meant by a field line velocity? In ideal MHD, a key field line velocity is V = c E B B 2 (46) This is the component of velocity perpendicular to B that the plasma is traveling at which is frozen into B The concept of a field line velocity is frought with peril. There is no way to distinguish one field line from another at different times There are an infinite number of field line velocities U is arbitrary (but in ideal MHD we usually pick V )
27 Eyink et al. (2013) argue that stochastic field line wandering in turbulence leads to Richardson diffusion that breaks down the frozen-in condition This result warrants further investigation by independent groups
28 Magnetic helicity measures the linkage of magnetic fields The magnetic field B can be written in terms of a vector potential B = A (47) while noting that A is not gauge invariant: A = A + φ The helicity of a magnetic field is given by H = A B dv (48) V A B should not be considered a helicity density because of gauge freedom! Helicity is approximately conserved during magnetic reconnection and topology changes Helicity can be injected into a system such as the solar corona. When too much builds up, helicity ends up being expelled through coronal mass ejections.
29 Helicity examples A single untwisted closed flux loop has H = 0 A single flux rope with a magnetic flux of Φ that twists around itself T times has a helicity of H = T Φ 2 Two interlinked untwisted flux loops with fluxes Φ 1 and Φ 2 have H = ±2Φ 1 Φ 2 where the sign depends on the sense of the linkedness
30 There are generalizations to allow for gauge-invariant definitions of helicity Berger & Field (1984) defined the relative magnetic helicity to be H = (A B A 0 B 0 ) dv (49) V where B 0 = A 0 is the potential field inside V with the same field outside of V (see also Finn & Antonsen 1985) In toroidal laboratory experiments, it is natural to consider the volume contained within conducting wall boundaries that are coincident with closed flux surfaces (i.e., the magnetic field along the wall is parallel to the boundary)
31 Helicity changes through parallel electric fields and helicity being added/removed to system The time evolution of magnetic helicity is given by dh dt = 2c E B dv + 2c A p E ds (50) V S where we choose A p = 0 and A p ds = 0 on S The first term represents helicity dissipation when E 0 But in ideal MHD, E = 0! The second term represents helicity fluxes in and out of system Flux emergence from the solar photosphere corresponds to helicity injection in the corona Helicity injection in laboratory experiments [e.g., the Steady Inductive Helicity Injected Torus (HIT-SI)]
32 The cross helicity H C measures the imbalance between interacting waves (important in MHD turbulence) The cross helicity is given by H C = In ideal MHD, the rate of change of H C is [( dh C 1 = ds dt S 2 V 2 + γ γ 1 V V B dv (51) ) ] p B V (V B) ρ (52) This vanishes when ds B = ds V = 0 along the boundary S or when the boundary conditions are periodic. Cross helicity is an ideal MHD invariant when this integral vanishes
33 Summary Ideal MHD conserves mass, momentum, and energy (as it darn well should!) In ideal MHD, the magnetic field is frozen-in to the plasma flow It is often useful to think in terms of a magnetic field line velocity, but there are caveats! Helicity and cross-helicity provide two topological constraints that are conserved in ideal MHD
Ideal 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 informationThe Virial Theorem, MHD Equilibria, and Force-Free Fields
The Virial Theorem, MHD Equilibria, and Force-Free Fields Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 10 12, 2014 These lecture notes are largely
More 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 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 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 informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
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 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 information13. ASTROPHYSICAL GAS DYNAMICS AND MHD Hydrodynamics
1 13. ASTROPHYSICAL GAS DYNAMICS AND MHD 13.1. Hydrodynamics Astrophysical fluids are complex, with a number of different components: neutral atoms and molecules, ions, dust grains (often charged), and
More informationTopological Methods in Fluid Dynamics
Topological Methods in Fluid Dynamics Gunnar Hornig Topologische Fluiddynamik Ruhr-Universität-Bochum IBZ, Februar 2002 Page 1 of 36 Collaborators: H. v. Bodecker, J. Kleimann, C. Mayer, E. Tassi, S.V.
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 informationMagnetic Reconnection in Laboratory, Astrophysical, and Space Plasmas
Magnetic Reconnection in Laboratory, Astrophysical, and Space Plasmas Nick Murphy Harvard-Smithsonian Center for Astrophysics namurphy@cfa.harvard.edu http://www.cfa.harvard.edu/ namurphy/ November 18,
More informationMacroscopic plasma description
Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion
More informationFundamentals of Magnetohydrodynamics (MHD)
Fundamentals of Magnetohydrodynamics (MHD) Thomas Neukirch School of Mathematics and Statistics University of St. Andrews STFC Advanced School U Dundee 2014 p.1/46 Motivation Solar Corona in EUV Want to
More informationCHAPTER 8 ENTROPY GENERATION AND TRANSPORT
CHAPTER 8 ENTROPY GENERATION AND TRANSPORT 8.1 CONVECTIVE FORM OF THE GIBBS EQUATION In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationMHD RELATED TO 2-FLUID THEORY, KINETIC THEORY AND MAGANETIC RECONNECTION
MHD RELATED TO 2-FLUID THEORY, KINETIC THEORY AND MAGANETIC RECONNECTION Marty Goldman University of Colorado Spring 2017 Physics 5150 Issues 2 How is MHD related to 2-fluid theory Level of MHD depends
More informationThe Divergence Theorem Stokes Theorem Applications of Vector Calculus. Calculus. Vector Calculus (III)
Calculus Vector Calculus (III) Outline 1 The Divergence Theorem 2 Stokes Theorem 3 Applications of Vector Calculus The Divergence Theorem (I) Recall that at the end of section 12.5, we had rewritten Green
More informationr r 1 r r 1 2 = q 1 p = qd and it points from the negative charge to the positive charge.
MP204, Important Equations page 1 Below is a list of important equations that we meet in our study of Electromagnetism in the MP204 module. For your exam, you are expected to understand all of these, and
More informationMHD turbulence in the solar corona and solar wind
MHD turbulence in the solar corona and solar wind Pablo Dmitruk Departamento de Física, FCEN, Universidad de Buenos Aires Motivations The role of MHD turbulence in several phenomena in space and solar
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 informationEECS 117. Lecture 17: Magnetic Forces/Torque, Faraday s Law. Prof. Niknejad. University of California, Berkeley
University of California, Berkeley EECS 117 Lecture 17 p. 1/? EECS 117 Lecture 17: Magnetic Forces/Torque, Faraday s Law Prof. Niknejad University of California, Berkeley University of California, Berkeley
More informationNovember 2, Monday. 17. Magnetic Energy Release
November, Monday 17. Magnetic Energy Release Magnetic Energy Release 1. Solar Energetic Phenomena. Energy Equation 3. Two Types of Magnetic Energy Release 4. Rapid Dissipation: Sweet s Mechanism 5. Petschek
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 informationComputational Astrophysics
Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture
More informationAccretion Disks I. High Energy Astrophysics: Accretion Disks I 1/60
Accretion Disks I References: Accretion Power in Astrophysics, J. Frank, A. King and D. Raine. High Energy Astrophysics, Vol. 2, M.S. Longair, Cambridge University Press Active Galactic Nuclei, J.H. Krolik,
More informationAA210A Fundamentals of Compressible Flow. Chapter 5 -The conservation equations
AA210A Fundamentals of Compressible Flow Chapter 5 -The conservation equations 1 5.1 Leibniz rule for differentiation of integrals Differentiation under the integral sign. According to the fundamental
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 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 informationToday in Physics 218: the classic conservation laws in electrodynamics
Today in Physics 28: the classic conservation laws in electrodynamics Poynting s theorem Energy conservation in electrodynamics The Maxwell stress tensor (which gets rather messy) Momentum conservation
More informationSpecial topic JPFR article Prospects of Research on Innovative Concepts in ITER Era contribution by M. Brown Section 5.2.2
Special topic JPFR article Prospects of Research on Innovative Concepts in ITER Era contribution by M. Brown Section 5.2.2 5.2.2 Dynamo and Reconnection Research: Overview: Spheromaks undergo a relaxation
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 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 informationE NCYCLOPEDIA OF A STRONOMY AND A STROPHYSICS
Magnetohydrodynamics Magnetohydrodynamics (or MHD for short) is the study of the interaction between a magnetic field and a plasma treated as a continuous medium (e.g. Cowling 1957, Roberts 1967, Priest
More information221A Miscellaneous Notes Continuity Equation
221A Miscellaneous Notes Continuity Equation 1 The Continuity Equation As I received questions about the midterm problems, I realized that some of you have a conceptual gap about the continuity equation.
More informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
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 informationAstronomy 404 October 18, 2013
Astronomy 404 October 18, 2013 Parker Wind Model Assumes an isothermal corona, simplified HSE Why does this model fail? Dynamic mass flow of particles from the corona, the system is not closed Re-write
More informationQuasi-separatrix layers and 3D reconnection diagnostics for linetied
Quasi-separatrix layers and 3D reconnection diagnostics for linetied tearing modes John M. Finn, LANL A. Steve Richardson, NRL CMSO, Oct 2011 Operated by Los Alamos National Security, LLC for the U.S.
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationBefore seeing some applications of vector calculus to Physics, we note that vector calculus is easy, because... There s only one Theorem!
16.10 Summary and Applications Before seeing some applications of vector calculus to Physics, we note that vector calculus is easy, because... There s only one Theorem! Green s, Stokes, and the Divergence
More informationLecture 3: 1. Lecture 3.
Lecture 3: 1 Lecture 3. Lecture 3: 2 Plan for today Summary of the key points of the last lecture. Review of vector and tensor products : the dot product (or inner product ) and the cross product (or vector
More informationCreation and destruction of magnetic fields
HAO/NCAR July 30 2007 Magnetic fields in the Universe Earth Magnetic field present for 3.5 10 9 years, much longer than Ohmic decay time ( 10 4 years) Strong variability on shorter time scales (10 3 years)
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationIntroduction to Electromagnetic Theory
Introduction to Electromagnetic Theory Lecture topics Laws of magnetism and electricity Meaning of Maxwell s equations Solution of Maxwell s equations Electromagnetic radiation: wave model James Clerk
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 informationMagnetic helicity of solar active regions
Magnetic helicity of solar active regions J. K. Thalmann Solar Group Seminar 12 Jul, 2011 Helicity of interlinked curves Number of enlacements of two interlinked curves l and l (Gauss, C. F., 1833, Werke,
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationNotes on Entropy Production in Multicomponent Fluids
Notes on Entropy Production in Multicomponent Fluids Robert F. Sekerka Updated January 2, 2001 from July 1993 Version Introduction We calculate the entropy production in a multicomponent fluid, allowing
More informationThe Navier-Stokes Equations
s University of New Hampshire February 22, 202 and equations describe the non-relativistic time evolution of mass and momentum in fluid substances. mass density field: ρ = ρ(t, x, y, z) velocity field:
More informationP. Démoulin and E. Pariat
Mem. S.A.It. Vol.?, 1 c SAIt 2004 Memorie della ÓÑÔÙØ Ò Ñ Ò Ø Ò Ö Ý Ò Ð ØÝ ÙÜ ÖÓÑ Ö Ó Ñ Ò ØÓ Ö Ñ P. Démoulin and E. Pariat Observatoire de Paris, section de Meudon, LESIA, UMR 8109 (CNRS), F-92195 Meudon
More informationTurbulent Magnetic Helicity Transport and the Rapid Growth of Large Scale Magnetic Fields
Turbulent Magnetic Helicity Transport and the Rapid Growth of Large Scale Magnetic Fields Jungyeon Cho Dmitry Shapovalov MWMF Madison, Wisconsin April 2012 The Large Scale Dynamo The accumulation of magnetic
More informationA Brief Revision of Vector Calculus and Maxwell s Equations
A Brief Revision of Vector Calculus and Maxwell s Equations Debapratim Ghosh Electronic Systems Group Department of Electrical Engineering Indian Institute of Technology Bombay e-mail: dghosh@ee.iitb.ac.in
More informationAsymmetric Magnetic Reconnection in Coronal Mass Ejection Current Sheets
Asymmetric Magnetic Reconnection in Coronal Mass Ejection Current Sheets Nicholas Murphy, 1 Mari Paz Miralles, 1 Crystal Pope, 1,2 John Raymond, 1 Kathy Reeves, 1 Dan Seaton, 3 & David Webb 4 1 Harvard-Smithsonian
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 informationPlasma Physics for Astrophysics
- ' ' * ' Plasma Physics for Astrophysics RUSSELL M. KULSRUD PRINCETON UNIVERSITY E;RESS '. ' PRINCETON AND OXFORD,, ', V. List of Figures Foreword by John N. Bahcall Preface Chapter 1. Introduction 1
More informationNavier-Stokes Equation: Principle of Conservation of Momentum
Navier-tokes Equation: Principle of Conservation of Momentum R. hankar ubramanian Department of Chemical and Biomolecular Engineering Clarkson University Newton formulated the principle of conservation
More informationBasics of MHD. Kandaswamy Subramanian a. Pune , India. a Inter-University Centre for Astronomy and Astrophysics,
Basics of MHD Kandaswamy Subramanian a a Inter-University Centre for Astronomy and Astrophysics, Pune 411 007, India. The magnetic Universe, Feb 16, 2015 p.0/27 Plan Magnetic fields in Astrophysics MHD
More informationMagnetic reconnection in coronal plasmas
UW, 28 May, 2010 p.1/17 Magnetic reconnection in coronal plasmas I.J.D Craig Department of Mathematics University of Waikato Hamilton New Zealand UW, 28 May, 2010 p.2/17 Why reconnection? Reconnection
More informationEuler equation and Navier-Stokes equation
Euler equation and Navier-Stokes equation WeiHan Hsiao a a Department of Physics, The University of Chicago E-mail: weihanhsiao@uchicago.edu ABSTRACT: This is the note prepared for the Kadanoff center
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 informationField line helicity as a tool for coronal physics
Field line helicity as a tool for coronal physics Anthony Yeates with Gunnar Hornig (Dundee), Marcus Page (Durham) Helicity Thinkshop, Tokyo, 20-Nov-2017 What is magnetic helicity? The average pairwise
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
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 informationOutline of Presentation. Magnetic Carpet Small-scale photospheric magnetic field of the quiet Sun. Evolution of Magnetic Carpet 12/07/2012
Outline of Presentation Karen Meyer 1 Duncan Mackay 1 Aad van Ballegooijen 2 Magnetic Carpet 2D Photospheric Model Non-Linear Force-Free Fields 3D Coronal Model Future Work Conclusions 1 University of
More information1. FUNDAMENTAL CONCEPTS AND MATH REVIEW
1. FUNDAMENTAL CONCEPTS AND MATH REVIEW 1.1. Introduction Here we provide for your reading pleasure a review of some of the math concepts used in part of this course. Most of this falls under the heading
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 informationFrom the last time, we ended with an expression for the energy equation. u = ρg u + (τ u) q (9.1)
Lecture 9 9. Administration None. 9. Continuation of energy equation From the last time, we ended with an expression for the energy equation ρ D (e + ) u = ρg u + (τ u) q (9.) Where ρg u changes in potential
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 informationWhere does Bernoulli's Equation come from?
Where does Bernoulli's Equation come from? Introduction By now, you have seen the following equation many times, using it to solve simple fluid problems. P ρ + v + gz = constant (along a streamline) This
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More information[variable] = units (or dimension) of variable.
Dimensional Analysis Zoe Wyatt wyatt.zoe@gmail.com with help from Emanuel Malek Understanding units usually makes physics much easier to understand. It also gives a good method of checking if an answer
More informationClass 4 : Maxwell s Equations for Electrostatics
Class 4 : Maxwell s Equations for Electrostatics Concept of charge density Maxwell s 1 st and 2 nd Equations Physical interpretation of divergence and curl How do we check whether a given vector field
More informationIntroduction to Magnetohydrodynamics (MHD)
Introduction to Magnetohydrodynamics (MHD) Tony Arber University of Warwick 4th SOLARNET Summer School on Solar MHD and Reconnection Aim Derivation of MHD equations from conservation laws Quasi-neutrality
More 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 informationブラックホール磁気圏での 磁気リコネクションの数値計算 熊本大学 小出眞路 RKKコンピュー 森野了悟 ターサービス(株) BHmag2012,名古屋大学,
RKK ( ) BHmag2012,, 2012.2.29 Outline Motivation and basis: Magnetic reconnection around astrophysical black holes Standard equations of resistive GRMHD Test calculations of resistive GRMHD A simulation
More informationPROBLEM SET. Heliophysics Summer School. July, 2013
PROBLEM SET Heliophysics Summer School July, 2013 Problem Set for Shocks and Particle Acceleration There is probably only time to attempt one or two of these questions. In the tutorial session discussion
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 informationLecture 3: Vectors. Any set of numbers that transform under a rotation the same way that a point in space does is called a vector.
Lecture 3: Vectors Any set of numbers that transform under a rotation the same way that a point in space does is called a vector i.e., A = λ A i ij j j In earlier courses, you may have learned that a vector
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 informationCreation and destruction of magnetic fields
HAO/NCAR July 20 2011 Magnetic fields in the Universe Earth Magnetic field present for 3.5 10 9 years, much longer than Ohmic decay time ( 10 4 years) Strong variability on shorter time scales (10 3 years)
More informationENGI 4430 Gauss & Stokes Theorems; Potentials Page 10.01
ENGI 443 Gauss & tokes heorems; Potentials Page.. Gauss Divergence heorem Let be a piecewise-smooth closed surface enclosing a volume in vector field. hen the net flux of F out of is F d F d, N 3 and let
More informationKeywords: 2015, IJARCSSE All Rights Reserved Page 315
Volume 5, Issue 3, March 2015 ISSN: 2277 128X International Journal of Advanced Research in Computer Science and Software Engineering Research Paper Available online at: www.ijarcsse.com Special Issue
More informationUnit-1 Electrostatics-1
1. Describe about Co-ordinate Systems. Co-ordinate Systems Unit-1 Electrostatics-1 In order to describe the spatial variations of the quantities, we require using appropriate coordinate system. A point
More informationSection 2: Lecture 1 Integral Form of the Conservation Equations for Compressible Flow
Section 2: Lecture 1 Integral Form of the Conservation Equations for Compressible Flow Anderson: Chapter 2 pp. 41-54 1 Equation of State: Section 1 Review p = R g T " > R g = R u M w - R u = 8314.4126
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 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 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 informationSOLAR MHD. S. S. Hasan Indian Institute of Astrophysics, Bangalore, India
SOLAR MHD S. S. Hasan (hasan@iiap.res.in) Indian Institute of Astrophysics, Bangalore, India Aim: The aim of the course is to provide a fundamental background in MHD to understand various processes that
More informationMagnetic helicity conservation in a solar active event
IUTAM Symposium: Helicity, Structures and Singularity in Fluid and Plasma Dynamics Venezia, Italy, April 11 th 2016 Magnetic helicity conservation in a solar active event É. Pariat 1, G. Valori 2, P. Démoulin
More informationCorresponding Author: Kandie K.Joseph. DOI: / Page
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 13, Issue 5 Ver. 1 (Sep. - Oct. 2017), PP 37-47 www.iosrjournals.org Solution of the Non-Linear Third Order Partial Differential
More informationLinearized Gravity Return to Linearized Field Equations
Physics 411 Lecture 28 Linearized Gravity Lecture 28 Physics 411 Classical Mechanics II November 7th, 2007 We have seen, in disguised form, the equations of linearized gravity. Now we will pick a gauge
More informationIntroduction to Plasma Physics
Introduction to Plasma Physics Hartmut Zohm Max-Planck-Institut für Plasmaphysik 85748 Garching DPG Advanced Physics School The Physics of ITER Bad Honnef, 22.09.2014 A simplistic view on a Fusion Power
More informationIII.3 Momentum balance: Euler and Navier Stokes equations
32 Fundamental equations of non-relativistic fluid dynamics.3 Momentum balance: Euler and Navier tokes equations For a closed system with total linear momentum P ~ with respect to a given reference frame
More informationLecture Notes in Astrophysical Fluid Dynamics
Lecture Notes in Astrophysical Fluid Dynamics Mattia Sormani November 19, 2017 1 Contents 1 Hydrodynamics 6 1.1 Introductory remarks.......................... 6 1.2 The state of a fluid...........................
More informationFluid Dynamics and Balance Equations for Reacting Flows
Fluid Dynamics and Balance Equations for Reacting Flows Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Balance Equations Basics: equations of continuum mechanics balance equations for mass and
More informationChapter 2 The Continuum Equations
Chapter 2 The Continuum Equations 2.1 The conservation of mass In solid mechanics and thermodynamics we keep track of a well defined fluid mass and this mass is usually trivially specified, as in the case
More informationLecture 3: The Navier-Stokes Equations: Topological aspects
Lecture 3: The Navier-Stokes Equations: Topological aspects September 9, 2015 1 Goal Topology is the branch of math wich studies shape-changing objects; objects which can transform one into another without
More informationComputational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Computational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture No. # 02 Conservation of Mass and Momentum: Continuity and
More information13. REDUCED MHD. Since the magnetic field is almost uniform and uni-directional, the field has one almost uniform component ( B z
13. REDUCED MHD One often encounters situations in which the magnetic field is strong and almost unidirectional. Since a constant field does not produce a current density, these fields are sometimes said
More information