Next Article in Journal
Solar Radiation Prediction Using a Novel Hybrid Model of ARMA and NARX
Next Article in Special Issue
Model Development and Transient Analysis of the HCPB BB BOP DEMO Configuration Using the Apros System Code
Previous Article in Journal
Accurate Demand Forecasting: A Flexible and Balanced Electric Power Production Big Data Virtualization Based on Photovoltaic Power Plant
Previous Article in Special Issue
MHD R&D Activities for Liquid Metal Blankets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Simulation of an Out-Vessel Loss of Coolant from the Breeder Primary Loop Due to Large Rupture of Tubes in a Primary Heat Exchanger in the DEMO WCLL Concept

by
Francesco Galleni
1,*,
Marigrazia Moscardini
1,
Andrea Pucciarelli
1,
Maria Teresa Porfiri
2 and
Nicola Forgione
1
1
Laboratory of Numerical Simulations for Nuclear Thermal-Hydraulics, Department of Civil and Industrial Engineering, University of Pisa, 56122 Pisa, Italy
2
ENEA-CR Frascati, Dipartimento Fusione e Tecnologie per la Sicurezza Nucleare, 00044 Roma, Italy
*
Author to whom correspondence should be addressed.
Energies 2021, 14(21), 6916; https://doi.org/10.3390/en14216916
Submission received: 17 August 2021 / Revised: 13 October 2021 / Accepted: 18 October 2021 / Published: 21 October 2021
(This article belongs to the Special Issue Thermal-Hydraulics in Nuclear Fusion Technology: R&D and Applications)

Abstract

:
This work presents a thermohydraulic analysis of a postulated accident involving the rupture of the breeder primary cooling loop inside a heat exchanger (once through steam generator). After the detection of the loss of pressure inside the primary loop, a plasma shutdown is actuated with a consequent plasma disruption, isolation of the secondary loop, and shutoff of the pumps in the primary; no other safety counteractions are postulated. The objective of the work is to analyze the pressurization of the primary and secondary sides to show that the accidental overpressure in the two sides of the steam generators is safely accommodated. Furthermore, the effect of the plasma disruption on the FW, in terms of temperatures, should be analyzed. Lastly, the time transients of the pressures and temperatures in the HX and BB for a time span of up to 36 h should be obtained to assess the effect of the decay heat over a long period. A full nodalization of the OTSG was realized together with a simplified nodalization of the whole PHTS BB loop. The code utilized was MELCOR for fusion version 1.8.6. The accident was simulated by activating a flow path which directly connected one section of the primary with the parallel section of the secondary side. It is shown here that the pressures and the temperatures inside the whole PHTS system remain below the safety thresholds for the whole transient.

1. Introduction

In the Roadmap to Fusion Electricity Horizon 2020, the European DEMOnstration Fusion Reactor (DEMO) is expected to be a nuclear fusion power plant with the aim of showing the feasibility of the production of electrical power through the conversion of around 2 GWth, generated continuously by the fusion reaction. The operational sequence is a pulsed operation, which consists of 11 pulses per day; each pulse comprises a burn time of 2 h (power pulse period; 100% of fusion power) and a dwell time of 10 min (1% of fusion power generated due to the decay heat) [1,2].
The work presented here aims to investigate the consequences of a loss of coolant accident in the primary heat transport system (PHTS) of the water-cooled lithium lead (WCLL) breeding blanket (BB) DEMO concept [2,3,4]. In particular, the objective of this work is the simulation of a rupture of pipes of the primary system into the PCS (power conversion system or secondary system) inside the once through steam generator (OTSG) which acts as interface between the two systems [5,6,7]. In normal operations, the energy transferred from the breeding zone (BZ) and from the first wall (FW) to PCS through steam generators is used to produce the main steam at condition suitable to feed the steam turbine.
The postulated accident foresees the rupture of the primary cooling loop inside the OTSG with a discharge of primary coolant in the secondary side.
In this work, a detailed nodalization of the primary and secondary system of the OTSG is presented together with a simplified nodalization of the whole BZ loop. A connection between the two sides of the OTSG is created in order to simulate the LOCA scenario, which comprises a mitigated plasma disruption, and the results of the consequent transients are analyzed and discussed.
The nodalization of the involved systems was created with MELCOR for fusion 1.8.6, and all the simulations were performed with the same code. The MELCOR code has been under development for fusion applications for many decades [8,9]. Chiefly because of its capabilities of assessing thermal–hydraulic transients of fusion reactor systems and the transport of radionuclides, MELCOR was chosen, together with other system thermal hydraulics codes, to be used to perform safety analyses for the ITER project and, consequently, for the DEMO project [10]. Especially in the last 10 years, MELCOR was widely used for accident analysis related to ITER reactor safety and in the preliminary design phase of the DEMO reactor ([11,12,13,14,15,16,17], among many others).

2. DEMO WCLL-BB PHTS System

2.1. General Parameters and Power Data

The thermodynamic cycle used as reference for the design of the WCLL BB PHTS DEMO reactor was mainly based on parameters similar to pressurized water reactors: the coolant was water at 15.5 MPa with inlet and outlet temperatures equal to 295 °C and 328 °C, respectively [3,7,8,9]. The main working parameters of the BB PHTS are given in Table 1, whilst the power parameters are reported in Table 2. Furthermore, Table 3 collects the mass flow rates of the FW and the BZ cooling systems. All the parameters and characteristics summarized here can be found in several technical and published studies [18].
The WCLL DCD BOP BB PHTS constitutes two independent primary systems:
  • The BZ primary system (BZ PHTS);
  • The FW primary system (FW PHTS).
