压水堆核动力系统瞬态热工水力特性分析仿真软件.docVIP

压水堆核动力系统瞬态热工水力特性分析仿真软件 第31卷 201O 第2期 年4月 核动力工程 NuclearPowerEngineering Vb1.31.NO.2 Apr.2010 文章编号:0258—0926(2010)02-0047—04 压水堆核动力系统瞬态热工水力 特性分析仿真软件 巫英伟,庄程军,苏光辉,秋穗正 (西安交通大学多相流国家重点实验室,西安,710049) 摘要:采用点堆中子动力学模型,两相漂移流蒸汽发生器模型,三区不平衡稳压器模型,主循环泵四象 限特性模型和非能动应急余热导出系统模型,并利用CompaqVisualFortran6.0语言开发了微机型压水反应堆 瞬态热工水力特性分析程序,并利用MicrosoftVisualStudio.NE借言实现输入参数的可视化,输出结果的实 时处理和动态显示.利用RELAP5程序对本瞬态安全分析软件进行了可靠性验证,结果表明,本软件求解精 度较高,速度快,界面新颖,功能完善,可操作性强.此外,利用本软件对秦山核电站事故瞬态工况下的热 工水力特性进行了分析,得出了一些具有工程价值的结论. 关键词:压水堆;核动力系统;瞬态特性;输入可视化;动态显示 中图分类号:TL33文献标识码:A 1前言 与大型轻水堆热工水力分析程序相比,微机 型安全分析程序具有编程工作量小,程序容易修 改和补充,可完善性好,分析各种瞬态构成所需 要准备的数据量少,运算灵活方便,计算费用低 廉,对用户的专门技术要求低等优点.随着我国 核能事业的发展,国内在大型系列程序的引进, 开发和应用方面做了大量工作,但由于种种原因, 微机型程序的引进,开发和利用工作进行得很少. 因此,从核动力系统安全分析的发展趋势来看, 不失时机地开发微机型程序是十分必要的. 本文采用比较完善的压水堆核动力系统瞬态 数学模型,选用适合求解刚性微分方程组的吉尔 方法编制微机型核动力系统瞬态热工水力特性分 析程序,利用MicrosoftVisualStudio.NET语言实 现程序的输入可视化及输出结果的实时处理和动 态显示,利用RELAP5程序对本文开发的瞬态安 全分析软件进行了可靠性验证. 2系统模型 2.1堆芯模型 堆芯热工水力计算采用单通道模型;用具有 6组缓发中子,考虑燃料多普勒效应和慢化剂密 收稿日期:2008-12-26;修回日期:2009-06—24 度反应性反馈的点堆中子动力学方程描述堆芯内 中子动力学行为;用一维传导模型求解棒状元件 燃料芯块和包壳的温度场【l】. 2.2稳压器模型 本文所采用的稳压器模型在早期建立的稳压 器模型基础【2J上,考虑了各区在不同条件下所 处的状态以及在稳压器内发生的所有重要的热工 水力现象.稳压器中的工质分为3个区域:波动 水区,主水区和蒸汽区,并做如下简化假设: (1)在同一时刻,3个区域具有相同的压力. (2)同一区域内,相同工质在同一时刻具有 同一热力学参数. (3)忽略稳压器向外散热. 基于质量,动量和能量守恒方程,可建立稳 压器的水位方程,压力方程,主水区质量方程, 波动水区质量方程,主水区焓方程和波动水区焓 方程等. 2.3蒸汽发生器模型 针对立式u型管束自然循环蒸汽发生器,采 用具有漂移流模型的一维两流体蒸汽发生器模 型,建立蒸汽发生器一,二次侧流体焓方程以及 二次侧压力方程,水位方程,二次侧再循环流量 方程等.具体的数学物理模型见文献【7】. 核动力工程V.o1.31.No.2.2010 2.4非能动应急余热排出系统模型 非能动应急堆芯余热排出系统的重要功能是 在反应堆失去正常的冷却手段,即反应堆冷却剂 系统.蒸汽发生器.二回路系统.汽轮机做功和冷凝 器向环境释热这一堆芯冷却链遭到破坏时,能保 证堆芯的余热排出,使堆芯的温度维持在规定的 安全限值之内.本文采用一维两流体二次侧非能 动应急堆芯余热排出系统模型计算分析该系统的 热工水力特性J. 2.5主循环泵模型 描述主循环泵的主要参数有泵的扬程,转矩, 体积流量和角速度.通常把由实验得出的这些主 要参数之间对应关系的曲线称为泵的四象限曲 线.然而,这种形状的曲线很难应用于程序分析 中,故其转换成一种较简单的类比曲线.类比曲 线是以扬程比和转矩比(真实值与额定值之比) 的形式画出的,是泵的转速比和容积流量比的函 数.该曲线以表格形式输入,因变量作为自变量 的函数由表格查找或线性内插获得,并应用到微 机型程序中. 3软件编制 3.1程序编制 在建立合理的数学物理模型基础上,采用 CompaqVisualFortran6.0程序设计语言开发了压 水堆核动力系统瞬态热工水力特性分析程序;程 序仿真流程如图1所示.程序完全采用模块化结 构设计,在计算方法的选取上,程序可随时判定 每一时刻微分方程组刚性的强弱,自动选取阿当