基于RELAP5与MELCOR联合分析方法的压水堆严重事故研究

针对严重事故的模拟研究,本文提出结合热工水力系统程序和严重事故一体化程序的分析方法,以典型三环路传统压水堆为对象,分别采用RELAP5和MELCOR程序建立模型,分析在全厂断电叠加汽动辅助给水泵失效事故下系统的瞬态响应。为了尽可能地利用RELAP5计算早期热工水力响应,同时保证严重事故计算结果的准确性,以MELCOR锆合金氧化模型开始工作温度的下限,即包壳温度达到1 100 K作为程序衔接准则并利用RELAP5的大编辑功能,提取所需计算结果导入MELCOR输入卡作为初始参数继续模拟。计算结果表明,数据连接过程整体保持了连续性,两种方法计算得出的主冷却剂系统压力、堆芯和稳压器水位、燃料包壳温度等参数的数值以及堆芯传热恶化和压力容器失效等现象的时序存在不同程度的差异,例如堆芯熔毁时间延后了约538 s。由于采用了RELAP5计算严重事故前的系统暂态响应,联合分析方法的计算结果比单独使用MELCOR分析的结果更加准确,该方法可以提高传统严重事故分析的可靠性。     

基于堆芯模型耦合过渡方法的严重事故模拟机开发研究

采取系统分析程序耦合过渡一体化严重事故（SA）分析程序的方法，对严重事故模拟机的开发进行研究。该方法首先使用系统分析程序计算事故早期响应，当满足耦合条件时，系统程序停止计算，切换至严重事故程序计算模拟事故中晚期。为实现切换时参数平滑过渡，以全范围模拟机常用程序RELAP5和严重事故程序MAAP4为例，主要分析了两程序热工水力模型重叠部分的堆芯区域的物理模型，选择传递了堆芯节点的芯块温度、包壳温度和堆芯功率。基于通用百万千瓦级压水堆小破口失水事故（SBLOCA）模型，使用该方法计算和SA程序单独计算进行对比验证。结果表明，过渡参数的选取是正确的，该系统分析程序耦合过渡SA程序的方法不仅能成功平滑地过渡参数，还保证了后续计算的准确性。      

SCDAP/RELAP5与MELCOR程序对堆芯损伤过程预测的比较

SCDAP／RELAP5与MELCOR程序是目前得到广泛使用的两个严重事故分析程序．它们在模拟堆芯溶化及压力容器下封头失效过程中采用了基于不同理论的计算模型。本文利用两个程序分别对秦山二期核电厂发生假想的全厂断电事故下的堆芯损伤过程进行预测．并对比分析了这2个严重事故分析程序的优点及相应的计算结果．     

非能动压水堆热工水力多尺度耦合计算分析研究

本文为了更加真实准确的模拟非能动核电机组复杂的热工水力工况,提高事故分析计算精度,开展了非能动压水堆热工水力多尺度耦合计算分析研究。首先应用热工水力系统分析程序RELAP5对AP1000机组进行系统建模,并开展冷却剂强迫流动完全丧失事故(全失流事故)的分析计算,得到堆芯相关热工水力参数。然后将RELAP5程序的计算结果作为边界条件,分别利用子通道程序COBRA-Ⅳ、计算流体力学程序FLUENT以及基于两个程序的耦合程序对AP1000堆芯组件进行建模,并分别开展全失流事故过程中堆芯热工水力分析计算。最终通过三个程序计算结果的对比,表明应用耦合程序开展堆芯热工水力分析的方法可行,建立的堆芯组件模型合理,计算结果更加接近真实情况,有效减少了单一程序计算的过度保守性。     

三维瞬态物理分析模型与RELAP5/mod3.2程序的耦合研究

在系统热工水力程序RELAP5/mod3.2的基础上,采用显式方法建立了堆芯三维时空中子动力学与一维热工水力计算的耦合模型,接入基于非线性迭代半解析节块法的三维瞬态物理分析模型(NLSANMT)后,形成了一个具有堆芯三维瞬态物理特性分析能力的系统计算程序NLSANMT/RELAP5(mod3.2).通过核动力反应堆温度反馈系数、堆芯功率分布参数的校算及单束控制棒失控抽出事故的模拟分析,验证了接口的正确性.验证结果表明,与RELAP5/mod3.2相比,所开发的NLSANMT/RELAP5(mod3.2)程序具有更强的堆芯物理瞬态分析能力.     

大功率非能动压水堆DVI管破叠加IRWST失效触发严重事故分析

