paul scherrer institut€¦ · aea-trs-1003 psi report 72 oecd-loft large break loca experiments:...

45
CHIOOO^IO 1 = 8 -LUI a << TSi *£. Paul Scherrer Institut Labor für Thermohydraulik OECD-LOFT Large Break LOCA Experiments: Phenomenology and Computer Code Analyses I. Brittain, S.N. Aksan m%&^y:>;. Paul Scherrer Institut Telefon 056/99 2111 Würenlingen und Villigen Telex 82 7414 psl ch

Upload: others

Post on 22-Aug-2021

0 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

CHIOOO^IO

1 = 8

-LUI a < <

TSi — *£.

Paul Scherrer Institut

Labor für Thermohydraulik

OECD-LOFT Large Break LOCA Experiments: Phenomenology and Computer Code Analyses

I. Brittain, S.N. Aksan

m%&^y:>;.

Paul Scherrer Institut Telefon 056/99 2111 Würenlingen und Villigen Telex 82 7414 psl ch

Page 2: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

AEA-TRS-1003 PSI Report 72

OECD-LOFT LARGE BREAK LOCA EXPERIMENTS:

PHENOMENOLOGY AND COMPUTER CODE ANALYSES

I Brittain* S N Aksan+

# AEA Technology + Paul Scherrer Institute Winfrith Technology Centre CH-5232 Dorchester Villigen PSI Dorset DT2 8DH Switzerland UK

Würenlingen und Vi l l i gen

August 1990

Page 3: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

(ü)

OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES

I Brittain, AEA Technology, Wlnfrith, UK

S N Aksan, Paul Schcner Institute, CH-5232 Villigen PSI, Switzerland

ABSTRACT

Large break LOCA data from LOFT are a very important part of the world database. This paper describes the two double-ended cold leg break tests LP-02-6 and LP-LB-1 carried out within ttie OECD-LOFT Programme. Tests in LOFT were the first to show the Importance of both bottonrup and top-down quenching during blowdown in removing stored energy from the fuel. These phenomena are discussed in detail, together with the related topics of the thermal performance of nuclear fuel and its simulation by electrical fuel rod simulators, and the accuracy of cladding external themocouples. The LOFT data are particularly important in the validation of integral thermal-hydraulics codes such as TRAC and RELAPS. Several OECD partner countries contributed analyses of the large break tests. Results of these analyses are summarised and some conclusions drawn.

Paper presented at the Open Forum on the Achievements of the OECD-LOFT Project, 9-11 May 1990, Madrid

Page 4: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

(iü)

CONTENTS

Page

SECTION 1 INTRODUCTION 1

2 LARGE BREAK LOCA TESTS LP-02-6 AND LP-LB-1 2

2.1 Test LP-02-6 2

2.2 Test LP-LB-1 3

3 EFFECTS OF PRIMARY COOLANT PUMP OPERATION 4

4 THE BLOWDOWN BOTTOM-UP CORE QUENCH 5

5 THE BLOWDOWN TOP-DOWN QUENCH 6

6 NUCLEAR FUEL VERSUS ELECTRICAL FUEL ROD 7 SIMULATORS: SIMULATION LIMIATIONS DURING BLOWDOWN

7 EFFECT OF CLADDING SURFACE THERMOCOUPLES ON BLOWDOWN 8 HEAT TRANSFER

8 REFLOODINCJkND BOIL-OFF: EXTERNAL CLADDING 13 THERMOCOUPLE EFFECT, AND NUCLEAR FUEL ROD AND ELECTRICAL HEATER ROD BEHAVIOUR

9 COMPUTER CODE ANALYSES 14

9.1 Analyses of LP-02-6 15

9.2 Analyses of LP-LB-1 17

10 GENERAL CONCLUSIONS 20

11 REFERENCES 20

Page 5: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

1

1 INTRODUCTION

The LOFT experimental PWR reactor facility was specially designed for thermal-hydraulic transient tests up to and including a double-ended (200Z) cold leg LOCA. The principal objectives at the beginning of the USNRC programme were to assess the adequacy of the engineered safety systems and to validate the computer codes used in safety analysis. At that time the emphasis was on large break LOCA. The first (low power) nuclear experiment L2-2 was a 200Z cold leg break and was carried out in December 1978. During preparations for the corresponding high power experiment, L2-3, the Three Mile Island (TMI) accident occurred. This led to a considerable change in the LOFT programme, with the emphasis switching to small break LOCA and plant transients. Large break test L2-3 was carried out in May 1979, but only one further test, L2-5, was carried out before the end of the NRC programme. It was realised that the LOFT tests were a very important part of the (small) worldwide integral test database for large breaks, and two further tests were therefore carried out within the international OECD-LOFT project. One of these tests, LP-02-6, was designed to be representative of US Appendix 'K' licensing limits. This test was carried out in October 1983. The second test, LF-LB-1, was based on UK licensing assumptions, in particular the "minimum safeguards'* ECCS injection assumption. This test was carried out in February 1984.

The LOFT facility has a number of strengths and weaknesses as far as large break LOCA tests are concerned. On the positive side the most Important point is that the facility is a real FWR using real PWR fuel. Tests can be carried out from full power with appropriate coolant conditions. An extensive Instrumentation system provides information on coolant flows, temperatures and densities, and clad temperatures throughout the core. A few centreline fuel thermocouples are also provided. Some weaknesses arise from the fact that LOFT is not an exact replica of a 4-loop FWR. Table 1 shows some of the scale factors.

TABLE 1

Parameter

Power (LB-1)

Coolant Volume

Upper Plenum, Volume

Vessel cross-section area

Fuel length

Number of SG loops

Scale to 4-

1

1

1

1

1

-loop plant

: 63

: 45

: 24

: 30

: 2.2

: 4

Figure 1 shows the LOFT system configured for large break LOCA experiments. A major atypicality is the short (1.68 m) core. For gravity-driven reflood the relative vertical positions of cold leg nozzles, core and downcomer will be important. Figure 2 shows a sketch of the 4-loop FWR vessel with the elevations of the LOFT vessel superimposed using the

Page 6: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

2

bottoms of the cold legs as the common point. It can be seen that the LOFT core sits approximately at the centre of the 4-loop PUR core and is therefore effectively lower down in its vessel than is the 4-loop FWR core. The figure also shows that the LOFT vessel upper head is relatively much smaller than In the 4-loop plant. A possible further problem is the use of externally-mounted thermocouples to measure clad temperatures, and this will be the subject of a later section of this report.

Tests LP-02-6 and LP-LB-1 will be described in Section 2 below, while Sections 3 to 8 will discuss in more detail the important phenomena seen in these tests and the relevance to large break LOCA in the full size plant. One of the primary objectives of the LOFT programme was to provide validation data for integral system codes such as TRAC and RE1AP5. Several participating countries carried out analyses of the large break tests using these codes, and the results and lessons learned from these analyses are summarised in Section 9. Finally, Section 10 gives some overall conclusions from the LOFT large break test programme.

2 LARGE BREAK LOCA TESTS LP-02-6 AND LP-LB-1

2.1 Test LP-02-6

This test was a double-ended cold leg break with conditions representative of US licensing limits, namely loss of offsite power coincident with LOCA and minimum US emergency core coolant (ECC) injection. The peak linear heat rate was 49 kW/m, corresponding to a core power of 46 MM. All the fuel pins in the central bundle except the outer rows were pre-pressurlsed with helium gas to 2.41 HPa, typical of beginnlng-of-life conditions. After pump trip the pumps coasted down normally until the speed reached 750 rpm, at which point the flywheels were disconnected to minimise possible damage to the rotors. At this speed the pumps produce a very low head. After flywheel disconnection the pump resistance will be different in the turbine mode. The main objectives of this test were to provide data for computer code assessment and to provide data on fuel clad ballooning under realistic reactor conditions.

As usual, the experiment was initiated by opening the quick-opening blowdown valves in the broken hot and cold legs. The reactor scrammed automatically at 0,1 seconds, when the pressure had fallen to 14.8 MPa. The subsequent sequence of main events is given In Table 2. Rapid voiding of the core led to dryout within 1 second, followed by near-adiabatic heat-up of the cladding to an indicted peak temperature of 1074K at 4.9 seconds. An important feature of this transient was the surge of liquid into the bottom of the core beginning at about 5 seconds and persisting for about 5 seconds. This phenomenon had been seen in the earlier tests L2-2 and L2-3, and occurs when pumped flow from the intact loop exceeds the vessel-side broken cold leg flow. A more-detailed description and analysis of this phase of the transient will be found in a later Section. The bottom-up surge of liquid produced considerable cooling and quenching of the core. Figure 3 shows traces from four clad thermocouples attached at various elevations to fuel pin 5G06 in the central bundle. The thermocouples at elevations of 11, 30 and 45 inches quenched In sequence over a period of about 4 seconds. However, the thermocouple at 62 inches did not quench, although it did indicate significant cooling. Overall the thermocouples indicated quenching over the lower 2/3 of

Page 7: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

3

the core in the 5-10 second interval. Again, this behaviour and the possible effects of the external clad thermocouples are discussed in «ore detail in a later section. The whole core was in dryout again by about 11 seconds. At about IS seconds a top-down liquid flow caused further quenching. Figure 3 shows a clear quench of the thermocouples at 62 and 45 inches. The thermocouple at 11 inches shows significant cooling, but does not actually quench. Overall, this second quench extended over about the top 1/3 of the core. Final dryout of the core occurred at about 20 seconds. The ECCS was activated normally, reflood began at about 30 seconds and the entire core was quenched at 56 seconds. Peak clad temperature in reflood was 840K.

A unique feature of the LOFT facility is the presence in the core of a number of self-powered neutron detectors (SPNDs). These are devices which respond to neutron and gamma fluxes, and, as the name implies, are normally used to measure neutron flux. However, the fluxes arriving at the detector are quite strongly dependent on the local fluid density. It is therefore possible to use the response of the SPND to infer how the fluid density is changing [l]. To extract actual fluid density from the SPND readings requires a mixture of calibration and theoretical calculation, and the resulting values will not be very accurate. Nevertheless the qualitative information from these devices is very valuable for confirming the picture built up from the thermocouple responses. Figure 4 shows the fluid density deduced from the SPND at the 11 Inch elevation in location 5D08 compared with the thermocouple trace from the same elevation in location 5606. It is very clearly seen that the period of quench during blowdown coincides precisely with a period of increased coolant density.

