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.

At least 217 records · Page 12

Fukushima Daiichi Unit 3 MELCOR Investigation

The Department of Energy sponsored a Fukushima Daiichi accident study as a collaboration between Sandia, Oak Ridge (ORNL), and Idaho National Laboratories. The purpose of the effort was to compile relevant data, reconstruct the accident progression using computer codes, assess the codes predictive capabilities, and to identify future data needs. The following summarizes MELCOR simulations performed at ORNL on Fukushima Daiichi Unit 3.

Fukushima Daiichi Unit 3↗

PSA 2025 Presentation: "Modeling and Sensitivity Analysis of a Generation IV Pebble Bed Reactor Using MELCOR 2.2"

Accompanying the advancement of reactor technologies is the need for computational modeling and simulation to predict their behavior under normal operating conditions and accident scenarios. New Generation IV reactor designs which employ non-conventional fuel have a particular need for modeling the behavior and release of radionuclides and other material from the fuel. In this work, MELCOR version 2.2, a system-level safety and accident scenario code developed by Sandia National Laboratories, was used to model a 200-MWth pebble bed modular reactor and calculate the inventories of circulating and deposited graphite, metal dust, and elemental components released from the fuel elements. A base case modeling the reactor under standard operating conditions was calculated using MELCOR and the inventories were extrapolated to 30 years of operation time using a logarithmic regression fit. A sensitivity analysis was also performed in which several key parameters for the base case model were modified to explore the effect of these changes on the inventories calculated by MELCOR. A set of transient scenario simulations for a depressurized loss of forced cooling (DLOFC) accident were also performed. The results of the sensitivity analysis and transient simulations are reported and discussed in relation to the modeling techniques used for this study.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Sodium Fire Collaborative Study Progress CNWG Fiscal Year 2021

This report discusses the progress on the collaboration between Sandia National Laboratories (Sandia) and Japan Atomic Energy Agency (JAEA) on the sodium fire research in fiscal year (FY) 2021 and is a continuation of the FY2020 progress report. In this report, we only report the changes made to the current sodium pool fire model in MELCOR. We modified and corrected many control functions to enhance the current model from the FY2020 report. This year’s enhancements relate to better agreement of the suspended aerosol measurement from JAEA’s F7 series tests. Staff from Sandia and JAEA conducted the validation studies of the sodium pool fire model in MELCOR. To validate this pool fire model with the latest enhancement, JAEA sodium pool fire experiments (F7-1 and F7-2) were used. The results of the calculation, including the code-to-code comparisons are discussed as well as suggestions for further model improvement. Finally, recommendations are made for new MELCOR simulations for FY2022.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Sodium Fire Collaborative Study Progress (CNWG Fiscal Year 2022)

This report discusses the progress on the collaboration between Sandia National Laboratories (Sandia) and Japan Atomic Energy Agency (JAEA) on the sodium fire research in fiscal year (FY) 2022 and is a continuation of the FY 2021 progress report. We only report the changes made to the current sodium pool fire model in MELCOR. We modified and corrected many control functions to enhance the fraction of oxygen consumed that reacts to form monoxide (FO2) parameter in the current model from the FY2021 report. This year's enhancements relate to better agreement of the suspended aerosol measurement from JAEA's F7 series tests. Staff from Sandia and JAEA conducted the validation studies of the sodium pool fire model in MELCOR. To validate this pool fire model with the latest enhancement, JAEA sodium pool fire experiments (F7-1 and F7-2) were used. The results of the calculation, including the code-to-code comparisons are discussed as well as suggestions for further model improvement. Finally, recommendations are made for new MELCOR simulations for FY 2023.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Heat up and potential failure of BWR upper internals during a severe accident

