Investigations of Electromagnetic Space-Charge Effects in Beam Physics

Similar documents
D zone schemes

Technical Vibration - text 2 - forced vibration, rotational vibration

Ans: In the rectangular loop with the assigned direction for i2: di L dt , (1) where (2) a) At t = 0, i1(t) = I1U(t) is applied and (1) becomes

USING LOWER CLASS WEIGHTS TO CORRECT AND CHECK THE NONLINEARITY OF BALANCES

ME 141. Engineering Mechanics

r r r r r EE334 Electromagnetic Theory I Todd Kaiser

LECTURE 5. is defined by the position vectors r, 1. and. The displacement vector (from P 1 to P 2 ) is defined through r and 1.

Circuits 24/08/2010. Question. Question. Practice Questions QV CV. Review Formula s RC R R R V IR ... Charging P IV I R ... E Pt.

() t. () t r () t or v. ( t) () () ( ) = ( ) or ( ) () () () t or dv () () Section 10.4 Motion in Space: Velocity and Acceleration

5.1-The Initial-Value Problems For Ordinary Differential Equations

Electric Potential. and Equipotentials

Two-Pion Exchange Currents in Photodisintegration of the Deuteron

Molecular Dynamics Simulations (Leach )

Motion on a Curve and Curvature

Solvability of nonlinear Klein-Gordon equation by Laplace Decomposition Method

Some algorthim for solving system of linear volterra integral equation of second kind by using MATLAB 7 ALAN JALAL ABD ALKADER


Faraday s Law. To be able to find. motional emf transformer and motional emf. Motional emf

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

f(x) dx with An integral having either an infinite limit of integration or an unbounded integrand is called improper. Here are two examples dx x x 2

Mark Scheme (Results) January 2008

Sharif University of Technology - CEDRA By: Professor Ali Meghdari

Field due to a collection of N discrete point charges: r is in the direction from

v T Pressure Extra Molecular Stresses Constitutive equations for Stress v t Observation: the stress tensor is symmetric

P441 Analytical Mechanics - I. Coupled Oscillators. c Alex R. Dzierba

Modeling of the thermodilution curve for coronary flow assessment. Carel Stasse TU/e BMTE07.22

ESS 265 Spring Quarter 2005 Kinetic Simulations

Solutions to Midterm Physics 201

PHYSICS 102. Intro PHYSICS-ELECTROMAGNETISM

Chapter 10. Simple Harmonic Motion and Elasticity. Goals for Chapter 10

A New Approach to the Numerical Solution of Fractional Order Optimal Control Problems

Fourier-Bessel Expansions with Arbitrary Radial Boundaries

Lecture 22 Electromagnetic Waves

Supplementary material for " Coherent and Tunable Terahertz Radiation from Graphene Surface Plasmon Polarirons Excited by Cyclotron Electron Beam "

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

1. Kinematics of Particles

Physics 1502: Lecture 2 Today s Agenda

ECE Microwave Engineering. Fall Prof. David R. Jackson Dept. of ECE. Notes 10. Waveguides Part 7: Transverse Equivalent Network (TEN)

Average & instantaneous velocity and acceleration Motion with constant acceleration

Answers to test yourself questions

Narrow-band Receiver Radio Architectures

Contraction Mapping Principle Approach to Differential Equations

Week 10: DTMC Applications Ranking Web Pages & Slotted ALOHA. Network Performance 10-1

Homework 5 for BST 631: Statistical Theory I Solutions, 09/21/2006

General Non-Arbitrage Model. I. Partial Differential Equation for Pricing A. Traded Underlying Security

Two-dimensional Effects on the CSR Interaction Forces for an Energy-Chirped Bunch. Rui Li, J. Bisognano, R. Legg, and R. Bosch

7 Wave Equation in Higher Dimensions

Ch.4 Motion in 2D. Ch.4 Motion in 2D

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

Control Volume Derivation