One of the main objectives of this test was to determine whether clad ballooning or rupture would occur in pressurised fuel. Figure 5 shows the maximum cladding teaperatures measured in the test plotted against the differential pressure across the clad. The hatched area on the figure is the region of predicted onset of fuel damage based on separate effects tests, taken from Reference 2. It can be seen that no fuel rod damage would be expected in LP-02-6 even allowing for some uncertainty in the measurement of clad temperatures. This was indeed the result of the test as shown by failure to detect fission product activity in the primary coolant and by post-test examination of the fuel rod bundle.

2.2 Test LP-LB-1

This test was a double-ended cold leg break with conditions relevant to UK licensing limits, naaely loss of offsite power coincident with LOCA and ECC injection rates determined by the UK "minimus safeguards" assumption. The values used were 70Z of the accumulator liquid volume and 50X of the pumped ECC injection used in LP-02-6. The peak linear heat rate was 51.7 kW/m, corresponding to a core power of 49.3 MW. The bundle containing pre-pressurised fuel pins used in LP-02-6 was replaced by a bundle of unpressurlsed pins for this test. One objective of the test was to maximise the fraction of the core that would not rewet during blowdown, thereby maximising the peak clad temperatures during reflood. Thus the primary coolant pumps were tripped and disconnected from their flywheels within the first second of the transient to prevent, or at least minimise, the upflow cooling/quenching phase during blowdown. Other stated objectives [3] of this test were:

Page 8: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

4

1 To provide further data to assist in the evaluation of the degree of conservatism provided by the assumptions used in Evaluation Model analyses.

2 To provide further data for BE code assessaent.

The accident Initiation and reactor trip were essentially the saae as In LP-02-6. Differences occurred fro« about 5 seconds onwards due to the different primary pump operation. The sequence of main events is given in Table 3. Figure 6 shows traces froa four clad thermocouples attached at various elevations to fuel pin 5G06 in the central bundle. Although there was clearly no massive cooling end quenching in the 5-10 second period there Is evidence of some enhanced heat transfer from the fuel. The thermocouple trace for the 11 inch elevation shows a cooling of ~ 50K between 5 and 7 seconds, while the thermocouple at 30 Inches shows a distinct change of slope. Thus, there was clearly an increase in forward convective flow at this time, but insufficient to cause significant cooling and quenching of the fuel. At about 13 seconds the core heat-up was stopped by a fallback of liquid froa the upper plenum. This downflow was sufficient to quench the thermocouple at the 62 Inch elevation almost immediately, and the one at 45 inches also quenched about 2 seconds later. The two lower thermocouples show evidence of enhanced cooling at this time, but only fall in temperature by a few tens of degrees. The maximum temperature during blowdown was 1261K, reached just before the top-down cooling began. Overall, approximately the top 1/3 of the central bundle was quenched by the top-down flow, while quenching was considerably •ore extensive in the lower power peripheral bundles. The entire core was in dryout again after 25 seconds. The ECCS injection turned over the peak clad temperature at about 27 seconds and began to quench the fuel at 34 seconds. Quenching proceeded froa top and bottom of the core, and was coaplete at 72 seconds. The peak clad temperature reached during reflood was 1257K.

3 EFFECTS OF PRIMARY COOLANT PIMP OPERATION

The LOFT large break LOCA experiments have shown that the early bottom-up quench Is a result of a flow reversal in the core inltated by a change in the fine balance of coolant Inflow and outflow from the dovncoaer. Review of the break flow [4] In all of the large break experiments shows that the break flows are dependent on primary system pressure and coolant temperature upstream of the break. These parameters were nearly the same for all the experiments. Therefore, the early bottom-up fuel cladding quench depends on hydraulic parameters within the reactor vessel and intact loop.

Three operational modes of the reactor coolant pumps were used in the LOFT large break experiments. These modes are (a) continuous puap operation, (b) early pump trip with typical pump coastdown, and (c) early pump trip with fast pump coastdown (decoupled flywheels). The early quench did not occur In the experiments (L2-5 and LP-LB-1) in which the pumps were tripped within 1 second of transient initiation and disengaged from the flywheels (Figure 7). The experimental results show that the early quench is a function of the pump operation mo4e, puap characteristics, and initial flow conditions In the intact loop.

The first two large break experiments, L2-2 and L2-3, were conducted with reactor coolant pumps running. The coolant mass flow rate in the cold

Page 9: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

5

leg reaalned alaost constant until 6 seconds during both experlaents. Due to this and the decreasing aass flow rate in the broken loop cold leg, the intact loop cold leg aass flow rate exceeded the broken loop cold leg aass flow rate. The difference In the flow rates resulted in an excess of 700 kg of water being delivered to the downcomer between 4 and 6 seconds. The excess aass of water resulted in the propagation of a density wave upward through the core during this tlae. In these experiaents the early quench das also the aost complete, extending through the entire core radially ana axially.

In experiment LP-02-6 the pumps were tripped at the beginning of the transient and allowed to coast down under the Influence of the flywheels. Figure 8 shows the aass flow rate in the intact loop cold leg. This mode of pump operation provided enough coolant and head to initiate the bottom-up quench, but the quench front did not propagate through the entire core [*]. This is also Influenced by the higher saturated steam flow as in experiment LP-FP-1.

In experiment L2-5 the reactor coolant pumps were also tripped at about 1 second but the flywheels were disconnected from the pumps resulting in a very fast pump coastdown. Consequently, there was no time interval early in the transient during which the Intact loop cold leg mass flow exceeded the broken loop cold leg mass flow. Therefore, not only did positive core flow fail to be re-established, but also core flow was sufficiently near stagnation that the fuel rods did not rewet. The peak cladding temperature measured during the L2-5 experiment is compared to the other large break experiments in Figure 7 and confirm the importance of the primary system coolant pumps in controlling the core flow during the large break LOCA. Experiment LP-LB-1 was performed with reactor coolant puaps in the same operation mode as in Experiment L2-5. This experiment also did not contain a bottom-up core quench, except at the very bottoa (2 Inches) of the core.

Experlaent LP-FP-1, in which the reactor coolant puaps were tripped early and disconnected froa the flywheels, contained an early quench. The quench occurred in this case because of the higher initial, or steady state, aass flow In the Intact loop. The additional aass flow caused the magnitude of the intact loop cold leg flow to be larger than the broken loop cold leg flow at the tlae of transition to saturated critical flow at the break.

4 THE BLOWDOWN BOTTOM-UP CORE QUENCH

The early core cooling during blowdown was observed for the first tlae during L2-2, the first LOFT nuclear experlaent. This behaviour was definitely different froa the expected core thermal behaviour. At that tlae according to the understanding of reactor system behaviour during large break LOCA the cladding temperature should rise rapidly after LOCA initiation and reach a aaxlaua during the blowdown phase. A precursory cooling due to droplet entralnaent froa the lower plenum would reduce the core heatup and finally the cladding temperatures would be quenched to the saturation teaperatures as result of core reflood with ECC water. This classic large break LOCA scenario was supported by experimental evidence froa facilities such as SEMISCALE and with code analysis.

Experiment L2-2 and then later other LOFT experiaents with similar boundary conditions have shown that the cladding temperature rises as expected

Page 10: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

6

but the temperature transient Is arrested after a few seconds and followed by a core-wide bottoa-up cladding quench. The cladding teaperature rises again later and reaches another aaxlaua during the reflood phase. The highest peak cladding temperatures are measured during the first heatup. The main reason for this unexpected core thermal behaviour is the hydraulic process within the primary system and the resulting very high upward mass flows in the core during the blowdown phase, within the first 10 seconds. This behaviour confirms the conservatism of the licensing models. Also the early quench provides a unique challenge to the best-estimate codes wherein the thermal-hydraulic phenomena are physically modelled. The thermal-hydraulic conditions that resulted in the early blowdown quench in LOFT are reviewed in some detail in References [4] to [ll] to establish the nature of the rapid core flow, at high system pressures, that resulted in quenching of the fuel rods. Although the core hydraulic conditions were not directly measured during the LOFT loss of coolant experiments, it is possible to obtain certain Insight into the core flow conditions that resulted In early cladding quench by using the data from the self-powered neutron detectors (SPNDs), lower plenum velocity and aoaentua flux measurements (used in LP-02-6 and LP-LB-1 tests, [l2]) and temperature measurements in the downcomer, core and upper plenum regions. In sumaary, the rapid cladding cooling was primarily a result of low-quality, high upward core flow at a time when the system pressure was still relatively high (7 HPa). The velocity of the density wave through the core is calculated to be approximately 1.8 m/s (and the core inlet mass flux is estimated to be about 615 to 1050 kg/a2 sec). This moving cooling wave causes a quench propagation of approximately 1.0 to 1.5 m/s in the core. The resulting measured peak cladding temperature response during the LOFT L2-2, L2-3 and LP-02-6 large break loss-of-coolant experiments is coapared In Figure 9. From these experiments, it can be seen that during the first 12 seconds the cladding is observed to quench rapidly due to core positive flow (Figure 10) and remain quenched for several seconds. More than 60% of the stored energy in the core was transferred froa the fuel rods by the high core Inlet flow causing a cladding rewet.

Even though correlation of the reactor vessel flow measurements, upper plenum thermocouple measurements and SPND measurements suggest that the early quench characteristics during LOFT experiments were hydraullcally controlled, the accuracy of the cladding temperature data has been questioned because of hypothesized selective cooling of the external surface thermocouples which aay cool the entire fuel rod due to a fin effect. To evaluate the accuracy of the measured LOFT cladding temperatures and to provide an Independent database for the rapid cooling transients under known high pressure hydraulic conditions, separate effect experiments were conducted and the effects of cladding external thermocouples on the early quench phenomenon were investigated. Additionally a special model to simulate the behaviour of the external thermocouples attached to the LOFT fuel rods was developed. Summary of this experimental and analytical work will be given in a later section.

5 THE BLOWDOWN TOP-DOWN QUENCH

The core cooling during blowdown contains two phenomena: the bottoa-up cooling, discussed In the previous sections, and the top-down cooling. Again, LOFT experiments were the first to show the top-down cooling phenomenon. Cladding temperatures measured during the LP-02-6 experiment, as shown in Figure 3, indicate a second quench in the upper part of the core which moved

Page 11: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

7

