Monte-Carlo finite elements gyrokinetic simulations of Alfven modes in tokamaks.
|
|
- Clifton Wilfred Wilkerson
- 5 years ago
- Views:
Transcription
1 Monte-Carlo finite elements gyrokinetic simulations of Alfven modes in tokamaks. A. Bottino 1 and the ORB5 team 1,2,3 1 Max Planck Institute for Plasma Physics, Germany, 2 Swiss Plasma Center, Switzerland, 3 University of Warwick, UK NUMKIN 2016 : International Workshop on Numerical Methods for Kinetic Equations with contributions from: A. Biancalani 1, J. Dominski 2, E. Sonnendruecker 1 and N. Tronko 1 in the framework of the Eurofusion ER projects: VeriGyro (E. Sonnendruecker) and NLED (F. Zonca)
2 Outline GOAL: Illustrate the progress made with the electromagnetic global gyrokinetic code ORB5. Going from turbulence to MHD: Kinetic effects on MHD instabilities in tokamaks Kinetic destabilisation of MHD-stable modes, energetic particle physics EM mictroturbulence Interaction between MHD and microturbulence
3 A problem to face... Trustable nonlinear electromagnetic gyrokinetic simulations require a robust self-consistent discrete GK model, obvious for electrostatic, not at all for electromagnetic 1) Alfvén dynamics opens a lot of new channels for energy transfer in the system. Interaction between Alfvén waves and energetic particles. The model and the code must know where the energy goes. 2) kinetic electrons (real mass ratio) do not forgive careless numerical discretisation (e.g. boundary conditions,...). On the plus side: consistent physics model when kinetic electrons are included, no artifacts as the electrostatic limit of the Alfven oscillations, a.k.a. Ω H mode.
4 The discrete ORB5 model from GK field theory Discrete electromagnetic gyrokinetic equations are suited for simulations if 1) preserve symmetries: conserved quantities (energy). 2) contain (only) relevant physics: approximation are needed, but must not break self-consistency. The ORB5 model is based on the following: Establish a GK Lagrangian. Discretise the Lagrangian. Construct discrete equations from the discrete Lagrangian. The resulting discrete Vlasov-Maxwell equations will keep the nice symmetry properties of the discrete Lagrangian.
5 Particle Lagrangian with time dependence in Hamiltonian STARTING POINT: Gyrocenter coordinate (Lie, Gauge,...) transformed low-frequency particle Lagrangian ( e ) L p c A + p b Ṙ + mc e µ θ H R, gyrocenter positions; µ mv 2 2B, magnetic moment; p mu e c J 0A, canonical parallel momentum; θ, gyroangle; B = A, background (static) magnetic field; J 0 gyroaverage operator, U parallel velocity. Coordinate transform method is rather general: choices can be made to arrange L p so that the symplectic part depends only on the background, while all the time varying fields are contained in the Hamiltonian H. [Hahm 1988, Brizard 2007, Miyato 2009, Scott 2015,...].
6 GK total Lagrangian contains all the needed physics Following [Sugama 2000], Lagrangian for particles AND fields is: L = dw 0 dv f(z 0, t 0 )L p (Z(Z 0, t 0 ; t), Ż(Z 0, t 0 ; t), t) sp + dv E2 B 2 8π Z (R, p, µ, θ); dw 2π/(m 2 )B dp dµ; B 2 = A 2 f (Z 0, t 0 ) is the distribution function for the species sp at an arbitrary initial time t 0. L p is the particle Lagrangian written in terms of the gyro-center coordinates Z. The second term is the Lagrangian for the perturbed electromagnetic fields.
7 Total Lagrangian with second order Hamiltonian L = sp L p = dwdv f(z, t)l p + ( e c A + p b) Ṙ + mc e µ θ H H = m U2 2 + µb + ej 0Φ mc2 2B 2 Φ 2 dv E2 B 2 8π H = H 0 + H 1 + H 2, 2 nd order in both fields: p mu e c J 0A H 0 p2 2m + µb H 1 e(j 0 Φ p mc J 0A ) ej 0 Ψ e 2 H 2 2mc 2 (J 0A ) 2 mc2 2B 2 Φ 2 Contains reduced MHD [Miyato 2009, Krommes 2013].
8 Total Lagrangian with second order Hamiltonian L = sp L p = dwdv f(z, t)l p + ( e c A + p b) Ṙ + mc e µ θ H H = m U2 2 + µb + ej 0Φ mc2 2B 2 Φ 2 dv E2 B 2 8π H = H 0 + H 1 + H 2, 2 nd order in both fields: p mu e c J 0A H 0 p2 2m + µb H 1 e(j 0 Φ p mc J 0A ) ej 0 Ψ e 2 H 2 2mc 2 (J 0A ) 2 mc2 2B 2 Φ 2 Contains reduced MHD [Miyato 2009, Krommes 2013].
9 Total Lagrangian with second order Hamiltonian L = sp L p = dwdv f(z, t)l p + ( e c A + p b) Ṙ + mc e µ θ H dv E2 B 2 8π H = p2 2m + µb + e(j 0Φ p mc J 0A ) + e2 2mc 2 (J 0A ) 2 mc2 2B 2 Φ 2 This is all we need from GK... from now on, field theory. In the context of field theory, this Lagrangian can be further approximated, without loosing self-consistency and energetic consistency of the final equations. Symmetry property of the Lagrangian will be automatically transferred to the equations.
10 Quasi-neutrality approximation L = sp L p = dwdv f(z, t)l p + ( e c A + p b) Ṙ + mc e µ θ H ( ) E 2 dv 8π B2 8π H = p2 2m + µb + e(j 0Φ p mc J 0A ) + e2 2mc 2 (J 0A ) 2 mc2 2B 2 Φ 2 Actually it is not directly a statement about quasi-neutrality, statement about relative magnitude of electric field energy and particle v E B kinetic energy.
11 L = sp L p = Quasi-neutrality approximation dwdv f(z, t)l p + ( e c A + p b) Ṙ + mc e µ θ H ( ) E 2 dv 8π B2 8π H = p2 2m + µb + e(j 0Φ p mc J 0A ) + e2 2mc 2 (J 0A ) 2 mc2 2B 2 Φ 2 dv E2 8π + dwdv f m 2 c 2 B 2 Φ 2 = 1 8π λ 2 d k BT e Debye length; ρ 2 4πne 2 S k BT emc 2 e 2 B 2 ( ) dv 1 + ρ2 S λ 2 Φ 2 d ion sound Larmor radius. In fusion plasmas : ρ 2 S λ 2 d = 4πnmc2 B 2 where v a is the Alfvén velocity, c speed of light. = c2 v 2 a 1
12 Linearised polarisation approximation L = sp (( e ) dvdw c A + p b Ṙ + mc e µ θ ) H f dv B2 8π Start from H = H 0 + H 1 + H 2 In the Lagrangian H 0 + H 1 only multiplies f : (H 0 + H 1 )f For H 2, f is replaced by an equilibrium distribution function f M independent of time: H 2 f M L = sp (( e ) dvdw c A + p b Ṙ + mc e µ θ ) H 0 H 1 f + sp dvdwh 2 f M This approximation will lead to linearised field equations. dv B2 8π
13 Discretized Lagrangian: Monte Carlo + finite elements L = sp (( e ) dvdw c A + p b Ṙ + mc e µ θ ) H 0 H 1 f + sp dvdwh 2 f M dv B2 8π If in this expression we discretize all the integrals in f via Monte-Carlo, discretize all the dynamic fields using finite elements, we will automatically get the discrete model of the ORB5 code. No need to introduce markers as physically meaningful quantities: Note: PIC from a Lagrangian is not new... [Lewis 1970,..]
14 Monte-Carlo approximation, operative definitions The principle of a Monte Carlo method is to approximate the expected value of a random variable by an average over a large number of samples. For our purposes a random variable X is a function that can take values in R which are distributed according to a probability density function. Vlasov-Maxwell GK: the random variable X is a phase space position, Z, distributed according to the initial particle phase-space density f. PIC language: each marker can then be seen as one realization (i.e. one random draw) of the random variable Z.
15 Monte-Carlo approximation, operative definitions Assume that a random variable X of probability density function f is available. Arbitrary smooth functions ψ of X define new random variables, and we can compute expected values of these using the definition: E(ψ(X )) = ψ(x)f (x) dx. The Monte Carlo method consists in approximating the expected value of random variables by a sample average. To approximate E(ψ(X )) we consider a sequence of independent random variables (X i ) i distributed like X and approximate E(ψ(X )) by the sample mean E(ψ(X )) M N = 1 N N ψ(x i ). i=1
16 Discretized Lagrangian: MC Lagrangian Go back to the general formulation of the Lagrangian: L = f sp (Z, t)l sp (Z, Ż)dZ + L F. sp The particle contribution to the Lagrangian is a sum of integrals with densities f sp. Each integral them can be replaced by an expected value of a random variable Z sp (t) of probability density f sp (, t): L = ( ) E L sp (Z sp (t), Ż sp (t)) + L F. sp MC method can be applied by replacing the expected value by a sample mean over a large number of samples being drawn independently according to the initial distributions f sp (, t 0 ): L MC = sp 1 N p L sp (z k (t), ż k (t)) + L F. N p k=1
17 Discrete Lagrangian: finite elements Basis functions Λ ν of the finite dimensional function space, all the functions in this space can be expressed as N g Φ h (x, t) = Φ µ (t)λ µ (x); A h (x, t) = A µ (t)λ µ (x) µ=1 N g µ=1 Replacing Φ by Φ h and A by A h in the Lagrangian: + sp L h,mc = sp 1 N p N p k=1 (( e ) w k c A(R k) + p,k b(r k ) Ṙk + mc e µ θ k k m U2 2 µ kb(r k ) ej 0 Φ h (R k ) ep,k mc J 0A h (R k ) dω mc2 2B 2 Φ h 2 f M + sp e2 dω 2mc 2 (J 0A h ) 2 f M ) dv 1 8π A h 2.
18 GK field and particle equations from variational principles Action functional I [Sugama 2000,...]: δi [Z, Φ, A ] = t2 t1 δl[z, Φ], A dt, Functional derivative of L with respect to a function ψ is defined for an arbitrary variation δψ in the same function space as ψ δl δψ δψ = d L(ψ + ɛδψ) L(ψ) L(ψ + ɛδψ) = lim. dɛ ɛ=0 ɛ 0 ɛ Particle equations of motion can obtained by taking the functional derivatives with respect to the particle phase space positions Z = (R, p, µ): δi δz = 0 δl δz = 0 as t 1 and t 2 are arbitrary.
19 GK field and particle equations from variational principles The equations of motion for the gyrocenters are Ṙ = (H 0 + H 1 ) B p B p = B B (H 0 + H 1 ) where we denote by B = B b. Polarization equation: Induction equation: c eb δi δφ = 0 δl δφ = 0 δi δa = 0 δl δa = 0 b (H 0 + H 1 )
20 Variational principle on the discrete Lagrangian, hints Field equations are obtained by taking the variation of the spline coefficients Φ µ (t) and A µ (t). Note: field equations for A µ (t) are difficult to solve properly. At the moment, ORB5 uses methods developed based on control variates by Mishchenko, Hatzky et al. Functional derivatives become partial derivatives [Bottino, Sonnendruecker JPP 2015; Tronko et al., to be submitted, 2016]. As a results, the weak form of the field equation (Galerkin) is automatically obtained. The combination of GK field theory, finite elements and Monte Carlo provides a self-consistent discrete model, with nice symmetry properties, perfect for code developping. Physics can be added by modifying the Lagrangian step by step.
21 Polarization equation from the discrete Lagrangian 0 = L = Φ ν sp N g Φ µ µ=1 sp 1 N p N p N g Φ h (x, t) = Φ µ (t)λ µ (x) µ=1 k=1 w k ( ej 0 Λ ν (R k ))+ sp dω mc2 B 2 Λ ν Λ µ = sp A µν Φ µ = b ν A µν = sp µ dω mc2 N B 2 g Λ ν ( Φ µ Λ µ ) 1 N p N p k=1 dω mc2 B 2 Λ ν Λ µ µ=1 w k (ej 0 Λ ν (R k )) b ν = sp 1 N p w k (ej 0 Λ ν (R k )). N p k=1
22 Discretized Lagrangian, comments No time discretisation has been performed, semi-discrete As the semi-discrete Lagrangian is still invariant with respect to a time translation, the corresponding Noether energy is exactly conserved. Time discretisation in the equation breaks the exact energy conservation. Time integrators: RK or symplectic Gauss-Legendre method [Krauss 2014] Where are the markers, the δ-functions the phase-space volumes associated to the markers??? Of course, the same exact discrete quations could be reconstructed in the traditional way. The PIC formalism is entirely reformulated as Monte Carlo, all the noise reduction methods of the MC community are available: control variates (δf PIC), importance sampling, dissipation...
23 ITPA EP-TAE linear benchmark Benchmark of gyrokinetic, kinetic MHD and gyrofluid codes for the linear calculation of energetic particle (EP) driven Toroidal Alfven Eigenmodes (TAE) dynamics [Köenies, IAEA FEC 2012]. Around 10 codes involved. Circular large aspect ratio plasma, ɛ = 0.1, β e = , ρ = 2/1854. Toroidal mode number n = 6. No energetic particle finite-larmor-radius retained (in this benchmark).
24 . NLED project: ITPA EP TAE linear benchmark. Equilibrium differences, quasi-neutrality was not enforced in ORB5 simulations (flat electron density profiles in ORB5), new simulations on the way... [A. Biancalani et al., PoP 2016]
25 NLED project: linear vs nonlinear EP finite-larmor-radius retained (T EP =500 kev). Wave-particle nonlinearities retained. Nonlinear structure modification observed. φ max (t)/φ max (t=0) 10 2 LIN NL t (Ω i ) x 10 4 [A. Biancalani et al., submitted to PPCF]
26 VeriGyro project: EM global benchmark. Global Cyclone base case. Circular equilibria. 2 Vlasov and 2 PIC codes involved. Linear β scan. T. Goerler, N. Tronko, W. Hornsby et al. PoP 2016 [T. Goerler, N. Tronko, W. Hornsby et al. PoP 2016]
27 VeriGyro project: EM global benchmark. 2 Vlasov and 2 PIC codes involved. Electromagnetic Global Cyclone case. Circular equilibria. Linear toroidal mode number scan. Difference at high n, due to long wavelength Poisson solver in ORB5. [T. Goerler, N. Tronko, W. Hornsby et al. PoP 2016]
28 Global ORB5 code: implementation of a field solver Implementation of a new field solver in order to account for short scales (fine radial structures) [Dominski, PhD 2016] Gyrokinetic quasi-neutrality equation for arbitrary wavelength solver with finite-elements, following a variation of [Z. Lin and W.W. Lee, PRE 1995] and [A. Mishchenko PoP 2005], q j dzδ(x + ρ x)δ φ q j f 0j = q s dzδ(x + ρ x)δf j, B 0 µ }{{} j {kin} }{{} linearized polarization-drift gyro-density j {kin} QNE for Padé approximant in a 2 species plasma qi 2 N i0 φ = q i (1 ρ 2 i ) dzδ(x + ρ x)δf j. m i Ω 2 0i j=i,e QNE for long wavelength approximation in a 2 species plasma qi 2 N i0 φ = q i dzδ(x + ρ x)δf j. m i Ω 2 0i j=i,e
29 Benchmark against the global version of the GENE code ω r,γ R/c s 2 ORB5 (arbitrary) Full line: growth rate 1.5 ORB5 (Pade) ORB5 (long wav. approx.) 1 GENE Dashed line: real frequency k θρ s [Dominski, PhD 2016] n s = 240 n θ = 1024 n ϕ = 512 n avg = 9 nptot = 16M φ θ n = 24 Fully kinetic n s = 800 Eigenmode envelope ORB5 (arbitrary) ORB5 (Pade) GENE ρ vol FWHM 0.3ρ i #part 32M (Padà ), 64M (arb)
30 Conclusions The combination of GK field theory, finite elements and Monte Carlo provides a perfect tool for constructing fully consistent (kinetic electrons) GK EM codes. The discrete equations have the same symmetry properties of the discrete Lagrangian. Physics can be add (or removed) just by changing the Hamiltonian. Most of the Monte Carlo noise reduction methods known in literature can, in principle, be applied.
31 L = sp L p = Quasi-neutrality approximation dwdv f(z, t)l p + ( e c A + p b) Ṙ + mc e µ θ H ( ) E 2 dv 8π B2 8π H = p2 2m + µb + e(j 0Φ p mc J 0A ) + e2 2mc 2 (J 0A ) 2 mc2 2B 2 Φ 2 dv E2 8π + dwdv f m 2 c 2 B 2 Φ 2 = 1 8π λ 2 d k BT e Debye length; ρ 2 4πne 2 S k BT emc 2 e 2 B 2 ( ) dv 1 + ρ2 S λ 2 Φ 2 d ion sound Larmor radius. Fusion plasmas : ρ 2 S λ 2 d = 4πnmc2 B 2 where v a is the Alfvén velocity, c speed of light. = c2 v 2 a 1
32 GK field and particle equations from variational principles Action functional I [Sugama 2000,...]: δi [Z, Φ, A ] = t2 t1 δl[z, Φ], A dt, Functional derivative of L with respect to a function ψ is defined for an arbitrary variation δψ in the same function space as ψ δl δψ δψ = d L(ψ + ɛδψ) L(ψ) L(ψ + ɛδψ) = lim. dɛ ɛ=0 ɛ 0 ɛ Particle equations of motion can obtained by taking the functional derivatives with respect to the particle phase space positions Z = (R, p, µ): δi δz = 0 δl δz = 0 as t 1 and t 2 are arbitrary.
33 GK field and particle equations from variational principles The equations of motion for the gyrocenters are Ṙ = (H 0 + H 1 ) B p B p = B B (H 0 + H 1 ) where we denote by B = B b. Polarization equation: Induction equation: c eb δi δφ = 0 δl δφ = 0 δi δa = 0 δl δa = 0 b (H 0 + H 1 )
34 Vlasov equation Distribution function is conserved along the particle trajectories f (Z(Z 0, t 0 ; t), t) = f (Z 0, t 0 ) yields the GK Vlasov equation by taking the time derivative d dt f (Z(Z 0, t 0 ; t), t) = dz f (Z, t) + t dt f (Z, t) = 0. Z The particle number conservation follows from Liouville s theorem for a time independent Jacobian: Z dz (B dt ) = 0 and the Vlasov equation can be written in the conservative form t ( 2π m 2 B f ) + Z ( 2π dz m 2 B dt f ) = 0. Full derivation, for example: [Miyato 2009].
35 Induction equations n 0 e 2 Σ sp mc 2 A + Σ sp e kbtn 0 4B 2 A 1 4π 2 A = Σ sp The p z formulation implies the appearance of large unphysical terms in the parallel Ampère s law. One can define the adiabatic part of the equlibrium distribution function as dw e p z c m J 0f f 0 = f M + f ad (1) and f ad ef M ( ) J0 Φ p z J 0 A k B Tc and an adiabatic current j ad j ad Σ sp e mc : = Σ sp e 2 mc 2 k B T (2) dwδ(r + ρ x)f ad p z dwδ(r + ρ x)f M p 2 z J 0 A (3)
36 Induction equations note that the simmetry of the Maxwellian in p z implies: ( pz ) 2 dwf M = n 0 k B T (4) m p z dwf M m = 0 (5) for electrons J 0 = 1: j ad e dwf ad p z = cm e 2 mc 2 k B T A n 0 k B T = e2 n 0 mc 2 A f = δf nonad + f M + f ad (6) n 0 e 2 Σ sp mc 2 A + Σ sp e kbtn 0 4B 2 A 1 4π 2 A = (7) Σ sp dw e ( p z c m J 0 δf nonad ef ) M k B Tc p za Σ sp dw p z e 2 f M m k B Tc 2 p zj 0 A
37 Energy theorem Thanks to its derivation from a Lagrangian density which does not directly depend on time, there is a conserved energy. In our case the following energy is conserved: E(t) = ( dω(h 0 + H 1 )f + sp dωh 2 f M ). (8) Let us verify this by direct computation. As H 0 and f M do not depend on time de dt (t) = ( dω(h 0 + H 1 ) f t + dω H 1 t f + dω H ) 2 t f M. sp (9) We first notice that H i t so that ( dω H 1 t f + = δh i δφ Φ t dω H ) 2 t f M = δl δφ Φ = 0. (10) t
38 Energy theorem On the other hand, we have for each species independently, denoting by H e = H 0 + H 1, dωh e f t = 0. This follows from the Hamiltonian structure of the Vlasov equation (Poisson brackets), but can also be verified as follows. As B = 0 and B / p = c/e b, we have (B H2 e f )+ (f (He 2 B )) = B f H2 e B H 2 f e c p p p p e f b H2 e. Integrating over phase space, the terms on the left hand side vanish. On the other hand (11) f b H 2 e = H 2 e (b f + (f b)). (12)
39 using the variational form of the polarisation equation Eq. (??) with δφ = Φ which reads Energy theorem Then, as dvdp dµ H 2 e (fb) = 0, using the GK Vlasov equation we get 2 ( dvdp dµb H f e t = dvdp dµ H 2 e c e b f B f H2 e (13) In the electrostatic, quasi-neutral limit that we consider, we have ) H 2 f e p p H 0 = m U2 2 + µb, H 1 = ej 0 Φ, H 2 = mc2 2B 2 Φ 2. (14) Then our conserved energy from Eq. (8) becomes: E(t) = ) dωf (m U2 2 + µb + ej mc 2 0Φ dωf M 2B 2 Φ 2 sp = ) dωf (m U2 2 + µb mc 2 + dωf M 2B 2 Φ 2 (15) sp
40 sp Energy theorem dωf M mc 2 B 2 Φ 2 = sp dωfej 0 Φ. (16) Note that the same relation can be obtained by multiplying the polarisation equation Eq. (??) by Φ and integrating (by part) over volume. Because of this last relation the energy can be written equivalently E = sp dωf(m U2 2 + µb ej 0Φ) E K + E F (17) with E F 1/2 sp dωefj0 Φ. The power balance equation, also called E B-thermal transfer equations is: de k dt (t) = de F (t) (18) dt
41 Energy theorem It can be verified, using the Euler-Lagrange equations, that de k dt (t) = dωfe (J 0 Φ) Ṙ0, (19) sp where R 0 represents the part of Eq.(??) which does not contain terms in Φ. This quantity can be compared to the time derivative of the field energy: de F dt (t) = d dt ( sp dωf 1 2 ej 0Φ In numerical simulations it is particularly useful to consider the following form of the power balance equation: ) (20) 1 de k E F dt (t) = 1 de F (t). (21) E F dt
Gyrokinetic Field Theory without Lie Transforms
ASDEX Upgrade Gyrokinetic Field Theory without Lie Transforms for realisable tokamak codes B. Scott Max Planck Institut für Plasmaphysik Euratom Association D-85748 Garching, Germany Lectures in Gyrokinetic
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 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 informationTotal-f Gyrokinetic Energetics and the Status of the FEFI Code
ASDEX Upgrade Total-f Gyrokinetic Energetics and the Status of the FEFI Code B. Scott Max Planck Institut für Plasmaphysik Euratom Association D-85748 Garching, Germany JAEA Rokkasho, Feb 2015 Outline
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 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 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 informationGeometric Gyrokinetic Theory and its Applications to Large-Scale Simulations of Magnetized Plasmas
Geometric Gyrokinetic Theory and its Applications to Large-Scale Simulations of Magnetized Plasmas Hong Qin Princeton Plasma Physics Laboratory, Princeton University CEA-EDF-INRIA School -- Numerical models
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 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 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 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 informationWhat place for mathematicians in plasma physics
What place for mathematicians in plasma physics Eric Sonnendrücker IRMA Université Louis Pasteur, Strasbourg projet CALVI INRIA Nancy Grand Est 15-19 September 2008 Eric Sonnendrücker (U. Strasbourg) Math
More informationBounce-averaged gyrokinetic simulations of trapped electron turbulence in elongated tokamak plasmas
Bounce-averaged gyrokinetic simulations of trapped electron turbulence in elongated tokamak plasmas Lei Qi a, Jaemin Kwon a, T. S. Hahm a,b and Sumin Yi a a National Fusion Research Institute (NFRI), Daejeon,
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 informationGlobal gyrokinetic modeling of geodesic acoustic modes and shear Alfvén instabilities in ASDEX Upgrade.
1 EX/P1-18 Global gyrokinetic modeling of geodesic acoustic modes and shear Alfvén instabilities in ASDEX Upgrade. A. Biancalani 1, A. Bottino 1, S. Briguglio 2, G.D. Conway 1, C. Di Troia 2, R. Kleiber
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 informationSummer College on Plasma Physics August Introduction to Nonlinear Gyrokinetic Theory
2052-24 Summer College on Plasma Physics 10-28 August 2009 Introduction to Nonlinear Gyrokinetic Theory T.S. Hahm Princeton Plasma Physics Laboratory Princeton University USA Introduction to Nonlinear
More informationTurbulence in Tokamak Plasmas
ASDEX Upgrade Turbulence in Tokamak Plasmas basic properties and typical results B. Scott Max Planck Institut für Plasmaphysik Euratom Association D-85748 Garching, Germany Uni Innsbruck, Nov 2011 Basics
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 informationGlobal gyrokinetic particle simulations with kinetic electrons
IOP PUBLISHING Plasma Phys. Control. Fusion 49 (2007) B163 B172 PLASMA PHYSICS AND CONTROLLED FUSION doi:10.1088/0741-3335/49/12b/s15 Global gyrokinetic particle simulations with kinetic electrons Z Lin,
More informationZ. Lin University of California, Irvine, CA 92697, USA. Supported by SciDAC GPS-TTBP, GSEP & CPES
GTC Framework Development and Application Z. Lin University of California, Irvine, CA 92697, USA and dgpsttbp GPS-TTBP Team Supported by SciDAC GPS-TTBP, GSEP & CPES GPS-TTBP Workshop on GTC Framework
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 informationGyrokinetic theory for particle transport in fusion plasmas
Gyrokinetic theory for particle transport in fusion plasmas Matteo Valerio Falessi 1,2, Fulvio Zonca 3 1 INFN - Sezione di Roma Tre, Via della Vasca Navale, 84 (00146) Roma (Roma), Italy 2 Dipartimento
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 informationGyrokinetic simulations including the centrifugal force in a strongly rotating tokamak plasma
Gyrokinetic simulations including the centrifugal force in a strongly rotating tokamak plasma F.J. Casson, A.G. Peeters, Y. Camenen, W.A. Hornsby, A.P. Snodin, D. Strintzi, G.Szepesi CCFE Turbsim, July
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 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 informationFinite-Orbit-Width Effect and the Radial Electric Field in Neoclassical Transport Phenomena
1 TH/P2-18 Finite-Orbit-Width Effect and the Radial Electric Field in Neoclassical Transport Phenomena S. Satake 1), M. Okamoto 1), N. Nakajima 1), H. Sugama 1), M. Yokoyama 1), and C. D. Beidler 2) 1)
More informationTokamak Edge Turbulence background theory and computation
ASDEX Upgrade Tokamak Edge Turbulence background theory and computation B. Scott Max Planck Institut für Plasmaphysik Euratom Association D-85748 Garching, Germany Krakow, Sep 2006 Outline Basic Concepts
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 informationParticle flux representations with FLR effects in the gyrokinetic model
Particle flux representations with FLR effects in the gyrokinetic model ジャイロ運動論モデルにおける有限ラーマ 半径効果を含む粒子フラックスの表現 N. Miyato 1), M. Yagi 1), B. Scott 2) 1) Japan Atomic Energy Agency 2) Max Planck Institut
More informationContinuum Edge Gyrokinetic Theory and Simulations 1
1 TH/P6-23 Continuum Edge Gyrokinetic Theory and Simulations 1 X.Q. Xu 1), K. Bodi 2), J. Candy 3), B. I. Cohen 1), R. H. Cohen 1), P. Colella 4), M. R. Dorr 1), J. A. Hittinger 1), G. D. Kerbel 1), S.
More informationElectric and Magnetic Forces in Lagrangian and Hamiltonian Formalism
Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Benjamin Hornberger 1/26/1 Phy 55, Classical Electrodynamics, Prof. Goldhaber Lecture notes from Oct. 26, 21 Lecture held by Prof. Weisberger
More informationVerification and validation of integrated simulation of energetic particles in fusion plasmas I: linear simulations
Page of 0 AUTHOR SUBMITTED MANUSCRIPT - NF- 0 0 0 Verification and validation of integrated simulation of energetic particles in fusion plasmas I: linear simulations S. Taimourzadeh, E. M. Bass, Y. Chen,
More informationBenchmark of gyrokinetic, kinetic MHD and gyrofluid codes for the linear calculation of fast particle driven TAE dynamics
ITR/P-34 Benchmark of gyrokinetic, kinetic MHD and gyrofluid codes for the linear calculation of fast particle driven TAE dynamics A.Könies,S.Briguglio 2,N.Gorelenkov 3,T.Fehér 4,M.Isaev 5,P.Lauber 4,A.Mishchenko,
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 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 informationL Aquila, Maggio 2002
Nonlinear saturation of Shear Alfvén Modes and energetic ion transports in Tokamak equilibria with hollow-q profiles G. Vlad, S. Briguglio, F. Zonca, G. Fogaccia Associazione Euratom-ENEA sulla Fusione,
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 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 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 informationVerification and validation of integrated simulation of energetic particles in fusion plasmas I: linear simulations
Verification and validation of integrated simulation of energetic particles in fusion plasmas I: linear simulations S. Taimourzadeh 1, E. M. Bass 2, Y. Chen 3, C. Collins 4, N. N. Gorelenkov 5, A. Könies
More informationarxiv: v1 [physics.plasm-ph] 3 Apr 2011
A comparison of Vlasov with drift kinetic and gyrokinetic theories arxiv:1104.0427v1 [physics.plasm-ph] 3 Apr 2011 H. Tasso 1, G. N. Throumoulopoulos 2 1 Max-Planck-Institut für Plasmaphysik, Euratom Association,
More informationCharged particle motion in external fields
Chapter 2 Charged particle motion in external fields A (fully ionized) plasma contains a very large number of particles. In general, their motion can only be studied statistically, taking appropriate averages.
More informationSimulation Study of Interaction between Energetic Ions and Alfvén Eigenmodes in LHD
1 Simulation Study of Interaction between Energetic Ions and Alfvén Eigenmodes in LHD Y. Todo 1), N. Nakajima 1), M. Osakabe 1), S. Yamamoto 2), D. A. Spong 3) 1) National Institute for Fusion Science,
More informationarxiv: v2 [physics.plasm-ph] 12 Jul 2011
July 10, 2011 Exact momentum conservation laws for the gyrokinetic Vlasov-Poisson equations Alain J. Brizard 1 and Natalia Tronko 2,3 1 Department of Chemistry and Physics, Saint Michael s College, Colchester,
More informationParticle-in-cell simulations of electron transport from plasma turbulence: recent progress in gyrokinetic particle simulations of turbulent plasmas
Institute of Physics Publishing Journal of Physics: Conference Series 16 (25 16 24 doi:1.188/1742-6596/16/1/2 SciDAC 25 Particle-in-cell simulations of electron transport from plasma turbulence: recent
More informationElectromagnetic theory of turbulent acceleration of parallel flow and momentum conservation. Abstract
Electromagnetic theory of turbulent acceleration of parallel flow and momentum conservation Shuitao Peng and Lu Wang * State Key Laboratory of Advanced Electromagnetic Engineering and Technology, School
More information0 Magnetically Confined Plasma
0 Magnetically Confined Plasma 0.1 Particle Motion in Prescribed Fields The equation of motion for species s (= e, i) is written as d v ( s m s dt = q s E + vs B). The motion in a constant magnetic field
More informationGlobal particle-in-cell simulations of Alfvénic modes
Max-Planck-Intitut für Plamaphyik, EURATOM Aociation WENDELSTEIN 7-X Global particle-in-cell imulation of Alfvénic mode Alexey Mihchenko, Axel Könie, Roman Hatzky Keyword: gyrokinetic, global particle-in-cell,
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 information(a) Write down the total Hamiltonian of this system, including the spin degree of freedom of the electron, but neglecting spin-orbit interactions.
1. Quantum Mechanics (Spring 2007) Consider a hydrogen atom in a weak uniform magnetic field B = Bê z. (a) Write down the total Hamiltonian of this system, including the spin degree of freedom of the electron,
More informationHPC Simulations of Magnetic Fusion Plasmas
HPC Simulations of Magnetic Fusion Plasmas L. Villard S. Brunner, A. Casati, W.A. Cooper, J.P. Graves, M. Jucker, S. Khoshagdam, X. Lapillonne, B.F. McMillan, P. Popovich, O. Sauter, T.M. Tran, T. Vernay
More informationNonlinear hydrid simulations of precessional Fishbone instability
Nonlinear hydrid simulations of precessional Fishbone instability M. Faganello 1, M. Idouakass 1, H. L. Berk 2, X. Garbet 3, S. Benkadda 1 1 Aix-Marseille University, CNRS, PIIM UMR 7345, Marseille, France
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 informationLectures on basic plasma physics: Hamiltonian mechanics of charged particle motion
Lectures on basic plasma physics: Hamiltonian mechanics of charged particle motion Department of applied physics, Aalto University March 8, 2016 Hamiltonian versus Newtonian mechanics Newtonian mechanics:
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 informationarxiv: v1 [physics.plasm-ph] 23 Jan 2019
Gauge-free electromagnetic rokinetic theory J. W. Burby,3 and A. J. Brizard 2 Courant Institute of Mathematical Sciences, New York, NY 2, USA 2 Department of Physics, Saint Michael s College, Colchester,
More informationGuiding Center Orbit Studies in a Tokamak Edge Geometry Employing Boozer and Cartesian Coordinates
Contrib. Plasma Phys. 48, No. -3, 4 8 (8) / DOI./ctpp.839 Guiding Center Orbit Studies in a Tokamak Edge Geometry Employing Boozer and Cartesian Coordinates Y. Nishimura,Y.Xiao,and Z. Lin Department of
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 informationHybrid Kinetic-MHD simulations Status and Updates
in NIMROD simulations Status and Updates Charlson C. Kim 1,2 Yasushi Todo 2 and the NIMROD Team 1. University of Washington, Seattle 2. National Institute for Fusion Science NIMROD Team Meeting Austin,
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 informationNonlinear Zonal Dynamics of Drift and Drift-Alfvén Turbulences in Tokamak Plasmas
Nonlinear Zonal Dynamics of Drift and Drift-Alfvén Turbulences in Toama Plasmas Liu Chen, Zhihong Lin, Roscoe B. White and Fulvio Zonca Department of Physics and Astronomy, University of California, Irvine,
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 informationElectromagnetic Turbulence Simulations with Kinetic Electrons from the the Summit Framework
19th IAEA Fusion Energy Conference Tuesday, October 15, 2002 Paper: TH/P1-13 Electromagnetic Turbulence Simulations with Kinetic Electrons from the the Summit Framework Scott Parker and Yang Chen University
More informationGeneral Relativity and Cosmology Mock exam
Physikalisches Institut Mock Exam Universität Bonn 29. June 2011 Theoretische Physik SS 2011 General Relativity and Cosmology Mock exam Priv. Doz. Dr. S. Förste Exercise 1: Overview Give short answers
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 informationA global collisionless PIC code in magnetic coordinates
A global collisionless PIC code in magnetic coordinates S. Jolliet a,, A. Bottino b, P. Angelino a, R. Hatzky c, T. M. Tran a, B. F. Mcmillan a, O. Sauter a, K. Appert a, Y. Idomura d, L. Villard a a Ecole
More informationCoarse-graining the electron distribution in turbulence simulations of tokamak plasmas
Coarse-graining the electron distribution in turbulence simulations of tokamak plasmas Yang Chen and Scott E. Parker University of Colorado at Boulder Gregory Rewoldt Princeton Plasma Physics Laboratory
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 informationHamiltonian formulation of reduced Vlasov-Maxwell equations
Hamiltonian formulation of reduced Vlasov-Maxell equations Cristel CHANDRE Centre de Physique Théorique CNRS, Marseille, France Contact: chandre@cpt.univ-mrs.fr importance of stability vs instability in
More informationSelf-consistent particle tracking in a simulation of the entropy mode in a Z pinch
Self-consistent particle tracking in a simulation of the entropy mode in a Z pinch K. Gustafson, I. Broemstrup, D. del-castillo-negrete, W. Dorland and M. Barnes Department of Physics, CSCAMM, University
More informationHamiltonian and Non-Hamiltonian Reductions of Charged Particle Dynamics: Diffusion and Self-Organization
NNP2017 11 th July 2017 Lawrence University Hamiltonian and Non-Hamiltonian Reductions of Charged Particle Dynamics: Diffusion and Self-Organization N. Sato and Z. Yoshida Graduate School of Frontier Sciences
More informationEffects of Alpha Particle Transport Driven by Alfvénic Instabilities on Proposed Burning Plasma Scenarios on ITER
Effects of Alpha Particle Transport Driven by Alfvénic Instabilities on Proposed Burning Plasma Scenarios on ITER G. Vlad, S. Briguglio, G. Fogaccia, F. Zonca Associazione Euratom-ENEA sulla Fusione, C.R.
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 informationModelling of Frequency Sweeping with the HAGIS code
Modelling of Frequency Sweeping with the HAGIS code S.D.Pinches 1 H.L.Berk 2, S.E.Sharapov 3, M.Gryaznavich 3 1 Max-Planck-Institut für Plasmaphysik, EURATOM Assoziation, Garching, Germany 2 Institute
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 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 informationProblems in Magnetostatics
Problems in Magnetostatics 8th February 27 Some of the later problems are quite challenging. This is characteristic of problems in magnetism. There are trivial problems and there are tough problems. Very
More informationPenning Traps. Contents. Plasma Physics Penning Traps AJW August 16, Introduction. Clasical picture. Radiation Damping.
Penning Traps Contents Introduction Clasical picture Radiation Damping Number density B and E fields used to increase time that an electron remains within a discharge: Penning, 936. Can now trap a particle
More informationAMSC 664 Final Report: Upgrade to the GSP Gyrokinetic Code
AMSC 664 Final Report: Upgrade to the GSP Gyrokinetic Code George Wilkie (gwilkie@umd.edu) Supervisor: William Dorland (bdorland@umd.edu) May 14, 2012 Abstract Simulations of turbulent plasma in a strong
More informationLecture 16 March 29, 2010
Lecture 16 March 29, 2010 We know Maxwell s equations the Lorentz force. Why more theory? Newton = = Hamiltonian = Quantum Mechanics Elegance! Beauty! Gauge Fields = Non-Abelian Gauge Theory = Stard Model
More informationHigh-order energy conserving, (discontinuous) finite-element algorithms for (gyro) kinetic simulations of plasmas
High-order energy conserving, (discontinuous) finite-element algorithms for (gyro) kinetic simulations of plasmas Ammar H. Hakim 1 G. Hammett 1 E. Shi 1,2 1 Princeton Plasma Physics Laboratory, Princeton,
More informationSynchro-Betatron Motion in Circular Accelerators
Outlines Synchro-Betatron Motion in Circular Accelerators Kevin Li March 30, 2011 Kevin Li Synchro-Betatron Motion 1/ 70 Outline of Part I Outlines Part I: and Model Introduction Part II: The Transverse
More informationGreg Hammett Imperial College, London & Princeton Plasma Physics Lab With major contributions from:
Greg Hammett Imperial College, London & Princeton Plasma Physics Lab With major contributions from: Steve Cowley (Imperial College) Bill Dorland (Imperial College) Eliot Quataert (Berkeley) LMS Durham
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 informationLienard-Wiechert for constant velocity
Problem 1. Lienard-Wiechert for constant velocity (a) For a particle moving with constant velocity v along the x axis show using Lorentz transformation that gauge potential from a point particle is A x
More informationTH/P6-14 Integrated particle simulation of neoclassical and turbulence physics in the tokamak pedestal/edge region using XGC a)
1 TH/P6-14 Integrated particle simulation of neoclassical and turbulence physics in the tokamak pedestal/edge region using XGC a) 1 Chang, C.S., 1 Ku, S., 2 Adams M., 3 D Azevedo, G., 4 Chen, Y., 5 Cummings,
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 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 informationSingle Particle Motion
Single Particle Motion C ontents Uniform E and B E = - guiding centers Definition of guiding center E gravitation Non Uniform B 'grad B' drift, B B Curvature drift Grad -B drift, B B invariance of µ. Magnetic
More informationGyrokinetic simulations of magnetic fusion plasmas
Gyrokinetic simulations of magnetic fusion plasmas Tutorial 1 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr
More informationMHD-particle simulations and collective alpha-particle transport: analysis of ITER scenarios and perspectives for integrated modelling
MHD-particle simulations and collective alpha-particle transport: analysis of ITER scenarios and perspectives for integrated modelling G. Vlad, S. Briguglio, G. Fogaccia, F. Zonca Associazione Euratom-ENEA
More informationMotion of Charged Particles in Fields
Chapter Motion of Charged Particles in Fields Plasmas are complicated because motions of electrons and ions are determined by the electric and magnetic fields but also change the fields by the currents
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 informationPhysics 5153 Classical Mechanics. Canonical Transformations-1
1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant
More informationActive and Fast Particle Driven Alfvén Eigenmodes in Alcator C-Mod
Active and Fast Particle Driven Alfvén Eigenmodes in Alcator C-Mod JUST DID IT. J A Snipes, N Basse, C Boswell, E Edlund, A Fasoli #, N N Gorelenkov, R S Granetz, L Lin, Y Lin, R Parker, M Porkolab, J
More informationA field theory approach to plasma self-organization
WILLIAM E. BOEING DEPARTMENT OF AERONAUTICS & ASTRONAUTICS A field theory approach to plasma self-organization Setthivoine You Jens von der Linden, E. Sander Lavine, Alexander Card, Evan G. Carroll, Manuel
More information