Exact canonical drift Hamiltonian formalism with pressure anisotropy and finite perturbed fields
|
|
- Marion Townsend
- 5 years ago
- Views:
Transcription
1 PHYSICS OF PLASMAS 14, Exact canonical drift Hailtonian foralis with pressure anisotropy and finite perturbed fields G. A. Cooper a epartent of Physics, University of the South, Sewanee, Tennessee 37383, USA and Ecole Polytechnique Fédérale de Lausanne, Centre de Recherches en Physique des Plasas, Association Eurato-Suisse, CH1015 Lausanne, Switzerland M. Jucker, W. A. Cooper, and J. P. Graves Ecole Polytechnique Fédérale de Lausanne, Centre de Recherches en Physique des Plasas, Association Eurato-Suisse, CH1015 Lausanne, Switzerland M. Yu. Isaev Nuclear Fusion Institute, RRC Kurchatov Institute, Moscow, Russia Received 10 August 2007; accepted 27 August 2007; published online 11 October 2007 A Hailtonian forulation of the guiding center drift orbits is extended to pressure anisotropy and field perturbations in axisyetric systes. The Boozer agnetic coordinates are shown to retain canonical properties in anisotropic pressure plasas with finite electrostatic perturbations and electroagnetic perturbed fields that solely affect the parallel coponent of the agnetic vector potential. The equations of otion developed in the Boozer coordinate frae are satisfied by direct verification of the drift velocities. A nuerical application illustrates the significance of retaining all second order ters Aerican Institute of Physics. OI: / I. INTROUCTION When driven by the neutral bea injection or ion cyclotron resonance heating required to achieve high teperatures in agnetic confineent systes, energetic particles behave like independent test particles whose otion can be accurately tracked using the guiding center drift approxiation. In addition, these energetic particles can generate a significant level of pressure anisotropy in the background equilibriu state. The Lagrangian for particle drift otion constitutes an invaluable tool in deterining canonical variables. 1,2 However, the canonical angular variables derived in Refs. 1 and 2 do not satisfy toroidal periodicity conditions. White and Zakharov have developed a drift Hailtonian forulation in generalized, nonstraight equilibriu agnetic field line coordinates that satisfies periodicity. 3 Yet, in practice this forulation appears to be highly cubersoe as ost guiding center orbit codes are ipleented in Boozer coordinates. The principle issue that lies in the transforation to Boozer coordinates, 4 where periodicity is guaranteed, usually involves neglecting higher order ters. 2,5 However, Wang 6 has developed a technique for axisyetric equilibriu fields in which the transforation to the Boozer coordinate frae retains all the relevant higher order ters neglected in previous treatents. The orbits can then be exactly identified with the drift velocity equation satisfying Liouville s theore, 7 V d e B + c B K B 1+ c 0, 1 B 2 where 0 KB. In particular, the corrections see to be of iportance when the plasa is at high and the particle has relatively large energy. In this article, we extend the ethod developed by Wang to deonstrate that the Boozer agnetic coordinates also retain their canonical properties in equilibriu fields sustained with anisotropic pressure and containing both finite electrostatic and nearly incopressible electroagnetic field perturbations. The Lagrangian of the syste is anipulated in Sec. II in order to define a set of canonical variables applicable to guiding center otion. In Sec. III, these variables for the basis of a Hailtonian forulation, obtaining equations of otion that are then transfored to a canonical Boozer coordinate frae. A nuerical application of the guiding center drifts is presented in Sec. IV, with a conversion to the notation of the VENUS code 8 outlined in the Appendix. Section V contains the conclusions and discussion. II. LAGRANGIAN ETERMINATION OF CANONICAL VARIABLES We begin by establishing, respectively, the covariant and contravariant representations of the equilibriu agnetic field and deriving its vector potential fro the latter, a Electronic ail: coopega0@sewanee.edu B = g + I + g,, X/2007/1410/102506/7/$ , Aerican Institute of Physics
2 Cooper et al. Phys. Plasas 14, B = + q = q 3 c,,p, 11 = =, where B is the equilibriu agnetic field,,, are the Boozer agnetic coordinates and +q. The basis for the difference in the representations of Eq. 2 between the anisotropic and isotropic pressure liits lies with the description of Apère s law. For anisotropic pressure plasas, Apère s law is H= 0 K, where HB while in the isotropic liit where =1, it is written as B= 0 j. The fields K and j correspond to the current density. Therefore, in the covariant Boozer representation, the agnetic field intensity H rather than its induction B is expressed in Eq. 2 under anisotropic pressure conditions. A valid interpretation is that the pressure anisotropy effectively odifies the pereability fro that of free space 0 to 0 /. We also find that the equilibriu vector potential of the agnetic field strength, A e, is. The drift approxiation of the oenta is defined as P = P B B + ea, where the parallel oentu P v, v is the particle velocity projection along the equilibriu field lines, while and e correspond to the particle ass and charge, respectively. The vector potential consists of equilibriu and perturbed coponents A=A e +A p, where A p is a finite perturbed function. We shall consider electroagnetic perturbations parallel to the equilibriu agnetic field only. Hence we define A p V,,,tB, 2,8 11 where =1 0 p p /B 2 is the pressure anisotropy paraeter with p and p representing the parallel and perpendicular pressures, respectively. 5 Thus the Lagrangian can be specified, noting that P /eb, L = P dx H =e c B + A e dx H, where c is defined as +V. By substituting for in the covariant representation of the agnetic field, we deterine B = h + g, g. We have defined, as,+q and h as I+gq, with the sybol denoting the derivative of a flux surface quantity with respect to. The Lagrangian becoes L e = cgd + c hd c g,d H e, 9 P = c g. 10 We know that the Lagrangian ust be of the for L = i p i dq i H in order to eet the conditions of a canonical coordinate syste. 1 Thus it is necessary to eliinate the d ter. This is facilitated by the introduction of the following: 6,,P P such that w,dw, d c = d +,d P,dP + d, L e = cgd c +,d + c h + P d P P,dP The equations of otion are invariant to the addition of a full tie differential ter to the Lagrangian. 1 Thus we introduce 6 S P L e w w,dw, ds = cgd c + c h + P P w w dwd By setting, Q,, 6 we deduce that the poloidal coponent of the oentu in the covariant representation is P = c h + c gq, +P Qw,dw. 17 We now reconsider and c as functions of the canonical variables P, P,, where c is an ignorable variable. Thus fro Eqs. 10 and 17, and setting gq+i +gi c gi c g 2,, we calculate that P = g, P = g q + QP, Q, + I g, = g c g Q P Q w dw, c = 1 cg, 21 P c = 1 cgq + QP, Q, + I P g 1 g, 22
3 Exact canonical drift Hailtonian foralis Phys. Plasas 14, c = 1 cg c g Q P Q d c dw. = e2 B 2 w 1 c g + g, 23 Recalling that =q c,, P, we also calculate the partial derivatives of the toroidal angle with respect to the canonical variables =,, P P =, + P,, P P = q, Q, + QP,, with the c differential of equal to the trivial negative unity. III. HAMILTONIAN FORMULATION OF RIFT ORBITS The Hailtonian associated with these canonical variables is actually H/e=e 2 2 B 2 /2+B/e+H e where the equations of otion are defined as dp = H e P,P, c, d = H e P P,, c, d c dp = H e, cp,p, = H e P P,, c. 27 The variation of H e with respect to the canonical variables is deterined using Eqs , dp = e2 B 2 1 c g + g, c g Q Q P w q + QP, Q,, dw 28 d = e2 B 2 1 c g + g,, 29 dp =, 30 q + QP, Q, e2 B 2 q + c I c g + I, + P,, 31 where the,, and ters are equal to 1 + 2, 1 + 2, and 1 + 2, respectively: 1 = e2 B 2 V, 2 = eb 2 + e B 1 = e2 B 2 2 = V, + eb2 eb 2 + e B +, 1 = e2 B 2 2 =. V, 2 +, Rather than following the oenta and the coordinate c which does not trivially satisfy periodicity, it is both ore convenient and physically intuitive to evaluate the particle otion in the Boozer coordinate frae and the parallel gyroradius. 1 By expanding the tie differential of P, P, and calculating the tie differential of =q we find that d = I g, d = e2 B 2 q + c I c g + g, I Furtherore we calculate the otion equation along by subtracting the tie differential of V,,,t fro that of c P, P,, d c = 1 cg q + I g, 40
4 Cooper et al. Phys. Plasas 14, FIG. 1. Color A bird s eye view of the torus showing the orbits of a trapped particle coencing at the vertical idplane. The blue trajectory corresponds to the standard forulation of guiding center otion, while the red orbit is integrated using Eqs with perturbed fields neglected. Under these conditions, the inclusion of full second order coponents solely affects the toroidal drift. This particular perspective gives a clear presentation of the toroidal displaceent of the two orbits. q + I g dv = 1 cg 1 + I V 2 V 2+ V t, 1 + g V 2 V 2 +, V 2 V 2 41 I V d = gv g V V g, 1+ cg g, V IV q ci + c g 2 V t. 2 42
5 Exact canonical drift Hailtonian foralis Phys. Plasas 14, FIG. 2. Color A full toroidal view illustrates another perspective of the toroidal displaceent of the two orbits. The blue trajectory corresponds to the standard forulation of the guiding center otion, while the red orbit is integrated using Eqs with perturbed fields neglected. Under these conditions, the inclusion of full second order coponents solely affects the toroidal drift. An evaluation of the projections of the drift velocity, V d, with respect to,,, and recovers precisely the Hailtonian forulation of the Boozer coordinate equations of otion calculated and expressed in Eqs. 29, 38, 39, and 42. The drift velocity that satisfies Liouville s theore is given by Eq It can be deonstrated that =V d, =V d, =V d, and =t +V d, where t = t V. ds = 0Js + 0Is, 43 v v d = e2 B 2 s + c 0 Is v 0Is s B s, v v 44 IV. APPLICATION WITH THE VENUS COE In order to ipleent Eqs. 29, 38, 39, and 42 into the VENUS code, 8,5 they ust first be converted to the notation eployed in the progra. When the required transforations are ade see the Appendix, the equations of otion becoe d = e2 B 2 e2 B 2 s + c 0 Js v 0Js s + B s v v c B s, 45 v
6 Cooper et al. Phys. Plasas 14, FIG. 3. Color A graph of the difference between the toroidal drifts with respect to c the toroidal drift using the conventional forulation of the orbit equation illustrates the deviation that arises when certain second order ters are ignored. The difference is calculated as c subtracted by n, where n is the toroidal drift otion calculated using Eq. 45 with perturbed fields neglected. V 0Js d = + IsV v s2 s + 0 Is V s + c 0 Is v v 2 0 Js V s + s + c 0 Js 2 V t + V V 2 v 2 v B s + c 2 v B s. 46 The ters containing B s, in soe for or another, were factored out because they are essentially those that have not previously been included in the VENUS code. These ters are all of second order and while previous treatents contained soe ters of this order, this forulation is significant as it contains all ters of O. In order to deonstrate the significance of retaining all second order ters, the equations of otion were stripped of all perturbations before being ipleented into the code. In the absence of perturbed fields, the full second order drift contributions to the guiding center orbits only alters the equation of otion for the toroidal drift. Our nuerical application with the VENUS code was considered for a 10 MeV trapped proton orbit at a high of 19%. In addition, the equilibriu agnetic field was obtained fro the VMEC code 12 with a central value of 5.6 T. We estiate that a 3.5 MeV particle in a 6.6 T field would recover the sae difference in orbits as that investigated. The toroidal displaceent is best observed fro the perspective of Fig. 1. The orbit described by the red curve has been integrated using Eqs while the blue curve represents the standard forulation of these equations of otion in which certain higher order ters are neglected. The axiu deviation between the two orbits occurs about halfway between the vertical idplane and the turning points of the orbits. A full toroidal view is presented in Fig. 2. Figure 3 is a graph of the difference of the two toroidal drifts with respect to the conventional toroidal drift in which certain second order ters are ignored. The iportance of including all second order ters is deonstrated as the percent difference plotted here reaches up to 7%. V. ISCUSSION The ethod developed by Wang 6 is extended to deonstrate that Boozer coordinates are also canonical for guiding center drift orbit otion in the presence of both arbitrary electrostatic perturbed fields and electroagnetic perturbations with vector potential parallel to the equilibriu agnetic field. The forulation is further extended to conditions where energetic particles drive an anisotropic pressure plasa. The equations of otion we have derived in our Hailtonian foralis are recovered exactly fro direct verification of the corresponding projections of the guiding center drift velocity. In addition, a nuerical application in the VENUS code was presented to illustrate the significance of retaining all second order ters. The introduction of S into our Lagrangian eliinates the d ter and establishes c and as our canonical coordinate variables. The Hailtonian derivation of these coordinates and their respective oenta develops a set of equations of otion that, while canonical, do not satisfy other criteria such as periodicity. Iportantly, when the required conversion to the Boozer reference frae is ade the equations of otion retain their canonical nature, facilitating the ipleentation of nuerical schees. It is to be noted that our forulation of the perturbed vector potential rests on the assuption ade extensively in existing codes of a relatively low ratio of kinetic to ag-
7 Exact canonical drift Hailtonian foralis Phys. Plasas 14, netic pressures,. 2,8 11 Furtherore, the forulation is calculated in axisyetry and thus is only applicable to tokaak and reverse field pinch conditions. ACKNOWLEGMENTS This research was partially sponsored by the Fonds Nationale Suisse de la Recherche Scientifique and Eurato. We thank r. S. P. Hirshan for the use of the VMEC code. APPENIX: CONVERSION TO VENUS NOTATION The following transforations were ade to Eqs. 29, 38, 39, and 42: g 0 I, I 0 J, g, B, gq + I gb 2. The VENUS code does not include as a coordinate variable. Instead the code is written in ters of s, where =s. The ter s is usually proportional to the plasa volue enclosed and varies fro 0 at the agnetic axis to unity at the plasavacuu interface. The following transforations were ade with representing any flux surface quantity, s, Also, s s. B B s s. Note that the ter becoes s where s s. Furtherore, our Jacobian in the VENUS notation becoes g = 0 J Iq/B 2. However the Jacobian of the VENUS code is in fact gv = 0 sjs sis/b 2, where sqs. Thus, v s. 1 R. G. Littlejohn, Phys. Fluids 28, R. B. White, Phys. Fluids B 2, R. B. White and L. E. Zakharov, Phys. Plasas 10, A. H. Boozer, Phys. Fluids 25, W. A. Cooper, J. P. Graves, M. Jucker, and M. Yu. Isaev, Phys. Plasas 13, S. Wang, Phys. Plasas 13, R. B. White, A. H. Boozer, and R. Hay, Phys. Fluids 25, O. Fischer, W. A. Cooper, M. Yu. Isaev, and L. Villard, Nucl. Fusion 42, R. B. White and M. S. Chance, Phys. Fluids 27, R. B. White and A. H. Boozer, Phys. Plasas 2, S.. Pinches, L. C. Appel, J. Candy, S. E. Sharapov, H. L. Berk,. Borba, B. N. Breizan, T. C. Hender, K. I. Hopcraft, G. T. A. Huysans, and W. Kerner, Coput. Phys. Coun. 111, S. P. Hirshan and J. C. Whitson, Phys. Fluids 26,
2 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 informationToroidal coupling of ideal magnetohydrodynamic instabilities in tokamak plasmas
Toroidal coupling of ideal agnetohydrodynaic instabilities in tokaak plasas C. C. Hegna, J. W. Connor, R. J. Hastie, and H. R. Wilson Citation: Phys. Plasas 3, 584 (996); doi: 0.063/.87886 View online:
More informationNew Schemes for Confinement of Fusion Products in Stellarators
New Schemes for Confinement of Fusion Products in Stellarators W.A. Cooper ), M.Yu. Isaev 1), M.F. Heyn 5), V.N. Kalyuzhnyj 3), S.V. Kasilov 3), W. Kernbichler 5), A.Yu. Kuyanov 1), M.I. Mikhailov 1),
More informationPHYSICS 110A : CLASSICAL MECHANICS MIDTERM EXAM #2
PHYSICS 110A : CLASSICAL MECHANICS MIDTERM EXAM #2 [1] Two blocks connected by a spring of spring constant k are free to slide frictionlessly along a horizontal surface, as shown in Fig. 1. The unstretched
More informationComparison of Charged Particle Tracking Methods for Non-Uniform Magnetic Fields. Hann-Shin Mao and Richard E. Wirz
42nd AIAA Plasadynaics and Lasers Conferencein conjunction with the8th Internati 27-30 June 20, Honolulu, Hawaii AIAA 20-3739 Coparison of Charged Particle Tracking Methods for Non-Unifor Magnetic
More informationPh 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 informationCHAPTER 4 TWO STANDARD SHORTCUTS USED TO TRANSFORM ELECTROMAGNETIC EQUATIONS 4.1 THE FREE-PARAMETER METHOD
CHAPTER 4 TWO STANDARD SHORTCUTS USED TO TRANSFORM ELECTROMAGNETIC EQUATIONS The last several chapters have explained how the standard rules for changing units apply to electroagnetic physical quantities.
More informationParticle interactions with spatially localized wavepackets
PSFC/JA-10-36 Particle interactions with spatially localized wavepackets Y. Koinis, a K. Hizanidis, a and A.K. Ra October 010 Plasa Science and Fusion Center, Massachusetts Institute of Technology Cabridge,
More informatione = n 1 ( ) 3 [ m 3] = n [ m 3] n
Magnetospheric Physics - Hoework Solutions, /7/4 7. Plasa definition Can a plasa be aintained at teperatures of T e K Hint: Calculate the density liit using the plasa paraeter and explain your result).
More informationEffects of an Inhomogeneous Magnetic Field (E =0)
Effects of an Inhoogeneous Magnetic Field (E =0 For soe purposes the otion of the guiding centers can be taken as a good approxiation of that of the particles. ut it ust be recognized that during the particle
More informationFinal Exam Classical Mechanics
Final Ea Classical Mechanics. Consider the otion in one diension of a article subjected to otential V= (where =constant). Use action-angle variables to find the eriod of the otion as a function of energ.
More informationScattering and bound states
Chapter Scattering and bound states In this chapter we give a review of quantu-echanical scattering theory. We focus on the relation between the scattering aplitude of a potential and its bound states
More informationMotion of Charges in Uniform E
Motion of Charges in Unifor E and Fields Assue an ionized gas is acted upon by a unifor (but possibly tie-dependent) electric field E, and a unifor, steady agnetic field. These fields are assued to be
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 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 informationRadiation Reaction in High-Intense Fields
Radiation Reaction in High-Intense Fields Keita Seto Extree Light Infrastructure Nuclear Physics (ELI-NP) / Horia Hulubei National Institute for R&D in Physics and Nuclear Engineering (IFIN-HH), 3 Reactorului
More information1 Graded problems. PHY 5246: Theoretical Dynamics, Fall November 23 rd, 2015 Assignment # 12, Solutions. Problem 1
PHY 546: Theoretical Dynaics, Fall 05 Noveber 3 rd, 05 Assignent #, Solutions Graded probles Proble.a) Given the -diensional syste we want to show that is a constant of the otion. Indeed,.b) dd dt Now
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 informationNumerical Studies of a Nonlinear Heat Equation with Square Root Reaction Term
Nuerical Studies of a Nonlinear Heat Equation with Square Root Reaction Ter Ron Bucire, 1 Karl McMurtry, 1 Ronald E. Micens 2 1 Matheatics Departent, Occidental College, Los Angeles, California 90041 2
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 informationExploration of Configurational Space for Quasi-isodynamic Stellarators with Poloidally Closed Contours of the Magnetic Field Strength
Exploration of Configurational Space for Quasi-isodynamic Stellarators with Poloidally Closed Contours of the Magnetic Field Strength V.R. Bovshuk 1, W.A. Cooper 2, M.I. Mikhailov 1, J. Nührenberg 3, V.D.
More informationData-Driven Imaging in Anisotropic Media
18 th World Conference on Non destructive Testing, 16- April 1, Durban, South Africa Data-Driven Iaging in Anisotropic Media Arno VOLKER 1 and Alan HUNTER 1 TNO Stieltjesweg 1, 6 AD, Delft, The Netherlands
More information7. Renormalization and universality in pionless EFT
Renoralization and universality in pionless EFT (last revised: October 6, 04) 7 7. Renoralization and universality in pionless EFT Recall the scales of nuclear forces fro Section 5: Pionless EFT is applicable
More informationPhysics 2107 Oscillations using Springs Experiment 2
PY07 Oscillations using Springs Experient Physics 07 Oscillations using Springs Experient Prelab Read the following bacground/setup and ensure you are failiar with the concepts and theory required for
More informationFour-vector, Dirac spinor representation and Lorentz Transformations
Available online at www.pelagiaresearchlibrary.co Advances in Applied Science Research, 2012, 3 (2):749-756 Four-vector, Dirac spinor representation and Lorentz Transforations S. B. Khasare 1, J. N. Rateke
More informationHORIZONTAL MOTION WITH RESISTANCE
DOING PHYSICS WITH MATLAB MECHANICS HORIZONTAL MOTION WITH RESISTANCE Ian Cooper School of Physics, Uniersity of Sydney ian.cooper@sydney.edu.au DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS ec_fr_b. This script
More informationQuantum Chemistry Exam 2 Take-home Solutions
Cheistry 60 Fall 07 Dr Jean M Standard Nae KEY Quantu Cheistry Exa Take-hoe Solutions 5) (0 points) In this proble, the nonlinear variation ethod will be used to deterine an approxiate solution for the
More informationVariation of Plasma Frequency with Applied Magnetic Fields in a Single Walled Carbon Nanotube
International Journal of Physics and Applications. ISSN 0974-3103 Volue 2, Nuber 1 (2010), pp. 39--43 International Research Publication House http://www.irphouse.co Variation of Plasa Frequency with Applied
More informationKinetic Theory of Gases: Elementary Ideas
Kinetic Theory of Gases: Eleentary Ideas 17th February 2010 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion
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 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 informationNewton's Laws. Lecture 2 Key Concepts. Newtonian mechanics and relation to Kepler's laws The Virial Theorem Tidal forces Collision physics
Lecture 2 Key Concepts Newtonian echanics and relation to Kepler's laws The Virial Theore Tidal forces Collision physics Newton's Laws 1) An object at rest will reain at rest and an object in otion will
More informationKinetic Theory of Gases: Elementary Ideas
Kinetic Theory of Gases: Eleentary Ideas 9th February 011 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion of
More informationDETECTION OF NONLINEARITY IN VIBRATIONAL SYSTEMS USING THE SECOND TIME DERIVATIVE OF ABSOLUTE ACCELERATION
DETECTION OF NONLINEARITY IN VIBRATIONAL SYSTEMS USING THE SECOND TIME DERIVATIVE OF ABSOLUTE ACCELERATION Masaki WAKUI 1 and Jun IYAMA and Tsuyoshi KOYAMA 3 ABSTRACT This paper shows a criteria to detect
More informationProblem Set II Solutions
Physics 31600 R. Wal Classical Mechanics Autun, 2002 Proble Set II Solutions 1) Let L(q, q; t) be a Lagrangian [where, as in class, q stans for (q 1,..., q n )]. Suppose we introuce new coorinates (Q 1
More information13 Harmonic oscillator revisited: Dirac s approach and introduction to Second Quantization
3 Haronic oscillator revisited: Dirac s approach and introduction to Second Quantization. Dirac cae up with a ore elegant way to solve the haronic oscillator proble. We will now study this approach. The
More informationSupplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators
Suppleentary Inforation for Design of Bending Multi-Layer Electroactive Polyer Actuators Bavani Balakrisnan, Alek Nacev, and Elisabeth Sela University of Maryland, College Park, Maryland 074 1 Analytical
More informationStern-Gerlach Experiment
Stern-Gerlach Experient HOE: The Physics of Bruce Harvey This is the experient that is said to prove that the electron has an intrinsic agnetic oent. Hydrogen like atos are projected in a bea through a
More informationP (t) = P (t = 0) + F t Conclusion: If we wait long enough, the velocity of an electron will diverge, which is obviously impossible and wrong.
4 Phys520.nb 2 Drude theory ~ Chapter in textbook 2.. The relaxation tie approxiation Here we treat electrons as a free ideal gas (classical) 2... Totally ignore interactions/scatterings Under a static
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 informationCharacteristics of Low-Temperature Plasmas Under Nonthermal Conditions A Short Summary
1 1 Characteristics of Low-Teperature Plasas Under Nontheral Conditions A Short Suary Alfred Rutscher 1.1 Introduction The concept of a plasa dates back to Languir (1928) and originates fro the fundaental
More information72. (30.2) Interaction between two parallel current carrying wires.
7. (3.) Interaction between two parallel current carrying wires. Two parallel wires carrying currents exert forces on each other. Each current produces a agnetic field in which the other current is placed.
More informationRelativity and Astrophysics Lecture 25 Terry Herter. Momenergy Momentum-energy 4-vector Magnitude & components Invariance Low velocity limit
Mo Mo Relativity and Astrophysics Lecture 5 Terry Herter Outline Mo Moentu- 4-vector Magnitude & coponents Invariance Low velocity liit Concept Suary Reading Spacetie Physics: Chapter 7 Hoework: (due Wed.
More informationProblem T1. Main sequence stars (11 points)
Proble T1. Main sequence stars 11 points Part. Lifetie of Sun points i..7 pts Since the Sun behaves as a perfectly black body it s total radiation power can be expressed fro the Stefan- Boltzann law as
More informationarxiv: v2 [hep-th] 16 Mar 2017
SLAC-PUB-6904 Angular Moentu Conservation Law in Light-Front Quantu Field Theory arxiv:70.07v [hep-th] 6 Mar 07 Kelly Yu-Ju Chiu and Stanley J. Brodsky SLAC National Accelerator Laboratory, Stanford University,
More informationModelling diabatic atmospheric boundary layer using a RANS-CFD code with a k-ε turbulence closure F. VENDEL
Modelling diabatic atospheric boundary layer using a RANS-CFD code with a k-ε turbulence closure F. VENDEL Florian Vendel 1, Guillevic Laaison 1, Lionel Soulhac 1, Ludovic Donnat 2, Olivier Duclaux 2,
More informationOn the PPN 1+2 Body Problem
On the PPN 1 Body Proble arxiv:gr-qc/971039v1 8 Dec 1997 D Şelaru, I Dobrescu Gravitational esearch Laboratory, 1-5 Mendeleev str, 70168 Bucharest, oania e-ail: dselaru@scou1ctro, idobrescu@scou1ctro June
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 informationInfluence lines for statically indeterminate structures. I. Basic concepts about the application of method of forces.
Influence lines for statically indeterinate structures I. Basic concepts about the application of ethod of forces. The plane frae structure given in Fig. is statically indeterinate or redundant with degree
More informationSOLUTIONS. PROBLEM 1. The Hamiltonian of the particle in the gravitational field can be written as, x 0, + U(x), U(x) =
SOLUTIONS PROBLEM 1. The Hailtonian of the particle in the gravitational field can be written as { Ĥ = ˆp2, x 0, + U(x), U(x) = (1) 2 gx, x > 0. The siplest estiate coes fro the uncertainty relation. If
More information2.141 Modeling and Simulation of Dynamic Systems Assignment #2
2.141 Modeling and Siulation of Dynaic Systes Assignent #2 Out: Wednesday Septeber 20, 2006 Due: Wednesday October 4, 2006 Proble 1 The sketch shows a highly siplified diagra of a dry-dock used in ship
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 informationENGI 3424 Engineering Mathematics Problem Set 1 Solutions (Sections 1.1 and 1.2)
ENGI 344 Engineering Matheatics Proble Set 1 Solutions (Sections 1.1 and 1.) 1. Find the general solution of the ordinary differential equation y 0 This ODE is not linear (due to the product y ). However,
More informationHEAT TRANSFER IN FERROFLUID IN CHANNEL WITH POROUS WALLS
HEAT TRANSFER IN FERROFLUID IN CHANNEL WITH POROUS WALLS T. Strek Institute of Applied Mechanics, Poznan University of Technology, Poland Abstract The viscous, two-diensional, incopressible and lainar
More informationBALLISTIC PENDULUM. EXPERIMENT: Measuring the Projectile Speed Consider a steel ball of mass
BALLISTIC PENDULUM INTRODUCTION: In this experient you will use the principles of conservation of oentu and energy to deterine the speed of a horizontally projected ball and use this speed to predict the
More informationIN A SENSE, every material is a composite, even if the
IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, VOL. 47, NO. 11, NOVEMBER 1999 2075 Magnetis fro Conductors and Enhanced Nonlinear Phenoena J. B. Pendry, A. J. Holden, D. J. Robbins, and W. J. Stewart,
More informationASSIGNMENT BOOKLET Bachelor s Degree Programme (B.Sc./B.A./B.Com.) MATHEMATICAL MODELLING
ASSIGNMENT BOOKLET Bachelor s Degree Prograe (B.Sc./B.A./B.Co.) MTE-14 MATHEMATICAL MODELLING Valid fro 1 st January, 18 to 1 st Deceber, 18 It is copulsory to subit the Assignent before filling in the
More informationClassical systems in equilibrium
35 Classical systes in equilibriu Ideal gas Distinguishable particles Here we assue that every particle can be labeled by an index i... and distinguished fro any other particle by its label if not by any
More informationBasic concept of dynamics 3 (Dynamics of a rigid body)
Vehicle Dynaics (Lecture 3-3) Basic concept of dynaics 3 (Dynaics of a rigid body) Oct. 1, 2015 김성수 Vehicle Dynaics Model q How to describe vehicle otion? Need Reference fraes and Coordinate systes 2 Equations
More informationA NEW APPROACH FOR CALCULATING AVERAGE CROSS SECTIONS IN THE UNRESOLVED ENERGY REGION
Nuclear Matheatical Coputational Sciences: A Century in Review, A Century Anew Gatlinburg, Tennessee, April 6-, 2003, on CD-ROM, Aerican Nuclear Society, LaGrange Park, IL (2003) A NEW APPROACH FOR CALCULATING
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 informationOscillatory Hydromagnetic Couette Flow in a Rotating System with Induced Magnetic Field *
CHAPTER-4 Oscillator Hdroagnetic Couette Flow in a Rotating Sste with Induced Magnetic Field * 4. Introduction Lainar flow within a channel or duct in the absence of agnetic field is a phenoenon which
More informationP235 Midterm Examination Prof. Cline
P235 Mier Exaination Prof. Cline THIS IS A CLOSED BOOK EXAMINATION. Do all parts of all four questions. Show all steps to get full credit. 7:00-10.00p, 30 October 2009 1:(20pts) Consider a rocket fired
More informationMass Spectrum and Decay Constants of Conventional Mesons within an Infrared Confinement Model
Mass Spectru and Decay Constants of Conventional Mesons within an Infrared Confineent Model Gurjav Ganbold (BLTP, JINR; IPT MAS (Mongolia)) in collaboration with: T. Gutsche (Tuebingen) M. A. Ivanov (Dubna)
More informationThe accelerated expansion of the universe is explained by quantum field theory.
The accelerated expansion of the universe is explained by quantu field theory. Abstract. Forulas describing interactions, in fact, use the liiting speed of inforation transfer, and not the speed of light.
More informationSolving initial value problems by residual power series method
Theoretical Matheatics & Applications, vol.3, no.1, 13, 199-1 ISSN: 179-9687 (print), 179-979 (online) Scienpress Ltd, 13 Solving initial value probles by residual power series ethod Mohaed H. Al-Sadi
More informationDefinition of Work, The basics
Physics 07 Lecture 16 Lecture 16 Chapter 11 (Work) v Eploy conservative and non-conservative forces v Relate force to potential energy v Use the concept of power (i.e., energy per tie) Chapter 1 v Define
More informationDESIGN OF MECHANICAL SYSTEMS HAVING MAXIMALLY FLAT RESPONSE AT LOW FREQUENCIES
DESIGN OF MECHANICAL SYSTEMS HAVING MAXIMALLY FLAT RESPONSE AT LOW FREQUENCIES V.Raachran, Ravi P.Raachran C.S.Gargour Departent of Electrical Coputer Engineering, Concordia University, Montreal, QC, CANADA,
More informationNuclear Physics (10 th lecture)
~Theta Nuclear Physics ( th lecture) Content Nuclear Collective Model: Rainwater approx. (reinder) Consequences of nuclear deforation o Rotational states High spin states and back bending o Vibrational
More informationarxiv: v1 [cond-mat.stat-mech] 22 Dec 2017
Notes on the Hybrid Monte Carlo Method arxiv:1712.08278v1 [cond-at.stat-ech] 22 Dec 2017 Jerey C. Paler, 1, Air Haji-Akbari, 2, Rakesh S. Singh, 2 Fausto Martelli, 3 Roberto Car, 3 Athanassios Z. Panagiotopoulos,
More informationAn Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period
An Approxiate Model for the Theoretical Prediction of the Velocity... 77 Central European Journal of Energetic Materials, 205, 2(), 77-88 ISSN 2353-843 An Approxiate Model for the Theoretical Prediction
More information1 (40) Gravitational Systems Two heavy spherical (radius 0.05R) objects are located at fixed positions along
(40) Gravitational Systes Two heavy spherical (radius 0.05) objects are located at fixed positions along 2M 2M 0 an axis in space. The first ass is centered at r = 0 and has a ass of 2M. The second ass
More informationCurrent-induced switching of a single-molecule magnet with an arbitrary oriented easy axis *
Materials Science-Poland, Vol. 5, No. 4, 007 Current-induced switching of a single-olecule agnet with an arbitrary oriented easy axis * M. MISIORNY 1 1, **, J. BARNAŚ 1 Departent of Physics, A. Mickiewicz
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 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 informationPhysical interpretation of the Riemann hypothesis
Physical interpretation of the Rieann hypothesis Ditry Pozdnyaov Faculty of Radiophysics and Coputer Technologies of Belarusian State University Nezavisiosty av4 Mins Belarus E-ail: pozdnyaov@tutby Keywords:
More informationGeneralized r-modes of the Maclaurin spheroids
PHYSICAL REVIEW D, VOLUME 59, 044009 Generalized r-odes of the Maclaurin spheroids Lee Lindblo Theoretical Astrophysics 130-33, California Institute of Technology, Pasadena, California 9115 Jaes R. Ipser
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 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 information760 S.K. El-Labany et al Vol. a concluding discussion is presented and a coparison with previous results is considered..basic equations and derivation
Vol No 7, July 003 cfl 003 Chin. Phys. Soc. 009-963/003/(07)/0759-06 Chinese Physics and IOP Publishing Ltd Modulational instability of a wealy relativistic ion acoustic wave in a war plasa with nontheral
More informationDevelopment of point source method and its practical significance
Water Science and Engineering, 9, (): 19-31 doi:1.388/j.issn.1674-37.9..3 http://kkb.hhu.edu.cn e-ail: wse@hhu.edu.cn Developent of point source ethod and its practical significance Bidya Sagar PANI* Civil
More informationThe path integral approach in the frame work of causal interpretation
Annales de la Fondation Louis de Broglie, Volue 28 no 1, 2003 1 The path integral approach in the frae work of causal interpretation M. Abolhasani 1,2 and M. Golshani 1,2 1 Institute for Studies in Theoretical
More informationRegular article Maximum radius of convergence perturbation theory: test calculations on Be, Ne, H 2 and HF
Theor Che Acc (2003) 110: 185 189 DOI 10.1007/s00214-003-0473-z Regular article Maxiu radius of convergence perturbation theory: test calculations on Be, Ne, H 2 and HF Kotaro Yokoyaa 1, Haruyuki Nakano
More informationOptical Properties of Plasmas of High-Z Elements
Forschungszentru Karlsruhe Techni und Uwelt Wissenschaftlishe Berichte FZK Optical Properties of Plasas of High-Z Eleents V.Tolach 1, G.Miloshevsy 1, H.Würz Project Kernfusion 1 Heat and Mass Transfer
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 informationEfficient Filter Banks And Interpolators
Efficient Filter Banks And Interpolators A. G. DEMPSTER AND N. P. MURPHY Departent of Electronic Systes University of Westinster 115 New Cavendish St, London W1M 8JS United Kingdo Abstract: - Graphical
More informationCelal S. Konor Release 1.1 (identical to 1.0) 3/21/08. 1-Hybrid isentropic-sigma vertical coordinate and governing equations in the free atmosphere
Celal S. Konor Release. (identical to.0) 3/2/08 -Hybrid isentropic-siga vertical coordinate governing equations in the free atosphere This section describes the equations in the free atosphere of the odel.
More informationAsymptotic equations for two-body correlations
Asyptotic equations for two-body correlations M. Fabre de la Ripelle Abstract. An asyptotic equation for two-body correlations is proposed for a large nubers of particles in the frae work of the Integro-Differential
More informationEnergy and Momentum: The Ballistic Pendulum
Physics Departent Handout -10 Energy and Moentu: The Ballistic Pendulu The ballistic pendulu, first described in the id-eighteenth century, applies principles of echanics to the proble of easuring the
More informationi ij j ( ) sin cos x y z x x x interchangeably.)
Tensor Operators Michael Fowler,2/3/12 Introduction: Cartesian Vectors and Tensors Physics is full of vectors: x, L, S and so on Classically, a (three-diensional) vector is defined by its properties under
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 information2.003 Engineering Dynamics Problem Set 2 Solutions
.003 Engineering Dynaics Proble Set Solutions This proble set is priarily eant to give the student practice in describing otion. This is the subject of kineatics. It is strongly recoended that you study
More informationMeasuring Temperature with a Silicon Diode
Measuring Teperature with a Silicon Diode Due to the high sensitivity, nearly linear response, and easy availability, we will use a 1N4148 diode for the teperature transducer in our easureents 10 Analysis
More informationA simple formula for the trapped fraction in tokamaks including the effect of triangularity
A simple formula for the trapped fraction in tokamaks including the effect of triangularity O. Sauter Centre de Recherches en Physique des plasmas, Ecole Polytechnique Fédérale de Lausanne, CRPP-EPFL,
More informationJOURNAL OF PHYSICAL AND CHEMICAL SCIENCES
JOURNAL OF PHYSIAL AND HEMIAL SIENES Journal hoepage: http://scienceq.org/journals/jps.php Review Open Access A Review of Siple Haronic Motion for Mass Spring Syste and Its Analogy to the Oscillations
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 informationthe static friction is replaced by kinetic friction. There is a net force F net = F push f k in the direction of F push.
the static friction is replaced by kinetic friction. There is a net force F net = F push f k in the direction of F push. Exaple of kinetic friction. Force diagra for kinetic friction. Again, we find that
More informationSensorless Control of Induction Motor Drive Using SVPWM - MRAS Speed Observer
Journal of Eerging Trends in Engineering and Applied Sciences (JETEAS) 2 (3): 509-513 Journal Scholarlink of Eerging Research Trends Institute in Engineering Journals, 2011 and Applied (ISSN: 2141-7016)
More informationNow multiply the left-hand-side by ω and the right-hand side by dδ/dt (recall ω= dδ/dt) to get:
Equal Area Criterion.0 Developent of equal area criterion As in previous notes, all powers are in per-unit. I want to show you the equal area criterion a little differently than the book does it. Let s
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 information