downwards and revetted the cladding at the 45 inch elevation at 17.5 seconds. The thermocouple measuring the highest temperature at the 30 inch elevation did not indicate the top-down quench. The fuel bundles closest to the intact loop hot leg, which is a source of water which drains into the reactor vessel, clearly indicate top-down quench. As shown in Reference [4], the top-down quench is multidimensional in contrast to the bottom-up quench which can be treated as one-dimensional as it rewets the centre fuel module and the peripheral modules at the same time.

Analogous top-down quench phenomena were detected in other LOFT large break experiments. Experiment L2-5 was performed with a rapid reactor coolant pump coastdown to prevent the bottom-up quench, however this operation did not prevent the top-down quench. The upper half of the centre fuel assembly was rewet between 12 and 15 seconds by liquid draining from the upper plenum and Intact loop hot leg and pressuriser into the core. The top-down quench was also present in the peripheral fuel modules.

The top-down quench phenomenon is concluded to be significant because it decreases cladding temperatures which reduces the time needed to quench the fuel cladding during core reflooding with the ECCS. Further, there are almost no separate effects data on this phenomenon, and it may be even more important in the full scale plant because of the availability of much larger amounts of water from the upper head.

6 NUCLEAR FUEL VERSUS ELECTRICAL FUEL ROD SIMULATORS: SIMULATION LIMITATIONS DURING BLOWDOWN

The typicality of the blowdown quench behaviour of a solid-type electrical heater rod relative to that of a nuclear fuel rod has been questioned because of the different thermal properties and lack of a siaulated fuel-pellet cladding gap. In this respect, LTSF experiments investigating the blowdown quench behaviour of a SEMISCALE solid-type heater rod with only internal thermocouples were used as a basis for evaluation of the model calculations. The details of these calculations using the RELAP4/MOD6 computer code are given in References [7] and [l3]. Having established the validity of the heat transfer models to calculate the initial cooldown rate of a quench, a series of RELAP4 calculations was performed to compare the initial cooldown rates of the nuclear fuel with gap, REBEKA cartridge-type electrical heater rod with gap, and SEMISCALE solid-type electrical heater rod for rapid cooling transients. These calculations under typical LTSF single-rod test conditions were performed by substituting the individual rods in the RELAP4/M0D6 LTSF model. The calculated results for one of the LTSF experimental case (Experiment no 12) are given in Figure 11. These calculations show that the REBEKA heater rod is expected to simulate the nuclear fuel rod behaviour very well for the conditions Investigated. The nuclear fuel rod with gap was calculated to cool approximately five times faster than the SEMISCALE solid-type heater rod. These comparison calculations [13] were carried out over a range of inlet flooding velocities (2 to 6 ri/s) and the predicted initial cooling rates are summarised In Figure 12. Also Indicated in this figure are the initial cladding cooling rates measured on nuclear fuel rods from the PBF Thermocouple Evaluation Experiment Series [14]. The limited number of nuclear fuel rod data suggest that Che calculated cooldown rates may be 10 to 20% too high. The results of those FBF experiments also indicate that t'.ie thermal decoupling of the

Page 12: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

8

cladding and fuel «as apparendy significant, allowing the cladding to quench rapidly during the blovdown phase. This thermal decoupling of fuel and cladding deaonscrates Che laporcance of in-pile experiments, and of out-of-plle experiments where the fuel-to-cladding gap is properly simulated.

Additional experiaents were conducted In LTSF using a REBEKA cartridge-type fuel rod simulator with gap, zircaloy cladding and thermal diffusivity much closer to the nuclear rod diffusivity. The experimental results without external thermocouples show very rapid cooling (150 to 200 K/s) and quench times (2 to 3 seconds) similar to the nuclear fuel rod data at 4 m/s Inlet flooding rates (Figure 12). A comparison of Che cladding temperature response of the REBEKA rod with external thermocouples and a nuclear fuel rod wich external thermocouples, where Che initial temperatures of the rods prior to quenching were about the same (900K) is shown In Figure 13. Similar results exist for rods without external thermocouples as mentioned above. The quench behaviour of the REBEKA rod is similar to that of a nuclear fuel rod, which is also confirmed by the results of calculations performed with the RELAP4/M0D6 code (Figure 11). Additionally Che bundle tests performed In the LTSF test facility produced valuable evidence with respect Co che behaviour of cartridge-type REBEKA and the solid- type FEBA fuel rod simulators (Figure 14). The analysis indicates that the solid-type heater rod temperatures are controlled by convective heat transfer at the cladding surface. This is due to che high thermal diffusivity of the rod allowing rod internal energy to be transferred rapidly to Che cladding. As the rod begins Co cool down in film boiling, internal rod energy is conducted Co Che cladding about as fast as the surface convective heac transfer removes Che energy; thus, Che cladding cools slowly since the rod's energy must be transferred before the cladding temperature Is low enough to allow surface quenching. This Is clearly seen in Figure 11, where the film boiling cooldown lasts for 8 seconds. The nuclear rod, by contrast, is Internally conduction limited by the greater thermal resistance of Che U02 fuel and the fuel-cladding gap thermal resistance. The inability of the nuclear rod to transfer rapidly rod internal energy to the cladding, together with a much smaller zircaloy cladding heat capacity, significantly changes the energy balance at the cladding surface, causing a more rapid cooldown during film boiling (see Figure 11). In comparison, at these high flow rates the nuclear rod cooling la controlled more by the cladding stored energy, while Che cooling of Che solid-type heater rod is controlled more by total rod internal energy.

7 EFFECT OF CLADDING SURFACE THERMOCOUPLES ON BLOVDOWN HEAT TRANSFER

Two differenc sets of experiments were conducced In Che LOFT Test SupporC Facility (LTSF) and another set of experiments with nuclear fuel in the Power Burst Facility (PBF) at Che Idaho National Engineering Laboratory (INEL) in order to Investigate the effects of cladding external thermocouples on the early quench phenomenon.

The LTSF quench experiments provide a simple geometrical configuration with well-quantified inlet hydraulics and the capability to maintain the system pressure as high as 7MPa. A detailed description of this test facility is given in References [l5] and [16]. The first series of tests in LTSF was performed using a single SEMISCALE fuel rod simulator which is representative of "solid" internally heated fuel rod simulators used in many light water reactor research projects. Four external thermocouples had been laser welded

Page 13: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

9

to the outer surface of the solid heater rod with junctions and attachments similar to LOFT cladding thermocouples (Figure IS). The heater rod also had four Internal cladding thermocouples to measure the rod's temperature response. In the same test facility, identical experiments were conducted with and without the cladding external thermocouples. The results of one of the 20 quench tests conducted in LTSF are shown in Figure 16, for the initial conditions given in the figure. The data presented were taken from cladding Internal and external thermocouples located at the heater rod hot spot. The time of coolant arrival at the thermocouple location is indicated by the rapid change in the test section gamma densitometer response. Thus, the quench times can be estimated with respect to coolant arrival. During the high pressure (7MPa) tests the heater rod with external thermocouples consistently quenched in about half the time required by the heater rod without surface thermocouples. It can also be seen from the data that the surface thermocouple is selectively cooled and quenches much sooner than the cladding, as indicated by the internal thermocouple reading. Further detail': of these experimental results can be found in References [15] and [17 ]. An analysis of the LTSF early quench tests without the external surface theisocouples was also performed [7], [13] and as a result of this analysis it was concluded that the solid-internal heater rods of the SEMISCALE design type cannot simulate the rapid quenching of a nuclear rod, due to the relatively high thermal dlffusivity. The typicality of the blowdown quench behaviour of a solid-type electrical heater rod relative to that of a nuclear fuel rod has been questioned because of these different thermal properties and lack of a simulated fuel-pellet cladding gap. In this respect, the LTSF experiments which investigated the blowdown qu nch behaviour of a SEMISCALE solid-type heater rod can be used as a basis for evaluation of code model calculations but not for direct extrapolation to reactor conditions [18]. Thus relating these LTSF high pressure quench test results to LOFT must be done with care, especially due to differences in the thermal response of the electrical heater rods.

Another separate effects experiments programme was conducted in the LTSF test facility with a nine-rod C3x3) bundle of electrical heater rods. A REBEKA cartridge-type heater rod ana a FEBA solid-type heater rod (similar to SEMISCALE heater rods) were each tested in the centre position in the nine-rod bundle, which provided a geometry and thermal-hydraulic environment more typical of a nuclear fuel rod cluster. The REBEKA heater rod has Zircaloy cladding and aluminium oxide pellet construction with a pellet-cladding gap to simulate the thermal characteristics of . nuclear fuel rod. The quench behaviour of the REBEKA heater rod has been compared with that of a LOFT nuclear fuel rod through calculations using the RELAF4/M0D6 computer code [13] and with that of the Power Burst Facility nuclear rod through experiment data evaluation [16]. Both the analytical and experimental data comparisons showed that the REBEKA heater rod provides good simulation of nuclear fuel rod thermal response for rapid cooling at high pressure in a blowdown phase transient. The REBEKA heater rod was tested with and without cladding external thermocouples. The main objectives of this experimental programme were to evaluate the effect of cladding external thermocouples on the early blowdown phase quench behaviour of a cartridge-type nuclear fuel rod simulator, to determine how accurately cladding external thermocouples measure cladding temperature during a blowdown phase quench, and to compare the high pressure quench behaviour of a cartridge-type heater with that of a solid-type heater rod under similar thermal-hydraulic conditions to those that occurred

Page 14: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

10

during the blovdovn phase of LOFT experiments. The experimental programme and the results of these tests are given in detail in [l6]. The results of the experimental programme indicate that cladding external thermocouples had a negligible effect on the cooldown rate and quench behaviour of a REBEKA cartridge-type heater rod under rapid (1 to 2 m/s) flooding conditions at high pressure (Figure r*). Also, the cladding external thermocouples are selectively cooled during the quenching process and do not accurately measure cladding temperature during this part of the transient. Since the REBEKA rod has been shown to simulate satisfactorily the thermal response of a nuclear rod [7], these results are considered applicable to LOFT nuclear fuel rods. Consequently, the value of LOFT external thermocouple data in validating computer models during quenching is somewhat limited. Further results also show that the quench behaviour of a FEBA solid-type heater rod is significantly different to that of a cartridge-type heater rod. The REBEKA rod quenched in less than 3 seconds from about 900K, whereas the FEBA heater rods experienced an extended period (10 seconds) of precursory cooling before quenching at about 700K (Figure 14). Since it has been shown that the REBEKA rod provides a good simulation of a nuclear rod, it is inferred that solid-type electrical heater rods do not provide a good simulation of the thermal response of nuclear fuel rods under high pressure, rapid cooling thermal-hydraulic conditions that occur during the blowdown phase of a large break loss-of-coolant accident in the LOFT experiments.

