專利名稱：一種預(yù)測(cè)腔體濾波器微放電閾值的方法
技術(shù)領(lǐng)域：
本發(fā)明涉及空間特殊效應(yīng)應(yīng)用領(lǐng)域，尤其涉及一種預(yù)測(cè)腔體濾波器微放電閾值的數(shù)值模擬方法。
背景技術(shù)：
微放電是真空環(huán)境下在微波器件內(nèi)由于二次電子發(fā)射產(chǎn)生的一種空間特殊效應(yīng)。 微放電發(fā)生時(shí)會(huì)帶來(lái)底噪增高、信號(hào)惡化等缺陷，甚至導(dǎo)致微波器件的永久損壞。為了在部件設(shè)計(jì)之初消除微放電隱患，提高微波部件的可靠性與使用壽命，必須預(yù)測(cè)微波部件的微放電閾值并加以避免。
尤其對(duì)于腔體濾波器而言，由于其電磁場(chǎng)分布復(fù)雜，簡(jiǎn)單的近似等效無(wú)法預(yù)測(cè)其微放電閾值。目前，預(yù)測(cè)同軸腔體濾波器與多階脊波導(dǎo)濾波器等復(fù)雜結(jié)構(gòu)腔體濾波器的微放電閾值主要通過反復(fù)進(jìn)行微放電實(shí)驗(yàn)得到，實(shí)驗(yàn)周期長(zhǎng)、研制經(jīng)費(fèi)高昂，嚴(yán)重阻礙了空間應(yīng)用腔體濾波器向更大功率容量的發(fā)展。
查閱專利數(shù)據(jù)庫(kù)與文獻(xiàn)，在數(shù)據(jù)庫(kù)中分別搜索“Multipactor”、“Multipactor Threshold”、“微放電”和“微放電閾值”,可以發(fā)現(xiàn)多篇關(guān)于微放電閾值的文章。其中，針對(duì)大功率微波部件，業(yè)內(nèi)傳統(tǒng)的微放電余量設(shè)計(jì)與閾值功率預(yù)測(cè)通?；跉W洲空間局提出的敏感曲線進(jìn)行(ESA-ESTEC, Space Engineering Multipacting Design and Test, vol. ECSS-20-01A, ESA Publication Division, the Netherlands, 2003.) 該曲線基于無(wú)窮大平行平板理想模型與均勻時(shí)諧電場(chǎng)變化假設(shè)，忽略磁場(chǎng)作用，結(jié)合微波部件最窄間距處物理尺寸與工作頻點(diǎn)，建立電子在平行平板間運(yùn)動(dòng)的解析求解方程并具體求解，通過假設(shè)功率增加至電子在最窄間距處運(yùn)動(dòng)與金屬平板發(fā)生碰撞時(shí)二次電子產(chǎn)額約等于I時(shí)將發(fā)生微放電，實(shí)現(xiàn)微波部件微放電閾值預(yù)測(cè)，與本發(fā)明方法存在本質(zhì)不同。同時(shí)該方法所作物理近似多，理論模型簡(jiǎn)單，預(yù)測(cè)精度有限，不適用于復(fù)雜結(jié)構(gòu)、強(qiáng)諧振結(jié)構(gòu)微波部件微放電閾值預(yù)測(cè)。
其他進(jìn)行大功率微波部件微放電閾值預(yù)測(cè)研究的主要研究機(jī)構(gòu)有瑞士的 Chalmers理工大學(xué),西班牙Valencia大學(xué),俄羅斯科學(xué)院應(yīng)用物理所(IAP)、通用物理所 (IGP)，美國(guó)Michigan大學(xué)。其中，涉及預(yù)測(cè)微放電閾值的研究主要有
Semenov與Sazontov等人提出的’MuLSl’軟件，P6rez與Vicente等人基于微波傳輸線研究的微放電三維數(shù)值模擬方法，以及Gusarova與Isaev等人針對(duì)非對(duì)稱射頻結(jié)構(gòu)進(jìn)行微放電數(shù)值模擬方法研究，均采用對(duì)單個(gè)初始電子的運(yùn)動(dòng)軌跡進(jìn)行追蹤，并記錄每次碰撞的二次電子產(chǎn)額，通過設(shè)置大量初始電子并對(duì)總的二次電子產(chǎn)額求平均，若平均值大于I 則認(rèn)為發(fā)生微放電，該方法與本發(fā)明方法理論模型不同，同時(shí)該方法存在一定的隨機(jī)偏差， 隨著仿真時(shí)間的改變與初始電子數(shù)的改變，微放電閾值預(yù)測(cè)值存在一定偏差。
Perez與Vicente等人基于同軸傳輸線系統(tǒng)提出了微放電閾值的預(yù)測(cè)方法,通過記錄微波部件中所有電子數(shù)目隨時(shí)間的變化進(jìn)行微放電判斷，該方法與本發(fā)明方法存在不同，同時(shí)該方法耗費(fèi)仿真時(shí)間長(zhǎng)，閾值判斷基于研究者經(jīng)驗(yàn)進(jìn)行，不能進(jìn)行自動(dòng)判斷。CN 102930100 A書明說(shuō)2/5頁(yè)
Frotanpour與Dadashzadeh等人基于橢圓波導(dǎo)提出的2維微放電數(shù)值模擬方法， 通過Monte Carlo數(shù)值計(jì)算方法與有效電子方法進(jìn)行微放電解析求解，建立單個(gè)初始電子的軌跡公式，同樣采用對(duì)單個(gè)初始電子每次碰撞二次電子產(chǎn)額求平均的方式進(jìn)行微放電閾值判定，與本發(fā)明方法存在不同。發(fā)明內(nèi)容
本發(fā)明的目的在于克服現(xiàn)有技術(shù)的上述不足，提供一種預(yù)測(cè)腔體濾波器微放電閾值的方法，該方法能夠簡(jiǎn)化腔體濾波器的微放電設(shè)計(jì)過程，縮短研制周期，大幅度降低研發(fā)成本，適合各種腔體濾波器的微放電設(shè)計(jì)。
本發(fā)明的上述目的是通過如下技術(shù)方案予以實(shí)現(xiàn)的
一種預(yù)測(cè)腔體濾波器微放電閾值的方法，包括下列步驟
a.設(shè)置輸入腔體濾波器的初始電磁信號(hào)及其頻率&、仿真時(shí)間T、時(shí)間步長(zhǎng)t0、總時(shí)間步Nt、輸入功率P0及其功率掃描范圍P1 P2, P1 ^ P0 ^ P2 ；
b.對(duì)腔體濾波器進(jìn)行三維幾何建模，設(shè)置腔體濾波器的金屬材料與金屬材料的二次電子發(fā)射理論模型，在腔體濾波器中空部分建立粒子模擬區(qū)域，將三維幾何模型與粒子模擬區(qū)域剖分成多個(gè)六面體網(wǎng)格，粒子隨機(jī)分布在粒子模擬區(qū)域；
c.對(duì)上述建立的模型按時(shí)間步迭代進(jìn)行微放電數(shù)值模擬，達(dá)到仿真時(shí)間T后結(jié)束數(shù)值模擬，然后轉(zhuǎn)入步驟d;
在每一個(gè)時(shí)間步η (I ^ n ^ Nt)進(jìn)行的微放電數(shù)值模擬包括
(i)求解時(shí)域有限差分離散形式的麥克斯韋方程得到每一個(gè)網(wǎng)格節(jié)點(diǎn)和網(wǎng)格線處的電場(chǎng)值和磁場(chǎng)值，求解時(shí)域差分離散形式的牛頓_洛倫茲方程得到每一個(gè)粒子的位移與動(dòng)量;
(ii)根據(jù)粒子的位移判斷粒子是否與腔體濾波器碰撞，若發(fā)生碰撞，根據(jù)碰撞時(shí)粒子的動(dòng)量與金屬材料的二次電子發(fā)射理論模型得到出射的二次粒子產(chǎn)額和二次粒子的能量與相位分布，累計(jì)所有發(fā)生碰撞的粒子數(shù)Mi (η)，所有出射的二次粒子的總數(shù)Ms (η)；M ( )
d.求解平均二次電子發(fā)射Ma(n)，其中=
權(quán)利要求
1.一種預(yù)測(cè)腔體濾波器微放電閾值的方法，其特征在于，包括下列步驟 a.設(shè)置輸入腔體濾波器的初始電磁信號(hào)及其頻率&、仿真時(shí)間T、時(shí)間步長(zhǎng)h、總時(shí)間步Nt、輸入功率P。及其功率掃描范圍P1 P2，P1 ^ P0 ^ P2 ; b.對(duì)腔體濾波器進(jìn)行三維幾何建模，設(shè)置腔體濾波器的金屬材料與金屬材料的二次電子發(fā)射理論模型，在腔體濾波器中空部分建立粒子模擬區(qū)域，將三維幾何模型與粒子模擬區(qū)域剖分成多個(gè)六面體網(wǎng)格，粒子隨機(jī)分布在粒子模擬區(qū)域； c.按時(shí)間步迭代進(jìn)行微放電數(shù)值模擬，達(dá)到仿真時(shí)間T后結(jié)束數(shù)值模擬，然后轉(zhuǎn)入步驟d; 在每一個(gè)時(shí)間步η (I ^ n ^ Nt)進(jìn)行的微放電數(shù)值模擬包括 (i)求解時(shí)域有限差分離散形式的麥克斯韋方程得到每一個(gè)網(wǎng)格節(jié)點(diǎn)和網(wǎng)格線處的電場(chǎng)值和磁場(chǎng)值，求解時(shí)域差分離散形式的牛頓-洛倫茲方程得到每一個(gè)粒子的位移與動(dòng)量; ( )根據(jù)粒子的位移判斷粒子是否與腔體濾波器碰撞，若發(fā)生碰撞，根據(jù)碰撞時(shí)粒子的動(dòng)量與金屬材料的二次電子發(fā)射理論模型得到出射的二次粒子產(chǎn)額和二次粒子的能量與相位分布，累計(jì)所有發(fā)生碰撞的粒子數(shù)M1 (η)，所有出射的二次粒子的總數(shù)Ms (η)； d.求解平均二次電子發(fā)射]\^(11)，其中
2.根據(jù)權(quán)利要求I所述的一種預(yù)測(cè)腔體濾波器微放電閾值的方法，其特征在于，所述第一間隔為3dB，所述第二間隔為X dB,0. I ^ O. 5。
全文摘要
本發(fā)明公開了一種預(yù)測(cè)腔體濾波器微放電閾值的方法，包括確定微放電數(shù)值模擬參數(shù)、建立腔體濾波器粒子模擬區(qū)域、進(jìn)行微放電數(shù)值模擬、記錄微放電數(shù)值模擬過程中隨時(shí)間變化的二次電子發(fā)射并進(jìn)行數(shù)值處理、微放電閾值功率自動(dòng)掃描與預(yù)測(cè)；其以腔體濾波器幾何模型與粒子模擬區(qū)域?yàn)榛A(chǔ)，結(jié)合金屬材料二次電子發(fā)射特性，利用粒子模擬技術(shù)，設(shè)計(jì)了腔體濾波器微放電數(shù)值模擬與閾值功率判斷方法，實(shí)現(xiàn)了腔體濾波器微放電閾值的準(zhǔn)確數(shù)值分析與預(yù)測(cè)，大大縮短了研制周期、節(jié)約了高昂的微放電實(shí)驗(yàn)費(fèi)用，非常適用于復(fù)雜結(jié)構(gòu)腔體濾波器微放電閾值預(yù)測(cè)與大功率微波部件微放電設(shè)計(jì)。
文檔編號(hào)G06F17/50GK102930100SQ201210429170
公開日2013年2月13日 申請(qǐng)日期2012年10月31日 優(yōu)先權(quán)日2012年10月31日
發(fā)明者李韻, 王新波, 崔萬(wàn)照 申請(qǐng)人:西安空間無(wú)線電技術(shù)研究所