1. 目录

[TOC]

2. 编写输入卡

output/input.txt 为输入卡，格式比relap5精简。

2.1. 标题-时间100-201

= simple pipe flow

******time card*****
*迭代输出100：GMRES内迭代残差输出时间; 牛顿外迭代残差输出时间
*时间离散(0一阶精度1二阶); 空间离散=0：迎风差分；=1：中心差分；=2：二阶迎风；*=3：quick；=4：三阶迎风；=5：fromm；=6：minmod；=7：muscl；=8：koren；=9：smart；=10：superbee；=11：harmonic；=12：ospre；=13：van albada
101   0   0

200   1 *时间设置组数 
201   0.0   10.0   0.01 *开始时间; 结束时间; 步长201;

2.2. 水力部件CCCXXXX(7位)

2.2.1. tmdpvol

tmdpvol *不写卡号和部件名称，直接写类型

1000100  1.13097334e-4  0.1  0.0 *面积; 长度; 倾角

*物性参数组数; 控制体参数标记（0p,Ul,Ug,α， 1T,x0， 2p,x0， 3p,T, 其中x0平衡静态含气率)
1000200  2  0
*物性参数组1、2：时间;控制体参数标记(以此为准); 压力; 液相内能(J/kg); 气相内能; 空泡份额
1000201  0.0  0  6.89e6   0.8061e6   2.58198e6  0.0
1000202  100.0  0  6.89e6   0.8061e6   2.58198e6  0.0

2.2.2. tmdpjun

tmdpjun

*上游控制体（不区分面，100组件的01控制体）; 下游控制体; 面积
1030101  10001  10501  1.13097334e-4

1030200  2  0 *物性组数; 接管参数标记（0液、汽流速，1液、汽流量）
1030201  0.0  0  1.1204  1.1204 *组1：时间; 接管参数标记(以此为准); 液、汽流速
1030202  100.0  0  1.1204  1.1204

2.2.3. sngljun

sngljun

1100101  10520  11501  1.13097334e-4 *上游控制体; 下游; 面积
*接管参数标记（0液、汽流速，1液、汽流量）; 液相速度; 汽相速度; 4正向阻力系数; 
        *5反向阻力系数; 过冷喷放系数; 两相喷放系数; 接管临界流、阻力系数模型选项ab
1100201  0  1.1204  1.1204  0.0  0.0  0.0  0.0  0
*接管选项ab:*a：是否启用临界流，=0：不启用；=1：启用
*b：接管阻力计算选项，=0：不计算阻力；=1：只计算突缩、突扩的形阻系数；
*=2：只计算输入卡给定的正反向阻力系数造成的局部阻力系数；=3：1、2都计算

2.2.4. pipe

pipe
1050000  20 *控制体数

1050100  1 *控制体信息组数
*组1：面积; 长度; 倾角; （水位模型标记？壁面粗糙度？）; 控制体编号
1050101  1.13097334e-4  0.075  90.0  0  20

1050200  2 *控制体物性参数组数
*组1/2：参数标记(0p,Ul,Ug,α， 1T,α， 2p,α， 3p,T); 
        *压力; 液相内能(J/kg); 气相内能; 空泡份额; 控制体编号
1050201  0  6.89e6   0.8061e6   2.58198e6  0.0  10
1050202  0  6.89e6   0.8061e6   2.58198e6  0.0  20

1050300  1 *接管物性参数组数
*组1：接管参数标记（0液/汽流速，1液/汽流量）; 液相速度/流量; 汽相速度/流量; 正向阻力系数; 5反向阻力系数; 
        *6过冷喷放系数; 两相喷放系数; 管临界流、阻力系数模型选项ab; 接管编号
1050301  0  1.1204  1.1204  0.0  0.0  1.0  1.0  0  19

2.3. 热构件(1CCCGXNN, 8位)

hstructure

*轴向控制体数量; 径向网格点数; 几何类型(1矩形体2圆柱3球形);
*左边界坐标(管内径); 再淹没标记1(0禁用模型1启用); 再淹没标记2
11050000   20   5   2   0.006   0   0

11050100   1 *网格格式标记
11050101   4   0.008 *网格数(网格点数-1);右边界坐标(管外径)

11050200   1 *径向材料成分组数？
11050201   5   4 *材料编号(201MMM00卡中的MMM，此处为20100500); 网格编号

11050300   1 *径向功率分布组数？
11050301   1.0   4 *径向功率因子; 网格编号

11050401   464.0 *初始温度

11050500   1 *左边界条件组数？
*边界条件类型(0绝热后两参数无，1定壁温206表温度-时间，2定热流203表热流-时间，3对流); 边界相邻水力控制体/通用表编号; 0; 周向网格长度; 热构件编号
11050501   3  10501  0  0.075   20
*类型W1=3：如果w2=300~399，换热系数-时间表（204表）；w3指定冷源温度，w3小于999，温度-时间表编号；w3大于10000，冷源是控制体，w3=控制体编号
       *如果w2=400~499，换热系数-温度表（205表），w3指定冷源温度，w3小于999，温度-时间表编号；w3大于10000，冷源是控制体，w3=控制体编号
       *如果w2=其它，w2=控制体编号，w3=0

11050600   1 *右边界条件组数？
11050601   0  0  0  0.075   20 *同上

11050700   1 *热源分布组数？
11050701   999   1.0   20 * 热源类型(999查总表202999XX); 热源乘子; 热构件编号

2.4. 总表table(CCCXXXNN, 8位)

• 包含了relap5中的热物性表和总表。
• CCC表示表的类型: 201 材料物性表（1-49导热系数，51-99体积热容），202 功率时间表，203 热流时间表，204 换热系数时间表，205 换热系数温度表，206 温度时间表。
• XXX为查表时所需编号。

*201XXXNN材料物性表，NN分两部分：01-49导热系数，51-99体积热容
table

20100500 4 *导热系数的组数

20100501 3. 14.88 *温度; 导热系数
20100502 373. 15.05
20100503 423. 15.56 
20100504 473. 16.09

20100550 4 *体积热容的组数

20100551 3. 3.81e6 *温度; 体积热容
20100552 373. 3.83e6                      
20100553 423. 3.91e6                      
20100554 473. 3.99e6

*202XXXNN功率时间表
table

20299900   2 *功率组数
20299901   0.0   3195.0 *时间; 功率
20299902   1000.0   3195.0

2.5. 结束

.end. *核对某行前五个字符

3. 输出

参考文件NUSOL_FI_JFNK.f90

100——input.txt
200——output.txt
101——velocity_gas.txt
102——velocity_liquid.txt
103——internal_energy_gas.txt
104——internal_energy_liq.txt
105——void_fraction.txt
106——pressure.txt
107——tempg.txt
108——tempf.txt
109——JFNK.txt
110——wucha.txt
111——satt.txt
112——sathg.txt
113——sathf.txt
150——1111.txt
114——fi.txt
115——quale.txt
116——mflow.txt
117——htemp.txt
218——218.txt

4. 杂项

1. 输入、输出均在output文件夹
2. 103 format(A200)，定义输出格式

5. 程序文件表

6. 变量表