The main components of the WCLL BB PHTS are indicated in the 3D CAD model (see Figure 1). Table 4 reviews the system-relevant data [2,7,19,20]. Details of the system architecture can be found in [18].
During normal operations, in pulse time, BZ and FW PHTSs transfer power to the PCS, through two once through steam generators (OTSGs) per system (i.e., four SGs in total). The BZ PHTS power is 1483 MW, and the FW PHTS power is 439.8 MW. A total of six main coolant pumps (MCPs) are installed to allow the circulation of the primary coolant (four pumps for BZ PHTS and two for FW PHTS). Each PHTS is equipped with a pressurizer (PRZ).

2.2. BZ Once through Steam Generator

In the BZ Primary System, each OTSG removes 742 MWth of thermal power, with a mass flow rate of coolant equal to 3831.2 kg/s. On the secondary side, the water is assumed to be at a pressure of 6.4 MPa, and the feedwater coolant inlet temperature is expected to be at 238 °C. The feedwater flow rate is imposed at 404 kg/s, in order to produce the same amount of superheated steam at 299 °C. It is important to underline that the OTSG is still in a design phase [21,22]. A simplified scheme of a generic OTSG is shown in Figure 2 and the main characteristics are summarized in Table 5.
The selected OTSG is characterized by 7569 tubes, with a length of 12.987 m.

2.3. BB and PHTS Integration

The WCLL DEMO 2018 blanket is composed of 16 sectors (22.5°) in the toroidal direction [4,23]. Each sector consists of three poloidal segments in the OB (outboard) blanket and two poloidal segments in the IB (inboard) blanket. The geometrical features and the inventory of the WCLL-BB are provided in Table 6.
The inner structure of each segment comprises a stack of about 100 toroidal–radial so-called breeding units (BUs). The BU is a sort of elementary cell repeated along the poloidal direction of each segment. It features the integrated FW–SW complex, the BZ, and the corresponding part of manifolds and BSS.
The coolant feeding pipes deliver cold water to FW and BZ systems, whilst the outlet pipes collect hot water and distribute it to the PHTS.
Considering the FW and BZ systems, the volume of water is estimated taking into account the BUs (assuming a rough value of 100 for each sector), the manifolds, and the feeding pipes. Furthermore, the velocities in the inlet feeding pipes and distributors were calculated assuming a density of 737.1 kg/m3 (15.5 MPa and 295 °C), while a density of 657.5 kg/m3 (15.5 MPa and 328 °C) was taken for the outlet pipes and collectors. The main geometrical data are provided in Table 6.

3. Postulated Accidental Scenario

In this section, the specifications of the accident—such as the event sequence and the main objectives—are reported and described in Table 7, Table 8, Table 9 and Table 10.

4. Nodalization of the WCLL-PHTS

Figure 3 shows the whole nodalization of the BZ loop. Since the loop is symmetrical with respect to the two OTSGs, only half of the loop is represented, without losing consistency. All the geometrical characteristics of the components of the nodalization were retrieved from the references presented above and summarized in the previous sections of this document. A list of all the main characteristics can be found in Table 7 and Table 8.
The whole nodalization can be roughly divided into two sections: the main one, more complex, which is the nodalization of the BZ OTSG, and the simpler one, which is the nodalization of the BB volumes; the two regions are highlighted in Figure 3, whilst a focus on the BZ OTSG is shown in Figure 4. The cold and hot rings provide a connection between the two regions. The pressurizer is directly connected to the hot ring.
The OTSG is partitioned in the primary and secondary side. Each side is divided into several control volumes (CVH in MELCOR), which correspond to the hemispherical top and bottom, as well as to the sections created by the support plates. The division into different sections is required to achieve a temperature gradient as close as possible to the actual one. All the control volumes are connected by flow paths (FL). CVH 10 and 11 (see Figure 4) are dummy volumes used only to impose the BCs to the secondary side at normal operations. The CHV and FL geometrical characteristics are summarized in Table 11 and Table 12, respectively.
In normal conditions, the primary and secondary side are connected only by heat structures (HSs), which allows heat transfer between the volumes but do not permit any passage of fluids. Between the primary and secondary side, there are a total of 17 equal heat structures. Each HS has a cylindrical geometry divided into four nodes with a total thickness of 0.8 mm and a height of 0.7825 m, which corresponds to the height of the connected volumes. The multiplicity of the HS is equal to 7569, which means that the HS is equivalent to 7569 heat structures of the same type, for a total surface area of 4795 m2. All HSs are made of stainless steel, and a convective condition is imposed on both sides of the HS.
The thermal power of the OTSG is simulated by means of two heat structures connected with the inboard and outboard sections of the BB regions. These two HSs provide 592 MW of thermal power to the outboard section (CVH 15) and 149 MW to the inboard section (CVH 14).
The accident is simulated by activating dummy flow paths which connect one volume of the primary side with the corresponding volume of the secondary side. The section area of this FL is equal to the rupture area (0.00312 m2 for nine tubes on the OTSG primary side). These FLs are located at different heights to allow sensitivity analysis on the location of the rupture.

4.1. Boundary Conditions

