Explicit time-reversible orbit integration in Particle In Cell codes with static homogeneous magnetic field

Patacchini, Leonardo; Hutchinson, I. H.

doi:10.1016/j.jcp.2008.12.021

Cited by 18 publications

(23 citation statements)

References 15 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…We note that all of these substeps are time symmetric, i.e., ϕ T −h = (ϕ T h ) −1 , and so on. We also note that the subflows (12) and (14) . Its determinant is given by det ∂y 1 ∂y 0 = det exp hΩ( q) .…”

Section: Splitting Methodsmentioning

confidence: 81%

Splitting methods for time integration of trajectories in combined electric and magnetic fields

et al. 2015

View full text Add to dashboard Cite

The equations of motion of a single particle subject to an arbitrary electric and a static magnetic field form a Poisson system. We present a second-order time integration method which preserves well the Poisson structure and compare it to commonly used algorithms, such as the Boris scheme. All the methods are represented in a general framework of splitting methods. We use the so-called φ functions, which give efficient ways for both analyzing and implementing the algorithms. Numerical experiments show an excellent long term stability for the new method considered.

show abstract

Section: Splitting Methodsmentioning

confidence: 81%

Splitting methods for time integration of trajectories in combined electric and magnetic fields

et al. 2015

View full text Add to dashboard Cite

show abstract

“…Code operation in this regime is deferred to future work. The n part particles representing ions are advanced in Cartesian coordinates using the Cyclotronic integration scheme [19], in the frame moving with velocity v ⊥ where E conv vanishes. This enables us to use longer time-steps as the strong convective acceleration need not be resolved.…”

Section: Code Operationmentioning

confidence: 99%

Spherical probes at ion saturation inE×Bfields

Patacchini

Hutchinson

2010

Plasma Phys. Control. Fusion

View full text Add to dashboard Cite

The ion saturation current to a spherical probe in the entire range of ion magnetization is computed with SCEPTIC3D, new three-dimensional version of the kinetic code SCEPTIC designed to study transverse plasma flows. Results are compared with prior two-dimensional calculations valid in the magnetic-free regime [I.H. Hutchinson PPCF 44:1953], and with recent semi-analytic solutions to the strongly magnetized transverse Mach probe problem [L. Patacchini and I.H. Hutchinson PRE ???]. At intermediate magnetization (ion Larmor radius close to the probe radius) the plasma density profiles show a complex three-dimensional structure that SCEPTIC3D can fully resolve, and contrary to intuition the ion current peaks provided the ion temperature is low enough. Our results are conveniently condensed in a single factor M c , function of ion temperature and magnetic field only, providing the theoretical calibration for a transverse Mach probe with four electrodes placed at 45 o to the magnetic field in a plane of flow and magnetic field.

show abstract

“…[15]. A uniform external magnetic field B and a uniform convective electric field E conv drive a background cross-field drift v ⊥ = E conv ×B/B 2 ; hence in Eqs (1,2) Φ is the dust-induced potential perturbation, and E = E conv −∇Φ is the total electric field acting on the ions and electrons.…”

Section: Methods For Force Calculations 21 Sceptic3dmentioning

confidence: 99%

Forces on a spherical conducting particle inE×Bfields

Patacchini¹,

Hutchinson²

2011

Plasma Phys. Control. Fusion

View full text Add to dashboard Cite

The forces acting on a spherical conducting particle in a transversely flowing magnetized plasma are calculated in the entire range of magnetization and Debye length, using the particle code SCEPTIC3D [Patacchini and Hutchinson, Plasma Phys. Control. Fusion 52 035005 (2010), Plasma Phys. Control. Fusion 53 025005 (2011)]. In short Debye length (i.e. high density) plasmas, both the ion-drag and Lorentz force arising from currents circulating inside the dust show strong components antiparallel to the convective electric field, suggesting that a free dust particle should gyrate faster than what predicted by its Larmor frequency. In intermediate to large Debye length conditions, by a downstream depletion effect already reported in unmagnetized strongly collisional regimes, the ion-drag in the direction of transverse flow can become negative. The internal Lorentz force however remains in the flow direction, and large enough in magnitude so that no spontaneous dust motion should occur.

show abstract

Explicit time-reversible orbit integration in Particle In Cell codes with static homogeneous magnetic field

Cited by 18 publications

References 15 publications

Splitting methods for time integration of trajectories in combined electric and magnetic fields

Splitting methods for time integration of trajectories in combined electric and magnetic fields

Spherical probes at ion saturation inE×Bfields

Forces on a spherical conducting particle inE×Bfields

Contact Info

Product

Resources

About