Mingrui YANG, Qizheng SUN, Chixu LUO, Donghao HE, Xiaojing LIU, and Tengfei ZHANG*
Fig. 1. Flowchart of the transient calculations with SCM in MORPHY
Fig. 2. Channel equivalent diagram
Fig. 3. The coupling method in MORPHY (a) OSSI method, (b) FPI method
Fig. 4. Layout of the TWIGL benchmark problem
Fig. 5. Results of the TWIGL A1 (a), A2 (b) problem
Fig. 6. Layout of the Dodds problem (a) Material layout, (b) Radial mesh generation
Fig. 7. Dodds benchmark problem relative power vs. time
Fig. 8. Layout of the NEACRP benchmark problem
Fig. 9. NEACRP benchmark problem relative power vs. time (a) NEACRP A1, (b) NEACRP A2, (c) NEACRP B1, (d) NEACRP B2
Fig. 10. NEACRP benchmark power distribution at peak power
时间 Time / s | DAISY | MORPHY Δt=5 ms | MORPHY Δt=20 ms | 误差Error / % |
---|
vs. Δt =5 ms | vs. Δt =20 ms |
---|
0.0 | 1.000 | 1.000 | 1.000 | 0.00 | 0.00 | 0.1 | 1.318 | 1.313 | 1.313 | -0.38 | -0.38 | 0.2 | 1.985 | 1.982 | 1.982 | -0.15 | -0.15 | 0.3 | 2.103 | 2.103 | 2.103 | 0.00 | 0.00 | 0.4 | 2.121 | 2.121 | 2.120 | 0.00 | -0.05 | 0.5 | 2.139 | 2.140 | 2.139 | 0.05 | 0.00 |
|
Table 1. Comparison of core relative powers for the TWIGL A1 problem
时间 Time / s | DAISY | MORPHY Δt=5 ms | MORPHY Δt=20 ms | 误差Erorr / % |
---|
vs. Δt= 5 ms | vs. Δt= 20 ms |
---|
0.0 | 1.000 | 1.000 | 1.000 | 0.00 | 0.00 | 0.1 | 2.088 | 2.089 | 2.087 | 0.05 | -0.05 | 0.2 | 2.106 | 2.107 | 2.106 | 0.05 | 0.00 | 0.3 | 2.124 | 2.125 | 2.123 | 0.05 | -0.05 | 0.4 | 2.143 | 2.143 | 2.141 | 0.00 | -0.09 | 0.5 | 2.161 | 2.161 | 2.160 | 0.00 | -0.05 |
|
Table 2. Comparison of core relative powers for the TWIGL A2 problem
程序Code | 临界硼浓度Critical boron concentration / 10-6 | 控制棒价值Control rod value / 10-5 |
---|
A1 | A2 | B1 | B2 | A1 | A2 | B1 | B2 |
---|
PARCS(Reference) | 561.26 | 1 158.85 | 1 248.21 | 1 185.55 | 827.70 | 89.30 | 829.70 | 98.40 | PANTHER(1993) | 567.70 | 1 160.60 | 1 254.60 | 1 189.40 | 821.80 | 89.50 | 831.30 | 99.10 | PANTHER(1997) | 561.20 | 1 156.63 | 1 247.98 | 1 183.83 | 821.80 | 89.50 | 831.00 | 99.10 | DAISY | 563.28 | 1 157.06 | 1 247.12 | 1 185.52 | 820.99 | 90.88 | 826.19 | 93.25 | MORPHY(S2) | 566.62 | 1 159.23 | 1 253.51 | 1 187.72 | 819.74 | 88.39 | 828.25 | 93.79 | 绝对误差Absolute error | 5.36 | 0.38 | 5.30 | 2.17 | 7.96 | 0.91 | 1.45 | 4.61 |
|
Table 3. Comparison of steady-state results for NEACRP problems
工况 Case | 时间步长 Time-step sizes / ms | 峰值功率 Power peak / % |
---|
A1 | 2.5 | 149.82 | A1 | 5 | 155.99 | A1 | 10 | 169.99 | A2 | 2.5 | 108.19 | A2 | 5 | 108.20 | A2 | 10 | 108.23 |
|
Table 4. Peak power at different time-step sizes for OSSI method
ΔtN / ms | ΔtTh / ms | 峰值功率 Power peak / % | 绝对误差 Absolute error |
---|
10 | 1 | 141.74 | — | 20 | 1 | 140.10 | -0.64 | 20 | 20 | 135.55 | -6.19 | 50 | 1 | 129.02 | -12.72 | 50 | 50 | 119.99 | -21.75 |
|
Table 5. Peak power of different thermal-hydraulic time-step sizes using FPI method in case A1
程序 Code | A1 | A2 |
---|
PMax / % | tMax / % | P5 / % | TM / ℃ | TC / ℃ | PMax / % | tMax / % | P5 / % | TM / ℃ | TC / ℃ |
---|
PARCS | 126.19 | 0.54 | 19.90 | 293.38 | 686.77 | 108.14 | 0.09 | 103.55 | 325.03 | 1 703.36 | MORPHY(S2) | 141.74 | 0.54 | 21.93 | 293.45 | 716.35 | 108.17 | 0.10 | 103.54 | 325.85 | 1 716.80 | 绝对误差Absolute error | 15.55 | 0 | 2.03 | 0.07 | 29.91 | 0.03 | 0.01 | 0.01 | 0.82 | 13.44 |
|
Table 6. Comparison of MORPHY A1 and A2 results with PARCS