In boiling water reactors, the steam dome, steam separators, and dryers above the core are comprised of approximately 100 tons of stainless steel. During a severe accident in which the coolant boils away and exothermic oxidation of zirconium occurs, gases (steam and hydrogen) are superheated in the core region and pass through the upper internals. Historically, the upper internals have been modeled using severe accident codes with relatively simple approximations. The upper internals are typically modeled in MELCOR as two lumped volumes with simplified heat transfer characteristics, with no structural integrity considerations, and with limited ability to oxidize, melt, and relocate. The potential for and the subsequent impact of the upper internals to heat up, oxidize, fail, and relocate during a severe accident was investigated. A higher fidelity representation of the shroud dome, steam separators, and steam driers was developed in MELCOR v1.8.6 by extending the core region upwards. This modeling effort entailed adding 45 additional core cells and control volumes, 98 flow paths, and numerous control functions. The model accounts for the mechanical loading and structural integrity, oxidation, melting, flow area blockage, and relocation of the various components. The results indicate that the upper internals can reach high temperatures during a severe accident; they are predicted to reach a high enough temperature such that they lose their structural integrity and relocate. The additional 100 tons of stainless steel debris influences the subsequent in-vessel and ex-vessel accident progression.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Insights on Fukushima Damage Progression based on PCV Inspections and Implications for Decommissioning Data Collection and Code Model Refinement [Slides]

Outline points are: Review what is known from experiments and how codes are modeling phenomena; Materials interactions are very important and key interactions will be identified and discussed; Chronology of damage progression roughly follows in order of increasing melting/liquefaction temperatures; Examine a plausible sequence to explain robotic visual examinations; Highlight MELCOR modeling observations; Highlight potential decommissioning phase data collection needs; and, Knowledge advance is iterative process of reconciling observations with code predictions, improving code models, and comparing to emerging new observations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Application of SCALE to Molten Salt Fueled Reactor Physics in Support of Severe Accident Analyses

As part of a US Nuclear Regulatory Commission–sponsored project to assess the modeling and simulation capabilities for accident progression, source term, and consequence analysis for advanced reactor technologies with SCALE and MELCOR, SCALE was used for the modeling and simulation of a molten salt-fueled reactor (MSR). SCALE capabilities for the modeling of MSR physics were demonstrated based on the Molten Salt Reactor Experiment (MSRE). Of primary interest were the determination of the system’s nuclide inventory, as well as the inventories in the various regions of the loop, considering that the fuel is continuously pumped through the system. This report contains discussions on the following: 1. Determination of the system-average fuel salt inventory considering fission gas removal in the off-gas system and noble metal removal through plating out at the heat exchanger using recent enhancements in SCALE’s depletion sequence TRITON, 2. Assessment of the nuclide spatial distribution throughout the loop using SCALE’s depletion solver ORIGEN, 3. Calculation of the core’s power profile, flux profile, temperature reactivity coefficients, and xenon reactivity using full-core calculations with SCALE’s Monte Carlo code KENO-VI. The results obtained with SCALE were post-processed to provide the MELCOR team with the core inventory and decay heat of the system, as well as the inventory and decay heat of individual regions in the loop, a zone-wise power profile, temperature feedback coefficients, and the xenon worth.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Investigation on the Core Bypass Flow in a Very High Temperature Reactor

Uncertainties associated with the core bypass flow are some of the key issues that directly influence the coolant mass flow distribution and magnitude, and thus the operational core temperature profiles, in the very high-temperature reactor (VHTR). Designers will attempt to configure the core geometry so the core cooling flow rate magnitude and distribution conform to the design values. The objective of this project is to study the bypass flow both experimentally and computationally. Researchers will develop experimental data using state-of-the-art particle image velocimetry in a small test facility. The team will attempt to obtain full field temperature distribution using racks of thermocouples. The experimental data are intended to benchmark computational fluid dynamics (CFD) codes by providing detailed information. These experimental data are urgently needed for validation of the CFD codes. The following are the project tasks: (1) Construct a small-scale bench-top experiment to resemble the bypass flow between the graphite blocks, varying parameters to address their impact on bypass flow. Wall roughness of the graphite block walls, spacing between the blocks, and temperature of the blocks are some of the parameters to be tested. (2) Perform CFD to evaluate pre- and post-test calculations and turbulence models, including sensitivity studies to achieve high accuracy. (3) Develop the state-of-the art large eddy simulation (LES) using appropriate subgrid modeling. (4) Develop models to be used in systems thermal hydraulics codes to account and estimate the bypass flows. These computer programs include, among others, RELAP3D, MELCOR, GAMMA, and GAS-NET. Actual core bypass flow rate may vary considerably from the design value. Although the uncertainty of the bypass flow rate is not known, some sources have stated that the bypass flow rates in the Fort St. Vrain reactor were between 8 and 25 percent of the total reactor mass flow rate. If bypass flow rates are on the high side, the quantity of cooling flow through the core may be considerably less than the nominal design value, causing some regions of the core to operate at temperatures in excess of the design values. These effects are postulated to lead to localized hot regions in the core that must be considered when evaluating the VHTR operational and accident scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Pruning of Discrete Dynamic Event Trees Using Density Peaks and Dynamic Time Warping

