Control of chaos in Hamiltonian systems

Size: px
Start display at page:

Download "Control of chaos in Hamiltonian systems"

Transcription

1 Control of chaos in Hamiltonian systems G. Ciraolo, C. Chandre, R. Lima, M. Vittot Centre de Physique Théorique CNRS, Marseille M. Pettini Osservatorio Astrofisico di Arcetri, Università di Firenze Ph. Ghendrih Association Euratom-CEA, CEA Cadarache Contact: Financial support: Association Euratom-CEA, CNRS, INFM References: J. Phys. A: Math. Gen. 37, 3589 (4) J. Phys. A: Math. Gen. 37, 6337 (4) Phys. Rev. E 69, 5613 (4) Nonlinearity 18 (5) Euratom Fréjus October, 4

2 Our strategy of control 1 original system dx (,) dt = X x t controlled system dx (,) c dt = X x t with X c = X + f chaotic dynamics regular dynamics control term Aim: find a small and apt control f > small modification of the system / great influence on the dynamics

3 Our strategy of control Aim: H = H +εv chaotic Hc = H + εv + f regular f = εv obvious and useless solution tailoring the control term f Requirements on f : small with respect to the perturbation εv > here, we require that f = O( ε ) localized in phase space > accessible region, fewer energy for the control with a certain shape > robustness other requirements?

4 Physical situations > electrostatic turbulence : E B drift motion y x ϕ θ V d x c c y dt y = E B= B B V x where E = V ( xy,, t) > magnetic field lines dψ H = dϕ θ dθ H =+ dϕ ψ H ( ψθ,, ϕ) = H ( ψ ) + εh ( ψ,, θ ϕ) 1 > atomic physics : atoms in external fields, traps > particle accelerators, free electron laser,...?

5 I. Global control of Hamiltonian systems : - method and illustration - application to edge plasma turbulence II. Localised control of Hamiltonian systems - method and illustration - application to magnetic field lines

6 In general Lie formalism In action-angle angle variables: H = H + V H( A, ϕ ) = H( A) + Vk( A) e k n where ( Αϕ, ) T n n i k ϕ LH given by LHV = { H, V} {, } = Poisson bracket L H V H V V H = A ϕ A ϕ n = iω k V k k ( A) e i k ϕ where H ω = A Γ pseudo inverse of L : L H Γ = H L H Γ V Vk( A) = n i ω k k ω k e i k ϕ R projector: R = 1 LH Γ RV V ( A) e = n k ω k= k i k ϕ

7 Global control Proposition 1 : H = H + V + f H + RV c and are canonically conjugate where f n ( 1) n = LΓ V ( nr + 1) V ( n + 1)! n= 1 V ε if is of order is of order f O ( ε ) > f = n= f n where f n is of order ε n H > If ω = is non-resonant and in many other situations, H + RV is integrable A

8 Control of chaotic transport in a model for E B drift in magnetized plasmas y B E( xyt,, ) B x V d x c c y where ( xy,, t) dt y = E B= E= V B B V x V is the Hamiltonian x, y are conjugate variables Algorithm for the computation of the control term f Numerical investigation: - of the effect of the control term (in practice f ) - of its robustness

9 Control of E B E B drift Electrostatic potential known on a spatio-temporal grid. Vxyt (,,) Vxyt= (,, ) 1) Expansion of iωk Vxyt (,,) Vk (, xy) e k t ) Computation of ΓV = ω k Vk (, x y) iωk e i ω k t 3) Computation of f 1 = Γ { V, V } where {, } = x y y x f( x,, y t = ) 4) Higher orders? ( recall: f = f n ). n= f n 1 = Γ n { V, f } n 1

10 Example: Vxyt= (,, ) ak Vxyt (,,) = sin ( π ( nx+ my) +ϕ 3/ nmk ωkt) ( n + m ) mnk,, Γ Vxyt (,, = ) ϕnmk random phases f3( xyt=,, ) f (,, xyt= )

11 The control term is a small modification of the potential Vxyt (,,) Vx (, y, t) + f (,, x y t)

12 Trajectories of particles > Poincaré surface of section a =.8 without control with control f > Projection on the x axis

13 Diffusion of test particles Diffusion coefficient without control D M 1 = lim xi( t) xi() t Mt i= 1 with control f > significant reduction of diffusion by f Probability Distribution Function of step sizes Step size distance between two successive sign reversals of the velocity without control with control f > the control quenches the large steps

14 Robustness of the control : a crucial requirement experimentally accessible H + ε V + f H + ε V + f domain of efficient control H + εv H + εv H + εv + f H + εv + f uncontrolled Hamiltonian controlled Hamiltonian, computed theoretically controlled Hamiltonian, implemented experimentally

15 Robustness 1 > Modification of the amplitude of the control term : f δ f a =.7 H c = H +ε V + δ f without control > optimal control for δ = 1 > increasing its amplitude does not improve control > decreasing the amplitude does still give a good control > reducing energy to control the system (δ =.5 reduction of 3%)

16 Robustness if the potential is known on a coarse grained grid: truncation of the Fourier series of the control term f a =.8 without control control f > efficient control with few Fourier modes: 1 modes reduction of 5% > simplification of the control term > reduction of the energy necessary for the control

17 I. Global control of Hamiltonian systems : - method and illustration - application to edge plasma turbulence II. Localised control of Hamiltonian systems - method and illustration - application to magnetic field lines

