RPV侧边破口条件下堆外蒸汽爆炸计算分析

针对实际过程中更有可能发生的压力容器（RPV）侧边破口条件开展蒸汽爆炸计算分析。根据经济合作与发展组织（OECD）发布的现象识别与重要度排序表（PIRT）,选取堆外蒸汽爆炸敏感性分析参数,使用MC3D软件建立三维局部破口和二维环状破口几何模型,对影响计算结果的重要参数（破口尺寸、堆坑水位、破口位置、触发条件、液柱碎化和液滴碎化模型）开展RPV侧边破口条件下敏感性分析,获得最恶劣计算工况条件。敏感性分析结果表明,在大破口失水事故（LBLOCA）工况下,当堆坑处于满水位、RPV发生二维侧边环状破口、接触堆坑侧壁面时触发蒸汽爆炸、采用CONST模型和Classical模型时,堆坑侧壁面的压力载荷计算结果最为保守,对堆坑和安全壳完整性威胁最大。      

最佳估算方法中不确定度评估关键问题分析

应用最佳估算+不确定度（BEPU）分析方法对核电厂进行事故分析或安全评审已成为国际发展趋势。本文对最佳估算分析中基于输入传递的统计类不确定度评估的流程进行了总结,并对其关键步骤进行了分析和研究。分析认为,评估流程可分为确定目标参数、确定重要输入参数及其分布、抽样、模型分析和目标参数分析5步,其中现象识别和重要度排序表（PIRT）是一种适用的重要输入参数确定方法,输入参数的分布需根据试验数据或专家判断确定;抽样方法上,可采用参数抽样或非参数抽样,后者可大幅减小抽样数量;不确定度评估所用模型须经过充分试验或分析证明其适用性;通过对目标参数进行统计,可获得不确定度范围及输入参数的敏感性。     

GOTHIC程序在海洋核动力平台反应堆舱安全分析领域的适用性评估

针对海洋核动力平台反应堆舱热工水力分析程序缺乏的现状，以一回路失水事故（LOCA）下反应堆舱压力响应为评价基准，基于安全壳现象识别与排序表（PIRT）分析方法，通过开展LOCA下反应堆舱热工水力现象识别、现象分级研究，建立了反应堆舱PIRT。通过开展GOTHIC程序模型验证矩阵与PIRT的匹配性分析，确认GOTHIC程序在海洋核动力平台反应堆舱热工水力分析领域的适用性。本文分析方法对其他安全分析程序在核电等领域的跨领域适用性评估具有一定参考价值。      

适用于核电厂BEPU的高效全局敏感性分析方法开发及应用

最佳估算加不确定性（BEPU）分析是IAEA推荐用于核电厂事故安全分析的方法,该方法中一个关键步骤为评估输入参数对目标输出的影响大小,即定量敏感性分析。传统BEPU分析中常使用基于线性或单调假设的局部敏感性分析方法,其难以适用于复杂的核反应堆系统,而全局敏感性分析则由于计算成本过高而难以在实际工程中应用。本研究中针对矩独立全局敏感性分析方法开展了优化研究,使用高阶模型表示、高斯求积公式等方法降低矩独立敏感性度量的计算成本,得到了一种高效的敏感性分析方法。使用了多个例题对优化方法的可靠性进行了验证,并将其应用于LOFT（loss of fluid test）大破口事故的敏感性分析。结果表明,该高效敏感性分析方法能准确识别核反应堆事故工况中的重要参数,并能对参数重要度进行定量排序。     

大型先进压水堆弹棒事故模拟及敏感性分析

