Search NASA⌕ Search

SEARCH · Search NASA

Results for “MELCOR Code”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

223 records · Page 13

Terry Turbopump Expanded Operating Band Modeling and Simulation Efforts in Fiscal Year 2020 - Progress Report

The Terry Turbine Expanded Operating Band Project is currently conducting testing at Texas A&M University as part of a revised experimental program meant to supplant previous full-scale testing plans under the headings of Milestone 5 and Milestone 6. In consultation with Sandia National Laboratories technical staff and with modeling and simulation support from the same, the hybrid Milestone 5&6 plan is moving forward with experiments aimed at addressing knowledge gaps regarding scale, working fluid, and turbopump self-regulation. Modeling and simulation efforts at Sandia National Laboratories in FY20 fell under the broad umbrella of Milestone 7 and consisted exclusively of MELCOR-related tasks aimed at: 1) Constructing/improving input models of Texas A&M University experiments, 2) Constructing a generic boiling water reactor input model according to best practices with systems-level Teny turbine capabilities, and 3) Adding code capability in order to leverage experimental data/findings, address bugs, and improve general code robustness Project impacts of the Covid-19 pandemic have fortunately been minimal thus far but are mentioned as necessary when discussing the hybrid Milestone 5&6 progress as well as the corresponding Milestone 7 modeling and simulation progress.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

MELCOR Model of the Spent Fuel Pool of Fukushima Dai-ichi Unit 4

Unit 4 of the Fukushima Dai-ichi Nuclear Power Plant suffered a hydrogen explosion at 6:00 am on March 15, 2011, exactly 3.64 days after the earthquake hit the plant and the off-site power was lost. The earthquake occurred on March 11 at 2:47 pm. Since the reactor of this Unit 4 was defueled on November 29, 2010, and all its fuel was stored in the spent fuel pool (SFP4), it was first believed that the explosion was caused by hydrogen generated by the spent fuel, in particular, by the recently discharged core. The hypothetical scenario was: power was lost, cooling to the SFP4 water was lost, pool water heated/boiled, water level decreased, fuel was uncovered, hot Zircaloy reacted with steam, hydrogen was generated and accumulated above the pool, and the explosion occurred. Recent analyses of the radioisotopes present in the water of the SFP4 and underwater video indicated that this scenario did not occur - the fuel in this pool was not damaged and was never uncovered the hydrogen of the explosion was apparently generated in Unit 3 and transported through exhaust ducts that shared the same chimney with Unit 4. This paper will try to answer the following questions: Could that hypothetical scenario in the SFP4 had occurred? Could the spent fuel in the SPF4 generate enough hydrogen to produce the explosion that occurred 3.64 days after the earthquake? Given the magnitude of the explosion, it was estimated that at least 150 kg of hydrogen had to be generated. As part of the investigations of this accident, MELCOR models of the SFP4 were prepared and a series of calculations were completed. The latest version of MELCOR, version 2.1 (Ref. 1), was employed in these calculations. The spent fuel pool option for BWR fuel was selected in MELCOR. The MELCOR model of the SFP4 consists of a total of 1535 fuel assemblies out of which 548 assemblies are from the core defueled on Nov. 29, 2010, 783 assemblies are older assemblies, and 204 are new/fresh assemblies. The total decay heat of the fuel in the pool was, at the time of the accident, 2.284 MWt, of which 1.872 MWt were from the 548 assemblies of the last core discharged and 0.412 MWt were from the older 783 assemblies. These decay heat values were calculated at Oak Ridge National Laboratory using the ORIGEN2.2 code (Ref. 2) - they agree with values reported elsewhere (Ref. 3). The pool dimensions are 9.9 m x 12.2 m x 11.8 m (height), and with the water level at 11.5 m, the pool volume is 1389 m3, of which only 1240 m3 is water, as some volume is taken by the fuel and by the fuel racks. The initial water temperature of the SFP4 was assumed to be 301 K. The fuel racks are made of an aluminum alloy but are modeled in MELCOR with stainless steel and B4C. MELCOR calculations were completed for different initial water levels: 11.5 m (pool almost full, water is only 0.3 m below the top rim), 4.4577 m (top of the racks), 4.2 m, and 4.026 m (top of the active fuel). A calculation was also completed for a rapid loss of water due to a leak at the bottom of the pool, with the fuel rapidly uncovered and oxidized in air. Results of these calculations are shown in the enclosed Table I. The calculation with the initial water level at 11.5 m (full pool) takes 11 days for the water to boil down to the top of the fuel racks, 11.5 days for the fuel to be uncovered, 14.65 days to generate 150 kg of hydrogen and 19 days for the pool to be completely dry. The calculation with the initial water level at 4.4577 m, takes 1.1 days to uncover the fuel and 4.17 days to generate 150 kg of hydrogen. The calculation with the initial water level at 4.02 m takes 3.63 days to generate 150 kg of hydrogen this is exactly the time when the actual explosion occurred in Unit 4. Finally, fuel oxidation in air after the pool drained the water in 20 minutes, generates only 10 kg of hydrogen this is because very little steam is available and Zircaloy (Zr) oxidation with the oxygen of the air does not generate hydrogen. MELCOR calculated water levels and hydrogen generated in the SFP4 as a function of time for initial water levels of 4.457 m, 4.2 m and 4.02 m are shown in Figs. 1 and 2. Water levels increase at the beginning due to the expansion of the water during the heat-up from 301 K to 373 K. Boiling occurs after the water temperature reaches 373 K. The total amount of hydrogen generated is ~2000 kg, this amount includes hydrogen generated from Zr, which is the largest amount (~1580 kg), from stainless steel (~360 kg), and from B4C (~60 kg). In theory, it is possible to generate up to 3.4 kg of hydrogen per assembly (from oxidation of Zr in the fuel cladding and box), or a total of 4,525 kg from the hot 1331 assemblies stored in the SFP4. The hydrogen generated from oxidation of steel and B4C will be additional. So the answers to the questions are YES according to these MELCOR calculations, enough hydrogen (150 kg) could be generated in the SFP4 3.64 days after the earthquake to produce ...

