Fundamental Approach to TRIGA SteadyState Thermal Hydraulic CHF

  • Slides: 31
Download presentation
Fundamental Approach to TRIGA Steady-State Thermal. Hydraulic CHF Analysis National Organization of Test, Research,

Fundamental Approach to TRIGA Steady-State Thermal. Hydraulic CHF Analysis National Organization of Test, Research, and Training Reactors (TRTR) Meeting Lincoln City, Oregon September 17 -20, 2007 Earl E. Feldman

Outline n Two-Step Process (Step 1: Flow; Step 2: CHF) n Flow – Coolant

Outline n Two-Step Process (Step 1: Flow; Step 2: CHF) n Flow – Coolant channel geometry of models – Computer codes (STAT & RELAP 5) – Nodal structure of RELAP 5 models used to determine flow – List representative parameters for two generic TRIGA reactors -- a hexagonal pitch TRIGA and a rectangular pitch TRIGA – Compare STAT and RELAP 5 flow results for a representative hexagonal pitch TRIGA reactor 2

Outline (continued) n Critical Heat Flux (CHF) – Bernath correlation – Groeneveld tables (1986,

Outline (continued) n Critical Heat Flux (CHF) – Bernath correlation – Groeneveld tables (1986, 1995, 2006) – Hall and Mudawar (Purdue) outlet correlation – PG-CHF (Czech Republic) correlations – Compare CHF correlations for representative TRIGA reactor conditions – Compare CHF power predictions for a representative hexagonal pitch TRIGA reactor n Suggested Approach to CHF n Conclusions 3

Geometric Model for Calculation of Coolant Flow Rates (Step 1) n The core flow

Geometric Model for Calculation of Coolant Flow Rates (Step 1) n The core flow area is divided into subchannels defined by the cusps between adjacent fuel rods. n Assume no mass exchange or heat transfer between adjacent subchannels, i. e, each subchannel behaves independently of its neighbors and can be analyzed separately. n Only potentially limiting subchannels need be considered. n Divide the length of the subchannel being analyzed into a series of horizontal layers or nodes. The 15 -inch (0. 381 -m) heated length was divided into 15 1 -inch layers. 4

Codes Being Used for Thermal-Hydraulic Analysis n STAT – GA-developed code with fixed geometry

Codes Being Used for Thermal-Hydraulic Analysis n STAT – GA-developed code with fixed geometry of one subchannel. – Custom made for TRIGA reactor hydraulics. – Steady state only. – No fuel rod temperature model – Has 2 CHF correlations • Bernath (1960) • Mc. Adams (1949) n RELAP 5 -3 D (Version 2. 3) – Current developer is the Idaho National Laboratory – General transient thermal-hydraulic neutronics reactor code. No fixed geometry. Uses a series of coolant nodes and junctions. Heat structures attached to coolant nodes represent solid regions, such as fuel rods. – Has 2 CHF correlation options • 1986 Groeneveld table • PG-CHF from the Czech Republic (~1994) 5

RELAP 5 Thermal-Hydraulic Model for Current Analysis 6

RELAP 5 Thermal-Hydraulic Model for Current Analysis 6

Representative TRIGA Generic Reactor Parameters (Not the Most Limiting Values for Safety Analysis) Reactor

Representative TRIGA Generic Reactor Parameters (Not the Most Limiting Values for Safety Analysis) Reactor Parameter Hexagonal Pitch Rectangular Pitch Hexagonal Rectangular conversion Flow area per rod, cm 2 5. 464 5. 532 Hydraulic diameter, mm 18. 64 19. 65 Rod (heated) diameter, mm 37. 34 35. 84 25 (77) 30 (86) 1. 68 1. 80 114. 8 (238. 6) 116. 9 (242. 4) Inlet K-loss 3. 58 1. 672 Exit K-loss 3. 0 0. 6 Reactor power, MW 2. 0 1. 0 Number of rods 100 90 Radial power factor (hot. rod) 1. 565 Power of hottest rod, k. W 30. 0 17. 4 Fuel element pitch Inlet temperature, C (F) Pressure (~mid-core), bars Saturation temperature, C (F) 7

Axial Power Shape for Hottest Rod of Hexagonal Pitch TRIGA 8

Axial Power Shape for Hottest Rod of Hexagonal Pitch TRIGA 8

Comparison of STAT and RELAP 5 Results n STAT void detachment fraction is assumed

Comparison of STAT and RELAP 5 Results n STAT void detachment fraction is assumed to be zero. n RELAP 5 fails to provide a stable (nonoscillatory) solution above 48 k. W/rod. 9

Representative TRIGA CHF Parameters for the Limiting Channel (Step 2) Nominal Conditions CHF Conditions

Representative TRIGA CHF Parameters for the Limiting Channel (Step 2) Nominal Conditions CHF Conditions Less than boiling Less than or at boiling Mass flux, kg/m 2 -s ~100 ~300 Velocity, cm/s (ft/s) ~10 (~1/3) ~30 (~1) ~1. 8 Mixed-mean coolant temperature Pressure, bar n Difficulty: Much of the published CHF measurements is focused on power reactors, which operate at high pressures and flow rates. However, TRIGA reactors operate at low pressures and at low (natural-convective) flow rates. 10

CHF Correlations Considered n Bernath (1960) – Used in STAT code along with Mc.

CHF Correlations Considered n Bernath (1960) – Used in STAT code along with Mc. Adams (1949) (STAT results indicate that for TRIGA reactors the Bernath correlation predicts lower CHF values than does the Mc. Adams correlation. ) n 1986 Groeneveld Table – RELAP 5 option n 1995 Groeneveld Table – Not available in RELAP 5 n 2006 Groeneveld Table – Not available in RELAP 5 n Hall and Mudawar (Purdue) – Proprietary 1998 collection of world’s CHF data in water. Has a simple correlation for subcooled boiling. for quality < − 0. 05 and G>300 kg/m 2 -s n PG-CHF (Czech Republic, ~1994) – RELAP 5 rod-bundle option, 4 flavors 11

Bernath Correlation (1960) n Based on low pressure subcooled measured data – 1956 Columbia

Bernath Correlation (1960) n Based on low pressure subcooled measured data – 1956 Columbia University data • Annulus formed by 27. 4 -mm (1. 08 -inch) diameter heater inside an unheated tube • 14 tests with approximate ranges of 2 to 4 bar, 80 to 110° C, 1800 to 9000 kg/m 2 -s (6 to 30 ft/s), and dh = 10. 6 to 14. 7 mm – 1949 Mc. Adams data – 0. 25” heater inside 0. 77” tube (dh = 13. 2 mm) n Checked by Bernath against several sets of independently measured data covering a wide range of parameters n Applicable to subcooled boiling; Limited applicability to low-pressure bulk boiling 12

Bernath CHF Correlation n CHF = CHF, pound centigrade units per hr-ft 2 (1

Bernath CHF Correlation n CHF = CHF, pound centigrade units per hr-ft 2 (1 p. c. u. = 1. 8 Btu) n film coefficient at CHF, p. c. u. /hr-ft 2 -C n TWBO = wall temperature at CHF, C n Tb = bulk coolant temperature, C n De = hydraulic diameter, ft n Di = diameter of the heated surface = heat perimeter / π, ft (In STAT code, diameter of fuel rod) n P = pressure, psia n V = coolant velocity, ft/s 13

1986, 1995, and 2006 Groeneveld CHF Look-Up Tables n CHFtable is a function of:

1986, 1995, and 2006 Groeneveld CHF Look-Up Tables n CHFtable is a function of: – pressure (k. Pa) – mass flux (kg/m 2 -s) – quality – Negative values are used to represent subcooled conditions n Based on water flowing inside an 8 mm diameter tube that is heated from the periphery n Linear interpolation used for values between table entries n Multiplicative factors for other geometries and conditions – CHFbundle = CHFtable × K 1 × K 2 × K 3 × K 4 × K 5 × K 6 × K 7 – 1986 has 6 factors. – Factors have changed after 1986. Later ones have 7 factors. – Some of the newer factors are tentative or not well defined – Most factors should be close to 1. 0 14

Groeneveld K 1 and K 2 Factors n 1986 K 1 (hydraulic diameter, dh)

Groeneveld K 1 and K 2 Factors n 1986 K 1 (hydraulic diameter, dh) n After 1986 K 1 (hydraulic diameter, dh) – For dh = 18. 64 mm (hexagonal pitch TRIGA): • 1986 => K 1=0. 79 • After 1986 => K 1=0. 66 – After 1986 / 1986 = 0. 83 n For K 2 (rod bundle factor) – After 1986 a tentative new relationship was suggested. – The 1986 relationship will be assumed to apply to all years. – It is K 2 = min[ 0. 8, 0. 8 × exp(-0. 5 × quality(1/3) ] – Therefore, K 2 = 0. 8 for subcooled regions and less for bulk boiling regions. 15

Groeneveld K 4 Factors n For K 4 (heated length factor) – It appears

Groeneveld K 4 Factors n For K 4 (heated length factor) – It appears that it has not been changed between 1986 and 2006. – The following is based on the RELAP 5 source code: • X = quality • L = heated distance from channel inlet to middle of node • D = heated diameter (i. e. , 4 × flow area / heated perimeter) • ρf and ρg are the densities of saturated liquid and vapor, respectively. • If X < 0, X = 0 • If L/D < 5, L/D = 5 • α = X / (X + ρg (1 − X) / ρf ) • K 4 = exp( D/L × exp( 2 × α ) ) – For X slightly greater than 0, K 4 increases rapidly with quality. This does not seem to affect the near limiting CHF powers for the generic hexagonal pitch TRIGA. 16

Errors Associated with 2006 Groeneveld Table* n For the region of the table of

Errors Associated with 2006 Groeneveld Table* n For the region of the table of interest for TRIGA reactors, the CHF values are not a result of direct measurement. These regions, Groeneveld* states, “represent calculated values based on selected prediction methods …” n In addition, Groeneveld* uses smoothing methods to eliminate discontinuities that are a result of scatter in the measured data. The paper provides RMS errors between the measured data and the smoothed entries in the table. For the direct substitution method being used in the current analysis, negative qualities in the measured regions of the table have an RMS error of 14. 74%. Positive quality regions have much higher RMS errors. * D. C. Groeneveld, J. Q. Shan, A. Z. Vasić, L. K. H. Leung, A. Durmayaz, J. Yang, S. C. Cheng, and A. Tanase, “The 2006 CHF look-up table, ” Nuclear Engineering and Design 237 (2007) 1909 -1922. 17

Hall & Mudawar (Purdue) CHF Outlet Correlation Symbol Variable Minimum Maximum D Hydraulic Diameter,

Hall & Mudawar (Purdue) CHF Outlet Correlation Symbol Variable Minimum Maximum D Hydraulic Diameter, mm 0. 25 15. 0 G Mass Flux, kg/s-m 2 300 30, 000 1 2000 -1. 00 -0. 05 Pressure, bar x 0 Quality hfg Latent heat of vaporization σ Surface tension ρf Density of saturated liquid ρg Density of saturated vapor 18

PG-CHF (Czech Republic) CHF Data n One of 2 CHF options built into RELAP

PG-CHF (Czech Republic) CHF Data n One of 2 CHF options built into RELAP 5. (The other is Groeneveld 1986. ) n Based on three separate experimental databases – one for tubes, one for rod bundles, and one for annuli. n For each geometry there are four PG-CHF forms called: “Basic, ” “Flux, ” “Geometry, ” and “Power” (It appears RELAP 5 produces obviously erroneous results for the “Basic, ” “Flux, ” and “Geometry” forms. ) n Rod bundle database – 153 test geometries – 7, 616 total points n Data ranges for rod bundles: – Pressure: 2. 8 to 187. 3 bar (TRIGA ~1. 8 bar) – Mass flux: 34. 1 to 7478 kg/m 2 -s – Quality: subcooled to 100% steam – Heated length: 0. 4 to 7. 0 m (TRIGA 0. 381 m) – Fuel rod diameter: 5 to 19. 05 mm (TRIGA ~37 mm) 19

CHF vs. Coolant Quality for 8 mm Diameter Tube 1. 8 bar, 300 kg/m

CHF vs. Coolant Quality for 8 mm Diameter Tube 1. 8 bar, 300 kg/m 2 -s 11. 5 C 38. 0 C 64. 4 C 90. 8 C 116. 9 C Coolant Temperature 20

CHF vs. Temperature for 19. 65 mm Diameter Tube 1. 8 bar, 300 kg/m

CHF vs. Temperature for 19. 65 mm Diameter Tube 1. 8 bar, 300 kg/m 2 -s (Rectangular Pitch TRIGA) 21

CHF Ratios for Hexagonal Pitch TRIGA Evaluated at Nominal Power, where Highest Power Rod

CHF Ratios for Hexagonal Pitch TRIGA Evaluated at Nominal Power, where Highest Power Rod is 30 k. W n CHR Ratio = local CHF prediction / local heat flux n Thermal-hydraulics code is shown in parentheses Directly from STAT Directly from RELAP 5 (1 Corresponds to 30 k. W/rod) 22

PG-CHF Ratios for Hexagonal Pitch TRIGA Evaluated at Nominal Power, where Highest Power Rod

PG-CHF Ratios for Hexagonal Pitch TRIGA Evaluated at Nominal Power, where Highest Power Rod is 30 k. W n RELAP 5 flow except for Bernath, which uses STAT flow 23

CHF Power Prediction of Hexagonal Pitch TRIGA Based on Groeneveld 2006 Table 24

CHF Power Prediction of Hexagonal Pitch TRIGA Based on Groeneveld 2006 Table 24

CHF Power Prediction of Hexagonal Pitch TRIGA Based on Bernath (1960) Correlation 25

CHF Power Prediction of Hexagonal Pitch TRIGA Based on Bernath (1960) Correlation 25

CHF Power Prediction of Hexagonal Pitch TRIGA Based on Purdue (Outlet) Correlation n Not

CHF Power Prediction of Hexagonal Pitch TRIGA Based on Purdue (Outlet) Correlation n Not valid because at CHF conditions the mass fluxes, G, is less than 300 kg/s-m 2 and the quality, X, is greater than -0. 05. For a CHF power of 50. 6 k. W, G is 265 kg/s-m 2 and X is -0. 02 at the limiting axial location. 26

CHF Power Prediction of Hexagonal Pitch TRIGA Based on PG-CHF (~1994) Correlations 27

CHF Power Prediction of Hexagonal Pitch TRIGA Based on PG-CHF (~1994) Correlations 27

Summary of CHF Results for Hexagonal Pitch TRIGA B – Flux & Power PG-CHF

Summary of CHF Results for Hexagonal Pitch TRIGA B – Flux & Power PG-CHF A B – Groeneveld 2006 B – RELAP 5/Bernath & B - Purdue A 28

Summary of CHF Results for Hexagonal Pitch TRIGA (continued) Flow STAT CHF Correlation Rod

Summary of CHF Results for Hexagonal Pitch TRIGA (continued) Flow STAT CHF Correlation Rod CHF Power, k. W A** Bernath B+ C++ 37. 1 52. 5 57. 5 Bernath 49. 6 50. 6 Purdue 48. 9 50. 6 Groeneveld 2006 62. 1 68. 9 RELAP 5 Groeneveld 1986 71. 9 CHF Ratio* A** B+ C++ 1. 24 1. 75 1. 69 1. 92 1. 63 1. 69 2. 07 2. 30 100. 3 2. 40 3. 30 PG-CHF, Basic 105. 9 124. 4 3. 53 4. 15 PG-CHF, Geometry 108. 9 129. 7 3. 63 4. 32 PG-CHF, Power or Flux 109. 2 128. 7 3. 64 4. 29 *1. 0 corresponds to 30 k. W for the highest power rod and 2. 0 MW for the reactor. **A (RELAP 5 Flow): CHF curve at maximum calculated flow per rod (0. 1394 kg/s, thin vertical black line A-A in the previous figure), where RELAP 5 flow begins to oscillate. +B (Extrapolated RELAP 5 Flow): Intersection of a CHF correlation curve and a reactor flow curve, as shown on the previous figure. ++C (Not Recommended): CHF based on calculated reactor power and flow at 30 k. W/rod. 29

Suggested Approach to CHF n Use the 2006 Groeneveld CHF table, with K 1

Suggested Approach to CHF n Use the 2006 Groeneveld CHF table, with K 1 (the newer one), K 2, and K 4, as provided above. n Evaluate the CHF table at the power that produces CHF, i. e. , CHF power = channel power. n Use RELAP 5, or other suitable code, to predict flow. If flow extrapolation is needed, be conservative. n NUREG-1537, Part 1, Appendix 14. 1, page 5 recommends minimum CHF ratios of at least 2. 0 for reactors with engineered cooling systems. TRIGA reactors with natural-convective primary flow do not have engineered cooling systems. A minimum CHF ratio is under discussion. 30

Conclusions n Flow Rate: – For the hexagonal pitch TRIGA reactor, the RELAP 5

Conclusions n Flow Rate: – For the hexagonal pitch TRIGA reactor, the RELAP 5 flow rate predictions are greater than the STAT predictions, especially at power levels approaching CHF conditions. n CHF – There is substantial uncertainty in the data. Correlation predictions differ greatly. – The 2006 Groeneveld table, with K 1, K 2, and K 4 as outlined above, is judged to be the best choice for TRIGA reactors. n For the hexagonal pitch TRIGA reactor: – The proposed 2006 Groeneveld CHF and RELAP 5 flow combination (column A of the previous table) predicts 62. 1 k. W/rod. – The traditional method of using the STAT code with the Bernath CHF correlation predicts 37. 1 k. W – Thus, in this example, the proposed method predicts the CHF power to be 67%, i. e. , (62. 1/37. 1 – 1) × 100%, greater. 31