Skip to Content.
Sympa Menu

opal - Re: [Opal] OPAL code for simulating RFQs

opal AT lists.psi.ch

Subject: The OPAL Discussion Forum

List archive

Re: [Opal] OPAL code for simulating RFQs


Chronological Thread 
  • From: Matt Easton <matt.easton AT pku.edu.cn>
  • To: <Opal AT lists.psi.ch>
  • Subject: Re: [Opal] OPAL code for simulating RFQs
  • Date: Mon, 17 Sep 2018 12:24:23 +0800

Forget that last part; it’s the same problem I had on Cygwin, with HDF5 installed without parallelisation. Fixed by `brew reinstall hdf5 --with-mpi`.

 

From: Matt Easton <matt.easton AT pku.edu.cn>
Date: Monday 17 September 2018 at 10:57
To: <Opal AT lists.psi.ch>
Subject: Re: [Opal] OPAL code for simulating RFQs

 

Hi there,

 

Is anyone from the OPAL development team attending Linac 2018 right now? I’d love to have a chance to talk through some of these issues in person.

 

I’m now trying to install OPAL on my MacBook, where I don’t have the luxury of the precompiled binaries. At the point of installing H5Hut, it doesn’t seem to be able to find my Open-MPI installation, even though everything is installed in /usr/local. I’m getting the following error:

 

                src/h5core/private/h5_types.h:38:9: error: unknown type name 'MPI_Comm'

 

Any idea how to fix this?

 

Many thanks for any help,

Matt


Dr Matt Easton 
董昊

Post-Doctoral Researcher 博士后研究员

Institute of Heavy Ion Physics 重离子研究所

Peking University 北京大学
matteaston.net/work
matt.easton AT pku.edu.cn
+86 185-1427-0582
北京市海淀区成府路201号加速器楼220

Sent from Outlook 2013

 

From: Matt Easton èŒ£æ˘… <matt.easton AT pku.edu.cn>
Organization: Peking University
北京大学
Date: Thursday 23 August 2018 at 08:33
To: <Opal AT lists.psi.ch>
Subject: RE: [Opal] OPAL code for simulating RFQs

 

Dear all,

 

I’m very new to using OPAL, so I’d appreciate a bit of guidance. In order to model an RFQ, I want to import a field map. How do I add this field map into the input file?

 

On the wiki, I can only find how to apply field maps as edge fields on a magnet. I want the field map to fill the space where I currently have a drift in the attached input file.

 

I’m sure I’ve made many other mistakes in the input file. Please do correct me—I’m eager to learn!

 

Many thanks for any help.

Matt


Dr Matt Easton
董昊

Post-Doctoral Researcher 博士后研究员

Institute of Heavy Ion Physics 重离子研究所

Peking University 北京大学
matteaston.net/work
matt.easton AT pku.edu.cn
+86 185-1427-0582
北京市海淀区成府路201号加速器楼220

Sent from Outlook 2016

 

From: srs0=ip/g=lf=lists.psi.ch=opal-owner AT srszone.org <srs0=ip/g=lf=lists.psi.ch=opal-owner AT srszone.org> On Behalf Of Adelmann Andreas (PSI)
Sent: Wednesday, August 22, 2018 11:42 PM
To: Daniel Winklehner <winklehn AT mit.edu>; Matt Easton
董昊 <matt.easton AT pku.edu.cn>
Cc: Opal AT lists.psi.ch
Subject: Re: [Opal] OPAL code for simulating RFQs

 

Dear all, let me clarify a few things:

 

 

Dear Matt,

 

To be honest, I haven’t pursued using OPAL for RFQ’s much after our last email exchange.

 

To answer some of your questions: I don’t think OPAL is related to Impact, but maybe it uses similar style input? Andreas can definitely comment on that.

 

OPAL is an electrostatic PIC and that is the only relation to Impact. 

 

 





You can definitely load a 3D electrostatic field in OPAL-T and give it a sinusoidal variation at your desired frequency. What is missing, in my opinion, to make it suitable for RFQ’s are the correct periodic boundary conditions in the longitudinal direction (they had been implemented previously, but don’t seem to be working at the moment).

 

