Constant of Motion for a One- Dimensional and nth-order Autonomous System and Its Relation to the Lagrangian and Hamiltonian

Similar documents
PROJECT PROGRESS REPORT (03/lfi?lfibr-~/15/1998):

Bulk Modulus Capacitor Load Cells

Start-up Noise in 3-D Self-AmpMed

Los Alamos. Nova 2. Solutions of the Noh Problem for Various Equations of State Using Lie Groups LA-13497

(4) How do you develop an optimal signal detection technique from the knowledge of

A Comparison of Least Squares Linear Regression and Measurement Error Modeling of WadCold Multipole Correlation in SSC Prototype Dipole Magnets

PROJECT PROGRESS REPORT (06/16/1998-9/15/1998):

CQNl_" RESPONSE TO 100% INTERNAL QUANTUM EFFICIENCY SILICON PHOTODIODES TO LOW ENERGY ELECTRONS AND IONS

GA A24166 SUPER-INTENSE QUASI-NEUTRAL PROTON BEAMS INTERACTING WITH PLASMA: A NUMERICAL INVESTIGATION

BASAL CAMBRIAN BASELINE GEOLOGICAL CHARACTERIZATION COMPLETED

D. L. Strenge PNL-SA A GENERAL ALGORITHM FOR RADIOACTIVE DECAY WITH BRANCHING AND LOSS FROM A MEDIUM DISCLAIMER. July 1995

LQSAlamos National Laboratory

sample-specific X-ray speckle contrast variation at absorption edges $ & ~ 0

Plasma Response Control Using Advanced Feedback Techniques

Three-Dimensional Silicon Photonic Crystals

Modeling Laser and e-beam Generated Plasma-Plume Experiments Using LASNEX

Tell uric prof i 1 es across the Darrough Known Geothermal Resource Area, Nevada. Harold Kaufniann. Open-file Report No.

Alamos National Laboratory is operated by the University of California for the United States Department of Energy under contract W-7405-ENG-36

8STATISTICAL ANALYSIS OF HIGH EXPLOSIVE DETONATION DATA. Beckman, Fernandez, Ramsay, and Wendelberger DRAFT 5/10/98 1.

GA A27857 IMPACT OF PLASMA RESPONSE ON RMP ELM SUPPRESSION IN DIII-D

Los Alamos IMPROVED INTRA-SPECIES COLLISION MODELS FOR PIC SIMULATIONS. Michael E. Jones, XPA Don S. Lemons, XPA & Bethel College Dan Winske, XPA

PLASMA MASS DENSITY, SPECIES MIX AND FLUCTUATION DIAGNOSTICS USING FAST ALFVEN WAVE

Quarterly Report April 1 - June 30, By: Shirley P. Dutton. Work Performed Under Contract No.: DE-FC22-95BC14936

SPHERICAL COMPRESSION OF A MAGNETIC FIELD. C. M. Fowler DISCLAIMER

Plutonium 239 Equivalency Calculations

e, )2 Crossing Angles In The Beam-Beam Interaction R.H. Siernann* II. PERTURBATION FORMALISM SLAC-PUB-7225

Applications of Pulse Shape Analysis to HPGe Gamma-Ray Detectors

GA A23736 EFFECTS OF CROSS-SECTION SHAPE ON L MODE AND H MODE ENERGY TRANSPORT

N. Tsoupas, E. Rodger, J. Claus, H.W. Foelsche, and P. Wanderer Brookhaven National Laboratory Associated Universities, Inc. Upton, New York 11973

4kU. Measurement of Storage Ring Motion at the Advanced Light Source. QSTt ERNESTORLANDO LAWRENCE BERKELEYNATIONAL LABORATORY

Design of Current Leads for the MICE Coupling Magnet

UNDERSTANDING TRANSPORT THROUGH DIMENSIONLESS PARAMETER SCALING EXPERIMENTS

AC dipole based optics measurement and correction at RHIC

Ray M. Stringfield Robert J. Kares J. R. Thomas, Jr.

