Transcript Huang_lattice_modelling

Lattice modeling for a storage ring with magnetic
field data
X. Huang, J. Safranek (SLAC)
Y. Li (BNL)
3/5/2012
3/5/2012
FLS2012, Ring WG, X. Huang
1
Motivation
• Discrepancy between original model and measurements.
• Understanding dynamic effects of rectangular gradient
dipoles.
x
X
s
Z
1. The ideal trajectory in RGD is not a circular arc.
2. Gradient varies with s-variable
3. Off-plane longitudinal field
• Understanding the sources of discrepancies in linear and
nonlinear characteristics between models and
measurements.
– Fringe field of dipoles
– Fringe field of quadrupoles
– Cross-talk of fields between adjacent magnets?
3/5/2012
FLS2012, Ring WG, X. Huang
2
The field-integration approach
An AT pass-method that transfer phase space coordinates from one end to the other
of a magnetic field region with Bx, By, Bz defined as function of (x, y, z).
Equation of motion when using z as free variable.
Coordinate transformation at the edges.
For dipoles, additional transformation is needed.
3/5/2012
FLS2012, Ring WG, X. Huang
3
Magnetic field in a standard SPEAR3 dipole
We started examining our lattice model from magnetic fields in magnets.
By at z=0
1.5
1.4
Tanabe
RADIA new
Measured 2007
By (normalized)
1
RADIA new
Measured 2001
fit meas 2001
measured 2007
0.5
0
-1
-0.5
0
Z (m)
1
0.8
0.6
0.4
0.2
0.5
1
We have coil, wire measurements.
Hall probe scans along Z in 2001, 2007.
Hall probe scans on X-Z plane in 2007.
0
-100
-50
0
50
100
x (mm)
150
200
250
300
0
-0.5
By (T)
By/Bymax
1.2
-1
-1.5
1
0.8
0.6
Z (m) 145D
3/5/2012
0.4
-0.05
0
0.05
0.1
0.15
0.25
0.2
X (m)
FLS2012, Ring WG, X. Huang
Hall probe x-z scan (2007)
4
0.3
An analytic dipole field model
An analytical field model can be built according to general field expansion to obtain
the full magnetic field distribution in the dipole. This also removes noise from field
measurements.
Take the dipole component as an example.
Note that the B0/B1 ratio is not constant in the fringe region.
3/5/2012
FLS2012, Ring WG, X. Huang
5
Energy calibration
How to calibrate the require bending field for a 3GeV beam?
+
2001 Z-scan
(Corbett & Tanabe, 2002)
(Yoon, et al, NIMA 2004). The virtual center
was held constant (392.35 mm).
This is how dipole magnets were positioned: adjust the dipole current
(converted to K-value) until the alignment requirement is met.
Following this procedure, the required field integral is calculated to be
(1) 1.86420 T-m, with a fixed virtual center, while the measured field integral is
1.86413 T-m for 587.6909 A (operating current since day 1 of SPEAR3).
(2) 1.86615 T-m, with the fitted field profile.
So the SPEAR3 beam energy may be lower than the nominal value by 0.1%.
Energy measurement at SPEAR3 confirmed the prediction with high precision.
3/5/2012
FLS2012, Ring WG, X. Huang
6
Effects of quadrupole fringe field
A general Hamiltonian (including longitudinal field variation) can be derived using a
proper magnetic field expansion(1).
H ( s) 
1 2
1
1
1
( Px  Py2 )  k ( s)( x 2  y 2 )  k ' ( s)( x 2  y 2 )( xPx  yPy )  k ' ' ( s)( x 4  y 4 )  O( X 6 )
2
2
4
12
J. Irwin, C.X. Wang
The leading correction for a hard-edge model is from the last two terms, which are
nonlinear(2).
The leading correction term from a soft fringe model is linear(3).
A perturbation approach
* El-Kareh; Forest; Bassetti & Biscari
** Lee-Whiting, Forest & Milutinovic, Irwin & Wang, Zimmermann
***Irwin & Wang (PAC’95), D. Zhou (IPAC10).
3/5/2012
FLS2012, Ring WG, X. Huang
7
The linear correction to quadrupole map
The correction map
J. Irwin, C.X. Wang, PAC95
The generating function for the correction map
f2 
I1
( xPx  yPy )
2
matrix
leading contribution
diag(e I1 , e I1 , e I1 , e I1 )
For a symmetric quadrupole, the entrance edge
has a reversed sign for I1
The tune changes are (always negative)
k02 L x I1
 x  
| |,
2 k0
k02 L y I1
 y  
| |
2 k0
For SPEAR3, quadrupole fringe fields cause tune changes of [-0.065, -0.059], in
agreement with the predictions by the above equation.
3/5/2012
FLS2012, Ring WG, X. Huang
8
The nonlinear correction
The generating function for the correction map (exit edge)
f4 
1
1
k0 ( x 3 Px  3xy 2 Px  y 3 Py  3x 2 yPy ) 
k skew ( x 3 Py  y 3 Px )
12(1   )
6(1   )
Forest & Milutinovic
The function for the entrance edge has an opposite sign.
F. Zimmermann derived the average Hamiltonian that include both edges.
Hard edge