In order to impose the required primary coolant mass flow rate of 3831.2 kg/s, a fluid velocity was imposed at FL12, since MELCOR does not allow to impose directly the mass flowrate. Furthermore, only for the initial steady condition, before the onset of the accident, the temperature on the bottom of the OTSG (CVH 5) was imposed to be 295 °C and the pressure was set at 15.5 MPa.
As mentioned in the previous section, two dummy CVHs (10 and 11) were used to impose the BCs at the secondary side. A pressure of 6.4 MPa and a feedwater coolant inlet temperature of 238 °C were imposed. The feedwater flow rate is kept at 404 kg/s by means of imposing an equivalent constant velocity FL, in order to produce the same amount of superheated steam at 299 °C.
The velocity (mass flow rate) imposed at FL12 was set to 0.0 at the FPS since the pumps were shut down due to the loss of offsite power. Furthermore, at the FPS, the FLs connecting the secondary side of the OTSG to the secondary loop are closed (thus acting as isolation valves), also interrupting the feed water flow in the secondary side. At the current stage, the design of the secondary loop is not completed.

4.2. Decay Heat

The decay heat was imposed using volumetric heat generation from heat structures connected with the outboard and inboard CVH. The total values and time trends, calculated using the most up-to-date volumetric nuclear heating data given in [24], are reported in Table 13.
However, it is important to note that the values reported in [24] did refer to a different configuration of the WCLL, namely, the one which was composed of modules and, therefore, had a different material ratio; in this work, the total nuclear heating was recalculated scaling the values with an estimation of the total volume obtained from the most recent configuration of the inboard and outboard segments, and this may represent a source of error. Furthermore, it is difficult to estimate the role of the FW in cooling the decay heat of the breeding unit under this particular scenario, since—as shown in the results—the FW loop is not damaged by the plasma disruption, and it can be assumed still functional. As a first approach in this work, only the decay nuclear heating relative to the BM back wall and the manifolds is assumed to have an effect on the damaged BZ cooling loops.

5. Numerical Analyses

In the simulation, the system was initially left to run for 1000 s, to reach a stable steady state in operational conditions, and then the FL connecting the volumes involved in the rupture was opened. Three seconds after the low-pressure signal, i.e., when the pressure in the loop fell below 12.4 MPa, the FPS was activated; plasma disruption occurred, the pumps were stopped, and the secondary loop was isolated.
In total, three different cases were run, at three different rupture locations (top, middle, and bottom). Table 14 summarizes the different cases.
The event sequence of the accident is reported in Table 15.

5.1. Pressure Evolution

Figure 5 and Figure 6 show the pressure trend after the rupture, at different time scales. The pressure started decreasing and then showed a sharp peak when the low-pressure signal was reached. This was due to the rapid shutdown of the pump and the isolation of the loop. The same behavior was seen in all the cases, but it is worth noting that, for Case 1, the pressure drop was slower and the peak remained significantly lower than the other two cases. This was due to the different rupture conditions; in Case 1, the primary side discharged in a region of the secondary, which was filled with superheated vapor, whereas, in Cases 2 and 3, it was discharged in regions containing liquid water. Figure 7 shows the time trend of the vapor quality calculated in CVH 85 for Case 1.
Figure 8 and Figure 9 show the pressure trend at the rupture locations of both the primary and the secondary side. In the secondary side of the OTSG, no pressure peaks were predicted, and the pressure increased to an equilibrium value by around 100 s in all the cases. Then, 200 s after the rupture, a steady state was reached, with the pressure lying between 12 and 13 MPa.
All pressures remained well below the threshold of 18 MPa, with the pressure peak for Cases 2 and 3 having a value of about 15 MPa, and that for Case 1 having a value of about 12.5 MPa.

5.2. Temperature Evolution in the BB

Figure 10, Figure 11 and Figure 12 show the time evolution of the temperature at the inlets and outlets of the BB, for the whole transient up to 36 h. Because of the effect of the decay heat, the temperatures did not reach a steady state, but tended to increase. The faster temperature increase after 104 s was due to the onset of weak natural circulation in the loop, with the water flow even changing direction several times in Cases 2 and 3, as shown in Figure 13. However, the increase could be considered slow throughout the whole transient, and all temperatures remained well below acceptable values.
As mentioned above, the difference between Case 1 and Cases 2 and 3 was due to a different height of the rupture in the primary side of OTSG for the three cases, with a jump of around 5 m between each case. The transient for Case 1 showed a different behavior from the very beginning, because, in this case, the primary discharged, through the rupture, into a region of the secondary which was filled with superheated vapor and not liquid water, as in Cases 2 and 3. This early difference consequently led to a radically different transient of the pressure and general behavior of the loop, compared with the other two cases.

6. Conclusions

This work presented a thermohydraulic analysis of a postulated accident involving the rupture of a breeder primary loop inside an HX (once through steam generator). After the detection of the loss of pressure inside the primary cooling loop, an FPS is actuated with a consequent plasma disruption, isolation of the secondary loop, and shutoff of the pumps in the primary; no other safety counteractions are postulated.
The objective of the work was to analyze the pressurization of the primary and secondary sides to show that the accidental overpressure in the two sides of the steam generators is safely accommodated. Furthermore, the effect of the plasma disruption on the FW, in terms of temperatures, should be analyzed. Lastly, the time transients of the pressures and temperatures in the HX and BB for a time span of up to 36 h should be obtained to assess the effect of the decay heat over a long period.
A full nodalization of the OTSG was realized together with a simplified nodalization of the whole PHTS BB loop. The code utilized was MELCOR for fusion version 1.8.6. The accident was simulated by activating a flow path which directly connected one section of the primary with the parallel section of the secondary side.
It was shown here that the pressures and the temperatures inside the whole PHTS system remained below the safety thresholds for the whole transient.
The only caveat of this analysis is due to the calculation of the effect of the decay heat, since the tables used for the estimation of this parameter referred to an old version of the WCLL-BB, as no other references are available on this issue.
It is believed that, for a fully comprehensive evaluation of the impact of the decay heat, a coupled calculation should be realized, connecting a nodalization of the FW PHTS loop and the PbLi loop with the nodalization of the BZ PHTS loop presented here.

