Vlasov - Poisson sytem of equations and Landau damping
|
|
- Harvey Williams
- 5 years ago
- Views:
Transcription
1 7/1/212 Kinetic theory Vlasov-Poisson syste 1 Vlasov - Poisson syte of equations and Landau daping.5 The kinetic approach to finding the dispersion relation of electrostatic waves uses the Vlasov equation to find the evolution of the distribution function in one diension: E j f f f v a t x v jt where f(x, is the distribution function, x is distance, v is velocity, t is tie and a is acceleration. Poisson 's equation is: de dx nq 1 f dv The Vlasov equation can be written where the ion density is assued to have a constant value n, and q and and have their usual eaning. Df Dt where D/Dt indicates a convective derivative. If we follow a group of particles, the vaue of f(x, reains constant on the trajectory in x,v space. For exaple, if a particle oves to x fro x-dx, and siultaneously to v fro v-d the value that f(x, will have at the new tie and place is the value of f(x-dx, v-d t -d. The following prescription ipleents the change and carries us ward one tie step: f( xvt ) f( x dxv dvt d where dx = v dt, dv = a dt, E is the electric field, and a = qe/. This procedure is repeated to find the evolution of f(x,. Cheng and Knorr (see reference) developed the following ethod that is accurate to second order. The ethod takes a half step in the dx direction, calculates E, takes a full step in the dv direction, and then a half step in the dx direction: dx 1. f tep1 ( xv ) f x v dvt dt 2. E 1 f v 2 tep1 ( xv ) d 3. f tep2 ( xv ) f tep1 ( xv dv) 4. f ( xvt ) f tep2 x dx v 2 where f tep denotes a teporary function. At each step, a cubic spline function is fit to f(x,v), and the spline function is used to find the values at x - dx/2 and v - dv by interpolation. The interpolation adds significantly to the tie needed to find a solution.
2 7/1/212 Kinetic theory Vlasov-Poisson syste 2 Diensionless equations: The following substitutions will ake the Vlasov equation and Poisson's equation diensionless: v v v t x xω p v t qe E nq 2 v t ω v t ω p f v p ε t f λ D ω p 2 where v t is the theral velocity and D is defined but not used. The diensionless equations are: f f v E t x de dx f dv f v x, t Diensionless Vlasov 1 Diensionless Poisson The x and v grids In x, our grid will span 12 Debye lengths: L 12 k Define the grid using 33 points: kax 32 k kax x L Δx x x k kax 1 Δx.375 The v grid will be -3 to + 3 in diensionless units: ax kax ax v 3 6 Δv v v Δv.188 ax 1 The v and x grids ust have the sae nuber of eleents the 2-d spline to work. The electrostatic wave The diensionless wave aplitude will be δ.25 =.3 is is the largest value that does not ake f(x,v) go negative (unphysical). fstart(x k,v ) is the density-odulated distribution function that we use at the start : fstart 1 δcos k 2πx k L e v 2 There is one wavelength in the siulation region. π f is noralized to unity in 1-d. Check noralization to unity: 2 v e Δv π Check ean squared-velocity: 2 v 2 v e Δv π
3 7/1/212 Kinetic theory Vlasov-Poisson syste 3 Plot of f(x,v): x is front to back and v is left to right fstart Poisson's equation: finding E fro f(x,v) Poisson equation finds de/dx at each x k : dedx 1 k fstart Δv k E is found by integrating de/dx, using Sipson's rule, which integrates across an interval using the average of the two derivative values at the ends of the interval. This ethod has second order accuracy. Estart k 1 kax By using the cosine function the density perturbation, we have a zero electric field at the left boundary. Now all the values k >, use Sipson's rule to integrate de/dx: Estart Estart.5 dedx dedx Δx k k1 k k1 Reset range of k: k kax dedx k.2 Estart k.2 Plot of initial E and de/dx values x k
4 7/1/212 Kinetic theory Vlasov-Poisson syste 4 Define the tie grid Our diensionless tie scale is defined by the plasa period. We will follow three oscillations which will put the ending tie at 3 x (2) = 18, approxiately. EndTie 18 For accuracy, we will use a tie step of.25. Δt.25 EndTie The nuber of iterations we will need is: jax ceil jax 72 Δt Periodic boundary conditions x jaxδt 18 If the particle oves to a value of x that is outside the siulations box, we ust use the odulo function od to put the x value into the box. This function will return the reainder after dividing by the length of the box L. If the value of x is less than zero, od will return a negative nuber that is outside our box. So we ust add L to the value of x bee using the od function: This returns a value x that is always between and L. od L x.5v ΔtL k Nuber density of particles in the siulation box We can find the nuber of particles in the box by integrating f(x,v) with dx and dv. The integral is pered by a double suation. We divide by L to find the nuber per unit length. nf ( ) kax 1 L kk 1.5 f f Δv kk kk1 Δx n( fstar For accuracy, we used Sispon's rule when integratring over x. It is not necessary to do this the v integral because the values of f(x,v) are nearly zero at the end points. Total energy of particles in the siulation box For the total kinetic energy W, we integrate (su).5 v 2 with the distribution function: W( f ) kax kk 1.5 f f.5 v kk kk1 2 Δv Δx W( fstar Total energy in the electric field This is found by suing.5 E 2 with dx, using Sipson's rule: Eenergy( Estar kax kk 1 Grid point ination.5.5 Estart kk 2 Estart kk1 2 Δx Eenergy( Estar The spline fitting routine needs to know the values of the x and v grid points, so we stack these vectors into a atrix with two coluns: xv augent( xv )
5 7/1/212 Kinetic theory Vlasov-Poisson syste 5 The progra loop j is tie index, k is distance index, is velocity index M f fstart f1 kaxax f2 kaxax E ax Esave Estart j jax SplineCoeffs ax k kax f1 k SplineCoeffs cspline( xvf) od L x.5v ΔtL k interpsplinecoeffsxvf v k kax ax dedx 1 k k 1 kax cspline( xvf1) f1 Δv k E E.5 dedx dedx Δx k k1 k k1 ax k kax f2 k SplineCoeffs f k Esave x k interpsplinecoeffsxvf1 v ax k kax augent( Esavef) cspline( xvf2) E Δt k od L x.5v ΔtL k interpsplinecoeffsxvf2 v augent( EsaveE) Initialize atrices f, f1, f2 Initialize vectors E, Esave. fit spline to f find f1 fit spline to f1 find de/dx find E k find f2 fit spline to f2 find the new f
6 7/1/212 Kinetic theory Vlasov-Poisson syste 6 Augent(Esave, E) is used to stack the E values at each tie step into an array Esave. Augent (Esave,f) is used to stack the final values f(x,v) into the atrix Esave that has the E values. The first thing we do below is unstack the E values and f(x,v) into separate atrices: This stateent puts the E(x, values into the atrix Ext: Ext subatrix( M kax jax 1) This stateent recovers the values of f(x,v) into a new atrix FinalF: FinalF subatrix( M kaxjax 2cols( M) 1) This subatrix has E(x,: Ext j jax E as a function of x at 3 ties Peak value of E as a fcn of tie Ext k Ext k Ext k Ext kax j x k jδt
7 7/1/212 Kinetic theory Vlasov-Poisson syste 7 Plot of E(x,. Ext Distance is along the horizontal axis and tie is on the vertical axis. The E field "seesaws" as tie increases. This plot is of f(x,v) at the end tie ( FinalF) Check here to see if f has gone negative: If so, the wave aplitude ust be reduced. in( FinalF)
8 7/1/212 Kinetic theory Vlasov-Poisson syste 8 Did the reduced wave energy show up as increased particle energy? The final E field is: Efinal Ext k kjax The change in the E field energy is: because jax is the end tie. Eenergy( Efinal) Eenergy( Estar where: Eenergy( Estar Eenergy( Efinal) The change in the particle energy is: W( FinalF) W( fstar where: W( fstar W( FinalF) 3.6 Energy is poorly conserved because we have used a rather coarse spacing in x, v and t in order to reduce the calculation tie. Try it: Change t fro.25 to.1 and notice that energy conservation is iproved. Modification of the distribution function The daping of the waves flattens the region of f(x,v) where v is near to the phase velocity of the wave. For a nicer plot, we will average the final value of F over all x. 1 faverage kax 1 1 k FinalF k e Copare this to: Fax v 2 π.1 faverage Fax A "bulge" in the distribution function is visible at a velocity near 2.5. Try it: Exaine the E( and calculate the wave frequency. The wavenuber is 2/L. Is the phase velocity /k near to the velocity where there is a bulge in the distribution function? v
9 7/1/212 Kinetic theory Vlasov-Poisson syste 9 k dependence of Landau daping The figures below show that Landau daping is ore severe shorter wavelengths: Daping L = 12 Daping L = Ext kax j 4.5 Ext kax j jt jt Try it: How does the observed daping copare with the daping fro theory? How heavily daped is the wave with L = 6? 4? Reference: C. Z. Cheng and G. Knorr, J. Cop. Phys. 22, 33 (1976).
Ph 20.3 Numerical Solution of Ordinary Differential Equations
Ph 20.3 Nuerical Solution of Ordinary Differential Equations Due: Week 5 -v20170314- This Assignent So far, your assignents have tried to failiarize you with the hardware and software in the Physics Coputing
More informationProjectile Motion with Air Resistance (Numerical Modeling, Euler s Method)
Projectile Motion with Air Resistance (Nuerical Modeling, Euler s Method) Theory Euler s ethod is a siple way to approxiate the solution of ordinary differential equations (ode s) nuerically. Specifically,
More informationReading from Young & Freedman: For this topic, read the introduction to chapter 25 and sections 25.1 to 25.3 & 25.6.
PHY10 Electricity Topic 6 (Lectures 9 & 10) Electric Current and Resistance n this topic, we will cover: 1) Current in a conductor ) Resistivity 3) Resistance 4) Oh s Law 5) The Drude Model of conduction
More informationCh 12: Variations on Backpropagation
Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith
More information2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all
Lecture 6 Introduction to kinetic theory of plasa waves Introduction to kinetic theory So far we have been odeling plasa dynaics using fluid equations. The assuption has been that the pressure can be either
More informationN-Point. DFTs of Two Length-N Real Sequences
Coputation of the DFT of In ost practical applications, sequences of interest are real In such cases, the syetry properties of the DFT given in Table 5. can be exploited to ake the DFT coputations ore
More informationPhysics 41 Homework #2 Chapter 16. fa. Here v is the speed of the wave. 16. The speed of a wave on a massless string would be infinite!
Physics 41 Hoewor # Chapter 1 Serway 7 th Conceptual: Q: 3,, 8, 11, 1, Probles P: 1, 3, 5, 9, 1, 5, 31, 35, 38, 4, 5, 57 Conceptual 3. (i) d=e, f, c, b, a (ii) Since, the saller the (the coefficient of
More informationANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER
IEPC 003-0034 ANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER A. Bober, M. Guelan Asher Space Research Institute, Technion-Israel Institute of Technology, 3000 Haifa, Israel
More informationDispersion. February 12, 2014
Dispersion February 1, 014 In aterials, the dielectric constant and pereability are actually frequency dependent. This does not affect our results for single frequency odes, but when we have a superposition
More informationNumerical Solution of the MRLW Equation Using Finite Difference Method. 1 Introduction
ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.1401 No.3,pp.355-361 Nuerical Solution of the MRLW Equation Using Finite Difference Method Pınar Keskin, Dursun Irk
More informationPhysics 41 HW Set 1 Chapter 15 Serway 7 th Edition
Physics HW Set Chapter 5 Serway 7 th Edition Conceptual Questions:, 3, 5,, 6, 9 Q53 You can take φ = π, or equally well, φ = π At t= 0, the particle is at its turning point on the negative side of equilibriu,
More informationA Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair
Proceedings of the 6th SEAS International Conference on Siulation, Modelling and Optiization, Lisbon, Portugal, Septeber -4, 006 0 A Siplified Analytical Approach for Efficiency Evaluation of the eaving
More informationAnalysis of Impulsive Natural Phenomena through Finite Difference Methods A MATLAB Computational Project-Based Learning
Analysis of Ipulsive Natural Phenoena through Finite Difference Methods A MATLAB Coputational Project-Based Learning Nicholas Kuia, Christopher Chariah, Mechatronics Engineering, Vaughn College of Aeronautics
More informationBlock designs and statistics
Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent
More informationOcean 420 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers
Ocean 40 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers 1. Hydrostatic Balance a) Set all of the levels on one of the coluns to the lowest possible density.
More informationPhysics 2210 Fall smartphysics 20 Conservation of Angular Momentum 21 Simple Harmonic Motion 11/23/2015
Physics 2210 Fall 2015 sartphysics 20 Conservation of Angular Moentu 21 Siple Haronic Motion 11/23/2015 Exa 4: sartphysics units 14-20 Midter Exa 2: Day: Fri Dec. 04, 2015 Tie: regular class tie Section
More informationPH 222-2C Fall Electromagnetic Oscillations and Alternating Current. Lectures 18-19
H - Fall 0 Electroagnetic Oscillations and Alternating urrent ectures 8-9 hapter 3 (Halliday/esnick/Walker, Fundaentals of hysics 8 th edition) hapter 3 Electroagnetic Oscillations and Alternating urrent
More informationChaotic Coupled Map Lattices
Chaotic Coupled Map Lattices Author: Dustin Keys Advisors: Dr. Robert Indik, Dr. Kevin Lin 1 Introduction When a syste of chaotic aps is coupled in a way that allows the to share inforation about each
More informationIn this chapter we will start the discussion on wave phenomena. We will study the following topics:
Chapter 16 Waves I In this chapter we will start the discussion on wave phenoena. We will study the following topics: Types of waves Aplitude, phase, frequency, period, propagation speed of a wave Mechanical
More informationIII.H Zeroth Order Hydrodynamics
III.H Zeroth Order Hydrodynaics As a first approxiation, we shall assue that in local equilibriu, the density f 1 at each point in space can be represented as in eq.iii.56, i.e. f 0 1 p, q, t = n q, t
More informationFigure 1: Equivalent electric (RC) circuit of a neurons membrane
Exercise: Leaky integrate and fire odel of neural spike generation This exercise investigates a siplified odel of how neurons spike in response to current inputs, one of the ost fundaental properties of
More informationFor a situation involving gravity near earth s surface, a = g = jg. Show. that for that case v 2 = v 0 2 g(y y 0 ).
Reading: Energy 1, 2. Key concepts: Scalar products, work, kinetic energy, work-energy theore; potential energy, total energy, conservation of echanical energy, equilibriu and turning points. 1.! In 1-D
More informationUSEFUL HINTS FOR SOLVING PHYSICS OLYMPIAD PROBLEMS. By: Ian Blokland, Augustana Campus, University of Alberta
1 USEFUL HINTS FOR SOLVING PHYSICS OLYMPIAD PROBLEMS By: Ian Bloland, Augustana Capus, University of Alberta For: Physics Olypiad Weeend, April 6, 008, UofA Introduction: Physicists often attept to solve
More informationFirst of all, because the base kets evolve according to the "wrong sign" Schrödinger equation (see pp ),
HW7.nb HW #7. Free particle path integral a) Propagator To siplify the notation, we write t t t, x x x and work in D. Since x i, p j i i j, we can just construct the 3D solution. First of all, because
More informationHee = ~ dxdy\jj+ (x) 'IJ+ (y) u (x- y) \jj (y) \jj (x), V, = ~ dx 'IJ+ (x) \jj (x) V (x), Hii = Z 2 ~ dx dy cp+ (x) cp+ (y) u (x- y) cp (y) cp (x),
SOVIET PHYSICS JETP VOLUME 14, NUMBER 4 APRIL, 1962 SHIFT OF ATOMIC ENERGY LEVELS IN A PLASMA L. E. PARGAMANIK Khar'kov State University Subitted to JETP editor February 16, 1961; resubitted June 19, 1961
More information( ). One set of terms has a ω in
Laptag Class Notes W. Gekelan Cold Plasa Dispersion relation Suer Let us go back to a single particle and see how it behaves in a high frequency electric field. We will use the force equation and Maxwell
More informationDesign of Sliding Mode Stabilizer for Wind Turbine Generator using Dynamic Compensation Observer Technique
Proceedings of the 6th WSES International Conference on Power Systes, Lisbon, Portugal, Septeber -4, 6 84 Design of Sliding Mode Stabilizer for Wind urbine Generator using Dynaic Copensation Observer echnique
More informationRECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS MEMBRANE
Proceedings of ICIPE rd International Conference on Inverse Probles in Engineering: Theory and Practice June -8, 999, Port Ludlow, Washington, USA : RECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS
More informationLecture 21. Interior Point Methods Setup and Algorithm
Lecture 21 Interior Point Methods In 1984, Kararkar introduced a new weakly polynoial tie algorith for solving LPs [Kar84a], [Kar84b]. His algorith was theoretically faster than the ellipsoid ethod and
More informationA note on the multiplication of sparse matrices
Cent. Eur. J. Cop. Sci. 41) 2014 1-11 DOI: 10.2478/s13537-014-0201-x Central European Journal of Coputer Science A note on the ultiplication of sparse atrices Research Article Keivan Borna 12, Sohrab Aboozarkhani
More informationQ5 We know that a mass at the end of a spring when displaced will perform simple m harmonic oscillations with a period given by T = 2!
Chapter 4.1 Q1 n oscillation is any otion in which the displaceent of a particle fro a fixed point keeps changing direction and there is a periodicity in the otion i.e. the otion repeats in soe way. In
More informationLecture 9 November 23, 2015
CSC244: Discrepancy Theory in Coputer Science Fall 25 Aleksandar Nikolov Lecture 9 Noveber 23, 25 Scribe: Nick Spooner Properties of γ 2 Recall that γ 2 (A) is defined for A R n as follows: γ 2 (A) = in{r(u)
More informationma x = -bv x + F rod.
Notes on Dynaical Systes Dynaics is the study of change. The priary ingredients of a dynaical syste are its state and its rule of change (also soeties called the dynaic). Dynaical systes can be continuous
More informationKinematics and dynamics, a computational approach
Kineatics and dynaics, a coputational approach We begin the discussion of nuerical approaches to echanics with the definition for the velocity r r ( t t) r ( t) v( t) li li or r( t t) r( t) v( t) t for
More informationExperimental Design For Model Discrimination And Precise Parameter Estimation In WDS Analysis
City University of New York (CUNY) CUNY Acadeic Works International Conference on Hydroinforatics 8-1-2014 Experiental Design For Model Discriination And Precise Paraeter Estiation In WDS Analysis Giovanna
More informationSeismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi.
Seisic Analysis of Structures by K Dutta, Civil Departent, II Delhi, New Delhi. Module 5: Response Spectru Method of Analysis Exercise Probles : 5.8. or the stick odel of a building shear frae shown in
More information2nd Workshop on Joints Modelling Dartington April 2009 Identification of Nonlinear Bolted Lap Joint Parameters using Force State Mapping
Identification of Nonlinear Bolted Lap Joint Paraeters using Force State Mapping International Journal of Solids and Structures, 44 (007) 8087 808 Hassan Jalali, Haed Ahadian and John E Mottershead _ Γ
More informationFinite fields. and we ve used it in various examples and homework problems. In these notes I will introduce more finite fields
Finite fields I talked in class about the field with two eleents F 2 = {, } and we ve used it in various eaples and hoework probles. In these notes I will introduce ore finite fields F p = {,,...,p } for
More informationPhysically Based Modeling CS Notes Spring 1997 Particle Collision and Contact
Physically Based Modeling CS 15-863 Notes Spring 1997 Particle Collision and Contact 1 Collisions with Springs Suppose we wanted to ipleent a particle siulator with a floor : a solid horizontal plane which
More informationMassachusetts Institute of Technology Quantum Mechanics I (8.04) Spring 2005 Solutions to Problem Set 4
Massachusetts Institute of Technology Quantu Mechanics I (8.04) Spring 2005 Solutions to Proble Set 4 By Kit Matan 1. X-ray production. (5 points) Calculate the short-wavelength liit for X-rays produced
More informationCHAPTER 15: Vibratory Motion
CHAPTER 15: Vibratory Motion courtesy of Richard White courtesy of Richard White 2.) 1.) Two glaring observations can be ade fro the graphic on the previous slide: 1.) The PROJECTION of a point on a circle
More informationTopic 5a Introduction to Curve Fitting & Linear Regression
/7/08 Course Instructor Dr. Rayond C. Rup Oice: A 337 Phone: (95) 747 6958 E ail: rcrup@utep.edu opic 5a Introduction to Curve Fitting & Linear Regression EE 4386/530 Coputational ethods in EE Outline
More informationPeriodic Motion is everywhere
Lecture 19 Goals: Chapter 14 Interrelate the physics and atheatics of oscillations. Draw and interpret oscillatory graphs. Learn the concepts of phase and phase constant. Understand and use energy conservation
More informationIntelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley osig 1 Winter Seester 2018 Lesson 6 27 February 2018 Outline Perceptrons and Support Vector achines Notation...2 Linear odels...3 Lines, Planes
More informationDepartment of Physics Preliminary Exam January 3 6, 2006
Departent of Physics Preliinary Exa January 3 6, 2006 Day 1: Classical Mechanics Tuesday, January 3, 2006 9:00 a.. 12:00 p.. Instructions: 1. Write the answer to each question on a separate sheet of paper.
More informationUsing a De-Convolution Window for Operating Modal Analysis
Using a De-Convolution Window for Operating Modal Analysis Brian Schwarz Vibrant Technology, Inc. Scotts Valley, CA Mark Richardson Vibrant Technology, Inc. Scotts Valley, CA Abstract Operating Modal Analysis
More informationPH 221-1D Spring Oscillations. Lectures Chapter 15 (Halliday/Resnick/Walker, Fundamentals of Physics 9 th edition)
PH 1-1D Spring 013 Oscillations Lectures 35-37 Chapter 15 (Halliday/Resnick/Walker, Fundaentals of Physics 9 th edition) 1 Chapter 15 Oscillations In this chapter we will cover the following topics: Displaceent,
More informationPhysics 201, Lecture 15
Physics 0, Lecture 5 Today s Topics q More on Linear Moentu And Collisions Elastic and Perfect Inelastic Collision (D) Two Diensional Elastic Collisions Exercise: Billiards Board Explosion q Multi-Particle
More informationExam 3 Solutions. 1. Which of the following statements is true about the LR circuit shown?
PHY49 Spring 5 Prof. Darin Acosta Prof. Paul Avery April 4, 5 PHY49, Spring 5 Exa Solutions. Which of the following stateents is true about the LR circuit shown? It is (): () Just after the switch is closed
More informationPH 221-2A Fall Waves - I. Lectures Chapter 16 (Halliday/Resnick/Walker, Fundamentals of Physics 9 th edition)
PH 1-A Fall 014 Waves - I Lectures 4-5 Chapter 16 (Halliday/Resnick/Walker, Fundaentals of Physics 9 th edition) 1 Chapter 16 Waves I In this chapter we will start the discussion on wave phenoena. We will
More informationSOLVING LITERAL EQUATIONS. Bundle 1: Safety & Process Skills
SOLVING LITERAL EQUATIONS Bundle 1: Safety & Process Skills Solving Literal Equations An equation is a atheatical sentence with an equal sign. The solution of an equation is a value for a variable that
More informationCHECKLIST. r r. Newton s Second Law. natural frequency ω o (rad.s -1 ) (Eq ) a03/p1/waves/waves doc 9:19 AM 29/03/05 1
PHYS12 Physics 1 FUNDAMENTALS Module 3 OSCILLATIONS & WAVES Text Physics by Hecht Chapter 1 OSCILLATIONS Sections: 1.5 1.6 Exaples: 1.6 1.7 1.8 1.9 CHECKLIST Haronic otion, periodic otion, siple haronic
More informationChapter 16 Solutions
Chapter 16 Solutions 16.1 Replace x by x vt = x 4.5t to get y = 6 [(x 4.5t) + 3] 16. y (c) y (c) y (c) 6 4 4 4 t = s t = 1 s t = 1.5 s 0 6 10 14 x 0 6 10 14 x 0 6 10 14 x y (c) y (c) 4 t =.5 s 4 t = 3
More informationON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD
PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY Physical and Matheatical Sciences 04,, p. 7 5 ON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD M a t h e a t i c s Yu. A. HAKOPIAN, R. Z. HOVHANNISYAN
More informationwhich proves the motion is simple harmonic. Now A = a 2 + b 2 = =
Worked out Exaples. The potential energy function for the force between two atos in a diatoic olecules can be expressed as follows: a U(x) = b x / x6 where a and b are positive constants and x is the distance
More informationLecture #8-3 Oscillations, Simple Harmonic Motion
Lecture #8-3 Oscillations Siple Haronic Motion So far we have considered two basic types of otion: translation and rotation. But these are not the only two types of otion we can observe in every day life.
More informationEFFECT OF MATERIAL PROPERTIES ON VIBRATIONS OF NONSYMMETRICAL AXIALLY LOADED THIN-WALLED EULER-BERNOULLI BEAMS
Matheatical and Coputational Applications, Vol. 5, No., pp. 96-07, 00. Association for Scientific Research EFFECT OF MATERIAL PROPERTIES ON VIBRATIONS OF NONSYMMETRICAL AXIALLY LOADED THIN-WALLED EULER-BERNOULLI
More informationGeneralized Rayleigh Wave Dispersion in a Covered Half-space Made of Viscoelastic Materials
Copyright 7 Tech Science Press CMC vol.53 no.4 pp.37-34 7 Generalized Rayleigh Wave Dispersion in a Covered Half-space Made of Viscoelastic Materials S.D. Akbarov and M. Negin 3 Abstract: Dispersion of
More information3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme
P-8 3D acoustic wave odeling with a tie-space doain dispersion-relation-based Finite-difference schee Yang Liu * and rinal K. Sen State Key Laboratory of Petroleu Resource and Prospecting (China University
More information26 Impulse and Momentum
6 Ipulse and Moentu First, a Few More Words on Work and Energy, for Coparison Purposes Iagine a gigantic air hockey table with a whole bunch of pucks of various asses, none of which experiences any friction
More informationAN APPLICATION OF CUBIC B-SPLINE FINITE ELEMENT METHOD FOR THE BURGERS EQUATION
Aksan, E..: An Applıcatıon of Cubıc B-Splıne Fınıte Eleent Method for... THERMAL SCIECE: Year 8, Vol., Suppl., pp. S95-S S95 A APPLICATIO OF CBIC B-SPLIE FIITE ELEMET METHOD FOR THE BRGERS EQATIO by Eine
More information6.2 Grid Search of Chi-Square Space
6.2 Grid Search of Chi-Square Space exaple data fro a Gaussian-shaped peak are given and plotted initial coefficient guesses are ade the basic grid search strateg is outlined an actual anual search is
More informatione-companion ONLY AVAILABLE IN ELECTRONIC FORM
OPERATIONS RESEARCH doi 10.1287/opre.1070.0427ec pp. ec1 ec5 e-copanion ONLY AVAILABLE IN ELECTRONIC FORM infors 07 INFORMS Electronic Copanion A Learning Approach for Interactive Marketing to a Custoer
More informationSimple Harmonic Motion
Reading: Chapter 15 Siple Haronic Motion Siple Haronic Motion Frequency f Period T T 1. f Siple haronic otion x ( t) x cos( t ). Aplitude x Phase Angular frequency Since the otion returns to its initial
More informationParticle-In-Cell Code for Two-Stream Instability
7// PIC codes Two-stream instability Particle-In-Cell Code for Two-Stream Instability This exercise uses the previously-developed PIC code to investigate the two-stream instability. The electron distribution
More informationChapter 6 1-D Continuous Groups
Chapter 6 1-D Continuous Groups Continuous groups consist of group eleents labelled by one or ore continuous variables, say a 1, a 2,, a r, where each variable has a well- defined range. This chapter explores:
More informationOBJECTIVES INTRODUCTION
M7 Chapter 3 Section 1 OBJECTIVES Suarize data using easures of central tendency, such as the ean, edian, ode, and idrange. Describe data using the easures of variation, such as the range, variance, and
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 informationTwo-loop evaluation of large Wilson loops with overlap fermions: the b-quark mass shift, and the quark-antiquark potential
arxiv:hep-lat/0509044v 4 Sep 2005 Two-loop evaluation of large Wilson loops with overlap ferions: the b-quark ass shift, and the quark-antiquark potential Departent of Physics, University of Cyprus, Nicosia
More informationWork, Energy and Momentum
Work, Energy and Moentu Work: When a body oves a distance d along straight line, while acted on by a constant force of agnitude F in the sae direction as the otion, the work done by the force is tered
More informationOSCILLATIONS AND WAVES
OSCILLATIONS AND WAVES OSCILLATION IS AN EXAMPLE OF PERIODIC MOTION No stories this tie, we are going to get straight to the topic. We say that an event is Periodic in nature when it repeats itself in
More informationPolygonal Designs: Existence and Construction
Polygonal Designs: Existence and Construction John Hegean Departent of Matheatics, Stanford University, Stanford, CA 9405 Jeff Langford Departent of Matheatics, Drake University, Des Moines, IA 5011 G
More informationReducing Vibration and Providing Robustness with Multi-Input Shapers
29 Aerican Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June -2, 29 WeA6.4 Reducing Vibration and Providing Robustness with Multi-Input Shapers Joshua Vaughan and Willia Singhose Abstract
More informationFlipping Physics Lecture Notes: Free Response Question #1 - AP Physics Exam Solutions
2015 FRQ #1 Free Response Question #1 - AP Physics 1-2015 Exa Solutions (a) First off, we know both blocks have a force of gravity acting downward on the. et s label the F & F. We also know there is a
More informationPhysics 139B Solutions to Homework Set 3 Fall 2009
Physics 139B Solutions to Hoework Set 3 Fall 009 1. Consider a particle of ass attached to a rigid assless rod of fixed length R whose other end is fixed at the origin. The rod is free to rotate about
More informationUnderstanding the coefficient of restitution (COR) using mass/spring systems
Understanding the coefficient of restitution (COR) using ass/spring systes Dr. David Kagan Departent of Physics California State University, Chico Chico, CA 9599-00 dkagan@csuchico.edu The coefficient
More informationIn this chapter we will study sound waves and concentrate on the following topics:
Chapter 17 Waves II In this chapter we will study sound waves and concentrate on the following topics: Speed of sound waves Relation between displaceent and pressure aplitude Interference of sound waves
More informationTime-of-flight Identification of Ions in CESR and ERL
Tie-of-flight Identification of Ions in CESR and ERL Eric Edwards Departent of Physics, University of Alabaa, Tuscaloosa, AL, 35486 (Dated: August 8, 2008) The accuulation of ion densities in the bea pipe
More informationUniform Approximation and Bernstein Polynomials with Coefficients in the Unit Interval
Unifor Approxiation and Bernstein Polynoials with Coefficients in the Unit Interval Weiang Qian and Marc D. Riedel Electrical and Coputer Engineering, University of Minnesota 200 Union St. S.E. Minneapolis,
More informationSIMPLE HARMONIC MOTION: NEWTON S LAW
SIMPLE HARMONIC MOTION: NEWTON S LAW siple not siple PRIOR READING: Main 1.1, 2.1 Taylor 5.1, 5.2 http://www.yoops.org/twocw/it/nr/rdonlyres/physics/8-012fall-2005/7cce46ac-405d-4652-a724-64f831e70388/0/chp_physi_pndul.jpg
More informationSlanted coupling of one-dimensional arrays of small tunnel junctions
JOURNAL OF APPLIED PHYSICS VOLUME 84, NUMBER 1 15 DECEMBER 1998 Slanted coupling of one-diensional arrays of sall tunnel junctions G. Y. Hu Departent of Physics and Astronoy, Louisiana State University,
More informationPoornima University, For any query, contact us at: , 18
AIEEE//Math S. No Questions Solutions Q. Lets cos (α + β) = and let sin (α + β) = 5, where α, β π, then tan α = 5 (a) 56 (b) 9 (c) 7 (d) 5 6 Sol: (a) cos (α + β) = 5 tan (α + β) = tan α = than (α + β +
More informationDonald Fussell. October 28, Computer Science Department The University of Texas at Austin. Point Masses and Force Fields.
s Vector Moving s and Coputer Science Departent The University of Texas at Austin October 28, 2014 s Vector Moving s Siple classical dynaics - point asses oved by forces Point asses can odel particles
More information= T. Oscillations and Waves. Example of an Oscillating System IB 12 IB 12
Oscillation: the vibration of an object Oscillations and Waves Eaple of an Oscillating Syste A ass oscillates on a horizontal spring without friction as shown below. At each position, analyze its displaceent,
More informationRandomized Recovery for Boolean Compressed Sensing
Randoized Recovery for Boolean Copressed Sensing Mitra Fatei and Martin Vetterli Laboratory of Audiovisual Counication École Polytechnique Fédéral de Lausanne (EPFL) Eail: {itra.fatei, artin.vetterli}@epfl.ch
More informationSimple and Compound Harmonic Motion
Siple Copound Haronic Motion Prelab: visit this site: http://en.wiipedia.org/wii/noral_odes Purpose To deterine the noral ode frequencies of two systes:. a single ass - two springs syste (Figure );. two
More informationVariations on Backpropagation
2 Variations on Backpropagation 2 Variations Heuristic Modifications Moentu Variable Learning Rate Standard Nuerical Optiization Conjugate Gradient Newton s Method (Levenberg-Marquardt) 2 2 Perforance
More informationProc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES
Proc. of the IEEE/OES Seventh Working Conference on Current Measureent Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Belinda Lipa Codar Ocean Sensors 15 La Sandra Way, Portola Valley, CA 98 blipa@pogo.co
More informationCOS 424: Interacting with Data. Written Exercises
COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well
More informationBootstrapping Dependent Data
Bootstrapping Dependent Data One of the key issues confronting bootstrap resapling approxiations is how to deal with dependent data. Consider a sequence fx t g n t= of dependent rando variables. Clearly
More informationThe Lagrangian Method vs. other methods (COMPARATIVE EXAMPLE)
The Lagrangian ethod vs. other ethods () This aterial written by Jozef HANC, jozef.hanc@tuke.sk Technical University, Kosice, Slovakia For Edwin Taylor s website http://www.eftaylor.co/ 6 January 003 The
More informationInteractive Markov Models of Evolutionary Algorithms
Cleveland State University EngagedScholarship@CSU Electrical Engineering & Coputer Science Faculty Publications Electrical Engineering & Coputer Science Departent 2015 Interactive Markov Models of Evolutionary
More informationAnswers to assigned problems from Chapter 1
Answers to assigned probles fro Chapter 1 1.7. a. A colun of ercury 1 in cross-sectional area and 0.001 in height has a volue of 0.001 and a ass of 0.001 1 595.1 kg. Then 1 Hg 0.001 1 595.1 kg 9.806 65
More informationNUMERICAL MODELLING OF THE TYRE/ROAD CONTACT
NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT PACS REFERENCE: 43.5.LJ Krister Larsson Departent of Applied Acoustics Chalers University of Technology SE-412 96 Sweden Tel: +46 ()31 772 22 Fax: +46 ()31
More informationPlasma-Wall Interaction: Sheath and Pre-sheath
Plasa-Wall Interaction: Sheath and Pre-sheath Under ost conditions, a very thin negative sheath appears in the vicinity of walls, due to accuulation of electrons on the wall. This is in turn necessitated
More informationNonmonotonic Networks. a. IRST, I Povo (Trento) Italy, b. Univ. of Trento, Physics Dept., I Povo (Trento) Italy
Storage Capacity and Dynaics of Nononotonic Networks Bruno Crespi a and Ignazio Lazzizzera b a. IRST, I-38050 Povo (Trento) Italy, b. Univ. of Trento, Physics Dept., I-38050 Povo (Trento) Italy INFN Gruppo
More informationCHAPTER 19: Single-Loop IMC Control
When I coplete this chapter, I want to be able to do the following. Recognize that other feedback algoriths are possible Understand the IMC structure and how it provides the essential control features
More informationm A 1 m mgd k m v ( C) AP Physics Multiple Choice Practice Oscillations
P Physics Multiple Choice Practice Oscillations. ass, attached to a horizontal assless spring with spring constant, is set into siple haronic otion. Its axiu displaceent fro its equilibriu position is.
More information12 Towards hydrodynamic equations J Nonlinear Dynamics II: Continuum Systems Lecture 12 Spring 2015
18.354J Nonlinear Dynaics II: Continuu Systes Lecture 12 Spring 2015 12 Towards hydrodynaic equations The previous classes focussed on the continuu description of static (tie-independent) elastic systes.
More informationThe Solution of One-Phase Inverse Stefan Problem. by Homotopy Analysis Method
Applied Matheatical Sciences, Vol. 8, 214, no. 53, 2635-2644 HIKARI Ltd, www.-hikari.co http://dx.doi.org/1.12988/as.214.43152 The Solution of One-Phase Inverse Stefan Proble by Hootopy Analysis Method
More information