relap5介绍..docx

第二章 RELAP5介绍2.1 RELAP5功能简介RELAP5是Idaho国家工程实验室（INEL）为核管会（NRC）开发的轻水堆瞬态分析程序，现已经成为核电厂分析器的基础。RELAP5几乎可以覆盖轻水堆核电厂所有热工水力瞬变和事故谱，其研究和分析领域可运用于轻水堆的设计基准事故以及瞬态和稳态的各种工况模拟。例如：LOCA(失水事故)、LOFT（失流事故）、给水丧失、意外失去厂外电源、以及汽轮机跳闸等等。图2.1.1 RELAP5的功能图RELAP5的部件模型包括管道、阀、泵、安注箱、汽轮机、环路、应急堆芯冷却器、蓄压箱、分支和控制系统部件等（见图2.1.2）。RELAP5还对截面突变、交叉流、壅塞流（choke flow）、混合水位跟踪、不凝气体迁移等现象建立了特殊分析模型（见图2.1.3）[16]。图2.1.2 RELAP5的组件模型图图2.1.3 RELAP5的特殊过程模型图RELAP5程序里还包括了一些简单的模型：热力平衡态模型、均相流模型以及无摩擦流动模型等。这些模型可供使用者单独选用或组合选用[7]。RELAP5以两相系统的非均相和不平衡态模型为基础，采用快速半隐式数值方法进行求解[7]。RELAP5的发展情况：自1967年RELAP5程序开始研制，经过20多年，出现了今天的RELAP5/MOD3版本。这个软件集中了人们在两相流理论研究、数值求解方法、计算机编程技巧和各种规模实验等方面取得的研究成果[7]。RELAP5已经由许多实例和各个国家相关的实验计划的实验数据进行了验证和比较，但目前它依旧在扩大和发展之中，以便能在新的先进压水堆中得到更好的应用。接近1995年末时,为了维持一个单一版本程序源的显示效率，美国核管会（NRC）及美国能源部门（DOE）都必须要额外的努力去克服一些相冲突的要求。这个单一版本程序源便被划分成两种版本,一个是NRC版,另一个是DOE版。DOE版本保持了原来程序所有的功能,再加上划分之后所提供的功能。目前见到的最新版本有RELAP5/MOD3.4和RELAP5/MOD4.0。2.1.2场方程图2.1.4 RELAP5中的场方程RELAP5程序要求求解八个一次因变量的场方程。这8个一次因变量是：液相比内能Uf，压力P，蒸汽体积份额(空泡份额)αg，汽相比内能Ug，液相速度Vf，汽相速度Vg，硼密度ρb和不凝气体质量含量Xn。场方程中还有二次因变量，它们分别为相温度Tg和Tf，相密度ρg和ρf，每一种不凝气体质量份额和饱和温度Ts [7]。RELAP5程序中的基本场方程包括两相动量守恒方程、两相质量连续性方程和两相能量守恒方程。这些方程以微分形式出现，自变量是时间和一维空间，因变量是按时间和体积求平均的量（见图2.1.4）[7]。质量守恒方程：汽相 （2-1）液相 （2-2）一般不考虑质量源或热阱，因此有Гg=-Гf 动量方程汽相（2-3）液相（2-4）RELAP5程序中的动量方程做了一些简化：假定各相的压力是相等的，假定界面的压力等于相压力（层流除外），不考虑应变项，忽略各相的粘性应力，忽略相界面上的动量储存，相界面上的力项包括压力和粘性应力，用可变面积动量流关系式来建立壁面力模型[7]。上面动量方程右边的各项为：压力梯度、体积力、壁面的摩擦力、界面上的质量转移引起的动量转移、界面的摩擦力、虚拟质量力。FWF和FWG项是壁面摩擦阻力项，这里把它们考虑成与相速度成正比的线性项，它们是单位体积内的摩擦面积、壁面摩擦系数和相速度大小的乘积。FWF和FIG项是界面摩擦阻力项，这里把它们考虑成与相对速度成正比的线性项，它们是单位体积内的摩擦面积、界面的摩擦系数以及相对速度的乘积[7]。上式利用RELAP5/MOD1中采用的形式对虚拟质量项进行处理。而在RELAP5/MOD3中，这项做了简化，特别是略去了对空间变化项。能量守恒方程汽相 （2-5）液相 （2-6）上面能量方程右边的各项为：膨胀或压缩功，内能变化率，Qwf，Qwg是每单位体积的壁面热交换量，Qif，Qig是界面上的热交换量，Гighf*，Гighg*是界面的质量迁移引起的能量迁移，，为壁面（热边界层）质量迁移引起的能量迁移，DISSg、DISSf是相能量消耗项[7]。 ， 2.1.3 状态方程在上面三大守恒方程中采用五个独立变量，它们是压力P 、蒸汽体积份额ag 、汽相比内能Ug 、液相比内能Uf和不凝气体质量含量xn 。其余的热力学变量(如温度、密度、局部压力、含汽量等)则可用这几个独立变量的函数来表示。为使方程线性化，除这些变量外, 还需要一些状态变量，表现为状态关系式，用以确定状态参数（见图2.1.5）。例如 图2.1.5 状态关系式图2.1.4 传热模型的分析在RELAP5/MOD1版本中，传热模型共分为