18 Localised control : channeling chaos by building barriers Proposition : There exists a control term f of order ε H c ( ϕ) ( A ϕh(, ϕ) ) ( A, ϕ) = H( A, ϕ) + f Ω + εγ has an invariant torus at A = εγ H(, ϕ). It is given by f ( ϕ) = H ( εγ H(, ϕ), ϕ). ϕ ϕ such that > expansion of H around A = A, eg.., A = : Ω( x) Remark: for d.o.f. barrier to diffusion for n d.o.f. effective barriers η η x Advantages: > Explicit expression for the control term: existence and regularity > Explicit expression for the invariant torus which has been created > Persistence of the created torus for arbitrarily large values of ε (provided f V )

19 Algorithm of construction of a localised control term p on an example: Hpxt (,,) = +εvxt (,) 1) translation of the momentum p by ω HpExt (,,, ) = ω p+ E+ε Vxt (, ) + ) computation of Γ V = ( Γ V, Γ V) 1 k1 Γ V ϕ = V k1, k ω k k 1 + k ikx e 1 + kt 1 ( ) (,, t) 3) computation of f = H εγ V, εγ V x ε f = ϕ x k1 k k ωk + k, 1 1 4) equation of the invariant torus p (, x t) = ω ε 1 x V kk 1 k1 ω k + k k, k 1 t kk e t ikx ( + kt) V 1 kk 1 e p ikx ( + kt) 1 A = (, pe) ϕ = (,) xt ω = ( ω,1)

20 Localised control: forced pendulum p Hc = +ε x x t ε cos x cos( x t) ( cos + cos( )) + Ω( p p (,) x t ) ω ω 1 without control ε =.5 with control Poincaré section 1 trajectory control 4% of V on 1% phase space control 4% of V on 1% phase space > uncontrolled case : no barrier for ε.8

21 Localised control: channelling chaos by building barriers Creation of barriers at different locations Poincaré section Poincaré section Creation of set of barriers M p Hc = +ε x i i x ( cos + cos( x t) ) + ε f(,) x t Ω( p p (,) t ) i= 1 > Trapping of particles ε =.5 Poincaré section

22 Localised control: magnetic field lines Magnetic field line dynamics in a toroidal geometry dϑ H =+ dϕ ψ dψ H = dϕ ϑ ϕ θ Nearly axisymmetric case: H( ψϑϕ,, ) = H( ψ ) + εh1( ψ, ϑ, ϕ) regular part: H ( ψ ) = dψ q( ψ) ε =.3 safety factor: 4 q( ψ ) = ψ ψ+ ψ ( )( ) magnetic perturbation: [ m m ϑ ] sin ( ) H m n m d m/ 1( ψϑϕ,, ) = ( 1) ψ cos( ϑ ϕ) m π( m m)

23 Localised control: magnetic field lines around ψ ψ : H ( ψϑϕ,, ) = H ( ψ ) + εh ( ψ, ϑ, ϕ) + ε f ( ϑ, ϕ) Ω( ψ ψ + εγ H ( ψ, θ, ϕ) ) c 1 θ 1 ε =.3 ε =.3 control 5% f V on 7% phase space ε =.3

24 Summary and outlook > Global control of Hamiltonian systems - application to electrostatic turbulence -robustness > Localised control of Hamiltonian systems - application to magnetic field lines - barrier to diffusion > control of chaos in area-preserving maps [cf. Poster] > application on a TWT [cf. Doveil s talk] > extension to more realistic models in fusion plasmas > other applications to atomic physics, particle accelerators, free electron laser..?

25 Robustness 3 > phase error in f ϕ = ϕ + γ ϕ a =.7 err nm nm nm without control control f > 5% reduction of the diffusion coefficient for 5 % error in the phases

Controlling chaotic transport in Hamiltonian systems

Controlling chaotic transport in Hamiltonian systems Controlling chaotic transport in Hamiltonian systems Guido Ciraolo Facoltà di Ingegneria, Università di Firenze via S. Marta, I-50129 Firenze, Italy Cristel Chandre, Ricardo Lima, Michel Vittot CPT-CNRS,

More information

Chaotic Motion in Problem of Dumbell Satellite

Chaotic Motion in Problem of Dumbell Satellite Int. J. Contemp. Math. Sciences, Vol. 6, 2011, no. 7, 299-307 Chaotic Motion in Problem of Dumbell Satellite Ayub Khan Department of Mathematics, Zakir Hussain College University of Delhi, Delhi, India

More information

arxiv:nlin/ v1 [nlin.cd] 17 Dec 2003

arxiv:nlin/ v1 [nlin.cd] 17 Dec 2003 Control of Hamiltonian chaos as a possible tool to control anomalous transport in fusion plasmas Guido Ciraolo Facoltà di Ingegneria, Università di Firenze, via S. Marta, I-529 Firenze, Italy, and I.N.F.M.

More information

Gyrokinetic simulations of magnetic fusion plasmas

Gyrokinetic simulations of magnetic fusion plasmas Gyrokinetic simulations of magnetic fusion plasmas Tutorial 2 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr

More information

EUROFUSION WPJET1-PR(16) CG Albert et al.

EUROFUSION WPJET1-PR(16) CG Albert et al. EUROFUSION WPJET1-PR(16) 15331 CG Albert et al. Hamiltonian approach for evaluation of toroidal torque from finite amplitude non-axisymmetric perturbations of a tokamak magnetic field in resonant transport

More information

Internal transport barriers in some Hamiltonian systems modeling the magnetic lines dynamics in tokamak