Three series of light water reactor fuel behaviour tests (Thermocouple Effects Test series TC-1, TC-3 and TC-4) were performed in the Power Burst Facility (PBF) at 1NEL specifically to evaluate the influence of cladding surface thermocouples 01 the thermal behaviour of nuclear fuel rods under LOCA conditions. A total of twelve tests were conducted. Each of these tests were performed with four LOFT-type fuel rods contained in individual flow shrouds while two of the rods were instrumented with four LOFT cladding external thermocouples, with junctions located near the high power region of the fuel rods. All four rods were also instrumented with internal thermocouples, the junctions being at the same axial level as the external thermocouples. Details of the experiment design, conduct and results are given in References [14], [19] and [20]. Evaluation of the measured temperature drop across the cladding indicates that the cladding surface thermocouples measured cladding surface peak temperatures during blowdown that were only slightly lover (20 to 30K) than the actual temperatures. However, the surface thermocouples influenced the cladding temperatures during the blowdown phase of the TC tests in two respects. The initial effect of the surface thermocouples is to delay the time of occurrence of CHF on the externally instrumented rods with respect to the time-to-CHF on the bare rods; consequently the externally instrumented fuel rods experience a longer period of high cladding surface heat transfer which in turn reduces the stored energy in the fuel rods at the time-of-CHF. In addition to the effect of delaying the time-to-CHF, the cladding surface thermocouples increase the surface heat transfer which in turn reduces the stored energy in the fuel rods during the heat-up phase due to steam cooling and further Influence the fuel rod thermal response through a "fin cooling" effect. Peak temperatures measured during blowdown were about 7SK lower for each second of delay in CHF, and the improved cladding heat transfer due to fin cooling reduced the peak temperatures during the blowdown phase of the PBF tests by 100 to 115K. The major contributor to the surface thermocouple effect during blowdown quench appears to be the reduced blowdown peak temperature due to the delay in

Page 15: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

11

time-to-CHF and Initial fin cooling, ie the lower cladding temperatures just prior to the blowdown quench. Garner [14] estimated that 35-58Z of this reduction is related to the delay in time-to-CHF with peak cladding temperatures educed by 74K for each second of delay, and 42-65Z of peak cladding temperature reduction is due to fin cooling. The PBF-TC test results also indicate that the effect of surface thermocouples on rod thermal response during blowdown quench (Figure 17) appears to be relatively small and, especially if the quality of the slug flow is very low and the rods are quenched extremely fast, the effect even d sappears. In addition, the effect of surface thermocouples on nuclear fuel rod thermal response during blowdown quench appears to decrease as the rod initial power decreases and, at low power, the effect again disappears.

As a further study of the thermal response during the blowdown phase, fuel centreline thermocouples were designed for placement in the LOFT Centre Fuel Modules (CFHs). These fuel centreline thermocouples were used In two LOFT large break LOCA experiments LP-02-6 (with early quench) and LP-LB-1 (no early quench). Some analysis of the results of fuel centreline thermocouples have been included in References [4], [21 ]. Figure 18 shows the response of the fuel centreline temperature to the measured cladding temperature at the same elevation for the LP-02-6 experiment. The fuel centreline temperature responds to the change in cladding temperature shortly after the blowdown quench at approximately 8 seconds and similarly after the reflood quench at approxlmtely 53 seconds. Figure 19 shows the fuel centerline temperature for two fuel rods which did not have externally mounted cladding thermocouples compared with the fuel centreline temperature shown in Figure 18. The INEL view is that these data show that the hydraulics cause complete fuel cladding quench and not just thermocouple quench or localized cladding quench. The UK takes the opposite view that these data in conjuction with pin heat conduction calculations actually confirm the UK belief that the external thermocouples quench 1-2 seconds earlier than the bare clad.

The LP-LB-1 experiment also contains some data showing the strong dependency of the fuel centreline temperature on the cladding temperature and heat transfer. In LP-LB-1 the early bottom-up quench phenomenon was suppressed. However, there was a partial top-down quench that occurred in the 10-30 seconds time interval and extended over approximately the top third of the core [21 ]. Figure 20 shows fuel centreline temperature at the 27 inch elevation for rods with and without cladding surface thermocouples. The centreline temperature behaviour indicates no quench in agreement with the measured cladding temperature. Figure 21 shows similar temperature data at the 43.8 inch elevation. A small early cooling occurs in this region as indicated by the cladding temperature in Figure 21 compared with that in Figure 20. The fuel centreline temperature is sufficiently sensitive to show even this small cooling. The data in Figures 20 and 21 during the final quench does show that the final quench occurs up to approximately 20 seconds earlier on fuel rods with thermocouples. These results are consistent with fuel rod results in PBF [14].

On effects of the LOFT external thermocouples on the blowdown quench phenomenon, the view of the UK differs and stronger conclusions are drawn. It is postulated that the external thermocouples are thermally decoupled from the clad and can Induce early quenching of the clad, rather than being true indicators of the clad temperature behaviour. Thus the UK view is that the

Page 16: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

12

