Fast particle-in-cell simulation method of calculating the multipactor thresholds of microwave devices based on their frequency-domain EM field solutions

Abstract

In order to compute the multipactor thresholds of microwave devices with high efficiency and precision,a novel fast particle-in-cell (PIC) method is proposed,which takes advantage of the frequency-domain (FD) electromagnetic field solver of CST Microwave Studio (MWS).At the initial stage of multipactor (when there are not many electrons in the device),the self-consistent field generated by the electrons is much smaller than the applied electromagnetic field. Therefore it can be ignored in calculating the multipactor threshold and this will significantly reduce the computation burden.During simulations of multipactor process,the FD field pre-calculated by CST MWS is converted into timedomain (TD) scaling with the square root of the input power.Then the electron motion is investigated by Boris algorithm.When the electrons hit the boundaries of the simulation region,where triangular facets from CST are used for discretization,the secondary electrons will be emitted.After a series of simulations with variable input powers,the multipactor threshold is determined according to time evolution of the electron number.The multipactor thresholds in a parallel plate and a coaxial transmission line are investigated,and used as relevant verifications.Compared with the CST Particle Studio (PS),the fast method obtains almost the same thresholds,while the computational efficiency is improved by more than one order of magnitude.Since the self-consistent field generated by the electrons is ignored in the fast method and it is considered in CST PS,the results validate that the self-consistent field can be ignored in calculating the multipactor threshold.Finally,taking for example a parallel plate transmission line and a stepped impedance transformer,we study the effect of the number of initial macro-particles on the calculation precision.When the initial particles are so few that they can hardly reflect the randomness of the multipactor process,a higher calculated value will be resulted in.With the increase of the number of initial macro-particles,the calculated multipactor threshold is lower and more accurate.It is convergent when the number reaches about 2000 for the parallel plate transmission line and 4000 for the stepped impedance transformer,respectively.Taking into account other microwave devices with more complex electromagnetic field distribution,in order to ensure precision,it is recommended to select the number of initial macro-particles to be 8000.In addition,although CST MWS is used to obtain the electromagnetic field and boundary information in this paper,of course,other electromagnetic softwares (such as HFSS) can also be adopted as an alternative.