Author Contributions

Conceptualization, F.G., M.M., A.P. and M.T.P.; methodology, F.G.; validation, F.G.; formal analysis, F.G.; investigation, F.G.; resources, N.F.; data curation, F.G.; writing—original draft preparation, F.G.; writing—review and editing, F.G., M.M., A.P., M.T.P. and N.F.; visualization, F.G.; supervision, M.T.P. and N.F.; project administration, M.T.P. and N.F.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Euratom research and training program 2014–2018 and 2019–2020, grant number 633053.

Acknowledgments

This work was carried out within the framework of the EUROfusion Consortium and received funding from the Euratom research and training program 2014–2018 and 2019–2020 under grant agreement No 633053. The views and opinions expressed herein do not necessarily reflect those of the European Commission.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ACPActivate corrosion product
BBBreeder blanket
BCBoundary condition
BOPBalance of plant
BSSBack supporting structure
BZBreeder zone
CADComputer=aided design
CVHControl volume hydrodynamics
DHRDecay heat removal
DCDDirect coupling design
DEMODEMOnstration power plant
FLFlow path
FPSFast plasma shutdown
FWFirst wall
HT packageHeat structure package
IB Inboard
LOCALoss of coolant accident
OB Outboard
OTSGOnce through steam generator
PCSPower conversion system
PHTSPrimary heat transport system
PbLiLead–lithium
PRZPressurizer
VVVacuum vessel
VVPSSVacuum vessel pressure suppression system
WCLLWater-cooled lithium lead
WPSAEWork Package Safety and Environment