Discrete Dynamic Event Trees (DDETs) are used in probabilistic risk assessment (PRA) to determine the risk associated with complex systems, such as nuclear power plants, where accident scenarios need to be realistically modeled. One of the most challenging aspects of the DDET method is that a large number of branches are produced for each initiating event. This can result in data that is difficult to organize and time consuming to analyze. As the number of uncertain parameters increases, so does the complexity of the event tree. An event tree often has nodes that do not contribute to the overall results and predictive capabilities; therefore, we suggest a pre-pruning algorithm that can work to remove these nodes, and consequently, reduce the complexity of the tree. Sandia National Laboratories (SNL) is currently working on a two year dynamic cyber-risk informed project to assess cyber failure in a nuclear power plant. This will be achieved by coupling SNL's dynamic scheduler known as Analysis of Dynamic Accident Progression Trees (ADAPT) with MELCOR, a reactor accident simulator, and a System for Circuit Evaluation and Prediction of Transient Radiation Effects (SCEPTRE), a cyber Emulytics code. This paper will discuss the preliminary results of a feasibility study for using density peaks (DP) clustering with dynamic time warping (DTW) to prune branches that are formed from large DDETs, like the ones that will be produced throughout SNL's project. Papers from the project that are submitted to this conference discuss the following: - The MELCOR model for a Pressurized Water Reactor that uses a digital instrumentation and control system. - Proposed generic cyber exploits for the above reactor. - Modifications to ADAPT to extend the framework for multiple simulators, incorporate importance measures into DDET analysis, as well as sort and reduce DDETs based on inputs to or outputs from the tree. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

ITER Port Interspace Pressure Calculations

The ITER ('the way' in Latin) project in Cadarache, France, will operate the largest tokamak ever built to demonstrate the physical and technical feasibility of nuclear fusion as an energy source. The ITER vacuum vessel is equipped with 50 access ports. Each port has an opening in the bio-shield that communicates with a dedicated port cell. During tokamak operation, the bio-shield opening must be closed with a concrete plug to shield the radiation coming from the plasma. This port plug separates the port cell into a port interspace (between vacuum vessel closure lid and port plug) on the inner side and a port cell on the outer side. This paper presents calculations of pressures and temperatures in the ITER port interspace after a double-ended guillotine break (DEGB) of a pipe in the tokamak cooling water system (TCWS). It is assumed that the DEGB occurs during the worst possible conditions-during water baking operation, with water at 523 K (250 deg. C) and 4.4 MPa. These conditions are more severe than during normal tokamak operation with the water at 398 K (125 deg. C) and 2 MPa. Two computer codes are employed in these calculations: RELAP5-3D Version 4.2.1 to calculate blowdown releases from the pipe break, and MELCOR, Version 1.8.6 to calculate pressures and temperatures in the port interspace. A sensitivity study has been performed to optimize the maximum peak pressure in the port interspace volume of ITER. The optimum case was for a labyrinth area (FL131) of 0.5 m{sup 2} and a relief panel (FL135) kept open after the first cycle. This case keeps the pressure in the port interspace under 150 kPa and the temperatures under 385 K. (authors)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modular Accident Analysis Program (MAAP) - MELCOR Crosswalk: Phase II Analyzing a Partially Recovered Accident Scenario

