Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
|
|
- Jasper Short
- 5 years ago
- Views:
Transcription
1 0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle, WA IEEE Pulsed Power and Plasma Sciences Conference June 2007 Albuquerque, NM
2 Abstract Simulations are performed using the full two-fluid plasma system. The twofluid plasma model is investigated for its capabilities of capturing physics that is lost with simpler fluid models such as Magnetohydrodynamics (MHD). The regime in between the full two-fluid equation system and MHD is explored by the application of asymptotic approximations to the two-fluid system. These asymptotic approximations involve ignoring electron inertia, setting the speed of light to infinity and ensuring charge neutrality. Applying all three approximations together gives Hall-MHD. The results obtained by applying the asymptotic approximations are compared to the two-fluid plasma model to determine what physics is lost by the application of each approximation. Simulations of electromagnetic plasma shock and collisionless reconnection will be presented to demonstrate the various physical effects in comparison to the results of the two-fluid model.
3 Motivation To study the two-fluid plasma model To study the physical effects when each of the asymptotic approximations are applied independently to the two-fluid plasma model To implement Hall-MHD and compare it to the two-fluid model The equation sets used are balance laws of the form, Q t + F(Q) = S(Q)
4 Two-Fluid Plasma Model Euler equations for ion and electron fluids in 1D for simplicity n i n i u i 0 n i u i n i u 2 i + p n i q i i m i (E x + v i B z w i B y ) n i v i n i u i v i n i q i m i (E y + w i B x u i B z ) n i w i n i u i w i n i q i m i (E z + u i B y v i B x ) e i t + (e i + p i )u i n n e x = i q i (E x u i + E y v i + E z w i ) n e u e 0 Mn e u e Mn e u 2 e + p e n eq e m i (E x + v e B z w e B y ) Mn e v e Mn e u e v e n eq e m i (E y + w e B x u e B z ) Mn e w e Mn e u e w e n eq e m i (E z + u e B y v e B x ) (e e + p e )u e n e q e (E x u e + E y v e + E z w e ) e e
5 The electromagnetic fields are described by Maxwell equations: B t + E = 0, E t c2 B = 1 ɛ 0 E = 1 ɛ 0 (q i n i + q e n e ) B = 0. s q s m s ρ s v s Writing Maxwell s equations in balance law form in 1D: E x 0 1 ε 0 (n i q i u i n e q e u e ) E y c 2 B z 1 ε 0 (n i q i v i n e q e v e ) E z + c 2 B y 1 ε = 0 (n i q i w i n e q e w e ). t B x x 0 0 B y E z 0 0 B z E y
6 The ion and electron energies e i and e e are M: electron-to-ion mass ratio n: number density u, v, w: velocity in x-,y-,z-directions p: pressure E: electric field B: magnetic field q i, q e : ion and electron charges c: speed of light e i = p i γ n i(u 2 i + vi 2 + wi 2 ) e e = p e γ Mn e(u 2 e + ve 2 + we) 2
7 Asymptotic Approximations to be applied to the Two-Fluid Model Each of the following approximations are applied independently. Charge neutrality, n i = n e No displacement currents, ε 0 0, c Ignore electron inertia, M 0 Applying all 3 of the above together gives Hall-MHD
8 Approximation 1: Apply Quasineutrality only, n i = n e The Two-Fluid equations remain unchanged except that there is only one continuity equation for both fluids. The dispersion relation for parallel propagation is ωpi 2 + ωpe 2 + ω 2 ω2 c 2 k 2 ω 2 c 2 = 0. (1) sik2 The dispersion relation for perpendicular propagation is where c si = T i mi. ωpi 2 + ωpe 2 + ω 2 ω2 c 2 k 2 ω 2 c 2 + ω 2 c 2 k 2 si k2 ω2 ci ω 2 c 2 = 0 (2) sik2
9 w k w 2 C 2 k w Left: Dispersion relation for the quasineutrality condition, ω Vs k. Right: Plot of the phase velocity squared, v 2 φ = ω2 k 2 Black: Parallel propagation, v φ = c si c for ω = 0 Vs ω Red: Perpendicular propagation similar to parallel with an additional dependence on ω ci. v 2 φ = c2 ω 2 ci +c2 si ω2 p c 2 (ω 2 ci +ωp2 ) for ω = 0, where ω2 p = ω 2 pi + ω2 pe
10 Approximation 2: Apply ε 0 0 only, implies c and no displacement currents The Electric field time derivative drops out of Ampere s law to give: B = µ 0 (n i qu i n e qu e ) (3) Quasineutrality is automatically implied through the math The dispersion relation for parallel propagation is which is just the ion acoustic wave. ω 2 = c 2 sk 2 (4) The dispersion relation for perpendicular propagation is where c 2 s = T i+t e m i +m e. ω 2 = c 2 sk 2 + 2ω ci ω ce (5)
11 w k w k w Left: Dispersion relation for the c condition, ω Vs k. Right: Plot of the phase velocity squared, v 2 φ = ω2 k 2 Black: Parallel propagation, ion acoustic wave Vs ω Red: Perpendicular propagation dispersion relation similar to parallel with an additional dependence on ω ci and a cutoff at ω = 2ω ce ω ci.
12 Approximation 3: Ignore electron inertia only, M 0 The time derivative and convective derivative terms drop out of the electron momentum equation to give Ohm s law: P e n e q = (E + u e B) (6) The massless electrons respond instantaneously to the ions in the direction of the magnetic field, so E =0.
13 The dispersion relation for parallel propagation for M 0 is ω 2 = c 2 sik 2. (7) The dispersion relation for perpendicular propagation is 1 = k2 c 2 ω 2 + ω2 pi ω 2 + ω2 pi ωci 2 k 2 T e 1 m i ω 2 ω4 pi ω 2 ci ω2 + 2k2 Te m i + ωci 2 + ω2 pi ω 2 k 2 c 2 si + ω2 ci + ω2 pi ( ω 4 pi ωci 2 ω2 pi ω2 ω 2 (8) ) where T e is electron temperature.
14 w w 2 C 2 k k w Left: Dispersion relation for the me mi 0 condition, ω Vs k Right: Plot of the phase velocity squared, v 2 φ = ω2 k 2 Black: Parallel propagation, ion sound wave Vs ω Red: Perpendicular propagation, asymptotes at the ion sound wave Blue: Perpendicular propagation, asymptotes at the speed of light
15 Applying all 3 asymptotic approximations gives Hall-MHD The dispersion relation for Hall-MHD is ( (ω 2 vak 2 [ω 2 ) 4 c 2 si + T ) e + va 2 k 2 ω 2 + m i = v2 A k2 k2 ω 2 n 0 eµ 0 ( c 2 si + T e m i ) v 2Ak 2 k2 ] ( ω 2 c 2 sik 2 T ) e k 2 m i (9) (10) where v A = B 0 n0 m i µ 0 is the Alfven speed.
16 Waves propagating parallel to the magnetic field are ω = ± 1 2 ω ciδ 2 i k 2 + where δ i is the ion skin depth. w 10 v 2 A k ω2 ci δ4 i k4, ω = c si k (11) k The Whistler wave (Green) grows quadratically without bound when δ i and ω ci are significant and Hall effects become important. If λ δ i >> 1 then the dispersion diagram sits near the origin and resembles Ideal MHD.
17 For perpendicular propagation, the dispersion relation gives ω = k va 2 + c2 s (12) where, c 2 s = T e+t i m i. This is the magnetosonic wave. The Whistler wave, the Alfven wave and the Magnetosonic wave need to be resolved in this system to capture the Hall effects.
18 Hall-MHD Equation system for cold electrons Written as Ideal MHD + source terms: ρ i t + (ρ iu i ) = 0 ρ i u i + (ρ i u i u i + P i I BB + B2 I) = 0 t µ 0 2µ 0 ε i t + ) [(ε i + P i + B2 u i (B u ] i) B = 1 2µ 0 µ 0 [ ] ( B) B B t + (u ib Bu i ) = n i eµ 0 eµ 2 0 ( ) B 2 n i ( B) Induction equation source term is tough to handle with shock problems such as the BrioWu shock. Need to either add artificial viscosity, or could use the electron energy equation with the P e source term in the induction equation.
19 BrioWu shock problem using Ideal MHD and the Two-Fluid Plasma Model 1 n Vs x n Vs x Two Fluid model MHD model number density number density x x Left: An MHD shock is initialized with a jump in number density, pressure and magnetic field Right: The Two-Fluid solution compared to the ideal MHD solution. Want to see how Hall-MHD and the other reduced fluid models compare to these results.
20 GEM challenge problem - Magnetic Reconnection Cannot simulate this in Ideal MHD because the fluid is frozen to the magnetic field and there are no non-ideal terms that allow the field lines to connect to each other The Hall term implies that the electrons are frozen to the magnetic field instead of the ions, so to get reconnection we need non-ideal terms like artificial viscosity, resistivity or an electron pressure gradient with the Cold Electron Hall-MHD model discussed. Two-Fluid GEM challenge is simulated, and the reduced fluid models will be compared to these results in the future
21 Two-Fluid Magnetic Reconnection initial conditions ( 1 ( n e = n i = n y ) ) sech2 λ (13) P e = 1 B 12µ 0n 2 e 0 (14) P i = 5 B 12µ 0n 2 i (15) 0 ( y B x = B 0 tanh + λ) ( ) ( ) 0 π 2πx πy cos sin (16) 10 L y B y = B π L x sin J ez = µ 0B 0 λ where λ = 0.5d i, d i = ( ) ( ) 2πx πy cos L x sech2 ( y λ ) c ω pi, L x = 25.6, L y L y L x L y (17) (18) = 12.8, with a grid resolution of
22 y y y x x x Two-Fluid Magnetic Reconnection: Left: Electron density is shown over characteristic transit times, t=0, 25, 50 Right: The blue dots represent the two-fluid reconnected flux, plot from [2]
23 Numerical Algorithms Used System to be solved: Q t + F x = S, where Q = Conserved variables, F = Flux, S = Source terms Finite Volume Method: High Resolution Wave Propagation method Finite Element Method: Runge-Kutta Discontinuous Galerkin method Both methods involve solving the Riemann problem across each cell interface
24 Finite Volume Method: High Resolution Wave Propagation Method First solve the hyperbolic, homogeneous equation using the 2 nd order update formula [3] Q n+1 i = Q n i t x t x Q t + F x = 0 [ A + Q i 1/2 + A Q i+1/2 ] ([ F] i+1/2 [ F] i 1/2 ) where the eigenvalues, eigenvectors and flux differences give: A ± Q i±1/2 : right- and left-going fluctuations at each cell interface, F] i±1/2 : correction fluxes at each cell interface to get to 2 nd order accuracy Limiters are applied to reduce the oscillations that form due to this higher order (Lax-Wendroff type scheme).
25 Source Term Handling Strang splitting is used. First solve ODE dq/dt = S over t/2 Solve homogeneous equation over t Again solve ODE over t/2 To advance the ODE dq/dt = S, a higher order Runge-Kutta solver is used with Strang splitting.
26 Finite Element Method: Discontinuous Galerkin Method Basic Scheme Apply polynomial basis functions to conserved quantities within each cell: Q(x) = Q r v r (x) r=0 where v r (x): basis functions locally defined within each cell Q r : expansion coefficients r: spatial order Balance law: t I i v r (x)qdx + I i v r (x) F x dx = I i v r (x)s.
27 Legendre polynomials chosen as basis functions (orthogonality property). The scheme becomes [4]: Q r t = F i+1/2 ( 1) r F i 1/2 x + 1 x 1 1 dp r (η) dη Fdη P r Sdη where F(η) F(Q(x(η), t)) and S(η) S(Q(x(η), t)). Interface fluxes computed using Riemann problem solver like with Wave Propagation method Integrals computed using Gaussian quadrature ODE dq r dt = RHS solved using Runge-Kutta time integration (RKDG scheme). Temporal order spatial order.
28 Limiters for RKDG scheme Coefficient of linear term modified using a min-mod limiter. Limiter limits order of algorithm in regions where it is applied. 2 types of limiters used: Characteristics-based limiters: Conserved variables transformed to characteristic variables and then apply limiter. Better accuracy. Component-based limiters: Directly apply limiters to conserved variables without transforming them. Faster, no Jacobian required For stability, CFL = 1/(2r 1), r: spatial order
29 Conclusions and Further Study The MHD shock problem is presented using Ideal MHD and the Two-Fluid Plasma Model The Magnetic Reconnection problem is presented using the Two-Fluid Model The dispersion relations are studied for each of the approximations that are independently applied to reduce the Two-Fluid Model. When all 3 approximations are applied together Hall-MHD is obtained. Hall-MHD is studied and will be implemented with either artificial viscosity or an electron pressure gradient and the electron energy equation. Hall-MHD and the other reduced fluid models will be used to simulate the MHD shock and Magnetic Reconnection problems and the results will be compared to Ideal MHD and the Two-Fluid Model.
30 References [1] D. Bale, R. J. LeVeque, S. Mitran, and J. A. Rossmanith. A wave-propagation method for conservation laws and balance laws with spatially varying flux functions. SIAM Journal of Scientific Computing, 24: , [2] M.A. Shay et al. Alfvenic collisionless magnetic reconnection and the hall term. Journal of Geophysical Research, 106(A3): , [3] A. Hakim, J. Loverich, and U. Shumlak. A high resolution wave propagation scheme for ideal Two-Fluid plasma equations. Journal of Computational Physics, [4] J. Loverich, A. Hakim, and U. Shumlak. A discontinous Galerkin method for ideal two-fluid plasma equations. Journal of Computational Physics, [5] J. Loverich and U. Shumlak. Approximate Riemann solver for the two-fluid plasma model. Journal of Computational Physics, 187: , [6] S. Ohsaki and S. Mahajan. Hall current and alfven wave. Physics of Plasmas, 11(3): , [7] P. Roe. Approximate Riemann solvers, parameter vectors, and difference schemes. Journal of Computational Physics, 43: , 1981.
A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1
A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1 B. Srinivasan 2, U. Shumlak Aerospace and Energetics Research Program University of Washington,
More informationA Study of 3-Dimensional Plasma Configurations using the Two-Fluid Plasma Model
A Study of 3-Dimensional Plasma Configurations using the Two-Fluid Plasma Model B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program University of Washington IEEE International Conference
More informationShock-capturing Schemes for a Collisionless Two-fluid Plasma Model
Shock-capturing Schemes for a Collisionless Two-fluid Plasma Model E. Alec Johnson Department of Mathematics, UW-Madison Presented on August 29, 2 at Sandia National Laboratories c 23, Kenneth Lang, Tufts
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 informationTwo-scale numerical solution of the electromagnetic two-fluid plasma-maxwell equations: Shock and soliton simulation
Mathematics and Computers in Simulation 76 (2007) 3 7 Two-scale numerical solution of the electromagnetic two-fluid plasma-maxwell equations: Shock and soliton simulation S. Baboolal a,, R. Bharuthram
More informationTransition From Single Fluid To Pure Electron MHD Regime Of Tearing Instability
Transition From Single Fluid To Pure Electron MHD Regime Of Tearing Instability V.V.Mirnov, C.C.Hegna, S.C.Prager APS DPP Meeting, October 27-31, 2003, Albuquerque NM Abstract In the most general case,
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 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 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 information20. Alfven waves. ([3], p ; [1], p ; Chen, Sec.4.18, p ) We have considered two types of waves in plasma:
Phys780: Plasma Physics Lecture 20. Alfven Waves. 1 20. Alfven waves ([3], p.233-239; [1], p.202-237; Chen, Sec.4.18, p.136-144) We have considered two types of waves in plasma: 1. electrostatic Langmuir
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 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 informationConstrained Transport Method for the Finite Volume Evolution Galerkin Schemes with Application in Astrophysics
Project work at the Department of Mathematics, TUHH Constrained Transport Method for the Finite Volume Evolution Galerkin Schemes with Application in Astrophysics Katja Baumbach April 4, 005 Supervisor:
More informationComputational Astrophysics
16 th Chris Engelbrecht Summer School, January 2005 3: 1 Computational Astrophysics Lecture 3: Magnetic fields Paul Ricker University of Illinois at Urbana-Champaign National Center for Supercomputing
More informationDispersive Media, Lecture 7 - Thomas Johnson 1. Waves in plasmas. T. Johnson
2017-02-14 Dispersive Media, Lecture 7 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasmas as a coupled system Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas
More information2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1. Waves in plasmas. T. Johnson
2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasma physics Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas Transverse waves
More informationFast Magnetic Reconnection in Fluid Models of (Pair) Plasma
Fast Magnetic Reconnection in Fluid Models of (Pair) Plasma E. Alec Johnson Department of Mathematics, UW-Madison Presented on September 10, 2009, Postdoctoral Research Symposium, Argonne National Laboratories.
More informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
More informationREPORT DOCUMENTATION PAGE
REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationUNCERTAINTY QUANTIFICATION OF THE GEM CHALLENGE MAGNETIC RECONNECTION PROBLEM USING THE MULTILEVEL MONTE CARLO METHOD
International Journal for Uncertainty Quantification, 5 (4): 327 339 (2015) UNCERTAINTY QUANTIFICATION OF THE GEM CHALLENGE MAGNETIC RECONNECTION PROBLEM USING THE MULTILEVEL MONTE CARLO METHOD Éder M.
More informationA high resolution wave propagation scheme for ideal Two-Fluid plasma equations
Journal of Computational Physics 29 (26) 48 442 www.elsevier.com/locate/jcp A high resolution wave propagation scheme for ideal Two-Fluid plasma equations A. Hakim *, J. Loverich, U. Shumlak Aerospace
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationRadiative & Magnetohydrodynamic Shocks
Chapter 4 Radiative & Magnetohydrodynamic Shocks I have been dealing, so far, with non-radiative shocks. Since, as we have seen, a shock raises the density and temperature of the gas, it is quite likely,
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 informationPhase Plane Behavior of Solitary Waves in Nonlinear Layered Media
Phase Plane Behavior of Solitary Waves in Nonlinear Layered Media Randall J. LeVeque 1 and Darryl H. Yong 2 1 Department of Applied Mathematics, University of Washington, Box 35242, Seattle, WA 98195-242
More informationSelf-Organization of Plasmas with Flows
Self-Organization of Plasmas with Flows ICNSP 2003/ 9/10 Graduate School of Frontier Sciences,, National Institute for Fusion Science R. NUMATA, Z. YOSHIDA, T. HAYASHI ICNSP 2003/ 9/10 p.1/14 Abstract
More informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More informationA Finite Volume Code for 1D Gas Dynamics
A Finite Volume Code for 1D Gas Dynamics Michael Lavell Department of Applied Mathematics and Statistics 1 Introduction A finite volume code is constructed to solve conservative systems, such as Euler
More informationPlasma waves in the fluid picture I
Plasma waves in the fluid picture I Langmuir oscillations and waves Ion-acoustic waves Debye length Ordinary electromagnetic waves General wave equation General dispersion equation Dielectric response
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationChapter 9 WAVES IN COLD MAGNETIZED PLASMA. 9.1 Introduction. 9.2 The Wave Equation
Chapter 9 WAVES IN COLD MAGNETIZED PLASMA 9.1 Introduction For this treatment, we will regard the plasma as a cold magnetofluid with an associated dielectric constant. We then derive a wave equation using
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 informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationVlasov-Maxwell Equations and Cold Plasma Waves
Astronomy 53 Spring 016) uening Bai Mar.,, 016 Vlasov-Maxwell Equations and Cold Plasma Waves The Vlasov-Maxwell equations Consider a plasma as a collection of N charged particles, each particle i has
More informationHigh Resolution Wave Propagation Schemes for Two-Fluid Plasma Simulations
High Resolution Wave Propagation Schemes for Two-Fluid Plasma Simulations Ammar H. Hakim A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy University
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 informationHeating and current drive: Radio Frequency
Heating and current drive: Radio Frequency Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 13 th February 2012 Dr Ben Dudson Magnetic Confinement Fusion (1 of 26)
More 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 informationThe importance of including XMHD physics in HED codes
The importance of including XMHD physics in HED codes Charles E. Seyler, Laboratory of Plasma Studies, School of Electrical and Computer Engineering, Cornell University Collaborators: Nat Hamlin (Cornell)
More informationMulti-Temperature, Thermal & Ion Fraction Effects over Wedge and Bluff Body Shapes in Hypervelocity Flow. Ward W. Vuillemot, Uri Shumlak
Multi-Temperature, Thermal & Ion Fraction Effects over Wedge and Bluff Body Shapes in ypervelocity Flow Ward W. Vuillemot, Uri Shumlak 1 University of Washington - Computational Fluid Dynamics Laboratory
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationCalculating Storm Surge and Other Coastal Hazards Using Geoclaw
Calculating Storm Surge and Other Coastal Hazards Using Geoclaw Kyle T. Mandli Department of Applied Mathematics University of Washington Seattle, WA, USA Modeling and Computations of Shallow-Water Coastal
More informationA Numerical Study of Magnetic Reconnection: A Central Scheme for Hall MHD
Proceedings of Symposia in Applied Mathematics A Numerical Study of Magnetic Reconnection: A Central Scheme for Hall MHD Xin Qian, Jorge Balbás, Amitava Bhattacharjee, and Hongang Yang Abstract. Over the
More informationAMath 574 February 11, 2011
AMath 574 February 11, 2011 Today: Entropy conditions and functions Lax-Wendroff theorem Wednesday February 23: Nonlinear systems Reading: Chapter 13 R.J. LeVeque, University of Washington AMath 574, February
More informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
More informationBasic plasma physics
Basic plasma physics SPAT PG Lectures Jonathan Eastwood 10-14 October 2016 Aims Provide new PhD students in SPAT and the SPC section with an overview of the most important principles in space plasma physics,
More informationSpace Physics. ELEC-E4520 (5 cr) Teacher: Esa Kallio Assistant: Markku Alho and Riku Järvinen. Aalto University School of Electrical Engineering
Space Physics ELEC-E4520 (5 cr) Teacher: Esa Kallio Assistant: Markku Alho and Riku Järvinen Aalto University School of Electrical Engineering The 6 th week: topics Last week: Examples of waves MHD: Examples
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 informationMagnetospheric Physics - Final Exam - Solutions 05/07/2008
Magnetospheric Physics - Final Exam - Solutions 5/7/8. Dipole magnetic field a Assume the magnetic field of the Earth to be dipolar. Consider a flux tube with a small quadratic cross section in the equatorial
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 informationEnergy Stable Discontinuous Galerkin Methods for Maxwell s Equations in Nonlinear Optical Media
Energy Stable Discontinuous Galerkin Methods for Maxwell s Equations in Nonlinear Optical Media Yingda Cheng Michigan State University Computational Aspects of Time Dependent Electromagnetic Wave Problems
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 informationA Three-Fluid Approach to Model Coupling of Solar Wind-Magnetosphere-Ionosphere- Thermosphere
A Three-Fluid Approach to Model Coupling of Solar Wind-Magnetosphere-Ionosphere- Thermosphere P. Song Center for Atmospheric Research University of Massachusetts Lowell V. M. Vasyliūnas Max-Planck-Institut
More informationLesson 3: MHD reconnec.on, MHD currents
Lesson3:MHDreconnec.on, MHDcurrents AGF 351 Op.calmethodsinauroralphysicsresearch UNIS,24. 25.11.2011 AnitaAikio UniversityofOulu Finland Photo:J.Jussila MHDbasics MHD cannot address discrete or single
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
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 informationGyrokinetics an efficient framework for studying turbulence and reconnection in magnetized plasmas
Frank Jenko Gyrokinetics an efficient framework for studying turbulence and reconnection in magnetized plasmas Max-Planck-Institut für Plasmaphysik, Garching Workshop on Vlasov-Maxwell Kinetics WPI, Vienna,
More informationMagnetohydrodynamic waves in a plasma
Department of Physics Seminar 1b Magnetohydrodynamic waves in a plasma Author: Janez Kokalj Advisor: prof. dr. Tomaž Gyergyek Petelinje, April 2016 Abstract Plasma can sustain different wave phenomena.
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More informationA Numerical Study of Magnetic Reconnection: A Central Scheme for Hall MHD
Proceedings of Symposia in Applied Mathematics A Numerical Study of Magnetic Reconnection: A Central Scheme for Hall MHD Xin Qian, Jorge Balbás, Amitava Bhattacharjee, and Hongang Yang Abstract. Over the
More informationNumerical simulation of cylindrical solitary waves in periodic media
Numerical simulation of cylindrical solitary waves in periodic media Manuel Quezada de Luna David I. Ketcheson June, 3 Abstract We study the behavior of nonlinear waves in a two-dimensional medium with
More informationEnergy-Conserving Numerical Simulations of Electron Holes in Two-Species Plasmas
Energy-Conserving Numerical Simulations of Electron Holes in Two-Species Plasmas Yingda Cheng Andrew J. Christlieb Xinghui Zhong March 18, 2014 Abstract In this paper, we apply our recently developed energy-conserving
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 informationFinite volumes for complex applications In this paper, we study finite-volume methods for balance laws. In particular, we focus on Godunov-type centra
Semi-discrete central schemes for balance laws. Application to the Broadwell model. Alexander Kurganov * *Department of Mathematics, Tulane University, 683 St. Charles Ave., New Orleans, LA 708, USA kurganov@math.tulane.edu
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationHybrid Simulations: Numerical Details and Current Applications
Hybrid Simulations: Numerical Details and Current Applications Dietmar Krauss-Varban and numerous collaborators Space Sciences Laboratory, UC Berkeley, USA Boulder, 07/25/2008 Content 1. Heliospheric/Space
More informationHybrid Simulation Method ISSS-10 Banff 2011
Hybrid Simulation Method ISSS-10 Banff 2011 David Burgess Astronomy Unit Queen Mary University of London With thanks to Dietmar Krauss-Varban Space Plasmas: From Sun to Earth Space Plasma Plasma is (mostly)
More informationTwo-layer shallow water system and its applications
Proceedings of Symposia in Applied Mathematics Two-layer shallow water system and its applications Jihwan Kim and Randall J. LeVeque Abstract. The multi-layer shallow water system is derived by depth averaging
More informationResistive MHD, reconnection and resistive tearing modes
DRAFT 1 Resistive MHD, reconnection and resistive tearing modes Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3NP, UK (This version is of 6 May 18 1. Introduction
More informationCold plasma waves. Waves in non-magnetized plasma Cold plasma dispersion equation Cold plasma wave modes
Cold plasma waves Waves in non-magnetized plasma Cold plasma dispersion equation Cold plasma wave modes EM wave propagation through and interaction with plasmas belong to central issues of plasma physics.
More informationDevelopment of Electromagnetic Solvers for Use with The Two-Fluid Plasma Algorithm
Development of Electromagnetic Solvers for Use with The Two-Fluid Plasma Algorithm Andree Susanto A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Aeronautics
More informationKinetic, Fluid & MHD Theories
Lecture 2 Kinetic, Fluid & MHD Theories The Vlasov equations are introduced as a starting point for both kinetic theory and fluid theory in a plasma. The equations of fluid theory are derived by taking
More informationSolution of time-dependent Boltzmann equation for electrons in non-thermal plasma
Solution of time-dependent Boltzmann equation for electrons in non-thermal plasma Z. Bonaventura, D. Trunec Department of Physical Electronics Faculty of Science Masaryk University Kotlářská 2, Brno, CZ-61137,
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
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 informationarxiv: v2 [math.na] 4 Mar 2014
Finite Difference Weighted Essentially Non-Oscillatory Schemes with Constrained Transport for Ideal Magnetohydrodynamics arxiv:1309.3344v [math.na] 4 Mar 014 Andrew J. Christlieb a, James A. Rossmanith
More informationStable Numerical Scheme for the Magnetic Induction Equation with Hall Effect
Stable Numerical Scheme for the Magnetic Induction Equation with Hall Effect Paolo Corti joint work with Siddhartha Mishra ETH Zurich, Seminar for Applied Mathematics 17-19th August 2011, Pro*Doc Retreat,
More informationThe Center for Astrophysical Thermonuclear Flashes. FLASH Hydrodynamics
The Center for Astrophysical Thermonuclear Flashes FLASH Hydrodynamics Jonathan Dursi (CITA), Alan Calder (FLASH) B. Fryxell, T. Linde, A. Mignone, G. Wiers Many others! Mar 23, 2005 An Advanced Simulation
More informationThe RAMSES code and related techniques 2- MHD solvers
The RAMSES code and related techniques 2- MHD solvers Outline - The ideal MHD equations - Godunov method for 1D MHD equations - Ideal MHD in multiple dimensions - Cell-centered variables: divergence B
More informationTarget Simulations. Roman Samulyak in collaboration with Y. Prykarpatskyy, T. Lu
Muon Collider/Neutrino Factory Collaboration Meeting May 26 28, CERN, Geneva U.S. Department of Energy Target Simulations Roman Samulyak in collaboration with Y. Prykarpatskyy, T. Lu Center for Data Intensive
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationWARPX - Washington Approximate Riemann Plasma Solver
WARPX - Washington Approximate Riemann Plasma Solver Current and Past Developers: S.Gilliam, A.Hakim, R.Lilly, N.Reddell, B.Srinivasan, A.Susanto Computational Plasma Dynamics Lab University of Washington,
More informationOverview of FRC-related modeling (July 2014-present)
Overview of FRC-related modeling (July 2014-present) Artan Qerushi AFRL-UCLA Basic Research Collaboration Workshop January 20th, 2015 AFTC PA Release# 15009, 16 Jan 2015 Artan Qerushi (AFRL) FRC modeling
More informationFinite Element Solver for Flux-Source Equations
Finite Element Solver for Flux-Source Equations Weston B. Lowrie A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Aeronautics Astronautics University
More informationMagnetic Reconnection: explosions in space and astrophysical plasma. J. F. Drake University of Maryland
Magnetic Reconnection: explosions in space and astrophysical plasma J. F. Drake University of Maryland Magnetic Energy Dissipation in the Universe The conversion of magnetic energy to heat and high speed
More informationSimulation results for magnetized plasmas
Chapter 4 Simulation results for magnetized plasmas In this chapter, we consider the dust charge fluctuation mode and lower hybrid wave damping in a magnetized plasma. Also, we consider plasma instabilities
More informationOn the Comparison of the Finite Volume and Discontinuous Galerkin Methods
Diploma Thesis Institute for Numerical Simulation, TUHH On the Comparison of the Finite Volume and Discontinuous Galerkin Methods Corrected version Katja Baumbach August 17, 2006 Supervisor: Prof. Dr.
More informationOn divergence-free reconstruction schemes for CED and MHD
On divergence-free reconstruction schemes for CED and MHD Praveen Chandrashekar praveen@math.tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore-560065, India
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationThe RAMSES code and related techniques 4. Source terms
The RAMSES code and related techniques 4. Source terms Outline - Optically thin radiative hydrodynamics - Relaxation towards the diffusion limit - Hydrodynamics with gravity source term - Relaxation towards
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 informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
More informationHydrodynamic modes of conducting liquid in random magnetic field
arxiv:1602.08543v1 [physics.plasm-ph] 27 Feb 2016 Hydrodynamic modes of conducting liquid in random magnetic field A. A. Stupka November 7, 2017 Oles Honchar Dnipropetrovs k National University, Gagarin
More informationProjection Dynamics in Godunov-Type Schemes
JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and
More informationAircraft Dynamics First order and Second order system
Aircraft Dynamics First order and Second order system Prepared by A.Kaviyarasu Assistant Professor Department of Aerospace Engineering Madras Institute Of Technology Chromepet, Chennai Aircraft dynamic
More informationElectromagneic Waves in a non- Maxwellian Dusty Plasma
Electromagneic Waves in a non- Maxwellian Dusty Plasma Nazish Rubab PhD student, KF University Graz IWF-OEAW Graz 26 January, 2011 Layout Dusty Plasma Non-Maxwellian Plasma Kinetic Alfven Waves Instability
More information