Numerical Analysis of the CIRCE-HERO PLOFA Scenarios

The present work deals with simulations carried out at the University of Pisa by using the System Thermal Hydraulics code RELAP5/Mod3.3 to support the experimental campaign conducted at the ENEA (Energia Nucleare ed Energie Alternative) Brasimone Research Centre on the CIRColazione Eutettico—Heavy liquid mEtal pRessurized water cOoled tubes (CIRCE-HERO) facility. CIRCE is an integral effect pool type facility dedicated to the study of innovative nuclear systems and cooled by heavy liquid metal, while HERO is a heat exchanger heavy liquid metal/ pressurized cooling water system hosted inside the CIRCE facility. Beside the H2020 project Multi-Purpose Hybrid Research Reactor for High-Tech Applications (MYRRHA) Research and Transmutation Endeavour (MYRTE), a series of experiments were performed with the CIRCE-HERO facility, for both nominal steady-state settings and accidental scenarios. In this framework, the RELAP5/Mod3.3 code was used to simulate the experimental tests assessing the heat losses of the facility and analyzing the thermal hydraulics phenomena occurring during the postulated Protected Loss Of Flow Accident (PLOFA). The modified version Mod. 3.3 of the source code RELAP5 was developed by the University of Pisa to include the updated thermo–physical properties and convective heat transfer correlations suitable for heavy liquid metals. After reproducing the facility through an accurate nodalization, boundary conditions were applied according to the experiments. Then, the PLOFA scenarios were reproduced by implementing the information obtained by the experimental campaign. Sensitivity analyses of the main parameters affecting the thermofluidynamics of the Lead-Bismuth Eutectic (LBE) were carried out. In the simulated scenario, the LBE mass flow rate strongly depends on the injected argon flow time trend. The numerical results are in agreement with the experimental data, however further investigations are planned to analyze the complex phenomena involved.


Introduction
The Lead-Cooled Fast Reactor (LFR) is one of the six nuclear reactor concepts selected by the Generation IV International Forum (GIF) to sustain the R&D of the next generation of nuclear energy systems [1]. The LFRs are cooled by Heavy Liquid Metals (HLM) such as lead or Lead-Bismuth Eutectic (LBE), employing a fast-neutron spectrum and a closed fuel cycle for efficient conversion of fertile uranium and management of actinides. The main advantages of using HLMs as coolant are [2]: • high thermal inertia; • very low absorption and scattering cross sections; • high boiling points allowing a low-pressure operation;

Experimental Facility
CIRCE is an integral effect pool-type facility, designed and constructed at the ENEA Brasimone research center. It consists of a cylindrical main vessel of AISI316L partially filled with nearly 70 tons of molten LBE covered with argon at about 0.2 bar [7,8]. The vessel is externally equipped with electrical heating cables installed at the bottom and on the lateral surface and operating in a temperature range of 200 • C-400 • C. The pool vessel parameters are summarized in Table 1. The HERO-loop test section is installed inside CIRCE from the top of the main vessel through a coupling flange; a schematic view of the CIRCE-HERO facility is reported in Figure 1. HERO-loop test section is installed inside CIRCE from the top of the main vessel through a coupling flange; a schematic view of the CIRCE-HERO facility is reported in Figure 1.  A Steam Generator Bayonet Tube (SGBT), the heat sink of the loop for heat removal. It is composed of 7 double wall bayonet tubes. Each tube represents in scale 1:1 the tube foreseen for the SGBT to be adopted in the Advanced Lead Fast Reactor European Demonstrator-super (ALFRED) [7]. This configuration guarantees the separation between the primary side (LBE) and the secondary side (steam-water) decreasing the probability of LBE-water interaction and facilitating the monitoring of leakages by pressurizing the annular separation region. The HERO-loop consists of:  [7]. This configuration guarantees the separation between the primary side (LBE) and the secondary side (steam-water) decreasing the probability of LBE-water interaction and facilitating the monitoring of leakages by pressurizing the annular separation region.
• An argon injection device is located at the inlet section of the riser. It circulates the LBE by imposing a fixed argon mass flow rate. • A dead volume, placed above the fitting volume and insulated from the pool, which encloses and maintains the wirings and electrical connections of the FPS.
A feeding conduit equipped with an LBE Venturi flow meter and a differential pressure transmitter allows the measuring of the LBE flowing upwards through the FPS region where it is heated up by the 37 electrical pins. Then the LBE passes through the fitting volume and the riser (at this point the flow is enhanced by the argon gas injection at the top of the fitting volume) reaching the separator. From the separator, the LBE enters the shell side (primary side) of the HERO-SGBT flowing downwards. Here, the LBE is cooled by the feedwater flowing in the 7 Bayonet Tubes (BTs) (1 central tube and 6 lateral tubes allocated in a hexagonal arrangement) to then be discharged in the bottom zone of the pool. Further details regarding the facility are reported in [7,8].