Internal transport barriers in some Hamiltonian systems modeling the magnetic lines dynamics in tokamak Journal of Physics: Conference Series Internal transport barriers in some Hamiltonian systems modeling the magnetic lines dynamics in tokamak To cite this article: D Constantinescu et al 2005 J. Phys.:

More information

Channeling chaotic transport in a wave-particle experiment

Channeling chaotic transport in a wave-particle experiment Channeling chaotic transport in a wave-particle experiment Alessandro Macor, Fabrice Doveil, Cristel Chandre, Guido Ciraolo, Ricardo Lima, Michael Vittot To cite this version: Alessandro Macor, Fabrice

More information

A Lie Algebra version of the Classical or Quantum Hamiltonian Perturbation Theory and Hamiltonian Control with Examples in Plasma Physics

A Lie Algebra version of the Classical or Quantum Hamiltonian Perturbation Theory and Hamiltonian Control with Examples in Plasma Physics A Lie Algebra version of the Classical or Quantum Hamiltonian Perturbation Theory and Hamiltonian Control with Examples in Plasma Physics Michel VITTOT Center for Theoretical Physics (CPT), CNRS - Luminy

More information

Hamiltonian formulation of reduced Vlasov-Maxwell equations

Hamiltonian formulation of reduced Vlasov-Maxwell equations Hamiltonian formulation of reduced Vlasov-Maxell equations Cristel CHANDRE Centre de Physique Théorique CNRS, Marseille, France Contact: chandre@cpt.univ-mrs.fr importance of stability vs instability in

More information

Enhancement of particle trapping in the wave-particle interaction

Enhancement of particle trapping in the wave-particle interaction Communications in onlinear Science and umerical Simulation 13 (2008) 660 665 www.elsevier.com/locate/cnsns Enhancement of particle trapping in the wave-particle interaction R. Bachelard a, *, A. Antoniazzi

More information

= 0. = q i., q i = E

= 0. = q i., q i = E Summary of the Above Newton s second law: d 2 r dt 2 = Φ( r) Complicated vector arithmetic & coordinate system dependence Lagrangian Formalism: L q i d dt ( L q i ) = 0 n second-order differential equations

More information

Chaotic motion. Phys 420/580 Lecture 10

Chaotic motion. Phys 420/580 Lecture 10 Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t

More information

Oscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution

Oscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ small angle approximation θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic

More information

Oscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution

Oscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic angular frequency

More information

Hamiltonian aspects of fluid dynamics

Hamiltonian aspects of fluid dynamics Hamiltonian aspects of fluid dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS 01/29/08, 01/31/08 Joris Vankerschaver (CDS) Hamiltonian aspects of fluid dynamics 01/29/08, 01/31/08 1 / 34 Outline

More information

Physics 106b: Lecture 7 25 January, 2018

Physics 106b: Lecture 7 25 January, 2018 Physics 106b: Lecture 7 25 January, 2018 Hamiltonian Chaos: Introduction Integrable Systems We start with systems that do not exhibit chaos, but instead have simple periodic motion (like the SHO) with

More information

An introduction to Birkhoff normal form

An introduction to Birkhoff normal form An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an

More information

CPT proposals for Master internships

