发明名称 一种感应电压叠加器次级磁绝缘电子鞘层边界的确定方法
摘要 本发明涉及一种感应电压叠加器次级磁绝缘电子鞘层边界的确定方法包括,(1)建立磁绝缘电子鞘层边界满足的数学物理方程,推导电子鞘层边界内/外侧相对论因子的通用解析表达式。(2)由实验测量阴极角向磁场空间分布,结合数值求解方法,获得通解表达式中的待定系数。(3)假定电子鞘层边界满足几种给定的分布形式(如余弦分布、二次抛物线分布、高斯分布等),确定各假定分布中的待定参数和计算误差,误差最小的分布视为电子鞘层边界的真实分布。相对于现有粒子模拟确定方法,本方法回避了目前国内三维粒子模拟软件及代码不成熟的技术问题,具有准确性高、耗时短的优点。
申请公布号 CN106096244A 申请公布日期 2016.11.09
申请号 CN201610388829.2 申请日期 2016.06.02
申请人 西北核技术研究所 发明人 魏浩;孙凤举;邱爱慈;呼义翔;梁天学;尹佳辉
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 西安智邦专利商标代理有限公司 61211 代理人 杨亚婷
主权项 一种感应电压叠加器次级磁绝缘电子鞘层边界的确定方法,其特征在于:包括以下步骤:1)建立IVA次级磁绝缘电子鞘层边界r<sub>s</sub>(θ)满足的数理方程,并明确数理方程的待确定常数,具体是:在磁绝缘电子鞘层边界r<sub>s</sub>(θ)上,磁绝缘电子区域的相对论因子γ<sub>i</sub>(r,θ)和真空区域相对论因子γ<sub>o</sub>(r,θ)自身及其法向导数均连续,即<maths num="0001"><math><![CDATA[<mrow><msub><mi>&gamma;</mi><mi>s</mi></msub><mo>=</mo><msub><mi>&gamma;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&gamma;</mi><mi>o</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>;</mo><mi>&theta;</mi><mo>)</mo></mrow><msub><mo>|</mo><mrow><mi>r</mi><mo>=</mo><msub><mi>r</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>a</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001007990020000014.GIF" wi="1245" he="104" /></maths><maths num="0002"><math><![CDATA[<mrow><mfrac><mrow><mo>&part;</mo><msub><mi>&gamma;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow><mrow><mo>&part;</mo><mi>r</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mo>&part;</mo><msub><mi>&gamma;</mi><mi>o</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow><mrow><mo>&part;</mo><mi>r</mi></mrow></mfrac><msub><mo>|</mo><mrow><mi>r</mi><mo>=</mo><msub><mi>r</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>b</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001007990020000011.GIF" wi="1006" he="143" /></maths>其中:γ<sub>s</sub>为电子鞘层边界r<sub>s</sub>(θ)上的相对论因子;γ<sub>i</sub>(r,θ)和γ<sub>o</sub>(r,θ)表达式如公式(2‑a)和(2‑b)所示:<maths num="0003"><math><![CDATA[<mrow><msub><mi>&gamma;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><mi>cosh</mi><mo>&lsqb;</mo><msub><mi>a</mi><mn>0</mn></msub><mi>l</mi><mi>n</mi><mrow><mo>(</mo><mfrac><mi>r</mi><msub><mi>r</mi><mi>c</mi></msub></mfrac><mo>)</mo></mrow><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><msub><mi>a</mi><mi>j</mi></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><msup><mrow><mo>(</mo><mfrac><mi>r</mi><msub><mi>r</mi><mi>c</mi></msub></mfrac><mo>)</mo></mrow><mi>j</mi></msup><mo>-</mo><msup><mrow><mo>(</mo><mfrac><msub><mi>r</mi><mi>c</mi></msub><mi>r</mi></mfrac><mo>)</mo></mrow><mi>j</mi></msup><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><mi>j</mi><mi>&theta;</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>-</mo><mi>a</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001007990020000012.GIF" wi="1558" he="175" /></maths><maths num="0004"><math><![CDATA[<mrow><msub><mi>&gamma;</mi><mi>o</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&gamma;</mi><mn>0</mn></msub><mo>+</mo><msub><mi>b</mi><mn>0</mn></msub><mi>l</mi><mi>n</mi><mrow><mo>(</mo><mfrac><mi>r</mi><msub><mi>r</mi><mi>a</mi></msub></mfrac><mo>)</mo></mrow><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><msub><mi>b</mi><mi>j</mi></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><msup><mrow><mo>(</mo><mfrac><mi>r</mi><msub><mi>r</mi><mi>a</mi></msub></mfrac><mo>)</mo></mrow><mi>j</mi></msup><mo>-</mo><msup><mrow><mo>(</mo><mfrac><msub><mi>r</mi><mi>a</mi></msub><mi>r</mi></mfrac><mo>)</mo></mrow><mi>j</mi></msup><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><mi>j</mi><mi>&theta;</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>-</mo><mi>b</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001007990020000013.GIF" wi="1566" he="159" /></maths>其中,r<sub>c</sub>为IVA次级阴极半径,r<sub>a</sub>为阳极半径;γ<sub>0</sub>=1+eV<sub>0</sub>/mc<sup>2</sup>为阳极相对论因子,V<sub>0</sub>为次级阴‑阳极间隙电压,m为电子质量,m≈9.11×10<sup>‑31</sup>kg,c为光速,c=3×10<sup>8</sup>m/s;n为正整数,表征阴极表面角向磁场B<sub>θ</sub>(r<sub>c</sub>,θ)角向分布的模数,n=0时,B<sub>θ</sub>(r<sub>c</sub>,θ)角向均匀分布,即一维轴对称磁绝缘;0≤j≤n,j为整数;a<sub>j</sub>、b<sub>j</sub>为待定常数,a<sub>j</sub>取决于阴极表面角向磁场B<sub>θ</sub>(r<sub>c</sub>,θ)空间分布,B<sub>θ</sub>(r<sub>c</sub>,θ)空间分布则取决于感应腔初级馈入电流空间分布;2)确定数理方程的待确定常数a<sub>j</sub>及b<sub>j</sub>2.1)a<sub>j</sub>的确定2.1.1)测量IVA次级阴极角向磁场B<sub>θ</sub>(r<sub>c</sub>,θ)空间分布,将不同角向位置的实测B<sub>θ</sub>(r<sub>c</sub>,θ)数值进行拟合,以获得阴极角向磁场B<sub>θ</sub>(r<sub>c</sub>,θ)满足分布B(θ);2.1.2)在[0,2π]区间内对阴极角向磁场B<sub>θ</sub>(r<sub>c</sub>,θ)满足的分布B(θ)作Fourier级数展开,得到系数a<sub>j</sub>;<maths num="0005"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>a</mi><mn>0</mn></msub><mo>=</mo><mfrac><mrow><mi>&kappa;</mi><mo>&CenterDot;</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mrow><mn>2</mn><mi>&pi;</mi></mrow></msubsup><mi>B</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>d</mi><mi>&theta;</mi></mrow><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>a</mi><mi>j</mi></msub><mo>=</mo><mfrac><mrow><mi>&kappa;</mi><mo>&CenterDot;</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mrow><mn>2</mn><mi>&pi;</mi></mrow></msubsup><mi>B</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>cos</mi><mrow><mo>(</mo><mrow><mi>j</mi><mi>&theta;</mi></mrow><mo>)</mo></mrow><mi>d</mi><mi>&theta;</mi></mrow><mrow><mn>2</mn><mi>n</mi><mi>&pi;</mi></mrow></mfrac></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>&le;</mo><mi>j</mi><mo>&le;</mo><mi>n</mi></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001007990020000021.GIF" wi="1302" he="443" /></maths>其中,κ为常数,<img file="FDA0001007990020000022.GIF" wi="262" he="150" />2.2)b<sub>j</sub>的确定2.2.1)根据IVA次级磁绝缘电子鞘层边界的物理内涵和可能分布,假定电子鞘层边界满足某种特定的函数分布f<sub>0</sub>(θ);2.2.2)离散化数理方程,求取b<sub>j</sub>与a<sub>j</sub>及假定电子鞘层边界f<sub>0</sub>(θ)的关系;2.2.3)最优化问题确定假定分布函数f<sub>0</sub>(θ)中的待定参数;<maths num="0006"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow></mtd><mtd><mrow><mi>&sigma;</mi><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mrow><mn>2</mn><mi>&pi;</mi></mrow></msubsup><mo>&lsqb;</mo><msub><mi>&gamma;</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mn>0</mn></msub><mo>(</mo><mi>&theta;</mi><mo>)</mo><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&gamma;</mi><mi>o</mi></msub><mo>(</mo><mrow><msub><mi>f</mi><mn>0</mn></msub><mo>,</mo><mi>&theta;</mi></mrow><mo>)</mo><mo>)</mo><msup><mo>&rsqb;</mo><mn>2</mn></msup><mi>d</mi><mi>&theta;</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001007990020000023.GIF" wi="1475" he="151" /></maths>f(θ)s.t.{r<sub>c</sub>≤f(θ)≤r<sub>a</sub>2.2.4)利用假定的电子鞘层边界函数分布替代建立的电子鞘层边界数理方程求解出b<sub>j</sub>;2.2.5)求取电子鞘层边界数理方程与假定的电子鞘层边界函数分布之间的误差;3)电子鞘层边界函数的确定3.1)根据IVA次级磁绝缘电子鞘层边界的物理内涵和可能分布,假定电子鞘层边界满足其它类型的函数分布f<sub>0</sub>(θ);3.2)重复步骤2.2.2‑2.2.5);3.3)比较多种假定的电子鞘层边界函数与所建立的电子鞘层边界数理方程之间的误差;3.4)选择误差最小的假定电子鞘层边界函数为最终的电子鞘层边界。
地址 710024 陕西省西安市69信箱