PLOFA Test
The test aims at analyzing the mixed convection and thermal stratification phenomena that take place into a pool filled with heavy liquid metal during the transition from forced convection (nominal conditions) to natural circulation (accidental scenario). The experiment involves a PLOFA scenario where a sudden degradation in the heat transfer capability occurs due to a sharp decrease in the flow rate. The test described in this work was selected to represent the reference test for the benchmark calculations of SESAME project. The test starts when the steady-state condition is approximately reached with an FPS power of 356 kW and forced circulation condition in the HERO-loop. The forced circulation of the LBE was provided by buoyancy effects and enhanced by argon gas injection. During the initial stage of the experiment, the argon flow rate was 2.75 Nl/s and the resulting LBE mass flow rate was about 34 Kg/s. Inside HERO SGBT the heat was removed by means of water flowing in the secondary side with a total mass flow rate of 0.264 Kg/s, an inlet temperature of about 336 • C and pressure equal to 175 bar. The transient was triggered after about half an hour, and the FPS power was reduced to a value of about 20 kW according to the curve shown in Figure 2.
imposing a fixed argon mass flow rate. • A dead volume, placed above the fitting volume and insulated from the pool, which encloses and maintains the wirings and electrical connections of the FPS.
A feeding conduit equipped with an LBE Venturi flow meter and a differential pressure transmitter allows the measuring of the LBE flowing upwards through the FPS region where it is heated up by the 37 electrical pins. Then the LBE passes through the fitting volume and the riser (at this point the flow is enhanced by the argon gas injection at the top of the fitting volume) reaching the separator. From the separator, the LBE enters the shell side (primary side) of the HERO-SGBT flowing downwards. Here, the LBE is cooled by the feedwater flowing in the 7 Bayonet Tubes (BTs) (1 central tube and 6 lateral tubes allocated in a hexagonal arrangement) to then be discharged in the bottom zone of the pool. Further details regarding the facility are reported in [7,8].

PLOFA Test
The test aims at analyzing the mixed convection and thermal stratification phenomena that take place into a pool filled with heavy liquid metal during the transition from forced convection (nominal conditions) to natural circulation (accidental scenario). The experiment involves a PLOFA scenario where a sudden degradation in the heat transfer capability occurs due to a sharp decrease in the flow rate. The test described in this work was selected to represent the reference test for the benchmark calculations of SESAME project. The test starts when the steady-state condition is approximately reached with an FPS power of 356 kW and forced circulation condition in the HERO-loop. The forced circulation of the LBE was provided by buoyancy effects and enhanced by argon gas injection. During the initial stage of the experiment, the argon flow rate was 2.75 Nl/s and the resulting LBE mass flow rate was about 34 kg/s. Inside HERO SGBT the heat was removed by means of water flowing in the secondary side with a total mass flow rate of 0.264 kg/s, an inlet temperature of about 336 °C and pressure equal to 175 bar. The transient was triggered after about half an hour, and the FPS power was reduced to a value of about 20 kW according to the curve shown in Figure 2.
The gas injection was decreased sharply from 2.75 Nl/s to 0 generating the corresponding LBE mass flow rate shown in Figure 3. After the transient, the LBE flow was induced only by the natural circulation in the loop and was oscillating around a value of 6 kg/s. The feedwater was reduced to 30% of the steady-state value in a time ramp of 2 s as reported inFigure 4.

RELAP5/Mod.3.3 Nodalization
The RELAP5 nodalization of the CIRCE-HERO facility is Figure 5. It is worth mentioning the presence of RELAP5 default structures called ANNULUS-used in this case to simulate the steam generator components-which are similar to pipes, but require that all the water must be in a film on the wall (no drops) in the annular mist flow regime (ANNULUS structure [5,6]). Details concerning the implemented components are reported in the Appendix A. The nodalization of the HERO-loop is designed according to the following configuration:

•
The feeding conduit is modelled with PIPE 20 which is connected to the downstream mixing zone of the FPS (PIPE 50) by means of a Single Junction, (SNGLJUN [5,6].).