PREDICTIVE MODELING OF PLASMA HALO EVOLUTION IN POST-THERMAL QUENCH DISRUPTING PLASMAS

Turbulent Scaling in Fluids. Robert Ecke, MST-10 Ning Li, MST-10 Shi-Yi Li, T-13 Yuanming Liu, MST-10

A Geological and Geophysical Assessment of the Royal Center Gas Storage Field in North-Central Indiana, a Joint NIPSCO, DOE & GRI Case Study

Undulator Interruption in

ust/ aphysics Division, Argonne National Laboratory, Argonne, Illinois 60439, USA

RECXWH2 W/o s3-1

A Distributed Radiator, Heavy Ion Driven Inertial Confinement Fusion Target with Realistic, Multibeam Illumination Geometry

Contractor Name and Address: Oxy USA, Inc. (Oxy), Midland, Texas OBJECTIVES

National Accelerator Laboratory

Fission-Fusion Neutron Source

RWM FEEDBACK STABILIZATION IN DIII D: EXPERIMENT-THEORY COMPARISONS AND IMPLICATIONS FOR ITER

Peak Reliability Delivering near real-time phase angle deltas using Inter-Control Center Communication Protocol (ICCP)

RESEARCH IN ACCELERATOR PHYSICS (THEORY)

Final Technical Report. Department of Energy. for

Diffractive Dijet Search with Roman Pots at CDF

OAK RIDGE NATIONAL LABORATORY

GA A26474 SYNERGY IN TWO-FREQUENCY FAST WAVE CYCLOTRON HARMONIC ABSORPTION IN DIII-D

Quench Propagation Velocity for Highly Stabilized Conductors

Magnetic Measurements of the Elliptical Multipole Wiggler Prototype

DE '! N0V ?

Colliding Crystalline Beams

National Accelerator Laboratory

GA A27235 EULERIAN SIMULATIONS OF NEOCLASSICAL FLOWS AND TRANSPORT IN THE TOKAMAK PLASMA EDGE AND OUTER CORE

Scaling between K+ and proton production in nucleus-nucleus collisions *

Optimization of NSLS-II Blade X-ray Beam Position Monitors: from Photoemission type to Diamond Detector. P. Ilinski

Multi-Scale Chemical Process Modeling with Bayesian Nonparametric Regression

Analysis of Shane Telescope Aberration and After Collimation

A Few-Group Delayed Neutron Model Based on a Consistent Set of Decay Constants. Joann M. Campbell Gregory D. Spriggs

Gordon G. Parker Mechanical Engineering Dept., Michigan Institute of Technology Houghton, MI

August 3,1999. Stiffness and Strength Properties for Basic Sandwich Material Core Types UCRL-JC B. Kim, R.M. Christensen.

National Accelerator Laboratory

Alex Dombos Michigan State University Nuclear and Particle Physics

Data Comparisons Y-12 West Tower Data

RUDOLPH J. HENNINGER, XHM PAUL J. MAUDLIN, T-3 ERIC N. HARSTAD, T-3

Survey, Alignment and Beam Stability at the Advanced Light Source

GA A22689 SLOW LINER FUSION

RICE COAL COMBUSTION: EFFECT OF PROCESS CONDITIONS ON CHAR REACTIVITY. Quarterly Technical Report Performance Period: 10/1/94 42/31/94 (Quarter #13)

Hyperon Particle Physics at JHF. R. E. Mischke

GA A27806 TURBULENCE BEHAVIOR AND TRANSPORT RESPONSE APPROACHING BURNING PLASMA RELEVANT PARAMETERS

Development of a High Intensity EBIT for Basic and Applied Science

+.V) eo(o) -2 so - sx. hngc)90w3-- Beam-beam collisions and crossing angles in RHIC*

GA A THERMAL ION ORBIT LOSS AND RADIAL ELECTRIC FIELD IN DIII-D by J.S. degrassie, J.A. BOEDO, B.A. GRIERSON, and R.J.