相对于传统堆型,大型非能动先进压水堆堆芯设计具有重大改变,这些改变对弹棒事故分析具有重要影响,进而影响反应堆的安全性。通过选取典型的四类工况（寿期初满功率、寿期初零功率、寿期末满功率和寿期末零功率）,利用中子动力学软件和燃料性能分析程序开展大型先进压水堆CAP1400的弹棒事故模拟计算,验证大型先进压水堆弹棒事故工况下的安全性,并针对弹棒事故分析关键输入参数开展敏感性分析。计算分析结果表明：大型先进压水堆发生弹棒事故时,其结果能够满足验收准则的要求,反应堆处于安全可控状态;弹棒事故分析中功率峰值对弹棒价值最敏感,事故分析结果对停堆反应性敏感性较小。     

安全壳卸压排放过程模化分析

针对大型非能动先进压水堆安全壳卸压排放过程中涉及的重要热工现象，采用系统性的关键现象识别及重要性分析方法，得到了大型非能动先进压水堆卸压排放过程中的现象过程识别与排序表（PIRT）。结果表明：排放管线及鼓泡器中对安全壳卸压排放过程影响程度较高的现象为临界和摩擦流、两相压降、几何尺寸及流动状态；乏燃料水池中对安全壳卸压排放过程影响程度较高的现象为冷凝、传热、几何尺寸、流体混合、不凝性气体及热分层。利用关键现象识别及重要性分析结果与现有缩放实验台架的搭建经验及研究结果，得到了安全壳卸压排放过程验证性试验装置搭建中应该遵循的相似准则，从而为安全壳卸压排放验证性试验装置的搭建提供设计基础和理论依据。     

基于HFR的核动力装置非能动系统敏感性分析方法研究

敏感性分析应用于反应堆非能动系统热工水力过程的不确定性分析和可靠性分析，能够定量识别对系统热工水力行为具有重要影响的不确定性输入参数。基于混合随机均衡-傅里叶幅度敏感性测试（HFR）方法，以某型核动力装置非能动余热排出试验系统作为算例进行全局敏感性分析研究，仿真结果证明了HFR方法的可行性与正确性。敏感性分析给出了系统输入参数重要度随时间的变化规律以及系统稳定运行时输入参数的重要度排序，分析结果有助于指导系统的设计优化及运行管理。      

应用DAKOTA程序耦合WGOTHIC程序进行安全壳压力响应敏感性分析

安全壳压力响应分析是验证非能动安全壳冷却系统(PCS)设计的重要内容,需考虑PCS的传热传质等各种现象的影响。本文应用DAKOTA程序耦合WGOTHIC程序对大型先进压水堆非能动安全壳压力响应进行敏感性分析,通过偏相关系数,定量评价了重要现象识别和排序表(PIRT)中各种现象对安全壳压力的影响程度。研究结果表明:质能释放现象、安全壳内初始环境条件、冷凝/蒸发现象显著影响安全壳压力。该研究结果为安全壳设计、安全分析和安全审评提供技术支持。     

基于抽样的敏感性分析方法在LBLOCA质能释放PIRT评级中的应用

基于随机抽样的非参量敏感性统计分析方法是一种有效的敏感性分析方法,通过计算热工水力分析程序多个抽样输入参数与输出参数之间的相关系数来评价各输入参数对输出参数影响的重要程度。通过耦合DAKOTA和WCOBRA/TRAC程序,开发了基于抽样的适用于非能动核电厂大破口失水事故质能释放的敏感性分析方法,该方法可全面定量评估各敏感性参数对计算结果的影响。计算结果表明：堆芯初始功率、燃耗、衰变热、安注箱初始水温、初始水体积、安注箱管道阻力系数、堆芯补水箱初始水温、喷放系数及破口阻力系数对破口质能释放具有显著影响。该分析结果可为大破口失水事故质能释放分析现象识别和重要度排序表评级提供定量依据。     

适用于核电厂BEPU的高效全局敏感性分析方法开发及应用