References

  1. Federici, G.; Bachmann, C.; Barucca, L.; Biel, W.; Boccaccini, L.; Brown, R.; Bustreo, C.; Ciattaglia, S.; Cismondi, F.; Coleman, M.; et al. DEMO design activity in Europe: Progress and updates. Fusion Eng. Des. 2018, 136, 729–741. [Google Scholar] [CrossRef]
  2. Barucca, L.; Bubelis, E.; Ciattaglia, S.; D’Alessandro, A.; Del Nevo, A.; Giannetti, F.; Hering, W.; Lorusso, P.; Martelli, E.; Moscato, I.; et al. Pre-conceptual design of EU DEMO balance of plant systems: Objectives and challenges. Fusion Eng. Des. 2021, 169, 112504. [Google Scholar] [CrossRef]
  3. Narcisi, V.; Ciurluini, C.; Giannetti, F.; Del Nevo, A. WCLL BB Phts Ddd (Direct Coupling Option with Small ESS). 2020. Report No.: IDM Ref. EFDA_D_2NURWJ. Available online: https://idm.euro-fusion.org/ (accessed on 13 October 2021).
  4. Del Nevo, A.; Arena, P.; Caruso, G.; Chiovaro, P.; Di Maio, P.A.; Eboli, M.; Edemetti, F.; Forgione, N.; Forte, R.; Froio, A.; et al. Recent progress in developing a feasible and integrated conceptual design of the WCLL BB in EUROfusion project. Fusion Eng. Des. 2019, 146, 1805–1809. [Google Scholar] [CrossRef] [Green Version]
  5. Zaupa, M.; Palma, M.D.; Del Nevo, A.; Moscato, I.; Tarallo, A.; Barucca, L. Preliminary Thermo-Mechanical Design of the Once Through Steam Generator and Molten Salt Intermediate Heat Exchanger for EU DEMO. IEEE Trans. Plasma Sci. 2020, 48, 1726–1732. [Google Scholar] [CrossRef]
  6. Martelli, E.; Giannetti, F.; Ciurluini, C.; Caruso, G.; Del Nevo, A. Thermal-hydraulic modeling and analyses of the water-cooled EU DEMO using RELAP5 system code. Fusion Eng. Des. 2019, 146, 1121–1125. [Google Scholar] [CrossRef]
  7. Ciurluini, C.; Giannetti, F.; Martelli, E.; Del Nevo, A.; Barucca, L.; Caruso, G. Analysis of the thermal-hydraulic behavior of the EU-DEMO WCLL breeding blanket cooling systems during a loss of flow accident. Fusion Eng. Des. 2021, 164, 112206. [Google Scholar] [CrossRef]
  8. Merrill, B.; Moore, R.; Polkinghorne, S.; Petti, D. Modifications to the MELCOR code for application in fusion accident analyses. Fusion Eng. Des. 2000, 51–52, 555–563. [Google Scholar] [CrossRef]
  9. Merrill, B.J. Recent Updates to the MELCOR 1.8. 2 Code for ITER Applications; Idaho National Laboratory (INL): Idaho Falls, ID, USA, 2007. [Google Scholar]
  10. Merrill, B.J.; Humrickhouse, P.; Moore, R.L. A recent version of MELCOR for fusion safety applications. Fusion Eng. Des. 2010, 85, 1479–1483. [Google Scholar] [CrossRef]
  11. D’Ovidio, G.; Martín-Fuertes, F. Accident analysis with MELCOR-fusion code for DONES lithium loop and accelerator. Fusion Eng. Des. 2019, 146, 473–477. [Google Scholar] [CrossRef]
  12. Grief, A.; Owen, S.; Murgatroyd, J.; Panayotov, D.; Merrill, B.; Humrickhouse, P.; Saunders, C. Qualification of MELCOR and RELAP5 models for EU HCLL TBS accident analyses. Fusion Eng. Des. 2017, 124, 1165–1170. [Google Scholar] [CrossRef]
  13. Panayotov, D.; Grief, A.; Merrill, B.J.; Humrickhouse, P.; Trow, M.; Dillistone, M.; Murgatroyd, J.T.; Owen, S.; Poitevin, Y.; Peers, K.; et al. Methodology for accident analyses of fusion breeder blankets and its application to helium-cooled pebble bed blanket. Fusion Eng. Des. 2015, 109–111, 1574–1580. [Google Scholar] [CrossRef]
  14. Nakamura, M.; Tobita, K.; Someya, Y.; Utoh, H.; Sakamoto, Y.; Gulden, W. Safety research on fusion DEMO in Japan: Toward development of safety strategy of a water-cooled DEMO. Fusion Eng. Des. 2016, 109–111, 1417–1421. [Google Scholar] [CrossRef]
  15. Pescarini, M.; Mascari, F.; Mostacci, D.; De Rosa, F.; Lombardo, C.; Giannetti, F. Analysis of unmitigated large break loss of coolant accidents using MELCOR code. J. Phys. Conf. Ser. 2017, 923, 012009. [Google Scholar] [CrossRef] [Green Version]
  16. D’Onorio, M.; Giannetti, F.; Porfiri, M.T.; Caruso, G. Preliminary safety analysis of an in-vessel LOCA for the EU-DEMO WCLL blanket concept. Fusion Eng. Des. 2020, 155, 111560. [Google Scholar] [CrossRef]
  17. D’Onorio, M.; Giannetti, F.; Porfiri, M.T.; Caruso, G. Preliminary sensitivity analysis for an ex-vessel LOCA without plasma shutdown for the EU DEMO WCLL blanket concept. Fusion Eng. Des. 2020, 158, 111745. [Google Scholar] [CrossRef]
  18. Narcisi, V.; Giannetti, F.; Del Nevo, A. WCLL BB PHTS Architecture Description & BOM (Direct Coupling Option with Small ESS). Report No.: IDM Ref. EFDA_D_2PC2N9. Available online: https://idm.euro-fusion.org/ (accessed on 13 October 2021).
  19. Tarallo, A. WCLL BB PHTS CAD Model Description (Direct Coupling Option with Small ESS). Report No.: IDM Ref. EFDA_D_2P9QFM. 2020. Available online: https://idm.euro-fusion.org/ (accessed on 13 October 2021).
  20. Tarallo, A. WCLL BB PHTS Drawings (Direct Coupling Option with Small ESS). Report No.: IDM Ref. EFDA_D_2P95DT. 2020. Available online: https://idm.euro-fusion.org/ (accessed on 13 October 2021).
  21. Tarallo, A. Preliminary Mechanical Design and Verification of WCLL BB PHTS Once through Steam Generator (OTSG). Report No.: IDM Ref EFDA_D_2MT7PL. 2020. Available online: https://idm.euro-fusion.org/ (accessed on 13 October 2021).
  22. Tarallo, A. WCLL OTSG CAD Model. Report No.: IDM Ref. UID 2NJQM8. 2020. Available online: https://idm.euro-fusion.org/ (accessed on 13 October 2021).
  23. Catanzaro, I.; Arena, P.; Basile, S.; Bongiovì, G.; Chiovaro, P.; Del Nevo, A.; Di Maio, P.A.; Forte, R.; Maione, I.A.; Vallone, E. Structural assessment of the EU-DEMO WCLL Central Outboard Blanket segment under normal and off-normal operating conditions. Fusion Eng. Des. 2021, 167, 112350. [Google Scholar] [CrossRef]
  24. Porfiri, M.T.; Mazzini, G. DEMO BB Safety Data List (SDL). Report No.: IDM Ref. EFDA_D_2MF8KU. 2018. Available online: https://idm.euro-fusion.org/ (accessed on 13 October 2021).
