Outline. 1 Why PIC Simulations Make Sense. 2 The Algorithm. 3 Examples Surface High Harmonics Generation. 4 Extensions Of The PIC Algorithm
|
|
- Cornelius Gibson
- 5 years ago
- Views:
Transcription
1 PIC Simulations an Introduction GRK 1203 Meeting February 12-15, 2008, Oelde
2 Outline 1 Simulations Make Sense 2 3 Surface High Harmonics Generation 4 Of PIC
3 Plasma Physics Is Complex Experiment real thing we are actually interested in Plasma diagnostics mostly very indirect Modication of certain parameters may require big eort Analytical ory Yields a profound physical understanding Relies on extensive approximations and assumptions Real experiments often too complex to be fully described by analytical theory Bridging the gap Simulations provide a bridge between experiment and theory
4 How To Simulate? Full kinetic description of a collisionless plasma: Vlasov equation Hydro-Codes p t f + m 1 + xf + F p 2 m pf = 0 Further assume thermal equilibrium Derive momenta equations, easier to solve But: Many interesting plasma phenomena are essentially non-thermal So we stick to the full kinetic description!
5 How To Simulate? Full kinetic description of a collisionless plasma: Vlasov equation Hydro-Codes p t f + m 1 + xf + F p 2 m pf = 0 Further assume thermal equilibrium Derive momenta equations, easier to solve But: Many interesting plasma phenomena are essentially non-thermal So we stick to the full kinetic description!
6 Vlasov PIC It is a challenge to represent f (x, p,t) numerically Direct Vlasov Direct sampling on a 6D Eulerian mesh Straightforward Naturally produces smooth results ly extremely intensive PIC Sampling using Macroparticles: f (x, p) = p W n S(x x p, p p p ) ly eective Intuitive picture of particle clouds
7 Basic Equations Maxwell E = 4πρ B = 0 E = B = 1 c Ḃ c j + 1 c Ė EOM ẋ p = pp m 1+p 2 p ṗ p = F p Fluid elements have to obey equation of motion EOM for macroparticles substitutes Vlasov Macroparticles single electrons/ions How to solve these equations numerically?
8 Basic Equations Maxwell E = 4πρ B = 0 E = B = 1 c Ḃ c j + 1 c Ė EOM ẋ p = pp m 1+p 2 p ṗ p = F p Fluid elements have to obey equation of motion EOM for macroparticles substitutes Vlasov Macroparticles single electrons/ions How to solve these equations numerically?
9 Particles In Cells Fields dened at grid points (cells) E c, B c, ρ c, j c Particles move freely in space x p, p p (taken from Birdsall, Langdon - Plasma Physics via Computer Simulation)
10 Cycle Integrate EOM F p for each particle p Interpolate to particle positions x p, p p for each particle p E c, B c for each cell c Interpolate to cell positions ρ c, j c for each cell c This computational cycle represents one time step Dierent realisations possible Integrate Maxwell Now let us take a walk through the cycle of our code VLPL!
11 Cycle Integrate EOM F p for each particle p Interpolate to particle positions x p, p p for each particle p E c, B c for each cell c Interpolate to cell positions ρ c, j c for each cell c This computational cycle represents one time step Dierent realisations possible Integrate Maxwell Now let us take a walk through the cycle of our code VLPL!
12 Integration Of EOM using the leap-frog scheme (taken from Birdsall, Langdon - Plasma Physics via Computer Simulation) Velocity (respectively momentum p) and position x are dened at staggered timesteps d p = F dt d x = p dt 1+p 2 pnew pold = t Fold xnew xold = p new t 1+p 2 new
13 Accuracy Of Leap-Frog Scheme Test the scheme by applying it to a non-relativistic harmonic potential F = mω 2 0 x Analysis yields: Amplitude is reproduced without mistake Phase error in third order: φ = ω 0 t (ω 0 t) 3 + What about higher order schemes? Big number of particles in simulation Higher order schemes use more storage usually not appropriate for PIC Conclusion simple leap-frog scheme is ideal for PIC
14 Lorentz Force Boris scheme Problem Magnetic force depends on p Leap frog scheme: p not given at the same time as F pnew pold t = q ( E + 1 c pnew + pold 2γ wherein γ is the time-centred relativistic energy B This can indeed be solved for pnew (Boris, 1970) In principle: apply half of the E-eld, then B-eld, then rest of the E-eld )
15 Maxwell Equations I the Yee mesh Use Ampère law to propagate E-eld: Ė = c B 4πj E n+1 E n = ( t) (c ˆ ) B n+1/2 4πj n+1/2, analogue Faraday for B-eld ˆ nite dierence version of curl operator staggered spatial lattice Yee mesh j i 2D geometry B y E x E z E y B z Bx k i B E z x E y B y B x E z 3D geometry j
16 Maxwell Equations II numerical dispersion free scheme k y h y 0 π (a) π k x h x k y h y 0 π (b) k x h x π taken from A Pukhov Straightforward nite dierence translation ( x f ( ) f x+ x/2 f x x/2 / x) of the curl operator ) Stable for (c t) 2 1/ (( x) 2 + ( y) 2 + ( z) 2 Numerical dispersion for short wavelengths (see Fig. (a)) Many problems in laser plasma physics: non-isotropic VLPL also implements a more sophisticated scheme Stable for c t x No articial dispersion in x-direction (see Fig. (b))
17 Maxwell Equations III divergence equations Substitute divergence equations ( E = 4πρ, B = 0 ) by the charge continuity equation ρ + j = 0 (1) = If divergence equations satised initially, they will automatically remain In VLPL, (1) is guaranteed by the proper (ρ, j) interpolation scheme
18 Interpolating ρ c Imagine macroparticles as massive cubes with uniform charge distribution edge length = grid step Distribute charge to the cells appropriately linear interpolation scheme Figure shows example in 1D taken from Birdsall / Langdon (a) shows the macroparticle inside the grid (b) shows the charge contained in grid cell x i as macroparticle passes by
19 Interpolating j c And Field Quantities Important for rigorous charge and energy conservation! Dene j c half timesteps shifted in relation to ρ c and at the cell boundaries Now follow the particle trajectories t j c = da dt vw p S p = da c 0 to perfectly fulll the continuity equation c xnew xold d xw p S p For optimal energy conservation, E-eld interpolation must be consistent with the current weighting
20 Input Parameters for a PIC run Initial conditions Particles Fields plasma or electron beams can be initialised inside the simulation box laser pulses may be initialised inside the box various shapes possible and readily implemented Boundary conditions Can be reecting, periodic, or absorbing both for particles and elds Particles or laser pulses may also enter the simulation box
21 Parallelisation For simple yet eective parallelisation, the simulation box is split in partitions A special boundary condition allows for data exchange between the partitions In the VLPL code, parallelisation is realised using MPI (message passing interface)
22 Diagnostics Because PIC produces huge amounts of data, it is usually not possible to store everything Before the simulation one must think of what to save: Time histories E.g. record the E- and B-elds at certain points in space Later one may Fourier-transform to obtain spectra Snapshots Movies Save all data at a certain time step Take a cross section and save periodically
23 Surface High Harmonics Generation Relativistic laser reected at an overdense plasma surface Reected light spectrum contains multiples of the laser frequency (harmonics) Power law spectrum found by Baeva, Gordienko, Pukhov Results conrmed rst by PIC simulation, then by experiment PIC simulation provides vivid insight to what's happening
24 Surface High Harmonics Generation Movie
25 Surface High Harmonics Generation Movie
26 Surface High Harmonics Generation Movie
27 Surface High Harmonics Generation Movie
28 Surface High Harmonics Generation Movie
29 Surface High Harmonics Generation Movie
30 Surface High Harmonics Generation Movie
31 Surface High Harmonics Generation Movie
32 Surface High Harmonics Generation Movie
33 Surface High Harmonics Generation Movie
34 Surface High Harmonics Generation Movie
35 Surface High Harmonics Generation Movie
36 Surface High Harmonics Generation Movie
37 Surface High Harmonics Generation Movie
38 Surface High Harmonics Generation Movie
39 Surface High Harmonics Generation Movie
40 Surface High Harmonics Generation Movie
41 Surface High Harmonics Generation Movie
42 Surface High Harmonics Generation Movie
43 Surface High Harmonics Generation Movie
44 Pictures and Simulation by Anupam Karmakar Electron beam propagating through a plasma provokes instability PIC brings not only marvellous pictures, but also deep insight into the physics E.g., collisions and temperature can be switched on or o at will to study their eects separately
45 Picture by A Pukhov Ultra-short, ultra-relativistic laser pulse propagating through tenuous plasma In the wakeeld a bubble emerges, wherein electrons are accelerated to very high energies Predicted by PIC, conrmed by experiments
46 Of PIC small modications Ionisation Using a Monte Carlo method, implemented for VLPL by A Karmakar Allows for simulation of laser interaction with high Z ions Collisions Using a Monte Carlo method, implemented for VLPL by A Karmakar Interesting e.g. for instability Lorentz transformation of initial conditions and results Changing the frame of reference allows to simulate some cases more eectively
47 Of PIC big modications Hybrid codes A part of the plasma is represented hydrodynamically, the non-thermal part by PIC Allows for simulation of non-thermal high density plasmas without resolving the plasma frequency Implemented in 1D by A Karmakar (One dimensional electromagnetic relativistic PIC-hydrodynamic hybrid simulation code H-VLPL (Hybrid Virtual Laser Plasma Lab), in press, Computer Physics Communication), 3D version planned Other interesting ideas: Quasi-static PIC codes Adaptive codes...
48 PIC simulation is an indispensable tool for modern plasma physics. algorithm is a concatenation of basically simple, traceable steps. Clever extensions allow to simulate even cases that are not covered by classical PIC.
49 For Further Reading Appendix For Further Reading C K Birdsall and A B Langdon. Plasma Physics via Computer Simulation. IOP Publishing, A Pukhov. Three-Dimensional Particle-in-Cell Simulations of Relativistic Laser-Plasma Interactions. Lecture Notes. Osaka University, 1999.
Particle-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 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 informationHybrid Simulation Method ISSS-10 Banff 2011
Hybrid Simulation Method ISSS-10 Banff 2011 David Burgess Astronomy Unit Queen Mary University of London With thanks to Dietmar Krauss-Varban Space Plasmas: From Sun to Earth Space Plasma Plasma is (mostly)
More 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 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 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 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 informationStudy of Laser Plasma Interactions Using an Eulerian Vlasov Code
PSFC/JA-04-6 Study of Laser Plasma Interactions Using an Eulerian Vlasov Code D. J. Strozzi, M. M. Shoucri*, and A. Bers March 2004 Plasma Science and Fusion Center Massachusetts Institute of Technology
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 informationBeam-Beam Simulations for e + e Colliders
Beam-Beam Simulations for e + e Colliders J. Rogers, Cornell University Motivation for simulations from observed effects Methods algorithms, simulation codes, and some results. Conclusions Brown Bag Seminar,
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More 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 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 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 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 informationMacroscopic plasma description
Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion
More informationSpace Plasma Physics Thomas Wiegelmann, 2012
Space Plasma Physics Thomas Wiegelmann, 2012 1. Basic Plasma Physics concepts 2. Overview about solar system plasmas Plasma Models 3. Single particle motion, Test particle model 4. Statistic description
More informationThree-dimensional relativistic particle-in-cell hybrid code based on an exponential integrator
Three-dimensional relativistic particle-in-cell hybrid code based on an exponential integrator T. Tückmantel, A. Pukhov Institut für Theoretische Physik I Heinrich-Heine Universität Düsseldorf 45, Düsseldorf,
More informationParticle-in-cell (PIC) simulation output for the temporal evolution of magnetic fields.
Type of file: pdf Title of file for HTML: Supplementary Information Description: Supplementary Figures and Supplementary Discussion. Type of file: MOV Title of file for HTML: Supplementary Movie 1 Description:
More informationReduced MHD. Nick Murphy. Harvard-Smithsonian Center for Astrophysics. Astronomy 253: Plasma Astrophysics. February 19, 2014
Reduced MHD Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 19, 2014 These lecture notes are largely based on Lectures in Magnetohydrodynamics by Dalton
More informationElectron Acceleration in Laser Plasma
Electron Acceleration in Laser Plasma Vojt ech Horn y IPP AV CR 4th December 2014, Praha Vojt ech Horn y (IPP AV CR) Electron Acceleration in Laser Plasma 4th December 2014, Praha 1 / 34 List of Contents
More informationLecture 13 Notes, Electromagnetic Theory I Dr. Christopher S. Baird University of Massachusetts Lowell
Lecture 13 Notes, Electromagnetic Theory I Dr. Christopher S. Baird University of Massachusetts Lowell 1. Static Equations and Faraday's Law - The two fundamental equations of electrostatics are shown
More informationApplicability of atomic collisional ionization cross sections in plasma environment
Applicability of atomic collisional ionization cross sections in plasma environment Viktoriia Isaenko,Voronezh State University, Russia CFEL Theory Division September 6, 2018 Abstract Secondary ionization
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationWeibel instability and filamentary structures of a relativistic electron beam in plasma
Mitglied der Helmholtz-Gemeinschaft 7 th Direct Drive and Fast Ignition Workshop Prague, 3-6 May, 009 Weibel instability and filamentary structures of a relativistic electron beam in plasma Anupam Karmakar
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 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 informationCurrent sheath formation in the plasma focus
Plasma Science and Applications (ICPSA 2013) International Journal of Modern Physics: Conference Series Vol. 32 (2014) 1460321 (8 pages) The Author DOI: 10.1142/S2010194514603214 Current sheath formation
More informationSummer College on Plasma Physics. 30 July - 24 August, The forming of a relativistic partially electromagnetic planar plasma shock
1856-31 2007 Summer College on Plasma Physics 30 July - 24 August, 2007 The forming of a M. E. Dieckmann Institut fuer Theoretische Physik IV, Ruhr-Universitaet, Bochum, Germany The forming of a The forming
More informationLaser trigged proton acceleration from ultrathin foil
Laser trigged proton acceleration from ultrathin foil A.V. Brantov 1, V. Yu. Bychenkov 1, D. V. Romanov 2, A. Maksimchuk 3 1 P. N. Lebedev Physics Institute RAS, Moscow 119991, Russia 2 All-Russia Research
More informationClassical Scattering
Classical Scattering Daniele Colosi Mathematical Physics Seminar Daniele Colosi (IMATE) Classical Scattering 27.03.09 1 / 38 Contents 1 Generalities 2 Classical particle scattering Scattering cross sections
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 informationTitle Surface Interaction under Oblique Intense. Author(s) Ruhl, H.; Sentoku, Y.; Mima, K.; Ta. Citation Physical Review Letters. 82(4) P.
Title Collimated Electron Jets by Surface Interaction under Oblique Intense I Author(s) Ruhl, H.; Sentoku, Y.; Mima, K.; Ta Citation Physical Review Letters. 82(4) P.74 Issue 1999-01-25 Date Text Version
More informationPHYS 1444 Section 004 Lecture #22
PHYS 1444 Section 004 Lecture #22 Monday, April 23, 2012 Dr. Extension of Ampere s Law Gauss Law of Magnetism Maxwell s Equations Production of Electromagnetic Waves Today s homework is #13, due 10pm,
More informationNOVEL SOLITONS? Classical and quantum solitons Solitons in bers 2D & 3D solitons in parametric ampliers Optical switching with parametric solitons Gap
SOLITONS: NOVEL AND ATOM LASERS PARAMPS Drummond Peter He Hao Kheruntsyan Karen of Queensland University September 18, 1998 Drummond, He, Kheruntsyan { Novel solitons: paramps and atom lasers NOVEL SOLITONS?
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 informationGuangye Chen, Luis Chacón,
JIFT workshop! Oct. 31, 2014 New Orleans, LA.! Guangye Chen, Luis Chacón, CoCoMANs team Los Alamos National Laboratory, Los Alamos, NM 87545, USA gchen@lanl.gov 1 Los Alamos National Laboratory Motivation
More 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 informationSource Free Surface x
Finite-dierence time-domain model for elastic waves in the ground Christoph T. Schroeder and Waymond R. Scott, Jr. School of Electrical and Computer Engineering Georgia Institute of Technology Atlanta,
More informationPHYSICS Computational Plasma Physics
PHYSICS 78 - Computational Plasma Physics INSTRUCTOR Dr. Earl Scime (escime@wvu.edu) 93-34, ext. 1437 Office hours: MW :30 3:30 and whenever door is open Rm 18 & 05 Hodges Hall Class: MWF 1:30-:0 Rm 334
More informationlimit of the time-step decreases as more resolutions are added requires the use of an eective multitime-stepping algorithm, that will maintain the req
Invited to the Session: "Wavelet and TLM Modeling Techniques" organized by Dr. W. J. R. Hoefer in ACES 2000 COMPUTATIONAL OPTIMIZATION OF MRTD HAAR-BASED ADAPTIVE SCHEMES USED FOR THE DESIGN OF RF PACKAGING
More informationCold plasma waves. Waves in non-magnetized plasma Cold plasma dispersion equation Cold plasma wave modes
Cold plasma waves Waves in non-magnetized plasma Cold plasma dispersion equation Cold plasma wave modes EM wave propagation through and interaction with plasmas belong to central issues of plasma physics.
More informationIntroduction Introduction
1 Introduction This book is an introduction to the theory of charged particle acceleration. It has two primary roles: 1.A unified, programmed summary of the principles underlying all charged particle accelerators.
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 informationComputational Methods in Plasma Physics
Computational Methods in Plasma Physics Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Purpose of Talk Describe use of numerical methods to solve simple problem in plasma
More 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 informationKinetic Plasma Simulations in Astrophysics. Lorenzo Sironi
Kinetic Plasma Simulations in Astrophysics Lorenzo Sironi Outline Plasma physics in astrophysics. The Vlasov-Maxwell system. Fully-kinetic particle-in-cell codes. 1. Electrostatic codes. 2. Electromagnetic
More informationESSENTIAL QUANTUM PHYSICS PETER LANDSHOFF. University of Cambridge ALLEN METHERELL. University of Central Florida GARETH REES. University of Cambridge
ESSENTIAL QUANTUM PHYSICS PETER LANDSHOFF University of Cambridge ALLEN METHERELL University of Central Florida GARETH REES University of Cambridge CAMBRIDGE UNIVERSITY PRESS Constants of quantum physics
More informationMethod of Perturbative-PIC Simulation for CSR Effect
Method of Perturbative-PIC Simulation for CSR Effect Jack J. Shi Department of Physics & Astronomy, University of Kansas OUTLINE Why Do We Want to Do This? Perturbation Expansion of the Retardation of
More informationAnisotropic electron distribution functions and the transition between the Weibel and the whistler instabilities
Anisotropic electron distribution functions and the transition between the Weibel and the whistler instabilities F. Pegoraro, L. Palodhi, F. Califano 5 th INTERNATIONAL CONFERENCE ON THE FRONTIERS OF PLASMA
More informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More informationIntroduction to electromagnetic theory
Chapter 1 Introduction to electromagnetic theory 1.1 Introduction Electromagnetism is a fundamental physical phenomena that is basic to many areas science and technology. This phenomenon is due to the
More informationAntennas and Propagation. Chapter 2: Basic Electromagnetic Analysis
Antennas and Propagation : Basic Electromagnetic Analysis Outline Vector Potentials, Wave Equation Far-field Radiation Duality/Reciprocity Transmission Lines Antennas and Propagation Slide 2 Antenna Theory
More informationIntro. Lecture 04: Numerical Methods for Particle and Distribution Methods: Introduction to the PIC Method *
Intro. Lecture 04: Numerical Methods for Particle and Distribution Methods: Introduction to the PIC Method * Prof. Steven M. Lund Physics and Astronomy Department Facility for Rare Isotope Beams (FRIB)
More informationarxiv: v1 [physics.comp-ph] 4 Feb 2015
arxiv:1502.01376v1 [physics.comp-ph] 4 Feb 2015 Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-fft Maxwell solver and a local charge conserving current deposit
More informationHybrid modeling of plasmas
Hybrid modeling of plasmas Mats Holmström Swedish Institute of Space Physics Kiruna, Sweden ENUMATH 2009 Uppsala, June 30 matsh@irf.se www.irf.se/~matsh/ Outline Background and motivation Space plasma
More informationCitation for published version (APA): Martinus, G. H. (1998). Proton-proton bremsstrahlung in a relativistic covariant model s.n.
University of Groningen Proton-proton bremsstrahlung in a relativistic covariant model Martinus, Gerard Henk IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you
More informationChapter 6. Maxwell s Equations, Macroscopic Electromagnetism, Conservation Laws
PHYS 532 Lecture 7 Page 1 Chapter 6. Maxwell s Equations, Macroscopic Electromagnetism, Conservation Laws 6.1 Maxwell Displacement Current and Maxwell Equations Differential equations for calculating fields
More informationPhysics 3312 Lecture 9 February 13, LAST TIME: Finished mirrors and aberrations, more on plane waves
Physics 331 Lecture 9 February 13, 019 LAST TIME: Finished mirrors and aberrations, more on plane waves Recall, Represents a plane wave having a propagation vector k that propagates in any direction with
More informationMaxwell s equations and EM waves. From previous Lecture Time dependent fields and Faraday s Law
Maxwell s equations and EM waves This Lecture More on Motional EMF and Faraday s law Displacement currents Maxwell s equations EM Waves From previous Lecture Time dependent fields and Faraday s Law 1 Radar
More informationMonte Carlo simulations of harmonic and anharmonic oscillators in discrete Euclidean time
Monte Carlo simulations of harmonic and anharmonic oscillators in discrete Euclidean time DESY Summer Student Programme, 214 Ronnie Rodgers University of Oxford, United Kingdom Laura Raes University of
More informationPage 1 of 9. Curriculum Map: Physics/Lab Course: Physics Sub-topic: Physics. Unit: Language of Physics Timeline: 2 Weeks Unit Description:
Curriculum Map: Physics/Lab Course: Physics Sub-topic: Physics Grade(s): 10 to 12 Course Course Textbooks, Workbooks, Materials Citations: Through theory and laboratory practices, the student will learn
More informationWhile the Gauss law forms for the static electric and steady magnetic field equations
Unit 2 Time-Varying Fields and Maxwell s Equations While the Gauss law forms for the static electric and steady magnetic field equations remain essentially unchanged for the case of time-varying fields,
More informationPHYSICS CURRICULUM. Unit 1: Measurement and Mathematics
Chariho Regional School District - Science Curriculum September, 2016 PHYSICS CURRICULUM Unit 1: Measurement and Mathematics OVERVIEW Summary Mathematics is an essential tool of physics. This unit will
More informationLinear and non-linear evolution of the gyroresonance instability in Cosmic Rays
Linear and non-linear evolution of the gyroresonance instability in Cosmic Rays DESY Summer Student Programme, 2016 Olga Lebiga Taras Shevchenko National University of Kyiv, Ukraine Supervisors Reinaldo
More informationBenchmarks in Computational Plasma Physics
Benchmarks in Computational Plasma Physics P. Londrillo INAF, Bologna, Italie S. Landi Università di Firenze, Italie What you compute when you do computations of the Vlasov equation? Overview A short review
More informationPrancing Through Quantum Fields
November 23, 2009 1 Introduction Disclaimer Review of Quantum Mechanics 2 Quantum Theory Of... Fields? Basic Philosophy 3 Field Quantization Classical Fields Field Quantization 4 Intuitive Field Theory
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 informationLecture 2. Introduction to plasma physics. Dr. Ashutosh Sharma
Preparation of the concerned sectors for educational and R&D activities related to the Hungarian ELI project Ion acceleration in plasmas Lecture 2. Introduction to plasma physics Dr. Ashutosh Sharma Zoltán
More informationIntegrated Modeling of Fast Ignition Experiments
Integrated Modeling of Fast Ignition Experiments Presented to: 9th International Fast Ignition Workshop Cambridge, MA November 3-5, 2006 R. P. J. Town AX-Division Lawrence Livermore National Laboratory
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 informationLectures on A coherent dual vector field theory for gravitation.
Lectures on coherent dual vector field theory for gravitation. The purpose of these lectures is to get more familiarized with gyrotation concepts and with its applications. Lecture : a word on the Maxwell
More informationJacob Trier Frederiksen, Niels Bohr Institute, Copenhagen, DK
Jacob Trier Frederiksen, Niels Bohr Institute, Copenhagen, DK Computational Astrophysics group Niels Bohr Institute Computational power History Present state at the Niels Bohr Institute Prospects Plasma
More informationElectrodynamics of Magnetized Rotators Anatoly Spitkovsky,, UC Berkeley
Electrodynamics of Magnetized Rotators Anatoly Spitkovsky,, UC Berkeley Magnetized rotators are ubiquitous: pulsars, AGN, GRBs (?) Rotation very efficient at long-term energy storage Extraction of rotational
More informationThe Larmor Formula (Chapters 18-19)
2017-02-28 Dispersive Media, Lecture 12 - Thomas Johnson 1 The Larmor Formula (Chapters 18-19) T. Johnson Outline Brief repetition of emission formula The emission from a single free particle - the Larmor
More informationSimulating experiments for ultra-intense laser-vacuum interaction
Simulating experiments for ultra-intense laser-vacuum interaction Nina Elkina LMU München, Germany March 11, 2011 Simulating experiments for ultra-intense laser-vacuum interaction March 11, 2011 1 / 23
More informationElectromagnetic Waves Retarded potentials 2. Energy and the Poynting vector 3. Wave equations for E and B 4. Plane EM waves in free space
Electromagnetic Waves 1 1. Retarded potentials 2. Energy and the Poynting vector 3. Wave equations for E and B 4. Plane EM waves in free space 1 Retarded Potentials For volume charge & current = 1 4πε
More informationCondensed Matter Physics Prof. G. Rangarajan Department of Physics Indian Institute of Technology, Madras
Condensed Matter Physics Prof. G. Rangarajan Department of Physics Indian Institute of Technology, Madras Lecture - 10 The Free Electron Theory of Metals - Electrical Conductivity (Refer Slide Time: 00:20)
More informationCHAPTER 32: ELECTROMAGNETIC WAVES
CHAPTER 32: ELECTROMAGNETIC WAVES For those of you who are interested, below are the differential, or point, form of the four Maxwell s equations we studied this semester. The version of Maxwell s equations
More informationResearch on Structure Formation of Plasmas Dominated by Multiple Hierarchical Dynamics
Chapter 4 Epoch Making Simulation Research on Structure Formation of Plasmas Dominated by Multiple Hierarchical Dynamics Group Representative Yasuaki Kishimoto Naka Fusion Research Establishment, Japan
More informationParticle in Cell method
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
More informationNumerical simulation methods for laser target interactions
Numerical simulation methods for laser target interactions Jiří Limpouch Czech Technical University in Prague Faculty of Nuclear Sciences & Physical Engineering Brehova 7, 115 19 Praha 1, Czechia jiri.limpouch@fjfi.cvut.cz
More informationLectures on basic plasma physics: Kinetic approach
Lectures on basic plasma physics: Kinetic approach Department of applied physics, Aalto University April 30, 2014 Motivation Layout 1 Motivation 2 Boltzmann equation (a nasty bastard) 3 Vlasov equation
More informationIntroduction to the School
Lucio Crivellari Instituto de Astrofísica de Canarias D.pto de Astrofísica, Universidad de La Laguna & INAF Osservatorio Astronomico di Trieste (Italy) Introduction to the School 10/11/17 1 Setting the
More informationPIC-MCC/Fluid Hybrid Model for Low Pressure Capacitively Coupled O 2 Plasma
PIC-MCC/Fluid Hybrid Model for Low Pressure Capacitively Coupled O 2 Plasma Kallol Bera a, Shahid Rauf a and Ken Collins a a Applied Materials, Inc. 974 E. Arques Ave., M/S 81517, Sunnyvale, CA 9485, USA
More informationFundamentals of wave kinetic theory
Fundamentals of wave kinetic theory Introduction to the subject Perturbation theory of electrostatic fluctuations Landau damping - mathematics Physics of Landau damping Unmagnetized plasma waves The plasma
More informationPreparation of the concerned sectors for educational and R&D activities related to the Hungarian ELI project
Preparation of the concerned sectors for educational and R&D activities related to the Hungarian ELI project Ion acceleration in plasmas Lecture 10. Collisionless shock wave acceleration in plasmas pas
More informationParticle in cell simulations
Particle in cell simulations Part III: Boundary conditions and parallelization Benoît Cerutti IPAG, CNRS, Université Grenoble Alpes, Grenoble, France. 1 Astrosim, Lyon, June 26 July 7, 2017. Plan of the
More informationS1: Particle Equations of Motion S1A: Introduction: The Lorentz Force Equation
S1: Particle Equations of Motion S1A: Introduction: The Lorentz Force Equation The Lorentz force equation of a charged particle is given by (MKS Units):... particle mass, charge... particle coordinate...
More informationブラックホール磁気圏での 磁気リコネクションの数値計算 熊本大学 小出眞路 RKKコンピュー 森野了悟 ターサービス(株) BHmag2012,名古屋大学,
RKK ( ) BHmag2012,, 2012.2.29 Outline Motivation and basis: Magnetic reconnection around astrophysical black holes Standard equations of resistive GRMHD Test calculations of resistive GRMHD A simulation
More informationDerivation of the General Propagation Equation
Derivation of the General Propagation Equation Phys 477/577: Ultrafast and Nonlinear Optics, F. Ö. Ilday, Bilkent University February 25, 26 1 1 Derivation of the Wave Equation from Maxwell s Equations
More informationFundamental Interactions (Forces) of Nature
Chapter 14 Fundamental Interactions (Forces) of Nature Interaction Gauge Boson Gauge Boson Mass Interaction Range (Force carrier) Strong Gluon 0 short-range (a few fm) Weak W ±, Z M W = 80.4 GeV/c 2 short-range
More informationPart VIII. Interaction with Solids
I with Part VIII I with Solids 214 / 273 vs. long pulse is I with Traditional i physics (ICF ns lasers): heating and creation of long scale-length plasmas Laser reflected at critical density surface Fast
More informationxkcd.com It IS about physics. It ALL is.
xkcd.com It IS about physics. It ALL is. Introduction to Space Plasmas The Plasma State What is a plasma? Basic plasma properties: Qualitative & Quantitative Examples of plasmas Single particle motion
More informationIntro. Lecture 04: Numerical Methods for Particle and Distribution Methods: Introduction to the PIC Method *
Intro. Lecture 04: Numerical Methods for Particle and Distribution Methods: Introduction to the PIC Method * Prof. Steven M. Lund Physics and Astronomy Department Facility for Rare Isotope Beams (FRIB)
More informationBaccalieu Collegiate. Physics Course Outline
Baccalieu Collegiate Physics 2204 Course Outline Course Content: Unit 1: Kinematics Motion is a common theme in our everyday lives: birds fly, babies crawl, and we, ourselves, seem to be in a constant
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 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 informationThermal and Ionization Aspects of Flows from Hot Stars: Observations and Theory
Thermal and Ionization Aspects of Flows from Hot Stars: Observations and Theory ASP Conference Series, Vol. XXX, YYYY Henny J.G.L.M. Lamers, and Arved Sapar, eds. Outer Wind Evolution of Instability-Generated
More informationPIC simulations of laser interactions with solid targets
PIC simulations of laser interactions with solid targets J. Limpouch, O. Klimo Czech Technical University in Prague, Faculty of Nuclear Sciences and Physical Engineering, Břehová 7, Praha 1, Czech Republic
More information