INTERMOLECULAR POTENTIAL FUNCTIONS AND HIGH RESOLUTION MOLECULAR SPECTROSCOPY OF WEAKLY BOUND COMPLEXES. Final Progress Report

GA A22677 THERMAL ANALYSIS AND TESTING FOR DIII D OHMIC HEATING COIL

Use of Gamma Rays from the Decay of 13.8-sec "Be to Calibrate a Germanium Gamma Ray Detector for Measurements up to 8 MeV

GA A24016 PHYSICS OF OFF-AXIS ELECTRON CYCLOTRON CURRENT DRIVE

The Radiological Hazard of Plutonium Isotopes and Specific Plutonium Mixtures

GA A26686 FAST ION EFFECTS DURING TEST BLANKET MODULE SIMULATION EXPERIMENTS IN DIII-D

12/16/95-3/15/96 PERIOD MULTI-PARAMETER ON-LINE COAL BULK ANALYSIS. 2, 1. Thermal Neutron Flux in Coal: New Coal Container Geometry

National Accelerator Laboratory

Valley-Fill Sandstones in the Kootenai Formation on the Crow Indian Reservation, South-Central Montana

EXPLOSIVE PARTICLES PARVLENE ENCAPSUTION. lac0 b Sando v a t. Normal Process Development OCTOBER DEVEZOPMENT D I V I S I O N DEZEMBER 1971

Role of Flow Shear in Enhanced Core Confinement

Measurement of Low Levels of Alpha in 99M0Product Solutions

GA A25853 FAST ION REDISTRIBUTION AND IMPLICATIONS FOR THE HYBRID REGIME

Manifestations of the axial anomaly in finite temperature QCD

STP-TS THERMOPHYSICAL PROPERTIES OF WORKING GASES USED IN WORKING GAS TURBINE APPLICATIONS

GA A23713 RECENT ECCD EXPERIMENTAL STUDIES ON DIII D

Bandgap Photons to GaInP2

IMPACT OF EDGE CURRENT DENSITY AND PRESSURE GRADIENT ON THE STABILITY OF DIII-D HIGH PERFORMANCE DISCHARGES

Experiment. The Pinhole Neutron. Pinex. c. c. sartain UCRL-ID November 21,1958

41--tsh.. J. Halperin M. P. Lietzke ORNL-2823 EFFECTIVE CADMIUM CUTOFF ENERGIES SUPPLEMENT OAK RIDGE NATIONAL LABORATORY. R. W.

GA A26057 DEMONSTRATION OF ITER OPERATIONAL SCENARIOS ON DIII-D

Direct Determination of the Stacking Order in Gd 2 O 3 Epi-Layers on GaAs

A NEW TARGET CONCEPT FOR PROTON ACCELERATOR DRIVEN BORON NEUTRON CAPTURE THERAPY APPLICATIONS* Brookhaven National Laboratory. P.O.

GA A27805 EXPANDING THE PHYSICS BASIS OF THE BASELINE Q=10 SCENRAIO TOWARD ITER CONDITIONS

Transcription:

in u a, IK u SSCL-Preprint-543 December 1993 Distribution Category: 414 G. Lopez Constant of Motion for a One- Dimensional and nth-order Autonomous System and Its Relation to the Lagrangian and Hamiltonian Superconducting Super Collider Laboratory APPROVED FOR RELEASE OR PUBLICATION - O.R. PATENTS". djoup BY.. (IDT * 1>ATE. %L$f. IS*

Disclaimer Notice This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government or any agency thereof, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. Superconducting Super Collider Laboratory is an equal opportunity employer.

DISCLAIMER Portions of this document may be illegible in electronic image products. Images are produced from the best available original document.

SSCL-Preprint-543 Constant of Motion for a One-Dimensional and nth-order Autonomous System and Its Relation to the Lagrangian and Hamiltonian* G. Lopez Superconducting Super Collider Laboratory^ 2550 Beckleymeade Avenue Dallas, Texas 75237 December 1993 * Submitted to Communications in Mathematical Physics. + Operated by the Universities Research Association, Inc., for the U.S. Department of Energy under Contract No. DE-AC35-89ER40486. DISTRIBUTION OF THIS DOCUMENT IS UNUMITH MASTER