CPT proposals for Master internships CPT proposals for Master internships Academic year 2018-2019 Warning(!!): Due to a limited funding, the CPT director will fix the number of financial rewardings up to a max of 6 internships. (This warning

More information

Theory for Neoclassical Toroidal Plasma Viscosity in a Toroidally Symmetric Torus. K. C. Shaing

Theory for Neoclassical Toroidal Plasma Viscosity in a Toroidally Symmetric Torus. K. C. Shaing Theory for Neoclassical Toroidal Plasma Viscosity in a Toroidally Symmetric Torus K. C. Shaing Plasma and Space Science Center, and ISAPS, National Cheng Kung University, Tainan, Taiwan 70101, Republic

More information

Hamiltonian and Non-Hamiltonian Reductions of Charged Particle Dynamics: Diffusion and Self-Organization

Hamiltonian and Non-Hamiltonian Reductions of Charged Particle Dynamics: Diffusion and Self-Organization NNP2017 11 th July 2017 Lawrence University Hamiltonian and Non-Hamiltonian Reductions of Charged Particle Dynamics: Diffusion and Self-Organization N. Sato and Z. Yoshida Graduate School of Frontier Sciences

More information

Chaotic motion. Phys 750 Lecture 9

Chaotic motion. Phys 750 Lecture 9 Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to

More information

Overview of Gyrokinetic Theory & Properties of ITG/TEM Instabilities

Overview of Gyrokinetic Theory & Properties of ITG/TEM Instabilities Overview of Gyrokinetic Theory & Properties of ITG/TEM Instabilities G. W. Hammett Princeton Plasma Physics Lab (PPPL) http://w3.pppl.gov/~hammett AST559: Plasma & Fluid Turbulence Dec. 5, 2011 (based

More information

Particle interactions with spatially localized wavepackets

Particle 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 information

Beams and magnetized plasmas

Beams and magnetized plasmas Beams and magnetized plasmas 1 Jean-Pierre BOEUF LAboratoire PLAsma et Conversion d Energie LAPLACE/ CNRS, Université Paul SABATIER, TOULOUSE Beams and magnetized plasmas 2 Outline Ion acceleration and

More information

Assignment 8. [η j, η k ] = J jk

Assignment 8. [η j, η k ] = J jk Assignment 8 Goldstein 9.8 Prove directly that the transformation is canonical and find a generating function. Q 1 = q 1, P 1 = p 1 p Q = p, P = q 1 q We can establish that the transformation is canonical

More information

arxiv: v1 [physics.plasm-ph] 15 Mar 2011

arxiv: v1 [physics.plasm-ph] 15 Mar 2011 Gyromap for a two-dimensional Hamiltonian fluid model derived from raginskii s closure for magnetized plasmas arxiv:1103900v1 physicsplasm-ph 15 Mar 011 O Izacard, 1 C Chandre, 1 E Tassi, 1 and G Ciraolo

More information

Chapter 29. Quantum Chaos

Chapter 29. Quantum Chaos Chapter 29 Quantum Chaos What happens to a Hamiltonian system that for classical mechanics is chaotic when we include a nonzero h? There is no problem in principle to answering this question: given a classical

More information

Mechanics Physics 151

Mechanics Physics 151 Mechanics Physics 151 Lecture 4 Continuous Systems and Fields (Chapter 13) What We Did Last Time Built Lagrangian formalism for continuous system Lagrangian L Lagrange s equation = L dxdydz Derived simple

More information

Oh, That. P. J. Morrison. Department of Physics and Institute for Fusion Studies The University of Texas at Austin

Oh, That. P. J. Morrison. Department of Physics and Institute for Fusion Studies The University of Texas at Austin Oh, That. P. J. Morrison Department of Physics and Institute for Fusion Studies The University of Texas at Austin morrison@physics.utexas.edu Overview 1. BGK Modes 2. Inverse Scattering 3. Torus Breakup

More information

Lectures on basic plasma physics: Hamiltonian mechanics of charged particle motion

Lectures on basic plasma physics: Hamiltonian mechanics of charged particle motion Lectures on basic plasma physics: Hamiltonian mechanics of charged particle motion Department of applied physics, Aalto University March 8, 2016 Hamiltonian versus Newtonian mechanics Newtonian mechanics:

More information

What place for mathematicians in plasma physics

What place for mathematicians in plasma physics What place for mathematicians in plasma physics Eric Sonnendrücker IRMA Université Louis Pasteur, Strasbourg projet CALVI INRIA Nancy Grand Est 15-19 September 2008 Eric Sonnendrücker (U. Strasbourg) Math

More information

Hamiltonian Dynamics

Hamiltonian Dynamics Hamiltonian Dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS Feb. 10, 2009 Joris Vankerschaver (CDS) Hamiltonian Dynamics Feb. 10, 2009 1 / 31 Outline 1. Introductory concepts; 2. Poisson brackets;

More information

Intermittent onset of turbulence and control of extreme events

Intermittent onset of turbulence and control of extreme events Intermittent onset of turbulence and control of extreme events Sergio Roberto Lopes, Paulo P. Galúzio Departamento de Física UFPR São Paulo 03 de Abril 2014 Lopes, S. R. (Física UFPR) Intermittent onset

More information

Ion plateau transport near the tokamak magnetic axis. Abstract

Ion plateau transport near the tokamak magnetic axis. Abstract Ion plateau transport near the tokamak magnetic axis K.C. Shaing and R.D. Hazeltine Institute for Fusion Studies, The University of Texas at Austin, Austin, Texas 78712 (December 2, 1997) Abstract Conventional

More information

Bursty Transport in Tokamaks with Internal Transport Barriers

Bursty Transport in Tokamaks with Internal Transport Barriers Bursty Transport in Tokamaks with Internal Transport Barriers S. Benkadda 1), O. Agullo 1), P. Beyer 1), N. Bian 1), P. H. Diamond 3), C. Figarella 1), X. Garbet 2), P. Ghendrih 2), V. Grandgirard 1),

More information

M2A2 Problem Sheet 3 - Hamiltonian Mechanics

M2A2 Problem Sheet 3 - Hamiltonian Mechanics MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,

More information

On long-term dynamics of slow-fast systems with passages through resonances

On long-term dynamics of slow-fast systems with passages through resonances On long-term dynamics of slow-fast systems with passages through resonances Anatoly Neishtadt Loughborough University, UK Space Research Institute, Moscow, Russia Based on joint works with A.Artemyev,

More information

Geometric Gyrokinetic Theory and its Applications to Large-Scale Simulations of Magnetized Plasmas

Geometric Gyrokinetic Theory and its Applications to Large-Scale Simulations of Magnetized Plasmas Geometric Gyrokinetic Theory and its Applications to Large-Scale Simulations of Magnetized Plasmas Hong Qin Princeton Plasma Physics Laboratory, Princeton University CEA-EDF-INRIA School -- Numerical models

More information

Physics 486 Discussion 5 Piecewise Potentials

Physics 486 Discussion 5 Piecewise Potentials Physics 486 Discussion 5 Piecewise Potentials Problem 1 : Infinite Potential Well Checkpoints 1 Consider the infinite well potential V(x) = 0 for 0 < x < 1 elsewhere. (a) First, think classically. Potential

More information

A statistical study of gyro-averaging effects in a reduced model of drift-wave transport. Abstract

A statistical study of gyro-averaging effects in a reduced model of drift-wave transport. Abstract A statistical study of gyro-averaging effects in a reduced model of drift-wave transport J. D. da Fonseca, 1, D. del-castillo-negrete, 2, I. M. Sokolov, 3, and I. L. Caldas 1, 1 Physics Institute, University