•
The FPS region is modelled with PIPES 50, 60, and 70 where the active length of the FPS (1000 mm active region) is modelled with PIPE 60. The FPS is modelled using a heat structure inside the pipe 60 with a source term representing the power generated by the 37 electrical pins. The upstream mixing zone of the FPS (PIPE 70) is connected to the fitting volume.

•
The fitting volume is represented via PIPE 80 and BRANCH 90 [5,6] and it is connected to the riser by means of a SNGLJUN.

•
At the inlet of the riser (PIPE 130) argon gas is injected to enhance the LBE circulation. This argon injection system is modelled with Time Dependent Volume 3, (TMDPVOL [5,6]).

•
After the riser, the separator is modelled by BRANCH 132. Furthermore, it is connected to the argon cover gas volume (TMDPVOL 151) through the BRANCH 150.

•
Connected to the separator there is the PIPE 172 which represents the primary side (LBE side) of the SGBT. LBE is released into the pool by means of PIPE 174.

RELAP5/Mod.3.3 Nodalization
The RELAP5 nodalization of the CIRCE-HERO facility is Figure 5. It is worth mentioning the presence of RELAP5 default structures called ANNULUS-used in this case to simulate the steam generator components-which are similar to pipes, but require that all the water must be in a film on the wall (no drops) in the annular mist flow regime (ANNULUS structure [5,6]). Details concerning the implemented components are reported in the Appendix A. The nodalization of the HERO-loop is designed according to the following configuration: • The feeding conduit is modelled with PIPE 20 which is connected to the downstream mixing zone of the FPS (PIPE 50) by means of a Single Junction, (SNGLJUN [5,6].).

•
The FPS region is modelled with PIPES 50, 60, and 70 where the active length of the FPS (1000 mm active region) is modelled with PIPE 60. The FPS is modelled using a heat structure inside the pipe 60 with a source term representing the power generated by the 37 electrical pins. The upstream mixing zone of the FPS (PIPE 70) is connected to the fitting volume.

•
The fitting volume is represented via PIPE 80 and BRANCH 90 [5,6] and it is connected to the riser by means of a SNGLJUN.

•
At the inlet of the riser (PIPE 130) argon gas is injected to enhance the LBE circulation. This argon injection system is modelled with Time Dependent Volume 3, (TMDPVOL [5,6]).

•
After the riser, the separator is modelled by BRANCH 132. Furthermore, it is connected to the argon cover gas volume (TMDPVOL 151) through the BRANCH 150.

RELAP5/Mod.3.3 Nodalization
The RELAP5 nodalization of the CIRCE-HERO facility is Figure 5. It is worth mentioning the presence of RELAP5 default structures called ANNULUS-used in this case to simulate the steam generator components-which are similar to pipes, but require that all the water must be in a film on the wall (no drops) in the annular mist flow regime (ANNULUS structure [5,6]). Details concerning the implemented components are reported in the Appendix A. The nodalization of the HERO-loop is designed according to the following configuration: • The feeding conduit is modelled with PIPE 20 which is connected to the downstream mixing zone of the FPS (PIPE 50) by means of a Single Junction, (SNGLJUN [5,6].).

•
The FPS region is modelled with PIPES 50, 60, and 70 where the active length of the FPS (1000 mm active region) is modelled with PIPE 60. The FPS is modelled using a heat structure inside the pipe 60 with a source term representing the power generated by the 37 electrical pins. The upstream mixing zone of the FPS (PIPE 70) is connected to the fitting volume.

•
The fitting volume is represented via PIPE 80 and BRANCH 90 [5,6] and it is connected to the riser by means of a SNGLJUN.

•
At the inlet of the riser (PIPE 130) argon gas is injected to enhance the LBE circulation. This argon injection system is modelled with Time Dependent Volume 3, (TMDPVOL [5,6]).

•
After the riser, the separator is modelled by BRANCH 132. Furthermore, it is connected to the argon cover gas volume (TMDPVOL 151) through the BRANCH 150.

Results and Discussions
In this section, the numerical results obtained by using the nodalization described in Section 3 and applying the boundary conditions of the experimental test reported in Section 2 are presented. According to the boundary conditions of the experimental test, a transient of 3600 s was simulated by imposing the PLOFA starting time at 1800 s. The FPS power was imposed according to Figure 2.
Similarly to the experiment, the required LBE mass flow rate is obtained by imposing the Argon mass flow rate (TDPV3 and J004 in Figure 5), as explained in Section 2. In Figure 6 the experimental LBE mass flow rate is reported with black diamonds, while the numerical value obtained by applying the nominal conditions is reported with a dotted light blue line. Due to the approximation introduced by the simplified numerical geometry and the assumptions adopted in terms of concentrated pressure drops, the obtained LBE mass flow rate overestimates the experimental value. The primary side of the SGBT is connected to the secondary side (steam-water) by means of a heat structure. The secondary side is modelled as follow: • The feedwater manifold is modelled with TMDPVOL 410.