ANALYSIS OF KINEMATICS AND KINETOSTATICS OF FOUR-BAR LINKAGE MECHANISM BASED ON GIVEN PROGRAM

Fundamental Vehicle Loads & Their Estimation

Physics 232 Exam I Feb. 13, 2006

Science Advertisement Intergovernmental Panel on Climate Change: The Physical Science Basis 2/3/2007 Physics 253

Reinforcement learning

Problems on transformer main dimensions and windings

Physics 232 Exam I Feb. 14, 2005

6. Gas dynamics. Ideal gases Speed of infinitesimal disturbances in still gas

Today - Lecture 13. Today s lecture continue with rotations, torque, Note that chapters 11, 12, 13 all involve rotations

Circular Motion. Radians. One revolution is equivalent to which is also equivalent to 2π radians. Therefore we can.

DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING FLUID MECHANICS III Solutions to Problem Sheet 3

Lecture-V Stochastic Processes and the Basic Term-Structure Equation 1 Stochastic Processes Any variable whose value changes over time in an uncertain

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 9: The High Beta Tokamak

M r. d 2. R t a M. Structural Mechanics Section. Exam CT5141 Theory of Elasticity Friday 31 October 2003, 9:00 12:00 hours. Problem 1 (3 points)

One of the common descriptions of curvilinear motion uses path variables, which are measurements made along the tangent t and normal n to the path of

Stress Analysis of Infinite Plate with Elliptical Hole

e t dt e t dt = lim e t dt T (1 e T ) = 1

MEEN 617 Handout #11 MODAL ANALYSIS OF MDOF Systems with VISCOUS DAMPING

Chapter 7. Interference

Lecture 5. Chapter 3. Electromagnetic Theory, Photons, and Light

QUESTION PAPER PHYSICS-PH

AQA Maths M2. Topic Questions from Papers. Circular Motion. Answers

IX.1.1 The Laplace Transform Definition 700. IX.1.2 Properties 701. IX.1.3 Examples 702. IX.1.4 Solution of IVP for ODEs 704

Physics 11b Lecture #11

The Production of Polarization

This immediately suggests an inverse-square law for a "piece" of current along the line.

TIME DELAY BASEDUNKNOWN INPUT OBSERVER DESIGN FOR NETWORK CONTROL SYSTEM

0 for t < 0 1 for t > 0

