Accurate representation of velocity space using truncated Hermite expansions.
|
|
- Cecily Wilkins
- 5 years ago
- Views:
Transcription
1 Accurate representation of velocity space using truncated Hermite expansions. Joseph Parker Oxford Centre for Collaborative Applied Mathematics Mathematical Institute, University of Oxford Wolfgang Pauli Institute, Vienna 3 rd March 01 Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 1 / 34
2 Model for ITG instabilities (from Pueschel et al., 010) Test problem for GENE, simplified gyrokinetic equations: 1D slab, (z, v ) adiabatic electrons B = Bz homogeneous and static drift-kinetic ions electrostatic: E = Φ quasineutral uniform density and temperature gradients in background perturbation about Maxwellian F 0 (v ) = π 1/ exp( v ) nonlinearity dropped May also derive directly from Vlasov equation and quasineutrality. Electrostatic limit of Belli & Hammett (005) model Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 / 34
3 Model for ITG instabilities (from Pueschel et al., 010) Linearized kinetic equation F 1 t + α iv F 1 z + [ω n + ω Ti ( v 1 )] ik y ΦF 0 + τ e α i v Φ z F 0 = 0 Quasineutrality Φ = F 1 dv. Parameters ω n = 1, normalized density gradient ω Ti = 10, normalized ion temperature gradient k y = 0.3, perpendicular wavenumber τ e = 1, species temperature ratio (T i = T e ) α i = 0.34 nondimensional constant: velocity and length scales Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 3 / 34
4 Wave-like solutions Leads to eigenvalue problem, F 1 = f(v ) exp ( i ( k z ωt )) L 0 f = ωf ( where L 0 f = α i kv f + [ω n + ω Ti v 1 )] k y ΦF 0 + τ e α i kv ΦF 0 L real = eigenvalues in complex conjugate pairs = no damping So actually want to solve, lim L ɛf = ωf, for L ɛ = L 0 + iɛc ɛ 0 where C is collisions, e.g. Lénard Bernstein or BGK. Breaks symmetry, so decay rate may tend to nonzero limit as ɛ 0. Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 4 / 34
5 0.4 Dispersion relation 0. growth rate, γ parallel wavenumber, k Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 5 / 34
6 Discretization Keep exp(ik z) behaviour, represent v on a uniform grid, Obtain linear ODE system f t = imf, Φ = N vf(v i ) i=1 M real N N matrix Matrix eigenvalue problem (quicker than running IVPs): iωf = f t = imf Growth rate is γ = max I(ω), for ω spectrum(m). Matrix M depends on: All parameters, particularly parallel wavenumber. Discretization and resolution of grid in v. Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 6 / 34
7 Growth rate for simple discretized system N = 10 N = 100 dispersion relation growth rate, γ parallel wavenumber, k Capture growth, but convergence slow with N. No decay. ɛ 0 at fixed N resolution. Need to restore collisions. Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 7 / 34
8 Distribution function in velocity space k =.75 k = 3.5 k = 3.75 f(v) v Fine scales develop as I(ω) 0. Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 8 / 34
9 Hermite expansion in velocity space Use expansion in asymmetric Hermite functions, f(v) = a m φ m (v), φ m (v) = H m(v) m m!, φ m(v) = F 0 (v)φ m (v) m=0 Bi-orthogonal polynomials with Maxwellian as weight function φ m φ n dv = δ mn Recurrence relation m + 1 vφ m = φ m+1 + m φ m 1 = particle streaming becomes mode coupling Represent velocity space scales, ( H m (v) cos v m mπ ), as m large m = fine scales Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 9 / 34
10 Relative Entropy and Free Energy Relative entropy is R[F F 0 ] = ( ) F F log F + F 0 dv F 0 F log F dv + U + function(density) Expansion F = F 0 + ɛf 1 gives (at leading order) R[F F 0 ] = ɛ F 1 F 0 dv = ɛ a m + O(ɛ 3 ) m=0 Define free energy of each Hermite mode, E m, E = E m, with E m = 1 a m m=0 cf. energy spectra in Fourier space for Navier Stokes turbulence. Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March / 34
11 Moment system Replace F 0 φ 0 and v with Hermite functions [ F 1 t + iα ikvf 1 + ik y Φ ω n φ 0 + ω ] T i φ + i kτ eα i Φ φ 1 = 0 Put F 1 = a m φ m (implicit sum over repeated m) [ a m t φ m + iα i kva m φ m + ik y Φ ω n φ 0 + ω ] T i φ + i kτ eα i Φ φ 1 = 0 Use recurrence relation on particle streaming ( ) a m m + 1 m t φ m + iα i ka m φ m+1 + φ m = 0 Gives infinite set of coupled algebraic equations for {a n }. Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March / 34
12 System of equations for expansion coefficients ( ) m + 1 m ωa m = a m+1 + a m 1 [ ] + k y ω n a 0 δ m0 + ω T i a 0 τ e a δ m + 0 δ m1 α i k }{{}}{{} driving Boltzmann response We have also used φ 0 1, so that, Φ = F 1 φ 0 dv = m=0 a m φ m φ 0 dv = a 0 Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 1 / 34
13 Matrix equation ω nk y 1 α i k τ e k yω Ti αi k N N a = ωa Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March / 34
14 Simple Truncation ω nk y 1 α i k τ e k yω Ti αi k a = ωa Assume a N+1 = 0 at point of truncation. Last row of matrix: N + 1 N a N+1 + a N 1 = ωa N Equivalent to discretization on Gauss Hermite points {v j } Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March / 34
15 10 0 Hermite spectra of eigenfunctions free energy, am k =, γ > 0 k = 6, γ = Hermite mode, m + 1 Growing mode: a m decay as m increases = a N+1 0 is okay (Not) decaying mode not resolved Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March / 34
16 Collisions F 1 t + Free energy should cascade to large m (fine scales in v). Free energy has nowhere to go when a N+1 = 0. Restore collisions: [ω n + ω Ti ( v 1 Desirable properties for C[F 1 ] )] ik y ΦF 0 + iα i kv F 1 + iτ e α i kv ΦF 0 = C[F 1 ] conserves mass, momentum and energy C[F 1 ] dv = vc[f 1 ] dv = v C[F 1 ] dv = 0 satisfies a linearized H theorem: dr dt = ɛ F 1 C[F 1 ] F 0 dv 0 represents small-angle collisions (contains v-derivatives) Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March / 34
17 Lénard Bernstein (1958) collisions Simple member of linearized Landau/Fokker Planck class: C [F 1 ] = ν [ vf ] F 1 v v collision frequency ν Hermite functions are eigenfunctions: C [a m φ m ] = νma m φ m = easy to implement in Hermite space Conserves mass, but not momentum or energy Satisfies a linearized R theorem: dr dt = ɛ ν m m a m 0 Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March / 34
18 Lénard Bernstein in Hermite space ω nk y 1 α i k τ e + 1 k yω Ti iν 1 αi k 1 iν 3 3 3iν N N inν a = ωa matrix now complex roots not in complex-conjugate pairs = can find negative growth rates (we can manually conserve momentum and energy) Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March / 34
19 Lénard Bernstein in Hermite space ω nk y 1 α i k τ e + 1 iν 1 k yω Ti αi 1 iν k 3 3 3iν N N inν a = ωa matrix now complex roots not in complex-conjugate pairs = can find negative growth rates (we can manually conserve momentum and energy) Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March / 34
20 Growth rate, Lenard Berstein collisions, N = growth rate, γ ν = 10 1 ν = 10 ν = 10 3 dispersion relation parallel wavenumber, k If collision frequency ν large enough to get k > 4 correct, then shape for k < 4 distorted. Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March / 34
21 Growth rate, Lenard Berstein collisions, ν = growth rate, γ N = 100 N = 00 N = 300 dispersion relation parallel wavenumber, k resolves dissipative scales expensive: 300 modes Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 0 / 34
22 Lenard Bernstein collisions, ν = 10 k = k = 6 free energy, am Hermite mode, m + 1 Appreciable damping along whole spectrum. Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 1 / 34
23 Hypercollisions Would prefer: low modes undamped high modes strongly damped... whatever the truncation point N Iterate Lénard Bernstein collisions: L[F 1 ] = ( vf ) F 1, C[F 1 ] = ν( N) n L n [F 1 ] v v (like hyperdiffusion: ( ) n in physical space) Hermite functions are still eigenfunctions: ( m ) n C [a m φ m ] = ν am φ m N ν sets the decay rate of the highest mode linearized R-theorem: dr ( m ) n dt = ɛ ν am 0 N m=0 n = 1 corresponds to Lénard Bernstein collisions. two parameters: n 6, ν 10 (robust to variation) Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 / 34
24 0. Growth rate with hypercollisions, ν = 10, n = 6 growth rate, γ simple truncation, N = 100 hypercollisions, N = 10 exact dispersion relation parallel wavenumber, k Captures decaying parts of spectrum. Excellent fit, fast convergence. Appropriate for nonlinear problems. Two parameters, n, ν, robust to variation. Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 3 / 34
25 10 0 Hermite spectra with hypercollisions, k = 10 1 free energy, am / N = 10 N = Hermite mode, m + 1 Low moments largely undamped Damping at high m, for any N. Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 4 / 34
26 10 0 Hermite spectra with hypercollisions, k = 6 free energy, am / N = 10 N = Hermite mode, m + 1 Low moments largely undamped Damping at high m, for any N. Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 5 / 34
27 Energy equations and theoretical spectra Equation for coefficients, valid for m 3, ( ) m + 1 m ωa m = a m+1 + a m 1 + driving + Boltzmann response Treat as a finite difference approximation in continuous m. Energy equation for E m = a m / (Zocco & Schekochihin, 011) E m + ( ) ( m ) n mem = ν Em. t m N For a mode with growth rate γ, ( E m = C exp γ ( m m γ m γ ) 1/ ( ) ) m n+1/, m c with the cutoffs, m γ = 1 8γ, m(n+1/) c = [ N n ] (n + 1/) ν Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 6 / 34
28 Energy equations and theoretical spectra Equation for coefficients, valid for m 3, ( ) m + 1 m ωa m = a m+1 + a m 1 + hypercollisions Treat as a finite difference approximation in continuous m. Energy equation for E m = a m / (Zocco & Schekochihin, 011) E m + ( ) ( m ) n mem = ν Em. t m N For a mode with growth rate γ, ( E m = C exp γ ( m m γ m γ ) 1/ ( ) ) m n+1/, m c with the cutoffs, m γ = 1 8γ, m(n+1/) c = [ N n ] (n + 1/) ν Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 6 / 34
29 Growing mode, k = collision dominated free energy, Em m 1/ behaviour theoretical spectrum numerical spectrum growth rate cutoff collisional cutoff growth rate dominated E m = Hermite mode, m + 1 ( C exp γ ( ) m 1/ ( ) ) m n+1/ m γ m γ m c Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 7 / 34
30 10 1 Decaying mode, k = 6 free energy, Em m 1/ growth rate dominated theoretical spectrum numerical spectrum growth rate cutoff collisional cutoff collision dominated E m = Hermite mode, m + 1 ( C exp γ ( ) m 1/ ( ) ) m n+1/ m γ m γ m c Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 8 / 34
31 How strong should collisions be? We can find the collision strength required for a given resolution. Write hypercollisions as, C[F 1 ] = ν( L) n [F 1 ] (i.e. remove N n, damping strength expressed just by ν.) Need to resolve collisional cutoff, ( ) (n + 1/) 1/(n+1/) N > m c = ν as ν 0 Need infinite resolution to resolve collisionless case. Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 9 / 34
32 free energy, am / Hermite spectra with different ν, k = 6, n =4, N = ν = ν = ν = Hermite mode, m + 1 Weaker collisions = finer scales in spectra Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March / 34
33 growth rate, γ Growth rate against collision strength, k = 6, n = N = 0 N = 50 N = 100 correct γ hypercollision strength, ν a range of ν give the correct growth rate range extends to smaller ν as N increases Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March / 34
34 Hypercollisions on other grids Easiest to implement in Hermite space: C[a m ] = ν(m/n) n a m C[a] = Da But may be used on any grid. Map function values f to Hermite space by a = Mf, collide and map back, C[f] = M 1 DMf Only need to calculate M 1 DM once. growth rate, γ Growth rate calculated on uniform grid, N = 16 simple truncation hypercollisions exact dispersion relation parallel wavenumber, k Example: hypercollisions implemented on a uniform grid. Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 3 / 34
35 Summary 1D model for ITG instability velocity space discretization = eigenvalue problem finite composition of normal modes = no Landau damping Hermite representation in decaying modes, have energy pile-up at small scales Damping with collisions Lénard Bernstein collisions finds damping, but requires O(100) terms hypercollisions excellent agreement with dispersion relation theoretical expression for eigenfunctions only 10 terms robust parameters easy to implement in Hermite space can use on any grid Vanishing collisions is different from collisionless. Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March / 34
36 References BELLI, E. A. & HAMMETT, G. W. (005) A numerical instability in an ADI algorithm for gyrokinetics. Computer Physics Communications 17 (), LANDAU, L. D. (1946) On the Vibrations of the Electronic Plasma. Journal of Physics-U.S.S.R. 10 (5). LÉNARD, A. & BERNSTEIN, I. B. (1958) Plasma oscillations with diffusion in velocity space. Physical Review 11 (5), PUESCHEL, M. J., DANNERT, T. & JENKO, F. (010) On the role of numerical dissipation in gyrokinetic Vlasov simulations of plasma microturbulence. Computer Physics Communications 181, VAN KAMPEN, N. G. (1955) On the Theory of Stationary Waves in Plasmas. Physica 1, ZOCCO, A. & SCHEKOCHIHIN, A. A. (011) Reduced fluid-kinetic equations for low-frequency dynamics, magnetic reconnection, and electron heating in low-beta plasmas. Physics of Plasmas 18 (10), Supported by Award No KUK-C from King Abdullah University of Science and Technology (KAUST). Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March / 34
37 Drift approximation in the Vlasov equation E + v B c Vlasov equation becomes = 0, v = v ẑ c Φ B, B = Bẑ, E = Φ. B f ( t + v ẑ c ) B Φ B f q Φ m z f = 0. v Perturb about stationary equilibrium, f(x, v, t) = F 0 (x, v) + ɛf 1 (x, v, t) Impose density and temperature gradients in x, [ ( F 0 x = 1 v + v )] ω n + ω Ti L vth 3 F 0 Assume a Fourier mode in y: y ik y. Integrate out perpendicular directions in v. ( F 1 [ω t + n + ω Ti v 1 )] F 1 ik y ΦF 0 + α i v z + τ Φ eα i v z F 0 = 0 Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01 1
38 Quasineutrality Poisson s equation Boltzmann electrons Φ = 4πq(n i n e ), ( ) qφ n e = n e exp Nondimensionalize (with ɛ = ρ s /L) ( ɛt e 4π n i q ρ s L Φ = 1 + ɛ ( ɛ Φ = 1 n ) e + ɛ n i ( Φ, potential for electrostatic perturbation Φ = T e ) F 1 dv n e exp (ɛφ). n i F 1 dv n ) e Φ n i F 1 dv. + O(ɛ ). Joseph Parker (University of Oxford) Velocity space using Hermite expansions 3 rd March 01
Kinetic damping in gyro-kinetic simulation and the role in multi-scale turbulence
2013 US-Japan JIFT workshop on New Aspects of Plasmas Kinetic Simulation NIFS, November 22-23, 2013 Kinetic damping in gyro-kinetic simulation and the role in multi-scale turbulence cf. Revisit for Landau
More informationGyrokinetic Large Eddy Simulations
Gyrokinetic Large Eddy Simulations A. Bañón Navarro 1, P. Morel 1, M. Albrecht-Marc 1, D. Carati 1, F. Merz 2, T. Görler 2, and F. Jenko 2 1 Laboratoire de Physique Statistique et des Plasmas Université
More informationSummer College on Plasma Physics. 30 July - 24 August, The particle-in-cell simulation method: Concept and limitations
1856-30 2007 Summer College on Plasma Physics 30 July - 24 August, 2007 The particle-in-cell M. E. Dieckmann Institut fuer Theoretische Physik IV, Ruhr-Universitaet, Bochum, Germany The particle-in-cell
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 informationMultiscale, multiphysics modeling of turbulent transport and heating in collisionless, magnetized plasmas
Multiscale, multiphysics modeling of turbulent transport and heating in collisionless, magnetized plasmas Michael Barnes Plasma Science & Fusion Center Massachusetts Institute of Technology Collaborators:
More informationCollisional damping of plasma waves on a pure electron plasma column
PHYSICS OF PLASMAS 14, 112110 2007 Collisional damping of plasma waves on a pure electron plasma column M. W. Anderson and T. M. O Neil Department of Physics, University of California at San Diego, La
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 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 informationFourier Hermite spectral representation for the Vlasov Poisson system in the weakly collisional limit
To appear in J. Plasma Phys. 1 Fourier Hermite spectral representation for the Vlasov Poisson system in the weakly collisional limit J O S E P H T. P A R K E R A N D P A U L J. D E L L A R OCIAM, Mathematical
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 informationIII.H Zeroth Order Hydrodynamics
III.H Zeroth Order Hydrodynamics As a first approximation, we shall assume that in local equilibrium, the density f 1 at each point in space can be represented as in eq.(iii.56), i.e. [ ] p m q, t)) f
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 informationTURBULENT TRANSPORT THEORY
ASDEX Upgrade Max-Planck-Institut für Plasmaphysik TURBULENT TRANSPORT THEORY C. Angioni GYRO, J. Candy and R.E. Waltz, GA The problem of Transport Transport is the physics subject which studies the physical
More informationLectures on basic plasma physics: Kinetic approach
Lectures on basic plasma physics: Kinetic approach Department of applied physics, Aalto University April 30, 2014 Motivation Layout 1 Motivation 2 Boltzmann equation (a nasty bastard) 3 Vlasov equation
More informationPHYSICS OF HOT DENSE PLASMAS
Chapter 6 PHYSICS OF HOT DENSE PLASMAS 10 26 10 24 Solar Center Electron density (e/cm 3 ) 10 22 10 20 10 18 10 16 10 14 10 12 High pressure arcs Chromosphere Discharge plasmas Solar interior Nd (nω) laserproduced
More informationFundamentals of wave kinetic theory
Fundamentals of wave kinetic theory Introduction to the subject Perturbation theory of electrostatic fluctuations Landau damping - mathematics Physics of Landau damping Unmagnetized plasma waves The plasma
More information= qe cos(kz ωt). ω sin(αt) This is further averaged over the distribution of initial ( ) α + ω. = 2m k P g(α) sin(αt) g(α) = g(0) + αg (0) + α2 2 g +
6 Landau Damping 61 Physical Picture of Landau Damping Consider a 1-dimensional electrostatic (longitudinal) wave with k E in the absence of magnetic field Taking v = ẑv and E = ẑe cos(kz ωt), the singleparticle
More informationCo-existence and interference of multiple modes in plasma turbulence: Some recent GENE results
Co-existence and interference of multiple modes in plasma turbulence: Some recent GENE results Frank Jenko IPP Garching, Germany University of Ulm, Germany Acknowledgements: F. Merz, T. Görler, D. Told,
More informationCoupled radius-energy turbulent transport of alpha particles
Coupled radius-energy turbulent transport of alpha particles George Wilkie, Matt Landreman, Ian Abel, William Dorland 24 July 2015 Plasma kinetics working group WPI, Vienna Wilkie (Maryland) Coupled transport
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 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 informationStructure formation. Yvonne Y. Y. Wong Max-Planck-Institut für Physik, München
Structure formation Yvonne Y. Y. Wong Max-Planck-Institut für Physik, München Structure formation... Random density fluctuations, grow via gravitational instability galaxies, clusters, etc. Initial perturbations
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 informationA numerical instability in an ADI algorithm for gyrokinetics
Computer Physics Communications 172 2005) 119 132 www.elsevier.com/locate/cpc A numerical instability in an ADI algorithm for gyrokinetics E.A. Belli,G.W.Hammett Princeton Plasma Physics Laboratory, Princeton,
More informationGyrokinetic Turbulence Simulations at High Plasma Beta
Gyrokinetic Turbulence Simulations at High Plasma Beta Moritz J. Pueschel Thanks to F. Jenko and M. Kammerer Ringberg Theory Meeting, Nov. 18, 2008 1 Motivation 2 3 The Beta Parameter Definition β β e
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 informationGyrokinetic Transport Driven by Energetic Particle Modes
Gyrokinetic Transport Driven by Energetic Particle Modes by Eric Bass (General Atomics) Collaborators: Ron Waltz, Ming Chu GSEP Workshop General Atomics August 10, 2009 Outline I. Background Alfvén (TAE/EPM)
More informationTurbulent Transport due to Kinetic Ballooning Modes in High-Beta Toroidal Plasmas
1 TH/P-3 Turbulent Transport due to Kinetic allooning Modes in High-eta Toroidal Plasmas A. Ishizawa 1, S. Maeyama, T.-H. Watanabe 1, H. Sugama 1 and N. Nakajima 1 1 National Institute for Fusion Science,
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 informationGlobal Nonlinear Simulations of Ion and Electron Turbulence Usintg a Particle-In-Cell Approach
Global Nonlinear Simulations of Ion and Electron Turbulence Usintg a Particle-In-Cell Approach S. Jolliet 1), B. F. McMillan 1), T. M. Tran 1), X. Lapillonne 1), L. Villard 1), A. Bottino 2), P. Angelino
More informationAST 553. Plasma Waves and Instabilities. Course Outline. (Dated: December 4, 2018)
AST 553. Plasma Waves and Instabilities Course Outline (Dated: December 4, 2018) I. INTRODUCTION Basic concepts Waves in plasmas as EM field oscillations Maxwell s equations, Gauss s laws as initial conditions
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 informationFine-Scale Zonal Flow Suppression of Electron Temperature Gradient Turbulence
Fine-Scale Zonal Flow Suppression of Electron Temperature Gradient Turbulence S.E. Parker, J.J. Kohut, Y. Chen, Z. Lin, F.L. Hinton and W.W. Lee Center for Integrated Plasma Studies, University of Colorado,
More informationComputer Physics Communications
Computer Physics Communications 181 010) 148 1437 Contents lists available at ScienceDirect Computer Physics Communications www.elsevier.com/locate/cpc On the role of numerical dissipation in gyrokinetic
More informationWater-bag reduced gyrokinetic model for the study of the spatial structure of ITG instability linear modes
Water-bag reduced gyrokinetic model for the study of the spatial structure of ITG instability linear modes PhD work directed by N. Besse with invaluable input from P. Bertrand, E. Gravier, P. Morel,R.
More informationPer Helander. Contributions from: R. Kleiber, A. Mishchenko, J. Nührenberg, P. Xanthopoulos. Wendelsteinstraße 1, Greifswald
Rotation and zonal flows in stellarators Per Helander Wendelsteinstraße 1, 17491 Greifswald Contributions from: R. Kleiber, A. Mishchenko, J. Nührenberg, P. Xanthopoulos What is a stellarator? In a tokamak
More informationPlasmas as fluids. S.M.Lea. January 2007
Plasmas as fluids S.M.Lea January 2007 So far we have considered a plasma as a set of non intereacting particles, each following its own path in the electric and magnetic fields. Now we want to consider
More informationThe Physics of Fluids and Plasmas
The Physics of Fluids and Plasmas An Introduction for Astrophysicists ARNAB RAI CHOUDHURI CAMBRIDGE UNIVERSITY PRESS Preface Acknowledgements xiii xvii Introduction 1 1. 3 1.1 Fluids and plasmas in the
More informationMacroscopic plasma description
Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion
More informationPlasma 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 informationGyrokinetic simulation of entropy cascade in two-dimensional electrostatic turbulence
J. Plasma Fusion Res. SERIES, Vol. 9 (2010) Gyrokinetic simulation of entropy cascade in two-dimensional electrostatic turbulence T. TATSUNO 1), M. ARNES 1,2,3), S. C. COWLEY 3,4), W. DORLAND 1), G. G.
More informationKinetic theory of ions in the magnetic presheath
Kinetic theory of ions in the magnetic presheath Alessandro Geraldini 1,2, Felix I. Parra 1,2, Fulvio Militello 2 1. Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, Oxford
More informationInternational Workshop on the Frontiers of Modern Plasma Physics July On the Nature of Plasma Core Turbulence.
1953-43 International Workshop on the Frontiers of Modern Plasma Physics 14-25 July 2008 On the Nature of Plasma Core Turbulence. F. Jenko Max-Planck Institute fuer Plasmaphysik Garching bei Munchen Germany
More informationEXAMINATION QUESTION PAPER
Faculty of Science and Technology EXAMINATION QUESTION PAPER Exam in: Fys-2009 Introduction to Plasma Physics Date: 20161202 Time: 09.00-13.00 Place: Åsgårdvegen 9 Approved aids: Karl Rottmann: Matematisk
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 informationStudies of Turbulence-driven FLOWs:
Studies of Turbulence-driven FLOWs: a) V ", V Competition in a Tube b) Revisiting Zonal Flow Saturation J.C. Li, P.H. Diamond, R. Hong, G. Tynan University of California San Diego, USA This material is
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 informationLecture 5: Kinetic theory of fluids
Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities
More informationA particle-in-cell method with adaptive phase-space remapping for kinetic plasmas
A particle-in-cell method with adaptive phase-space remapping for kinetic plasmas Bei Wang 1 Greg Miller 2 Phil Colella 3 1 Princeton Institute of Computational Science and Engineering Princeton University
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationGyrokinetic Simulations of Tokamak Microturbulence
Gyrokinetic Simulations of Tokamak Microturbulence W Dorland, Imperial College, London With key contributions from: S C Cowley F Jenko G W Hammett D Mikkelsen B N Rogers C Bourdelle W M Nevins D W Ross
More informationKinetic relaxation models for reacting gas mixtures
Kinetic relaxation models for reacting gas mixtures M. Groppi Department of Mathematics and Computer Science University of Parma - ITALY Main collaborators: Giampiero Spiga, Giuseppe Stracquadanio, Univ.
More informationNumKin, Strasbourg, October 17 th, 2016
F. Palermo 1 A.Biancalani 1, C.Angioni 1, F.Zonca 2, A.Bottino 1, B.Scott 1, G.D.Conway 1, E.Poli 1 1 Max Planck Institut für Plasmaphysik, Garching, Germany 2 ENEA C. R. Frascati - Via E. Fermi 45, CP
More informationKinetic Models for Granular Flow
Journal of Statistical Physics, Vol. 97, Nos. 12, 1999 Kinetic Models for Granular Flow J. Javier Brey, 1 James W. Dufty, 2 and Andre s Santos 3 Received September 28, 1998; final March 29, 1999 The generalization
More informationKinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles
Kinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3. The 8th International Conference on Computational
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 informationNONLINEAR MHD WAVES THE INTERESTING INFLUENCE OF FIREHOSE AND MIRROR IN ASTROPHYSICAL PLASMAS. Jono Squire (Caltech) UCLA April 2017
NONLINEAR MHD WAVES THE INTERESTING INFLUENCE OF FIREHOSE AND MIRROR IN ASTROPHYSICAL PLASMAS Jono Squire (Caltech) UCLA April 2017 Along with: E. Quataert, A. Schekochihin, M. Kunz, S. Bale, C. Chen,
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 informationBernstein-Greene-Kruskal (BGK) Modes in a Three Dimensional Unmagnetized Plasma
Bernstein-Greene-Kruskal (BGK) Modes in a Three Dimensional Unmagnetized Plasma C. S. Ng & A. Bhattacharjee Space Science Center University of New Hampshire Space Science Center Seminar, March 9, 2005
More informationKinetic Alfvén waves in space plasmas
Kinetic Alfvén waves in space plasmas Yuriy Voitenko Belgian Institute for Space Aeronomy, Brussels, Belgium Solar-Terrestrial Center of Excellence, Space Pole, Belgium Recent results obtained in collaboration
More informationPHYSICS Computational Plasma Physics
PHYSICS 78 - Computational Plasma Physics INSTRUCTOR Dr. Earl Scime (escime@wvu.edu) 93-34, ext. 1437 Office hours: MW :30 3:30 and whenever door is open Rm 18 & 05 Hodges Hall Class: MWF 1:30-:0 Rm 334
More informationAPPENDIX Z. USEFUL FORMULAS 1. Appendix Z. Useful Formulas. DRAFT 13:41 June 30, 2006 c J.D Callen, Fundamentals of Plasma Physics
APPENDIX Z. USEFUL FORMULAS 1 Appendix Z Useful Formulas APPENDIX Z. USEFUL FORMULAS 2 Key Vector Relations A B = B A, A B = B A, A A = 0, A B C) = A B) C A B C) = B A C) C A B), bac-cab rule A B) C D)
More informationThe Weibel Instability in Collisionless Relativistic Shocks
The Weibel Instability in Collisionless Relativistic Shocks Tanim Islam University of Virginia The Weibel instability[1], purely electromagnetic in nature, has been used to explain the strong magnetic
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 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 informationA Numerical Method for Parallel Particle Motions in Gyrokinetic Vlasov Simulations )
A Numerical Method for Parallel Particle Motions in Gyrokinetic Vlasov Simulations ) Shinya MAEYAMA, Akihiro ISHIZAWA 1), Tomo-Hiko WATANABE 1), Noriyoshi NAKAJIMA 1), Shunji TSUJI-IIO and Hiroaki TSUTSUI
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 informationarxiv:astro-ph/ v2 19 Nov 2002
Astronomy & Astrophysics manuscript no. A AMRFF February, 8 (DOI: will be inserted by hand later) Jeans gravitational instability and nonextensive kinetic theory J. A. S. Lima, R. Silva, J. Santos Departamento
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 informationLectures on basic plasma physics: Introduction
Lectures on basic plasma physics: Introduction Department of applied physics, Aalto University Compiled: January 13, 2016 Definition of a plasma Layout 1 Definition of a plasma 2 Basic plasma parameters
More informationIon vs. Electron Heating by Plasma Turbulence
WPI, Vienna 30 July 2018 Ion vs. Electron Heating by Plasma Turbulence (an accretion-disc problem that opened the plasmagates ) Yohei Kawazura, Michael Barnes & Alex Schekochihin (Oxford) with thanks to
More informationAnomalous transport of particles in Plasma physics
Anomalous transport of particles in Plasma physics L. Cesbron a, A. Mellet b,1, K. Trivisa b, a École Normale Supérieure de Cachan Campus de Ker Lann 35170 Bruz rance. b Department of Mathematics, University
More informationFrom Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray
From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray C. David Levermore Department of Mathematics and Institute for Physical Science and Technology University of Maryland, College Park
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 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 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 informationKinetic Theory of the Presheath and the Bohm Criterion
UW-CPTC 10-3 Kinetic Theory of the Presheath and the Bohm Criterion S D Baalrud and C C Hegna Department of Engineering Physics, University of Wisconsin-Madison, 1500 Engineering Drive, Madison, WI 53706,
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationProperties of freely decaying and driven turbulence of fusion plasmas using gyrokinetic particle simulation
J. Plasma Fusion Res. SERIES, Vol. 9 () Properties of freely decaying and driven turbulence of fusion plasmas using gyrokinetic particle simulation R. Ganesh Institute for Plasma Research, Bhat Village,
More informationEntropy evolution and dissipation in collisionless particle-in-cell gyrokinetic simulations
Max-Planck-Insititut für Plasmaphysik Entropy evolution and dissipation in collisionless particle-in-cell gyrokinetic simulations A. Bottino Objectives Develop a numerical tool able to reproduce and predict
More informationIntroduction. Chapter Plasma: definitions
Chapter 1 Introduction 1.1 Plasma: definitions A plasma is a quasi-neutral gas of charged and neutral particles which exhibits collective behaviour. An equivalent, alternative definition: A plasma is a
More informationGyrokinetic simulations with GYSELA: Main current issues in physics & numerics
Gyrokinetic simulations with GYSELA: Main current issues in physics & numerics Y. Sarazin, Y. Asahi 2, N. Bouzat, G. Dif-Pradalier, P. Donnel, C. Ehrlacher, C. Emeriau 3, X. Garbet, Ph. Ghendrih, V. Grandgirard,
More informationOptimal design of 2-D and 3-D shaping for linear ITG stability*
Optimal design of 2-D and 3-D shaping for linear ITG stability* Mordechai N. Rorvig1, in collaboration with Chris C. Hegna1, Harry E. Mynick2, Pavlos Xanthopoulos3, and M. J. Pueschel1 1 University of
More informationThe gyrokinetic turbulence code GENE - Numerics and applications
Contributors: T. Dannert (1), F. Jenko (1),F. Merz (1), D. Told (1), X. Lapillonne (2), S. Brunner (2), and others T. Görler (1) The gyrokinetic turbulence code GENE - Numerics and applications (1) Max-Planck-Institut
More informationComputational Methods in Plasma Physics
Computational Methods in Plasma Physics Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Purpose of Talk Describe use of numerical methods to solve simple problem in plasma
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 informationFundamentals of Plasma Physics Transport in weakly ionized plasmas
Fundamentals of Plasma Physics Transport in weakly ionized plasmas APPLAuSE Instituto Superior Técnico Instituto de Plasmas e Fusão Nuclear Luís L Alves (based on Vasco Guerra s original slides) 1 As perguntas
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 informationBursty Transport in Tokamaks with Internal Transport Barriers
Bursty Transport in Tokamaks with Internal Transport Barriers S. Benkadda 1), O. Agullo 1), P. Beyer 1), N. Bian 1), P. H. Diamond 3), C. Figarella 1), X. Garbet 2), P. Ghendrih 2), V. Grandgirard 1),
More informationFrequency spectra at large wavenumbers in two-dimensional Hasegawa-Wakatani turbulence
Frequency spectra at large wavenumbers in two-dimensional Hasegawa-Wakatani turbulence Juhyung Kim and P. W. Terry Department of Physics, University of Wisconsin-Madison October 30th, 2012 2012 APS-DPP
More informationLandau Damping Simulation Models
Landau Damping Simulation Models Hua-sheng XIE (u) huashengxie@gmail.com) Department of Physics, Institute for Fusion Theory and Simulation, Zhejiang University, Hangzhou 310027, P.R.China Oct. 9, 2013
More informationTowards Multiscale Gyrokinetic Simulations of ITER-like Plasmas
Frank Jenko Max-Planck-Institut für Plasmaphysik, Garching Universität Ulm Towards Multiscale Gyrokinetic Simulations of ITER-like Plasmas 23 rd IAEA Fusion Energy Conference 11-16 October 2010, Daejeon,
More information12. MHD Approximation.
Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal
More informationGlobal particle-in-cell simulations of Alfvénic modes
Global particle-in-cell simulations of Alfvénic modes A. Mishchenko, R. Hatzky and A. Könies Max-Planck-Institut für Plasmaphysik, EURATOM-Association, D-749 Greifswald, Germany Rechenzentrum der Max-Planck-Gesellschaft
More informationThe Euler Equation of Gas-Dynamics
The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u
More informationBenchmarks in Computational Plasma Physics
Benchmarks in Computational Plasma Physics P. Londrillo INAF, Bologna, Italie S. Landi Università di Firenze, Italie What you compute when you do computations of the Vlasov equation? Overview A short review
More informationGTC Simulation of Turbulence and Transport in Tokamak Plasmas
GTC Simulation of Turbulence and Transport in Tokamak Plasmas Z. Lin University it of California, i Irvine, CA 92697, USA and GPS-TTBP Team Supported by SciDAC GPS-TTBP, GSEP & CPES Motivation First-principles
More informationOn the Boltzmann equation: global solutions in one spatial dimension
On the Boltzmann equation: global solutions in one spatial dimension Department of Mathematics & Statistics Colloque de mathématiques de Montréal Centre de Recherches Mathématiques November 11, 2005 Collaborators
More informationUne décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion
Une décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion Anaïs Crestetto 1, Nicolas Crouseilles 2 et Mohammed Lemou 3 La Tremblade, Congrès SMAI 2017 5
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 information