Figure 1. Overview of the DEMO WCLL BB PHTS.
Figure 1. Overview of the DEMO WCLL BB PHTS.
Energies 14 06916 g001
Figure 2. Conceptual scheme for once through steam generator [21].
Figure 2. Conceptual scheme for once through steam generator [21].
Energies 14 06916 g002
Figure 3. MELCOR nodalization of PHTS loop.
Figure 3. MELCOR nodalization of PHTS loop.
Energies 14 06916 g003
Figure 4. Detail of BZ OTSG nodalization.
Figure 4. Detail of BZ OTSG nodalization.
Energies 14 06916 g004
Figure 5. Pressure trend after the rupture: primary side, 1000–1030 s.
Figure 5. Pressure trend after the rupture: primary side, 1000–1030 s.
Energies 14 06916 g005
Figure 6. Pressure trend after the rupture: primary side.
Figure 6. Pressure trend after the rupture: primary side.
Energies 14 06916 g006
Figure 7. Vapor quality in volume 85 for Case 1.
Figure 7. Vapor quality in volume 85 for Case 1.
Energies 14 06916 g007
Figure 8. Pressure trend after the rupture: primary and secondary sides, 1000–1200 s.
Figure 8. Pressure trend after the rupture: primary and secondary sides, 1000–1200 s.
Energies 14 06916 g008
Figure 9. Pressure trend after the rupture: primary and secondary sides.
Figure 9. Pressure trend after the rupture: primary and secondary sides.
Energies 14 06916 g009
Figure 10. Temperature evolution in the BB: Case 3.
Figure 10. Temperature evolution in the BB: Case 3.
Energies 14 06916 g010
Figure 11. Temperature evolution in the BB: Case 2.
Figure 11. Temperature evolution in the BB: Case 2.
Energies 14 06916 g011
Figure 12. Temperature evolution in the BB: Case 1.
Figure 12. Temperature evolution in the BB: Case 1.
Energies 14 06916 g012
Figure 13. Water mass flow rate in the loop: 104–105 s.
Figure 13. Water mass flow rate in the loop: 104–105 s.
Energies 14 06916 g013
Table 1. WCLL DCD BOP BB cooling system parameters [3].
Table 1. WCLL DCD BOP BB cooling system parameters [3].
DescriptionUnitsParameters
Typology of coolant Water
PressureMPa15.5
Temperature range°C295–328
Coolant density (average)kg/m3701.3
Design pressureMPa17.8
Design temperature°C345
Total flow ratekg/s9936.0
FW flow rate
Inboard (IB) blanket
Outboard (OB) blanket
kg/s17.8
35.5
BZ flow rate
IB blanket
OB blanket
kg/s48.2
127.5
Table 2. DEMO and WCLL DCD BOP BB power balance [3].
Table 2. DEMO and WCLL DCD BOP BB power balance [3].
DescriptionUnitWCLL 2018 Design
Total nuclear heatingMW1650.3
Total FW Heat FluxMW272.7
Neutron Wall LoadMW167.0
Total FW powerMW439.8
Total BZ powerMW1483.2
Total powerMW1923.2
Table 3. WCLL DCD BOP BB power and coolant flow rates (TD 295–328 °C) [3].
Table 3. WCLL DCD BOP BB power and coolant flow rates (TD 295–328 °C) [3].
DescriptionPower (MW)Mass Flow Rate (kg/s)
----TotalFWBZTotalFWBZ
IB segment12.7803.4499.33166.017.848.2
IB sector25.5606.89818.662132.035.696.4
Total
IB sectors (16)
409.960110.368298.5922112.0569.61542.4
OB segment31.5456.86724.678163.035.5127.5
OB sector94.63520.60174.034489.0106.5382.5
Total
OB sectors (16)
1514.160329.6161184.547824.1704.6120.
Total reactor1923.120439.9841483.1369936.02273.67662.4
Table 4. System-relevant parameters.
Table 4. System-relevant parameters.
WCLL DCD BOP BB PHTS Design and Operating ParametersBZFW
General
Thermal power (MW)1483439.8
Operating pressure (MPa)15.515.5
Reactor vessel inlet temperature (°C)295295
Reactor vessel outlet temperature (°C)328328
Overall volume (m3)~563~159
Overall PHTS piping length (km)~3.2~3.7
Number of loops22
Loop data
Piping length (km)~1.25~1.52
Hot/cold manifolds per loop and size8/8
DN-150, DN-200, DN-350
8/8
DN-100, DN-125, DN-200
Hot/cold legs per loop and size1/2
DN-850, DN-650
1/1
DN-500, DN-500
Hot/cold ring header per loop and size1/2
DN-650
1/1
DN-350
Pump
Number of pumps per loop21
TypeCentrifugal, vertical single-stage (RSR)Centrifugal, vertical single-stage (RSR)
Effective pump power to coolant (MW)3.031.79
Pressurizer
Number of units11
Total volume (m3)101.432.5
Liquid volume (m3)44.816.3
Heat Exchanger (Steam Generator)
Number of units22
Steam generator power (MWth/unit)742219.9
TypeOTSGOTSG
Heat transfer area (m2/unit)49031423
Feedwater temperature (°C)238238
Exit steam pressure (MPa)6.46.4
Steam flow per SG (kg/s)404119.9
Flow rate per SG (kg/s)3831.21136
Table 5. Main features of an OTSG.
Table 5. Main features of an OTSG.
SG PowerMWth742
Primary side pressureMPa15.5
Primary side water Tin°C328
Primary side water Tout°C295
Secondary side pressureMPa6.41
Secondary side water Tin°C238
Secondary side water Tout°C299
No. of tubes--7569
Tube ODmm15.88
Tube thicknessmm0.864
Tube lengthm12.987
Tubesheet lattice--Square
Tubesheet p/D--1.28
Heat transfer aream24903
V water tubesm320
Dext vesselm2.9
Table 6. WCLL BB In-Vessel inventory.
Table 6. WCLL BB In-Vessel inventory.
DescriptionPipe SizeOD
(mm)
Thick.
(mm)
MF Total
(kg/s)
v
(m/s)
L
(m)
L Tot
(m)
V
(m3)
H2O Mass Inventory (kg)
Inlet IB FW distributorDN-100114.38.835.66.617.6281.22.071522.21
Inlet IB FW LIBS/RIBSDN-100114.38.817.83.35.8186.51.371009.43
Outlet IB FW LIBS/RIBSDN-100114.38.817.83.76.4203.91.50984.62
Outlet IB FW collectorDN-100114.38.835.67.417.6281.22.071357.82
Inlet IB BZ distributorDN-200219.117.596.44.917.5280.37.465500.35
Inlet IB BZ LIBS/RIBSDN-150168.312.748.24.14.0128.12.051513.83
Outlet IB BZ LIBS/RIBSDN-150168.312.748.24.64.7149.12.391572.12
Outlet IB BZ collectorDN-200219.117.596.45.517.5279.97.454898.97
Inlet OB FW distributorDN-200219.117.5106.55.411.3180.84.813548.13
Inlet OB FW COBSDN-125139.711.035.54.47.9126.41.381013.89
Inlet OB FW LOBS/ROBSDN-125139.711.035.54.49.2293.93.202357.06
Outlet OB FW COBSDN-125139.711.035.55.08.3133.11.45952.47
Outlet OB FW LOBS/ROBSDN-125139.711.035.55.09.4301.63.282157.42
Outlet OB FW collectorDN-200219.117.5106.56.111.5183.34.883208.51
Inlet OB BZ distributorDN-350355.628.0382.47.410.6170.312.008847.80
Inlet OB BZ COBSDN-200219.117.5127.56.57.0111.52.972187.86
Inlet OB BZ LOBS/ROBSDN-200219.117.5127.56.57.4236.16.294633.48
Outlet OB BZ COBSDN-200219.117.5127.57.37.7122.73.272147.62
Outlet OB BZ LOBS/ROBSDN-200219.117.5127.57.38.4267.97.134688.31
Outlet OB BZ collectorDN-350355.628.0382.48.211.0176.312.438172.15
Total 4094.189.4462,274.06
Table 7. Accidental scenario general specifications.
Table 7. Accidental scenario general specifications.
ParameterSpecification
Name of eventLBO3 out-vessel loss of coolant from the breeder primary loop due to large rupture of tubes in a primary HX in WCLL concept
CategoryAccident
ObjectivesShow that accidental overpressure in the secondary loop of PHTS will be safely accommodated.
Show that post-accident cooling of the decay heat removal system of the VV is sufficient to remove decay heat during the 32 h of the offsite loss of power.
Show that radioactive releases, if any, are adequately confined.
Scope of analysisIntegrated breeder blanket PHTS thermohydraulic analysis.
VV, VVPSS(W) pressure transient analysis if melting occurs due to mitigated plasma disruption.
Analysis of ACPs and tritium transport in the containment volumes as a consequence of the plasma disruption, if melting occurs.
Acceptance criteriaMaximum Eurofer temperature < 550 °C
Confinement integrity: BB module safety assessment pressure < 18 MPa
Confinement integrity: VV safety assessment pressure < 200 kPa
Table 8. Accidental scenario main events sequence.
Table 8. Accidental scenario main events sequence.
ParameterSpecification
Definition of initiating eventBreak in the primary cooling loop side (break size 0.0028 m2, corresponding to the rupture of 9 tubes of the primary side) inside the HX toward the secondary side.
Possible transient sequence Fusion power is terminated by loss of pressure (−20% of operating pressure) in the BB cooling loop inside the HX. The initiating event is followed by a fast plasma shutdown (FPS) actuated 3 s from the low signal, which leads to a mitigated plasma disruption for 5 ms.
The disruption could cause failure of the FW cooling pipes in a BB module if temperature melting is reached in FW Eurofer. In such a case, the break flow area to be considered is reported in [SDL19], chapter 2.1.
The VV decay heat removal cooling loops will cool down the in-vessel components post accident.
Ingress of coolant and radioactive inventories (tritium, dust, and suspended products) will be mobilized.
The rupture discs toward the VVPSS(W) open upon reaching the set VV pressure point.
Mobilized radioactivity is transported into the VVPSS(W).
After the coolant inventory is lost, the FW/breeder blanket modules will be cooled by steam convection and thermal conduction/radiation to the VV.
The DV components are accounted for in this analysis as heat structures at the initial temperature and without cooling during the transient.
Aggravating failures None.
Loss of power A loss of offsite power occurs at the same time of the plasma disruption.
Table 9. Accidental scenario system assumptions.
Table 9. Accidental scenario system assumptions.
ParameterSpecification
Process system assumptionsVV decay heat removal (DHR) will remove the decay heat. The temperature of the VV is maintained at 40 °C by DHR.
Safety systems assumptionsVV pressure limit is 0.2 MPa.
VVPSS(W) rupture discs open upon reaching the VV pressure set point.
Source termTritium, dust in the VV, ACP products in the PHTS, tritium and activated products in the breeder materials and/or purge gas.
Table 10. Expected results from the analysis of the accidental scenario.
Table 10. Expected results from the analysis of the accidental scenario.
Parameter Specification
General The output locations, parameters, and time trace should show the results of safety analyses to support objectives and purposes.
The time span should cover until the transients are stabilized.
FW temperature Transient curve for FW and BB module temperature of the failed loop.
Confinement response Pressures and atmosphere temperature in BB cooling loop, HX.
Pressures and atmosphere temperature in VV and VVPSS(W) only in case of FW Eurofer structure melting.
Cooling system conditions Water break flow versus time in HX and in VV and VVPSS(W) only in case of FW Eurofer structure melting.
Water inventory inside the affected FW/breeder blanket cooling system.
Radioactive transport Transient curves for tritium, dust, and sputtering concentration (airborne, deposited) in VVPSS(W) if FW structure melts.
Bookkeeping of mobilized tritium, dust, and sputtering products if FW structure melts.
Table 11. CVH number and characteristics.
Table 11. CVH number and characteristics.
Nodalization
CVH
Volume Water
(m3)
Height/Length
(m)
Section
(m2)
Description
Primary side
13.221.154.18Top hemisphere
20.720.611.18Top tubesheet
31–4615.360.78251.18Shroud tubes primary
40.720.611.18Bottom tubesheet
53.221.154.18Bottom hemisphere
Secondary side
612.657.631.66Shell boiler
71–8634.0912.702.68Shroud secondary
81.810.306.03Steam chamber
98.425.081.66Shell SH
Rings and BB
1236.04146.50.246Hot ring
1336.04146.50.246Cold ring
167.7317.90.027Cold ring distributors IB
1814.6713.10.07Cold ring distributors OB
177.9918.50.027Hot ring collectors IB
1915.914.20.07Hot ring collectors OB
1473.714.7-BB inboard
15121.61612.3-BB outboard
Table 12. FL number and characteristics.
Table 12. FL number and characteristics.
Nodalization
FL
From
CVH
To
CVH
From Height
(m)
To Height
(m)
Section
(m2)
Hydraulic Diameter
(m)
Length
(m)
Primary side
FL001001214.2516.52720.41170.72445
FL002001215.373215.37321.18181.22670
FL0030023114.763614.76361.18181.22670
FL03100 TO FL04600(31)(46)13.981113.94911.18180.01410.032
FL004004641.76361.76361.18181.22670
FL00500451.1541.1541.18181.22670
Secondary side
FL00600106990.1986390.08890
FL07100 TO FL08600(71)(86)13.949113.98110.894660.0126360.032
FL007006712.22.24.3303490.845790
FL0080086814.463614.46362.68460.02790
Rings and BB
FL01200513020.2460.5623
FL0130013181.7241.7240.620.310.6
FL014001815−7.746−7.7460.5640.184110.6
FL015001519−7.746−7.7460.5640.184110.6
FL0160019123.7843.7840.620.310.6
FL0170013161.7241.7240.1560.18417.5
FL018001614−4.446−4.4460.2550.14294
FL019001417−4.446−4.4460.2550.14294
FL0200017123.7843.7840.1560.18417.5
Accident
FL10101338412.712.70.00280.01410
FL1010239787.97.90.00280.01410
FL1010345722.22.20.00280.01410
Table 13. Decay heat densities (MW/m3) integrated into the WCLL full reactor [24].
Table 13. Decay heat densities (MW/m3) integrated into the WCLL full reactor [24].
Entire Reactor
Name of ZoneNuclear HeatingCooling Time
MW/m3MW/m3
First wall (FW) 0 s1 s1 h1 day1 week
W2.18E+014.95E−014.93E-014.63E-012.27E-015.66E-03
Eurofer7.14E+001.76E-011.76E-011.01E-011.34E-026.91E-03
Breeder module (BM)
BM caps and lateral walls1.81E+001.08E-021.08E-027.32E-038.82E-042.57E-04
BM material mixture1.17E+001.23E-029.35E-033.60E-031.07E-039.01E-04
BM backwall1.06E-011.50E-031.50E-031.02E-031.20E-044.24E-05
BM back support/manifold5.64E-024.31E-044.28E-042.90E-043.45E-051.32E-05
Sum (MW/m3)1.65E+032.21E+011.91E+019.68E+002.26E+001.25E+00
Table 14. Summary of simulation parameters.
Table 14. Summary of simulation parameters.
Rupture LocationNumber of PipesRupture Size (m2)FL Involved
Case 1Top OTSG90.0028FL10101
Case 2Middle OTSG90.0028FL10102
Case 3Bottom OTSG90.0028FL10103
Table 15. Event sequence.
Table 15. Event sequence.
Event Time (s)
Steady state 0–1000
Low-pressure signalCase 11016.5
Case 21008.5
Case 31007.5
FPSCase 11019.5
Case 21011.5
Case 31010.5
End of the simulation 150,000
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Galleni, F.; Moscardini, M.; Pucciarelli, A.; Porfiri, M.T.; Forgione, N. Numerical Simulation of an Out-Vessel Loss of Coolant from the Breeder Primary Loop Due to Large Rupture of Tubes in a Primary Heat Exchanger in the DEMO WCLL Concept. Energies 2021, 14, 6916. https://doi.org/10.3390/en14216916

AMA Style

Galleni F, Moscardini M, Pucciarelli A, Porfiri MT, Forgione N. Numerical Simulation of an Out-Vessel Loss of Coolant from the Breeder Primary Loop Due to Large Rupture of Tubes in a Primary Heat Exchanger in the DEMO WCLL Concept. Energies. 2021; 14(21):6916. https://doi.org/10.3390/en14216916

Chicago/Turabian Style

Galleni, Francesco, Marigrazia Moscardini, Andrea Pucciarelli, Maria Teresa Porfiri, and Nicola Forgione. 2021. "Numerical Simulation of an Out-Vessel Loss of Coolant from the Breeder Primary Loop Due to Large Rupture of Tubes in a Primary Heat Exchanger in the DEMO WCLL Concept" Energies 14, no. 21: 6916. https://doi.org/10.3390/en14216916

APA Style

Galleni, F., Moscardini, M., Pucciarelli, A., Porfiri, M. T., & Forgione, N. (2021). Numerical Simulation of an Out-Vessel Loss of Coolant from the Breeder Primary Loop Due to Large Rupture of Tubes in a Primary Heat Exchanger in the DEMO WCLL Concept. Energies, 14(21), 6916. https://doi.org/10.3390/en14216916

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop