Study of magnetic core reset by direct current for induction voltage adder
-
摘要: 根据感应电压叠加器感应腔的工作需求，针对感应腔磁芯的高效可靠应用，一般在感应腔工作前都要对磁芯进行去磁，以实现磁芯的利用最大化，减少磁性材料用量，降低设备体积和造价。通过对磁性材料性能的研究分析，得出了相对于传统的脉冲去磁方式，直流去磁具有电压低、电流小、易绝缘、电极烧蚀小、对变压器油无污染、成本低等优点的结论。根据去磁电流引入位置的不同，研究了两种使用直流对感应腔磁芯进行去磁的方法，并对这两种直流去磁的方法进行了深入地分析和模拟计算，验证了直流去磁的可行性和优越性。Abstract: Because of the working requirement of induction cell and efficient application of magnetic core, the magnetic core must be reset before the induction cell works. Reset can maximize the utilization of magnetic cores, reduce magnetic material consumption, and reduce the volume and cost of equipment. In this paper, the properties of magnetic materials are studied and analyzed. It is found that compared with the pulsed reset method, the direct current reset method has many advantages, such as low voltage, low current, little electrode erosion and low cost. According to the different import position of reset current, two direct current reset methods have been designed. Based on the theoretical analysis and numerical results, the implementation plans have been formulated.
-
Key words:
- IVA /
- induction cell /
- magnetic core /
- magnetic core reset
-
自福岛核电站事故发生以来，反应堆的安全性越来越受到公众的关注，成为影响核能发展的关键因素。采用非能动系统来提高反应堆的固有安全性已经在新的反应堆设计中得到了广泛的应用。其中，利用依靠回路内冷热段内冷却剂的密度差来产生驱动压头带走堆芯余热的自然循环成为主要的研究方向。目前，国际上对自然循环的研究以试验和数值模拟为主。试验研究方面，杨祖毛等[1]对单相自然循环进行了大量的试验，研究了稳态自然循环流量与热源功率、进出口温差的关系；黄彦平[2]通过试验研究了自然循环能力的影响因素等。但试验研究大多针对具体的回路及系统，不具有良好的移植性，且需要大量的资金投入。数值模拟方面，国内外通常使用一些系统分析程序，如RELAP5，RETRAN-02等大型商用软件，对于某些反应堆，也会采用独立开发的专用程序，如田文喜等[3]开发了CARR堆的自然循环能力计算程序，黄洪文等[4]开发了池式研究堆的自然循环能力分析程序，Kaminaga等[5]开发了JRR-3M自然循环能力程序。针对JRR-3M类型的研究堆，稳定自然循环形成前存在冷却剂流动方向反转的过程，有必要开展瞬态安全计算分析。本文利用RELAP5/Mod 3.4对JRR-3M进行了细致的堆芯建模，研究了不同事故条件下自然循环建立前的流量反转过程瞬态热工特性，并得到了瞬态时自然循环的载热极限，同时验证了计算模型和计算方法的可行性。
1. 计算模型
1.1 系统概述
RELAP程序是美国爱达荷国家工程实验室开发、美国核管理委员会批准的用于工程审评的大型瞬态热工水力计算程序，是一维瞬态、两流体六方程流体力学及点堆动力学模型[6]，现已广泛用于压水堆的事故模拟分析。JRR-3M是一座重水慢化、轻水冷却的池式研究堆，在正常运行工况下，冷却剂在两台主泵的驱动下从堆芯上部水池进入堆芯，向下流经堆芯带走裂变热后流入下部的下联箱，与换热器内二次侧的冷却水进行热交换冷却后重新流入主泵。在自然循环阶段，下联箱与水池连接的自然循环阀开启，池水在自然循环压头的驱动下进入下联箱，向上流过堆芯带走衰变热后返回水池，形成自然循环回路。应急冷却系统由两台由不间断电源供电的辅助泵组成，它们的额定流量为75.00 kg·s-1，最多能持续工作3 h。在发生场外电源丧失事故时，一台辅助泵投入运行，带走堆芯热量。
在一回路失去强迫流动后，存在一个由于冷却剂流动方向反转所导致的零流量时刻，此时由于没有冷却剂的流动带走衰变热，堆芯温度迅速上升，该时刻堆芯最有可能被烧毁。基于RELAP5/Mod 3.4程序，对JRR-3M进行建模，反应堆节点图如图 1所示。
其中，考虑到堆水池中存在明显的热分层现象以及内部循环，用两个管道(100，102)来模拟堆芯上部水池，两根管间用横向接管连接，最上方的控制体为空气层，时间控制体130表示空气的温度压力边界，取35 ℃和1.013×105 Pa。104和106是两种不同燃料组件的平均通道，108表示热通道，110和112表示旁路流道，堆芯流道均划分为15个控制体。114表示下联箱，116是环绕堆芯的水池，它们之间用自然循环阀连接。正常运行时从堆芯流出的冷却剂在热交换器118中与202中的冷却水进行冷却热交换。冷却剂由主泵流出后在分支部件126汇合后流回堆芯上部水池。
对于堆芯内功率的分布取该池式反应堆的物理设计结果，其功率衰变曲线采用中国行业标准EJ/T 745-2001[7]推荐公式计算给出, 即
P(τ0)=G(t)4∑i=1∫τ00Pi(t)Qif(t)dt (1) Q(t)=1.0+τ0.40φ(3.24×10−6+5.23×10−10t) (2) 式中：P(τ0)是τ0时刻的裂变总功率；G(t)是考虑裂变产物中子俘获的校正因子；Pi(t)是核素i在t时刻的裂变功率；Qi为核素i一次裂变可回收能量(MeV)；f(t)为可裂变核素的衰变热功率函数，它的值由标准附表给出(略)；φ为一个初始裂变原子的裂变数；i取1，2，3，4分别代表235U热中子裂变、239Pu热中子裂变、238U快中子裂变和241Pu热中子裂变对衰变热功率的贡献；t为停堆后时间(s)。
1.2 稳态参数
将系统调节至稳态工况是进行瞬态计算的前提条件。在自然循环的瞬态计算中，系统各个参数应在某个预计工况下稳态运行，随后关闭主泵打开自然循环阀进行瞬态计算，而这个稳态的工况就决定了事故的初始状态。通过调节回路内的物理参数，使得稳态工况下的热工参数与设计的热工参数一致且达到稳定。表 1为该池式研究堆的热工设计参数。
表 1 池式研究堆主要设计参数Table 1. Main design parameters of pool-type research reactortotal power/MW total mass flow rate/(kg·s-1) inlet temperature/℃ outlet temperature/℃ core inlet pressure/MPa 20.00 661.91 35.00 42.19 1.69 1.3 安全准则
按照该堆设计确定的安全标准为：
(1) 使用Sudo临界热流密度经验公式[8]，计算的最小偏离泡核沸腾比(MDNBR)应大于限值1.5。
(2) 燃料元件芯体最高温度低于250 ℃。
2. 计算结果与分析
2.1 无应急冷却的场外电源丧失事故分析
反应堆在满功率正常运行了100 s后发生外电源丧失事故，主泵由于失去电力供应立即失电惰转，二次侧热阱立即丧失，自然循环阀打开，1 s后控制棒落下反应堆进行停堆，此时堆芯所产生的热量主要依靠主泵惰转和自然循环带走。整个过程的事件序列列于表 2。
表 2 无应急冷却系统事件序列Table 2. Accident process without emergency cooling systemevent start time event start time main pump loss of off-site power supply 100 s control rods drop 101 s secondary side loss of hot trap 100 s finish main pump coasting 220 s open the natural circulation valve 100 s reverse coolant flow 221 s 图 2为事故分析的结果曲线。
从图 2(a)和图 2(b)中可以看出，在0~100 s内，反应堆满功率稳态运行，堆芯流量保持不变。在100 s时刻，由于主泵失电惰转流量迅速减少，在t=220 s时主泵惰转停止，堆芯流量减少至0。此时自然循环阀已经打开，在自然循环驱动压头的作用下，由于热通道中热流密度最大，故最先产生流量反转，随后堆芯流量反转，冷却剂由下而上地通过堆芯，自然循环迅速建立。
从图 2(c)看出，在稳态运行时，燃料芯块最热温度为108.27 ℃，包壳最高温度为99.37 ℃, 在101 s时控制棒开始下落后，由于堆芯功率的呈指数式降低，燃料板和冷却剂的温度也随之急剧下降。随后，由于冷却剂流量不足以完全带走燃料所产生的热量，热构件和冷却剂的温度逐渐上升。在220 s时冷却剂流量为0，燃料外表面与冷却剂的对流换热性能最差，燃料板的温度达到峰值，燃料芯块的峰值温度为114.30 ℃, 包壳峰值温度为114.07 ℃。此后，由于自然循环流量的增大，对流换热作用增强，燃料板和冷却剂的温度开始下降。在稳定的自然循环建立后，燃料板和冷却剂的温度由于堆芯功率的缓慢衰减而缓慢降低。
图 2(d)是整个过程中的MDNBR变化曲线。从图中可以看出，热管处的MDNBR在稳态运行时为2.66，在反应堆停堆后由于燃料板热流密度的急剧减少而迅速增加，随后随着流量的降低而降低，在流量反转时刻为最小值，最小值为0.31，远低于安全限值1.50。
2.2 有应急冷却的场外电源丧失事故分析
当该池式研究堆在满功率正常运行了100 s后发生外电源丧失事故，主泵失电惰转，当冷却剂流量低于75.00 kg·s-1时，一台辅助泵投入运行。当t=3 188 s时，衰变功率已降至200 kW，辅助泵停止运行，在辅助泵停止运行前5 s打开自然循环阀。事故过程中的事件时间序列如表 3所示。
表 3 有应急冷却事件序列Table 3. Accident process with emergency cooling systemevent start time event start time main pump loss of off-site power supply 100 s close auxiliary pump 3 188 s control rods drop 101 s auxiliary pump finish coast 3 198 s input auxiliary pump 110 s coolant flow reverse 3 199 s open the natural circulation valve 3 183 s 图 3为事故分析的结果曲线。
从图 3(a)和图 3(b)中可以看出，在100 s以前，反应堆满功率稳态运行，堆芯流量不变。100 s后主泵失电惰转，流量降低。当流量降低至75.00 kg·s-1时，辅助泵投入运行，堆芯流量保持不变。当t=3 183 s时，自然循环阀打开，从图中可以发现堆芯流量平台有一个下降。这时由于此时堆水池与反应堆一回路连接，具有了一定的流量，导致流过堆芯的流量减少。在t=3 188 s时关闭辅助泵，堆芯流量逐渐降低至0，随后由于自然循环的建立，流量发生反转，最后趋向稳定。
通过比较图 2(c)和图 3(c)可以发现，在投入应急冷却系统的事故工况下，由于堆芯拥有一定程度的强迫流量，这一部分的流量带走了部分衰变热，有效地降低了燃料板及冷却剂的温度，相较于未投入应急冷却系统的条件下，燃料芯块峰值温度降低了46.10 ℃，堆芯材料不会发生熔毁。同时可以发现冷却剂的温度和燃料板温度的变化趋势相似，但有一定的滞后。
从图 3(d)可以看出，有应急冷却时事故过程中的MDNBR始终大于安全限值，不会发生偏离核态沸腾，池式研究堆始终处于安全状态。
2.3 最大载热能力计算
为进一步探究池式研究堆自然循环瞬态过程中的最大载热能力，在堆芯功率衰减至不同功率水平时关闭辅助泵并打开自然循环阀，不同功率下自然循环过程中的最高燃料温度和MDNBR如图 4所示。
从图 4可以明显看出，此时反应堆的安全性主要受到来自MDNBR的限制，因此，池式研究堆最大瞬态自然循环能力约为590 kW。
3. 结论
本文利用系统性分析程序RELAP5/Mod 3.4对JRR-3M池式研究堆进行了建模，计算了在无应急冷却和有应急冷却情况下的丧失场外电源事故。RELAP5很好地模拟了堆芯流量的反转过程以及自然循环的建立过程。针对这两种事故工况，计算得到了堆芯各个流道的温度、压力、冷却剂流量等热工水力参数。计算结果表明：若未投入应急冷却系统，反应堆在失去场外电源后热管内MDNBR为0.31，小于1.5的安全限值，会发生偏离核态沸腾，存在烧毁堆芯的风险。应急冷却系统能有效地降低事故后燃料板及冷却剂的温度，提高了系统的安全性。在瞬态条件下，池式研究堆的最大自然循环能力为590 kW，大于设计值200 kW，表明该堆具有良好的固有安全性。
-
图 9 感应腔模拟结构及去磁磁场分布图
Figure 9. Simulation structure of induction cell and distribution of demagnetized magnetic field
图 10 不同位置磁场强度的分布曲线
Figure 10. Magnetic field intensity curves for the three lines in fig. 9
表 1 不同励磁脉冲宽度下铁基非晶磁芯的Hc
Table 1. Coercive force (Hc) under different pulse width of current
pulse width of current/μs Hc/(A·m－1) DC 4 20 000 10 120 220 27 300 3 500 
下载: 导出CSV
-
[1] 刘锡三. 高功率脉冲技术[M]. 北京: 国防工业出版社. 2007: 468-472.Liu Xisan. High pulsed power technology. Beijing: National Defense Industry Press. 2007: 468-472 [2] 丛培天, 张国伟, 吴撼宇, 等. 3 MV感应电压叠加器的磁感应腔研制[J]. 强激光与粒子束, 2011, 23(2): 563-568. http://www.hplpb.com.cn/article/id/5012Cong Peitian, Zhang Guowei, Wu Hanyu, et al. Design of inductive cavity for 3MV induction voltage adder. High Power Laser and Particle Beams, 2011, 23(2): 563-568 http://www.hplpb.com.cn/article/id/5012 [3] Smith I D. Induction voltage adders and the induction accelerator family[J]. Physical Review Special Topics—Accelerators and Beams, 2004, 7: 064801. doi: 10.1103/PhysRevSTAB.7.064801 [4] Ramirez J J, Prestwich K R, Burgess E L, et al. The Hermes-Ⅲ program[C]//Proceedings of 6th International Pulsed Power Conference. 1987: 294-299. [5] Johnson D L, Ramirez J J, Huddle C W, et al. Hermes-Ⅲ prototype cavity tests[C]//Proceedings of 6th International Pulsed Power Conference. 1987: 482-485. [6] Weidenheimer D, Corcoran P, Altes R, et al. Design of a driver for the Cygnus X-ray source[C]//Proceedings of 13th International Pulsed Power Conference. 2001: 591-595. [7] 张国伟, 丛培天, 乔开来, 等. 脉冲磁化条件下非晶磁芯的损耗特性[J]. 强激光与粒子束, 2012, 24(5): 1247-1250. doi: 10.3788/HPLPB20122405.1247Zhang Guowei, Cong Peitian, Qiao Kailai, et al. Loss characteristics of glassy alloy magnetic core under pulsed magnetization. High Power Laser and Particle Beams, 2012, 24(5): 1247-1250 doi: 10.3788/HPLPB20122405.1247 [8] 王志国, 孙凤举, 邱爱慈, 等. 快放电直线型变压器驱动源磁芯脉冲损耗特性[J]. 强激光与粒子束, 2012, 24(3): 524-528. http://www.hplpb.com.cn/article/id/6768Wang Zhiguo, Sun Fengju, Qiu Aici, et al. Loss characteristics of FLTD magnetic cores under fast pulsed voltage. High Power Laser and Particle Beams, 2012, 24(3): 524-528 http://www.hplpb.com.cn/article/id/6768 [9] 张瑜, 王仕文, 刘金亮, 等. 铁基非晶磁环的高频响应和ns级脉冲响应特性[J]. 强激光与粒子束, 2011, 23(7): 1979-1984. http://www.hplpb.com.cn/article/id/5323Zhang Yu, Wang Shiwen, Liu Jinliang, et al. High frequency response and ns range pulse response of iron based amorphous cores. High Power Laser and Particle Beams, 2011, 23(7): 1979-1984 http://www.hplpb.com.cn/article/id/5323 [10] Barrett D M. Core reset considerations in magnetic pulse compression network[C]//Proceedings of 10th International Pulsed Power Conference. 1995: 1160-1165. [11] 许军, 李坤. 电接触的接触电阻研究[J]. 电工材料, 2011(1): 10-13. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHJ201101003.htmXu jun, Li kun. The research on resistance of electrical contact. Electrical Engineering Materials, 2011(1): 10-13 https://www.cnki.com.cn/Article/CJFDTOTAL-JGHJ201101003.htm [12] 吴楠, 庄劲武, 江壮贤, 等. 螺栓紧固铜排间接触电阻研究[J]. 海军工程大学学报, 2016, 28(6): 28-31. https://www.cnki.com.cn/Article/CJFDTOTAL-HJGX201606006.htmWu Nan, Zhuang Jinwu, Jiang Zhuangxian, et al. Contact resistance between bolt fastening copper platoons. Journal of Naval University of Engineering, 2016, 28(6): 28-30 https://www.cnki.com.cn/Article/CJFDTOTAL-HJGX201606006.htm 期刊类型引用(3)
1. 何进，赵士龙，郭凌伟，华有杰，叶仁广，徐时清. 反射式体布拉格光栅设计及其在激光器中的应用研究. 光电子·激光. 2024(02): 122-127 . 
百度学术2. 李奎奎，张伟桥，董风鑫，张建心，周旭彦. 高功率半导体激光叠阵光纤耦合系统与光束整形设计. 中国激光. 2024(17): 98-103 . 百度学术3. 任慧娟. 激光惯性约束下多帧影视图像视觉传达模型构建. 激光杂志. 2023(09): 103-108 . 百度学术其他类型引用(1)
-
下载:
点击查看大图
计量
- 文章访问数: 1411
- HTML全文浏览量: 334
- PDF下载量: 124
- 被引次数: 4
