Analysis of high-altitude electromagnetic pulse effect on wireless communication network from hierarchical perspective
-
摘要: 从分层视角研究无线通信网络高空电磁脉冲效应，提出了一种具有层级属性的网络级效应评估方法。设计研制了一种无线通信网络效应试验仿真平台，网络使用典型全连通拓扑结构，网络节点使用超外差式收发机，网络运行在超短波频段。开展了小型网络脉冲辐照效应试验和效应机理分析，网络仿真研究了网络规模不同条件下节点失效概率对网络性能的影响。得到以下结论：（1）设备级效应对低层级效应具备一定冗余容错能力,这取决于效应评价指标的选择和损伤组件的功能属性；（2）网络级效应对设备级效应具有潜在放大效果，个别节点失效可能造成网络性能严重下降甚至瘫痪。Abstract: To study the representation of high altitude electromagnetic pulse effect in communication network from a hierarchical perspective, a network-level effect assessment method with hierarchical feature was proposed. A network effect test and simulation platform for wireless network was designed and developed. The platform employs a typically fully-connected network, in which the communication node is designed based on the framework of superheterodyne transceivers with operating on very high frequency ranges. The effect of pulse irradiation test on this network platform was tested and its effect mechanism was analyzed. The influence of node failure probability on network performance under different network scales was studied by using OPNET software. The following conclusions are drawn: (1) The equipment-level effect has a certain redundant fault-tolerance for the lower-level damage effect, which depends on the selection of effect assessment indicators and the function properties of the damaged components; (2) The network-level effect has a potential effect on the equipment-level effect, and the failure of individual nodes may cause serious degradation or even paralysis of the whole network.
-
自福岛核电站事故发生以来，反应堆的安全性越来越受到公众的关注，成为影响核能发展的关键因素。采用非能动系统来提高反应堆的固有安全性已经在新的反应堆设计中得到了广泛的应用。其中，利用依靠回路内冷热段内冷却剂的密度差来产生驱动压头带走堆芯余热的自然循环成为主要的研究方向。目前，国际上对自然循环的研究以试验和数值模拟为主。试验研究方面，杨祖毛等[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，表明该堆具有良好的固有安全性。
-
表 1 基于层级关系的效应现象列表
Table 1. List of effect phenomena based on hierarchy
layer effect（indicator） effect type network throughput declined rapidly performance degradation equipment receive sensitivity declined by 30～50 dB performance degradation component LNA grain declined by 30～50 dB damage device the internal triodes of LNA have different degrees of breakdown damage 
下载: 导出CSV
-
[1] 毛从光, 程引会, 谢彦召. 高空电磁脉冲技术基础[M]. 北京: 科学出版社, 2019Mao Congguang, Cheng Yinhui, Xie Yanzhao. Fundamentals of high altitude electromagnetic pulse technology[M]. Beijing: Science Press, 2019 [2] 周璧华, 陈彬, 石立华. 电磁脉冲及其工程防护[M]. 北京: 国防工业出版社, 2003Zhou Bihua, Chen Bin, Shi Lihua. EMP and EMP protection[M]. Beijing: National Defense Industry Press, 2003 [3] Giri D V, Hoad R, Sabath F. High-power electromagnetic effects on electronic systems[M]. Norwood: Artech House, 2020: 20-30. [4] Foster J S Jr, Gjelde E, Graham W R, et al. Report of the commission to assess the threat to the United States from electromagnetic pulse (EMP) attack: critical national infrastructures[R]. EMP Commission USA, 2008. [5] GJB538-88, 半导体器件电磁脉冲损伤阈值试验方法[S]GJB538-88, 半导体器件电磁脉冲损伤阈值试验方法[S]. (GJB538-88, Test method for electromagnetic pulse damage threshold of semiconductor devices[S] [6] 汪海洋, 周翼鸿, 李家胤, 等. 低噪声放大器有意电磁干扰效应(英文)[J]. 强激光与粒子束, 2011, 23(11):2865-2871. (Wang Haiyang, Zhou Yihong, Li Jiayin, et al. LNA malfunctions under intentional electromagnetic interference[J]. High Power Laser and Particle Beams, 2011, 23(11): 2865-2871 doi: 10.3788/HPLPB20112311.2865 [7] 任兴荣. 半导体器件的电磁损伤效应与机理研究[D]. 西安: 西安电子科技大学, 2014Ren Xingrong. Research on the electromagnetic damage effects and mechanisms of semiconductor devices[D]. Xi’an: Xidian University, 2014 [8] 王洋. 接收机前端低噪声放大器设计及其电磁脉冲防护研究[D]. 上海: 上海交通大学, 2011Wang Yang. Design and electromagnetic pulse protection research of receiver front-end LNA[D]. Shanghai: Shanghai Jiao Tong University, 2011 [9] Camp M, Gerth H, Garbe H, et al. Predicting the breakdown behavior of microcontrollers under EMP/UWB impact using a statistical analysis[J]. IEEE Transactions on Electromagnetic Compatibility, 2004, 46(3): 368-379. doi: 10.1109/TEMC.2004.831816 [10] Nitsch D, Camp M, Sabath F, et al. Susceptibility of some electronic equipment to HPEM threats[J]. IEEE Transactions on Electromagnetic Compatibility, 2004, 46(3): 380-389. doi: 10.1109/TEMC.2004.831842 [11] Coburn W O, Nguyen E, Reyzer R J, et al. High-altitude electromagnetic pulse survivability assessment of the Harris RF-3200 transceiver[R]. AD-A258 347, Adelphi, MD: Harry Diamond Laboratories, 1992. [12] Hoad R, Carter N J, Herke D, et al. Trends in EM susceptibility of IT equipment[J]. IEEE Transactions on Electromagnetic Compatibility, 2004, 46(3): 390-395. doi: 10.1109/TEMC.2004.831815 [13] Hoad R, Carter N J, Herke D, et al. An investigation into radiated susceptibility of IT networks[C]//Proceedings of the International Symposium on Electromagnetic Compatibility. Eindhoven, The Netherland: QinetiQ, Electromagnetic Compatibility Group, 2004. [14] IEC TR 61000-1-3, Electromagnetic compatibility (EMC) –part 1-3: general – the effects of high-altitude EMP (HEMP) on civil equipment and systems[S]. [15] IEC 61000-5-9, Electromagnetic compatibility (EMC) - part 5-9: installation and mitigation guidelines – system-level susceptibility assessments for HEMP and HPEM[S]. [16] Frank H, Frisch I. Analysis and design of survivable networks[J]. IEEE Transactions on Communication Technology, 1970, 18(5): 501-519. doi: 10.1109/TCOM.1970.1090419 [17] 魏福林. 野战地域通信网拓扑层抗毁性研究[D]. 郑州: 解放军信息工程大学, 2006Wei Fulin. Research on invulnerability based on topological structure of regional communication networks[D]. Zhengzhou: Information Engineering University, 2006 [18] Mao Congguang, Canavero F. System-level vulnerability assessment for EME: from fault tree analysis to Bayesian networks—part I: methodology framework[J]. IEEE Transactions on Electromagnetic Compatibility, 2016, 58(1): 180-187. doi: 10.1109/TEMC.2015.2484067 [19] Mao Congguang, Canavero F G, Cui Zhitong, et al. System-level vulnerability assessment for EME: from fault tree analysis to Bayesian networks-part II: illustration to microcontroller system[J]. IEEE Transactions on Electromagnetic Compatibility, 2016, 58(1): 188-196. doi: 10.1109/TEMC.2015.2502591 期刊类型引用(3)
1. 何进，赵士龙，郭凌伟，华有杰，叶仁广，徐时清. 反射式体布拉格光栅设计及其在激光器中的应用研究. 光电子·激光. 2024(02): 122-127 . 
百度学术2. 李奎奎，张伟桥，董风鑫，张建心，周旭彦. 高功率半导体激光叠阵光纤耦合系统与光束整形设计. 中国激光. 2024(17): 98-103 . 百度学术3. 任慧娟. 激光惯性约束下多帧影视图像视觉传达模型构建. 激光杂志. 2023(09): 103-108 . 百度学术其他类型引用(1)
-
下载:
点击查看大图
计量
- 文章访问数: 1171
- HTML全文浏览量: 511
- PDF下载量: 49
- 被引次数: 4