LOFT blowdown quenches vere not as extensive as the measurements of the external thermocouples suggest, but some quenching of the rods would have occurred. A study has been carried out at Winfrith, and the development of a model to simulate the behaviour of the external thermocouples attached to the LOFT fuel rods during the blowdown phase of a large break LOCA has been undertaken [25], [26 J, [27] and [28J. To establish the model and determine the thermal coupling between the thermocouple and fuel rod, extensive use was made of the data of the two series of experiments (mentioned above) performed in the LTSF. The investigation of the LTSF experiments, and development of the model to simulate the behaviour of an externally mounted thermocouple and its application to a LOFT large break LOCA were performed in the following steps:

1 As the starting basis, the experiments with solid fuel rod simulators, in which there were no external thermocouples, were used to evaluate surface heat fluxes (hence heat transfer coefficients) from the rods to the coolant. It was assumed that the surface heat fluxes evaluted from the experimental series with no external instrumentation would apply to heater rods with external thermocouples under an equivalent set of independent variables, ie pressure, coolant flow rate, quality etc.

2 A 2-D finite element heat conduction model of the fuel rod, external thermocouple, and fuel rod-thermocouple weld was constructed using the general purpose heat conduction code TAU.

3 A 3-D model was constructed in order to improve on the calculated thermocouple response because of the fact that the thermocouples are only welded on the fuel rod intermittently along their lengths.

4 The thermal resistance of the weld between the thermocouple and rod, as the unknown parameter, was determined using the LTSF experiment data.

5 The 3-D fuel rod-external thermocouple model was then used to simulate the best available approximation of the cooling and quench behaviour during the blowdown phase of a LOFT large break transient (Figure 22).

As a result of the study, it is shown that it is possible to construct a 3-D model for a heater or nuclear fuel rod with simulation of external thermocouples. The calculations performed using this model indicate that a rapid cooldown and quench of the external thermocouple can occur while the cladding of the nuclear fuel remains dry, as a result of the poor thermal coupling that exists between the external thermocouple and the cladding. It la concluded that the behaviour of the external thermocouples can be predicted by a code using such a model, and under rapidly changing transient conditions the thermocouple cannot be relied upon to represent the cladding temperature of the rod accurately. In the particular case of LP-02-6 it is concluded that the thermocouples in the high temperature regions will quench 1-2 seconds before bare cladding quenches. This has an important consequence on the stored energy in the core at the end of blowdown.

Page 17: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

13

8 RELOODING AND BOIL-OFF: EXTERNAL CLADDING THERMOCOUPLE EFFECT, AND NUCLEAR FUEL ROD AND ELECTRICAL HEATER ROD BEHAVIOUR

Large break experiments in LOFT were intended to validate the performance of the emergency core cooling systems for the design basis loss-of-coolant accident. As discussed in previous sections the L2-2, L2-3 and LP-02-6 experiments showed that about 60% of the initial steady state stored energy is transferred to the primary coolant prior to emergency core coolant delivery to the core. The final core quenches are primarily due to accumulator fluid delivery. The characteristics of the relatively rapid (10 to IS cm/s) core quenching for the L2-3 and L2-S experiments are compared in Figure 23. The core reflood behaviour was very similar during both experiments even though a significant difference in Initial stored energy and cladding temperatures existed. All of the other LOFT large break loss-of-coolant experiments (L2-2, LP-LB-1, LP-02-6), showed similar type of rapid quenching behaviour during the core «flooding phase. A significant observation is that the reflooding rates always exceeded the 2.5 cm/s Evaluation Model licensing regulation limitation. Additionally the LOFT large break loss-of- coolant experiments showed that the biowdown hydraulics and heat transfer (early blowdown quench) is more important in removing the initial fuel rod stored energy than reflood heat transfer. The relatively rapid reflooding observed during the LOFT large break loss-of-coolant experiments was questioned because of the fin effect of external thermocouples. As a consequence, separate effect bundle reflooding experiments using electrical heater rods instrumented with both external and internal cladding thermocouples were conducted in the NEPTUN test facility [29] in Switzerland at the Paul Scherrer Institute (formerly EIR). The NEPTUN-I series of experiments were performed with five central heater rods instrumented with both external and internal cladding thermocouples and the NEPTUN-II series of experiments were performed only with internal cladding thermocouples. The results and comparison of experimental data from these two experiment series [30], [31 ], indicated that electrical heater rods Instrumented with LOFT external thermocouples experience preferential cooling during reflooding compared to heater rods with internal embedded cladding thermocouples (Figure 24). The effect Is reduced with higher reflooding rates (eg 15 cm/s). During the precursory cooling before the quench, the rods with external thermocouples show comparable temperature histories to the rods with internal thermocouples, but heater rods with external thermocouples quench at higher temperatures and earlier than other heater rods. An overall comparison between repeat experiments NEPTUN-I (five central rods equipped with external thermocouples) and NEPTUN-II (all thermocouples embedded in the cladding of the heater rods) is shown In Figure 25. Differences between simllai experiments are small especially during precursory cooling.

The ability of electrical heater rods to duplicate nuclear fuel rod thermal response during reflooding was also questioned because of the large differences in electrical and nuclear fuel rod thermal properties. In this regard the Halden Project Test Programme Instrumented Fuel Assembly 511 (IFA 511) In Norway, in the Halden Research Reactor, was designed to evaluate systematically the ability of electrical heater rods (SEMISCALE solid-type) to simulate the response of nuclear fuel rods during the heatup and reflood phases of a -arge break loss-of-coolant accident. The experiment rods were also instrumented with both external and Internal cladding thermocouples to determine if external thermocouples provide an accurate measurement of

Page 18: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

14

cladding temperature. The experiments consisted of one series with nuclear fuel rods [32] and two other series with electric heater rods, SEMISCALE solid-type rods [33], [34] and REBEKA heater rods with cladding gap [35]. In these experiments the nuclear fuel rods were quenched substantially earlier (four times faster) than solid-type electrical heater rods (Figure 26). The REBEKA heater rods closely simulated the actual nuclear fuel rod behaviour under «flooding conditions. Also, the electrical rod, unlike the nuclear fuel rod, is characterised by a well-defined quench. Experimental data also show that the response of the external thermocouples was significantly different than the comparative internal cladding thermocouples during reflooding at about 7 cm/a flooding rries. The indicated temperature of the external thermocouple was at least 50K less than that indicated by the internal thermocouples throughout reflood and the external thermocouples indicated quench 20 seconds earlier. The different thermal behaviour indicated by the external thermocouples was primarily caused by the fin cooling effects. Additional experiments were performed in the FEBA test facility (at KFK Karlsruhe, Federal Republic of Germany) within the SEFLEX test programme in order to quantify the influence of the design of different fuel rod simulators on cladding temperature transients under reflood conditions. These experiments were done by employing a solid-type FEBA heater rod bundle and a corresponding bundle of REBEKA fuel rod simulators with gap. The experimental data [36] indicate that the reflooding behaviour between the two bundles consisting of 5x5 FEBA and 5x5 REBEKA rods is significantly different (Figure 27). At an inlet flooding velocity of 3.8 ca/s, the influence of the rod design on the peak cladding temperature is around 100K. lower for REBEKA rods and the quench times are about 100 seconds shorter. The reason for the lower cladding temperatures and the faster quench front progression for the REBEKA rod bundles are the lower heat capacity of the zircaloy cladding and the pronounced decoupling of the cladding from the heat source due to the cladding gap and the Influence that this has on the channel hydraulic conditions.

Core uncovery (boil-off) experiments were conducted in the NEPTUN experiment facility which has already been mentioned above. The results of NEPTUN boil-off experiments at 5 bar [37] showed that the cladding thermocouples do not cause a significant cooling effect on the rods to which they are attached. The dry-out times of the Internal and external cladding thermocouples were within 10 seconds of each other at any axial elevation for all rods In the bundle. The cladding external thermocouples measure the cladding temperatures that would have been measured in their absence within 0 to - 20K (Figure 28). The experimental data from the IFA 511 experiments for the heet-up phase at low pressures showed that the response of the external and Internal thermocouples was nearly Identical through heat-up until temperatures exceeded 700K. However, after about 700K, the cladding surface temperature measured by the external thermocouple was lower than that measured by internal thermocouples and the difference increases thereafter (Figure 26). The measured cladding peak temperature was 25 to 40K less, for both electrical heater and nuclear fuel rods. These results confirm the findings of the NEPTUN boil-off experiments.

9 COMPUTER CODE ANALYSES

Computer code analyses of LOFT tests played a specific role in pre-test planning and in post-test analyses of the LB experiments, and continue to play

Page 19: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

15

a wider, more general role within code assessment programmes. Planning calculations and pre-test predictions were carried out using an early version of the TRAC code, TRAC-PD2/M0D1, for LP-02-6 at INEL [38,39] and for LP-LB-1 at Winfrith [40] and INEL [41]. However, because this version of the code is no longer used, these analyses will not be discussed further in this paper.

As part of the formal Project agreement, partner countries generally undertook to produce both an Experiment Analysis and Summary Report (EASR) and a Comparison Report of code analyses for each test. In the event, for LP-02-6 no EASR was produced, but Switzerland produced a Comparison Report [42]. For LP-LB-1 the UK produced an EASR [43] and a Comparison Report [44]. A brief summary of the analyses and findings is given below.

It is worthwhile noting at this point that the LOFT results themselves do not tell us what would happen in the full-size plant and therefore they cannot be used directly in the Licensing process. Best-estimate computer codes must be used to extrapolate to plant conditions, taking into account scaling and other factors inevitably introduced by the experimental procedures Involved. The Importance of the LOFT data and analyses is in code validation and the quantification of uncertainties, and the data will continue to be used for this purpose for many years.

9.1 Analyses of LP-02-6

Table 4 lists the pre- and post-test analyses which have been reported. The LANL calculations are reported in Reference 45. All the TRAC calculations use a 3-dlmensional representation of the reactor vessel together with a 1-dimensional representation of all other components. The RELAP5 and DRUFAN calculations perforce use only 1-dimensional representations. In these cases 3-dimensional effects are simulated to some

TABLE 4

Computer Codes Used For Analyses of LP-02-6

Computer Code

TRAC-PD2/M0D1 TRAC-PD2/M0D1 TRAC-PD2/M0D1 TRAC-PF1/M0D1 TRAC-PF1/M0D1 RELAP5/M0D2 RELAP5/MOD2 DRUFAN-2

Participant

USA-INEL USA-LANL Spain-ETS USA-LANL UK-Winfrith Japan-JAERI Switzerland-EIR Germany-GRS

Type

Pre-Test Pre-Test Post-Test Post-Test Post-Test Post-Test Post-Test Post-Test

extent. For example, downcomer and core are both subdivided into two parallel flow paths, with various degrees of cross-connection. Other modelling assumptions also show some variation, in particular the assumed total core bypass flow varied from 4% (Japan) to 10% (UK).

Page 20: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

16

All the code calculations give a similar overall picture for the early hydraulic behaviour during the test, in quite good agreement with the data. That is to say all codes predict rapid core flow reversal accompanied by core voiding, a return to positive flow at about S seconds and some top-down flow at about 15 seconds. The length of the blowdown phase is well-predicted showing that the critical flow calculations are probably adequate. There are differences in the detailed behaviour and hence deviations from the data. The TRAC calculations tend to show periods of complete core voiding whereas the RELAP5 calculations tend to show some liquid remaining in the core at all times. However the biggest differences, and the biggest deviations from the data, are in the core thermal behaviour. It will be recalled that the clad thermocouples indicated an extensive bottom-up quench beginning at 5 seconds, covering at least 2/3 of the core. Host of the calculations Indicate a cooling during this period, but insufficient to reach the temperature Tmln at which the codes would switch to a wetted surface heat transfer mode. The subsequent temperature histories are thereby in error, with temperatures and stored energies too high. Figures 29 and 30 are example comparisons of clad temperature predictions at the 31 inch elevation illustrating these points. These errors have a feedback effect on the hydraulic behaviour in reflood. In complete contrast the version of RELAP5 used by Japan and the version of TRAC used by Spain were modified to use very high values of Tmin, completely inconsistent with steady-state film boiling data and much higher than the homogeneous nucleatlon temperature for water, during the blowdown phase. These modified codes thereby more accurately reproduced the experimental thermocouple data during blowdown, but the calculations then show too little stored energy in the reflood phase, which leads to early final quenching of the core. [It should be noted that there is considerable doubt over whether the external clad thermocouples used in LOFT give an accurate picture of the clad temperature transient during blowdown. This topic is discussed in more detail elsewhere in this paper]. The calculations predict some top-down flow after about 15 seconds, but again the thermal response of the core in poorly predicted. In general, too little heat is removed in this phase also. The experimental data show that this top-down flow was highly asymmetric, with most flow being close to the intact hot leg. The 1-dimenslonal codes of course cannot model this. The 3-dimensional TRAC calculations do show the effect qualitatively, but are quantitatively poor. Additionally, recent work carried out at FSI [23] presents some improvements in predicting the LP-02-6 blowdown phase quench by changing the selection logic related to nucleate, transition and film boiling in RELAP5/HOD2. It should be noted that the frozen verlon of the code predicts no blowdown phase quench of the central high power regions under the same hydraulic conditions.

In refill and reflood the calculations show greater variation in behaviour. Most calculations agree with the data that very little liquid Is bypassed directly to the broken cold leg. However, all fall to predict the Intact cold leg density variations seen in the data. Presumably, this is an indication of poor modelling of condensation. The calculations are generally In agreement that very little liquid enters the core until the accumulator empties. All the RELAP5 calculations predict complete quenching of the core within a few seconds of each other, regardless of the aaount of stored energy being removed by the reflood. The TRAC calculations which failed to remove sufficient stored energy during blowdown show oscillatory behaviour during reflood, and a delayed final quench. The RELAP5 calculation frrm Japan, which was the only calculation to reproduce all the thermocouple data during

Page 21: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

17

blowdown, «as the only calculation to predict the reflood quench earlier than the data, with an atypical predicted core liquid content of only about 30Z. The German DRUFAN calculation did not extend Into the reflood phase.

9.2 Analyses of LP-LB-1

Table 5 lists the pre- and post-test analyses which have been reported.

TABLE 5

Computer Codes Used For Analyses of LP-LB-1

Computer Code

TRAC-PD2/M0D1 TRAC-PD2/M0D1 THAC-PF1/M0D1 RELAP5/M001 RELAP5/MOD2 RELAF5/H0D2 DRUFAN/VLUT

Participant

USA-INEL USA-1NEL UK-Vlnfrlth Italy-OofBologna Flnland-VTT Switzerland-EIR Germany-GRS

Type

Pre-Te.it Post-Test Post-Test Po»t-Test Post-Test Post-Test Post-Test

The modelling details are essentially the same »s noted for LP-02-6 above. The one major difference Is that the German submission Included a reflood phase calculation, made using the FLUT 1-dimensional code. The FLUT model is unique in having a single channel representation of the downcomer.

All the code calculations give satisfactory agreement with the data for flows in the Intact and broken loops during the blowdown phase. Indeed all the flow rates are within the experimental uncertainty bars. However, this does not mean that the core flow predictions are similar. All the 1-dlnenslonal calculations show some upward flow of liquid into the core after about 5-7 seconds. The TRAC calculations do not show this flow, in agreement with evidence from thermocouples at the bottom of the core. There Is wide variation in the predicted core thermal response during this phase. In one of the 1-D calculations a significant fraction of the cladding remained wet during the whole of the blowdown period. In another a large quantity of liquid was predicted to flow up through the core, but this was accompanied by a relatively small amount of heat transfer. In general, as in LP-02-6, the TRAC calculations show more rapid voiding of the core than do RELAP5 calculations.

The data show a top-down flow of liquid through the core central bundle beginning at about 12 seconds. In the peripheral bundles there is evidence of more than one top-down flow in the period up to 16 seconds. Only one of the 1-D calculations predicts any downflow of liquid at the top of the core during blowdown, and then not until after 16 seconds. They all however predict downflow at the bottom of the core a» a result of erroneously predicting upflow earlier in the transient. TRAC-PFl/MOOl correctly predicts the downflow In the central bundle, accompanied by a reasonably accurate prediction of the quenching behaviour. However, In the TRAC calculation the

Page 22: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

18

fuel rapidly dries out again, whereas the thermocouple data show an extended period of quench and a auch greater removal of stored energy. In the peripheral bundles the predictions are considerably worse. The experimental data indicate extensive quenching, with evidence of asymmetry in flows. (As might be expected, more cooling and quenching is seen in the peripheral fuel in the intact hot leg quadrant than in the opposite quadrant). In the TRAC calculations there is essentially no asymmetry and no quenching. The thermal behaviour predicted by the other codes is significantly different from the TIAC predictions. At the very bottom of the core both RELAP and DRUFAN/FLUT show virtually no heat-up. The remaining lower half of the core is reasonably well predicted by DRUFAN/FLUT, while RELAP shows generally too much heat removal and no heat-up after about 5 seconds. At the higher elevations DRDFAN/FLUT shows no top-down cooling and therefore predictions are generally poor, while RELAP tends to show no heat-up phase and predictions are poor in the opposite direction.

Figures 31 and 32 illustrate the differences In behaviour between the central bundle and peripheral bundles, and show typical examples of differences in the code predictions. Overall, the objective of maximising the core temperatures at the beginning of ref111/reflood could be said to have been partially achieved. Most of the central bundle saw little cooling during blowdown, however this is offset by the extensive cooling of the peripheral bundles, [it is difficult to see how the experiment could have been carried out differently in such a way as to reduce further the blowdown cooling].

The refill and reflood phases show a wider divergence of predicted hydraulic behaviours. The accumulator liquid flow period is from about 17 to 46 seconds, and all calculations predict this reasonably well. The experimental data clearly show that there is no direct bypassing of subcooled liquid during this period. The 3-D TRAC calculations correctly predict this. However, all the 1-D calculations predict some degree of direct bypass, In spite of using split downcwr models in most cases. The single downcoaer model In FLUT inevitably predicts direct bypass, but also predicts liquid penetration down the downcomer by having countercurrent flow of steam and subcooled liquid, which is quite unlike the actual situation In the downcoaer. Refill depends on the degree of voiding during blowdown, which varies markedly from calculation to calculation. For example, the mlniaua liquid fractions in the lower plenum and downcoaer range from 0.05 to O.SS and froa 0.15 to 0.4 respectively. It therefore alght be expected that there would be a spread In the predicted onset of reflood. However, this is not the case, and aost calculations predict the onset of reflood within a few seconds of the aeasured tlae of 33 seconds. The early reflood proceeds with a series of fluid oscillations between the core and downcoaer, extending to the intact cold leg. These appear to be driven by a coablnation of core heat transfer and downcomer/cold leg condensation. The code calculations generally predict the oscillations in a qualitative, but not quantitative, sense. One of the reasons why quantitative predictions are not obtained is the poor prediction of core stored energy at the end of blowdown. At about 46 seconds the accumulator Is eapty of liquid and nitrogen begins to flow into the primary circuit. This has the effect of pressurising the cold leg and upper downcoaer, due to the Increased volumetric flow and the reduced condensation rate. TRAC models the nitrogen flow but appears to over-predict the pressure effect somewhat. The other codes do not model Che nitrogen flow, but do fee an effect froa the reduced condensation rate. Most calculations are

Page 23: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

19

qualitatively correct in predicting a large surge of liquid Into the core at this time. However, one of the 1-0 calculations does not predict this surge because the lower plenum is highly voided at this time, and the surge of liquid merely fills it up. The surge of liquid into the core is seen quite clearly in the SPMD responses. The sudden increase in steam production in the core produces a pressure Increase which transiently drives liquid out of the core and up the downcomer. There is a period of Increased loss of fluid from the broken cold leg lasting about 10 seconds and including direct bypassing of subcooled liquid. The TRAC calculation predicts such an increase, but probably over-estimates the amount of liquid lost. The 1-D calculations tend not to predict an increase in break flow. Although the surge of liquid at about 46 seconds cause considerable cooling and quenching, unlike other LOFT LB transients It does not completely quench the core. Thus there follows a period of about 25 seconds of 'classical' reflood under pumped ECCS conditions, easily the most extended in the LOFT LB test series. This phase of the reflood has fairly regular oscillations of fluid between the core and the downcomer. The code calculations generally predict this oscillatory phase, at least qualitatively. However, there is a wide divergence in the predictions of core average liquid content, in the interval 60-80 seconds the predicted average liquid fraction in the core is about 0.60 for TRAC, 0.25 for FLUT and as low as 0.20 for RELAP5/M002. The TRAC calculation also shows a tendency for the vessel liquid content to decrease with time after about 50 seconds.

The thermal behaviour during refill/reflood is, not unexpectedly, closely linked to the hydraulic behaviour. In the experimental data it is clear that there is a period of steam cooling beginning at about 25 seconds which turns over the cladding temperature from its peak of 1260K. Little quenching occurs until about 50 seconds in either the central or peripheral bundles. In the central bundle at the lower levels the cooling is quite well-predicted by all the codes. At mid-elevations TRAC tends to over-predict somewhat, while the other codes are generally good. At 49 inches and above the TRAC temperatures are too high because too little heat was removed during the earlier top-down flow, while most RELAP calculations show temperatures which are too low because there was virtually no heat-up phase earlier. The FLUT calculation showed no top-down cooling in the earlier phase but still shows a quench (from a much higher temperature) at about the right time. The overall quenching behaviour predicted by the codes shows wide variations. The TRAC T B i n correlation does not allow quenching until much later and from much lower temperatures than those seen In the data, [it should be noted that this is qualitatively consistent with the evidence from the NEPTUN facility on the effect of the LOFT external thermocouples on reflood quenching as mentioned In an earlier section of this paper]. The other calculations show a range of apparent quench temperatures ranging from well above the data to cases where no distinct quench temperature can be identified. The information available for the peripheral bundles is more limited. TRAC consistently over-predicts the temperatures and does not show the marked asymmetry seen in the data. The FLUT calculations are similar to those for pins of the same rating placed In the central bundle, and are in reasonable agreement with the data. The Halted RELAP data suggests that no or little heaL-up was predicted In the peripheral bundles due to the excess liquid remaining in the core throughout the transient.

Page 24: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

20

10 GENERAL CONCLUSIONS

The LOFT programme has generated a very large aaount of Information, much of it of a unique nature. It is not practical in a review paper of this type to list here all the detailed conclusions which have arisen, though many of them are implicit in the main text. However, it is felt that it is worthwhile to draw some general conclusions:

1 The worldwide database for Integral large break tests is very limited. The LOFT facility has been a major contributor to this database of high quality data.

2 The 0ECD-L0FT tests LP-02-6 and LP-LB-1 essentially achieved all their defined objectives, thereby adding significantly to the world database.

3 The hydraulic behaviour seen in LOFT is now well-understood. However, there remain uncertainties in our understanding of the core thermal behaviour.

4 Several countries collaborated in the computer code analyses of the tests, mostly using the TRAC and RELAPS codes. These analyses show that the hydraulic behaviour in the loops is reasonably well-predicted, nevertheless there is a considerable divergence in the predicted vessel hydraulics. This in turn leads to a wide divergence in the predicted core thermal behaviour, with none of the analyses showing good agreement with the data.

11 REFERENCES

[l] Adams, J.P.: "Monitoring PWR Reactor Vessel Liquid Level with SPNDs during LOCAs", Proceedings of the Regional Meeting on Advances in Reactor Physics and Core Thermal Hydraulics, Kiamesha Lake, New York, 22-24 September 1982.

[2] Olsen, C S . : "Burst Strength of Zircaloy Tubing with Embedded Thermocouple", LTR LO-14-82-081, April 1982.

[3] Birchley, J.C.: "OECD-LOFT Project Experiment Specification Document for Cold Leg Break Experiment LP-LB-1", OECD L0FT-T-3502, December 1983.

[4] Modro, S.M., Aksan, S N., Berta, V.T., Wahba, A.B: "Review of LOFT Large Break Experiments", NUREG/IA-0028, October 1989, (OECD-LOFT-T-3900 Topical Report, August 1988).

[5] Adams, J.P. and Berta, V.T.: "Response of LOFT SPNDs to Reactor Coolant Density Variations During a LOCA Simulation", Transactions of the American Nuclear Society, 16, 334, Winter Meeting 1980.

[6] Adams, J.P. and Berta, V.T.: "Monitoring Reactor Vessel Liquid Level with a Vertical String of SPNDs", ANS Winter Meeting, San Francisco, California, 1983.

Page 25: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

21

7] Aksan, S.N.: "Evaluation of Analytical Capability to Predict Cladding Quench", EGG-LOFT-5555, August 1982.

8] Lin, J.C.: "Post Test Analysis of LOFT Loss-of-Coolant Experiment L2-3", ECG-LOFT-5075, March 1980.

