Control of chaos in Hamiltonian systems
|
|
- Allen Parks
- 6 years ago
- Views:
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 Guido Ciraolo Facoltà di Ingegneria, Università di Firenze via S. Marta, I-50129 Firenze, Italy Cristel Chandre, Ricardo Lima, Michel Vittot CPT-CNRS,
More informationChaotic 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 informationarxiv: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 informationGyrokinetic simulations of magnetic fusion plasmas
Gyrokinetic simulations of magnetic fusion plasmas Tutorial 2 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr
More informationEUROFUSION 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 informationInternal 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 informationChanneling 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 informationA 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 informationHamiltonian 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 informationEnhancement 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
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 informationChaotic 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 informationOscillatory 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 informationOscillatory 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 informationHamiltonian 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 informationPhysics 106b: Lecture 7 25 January, 2018
Physics 106b: Lecture 7 25 January, 2018 Hamiltonian Chaos: Introduction Integrable Systems We start with systems that do not exhibit chaos, but instead have simple periodic motion (like the SHO) with
More informationAn 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 informationCPT 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 informationTheory 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 informationHamiltonian 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 informationChaotic 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 informationOverview 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 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 informationBeams 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 informationAssignment 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 informationarxiv: 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 informationChapter 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 informationMechanics 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 informationOh, 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 informationLectures on basic plasma physics: Hamiltonian mechanics of charged particle motion
Lectures on basic plasma physics: Hamiltonian mechanics of charged particle motion Department of applied physics, Aalto University March 8, 2016 Hamiltonian versus Newtonian mechanics Newtonian mechanics:
More informationWhat 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 informationHamiltonian 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 informationIntermittent 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 informationIon 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 informationBursty 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 informationM2A2 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 informationOn 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 informationGeometric 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 informationPhysics 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 informationA 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 informationPlasma 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 informationChapter 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 informationElectron 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 informationScattering 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 informationNon-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 informationGEOMETRIC 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 informationNotes 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 informationA 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 informationLectures 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 informationSome 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 informationDOPPLER RESONANCE EFFECT ON ROTATIONAL DRIVE BY ION CYCLOTRON MINORITY HEATING
DOPPLER RESONANCE EFFECT ON ROTATIONAL DRIVE BY ION CYCLOTRON MINORITY HEATING V.S. Chan, S.C. Chiu, Y.A. Omelchenko General Atomics, San Diego, CA, U.S.A. 43rd Annual APS Division of Plasma Physics Meeting
More informationStability 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 informationThe Linear Theory of Tearing Modes in periodic, cyindrical plasmas. Cary Forest University of Wisconsin
The Linear Theory of Tearing Modes in periodic, cyindrical plasmas Cary Forest University of Wisconsin 1 Resistive MHD E + v B = ηj (no energy principle) Role of resistivity No frozen flux, B can tear
More informationNon-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 informationTURBULENT 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 informationActa 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 informationAction-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 informationFlow, 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 informationarxiv: 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 informationTHE 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 informationLecture 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 informationFormation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )
Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University
More informationBACKGROUND 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 informationGyrokinetic 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 informationNonlinear 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 information2 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 informationRecursion 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 informationHeating 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 informationKAM 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 informationPage 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 informationImpact 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 informationAnalysis 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 informationModeling 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 informationLecture 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 informationDynamics of charged particles in spatially chaotic magnetic fields
PSFC/JA-1-38 Dynamics of charged particles in spatially chaotic magnetic fields Abhay K. Ram and Brahmananda Dasgupta a October 21 Plasma Science and Fusion Center, Massachusetts Institute of Technology
More informationWaves 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 informationNotes 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 informationSMR/ 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 information1 (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 informationChaos 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 informationSample 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 informationCalculation 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 informationOrbits, 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 informationEffect 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 informationCHAPTER 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 informationFAST 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 informationTheory 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 informationNumerical 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 informationINTERACTION 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 informationEigenmodes 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 informationLinear 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 information9.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 informationTwo 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 informationarxiv: 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 informationSimulation 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 informationSuperposition 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 informationPhysics 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 informationTowards 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