一种确定非均匀场分布多载波微放电二次电子数目的方法
【专利摘要】一种确定非均匀场分布多载波微放电二次电子数目的方法，首先对非均匀电场分布区域的横截面进行等分，对每个细分区域进行多载波二次电子统计累积，获得每个细分区域的电子数目后，以二次电子的扩散特性对不同的区间进行二次电子数目的重新分配，以分配后的二次电子数目为初始数目对每个区间再进行二次电子数目统计累积，从而获得非均匀场分布时多载波微放电二次电子数目。
【专利说明】-种确定非均匀场分布多载波微放电二次电子数目的方法
【技术领域】
[0001] 本发明涉及一种确定非均匀场分布多载波微放电二次电子数目的方法，属于微波 器件微放电领域。
【背景技术】
[0002] 目前大多数卫星都工作在多载波模式，通信卫星收发系统双工器的滤波器是信号 传输的公共通道，并且处于高谐振状态，存在多载波信号激励下微放电的分析与设计问题。
[0003] 多载波微放电分析不同于单载波情形最重要的区别在于，激励信号的幅度是随着 初始相位的不同随时间快速变化的。空间微波部件多载波微放电阈值分析需要确定能够激 励微放电的最小功率及其对应的相位分布。为了确定该相位分布需要进行大量的计算，数 值模拟方法由于计算速度慢无法解决这一问题。因此国际上提出采用概率统计方法进行多 载波微放电分析，大大节省了时间。
[0004] 但目前只能针对平板结构、电场均匀分布的情形进行计算，对于非均匀场条件下 无法考虑不同电子的扩散。因此对于非均匀场分布的微波部件无法计算。
【发明内容】
[0005] 本发明的技术解决问题是：克服现有技术的不足，提供了一种确定非均匀场分布 多载波微放电二次电子数目的方法，可以对上下两极板形状相同，对应的电场沿横截面的 一个方向成非均匀分布，而另外一个方向均匀分布的微波部件获得多载波合成信号激励条 件下的二次电子数目特性。
[0006] 本发明的技术解决方案是：一种确定非均匀场分布多载波微放电二次电子数目的 方法，步骤如下： η
[0007] ⑴令多载波信号为卬卜乙匕奶^^^/:广抑山其中心心…忑为载波频 率，奶，於，…，％为载波相位，Vi，V2，…，Vn为载波幅值，η为载波路数，t为仿真时间， 0 < t < T，T为预设的仿真时长；根据微波部件材料确定部件表面的二次电子发射特性；
[0008] (2)对待分析的微波部件施加外加电场，并在与外加电场垂直的方向对待分析的 微波部件横截面进行Num等分，每个等分区域记为Py i = 1. . . Num，在与外加电场垂直的方 向，任意两个等分区域中心之间的距离为Dij，i, j = 1. . . Num ;
[0009] (3)对预设的仿真时长T进行A等分，每一等分的时间长度记为扩散步长ΔΤ，把 每一个Δ T进行M2等分，等分后的时间步长记为累计步长Δ t，即T = 1^* Δ T = MfM# Δ t ;
[0010] ⑷通过CST软件或解析计算得到当多载波信号为VOi^AT+n^At)时，步骤（2) 中得到的微波部件各等分区域Pi的平均电压Vi (11^* Δ T+m2* Δ t)，i = 1. · · Num，!!^ = 0· · · Μρ m2 = 0· · · M2 ;
[0011] (5)令!^ = 0,每个等分区域Pi的初始电子数目为乂，每个等分区域Pi临时电子 数目为Xi，i = 1. . . Num，仿真初始时刻令Xi = NQ, i = 1. . . Num ;
[0012] (6)在叫+八!'到（mi+l)*AT时间段内，以At为步长，步骤（4)中得到的平均电 压ViOi^AT+n^At)为激励，以发生微放电微波部件上下极板之间的距离为L的平行板 为对象，Xi为初始电子数目，根据步骤（1)提供的二次电子发射特性，采用二次电子概率 统计方法计算步骤（2)中得到的微波部件各等分区域？ 1的电子数目NiOi^AT+n^At)，i =1. . .Num, m2 = 0. . .M2，最终得到（i^+lhAT时亥lj，微波部件各等分区域Pi的电子数目 ^((m^D^AT)；
[0013] (7)对步骤（6)得到的微波部件各等分区域Pi的电子数目NiUnii+lhAT)进行重 新分配获得((川i +1ΓΔ?,并令4 = ((川i +1)*ΔΓ);
[0014] (8)令叫=π^+1,判断叫是否等于％，若叫不等于％，则执行步骤（6)，若叫等于 Mi，则结束计算，获得电子数目随时间变化曲线队(〇,0彡t彡T。
[0015] 所述步骤（7)中对步骤（6)得到的微波部件各等分区域Pi的电子数目 队（(mi+l) * Λ T)进行重新分配获得#,.((岬+1)*ΔΓ)，并令4 = +1)*ΔΓ)，具体步骤如 下：
[0016] (i)计算 Δ 队』=队（(mi+Ι) * Δ T) *p (Du)，i = 1. · · Num，j = 1. · · Num，j 尹 i，其中 Λ &为从第i个等分区域Pi向第j个等分区域&扩散的电子数目，p (Dy为电子扩散因 子； Num Num _7] (ii)计算 Σ^+Σ^,:，将调整 M i'=i 后的电子数目々,((叫+ 1广ΔΓ)作为下一次迭代的初始电子数目进行后续计算，即 Zf=i$((，Wl+l)*AD。
[0018] 所述电子扩散因子ρΦμ)由公式：
[0019] p{Dij) = e-D-'/2/}nbin
[0020] 给出；其中，β = 0. 8, 为在与外加电场垂直的方向，任意两个等分区域中心之 间的距离，n = AT*f\。
[0021] 所述载波路数n>2 ;所述等分数目M，2 ;所述At〈lAlOO*fn);所述初始电子数目 N。为大于0的整数。
[0022] 本发明与现有技术相比的有益效果为：
[0023] (1)本发明方法以概率统计计算方法为基础，把非均匀区域进行细分，考虑了不 同区域之间的电子扩散，有效解决了非均匀场条件下多载波微放电二次电子数目计算的难 题，计算速度快，为多载波微放电分析提供了有效手段；
[0024] (2)本发明提供了确定多载波激励条件下电子数目的方法，不借助于数值分析过 程，不受计算机硬件导致的粒子数目的限制，对于单个合成周期内放电的情形能够计算，有 效克服了数值分析方法的缺陷；
[0025] (3)本发明通过把电场非均匀的区域进行等分，对每个等分的区域采用概率统计 方法进行处理，同时考虑不同区域之间电子的扩散，通过等分实现了非均匀场分布条件下 二次电子的计算，有效解决了非均匀场分布条件下微放电最坏状态分析中电子数目确定的 问题，有着实际应用价值。
【专利附图】
【附图说明】
[0026] 图1为本发明流程图；
[0027] 图2为本发明实施例中的仿真对象示意图；
[0028] 图3为本发明实施例中电场非均匀区域的分解图；
[0029] 图4为本发明实施例中多载波信号激励时获得的二次电子数目。
【具体实施方式】
[0030] 下面结合附图对本发明的【具体实施方式】进行详细的阐述。
[0031] 如图1所示为本发明的流程图，由图1可知，本发明提供的一种确定非均匀场分布 多载波微放电二次电子数目的方法，主要步骤如下： η
[0032] (1)令多载波信号为ΖΟΣ匕cos(2：l!#// i + 6)，其中f\，f2，…，4为载波 Μ 频率，仍，外，…，％为载波相位，Vi，V2，…，Vn为载波幅值，η为载波路数，n>2;t为仿 真时间，0 < t < T ;T为预设的仿真时长；同时令待分析微波部件发生微放电的上下极板之 间的距离为L，设定部件表面的二次电子发射特性，采用Vaughan模型设定部件表面的二次 电子发射特性；对于银而言二次电子发射特性设定四个值洱=30eV，E max = 165eV，E2 = 5000eV，δ max = 2. 22其中Ei为二次电子发射系数等于1对应的低端能量值，Emax为最大电 子发射系数对应的能量值，E 2为二次电子发射系数等于1对应的上端能量值，为最大电子发 射系数；
[0033] (2)对待分析的微波部件施加外加电场，并在与外加电场垂直的方向对器件横截 面进行Num等分，每个等分区域记为Py i = 1. . . Num，在与外加电场垂直的方向，任意两个 等分区域Pi中心之间的距离为Dm i，j = 1. . . Num ;
[0034] (3)对预设的仿真时长T进行A等分，Mi>2,每一等分的时间长度记为扩散步长 ΛΤ，把每一个ΛΤ进行M2等分，等分后的时间步长记为累计步长At，At〈lAl00*fn)，即 T = M# Δ T = M^Mg* Δ t ;
[0035] (4)通过通用电磁仿真软件CST或解析计算获得微波部件对应于 V(mi*AT+m2*At)时，步骤⑵中得到的微波部件各等分区域Pi的平均电压 (π^* Δ T+m2* Δ t)，i = 1. · · Num，1? = 0· · · Ml，m2 = 0· · · M2,在 CST 软件中米用 PSstudio η 工作室，建立待分析部件的040模型，设定激励信号为》/(〇 =乙【/,叫彳2$足：1!/广<+^/!)，设 Μ 定总的仿真时间为Τ，在每个等分区域Pi的中央加入电压监视器，并进行时域仿真，可记录 每个等分区域Pi的平均电压Vi (nii* Δ T+m2* Δ t);
[0036] (5)令叫=0,每个等分区域Pi的初始电子数目为NQ，NQ为大于0的整数，每个等 分区域 Pi 临时电子数目为 Xi，i = 1··· Num，X。= N。，i = 1··· Num ;
[0037] (6)在-八!1到（n^+lhAT时间段内，以At为步长，平均电压ViOi^AT+n^At) 为激励，上下极板之间的距离为L的平行板为对象，\为初始电子数目，根据步骤（1)提供 的二次电子发射特性，采用二次电子概率统计方法计算步骤（2)中得到的微波部件各等分 区域的电子数目队(nii* Δ T+m2* Δ t)，i = 1. · · Num，m2 = 0· · · M2 ;最终得到（mi+Ι) * Δ T 时 亥lj，微波部件各等分区域Pi的电子数目NjOiii+lhAT);所述二次电子概率统计方法详见 "Multipactor theory for multicarrier signals, Phys. Plasmas, vol. 18, no. 3, p. 0321 05, 2011.第 7-8 页公式（27)-公式（32);
[0038] (7)对每个等分区域？1的电子数目NjOiii+lhAT)进行重新分配获得 々,.((^ + 1)1:0,具体步骤如下：
[0039] (7. 1)计算 Δ N。. = & ((mdl) * Δ T) *p (Du)，（i = 1. · · Num)，（ j = 1. . . Num, j尹i)，△ \代表第i个等分区域Pi向第j个等分区域Pj扩散的电子数目，其中 ρ(Ο^) = β-°·'/2βη1--η ;
[0040] 其中，β = 0. 8, 在与外加电场垂直的方向，为任意两个等分区域中心之间的距 离，η为二次电子和器件金属表面碰撞的次数；η = Λ T*f\ ; Num Num
[0041 ] (7. 2)取/Wi +1：) *ΔΓ)=顺叫 +1) * ΔΓ) - Σ W + Σ W，将调整后的电子数目 j=i 7=1 Μ 作为下一次迭代的初始电子数目进行后续计算，即l = M((叫+1)"ΔΓ);
[0042] (8)令叫=1^+1，判断1?是否等于％，若1?不等于％，则执行步骤（6)，若1?等于 Mi，则结束计算，获得电子数目随时间变化曲线队(〇,0彡t彡T。
[0043] 实施例
[0044] (1)如图2所示为本发明实施例中的仿真对象示意图，针对矩形波导，发生微放电 的两极板为y = 〇mm和y = 〇. 43mm, X方向为宽边,长度为58. 17mm, y方向为窄边,长度为 0. 43mm，采用η = 6载波进行实施，6路载波频率分别为= 3. 57GHz，f2 = 3. 67GHz，f3 = 3. 77GHz，f4 = 3. 87GHz，f5 = 3. 97GHz 和 f6 = 4. 07GHz，载波幅度 ν! = V2 = V3 = V4 = V5 = 10V，6路载波的初始相位分别为奶=0",灼=0",约=〇(l，灼=O'1，灼=O' 1，％ = 〇"， 预设的仿真时长T = 40ns，求和获得多载波合成信号，微波部件发生微放电的上下极板之 间距离为0. 43mm，部件表面为银；
[0045] (2)对待分析的微波部件施加外加电场，如图3所示为本发明实施例中电场非均 匀区域的分解图，由图3可知，外加电场方向为y方向，在与外加电场方向垂直的X方向进 行细分，把矩形波导横截面进行Num = 11等分；
[0046] (3)对预设的仿真时长T进行A = 4等分，每一等分的时间长度记为扩散步长ΔΤ =10ns，把每一个Δ T进行M2 = 200等分，等分后的时间步长记为累计步长Δ t = 0. 05ns ;
[0047] (4)通过解析计算得到当多载波信号为VOi^AT+n^At)时，步骤⑵中得到的 微波部件各等分区域Pi的平均电压Vi (m# Δ T+m2* Δ t)，i = 1. · · Num，1? = 0· · · Mp m2 = 0. . . M2 ；
[0048] (5)令叫=0,每个等分区域Pi的初始电子数目为NQ = 20,每个等分区域Pi临时 电子数目为Xi，i = 1…Num，仿真初始时刻令X。= N。= 20, i = 1··· Num ;
[0049] (6)在叫水八!'到（n^+lhAT时间段内，以At为步长，步骤（4)中得到的平均电 压Vi (mi* Λ T+m2* At)为激励，以发生微放电上下极板之间的距离为L的平行板为对象，Xi 为初始电子数目，根据步骤（1)提供的二次电子发射特性，采用二次电子概率统计方法计 算步骤⑵中得到的微波部件各等分区域Pi的电子数目队(nii* Δ T+m2* Δ t)，i = 1. . . Num, m2 = 0· · · M2 ;
[0050] (7)对步骤（6)得到的微波部件各等分区域Pi的电子数目^((π^ΙΚΛΤ)进行重 新分配获得々,.((％+1)*ΔΠ ;
[0051] (7. 1)计算 Δ \ =队（(ii^+l) * Δ T) *p (D。·)，（i = 1. · · Num)，（j = 1. · · Num, j 尹 i)， 其中Λ &为从第i个等分区域Pi向第j个等分区域&扩散的电子数目，p(Dj为电子扩 散因子；其中户(1\.) = e_Ar/2/?"lM，其中 β = 〇· 8, n = AT*f\ = 35. 7 ; Num Num
[0052] (7.2)计算死((岬+1)*从)= ]^伽1+1)：^71)-艺應</ +艺餅力_,将调整 i'=i j幸i j封 后的电子数目M ((M + ip △「）作为下一次迭代的初始电子数目进行后续计算，即 ((州 Ι+1)*Δ；Γ);
[0053] (8)令1? = 1^+1，判断1?是否等于％，若1?不等于％，则执行步骤（6)，若1?等于 Mi，则结束计算，获得电子数目随时间变化曲线队⑴，0 < t < T，如图4所示为本发明实施 例中多载波信号激励时获得的二次电子数目，图4中X轴为仿真时间，y轴为电子数目，可 以看出电子数目随着多载波信号波动，并呈增长趋势，将会产生微放电。
[0054] 本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
【权利要求】
1. 一种确定非均匀场分布多载波微放电二次电子数目的方法，其特征在于步骤如下： ⑴令多载波信号为HO = ，其中f\，f2，…，fn为载波频率， Μ 约，朽，…，％为载波相位，Vi，V2，…，Vn为载波幅值，η为载波路数，t为仿真时间， 0 < t < T，T为预设的仿真时长；根据微波部件材料确定部件表面的二次电子发射特性； (2) 对待分析的微波部件施加外加电场，并在与外加电场垂直的方向对待分析的微波 部件横截面进行Num等分，每个等分区域记为Pi, i = 1. . .Num，在与外加电场垂直的方向， 任意两个等分区域中心之间的距离为D# i，j = 1. . . Num ; (3) 对预设的仿真时长T进行A等分，每一等分的时间长度记为扩散步长ΔΤ，把每一 个Δ T进行M2等分，等分后的时间步长记为累计步长Δ t，即T = 1^* Δ T = MfM# Δ t ; (4) 通过CST软件或解析计算得到当多载波信号为VOi^AT+n^At)时，步骤（2)中 得到的微波部件各等分区域Pi的平均电压Vi (11^* Δ T+m2* Δ t)，i = 1. · · Num，1? = 0· · · Μρ m2 = 0. . . M2 ； (5) 令!^ = 0,每个等分区域Ρ,的初始电子数目为队，每个等分区域？，临时电子数目 为 Xi，i = 1. . . Num，仿真初始时刻令 Xi = NQ, i = 1. . . Num ; (6) 在^八!"到（mi+l)*AT时间段内，以At为步长，步骤（4)中得到的平均电压 ViOi^AT+n^At)为激励，以发生微放电微波部件上下极板之间的距离为L的平行板为 对象，\为初始电子数目，根据步骤（1)提供的二次电子发射特性，采用二次电子概率统 计方法计算步骤（2)中得到的微波部件各等分区域？ 1的电子数目NiOi^AT+n^At)，i =1. . .Num, m2 = 0. . .M2，最终得到（i^+lhAT时亥lj，微波部件各等分区域Pi的电子数目 ^((m^D^AT)； (7) 对步骤（6)得到的微波部件各等分区域Pi的电子数目NjOiii+lhAT)进行重新分 配获得&.((% +丨)*ΔΓ)，并令不=々/((川1 +丨)* ΔΓ); (8) 令叫=1^+1，判断叫是否等于％，若叫不等于％，则执行步骤^，若叫等于札， 则结束计算，获得电子数目随时间变化曲线队(〇,0 < t < T。
2. 根据权利要求1所述的一种确定非均勻场分布多载波微放电二次电子数目的方 法，其特征在于：所述步骤（7)中对步骤（6)得到的微波部件各等分区域Pi的电子数目 NjOiii+lhAT)进行重新分配获得;?,((^+1)=^:0 ,并令义,. = &((/?,+1)*ΔΠ ,具体步骤 如下： (i)计算 Δ \ =队（(ii^+l) * Δ T) *p (Dy，i = 1. · · Num, j = 1. · · Num, j 关 i,其中 Δ Νυ 为从第i个等分区域Pi向第j个等分区域h扩散的电子数目，P (Dj为电子扩散因子； Num Num 出）计算灰((讲1+1)*么71)=%((所 1+1)*么71)-艺从^+￡从^,将调整后的电子数目 i=i 况((Wl +1) * ΔΓ)作为下一次迭代的初始电子数目进行后续计算，即A =灰((叫+1) * ΔΓ)。
3. 根据权利要求2所述的一种确定非均勻场分布多载波微放电二次电子数目的方法， 其特征在于：所述电子扩散因子P(Dj由公式 ： p(D") =，W 给出；其中，β =0.8, 为在与外加电场垂直的方向，任意两个等分区域中心之间的 距离，n = AT*f\。
4.根据权利要求1所述的一种确定非均勻场分布多载波微放电二次电子数目的方法， 其特征在于：所述载波路数n>2 ;所述等分数目MP2 ;所述△Κ?ΛΚΚ^?；);所述初始电子 数目队为大于0的整数。
【文档编号】G06F19/00GK104091050SQ201410286017
【公开日】2014年10月8日 申请日期:2014年6月24日 优先权日:2014年6月24日
【发明者】王新波, 张洪太, 李韵, 崔万照, 冉立新 申请人:西安空间无线电技术研究所