9] Grush, W.H., Lin, J.C., Linebarger, S.H., Nelson, R.A.Jr.: "Results and Predictions of Scaled Nuclear Large Break Loss of Coolant Experiments", ASME Winter Meeting, Chicago, Illinois, (80-WA/HT-49), 1980.

10] Nelson, R.A.: "Forced-Convective Post-CHF Heat Transfer and Quenching", Transactions of the ASME vol 104, 1982.

11] Adaas, J.P., Condie, K.G., Batt, D.L.: "Quick-Look Report on 0EC0/L0FT Experiment LP-02-6", 0ECD/L0FT-T-3404, October 1983.

12] Coddington, P., Gill, C : "TRAC-PF1/M0D1 Calculations of LOFT Experiment LP-02-6", AEEW - R 2464, AEE Winfrith, August 1987.

13] Aksan, S.N., Toloan, E.L., Nelson, R.A.: "Application of Analytical Capability to Predict Rapid Cladding Cooling and Quench during the Blowdown Phase of a Large Break Loss-of-Coolant Accident", Proceedings of the Second International Topical Meeting on Nuclear Reactor Thermal-Hydraulics (Vol 3), Santa Barbara, California, USA, 11-14 January 1983.

14] Garner, R.W., MacDonald, P.E.: "Power Burst Facility Thermocouple Effects Test Results Report, Test Series TC-1, TC-2 and TC-3", NUREG/CR-2665 (EGG-2190), June 1982.

