Computational Modeling of Fully Ionized Magnetized Plasmas Using the Fluid Approximation
|
|
- Christopher Gallagher
- 5 years ago
- Views:
Transcription
1 Computational Modeling of Fully Ionized Magnetized Plasmas Using the Fluid Approximation Dalton D. Schnack Center for Energy and Space Science Center for Magnetic Self-Organization in Space and Laboratory Plasmas Science Applications International Corp. San Diego, CA USA
2 Co-conspirators D. C. Barnes, Univ. Colorado, Boulder D. P. Brennan, MIT and GA C. Hegna, Univ. Wisconsin, Madison E. Held, Utah State Univ., Logan C. Kim, Univ. Washington, Seattle S. E. Kruger, TechX Corp., Boulder A. Y. Pankin, SAIC, San Diego C. R. Sovinec, Univ. Wisconsin, Madison + The NIMROD Team
3 Acknowledgements E. Belova A. Glasser S. Jardin J-N. Leboeuf J. Ramos L. Sugiyama CEMM Pioneering work: Theory J. Callen Multi-level modeling G. Fu W. Park S. Parker H. Strauss
4 Special Acknowledgement Prof. John Killeen Pioneer in computational physics Proponent of using implicit methods and primitive equations Founding Director of NERSC (CTRCC) Leader and mentor
5 The Problem Compute the low frequency dynamics of hot magnetized plasmas in realistic geometry in the presence of high frequency oscillations Incorporate the effects of lowest order kinetic corrections to the usual MHD equations Develop accurate and efficient algorithms that enable these goals
6 Modeling Magnetized Plasmas Plasma kinetic equation d dt f " (x,v,t) = #f " + v$ %f " + q " #t Maxwell s equations m " "B "t = #$ % E $ & E = ' q ( E + v & B)$ % v f " = ( C ",' ( f ", f ' ) ' q = * q ( ) f ( d 3 v $ % B # "E "t = µ 0J J = * q ( ) f ( vd 3 v ( ( Contains all information about plasma dynamics Impossible to solve analytically except in special cases Impractical for low frequencies, global geometry '
7 Fluid equations defined by taking moments of distribution function Define moments of distribution function M n (x, t) = # $ f (x, v, t)v n dv "# Knowledge of N moments allows (in principle) reconstruction of f at N points in velocity space N moments of plasma kinetic equation => N fluid equations satisfied by M N+1 Each additional moment equation yields more information about velocity distribution Must truncate moment equation hierarchy Approximate solution of kinetic equation
8 Closure of Moment Equations Use low-order truncation and closures Need to express high-order moments in terms of low-order moments q = q[ n,t,... ], " = "( p, V,...) Must be obtained from approximate solution of kinetic equation Analytical Numerical There is no general agreement on the closure of the moment equations for hot, magnetized plasmas!
9 Two-fluid Equations (m e ~ 0, n e ~ n i ) Lowest order moments for ions and electrons: "n "t = #$ % nv e = #$ % nv i = #$p i + ne( E + V i & B) # $ % ' i + R mn dv i dt 0 = #$p e # ne E + V e & B + Energy equation (+??) ( ) # $ % ' e # R ( ) J = ne V i # V e + Maxwell s equations (V 2 /c 2 <<1) + Closure expressions
10 Two-Fluid Equations Present Challenges for Computation Extreme separation of time scales { " A Alfvén transit time < { " S Sound transit time << " evol MHD evolution time { << { " R Extreme separation of spatial scales Internal boundary layers, localized and extended along magnetic field lines δ /L ~ S -α << 1 for S >> 1 Extreme anisotropy E.g., accurate treatment of B., χ /χ ~ 10 10, etc. Parasitic modes S = τ R /τ A Resistive diffusion time High frequency modes inherent in the formulation that may affect the low frequency dynamics
11 Dealing with Parasitic Modes The fundamental problem of computational MHD: Compute low frequency dynamics in presence of high frequency parasitic modes Reduction of mathematical model Eliminate parasitic modes analytically Example: V = 0 eliminates sound waves Strong toroidal field allows elimination of fast waves from MHD model Primitive equations and strongly implicit methods No analytic reduction of equations Use algorithms that allow very large time steps (CFL ~ ) Will concentrate on the second approach
12 There are Different Fluid Models Within fluid formulation, different terms are important in different parameter regimes Leads to different fluid models of plasmas MHD Hall MHD Drift MHD Transport Models distinguished by degree of force balance Obtained by non-dimensionalizing equations and systematically ordering small parameters: " = # i / L << 1, $ = % / & i, ' = V /V thi
13 Non-dimensional Equations Continuity: Ion momentum: Electron momentum: Pre-Maxwell: Orderings: Normalizations: "% #V i #t " #n #t = $%&'( nv i = $%&'( nv e +% 2 &V i ( 'V i = $ 1 n & * 'p i + ) i0 -, '( ) i / +% ( E+ V i 0 B), + p 0. %E= $%V e 0 B$ 1 n & * 'p e + ) e0 -, '( ) e / + p 0. " #B #t = $%&'0 E, J =%'0 B, J = n ( V i $ V e ) " = 1 2, % = V 0, & = 3 i V thi L << 1 time flow length 2 E 0 =V 0 B 0, J 0 = n 0 ev 0, p 0 = mn 0 V thi
14 Equation of Motion and Generalized Ohm s Law Adding and subtracting ion and electron equations: "J # B $ 1 n %&p * = n '" (V i +" 2 -, %V i ) &V i / $ (t. n % 0 i0 & ) 0 i p "Equilibrium" forces Dynamical response "( E + V i # B) = " 1 n J # B $% 1 ) n &p e + ' e0, + & ( ' e. * p Ideal MHD 2- fluid and FLR effects V B and J B enter at same order in ξ
15 Stress Tensor Scaling " = " + "^ + " # " = bb# " "^ = (I $ b) # " " % = (I & bb) # " Component Scaling Remarks " / p 0 (Braginskii) #$ /(% / &) D iverges for low collisionality ( % / & ~ $ 2 ) " / p 0 (Neo-classical) " ' / p 0 (# /$)(% / &) #$(% / &) "^ / p 0 (Gyro-viscosity) #$ O(#$) at low collisionality R emains in scale Vanishingly small at low collisionality Ignore Independent of collisionality Not dissipative Important FLR effects
16 Different Orderings Yield Different Fluid Models "J # B $ 1 n %&p * = n '" (V i + " 2 -, %V i ) &V i / $ (t. n % 0 i0 & ) 0 i p " Equilibrium" forces Dynamical response "( E + V i # B) = " 1 n J # B $% 1 ) n &p e + ' e0, + & ( ' e. * p 0 - Ideal MHD fluid and FLR effects Model V Force Balance Ohm s Law J $ B = n dv i dt E + V i $ B = Hall MHD V th /" # ci MHD V th "# ci Drift MHD "V th " 2 # ci + 1 n "2 (%p + % & ' gv ) + O(" 3 ) ( J $ B = " n dv i + * + %p- ) dt, + " 2 % & ' i gv + O(" 4 )..%p + J $ B = " 2 ( n dv i gv+ * + % & ' ) dt i -, + O(" 4 ) 1 n J $ B + O(" 2 ) E + V i $ B = 1 n J 12 $ 3 B." 1 n %p e O(") = O(") E + V i $ B = 1 ( n J $ B. %p e)
17 The Standard Drift Model In MHD, V i = E B/ B 2 = V E In drift ordering, V i = V E + V * + O(δ 2 ) Write drift equations in terms of V E : Ohm s Law $ E = " V E + V *i " 1 n J ' & #) * B " 1 % ( n +p e + O(, 2 ), ( ) = "V E * B " 1 n + p e + 1 n "+ # p + J * B O(,2 ), O(, 2 ) = "V E * B " 1 n + p e Valid only for slight deviations from equilibrium & " 2 ( ( n d ( dt ' Equation of Motion ( ) + n dv *i V i + V E ) gv + + # $ % dt i ( Vi ) + = * GVC, #p + J - B + O(" 4 ) Gyro-viscous cancellation gives simplified equations Exact form uncertain Only applicable to slight deviations from equilibrium We ignore for general application
18 Extended MHD Model Mn dv dt = "#p + J $ B " # % & i " # % & gvi E = "V $ B + 1 ( ne J $ B " #p e " # % & e ) + 'J (B (t = "# $ E, µ 0J = # $ B + Continuity and Energy equations + Closure expressions Encompasses Hall, MHD, and Drift models Terms can be selected by the user GV cancellation not explicitly implemented
19 Extended MHD Properties Dispersion Contains all MHD modes (ω 2 ~ k 2 ) Introduces dispersive modes (ω 2 ~ k 4 ) Electrons (whistlers) Ions + electrons (kinetic Alfvén wave) Ions only ( gyro-viscous waves) If extended MHD just produced more troublesome parasitic modes, who cares? However.. Stability Drift stabilization at moderate to high k Neo-classical de-stabilization of magnetic islands ++++?
20 Dispersion in Extended MHD Mode Origin Wave Equation Dispersion Comments Whistler J " B in Ohm # 2 B & #t 2 = $ V 2 ) A ( ' % + * 2 ( b, -) 2-2 B. 2 2 = V Ak / 0 ik 2 ( ) finite k electron response KAW - p e in Ohm # 2 B #t 2 = & V AV th* ) ( + ' % * 2 ( b, -) 2 [ ] - " bb, - " B = V 1 Ak 1+ 0s k 23 7 ( ) finite k, k 7 ion and electron response Parallel ion GV Perp. ion GV 8 4 term in -, 9 GV 0 #2 V 7 #t 2 = $ V7 8 3 term in -, 9 GV 0 #V 7 #t 1. L± = V A k ±1+ 1+ / 4 3 ( 0 i k ) / 5 1. R± = V A k ±1$ 1+ / 4 3 ( 0 i k ) / 5 = $ V :/ = V Ak / 16 0 ik 7 2 ( ) finite k ion response finite k 7 ion response Notation: 0 i = V th i / % is the ion gyro-radius; V th* = T e / m i ; 0 s = V th* / %; 8 4 = nt i / 2%; 8 3 = 28 4 ω 2 ~ k 4 Δt ~ Δx 2 Requires implicit methods
21 x Stability: Gravitational Interchange K. V. Roberts and J. B. Taylor, Phys. Rev. Letters 8, 197 (1962). "# "t + $ % #V = 0 ρ k B g y (" # $ ) = % & x 0 ' 0 # dv ' * B2 = &$ p + dt ) ( 2µ, + #g & $ % E = &V. B + M / #e # dv dt + $p 2 1 i & #g + $ % Assume electrostatic: " # E = 0 $ Gyro-viscosity: ( ) ( ikv x + & 0 ' 0 k 2 V y (" # $ ) = % & y ( 0 ' 0 ) ( ikv y % & 0 ' 0 k 2 V x " % V + 1 & " # dv dt ' 1 &( 2 "( # " % ) Extended MHD = 0
22 G-mode stabilization Dispersion Relation Solution Stabilizing Wave Number MHD (" = 0, # = 0) $ 2 + g% = 0 $ = i g% None 2-Fluid (" = 0, # = 1) Gyro-Viscosity (" = 1, # = 0) $ 2 & gk ' 0 $ + g% = 0 $ 2 &. 0 %k$ + g% = 0 2$ = gk ' 0 ± 2$ =. 0 %k ( gk + * - ), ' 0 2 & 4g% ± (. 0 %k) 2 & 4g% k 2 > 4%' 0 2 g k 2 > 4g. 0 2 % Full Extended MHD (" =1, # =1) $ 2 ( & gk + * +. 0 %k - $ + g% = 0 )' 0, 2$ = gk ' %k ± ( gk + * +. 0 %k - )' 0, 2 & 4g% k 2 > 4g% ( g + * +. 0 %- )' 0, 2
23 Form of the Gyro-viscous Stress (Hooke s Law for a Magnetized Plasma) Braginskii: " ^ = " gv = p ( 4# b$ W )% I + 3bb [ ( ) + transpose] W = "V + "V T # 2 3 I"$ V Suggested modifications for consistency (Mikhailovskii and Tsypin, Hazeltine and Meiss, Simakov and Catto, Ramos) involve adding term proportional to the ion heat rate of strain: "^q = 2 [ 5# b $W q % ( I + 3bb) + transpose] W q = &q i + &q i T ' 2 3 I& % q i Implicit numerical treatment difficult What is the effect of this term on dispersion and stability? Does it introduce new normal modes? Does it alter stability properties?
24 Ion Heat Stress Has Little Effect on Important Dynamics [ ( ) + transpose] " ^ q = 2 5# b$ W q % I + 3bb W q = &q i + &q i T ' 2 3 I&%q i q = '( & T '( ) & ) T '(^ b$ & ) T " 2 2 =C s k 2 % 1+ f (# ) $i k &' " 0 #V #t = $%p$ %& ' ^ q #p #t = $(p 0 %& V ( ) 2 ( )* f (0) = 0 f (+ / 2) = 1 Dispersive effect on compressional waves, but Negligible effect on g-mode stability Simplification: ignore these terms (for now!)
25 Careful Computational Approach is Required Spatial approximation Must capture anisotropy and global geometry Flux aligned grids High order finite elements Temporal approximation Must compute for long times Require implicit methods Semi-implicit methods have proven useful
26 Solenoidal Constraint Faraday: "B "t = #$ % E & " "t $ ' B = 0 Depends on " # "$ = 0 Different discrete approximations Modified wave system Projection Diffusion "U "t + # $ F = R# $ B B " = B + #$ # 2 $ = %# & B "B "t = #$ % E + &$$ ' B, " "t $ ' B = $ '&$$ ' B Grid properties " a # " b $ % 0 E.g., staggered grid, dual mesh
27 Galerkin Methods Finite differences and finite volumes minimize error locally Based on Taylor series expansion Galerkin methods minimize weighted error Based on expansion in basis functions Solve weak form of problem "u(x,t) "t % = Lu(x,t) # v ' "u &"t $ Lu ( + * dx = 0 ) Minimize error by expansion in basis functions and determining coefficients
28 Galerkin Discrete Approximation M ij du j dt = L ij u j M ij = $ dv" i # j Mass Matrix L ij = Solution generally requires inverting the mass matrix, even for explicit methods Different basis functions give different methods Usually: β i = α i α i =exp(ikx) => Fourier spectral methods α i =localized polynomial => finite element methods $ dv" i L# j Response Matrix
29 Finite Elements Project onto basis of locally defined polynomials of degree p Polynomials of degree p can converge as fast as h p+1 Integrate by parts: dv j ' " i " j dv = #' " i $% &(" j V j ) dv = ' $" i % &(" j V j )dv dt # ' " jn % &(" j V j )ds S Evaluate & only Boundary condition Simplifies implementation of complex closure relations Natural implementation of boundary conditions Automatically preserves self-adjointness Works well with arbitrary grid shapes
30 Three Examples of Favorable Properties of High Order Elements Z ln[ div(b) / b ] h 1 h 2 h 3 bilinear biquadratic bicubic 7 sim. data & fit 6 analytic W c w (cm) R ln(h) ! /! perp Grid used for ELM studies Non-uniform meshes retain high-order convergence rate Magnetic divergence constraint Scalings show expected convergence rates Critical island width for temperature flattening Dealing with extreme anisotropy Agreement on scaling
31 Multiple Time Scales (Parasitic Waves) MHD contains widely separated time scales (eigenvalues) "u "t = #u { = Fu { + Su { Full MHD operator Fast time scales: Alfvén waves, soundwaves, etc (parasitic waves) Slow time scales: Resistive instabilities, island evolution, (interesting physics) Parasitic waves are properties of the physics problem but are not the dynamics of interest Treat only fast part of operator implicitly to avoid time step restriction u n+1 " u n #t = Fu n+1 + Su n Precise decomposition of Ω for complex nonlinear system is often difficult or impractical to achieve algebraicly
32 Dealing with Parasitic Waves Original idea from André Robert (1971) Gravity waves in climate modeling F and Ω are often known, but an expression for S is difficult to achieve Ω: full MHD operator F: linearized MHD operator Use operator splitting: u n+1 " u n #t = Fu n+1 + ($ " F)u n = $u n % + #tf un+1 " u n ( ' * & #t ) Expression for S not needed " = F + S # S = " $ F
33 Semi-Implicit Method u n+1 " u n #t = Fu n+1 + ($ " F)u n = $u n % + #tf un+1 " u n ( ' * & #t ) Recognize that the operator F is completely arbitrary!! (I " #tg { )u n+1 = (I 1 " 4 #t$)u 243 n " #tgu 12 3 n S. I. operator Explicit S. I. operator G can be chosen for accuracy and ease of inversion G should be easier to invert than F (or Ω, e.g., toroidal coupling) G should approximate F for modes of interest Some choices are better than others! Has proven to be very useful for resistive and extended MHD Used for spheromak, RFP, tokamak, and solar corona modeling
34 SI Operator for MHD "B "t = # $ ( V $ B %&J ) "' = %# ( 'V "t "p = %# ( pv % () %1) p# ( V "t 2 ' "V "t = %'V ( #V % #p + J $ B + *+t 4, 2 # $ # $. "V "t $ B / 4 01 $ B Alfvén waves 5 + #)P 0 # ( "V "t 37 Sound waves 6 Ideal MHD operator (Lerbinger and Luciani) Anisotropic, self-adjoint Avoids implicit toroidal coupling (great simplification) Accurate linear results for CFL ~ (=> Condition number ~ 10 10!!) kv"t <1
35 Semi-Implicit Leap Frog t Variables staggered at different time levels SI operator on velocity "V =V j #V j#1 "p = p j+1/2 # p j#1/2 Δt t j+1/2 "V "t = #$p j#1/ 2 & +%"ts ( "V ' "t ) + * V p Δt t j t j-1/2 t j-1 V j = V + "V "p "t = #$P 0 %&V j p j+1/2 = p j#1/2 + "p
36 Extended MHD Time Advance Implicit leap-frog (also used in MHD) Maintains numerical stability without introducing numerical dissipation MHD advance unchanged (semi-implicit self-adjoint operators) Need to invert non-self-adjoint operators at each step for dispersive modes Requires high performance parallel linear algebra software
37 Implicit Leap Frog for Extended MHD % ( m i n j+1/ 2 ' "V + 1 "t 2 V j # $"V + 1 * ' j "V # $V * + "tl j+1/ 2 ("V) ' $ #, i ("V) = * & Implicit advection ) SI MHD Includes ALL stresses "n "t V j+1 # $"n = %$ # ( V j+1 # n j+1/ 2 ) 3n &"T # ( 2 ' "t ( ) J j+1/ 2 - B j+1/ 2 + m i n j+1/ 2 V j # $V j + $p j+1/ 2 + $ #, i V j V # j+1 ) $ %"T # % $ q # ("T # ) = * Anisotropic thermal conduction, 3n 2 V # j+1 $ %T # j+1/ 2, nt# j+1/ 2 % $ V# j+1, % $ q# T # j+1/ 2 ( ) ( ) + Q # j+1/ 2 "B "t V j+1 # $"B $ % 1 ne J j+1/ 2 % "B + "J % B j+1/ $ % &"J = Implicit HALL term Implicit resistive term ( ' $ % 1 ( ne J j+1/ 2 % B j+1/ 2 ' $p e ) ' V j+1 % B j+1/ 2 + &J j+1/ 2 + * - ), Momentum Continuity Energy Maxwell/Ohm
38 Nonlinear ELM Evolution Anisotropic thermal conduction ELM interaction with wall 70 kj lost in 60 µsec 2-fluid and gyro-viscosity have little effect on linear properties
39 Two-fluid Reconnection GEM Problem 2-D slab η = Good agreement with many other calculations Computed with same code used for tokamaks, spheromaks, RFPs
40 Heuristic Closure Captures Essential Neoclassical Physics Neo-classical theory gives flux surface average Local form for stress tensor forces: (Gianakon et.al., Phys. Plasmas 9, 536 (2002) "# $ % = & % µ % B 2 V % # e ' ( ) 2 e ' B % # e ' Valid for both ion and electrons Energy conserving and entropy producing Gives: bootstrap current neoclassical resistivity polarization current enhancement
41 Beyond Extended MHD: Parallel Kinetic closures Parallel closures for q and Π derived using Chapman-Enskog-like approach. Non-local; requires integration along perturbed field lines. General closures map continuously from collisional to nearly collisionless regime. General q closure predicts collisional response for heat flow inside magnetic island. As plasma becomes moderately collisional (T > 50 ev), general closure predicts correct flux limited response. Incorporated into global extended MHD algorithms. Thermal diffusivity as function of T showing T 5/2 response of Braginskii and general closure.
42 Beyond Extended MHD: Kinetic Minority ions species affects bulk evolution: Minority Species n h << n 0, " h ~ " 0 Mn dv dt = J " B # $ { % & Bulk Plasma # $ % & 12 3 h Hot Minority Ion Species "# h = % M (v $ V h )(v $ V h )"f (x,v)d 3 v δf determined by kinetic particle simulation in evolving fields Demonstrated transition from internal kink to fishbone Benchmark of three codes
43 Constraints on Modeling Balance of algorithm performance and problem requirements with available cycles N " Q 12 #t 3 Algorithms Algorithms: N - # of meshpoints for each dimension α - # of dimensions transport 3 (spatial) fluid 5-6 kinetic (spatial + velocity) Q - code-algorithm requirements (Tflop / meshpoint / timestep) Δt - time step (seconds) = 3 $10 7 %P CT { Constraints Constraints: P - peak hardware performance (Tflop/sec) ε - hardware efficiency εp - delivered sustained performance T - problem time duration (seconds) C - # of cases / year 1 case / week ==> C ~ 50
44 The Future???? N " Q = 3 $10 7 %P 12 #t 3 { Algorithms CT Constraints Assumptions: Performance is delivered Implicit algorithm Q ind. of Δt (!!) Requirements: At least 3-D physics required Required problem time: 1 msec - 1 sec Conclusions: 3-D (i.e., fluid) calculations for times of ~ 10 msec within reach Longer times require next generation computers (or better algorithms) Higher dimensional (kinetic) long time calculations unrealistic Integrated kinetic effects must come through low dimensionality fluid closures
45 Summary Fluid models are an approximation to the plasma kinetic equation, but are required for modeling low frequency response of hot, magnetized plasmas with global geometry Direct kinetic calculations are impractical Primitive equations and implicit methods have proven successful in modeling a variety of plasmas Implicit methods are required for handling the dispersive terms of MHD. An understanding of the dispersive characteristics of discretized equations needed. Kinetic effects must be captured through fluid closures Best form of fluid equations still unknown Often problem dependent Next step is direct coupling of kinetic/fluid/transport models
Computational modeling of fully ionized magnetized plasmas using the fluid approximation a
PHYSICS OF PLASMAS 13, 058103 2006 Computational modeling of fully ionized magnetized plasmas using the fluid approximation a D. D. Schnack b Center for Energy Space Science, Science Applications International
More informationComparison of Kinetic and Extended MHD Models for the Ion Temperature Gradient Instability in Slab Geometry
Comparison of Kinetic and Extended MHD Models for the Ion Temperature Gradient Instability in Slab Geometry D. D. Schnack University of Wisconsin Madison Jianhua Cheng, S. E. Parker University of Colorado
More informationModeling of ELM Dynamics for ITER
Modeling of ELM Dynamics for ITER A.Y. PANKIN 1, G. BATEMAN 1, D.P. BRENNAN 2, A.H. KRITZ 1, S. KRUGER 3, P.B. SNYDER 4 and the NIMROD team 1 Lehigh University, 16 Memorial Drive East, Bethlehem, PA 18015
More informationGA A25229 FINAL REPORT ON THE OFES ELM MILESTONE FOR FY05
GA A25229 ELM MILESTONE FOR FY05 by D.P. BRENNAN, E.D. HELD, S.E. KRUGER, A.Y. PANKIN, D.D. SCHNACK, AND C.R. SOVINEC APRIL 2006 DISCLAIMER This report was prepared as an account of work sponsored by an
More informationA 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 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 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 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 informationMAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT
MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT ABSTRACT A. G. Tarditi and J. V. Shebalin Advanced Space Propulsion Laboratory NASA Johnson Space Center Houston, TX
More informationGyroviscosity in the NIMROD code. A.Y. Pankin D. D. Schnack C. Sovinec S. E. Kruger
Gyroviscosity in the NIMROD code A.Y. Pankin D. D. Schnack C. Sovinec S. E. Kruger Outline Introduction Two possible implementations Testing GV in the NIMROD code Convergence Summary ε = ω / Ω i ξ = V
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
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,
More informationComputations with Discontinuous Basis Functions
Computations with Discontinuous Basis Functions Carl Sovinec University of Wisconsin-Madison NIMROD Team Meeting November 12, 2011 Salt Lake City, Utah Motivation The objective of this work is to make
More informationIntermediate Nonlinear Development of a Line-tied g-mode
Intermediate Nonlinear Development of a Line-tied g-mode Ping Zhu University of Wisconsin-Madison In collaboration with C. C. Hegna and C. R. Sovinec (UW-Madison) A. Bhattacharjee and K. Germaschewski
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 informationFLUID MODELS OF MAGNETIZED PLASMAS or BEYOND MHD or PLASMA PHYSICS IN 100 MINUTES!
FLUID MODELS OF MAGNETIZED PLASMAS or BEYOND MHD or PLASMA PHYSICS IN 100 MINUTES! Dalton D. Schnack April 18, 2006 1 Introduction You are now all experts in MHD. As you know, ideal MHD describes the dynamics
More informationGyrokinetic Simulations of Tearing Instability
Gyrokinetic Simulations of Tearing Instability July 6, 2009 R. NUMATA A,, W. Dorland A, N. F. Loureiro B, B. N. Rogers C, A. A. Schekochihin D, T. Tatsuno A rnumata@umd.edu A) Center for Multiscale Plasma
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 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 informationHybrid Kinetic-MHD simulations with NIMROD
simulations with NIMROD Charlson C. Kim and the NIMROD Team Plasma Science and Innovation Center University of Washington CSCADS Workshop July 18, 2011 NIMROD C.R. Sovinec, JCP, 195, 2004 parallel, 3-D,
More informationSteady-state Force Balance in the DEBS Code
Steady-state Force Balance in the DEBS Code D. D. Schnack Department of Physics University of Wisconsin - Madison Madison, WI 53706 November 10, 2013 The DEBS code advances the induction equation, the
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 informationGyrokinetic simulations of magnetic fusion plasmas
Gyrokinetic simulations of magnetic fusion plasmas Tutorial 2 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr
More information(a) (b) (c) (d) (e) (f) r (minor radius) time. time. Soft X-ray. T_e contours (ECE) r (minor radius) time time
Studies of Spherical Tori, Stellarators and Anisotropic Pressure with M3D 1 L.E. Sugiyama 1), W. Park 2), H.R. Strauss 3), S.R. Hudson 2), D. Stutman 4), X-Z. Tang 2) 1) Massachusetts Institute of Technology,
More informationExponential Growth of Nonlinear Ballooning Instability. Abstract
Exponential Growth of Nonlinear Ballooning Instability P. Zhu, C. C. Hegna, and C. R. Sovinec Center for Plasma Theory and Computation University of Wisconsin-Madison Madison, WI 53706, USA Abstract Recent
More informationNIMROD Boundary Conditions
NIMROD Boundary Conditions Carl Sovinec University of Wisconsin-Madison and the Two-Fluid and Transport Group Plasma Science and Innovation Center Annual Meeting Seattle Washington, August 13-14, 2007
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 informationSemi-implicit Treatment of the Hall Effect in NIMROD Simulations
Semi-implicit Treatment of the Hall Effect in NIMROD Simulations H. Tian and C. R. Soinec Department of Engineering Physics, Uniersity of Wisconsin-Madison Madison, WI 5376 Presented at the 45th Annual
More informationTheory for Neoclassical Toroidal Plasma Viscosity in a Toroidally Symmetric Torus. K. C. Shaing
Theory for Neoclassical Toroidal Plasma Viscosity in a Toroidally Symmetric Torus K. C. Shaing Plasma and Space Science Center, and ISAPS, National Cheng Kung University, Tainan, Taiwan 70101, Republic
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 informationExponential Growth and Filamentary Structure of Nonlinear Ballooning Instability 1
Exponential Growth and Filamentary Structure of Nonlinear Ballooning Instability 1 Ping Zhu in collaboration with C. C. Hegna and C. R. Sovinec University of Wisconsin-Madison Sherwood Conference Denver,
More informationSawteeth in Tokamaks and their relation to other Two-Fluid Reconnection Phenomena
Sawteeth in Tokamaks and their relation to other Two-Fluid Reconnection Phenomena S. C. Jardin 1, N. Ferraro 2, J. Chen 1, et al 1 Princeton Plasma Physics Laboratory 2 General Atomics Supported by the
More informationPLASMA-NEUTRAL MODELING IN NIMROD. Uri Shumlak*, Sina Taheri*, Jacob King** *University of Washington **Tech-X Corporation April 2016
PLASMA-NEUTRAL MODELING IN NIMROD Uri Shumlak*, Sina Taheri*, Jacob King** *University of Washington **Tech-X Corporation April 2016 Plasma-Neutral Model Physical Model is derived by E. Meier and U. Shumlak*
More informationTokamak Fusion Basics and the MHD Equations
MHD Simulations for Fusion Applications Lecture 1 Tokamak Fusion Basics and the MHD Equations Stephen C. Jardin Princeton Plasma Physics Laboratory CEMRACS 1 Marseille, France July 19, 21 1 Fusion Powers
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 informationNIMROD Project Overview
NIMROD Project Overview Christopher Carey - Univ. Wisconsin NIMROD Team www.nimrodteam.org CScADS Workshop July 23, 2007 Project Overview NIMROD models the macroscopic dynamics of magnetized plasmas by
More informationModels for Global Plasma Dynamics
Models for Global Plasma Dynamics F.L. Waelbroeck Institute for Fusion Studies, The University of Texas at Austin International ITER Summer School June 2010 Outline 1 Models for Long-Wavelength Plasma
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More informationHeat Transport in a Stochastic Magnetic Field. John Sarff Physics Dept, UW-Madison
Heat Transport in a Stochastic Magnetic Field John Sarff Physics Dept, UW-Madison CMPD & CMSO Winter School UCLA Jan 5-10, 2009 Magnetic perturbations can destroy the nested-surface topology desired for
More information36. TURBULENCE. Patriotism is the last refuge of a scoundrel. - Samuel Johnson
36. TURBULENCE Patriotism is the last refuge of a scoundrel. - Samuel Johnson Suppose you set up an experiment in which you can control all the mean parameters. An example might be steady flow through
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 informationMHD Linear Stability Analysis Using a Full Wave Code
US-Japan JIFT Workshop on Progress of Extended MHD Models NIFS, Toki,Japan 2007/03/27 MHD Linear Stability Analysis Using a Full Wave Code T. Akutsu and A. Fukuyama Department of Nuclear Engineering, Kyoto
More 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 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 informationFundamentals of Magnetic Island Theory in Tokamaks
Fundamentals of Magnetic Island Theory in Tokamaks Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Austin, TX, USA Talk available at http://farside.ph.utexas.edu/talks/talks.html
More informationGyrokinetic Theory and Dynamics of the Tokamak Edge
ASDEX Upgrade Gyrokinetic Theory and Dynamics of the Tokamak Edge B. Scott Max Planck Institut für Plasmaphysik D-85748 Garching, Germany PET-15, Sep 2015 these slides: basic processes in the dynamics
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 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 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 informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel
More informationPredictive Simulation of Global Instabilities in Tokamaks
Predictive Simulation of Global Instabilities in Tokamaks Stephen C. Jardin Princeton Plasma Physics Laboratory Fourth ITER International Summer School IISS 2010 Institute of Fusion Studies University
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 informationProgress in the Plasma Science and Innovation Center
1 THP/2-02 Progress in the Plasma Science and Innovation Center U. Shumlak 1, C. Akcay 1, A. H. Glasser 1, C. J. Hansen 1, E. D. Held 2, T. R. Jarboe 1, J.-Y. Ji 2, C. Kim 1, W. Lowrie 1, V. S. Lukin 4,
More informationElectron parallel transport for arbitrary collisionality
Utah State University DigitalCommons@USU Publications Atmospheric Imaging Laboratory 11-2017 Electron parallel transport for arbitrary collisionality Jeong-Young Ji Utah State University, j.ji@usu.edu
More informationMicrotearing Simulations in the Madison Symmetric Torus
Microtearing Simulations in the Madison Symmetric Torus D. Carmody, P.W. Terry, M.J. Pueschel - University of Wisconsin - Madison dcarmody@wisc.edu APS DPP 22 Overview PPCD discharges in MST have lower
More informationConservation Laws in Ideal MHD
Conservation Laws in Ideal MHD Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 3, 2016 These lecture notes are largely based on Plasma Physics for Astrophysics
More information14. Energy transport.
Phys780: Plasma Physics Lecture 14. Energy transport. 1 14. Energy transport. Chapman-Enskog theory. ([8], p.51-75) We derive macroscopic properties of plasma by calculating moments of the kinetic equation
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 informationSimulations of Sawteeth in CTH. Nicholas Roberds August 15, 2015
Simulations of Sawteeth in CTH Nicholas Roberds August 15, 2015 Outline Problem Description Simulations of a small tokamak Simulations of CTH 2 Sawtoothing Sawtoothing is a phenomenon that is seen in all
More informationHybrid Kinetic-MHD simulations with NIMROD
in NIMROD simulations with NIMROD Charlson C. Kim 1 Dylan P. Brennan 2 Yasushi Todo 3 and the NIMROD Team 1 University of Washington, Seattle 2 University of Tulsa 3 NIFS, Toki-Japan December 2&3, 2011
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 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 informationINTRODUCTION TO GYROKINETIC AND FLUID SIMULATIONS OF PLASMA TURBULENCE AND OPPORTUNITES FOR ADVANCED FUSION SIMULATIONS
INTRODUCTION TO GYROKINETIC AND FLUID SIMULATIONS OF PLASMA TURBULENCE AND OPPORTUNITES FOR ADVANCED FUSION SIMULATIONS G.W. Hammett, Princeton Plasma Physics Lab w3.pppl.gov/ hammett Fusion Simulation
More informationMeasuring from electron temperature fluctuations in the Tokamak Fusion Test Reactor
PHYSICS OF PLASMAS VOLUME 5, NUMBER FEBRUARY 1998 Measuring from electron temperature fluctuations in the Tokamak Fusion Test Reactor C. Ren, a) J. D. Callen, T. A. Gianakon, and C. C. Hegna University
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 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 informationOverview of Gyrokinetic Theory & Properties of ITG/TEM Instabilities
Overview of Gyrokinetic Theory & Properties of ITG/TEM Instabilities G. W. Hammett Princeton Plasma Physics Lab (PPPL) http://w3.pppl.gov/~hammett AST559: Plasma & Fluid Turbulence Dec. 5, 2011 (based
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 informationStability of a plasma confined in a dipole field
PHYSICS OF PLASMAS VOLUME 5, NUMBER 10 OCTOBER 1998 Stability of a plasma confined in a dipole field Plasma Fusion Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 Received
More information3D hybrid-kinetic turbulence and phase-space cascades
3D hybrid-kinetic turbulence and phase-space cascades ( in a β = 1 plasma ) Silvio Sergio Cerri Department of Astrophysical Sciences, Princeton University, USA 11th Plasma Kinetics Working Meeting WPI
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 informationPlasma instabilities. Dr Ben Dudson, University of York 1 / 37
Plasma instabilities Dr Ben Dudson, University of York 1 / 37 Previously... Plasma configurations and equilibrium Linear machines, and Stellarators Ideal MHD and the Grad-Shafranov equation Collisional
More informationImposed Dynamo Current Drive
Imposed Dynamo Current Drive by T. R. Jarboe, C. Akcay, C. J. Hansen, A. C. Hossack, G. J. Marklin, K. Morgan, B. A. Nelson, D. A. Sutherland, B. S. Victor University of Washington, Seattle, WA 98195,
More informationNIMROD simulations of dynamo experiments in cylindrical and spherical geometries. Dalton Schnack, Ivan Khalzov, Fatima Ebrahimi, Cary Forest,
NIMROD simulations of dynamo experiments in cylindrical and spherical geometries Dalton Schnack, Ivan Khalzov, Fatima Ebrahimi, Cary Forest, 1 Introduction Two experiments, Madison Plasma Couette Experiment
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 information35. RESISTIVE INSTABILITIES: CLOSING REMARKS
35. RESISTIVE INSTABILITIES: CLOSING REMARKS In Section 34 we presented a heuristic discussion of the tearing mode. The tearing mode is centered about x = 0, where F x that B 0 ( ) = k! B = 0. Note that
More informationEnergetic-Ion-Driven MHD Instab. & Transport: Simulation Methods, V&V and Predictions
Energetic-Ion-Driven MHD Instab. & Transport: Simulation Methods, V&V and Predictions 7th APTWG Intl. Conference 5-8 June 2017 Nagoya Univ., Nagoya, Japan Andreas Bierwage, Yasushi Todo 14.1MeV 10 kev
More informationAMSC 663 Project Proposal: Upgrade to the GSP Gyrokinetic Code
AMSC 663 Project Proposal: Upgrade to the GSP Gyrokinetic Code George Wilkie (gwilkie@umd.edu) Supervisor: William Dorland (bdorland@umd.edu) October 11, 2011 Abstract Simulations of turbulent plasma in
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 informationGA A27235 EULERIAN SIMULATIONS OF NEOCLASSICAL FLOWS AND TRANSPORT IN THE TOKAMAK PLASMA EDGE AND OUTER CORE
GA A27235 EULERIAN SIMULATIONS OF NEOCLASSICAL FLOWS AND TRANSPORT IN THE TOKAMAK PLASMA EDGE AND OUTER CORE by E.A. BELLI, J.A. BOEDO, J. CANDY, R.H. COHEN, P. COLELLA, M.A. DORF, M.R. DORR, J.A. HITTINGER,
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 informationExtended MHD simulation of Rayleigh-Taylor/Kelvin-Helmholtz instability
Extended MHD simulation of Rayleigh-Taylor/Kelvin-Helmholtz instability R. Goto (a), T. Hatori (a), H. Miura (a, b), A. Ito (a, b) and M. Sato (b) (a) The Graduate University for Advanced Studies (Sokendai)
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 informationIssues in Neoclassical Tearing Mode Theory
Issues in Neoclassical Tearing Mode Theory Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Austin, TX Tearing Mode Stability in Tokamaks According to standard (single-fluid)
More 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 informationCollisions and transport phenomena
Collisions and transport phenomena Collisions in partly and fully ionized plasmas Typical collision parameters Conductivity and transport coefficients Conductivity tensor Formation of the ionosphere and
More informationMagnetic Self-Organization in the RFP
Magnetic Self-Organization in the RFP Prof. John Sarff University of Wisconsin-Madison Joint ICTP-IAEA College on Plasma Physics ICTP, Trieste, Italy Nov 7-18, 2016 The RFP plasma exhibits a fascinating
More informationGuangye Chen, Luis Chacón,
JIFT workshop! Oct. 31, 2014 New Orleans, LA.! Guangye Chen, Luis Chacón, CoCoMANs team Los Alamos National Laboratory, Los Alamos, NM 87545, USA gchen@lanl.gov 1 Los Alamos National Laboratory Motivation
More informationInnovative Concepts Workshop Austin, Texas February 13-15, 2006
Don Spong Oak Ridge National Laboratory Acknowledgements: Jeff Harris, Hideo Sugama, Shin Nishimura, Andrew Ware, Steve Hirshman, Wayne Houlberg, Jim Lyon Innovative Concepts Workshop Austin, Texas February
More informationHybrid Kinetic-MHD simulations with NIMROD
simulations with NIMROD 1 Yasushi Todo 2, Dylan P. Brennan 3, Kwang-Il You 4, Jae-Chun Seol 4 and the NIMROD Team 1 University of Washington, Seattle 2 NIFS, Toki-Japan 3 University of Tulsa 4 NFRI, Daejeon-Korea
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 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 informationDerivation of dynamo current drive in a closed current volume and stable current sustainment in the HIT SI experiment
Derivation of dynamo current drive and stable current sustainment in the HIT SI experiment 1 Derivation of dynamo current drive in a closed current volume and stable current sustainment in the HIT SI experiment
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 informationTHE DIII D PROGRAM THREE-YEAR PLAN
THE PROGRAM THREE-YEAR PLAN by T.S. Taylor Presented to Program Advisory Committee Meeting January 2 21, 2 3 /TST/wj PURPOSE OF TALK Show that the program plan is appropriate to meet the goals and is well-aligned
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 informationNIMROD Team Meeting Minutes August 12-14, 2014, Logan, Utah The team meeting began with an informal discussion of features that would be beneficial
NIMROD Team Meeting Minutes August 12-14, 2014, Logan, Utah The team meeting began with an informal discussion of features that would be beneficial for applications and for development. Participants were
More informationNIMROD FROM THE CUSTOMER S PERSPECTIVE MING CHU. General Atomics. Nimrod Project Review Meeting July 21 22, 1997
NIMROD FROM THE CUSTOMER S PERSPECTIVE MING CHU General Atomics Nimrod Project Review Meeting July 21 22, 1997 Work supported by the U.S. Department of Energy under Grant DE-FG03-95ER54309 and Contract
More informationAccurate representation of velocity space using truncated Hermite expansions.
Accurate representation of velocity space using truncated Hermite expansions. Joseph Parker Oxford Centre for Collaborative Applied Mathematics Mathematical Institute, University of Oxford Wolfgang Pauli
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 information1 Energy dissipation in astrophysical plasmas
1 1 Energy dissipation in astrophysical plasmas The following presentation should give a summary of possible mechanisms, that can give rise to temperatures in astrophysical plasmas. It will be classified
More information