Abstract A constant of motion is defined for a one-dimensional and nth-differential-order autonomous system. A generalization of the Legendre transformation is given that allows one to obtain a relation among the constant of motion, the Lagrangian, and the Hamiltonian. The approach is used to obtain the constant of motion associated with the nonrelativistic third-differential-order Abraham-Lorentz radiation damping equation. I. INTRODUCTION Finding the constant of motion of a dynamical system is one of the best approaches to understanding some of the system's local and global characteristics [1]. In addition, the constants of motion are of considerable importance in physics [2,3], since they are closely connected with the concept of ''energy" and invariants of motion [4]. The constant of motion approach has been used to study the relationship among this constant, the Lagrangian, and the Hamiltonian for one-dimensional and second-differential-order dynamical systems [5,6]. In this paper, the study of the constant of motion and its relation with the Lagrangian and the Hamiltonian is extended to one-dimensional and nth-differential-order autonomous systems. The approach is applied to the nonrelativistic third-differential-order Abraham-Lorentz radiation damping problem. II. CONSTANT OF MOTION A one-dimensional autonomous system of nth-order of differentiation is described by the equation ^ = F(,,x«1 ),...,^- 1 )), (1) where t represents the time, x represents the position, x^*\ i = l,...,n 1 is the ithdifferentiation with respect to time (x^ = d l x/dt l ), and the function F does not depend explicitly on time. Using the definitions & = xw, i = l,...,n (2a) 2

and Zo = x. (26) Eq. (1) can be written as the following dynamical system: and A constant of motion of this system is a function, <*&- i = &, z = 1 n-1 (3a) dt ^ = F(x.6,...,^-i). (36) tf = A-(x, i,..., n -i), (4) such that its total differentiation with respect to time is null (dk/dt 0), i.e., the function K must be a solution of the following partial differential equation: E. uk _, ok,_. 6«r - + Fs,6,...,e B _ 1 w = o. (5) This equation can be solved by the characteristics method [7], where the equations for the characteristics are given by Finding n 1 characteristics curves, C, = C,(x, i,..., n _i), i l,...,n, the general solution of Eq. (6) is given by any arbitrary function of the characteristics, K = K(C u...,c n ). (7) 3