Following the conclusion of the first phase of the crosswalk analysis, one of the key unanswered questions was whether or not the deviations found would persist during a partially recovered accident scenario, similar to the one that occurred in TMI - 2. In particular this analysis aims to compare the impact of core degradation morphology on quenching models inherent within the two codes and the coolability of debris during partially recovered accidents. A primary motivation for this study is the development of insights into how uncertainties in core damage progression models impact the ability to assess the potential for recovery of a degraded core. These quench and core recovery models are of the most interest when there is a significant amount of core damage, but intact and degraded fuel still remain in the cor e region or the lower plenum. Accordingly this analysis presents a spectrum of partially recovered accident scenarios by varying both water injection timing and rate to highlight the impact of core degradation phenomena on recovered accident scenarios. This analysis uses the newly released MELCOR 2.2 rev. 966 5 and MAAP5, Version 5.04. These code versions, which incorporate a significant number of modifications that have been driven by analyses and forensic evidence obtained from the Fukushima - Daiichi reactor site.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

ITER Port Interspace Pressure Calculations

The ITER Vacuum Vessel (VV) is equipped with 54 access ports. Each of these ports has an opening in the bioshield that communicates with a dedicated port cell. During Tokamak operation, the bioshield opening must be closed with a concrete plug to shield the radiation coming from the plasma. This port plug separates the port cell into a Port Interspace (between VV closure lid and Port Plug) on the inner side and the Port Cell on the outer side. This paper presents calculations of pressures and temperatures in the ITER (Ref. 1) Port Interspace after a double-ended guillotine break (DEGB) of a pipe of the Tokamak Cooling Water System (TCWS) with high temperature water. It is assumed that this DEGB occurs during the worst possible conditions, which are during water baking operation, with water at a temperature of 523 K (250 C) and at a pressure of 4.4 MPa. These conditions are more severe than during normal Tokamak operation, with the water at 398 K (125 C) and 2 MPa. Two computer codes are employed in these calculations: RELAP5-3D Version 4.2.1 (Ref. 2) to calculate the blowdown releases from the pipe break, and MELCOR, Version 1.8.6 (Ref. 3) to calculate the pressures and temperatures in the Port Interspace. A sensitivity study has been performed to optimize some flow areas.

Interspace ITER Vacuum vessel↗

System code evaluation of near-term accident tolerant claddings during pressurized water reactor station blackout accidents

Following the Fukushima Daiichi nuclear accident in 2011, researches on Accident-tolerant fuels (ATFs) are currently of high interest in not only the nuclear industry but also governmental and international organizations. In this work, a quantitative evaluation of the performance of monolithic FeCrAl cladding and Cr-coated Zircaloy cladding has been performed for Pressurized Water Reactor (PWR) Station Blackout (SBO) accidents. A generic PWR model has been built in system thermal-hydraulics code TRACE based on the Surry Nuclear Power Station with counter-current natural circulation modelling capability for hotleg and steam generator U-tube components during the accidents. The base model results are then compared to MELCOR and RELAP simulations to verify the system component implementation in TRACE. Two PWR SBO scenarios were investigated, including: short-term SBO and long-term SBO with early reactor coolant pump (RCP) seal failure. These scenarios are defined to be very similar to the accidents studied in the State-of-the-Art Reactor Consequence Analysis (SOARCA) project. TRACE code is modified to reflect the oxidation kinetics of FeCrAl and Cr-coating. Larson-Miller creep rupture model is also implemented in TARCE using its built-in control systems to simulate the creep rupture of hotlegs. Additionally, the comparison between the TRACE models with and without the counter-current flow modeling resulted in significant difference when comparing ATF cladding to Zircaloy for short term SBO, while it marginal impacted the performance during long term SBO with RCP seal failure. For short term SBO, both ATF cladding underwent hot leg creep rupture ~20 min after Zircaloy cladding. While Zircaloy and Cr-coated cladding had generated significant amount of hydrogen gas (>10 kg) before the creep rupture event, FeCrAl cladding had only generated <0.5 kg of hydrogen gas. For long term SBO with RCP seal failure, significant hydrogen generation and fuel melting was predicted before hot leg creep rupture for the ATF cladding while providing only 10–20 min additional coping time compared to Zircaloy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

