Explicit integrators for the magnetized equations of motion in Particle in Cell codes
A new explicit time-reversible orbit integrator for the equations of motion in a static homogeneous magnetic field – called Cyclotronic integrator – is presented. Like Spreiter and Walter’s Taylor expansion algorithm, for sufficiently weak electric field gradients this second order method does not r...
| Main Authors: | , |
|---|---|
| Other Authors: | |
| Format: | Article |
| Language: | en_US |
| Published: |
Academic Press
2010
|
| Online Access: | http://hdl.handle.net/1721.1/58724 |
| _version_ | 1811082639580332032 |
|---|---|
| author | Patacchini, Leonardo Hutchinson, Ian Horner |
| author2 | Massachusetts Institute of Technology. Department of Nuclear Science and Engineering |
| author_facet | Massachusetts Institute of Technology. Department of Nuclear Science and Engineering Patacchini, Leonardo Hutchinson, Ian Horner |
| author_sort | Patacchini, Leonardo |
| collection | MIT |
| description | A new explicit time-reversible orbit integrator for the equations of motion in a static homogeneous magnetic field – called Cyclotronic integrator – is presented. Like Spreiter and Walter’s Taylor expansion algorithm, for sufficiently weak electric field gradients this second order method does not require a fine resolution of the Larmor motion; it has however the essential advantage of being symplectic, hence time-reversible. The Cyclotronic integrator is only subject to a linear stability constraint ([OmegaDelta t] < pi, [Omega] being the Larmor angular frequency), and is therefore particularly suitable to electrostatic Particle In Cell codes with uniform magnetic field where [Omega]is larger than any other characteristic frequency, yet a resolution of the particles’ gyromotion is required. Application examples and a detailed comparison with the well-known (time-reversible) Boris algorithm are presented; it is in particular shown that implementation of the Cyclotronic integrator in the kinetic codes SCEPTIC and Democritus can reduce the cost of orbit integration by up to a factor of ten. |
| first_indexed | 2024-09-23T12:06:36Z |
| format | Article |
| id | mit-1721.1/58724 |
| institution | Massachusetts Institute of Technology |
| language | en_US |
| last_indexed | 2024-09-23T12:06:36Z |
| publishDate | 2010 |
| publisher | Academic Press |
| record_format | dspace |
| spelling | mit-1721.1/587242023-02-26T02:06:24Z Explicit integrators for the magnetized equations of motion in Particle in Cell codes Explicit time-reversible orbit integration in Particle In Cell codes with static homogeneous magnetic field Patacchini, Leonardo Hutchinson, Ian Horner Massachusetts Institute of Technology. Department of Nuclear Science and Engineering Massachusetts Institute of Technology. Plasma Science and Fusion Center Hutchinson, Ian H. Hutchinson, Ian H. Patacchini, Leonardo A new explicit time-reversible orbit integrator for the equations of motion in a static homogeneous magnetic field – called Cyclotronic integrator – is presented. Like Spreiter and Walter’s Taylor expansion algorithm, for sufficiently weak electric field gradients this second order method does not require a fine resolution of the Larmor motion; it has however the essential advantage of being symplectic, hence time-reversible. The Cyclotronic integrator is only subject to a linear stability constraint ([OmegaDelta t] < pi, [Omega] being the Larmor angular frequency), and is therefore particularly suitable to electrostatic Particle In Cell codes with uniform magnetic field where [Omega]is larger than any other characteristic frequency, yet a resolution of the particles’ gyromotion is required. Application examples and a detailed comparison with the well-known (time-reversible) Boris algorithm are presented; it is in particular shown that implementation of the Cyclotronic integrator in the kinetic codes SCEPTIC and Democritus can reduce the cost of orbit integration by up to a factor of ten. National Science Foundation (U.S.) and United States. Dept. of Energy (DE-FG02- 06ER54891) United States. Dept. of Energy (DE-FC02-99ER54512) 2010-09-28T12:46:50Z 2010-09-28T12:46:50Z 2008-12 2008-10 Article http://purl.org/eprint/type/SubmittedJournalArticle 0021-9991 http://hdl.handle.net/1721.1/58724 Patacchini, L., and I.H. Hutchinson. “Explicit time-reversible orbit integration in Particle In Cell codes with static homogeneous magnetic field.” Journal of Computational Physics 228.7 (2009): 2604-2615. © 2009 Elsevier Inc. en_US http://dx.doi.org/10.1016/j.jcp.2008.12.021 Journal of Computational Physics Article is made available in accordance with the publisher's policy and may be subject to US copyright law. Please refer to the publisher's site for terms of use. application/pdf Academic Press Assistant to Ian Hutchinson |
| spellingShingle | Patacchini, Leonardo Hutchinson, Ian Horner Explicit integrators for the magnetized equations of motion in Particle in Cell codes |
| title | Explicit integrators for the magnetized equations of motion in Particle in Cell codes |
| title_full | Explicit integrators for the magnetized equations of motion in Particle in Cell codes |
| title_fullStr | Explicit integrators for the magnetized equations of motion in Particle in Cell codes |
| title_full_unstemmed | Explicit integrators for the magnetized equations of motion in Particle in Cell codes |
| title_short | Explicit integrators for the magnetized equations of motion in Particle in Cell codes |
| title_sort | explicit integrators for the magnetized equations of motion in particle in cell codes |
| url | http://hdl.handle.net/1721.1/58724 |
| work_keys_str_mv | AT patacchinileonardo explicitintegratorsforthemagnetizedequationsofmotioninparticleincellcodes AT hutchinsonianhorner explicitintegratorsforthemagnetizedequationsofmotioninparticleincellcodes AT patacchinileonardo explicittimereversibleorbitintegrationinparticleincellcodeswithstatichomogeneousmagneticfield AT hutchinsonianhorner explicittimereversibleorbitintegrationinparticleincellcodeswithstatichomogeneousmagneticfield |