08 HYDROGEN↗

FY24: Progress Report on the MELCOR Modeling of the Liquid Salt Test Loop

This report outlines the activities conducted in FY24 focused on updating the liquid salt test loop (LSTL) model through the integration of a new test section featuring 16 radial tubes coupled to a filter section. Some comparative analyses of the updated model with existing experimental data for the LSTL was made and benchmarked against other computational tools, such as the ORNL code, SAM. These actions are part of a comprehensive validation effort and promote collaboration among laboratories participating in the Molten Salt Reactor (MSR) campaign.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Integral Reactor Containment Condensation Model and Experimental Validation

This NEUP funded project, NEUP 12-3630, is for experimental, numerical and analytical studies on high-pressure steam condensation phenomena in a steel containment vessel connected to a water cooling tank, carried out at Oregon State University (OrSU) and the University of Wisconsin at Madison (UW-Madison). In the three years of investigation duration, following the original proposal, the planned tasks have been completed: (1) Performed a scaling study for the full pressure test facility applicable to the reference design for the condensation heat transfer process during design basis accidents (DBAs), modified the existing test facility to route the steady-state secondary steam flow into the high pressure containment for controllable condensation tests, and extended the operations at negative gage pressure conditions (OrSU). (2) Conducted a series of DBA and quasi-steady experiments using the full pressure test facility to provide a reliable high pressure condensation database (OrSU). (3) Analyzed experimental data and evaluated condensation model for the experimental conditions, and predicted the prototypic containment performance under accidental conditions (UW-Madison). A film flow model was developed for the scaling analysis, and the results suggest that the 1/3 scaled test facility covers large portion of laminar film flow, leading to a lower average heat transfer coefficient comparing to the prototypic value. Although it is conservative in reactor safety analysis, the significant reduction of heat transfer coefficient (50%) could under estimate the prototypic condensation heat transfer rate, resulting in inaccurate prediction of the decay heat removal capability. Further investigation is thus needed to quantify the scaling distortion for safety analysis code validation. Experimental investigations were performed in the existing MASLWR test facility at OrST with minor modifications. A total of 13 containment condensation tests were conducted for pressure ranging from 4 to 21 bar with three different static inventories of non-condensable gas. Condensation and heat transfer rates were evaluated employing several methods, notably from measured temperature gradients in the HTP as well as measured condensate formation rates. A detailed mass and energy accounting was used to assess the various measurement methods and to support simplifying assumptions required for the analysis. Condensation heat fluxes and heat transfer coefficients are calculated and presented as a function of pressure to satisfy the objectives of this investigation. The major conclusions for those tests are summarized below: (1) In the steam blow-down tests, the initial condensation heat transfer process involves the heating-up of the containment heat transfer plate. An inverse heat conduction model was developed to capture the rapid transient transfer characteristics, and the analysis method is applicable to SMR safety analysis. (2) The average condensation heat transfer coefficients for different pressure conditions and non-condensable gas mass fractions were obtained from the integral test facility, through the measurements of the heat conduction rate across the containment heat transfer plate, and from the water condensation rates measurement based on the total energy balance equation. 15 (3) The test results using the measured HTP wall temperatures are considerably lower than popular condensation models would predict mainly due to the side wall conduction effects in the existing MASLWR integral test facility. The data revealed the detailed heat transfer characteristics of the model containment, important to the SMR safety analysis and the validation of associated evaluation model. However this approach, unlike separate effect tests, cannot isolate the condensation heat transfer coefficient over the containment wall, and therefore is not suitable for the assessment of the condensation heat transfer coefficient against system pressure and noncondensable gas mass fraction. (4) The average condensation heat transfer coefficients measured from the water condensation rates through energy balance analysis are appropriate, however, with considerable uncertainties due to the heat loss and temperature distribution on the containment wall. With the consideration of the side wall conduction effects, the results indicate that the measured heat transfer coefficients in the tests is about 20% lower than the prediction of Dehbi’s correlation, mainly due to the side wall conduction effects. The investigation also indicates an increase in the condensation heat transfer coefficient at high containment pressure conditions, but the uncertainties invoked with this method appear to be substantial. (5) Non-condensable gas in the tests has little effects on the condensation heat transfer at high elevation measurement ports. It does affect the bottom measurements near the water level position. The results suggest that the heavier non-condensable gas is accumulated in the lower portion of the containment due to stratification in the narrow containment space. The overall effects of the non-condensable gas on the heat transfer process should thus be negligible for tall containments of narrow condensation spaces in most SMR designs. Therefore, the previous correlations with noncondensable gas effects are not appropriate to those small SMR containments due to the very poor mixing of steam and non-condensable gas. The MELCOR simulation results agree with the experimental data reasonably well. However, it is observed that the MELCOR overpredicts the heat flux for all analyzed tests. The MELCOR predicts that the heat fluxes for CCT’s approximately range from 30 to 45 kW/m2 whereas the experimental data (averaged) ranges from about 25 to 40 kW/m2. This may be due to the limited availability of liquid film models included in MELCOR. Also, it is believed that due to complex test geometry, measured temperature gradients across the heat transfer plate may have been underestimated and thus the heat flux had been underestimated. The MELCOR model predicts a film thickness on the order of 100 microns, which agrees very well with film flow model developed in this study for scaling analysis. However, the expected differences in film thicknesses for near vacuum and near atmospheric test conditions are not significant. Further study on the behavior of condensate film is expected to refine the simulation results. Possible refinements include but are not limited to, the followings: CFD simulation focusing on the liquid film behavior and benchmarking with experimental analyses for simpler geometries. 16 1 INTRODUCTION This NEUP funded project, NEUP 12-3630, is for experimental, numerical and analytical studies on high-pressure steam condensation phenomena in a steel containment vessel connected to a water cooling tank, carried out at Oregon State University (OrSU) and the University of Wisconsin at Madison (UW-Madison). The experimental results are employed to validate the containment condensation model in reactor containment system safety analysis code for integral SMRs. Such a containment condensation model is important to demonstrate the adequate cooling. In the three years of investigation, following the original proposal, the following planned tasks have been completed: (1) Performed a scaling study for the full pressure test facility applicable to the reference design for the condensation heat transfer process during design basis accidents (DBAs), modified the existing test facility to route the steady-state secondary steam flow into the high pressure containment for controllable condensation tests, and extended the operations at negative gage pressure conditions (OrSU). (2) Conducted a series of DBA and quasi-steady experiments using the full pressure test facility to provide a reliable high pressure condensation database (OrSU). (3) Analyzed experimental data and evaluated condensation model for the experimental conditions, and predicted the prototypic containment performance under accidental conditions (UW-Madison). The results are applicable to integral Small Modular Reactor (SMR) designs, including NuScale, mPower, Westinghouse SMR, Holtec-160 and other integral reactors with small containments of relatively high pressures under accidental conditions. Testing has been conducted at the OrSU laboratory in the existing MASLWR (Multi-Application Small Light Water Reactor) integral test facility sponsored by the US Department of Energy. Its highpressure stainless steel containment model (~2 MPa) is scaled to the NuScale SMR currently under development at NuScale Power, Inc.. Minor modifications to the model containment have been made to control the non-condensable gas fraction and to utilize the secondary loop stable steam flow for condensation testing. UW-Madison has developed a containment condensation model, which leveraged previous validated containment heat transfer work carried out at UW-Madison, and extended the range of applicability of the model to integral SMR designs that utilize containment vessels of high heat transfer efficiencies. In this final report, the research background and literature survey are presented in Chapter 2 and 3, respectively. The test facility description and modifications are summarized in Chapter 4, and the scaling analysis is introduced in Chapter 5. The tests description, procedures, and data analysis are presented in Chapter 6, while the numerical modeling is presented in Chapter 7, followed by a conclusion section in Chapter 8.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Interim MELCOR Simulation of the Fukushima Daiichi Unit 2 Accident Reactor Core Isolation Cooling Operation

