Boiling water reactors use the Pressure Suppression Pool (PSP) to relieve the containment pressure in case of an accident. During the event of a Loss of Coolant Accident (LOCA), drywell air and steam are injected into the PSP through blowdown pipes. This may lead to thermal stratification, which is a relevant safety issue as it leads to higher water surface temperatures than in mixed conditions and thus, to higher containment pressures. The Effective Heat (EHS) and Momentum (EMS) Source models were previously introduced to predict the effect of small-scale direct contact condensation phenomena on the large-scale pool water circulation. In this paper, the EHS/EMS models are extended by adding the effect of non-condensable gases on the chugging regime. The EHS/EMS models are implemented in the GOTHIC code to model a full-scale Nordic BWR containment under different LOCA scenarios. The results show that thermal stratification can be developed in the PSP.

Nordic boiling water reactors (BWR) employ ex-vessel debris coolability as a severe accident management (SAM) strategy. Melt release mode into a deep water pool located in the lower drywell is a major source of uncertainty for success of such strategy. The main goal of this work is to identify the major contributors to the uncertainty in prediction of vessel failure mode, timing and melt release conditions in Nordic BWR using MELCOR severe accident analysis code. There is a forest of control rod guide tubes (CRGTs) and instrumentation guide tubes (IGTs) in the lower head of a BWR. Failure of the penetrations is considered in this work along with the gross creep rupture of the vessel lower head. Modelling of penetrations failure in MELCOR code is based on parametric models, allowing a user to select different assumptions, such as temperature threshold for penetrations failure, modelling of penetrations assembly and respective failure modes, mode of solid and liquid debris ejection from the vessel. In this work we perform sensitivity analysis to the MELCOR modelling options and sensitivity parameters in unmitigated station blackout scenario (SBO). Results of analysis suggest that (i) vessel breach due to penetration failure is predicted to occur before creep-rupture failure of the vessel lower head, (ii) the mode of debris ejection from the vessel has the dominant effect on the likelihood of creep-rupture failure of the vessel lower head and debris ejection rate from the vessel.

Nordic boiling water reactors (BWR) employ ex-vessel debris coolability as a severe accident management (SAM) strategy. Melt release mode into a deep water pool located in the lower drywell is a major source of uncertainty for success of such strategy. The main goal of this work is to identify the major contributors to the uncertainty in prediction of vessel failure mode, timing and melt release conditions in Nordic BWR using MELCOR severe accident analysis code. There is a forest of control rod guide tubes (CRGTs) and instrumentation guide tubes (IGTs) in the lower head of a BWR. Failure of the penetrations is considered in this work along with the gross creep rupture of the vessel lower head. Modelling of penetrations failure in MELCOR code is based on parametric models, allowing a user to select different assumptions, such as temperature threshold for penetrations failure, modelling of penetrations assembly and respective failure modes, mode of solid and liquid debris ejection from the vessel. In this work we perform sensitivity analysis to the MELCOR modelling options and sensitivity parameters in unmitigated station blackout scenario (SBO). Results of analysis suggest that (i) vessel breach due to penetration failure is predicted to occur before creep-rupture failure of the vessel lower head, (ii) the mode of debris ejection from the vessel has the dominant effect on the likelihood of creep-rupture failure of the vessel lower head and debris ejection rate from the vessel.

Nordic Boiling Water Reactors (BWRs) rely on the flooding of the lower drywell (LDW) as a severe accident management (SAM) strategy. The termination of a SA is achieved by fragmenting and quenching of the melt released from the vessel. Success of SAM strategy depends on melt release and water pool conditions. The characteristics of the melt release are the major source of uncertainty in quantification of the risk of SAM failure. Vessel failure and melt release modes are subject to aleatory and epistemic uncertainties at the in-vessel accident progression stage. In this work we focus on predicting the properties of debris relocated to the lower plenum using MELCOR code. We address the effect of epistemic uncertainty in modeling parameters and models in the MELCOR code in different severe accident scenarios on main characteristics of the in-vessel accident progression in Nordic BWRs. Sensitivity analysis is performed to rank the importance of MELCOR modelling parameters and the effect of different MELCOR models is addressed by using different versions of the code. The results provide valuable insights regarding the effect of MELCOR models, modelling parameters and sensitivity coefficients on code predictions.

Nordic Boiling Water Reactors (BWRs) employ ex-vessel debris coolability as a severe accident management strategy (SAM). Core melt is released into a deep pool of water where formation of noncoolable debris bed and ex-vessel steam explosion can pose credible threats to containment integrity. Success of the strategy depends on the scenario of melt release from the vessel that determines the melt-coolant interaction phenomena. The melt release conditions are determined by the in-vessel phase of severe accident progression. Specifically, properties of debris relocated into the lower plenum have influence on the vessel failure and melt release mode. In this work we use MELCOR code for prediction of the relocated debris. Over the years, many code modifications have been made to improve prediction of severe accident progression in light-water reactors. The main objective of this work is to evaluate the effect of models and best practices in different versions of MELCOR code on the in-vessel phase of different accident progression scenarios in Nordic BWR. The results of the analysis show that the MELCOR code versions 1.86 and 2.1 generate qualitatively similar results. Significant discrepancy in the timing of the core support failure and relocated debris mass in the MELCOR 2.2 compared to the MELCOR 1.86 and 2.1 has been found for a domain of scenarios with delayed time of depressurization. The discrepancies in the results can be explained by the changes in the modeling of degradation of the core components and changes in the Lipinski dryout model in MELCOR 2.2.

