Particle in Cell method
|
|
- Louise Lloyd
- 5 years ago
- Views:
Transcription
1 Particle in Cell method Birdsall and Langdon: Plasma Physics via Computer Simulation Dawson: Particle simulation of plasmas Hockney and Eastwood: Computer Simulations using Particles we start with an electrostatic 1d-1v collissionless plasma: Lecture by G. Lapenta Vlasov equation for distribution function f s : f s t + v f s x + q se m s f s v =0 q s, m s : charge and mass of the species Poisson equation: 2 x 2 = charge density follows from distribution function (x, t) = f s (x, v, t)dv s
2 moving finite elements or superparticles: f p f s (x, v, t) = p f p (x, v, t) Ansatz: tensor product of shape functions: f p (x, v, t) =N p S x (x x p (t))s v (v v p (t)) S x, S v : shape functions N p : number of real particles in the superparticle Assumptions: 1. compact support 2. normalisation: = x or v 3. symmetric shapes: S ( p)d =1 S ( p) =S ( p )
3 Selection of particle shape S v (v v p )= (v v p ) b-splines first b-spline: flat-top function b 0 ( ) b 0 ( )= 1 if < 1/2 0 otherwise subsequent b-splines: b l ( )= d b 0 ( )b l 1 ( )
4
5 spatial shape function S x (x x p )=b l x x p p: size of superparticle p Remarks: few PIC codes use splines of order 1 but most use order 0 Vlasov equation looks linear, but is nonlinear due to E-field. Equation for superparticle: f p t + v f p x + qe m f p v =0 The arbitrary functional form chosen for the elements does not satisfy exactly the Vlasov equation. The usual procedure of the finite element method is to require that the moments of the equations be satisfied. Notation:... := dx dv
6 Moment 0 f p t + v f p x + qe m f p v =0 second and third term = 0 = N p = const Moment 1 x xf p t + xv f p x + qxe m f p v =0 last term = 0
7 consider first term: xf p = N p S v (v v p )dv xs x (x x p )dx = N p (x + x p )S x (x )dx = N p x p using symmetry of S x consider second term: Z Z vdv dx = Z vf p dxdv = hvf p i hvf p i = N p Z = N p Z Z vs v (v v p )dv S x (x (v 0 + v p )S v (v 0 )dv 0 x p )dx thus dx p = N p v p using symmetry of S v dt = v p
8 Moment 1 v vf p t + v 2 f p x + qev m f p v =0 second term = 0 first term has already been computed: vf p = N p v p remaining term: qs E dx v m s f p v dv = qs E m s dx f p dv = q s E f s m s This defines average electric field E p acting on computational particle: q s E q s f s = N p E p m s m s with E p = = S v (v S x (x v p )dv x p )E(x)dx S x (x x p )E(x)dx thus dv p dt = q s m s E p
9 Equations of motion Field equations: finite volume approach dn p dt dx p dt dv p dt = 0 = v p = q s m s E p cell centres: x i cell faces: x i+1/2 cell-averaged values i of potential Poisson equation: i+1 2 i + i 1 ( x) 2 = q i with cell-averaged density: i = 1 x x i+1/2 x i 1/2 (x)dx
10 recall definition of b-spline of order 0 xi+1/2 (x)dx = x i 1/2 b 0 x x i x (x)dx recall definition of the density xi+1/2 x i 1/2 (x)dx = p b 0 x x i x S x (x x p )dx define interpolation function W (x i x p )= S x (x x p )b 0 x x i x dx and get: q p x W (x i x p ) with q p = q s N p i = p
11 very simple expression, when we choose: p = x =) xi x p W (x i x p )=b l+1 x Electric field from central di erences E i = i+1 i 1 2 x Efieldisconstantineachcell: E(x) = X i E i b 0 xi x p x recall definition of E p : E p = X i E i Z b 0 x x xi S x (x x i ) and recall definition of interpolation function E p = X i E i W (x x i )
12 Simplest symplectic integrator: leap-frog xp n+1 = xp n + tvp n+1/2 v n+3/2 p = v n+1/2 p + te p (x n p ) second order accurate initial step: Euler forward v 1/2 = v 0 p + te p (x 0 p) Stability c t < x pe < 2 x < De
13 Elektromagnetic PIC FDTD and Yee mesh (1966) Particle push: Boris scheme (1970) conservative charge deposition (1991, 2001) Literature: inspired by lectures by A. Spitkovsky J. Villasenor, Bunemann T. Zh. Esirkepov T. Umeda
14 d dt "mv = q(e + v c # B) FDTD: second order in space and time Yee mesh: div B = 0
15 Equations of motion: leap-frog: looks implicit
16 Solution: Boris (1970) explicit with t t = q t 2 t m Bt =2arctan(t t )
17 Proof: We know: We want to proof: v + v = q t 2m (v+ + v ) B v + v = v s v = v + v t, t = q t 2m B, s = 2t 1+t 2 thus: v + v = v s +(v t) s v s = v B q t 2 2m 1+t 2 = v + B q t 2m = (v + v ) B q t 2m 2 1+t 2 +(v+ v ) 2 1+t t 2 +(v+ v ) 1 1+t 2
18 (v t) s = (v s) t 2 q t 1 = 2m 1+t 2 [(v+ v ) B] B + q t 2m = v + v = (v + v ) 1 1+t 2 + q t 2m t 2 B [(v+ v ) B] 1 1+t 2 (v+ v ) B (v + v ) B = q t 2m [(v+ + v ) B] B =: C B B (C B) =CB 2 BB C = q t 2m [(v+ + v ) B]B 2 = t 2 (v + v )= q t 2 B 2 [(v + + v ) B] q t 2m 2m t 2 = q t 2 B 2 = 2m
19 PIC Vlasov I) x n+1/2 x n 1/2 t = v n ˆf n+1 (x n+1/2, v n )= x ( t) f n (x n 1/2, v n ) II) j n = v n S(x ) x = x n+1/2 + x n 1/2 2 = x n + O( t 2 ) j n = v n f (x n, v n ) f (x n, v n )= f n (x n 1/2, v n )+ˆf n (x n+1/2, v n ) 2 = f n (x n, n v )+O( t 2 ) III) E n+1/2 E n 1/2 t = B n j n IV) B n+1 B n t = E n+1/2 V) Boris Boris v n+1 B t v n = Bn+1 + B n 2 = q m E n+1/2 (x n+1/2 )+ v n+1 + v n = B n+1/2 + O( t 2 ) 2 B (x n+1/2 ) f n+1 (x n+1/2, v n+1 )= v ( t)ˆf (x n+1/2, v n )
20 Charge and current deposition B/t = c( E) = B/t =0 Yee grid "# $ E "t = c# $ (# % B) & 4'# $ J "# "# "t = $% & J continuity equation is again an initial condition
21 ) =) projection method =) parabolic =) hyperbolic
22 If one uses standard weighting, charge is not conserved. Villasenor, Bunemann (1992): count current through faces Problem: lots of if -statements current deposition as expensive as particle push optimized deposit by Umeda (2003): zig-zag curve Esirkepov (2001) arbitrary form factor, easy and precise used in PSC (Ruhl) and EPOCH (UK code based on PSC)
23 Esirkepov: E n+1 E n = + B n+1/2 J n+1/2, dt B n+1/2 B n 1/2 = E n, dt + E n = ρ n, B n+1/2 = 0, u n+1/2 α x n+1 α u n 1/2 α dt x n α dt = un+1/2 α γα n+1/2, γ α = ( 1 + (u α ) 2) 1/2. = 2π q α m e m α e ( E n( x n α,t) + un α B n( ) x n α γ,t), α
24 E n = ( E 1 i,j+1/2,k+1/2,e2 i+1/2,j,k+1/2,e3 i+1/2,j+1/2,k) n, B n+1/2 = ( B 1 i+1/2,j,k,b2 i,j+1/2,k,b3 i,j,k+1/2) n+1/2, ρ n = ρ n i+1/2,j+1/2,k+1/2, J n+1/2 = ( Ji,j+1/2,k+1/2 1, J i+1/2,j,k+1/2 2, J 3 ( ) ( + fi+1,j,k f i,j,k f i,j,k = dx ( fi,j,k f i 1,j,k f i,j,k = ( dx, f i,j+1,k f i,j,k dy, f i,j,k f i,j 1,k dy i+1/2,j+1/2,k) n+1/2, ± ±, f i,j,k+1 f i,j,k dz, f i,j,k f i,j,k 1 dz ), ). = + + = 0, + = + =,
25 dropping indices ρ n+1 ρ n + + J n+1/2 = 0, dt B n+1/2 B n 1/2 = 0. dt + E = ρ and B = 0 are initial conditions continuity equation ρ n+1 i+1/2,j+1/2,k+1/2 ρn i+1/2,j+1/2,k+1/2 dt + J 2 i+1/2,j+1,k+1/2 J 2 i+1/2,j,k+1/2 dy + J 1 i+1,j+1/2,k+1/2 J 1 i,j +1/2,k+1/2 dx + J 3 i+1/2,j+1/2,k+1 J 3 i+1/2,j+1/2,k dz ± = 0.
26 charge density ρ i,j,k = α Q α S i,j,k (x α,y α,z α ), S i,j,k (x α,y α,z α ) = S(X i x α,y j y α,z k z α ), = X i,y j,z k denote coordinates of the grid, (x α,y α,z α ) re the form-factor can be interpreted as a charge dens density decomposition: due to linearity of continuity equation, consider single superparticle vector W defined by J 1 i+1,j,k J 1 i,j,k = Qdx dt W 1 i,j,k, J 2 i,j+1,k J 2 i,j,k = Qdy dt W 2 i,j,k, J 3 i,j,k+1 J 3 i,j,k = Qdz dt W 3 i,j,k. continuity equation: W 1 + W 2 + W 3 = S(x + x,y + y,z+ z) S(x, y, z).
27 Here ( x, y, z) is 3-dimensional shift of the quasi-particle due to the motion. Shift of the quasi-particle generates eight functions S(x, y, z), S(x + x, y, z), S(x, y + y, z), S(x, y, z + z), S(x + x,y + y, z), S(x + x,y,z+ z), S(x, y + y,z+ z), S(x + x,y + y,z+ z). must increase proportionally. (2) The quasi-particle trajectory over time step is treated as a straight lin decompose any three-dimensional shift of the form-factor S(x, y, z) into three one-dimensional shifts: S(x + x,y + + y,z+ + z) + S(x, y, z) = S(x + x,y,z) S(x, y, z) We assume + S(x + x,y + y,z) S(x + x,y,z) + S(x + x,y + y,z+ z) S(x + x,y + y, z). (1) The vector (Wi,j,k 1,W2 i,j,k,w3 i,j,k ) is a linear combination of eight functions, Eq. (21). (2) W 1 + W 2 + W 3 = S(x + x,y + y,z+ z) S(x, y, z). (2) The sum of components of the vector is equal to the finite difference of form-factors, Eq. (20). (3) If one of the shifts x, y, z is zero, the corresponding component of W is also zero: x = 0 W 1 = 0, y = 0 W 2 = 0, z = 0 W 3 = 0. (4) If S(x, y, z) is symmetrical with respect to the permutation of (x, y), S(x, y, z) = S(y, x, z),and x = y, then W 1 = W 2. The same property is assumed for symmetries with respect to permutations of pairs (x, z) and (y, z). Only one linear combination of eight functions,obeys the conditions (1) (4):
28 Solution: W 1 = 1 3 S(x + x,y + y,z+ z) 1 3 S(x, y + y,z+ z) S(x + x,y,z + z) 1 6S(x, y, z + z) S(x + x,y + y,z) 1 6S(x, y + y,z) S(x + x,y,z) 1 3S(x, y, z), W 2 = 1 3 S(x + x,y + y,z+ z) 1 3 S(x + x,y,z + z) S(x, y + y,z+ z) 1 6 S(x, y, z + z) S(x + x,y + y,z) 1 6S(x + x,y,z) S(x, y + y,z) 1 3S(x, y, z), W 3 = 1 3 S(x + x,y + y,z+ z) 1 3 S(x + x,y + y,z) S(x, y + y,z+ z) 1 6 S(x, y + y,z) S(x + x,y,z + z) 1 6S(x + x,y,z) S(x, y, z + z) 1 3 S(x, y, z).
Guangye Chen, Luis Chacón,
JIFT workshop! Oct. 31, 2014 New Orleans, LA.! Guangye Chen, Luis Chacón, CoCoMANs team Los Alamos National Laboratory, Los Alamos, NM 87545, USA gchen@lanl.gov 1 Los Alamos National Laboratory Motivation
More informationSimulation Techniques for HED I: Particle-in-Cell Methods
2015 HED Summer School Simulation Techniques for HED I: Particle-in-Cell Methods Frank S. Tsung UCLA Special Thanks Viktor Decyk Peicheng Yu < - Slides Michael Meyers Thamine Dalichaouch < - Demos Outline
More informationLecture 4: The particle equations (1)
Lecture 4: The particle equations (1) Presenter: Mark Eric Dieckmann Department of Science and Technology (ITN), Linköping University, Sweden July 17, 2014 Overview We have previously discussed the leapfrog
More informationEXASCALE COMPUTING. Implementation of a 2D Electrostatic Particle in Cell algorithm in UniÞed Parallel C with dynamic load-balancing
ExaScience Lab Intel Labs Europe EXASCALE COMPUTING Implementation of a 2D Electrostatic Particle in Cell algorithm in UniÞed Parallel C with dynamic load-balancing B. Verleye P. Henry R. Wuyts G. Lapenta
More information1/30. Rigid Body Rotations. Dave Frank
. 1/3 Rigid Body Rotations Dave Frank A Point Particle and Fundamental Quantities z 2/3 m v ω r y x Angular Velocity v = dr dt = ω r Kinetic Energy K = 1 2 mv2 Momentum p = mv Rigid Bodies We treat a rigid
More informationParticle-in-Cell Codes for plasma-based particle acceleration
A. Pukhov Institute for Theoretical Physics I University of Dusseldorf, Germany Particle-in-Cell Codes for plasma-based particle acceleration Outline Relativistic plasmas, acceleration and the simulation
More informationNotes for Lecture 10
February 2, 26 Notes for Lecture Introduction to grid-based methods for solving Poisson and Laplace Equations Finite difference methods The basis for most grid-based finite difference methods is the Taylor
More informationPractical in Numerical Astronomy, SS 2012 LECTURE 9
Practical in Numerical Astronomy, SS 01 Elliptic partial differential equations. Poisson solvers. LECTURE 9 1. Gravity force and the equations of hydrodynamics. Poisson equation versus Poisson integral.
More informationImplicit-Explicit Time Integration of a High-Order Particle-in-Cell Method with Hyperbolic Divergence Cleaning.
Implicit-Explicit Time Integration of a High-Order Particle-in-Cell Method with Hyperbolic Divergence Cleaning. G.B. Jacobs a J.S. Hesthaven b a Department of Aerospace Engineering & Engineering Mechanics,
More informationTight-Focusing of Short Intense Laser Beams in Particle-in-Cell Simulations of Laser-Plasma Interaction
28/03/2017, CTU in Prague Tight-Focusing of Short Intense Laser Beams in Particle-in-Cell Simulations of Laser-Plasma Interaction Bc. Petr Valenta (petr.valenta@eli-beams.eu) Supervisors: doc. Ing. Ondrej
More informationPartititioned Methods for Multifield Problems
Partititioned Methods for Multifield Problems Joachim Rang, 22.7.215 22.7.215 Joachim Rang Partititioned Methods for Multifield Problems Seite 1 Non-matching matches usually the fluid and the structure
More informationTight-Focusing of Short Intense Laser Pulses in Particle-in-Cell Simulations of Laser-Plasma Interaction
16/05/2017, CTU in Prague Tight-Focusing of Short Intense Laser Pulses in Particle-in-Cell Simulations of Laser-Plasma Interaction Bc. Petr Valenta (petr.valenta@eli-beams.eu) Supervisors: doc. Ing. Ondrej
More informationComputational Relativistic Electrodynamics: New Algorithms, Parallel Software, and Applications to Accelerator Design
Computational Relativistic Electrodynamics: New Algorithms, Parallel Software, and Applications to Accelerator Design A Dissertation Presented by Kwang Min Yu to The Graduate School in Partial Fulfillment
More informationExponential Recursion for Multi-Scale Problems in Electromagnetics
Exponential Recursion for Multi-Scale Problems in Electromagnetics Matthew F. Causley Department of Mathematics Kettering University Computational Aspects of Time Dependent Electromagnetic Wave Problems
More informationAdaptive semi-lagrangian schemes for transport
for transport (how to predict accurate grids?) Martin Campos Pinto CNRS & University of Strasbourg, France joint work Albert Cohen (Paris 6), Michel Mehrenberger and Eric Sonnendrücker (Strasbourg) MAMCDP
More informationMA 201: Partial Differential Equations Lecture - 2
MA 201: Partial Differential Equations Lecture - 2 Linear First-Order PDEs For a PDE f(x,y,z,p,q) = 0, a solution of the type F(x,y,z,a,b) = 0 (1) which contains two arbitrary constants a and b is said
More informationPhysics 6303 Lecture 2 August 22, 2018
Physics 6303 Lecture 2 August 22, 2018 LAST TIME: Coordinate system construction, covariant and contravariant vector components, basics vector review, gradient, divergence, curl, and Laplacian operators
More informationOutline. 1 Why PIC Simulations Make Sense. 2 The Algorithm. 3 Examples Surface High Harmonics Generation. 4 Extensions Of The PIC Algorithm
PIC Simulations an Introduction GRK 1203 Meeting February 12-15, 2008, Oelde Outline 1 Simulations Make Sense 2 3 Surface High Harmonics Generation 4 Of PIC Plasma Physics Is Complex Experiment real thing
More informationUne décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion
Une décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion Anaïs Crestetto 1, Nicolas Crouseilles 2 et Mohammed Lemou 3 Rennes, 14ème Journée de l équipe
More information25. Chain Rule. Now, f is a function of t only. Expand by multiplication:
25. Chain Rule The Chain Rule is present in all differentiation. If z = f(x, y) represents a two-variable function, then it is plausible to consider the cases when x and y may be functions of other variable(s).
More informationEOS 250 Geophysical Fields and Fluxes Volume Integrals
EOS 25 Geophysical Fields and Fluxes olume Integrals February 17, 21 Overview These notes cover the following: Quantities defined per unit volume : volume integrals olume Integrals We have seen in the
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationMATHS 267 Answers to Stokes Practice Dr. Jones
MATH 267 Answers to tokes Practice Dr. Jones 1. Calculate the flux F d where is the hemisphere x2 + y 2 + z 2 1, z > and F (xz + e y2, yz, z 2 + 1). Note: the surface is open (doesn t include any of the
More informationMicro-macro methods for Boltzmann-BGK-like equations in the diffusion scaling
Micro-macro methods for Boltzmann-BGK-like equations in the diffusion scaling Anaïs Crestetto 1, Nicolas Crouseilles 2, Giacomo Dimarco 3 et Mohammed Lemou 4 Saint-Malo, 14 décembre 2017 1 Université de
More informationAn asymptotic-preserving micro-macro scheme for Vlasov-BGK-like equations in the diffusion scaling
An asymptotic-preserving micro-macro scheme for Vlasov-BGK-like equations in the diffusion scaling Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3 Saint-Malo 13 December 2016 1 Université
More informationPRACTICE PROBLEMS. Please let me know if you find any mistakes in the text so that i can fix them. 1. Mixed partial derivatives.
PRACTICE PROBLEMS Please let me know if you find any mistakes in the text so that i can fix them. 1.1. Let Show that f is C 1 and yet How is that possible? 1. Mixed partial derivatives f(x, y) = {xy x
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More informationE[X n ]= dn dt n M X(t). ). What is the mgf? Solution. Found this the other day in the Kernel matching exercise: 1 M X (t) =
Chapter 7 Generating functions Definition 7.. Let X be a random variable. The moment generating function is given by M X (t) =E[e tx ], provided that the expectation exists for t in some neighborhood of
More informationSummer College on Plasma Physics. 30 July - 24 August, The particle-in-cell simulation method: Concept and limitations
1856-30 2007 Summer College on Plasma Physics 30 July - 24 August, 2007 The particle-in-cell M. E. Dieckmann Institut fuer Theoretische Physik IV, Ruhr-Universitaet, Bochum, Germany The particle-in-cell
More informationLandau Damping Simulation Models
Landau Damping Simulation Models Hua-sheng XIE (u) huashengxie@gmail.com) Department of Physics, Institute for Fusion Theory and Simulation, Zhejiang University, Hangzhou 310027, P.R.China Oct. 9, 2013
More 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 informationMath 31CH - Spring Final Exam
Math 3H - Spring 24 - Final Exam Problem. The parabolic cylinder y = x 2 (aligned along the z-axis) is cut by the planes y =, z = and z = y. Find the volume of the solid thus obtained. Solution:We calculate
More informationModeling of planar plasma diode
International Letters of Chemistry, Physics and Astronomy Online: 2013-05-03 ISSN: 2299-3843, Vol. 13, pp 220-242 doi:10.18052/www.scipress.com/ilcpa.13.220 2013 SciPress Ltd., Switzerland Modeling of
More informationKinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles
Kinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3. Workshop Asymptotic-Preserving schemes, Porquerolles.
More informationExplicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems. Abstract
Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems Jianyuan Xiao, 1, 2 Hong Qin, 1, 3, Jian Liu, 1, 2 Yang He, 1, 2 Ruili Zhang, 1, 2 and Yajuan Sun 4 arxiv:1510.06972v2
More informationA charge preserving scheme for the numerical resolution of the Vlasov-Ampère equations.
A charge preserving scheme for the numerical resolution of the Vlasov-Ampère equations. Nicolas Crouseilles Thomas Respaud January, Abstract In this report, a charge preserving numerical resolution of
More informationSimulation of Coulomb Collisions in Plasma Accelerators for Space Applications
Simulation of Coulomb Collisions in Plasma Accelerators for Space Applications D. D Andrea 1, W.Maschek 1 and R. Schneider 2 Vienna, May 6 th 2009 1 Institut for Institute for Nuclear and Energy Technologies
More informationA fully implicit, exactly conserving algorithm for multidimensional particle-in-cell kinetic simulations
A fully implicit, exactly conserving algorithm for multidimensional particle-in-cell kinetic simulations L. Chacón Applied Mathematics and Plasma Physics Group Theoretical Division Los Alamos National
More informationKinetic Plasma Simulations. Anatoly Spitkovsky (Princeton)
Kinetic Plasma Simulations Anatoly Spitkovsky (Princeton) Contents Plasma physics on computers How PIC works Electrostatic codes Charge assignment and shape factors Discretization effects Electromagnetic
More informationLECTURE 10: THE PARALLEL TRANSPORT
LECTURE 10: THE PARALLEL TRANSPORT 1. The parallel transport We shall start with the geometric meaning of linear connections. Suppose M is a smooth manifold with a linear connection. Let γ : [a, b] M be
More informationKinetic Solvers with Adaptive Mesh in Phase Space for Low- Temperature Plasmas
Kinetic Solvers with Adaptive Mesh in Phase Space for Low- Temperature Plasmas Vladimir Kolobov, a,b,1 Robert Arslanbekov a and Dmitry Levko a a CFD Research Corporation, Huntsville, AL 35806, USA b The
More informationHere are some solutions to the sample problems assigned for Chapter 4. Solution: Consider a function of 3 (independent) variables. treating as real.
Lecture 11 Appendix B: Some sample problems from Boas Here are some solutions to the sample problems assigned for Chapter 4. 4.1: 3 Solution: Consider a function of 3 (independent) variables,, ln z u v
More informationMATH 0350 PRACTICE FINAL FALL 2017 SAMUEL S. WATSON. a c. b c.
MATH 35 PRACTICE FINAL FALL 17 SAMUEL S. WATSON Problem 1 Verify that if a and b are nonzero vectors, the vector c = a b + b a bisects the angle between a and b. The cosine of the angle between a and c
More informationES.182A Topic 45 Notes Jeremy Orloff
E.8A Topic 45 Notes Jeremy Orloff 45 More surface integrals; divergence theorem Note: Much of these notes are taken directly from the upplementary Notes V0 by Arthur Mattuck. 45. Closed urfaces A closed
More informationCapacitors. Lecture 10. Chapter 26. My Capacitance is limited. PHYS.1440 Lecture 10 Danylov. Department of Physics and Applied Physics
Lecture 10 Chapter 26 Capacitors My Capacitance is limited Course website: http://faculty.uml.edu/andriy_danylov/teaching/physicsii Today we are going to discuss: Chapter 26: Section 26.2 The Geometry
More informationLecture Notes Special: Using Neural Networks for Mean-Squared Estimation. So: How can we evaluate E(X Y )? What is a neural network? How is it useful?
Lecture Notes Special: Using Neural Networks for Mean-Squared Estimation The Story So Far... So: How can we evaluate E(X Y )? What is a neural network? How is it useful? EE 278: Using Neural Networks for
More informationLecture 5: Kinetic theory of fluids
Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationPart 8: Rigid Body Dynamics
Document that contains homework problems. Comment out the solutions when printing off for students. Part 8: Rigid Body Dynamics Problem 1. Inertia review Find the moment of inertia for a thin uniform rod
More informationKinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles
Kinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3. The 8th International Conference on Computational
More informationPractice Qualifying Exam Questions, Differentiable Manifolds, Fall, 2009.
Practice Qualifying Exam Questions, Differentiable Manifolds, Fall, 2009. Solutions (1) Let Γ be a discrete group acting on a manifold M. (a) Define what it means for Γ to act freely. Solution: Γ acts
More informationLECTURE 16 GAUSS QUADRATURE In general for Newton-Cotes (equispaced interpolation points/ data points/ integration points/ nodes).
CE 025 - Lecture 6 LECTURE 6 GAUSS QUADRATURE In general for ewton-cotes (equispaced interpolation points/ data points/ integration points/ nodes). x E x S fx dx hw' o f o + w' f + + w' f + E 84 f 0 f
More informationDevelopment of 2D particle-in-cell code to simulate high current, low energy beam in a beam transport system
PRAMANA c Indian Academy of Sciences Vol. 9, No. journal of October 7 physics pp. 551 5 Development of D particle-in-cell code to simulate high current, low energy beam in a beam transport system S C L
More informationMultiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015
Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction
More information18.02 Multivariable Calculus Fall 2007
MIT OpenourseWare http://ocw.mit.edu 18.02 Multivariable alculus Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 18.02 Lecture 30. Tue, Nov
More informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF MATHEMATICS ACADEMIC YEAR / EVEN SEMESTER QUESTION BANK
KINGS COLLEGE OF ENGINEERING MA5-NUMERICAL METHODS DEPARTMENT OF MATHEMATICS ACADEMIC YEAR 00-0 / EVEN SEMESTER QUESTION BANK SUBJECT NAME: NUMERICAL METHODS YEAR/SEM: II / IV UNIT - I SOLUTION OF EQUATIONS
More informationSimulations of Plasmas with Electrostatic PIC Models Using the Finite Element Method
Brazilian Journal of Physics, vol. 33, no., June, 003 411 Simulations of Plasmas with Electrostatic PIC Models Using the Finite Element Method A. C. J. Paes, N. M. Abe, V. A. Serrão, and A. Passaro Instituto
More information(1 + 2y)y = x. ( x. The right-hand side is a standard integral, so in the end we have the implicit solution. y(x) + y 2 (x) = x2 2 +C.
Midterm 1 33B-1 015 October 1 Find the exact solution of the initial value problem. Indicate the interval of existence. y = x, y( 1) = 0. 1 + y Solution. We observe that the equation is separable, and
More informationOutline: Phenomenology of Wall Bounded Turbulence: Minimal Models First Story: Large Re Neutrally-Stratified (NS) Channel Flow.
Victor S. L vov: Israeli-Italian March 2007 Meeting Phenomenology of Wall Bounded Turbulence: Minimal Models In collaboration with : Anna Pomyalov, Itamar Procaccia, Oleksii Rudenko (WIS), Said Elgobashi
More informationFinite Difference Methods for
CE 601: Numerical Methods Lecture 33 Finite Difference Methods for PDEs Course Coordinator: Course Coordinator: Dr. Suresh A. Kartha, Associate Professor, Department of Civil Engineering, IIT Guwahati.
More informationScientific Computing
2301678 Scientific Computing Chapter 2 Interpolation and Approximation Paisan Nakmahachalasint Paisan.N@chula.ac.th Chapter 2 Interpolation and Approximation p. 1/66 Contents 1. Polynomial interpolation
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More informationMULTIVARIABLE INTEGRATION
MULTIVARIABLE INTEGRATION (PLANE & CYLINDRICAL POLAR COORDINATES) PLANE POLAR COORDINATES Question 1 The finite region on the x-y plane satisfies 1 x + y 4, y 0. Find, in terms of π, the value of I. I
More informationThe Accelerator Hamiltonian in a Straight Coordinate System
Hamiltoninan Dynamics for Particle Accelerators, Lecture 2 The Accelerator Hamiltonian in a Straight Coordinate System Andy Wolski University of Liverpool, and the Cockcroft Institute, Daresbury, UK. Given
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationWilliam P. Thurston. The Geometry and Topology of Three-Manifolds
William P. Thurston The Geometry and Topology of Three-Manifolds Electronic version 1.1 - March 00 http://www.msri.org/publications/books/gt3m/ This is an electronic edition of the 1980 notes distributed
More informationMathematical Tripos Part IA Lent Term Example Sheet 1. Calculate its tangent vector dr/du at each point and hence find its total length.
Mathematical Tripos Part IA Lent Term 205 ector Calculus Prof B C Allanach Example Sheet Sketch the curve in the plane given parametrically by r(u) = ( x(u), y(u) ) = ( a cos 3 u, a sin 3 u ) with 0 u
More informationNumerical tensor methods and their applications
Numerical tensor methods and their applications 8 May 2013 All lectures 4 lectures, 2 May, 08:00-10:00: Introduction: ideas, matrix results, history. 7 May, 08:00-10:00: Novel tensor formats (TT, HT, QTT).
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationUnderstanding the pulsar magnetosphere through first-principle simulations
Understanding the pulsar magnetosphere through first-principle simulations Alexander Y. Chen In collaboration with: Andrei Beloborodov Rui Hu The Many Faces of Neutron Stars August 25, 2015 Pulsars: Rotating
More information21-256: Partial differentiation
21-256: Partial differentiation Clive Newstead, Thursday 5th June 2014 This is a summary of the important results about partial derivatives and the chain rule that you should know. Partial derivatives
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationHigh Order Semi-Lagrangian WENO scheme for Vlasov Equations
High Order WENO scheme for Equations Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Andrew Christlieb Supported by AFOSR. Computational Mathematics Seminar, UC Boulder
More informationSummary of time independent electrodynamics
hapter 10 Summary of time independent electrodynamics 10.1 Electrostatics Physical law oulomb s law charges as origin of electric field Superposition principle ector of the electric field E(x) in vacuum
More informationEEL 5544 Noise in Linear Systems Lecture 30. X (s) = E [ e sx] f X (x)e sx dx. Moments can be found from the Laplace transform as
L30-1 EEL 5544 Noise in Linear Systems Lecture 30 OTHER TRANSFORMS For a continuous, nonnegative RV X, the Laplace transform of X is X (s) = E [ e sx] = 0 f X (x)e sx dx. For a nonnegative RV, the Laplace
More informationIn this section, mathematical description of the motion of fluid elements moving in a flow field is
Jun. 05, 015 Chapter 6. Differential Analysis of Fluid Flow 6.1 Fluid Element Kinematics In this section, mathematical description of the motion of fluid elements moving in a flow field is given. A small
More informationMultiple Random Variables
Multiple Random Variables Joint Probability Density Let X and Y be two random variables. Their joint distribution function is F ( XY x, y) P X x Y y. F XY ( ) 1, < x
More informationJUST THE MATHS UNIT NUMBER PARTIAL DIFFERENTIATION 1 (Partial derivatives of the first order) A.J.Hobson
JUST THE MATHS UNIT NUMBER 14.1 PARTIAL DIFFERENTIATION 1 (Partial derivatives of the first order) by A.J.Hobson 14.1.1 Functions of several variables 14.1.2 The definition of a partial derivative 14.1.3
More informationBrad Gibson Centre for Astrophysics & Supercomputing Swinburne University
N-body Simulations: Tree and Mesh Approaches Brad Gibson Centre for Astrophysics & Supercomputing Swinburne University Demystifying the Jargon Trees Collisionless Boltzmann Equation Symplectic Time Integration
More informationSimulation of heating in space plasma
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research, 2015, 6(6):103-113 ISSN: 0976-8610 CODEN (USA): AASRFC Simulation of heating in space plasma 1* Umar S. A., 2 Ibrahim
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationPIC solvers for Intense Beams in Synchrotrons
PIC solvers for Intense Beams in Synchrotrons Oliver Boine-Franenheim TU Darmstadt and GSI Contents o Application and challenges for beam simulations with space charge o D electrostatic Particle-In-Cell
More informationMoments conservation in adaptive Vlasov solver
12 Moments conservation in adaptive Vlasov solver M. Gutnic a,c, M. Haefele b,c and E. Sonnendrücker a,c a IRMA, Université Louis Pasteur, Strasbourg, France. b LSIIT, Université Louis Pasteur, Strasbourg,
More informationThe Conjugate Gradient Method
The Conjugate Gradient Method Jason E. Hicken Aerospace Design Lab Department of Aeronautics & Astronautics Stanford University 14 July 2011 Lecture Objectives describe when CG can be used to solve Ax
More informationLecture 5: Moment generating functions
Lecture 5: Moment generating functions Definition 2.3.6. The moment generating function (mgf) of a random variable X is { x e tx f M X (t) = E(e tx X (x) if X has a pmf ) = etx f X (x)dx if X has a pdf
More informationAdaptive simulation of Vlasov equations in arbitrary dimension using interpolatory hierarchical bases
Adaptive simulation of Vlasov equations in arbitrary dimension using interpolatory hierarchical bases Erwan Deriaz Fusion Plasma Team Institut Jean Lamour (Nancy), CNRS / Université de Lorraine VLASIX
More informationParsek2D: An Implicit Parallel Particle-in-Cell Code
NUMERICAL MODELING OF SPACE PLASMA FLOWS: ASTRONUM-2008 ASP Conference Series, Vol. 406, c 2009 Nikolai V. Pogorelov, Edouard Audit, Phillip Colella, and Gary P. Zank, eds. Parsek2D: An Implicit Parallel
More informationOn the cosmic no-hair conjecture in the Einstein-Vlasov setting
On the cosmic no-hair conjecture in the Einstein-Vlasov setting KTH, Royal Institute of Technology, Stockholm Recent Advances in Mathematical General Relativity Institut Henri Poincaré, Paris September
More informationanalysis for transport equations and applications
Multi-scale analysis for transport equations and applications Mihaï BOSTAN, Aurélie FINOT University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Numerical methods for kinetic equations Strasbourg
More informationComputational Physics
Interpolation, Extrapolation & Polynomial Approximation Lectures based on course notes by Pablo Laguna and Kostas Kokkotas revamped by Deirdre Shoemaker Spring 2014 Introduction In many cases, a function
More informationSDS 321: Introduction to Probability and Statistics
SDS 321: Introduction to Probability and Statistics Lecture 17: Continuous random variables: conditional PDF Purnamrita Sarkar Department of Statistics and Data Science The University of Texas at Austin
More informationVerification and Convergence Properties of a Particle-In-Cell Code
Verification and Convergence Properties of a Particle-In-Cell Code B. J. Winjum 1, J. J. Tu 2, S. S. Su 3, V. K. Decyk 1, W. B. Mori 1 1 UCLA 2 Fudan U., China 3 U. Michigan 2014 Anomalous Absorption Conference
More informationCH.3. COMPATIBILITY EQUATIONS. Multimedia Course on Continuum Mechanics
CH.3. COMPATIBILITY EQUATIONS Multimedia Course on Continuum Mechanics Overview Introduction Lecture 1 Compatibility Conditions Lecture Compatibility Equations of a Potential Vector Field Lecture 3 Compatibility
More informationChapter 2. Deriving the Vlasov Equation From the Klimontovich Equation 19. Deriving the Vlasov Equation From the Klimontovich Equation
Chapter 2. Deriving the Vlasov Equation From the Klimontovich Equation 19 Chapter 2. Deriving the Vlasov Equation From the Klimontovich Equation Topics or concepts to learn in Chapter 2: 1. The microscopic
More informationPhysical models for plasmas II
Physical models for plasmas II Dr. L. Conde Dr. José M. Donoso Departamento de Física Aplicada. E.T.S. Ingenieros Aeronáuticos Universidad Politécnica de Madrid Physical models,... Plasma Kinetic Theory
More informationPhysics 325: General Relativity Spring Final Review Problem Set
Physics 325: General Relativity Spring 2012 Final Review Problem Set Date: Friday 4 May 2012 Instructions: This is the third of three review problem sets in Physics 325. It will count for twice as much
More information