The OPAL Discussion Forum

["Calvo Portela, Pedro" <Pedro.Calvo AT ciemat.es>]

Dear Ezgi

Could you please provide me the input file, as well as the fieldmap and the distribution you are using? So I can try to run your simulation and find the specific error.

Regarding the gnu.out and eb.out files you need "OPTION INFO=TRUE to be activated.

Best regards

Pedro

-----------------------------------------------

Dr. Pedro Calvo Portela
Particle Accelerator Unit

Technology Department

CIEMAT - Av. Complutense 40, 28040 Madrid (Spain)
Tel: (+34) 914962542 / ext. 362542

---

De: opal-request AT lists.psi.ch <opal-request AT lists.psi.ch> en nombre de Ezgi Sunar <ezgi.sunar AT mnt.com.tr>
Enviado: martes, 18 de junio de 2024 13:41
Para: Adelmann Andreas; opal AT lists.psi.ch
Cc: Milo Conard; Uğur Bozluolcay; Ezgi Sunar
Asunto: RE: [Opal] Problem with Cyclotron Tune and Phase Calculation

 

Dear Andy,

 

Thank you so much for your prompt reply. The insightful explanation regarding the repository page helped me a lot. However, I need your support one more time. 

I've thoroughly explored the repository, but I'm stuck finding information on how to calculate orbits in this specific mode. Here's what I've encountered during my trials:

1. Single Particle Distribution: When running in single-particle mode, I utilize a "one particle distribution" with the "fromfile" type.
   This approach, according to the documentation, should generate output files named "gnu.out" and "eb.out." However, these files remain empty after the simulation.
   I can successfully generate a "track-Orbit.dat" file. Unfortunately, this file only shows the positions of the reference particle at only one orbit.
2. Gaussmatched Distribution: Switching to the "gaussmatched" distribution offers valuable insights into energy increase per orbit using Gordon Units. It generates a "tunes.dat" file containing data like energy increase, average radius, momentum, nu_R, and nu_Z, useful for the tune diagram, but lacking the detailed information needed for a full understanding of the orbits themselves.
3. Two-Particle Test: I also attempted to use two particles for Static Equilibrium Orbit (SEO) information. Unfortunately, this approach didn't yield significant insights.
4. Tuning Script Error: Additionally, I attempted running the "tuning.sh" script for Static Equilibrium Orbit (SEO) calculations using the cyclotron example from the OPAL wiki, but the script returned an "invalid file descriptor" error.