Additional soft edge contribution.
2 is fringe length.
6 I1
k0
and tune dependence on amplitude (only showing hard edge contribution below)
This agree with tracking quite well.
3/5/2012
FLS2012, Ring WG, X. Huang
9
An AT quadrupole passmethod with fringe field
Both linear and nonlinear effects are considered in the new quadrupole passmethod.
Forest & Milutinovic pointed out the skew quadrupole part corresponds to a ‘kick map’!
A normal quadrupole can thus be modeled by a pair of pi/4 rotation and a kick map.
This is the basis for the nonlinear part of the new AT quadrupole pass method.
-0.0122
-0.0122
The new quad passmethod agree
very well with the field-integration
method.
xpf (rad)
-0.0122
-0.0123
-0.0123
-0.0123
-0.0123
quadpass
quadpass+matrix
new quad pass
field pass
-0.0123
0.02
3/5/2012
0.02
FLS2012, Ring WG, X. Huang
0.02
0.02
0.02
xi (m)
0.02
0.02
10
The SPEAR3 quadrupole field profile
The analytical quadrupole field map for SPEAR3 magnet was based on magnet
modeling. Simulated field is converted to an analytic form.
1
' ' ( z )(3x 2 y  y 3 )]
12
1
B y  B1[ x( z )  ' ' ( z )(x 3  3xy 2 )]
12
Bz  sgn( z ) B1' ( z ) xy
1.5
Magnetic field
B1f (normalized)
Bx  B1[ y( z ) 
1
0.5
0
0
0.1
all quads at 72 A in modeling
I1a 
18
16
60Q
50Q
34Q
15Q
B1 (T/m)
14
12
0.2
0.3
0.4
Z (m) for 60Q
0.5
0.6
0.7
I1
6 I1
 0.61103 m 2 ,  
 0.060 m
k0
k0
10
8
6
4
2
0
-0.4
-0.3
-0.2
-0.1
0
0.1
0.2
Z (m)
0.3
0.4
0.5
0.6
All SPEAR3 quads have identical fringe profile.
3/5/2012
FLS2012, Ring WG, X. Huang
11
Comparison of models to measurement
The model is based on a calibrated experimental lattice with all IDs open (4/6/2009).
Parameter
Measure
d
All field model
Bend field,
quad i2k
with fringe
Bend field,
quad i2K
i2k old AT
model
Tune x
14.106
14.146
14.150
14.215
14.190
Tune y
6.177
6.119
6.121
6.180
6.431
Chrom x
1.7
-0.54
-0.53
-0.44
-0.60
Chrom y
2
0.89
0.90
0.73
1.90
Effect of the predicted -0.1% beam energy shift is not included, which change the tunes by
[0.023, -0.004] for [nux, nuy].
Dipole field is given by the field profile and alignment requirement.
Drift lengths neighboring to dipoles adjusted according to measured rf frequency.
Strengths of quadrupoles and sextupoles are derived from operating currents and
measured excitation curves.
No adjustment of any magnet strength!
The tune differences are [0.067 -0.060] between the best model and the measurement, a
big improvement from the original model.
3/5/2012
FLS2012, Ring WG, X. Huang
12
Beta beat and correction
0.1
0.2
0.15
0.1
 y/y
 x/x
0.05
0
0
-0.05
-0.05
-0.1
0
0.05
50
100
150
s (m)
no correction
corrected
-0.1
200
-0.15
0
250
50
100
150
s (m)
200
250
Beta beat is relative to the ideal lattice.
“No correction” is for “bend field + quad fringe”
“corrected” is after the quadrupole strength is adjusted to reduce beta beat (LOCO).
Possible causes of optics difference between measurement and un-adjusted
model:
(1) Interference of magnetic fields between neighboring magnets.
(2) Magnet calibration errors.
3/5/2012
FLS2012, Ring WG, X. Huang
13
The tune map
Chromaticities are corrected with SF/SD to obtain [1.65, 2.18].
Tunes are obtained by tracking 256 turns.
old model vs. new model
0.25
0.25
0.24
0.24
0.23
0.23
0.22
0.22
y
y
measured vs. new model
0.21
0.21
0.2
0.2
0.19
0.19
0.18
0.18
0.1
0.15
x
0.2
0.1
0.15
x
0.2
“new model” = field model for bend + quad fringe.
This model agrees with measurement better.
3/5/2012
FLS2012, Ring WG, X. Huang
14
High-order chromaticities
low tune [0.106, 0.177] LE
0.17
0.16
x-Cx  p/p
0.15
measured
old model
new model
0.14
0.13
0.12
0.11
0.1
-0.04
-0.02
0
 p/p
0.02
measured
New
model
chrx0
1.725*
1.725
1.647**
chrx1
28.8
28.4
30.6
chrx2
-569
-545
-557
chry0
2.081*
2.081
2.181**
chry1
16.7
30.3
28.9***
chry2
-207
16
-213
* Model chromaticities adjusted to match
measured values.
** model chromaticities adjusted, but not yet
completely on target.
*** improvement from old model.
0.23
y-Cy  p/p
Old
model
0.04
low tune [0.106, 0.177] LE
measured
old model
new model
0.22
Low
tune
0.21
0.2
0.19
0.18
0.17
-0.04
3/5/2012
-0.02
0
 p/p
0.02
0.04
FLS2012, Ring WG, X. Huang
15
Progress toward fast tracking for dipoles
(1) Extract a Lie map from the Taylor map obtained from the fieldintegration method (Yongjun Li). A map may also be obtained with
COSYInfinity.
(2) Split the f3 and f4 polynomials into individual terms for tracking (f4
terms are altered by splitting f3), ignore higher order polynomials.
Monomial maps have exact solutions (A. Chao, Lie Algrebra Notes).
An AT passmethod is written to track f3 (35 terms) and f4 (70 terms) maps (f2 is
supplied by a matrix).
3/5/2012
FLS2012, Ring WG, X. Huang
16
Comparison of Map-pass to field-pass
For comparison, the second order transport map is extracted with AT for the
SPEAR3 dipole, using the field-pass or the map-pass.
T1ij from the field pass
T1ij from the map pass
All transverse-only elements agree well (for T2ij, T3ij T4ij, too).
The discrepancy for the momentum-related elements may be caused by an problem
in the field-pass used for map extraction (different from the one compared to here).
The map-pass provides a symplectic tracking solution to the dipole model.
3/5/2012
FLS2012, Ring WG, X. Huang
17
Summary and Discussion
• We built a lattice model from magnetic field measurements and
alignment requirements and compared the linear optics to beam based
measurements.
– Improved: tunes, betatron functions.
– But: still up to 15% maximum beta beat (vertical)
• After optics and chromaticity corrections, nonlinear parameters from
the model are compared to beam based measurements.
– Improved: 2nd order vertical chromaticity, tune map.
– But: the tune map is still slightly different from measurement.
•
We have developed fast symplectic method to represent high order effects.
• An accurate model may be crucial for a smooth commissioning of a
new machine and for dynamic aperture optimization of existing
machines.
– Quadrupole fringe field effect (tune shifts and beta beat) would be larger for a large
ring (with more quads).
– Magnetic field based lattice can be used as a “reference” model.
More efforts are need to understand the discrepancies between model and measurements.
3/5/2012
FLS2012, Ring WG, X. Huang
18
This slide is left blank
3/5/2012
FLS2012, Ring WG, X. Huang
19
The dipole field map
B0, B1, B2
from X-Z scan
5
2007 xz scan
-5
B0
B1
B2
-10
B
0,1,2
, (T, T/m, T/m2)
0
-15
By(at y=0,z=0) = -1.233257 + 3.143436*x -0.324508*x^2
ByL = -1.857103 + 4.662405*x -0.931245*x^2
-20
0.4
0.5
0.6
0.7
Z (m)
0.8
0.9
1
Coil measurement gives
ByL (T m)= -1.8506 +4.6081 x - 1.2632 x^2
Note that the dipole/quadrupole ratio is constant (392.35 mm) in the magnet body, but
varying in the fringe.
The integrated quadrupole component is actually 2% weaker than the present model.
(The coil measurement gives an average ratio of 399.8+-1.7 mm)
3/5/2012
FLS2012, Ring WG, X. Huang
20
The linear correction to quadrupole map
A perturbation approach
Hard-edge model, for exit edge
Perturbation term
The map
J. Irwin, C.X. Wang, PAC95
The generating function for the correction map (only leading contribution is shown)
f2 
I1
( xPx  yPy )
2
matrix
diag(e I1 , e I1 , e I1 , e I1 )
For a symmetric quadrupole, the entrance edge
has a reversed sign for I1
The tune change would be (always negative)
3/5/2012
k02 L y I1
k02 L x I1
 x  
| |,  y  
| |
2 k0
2 k0
FLS2012, Ring WG, X. Huang
21
Verification of the quad fringe pass method
-0.0122
yi=0.005 m
-0.0122
xpf (rad)
-0.0122
-0.0123
-0.0123
-0.0123
-0.0123
quadpass
quadpass+matrix
new quad pass
field pass
-0.0123
0.02
0.02
0.02
0.02
0.02
xi (m)
0.02
0.02
Quadpass: quad transfer matrix
Quadpass+matrix: quad transfer matrix + linear edge
transfer matrix.
New quad pass: with linear and nonlinear corrrection.
Field pass: integration through magnetic field.
With the (quad+matrix) part subtracted.
Zimmerman result is from the average
Hamiltonian H1+2
3/5/2012
FLS2012, Ring WG, X. Huang
22
A pass method for magnetic field in AT
(1) Coordinate transformation at the entrance and exit of the magnets
(2) Integration of the Lorentz equation in the body of magnets.
Can we study beam dynamics with such a pass method?
With an accurate magnetic field model, we can reproduce reality in simulation.
Integration is slow and non-symplectic, not good for dynamic aperture tracking.
But it should be good for linear and nonlinear parameter evaluation.
7/20/2011
23