应用MELCOR 2.1程序,建立了大功率非能动压水堆核电厂主要回路系统及安全壳的热工水力模型,并以直接注水管线破口叠加内置换料水箱失效触发严重事故为对象进行了独立计算。计算结果与MAAP 4.04程序计算结果趋势一致,分析表明:MELCOR 2.1新版本对严重事故计算合理可信;部分非能动安全设施的启动有效地降低了主回路系统压力,防止高压熔堆,缓解了堆芯熔化进程,从而验证了非能动安全设施的有效性。     

压水堆核电站堆芯物理/热工水力耦合特性研究

采用RELAP5-HD作为堆芯耦合计算程序,以秦山核电二期工程反应堆堆芯为研究对象,建立堆芯活性区的物理/热工水力耦合模型,在此基础上进行了稳态计算和掉棒事故仿真研究。结果表明,使用RELAP5-HD计算得到的结果与电厂实测值符合较好,获得的掉棒事故参数曲线能准确反映事故工况下的参数变化趋势。稳态和事故工况的计算结果均符合堆芯物理/热工水力反馈效应的理论分析,证实了所建立的堆芯耦合模型的准确性,为下一步进行核电站系统的仿真分析提供基础。     

RELAP5/FLUENT耦合程序的开发

热工水力数值模拟是反应堆系统设计和安全分析的重要内容,以RELAP5为代表的系统程序可对瞬态或事故工况进行快速分析,同时以FLUENT为代表的计算流体动力学(CFD)程序对堆芯局部三维现象的分析也越来越重要。为综合利用两者的优点,以RELAP5/FLUENT为基础,利用对RELAP5程序源代码的二次开发和FLUENT的用户自定义函数(UDF)进行编程,开发了RELAP5/FLUENT耦合程序。利用flibe熔盐在水平圆管流动问题验证了程序耦合的正确性;针对2 MW熔盐堆进行了稳态模拟,耦合程序能详细分析熔盐堆的热工水力行为;模拟了2 MW熔盐堆功率突变的瞬态热工水力行为,相对于单独的RELAP5,耦合程序能更好地揭示熔盐堆系统和堆芯的三维物理现象。该耦合程序可用于解决熔盐堆热工水力分析中存在的显著三维混合现象的问题。     

用PARCS/TRACE/ROBIN程序系统研究秦山二期弹棒事故

利用美国核管制委员会(US NRC)堆芯三维中子动力学软件PARCS、热工水力软件TRACE、辅助建模软件SNAP以及具有国内自主知识产权的压水堆燃料组件计算软件RONBIN,建立了秦山二期两环路压水堆物理模型和热工水力系统模型,进行弹棒事故模拟计算,得出合理的计算结果。AFA 3G燃料组件的两维中子输运计算由ROBIN程序完成,生成的宏观中子截面参数被传递给PARCS程序作为输入。然后由PARCS程序进行堆芯三维弹棒模拟计算,得到事故过程中的核功率变化趋势。最后将反应堆功率瞬态数据输入TRACE热工水力系统模型计算系统压力响应以及燃料包壳和芯块温度。本文通过使用与设计单位完全不同的软件体系,独立地验证了该堆型在弹棒事故下的安全性。     

基于修改后RELAP5/MOD3.2的氟盐高温堆实验回路分析

大型热工流体整体效应系统实验(CIET)台架是为模拟氟盐冷却高温堆(FHR)热工水力响应而设计的实验回路,采用DOWTHERM A模拟氟盐作为冷却剂。通过在RELAP5/MOD3.2程序中加入DOWTHERM A物性参数以及传热关系式,计算FHR实验回路CIET在两种工况下的热工水力行为,并与实验结果进行对比,计算工况包括强迫循环条件与自然循环条件。计算结果表明:在强迫循环条件下,堆芯热量主要靠盘管式空气换热器(CTAH)排出,堆芯进出口冷却剂温度及CTAH出口冷却剂温度与实验值符合良好,CTAH进口冷却剂温度与实验值有些微偏差;在自然循环工况中,堆芯热量主要通过DHX与堆芯辅助冷却系统(DRACS)回路的换热带走,DHX及DRACS的流量与实验值接近,相对误差在10%左右,验证了修正后RELAP5/MOD3.2的正确性。     

Analysis of pressurizer level control system using MAAP4 code

