Introduction to Kinetic Simulation of Magnetized Plasma
|
|
- Gavin Barnett
- 5 years ago
- Views:
Transcription
1 Introduction to Kinetic Simulation of Magnetized Plasma Jae-Min Kwon National Fusion Research Institute, Korea 018 EASW8 July 30 Aug 3, 018 1
2 Outline Introduction to kinetic plasma model Very brief on essential things to understand kinetic simulation Reduced kinetic models for magnetized plasma 5D gyrokinetic, 4D bounce-averaged kinetic, 3D fluid with kinetic closure Numerical methods for kinetic simulation of magnetized plasma Particle-in-Cell methods, and related numerical issues Help students understand the basic idea behind the models and read related literatures for further studies
3 Introduction to Kinetic Plasma Model 3
4 Kinetic Plasma Model The most general description of physical system with many particles Each particle satisfies the following equations of motion dd xx dddd = vv, dd vv mm dddd = FF xx, tt = qqee + vv cc BB Number of particles in dd 3 xxdd 3 vv ff xx, vv, tt dd 3 xxxx 3 vv 4
5 Klimontovich Equation Exact description of classical particles interacting with self-consistent electromagnetic forces NN FF xx, vv, tt = δδ xx XX pp tt pp=1 δδ vv VV pp tt dd dddd XX pp tt = VV pp tt, mm ss dd dddd VV pp tt = qqee mm XX pp tt, tt + qq cc VV pp tt BB mm XX pp tt, tt Then, FF satisfies FF + vv FF + qq mm EEmm + vv cc BBmm vv FF = 0 Note that this equation contains whole spatio-temporal scales all the way down to particle distances. 5
6 Klimontovich Equation Since we are not interested in physical phenomena occurring in super micro-scales (actually, the equation itself is not valid in such scales), We separate quantities into two scales i.e. smooth part in large scale and non-smooth part in small scale, and We keep only the smooth part in left hand side and through out all the remaining into the right hand side and call them collision FF = ff + δδδδ, EE mm = EE + δδδδ, BB mm = BB + δδδδ ff + vv ff + qq mm EE + vv cc BB vvff = qq mm δδδδ + vv cc δδδδ vvδδδδ CC ρρ = qq ss ff ss dd vv ss ȷȷ = qq ss ff ss vvdd vv ss BB = 0 BB = EE EE = ρρ εε 0 EE = BB μμ 0 ȷȷ 6
7 Collision Operator Collision term can be derived by following BBGKY hierarchy and truncating higher order interaction terms (Introduction to Plasma Theory, Nicholson) CC ff = vv AAff vv + 1 vv vv : BBff vv AA vv, tt 8ππnn 0ee 4 llllλ mm ee BB vv, tt 4ππnn 0ee 4 llllλ mm ee vv dd vv ff vv, tt vv vv vv vv dd vv vv vv ff vv, tt For high temperature plasmas, the collision frequency becomes very small νν cccccccc nn/tt 3/ In tokamak plasma νν 1 cccccccc = 10 ~ 100 mmmm collisionless 7
8 Kinetic Phase Mixing vv + vv ff = 00 vv xx xx ff ff vv vv 8
9 Collisionless No Collision ff = CC(ff, ff) ff ff ee vv vv vv As Δvv 0, CC ff, ff ~ ff vv increases faster than streaming term Collision becomes important as fine scale structures are developed in velocity space 9
10 Kinetic Phase Mixing vv xx 10
11 Kinetic Plasma Model + vv + FF mm vvff = CC ff, ff Streaming motions of particles + Mutual interactions in large scales Tend to drive complicated phase space structures Mutual interactions in very short scales ~ Random collision Tend to erase fine scale phase space structures Weaker for small density and high temperature system. Stronger for smaller velocity space scales ff ee vv vv vv 11
12 Kinetic vs Fluid Shifted Gaussian distribution function ff xx, vv, tt = ππ 3/ dd vvff nn( xx) TT xx mm = nn 3/ mm vv uu(xx) exp[ ] TT(xx) Navier-Stokes (Fluid) Equation nn uu + uu uu = pp + μμ uu pp = nnnn uu vv tt = TT mm vv 1
13 Spatio-Temporal Scales of Fusion Plasma S.Either et al, IBM J. RES. & DEV. Vol
14 Turbulence in Fusion Device Fusion plasma confined by external magnetic field strongly magnetized plasma ρρ RR 14
15 Turbulence in Fusion Device In Tokamak, anomalous heat and particle transport driven by micro-scale fluctuations with relatively low frequency ff 300 kkkkkk and kk ρρ ii 1 Collisionless plasma νν cc /ff 1 MIR measured nn ee fluctuation on KSTAR L-mode (J.A. Lee et al, PoP 5, 0513(018)) ECEI measured TT ee fluctuation with ff 100 kkkkkk (M.J. Choi, 018) 15
16 Kinetic Plasma Model: Problem Size Problem size for KSTAR plasma Number of grids: NN xx NN yy NN zz NN vvxx NN vvyy NN vvzz ~ Electron-Ion mass ratio ~ 1:3600 time scale disparity ~ 100 Even for limited spatio-temporal scales, simulations based on brute-force approaches are practically impossible Reduced models are essential 16
17 Model Hierarchy for Magnetized Plasma 6D Boltzmann 5D Gyrokinetic (4D Kinetic) ff ss + vv ff ss + qq ss mm ss ρρ = qq ss ff ss ddvv ss EE + vv BB vv ff ss = CC ff ss jj = qq ss ff ss vvddvv ss 3D Fluid (Gyro-Landau-Fluid, MHD) BB = 00 BB = EE EE = ρρ εε 00 EE = BB μμ 00 jj 1D Diffusion 17
18 18 5D Gyrokinetic Model
19 Gyro Motion in Magnetized Plasma Electron Magnetic field μμ vv ΩΩ ii Adiabatic invariant of motion for time scales slower than ΩΩ ii 11 19
20 Basic Idea of Gyrokinetic Model Gyrokinetic orderings Small fluctuation: δδδδ ff 0 ~ eeeeee TT ~ δδδδ BB 0 1 Frieman, Chen, Phys. Fluids 5, 50(198) Hahm, Lee et al, Phys. Fluids 31, 1940(1988) Hahm, Phys. Fluids 31, 670(1988) Brizard, Hahm, Rev. Mod. Phys. 79, 41(007) Low frequency: ωω Ω ii 1 Fast MHD waves and cyclotron waves are ruled out (high freq. GK; Kolesnikov et al, Phys. Plasmas 14, 07506(007)) Anisotropic fluctuation: kk kk 1, kk ρρ ii ~ 1 Mild non-uniformity in plasma profiles, background magnetic field: ρρ ii LL TT,nn 1 Low beta: ββ 1 Free energy to drive turbulence (GK with strong gradient; Shear Alfven wave only (GK with Compressional Alfven; Brizard, Hahm 07 ) Hahm et al, Phys. Plasmas 16, 0305(009)) 0
21 Basic Idea of Gyrokinetic Model Guiding center transformation particle space xx, vv guiding center space XX, vv, μμ, θθ θθ: gyro-angle average out XX = xx ρρ ρρ = bb vv Ω Ω = eebb 0 mmmm ρρ θθ + vv = bb vv μμ = vv BB XX xx vv = vv bb + vv ee 1
22 Basic Idea of Gyrokinetic Model Schematics of guiding center transformation in Gyrokinetic model Particle Space ( xx, vv) Vlasov Equation (6D) Reduced Maxwell Equation (3D) Guiding Center Space ( XX, vv, μμ) GK Equation (5D) Recued Maxwell Equation (3D) Solve Vlasov equation in guiding center space and evaluate sources (nn ss, jj ss ) Transform sources (nn ss, jj ss ) to particle space Solve recued Maxwell equations to obtain EM fields Transform EM fields to guiding center space
23 Gyrokinetic Vlasov Equation for Low-ββ Plasma Transform original 6D Vlasov equation in particle space into guiding center space Take gyro-angle average to remove θθ reduction to 5D ff(xx, vv, μμ, tt), retaining only slow time scales Δtt 1/Ω ii ff + vv bb + μμ bb BB BB + cc bb δδδδ ff BB 0 XX + qq mm bb μμ BB bb δδδδ 1 cc δδaa ff vv = 0 δδδδ = δδδδ vv cc δδaa bb = bb + vv BB bb bb bb = gyro-phase averaged fluctuations 3
24 Gyrokinetic Model Simple View Gyrokinetic description of magnetized plasmas Motion of charged particle d x = v dt d q q v = δφ + v ( B0 + δb) dt m mc Motion of charged ring centered at XX ddxx dddd = vv bb + μμ bb BB BB + cc bb δδδδ BB 0 ddvv dddd = qq bb mm μμ BB + bb δδδδ + 1 cc δδaa 4
25 Gyrokinetic Model Simple View FF BB drift motion of charged particle drift motion of gyro-center in FF BB direction FF BB 5
26 Gyrokinetic Model Simple View Grad-B + Curvature drift dd XX dddd vv bb + δδaa bb cccc + mmvv + μμμμ bb BB qqbb + cc BB 0 bb δδδδ ExB drift Parallel motion along perturbed magnetic field ddvv dddd qq mm μμ bb BB + bb δδδδ + 1 cc mirror force δδaa Parallel E-field GK equations of motion are nothing but a combination of familiar drift motions ensuring phase space volume conservation and making them Hamiltonian flows BB ff + XX dd XX dddd BB ff + vv ddvv dddd BB ff = 0 6
27 Gyrokinetic Model Simple View Poisson equation with enhanced polarization shielding (1 + ρρ ii λλ ) δδδδ xx, tt = 4ππ DDDD ss qq ss NN ss Density from charged rings Additional shielding by polarization charges carried by charged rings Significantly enhanced compared to Debye shielding Ampere equation without displacement current, also for δδ AA = bbδδaa for low-ββ δδbb bbδδaa 4ππ cc SS JJ ss xx ii=1 NNgg δδaa = 4ππ cc ss JJ ss Parallel current carried by charged rings XX BB 0 7
28 8 Gyrokinetic Model Simple View Gyro-averaged potentials δδδδ, δδaa felt by charged ring = + + = + N l l i i i X N d X x x X dx d )) ( ( 1 )) ( ( 1 ) ( ) ) ( ( 1 θ ρ δφ θ θ ρ δφ π δφ θ ρ δ θ π δφ π π X ) ( l ρ i θ or in Fourier space (as is often done in continuum codes) ( ) ( ) ( ) ( ) + Ω = = = + = + k d e v k J k dk e d e k d dk e k d X x x X dx d X ik i X ik ik X ik i i i i cos 3 0 ) ( ) ˆ( 1 ) ˆ( 1 1 ) ˆ( 1 1 )) ( ( 1 ) ( ) ) ( ( 1 δφ π θ δφ π π θ δφ π π θ θ ρ δφ π δφ θ ρ δ θ π δφ π θ ρ π θ ρ π π Integration can be approximated by a few points sum
29 Reduced Problem Size Problem size for KSTAR plasma Number of grids: NN xx NN yy NN zz NN vvxx NN vvyy NN vvzz ~ Number of grids: NN xx NN yy NN zz NN vv NN μμ ~ Electron-Ion mass ratio ~ 1:3600 Time scale disparity ~ 100 Fluid Model? 9
30 4D Bounce-Averaged Kinetic Model 30
31 Bounce Motion in Magnetized Plasma vv FF = μμ BB Trapped particles vv vv = 1 BB BB mmmmmm Passing particles Passing particles vv Passing and trapped particles behave very differently Passing particles move freely along magnetic field line For slow perturbation, fluid model works well Trapped particles show non-trivial responses to slow perturbations Hard to capture in fluid model 31
32 Bounce Motion in Magnetized Plasma JJ 11 = PP dd ll = mmvv 11 qqqqqq RRRRRR = ππππvv ΩΩ μμ JJ = PP dddd mmvv zz dddd = mmvv zz dddd = ππππvv zz, vv zz = vv zz cccccc ωω bb tt ωω bb 3 Figure from Geomagnetism, Nathani Basavaiah
33 Bounce Motion in Magnetized Plasma Trapped electrons in fusion device Motions along banana are very fast detailed position is negligible Toroidal precession motions are slow comparable with ion transit motions Trapped electron bounce-centers behave like ions resonate with ion scale turbulence 33
34 Bounce-Averaged Kinetic Model Bounce-center coordinate Radial position of bounce center: ββ = rr Toroidal angle of bounce center at outer mid-plane: αα = ζζ 1 sin θθ/ Bounce phase: Ψ = ππ + ssssss vv sin κκ pp Ignorable variable, averaged out Bounce invariant: II qqrr 0 mmεεεεbb 0 κκ pp Further reduction of 5D gyrokinetic equation to 4D bounce-averaged kinetic model Fong, Hahm, Phys. Plasmas 6, 188 (1999) Qi, Kwon, Hahm, Jo, Phys. Plasmas 3, 06513(016) Kwon, Qi, Yi, Hahm, Comput. Phys. Commun. 177, 775(017) + ππ/ αα ββ Ψ 34
35 Bounce-Averaged Kinetic Model Hamiltonian and equations of motion: HH ββ, αα, ΙΙ, μμ = μμbb ΔΔ εε + κκ qq εε + εεεεbb 0 qqrr 0 mm II + qq ss φφ bb ddββ = cc φφ bb dddd αα ddαα dddd = cc ee HH 0 cc φφ bb ββ ββ dddd dddd = 0 dddd dddd = 0 Bounce averaged kinetic equation: dddd dddd = + ddββ dddd ββ + ddαα dddd αα = 0 35
36 Bounce-Averaged Kinetic Model Numerical calculation of bounce average φφ bb φφ dddd ll dddd ll = 1 TT φφ dddd ll Approximation of bounce orbit by unperturbed guiding center motion 36
37 Reduced Problem Size Problem size for KSTAR plasma Number of grids: NN xx NN yy NN zz NN vvxx NN vvyy NN vvzz ~ Number of grids: NN xx NN yy NN zz NN vv NN μμ ~ 10 9 Electron-Ion mass ratio ~ 1:3600 time scale disparity ~ 100 Δtt CCCCCC Δtt for ions only 37
38 Fluid Model with Kinetic Closure 38
39 Kinetic Phase Mixing vv + vv ff = 0 vv xx xx ff δδδδ(0) δδδδ tt < δδδδ(0) ff vv Decay rate? vv 39
40 Kinetic Closure: How to mimic kinetic process? Fluid Model For neutral gas + (Fick s law) VVVV = 0 with particle flux Γ = nnnn = DD Then, we have = DD xx nn Let s assume a stationary solution nn xx, tt = nn 0 and put a perturbation at tt = 0 as δδδδ xx, 0 = nn 1 ee iiiiii. Then, the solution of the equation becomes δδδδ xx, tt = δδδδ xx, 0 ee DDkk tt 40
41 Kinetic Closure: How to mimic kinetic process? Kinetic Model For neutral gas + vv = 0 With a stationary solution, ff 0 vv = nn 0 ππvv tt ee vv vv tt + (note that ddddff0 = nn 0 ) Let s put a perturbation at tt = 0, ff 1 xx, vv = nn 1ee iiiiii ππvv tt Then the solution of the kinetic equation becomes ee vv vv tt + (note that ddddff1 = nn 1 ee iiiiii ) ff xx, vv, tt = nn 0 + nn 1 ee iiii xx vvvv 1 ππvv tt ee vv /vv tt 41
42 Kinetic Closure: How to mimic kinetic process? Kinetic Model ff xx, vv, tt = nn 0 + nn 1 ee iiii xx vvvv 1 ππvv tt ee vv /vv tt The density evolution becomes nn xx, tt = dddddd = nn 0 + nn 1 ee iiiiii ππvv ddddee iiiiiiii ee vv vv tt = nn 0 + nn 1 ee iiiiii ee kk vv tt tt / tt i.e. δδδδ xx, tt = δδδδ xx, 0 ee kk vv tt tt / Kinetic vs Fluid kinetic ee kk vv tt tt / fluid ee DDkk tt 4
43 Kinetic Closure: How to mimic kinetic process? kinetic ee kk vv tt tt / dddd ee kk vv tt tt / = kk vv tt ππ fluid ee DDkk tt dddd ee DDkktt = DDkk tt Match two time responses ~ match the two areas DD = ππ vv tt kk Hammett et al, Phy Rev Lett 64, 3019(1990) Hammett et al, Phys Fluids B 4, 05(199) 43
44 Kinetic Closure: How to mimic kinetic process? With DD kk = ππ vv tt kk, the particle flux can be written as ΓΓ kk = DD kk iiiinn kk In real space, Γ = 1 ππ + ddddee iiiiii Γ kk = 1 ππ + ddddee iiiixx iiii nn kk kk Using delta function identities 1 = kk + δδ kkkkk ddddd, δδ kkxx = lim εε 0 1 ππ εε εε +kk xx Γ = vv tt ππ 3/ 0 ddxx nn xx+xx nn xx xx xx ( c.f. conventional closure Γ = DD ) i.e. DD kk is a non-local integral operator 44
45 45 Hammett et al, Phys Fluids B 4, 05(199) Beer, Hammett, Phys. Plasmas 3, 4046(1996) P. Snyder, Ph.D. Thesis (1999)
46 Numerical Methods for Kinetic Plasma Simulation Continuum (Eulerian) Method Discretize 5D/4D phase space, and apply FDM, FVM, FEM Computationally expensive, but enable high quality simulation Particle-in-Cell (Lagrangian) Method Computationally cheap (relative to continuum method) Noise issues Semi-Lagrangian Method 46
47 Particle-in-Cell Method for Kinetic Simulation C.K. Birdsall and A.B. Langdon, Plasma Physics via Computer Simulation, McGraw-Hill, 1985 R.W. Hockney and J.W. Eastwood, Computer Simulation using Particles, IPP, 1988 W.W. Lee, J. Comput. Phys. 7, 43 (1987) 47
48 From Klimontovich Equation to PIC FF + vv FF + qq mm EEmm + vv cc BBmm vv FF = 0 FF xx, vv, tt = NN pp=1 δδ xx XX pp tt δδ vv VV pp tt dd dddd XX pp tt = VV pp tt mm ss dd dddd VV pp tt = qqee mm XX pp tt, tt + qq cc VV pp tt BB mm XX pp tt, tt NN We want simulate NN-particle system with NN ss NN, GG = ss pp=1 δδ xx XXpp tt δδ vv VV pp tt NN ss FF WW xx, vv GG = WW pp δδ xx XX pp tt δδ vv VV pp tt pp=1 {WW pp } depends on marker particle loading scheme i.e. how to set {XX pp 0, VV pp (0)} For example, if GG xx, vv, tt = 0 FF(xx, vv, tt = 0) WW pp = NN NN ss More sophisticated schemes to minimize loading noise: quite starting scheme, optimal loading scheme etc. (J. Denavit and J.M. Walsh, Plasma Phys. Control. Fusion 6, 09 (1981))
49 Poisson Equation: φφ = dddddd φφ xx, tt = jj φφ jj tt SS jj xx = jj φφ jj tt SS xx xx jj xx jj xx xx If we write electrostatic potential as φφ xx, tt = jj φφ jj tt SS jj xx using a set of basis function SS jj xx, the Poisson equation becomes jj φφ jj tt SS jj xx = 4ππππ dddd FF xx, vv, tt = 4ππππ dddd WW pp δδ xx XX pp tt NN pp=0 δδ vv VV pp tt P.M. Prenter, Splines and Variational Methods, John Wiley & Sons
50 Poisson Equation: φφ = dddddd jj φφ jj tt SS jj xx = 4ππππ dddd FF xx, vv, tt = 4ππππ dddd WW pp δδ xx XX pp tt NN pp=0 δδ vv VV pp tt Appling dddd SS ii (xx) on both sides, jj dddd SS ii xx SS jj xx φφ jj (tt) = 4ππππ WW pp SS ii pp XX pp tt XX pp tt Note that if S is spline function, it satisfies SS xx yy = SS yy xx xx jj 1 xx jj xx jj+1 WW pp SS ii XX pp tt = WW pp SS XXpp tt (xx ii ) = dddd WW pp SS XXpp tt (xx ii )δδ xx XX pp tt pp pp pp δδ vv VV pp tt i.e. the distribution function can be interpreted as a set of particles with spatial shape SS xx P.M. Prenter, Splines and Variational Methods, John Wiley & Sons
51 Collision in PIC Simulation Fokker-Planck equation ff vv, tt = αα vv ff vv, tt + [ββ vv ff vv, tt ] vv It can be interpreted as a distribution function for particles with characteristic equation ddvv dddd = αα vv, tt + ββ vv, tt ξξ tt, ξξ tt ξξ tt = δδ tt tt Due to the drag and diffusion on the right hand side, particle trajectories are scattered as Δvv = ααδtt + 3 RR 0.5 ββδtt RR is a random number in [0, 1] J. Grasman and O.A. van Herwaarden, Asymptotic Methods for the Fokker-Planck Equation and the Exit Problem in Applications, Springer
52 Example: nd order Runge-Kutta dddd dddd = ZZ zz, Φ Start Calculate charge density using zz pp nn Solve Poisson equation for Φ nn Initialize Poisson solver Initialize zz pp 0 xx pp 0, vv pp 0 Push Particles zz pp nn+ 1 = zz pp nn + Δtt ZZ( zz pp nn, Φ nn ) Calculate charge density using zz pp nn+ 1 Diagnosis End Solve Poisson equation for Φ nn+1 Push Particles zz pp nn+1 = zz pp nn + Δtt ZZ( zz pp nn+ 1, Φ nn+1 )
53 Numerical Issues in PIC Simulation 53
54 NN-Particle System Temperature
55 NN-Particle System Fourier representation of potential and distribution function ΦΦ xx = ee iiiiii ΦΦ kk, kk ΦΦ kk = 1 LL ddddee iiiiii ΦΦ(xx) ff xx, vv, tt = ww pp δδ xx xx pp δδ(vv vv pp ), pp ff kk tt = 1 LL ddddee iiiiii ff xx, vv, tt = 1 LL pp ww pp ee iiiixx ppδδ vv vv pp Poisson equation in Fourier space ΦΦ = 4ππππ = 4ππππ dddddd kk ΦΦ kk = 4ππππ LL pp ww pp ee iiiixx pp Apply a filter SS kk on RHS ΦΦ kk = kk ΦΦ kk = 4ππee LL SS kk ww pp ee iiiixx pp pp 4ππππ LLkk SS kk ww pp ww ppp ee iiii(xx pp xx pp ) pp ppp
56 NN-Particle System Φ kk = 4ππππ LLkk SS kk ww pp ww ppp ee iiii(xx pp xx pp ) pp ppp For randomly scattered particles and NN 1, we can simplify RHS Φ kk 4ππee LLkk SS kk ww pp = pp 4ππee LLkk SS kk ww NN ee Φ kk TT 4ππee LLTTkk SS kk ww NN = 4ππππee TT 1 nn LL kk 4 SS kk ww NN = SS kk ww λλ 4 DD kk 4 NN ffffffffffffffffffffff eeeeeeeeeeee ttttttttttttt eeeeeeeeeeee ~ 1 NN
57 Dielectric Response Function and Fluctuation Dissipation Theorem Suppose we put an external test charge ρρ eeeeee xx, tt = ρρ eeeeee kk, ωω exp[ii kkkk ωωωω ] kk,ωω Then, plasma will respond to this external perturbation and generate density fluctuation ρρ(kk, ωω) If we write the total charge density as ρρ tttttt (kk, ww), dielectric response function is defined as ρρ tttttt kk, ωω = ρρ eeeeee kk, ωω + ρρ kk, ωω ρρ eeeeee kk, ωω εε kk, ωω, εε kk, ωω ρρ tttttt kk, ωω = ρρ eeeeee kk, ωω Fluctuation Dissipation Theorem tells us how system responds to small fluctuations and provides useful information for fluctuation spectrum LL ππ δδδδ ωω, kk = TT ωω IIII 1 εε(kk, ωω) Lower bound of fluctuation (or noise) for N-body system
58 Fluctuation and Dissipation Theorem: Implication for PIC simulation PIC simulation employs finite number of particles, actually far fewer than the real number of particles (e.g. in medium fusion device NN rrrrrrrr ~ 10 1 NN ssssss = 10 9 ~ 10 1 ) Fluctuation (or noise) from simulation should be much bigger than physical ones However, we have numerical tools to control too big fluctuation (or noise): Number of particles Particle shape function SS kk Spatio-temporal discretization Δxx, Δtt ee Φ kk TT SS kk ww λλ 4 DD kk 4 NN The strategy is to employ these to suppress noise to a level not to affect key physics
59 Noises NN rrrrrrrr NN ssssssssssssssssssss Noises particle shape, filtering, smoothing, errors etc. Noises
60 Summary Reduced kinetic models for magnetized plasma 5D gyrokinetic, 4D bounce-averaged kinetic, 3D fluid with kinetic closure These enable us to perform kinetic simulation with reasonable computing cost Numerical methods for kinetic simulation of magnetized plasma Particle-in-Cell method Fluctuation dissipation theorem to understand discrete particle noise Verification and benchmark test are necessary, and experimental validation should be the next step very active area of on-going research in fusion 60
Lecture 2: Plasma particles with E and B fields
Lecture 2: Plasma particles with E and B fields Today s Menu Magnetized plasma & Larmor radius Plasma s diamagnetism Charged particle in a multitude of EM fields: drift motion ExB drift, gradient drift,
More information(1) Introduction: a new basis set
() Introduction: a new basis set In scattering, we are solving the S eq. for arbitrary VV in integral form We look for solutions to unbound states: certain boundary conditions (EE > 0, plane and spherical
More information(1) Correspondence of the density matrix to traditional method
(1) Correspondence of the density matrix to traditional method New method (with the density matrix) Traditional method (from thermal physics courses) ZZ = TTTT ρρ = EE ρρ EE = dddd xx ρρ xx ii FF = UU
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 informationDependences of Critical Rotational Shear in DIII-D QH-mode Discharges
Dependences of Critical Rotational Shear in DIII-D QH-mode Discharges by T.M. Wilks 1 with K.H. Burrell 2, Xi. Chen 2, A. Garofalo 2, R.J. Groebner 2, P. Diamond 3, Z. Guo 3, and J.W. Hughes 1 1 MIT 2
More informationWork, Energy, and Power. Chapter 6 of Essential University Physics, Richard Wolfson, 3 rd Edition
Work, Energy, and Power Chapter 6 of Essential University Physics, Richard Wolfson, 3 rd Edition 1 With the knowledge we got so far, we can handle the situation on the left but not the one on the right.
More informationRevision : Thermodynamics
Revision : Thermodynamics Formula sheet Formula sheet Formula sheet Thermodynamics key facts (1/9) Heat is an energy [measured in JJ] which flows from high to low temperature When two bodies are in thermal
More informationCHAPTER 5 Wave Properties of Matter and Quantum Mechanics I
CHAPTER 5 Wave Properties of Matter and Quantum Mechanics I 1 5.1 X-Ray Scattering 5.2 De Broglie Waves 5.3 Electron Scattering 5.4 Wave Motion 5.5 Waves or Particles 5.6 Uncertainty Principle Topics 5.7
More informationRadiation. Lecture40: Electromagnetic Theory. Professor D. K. Ghosh, Physics Department, I.I.T., Bombay
Radiation Zone Approximation We had seen that the expression for the vector potential for a localized cuent distribution is given by AA (xx, tt) = μμ 4ππ ee iiiiii dd xx eeiiii xx xx xx xx JJ (xx ) In
More informationDiffusive DE & DM Diffusve DE and DM. Eduardo Guendelman With my student David Benisty, Joseph Katz Memorial Conference, May 23, 2017
Diffusive DE & DM Diffusve DE and DM Eduardo Guendelman With my student David Benisty, Joseph Katz Memorial Conference, May 23, 2017 Main problems in cosmology The vacuum energy behaves as the Λ term in
More informationElastic light scattering
Elastic light scattering 1. Introduction Elastic light scattering in quantum mechanics Elastic scattering is described in quantum mechanics by the Kramers Heisenberg formula for the differential cross
More informationA Posteriori Error Estimates For Discontinuous Galerkin Methods Using Non-polynomial Basis Functions
Lin Lin A Posteriori DG using Non-Polynomial Basis 1 A Posteriori Error Estimates For Discontinuous Galerkin Methods Using Non-polynomial Basis Functions Lin Lin Department of Mathematics, UC Berkeley;
More informationReview for Exam Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa
Review for Exam2 11. 13. 2015 Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa Assistant Research Scientist IIHR-Hydroscience & Engineering, University
More informationAngular Momentum, Electromagnetic Waves
Angular Momentum, Electromagnetic Waves Lecture33: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay As before, we keep in view the four Maxwell s equations for all our discussions.
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 informationQuantum Mechanics. An essential theory to understand properties of matter and light. Chemical Electronic Magnetic Thermal Optical Etc.
Quantum Mechanics An essential theory to understand properties of matter and light. Chemical Electronic Magnetic Thermal Optical Etc. Fall 2018 Prof. Sergio B. Mendes 1 CHAPTER 3 Experimental Basis of
More informationInteraction with matter
Interaction with matter accelerated motion: ss = bb 2 tt2 tt = 2 ss bb vv = vv 0 bb tt = vv 0 2 ss bb EE = 1 2 mmvv2 dddd dddd = mm vv 0 2 ss bb 1 bb eeeeeeeeeeee llllllll bbbbbbbbbbbbbb dddddddddddddddd
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 informationCharge carrier density in metals and semiconductors
Charge carrier density in metals and semiconductors 1. Introduction The Hall Effect Particles must overlap for the permutation symmetry to be relevant. We saw examples of this in the exchange energy in
More informationReview for Exam Hyunse Yoon, Ph.D. Assistant Research Scientist IIHR-Hydroscience & Engineering University of Iowa
57:020 Fluids Mechanics Fall2013 1 Review for Exam3 12. 11. 2013 Hyunse Yoon, Ph.D. Assistant Research Scientist IIHR-Hydroscience & Engineering University of Iowa 57:020 Fluids Mechanics Fall2013 2 Chapter
More informationPhysics 371 Spring 2017 Prof. Anlage Review
Physics 71 Spring 2017 Prof. Anlage Review Special Relativity Inertial vs. non-inertial reference frames Galilean relativity: Galilean transformation for relative motion along the xx xx direction: xx =
More informationRotational Motion. Chapter 10 of Essential University Physics, Richard Wolfson, 3 rd Edition
Rotational Motion Chapter 10 of Essential University Physics, Richard Wolfson, 3 rd Edition 1 We ll look for a way to describe the combined (rotational) motion 2 Angle Measurements θθ ss rr rrrrrrrrrrrrrr
More informationGradient expansion formalism for generic spin torques
Gradient expansion formalism for generic spin torques Atsuo Shitade RIKEN Center for Emergent Matter Science Atsuo Shitade, arxiv:1708.03424. Outline 1. Spintronics a. Magnetoresistance and spin torques
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 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 informationMathematics Ext 2. HSC 2014 Solutions. Suite 403, 410 Elizabeth St, Surry Hills NSW 2010 keystoneeducation.com.
Mathematics Ext HSC 4 Solutions Suite 43, 4 Elizabeth St, Surry Hills NSW info@keystoneeducation.com.au keystoneeducation.com.au Mathematics Extension : HSC 4 Solutions Contents Multiple Choice... 3 Question...
More informationWave Motion. Chapter 14 of Essential University Physics, Richard Wolfson, 3 rd Edition
Wave Motion Chapter 14 of Essential University Physics, Richard Wolfson, 3 rd Edition 1 Waves: propagation of energy, not particles 2 Longitudinal Waves: disturbance is along the direction of wave propagation
More informationVariations. ECE 6540, Lecture 02 Multivariate Random Variables & Linear Algebra
Variations ECE 6540, Lecture 02 Multivariate Random Variables & Linear Algebra Last Time Probability Density Functions Normal Distribution Expectation / Expectation of a function Independence Uncorrelated
More informationReview for Exam Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa
Review for Exam3 12. 9. 2015 Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa Assistant Research Scientist IIHR-Hydroscience & Engineering, University
More informationIon Heating During Local Helicity Injection Plasma Startup in the Pegasus ST
Ion Heating During Local Helicity Injection Plasma Startup in the Pegasus ST M.G. Burke, J.L. Barr, M.W. Bongard, R.J. Fonck, E.T. Hinson, J.M. Perry, J.A. Reusch 57 th Annual Meeting of the APS-DPP University
More informationChemical Engineering 412
Chemical Engineering 412 Introductory Nuclear Engineering Lecture 12 Radiation/Matter Interactions II 1 Neutron Flux The collisions of neutrons of all energies is given by FF = ΣΣ ii 0 EE φφ EE dddd All
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 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 informationCHAPTER 4 Structure of the Atom
CHAPTER 4 Structure of the Atom Fall 2018 Prof. Sergio B. Mendes 1 Topics 4.1 The Atomic Models of Thomson and Rutherford 4.2 Rutherford Scattering 4.3 The Classic Atomic Model 4.4 The Bohr Model of the
More informationLecture No. 5. For all weighted residual methods. For all (Bubnov) Galerkin methods. Summary of Conventional Galerkin Method
Lecture No. 5 LL(uu) pp(xx) = 0 in ΩΩ SS EE (uu) = gg EE on ΓΓ EE SS NN (uu) = gg NN on ΓΓ NN For all weighted residual methods NN uu aaaaaa = uu BB + αα ii φφ ii For all (Bubnov) Galerkin methods ii=1
More informationLecture 3 Transport in Semiconductors
EE 471: Transport Phenomena in Solid State Devices Spring 2018 Lecture 3 Transport in Semiconductors Bryan Ackland Department of Electrical and Computer Engineering Stevens Institute of Technology Hoboken,
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 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 informationCHAPTER 2 Special Theory of Relativity
CHAPTER 2 Special Theory of Relativity Fall 2018 Prof. Sergio B. Mendes 1 Topics 2.0 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 2.10 2.11 2.12 2.13 2.14 Inertial Frames of Reference Conceptual and Experimental
More informationPhotons in the universe. Indian Institute of Technology Ropar
Photons in the universe Photons in the universe Element production on the sun Spectral lines of hydrogen absorption spectrum absorption hydrogen gas Hydrogen emission spectrum Element production on the
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 informationEquation Sheet, Phys 1321 (Exam II), University of Houston, Fall 2016 Instructor: Dr. W. P. Su
vv (tt) = ddrr (tt) dddd vv aaaaaa = Δrr Δtt aa (tt) = ddvv (tt) dddd aa aaaaaa = Δvv Δtt Equation Sheet, Phys 1321 (Exam II), University of Houston, Fall 2016 Instructor: Dr. W. P. Su AAAAAA. ssssssssss
More informationTECHNICAL NOTE AUTOMATIC GENERATION OF POINT SPRING SUPPORTS BASED ON DEFINED SOIL PROFILES AND COLUMN-FOOTING PROPERTIES
COMPUTERS AND STRUCTURES, INC., FEBRUARY 2016 TECHNICAL NOTE AUTOMATIC GENERATION OF POINT SPRING SUPPORTS BASED ON DEFINED SOIL PROFILES AND COLUMN-FOOTING PROPERTIES Introduction This technical note
More informationPHL424: Nuclear fusion
PHL424: Nuclear fusion Hot Fusion 5 10 15 5 10 8 projectiles on target compound nuclei 1 atom Hot fusion (1961 1974) successful up to element 106 (Seaborgium) Coulomb barrier V C between projectile and
More informationApproximate Second Order Algorithms. Seo Taek Kong, Nithin Tangellamudi, Zhikai Guo
Approximate Second Order Algorithms Seo Taek Kong, Nithin Tangellamudi, Zhikai Guo Why Second Order Algorithms? Invariant under affine transformations e.g. stretching a function preserves the convergence
More informationLise Meitner, Otto Hahn. Nuclear Fission Hans-Jürgen Wollersheim
Lise Meitner, Otto Hahn Nuclear Fission Hans-Jürgen Wollersheim Details of the 252 Cf decay α s: 96.9% SF: 3.1% T 1/2 = 2.647 a Q α = 6.217 MeV E α = 6.118 MeV α α α α α-decay of 252 Cf Mass data: nucleardata.nuclear.lu.se/database/masses/
More informationDispersive Media, Lecture 7 - Thomas Johnson 1. Waves in plasmas. T. Johnson
2017-02-14 Dispersive Media, Lecture 7 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasmas as a coupled system Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas
More information2.4 Error Analysis for Iterative Methods
2.4 Error Analysis for Iterative Methods 1 Definition 2.7. Order of Convergence Suppose {pp nn } nn=0 is a sequence that converges to pp with pp nn pp for all nn. If positive constants λλ and αα exist
More informationSecondary 3H Unit = 1 = 7. Lesson 3.3 Worksheet. Simplify: Lesson 3.6 Worksheet
Secondary H Unit Lesson Worksheet Simplify: mm + 2 mm 2 4 mm+6 mm + 2 mm 2 mm 20 mm+4 5 2 9+20 2 0+25 4 +2 2 + 2 8 2 6 5. 2 yy 2 + yy 6. +2 + 5 2 2 2 0 Lesson 6 Worksheet List all asymptotes, holes and
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 informationDressing up for length gauge: Aspects of a debate in quantum optics
Dressing up for length gauge: Aspects of a debate in quantum optics Rainer Dick Department of Physics & Engineering Physics University of Saskatchewan rainer.dick@usask.ca 1 Agenda: Attosecond spectroscopy
More informationPHL424: Nuclear Shell Model. Indian Institute of Technology Ropar
PHL424: Nuclear Shell Model Themes and challenges in modern science Complexity out of simplicity Microscopic How the world, with all its apparent complexity and diversity can be constructed out of a few
More informationSECTION 5: CAPACITANCE & INDUCTANCE. ENGR 201 Electrical Fundamentals I
SECTION 5: CAPACITANCE & INDUCTANCE ENGR 201 Electrical Fundamentals I 2 Fluid Capacitor Fluid Capacitor 3 Consider the following device: Two rigid hemispherical shells Separated by an impermeable elastic
More informationClassical RSA algorithm
Classical RSA algorithm We need to discuss some mathematics (number theory) first Modulo-NN arithmetic (modular arithmetic, clock arithmetic) 9 (mod 7) 4 3 5 (mod 7) congruent (I will also use = instead
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 informationLast Name _Piatoles_ Given Name Americo ID Number
Last Name _Piatoles_ Given Name Americo ID Number 20170908 Question n. 1 The "C-V curve" method can be used to test a MEMS in the electromechanical characterization phase. Describe how this procedure is
More informationMath 171 Spring 2017 Final Exam. Problem Worth
Math 171 Spring 2017 Final Exam Problem 1 2 3 4 5 6 7 8 9 10 11 Worth 9 6 6 5 9 8 5 8 8 8 10 12 13 14 15 16 17 18 19 20 21 22 Total 8 5 5 6 6 8 6 6 6 6 6 150 Last Name: First Name: Student ID: Section:
More informationThe Bose Einstein quantum statistics
Page 1 The Bose Einstein quantum statistics 1. Introduction Quantized lattice vibrations Thermal lattice vibrations in a solid are sorted in classical mechanics in normal modes, special oscillation patterns
More informationOn The Cauchy Problem For Some Parabolic Fractional Partial Differential Equations With Time Delays
Journal of Mathematics and System Science 6 (216) 194-199 doi: 1.17265/2159-5291/216.5.3 D DAVID PUBLISHING On The Cauchy Problem For Some Parabolic Fractional Partial Differential Equations With Time
More informationAdvances in stellarator gyrokinetics
Advances in stellarator gyrokinetics Per Helander and T. Bird, F. Jenko, R. Kleiber, G.G. Plunk, J.H.E. Proll, J. Riemann, P. Xanthopoulos 1 Background Wendelstein 7-X will start experiments in 2015 optimised
More informationProf. Dr. Rishi Raj Design of an Impulse Turbine Blades Hasan-1
Prof. Dr. Rishi Raj Design of an Impulse Turbine Blades Hasan-1 The main purpose of this project, design of an impulse turbine is to understand the concept of turbine blades by defining and designing the
More informationGrover s algorithm. We want to find aa. Search in an unordered database. QC oracle (as usual) Usual trick
Grover s algorithm Search in an unordered database Example: phonebook, need to find a person from a phone number Actually, something else, like hard (e.g., NP-complete) problem 0, xx aa Black box ff xx
More informationYang-Hwan Ahn Based on arxiv:
Yang-Hwan Ahn (CTPU@IBS) Based on arxiv: 1611.08359 1 Introduction Now that the Higgs boson has been discovered at 126 GeV, assuming that it is indeed exactly the one predicted by the SM, there are several
More informationSaturation rules for ETG transport in quasilinear transport models
Saturation rules for ETG transport in quasilinear transport models J. Citrin 1, C. Bourdelle 2, N. Bonanomi 3, T. Goerler 4, P. Mantica 3 1 FOM Institute DIFFER, PO Box 6336, 5600 HH, Eindhoven, The Netherlands
More informationChapter 22 : Electric potential
Chapter 22 : Electric potential What is electric potential? How does it relate to potential energy? How does it relate to electric field? Some simple applications What does it mean when it says 1.5 Volts
More informationModeling of a non-physical fish barrier
University of Massachusetts - Amherst ScholarWorks@UMass Amherst International Conference on Engineering and Ecohydrology for Fish Passage International Conference on Engineering and Ecohydrology for Fish
More informationABSTRACT OF THE DISSERTATION
A Study of an N Molecule-Quantized Radiation Field-Hamiltonian BY MICHAEL THOMAS TAVIS DOCTOR OF PHILOSOPHY, GRADUATE PROGRAM IN PHYSICS UNIVERSITY OF CALIFORNIA, RIVERSIDE, DECEMBER 1968 PROFESSOR FREDERICK
More informationControl of Mobile Robots
Control of Mobile Robots Regulation and trajectory tracking Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Organization and
More informationLecture 3. STAT161/261 Introduction to Pattern Recognition and Machine Learning Spring 2018 Prof. Allie Fletcher
Lecture 3 STAT161/261 Introduction to Pattern Recognition and Machine Learning Spring 2018 Prof. Allie Fletcher Previous lectures What is machine learning? Objectives of machine learning Supervised and
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 informationPHY103A: Lecture # 4
Semester II, 2017-18 Department of Physics, IIT Kanpur PHY103A: Lecture # 4 (Text Book: Intro to Electrodynamics by Griffiths, 3 rd Ed.) Anand Kumar Jha 10-Jan-2018 Notes The Solutions to HW # 1 have been
More informationJasmin Smajic1, Christian Hafner2, Jürg Leuthold2, March 23, 2015
Jasmin Smajic, Christian Hafner 2, Jürg Leuthold 2, March 23, 205 Time Domain Finite Element Method (TD FEM): Continuous and Discontinuous Galerkin (DG-FEM) HSR - University of Applied Sciences of Eastern
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 informationECE 6540, Lecture 06 Sufficient Statistics & Complete Statistics Variations
ECE 6540, Lecture 06 Sufficient Statistics & Complete Statistics Variations Last Time Minimum Variance Unbiased Estimators Sufficient Statistics Proving t = T(x) is sufficient Neyman-Fischer Factorization
More informationCHEMISTRY Spring, 2018 QUIZ 1 February 16, NAME: HANS ZIMMER Score /20
QUIZ 1 February 16, 018 NAME: HANS ZIMMER Score /0 [Numbers without decimal points are to be considered infinitely precise. Show reasonable significant figures and proper units. In particular, use generally
More informationCold atoms in optical lattices
Cold atoms in optical lattices www.lens.unifi.it Tarruel, Nature Esslinger group Optical lattices the big picture We have a textbook model, which is basically exact, describing how a large collection of
More information7.3 The Jacobi and Gauss-Seidel Iterative Methods
7.3 The Jacobi and Gauss-Seidel Iterative Methods 1 The Jacobi Method Two assumptions made on Jacobi Method: 1.The system given by aa 11 xx 1 + aa 12 xx 2 + aa 1nn xx nn = bb 1 aa 21 xx 1 + aa 22 xx 2
More informationAcceleration to higher energies
Acceleration to higher energies While terminal voltages of 20 MV provide sufficient beam energy for nuclear structure research, most applications nowadays require beam energies > 1 GeV How do we attain
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 informationMaterials & Advanced Manufacturing (M&AM)
Modeling of Shear Thickening Fluids for Analysis of Energy Absorption Under Impulse Loading Alyssa Bennett (University of Michigan) Nick Vlahopoulos, PhD (University of Michigan) Weiran Jiang, PhD (Research
More informationChapter 5: Spectral Domain From: The Handbook of Spatial Statistics. Dr. Montserrat Fuentes and Dr. Brian Reich Prepared by: Amanda Bell
Chapter 5: Spectral Domain From: The Handbook of Spatial Statistics Dr. Montserrat Fuentes and Dr. Brian Reich Prepared by: Amanda Bell Background Benefits of Spectral Analysis Type of data Basic Idea
More informationSECTION 7: FAULT ANALYSIS. ESE 470 Energy Distribution Systems
SECTION 7: FAULT ANALYSIS ESE 470 Energy Distribution Systems 2 Introduction Power System Faults 3 Faults in three-phase power systems are short circuits Line-to-ground Line-to-line Result in the flow
More informationFermi Surfaces and their Geometries
Fermi Surfaces and their Geometries Didier Ndengeyintwali Physics Department, Drexel University, Philadelphia, Pennsylvania 19104, USA (Dated: May 17, 2010) 1. Introduction The Pauli exclusion principle
More informationWorksheets for GCSE Mathematics. Algebraic Expressions. Mr Black 's Maths Resources for Teachers GCSE 1-9. Algebra
Worksheets for GCSE Mathematics Algebraic Expressions Mr Black 's Maths Resources for Teachers GCSE 1-9 Algebra Algebraic Expressions Worksheets Contents Differentiated Independent Learning Worksheets
More informationThe impact of hot charge carrier mobility on photocurrent losses
Supplementary Information for: The impact of hot charge carrier mobility on photocurrent losses in polymer-based solar cells Bronson Philippa 1, Martin Stolterfoht 2, Paul L. Burn 2, Gytis Juška 3, Paul
More information14- Hardening Soil Model with Small Strain Stiffness - PLAXIS
14- Hardening Soil Model with Small Strain Stiffness - PLAXIS This model is the Hardening Soil Model with Small Strain Stiffness as presented in PLAXIS. The model is developed using the user-defined material
More informationCoarse-Graining via the Mori-Zwanzig Formalism
Coarse-Graining via the Mori-Zwanzig Formalism George Em Karniadakis & Zhen Li http://www.pnnl.gov/computing/cm4/ Supported by DOE ASCR Mori-Zwanzig can be used to: 1. Derive coarse-grained equations 2.
More informationEstimate by the L 2 Norm of a Parameter Poisson Intensity Discontinuous
Research Journal of Mathematics and Statistics 6: -5, 24 ISSN: 242-224, e-issn: 24-755 Maxwell Scientific Organization, 24 Submied: September 8, 23 Accepted: November 23, 23 Published: February 25, 24
More informationAnalysis of uncertainty sources in DNS of a turbulent mixing layer using Nek5000
Analysis of uncertainty sources in DNS of a turbulent mixing layer using Nek5000 Juan D. Colmenares F. 1 and Svetlana V. Poroseva 2 The University of New Mexico, Albuquerque, NM 87131 Yulia T. Peet 3 Arizona
More informationPhysics 141 Second Mid-Term Examination Spring 2015 March 31, 2015
Physics 141 Second Mid-Term Examination Spring 2015 March 31, 2015 Your Number Solutions LAST Name (print) FIRST Name (print) Signature: UIN #: Your Section: Barkan 9 am Barkan 11 am Goeckner 3 pm PROBLEM
More informationProblem 3.1 (Verdeyen 5.13) First, I calculate the ABCD matrix for beam traveling through the lens and space.
Problem 3. (Verdeyen 5.3) First, I calculate the ABCD matrix for beam traveling through the lens and space. T = dd 0 0 dd 2 ff 0 = dd 2 dd ff 2 + dd ( dd 2 ff ) dd ff ff Aording to ABCD law, we can have
More informationGravitation. Chapter 8 of Essential University Physics, Richard Wolfson, 3 rd Edition
Gravitation Chapter 8 of Essential University Physics, Richard Wolfson, 3 rd Edition 1 What you are about to learn: Newton's law of universal gravitation About motion in circular and other orbits How to
More informationNeoclassical transport
Neoclassical transport Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 28 th January 2013 Dr Ben Dudson Magnetic Confinement Fusion (1 of 19) Last time Toroidal devices
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 informationIntegrated Attitude Determination and Control System via Magnetic Measurements and Actuation
Proceedings of the th WSEAS International Conference on Automatic Control, Modelling and Simulation Integrated Attitude Determination and Control System via Magnetic Measurements and Actuation MOHAMMAD
More informationModule 7 (Lecture 25) RETAINING WALLS
Module 7 (Lecture 25) RETAINING WALLS Topics Check for Bearing Capacity Failure Example Factor of Safety Against Overturning Factor of Safety Against Sliding Factor of Safety Against Bearing Capacity Failure
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 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 informationSECTION 5: POWER FLOW. ESE 470 Energy Distribution Systems
SECTION 5: POWER FLOW ESE 470 Energy Distribution Systems 2 Introduction Nodal Analysis 3 Consider the following circuit Three voltage sources VV sss, VV sss, VV sss Generic branch impedances Could be
More informationSupport Vector Machines. CSE 4309 Machine Learning Vassilis Athitsos Computer Science and Engineering Department University of Texas at Arlington
Support Vector Machines CSE 4309 Machine Learning Vassilis Athitsos Computer Science and Engineering Department University of Texas at Arlington 1 A Linearly Separable Problem Consider the binary classification
More information