# Presentation: 061896

## BENCHMARKING OF SIMULATION CODES

**TRACK (ANL) AND ASTRA (DESY)**
**FOR THE FNAL HIGH-INTENSITY PROTON SOURCE**
**(Paper TUP088, LINAC06, Knoxville, TN)**
- J.-P. Carneiro
- Accelerator Physics Department, FNAL
- Argonne National Laboratory
- 08-SEPT-06

## 	• From 1997 2004 I worked on beam dynamics

***	***• From 1997 2004 I worked on beam dynamics
- for electron photo-injectors (FNAL, DESY) and XFELs.
- Beam dynamics codes : ASTRA, PARMELA, HOMDYN (Italy)
- 	 • In Jan. 2005 : Moved to TEVATRON Dept.
-  	 • In Dec. 2005 : Moved to HINS Group.
- B. Foster requested benchmark TRACK  with another code.
- I selected ASTRA (can do e-, e+, proton, H-)

***A**** ****S****pace charge ****TR****acking ****A****lgorithm *
*(free executable, PC, UNIX, LINUX)*
*(from K. Floettmann, DESY Hamburg)*

*And .....*

**SOME HISTORY**

## ASTRA 

**ASTRA**

**TRACK**

- • Integration of Equation of motion by Runge-Kutta Method of Forth Order

- • id but BINARY* *

- • Can read any arbitrary ASCII input distribution

-  • id

- • Can handle 3D E fields

- • id

- • _HANDLES ONLY 2D B fields_

* *• can read only 3D B fields

- • Can handle 3D space charge
- (Poisson Equation in Free Space)

-  • 3D space charge with Image    Charge

**BRIEF CODE COMPARISON**

## Benchmarking Philosophy

- • Take design of High Intensity Proton Source from P. Ostroumov
-    done with TRACK at 45 mA
-   674 m. exit RFQ to end last cryomodule, 454 Cav., 55 Sol., 117 Quads
-  18 RT-TSR + 18 SSR1 + 33 SSR2 + 42 TSR + 56 SBETA + 287 BETA1
-   Sol MEBT to SSR, Quads from SBETA to BETA1

- • Tune Lattice for Zero Current

- • Benchmark TRACK/ASTRA for Zero Current (10k p.)

- • Get back to 45 mA lattice and benchmark TRACK/ASTRA
- with 3D space charge routine (200k p.)

## Lattice 14-APR-06 tuned for zero current H-minus beam - HINS

- • Lattice tuned for 45 mA beam with zero current beam

- • Lattice tuned for zero current beam with zero current beam

- Bunch Length [Deg]

- J.-P. Carneiro						08-SEPTEMBER-06

## • Lattice tuned for 45 mA beam with zero current beam

- • Lattice tuned for 45 mA beam with zero current beam

- • Lattice tuned for zero current beam with zero current beam

- X/Y max [cm]
- X/Y rms  [cm]
- X=blue Y=red

- Lattice 14-APR-06 tuned for zero current H-minus beam - HINS

- J.-P. Carneiro						08-SEPTEMBER-06

## ENERGY

- TRACK / ASTRA Zero Current Lattice 10k part.

**ENERGY**

- J.-P. Carneiro						08-SEPTEMBER-06

## ENERGY SPREAD

- TRACK / ASTRA Zero Current Lattice 10k part.

**ENERGY SPREAD**

- J.-P. Carneiro						08-SEPTEMBER-06

## BUNCH LENGTH

- TRACK / ASTRA Zero Current Lattice 10k part.

**BUNCH LENGTH**

- J.-P. Carneiro						08-SEPTEMBER-06

## LONGITUDINAL EMITTANCE

- J.-P. Carneiro						08-SEPTEMBER-06

- TRACK / ASTRA Zero Current Lattice 10k part.

**LONGITUDINAL EMITTANCE**

## ENVELOPE

- TRACK / ASTRA Zero Current Lattice 10k part.

- J.-P. Carneiro						08-SEPTEMBER-06

**ENVELOPE**

## End SSR section

- Mismatched TRACK / ASTRA

- J.-P. Carneiro						08-SEPTEMBER-06

- TRACK / ASTRA Zero Current Lattice 10k part.

**2D vs 3D B FIELDS ?**

## TRANSVERSE EMITTANCE

- TRACK / ASTRA Zero Current Lattice 10k part.

- J.-P. Carneiro						08-SEPTEMBER-06

**TRANSVERSE EMITTANCE**

## TRACK / ASTRA Zero Current Lattice 10k part.

- J.-P. Carneiro						08-SEPTEMBER-06

**~ 4% difference**

## TRACK / ASTRA 3D SPACE CHARGE ROUTINE UNIF. 3D ELLIPSOID 500k

-  FNAL beams_doc_2395

- J.-P. Carneiro						08-SEPTEMBER-06

## EXCELLENT AGREEMENT TRACK/ASTRA/ANALYTICAL

- J.-P. Carneiro						08-SEPTEMBER-06

**EXCELLENT AGREEMENT TRACK/ASTRA/ANALYTICAL**

- TRACK / ASTRA 3D SPACE CHARGE ROUTINE UNIF. 3D ELLIPSOID

## ENERGY

- TRACK / ASTRA 45 mA 200k part.

- J.-P. Carneiro						08-SEPTEMBER-06

**ENERGY**

## ENERGY SPREAD

- J.-P. Carneiro						08-SEPTEMBER-06

- TRACK / ASTRA 45 mA 200k part.

**ENERGY SPREAD**

## BUNCH LENGTH

- J.-P. Carneiro						08-SEPTEMBER-06

- TRACK / ASTRA 45 mA 200k part.

**BUNCH LENGTH**

## LONGITUDINAL EMITTANCE

- J.-P. Carneiro						08-SEPTEMBER-06

- TRACK / ASTRA 45 mA 200k part.

**LONGITUDINAL EMITTANCE**

## ENVELOPE

- J.-P. Carneiro						08-SEPTEMBER-06

- TRACK / ASTRA 45 mA 200k part.

**ENVELOPE**

## TRANSVERSE EMITTANCE

- J.-P. Carneiro						08-SEPTEMBER-06

- TRACK / ASTRA 45 mA 200k part.

**TRANSVERSE EMITTANCE**

## J.-P. Carneiro						08-SEPTEMBER-06

- TRACK / ASTRA 45 mA 200k part.

**~ 5% difference**

**~ 5% difference**

**~ 7% difference**

**~ 5% difference**

## J.-P. Carneiro						08-SEPTEMBER-06

- CONCLUSION BENCHMARKING TRACK / ASTRA

- Overall Good Agreement Between TRACK/ASTRA for H-minus beam
- Zero Current and 45 mA beam with 3D space charge routine
- 3D B Fields in ASTRA or 2D in TRACK should better the agreement
- First H-minus simulations for ASTRA
_**New ASTRA**_ : includes CSR + Dipoles (available soon)