This is correct !

 





For testing purposes, you could inject a beam that has length 2-3 times of an RF period. And select only the particles from the middle during analysis. This is a bit of a crude workaround. Another question is particle termination on the RFQ vanes. I think all the code is there to even load a meshed (vtk) version of the actual vanes and use those for termination, but I don’t know if it is implemented for use in the T tracker (I used it for the central region of our cyclotron in the OPAL-cycl flavor).

 

 

Yes the code is there however OPAL-t can not access the routines for the moment.

 





For our project, I have been working with a student to create a python module that works seamlessly with another code, called WARP, and we have seen some success there. It is close to the point where I would feel comfortable if other people wanted to use it. I can definitely let you know how that endeavor progresses.

 

That being said, I would still like to see OPAL being able to handle RFQ’s as well. Benchmarking of the two codes against each other and against known codes like TOUTATIS would be really helpful. 

 

I agree it would be very helpful, however someone needs to take up the task and 

 

 - fix the boundary conditions

 - add the geometry to OPAL-t

 

 

Cheers Andreas 

 

 





 

Hope this helps a little bit,

   Daniel

 

From: Matt Easton 董昊 <matt.easton AT pku.edu.cn> 
Sent: Tuesday, August 21, 2018 2:46 AM
To: Daniel Winklehner <
winklehn AT mit.edu>
Cc: 'Adelmann Andreas (PSI)' <
andreas.adelmann AT psi.ch>
Subject: RE: OPAL code for simulating RFQs

 

Dear Daniel,

 

What is the latest status of simulating an RFQ using OPAL? I know it was a long time ago that I first asked about this, but with one thing and another it would be useful for me now to be able to import a field map for an RFQ.

 

In my fiddling around with OPAL and trying to learn from the Wiki, I think it said that OPAL-t is related in some way to Impact-T. Is that right? Right now I can model my RFQ in Impact-T by importing a 3D electrostatic field map and applying a sinusoidal variation at the RF frequency. I’d love to be able to do the same thing in OPAL.

 

I’ve set up a crude simulation with a single drift and a square aperture for the length of the RFQ, so all that is missing is the field. Fromhttps://gitlab.psi.ch/OPAL/Manual-2.0/wikis/fieldmaps it seems that such 3D dynamic field maps are supported, but I think these are applied only as edge fields for a magnet. How can I use a field map by itself, or apply it to a drift space?

 

Thanks for your help. If it would be better to send this on the OPAL support mailing list, I’m happy to do that too.

 

Regards,

Matt


Dr Matt Easton 
董昊

Post-Doctoral Researcher 博士后研究

Institute of Heavy Ion Physics 重离子研究所

Peking University 北京大学
matteaston.net/work
matt.easton AT pku.edu.cn
+86 185-1427-0582
北京市海淀区成府路201号加速器楼220

Sent from Outlook 2016

 

From: Daniel Winklehner <winklehn AT mit.edu> 
Sent: Friday, January 19, 2018 9:52 AM
To: Matt Easton 
董昊 <matt.easton AT pku.edu.cn>; 'Adelmann Andreas (PSI)' <andreas.adelmann AT psi.ch>
Cc: 'Suzanne Sheehy' <
suzie.sheehy AT physics.ox.ac.uk>
Subject: RE: OPAL code for simulating RFQs

 

Dear Matt,

 

Apologies for taking so long to respond. Unfortunately, currently, I would not recommend using OPAL for RFQ’s.

I have done some preliminary studies with some success, but there are a few issues that need to be solved in order for it to become widely usable. As Andreas said, the idea is to import a field map from a BEM or FEM program and track the particles in OPAL. This requires correct particle termination and periodic boundary conditions, both of which are currently a work in progress. These issues are being addressed, or will be in the near future, though and I will be happy to let you know about the new results when the time comes.

 

Best,

   Daniel

 

