一种感应电压叠加器次级磁绝缘电子鞘层边界的确定方法
【专利摘要】本发明涉及一种感应电压叠加器次级磁绝缘电子鞘层边界的确定方法包括，(1)建立磁绝缘电子鞘层边界满足的数学物理方程，推导电子鞘层边界内/外侧相对论因子的通用解析表达式。(2)由实验测量阴极角向磁场空间分布，结合数值求解方法，获得通解表达式中的待定系数。(3)假定电子鞘层边界满足几种给定的分布形式(如余弦分布、二次抛物线分布、高斯分布等)，确定各假定分布中的待定参数和计算误差，误差最小的分布视为电子鞘层边界的真实分布。相对于现有粒子模拟确定方法，本方法回避了目前国内三维粒子模拟软件及代码不成熟的技术问题，具有准确性高、耗时短的优点。
【专利说明】
一种感应电压叠加器次级磁绝缘电子鞘层边界的确定方法
技术领域
[0001] 本发明提出一种感应电压叠加器次级磁绝缘电子鞘层边界的确定方法，在感应电 压叠加器、真空磁绝缘传输线等脉冲功率装置中具有重要应用。
【背景技术】
[0002] 感应电压叠加器（induction voltage adders，IVA)是一种强流脉冲功率驱动源 技术，可产生电压几 MV~几十MV，电流几十~数百kA的高功率电脉冲，在γ射线辐射效应模 拟、高能脉冲X射线闪光照相等领域具有重要应用。从本质上讲，IVA可视为多个变比1:1的 脉冲变压器初级并联，次级线圈串联。每级感应腔相当于脉冲变压器的初级线圈，一根贯串 各级感应腔轴心的中心内筒视为次级线圈。为了提高功率传输密度和输出电压，IVA次级通 常工作于磁绝缘模式。
[0003] 脉冲电压通常是在IVA感应腔侧壁单端口或多端口馈入，导致馈入脉冲电流在感 应腔初级角向非均匀分布。特别地，当感应腔仅有一路脉冲电压单端口注入时，脉冲电流空 间分布不均匀程度尤为严重。感应腔初级馈入电流角向非均匀分布，相当于脉冲变压器非 均匀激磁，将影响IVA次级磁绝缘性能，如引起次级阴极表面磁绝缘电子非均匀发射、阴一 阳电极间电/磁场非均匀分布、电子鞘层边界(磁绝缘电子区域和真空区域分界面)偏心等 非对称特性。目前的磁绝缘模型（1975年Creedon提出的层流理论，1979年Mendel提出的任 意动量理论，1990年Lawconnell提出的通用Laminar理论，2006年Ottinger基于粒子模拟的 修订Mendel理论，以及2006年Stygar提出的考虑磁绝缘电子弹性碰撞理论)均用来描述一 维、轴对称磁绝缘问题，无法解决因馈入电流非均匀分布引起的IVA次级非对称磁绝缘问 题。与通常的一维、轴对称磁绝缘不同，IVA次级磁绝缘电子鞘层边界r s(0)不再是半径固定 的同心圆(rs(0) =常数），而是随角向位置呈某种分布rs(0)=f(0)。
[0004] 确定磁绝缘电子鞘层边界是评估IVA馈入电流非均匀分布、研究IVA次级功率传输 特性的关键环节。电子鞘层边界限定了磁绝缘空间电子活动区域，实际上可视作磁绝缘传 输线的虚阴极，电子鞘层边界(虚阴极)和阳极构成的传输线阻抗才是磁绝缘传输线的实际 运行阻抗或流阻抗。电子鞘层边界和阴极之间电子运动形成了磁绝缘鞘层电流，多数负载 情况下，鞘层电流是无用的。IVA次级磁绝缘线也一直着力降低鞘层电流。确定电子鞘层边 界就能够确定磁绝缘鞘层电流大小和电荷质心位置，进而确定磁绝缘运行阻抗或流阻抗。
[0005] 馈入电流非均匀分布引起IVA次级电子鞘层偏心，靠近馈入点的角向位置，电子鞘 层较薄;远离馈入点的位置，电子鞘层较厚。极端情况下，鞘层较厚位置的磁绝缘电子可能 达到阳极，造成磁绝缘阴-阳间隙短路。为了定量表征注入电流非均匀分布对磁绝缘电子鞘 层偏心程度和电流损失之间的影响程度，也必须确定电子鞘层边界。
[0006] 由于磁绝缘电子鞘层边界无法直接测量，目前国内外通常采用粒子模拟(PIC)方 法，通过跟踪大量磁绝缘电子的运动轨迹来确定电子鞘层边界，但粒子模拟计算量大、耗时 长、随机性强，其计算精度取决于网格剖分的疏密。另外，对于类似IVA次级非对称磁绝缘问 题，电子鞘层边界的确定需要采用三维PIC模拟，这对PIC模拟软件及代码提出了较高要求。 国外的大型三维粒子模拟软件如LSP，一直对我国禁运，国内的三维PIC模拟软件和代码尚 处于研制和测试阶段，距成熟应用尚有一定差距。因此，有必要探索一种实验间接测量和物 理模型相结合的新型磁绝缘电子鞘层确定方法。
【发明内容】
[0007] 为了准确评估馈入电流非均匀分布对IVA次级磁绝缘特性(电子鞘层偏心、电流损 失等)的影响，本发明提出一种实验间接测量和物理模型相结合的IVA次级磁绝缘电子鞘层 边界的确定方法。
[0008] 本发明的技术解决方案如下：
[0009] 本发明提供一种感应电压叠加器次级磁绝缘电子鞘层边界的确定方法，其特殊之 处在于:包括以下步骤：
[0010] 1)建立IVA次级磁绝缘电子鞘层边界rs(0)满足的数理方程，并明确数理方程的待 确定常数，具体是：
[0011] 在磁绝缘电子鞘层边界以0)上，磁绝缘电子区域的相对论因子γι(Γ，θ)和真空区 域相对论因子γ 〇(r，Θ)自身及其法向导数均连续，即
[0014] 其中：ys为电子鞘层边界^(0)上的相对论因子；
[0015] 丫4^0)和7。(^0)表达式如公式(2-&)和(2-13)所示 ：
[0018]其中，为IVA次级阴极半径，ra为阳极半径；
[0019] γ 〇= Ι+eVo/mc2为阳极相对论因子，Vo为次级阴-阳极间隙电压，m为电子质量，m~ 9.11X10-31kg，c 为光速，c = 3X108m/s;
[0020] η为正整数，表征阴极表面角向磁场Βθ (rc，Θ)角向分布的模数，η = 0时，Βθ (rc，Θ)角 向均匀分布，即一维轴对称磁绝缘；
[0021] 0彡j彡n，j为整数；
[0022] a」、b伪待定常数，a煎决于阴极表面角向磁场Be(rc，0)空间分布，Βθ(η，θ)空间分 布则取决于感应腔初级馈入电流空间分布；
[0023] 2)确定数理方程的待确定常数^及、
[0024] 2.1)aj 的确定
[0025] 2.1.1)测量IVA次级阴极角向磁场Be(rc，Θ)空间分布，将不同角向位置的实测Be (η，θ)数值进行拟合，以获得阴极角向磁场说(^，0)满足分布Β(θ);
[0026] 2.1.2)在[0，2jt ]区间内对阴极角向磁场Be (r c，Θ)满足的分布Β (Θ)作Four i er级数 展开，得到系数
[0029] 2.2)bj 的确定
[0030] 2.2.1)根据IVA次级磁绝缘电子鞘层边界的物理内涵和可能分布，假定电子鞘层 边界满足某种特定的函数分布f〇(9);
[0031] 2.2.2)离散化数理方程，求取匕与叫及假定电子鞘层边界?·〇(θ)的关系；
[0032] 2.2.3)最优化问题确定假定分布函数fo(0)中的待定参数；
[0034] f(0) s.t. {rc^f(9)^ra
[0035] 2.2.4)利用假定的电子鞘层边界函数分布替代建立的电子鞘层边界数理方程求 解出bj;
[0036] 2.2.5)求取电子鞘层边界数理方程与假定的电子鞘层边界函数分布之间的误差；
[0037] 3)电子鞘层边界函数的确定
[0038] 3.1)根据IVA次级磁绝缘电子鞘层边界的物理内涵和可能分布，假定电子鞘层边 界满足其它类型的函数分布fo(9);
[0039] 3.2)重复步骤2.2.2-2.2.5);
[0040] 3.3)比较多种假定的电子鞘层边界函数与所建立的电子鞘层边界数理方程之间 的误差；
[0041] 3.4)选择误差最小的假定电子鞘层边界函数为最终的电子鞘层边界。
[0042] 上述2.2.2)具体为：
[0045]将公式(4)代入方程(lb)，得到方程(5)
[0047]将方程(5)等号两边均写成级数求和形成
[0051]因方程(6)对于[0,2π]内任意Θ均成立，因此将Θ在[0,231]取(n+1)个离散点
[0053] 由于对于任意一个均满足方程(6)，可以得到
[0054]
[0055] 将公式(9)写成矩阵方程的形式，BP
[0056] Na=Mb (10)
[0057] 其中，a=(a〇，ai,…，an)T和b = (b〇，bi,…，bn)T为n+1阶向量；
[0058] N、M为(n+1) X (n+1)阶矩阵，其元素分别为：
[0059]
[0060]上述步骤2.1.2)是采用正交的三角函数系将由实验数据拟合的分布Β(θ)，在[0,2 π ]区间内作Four i er级数展开，得到系数aj的。
[0061 ]上述步骤2.2.3)的最优化问题采用全局搜索解法进行数值求解，以确定假定分布 函数中的待定常数，使各假定分布函数本身的计算误差最小；
[0062]所述全局搜索解法包括遗传算法、模拟退化方法或蚁群算法。
[0063]由于数值求解方法2.2.2和2.2.3是由数理方程lb推导而来，因此计算出的假定分 布f〇( θ)是否满足数理方程la需要进行校核；
[0064]另外，步骤2.2.3中求解假定分布函数中的待定常数时，实际上将上述定义计算误 差σ作为最优化问题的目标函数，采取全局搜索计算方法，寻找待定常数，使目标函数〇最 小；
[0065]为了定量描述这种电子鞘层边界确定方法的精度，定义计算误差为
[0067] 计算误差实际上表征在假定电子鞘层边界分布函数fQ(0)上，磁绝缘电子区域γι (?·〇(θ)，θ)和真空区域相对论因子YQ(fQ(0)，0)之间的误差。
[0068] 本发明与现有技术相比，有益效果是：
[0069] 1、本发明提出一种确定感应电压叠加器次级磁绝缘电子鞘层边界的方法，通过实 验间接测量阴极角向磁场空间分布，结合数学物理方程的数值算法，可确定非对称磁绝缘 (馈入电流非均匀分布引起的IVA次级磁绝缘、长距离磁绝缘线内筒偏心等）的电子鞘层边 界。相对于现有粒子模拟确定方法，本方法回避了目前国内三维粒子模拟软件及代码不成 熟的技术问题。本方法的准确性高、耗时短。
[0070] 2、本发明提出的方法，将确定电子鞘层边界问题转换为最优化求解问题。假定电 子鞘层满足某种分布，采用目前成熟的最优化问题数值求解算法，确定假定分布中的待定 参数，使目标函数(计算误差)最小，并给出数值计算误差。
[0071] 3、本发明提出的方法，可以界定IVA次级非对称磁绝缘空间电子的活动区域，获得 磁绝缘电子鞘层边界的偏心程度，对于研究非对称磁绝缘特性具有重要应用价值。
【附图说明】
[0072]图1为IVA次级非对称磁绝缘示意图。
[0073]图2由实验测量拟合得到Be(，Θ)分布。
[0074]图3a为三种形式的磁绝缘电子鞘层边界。
[0075]图3b极坐标下三种电子鞘层边界分布。
[0076] 图4磁绝缘空间电子二维分布。
[0077 ]图5阴极r。、阳极r4P电子鞘层边界处1^ (Θ)角向磁场
[0078] r。为IVA次级中心内筒的阴极半径，ra为次级导体的阳极半径。rc；〈r〈 rs (Θ)为磁绝 缘电子鞘层区域，rs(0)〈r〈ra为真空区域(无磁绝缘空间电子）， rs(0)为待确定的电子鞘层 边界。
【具体实施方式】
[0079] 本发明提出一种感应电压叠加器次级磁绝缘电子鞘层边界的确定方法，核心思想 是：（1)建立磁绝缘电子鞘层边界满足的数学物理方程，推导电子鞘层边界内/外侧相对论 因子的通用解析表达式。（2)由实验测量阴极角向磁场空间分布，结合数值求解方法，获得 通解表达式中的待定系数。（3)假定电子鞘层边界满足几种给定的分布形式(如余弦分布、 二次抛物线分布、高斯分布等），确定各假定分布中的待定参数和计算误差，误差最小的分 布视为电子鞘层边界的真实分布。
[0080] 以下结合附图，举例对本发明进行说明。
[0081] 感应电压叠加器给定参数：阴极半径rc = 0. lm，阳极半径^ = 0.2m，阴-阳极电压Vo = 4MV，γ Q=9<JVA次级工作在最小磁绝缘电流工作点。
[0082] 该参数下次级磁绝缘电子鞘层边界的确定方法主要包括以下步骤：
[0083] 1)建立IVA次级磁绝缘电子鞘层边界rs(0)满足的数理方程，并明确数理方程的待 确定常数，具体是：
[0084] 在磁绝缘电子鞘层边界^(0)上，磁绝缘电子区域的相对论因子Yl(r，0)和真空区 域相对论因子γ 〇(r，Θ)自身及其法向导数均连续，即
[0087] 其中：ys为电子鞘层边界^(0)上的相对论因子。
[0088] 丫410)和丫。(10)表达式如公式(2-&)和(2-13)所示。
[0091 ] 其中，r。为IVA次级阴极半径，ra为阳极半径。γ 〇=l+eV〇/mc2为阳极相对论因子，V〇 为次级阴-阳极间隙电压，m为电子质量，m~9.11 X l(T31kg，c为光速，c = 3X 108m/s。]!为非 负整数，表征阴极表面角向磁场Βθ (r。，Θ)角向分布的模数，η = 0时，Βθ (r。，Θ)角向均勾分布， 即一维轴对称磁绝缘。0<j<n，j为整数为待定常数。a诹决于阴极角向磁场Βθ(η，θ) 空间分布，Β θ(η，θ)空间分布则取决于感应腔初级馈入电流空间分布。
[0092] 2)确定数理方程的待确定常数a」及bj
[0093] 2.1)aj 的确定
[0094] 在感应电压叠加器次级内筒上，在[0,2π]区间内均匀布置多只B-dot探头，测量不 同角向位置Be(rc，0)数值。将不同角向位置的实测Be(r c，0)值进行数学插值拟合，经数据拟 合其满足余弦分布
[0095] B(0)=Bc〇+Bc〇 · δ . c〇s(0) (14)
[0096] 其中，BCC1 = 0.131T，Be(rc，0)角向不均匀系数δ为25%。
[0097]拟合的阴极角向磁场Β( Θ)分布如图2所不。
[0098] 采用正交的三角函数系cos(j0)，将Β(θ)在[0,2jt]区间内作Fourier级数展开，
[0100]由公式（3)确定系数&。= 7.7045;&1 = 0.9631;&」=0(]·彡 2).
[0101] 2.2)bj的确定及电子鞘层边界函数的确定
[0102] 假定IVA次级磁绝缘电子鞘层边界满足以下三种分布形式
[0103] (1)余弦分布
[0104] rsi(9) = -picos(9)+qi (15)
[0105] (2)高斯分布
[0106] rS2(9) =p2exp(-(9-ji)2/w22)+q2 (16)
[0107] (3)二次抛物线分布
[0108] rS3(9) = -p3(9-3r)2+q3 (17)
[0109] 2.3)最优化问题数值方法确定假定分布的待定参数
[0110] 采用最优化问题全局搜索解法，分别确定确定以下三个最优化问题（18-20)的待 定参数口1、91、口2、92、'?2、口3、和93。二个最优化问题的目标函数都是计算误差最小，约束条件 分别是r si(9)s.t.、rsi(9)s.t.和rsi(9)s.t.。
[0117]本实施例中，采用Matlab优化工具箱，确定了口1、91、口2、92、￥2、口3、和93取值分别为 0.01、0.121、0.023、0.112、1.4U0.002和0.126。三种分布形式下的计算误差分别为 2.7084、0.0682和9.6110。高斯分布的计算误差最小，因此，可认为此时磁绝缘电子鞘层边 界近似呈高斯分布形式。确定三种分布形式下电子鞘层边界如图3所示。
[0118]当感应腔初级馈入电流角向均匀性改变，即阴极角向磁场Βθ(η，θ)不均匀系数δ改 变时，系数W(j多1)也会随之改变。由上述方法得到余弦分布、高斯分布和二次抛物线三种 分布形式中待定参数和计算误差如表1所示。
[0119]表1表明，即使阴极角向磁场Be(rc，0)不均匀系数改变，假定的三种分布形式中，高 斯分布的计算误差均最小，因此认为IVA次级磁绝缘电子鞘层边界呈高斯分布。表1可以看 出，当阴极角向磁场Be(rc，0)分布比较均勾时(δ = 1% )，三种分布的计算结果较为接近，鞘 层边界rs(0)~〇.119m，即电子鞘层边界与中心内导体近似同轴心，此数值与一维、对称磁 绝缘Creedon理论计算的电子鞘层边界相同。
[0120]表1阴极角向磁场Be(rc，0)不均匀系数改变时，三种形式的磁绝缘电子鞘层边界及 计算误差。
[0121]
[0122] 采用上述方法确定电子鞘层边界后，IVA次级阴-阳极间隙(rc〈r〈ra)整个区域的磁 绝缘电参数将唯一确定。图4为磁绝缘空间电子的二维分布，在阴极角向磁场Β θ(^，θ)较大 的角度(〇°)，电子鞘层较薄，电子密度ne(r，0)偏大。图5为阴极r。、阳极r a和鞘层边界rs(0) 处角向磁场分布。阴极角向磁场说(^，0)实际上是实验测量值的拟合曲线，满足余弦分布。 采用本发明方法确定电子鞘层边界后，能够获得了电子鞘层边界Be(r s，0)和阳极Be(ra，0)。
[0123] 将图5所示电子鞘层边界rs(0)处角向磁场Be(rs，0)作环路积分，即为阳极电流
[0125]由式(21)计算磁绝缘阳极电流13=136.31^\<^表示真空中的磁导率，μ = 4πχ1〇 -7。
【主权项】
1. 一种感应电压叠加器次级磁绝缘电子銷层边界的确定方法，其特征在于:包括W下 步骤： 1) 建立IVA次级磁绝缘电子銷层边界。（0)满足的数理方程，并明确数理方程的待确定 常数，具体是： 在磁绝缘电子銷层边界。（0)上，磁绝缘电子区域的相对论因子yi(r，0)和真空区域相 对论因子丫。(r，Θ)自身及其法向导数均连续，即其中：丫 S为电子銷层边界。(Θ)上的相对论因子； 丫 i(r，目)和丫。(r，目)表达式如公式(2-a)和(2-b)所示：其中，。为IVA次级阴极半径，。为阳极半径； 丫 〇=l+eV〇/mc2为阳极相对论因子，Vo为次级阴-阳极间隙电压，m为电子质量，m>9.11 Xl〇-3ikg，c 为光速，c = 3Xl〇8m/s; η为正整数，表征阴极表面角向磁场Βθ(。，目）角向分布的模数，n = 0时，Βθ(。，目）角向均 匀分布，即一维轴对称磁绝缘； 〇《j《n，j为整数； aj、bj为待定常数，aj取决于阴极表面角向磁场Be(rc，0)空间分布，Be(rc，0)空间分布则 取决于感应腔初级馈入电流空间分布； 2) 确定数理方程的待确定常数^及6^ 2.1) ?的确定 2.1.1) 测量^4次级阴极角向磁场80(。，0)空间分布，将不同角向位置的实测80(。，0) 数值进行拟合，W获得阴极角向磁场Βθ0-ε，Θ)满足分布Β(θ); 2.1.2) 在[〇,231]区间内对阴极角向磁场80(。，目)满足的分布8(目）作化11'16'级数展开， 得到系数aj;其中，Κ为常数2.2)bj的确定 2.2.1) 根据IVA次级磁绝缘电子銷层边界的物理内涵和可能分布，假定电子銷层边界 满足某种特定的函数分布时(目）； 2.2.2) 离散化数理方程，求取6^与曰^及假定电子銷层边界时(0)的关系； 2.2.3) 最优化问题确定假定分布函数?·〇(θ)中的待定参数；f(白）s.t. {rc《f(白） 2.2.4) 利用假定的电子銷层边界函数分布替代建立的电子銷层边界数理方程求解出 bj; 2.2.5) 求取电子銷层边界数理方程与假定的电子銷层边界函数分布之间的误差； 3)电子銷层边界函数的确定 3.1) 根据IVA次级磁绝缘电子銷层边界的物理内涵和可能分布，假定电子銷层边界满 足其它类型的函数分布时(目）； 3.2) 重复步骤 2.2.2-2.2.5); 3.3) 比较多种假定的电子銷层边界函数与所建立的电子銷层边界数理方程之间的误 差； 3.4) 选择误差最小的假定电子銷层边界函数为最终的电子銷层边界。2.根据权利要求1所述的感应电压叠加器次级磁绝缘电子銷层边界的确定方法，其特 征在于： 所述2.2.2)具体为：其中，ε姻下式3. 根据权利要求1或2所述的感应电压叠加器次级磁绝缘电子銷层边界的确定方法，其 特征在于： 所述步骤2.1.2)是采用正交的Ξ角函数系将由实验数据拟合的分布Β(θ)，在[〇,2π]区 间内作化urier级数展开，得到系数aj的。4. 根据权利要求1或2所述的感应电压叠加器次级磁绝缘电子銷层边界的确定方法，其 特征在于： 步骤(2.2.3)的最优化问题采用全局捜索解法进行数值求解，W确定假定分布函数中 的待定常数，使各假定分布函数本身的计算误差最小； 所述全局捜索解法包括遗传算法、模拟退化方法或蚁群算法。5. 根据权利要求1或2所述的感应电压叠加器次级磁绝缘电子銷层边界的确定方法，其 特征在于： 由于数值求解方法2.2.2和2.2.3是由数理方程lb推导而来，因此计算出的假定分布fo (9)是否满足数理方程la需要进行校核； 另外，步骤2.2.3中求解假定分布函数中的待定常数时，实际上将上述定义计算误差σ 作为最优化问题的目标函数，采取全局捜索计算方法，寻找待定常数，使目标函数σ最小； 为了定量描述运种电子銷层边界确定方法的精度，定义计算误差为计算误差实际上表征在假定电子銷层边界分布函数fo(目）上，磁绝缘电子区域丫 i(f〇 (θ)，θ)和真空区域相对论因子γ〇(?·〇(θ)，θ)之间的误差。
【文档编号】G06F19/00GK106096244SQ201610388829
【公开日】2016年11月9日
【申请日】2016年6月2日 公开号201610388829.2, CN 106096244 A, CN 106096244A, CN 201610388829, CN-A-106096244, CN106096244 A, CN106096244A, CN201610388829, CN201610388829.2
【发明人】魏浩, 孙凤举, 邱爱慈, 呼义翔, 梁天学, 尹佳辉
【申请人】西北核技术研究所