最佳估算加不确定性(BEPU)分析是IAEA推荐用于核电厂事故安全分析的方法,该方法中一个关键步骤为评估输入参数对目标输出的影响大小,即定量敏感性分析。传统BEPU分析中常使用基于线性或单调假设的局部敏感性分析方法,其难以适用于复杂的核反应堆系统,而全局敏感性分析则由于计算成本过高而难以在实际工程中应用。本研究中针对矩独立全局敏感性分析方法开展了优化研究,使用高阶模型表示、高斯求积公式等方法降低矩独立敏感性度量的计算成本,得到了一种高效的敏感性分析方法。使用了多个例题对优化方法的可靠性进行了验证,并将其应用于LOFT(loss-of-fluid test)大破口事故的敏感性分析。结果表明,该高效敏感性分析方法能准确识别核反应堆事故工况中的重要参数,并能对参数重要度进行定量排序。     

Iris small break loca phenomena identification and ranking table (PIRT)

The international reactor innovative and secure (IRIS) is a modular pressurized water reactor with an integral configuration (all primary system components – reactor core, internals, pumps, steam generators, pressurizer, and control rod drive mechanisms – are inside the reactor vessel). The IRIS plant conceptual design was completed in 2001 and the preliminary design is currently underway. The pre-application licensing process with the United States Nuclear Regulatory Commission (USNRC) started in October 2002.The first line of defense in IRIS is to eliminate event initiators that could potentially lead to core damage. If it is not possible to eliminate certain accidents altogether, then the design inherently reduces their consequences and/or decreases their probability of occurring. One of the most obvious advantages of the IRIS Safety-by-Design™ approach is the elimination of large break loss-of-coolant accidents (LBLOCAs), since no large primary penetrations of the reactor vessel or large loop piping exist.While the IRIS Safety-by-Design™ approach is a logical step in the effort to produce advanced reactors, the desired advances in safety must still be demonstrated in the licensing arena. With the elimination of LBLOCA, an important next consideration is to show the IRIS design fulfills the promise of increased safety also for small break LOCAs (SBLOCAs). Accordingly, the SBLOCA phenomena identification and ranking table (PIRT) project was established. The primary objective of the IRIS SBLOCA PIRT project was to identify the relative importance of phenomena in the IRIS response to SBLOCAs. This relative importance, coupled with the current relative state of knowledge for the phenomena, provides a framework for the planning of the continued experimental and analytical efforts.To satisfy the SBLOCA PIRT project objectives, Westinghouse organized an expert panel whose members were carefully selected to insure that the PIRT results reflect internationally recognized experience in reactor safety analysis, and were not biased by program preconceptions internal to the IRIS program.The SBLOCA PIRT Panel concluded that continued experimental data and analytical tool development in the following areas, in decreasing level of significance, are perceived as important with respect to satisfying the safety analysis and licensing objectives of the IRIS program: (1) steam generator; (2) pressure suppression system, containment dry well and their interactions; (3) emergency heat removal system; (4) core, long-term gravity makeup system, automatic depressurization system, and pressurizer; (5) direct vessel injection system and reactor vessel cavity.     

Assessment study about the risk of Wilks’ formula for uncertainty quantification of design extension condition scenarios in prototype Gen-IV sodium fast reactor

Sensitivity analysis and uncertainty quantification using Wilks’ formula and Monte Carlo for Unprotected Loss of Flow (ULOF) and Unprotected Transient OverPower (UTOP) accidents of prototype Gen-IV sodium-cooled fast reactor were performed. Multi-dimensional analysis for reactor safety for liquid metal reactors code calculations were conducted while simultaneously varying the values of all uncertain parameters according to their distribution using parallel computing platform integrated for uncertainty and sensitivity analysis to obtain uncertainty bands for Figures of Merit (FOM) – coolant, fuel centerline, and cladding temperature at the hottest fuel rod. To specify the uncertainty range of the parameters for each accident scenario, literature survey and expert judgments were consulted. By the sensitivity analysis, the importance ranking of 25 parameters in model identification and ranking table based on phenomena identification and ranking table was identified. Through Monte Carlo calculation, 95% upper limit and 95% confidence level were obtained, and about 2% and 5% under-prediction (risk) of FOM of ULOF and UTOP accidents using Wilks’ formula were confirmed, respectively.     

