1 / 9

Positron Capture Simulations: Runge-Kutta vs Boris

Positron Capture Simulations: Runge-Kutta vs Boris. Leo Jenner, Daresbury Labs. Outline …. Positron capture optics What is astra? Spin tracking… T-BMT equation Runge-Kutta integration method Boris’ method. Capture Optics Design. V.S.Kashikhin. z. z.

saxon
Download Presentation

Positron Capture Simulations: Runge-Kutta vs Boris

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Positron Capture Simulations: Runge-Kutta vs Boris Leo Jenner, Daresbury Labs Spin Tracking Meeting

  2. Outline… • Positron capture optics • What is astra? • Spin tracking… • T-BMT equation • Runge-Kutta integration method • Boris’ method Spin Tracking Meeting

  3. Capture Optics Design V.S.Kashikhin z z Two coil magnet gives focussing solenoid field Spin Tracking Meeting

  4. ASTRA Klaus Floettmann Spin Tracking Meeting

  5. T-BMT Equation Lorentz Factor Differential Equation… Gyromagnetic anomaly Spin Tracking Meeting

  6. Let an initial value problem be specified as follows. • Then, the RK4 method for this problem is given by the following equation: • where • Thus, the next value (yn+1) is determined by the present value (yn) plus the • product of the size of the interval (h) and an estimated slope. The slope is a weighted average of slopes: • k1 is the slope at the beginning of the interval; • k2 is the slope at the midpoint of the interval, using slope k1 to determine the value of y at the point tn + h/2 using Euler's method; • k3 is again the slope at the midpoint, but now using the slope k2 to determine the y-value; • k4 is the slope at the end of the interval, with its y-value determined using k3. Runge-Kutta Integrator Scheme Andriy Ushakov Spin Tracking Meeting

  7. Boris’ Numerical Integration Scheme Jeff Smith • Advance w by vector term b a half-step • Advance w by M a full step • Advance w by b a half-step Spin Tracking Meeting

  8. R-K vs Boris Boris: 2nd order, good at preserving conserving quantities, 1 calculation per step R-K: 4th order, 4 calculations per step (speed!) Will try both! Spin Tracking Meeting

  9. Conclusions… • ASTRA can be used to simulate the motion of particles/bunches in the capture optics section • …but doesn’t have spin tracking included… • Spin tracking subroutine using Boris integrator underway • Spin tracking subroutine using Runge-Kutta integrator underway – cross-check • Merlin? Spin Tracking Meeting

More Related