United States Department of Energy severe accident research following the Fukushima Daiichi accidents

The U.S. Department of Energy (DOE) has played a major role in the U.S. response to the events at Fukushima Daiichi. During the first several weeks following the accident, U.S. assistance efforts were guided by results from a significant and diverse set of analyses. In the months that followed, a coordinated analysis activity aimed at gaining a more thorough understanding of the accident sequence was completed using laboratory-developed, system-level best-estimate accident analysis codes, while a parallel analysis was conducted by U.S. industry. A comparison of predictions for Unit 1 from these two studies indicated significant differences between MAAP and MELCOR results for key plant parameters, such as in-core hydrogen production. On that basis, a crosswalk was completed to determine the key modeling variations that led to these differences. In parallel with these activities, it became clear that there was a need to perform a technology gap evaluation on accident-tolerant components and severe accident analysis methodologies with the goal of identifying any data and/or knowledge gaps that may exist given the current state of light water reactor (LWR) severe accident research and augmented by insights from Fukushima. In addition, there is growing international recognition that data from Fukushima could significantly reduce uncertainties related to severe accident progression, particularly for boiling water reactors. On these bases, a group of U. S. experts in LWR safety and plant operations was convened by the DOE Office of Nuclear Energy (DOE-NE) to complete technology gap analysis and Fukushima forensics data needs identification activities. The results from these activities were used as the basis for refining DOE-NE's severe accident research and development (R&D) plan. Finally, this paper provides a high-level review of DOE-sponsored R&D efforts in these areas, including planned activities on accident-tolerant components and accident analysis methods.

61 RADIATION PROTECTION AND DOSIMETRY↗

Response Technical Tools User Guide [Manual]

Developed by the U.S. Nuclear Regulatory Commission for assessment of reactor accidents, Response Technical Manual (RTM) is a paper-based report which contains simple methods for estimating possible accident scenarios and relevant consequences for different kinds of radiological events. Based on RTM, a software called Response Technical Tools (RTT) was developed by Sandia National Laboratories to convert the paper manual into an automated and easy-to-use code. In particular, the RTT focuses on the nuclear power plant severe accidents and is informed by state-of-the-art analyses and software programs, such as MELCOR and MAAP. RTT evaluations can be used to track and predict, at a very coarse level, the progression of a severe accident in nuclear power plant. The RTT allows a user to track the progress of an accident in a nuclear power plant from the point of initiation through the point of containment breach and release to the environment.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Integrating static PRA information with risk informed safety margin characterization (RISMC) simulation methods

