relap5英文手册阅读笔记.docx

2 基本规则relap5是专门计算单相或两相流动系统的程序之一。在选择一种计算程序之前，要考虑该程序是否有能力处理和分析要解决的问题，对计算机有什么特殊要求，是否经济，以及怎样更好地使用该程序。2.1 relap5的功能 relap5用于轻水热工水力系统。最初是为了计算压水堆破口失水事故而编写的，经过不断地发展，现在可用于计算许多热工水力系统的瞬态过程。如：破口失水事故、运行瞬态。可用于生产堆、实验堆和反应堆模拟器。还可以用于沸水堆。relap5用六个非平衡、非均匀的方程模拟两流体系统。同时可模拟有硼和不凝气体存在的情况。中子动力学用点堆动力学方程。联合场方程和点堆方程可模拟反馈效应。可以模拟板材和流体间的传热过程。控制系统可以模拟设备的控制器，部件模型可以模拟核电厂配套设备（如汽轮机、泵、冷凝器等）。2.1.2 relap5主要用于计算运行和事故瞬态中流体的行为特性。计算出的流体瞬态压力原则上可用于部件的结构载荷分析（如：水锤、汽锤），还可以计算声波（压力波）在管道和设备内的传播。但relap5很少用于结构分析，因此，要谨慎对待计算结果，并且尽量避免把relap5应用于这些场合，因为relap5不是计算上述问题的最佳程序。2.2 热工水力分析的描述用relap5解决问题的基本步骤：（1）收集和整理描述边界条件和初始条件的详细信息。有些信息需要由专门的部门和机构提供。（2）详细描述问题并将问题节点化。需要成立专门小组来评价模型的合理性。（3）输入问题。将初始条件和边界条件输入到一个计算机文件中。（4）检查对问题的描述和输入。检查初始信息的来源是否合理，导出量的计算是否，对模型的简化是否合理。（5）运行程序，分析计算结果。2.2.1 收集和整理信息relap5的输入信息可分为四个方面：水力学、热构件、控制系统、中子动力学。水力学信息包括：（1）流动面积（2）流动长度（3）竖直取向（4）足以计算水力直径的几何条件（5）壁面粗糙度（6）计算流动损失的信息。例如：弯管几何形状、面积突扩形状、阀的几何形状、阀的额定流量、设备启动试验数据等。（7）设备初始条件（8）泵的参数热构件信息包括：（1）材料的厚度（2）材料的质量（3）管长（4）内外径（5）材料类型（6）材料的热物性随温度的变化规律。如：热导率、密度、比热等。（7）加热功率（如有热源项）（8）热源分布（9）初始温度分布控制系统信息包括：（1）控制系统方块图（2）控制系统与被控对象的关系（3）控制器特征（4）滤波器特征（5）设定值（6）增益（7）饱和限（8）滞后（9）被控设备特性（10）阀的行程速度（11）最大（最小）泵的转速（12）最大（最小）循环速度中子动力学信息包括：（1）初始反应性（2）照射数据（3）缓发中子份额（4）裂变产物产额（5）锕系元素产额（6）反应性特征2.2.2 详细描述问题并将问题节点化 将需要建立模型的系统隔离开来，明确其边界。系统的模型由若干控制体（节点）组成。建立各个控制体的过程就是将问题节点化的过程。建好的模型应由其他人检查以确保其正确性。2.2.2.1 确定系统边界2.2.2.2 将模型节点化 将模型划分为若干控制体，这些控制体主要是一些有入口接点和出口接点的流管。控制体由接管（junction）连接。程序计算控制体中心的平均流动特性和接点处的流动矢量特性。 节点的尺寸受计算数值稳定性、运行时间、空间收敛性的影响。一般说来在结果满足空间收敛性的前提下节点的尺寸越大越好。因为节点的尺寸直接影响运行时间；节点越小，最大时间步长越小，这是为了保证数值稳定性。模型中除了流管外，还有一些静止的构件，这些构件可以储存热量，有的甚至还含有内热源。程序可以模拟这些构件与流体之间的一维热传递。模型节点的划分非常重要，因此在考虑到所研究问题的特殊性的情况下要充分借鉴前人的研究成果。结果应由检查人员认真检查。划分节点的过程与计算模型各部分参数的过程同步。因此节点划分完成后，模型各部分的参数也就已知了。2.2.3 获得边界条件和初始条件2.2.3.1 输入数据程序可以检查输入错误，但不能仅依靠程序本身去检查，因为许多输入是程序无法检查的。程序能够检查模型环路的高度是否封闭，即组成环路的各个部件的高度变化之和是否为零。如果不为零，将会自动加上一个虚拟的泵，泵的压头等于环路高度不封闭而产生的静压差。通过程序的初始化来检查这种不匹配。2.2.3.2 系统的稳态系统的稳态计算非常重要，它是瞬态计算的基础。模型的稳态被设定为实际系统的初始工况。程序含有“稳态”选项，它使系统通过减少材料热容的方式很快达到设定的稳态点。2.2.4运行程序分析结果2.2.4.1运行程序运行程序时可能出现的问题：（1）出现不能被解释的错误。例如结果提示某个节点的压力或温度不符合实际或程序不收敛（converge）。这些问题有事可通过减小最大时间步长为几个时间步长来解决。（2）结果不趋同（convergent）。计算结果与前人的结果不同。2.2.4.2 解释结果目的：（1）检查程序是否收敛； （2）检查是否存在非自然规律（nonphysical）的结果。 （3）检查结果是否符合常理。 （4）结果的每一个部分应该都能得到合理的解释。2.3 relap5的基本建模元件实际系统有管道、阀门、泵和其他设备组成，类似地，relap5模型也是如此。relap5的基本建模块可分为四类：热工水力、热构件、trips和控制变量。热工水力部件用于模拟流体通道和流体处理设备。热构件用于模拟材料块和材料块与流体之间的相互作用。trips用于模拟各种设备的触发信号（如在特定时间开启泵；阀门在压力高于某一值时打开，低于某一值时关闭等）。控制变量用于模拟控制系统（如比例、积分、微分控制器；超前、滞后控制器）和集总节点（lumped-node）系统。热工水力建模块：单一控制体 ：代表一部分流管。管型（环型）：代表一段管道。内部可以有1到100个控制体，内部控制体由内部接管相连。 环型：一种特殊的管型，特点与管型相同，专门用于模拟环形流动通道。分支：代表流管的接合处。最多可规定10个接管。分离器：特殊的分支。用于模拟蒸汽发生器中的汽水分离器。混合喷射器：特殊的分支。用于模拟喷射泵。透平：特殊的分支。用于模拟汽轮机。eccmix：单一接管：用于连接两个部件。时间相关控制体：用于指定系统模型的边界条件。时间相关接管：连接两个部件，同时指定接管边界条件。阀：一种特殊的接管部件。用于模拟六种阀。即：制止阀、触发阀、惯性阀、马达阀、伺服阀、释放阀。泵：用于模拟离心泵。蓄压箱：用于模拟压水堆的蓄压箱，包括容器和管线。控制变量建模块：和差：对变量进行加减运算。乘：对变量进行乘法运算。除：对变量进行除法运算。求导：变量对时间求导。积分：变量对时间积分。函数：对变量定义一个表格函数关系。标准函数：对指定的变量进行以下运算（绝对值、平方根、指数、自然对数、正弦、余弦、正切、反正切、最小值、最大值）延迟：代表指定变量的时间延迟因子。单位trip：在指定时刻为真。也可以定义为互补函数。整数次幂：求变量的某一整数次幂。实数次幂：求变量的某一实数次幂。变量次幂：求变量的某一变量次幂。比例积分：定义一个比例积分控制器。滞后：定义一个滞后控制器。超前滞后：定义一个超前滞后控制器。常量：定义一个和其他控制变量一起使用的常量。轴：定义与发电机连接的轴的特性。泵控制：定义一个泵控制器（主要用于稳态分析）。蒸汽控制：定义一个蒸汽流动控制器（主要用于稳态分析）。给水控制：定义一个给水流动控制器（主要用于稳态分析）。2.3.1 热工水力组2.3.2 热构件热构件用于模拟系统的结构材料。如外部包有隔热层内部镀有不锈钢的铸铁管道。也可以模拟由内热源的构件，如核燃料和电加热元件。热构件的内部可以在一个方向上划分节点，以求得内部的温度分布。程序默认热传递方向与流动方向正交（垂直）。2.3.3 trips和控制变量trip可以在瞬态计算时指定动作。与控制变量共同使用时可以大大扩展程序的功能。trip可用于以下部分：时间相关控制体、时间相关接管、泵部件、阀部件、分支、蓄压箱还有描述反应堆动力学特性、热构件特性的表格。控制变量的作用可归结为三点：（1）模拟设备的控制系统（2）构造集总节点参数（3）增加热工水力系统和热构件的边界条件的维数。2.3.3.1 模拟设备的控制系统2.3.3.2 模拟集总节点系统2.3.3.3 增强模型的边界条件2.4 relap5基本建模指南2.4.1模拟系统的流动通道2.4.2 模拟系统的热构件2.4.3 trip和控制变量 变量trip前加负号表示“非”。3 总体规则3.1 标准步骤标准步骤包括：输入准备，调试模型的输入，执行问题，输出解释。3.1.1 输入准备31.1.1 模型记录和质量保证 要准备一个工作簿。其中记录建模的所有信息，包括参考资料、基本假设、计算数据等。包括最后的总的节点图。节点的编号要有规律。比如一个三回路的压水堆的第一个回路可以用100-199编号，二回路用200-299编号，三回路用300-399编号。并且要注意标号的对称性，例如120、220、320代表不同回路的同一部件。以上措施均有助于检查错误。部件号共9位，结构为xxxcc0000。xxx为部件号，100-999，应避免使用1到99为部件号，这样可能会导致错误。cc为部件的内部元件号。0000为尾部零，尾部零是为了提供一种定义多维部件的方式，尽管多维部件到目前为止还没有被加入到程序中。工作簿的工作表中记录所有输入信息，记录顺序应与输入顺序形同。对任意一个部件，记录顺序为水力学数据，热构件数据，控制变量数据。对任一部件，工作簿的记录为（1）建模数据来源（例如图形、报告等）（2）假设（3）对原始数据的计算步骤（4）最终要输入的参数的值。工作簿中的信息直接用于输入计算机文件来建模，输入后将计算机文件与工作簿对照已以检查是否有错。输入后由另外的人员对整个模型的检查对保证模型的正确性是非常重要的，虽然检查过程非常枯燥和乏味。检查者应认识到这项工作的重要性，并有积极的态度。3.1.1.2 输入卡的安排 每行数据的开头必须有一个卡号，其后是数据列，数据之间用一个或几个空格分隔。如果一个卡号被重复，程序将把它作为代替卡，并用最后一个卡的信息。 一个完整的输入应包含注释信息。注释信息在*或$之后，系统自动把星号和美元符之后的内容当做注释内容。 推荐的输入顺序为：题目，工作控制，时间步长控制，小编辑要求，trip说明，水力学部件，热构件，输入的数据表，控制变量，反应堆动力学说明。一般应按照卡号由小到大的顺序输入。3.1.1.3 规则 热构件号和控制变量号最好与相关的水力学控制体的部件号相同。 系统正常运行时，接管的方向最好与流动方向相同。 程序默认的热传递方向：由构件到流体为正，由流体到构件为负。例如，蒸汽发生器传热管的一次侧的热传递方向为负，二次侧为正。 建议将热构件的左表面作为内侧，右表面作为外侧。3.1.2 模型输入调试 通常一个模型会产生一些错误，有的错误还会进一步引起其他错误。因此显示的错误信息一般要多于实际的错误数。一些错误让人难以理解。比较有效的方法是先改正一些比较明显的错误，然后再次运行程序。错误信息的前面有8个星号“*”。因此，使用者应尽量避免使用一串星号来分隔输入文件。所有错误移除后显示“input processing completed successfully.”检错过程实际上分为几个阶段。有些情况下，检查到某种错误后检查过程会自动停止，此时显示的错误是已检查过的部分出现的错误，直至这些错误被移除后检错过程才能继续，新的错误会继续被发现。因此，有时，错误并不是一直减少的。检错过程包括检查环路的高度是否封闭，如果系统任一个环路的高度差大于0.0001米，将会显示错误信息。有一个虚假的错误曾被碰到过。在输入处理结束后显示程序执行错误，算法错误或算法溢出。这通常是因为错误的对不凝气体进行说明造成的。使用者在初始条件控制字中指定了不正确的值。控制字的4、5、6字是留给有不凝气体存在的情况，如果使用了这些字而在110卡中有没有指定不凝气体，就会出现错误。3.1.3 执行问题如果代码成功通过输入处理阶段（所有错误都被移除），将会生成初始时间编辑。如果选择“run”，程序接着从指定的初始状况开始执行。对new问题，初始时间是0时刻，对restart问题，初始时间为重新编辑时间。3.1.3.1 时间步长与编辑选择问题的执行受201-299时间步长控制卡的控制。每一个时间间隔都有一个最小时间步长和一个最大时间步长。程序按最大时间步长开始执行，时间步长根据一些测试自动减小。材料的库兰特限制不能被违背。在每一个计算元件，质量、流体特性和推断的错误将被监测，如果错误超出了内部预先设定的限制，时间步长将会减小。大编辑输出表明准则和模型区域造成时间步长的减小。时间步长通过不断被2除的方式减小，直至错误达到可以接受的限制范围内，或者达到最小时间步长，或者遇到错误。好像开始就设定一个较小的最大时间步长对计算是有利的。不可能推荐一个通用的最大或最小时间步长，时间步长的设定要根据模型的特点，要解决的问题和一些时间步长对计算的影响的研究来完成。但一个非正式的建议是：最小时间步长用110-7秒，最大时间步长用库兰特限制（但不要大于0.2秒）。时间步长控制中ssdtt通常设定为00003（或3）。小编辑频率、大编辑频率和重新启动编辑频率也在时间步长控制中设定。这些值是最大时间步长的整数乘子。例如，最大时间步长为0.1秒，小编辑频率为10，大编辑频率为100，重新启动编辑频率为200。则程序每1秒生成一个小编辑，每10秒生成一个大编辑，每20秒生成一个重新启动编辑。小编辑频率的选择非常重要，因为他代表了输出文件数据点的频率，这些点之间的数据点将不能得到。所有计算数据（如压力、温度、空泡份额、流速等）都可以在输出文件中得到，但并不是所有的计算数据都要在小编辑要求中指定，小编辑要求中指定的只是打印输出的参数。3.1.3.2 稳态、瞬态、strip方式计算方式在100卡中设置。在稳态计算中，为了减少运行时间，所有热构件的热容被人为的减小。当内部测试到参数的变化率已很小时，计算停止。一般情况下，稳态方式不推荐使用。实际上，稳态计算通过在瞬态计算模式下控制边界条件的方式进行。通过这种方式，再人为的减少热容，可以得到满意的稳态计算。在瞬态计算前，再把热容恢复到实际值即可。strip模式用于在已有的计算结果中提取特定的数据通道。3.1.3.3瞬态执行失败引起瞬态执行失败的一个常见原因是指定的初始条件和边界条件不合适。例如，一个冷却剂丧失事故的瞬态计算是从计算好的稳态状况开始的，计算刚开始就遇到失败。此时应该检查初始条件和边界条件是否合适。如果没问题，下一步就要减小最大时间步长。3.1.3.4瞬态执行的正常结束瞬态执行通常以下列三种方式结束：（1）问题时间达到在201-299卡中设定的结束时间。（2）花费的计算机时间达到在105卡中指定的限制。（3）被trip终止，正如在600卡中指定的一样。因此，使用者可以通过问题时间、计算机时间和计算过程中的事件（如压力超过限定值）灵活地结束计算。使用者确保以以上三种方式之一正常结束计算是很重要的，否则会出现不正常的结束，所有计算数据都可能丢失。3.1.4 输出有两种形式的计算结果输出：打印输出、重新启动/图表输出。3.1.4.1 打印输出输出文件的开头是输入数据列表，然后是输入处理过程。输入处理中会再次出现输入数据，这些数据是程序真正引用的数据。如果出现重复的卡号，用最后一个卡上的数据。输出文件中包含初始大编辑，接着是使用者要求的大编辑和程序生成的大编辑。使用者在301-399卡中要求的小编辑分散在大编辑中间。另外还有警告信息，里面包含一些不重要的计算困难。再启动编辑由再启动号标出，在大编辑数据之后。以下是一个大编辑输出的例子：1 tot代表初始大编辑，edit代表上一个大编辑；程序以最大时间步长（要求时间步长）计算，计算过程中时间步长没有减小，510-3秒。库兰特限制为6.3910-2秒，花费的cpu时间为68.6655秒，质量减少为-0.207kg,总质量为514348kg质量误差为10-7数量级。2 这个区域显示所有trip的当前状态，-1代表为假。正值代表为真。所显示的正值为该trip最后一次由假变为真的时刻。3 该区域显示所有水力学控制体的当前状态。分3a、3b、3c、3d四个区域。控制体号在左边缘显示。3a中显示的是流体的状态。3b中，rho-mix为两相混合物的密度。mixing cup quality为热力学含气率 static quality为截面含气率 noncondensable quality 为不凝气体份额。3c中显示控制体输运状况和特性。tot. ht.inp是热构件输入到控制体流体的总热流量（w），第二列是热构件输入到蒸汽中的热流量，二者相差一个输入到液体中的热流量。vapor-gen水蒸气的产生速率。水蒸气产生总速率是相间质量交换产生率与壁面沸腾产生率之和。wall-flashing代表壁面沸腾蒸汽产生率。上述二者之差为相间质量交换率。liquid- and vapor-side heat transfer coefficients液相气相传热系数。相间传热假设相间界面是饱和的。接下来控制体质量流量，各相的雷诺数，流型。bby=bubbly（泡状流）, hst=horizontal stratified（水平分层流）, anm=annularmist（环状雾状流）, slg=slug（弹状流）, ian=inverted annular（反向环状流）, isl=inverted slug（反向弹状流）, actb=churn turbulent（搅浑流）。3d 中带有reduce-列表示控制体导致时间步长减小。 4 4a该区域显示水力学接管的状态。接管的面积可能随时间变化（阀接管）。接管真实的面积等于接管面积与喉部比例的乘积。最后三列分别表示最后一个时间步长、最后一次大编辑和从开始计算到现在接管是否出现壅塞现象，0表示没有出现，1表示出现一个壅塞现象。4b开始两列为接管内的液相和气相份额。以“f”开头的表示接管的摩擦压降 ；“fij” 表示相间拖拽。“fwalfj” “fwalgj”表示作用在液相和气相的壁面拖拽。壁面拖拽指的是通过接管相连的两个相邻控制体的中心之间的流体。“fjunf” 和 “fjunr”表示用输入的流动损失系数计算的流动损失。“formfj” 和 “formgj”表示程序计算的流动损失（如面积突然变化引起的损失）。 5 显示热构件的当前状态。5a 实际是临界热流密度可能不会打印出来。如果没有出现沸腾，则临界热流密度不被计算，打印0.00000，例如单相强迫对流。下一列显示使用者指定的临界热流密度乘子。下一列是热传递方式号。int.-heatsource表示内热源。“conv+rad-source”表示热构件的总的热平衡。最后一列表示热构件的体积平均温度。 5b 显示当前热构件的节点温度。第一列表示左边界温度，最后一列表示右边界温度， 6 显示所有控制变量的当前值。控制变量号、变量名、变量类型、变量的当前值。 7 如果这个大编辑也是再启动编辑，则会打印出这个区域。再启动号是103卡中所要求的。block number没有被用到。3.1.4.2 再启动/绘图文件绘图点频率（小编辑频率）应随着问题的进行而改变。当参数变化较快时应采用较大频率，到达静态状态时可使用较小的频率。使用多重重新启动来完成计算。例如一个问题运行时间0-10秒。如果在10秒前也计算结束，接着就要重新运行（比如从10秒到30秒），以此类推。relap5允许在重新启动点改动模型，改动后的模型从该重新启动点开始往后一直生效，除非在此后的另一个再启动点再次改动模型。3.2 计算节点和网孔的大小 以下介绍如何选择水力学元件的节点大小和热构件的网孔大小。 从经济角度讲，水力学元件数和热构件数应尽可能少，因为计算机的运行时间几乎完全取决于模型的水力学元件数。热构件数也随水力学元件数的增加而增加。因此，减少水力学元件数对计算的经济性非常有利。此外，减少热构件内部的网孔数也有利于提高计算的经济性。但减少模型的trip数和控制变量数对模型的经济性贡献非常小。 较小的水力学元件尺寸会减小时间步长，因为时间步长受库兰特限制的约束。 但减小模型的元件数目（即增大元件的尺寸）要考虑到一些重要的现象不被遗漏。例如，大破口和小破口的现象有很大不同，因此，两种模型也有很大不同。3.2.1 水力学元件的大小 通常在模型的产生重要现象的区域（如反应堆压力容器、稳压器、蒸汽发生器）采用较小的节点尺寸。在一些不太重要的区域（如主蒸汽管线、给水序列）采用较大的节点尺寸。 考虑到协调好小元件和大元件之间的过渡，相邻的控制体的尺寸变化不应超过某一限制。3.2.2 热构件网孔的大小由于燃料元件气隙的比热较小，为了避免计算上的困难，通常在气隙内部不设置节点。3.3 通用选项的选择 所谓通用选项，是指对所用部件来说都是相同的选项。主要有：控制体选项、接管选项、初始条件和边界条件。3.3.1 控制体相关选项 控制体相关选项通过控制体控制标识tlpvbfe来设定。一般推荐的输入是省缺形式，即0000000。 t标识指出热分层模型是否使用。热分层模型只在竖直容器且内部不同层的流体温度相差分明（热流体在上冷流体在下）时才使用。 l 标识指出level模型是否使用。level模型只用在level相差分明的竖直管道或容器中（蒸汽在上水在下）。 p 标识指出是否使用水程序包，一般情况下使用。只有在因重复的水程序包造成计算困难时才不用。对tmdpvol, separatr, jetmixer,turbine, pump, and accum 部件，此标识不用，且必须输入0。 v 标识指出是否应用竖直分层模型。一般使用。只有在不需要竖直分层行为时才不使用。对tmdpvol,separatr, jetmixer, eccmix, turbine, pump, and accum 部件,此标志不用，且必须输入0。b 标识指出相间摩擦模型。一般使用管内相间摩擦模型。在有棒束结构的部分（如蒸汽发生器二次侧、堆芯）使用棒束相间摩察模型。separatr, jetmixer, eccmix, turbine, pump, and accum不使用此标志，必须输入0。f 标识指出是否计算壁面摩察。一般要计算。只有在一些不需要考虑壁面摩察的场合才不计算。对separatr and pump部件，不管使用者设定值是多少，此标识自动置为1（不计算壁面摩察）。e 标识指定是否使用相平衡选项。不平衡指液相和气相温度不相等。平衡指两相温度相等。通常选定为不平衡态。对accum 部件，此标识必须为0。3.3.2 接管相关选项 接管相关选项在接管控制标识efvcahs中设定，但时间相关接管没有此控制标识。 e 标识指定是否使用能量修正选项。0不使用，1使用。1适用于以下场合：大的扩张段；核岛内同时包含主回路和一些低压系统（如反应堆安全壳系统）的模型。 f 标识指定是否适用反向流动限制模型（ccfl，countercurrent flow limiting）。0不使用，1使用。一般用0。在希望出现ccfl现象的场合用1。对于和以下部件separatr, jetmixer, eccmix, turbine, and accum相接的接管，此标识不用，且需输入0。 v 标识指定是否使用horizontal stratification vapor pullthrough/liquid entrainment 模型。 c 标识指定是否使用壅塞（临界流）模型。0用，1不用。除非发生壅塞现象，一般应设定为1。 a 指定面积变化情况。0面积平滑变化，1面积突然变化。对于马达阀和伺服阀，可以是0或1。但0情况下，必须输入阀的cv表，如果没有输入cv表，则必须是1。对于其他阀，必须输入1。程序能自动计算面积突变造成的流动损失，这些损失可加到任何其他损失（如弯管损失）上。 h 标记指定用于接管的相速度假设。0表示非均相，液相和气相速度不相等。一般推荐设定0。2表示均相，两相速度相等。 f 标识设定动量流量选项。此标识可用来关掉横流的动量流量。以前的横流模型都是自动关闭动量通量的。3.3.3 初始条件选项 使用者需要为水力学控制体内的流体状态、水力学接管内的流动状态、热构件节点温度、控制变量状态设定初始条件。可选择是否为trip设定初始条件。初始条件不必设定的非常详细，有事过于详细反而不好。建好模型后，建议只给定一些粗略的初始条件，通过程序本身的计算来得到稳态初始条件。例如，所有流体的初始温度都可以设置为冷段温度；初始压力设置为冷段压力；所有初始流速可设为0。所有的热构件和控制变量可以自我初始化。3.3.3.1 控制体流体状态初始化 用流体状态控制字“ bt.”来指定流体的初始状态。指定流体，b指定是否含硼，t指定以下输入的2到5个字是如何被程序解释的。 建议用省缺形式（汽/水）或120-129流体系统控制卡来设定流体，而不用来对所有元件单独指定。因此令=0，或直接省略。 t=0-3指定单组份，4-6指定双组份。尽量避免使用0和6，因为这两种情况下的输入很复杂。建议一般使用2（平衡态压力和含汽率）和3（平衡态压力和温度）。2应用于蒸汽区域（如蒸汽管线和蒸汽腔室）和两相区域（如稳压器相间界面、锅炉、沸水堆堆芯），其实也就是饱和区。3应用于过冷区（如冷管段、热管段、压水堆堆芯）。要包含硼或不凝气体，没必要对系统的所有元件都设定（b=1表示含硼，t=4,5,6指定不凝气体）。当元件的硼浓度或不凝气体含量不为零才需要设定，在系统的任何一个元件设定后，自动扩散到系统的所有元件。3.3.3.2 接管流动初始化 用控制字（0为流速，1为质量流速）来对接管流动进行初始化。需要注意接管初始条件的设定要和该接管上游元件的流体状态和上游接管的初始条件一致。3.3.3.3 热构件初始化 用稳态初始化标识（热构件的第一个输入卡的第4个字）来设定热构件的初始状态。0表示需要输入节点的初始温度，1表示程序根据与热构件相邻的控制体的流体状态和稳态热传递计算初始节点温度。但即使是1，使用者也要输入各个节点的初始温度，这些温度作为稳态温度的起点。对一个新的模型建议用1。对于再启动问题，用0。3.3.3.4 控制变量初始化 用初始值字（控制变量第一个输入卡的第5个字）设定控制变量的初始状态。0表示用第4个字的值作为控制变量的初始值。1表示计算得到初始值。 注意与水力学参数有关的控制变量都用国际单位，尽管开始控制卡上输入的是英制单位。另外，注意控制变量的求解顺序是在最后（在水力学部件、热构件、trip之后）且按数字顺序。3.3.3.5 trip初始化 trip的初始化是选择性的。用timeof（变量trip语句的第8个子，逻辑trip语句的第5个字）进行初始化。-1表示假，一个非负的数表示真。对于大多数new问题和restart问题，建议输入0。3.3.4 边界条件选项3.3.4.1 流体状态边界条件 流体状态边界条件（压力、温度、含气率或空泡份额）在时间相关控制体中设定。这个名字是不准确的，因为最初流体的状态只能设定为时间的函数，而现在的时间相关控制体能将流体的状态设定为任何参量的函数。 时间相关控制体可以指定模型的某一点处的流体状态。应该将时间相关控制体当做一个无限大的流体源或汇聚流体的池子，它的状态可以根据要求而不变或改变，但流入和流出的状态是有差异的。3.3.4.2 流体流动边界条件 流体流动边界条件（流速、质量流量）在时间相关接管中设定。这个名字同样是不准确的，原因同上。时间相关接管的设定要考虑上游的流体的状态。当上游的状态由时间相关控制体设定并保持恒定时，时间相关接管可以根据需要设定流速或质量流量边界条件。但如果上游流体的状态随过程的进行而改变时，就要注意一些问题。例如，上游流体状态由时间相关控制体控制，并且有变化时，时间相关接管设定的参数（如流速）可以保持不变，但没设定的参数（如质量流量）会随上游流体状态的改变而改变，因为流体的密度在改变。 不建议将时间相关接管用于主冷却剂系统的某一出口或是冷却剂系统内部。因为上游流体的状态可能发生急剧变化而给计算带来困难。3.3.4.3 热构件边界条件3.4 特殊模型的应用3.4.1 破口模型 对于上面的小破口，单一接管950应选择面积突变选项，引用壅塞模型，接管初始速度为0。时间相关控制体960模拟安全壳。960的压力可以是定值。或者是一个变化值（如是泄漏率的积分的函数）。对于下面的双端剪切断裂，连接了两个时间相关控制体，因为一个时间相关控制体只能连一个接管。同小破口一样，要用面积突变和壅塞选项。950的初始质量流量为正，大小与进入110的质量流量相同。960的初始质量流量为负，大小与950相同。还可以用一个触发阀来模拟破口。小破口时安全壳内的状况无关紧要，但大破口时就要指定安全壳内的状况，因为此时破口处可能发生从安全壳到系统内的反向流动。可以通过安全壳内的质量平衡和热平衡来模拟安全壳内的状况。relap5对临界流的计算有很大的不确定性，因此有些时候要对结果进行评价。或者通过实验，或者通过临界流标准关系式。3.4.2 硼模型 使用硼模型要令b=1，并在系统的一个或几个控制体内设定初始硼浓度。硼模型具有跟踪整个系统内个部分硼浓度的能力。这是为了计算堆芯硼浓度以得到反应性反馈效应。硼模型用于使用反应堆动力学模型计算堆芯功率，并且堆芯内硼浓度不断变化的情况。relap5假定硼只存在于液体中，并且随液体迁移至系统各部分。硼浓度作为反应性的一部分，可通过两种途径用于反应堆动力学堆芯功率的计算：（1）表4反应堆动力学选项（2）通过表格或控制变量引入一个单独的反应性。前者生成一个四维表格来描述与之相关的反应性效应，即液体密度、液体温度、燃料温度、硼。后者认为硼是一个独立的反应性部件。3.4.3 不凝气体模型 不凝气体的设定在110卡中指定气体类型，并在一个或多个控制体初始条件控制卡中设定初始不凝气体份额。多种不凝气体混合物的设定在110卡中指定多种气体类型，并在115卡中设置各气体份额。一个问题中只能用一种不凝气体（或一种混合物），且如果使用安注箱部件，不凝气体必须是氮气（或含氮气的混合物）。 不凝气体模型假定气体被蒸汽携带。并且气体和蒸汽达到热平衡状态（即二者之间没有热交换）。计算的压力为气/汽混合物的总压力，蒸汽分压力可作为一个标准输出变量而获得。 不凝气体模型的使用受到一些限制。经验表明在“bt”中令t=6设定不凝气体控制体初始条件时会遇到一些麻烦，因此建议t=4。进一步经验表明，程序在气汽混合物产生或消失时（即很小的空泡份额）可能会遇到一些计算问题。解决方法是手动（例如在再启动时）加入或消除水蒸气。3.4.4 再淹没模型 通过在燃料元件通用热构件卡（1cccg000）中设定非零的再淹没条件标识来使用再淹没模型。在以下情况可以设定再淹没初始化条件：（1）相邻控制体是空的；（2）热构件烧干；（3）使用者设定的trip。 再淹没只出现在大破口事故中。由于relap5主要是一个小破口事故分析工具，所以目前关于再淹没的结果评估很少，能够给予的建议也很少。一些实验表明relap5的时间步长控制特性可能不足以解决再淹没问题。3.4.5 横流接管模型relap5是一个一维的编码，所有控制体只有一个入口和一个出口。横流接管模型用于模拟一些多维的现象。正常的接管只能连在元件的表面，而横流接管模型能连在元件的中心。当横流接管用于连接竖直部件和水平部件时，建议竖直部件的高度与水平部件的直径相一致（但高度不要小于1英尺（ft）。使用横流接管模型要用扩大连接编码。其格式为：cccvv000n。ccc为部件号，vv为控制体号，n为表面号。扩大连接编码假定每个控制体有六个表面，即三个坐标轴方向都可以有进口和出口。单一控制体时vv为01，管型中vv从01开始直到最后一个内部控制体号为止。为1和2分别为正常x轴方向入口和出口表面；n为3和4表示y轴方向入口和出口表面；n为5和6表示z轴方向入口和出口表面。因此当n为1,2是为正常接管连接，n为3到6表明为横流接管连接。尽管横流模型可以模拟一些多维流动现象，但横流模型并不是一个完全的多维模型。因此，在横流接管内流速与正常接管内流速相比大小相当甚至更大时，建议不要使用横流接管模型。稳压器波动管与热管段的连接是典型的应用横流接管的例子。 横流接管不能与管型内部的控制体单元相连，而只能与一个单独的控制体相连。3.4.6 反向流动限制（ccfl）模型 ccfl模型是relap5 新引入的模型，以前的版本没有。3.4.7 控制系统模型 控制变量可用于：进行有用的副计算，指定复杂的边界条件，模拟实际的控制系统。3.4.7.1 有用的副计算 例子：追踪蒸汽发生器二次侧总质量，积分注入流量，积分破口流量，计算蒸汽发生器总的换热速率。3.4.7.2 指定复杂的边界条件 例子：模拟上充与下泄系统。3.4.7.3 模拟实际控制系统 例子：模拟稳压器液位控制系统3.4.8 水位跟踪模型3.4.9 热分层模型 4 具体规则4.1 问题控制选择4.1.1 格式考虑第一个卡为问题的标题，以“=”开头。建议为输入卡加上注释。如果一个数据分布在两个卡上，则在第二个卡的开头加一个“+”，第二个卡不需要卡号。一个数据区域必须在一个卡内完成。建议（但不是必须）各卡以卡号从小到大顺序排列。以句号“.”或前斜线“/”开头的卡为结束卡。建议加上注释，如 . *end of input。4.1.2 问题类型100卡 relap5有三种计算类型：new、restart、strip。 new类型要求一个完整的数据输入流，所有的特性都要设定初始条件。 restart类型是在上次计算的基础上继续计算。可以是简单地延续上次的计算，也可以在再启动点对模型进行一些改动后再次计算。 strip类型是对输出数据的后处理过程。可以选择一些关心的参数绘图输出，其他的参数不输出，减少数据数量。对new类型有transnt（瞬态） 和 stdy-st（稳态）两个选项。但建议不要用stdy-st（稳态）选项。4.1.3 输入检查/运行101卡该卡是可选择的，如果省略则默认为运行（run）模式。建议用运行（run）模式。4.1.4 单位选择102卡 该卡是可选择的，如果省略则默认为si单位制。不管用户输入的是什么单位制，relap5内部计算和数据在输出文件的储存都使用si单位。所以为了避免一些错误，建议使用si单位。4.1.5 再启动控制103、104卡 103卡指定再启动号，这个号在上次计算的打印输出文件中再启动编辑中显示。如果需要，104卡可以用于阻止再启动/绘图文件的写入。再启动输入板的输入非常简短。只包括问题控制卡、时间步长控制卡和指定改变信息的卡。但注意，即使只改变某个部件的一个变量，与该部件相关的所有输入卡都要在再启动输入板中重新输入一遍。再启动输入时，用户要十分小心地指定初始条件。这在再启动点的大编辑数据中显示（因此建议大编辑频率与再启动编辑频率应相等）。4.1.6 cpu时间控制105卡 105卡可以在内部结束一个计算，这是通过达到cpu使用时间限制来实现的。建议使用这个卡来正常结束计算。如果通过外部来结束一个计算（比如，花费的计算时间达到了在外部工作卡上设定的最大值），一些输出数据可能会丢失。 字2应比字1大，字3是分配给cpu的时间。 cpu所需时间的估算方法：线性估计，与在同一计算机上运行的相似问题比较，（1）与问题中水力学单元的个数成正比；（2）与问题时间成正比；（3）与时间步长成反比。4.1.7 不凝气体类型110和115卡4.1.8 水力学系统定义12x卡 该卡用于指定用于水力学的流体类型：轻水、重水、氢气。该卡是选择性的，如果没有则默认为轻水。4.1.9 自我初始化选项140-147卡 这些卡是选择性的，不需要的。4.2 时间步长控制200-299卡4.3 小编辑和扩大编辑/绘图变量要求4.3.1 小编辑要求301-399卡对大多数数据，字1代表数据的类型，字2代表数据的位置。例如，管型部件125的第3个单元的压力可表示为：p 125030000。两点注意（1）部件量（component quantities）的数据位置是三位部件号ccc，而不是九位单元号cccnn0000。如，泵255的转速为：pmpvel255，而不是像上面例子那样：pmpvel225010000。（2）接管量的数据位置，如果只有一个接管与部件ccc相连(single-junction, valve, time-dependentjunction, and accumulator components include only one junction)，接管位置表示为ccc000000。如果有多个接管与部件ccc相连pipe/annulus, branch/separator/jet- mixer/turbine/emergency core cooling (ecc) mixer, pump, and multiple-junction components may include more than one junction，则接管位置为cccnn0000，nn为该部件内部的接管号。建议在一个新的输入板（模型）中包含一些小编辑要求来描述整个模型的行为。小编辑要求可以包含模型的重要部分的有代表性的参量如压力、温度、流速、速度，还应该包含一些重要控制变量的当前值。对于一个已经比较成熟的输入板（模型），小编辑要求只需包含那些表明瞬态过程的量，如堆芯功率、破口流动、燃料温度等。4.3.2 扩大编辑/绘图变量要求4.4 trips用户要注意trip值是在每一个时间步长，按照trip号由小到大的顺序被决定的。如果引用的是较小号的trip，则它的值为当前时间步长的值，如果引用的是较大号的trip，它的值为上一个时间步长的值。4.4.1 变量trips 401-599卡可以将两个参数进行比较，也可将讲一个参数和一个常量进行比较。例子：505 p 140010000 gt p 145010000 50. n -1. *delta p是否14001的压力比14501的压力大50 psia？n代表不锁定，即在每一个时间步长都问此问题（检查505的值）。于此不同，l代表锁定，每个时间步长都检查505的状态，直至它的值为真，此后不再检查，保持为真，不再变化。-1表示该trip的初始值为假，如果输入0或一个正数，表示初始值为真，如果省略，默认初始值为假。对泵部件，如果开始时没有用trip，后来想在再启动时加上trip，则会非常麻烦。因为仅仅为了加一个trip，却需要在再启动时把泵的全部输入卡（可能上百个）再次输入一遍。为此，可在第一次输入泵部件时就加一个dummy trip，它的值恒为真或恒为假，这样，再启动输入时只需改变trip的内容就行了。例如：507 time 0 lt null 0 1.e6 l 0. *always true508 time 0 gt null 0 1.e6 n -1. *always false507的初始值为真，且其后一直为真。508的初始值为假，且其后一直为假。（这里假设问题时间不可能超过1000000秒）4.4.2 逻辑trips 601-799卡例子：507 time 0 gt timeof 605 0.5 l -1.当前的时间比605最后一次改变的时间大0.5秒？4.4.3 用trip结束一个计算600卡501 mflowj 150000000 lt null 0 100. l -1600 501当接管150的流量小于100时，结束计算。4.5 交互式变量801-999卡当用户用npa（nuclear plant analyzer）软件作界面运行relap5时才用得到交互式变量。4.6 水力学部件 anticipated transient without scram (atws)ccc0000形式的卡说明水力部件的类型，建议将它放在每个部件的最前面。4.6.1 单一控制体部件单一控制体是relap5中最基本的水力学单元。管型可看作是一系列由单一接管连接的单一控制体组成。分支可看作可以连接一个或多个单一接管的单一控制体。描述其他部件的基本控制体的几何条件和初始条件的输入数据与单一控制体的输入数据相同。为了避免不匹配情况的发生，建议流动面积、长度、容积三者中有一个量输入0。方位角是为了三维建模而设置的，目前还没有用，建议输入0。控制体高度变化的绝对值与控制体长度的绝对值的比值算出的角度用于选择流型图。当该角大于45度时用竖直管流型图，小于45度时用水平管流型图。控制体的高度变化必须小于等于控制体长度sin（倾角）。对于绝对粗糙度，建议commercial steel finishes用0.0000457 m (0.00015 ft)，drawn tubing