•
The feedwater flow rate is supplied to the system via the TMDPJUN (Time Dependent Junction, [5,6]) connecting TMDPVOL 410 to the BRANCH 418, which models the feedwater inlet header to the BTs. Concerning the heat structures, the active length of the FPS region is modelled as a power source embedded inside the thermal structure of PIPE 60. The power source is defined by a power time trend card and the Ushakov heat transfer correlation [10] is used for the triangular lattice fuel bundle region. Furthermore, the following heat structures and wall boundary conditions were applied:

•
Heat structures between the primary and the secondary side of the In order to simulate the axial thermal conduction inside the pool, heat structures were applied according to the arrows reported in Figure 5 and connecting PIPES 203-213, 213-223, 223-253 and PIPE 253 with BRANCH 260.

Results and Discussions
In this section, the numerical results obtained by using the nodalization described in Section 3 and applying the boundary conditions of the experimental test reported in Section 2 are presented. According to the boundary conditions of the experimental test, a transient of 3600 s was simulated by imposing the PLOFA starting time at 1800 s. The FPS power was imposed according to Figure 2.
Similarly to the experiment, the required LBE mass flow rate is obtained by imposing the Argon mass flow rate (TDPV3 and J004 in Figure 5), as explained in Section 2. In Figure 6 the experimental LBE mass flow rate is reported with black diamonds, while the numerical value obtained by applying the nominal conditions is reported with a dotted light blue line. Due to the approximation introduced by the simplified numerical geometry and the assumptions adopted in terms of concentrated pressure drops, the obtained LBE mass flow rate overestimates the experimental value. Therefore, in order to assure the same initial conditions of the experiments, the LBE mass flow was reduced by following two different approaches. As a first approach, the pressure drops inside the riser (PIPE 130) were increased, the resulting LBE mass flow is reported with the dash-dotted orange line. As a second approach, the LBE mass flow was directly imposed as an input parameter according to the dashed purple line shown in Figure 6. In these simulations, mass flow rates are applied by using the component TMDPJUN, which joins a TMDPVOL to the corresponding hydrodynamic volume. The TMDPVOL imposes the desired temperature, pressure and static quality of the injected fluid, while the TMDPJUN applies the chosen mass flow rate (liquid and/or gas phase) as a function of the time. Thus, under both forced circulation during normal operation and natural circulation after the PLOFA, the thermal behavior of the system is simulated. The initial temperature of both the whole HERO-loop and the bottom part of the pool (PIPE 253-BRANCH 260) was set to 420 • C. According to the experiment, a temperature equal to 475 • C was initially applied in the rest of the pool [7,8]. For the secondary side of the SGBT the inlet temperature of the water was set equal to 336 • C and the related mass flow rate was implemented according to   After the PLOFA, the low LBE mass flow could generate a different flow path, avoiding LBE hitting the lower plenum after leaving the HERO-loop and flowing directly to the inlet of the FPS. This behavior joined to the applied combination "LBE mass flow rate-FPS power", which determines the same temperature at the end of the simulation. In this work results in terms of the temperatures reached inside the FPS and in the outlet region of HERO-SGBT are also reported and analyzed. In Figure 7 the temperature trends measured during the experiment at the inlet and outlet of the FPS are reported and compared to the related numerical values obtained by implementing the different trends of LBE mass flow rate reported in Figure 6. In the same way, Figure 8 shows the comparison between experimental and numerical results of the HERO-SGBT outlet temperature. In the legends of both Figures 7 and 8 the experimental data are nominated with "Exp.", the numerical analysis carried out under nominal condition with "Nom.", while numerical cases with increased pressure drops and LBE mass flow rate imposed as input are reported with "PressDrop" and "ImpInp", respectively.    According to the obtained numerical results, by approaching the experimental LBE mass flow rate the FPS temperatures evaluated numerically become closer to the related experimental results. Before the PLOFA starting time, the inlet and outlet temperature slightly increase in the experiments while both decrease in the numerical analyses. Furthermore, compared to the nominal case, when a lower LBE mass flow rate is applied a more effective cooling is performed by the water mass flow rate, which is kept constant. As a result, a lower LBE temperature is evaluated. After the transient, the experimental and numerical analyses behave similarly from a qualitative point of view for: both the nominal case and the numerical analysis with increased pressure drops in the riser. The temperature decreases in both cases, however, a difference of about 10% is found at the end of the simulation between experimental and numerical results. On the contrary, the third numerical analysis, named 'ImpInt', generates a temperature increase after the PLOFA, converging to the FPS temperature of the experiment, with a difference of about 0.1% at the end of the simulation. This behavior is determined by the increase in the temperature at the HERO-loop outlet (see Figure 7) due to the combination "LBE mass flow rate-FPS power" occurring soon after the PLOFA. Indeed, LBE mass flow rate decreases up to 6 Kg/s according to the experiments, while in the other two simulations the value stays around 10 Kg/s. Under this condition, the FPS power (20 kW) easily warms up the LBE increasing the temperature at the HERO-loop outlet and as a consequence at the FPS inlet. In all cases, the differences occurring before the PLOFA starting time could be due to the continuous reduction of the FPS inlet temperature, which should be instead almost constant. This reduction is generated by the same temperature reduction of the LBE at the outlet of HERO-SGBT (see Figure 8). This entails that when LBE comes out the SGBT it enters almost directly into the FPS maintaining the same temperature. Looking at the CFD (Computational Fluid Dynamics) results reported in the literature [17], a longer and more complex path is made by LBE from the outlet of HERO-SGBT to the inlet of the FPS. Indeed, according to [17], soon after entering in the pool, the LBE hits the lower plenum rising the pool in the opposite direction up to the outlet zone of the SGBT to then conclude the path flowing down to the inlet of the FPS. During this route, the LBE has more time to warm up, keeping its FPS inlet temperature almost constant at 420 • C. Therefore, as a further study, the volume of the lower plenum was intentionally increased by one order of magnitude to let the LBE warm sufficiently.   After the PLOFA, the low LBE mass flow could generate a different flow path, avoiding LBE hitting the lower plenum after leaving the HERO-loop and flowing directly to the inlet of the FPS. This behavior joined to the applied combination "LBE mass flow rate-FPS power", which determines the same temperature at the end of the simulation. After the PLOFA, the low LBE mass flow could generate a different flow path, avoiding LBE hitting the lower plenum after leaving the HERO-loop and flowing directly to the inlet of the FPS. This behavior joined to the applied combination "LBE mass flow rate-FPS power", which determines the same temperature at the end of the simulation. Figure 9 shows the comparison between the experimental FPS temperatures and the related numerical values obtained by increasing the volume of the pool lower plenum. Here, the case with LBE mass flow rate imposed as input parameter was considered. A good agreement between numerical and experimental results during the entire analysis is obtained for both quantitative and qualitative study. This reveals the importance of a coupled CFD-STH analysis to include the exact path made by the LBE and evaluate the resulting temperature. While FPS experimental temperature decreases after PLOFA, numerical results remain almost constant, generating a temperature difference of about 3% at 3600 s. This behavior is generated by the approximation introduced by increasing the volume of the lower plenum. Figure 10 shows the variation of the temperature trend generated in the lower plenum by increasing its volume by one order of magnitude starting from about 1 m 3 . When the volume increases the temperature remains almost constant generating a different outlet temperature in the FPS. path made by the LBE and evaluate the resulting temperature. While FPS experimental temperature decreases after PLOFA, numerical results remain almost constant, generating a temperature difference of about 3% at 3600 s. This behavior is generated by the approximation introduced by increasing the volume of the lower plenum. Figure 10 shows the variation of the temperature trend generated in the lower plenum by increasing its volume by one order of magnitude starting from about 1 m 3 . When the volume increases the temperature remains almost constant generating a different outlet temperature in the FPS.  The facility represents a complicated thermal-hydraulic system which involves two-phase flow, non-condensable gas and liquid metal flow conditions. In this framework, a stand-alone approach could not be sufficient for addressing all the involved phenomena, while a coupled CFD-STH application may represent the best option to obtain a suitable representation of the thermo-hydraulic behavior of this system. Thus, further CFD-STH coupled analyses were recently performed obtaining promising results in predicting the observed phenomena, details are reported in [19]. difference of about 3% at 3600 s. This behavior is generated by the approximation introduced by increasing the volume of the lower plenum. Figure 10 shows the variation of the temperature trend generated in the lower plenum by increasing its volume by one order of magnitude starting from about 1 m 3 . When the volume increases the temperature remains almost constant generating a different outlet temperature in the FPS.  The facility represents a complicated thermal-hydraulic system which involves two-phase flow, non-condensable gas and liquid metal flow conditions. In this framework, a stand-alone approach could not be sufficient for addressing all the involved phenomena, while a coupled CFD-STH application may represent the best option to obtain a suitable representation of the thermo-hydraulic behavior of this system. Thus, further CFD-STH coupled analyses were recently performed obtaining promising results in predicting the observed phenomena, details are reported in [19]. The facility represents a complicated thermal-hydraulic system which involves two-phase flow, non-condensable gas and liquid metal flow conditions. In this framework, a stand-alone approach could not be sufficient for addressing all the involved phenomena, while a coupled CFD-STH application may represent the best option to obtain a suitable representation of the thermo-hydraulic behavior of this system. Thus, further CFD-STH coupled analyses were recently performed obtaining promising results in predicting the observed phenomena, details are reported in [19].