From: Matt Easton 董昊 [mailto:matt.easton AT pku.edu.cn] 
Sent: Monday, January 15, 2018 9:14 PM
To: 'Adelmann Andreas (PSI)' <
andreas.adelmann AT psi.ch>
Cc: 'Suzanne Sheehy' <
suzie.sheehy AT physics.ox.ac.uk>; Daniel Winklehner <winklehn AT mit.edu>
Subject: RE: OPAL code for simulating RFQs

 

Dear Andreas,

 

Thanks for your response, and thanks to Suzie for putting us in touch. I’d be very interested to hear the status of the RFQ modelling effort, thanks for connecting me with Daniel.

 

Warm regards,

Matt


Dr Matt Easton 
董昊

Post-Doctoral Researcher 博士后研究

Institute of Heavy Ion Physics 重离子研究所

Peking University 北京大学
matteaston.net/work
matt.easton AT pku.edu.cn
+86 185-1427-0582
北京市海淀区成府路201号加速器楼220

Sent from Outlook 2016

 

From: Adelmann Andreas (PSI) [mailto:andreas.adelmann AT psi.ch] 
Sent: Monday, January 15, 2018 2:09 PM
To: Matt Easton 
董昊 <matt.easton AT pku.edu.cn>
Cc: Suzanne Sheehy <
suzie.sheehy AT physics.ox.ac.uk>; Daniel Winklehner <Winklehn AT mit.edu>
Subject: Re: OPAL code for simulating RFQs

 

Hello Matt yes indeed we have particle-matter integration integration in OPAL (https://gitlab.psi.ch/OPAL/src/wikis/home).

 

Daniel (in cc) can update you on our RFQ modelling effort w.r.t the IsoDAR RFQ injection project. If I remember correctly, we 

aim to import field maps from RFQ-Design programs and the track in OPAL.

 

In case you consider adding a specific RFQ Element in OPAL please do so, OPAL is open source easy to extend. 

 

Thanks for the interest 

 

Andreas 

 

 

------
Dr. sc. math. Andreas (Andy) Adelmann
Staff Scientist
Paul Scherrer Institut WBBA/219 CH-5232 Villigen PSI
Phone Office: xx41 56 310 42 33 Fax: xx41 56 310 31 91
Phone Home: xx41 62 891 91 44
-------------------------------------------------------
Friday: ETH CAB F37 or HPK G 28   +41 44 632 75 22
============================================
The more exotic, the more abstract the knowledge, 
the more profound will be its consequences.
Leon Lederman 
============================================ 

 

On 15 Jan 2018, at 06:00, Suzanne Sheehy <suzie.sheehy AT physics.ox.ac.uk> wrote:

 

HI Matt,  

Best to chat to Andreas Adelmann about those aspects, CC’d in.

I should be at IPAC this year as I’m planning to go down to Fermilab after for the IOTA meeting. Hope we can grab a drink there?

Cheers,

Suzie

On 9 Jan 2018, at 03:07, Matt Easton 董昊 <matt.easton AT pku.edu.cn> wrote:

 

Hi Suzie,

 

How are you doing? It’s been a while since we’ve been in touch. Are you going to be at IPAC this year?

 

I’ve been working on adding some particle–particle / particle–matter interactions (currently in quite a basic form) to the tracking code Impact-T, to which we’ve already added RFQ and DTL tracking modules. But then I saw a recent paper about OPAL, that it already has proper Monte Carlo simulations incorporated into it. And then I saw your name on the list of people that developed it!

 

Would you be able to tell me a bit more about OPAL, or point me in the direction of the right person to talk to? I’m interested in tracking high intensity beams through an RFQ, and processes such as stripping electrons from ions and dissociation of compound ions into smaller particles, etc.

 

Thanks a lot,

Matt


Dr Matt Easton 
董昊

Post-Doctoral Researcher 博士后研究

Institute of Heavy Ion Physics 重离子研究所

Peking University 北京大学
matteaston.net/work
matt.easton AT pku.edu.cn
+86 185-1427-0582
北京市海淀区成府路201号加速器楼220

Sent from Outlook 2016

 




Archive powered by MHonArc 2.6.19.

Top of Page