This work is motivated by the problem of restraining temperature escalation inside a porous heat-releasing media submerged in a pool of liquid coolant. When coolant temperature reaches saturation, boiling begins in the bulk of the porous bed, with void generation rate determined by the heating power. Amount of void determines hydrostatic pressure difference that drives natural circulation of two-phase flow through the porous material. At a certain critical value of the heat release rate, the driving head cannot overcome drag of the two phase porous media flow, which results in complete evaporation of coolant in some zone. Temperature of material in the dry zone increases significantly due to deterioration of heat exchange with single phase vapor flow in comparison with boiling flow heat transfer. The paper considers the problem of determining the critical conditions for onset of dryout in a heat-releasing porous bed of an arbitrary shape. The well-known one-dimensional problem for a flat top-flooded bed is revisited, and the functional form of the dryout boundary (expressed as the dryout heat flux, DHF) is derived using non-dimensional parameters. Asymptotic behavior of the solution is analyzed, and, by the method of asymptotic interpolation, a surrogate model is proposed consisting of three single-argument, non-dimensional functions. It is shown that such a model provides acceptable accuracy even in the cases where complete similarity of solutions is not achieved. The results obtained provide important insights into the physics of the problem, reduce the number of free parameters, and enable fast evaluation of dryout conditions without the need of numerical solution of algebraic equations involved in the exact formulation. The ultimate goal of the surrogate model development, i.e. its application to multidimensional configurations, is discussed.

The current paper describes the Computational Fluid-Dynamics (CFD) modelling of Heavy Liquid Metal (HLM) flows in a pool configuration and in particular how this is approached within the Horizon 2020 SESAME project. SESAME's work package structure, based on a systematic approach of redundancy and diversification, is explained along with its motivation. The main achievements obtained and the main lessons learned during the project are illustrated. The paper focuses on the strong coupling between the experimental activities and CFD simulations performed within the SESAME project. Two different HLM fluids are investigated: pure lead and Lead-Bismuth Eutectic. The objective is to make CFD a valid instrument used during the design of safe and innovative Gen-IV nuclear plants. Some effort has also been devoted to Proper Orthogonal Decomposition with Galerkin projection modelling (POD-Galerkin), a reduced order model suited for Uncertainty Quantification that operates by post-processing CFD results. Assessment of Uncertainty highly improves the reliability of CFD simulations. Dedicated experimental campaigns on heavily instrumented facilities have been conceived with the specific objective to build a series of datasets suited for the calibration and validation of the CFD modelling. In pool configuration, the attention is focused on the balance between conductive and convective heat transfer phenomena, on transient test-cases representative of incidental scenarios and on the possible occurrence of solidification phenomena. Four test sections have been selected to generate the datasets: (i) the CIRCE facility from ENEA, (ii) the TALL-3D pool test section from KTH, (iii) the TALL-3D Solidification Test Section (STS) from KTH and (iv) the SESAME Stand facility from CVR. While CIRCE and TALL-3D were existing facilities, the STS and SESAME Stand facility have been conceived, built and operated within the project, heavily relying on the use of CFD support. Care has been taken to ensure that almost all tasks were performed by at least two partners. Specific examples are given on how this strategy has allowed to uncover flaws and overcome pitfalls. Furthermore, an overview of the performed work and the achieved results is presented, as well as remaining or new uncovered issues. Finally, the paper is concluded with a description of one of the main goals of the SESAME project: the construction of the Gen-IV ALFRED CFD model and an investigation of its general circulation.

Spargers are multi-hole injection pipes used in Boiling Water Reactors (BWR) and Generation III/III+ Pressurized Water Reactors (PWR) to condense steam in large water pools. During the steam injection, high pool surface temperatures induced by thermal stratification can lead to higher containment pressures compared with completely mixed pool conditions, the former posing a threat for plant safety. The Effective Heat Source (EHS) and Effective Momentum Source (EMS) models were previously developed and validated for the modelling of a steam injection through blowdown pipes. The goal of this paper is to extend the EHS/EMS model capabilities towards steam injection through multi-hole spargers. The models are implemented in ANSYS Fluent 17.0 Computational Fluid Dynamics (CFD) code and calibrated against the spargers experiments performed in the PPOOLEX and PANDA facilities, analysed by the authors in Gallego-Marcos et al. (2018b). CFD modelling guidelines are established for the adequate simulation of the pool behaviour. A new correlation is proposed to model the turbulent production and dissipation caused by buoyancy. Sensitivity studies addressing the effect of different assumptions on the effective momentum magnitude, profile, angle and turbulence are presented. Calibration of the effective momentum showed an inverse proportionality to the sub-cooling. Differences between the effective momentum calibrated for PPOOLEX and PANDA are discussed. Analysis of the calculated flow above the cold stratified layer showed that the erosion of the layer is induced by the action of turbulence rather than mean shear flow.