The overall objective of the project was to develop a computationally feasible and user-friendly mechanized process to integrate traditional probabilistic risk assessment (PRA) and dynamic PRA (DPRA) results. Starting with the systematic identification of items in an existing PRA that need dynamic augmentation, the project used a generic 4-loop pressurized reactor (PWR) and 3-loop PWR as example plants. Station blackout (SBO) and large break loss of coolant accident SBLOCA) were selected as the example initiating events. Using the traditional event-tree (ET)/fault-tree (FT) methodology augmented by dynamic evet tree approach, the potential consequences of the initiating events were simulated with RELAP-3D and MELCOR/RASCAL codes to cover Level 1 through Level 3 of PRA. RAVEN and ADAPT software were used to generate Level 1 simulations with RELAP-3D and Level 2/3 simulations with MELCOR (Level 2)/RASCAL (Level 3), respectively. Example branching conditions (BCs) for SBO included AC power recovery time, valve repair failure time, reactor coolant pump leak time/break size and emergency power supply duration to a total of 9. Example BCs for LOCA included off-site power recovery time, diesel generator power recovery time, auxiliary feed water system operation time, safety relief valve failure to open upon demand, reactor coolant pump seal break time and size to a total of 21. Each RELAP-3D simulation (9,587 scenarios) was labelled OK or Core Damage based on the maximum allowed peak clad temperature (2,100oF). Each MELCOR simulation (4610 scenarios) was labeled as Bin over 10rem or Bin 0-10rem based on the dose at the site boundary. The scenarios were clustered based on the criteria above using the mean shift methodology. Classical PRA (CPRA) and DPRA results were compared to identify the ET sequences that need DPRA augmentation. Several approaches were proposed for the incorporation of these sequences into CPRA using clustering with the mean shift methodology, restructuring the CPRA ETs by adding new BCs/sequences, and using the concept of a limit surface. Procedures for decision making regarding the possible consequences of an initiating event (e.g. core damage or not, site evacuation or not) were developed using a convolutional neural network (CNN), a recurrent neural network (RNN) and a transformer neural network (TNN). The project has led to two PhD degrees, three archival journal papers and five refereed conference proceedings.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Preliminary Assessment of the Dose near the Shield Plug at Fukushima Daiichi Unit 2 (1F2)

Since the accident at Fukushima Daiichi on March 11, 2011, caused by the Great East Japan earthquake and subsequent tsunami, there have been significant coordinated international efforts to better understand the accident progression, phenomenology and consequences. Of particular note is the Organisation for Economic Cooperation and Development Nuclear Energy Agency's (OECD/NEA) Benchmark Study of the Accident at the Fukushima Daiichi Nuclear Power Plant (BSAF Project). Phase 1 of this ongoing effort focused on modeling the modeling of the three beyond design basis accidents at the site for the first six days of the accident. Phase 2 of the effort is a forensic effort and focuses on extending the accident simulations out to 21 days and assessing the consequences of the accident. This effort aims to inform the government of Japan and The Tokyo Electric Power Company (TEPCO)'s ongoing effort to decommission and decontaminate the site. One of the key questions posed by TEPCO in its efforts to decommission the site is the location of released radionuclides (RNs) that result in anomalously high dose rate measurements on and near the shield plugs of all three units. In this summary, we provide insight into the possible location of the source resulting in this high dose. Previous MELCOR accident progression analyses have predicted leakage of the containment through the drywell head, directly below the shield plug. Using this plausible accident scenario, MELCOR is able to predict the radionuclide source to the region below the shield plug and possible deposition amounts within the gaps of the shield plug. Using this predicted source and the transport and shielding code MCNP6 the dose in and around the shield plug of 1F2 was assessed. Unit 2 was chosen since the refueling bay is still intact and therefore a simpler initial problem compared to Units 1 and 3. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Assessment of Microreactor Safety Analysis Challenges and Recommendations for Utilization of the Comprehensive Reactor Analysis Bundle

To enable the broad deployment of microreactors in fundamentally new application regimes (i.e., mobile and autonomous operations), their safety must be indisputable in terms of possessing inherent resistance to severe offsite dose consequences. Therefore, mechanistic beyond-design-basis event source term calculations that demonstrate a sufficiently large margin of safety will be required to accommodate these new application regimes, which have no history of commercial regulation. Even for traditional reactor operation configurations, safety analysis expertise and familiarity with accident phenomena and conditions in microreactors—specifically those with heat pipe primary cooling arrangements—are lacking compared with other advanced reactor concepts and small modular reactors. Recently, modeling and simulation tools to account for unique heat pipe design aspects have been developed by Sandia National Laboratories with MELCOR and by the US Department of Energy’s (DOE’s) Office of Nuclear Energy Advanced Modeling and Simulation Program with BlueCRAB. However, further demonstration and assessment of potential knowledge gaps are needed to support these codes’ broad usage by the microreactor community. Through the DOE Microreactor Program, an initial assessment of these two tools and guidance on how an evaluation model could be constructed was performed and is reported herein. Moreover, a proposed approach for demonstrating an evaluation model using these two tools is outlined.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