Conclusions
This work aimed at evaluating the phenomena occurring in HLM pool-type systems as a result of a Protected Loss Of Flow Accident (PLOFA) by using the last updated version of the STH code RELAP5/mod3.3 made by the University of Pisa (UNIPI) to take into account liquid metals as working fluids.
For this purpose, an experimental campaign launched in the framework of different projects involved in HORIZON 2020 was considered. The numerical analyses were performed at UNIPI, while the experimental campaign was carried out at the ENEA Brasimone Research Centre.
The CIRCE-HERO test facility was simulated through a full detailed model by adopting the most convenient nodalization. Standalone simulations were performed to simulate the PLOFA scenario by setting up appropriate boundary conditions to reproduce the experimental trends. The PLOFA scenario is characterized by a sudden degradation in the heat transfer capability due to the sharp decrease in the 1.
Time-dependent volume (TMDPVOL) is a hydrodynamic boundary volume used to specify a fluid state boundary condition. Quantities such as pressure, liquid temperature, vapour temperature, void fraction, and quality can be set as boundary conditions. The TMDPVOL provides the user with a mechanism for absolutely defining the fluid condition at a point in the model. The user should consider that a TMDPVOL acts as an infinite fluid source or sink. Its conditions remain unchanged (or vary) as requested, but are invariant with inflow or outflow [5,6].