Data, a brief description of key boundary conditions, and results of Sandia National Laboratories’ ongoing MELCOR analysis of the Fukushima Unit 2 accident are given for the reactor core isolation cooling (RCIC) system. Important assumptions and related boundary conditions in the current analysis additional to or different than what was assumed/imposed in the work of SAND2012-6173 are identified. This work is for the U.S. Department of Energy’s Nuclear Energy University Programs fiscal year 2014 Reactor Safety Technologies Research and Development Program RC-7: RCIC Performance under Severe Accident Conditions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Overview of recent SCALE activities for Non-LWR inventory and decay heat analysis

In 2019, the US Nuclear Regulatory Commission initiated a project for the development and assessments of non-light-water reactor (non-LWR) accident progression using the SCALE and MELCOR simulation tools. SCALE simulations are used to generate nuclide inventories, full-core power distributions, decay heat, and kinetics parameters to initialize MELCOR simulations of severe accident scenarios. Five non-LWR concepts were studied: high-temperature gas-cooled reactor (HTGR), heat pipe reactor (HPR), high-temperature fluoride salt-cooled reactor (FHR), molten salt-fueled reactor (MSR), and sodium-cooled fast reactor (SFR). This paper summarizes the SCALE results obtained in 2021 for the first three non-LWR concepts, compares characteristics and results to common LWRs, and provides the strategy for the analysis of the remaining two non-LWRs. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