My understanding (please correct me if I'm wrong) is that the workflow involves running a single-particle simulation first. This analyzes the orbits at increasing energy levels. Subsequently, a two-particle simulation would be used with the aim of obtaining accurate tune calculations. However, even if this is the correct approach, I'm unsure how to use the output from the single-particle simulation as input for the two-particle one.

I've attached all the relevant output files for your reference.

Thank you once again for your continued support. I greatly appreciate your assistance in navigating this challenge.

Sincerely,
Ezgi

 

From: Adelmann Andreas <andreas.adelmann AT psi.ch>
Sent: Wednesday, May 29, 2024 9:54 AM
To: Ezgi Sunar <ezgi.sunar AT mnt.com.tr>
Cc: opal AT lists.psi.ch; Milo Conard <milo AT pac.be>; Uğur Bozluolcay <ugur.bozluolcay AT bozlu.com.tr>; Ezgi Sunar <sunarezgi AT gmail.com>
Subject: Re: [Opal] Problem with Cyclotron Tune and Phase Calculation

 

Hi Sunar

 

1. How and where can I define the initial and final energies in the input file?

 

You specify the injection energy in the input file, check line 44 in 

https://gitlab.psi.ch/OPAL/documentation/manual/blob/master/examples/Injector2.in for the PSI Injector example.

P0 will then be used on line 282 on the BEAM command

2. Which type of cyclotron should I run for this application? I suppose it would be CARBONCYCL. If so, I have some difficulties with the necessary structures of the data files; like BFIELD.

 

YES, CARBONCYCL and the detailed description can be found in the manual section 7.4.1.

3. How and where can I define and calculate the orbit coordinates according to the given energies in the input file? I found that OPAL uses the external program FIXPO for calculating the initial orbit coordinates based on our desired energy range. If this is correct, I cannot find how and where to set these parameters. 

 

FIXPO is used to calculate the reference solution only. This is only for making sure the OPAL tune calculation is correct. 

 

To me your particle (initial condition) starts out of the field map, please check.

 

Cheers Andy 
------
Dr. sc. math. Andreas (Andy) Adelmann
Paul Scherrer Institut OHSA/D17 CH-5232 Villigen PSI
Phone Office: xx41 56 310 42 33 Fax: xx41 56 310 31 91
Zoom ID: 470-582-4086 Password: AdA
Zoom Link: https://ethz.zoom.us/j/4705824086?pwd=dFcvT1pMMGY0bHg0dTNncUNZZTJkZz09

-------------------------------------------------------
Friday: ETH HPK G 28   +41 44 633 3076
============================================
The more exotic, the more abstract the knowledge, 
the more profound will be its consequences.
Leon Lederman 
============================================

On 28 May 2024, at 16:16, Ezgi Sunar <ezgi.sunar AT mnt.com.tr> wrote:

 

Hello again,

 

I apologize for the oversight in my previous email. The error message encountered during the run is provided in the attachment.

Thanks for your time.

 

Best Regards,

Ezgi

 

From: Ezgi Sunar 
Sent: Tuesday, May 28, 2024 5:02 PM
To: opal AT lists.psi.ch
Cc: Milo Conard <milo AT pac.be>; Uğur Bozluolcay <ugur.bozluolcay AT bozlu.com.tr>; Ezgi Sunar <sunarezgi AT gmail.com>
Subject: Problem with Cyclotron Tune and Phase Calculation

 

Dear Sir/Madam,

 

I am working on the design of a low energy cyclotron (for industrial applications).

The first result I’d like to get is a simple tune calculation in order to optimise the magnetic field

To achieve this, I have modified the cyclotron example file from the OPAL wiki page to suit my design parameters. Here is the relevant section of the input file:

 

Name_Cycl: CYCLOTRON, TYPE=CARBONCYCL, CYHARMON=4, FMLOWE=0.5, FMHIGHE=3.0, PHIINIT=0.0, PRINIT=pr0, RINIT=r0, SYMMETRY=4.0, RFFREQ=f1, FMAPFN="bfield.dat", BSCALE=1E-3;

L1: LINE = (Name_Cycl);

Dist1: DISTRIBUTION, TYPE=FILENAME, NSTEPS=100, SECTOR=FALSE, FNAME="dist1.dat";

Fs1: FIELDSOLVER, FSTYPE=NONE, MX=64, MY=64, MT=64, PARFFTX=true, PARFFTY=true, PARFFTT=false, BCFFTX=open, BCFFTY=open, BCFFTT=open;

Beam1: BEAM, PARTICLE=PROTON, PC=P0, NPART=2, BCURRENT=1.0E-6, BFREQ=f1;

SELECT, LINE=L1;

TRACK, LINE=L1, BEAM=Beam1, MAXSTEPS=nstep*turns, STEPSPERTURN=nstep;

RUN, METHOD="CYCLOTRON-T", BEAM=Beam1, FIELDSOLVER=Fs1, DISTRIBUTION=Dist1;

 

As part of the tune calculation mode, I am using data distribution for two particles: one as the reference particle with all data set to zero and another one off-center. The `dist1.dat` file includes the following information:

 

2

0.0         0.0           0.0           0.0          0.0           0.00

0.005       0.001         0.005         0.001        0.005         0.001

 

Additionally, the B-field data is created with:

- Grid points along azimuth (ntet): 239

- Grid points along radius (nrad): 4

 

4

4

0.0

1.5

239

48

11296.6  11317.7  11327.6  11337.9  11354.2

11367.0  11378.0  11393.0  11412.4  11435.4

11453.5  11472.2  11495.2  11523.0  11551.7 ......

 

Upon running this configuration, I encountered a few issues and a failure during execution. I would appreciate your guidance on the following points:

 

1. How and where can I define the initial and final energies in the input file?

2. Which type of cyclotron should I run for this application? I suppose it would be CARBONCYCL. If so, I have some difficulties with the necessary structures of the data files; like BFIELD.

3. How and where can I define and calculate the orbit coordinates according to the given energies in the input file? I found that OPAL uses the external program FIXPO for calculating the initial orbit coordinates based on our desired energy range. If this is correct, I cannot find how and where to set these parameters.

 

I have attached the error message encountered during the run for your reference.

 

Thank you for your assistance.

 

Best regards,  

Ezgi Sunar