The code initialization effort has been troubling code users for decades for system transient and severe accident analyses using codes such as RETRAN, MAAP4, MAAP5 and MELCOR. The purpose of this work is to demonstrate an approach that could be considered a generic method to address the code initialization problem. This was demonstrated by developing a pressurizer level control model and temperature dependent level control logic in MAAP4 without re-compiling with the source code. The method would enhance the simulation capability and accuracy of a severe accident analysis by transient and severe accident analyses codes. The demonstration case used MAAP4 to show that the adopted proportional-integral controller with the temperature dependent level control logic would reduce its code steady state errors to zero. The subsequent transient response would become more realistic. The proposed method provides a convenient and exemplified approach for code initialization which is applicable to the next generation of codes that couple with the balance of plant models. These codes include the MAAP5 code and others future codes that could simulate the whole plant by a single and elaborate plant model with exhausting component and phenomenological models.     

Simulation of the small modular reactor severe accident scenario response to SBO using MELCOR code

Advanced small modular reactors (SMRs) use different design in the systems, structures, components from large reactors for achieving a high level of safety and reliability. In present work, the SMRs severe accident caused by the station blackout (SBO) was modeled and analyzed using MELCOR code, and the simulation of the accident scenario response to SBO was conducted. Based on the steady state calculation, which agrees well with designed values, we introduced the SBO accident for transient calculation. First, the case of the SBO accident without the passive core cooling system (PXS) was calculated. The progression and scenario in the reactor pressure vessel (RPV) and the containment were simulated and analyzed, including the transient response, cooling capacity and thermal-hydraulic characteristics and so on. The station black-out transient in the SMR can be simulated accurately, and the main failure model in the accident process can be concluded. Then three other cases of the SBO accident with different passive safety systems (core makeup tank (CMT), accumulator (ACC), passive residual heat removal system heat exchanger (PRHR HX), automatic depressurization system (ADS)) of the PXS were calculated respectively, and the results for different passive safety systems were compared. The passive core cooling system can not only provide water to the primary coolant system, but also take away the reactor decay residual heat. So in a station black-out transient, we can get more time for restoring AC power, and effectively prevent the accidents such as Fukushima.     

严重事故分析程序

严重事故过程可分为堆芯解体、压力容器熔穿、安全壳失效三阶段，裂变产物随之释放。严重事故分析程序有两类，系统性程序与机理性程序。系统性程序能计算完整的事故序列，机理性程序偏重事故进程的局部细节。主要介绍系统性程序STCP、MELCOR、ASTEC与机理性程序RELAP5／SCDAP、VICTORIA、CONTAIN。     

大破口触发的严重事故分析及缓解措施研究

采用严重事故一体化分析程序MELCOR,对国产先进压水堆核电厂进行系统建模,选取大破口触发的严重事故进行校核计算研究,获得了严重事故工况下核电厂关键参数的瞬态特性和非能动系统响应特性,并与安全分析报告中MAAP的计算结果进行了对比分析。结果表明:虽然校核计算结果与安全分析报告中的结果存在一定差异,但总体上事故序列和主要参数的变化趋势吻合良好,并且都能够在严重事故情况下保持压力容器和安全壳的完整性,放射性裂变产物释放量极低,缓解措施的设计能够有效缓解事故进程,满足核电厂的安全要求。     

Simulation of the chernobyl accident

An analysis of the April 26, 1986 accident at the Chernobyl-4 nuclear power plant in the Soviet Union is presented. The peak calculated core power during the accident was 550 000 MW_t. The analysis provides insights that further understanding of the plant behavior during the accident. The plant was modeled with the RELAP5/MOD2 computer code using information available in the open literature. RELAP5/MOD2 is an advanced computer code designed for best-estimate thermal-hydraulic analysis of transients in light water reactors. The Chernobyl-4 model included the reactor kinetics effects of fuel temperature, graphite temperature, core average void fraction, and automatic regulator control rod position. Preliminary calculations indicated the effects of recirculation pump coast down during performance of a test at the plant were not sufficient to initiate a reactor kinetics-driven power excursion. Another mechanism, or “trigger” is required. The accident simulation assumed the trigger was recirculation pump performance degradation caused by the onset of pump cavitation. Fuel disintegration caused by the power excursion probably led to rupture of pressure tubes. To further characterize the response of the Chernobyl-4 plant during severe accidents, simulations of an extended station blackout sequence with failure of all feedwater are also presented. For those simulations, RELAP5/MOD2 and SCDAP/MOD1 (an advanced best-estimate computer code for the prediction of reactor core behavior during a severe accident) were used. The simulations indicated that fuel rod melting was delayed significantly because the graphite acted as a heat sink.