2.
A Time-dependent junction (TMDPJUN) component allows the user to impose a flow boundary condition on a model. It is possible to specify the flow condition as either a volumetric or mass flow rate. An example of this capability is the specification of an injection flow as a function of the coolant system pressure. TMDPJUN capabilities include varying the flow condition in any manner and as a function of any problem variable the user desires [5,6].

3.
A Single-Volume component (SNGLVOL) is the basic hydrodynamic cell unit. The flow area, length, and volume of the cell must be defined by the user to describe the geometry in the input file. The input flow area determines the flow velocity, the input length affects the calculated frictional pressure drop, and the input volume contributes to the overall fluid system volume [5,6].

4.
A Single-Junction component (SNGLJUN) is the basic hydrodynamic flow unit. Each junction is used to connect hydrodynamic volumes, which are characterized by an inlet and outlet face [5,6].

5.
A Pipe (PIPE) component is a hydrodynamic volume component. As a function of the length it can be set as series combination of single-volumes joined by single-junctions. The advantage of the pipe over the separate single components is primarily one of input efficiency, reducing significantly the number of data cards to include. By definition, the pipe component has only internal junctions associated with it. Any connections to the ends of a pipe must be made with external junctions. Furthermore, it is possible to connect external junctions to any face of internal pipe cells and any face of pipe cells at the ends of a pipe [5,6]. 6.
An Annulus component (ANNULUS) is identical to the pipe component except that an annular flow regime map is used. An annulus must be specified as a vertical component [5,6].

7.
A Branch component (BRANCH) is a single-volume component that may have single-junctions appended [X].
The heat structures are used to represent metal structures such as vessel walls, steam generator tubes, fuel rods, and reactor vessel internals in a facility. Each heat structure is defined to connect thermally a "left" and a "right" side. Each side of a heat structure may be connected to at most one hydrodynamic volume. However, more than one heat structure may be connected to the same hydrodynamic volume. The geometry of the heat structure can be chosen as well as the desired boundary condition (surface temperature or surface heat flux). Adiabatic, conductive and convective boundary conditions can be activated. The default convection and boiling correlations were derived mainly based on data from internal vertical pipe flow [5,6].