Motion. ( (3 dim) ( (1 dim) dt. Equations of Motion (Constant Acceleration) Newton s Law and Weight. Magnitude of the Frictional Force

Picking Coordinate Axes

Physics 201, Lecture 5

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

ÖRNEK 1: THE LINEAR IMPULSE-MOMENTUM RELATION Calculate the linear momentum of a particle of mass m=10 kg which has a. kg m s

Particle Horizons in Cosmology

SOLUTIONS TO CONCEPTS CHAPTER 11

PHYS PRACTICE EXAM 2

INTEGRALS. Exercise 1. Let f : [a, b] R be bounded, and let P and Q be partitions of [a, b]. Prove that if P Q then U(P ) U(Q) and L(P ) L(Q).

Orthotropic Materials

( ) ( ) Physics 111. Lecture 13 (Walker: Ch ) Connected Objects Circular Motion Centripetal Acceleration Centripetal Force Sept.

(A) 6.32 (B) 9.49 (C) (D) (E) 18.97

Chapter 25: Current, Resistance and Electromotive Force. ~10-4 m/s Typical speeds ~ 10 6 m/s

U>, and is negative. Electric Potential Energy

ρ θ φ δ δ θ δ φ δ φ π δ φ π δ φ π

International Journal of Pure and Applied Sciences and Technology

An Open cycle and Closed cycle Gas Turbine Engines. Methods to improve the performance of simple gas turbine plants

PHYSICS 1210 Exam 1 University of Wyoming 14 February points

Lecture 3. Electrostatics

Chapter 1. Model Theory

Chapter 28 Sources of Magnetic Field

Fig. 1S. The antenna construction: (a) main geometrical parameters, (b) the wire support pillar and (c) the console link between wire and coaxial

Transcription:

Invesigions of Elecogneic Spce-Chge Effecs in Be Physics Chong Shik Pk Depen of Physics Indin Univesiy Indin Univesiy Cycloon Fciliy @ effenson Lb Sep. 9

Ouline Moivions Exising Spce-Chge Modeling Spce-Chge Effecs in RF Phooinecos Spce-Chge Effecs wih Tnsvese Cuens Suy Fuue Plns 9//9 effeson Lb

Moivions Fuue cceleo-bsed expeiens dend high-bighness elecon bes fo phooinecos The elecon be hs low enegy so spce-chge foces cn be ipon elive o exenl gneic nd f fields Thee e wo in chllenges fo siulions of highbighness phooinecos: - Resoluion of sll lengh/ie scle spce-chge fields elive o long lengh/ie scles of ineco e.g. - ps bunch lenghs fo.3-.8 GH - Reovl of unphysicl siulion effecs such s nueicl gid dispesion nd nueicl Cheenkov effecs in FDTD ehods 9//9 effeson Lb 3

Exising Spce Chge Modeling 9//9 effeson Lb 4

Exising Elecosic Algoih SCHEFF (Spce Chge EFFec) - beks up he co-picles ino se of nnul ings - clcules he elecosic spce-chge foces in he be es fe nd Loen nsfos o he lb fe PARMELA (Phse And Rdil Moion in Elecon Line Acceleos) - Wokhose of phooineco design codes - iposes exenl f-fields copued fo SUPERFISH nd exenl gneic fields copued fo POISSON - cnno clcule wkefields self-consisenly ASTRA (A Spce chge TRcking Algoih) 9//9 effeson Lb 5

Exising Elecogneic Algoihs Yee/PIC lgoih (FDTD ehod) - solves Mxwell s equions on he wo ineleved E nd B gids MAFIA Nueicl Dispesion Nueicl Cheenkov Rdiion TREDI (Thee Diensionl Inecos) - Liend-Wieche Poenils - no conducing boundies 9//9 effeson Lb 6

Ohe EM SC Mehod Mode nlysis nd Seies Expnsions (Slh e l) - solves wve equions using seies expnsions nd Fouie nsfoion in nol odes. - clcules he spce-chge fields o biy ccucy fo given be chge nd cuen densiies - needs sufficien oun of eigenodes nd Fouie odes 9//9 effeson Lb 7

Spce Chge Effecs in RF Phooineco 9//9 effeson Lb 8

Elecogneic Spce-Chge Poenils nd Fields The elions of EM fields nd poenils A E φ B A Loen Guge A + c Wve Equions c φ φ ρ ε A µ effeson Lb ρ E B ε B E E B µ + c 9//9 *PRST-AB Mk Hess e l 9

Tie-Dependen Geen s Funcions Fo he specil cse of cuens in he xil diecion in n pipe wih chode he poenils e given by Tie Dependen Geen s Funcion: ( ) ( ) ( ) G d d ; 3 ρ ε φ φ ( ) ( ) ( ) G d d A A ; 3 µ chode wll side boundy A A φ Boundy Condiions: 9//9 ( ) ( ) G G c A δ δ φ ( ) ( ) ( ) ( ) ( ) ( ) [ ] + + + λ θ λ λ θ λ π θ θ φ n n n i n n n A k k e c G G chode A wll side A boundy G G G φ effeson Lb

Elecogneic Spce-Chge Fields Abiy pipe coss-secion wih longiudinl cuens Recngul pipe coss-secion wih biy cuens Cicul pipe coss-secion wih longiudinl cuens E c 3 ( ) ψ ( ) ( )[ Γ Γ+ ] ( ) n ψ * n ρ d d ε n Fo cicul pipe ψ n Γ ± ( ) ( λ ) θ ( λ ) n ( ) ( ) ± / ± ( n ) e ( ) iθ π + n whee c ± 9//9 λ ± effeson Lb

Nueicl Ipleenion IRPSS (Indin RF Phoochode Souce Siulo) IRPSS is 3-D elecogneic picle/slice pushing lgoih Hndles ellic boundy condiions self-consisenly (boundy will consis of chode cicul side wlls) Uses elecogneic Geen s funcion folis fo solving fields Tie-dependen Geen s funcion foulion hs he dvngeous popey h elecogneic wves fo igh elecon bunches Specify ecoy nd chge densiy of he be Copue E nd B due o spcechge Siule be dynics of es picles effeson Lb 9//9

Bunched Be Model The uli-sliced bunch odel is inoduced o genelie he finie sie bunch lengh. We se up unifoly spced (in ie) eqully chged slices which fo one coplee bunch. Needs n enough nube of slices Equion fo chge densiy (Cicled) Q ρ( ) ( ) Θ b i i b i δ π b i [ ( ) ] 9//9 effeson Lb 3 i

Eigenode Suions The equied eigenode nubes cn be esied fo he expnsion of he chge densiy of he be. The nube of nsvese eigenodes necessy fo ccuely deeining he fields is invesely popoionl o he nsvese sie of he be. [ ] ) ( ) ( ) ( ) ( ) ( Θ i n i n i n i i i b i b i b Q ψ ρ δ π ρ 9//9 4 b /.5 b /. effeson Lb

Eigenode Suion (con d) Fo he expnsion of he be densiy funcion we cn ge he inequliy b M >> π In ode o odel he fields wihin % ccucy i is necessy o su ove les odes coesponding o BNL.856 GH Phoochode gun. 9//9 effeson Lb 5

Nueicl Tie Inegion In he E clculions he longiudinl field songly depends on he be ecoy. The oscillion peiods of he Geen s funcions e deeined by he nsvese eigenode nube M nd he inegion sep sie '. A slle ' educes he pliude of he oscillion. c τ b <<. 9//9 effeson Lb 6

Benchk Modeling Q v -Q v Q v IRPSS siulion of disk bunch of chge eied ie fo he chode sufce oving unifoly wih speed v Anlyicl odel of wo disks of chge oving unifoly in opposie diecions fo ll ie nd inesecing 9//9 effeson Lb 7

Benchk Copison E vs. E vs. Blue: Benchk Red: IRPSS effeson Lb 9//9 8

Siuling Be Dynics wih ANL AWA.3 GH RF Phooineco The f fields wihin AWA gun e ppoxiely E E cos( k Ek E sin( k )sin( ω + ϕ) Ek Bθ cos( k )cos( ω + ϕ) c The longiudinl f fields ke be ecoies )sin( ω + ϕ ) ff E ϕ.3 GH 5 MV/ 65 deg pc nc. ps Qb b lse effeson Lb 9//9 9

Tnsvese Spce-Chge Foces Lge discepncies fo elie ie Sll diffeences fo le ie F / Q E + β ( cb) c/.5 c/.4 effeson Lb 9//9

Longiudinl Spce-Chge Foces The diffeence is slle hn nsvese cse Fo longe bunch lengh his will be incesed c/.5 c/.4 effeson Lb 9//9

E Spce-Chge Fields A key esul is h s he lse pulse lengh is incesed he discepncy beween elecosic nd elecogneic SC fields is incesed he chode when he bck of he bunch is eied. The sie of he discepncy lso depends on he be dius. Quliively s he be dius is incesed i.e. he be becoes oe pncke like nd he discepncy becoes slle. 9//9 effeson Lb

Be Dynics w/o Spce-Chge Effecs Wih only f fields picles e dilly focused ne he chode P /c vs. phsespce. c.3 c.3 c 9//9 effeson Lb 3

Be Dynics wih Spce-Chge Effecs.3 c Q pc.3 c effeson Lb 9//9 4

Be Dynics wih Spce-Chge Effecs.3 c.3 c Q 5 pc effeson Lb 9//9 5

Be Loss Mesueens By vying he oun of lse powe P (popoionl o bunch chge) We esue he be cuen I I P Blue: Elecosic Red: Elecogneic effeson Lb 9//9 6

Expeienl Mesueens Pefoed n expeienl be loss esueen on he.3 GH f gun he ANL AWA expeien. Below e plos of esued be chge vs. esued lse pulse inensiy. If no be loss ws o occu hen he plo should be line wih unifo slope. Howeve ciicl bunch chge. i.e. E (f) E (ciicl) fo fixed lse pulse lengh nd dius one would expec be loss o occu nd educion in he slope of he cuve. b 3.4 psec effeson Lb b.4 psec 9//9 7

Spce Chge Effecs wih Tnsvese Cuens 9//9 effeson Lb 8

Elecogneic Field fo Cicully Syeic Souces Genelie he exc folis fo he SC fields of cylindiclly syeic be in cicul conducing pipe Include he effec of he nsvese cuens Consuc elecogneic SC fields using he iedependen Geen s funcion ehod in he cylindicl conducing boundy condiions Cn odel he high SC doined syses such s highpowe icowve souces Cope o Elecosic (ES) esul which is fequenly used o odel high-powe icowve souces such s klyson 9//9 effeson Lb 9

Expnsions of Chge nd Cuen Densiies Be souce nd he syse e cylindiclly syeic* 9//9 3 ( ) ( ) ( ) ( ) θ θ ρ ρ ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) d d d d θ θ ρ ρ effeson Lb *Subied o IEEE Tnscion equion : Coninuiy + + ρ

TM Mode Spce-Chge Fields 9//9 3 EM fields geneed by nd ( ) E E ( ) G d d E µ ε ρ ( ) ( ) λ λ θ c G effeson Lb E c B E c B B E E E E B E B B + + + + Mxwell'sEquions : µ µ ε ρ θ θ θ θ

Nueicl Exple Siule dilly behing be Rdil be cuen bu no longiudinl be cuen Cope EM SC fields wih ES SC fields ρ ( ) Q π ( ) ( ) b ( ) θ ( b ) b ( ) [ θ ( L + ) θ ( L ) ] b ( ) + δ δ cos ( ω ) > effeson Lb 9//9 3

Nueicl Requieens Siil o he phooineco odeling Less nube of eigenodes due o lge be dius effeson Lb 9//9 33

Tnsvese Spce-Chge Fields wih Tnsvese Cuens ES fields e lso found using Geen s funcion ehod The be oscillion ss wih he iniil be dius b /.5 9.8 c ω c/./c. effeson Lb 9//9 34

Tnsvese Spce-Chge Fields wih Tnsvese Cuens (con d) Blue: Elecogneic Red: Elecosic.5T.5T.75T T effeson Lb 9//9 35

Longiudinl Spce-Chge Fields wih Tnsvese Cuens The effecs of he SC E fields e ipon ne he be edge /. /.5 effeson Lb 9//9 36

Suy Developed elecogneic spce-chge odels of elecon bes in he pesence of he conducing boundies Developed novel copuionl code IRPSS o copue he SC fields nueiclly Is cpble of siuling bes wih biily sll bunch Lenghs since i uses Geen s funcion ppoch Siuled he be dynics of he be ne he chode in he f phooineco Exended Geen s funcion ehods by including he nsvese cuens Invesiged he elecogneic SC fields fo dilly behing be oscillion effeson Lb 9//9 37

Fuue Plns Exend IRPSS code by including he effecs of iis(es) o disconinuiies of he cviy Ipove he code o self-consisenly clcule he ecoies due o boh he exenl nd SC fields Sudy how he be dynics e ffeced due o he SC fields in he designs of gneic focusing schees fo eince copensions Include n biy be cuens such s iuhlly vying cuens 9//9 effeson Lb 38

Thnk You! 9//9 effeson Lb 39