15] Gottula, R.C., and Good, J.A.: "The Effect of Cladding Surface Thermocouples on the Quench Behaviour of an Electrical Heater Rod", EG&G Idaho Report LO-00-80-115, 1982.

16] Gottula, R.C.: "Effects of Cladding Surface Thermocouples and Electrical Heater Rod Design on Quench Behaviour", NUREG/CR-2691, EGG-2186, February 1984.

17] Gottula, R.C., Carbonneau, M.L., Tolman, E.L.: " An Assessment of LOFT Fuel Rod Quench Behaviour based on Electric Rod Quench Tests", Proceedings of the International Symposium on Fuel Rod Simulators-* Development and Application, Gatlinburg, Tennessee, USA, October 1980, (Conf-801091, May 1981).

18] Aksan, S.N.: "Investigations on Rapid Cladding Cooling and Quench during the Blowdown Phase of a Large Break L0CA using RELAP5/MOD2", Proceedings of the Fourth International Topical Meeting on Nuclear Reactor Thermal-Hydraulics, NURETH-4, Vol 1, pp 214-220, Karlsruhe, FRG, 10-13 October 1989.

19] Yackle, T.R.: "An Assessment of the Influence of Surface Thermocouples on the Behaviour of Nuclear Fuel Rods during a Large Break L0CA", The 8th Water Reactor Safety Research Information Meeting, Gaithersburg, Maryland, USA, October 1980,

Page 26: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

22

Garner, R.W.: "PBF Thermocouple Effects Tests, Test Series TC-4, Quick Look Report", EGG-TFBP-5465, August 1981.

Adams, J.P., Birchley, J.C.: "Quick Look Report on OECD-LOFT Experiment LP-LB-1", OECD L0FT-T-3504, February 1984.

Aksan, S.M.: "A Review of Large Break Loss-of-Coolant Accident Slowdown Quench and Effect of External Thermocouples", Paper presented at the 3rd International Code Assessment and Application Programs« (ICAP), Specialist Meeting, Grenoble, France, 1-4 March 1988.

Aksan, S.N., Analytis, G.Th., LUbbesmeyer, D.: "Switzerland's Code Assessment Activities in Support of the International Code Assessment Programme (ICAP)", Paper presented at the 16th Hater Reactor Safety Information Meeting, Gaithersburg, Maryland, USA, 24-27 October 198t.

Ransom, V.H., et al: "RELAPS/M0D2 Code Manual Volume 1: Code Structure, System Models and Solution Methods", NUREG/CR-4312 (ECG-2396) Rev 1, March 1987.

Brlttain, I., Coney, M.W.E., "TRAC and RELAP5 Code Development within the UK", Paper presented at the 16th Water Reactor Safety Information Meeting, Gaithersburg, Maryland, USA, 24-27 October 1988.

Coddington, P., Gill, CR., " Development of an External Thermocouple Model to Analyse the LOFT Large Break Experiments L2-2, L2-3 and L2-6", Proceedings of the International ENS/ANS Conference on Thermal Reactor Safety (NUCSAFE 88), Avignon, France, 2-7 October 1988.

Gill, C.R., Coddington, P.: "Analysis of Heat Transfer from Fuel Rods with Externally Attached Thermocouples", AEEW - R 2352, AEE Winfrith, May 1988.

Coddington, P.: "Implementation of an External Thermocouple Model into TRAC-PF1/M0D2", AEEW - R 2547, AEE Winfrith, August 1989.

Grutter, R., Stlerli, F., Aksan, S.N., Varadl, G.: "NEPTUN Bundle Reflooding Experiments: Test Facility Description", EIR Report No 386, WUrenlingen, March 1980.

Aksan, S.N.: "NEPTUN Bundle Boil-off and Reflood Experimental Programme Results", The 10th Water Reactor Safety Research Information Meeting, Gaithersburg, Maryland, USA, NUREG/CP-0041, Vol 1 (1982).

Richner, M.: "Comparison of NEPTUN-I and NEPTUN-II Refloodlng Experiments: External Thermocouple Effect", (In German), EIR Internal Report (TM-32-87-06), WUrenlingen, May 1987.

Vltanzac, C , et al: "Blowdown/Reflood Tests With Nuclear Heated Rods (IFA-511.2)", OECD Halden Reactor Project, HPR-248, May 1986.

Johnsen, T., Vltanza, C : "Blowdown/Reflood Tests with SEMISCALE Heaters (IFA-511.3, Data Collection)", OECD Halden Reactor Project HWR-17, May 1981.

Page 27: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

23

[34] Johnsen, T., Vltanza, C : "Results of LOCA Tests with SEMISCALE heaters (IFA-511.5)", OECD Halden Reactor Project HPR-313, May 1984.

[35] Vltanza, C , Johnsen, T.: "Results of Blowdown/Reflood Tests with REBEKA Electric Simulators (IFA-S11.4, July 2)", OECD Halden Reactor project, HWR-85, May 1983.

[36] Ihle, P., Rust, K.: "Reflood Experiments Using Full Length Bundles of Rods with Zircaloy Claddings and Alumina Pellets (Results of the SEFLEX Programme)", Proceeding of Third International Topical Meeting on Reactor Thermal Hydraulics (Vol 1), Newport, Rhode Island, USA, 15-18 October 1985.

[37] Aksan, S.N., Stierli, F., Analytis, G.Th.: "Boil-off Experiments with the EIR-NEPTUN Facility: Analysis and Code Assessment Overview Report", SIR Report No 623, September 1987.

[38] Condie, K.G., et al: "Best Estimate Predictions for the OECD-LOFT Project Large Cold Leg Break Experiment LP-02-6", OECD LOFT-T-3403, September 1983.

[39] Demmie, P.N., et al: Planning Calculations for the OECD LOFT Project Urge Cold Leg Break Experiment LP-02-6", OECD LOFT-T-3406, September 1983.

[40] Brittain, I., and O'Mahoney, R.: "Planning Calculations using TRAC-PD2/M0D1 for the OECD LOFT Project Large Cold Leg Break Experiment LP-LB-1", AEEW - M 2129.

[41] Birchley, J.C., Newman, N.: "Best Estimate Prediction of OECD LOFT Project Large Cold Leg Break Experiment LP-LB-1", OECD LOFT-T-3503, January 1984.

[42] Guntay, S.: A Comparative Analysis on OECD-LOFT Experiment LP-02-6", OECD L0FT-T-3410, November 1987.

[43] Coddington, P.: "Analysis of LOFT Experiment LP-LB-1 Using the TRAC-PF1/M0D1 Code", AEEW - R 2039, January 1986.

[44] Coddington, P.: "OECD-LOFT LP-LB-1 Comparison Report", AEEW - R 2478, February 1989.

[45] Knight, T.D.: "TRAC Analysis of LOFT LP-02-6", OECD LOFTJT-3409, October 1985.

Page 28: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

24

TABLE 2

Chronology of Events for Experiment LP-02-6

Event

Blowdown valves opened

Reactor scrammed

Primary coolant pumps tripped

Cladding temperatures initially deviated from saturation

End of subcooled break flow

Maximum cladding temperature reached (blowdown)

Bottom-up core rewet initiated

Bottom-up core rewet complete

Partial core top-down quench initiated

Pressurizer emptied

Primary coolant pumps disconnected from flywheels

Accumulator Injection initiated

Partial core top-down quench complete

High-pressure injection initiated

Lower plenum refill complete

Low-pressure injection initiated

Maximum cladding temperature reached (reflood)

Accumulator Injection complete

Core quench complete

Core reflood complete

Time Seconds

0

0.1 ± 0.01

0.8 ± 0.01

0.9 ± 0.01

4.0 ± 0.5

4.9 ± 0.2

5.2 ± 0.2

9.1 ± 0.2

14.8 ± 0.2

15.5 ± 0.5

16.5 ± 0.01

17.5 ± 0.5

18.6 ± 0.2

21.8 ± 0.01

30.7 ± 0.2

34.8 1 0.01

41.0 ± 0.2

57 ± 5

56 ± 0.2

59 ± 1

Page 29: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

25

TABLE 3

Chronology of Events for Experiment LP-LB-1

Event

Blowdown valves opened

Reactor scrammed

Primary coolant pumps tripped

Primary coolant pumps disconnected from flywheels

Cladding temperatures initially deviated from saturation

End of subcooled break flow (broken hot leg) (broken cold leg)

Maximum cladding temperature reached (blowdown)

Partial core top-down quench initiated

Pressurizer emptied

Accumulator injection initiated

Partial core top-down quench complete

Maximum cladding temperature reached (reflood)

Lower plenum refill complete

Low-pressure injection initiated

Accumulator empty

Accumulator injection complete

Core reflood complete

Core quench complete

Time Seconds

0

0.13 + 0.01

0.24 ± 0.01

0.63 ± 0.01

0.8 ± 0.1

1.0 ± 0.1 3.5 ± 0.1

12.9 ± 0.5

13 ±0.5

15 ± 1

17.5 ± 0.05

25.5 ± 0.5

26.8 ± 0.5

34.5 ± 0.5

32 ± 1

40 ± 1

46 ± 2

50 ± 2

72 ± 1

Page 30: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

26

Figure 1 Axonometrie Projection of LOFT System Configuration for Large Break LOCA

VVVVvN

LfiPT

TIP .e u«

•r.t .p «•«

-t.TT.n »f oK

.5.Tr.M «P UP

Figure 2 FUR Pressure Vessel and LOFT Elevations

Page 31: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

27

en

1300

1200

1100

1000 -

Figure 3

TIME, SECONDS

LF-L2-6 Cladding Temperacures on Rod SG06

OS

E a. 5 H

10 15 10 IS 30 35 40 4S JO

TIME, SECONDS

1000

100

