Reminders: Show your work! As appropriate, include references on your submitted version. Write legibly!
|
|
- Gwen Gray
- 5 years ago
- Views:
Transcription
1 Phys Computer Simulation of Plasmas Homework # 4 (Project #1) Due Wednesday, October 22, 2014 Reminders: Show your work! As appropriate, include references on your submitted version. Write legibly! Instructions - Do not do every problem!!! Choose one problem to do and talk to Paul about it. No problem is to be done by more than one student in the class. Problems will be reserved on a first come, first served basis. Treat the provided instructions as a guideline, but be open to taking the problem in a different direction (in consultation with Paul). This one problem is expected to take you approximately 1.5 times a usual homework set. Therefore, you should start early (!). This project is weighted higher than the previous homeworks, so be clear in your responses and take pride in your work - this is a great opportunity to learn some really cool physics! 1. Test Particle Simulations of Velocity Shear due to Spatially Varying Electric Fields - This problem explores systems with a non-uniform electric field in a uniform magnetic field; the goal is to show that finite Larmor radius effects self-consistently introduce a density gradient. This is important because the density gradient can excite drift waves and instabilities, which may be important at the Earth s magnetopause and in the magnetotail plasma sheet. Consider an electric field of the form ( x E(r) = E 0 [2 + tanh ˆx, L)] where E 0 and L are constants, with a magnetic field where B 0 is a constant. B(r) = B 0 ẑ, (a) Predict the direction and magnitude (in real units, not normalized) of the E B drift for particles immersed in this field far from x = 0 (for both positive and negative x). (b) In non-normalized units, calculate the Larmor radius of ions for arbitrary (nonrelativistic) perpendicular speeds (ignoring the presence of the electric field) in terms of B 0. (c) Normalize the fields so that L = 2, B 0 = 1, and E 0 = Add an option to the test particle code that was provided in class to update this combination of electric and magnetic fields. (d) Modify the code to allow an option for multiple particles instead of just a single one. To do so, put a loop around the existing code and make appropriate changes to the way the code does output. Put the initial conditions inside the loop over particles so they can be different for each particle. (e) Perform simulations for a particle distribution with initial positions between 10 < x < 10 in steps of 1 and multiple initial velocities v y between -2.5 and 2.5 in steps of (Thus, there will be 441 particles total.) Confirm this range of speeds includes high enough values that there are particles with Larmor radius larger than L. Calculate the density by making a histogram of the final x position of particles. Show that the density near the gradient in electric field is a function of x.
2 2. Test Particle Simulations of a Magnetic Mirror/Converging Magnetic Fields - An example of a magnetic mirror field in cylindrical (r, θ, z) coordinates is with where B 0 and L are constants. B(r) = 1 ψ θ, (1) 2π ( ) ψ(r) = B 0 πr z2, L 2 (a) By dotting both sides of Eq. (1) with ψ, show that magnetic fields lines are confined to surfaces of constant ψ. (b) Plot a collection of lines of constant ψ in the r z plane for a collection of initial r values at fixed z (using the contour function in idl). It should reveal that the magnetic fields have the right structure to produce magnetic mirroring. (c) Calculate the magnetic field in rectangular coordinates (you can leave it in terms of r and θ, but have the components be ˆx, ŷ, and ẑ). (d) Normalize this magnetic field so that z is in units of L and the strength of the magnetic field for z = 0 is 1. Add this magnetic field configuration (with no electric field) to the test particle code given out in class. (e) Modify the test particle code to allow an option for multiple particles instead of just a single one. To do so, put a loop around the existing code and make appropriate changes to the way the code does output. Put the initial conditions inside the loop over particles so they can be different for each particle. (f) Run the simulation for electrons (q/m < 0) that start in the equatorial plane (z = 0) at (x, y, z) = (2, 0, 0) with a distribution of initial parallel v 0 and perpendicular v 0 velocities in steps of 0.1 from 0.5 < v y0 < 0.5 and 0.5 < v z < 0.5. (Thus, there will be 121 particles - note, you should print out data less often by increasing skip to 100 or so to avoid printing too much data.) Show the initial distribution of particles on a plot of v vs. v. Repeat for the final velocities of the same particles after enough time stepping that the particles bounce (200,000 time steps per particle). (g) Modify the idl code to read in the data for multiple particles. At each time step, use the µ analysis portion of the code to find the maximum and minimum of z using the max(z) and min(z) commands. Then, after this has been calculated for every particle, plot the initial speeds v 0 vs. v 0 for all particles that did not go outside z = 1. Repeat for the particles that did escape to z 1. What do you find for the distribution of particles that did not escape? (h) The expected result is that particles initially within a cone in v 0 v 0 space escape. (See, any plasma textbook, such as Bellan s Fundamentals of Plasma Physics, for the theory.) Analytically calculate the loss cone angle for this system and compare it to the numerical results you obtained. (i) How does v 0 compare to v 0 for particles that escape past z = 1? Why would this be useful to the operation of scanning electron microscopes?
3 3. Test Particle Simulations of RF Capacitively Coupled Plasmas (Fermi Acceleration) - Consider a 1D parallel plate capacitor with a plasma in the interior 0 < x < d. When the plate is charged, a thin sheath is created near the plate. The electric field at the capacitor plates is made to oscillate in the radio frequency (RF) range. This makes the thickness of the sheath s change in time. Here, we consider the single particle dynamics of electrons in this system. (a) The simplest way to estimate the sheath thickness s is to assume that ions have a uniform density n i and are stationary on relevant time scales and that there are no electrons present. (Note, this is not a very accurate model...) Integrate the 1D Poisson equation (analytically) and show that the associated electric field is E x (x) = en i ɛ 0 x + E 0, where e is the proton charge and E 0 is the electric field at x = 0. Thus, the electric field goes from E 0 (which is negative) at the surface linearly up to zero; the thickness of the sheath s is therefore s = ɛ 0 E 0 /en i. Thus, the electric field can be written as [ E x (x) = E 0 1 x ] s for 0 < x < s. Similarly, the sheath at the other plate has an electric field of [ E x (x) = E 0 1 d x ] s for d s < x < d. The minus sign enforces that the two sheaths are out of phase. The electric field for s < x < d s is zero. Now, the charge on the capacitor is made to oscillate in time. We model this as s evolving in time as s = s 0 (1 + cos ωt)/2, where ω is the frequency and s 0 is the maximum sheath thickness. (b) Normalize this electric field so that E 0 = 1 and d = 1 and code it into the test particle code given in class (with no magnetic field). Add a feature to the code where the simulation will stop if the particle goes outside 0 < x < d. Let s 0 = d/4. (c) Put an electron at x = d/2, and give it some initial speed in the x direction. Show that as you vary the initial speed, there are some values for which the particle rebounds at a lower speed and others where it rebounds with a higher speed. By plotting the sheath size as a function of time when the particle hits, determine what causes the particle to get a kick. (This is an example of what is called Fermi acceleration, which is important in many settings in plasma physics and astrophysics!) (d) For a particle that gained energy from the interaction, put particles in different initial positions and with different initial speeds that allow the particle to interact with the same sheath in the same way. Calculate the energy gain as a function of initial speed. You should find that the energy gain is proportional to a power of the velocity; what power is it? This specifies the order of the Fermi acceleration. (e) Modify the test particle code to allow an option for multiple particles instead of just a single one. To do so, put a loop around the existing code and make appropriate changes to the way the code does output. Put the initial conditions inside the loop over particles so they can be different for each particle. By plotting trajectories of multiple particles, find one that interacts sequentially with multiple sheaths. Electrons of this type gain a lot of energy, which may help generate the plasma!
4 4. Test Particle Simulations in a Reversed Field (Speiser Orbits) - There has been effort to understand aspects of dynamics in the Earth s magnetotail from single particle motion. The magnetotail consists of a current sheet with current in the y direction. The Earth s dipole magnetic field effectively adds a normal magnetic field in the z direction. We investigate single particle dynamics in simplified versions of these fields. Consider a reversing magnetic field of the form ( z ) B(z) = B 0 tanh ˆx, w where B 0 and w are constants. This magnetic field is associated with a current sheet in the y direction (it is called a Harris sheet ), which itself is associated with an electric field in the y direction. In a steady state, the electric field is uniform: where E 0 is a constant. E = E 0 ŷ, (a) Predict what electrons (q/m < 0) positioned at z = ±3w and released from rest would do in these fields. (b) Normalize the fields so that B 0 = 1 and E 0 = 0.1 and lengths are measured in terms of w. Code these electric and magnetic fields into the test particle code provided in class. (c) Run the simulation for a particle starting at rest at (x, y, z) = (0, 0, 3w). Describe the motion for z > 1 [make a quantitative comparison to part (a)!]. What qualitatively changes for z < 1? What is happening physically? (d) How does the kinetic energy of the particle change in time for z < 1? Is this perpendicular or parallel energy? Based on the physical situation, calculate the rate of change analytically and compare it to your numerical results (quantitatively!). (e) Interpret your results in terms of the magnetic moment µ (using the provided idl code). How much does µ vary for z > 1? Where does µ conservation strongly break down? How does this compare to expectations? (f) Now add a small normal magnetic field, ( z ) B(z) = B 0 tanh ˆx + B 1 ẑ, w where B 1 is small compared to B 0. Perform simulations for a few different (small) values of B 1. For what strength does B 1 begin to qualitatively change the trajectories and energy gains? What happens after this change? What do you conclude the impact is on particle acceleration due to the presence of a normal magnetic field? Interpret the results in terms of magnetic moment conservation. This is important for understanding the onset of the tearing instability (the linear phase of magnetic reconnection) in the magnetotail. It also suggests that the current sheet immediately around a reconnection site is not expected to give rise to much particle acceleration (for the case of anti-parallel magnetic fields). (g) For the case with the normal field B 1 that changed the energy gains, add in a so-called guide field in the y direction: ( z ) B(z) = B 0 tanh ˆx + B 1 ẑ + B g ŷ, w where B g is a constant. How does the particle energization change with a guide field of 1? Is the energization perpendicular or parallel to the magnetic field?
5 5. Test Particle Simulations in Oscillating Electric Fields with Magnetic Fields (Chaotic Trajectories) - We have already studied wave-particle interactions for particles in a traveling sinusoidal electric field; now we add in the effect of a uniform magnetic field perpendicular to the electric field. Consider an electric field of the form E(r, t) = E 0 cos(kx ωt)ˆx, where E 0, k, and ω are constants, and a uniform magnetic field B(r, t) = B 0 ẑ, where B 0 is a constant. The behavior of single particles in these fields is surprisingly rich considering the simplicity of the fields. (a) Using Newton s second law, find the equation of motion describing the evolution of the x coordinate of the particles. Normalize the result with distances to 1/k and times to 1/Ω = m/qb, showing that the result is d 2 x dt 2 + x = α cos(x νt ), where primes denote normalized quantities and ν = ω/ω. What is α in terms of the physical parameters of the system? This equation has no exact solution - numerical simulations can help! (b) Code the normalized electric and magnetic fields into the test particle stepper given in class. (c) In the limit of ν = 0 and small α = 0.5, plot a phase space portrait (v x vs. x, as we did in Homework 2) for a number of particles starting from rest at different initial x positions (-15 to 15 in increments of 5). Repeat for large α 43 with initial x positions (-25 to 25 in increments of 5) starting from rest and for particles initially at x = 0 but with v x ranging from 0 to 10 in increments of 2. (d) Now, let ν be non-zero, but small (ν = 0.5). Repeat part (c). How do the phase space portraits differ? (e) What happens to the closed orbits in the high α, non-zero ν case? You should find that they get smaller as they move to the right. It turns out this is key for this process producing irreversible heating - as the closed orbits shrink, particles near the separatrix gain energy. Interestingly, when we studied just the traveling electric field with no magnetic field, we showed that the particles gain energy, but the gain is reversible. It is surprising that one can get irreversible heating merely by adding in a uniform magnetic field! (f) Look at the particle trajectories for some of the particles for some of the parameters you have used. In particular, fix ν and start with small α, and gradually increase it. Notice that the trajectories are quasi-periodic for small α, but then undergo a period doubling route to chaos. Explore this transition.
6 6. Test Particle Simulations in Stochastic Magnetic Fields - Magnetic fields in tokamaks can become stochastic, meaning they become space filling. It has been proposed that stochastic fields enhance transport, which is usually not desirable for tokamak operation. Here, we investigate this, but we do so in a simplified rectangular geometry. Consider a magnetic field of the form ( ) ( ) ( y 2πx 2πx B(r) = B g ẑ + B 0 tanh ˆx + B 10 cos ŷ + B 11 cos + 2πz ) ŷ L y /2 L x L x L z on a rectangular periodic domain of size L x L y L z. (a) Write a field line tracing code (normalize the magnetic fields to B 0 ). Note that dx dz = B x B z and dy dz = B y B z are first order ODEs, so you can solve them using Euler, Runge-Kutta, etc. Assume periodic boundary conditions, so that if the field line leaves one side of the domain it comes in at the other side. (If a coordinate goes outside, say, x = L x, then subtract L x from it; do this for each of the six faces of the domain.) (b) Carry out the field line tracing for a number (about 25) of magnetic field lines at equally spaced y values for a domain of size with a step size of dz = (There is nothing special about these parameters - just choosing something!) A common way to plot the results is with a Poincaré plot - every time the magnetic field passes through z = 0, make a point on a plot. After many traversals of the domain in z, a pattern emerges. Make this plot for B g = 2, B 10 = 0.05, B 11 = 0.0 and B g = 2, B 10 = 0.0, B 11 = (c) For these two parameters, you should see an island -type shape in the Poincaré plot. Show (numerically, not analytically) that the half thickness of the island is approximately equal to B Also show that the case with B 10 0 has the island centered at y = 0, while the case with B 11 0 has the island centered at y 11 = (L y /2) tanh 1 (B g L x /L z ) (This is the q = 1 rational surface. ) (d) Now make Poincaré plots when both modes are present. Repeat part (b) for the following parameters: (i) B 0 = 2, B 10 = B 11 = 0.05, (ii) B 0 = 2, B 10 = B 11 = 0.15, (iii) B 0 = 2, B 10 = B 11 = 0.3, and (iv) B 0 = 2, B 10 = B 11 = 0.5. (e) You should have found that the field lines are mostly regular for some parameters and space filling over some range of y for other parameters. The condition that stochastic fields dominate is called the Chirikov condition. This happens when the islands overlap; for the case here where B 10 = B 11, this occurs when the width B10 exceeds y 11. Show analytically that this is not too different from your results. (f) Now, code these magnetic fields (with zero electric field) into the test particle code and implement an option for periodic boundary conditions in the code. (g) For a good number of particles, start them at various locations in the domain with a speed of 0.1 in various directions for case (i) and (iii) from part (d). (You may want to modify the test particle code to allow for multiple particles by putting a loop around it and changing how the code does output.) What fraction of them travel a significant distance in y for each case? What does this imply about transport in a tokamak with stochastic fields? For a particle that travels a significant distance in y, how far has it traveled in z? What does this imply about transport in a tokamak with stochastic fields?
7 7. Test Particle Simulations of the Ponderomotive Force in Laser-Plasma Interactions - A charged particle interacting with a wave made up of two high-frequency components of comparable frequency feels an effective force at their beat frequency; this is the ponderomotive force. We investigate this physics here. Consider an electric field of the form E(r, t) = E 0 cos(k 1 x ω 1 t)ẑ + E 0 cos(k 2 x ω 2 t + φ)ẑ, where E 0, k 1, k 2, ω 1, ω 2 and φ are constants. (a) If these are electromagnetic (light) waves, what is the appropriate form for the magnetic field? (b) Normalize the equations so that E 0 = B 0 = 1 and choose k 1 = k 2 = 2π. Code these electric and magnetic fields into the test particle code provided in class. (c) Run a test particle simulation for an electron starting from rest. Choose ω 1 /k 1 = 1 and φ = 0.05, and make the frequencies of the two components differ by a small amount: ω 1 = ω On a time scale of the waves ω 1, show that the particle gains and loses energy in the z direction quasi-periodically; this is called the quiver motion and does not yield any net acceleration of the particles. (d) Show that on a time scale of the beat frequency between the two waves, you get unidirectional motion in the x direction. This is due to the ponderomotive force. Qualitatively, how does the x position change in time? Confirm the relation to the beat frequency by choosing a different ω 1 ω 2. (e) Repeat part (c) for different field strengths; use twice, four times, and eight times as much. How does the peak in v x change with field strength? (f) Repeat part (c) for different frequencies. How does the peak v x change with frequency? (g) Investigate the theory of the ponderomotive force and confirm that the numerical results are consistent with the theory. The ponderomotive force is a very effective way to accelerate particles - an off-the-shelf laser can accelerate electrons to hundreds of MeV! 8. Devise Your Own Project - Using any of the tools we have discussed so far (ODE solvers, test particle simulations, diffusion and parabolic PDE solvers, relaxation of elliptic PDEs, Grad-Shafranov, convection solvers and hyperbolic PDEs) devise your own project. Obviously, you will have to talk to me about it and ensure it is both challenging enough and not too challenging.
Single particle motion
Single particle motion Plasma is a collection of a very large number of charged particles moving in, and giving rise to, electromagnetic fields. Before going to the statistical descriptions, let us learn
More informationSingle particle motion and trapped particles
Single particle motion and trapped particles Gyromotion of ions and electrons Drifts in electric fields Inhomogeneous magnetic fields Magnetic and general drift motions Trapped magnetospheric particles
More informationSingle Particle Motion
Single Particle Motion C ontents Uniform E and B E = - guiding centers Definition of guiding center E gravitation Non Uniform B 'grad B' drift, B B Curvature drift Grad -B drift, B B invariance of µ. Magnetic
More informationMotion of Charged Particles in Fields
Chapter Motion of Charged Particles in Fields Plasmas are complicated because motions of electrons and ions are determined by the electric and magnetic fields but also change the fields by the currents
More informationcos 6 λ m sin 2 λ m Mirror Point latitude Equatorial Pitch Angle Figure 5.1: Mirror point latitude as function of equatorial pitch angle.
Chapter 5 The Inner Magnetosphere 5.1 Trapped Particles The motion of trapped particles in the inner magnetosphere is a combination of gyro motion, bounce motion, and gradient and curvature drifts. In
More informationMechanisms for particle heating in flares
Mechanisms for particle heating in flares J. F. Drake University of Maryland J. T. Dahlin University of Maryland M. Swisdak University of Maryland C. Haggerty University of Delaware M. A. Shay University
More informationPlasma Physics Prof. V. K. Tripathi Department of Physics Indian Institute of Technology, Delhi
Plasma Physics Prof. V. K. Tripathi Department of Physics Indian Institute of Technology, Delhi Module No. # 01 Lecture No. # 22 Adiabatic Invariance of Magnetic Moment and Mirror Confinement Today, we
More informationFormation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )
Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University
More informationProblems in Magnetostatics
Problems in Magnetostatics 8th February 27 Some of the later problems are quite challenging. This is characteristic of problems in magnetism. There are trivial problems and there are tough problems. Very
More informationReport submitted to Prof. P. Shipman for Math 540, Fall 2009
Dynamics at the Horsetooth Volume 1, 009. Three-Wave Interactions of Spin Waves Aaron Hagerstrom Department of Physics Colorado State University aaronhag@rams.colostate.edu Report submitted to Prof. P.
More informationPoint Vortex Dynamics in Two Dimensions
Spring School on Fluid Mechanics and Geophysics of Environmental Hazards 9 April to May, 9 Point Vortex Dynamics in Two Dimensions Ruth Musgrave, Mostafa Moghaddami, Victor Avsarkisov, Ruoqian Wang, Wei
More informationDOPPLER RESONANCE EFFECT ON ROTATIONAL DRIVE BY ION CYCLOTRON MINORITY HEATING
DOPPLER RESONANCE EFFECT ON ROTATIONAL DRIVE BY ION CYCLOTRON MINORITY HEATING V.S. Chan, S.C. Chiu, Y.A. Omelchenko General Atomics, San Diego, CA, U.S.A. 43rd Annual APS Division of Plasma Physics Meeting
More information(a) Write down the total Hamiltonian of this system, including the spin degree of freedom of the electron, but neglecting spin-orbit interactions.
1. Quantum Mechanics (Spring 2007) Consider a hydrogen atom in a weak uniform magnetic field B = Bê z. (a) Write down the total Hamiltonian of this system, including the spin degree of freedom of the electron,
More informationPlanetary Magnetospheres: Homework Problems
Planetary Magnetospheres: Homework Problems s will be posted online at http://www.ucl.ac.uk/ ucapnac 1. In classical electromagnetic theory, the magnetic moment µ L associated with a circular current loop
More informationIon Dynamics in Magnetic Reconnection Region
Ion Dynamics in Magnetic Reconnection Region Nehpreet Kaur Walia Guru Nanak Dev University, Punjab, India Masahiro Hoshino, Takanobu Amano The University of Tokyo, Tokyo, Japan The Maxwellian ion distribution
More informationAsymmetric Magnetic Reconnection in Coronal Mass Ejection Current Sheets
Asymmetric Magnetic Reconnection in Coronal Mass Ejection Current Sheets Nicholas Murphy, 1 Mari Paz Miralles, 1 Crystal Pope, 1,2 John Raymond, 1 Kathy Reeves, 1 Dan Seaton, 3 & David Webb 4 1 Harvard-Smithsonian
More information(Refer Slide Time: 01:17)
Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay Lecture No. 7 Heat Conduction 4 Today we are going to look at some one dimensional
More informationObservations on the ponderomotive force
Observations on the ponderomotive force D.A. Burton a, R.A. Cairns b, B. Ersfeld c, A. Noble c, S. Yoffe c, and D.A. Jaroszynski c a University of Lancaster, Physics Department, Lancaster LA1 4YB, UK b
More informationPhysics 106b: Lecture 7 25 January, 2018
Physics 106b: Lecture 7 25 January, 2018 Hamiltonian Chaos: Introduction Integrable Systems We start with systems that do not exhibit chaos, but instead have simple periodic motion (like the SHO) with
More informationPhys102 Lecture 16/17 Magnetic fields
Phys102 Lecture 16/17 Magnetic fields Key Points Electric Currents Produce Magnetic Fields Force on an Electric Current in a Magnetic Field; Definition of B Force on an Electric Charge Moving in a Magnetic
More informationM01M.1 Massive Spring
Part I Mechanics M01M.1 Massive Spring M01M.1 Massive Spring A spring has spring constant K, unstretched length L, and mass per unit length ρ. The spring is suspended vertically from one end in a constant
More informationLectures on basic plasma physics: Hamiltonian mechanics of charged particle motion
Lectures on basic plasma physics: Hamiltonian mechanics of charged particle motion Department of applied physics, Aalto University March 8, 2016 Hamiltonian versus Newtonian mechanics Newtonian mechanics:
More information3. In the adjacent figure, E 1 = 6.0 V, E 2 = 5.5 V, E 3 = 2.0 V, R 1 = 1W, and R 2 = 6W. All batteries are ideal. Find the current in resistor R 1.
1. A cylindrical copper rod of length L and cross-sectional area A is re-formed to twice its original length with no change in volume. If the resistance between its ends was originally R, what is it now?
More informationPlasma Physics Prof. V. K. Tripathi Department of Physics Indian Institute of Technology, Delhi
Plasma Physics Prof. V. K. Tripathi Department of Physics Indian Institute of Technology, Delhi Lecture No. # 09 Electromagnetic Wave Propagation Inhomogeneous Plasma (Refer Slide Time: 00:33) Today, I
More informationThe Linear Theory of Tearing Modes in periodic, cyindrical plasmas. Cary Forest University of Wisconsin
The Linear Theory of Tearing Modes in periodic, cyindrical plasmas Cary Forest University of Wisconsin 1 Resistive MHD E + v B = ηj (no energy principle) Role of resistivity No frozen flux, B can tear
More informationPlasma spectroscopy when there is magnetic reconnection associated with Rayleigh-Taylor instability in the Caltech spheromak jet experiment
Plasma spectroscopy when there is magnetic reconnection associated with Rayleigh-Taylor instability in the Caltech spheromak jet experiment KB Chai Korea Atomic Energy Research Institute/Caltech Paul M.
More informationDynamics of charged particles in spatially chaotic magnetic fields
PSFC/JA-1-38 Dynamics of charged particles in spatially chaotic magnetic fields Abhay K. Ram and Brahmananda Dasgupta a October 21 Plasma Science and Fusion Center, Massachusetts Institute of Technology
More informationPhysics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top
Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem
More informationPotential/density pairs and Gauss s law
Potential/density pairs and Gauss s law We showed last time that the motion of a particle in a cluster will evolve gradually, on the relaxation time scale. This time, however, is much longer than the typical
More informationCoordinate systems and vectors in three spatial dimensions
PHYS2796 Introduction to Modern Physics (Spring 2015) Notes on Mathematics Prerequisites Jim Napolitano, Department of Physics, Temple University January 7, 2015 This is a brief summary of material on
More informationPhysics of fusion power. Lecture 13 : Diffusion equation / transport
Physics of fusion power Lecture 13 : Diffusion equation / transport Many body problem The plasma has some 10 22 particles. No description is possible that allows for the determination of position and velocity
More informationE&M. 1 Capacitors. January 2009
E&M January 2009 1 Capacitors Consider a spherical capacitor which has the space between its plates filled with a dielectric of permittivity ɛ. The inner sphere has radius r 1 and the outer sphere has
More informationUNIVERSITY OF MISSOURI-COLUMBIA PHYSICS DEPARTMENT. PART I Qualifying Examination. August 20, 2013, 5:00 p.m. to 8:00 p.m.
UNIVERSITY OF MISSOURI-COLUMBIA PHYSICS DEPARTMENT PART I Qualifying Examination August 20, 2013, 5:00 p.m. to 8:00 p.m. Instructions: The only material you are allowed in the examination room is a writing
More informationPhysics 221A Fall 1996 Notes 13 Spins in Magnetic Fields
Physics 221A Fall 1996 Notes 13 Spins in Magnetic Fields A nice illustration of rotation operator methods which is also important physically is the problem of spins in magnetic fields. The earliest experiments
More informationPlasmas as fluids. S.M.Lea. January 2007
Plasmas as fluids S.M.Lea January 2007 So far we have considered a plasma as a set of non intereacting particles, each following its own path in the electric and magnetic fields. Now we want to consider
More informationElectromagnetic Theory (Hecht Ch. 3)
Phys 531 Lecture 2 30 August 2005 Electromagnetic Theory (Hecht Ch. 3) Last time, talked about waves in general wave equation: 2 ψ(r, t) = 1 v 2 2 ψ t 2 ψ = amplitude of disturbance of medium For light,
More informationIon trap. 3D Quadrupole Trap. -Mass spectrometer -Ion trap (Quadrupol Ion Storage, Quistor) 18. April
Ion trap 3D Quadrupole Trap -Mass spectrometer -Ion trap (Quadrupol Ion Storage, Quistor) 18. April 011 1 The trajectories of the ions can be given: They form Mathieu differential equations: with r + Ion
More informationElectricity & Magnetism Qualifier
Electricity & Magnetism Qualifier For each problem state what system of units you are using. 1. Imagine that a spherical balloon is being filled with a charged gas in such a way that the rate of charge
More informationOn the limiting behaviour of regularizations of the Euler Equations with vortex sheet initial data
On the limiting behaviour of regularizations of the Euler Equations with vortex sheet initial data Monika Nitsche Department of Mathematics and Statistics University of New Mexico Collaborators: Darryl
More informationChapter 28. Magnetic Fields. Copyright 2014 John Wiley & Sons, Inc. All rights reserved.
Chapter 28 Magnetic Fields Copyright 28-1 Magnetic Fields and the Definition of B The Definition of B The Field. We can define a magnetic field B to be a vector quantity that exists when it exerts a force
More informationTwo Fluid Dynamo and Edge-Resonant m=0 Tearing Instability in Reversed Field Pinch
1 Two Fluid Dynamo and Edge-Resonant m= Tearing Instability in Reversed Field Pinch V.V. Mirnov 1), C.C.Hegna 1), S.C. Prager 1), C.R.Sovinec 1), and H.Tian 1) 1) The University of Wisconsin-Madison, Madison,
More informationControl of Neo-classical tearing mode (NTM) in advanced scenarios
FIRST CHENGDU THEORY FESTIVAL Control of Neo-classical tearing mode (NTM) in advanced scenarios Zheng-Xiong Wang Dalian University of Technology (DLUT) Dalian, China Chengdu, China, 28 Aug, 2018 Outline
More informationHighlights from (3D) Modeling of Tokamak Disruptions
Highlights from (3D) Modeling of Tokamak Disruptions Presented by V.A. Izzo With major contributions from S.E. Kruger, H.R. Strauss, R. Paccagnella, MHD Control Workshop 2010 Madison, WI ..onset of rapidly
More informationMarching on the BL equations
Marching on the BL equations Harvey S. H. Lam February 10, 2004 Abstract The assumption is that you are competent in Matlab or Mathematica. White s 4-7 starting on page 275 shows us that all generic boundary
More informationColumbia University Department of Physics QUALIFYING EXAMINATION
Columbia University Department of Physics QUALIFYING EXAMINATION Monday, January 14, 2013 3:10PM to 5:10PM Classical Physics Section 2. Electricity, Magnetism & Electrodynamics Two hours are permitted
More informationSingle Particle Motion in a Magnetized Plasma
Single Particle Motion in a Magnetized Plasma Aurora observed from the Space Shuttle Bounce Motion At Earth, pitch angles are defined by the velocity direction of particles at the magnetic equator, therefore:
More informationApplications of Ampere s Law
Applications of Ampere s Law In electrostatics, the electric field due to any known charge distribution ρ(x, y, z) may alwaysbeobtainedfromthecoulomblaw it sauniversal tool buttheactualcalculation is often
More informationPhysics 102: Magnetic Fields
Physics 102: Magnetic Fields Assist. Prof. Dr. Ali Övgün EMU Physics Department www.aovgun.com Electric Field & Magnetic Field Electric forces acting at a distance through electric field. Vector field,
More informationModule 3: Electromagnetism
Module 3: Electromagnetism Lecture - Magnetic Field Objectives In this lecture you will learn the following Electric current is the source of magnetic field. When a charged particle is placed in an electromagnetic
More informationESG Fall 2006 Final Exam
ESG 8.022 Fall 2006 Final Exam Instructor: Michael Shaw Tuesday, December 19th, 1:30PM to 4:30PM Instructions Show work on all problems. Partial credit cannot be granted without adequate progress. Please
More informationThree Loose Ends: Edge Focusing; Chromaticity; Beam Rigidity.
Linear Dynamics, Lecture 5 Three Loose Ends: Edge Focusing; Chromaticity; Beam Rigidity. Andy Wolski University of Liverpool, and the Cockcroft Institute, Daresbury, UK. November, 2012 What we Learned
More information10/24/2012 PHY 102. (FAWOLE O.G.) Good day. Here we go..
Good day. Here we go.. 1 PHY102- GENERAL PHYSICS II Text Book: Fundamentals of Physics Authors: Halliday, Resnick & Walker Edition: 8 th Extended Lecture Schedule TOPICS: Dates Ch. 28 Magnetic Fields 12
More informationElectromagnetic Waves
Electromagnetic Waves Our discussion on dynamic electromagnetic field is incomplete. I H E An AC current induces a magnetic field, which is also AC and thus induces an AC electric field. H dl Edl J ds
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 informationPhysics 351 Monday, April 3, 2017
Physics 351 Monday, April 3, 2017 This weekend you read Chapter 11 (coupled oscillators, normal modes, etc.), but it will take us another day or two to finish Chapter 10 in class: Euler angles; Lagrangian
More informationTHE UNIVERSITY OF SYDNEY FACULTY OF SCIENCE INTERMEDIATE PHYSICS PHYS 2912 PHYSICS 2B (ADVANCED) ALL QUESTIONS HAVE THE VALUE SHOWN
CC0936 THE UNIVERSITY OF SYDNEY FACULTY OF SCIENCE INTERMEDIATE PHYSICS PHYS 91 PHYSICS B (ADVANCED) SEMESTER, 014 TIME ALLOWED: 3 HOURS ALL QUESTIONS HAVE THE VALUE SHOWN INSTRUCTIONS: This paper consists
More informationPHYS4210 Electromagnetic Theory Spring Final Exam Wednesday, 6 May 2009
Name: PHYS4210 Electromagnetic Theory Spring 2009 Final Exam Wednesday, 6 May 2009 This exam has two parts. Part I has 20 multiple choice questions, worth two points each. Part II consists of six relatively
More informationChapter 11. Radiation. How accelerating charges and changing currents produce electromagnetic waves, how they radiate.
Chapter 11. Radiation How accelerating charges and changing currents produce electromagnetic waves, how they radiate. 11.1.1 What is Radiation? Assume a radiation source is localized near the origin. Total
More informationfree space (vacuum) permittivity [ F/m]
Electrostatic Fields Electrostatic fields are static (time-invariant) electric fields produced by static (stationary) charge distributions. The mathematical definition of the electrostatic field is derived
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 21 Quantum Mechanics in Three Dimensions Lecture 21 Physics 342 Quantum Mechanics I Monday, March 22nd, 21 We are used to the temporal separation that gives, for example, the timeindependent
More informationFigure 1: Grad, Div, Curl, Laplacian in Cartesian, cylindrical, and spherical coordinates. Here ψ is a scalar function and A is a vector field.
Figure 1: Grad, Div, Curl, Laplacian in Cartesian, cylindrical, and spherical coordinates. Here ψ is a scalar function and A is a vector field. Figure 2: Vector and integral identities. Here ψ is a scalar
More informationDate: 1 April (1) The only reference material you may use is one 8½x11 crib sheet and a calculator.
PH1140: Oscillations and Waves Name: Solutions Conference: Date: 1 April 2005 EXAM #1: D2005 INSTRUCTIONS: (1) The only reference material you may use is one 8½x11 crib sheet and a calculator. (2) Show
More informationPlasma Processes. m v = ee. (2)
Plasma Processes In the preceding few lectures, we ve focused on specific microphysical processes. In doing so, we have ignored the effect of other matter. In fact, we ve implicitly or explicitly assumed
More informationUniversity of Illinois at Chicago Department of Physics
University of Illinois at Chicago Department of Physics Electromagnetism Qualifying Examination January 4, 2017 9.00 am - 12.00 pm Full credit can be achieved from completely correct answers to 4 questions.
More informationFigure 1: three bodies orbiting each other
PHAS2443 Practical Mathematics II Orbiting Bodies: The Dance of the Stars Introduction When one thinks of planetary motion a similar image of our solar system comes up where one massive body is being orbited
More informationStudent number: Question # Mark Maximum Mark. Multiple Choice 20
Name: Student number: Academic Honesty: Cheating in an examination includes the following: 1. the unauthorized sharing of material such as textbooks during an open book examination; 2. concealing information
More informationPhysics 240 Fall 2003: Final Exam. Please print your name: Please list your discussion section number: Please list your discussion instructor:
Physics 40 Fall 003: Final Exam Please print your name: Please list your discussion section number: Please list your discussion instructor: Form #1 Instructions 1. Fill in your name above. This will be
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 informationDate: 31 March (1) The only reference material you may use is one 8½x11 crib sheet and a calculator.
PH1140: Oscillations and Waves Name: SOLUTIONS AT END Conference: Date: 31 March 2005 EXAM #1: D2006 INSTRUCTIONS: (1) The only reference material you may use is one 8½x11 crib sheet and a calculator.
More informationPlasma Astrophysics Chapter 1: Basic Concepts of Plasma. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University
Plasma Astrophysics Chapter 1: Basic Concepts of Plasma Yosuke Mizuno Institute of Astronomy National Tsing-Hua University What is a Plasma? A plasma is a quasi-neutral gas consisting of positive and negative
More informationMeasurement of magnetic eld line stochasticity in nonlinearly evolving, Department of Engineering Physics,
Measurement of magnetic eld line stochasticity in nonlinearly evolving, nonequilibrium plasmas Y. Nishimura, J. D. Callen, and C. C. Hegna Department of Engineering Physics, University of Wisconsin-Madison,
More informationPhysics 2212 GH Quiz #4 Solutions Spring 2016
Physics 2212 GH Quiz #4 Solutions Spring 2016 I. (18 points) A bar (mass m, length L) is connected to two frictionless vertical conducting rails with loops of wire, in the presence of a uniform magnetic
More informationLecture 2: Plasma particles with E and B fields
Lecture 2: Plasma particles with E and B fields Today s Menu Magnetized plasma & Larmor radius Plasma s diamagnetism Charged particle in a multitude of EM fields: drift motion ExB drift, gradient drift,
More informationPlasma Physics Prof. Vijayshri School of Sciences, IGNOU. Lecture No. # 38 Diffusion in Plasmas
Plasma Physics Prof. Vijayshri School of Sciences, IGNOU Lecture No. # 38 Diffusion in Plasmas In today s lecture, we will be taking up the topic diffusion in plasmas. Diffusion, why do you need to study
More informationECE 3209 Electromagnetic Fields Final Exam Example. University of Virginia Solutions
ECE 3209 Electromagnetic Fields Final Exam Example University of Virginia Solutions (print name above) This exam is closed book and closed notes. Please perform all work on the exam sheets in a neat and
More informationRevision History Date Version Description Author 03/06/ Initial Version Sergio Miguel Martin. Table of Contents
Icarus Physics Engine Manual 1 - Basic Engine Operation Version 1.0 Revision History Date Version Description Author 03/06/2010 1.0 Initial Version Sergio Miguel Martin Table of Contents 1. Introduction
More informationMAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT
MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT ABSTRACT A. G. Tarditi and J. V. Shebalin Advanced Space Propulsion Laboratory NASA Johnson Space Center Houston, TX
More information8.2.2 Rudiments of the acceleration of particles
430 The solar wind in the Universe intergalactic magnetic fields that these fields should not perturb them. Their arrival directions should thus point back to their sources in the sky, which does not appear
More informationProblem set 3. Electromagnetic waves
Second Year Electromagnetism Michaelmas Term 2017 Caroline Terquem Problem set 3 Electromagnetic waves Problem 1: Poynting vector and resistance heating This problem is not about waves but is useful to
More informationDiffusion equation, flux, diffusion coefficient scaling. Diffusion in fully ionized plasma vs. weakly ionized plasma. n => Coulomb collision frequency
Last Time Diffusion in plasma: the main reason why we need to control it (i.e. using magnetic field) Diffusion equation, flux, diffusion coefficient scaling o o t nx,t Dn D2 nx,t o D ~ L 2 T Diffusion
More informationSimulation of alpha particle current drive and heating in spherical tokamaks
Simulation of alpha particle current drive and heating in spherical tokamaks R. Farengo 1, M. Zarco 1, H. E. Ferrari 1, 1 Centro Atómico Bariloche and Instituto Balseiro, Argentina. Consejo Nacional de
More informationMagnetostatics: Part 1
Magnetostatics: Part 1 We present magnetostatics in comparison with electrostatics. Sources of the fields: Electric field E: Coulomb s law. Magnetic field B: Biot-Savart law. Charge Current (moving charge)
More informationSawteeth in Tokamaks and their relation to other Two-Fluid Reconnection Phenomena
Sawteeth in Tokamaks and their relation to other Two-Fluid Reconnection Phenomena S. C. Jardin 1, N. Ferraro 2, J. Chen 1, et al 1 Princeton Plasma Physics Laboratory 2 General Atomics Supported by the
More informationGuiding Center Orbit Studies in a Tokamak Edge Geometry Employing Boozer and Cartesian Coordinates
Contrib. Plasma Phys. 48, No. -3, 4 8 (8) / DOI./ctpp.839 Guiding Center Orbit Studies in a Tokamak Edge Geometry Employing Boozer and Cartesian Coordinates Y. Nishimura,Y.Xiao,and Z. Lin Department of
More informationMagnetism II. Physics 2415 Lecture 15. Michael Fowler, UVa
Magnetism II Physics 2415 Lecture 15 Michael Fowler, UVa Today s Topics Force on a charged particle moving in a magnetic field Path of a charged particle moving in a magnetic field Torque on a current
More informationExam 2 Solutions. Applying the junction rule: i 1 Applying the loop rule to the left loop (LL), right loop (RL), and the full loop (FL) gives:
PHY61 Eam Solutions 1. [8 points] In the circuit shown, the resistance R 1 = 1Ω. The batter voltages are identical: ε1 = ε = ε3 = 1 V. What is the current (in amps) flowing through the middle branch from
More informationLet s consider nonrelativistic electrons. A given electron follows Newton s law. m v = ee. (2)
Plasma Processes Initial questions: We see all objects through a medium, which could be interplanetary, interstellar, or intergalactic. How does this medium affect photons? What information can we obtain?
More informationSolution Set Eight. 1 Problem #1: Toroidal Electromagnet with Gap Problem #4: Self-Inductance of a Long Solenoid. 9
: Solution Set Eight Northwestern University, Electrodynamics I Wednesday, March 9, 6 Contents Problem #: Toroidal Electromagnet with Gap. Problem #: Electromagnetic Momentum. 3 3 Problem #3: Type I Superconductor.
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationOverview of Experiments for Magnetic Torque
Overview of Experiments for Magnetic Torque General Description of Apparatus The Magnetic Torque instrument consists of a pair of Helmholtz like coils with a brass air bearing mounted in the middle. (The
More informationGyrokinetic simulations of magnetic fusion plasmas
Gyrokinetic simulations of magnetic fusion plasmas Tutorial 2 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr
More informationPhysics 112. Study Notes for Exam II
Chapter 20 Electric Forces and Fields Physics 112 Study Notes for Exam II 4. Electric Field Fields of + and point charges 5. Both fields and forces obey (vector) superposition Example 20.5; Figure 20.29
More informationClass 4: More Pendulum results
Class 4: More Pendulum results The pendulum is a mechanical problem with a long and interesting history, from Galileo s first ansatz that the period was independent of the amplitude based on watching priests
More informationGraduate Written Examination Spring 2014 Part I Thursday, January 16th, :00am to 1:00pm
Graduate Written Examination Spring 2014 Part I Thursday, January 16th, 2014 9:00am to 1:00pm University of Minnesota School of Physics and Astronomy Examination Instructions Part 1 of this exam consists
More informationWhere k = 1. The electric field produced by a point charge is given by
Ch 21 review: 1. Electric charge: Electric charge is a property of a matter. There are two kinds of charges, positive and negative. Charges of the same sign repel each other. Charges of opposite sign attract.
More informationPart III. Interaction with Single Electrons - Plane Wave Orbits
Part III - Orbits 52 / 115 3 Motion of an Electron in an Electromagnetic 53 / 115 Single electron motion in EM plane wave Electron momentum in electromagnetic wave with fields E and B given by Lorentz
More informationAdiabatic Approximation
Adiabatic Approximation The reaction of a system to a time-dependent perturbation depends in detail on the time scale of the perturbation. Consider, for example, an ideal pendulum, with no friction or
More informationPlasma Physics Prof. V. K. Tripathi Department of Physics Indian Institute of Technology, Delhi
Plasma Physics Prof. V. K. Tripathi Department of Physics Indian Institute of Technology, Delhi Module No. # 01 Lecture No. # 02 Plasma Response to Fields: Fluid Equations Well, friends today I am going
More informationMATH 308 COURSE SUMMARY
MATH 308 COURSE SUMMARY Approximately a third of the exam cover the material from the first two midterms, that is, chapter 6 and the first six sections of chapter 7. The rest of the exam will cover the
More informationMagnetic Reconnection: Recent Developments and Future Challenges
Magnetic Reconnection: Recent Developments and Future Challenges A. Bhattacharjee Center for Integrated Computation and Analysis of Reconnection and Turbulence (CICART) Space Science Center, University
More information