III. RELATION WITH THE LAGRANGIAN Define the operator ft 1 " as and the generalized Legendre transformation as Q, = ^ ( _ i r -. e ( g ) i=l where is the Lagrangian associated with Eq. (1) and G(x, i,..., n -i) is any arbitrary function. Applying the total time differentiation operator, to Eq. (9) and making some rearrangements, the result is ( n-l 14^v <f" dl dl *' rft"d^ ax 6 (-iy +1 dg, x This expression brings about the relation among the constant of motion, the Legendre transformation, and the Lagrangian. If L is the Lagrangian of the system (1) (meaning that the expression inside the braces of Eq. (11) is null), the function G must be a constant of motion of the system (G = A'). On the other hand, if G is a constant of motion, then L represents the Lagrangian of the system ( i ^ 0). IV. RELATION WITH THE HAMILTONIAN Defining the generalized momenta u M, y. = 1,..., n 1, as dl -; (!. i,..., n _i) =,(JL = l,...,n-l, (12) the Legendre transformation (9) becomes ^EB^H^'-I, (i=i t=i as)

where uijj 1 represents the total (/J. i)th differentiation with respect to the time of the new variable u: M, and function H is given by H{x.*)\,...,*) n -i) = K[x,^i{x.~x!i...,uJ n -\),...,f _i(:r,u;i,...,il> n ~\)j (14) Of course, the jacobian of the transformation (tow must be different from zero. Making the differential variation in Eq. (13) and rearranging terms, the generalized Hamiltonian equations follow: dh dx 71 1 II M=lt=l r 9 ^ l ),,(/*-!) ^«^ ~ d x ~ + ^ ~d~x <9x (15a) and ^-EIB- 1 '-' (M-i) dum WA a '"a'^a.a = l...,n-1. (156) V. SPECIAL CASES Clearly, Eqs. (9), (11), (13), and (15) are reduced themselves to the well-known expressions for a second-order-differential autonomous system (n = 2). For a third-orderdifferential dynamical system, the Legendre transformation (9) represents a second-order partial differential equation of variable coefficients: / d\dl dl. r, [ a - V Jt)^JrV^- L = K^ (16a) where v = i is the velocity, a = (2 is the acceleration, and the total time differentiation operator is given by The solution of Eq. (16) is not trivial and depends on the dynamics, which is determined by the equation z ( 3 ) = F(x,v,a). (17) This fact contrasts with the second-order case of Reference 6, where all the dynamics is kept in the constant of motion K. 5

VI. ABRAHAM-LORENTZ RADIATION DAMPING The third-order Abraham-Lorentz equation for the electromagnetic radiation of a charged particle [8] that is moving freely in one-dimensional space can be written as x(2) _ TX(3) = 0 _ ( 18a j where the parameter r characterizes the dissipation of energy and is given by 2 e 2 3 mc* e being the charge of the particle, m its mass, and c the speed of light. The dynamical system associated with Eq. (18) is expressed by the equations z ( 1 ) = v, (19a) and x^ = a, (196) x (3) = a/r, (19c) where v represents the velocity, and a is the acceleration of the charge. The constant of motion satisfies the following partial differential equation: dk OK adk n ox dv T da It is not difficult to find from Eqs. (6) and (19) the following two characteristic curves: C\ = v ra (21a) and Therefore, the general solution of Eq. (20) is C 2 = r \v + (v -ra)log(ra) - r. (216) K(x.v,a) = G(C 1,C 2 ), (22) 6

where G is an arbitrary function of the characteristic curves. Selecting this function of the form G{Ci.C 2 ) = \mcl + Tg(C 2 ). (23) one is allowed to recover the usual expression for the nonrelativistic energy when the damping factor r goes to zero. (Clearly, the arbitrary function g must remain finite for r equal to zero.) Thus, the constant of motion for this system can be written as K(x. v. a) = -mv~ rmva -r~ma~ 4- rg(rv + r(v ra) log(ra) x). (24) The function g is arbitrary and may be determined from experimental analysis on the radiating system. Substituting this expression in Eq. (16) and solving the resulting partial differential equation, the Lagrangian should be obtained. However, to do this requires much more elaborate and complex analysis that will not be performed here. VII. CONCLUSIONS A constant of motion was defined for one-dimensional and nth-differential-order autonomous systems. By generalizing the Legendre transformation on these systems, a relation is obtained between the Lagrangian and the constant of motion. The relation between the constant of motion and the Hamiltonian also follows from the generalized Legendre transformation. Finally, for the nonrelativistic Abraham-Lorentz radiation damping system, an explicit constant of motion is deduced using the above approach. i

REFERENCES 1. E. A. Jackson. Perspectives of Nonlinear Dynamics, Chapter 2, Cambridge University Press. 1989. 2. L. D. Landau and E. M. Lifshitz, Mechanics, Pergamon Press, 1960. 3. H. Goldstein. Classical Mechanics. Addison-Wesley Press, 1950. 4. E. D. Courant and H. S. Snyder. Ann. of Phys., 3, 1 (1958). 5. G. Lopez and J. I. Hernandez, Ann. of Phys., 193. No. 1, 1 (1989). 6. G. Lopez, "Hamiltonians and Lagrangians for One-Dimensional Autonomous Systems and Dissipative Systems." Report SSCL-472. June 1991. 7. F. John. Partial Differential Equations. Springer-Verlag, 1974. 8. J. D. Jackson, Classical Electrodynamics, Chapter 17, John Wiley & Sons, 1975. 8