(00

10

Figure 4 LP-02-6 Clad Temperacure and SPND Fluid Density at 11 inches

s

s 400 £

200

Page 32: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

28

H öS

I td tu

s

1200

1000

800

600

-BOO

-JOOO

ARROWS INOICATE TIME FLOW

400 -5

A-J

-500

10

DIFFERENTIAL PRESSURE (MPa)

Figure 5 LP-02-6 Cladding Pressure Differencial and Ballooning Performance

Id OS

! w 0«

5

U00

»00

1100

1000

MO

•00

700

100

soo

400

• j\_fJ'^ p^ / 11 i n ^ ^ Y .

/ ' ' ^ \ \> ' 1 / ! 45 in .

/ i ! / 1 ! f> ,'' 62 in-

t 1 1 * ( ' I k4 / \

30 in.

'V:

*. \

\

" i *

,'i

i 10 10 «0 SO 70

TINE, SECONDS

Figure 6 LP-LB-1 Cladding Tenperacure* on Rod 5G06

Page 33: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

29

WOO

1200

T r T r

5«5 1000

2 «00 H PL«

s «00

400

Cor« hot loot LP-LB-1

LP-02-0 L2-5

-I L.

•10 0 10 20 30 40 80 «0 70 «0

TIME, SECONDS

Figure 7 Cladding Temperature Comparison

CO V)

400

350

300

250

200

150

100

50

0

-50

• • ' • 1 1 , . . . . , , p

12-5 LP-02-6

h WM \ \

\jiA VA wA

'^/\ l S-«W\ A 4i L *?SSW\ d ml • ^AXw^(HJwJWf

^•'«iJ'Waw 3H

' i i i ' i

-

) -

••

k .i *rA:

VAF-WU ircpy«

10 15 20 25 30

TIME, SECONDS

Figure 8 Coolant Mass Flow Ratas in the Intact Loop Cold Leg During Experiments L2-5 and LP-02-6

Page 34: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

30

1100

Cd CS

B

!

TIME, SECONDS

Figure 9 Peak Cladding Temperatures for LOFT Experiments L2-2, L2-3 and LP-02-6

400

-600

TIME, SECONDS

Figure 10 Difference Between Intact and Broken Cold Leg Mass Flows for LOFT Experiment LP-02-6

Page 35: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

31

1100

H

ei

a,

s

—i I i — Nuclear rod with gap

RE35KA rod with gap _ —- Solid U02 rod ZR clad

Solid U02 rod SS clad • — Semiscale heater roc — — LTS? dala

\ *

8 10 12 14 TIME, SECONDS

Figure 11 Calculated Nuclear Fuel, Cartridge-Type and Solid-Type Heater Rod Responses for LTSF Single Rod Blowdown Quench Tests

v. a H

S (9 Z

I u o z M o o

3 4 M H M

250

200

150

100

50

" I 1 - -T" " i LXSF data (Semiscale rod)

— Calculated nuclear roo . response for LTSF • ' test conditions y

y

/ 0

/ /

/ J /

/ *

/ /

i Nuclear data from PBF Test TC-3

^^* ^ '

^ 9 * * ^

i i f

Initial cladding temperature = 1025 K

Inlet quality = 0 Pressure = 7.0 MPa

8 INLET FLOW VELOCITY, m/s

Figure 12 Comparison of Measured SEMISCALE Heater Rod Initial Cooling Rates from LTSF Single Rod Quench Tests with Calculated Nuclear Fuel Rod Initial Cooling Rates

Page 36: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

32

liO

TIME, SECONDS

Figure 13 Blowdown Phase Quench Behaviour for a REBEXA Heater Rod (from LTSF Bundle Tests) and a Nuclear Fuel Rod (fron PBF-TC tests) with Cladding External Thermocouples

1000

900

800

700

600

900

InUl YcUdly • I M / f W»( auoBiv - 0Z O - 0J5 x 10* Ib/hr II*

FEBA ooo ooo ooo

SESEKA

_l_

0.0 2.5 5.0 7.5

TIME, SECONDS

10.0 12.5 15.0

Figure 14 LTSF Bundle Test: Conparison of REBEKA and FEBA Electrical Heater Rods Under Blowdown Quench Conditions

Page 37: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

33

OUTWIT SCEMCXT MCASMMC JUNCTOR

Figure 15 LOFT Fuel Rod Cladding External Thermocouple Atcachmenc

1200 1 1 i I i i Rod with surlaca and Internal thermoeouplts

Surface mennocsuple Internal ttieonocouole

Oensäonieter

Bare rod wich Internal

thermocouple

\

\

Pressurt a 7 MPa Flooding rale a 1.8 m/j Cladding iimoeraiure = 1175 K Inlet duality a 0

1 : • !_ 14 16 18 20

TIME, SECONDS

Figure 16 LTSF Single Rod Quench Test Thermocouple Response for a SEMISCALE Heater Rod

Page 38: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

34

1100

600 --"'.!„

Bor« TC4-G rod Inil. TC-t-E rod

500 10 U

TIME, SECONDS

Figure 17 Comparison Between Embedded and External Cladding Thermocouple Responses During Blowdown Quench in FBF-TC Tests

tu

s

2000

1500

1000

500

- TC-5007-27 -- TE-6D07-OZ7

•50 60 100 160

TIME, SECONDS

Figure 18 Fuel Centreline and Cladding Temperature at the 27 inch Elevation on Fuel Rod 5D07 during Experiment LP-02-6

Page 39: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

35

2000

1500

et

I H tu

s 1000 -

500

TC-5007-27 (with daddmg TCJ TC-6O0S-27 (wHhogl cladding TO TC-5OT0-27 (without cladding TO

J_ _l_

10 20 30 40

TIME, SECONDS

Figure 19 Fuel Centreline Temperatures at the 27 inch Elevation during Experiment LP-02-6

fid

s

2000

1500

1000

500

•60

1 i i TC-5F09-027 (cantartin* without TO TC 6007-027 (centarbie with TO TE-5D07-027 (cladding TO

_L

50 100 150

TIME, SECONDS

200 250 300

Figure 20 Fuel Rod Centreline Tenperature With and and Without Cladding Surface Thermocouple and Cladding Temperature Measurement During Experiment LP-LB-1 at the 27 inch Elevation

Page 40: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

36

2 M

s

2000

1S00

1000

500

-SO

1 1 1 TC-SJ07-43.8 (cmuaifcw without TO TC-Ct07-43J (cMtwfcM with TO TE-5107-43." KJmtfna TO

_1_

50 100

TIME. SECONDS

150 200 250 300

Figure 21 Fuel Rod Centreline Temperature With and Without Cladding Outer Surface Thermocouple and Cladding Temperature Measurement at the 43.8 inch Elevation

1100

Figure 22 TAU Calculation of LOFT Experiment LP-02-6 ac 21 inch Elevation

Page 41: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

37

s u 2 M

• z o HI H <

LEV

Cd

• J

%

70

60

SO

40

30

20

10

0

-

I I

us » 1

1 1

" T T - ( Ml

Art M / - /

H r -* w ^ j / # ™ *

j —J *—* H>

/ 1 — * HI

1 >~i

1 H

it—i

1H

-

_

_

-

-

10 20 30 40 SO 60 70

TIME THAT LIQUID LEVEL DETECTORS INDICATE ECC WATER ENTRY

Figure 23 Comparison of Liquid Level Versus Time for Final Reflood for LOFT Experiments L2-3 and L2-5

1000 1 1 LOFT Rod. laval 4 (extern*! TO LOFT Rod. Iivil 6 (external TO NEPTUN Rod. level 4 (embedded TO NEPTUN Rod. level 5 (embedded TO-

100 150

TIME, SECONDS

250

Figur« 24 Cladding Temperature* During NEPTUN-I and NEPTUN-II Experiments, With and Without External Cladding Surface Thermocouple*

Page 42: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

38

o o

PE

RA

TU

RE

,

1000

SOO

800

700

600

500

400

300

Level 2 Level 4

_i_

Lev.l 6

50 100 150

TIME, SECONDS

200 250

Figure 25 Cladding Temperatures During NEPTUN-I and NEPTUN-II Experiments, With and Without External Cladding Surface Thermocouples

w at

e a,

s

' Nucletr rod

' Electric rod

j^--Riflood -

_L 50 100 150 200 250 300 350 400

TIME, SECONDS

Figure 26 Cladding Temperatures of Electrical Heater Rods and Nuclear Fuel Rods During Quench Experiments in Halden

Page 43: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

39

noo

u o Id BS

I

100 200 300 400

TIME, SECONDS »oo

(00

400

200

(00

0 « 0 200 300 400 (00

TIME. SECONDS

Figure 27 Cladding Temperatures Measured at Two Different Axial Levels in FEBA and REBEKA Rod Bundles

I

suu

600

400

ZOO

n

i i

f U m Canttr LOFT rod* — . , _ j _

LOFT Ihwmocoupl*

,

i i

i i

Af *4U4££

V

-

29 Inch axial «Itviiion

i i

200 4 0 0 000

TIME, SECONDS

800 1000

Figure 28 Cladding Temperatures Measured in NEPTUN Experiment With External Cladding Surface (LOFT) Thermocouple end Vlth Embedded Thermocouples

Page 44: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

40

u OS 1 w ft.

§

»JO«

U N

1000

20 40 (0 10

TIME, SECONDS

100

Figure 29 LF-02-6 Comparison TRAC Calculations Centre Bundle Fuel Surface Temperatures at 31 inch Elevation

1200

I a.

20 40 U 10

TIME, SECONDS

120

Figure 30 LP-02-6 Cooparison RELAP and DRUFAN Calculations Centre Bundle Fuel Surface Temperatures at 31 inch Elevation

Page 45: Paul Scherrer Institut€¦ · AEA-TRS-1003 PSI Report 72 OECD-LOFT LARGE BREAK LOCA EXPERIMENTS: PHENOMENOLOGY AND COMPUTER CODE ANALYSES I Brittain* S N Aksan+ # AEA Technology

41

Cd OS

1 Cd

5

40 «0

TIME, SECONDS

10 100

Figure 31 LP-LB-1 Comparison. Fuel Rod Cladding Temperatures i n Central Fuel Bundle (5) Temperature at Elevation 43.8 inches

et

I

1200

1100

1000

M0

too

TOO

(00

SOO

400

1 i i i

} f\ f «V"*-% \ \ M \ III / I -• -̂ i n / $*>* [| A

falA? *^''*>'~~ '—1

'

— i i i * — — - > . —

20 40 M

TIME, SECONDS

10 leo

Figure 32 LF-LB-1 Comparison. Fuel Rod Cladding Temperatures in Feripheral Fuel Bundles 6, 4 Temperature at Elevation 39 inches