More information

Plasma turbulence measured by fast sweep reflectometry on TORE SUPRA

Plasma turbulence measured by fast sweep reflectometry on TORE SUPRA Plasma turbulence measured by fast sweep reflectometry on TORE SUPRA F. Clairet &, L. Vermare &, S. Heuraux, G. Leclert # & Association Euratom-CEA sur la fusion, DSM/DRFC/SCCP C.E. Cadarache, 8 Saint-Paul-lès-Durance,

More information

Chapter 2 Heisenberg s Matrix Mechanics

Chapter 2 Heisenberg s Matrix Mechanics Chapter 2 Heisenberg s Matrix Mechanics Abstract The quantum selection rule and its generalizations are capable of predicting energies of the stationary orbits; however they should be obtained in a more

More information

Electron Cyclotron Current Drive in Spherical Tokamaks with Application to ITER

Electron Cyclotron Current Drive in Spherical Tokamaks with Application to ITER 1 TH/P3-16 Electron Cyclotron Current Drive in Spherical Tokamaks with Application to ITER A.K. Ram 1), Y. Kominis 2), K. Hizanidis 2), J. Decker 3) 1) Plasma Science & Fusion Center, M.I.T., Cambridge,

More information

Scattering of ECRF waves by edge density fluctuations and blobs

Scattering of ECRF waves by edge density fluctuations and blobs PSFC/JA-14-7 Scattering of ECRF waves by edge density fluctuations and blobs A. K. Ram and K. Hizanidis a June 2014 Plasma Science and Fusion Center, Massachusetts Institute of Technology Cambridge, MA

More information

Non-linear dynamics! in particle accelerators!

Non-linear dynamics! in particle accelerators! Non-linear dynamics! in particle accelerators! Yannis PAPAPHILIPPOU Accelerator and Beam Physics group Beams Department CERN" Cockroft Institute Lecture Courses Daresbury, UK 16-19 September 2013 1! Contents

More information

GEOMETRIC QUANTIZATION

GEOMETRIC QUANTIZATION GEOMETRIC QUANTIZATION 1. The basic idea The setting of the Hamiltonian version of classical (Newtonian) mechanics is the phase space (position and momentum), which is a symplectic manifold. The typical

More information

Notes for Lecture 10

Notes for Lecture 10 February 2, 26 Notes for Lecture Introduction to grid-based methods for solving Poisson and Laplace Equations Finite difference methods The basis for most grid-based finite difference methods is the Taylor

More information

A path integral approach to the Langevin equation

