发明名称 一种计算等离子体安全因子分布的方法
摘要 本发明属于核聚变等离子体物理技术领域,具体涉及一种计算等离子体安全因子分布的方法。包括如下步骤:将测量得到的多路密度弦积分数据和弦积分的法拉第旋转角测量数据进行平滑处理;用多路密度弦积分数据进行差值拟合成多于30个空间点的数据,计算出密度分布的数据;根据测量通道与被测量等离子体截面的位置关系,将等离子体电流密度分布j(r)的形式采用峰化因子ν的大环径比圆截面装置的物理近似模型来表示;利用安培定律得到等离子体电流产生的磁场B<sub>p</sub>(r)分布;通过Abel数值反演和理论值与测量值不断趋近方式计算出电流峰化系数ν,并使理论值与测量值误差最小;将峰化因子代入公式得到电流密度的分布和安全因子的分布。本发明提高了测量数据的利用率。
申请公布号 CN104750966A 申请公布日期 2015.07.01
申请号 CN201310754156.4 申请日期 2013.12.31
申请人 核工业西南物理研究院 发明人 焦一鸣
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 核工业专利中心 11007 代理人 高尚梅
主权项 一种计算等离子体安全因子分布的方法,其特征在于:该方法包括如下步骤:步骤1.将测量得到的多路密度弦积分数据和弦积分的法拉第旋转角测量数据进行平滑处理;步骤2.用多路密度弦积分数据进行差值拟合成多于30个空间点的数据,并计算出密度分布的数据;步骤3.根据测量通道与被测量等离子体截面的位置关系,将等离子体电流密度分布j(r)的形式采用峰化因子ν的大环径比圆截面装置的物理近似模型来表示,表达式如下:j(r)=j<sub>0</sub>(1‑r<sup>2</sup>/a<sup>2</sup>)<sup>v</sup>                  (1)其中ν为电流峰化因子;j<sub>0</sub>为等离子体芯部的电流密度值,r为等离子体任意一点的半径,a为等离子体的小半径;利用安培定律得到等离子体电流产生的磁场B<sub>p</sub>(r)分布<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>B</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>&mu;</mi><mn>0</mn></msub><msub><mi>j</mi><mn>0</mn></msub><msup><mi>a</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><mrow><mo>(</mo><mi>v</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>&times;</mo><mfrac><mrow><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>r</mi><mn>2</mn></msup><mo>/</mo><msup><mi>a</mi><mn>2</mn></msup><mo>)</mo></mrow><mrow><mi>v</mi><mo>+</mo><mn>1</mn></mrow></msup></mrow><mi>r</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000451437060000011.GIF" wi="1436" he="160" /></maths>其中μ<sub>0</sub>为真空中的磁常数;将(2)代入偏振仪测量的法拉第旋转角<img file="FDA0000451437060000012.GIF" wi="53" he="62" />的公式<img file="FDA0000451437060000013.GIF" wi="1429" he="175" />其中B<sub>//</sub>=B<sub>p</sub>cosξ,λ,c,ε<sub>0</sub>,m<sub>e</sub>分别为探测束波长,光速,真空介电常数和电子质量,Z<sub>1</sub>和Z<sub>2</sub>为探测束通过等离子体的起始位置,<img file="FDA0000451437060000014.GIF" wi="50" he="57" />为法拉第旋转角测量值,n<sub>e</sub>为密度分布;通过Abel数值反演和理论值与测量值不断趋近方式计算出电流峰化系数v,并使理论值与测量值的误差最小;步骤4.将峰化因子代入公式(1)即得到电流密度的分布,将峰化因子代入公式(2)并利用安全因子q(r)的定义:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>q</mi><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>rB</mi><mi>t</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>R</mi><mn>0</mn></msub><msub><mi>B</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000451437060000021.GIF" wi="1289" he="159" /></maths>式中B<sub>t</sub>(r)为等离子体环向场的分布,R<sub>0</sub>为等离子体的大半径;即得到安全因子的分布。
地址 610041 四川省成都市武侯区二环路南三段三号