Axisymmetric charge-conservative electromagnetic particle simulation algorithm on unstructured grids: Application to microwave vacuum electronic devices

JOURNAL OF COMPUTATIONAL PHYSICS

卷:346

Article

Na, Dong-Yeop^1,2  Omelchenko, Yuri A.³  Moon, Haksu⁴  Borges, Ben-Hur V.⁵  Teixeira, Fernando L.^1,2 

[1] Ohio State Univ, ElectroSci Lab, Columbus, OH 43212 USA

[2] Ohio State Univ, Dept Elect & Comp Engn, Columbus, OH 43212 USA

[3] Trinum Res Inc, San Diego, CA 92126 USA

[4] Intel Corp, Hillsboro, OR 97124 USA

[5] Univ Sao Paulo, Dept Elect & Comp Engn, BR-13566590 Sao Carlos, SP, Brazil

关键词: Electromagnetic simulation;    Particle-in-cell;    Unstructured meshes;    Vacuum electronic devices;    Backward-wave oscillator;    Cerenkov radiation;   

DOI  :  10.1016/j.jcp.2017.06.016

来源: Elsevier

【 摘 要 】

We present a charge-conservative electromagnetic particle-in-cell (EM-PIC) algorithm optimized for the analysis of vacuum electronic devices (VEDs) with cylindrical symmetry (axisymmetry). We exploit the axisymmetry present in the device geometry, fields, and sources to reduce the dimensionality of the problem from 3D to 2D. Further, we employ 'transformation optics' principles to map the original problem in polar coordinates with metric tensor diag(1,rho(2), 1) to an equivalent problem on a Cartesian metric tensor diag(1, 1, 1) with an effective (artificial) inhomogeneous medium introduced. The resulting problem in the meridian (rho z) plane is discretized using an unstructured 2D mesh considering TE phi-polarized fields. Electromagnetic field and source (node-based charges and edge-based currents) variables are expressed as differential forms of various degrees, and discretized using Whitney forms. Using leapfrog time integration, we obtain a mixed epsilon - B finite-element time-domain scheme for the full-discrete Maxwell's equations. We achieve a local and explicit time update for the field equations by employing the sparse approximate inverse (SPAI) algorithm. Interpolating field values to particles' positions for solving Newton-Lorentz equations of motion is also done via Whitney forms. Particles are advanced using the Boris algorithm with relativistic correction. A recently introduced charge-conserving scatter scheme tailored for 2D unstructured grids is used in the scatter step. The algorithm is validated considering cylindrical cavity and space-charge-limited cylindrical diode problems. We use the algorithm to investigate the physical performance of VEDs designed to harness particle bunching effects arising from the coherent (resonance) Cerenkov electron beam interactions within micro-machined slow wave structures. (C) 2017 Elsevier Inc. All rights reserved.

【 授权许可】

Free