A path integral approach to the Langevin equation A path integral approach to the Langevin equation - Ashok Das Reference: A path integral approach to the Langevin equation, A. Das, S. Panda and J. R. L. Santos, arxiv:1411.0256 (to be published in Int.

More information

Lectures on Dynamical Systems. Anatoly Neishtadt

Lectures on Dynamical Systems. Anatoly Neishtadt Lectures on Dynamical Systems Anatoly Neishtadt Lectures for Mathematics Access Grid Instruction and Collaboration (MAGIC) consortium, Loughborough University, 2007 Part 3 LECTURE 14 NORMAL FORMS Resonances

More information

Some Collision solutions of the rectilinear periodically forced Kepler problem

Some Collision solutions of the rectilinear periodically forced Kepler problem Advanced Nonlinear Studies 1 (2001), xxx xxx Some Collision solutions of the rectilinear periodically forced Kepler problem Lei Zhao Johann Bernoulli Institute for Mathematics and Computer Science University

More information

DOPPLER RESONANCE EFFECT ON ROTATIONAL DRIVE BY ION CYCLOTRON MINORITY HEATING

DOPPLER RESONANCE EFFECT ON ROTATIONAL DRIVE BY ION CYCLOTRON MINORITY HEATING DOPPLER RESONANCE EFFECT ON ROTATIONAL DRIVE BY ION CYCLOTRON MINORITY HEATING V.S. Chan, S.C. Chiu, Y.A. Omelchenko General Atomics, San Diego, CA, U.S.A. 43rd Annual APS Division of Plasma Physics Meeting

More information

Stability and instability of solitons in inhomogeneous media

Stability and instability of solitons in inhomogeneous media Stability and instability of solitons in inhomogeneous media Yonatan Sivan, Tel Aviv University, Israel now at Purdue University, USA G. Fibich, Tel Aviv University, Israel M. Weinstein, Columbia University,

More information

Александр Баланов

Александр Баланов a.balanov@lboro.ac.uk Александр Баланов a.balanov@lboro.ac.uk Collaborators M.T. Greenaway (University of Nottingham, UK) T.M. Fromhold (University of Nottingham, UK) A.O. Selskii (Saratov State University,

More information

The Linear Theory of Tearing Modes in periodic, cyindrical plasmas. Cary Forest University of Wisconsin

The Linear Theory of Tearing Modes in periodic, cyindrical plasmas. Cary Forest University of Wisconsin The Linear Theory of Tearing Modes in periodic, cyindrical plasmas Cary Forest University of Wisconsin 1 Resistive MHD E + v B = ηj (no energy principle) Role of resistivity No frozen flux, B can tear

More information

Non-Twist Field Line Mappings for Tokamaks with Reversed Magnetic Shear

Non-Twist Field Line Mappings for Tokamaks with Reversed Magnetic Shear Brazilian Journal of Physics, vol. 34, no. 4B, December, 2004 1759 Non-Twist Field Line Mappings for Tokamaks with Reversed Magnetic Shear M. Roberto 1, E. C. da Silva 2, I. L. Caldas 2, and R. L. Viana

More information

TURBULENT TRANSPORT THEORY

TURBULENT TRANSPORT THEORY ASDEX Upgrade Max-Planck-Institut für Plasmaphysik TURBULENT TRANSPORT THEORY C. Angioni GYRO, J. Candy and R.E. Waltz, GA The problem of Transport Transport is the physics subject which studies the physical

More information

Acta Polytechnica 53(4): , 2013 Czech Technical University in Prague, 2013 available online at

Acta Polytechnica 53(4): , 2013 Czech Technical University in Prague, 2013 available online at Acta Polytechnica 53(4):359 364, 2013 Czech Technical University in Prague, 2013 available online at http://ctn.cvut.cz/ap/ PARTICLE DYNAMICS IN THE RMP ERGODIC LAYER UNDER THE INFLUENCE OF EDGE PLASMA

More information

Action-Angle Variables and KAM-Theory in General Relativity

Action-Angle Variables and KAM-Theory in General Relativity Action-Angle Variables and KAM-Theory in General Relativity Daniela Kunst, Volker Perlick, Claus Lämmerzahl Center of Space Technology and Microgravity University of Bremen, Germany Workshop in Oldenburg

More information

Flow, current, & electric field in omnigenous stellarators

Flow, current, & electric field in omnigenous stellarators Flow, current, & electric field in omnigenous stellarators Supported by U.S. D.o.E. Matt Landreman with Peter J Catto MIT Plasma Science & Fusion Center Oral 2O4 Sherwood Fusion Theory Meeting Tuesday

More information

arxiv: v7 [quant-ph] 22 Aug 2017

arxiv: v7 [quant-ph] 22 Aug 2017 Quantum Mechanics with a non-zero quantum correlation time Jean-Philippe Bouchaud 1 1 Capital Fund Management, rue de l Université, 75007 Paris, France. (Dated: October 8, 018) arxiv:170.00771v7 [quant-ph]

More information

THE POSSIBILITY OF PRECISE MEASUREMENT OF ABSOLUTE ENERGY OF THE ELECTRON BEAM BY MEANS OF RESONANCE ABSORPTION METHOD R.A.

THE POSSIBILITY OF PRECISE MEASUREMENT OF ABSOLUTE ENERGY OF THE ELECTRON BEAM BY MEANS OF RESONANCE ABSORPTION METHOD R.A. THE POSSIBILITY OF PRECISE MEASUREMENT OF ABSOLUTE ENERGY OF THE ELECTRON BEAM BY MEANS OF RESONANCE ABSORPTION METHOD R.A. Melikian Yerevan Physics Institute, Yerevan, Armenia Abstract In this report

More information

Lecture 8 Phase Space, Part 2. 1 Surfaces of section. MATH-GA Mechanics

Lecture 8 Phase Space, Part 2. 1 Surfaces of section. MATH-GA Mechanics Lecture 8 Phase Space, Part 2 MATH-GA 2710.001 Mechanics 1 Surfaces of section Thus far, we have highlighted the value of phase portraits, and seen that valuable information can be extracted by looking

More information

Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )

Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University

More information

BACKGROUND IN SYMPLECTIC GEOMETRY

BACKGROUND IN SYMPLECTIC GEOMETRY BACKGROUND IN SYMPLECTIC GEOMETRY NILAY KUMAR Today I want to introduce some of the symplectic structure underlying classical mechanics. The key idea is actually quite old and in its various formulations

More information

Gyrokinetic simulations of magnetic fusion plasmas

Gyrokinetic simulations of magnetic fusion plasmas Gyrokinetic simulations of magnetic fusion plasmas Tutorial 1 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr

More information

Nonlinear Dynamics: Synchronisation

Nonlinear Dynamics: Synchronisation Nonlinear Dynamics: Synchronisation Bristol Centre for Complexity Sciences Ian Ross BRIDGE, School of Geographical Sciences, University of Bristol October 19, 2007 1 / 16 I: Introduction 2 / 16 I: Fireflies

More information

2 Canonical quantization

2 Canonical quantization Phys540.nb 7 Canonical quantization.1. Lagrangian mechanics and canonical quantization Q: How do we quantize a general system?.1.1.lagrangian Lagrangian mechanics is a reformulation of classical mechanics.

More information

Recursion Systems and Recursion Operators for the Soliton Equations Related to Rational Linear Problem with Reductions

Recursion Systems and Recursion Operators for the Soliton Equations Related to Rational Linear Problem with Reductions GMV The s Systems and for the Soliton Equations Related to Rational Linear Problem with Reductions Department of Mathematics & Applied Mathematics University of Cape Town XIV th International Conference

More information

Heating and current drive: Radio Frequency

Heating and current drive: Radio Frequency Heating and current drive: Radio Frequency Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 13 th February 2012 Dr Ben Dudson Magnetic Confinement Fusion (1 of 26)

More information

KAM for NLS with harmonic potential

KAM for NLS with harmonic potential Université de Nantes 3rd Meeting of the GDR Quantum Dynamics MAPMO, Orléans, 2-4 February 2011. (Joint work with Benoît Grébert) Introduction The equation : We consider the nonlinear Schrödinger equation

More information

Page 684. Lecture 40: Coordinate Transformations: Time Transformations Date Revised: 2009/02/02 Date Given: 2009/02/02

Page 684. Lecture 40: Coordinate Transformations: Time Transformations Date Revised: 2009/02/02 Date Given: 2009/02/02 Page 684 Lecture 40: Coordinate Transformations: Time Transformations Date Revised: 2009/02/02 Date Given: 2009/02/02 Time Transformations Section 12.5 Symmetries: Time Transformations Page 685 Time Translation

More information

Impact of diverted geometry on turbulence and transport barrier formation in 3D global simulations of tokamak edge plasma

Impact of diverted geometry on turbulence and transport barrier formation in 3D global simulations of tokamak edge plasma 1 Impact of diverted geometry on turbulence and transport barrier formation in 3D global simulations of tokamak edge plasma D. Galassi, P. Tamain, H. Bufferand, C. Baudoin, G. Ciraolo, N. Fedorczak, Ph.

More information

Analysis of Nonlinear Dynamics by Square Matrix Method

Analysis of Nonlinear Dynamics by Square Matrix Method Analysis of Nonlinear Dynamics by Square Matrix Method Li Hua Yu Brookhaven National Laboratory NOCE, Arcidosso, Sep. 2017 Write one turn map of Taylor expansion as square matrix Simplest example of nonlinear

More information

Modeling of Transport Barrier Based on Drift Alfvén Ballooning Mode Transport Model

Modeling of Transport Barrier Based on Drift Alfvén Ballooning Mode Transport Model 9th IAEA TM on H-mode Physics and Transport Barriers Catamaran Resort Hotel, San Diego 3-9-5 Modeling of Transport Barrier Based on Drift Alfvén Ballooning Mode Transport Model A. Fukuyama, M. Uchida and

More information

Lecture 11 : Overview

Lecture 11 : Overview Lecture 11 : Overview Error in Assignment 3 : In Eq. 1, Hamiltonian should be H = p2 r 2m + p2 ϕ 2mr + (p z ea z ) 2 2 2m + eφ (1) Error in lecture 10, slide 7, Eq. (21). Should be S(q, α, t) m Q = β =

More information

Dynamics of charged particles in spatially chaotic magnetic fields

Dynamics of charged particles in spatially chaotic magnetic fields PSFC/JA-1-38 Dynamics of charged particles in spatially chaotic magnetic fields Abhay K. Ram and Brahmananda Dasgupta a October 21 Plasma Science and Fusion Center, Massachusetts Institute of Technology

More information

Waves and the Schroedinger Equation

Waves and the Schroedinger Equation Waves and the Schroedinger Equation 5 april 010 1 The Wave Equation We have seen from previous discussions that the wave-particle duality of matter requires we describe entities through some wave-form

More information

Notes on x-ray scattering - M. Le Tacon, B. Keimer (06/2015)

Notes on x-ray scattering - M. Le Tacon, B. Keimer (06/2015) Notes on x-ray scattering - M. Le Tacon, B. Keimer (06/2015) Interaction of x-ray with matter: - Photoelectric absorption - Elastic (coherent) scattering (Thomson Scattering) - Inelastic (incoherent) scattering

More information

SMR/ Summer College on Plasma Physics. 30 July - 24 August, Introduction to Magnetic Island Theory.

SMR/ Summer College on Plasma Physics. 30 July - 24 August, Introduction to Magnetic Island Theory. SMR/1856-1 2007 Summer College on Plasma Physics 30 July - 24 August, 2007 Introduction to Magnetic Island Theory. R. Fitzpatrick Inst. for Fusion Studies University of Texas at Austin USA Introduction

More information

1 (2n)! (-1)n (θ) 2n

1 (2n)! (-1)n (θ) 2n Complex Numbers and Algebra The real numbers are complete for the operations addition, subtraction, multiplication, and division, or more suggestively, for the operations of addition and multiplication

More information

Chaos in Hamiltonian systems

Chaos in Hamiltonian systems Chaos in Hamiltonian systems Teemu Laakso April 26, 2013 Course material: Chapter 7 from Ott 1993/2002, Chaos in Dynamical Systems, Cambridge http://matriisi.ee.tut.fi/courses/mat-35006 Useful reading:

More information

Sample Solutions of Assignment 9 for MAT3270B

Sample Solutions of Assignment 9 for MAT3270B Sample Solutions of Assignment 9 for MAT370B. For the following ODEs, find the eigenvalues and eigenvectors, and classify the critical point 0,0 type and determine whether it is stable, asymptotically

More information

Calculation of alpha particle redistribution in sawteeth using experimentally reconstructed displacement eigenfunctions

Calculation of alpha particle redistribution in sawteeth using experimentally reconstructed displacement eigenfunctions Calculation of alpha particle redistribution in sawteeth using experimentally reconstructed displacement eigenfunctions R. Farengo, H. E. Ferrari,2, M.-C. Firpo 3, P. L. Garcia-Martinez 2,3, A. F. Lifschitz

More information

Orbits, Integrals, and Chaos

Orbits, Integrals, and Chaos Chapter 7 Orbits, Integrals, and Chaos In n space dimensions, some orbits can be formally decomposed into n independent periodic motions. These are the regular orbits; they may be represented as winding

More information

Effect of ideal kink instabilities on particle redistribution

Effect of ideal kink instabilities on particle redistribution Effect of ideal kink instabilities on particle redistribution H. E. Ferrari1,2,R. Farengo1, P. L. Garcia-Martinez2, M.-C. Firpo3, A. F. Lifschitz4 1 Comisión Nacional de Energía Atómica, Centro Atomico

More information

CHAPTER 8 The Quantum Theory of Motion

CHAPTER 8 The Quantum Theory of Motion I. Translational motion. CHAPTER 8 The Quantum Theory of Motion A. Single particle in free space, 1-D. 1. Schrodinger eqn H ψ = Eψ! 2 2m d 2 dx 2 ψ = Eψ ; no boundary conditions 2. General solution: ψ

More information

FAST SYMPLECTIC MAPPING, QUASI-INVARIANTS, AND LONG-TERM STABILITY IN THE LHC

FAST SYMPLECTIC MAPPING, QUASI-INVARIANTS, AND LONG-TERM STABILITY IN THE LHC SLAC-PUB-96-7091 January 1996 FAST SYMPLECTIC MAPPING, QUASI-INVARIANTS, AND LONG-TERM STABILITY IN THE LHC ROBERT L. WARNOCK and J. SCOTT BERG Stanford Linear Accelerator Center, Stanford University Stanford,

More information

Theory of Adiabatic Invariants A SOCRATES Lecture Course at the Physics Department, University of Marburg, Germany, February 2004

Theory of Adiabatic Invariants A SOCRATES Lecture Course at the Physics Department, University of Marburg, Germany, February 2004 Preprint CAMTP/03-8 August 2003 Theory of Adiabatic Invariants A SOCRATES Lecture Course at the Physics Department, University of Marburg, Germany, February 2004 Marko Robnik CAMTP - Center for Applied

More information

Numerical methods for eigenvalue problems in the description of drift instabilities in the plasma edge

Numerical methods for eigenvalue problems in the description of drift instabilities in the plasma edge Numerical methods for eigenvalue problems in the description of drift instabilities in the plasma edge Dominik Löchel Supervisors: M. Hochbruck and M. Tokar Graduate school Dynamic of hot plasmas Mathematisches

More information

INTERACTION OF AN EXTERNAL ROTATING MAGNETIC FIELD WITH THE PLASMA TEARING MODE SURROUNDED BY A RESISTIVE WALL

INTERACTION OF AN EXTERNAL ROTATING MAGNETIC FIELD WITH THE PLASMA TEARING MODE SURROUNDED BY A RESISTIVE WALL INTERACTION OF AN EXTERNAL ROTATING MAGNETIC FIELD WITH THE PLASMA TEARING MODE SURROUNDED BY A RESISTIVE WALL S.C. GUO* and M.S. CHU GENERAL ATOMICS *Permanent Address: Consorzio RFX, Padova, Italy **The

More information

Eigenmodes for coupled harmonic vibrations. Algebraic Method for Harmonic Oscillator.

Eigenmodes for coupled harmonic vibrations. Algebraic Method for Harmonic Oscillator. PHYS208 spring 2008 Eigenmodes for coupled harmonic vibrations. Algebraic Method for Harmonic Oscillator. 07.02.2008 Adapted from the text Light - Atom Interaction PHYS261 autumn 2007 Go to list of topics

More information

Linear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems

Linear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems Linear ODEs p. 1 Linear ODEs Existence of solutions to linear IVPs Resolvent matrix Autonomous linear systems Linear ODEs Definition (Linear ODE) A linear ODE is a differential equation taking the form

More information

9.3: Separable Equations

9.3: Separable Equations 9.3: Separable Equations An equation is separable if one can move terms so that each side of the equation only contains 1 variable. Consider the 1st order equation = F (x, y). dx When F (x, y) = f (x)g(y),

More information

Two models for the parametric forcing of a nonlinear oscillator

Two models for the parametric forcing of a nonlinear oscillator Nonlinear Dyn (007) 50:147 160 DOI 10.1007/s11071-006-9148-3 ORIGINAL ARTICLE Two models for the parametric forcing of a nonlinear oscillator Nazha Abouhazim Mohamed Belhaq Richard H. Rand Received: 3

More information

arxiv: v1 [nlin.si] 11 Mar 2015

arxiv: v1 [nlin.si] 11 Mar 2015 Expanding the simple pendulum s rotation solution in action-angle variables Martin Lara a,1,, Sebastián Ferrer b,2 arxiv:1503.03358v1 [nlin.si] 11 Mar 2015 Abstract a Rua Talim, 330, 12231-280 São José

More information

Simulation study on the nonlinear EMIC waves

Simulation study on the nonlinear EMIC waves SH21B-2210 Simulation study on the nonlinear EMIC waves Kicheol Rha 1*, Chang-Mo Ryu 1 and Peter H Yoon 2 * lancelot@postech.ac.kr 1 Department of Physics, Pohang University of Science and Technology,

More information

Superposition of electromagnetic waves

Superposition of electromagnetic waves Superposition of electromagnetic waves February 9, So far we have looked at properties of monochromatic plane waves. A more complete picture is found by looking at superpositions of many frequencies. Many

More information

Physics 215b: Problem Set 5

Physics 215b: Problem Set 5 Physics 25b: Problem Set 5 Prof. Matthew Fisher Solutions prepared by: James Sully April 3, 203 Please let me know if you encounter any typos in the solutions. Problem 20 Let us write the wavefunction

More information

Towards stability results for planetary problems with more than three bodies

Towards stability results for planetary problems with more than three bodies Towards stability results for planetary problems with more than three bodies Ugo Locatelli [a] and Marco Sansottera [b] [a] Math. Dep. of Università degli Studi di Roma Tor Vergata [b] Math. Dep. of Università

More information