Coupled calculation of the radiological release and the thermal-hydraulic behaviour of a 3-loop PWR after a SGTR by means of the code RELAP5

To enable a more realistic and accurate calculation of the radiological consequences of a steam generator tube rupture (SGTR), a fission product transport model was developed. As the radiological releases strongly depend on the thermal-hydraulic transient, the model was included in the RELAP5 input decks of the Belgian Nuclear Power Plants. This enables the coupled calculation of the thermal-hydraulic transient and the radiological release. The fission product transport model tracks the concentration of the fission products in the primary circuit, in each of the SGs as well as in the condenser. This leads to a system of six coupled, first order ordinary differential equations with time dependent coefficients. Flashing, scrubbing, atomisation and dry out of the break flow are accounted for. Coupling with the thermal-hydraulic calculation and correct modelling of the break position enables an accurate calculation of the mixture level above the break. Pre- and post-accident spiking in the primary circuit are introduced. The transport times in the FW-system and the SG blowdown system are also taken into account, as is the decontaminating effect of the primary make-up system and of the SG blowdown system. Physical input parameters such as the partition coefficients, half life times and spiking coefficients are explicitly introduced so that the same model can be used for iodine, caesium and noble gases.     

Scaling methodology for a reduced-height reduced-pressure integral test facility to investigate direct vessel injection line break SBLOCA

A scaling methodology for a small-scale integral test facility was investigated in order to analyze thermal-hydraulic phenomena during a DVI (direct vessel injection) line SBLOCA (small break loss-of-coolant accident) in an APR1400 (advanced power reactor 1400 MWe) pressurized water reactor. The test facility SNUF (Seoul National University Facility) was utilized as a reduced-height and reduced-pressure integral test loop. To determine suitable test conditions for simulating the prototype in the SNUF experiment, the energy scaling methodology was propose to scale the coolant mass inventory and the thermal power for a reduced-pressure condition. The energy scaling methodology was validated with a system code (MARS) analysis for an ideally scaled-down SNUF model and that predicted a reasonable transient of pressure and coolant inventory when compared to the prototype model. For the actually constructed SNUF, the effect of scaling distortions in the test facility's thermal power and the loop geometry was analytically investigated. To overcome the limitation of the thermal power supply in the facility, the convective heat transfer between primary and secondary systems at the steam generator U-tubes was excluded and a modified power curve was applied for simulating the core decay heat. From the code analysis results for the actual SNUF model, the application of the modified power curve did not affect the major events occurring during the transient condition. The results revealed that the scaling distortion in the actual SNUF geometry also did not strongly disturb significant thermal-hydraulic phenomena such as the downcomer seal clearing. Thus, with an adoption of the energy scaling methodology, the thermal-hydraulic phenomena observed in the SNUF experiment can be properly utilized in a safety analysis for a DVI line break SBLOCA in the APR1400.     

A numerical investigation of SBLOCA scenario in nuclear heating reactor

Small break loss of coolant accident (SBLOCA) is one of the most important severe accidents in nuclear heating reactor. Nuclear heating reactor designed by Tsinghua University, whose primary loop is integrated layout and designed without main pump. The initial water volume in the reactor vessel is important to determine whether the reactor will be cooled or not as no safety injection system is designed for coolant makeup during the whole scenario. This paper simulates SBLOCA in nuclear heating reactor based on RELAP5. Transient behavior of relevant thermal parameters is specifically analyzed. Moreover, investigation also has been made on SBLOCA scenario based on different residual heat removal correlations and found the long-term residual heat removal capacity is decisive in determining the loss of coolant. The mathematical form of residual heat removal correlation is specifically deducted and can be widely applied to different situations. The envelope line that differentiates the region whether the core is safe or not under different maximum PRHRS capacity is also given.