The U.S. Nuclear Regulatory Commission initiated the state-of-the-art reactor consequence analyses (SOARCA) project to develop realistic estimates of the offsite radiological health consequences for potential severe reactor accidents. The SOARCA analysis of an ice condenser containment plant was performed because its relatively low design pressure and its reliance on igniters make it potentially susceptible to early containment failure from hydrogen combustion during a severe accident. The focus was on station blackout accident scenarios where all alternating current power is lost. Accident progression calculations used the MELCOR computer code and offsite consequence analyses were performed with MACCS. The analysis included more than 500 MELCOR and MACCS simulations to account for uncertainty in important accident progression and offsite consequence input parameters. Consequences from severe nuclear power plant accidents modeled in SOARCA are smaller than previously calculated. The delayed releases calculated provide more time for emergency response actions. The results show that early containment failure is very unlikely, even without successful use of igniters. The modeled behavior of safety valves is very important to this conclusion, but there is sparse data and a lack of established expert consensus on the failure rates under severe accident conditions. Even for scenarios resulting in early containment failure, the calculated individual latent fatal cancer risks are very small. Early and latent-cancer fatality risks are one focus of this paper. Regression results showing the most influential parameters are also discussed.
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.
Since the accidents at Fukushima-Daiichi, Sandia National Laboratories has been modeling these accident scenarios using the severe accident analysis code, MELCOR. MELCOR is a widely used computer code developed at Sandia National Laboratories since ~1982 for the U.S. Nuclear Regulatory Commission. Insights from the modeling of these accidents is being used to better inform future code development and potentially improved accident management. To date, our necessity to better capture in-vessel thermal-hydraulic and ex-vessel melt coolability and concrete interactions has led to the implementation of new models. The most recent analyses, presented in this paper, have been in support of the of the Organization for Economic Cooperation and Development Nuclear Energy Agency’s (OECD/NEA) Benchmark Study of the Accident at the Fukushima Daiichi Nuclear Power Station (BSAF) Project. The goal of this project is to accurately capture the source term from all three releases and then model the atmospheric dispersion. In order to do this, a forensic approach is being used in which available plant data and release timings is being used to inform the modeled MELCOR accident scenario. For example, containment failures, core slumping events and lower head failure timings are all enforced parameters in these analyses. This approach is fundamentally different from a blind code assessment analysis often used in standard problem exercises. The timings of these events are informed by representative spikes or decreases in plant data. The combination of improvements to the MELCOR source code resulting from analysis previous accident analysis and this forensic approach has allowed Sandia to generate representative and plausible source terms for all three accidents at Fukushima Daiichi out to three weeks after the accident to capture both early and late releases. In particular, using the source terms developed by MELCOR, the MACCS software code, which models atmospheric dispersion and deposition, we are able to reasonably capture the deposition of radionuclides to the northwest of the reactor site.
In this study, risk-significant pressurized-water reactor severe accident sequences are examined using MELCOR 1.8.5 to explore the range of fission product releases to the reactor containment building. Advances in the understanding of fission product release and transport behavior and severe accident progression are used to render best estimate analyses of selected accident sequences. Particular emphasis is placed on estimating the effects of high fuel burnup in contrast with low burnup on fission product releases to the containment. Supporting this emphasis, recent data available on fission product release from high-burnup (HBU) fuel from the French VERCOR project are used in this study. The results of these analyses are treated as samples from a population of accident sequences in order to employ approximate order statistics characterization of the results. These trends and tendencies are then compared to the NUREG-1465 alternative source term prescription used today for regulatory applications. In general, greater differences are observed between the state-of-the-art calculations for either HBU or low-burnup (LBU) fuel and the NUREG-1465 containment release fractions than exist between HBU and LBU release fractions. Current analyses suggest that retention of fission products within the vessel and the reactor coolant system (RCS) are greater than contemplated in the NUREG-1465 prescription, and that, overall, release fractions to the containment are therefore lower across the board in the present analyses than suggested in NUREG-1465. The decreased volatility of Cs 2 MoO 4 compared to CsI or CsOH increases the predicted RCS retention of cesium, and as a result, cesium and iodine do not follow identical behaviors with respect to distribution among vessel, RCS, and containment. With respect to the regulatory alternative source term, greater differences are observed between the NUREG-1465 prescription and both HBU and LBU predictions than exist between HBU and LBU analyses. Additionally, current analyses suggest that the NUREG-1465 release fractions are conservative by about a factor of 2 in terms of release fractions and that release durations for in-vessel and late in-vessel release periods are in fact longer than the NUREG-1465 durations. It is currently planned that a subsequent report will further characterize these results using more refined statistical methods, permitting a more precise reformulation of the NUREG-1465 alternative source term for both LBU and HBU fuels, with the most important finding being that the NUREG-1465 formula appears to embody significant conservatism compared to current best-estimate analyses. ACKNOWLEDGEMENTS This work was supported by the United States Nuclear Regulatory Commission, Office of Nuclear Regulatory Research. The authors would like to thank Dr. Ian Gauld and Dr. Germina Ilas, of Oak Ridge National Laboratory, for their contributions to this work. In addition to development of core fission product inventory and decay heat information for use in MELCOR models, their insights related to fuel management practices and resulting effects on spatial distribution of fission products in the core was instrumental in completion of our work.
Sandia National Laboratories (SNL) has conducted an uncertainty analysis (UA) on the Fukushima Daiichi unit (1F1) accident progression with the MELCOR code. The model used was developed for a previous accident reconstruction investigation jointly sponsored by the US Department of Energy (DOE) and Nuclear Regulatory Commission (NRC). That study focused on reconstructing the accident progressions, as postulated by the limited plant data. This work was focused evaluation of uncertainty in core damage progression behavior and its effect on key figures-of-merit (e.g., hydrogen production, reactor damage state, fraction of intact fuel, vessel lower head failure). The primary intent of this study was to characterize the range of predicted damage states in the 1F1 reactor considering state of knowledge uncertainties associated with MELCOR modeling of core damage progression and to generate information that may be useful in informing the decommissioning activities that will be employed to defuel the damaged reactors at the Fukushima Daiichi Nuclear Power Plant. Additionally, core damage progression variability inherent in MELCOR modeling numerics is investigated.
Efforts are being pursued to develop and qualify a system-level model of a reactor core isolation (RCIC) steam-turbine-driven pump. The model is being developed with the intent of employing it to inform the design of experimental configurations for full-scale RCIC testing. The model is expected to be especially valuable in sizing equipment needed in the testing. An additional intent is to use the model in understanding more fully how RCIC apparently managed to operate far removed from its design envelope in the Fukushima Daiichi Unit 2 accident. RCIC modeling is proceeding along two avenues that are expected to complement each other well. The first avenue is the continued development of the system-level RCIC model that will serve in simulating a full reactor system or full experimental configuration of which a RCIC system is part. The model reasonably represents a RCIC system today, especially given design operating conditions, but lacks specifics that are likely important in representing the off-design conditions a RCIC system might experience in an emergency situation such as a loss of all electrical power. A known specific lacking in the system model, for example, is the efficiency at which a flashing slug of water (as opposed to a concentrated jet of steam) could propel the rotating drive wheel of a RCIC turbine. To address this specific, the second avenue is being pursued wherein computational fluid dynamics (CFD) analyses of such a jet are being carried out. The results of the CFD analyses will thus complement and inform the system modeling. The system modeling will, in turn, complement the CFD analysis by providing the system information needed to impose appropriate boundary conditions on the CFD simulations. The system model will be used to inform the selection of configurations and equipment best suitable of supporting planned RCIC experimental testing. Preliminary investigations with the RCIC model indicate that liquid water ingestion by the turbine decreases the developed turbine torque; the RCIC speed then slows, and thus the pump flow rate to the RPV decreases. Subsequently, RPV water level decreases due to continued boiling and the liquid fraction flowing to the RCIC decreases, thereby accelerating the RCIC and refilling the RPV. The feedback cycle then repeats itself and/or reaches a quasi-steady equilibrium condition. In other words, the water carry-over is limited by cyclic RCIC performance degradation, and hence the system becomes self-regulating. The indications achieved to date with the system model are more qualitative than quantitative. The avenues being pursued to increase the fidelity of the model are expected to add quantitative realism. The end product will be generic in the sense that the RCIC model will be incorporable within the larger reactor coolant system model of any nuclear power plant or experimental configuration.
This report discusses some of the findings from a Fukushima analysis that relays deep insight into critical operating systems such as the RCIC cooling system.
A methodology for using the MELCOR code with the Latin Hypercube Sampling method was developed to estimate uncertainty in various predicted quantities such as hydrogen generation or release of fission products under severe accident conditions. In this case, the emphasis was on estimating the range of hydrogen sources in station blackout conditions in the Sequoyah Ice Condenser plant, taking into account uncertainties in the modeled physics known to affect hydrogen generation. The method uses user-specified likelihood distributions for uncertain model parameters, which may include uncertainties of a stochastic nature, to produce a collection of code calculations, or realizations, characterizing the range of possible outcomes. Forty MELCOR code realizations of Sequoyah were conducted that included 10 uncertain parameters, producing a range of in-vessel hydrogen quantities. The range of total hydrogen produced was approximately 583kg <U+F0B1> 131kg. Sensitivity analyses revealed expected trends with respected to the parameters of greatest importance, however, considerable scatter in results when plotted against any of the uncertain parameters was observed, with no parameter manifesting dominant effects on hydrogen generation. It is concluded that, with respect to the physics parameters investigated, in order to further reduce predicted hydrogen uncertainty, it would be necessary to reduce all physics parameter uncertainties similarly, bearing in mind that some parameters are inherently uncertain within a range. It is suspected that some residual uncertainty associated with modeling complex, coupled and synergistic phenomena, is an inherent aspect of complex systems and cannot be reduced to point value estimates. The probabilistic analyses such as the one demonstrated in this work are important to properly characterize response of complex systems such as severe accident progression in nuclear power plants.
This paper describes the knowledge advancements from the uncertainty analysis for the State-of- the-Art Reactor Consequence Analyses (SOARCA) unmitigated long-term station blackout accident scenario at the Peach Bottom Atomic Power Station. This work assessed key MELCOR and MELCOR Accident Consequence Code System, Version 2 (MACCS2) modeling uncertainties in an integrated fashion to quantify the relative importance of each uncertain input on potential accident progression, radiological releases, and off-site consequences. This quantitative uncertainty analysis provides measures of the effects on consequences, of each of the selected uncertain parameters both individually and in interaction with other parameters. The results measure the model response (e.g., variance in the output) to uncertainty in the selected input. Investigation into the important uncertain parameters in turn yields insights into important phenomena for accident progression and off-site consequences. This uncertainty analysis confirmed the known importance of some parameters, such as failure rate of the Safety Relief Valve in accident progression modeling and the dry deposition velocity in off-site consequence modeling. The analysis also revealed some new insights, such as dependent effect of cesium chemical form for different accident progressions. (auth)
In response to the accident at the Fukushima Daiichi nuclear power station in Japan, the U.S. Nuclear Regulatory Commission and US. Department of Energy agreed to jointly sponsor an accident reconstruction study as a means of assessing the severe accident modeling capability of the MELCOR code and developing an understanding of the likely accident progression. Objectives of the project included reconstruction of the accident progressions using computer models and accident data, and validation of MELCOR and the Fukushima models against plant data. In this study Sandia National Laboratories developed MELCOR 2.1 models of Fukushima Daiichi Units 1 (IFI), 2, and 3 as well as the Unit 4 spent fuel pool. This paper reports on the analysis of the 1F1 accident. Details are presented on the modeled accident progression, hypothesized mode of failures in the reactor pressure vessel (RPV) and containment pressure boundary, and release of fission products to the environment. The MELCOR-predicted RPV and containment pressure trends compare well with available measured pressures. Conditions leading up to the observed explosion of the reactor building are postulated based on this analysis where drywell head flange leakage is thought to have led to accumulation of flammable gases in the refueling bay. The favorable comparison of the results from the analyses with the data from the plant provides additional confidence in MELCOR to reliably predict real-world accident progression. The modeling effort has also provided insights into future data needs for both model development and validation.
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.
Representative accident source terms patterned after the NUREG-1465 Source Term have been developed for high burnup fuel in BWRs and PWRs and for MOX fuel in a PWR with an ice-condenser containment. These source terms have been derived using nonparametric order statistics to develop distributions for the timing of radionuclide release during four accident phases and for release fractions of nine chemical classes of radionuclides as calculated with the MELCOR 1.8.5 accident analysis computer code. The accident phases are those defined in the NUREG-1465 Source Term - gap release, in-vessel release, ex-vessel release, and late in-vessel release. Important differences among the accident source terms derived here and the NUREG-1465 Source Term are not attributable to either fuel burnup or use of MOX fuel. Rather, differences among the source terms are due predominantly to improved understanding of the physics of core meltdown accidents. Heat losses from the degrading reactor core prolong the process of in-vessel release of radionuclides. Improved understanding of the chemistries of tellurium and cesium under reactor accidents changes the predicted behavior characteristics of these radioactive elements relative to what was assumed in the derivation of the NUREG-1465 Source Term. An additional radionuclide chemical class has been defined to account for release of cesium as cesium molybdate which enhances molybdenum release relative to other metallic fission products.
A fission product release and transport model for High Temperature Gas cooled Reactors (HTGRs) is being developed for the MELCOR code. HTGRs use fuel in the form of TRISO coated fuel particles embedded in a graphitized matrix. The HTGR fission product model for MELCOR is being developed to calculate the released amounts and distribution offission products during normal operation and during accidents. The fission product release and transport model considers the important phenomena for fission product behavior in HTGRs, including the recoil and release offission products from the fuel kernel, transport through the coating layers, transport through the surrounding fuel matrix, release into circulating helium coolant, settling and plate-out on structural surfaces, adsorption by graphite dust in the primary system, and resuspension. The fraction of failed particles versus time is input by a particle failure fraction response surface of particle failure fraction as a function offuel temperature, and potentially, fuel burn-up. Fission product release from the fuel kernel and transport through the particle coating layers is calculated using diffusion-based release models. The models account for fission product release from uranium contamination in the graphitized matrix, and adsorption of fission products in the reactor system. The dust and its distribution can be determined from either MELCOR calculations of the reactor system during normal operation, or provided by other sources as input. The distribution of fission products is then normalized using the OR1GEN inventory to provide initial conditions for accident calculations. For the initial releases during an accident, the existing MELCOR aerosol transport models, with appropriate modifications, are being explored for calculating dust and fission product transport in the reactor system and in the confinement. For the delayed releases during the accident, which occur over many hours, and even days fission product release is calculated by combining the diffusion-based release rate with the failure fraction response surface input via a convolution integral. The decay of fission products is also included in the modeling.
Appendix A-5 of Draft Regulatory Guide DG-1199 'Alternative Radiological Source Term for Evaluating Design Basis Accidents at Nuclear Power Reactors' provides guidance - applicable to RADTRAD MSIV leakage models - for scaling containment aerosol concentration to the expected steam dome concentration in order to preserve the simplified use of the Accident Source Term (AST) in assessing containment performance under assumed design basis accident (DBA) conditions. In this study Economic and Safe Boiling Water Reactor (ESBWR) and Advanced Boiling Water Reactor (ABWR) RADTRAD models are developed using the DG-1199, Appendix A-5 guidance. The models were run using RADTRAD v3.03. Low Population Zone (LPZ), control room (CR), and worst-case 2-hr Exclusion Area Boundary (EAB) doses were calculated and compared to the relevant accident dose criteria in 10 CFR 50.67. For the ESBWR, the dose results were all lower than the MSIV leakage doses calculated by General Electric/Hitachi (GEH) in their licensing technical report. There are no comparable ABWR MSIV leakage doses, however, it should be noted that the ABWR doses are lower than the ESBWR doses. In addition, sensitivity cases were evaluated to ascertain the influence/importance of key input parameters/features of the models.
The Oak Ridge National Laboratory computer code, ORIGEN2.2 (CCC-371, 2002), was used to obtain the elemental composition of irradiated low-enriched uranium (LEU)/mixed-oxide (MOX) pressurized-water reactor fuel assemblies. Described in this report are the input parameters for the ORIGEN2.2 calculations. The rationale for performing the ORIGEN2.2 calculation was to generate inventories to be used to populate MELCOR radionuclide classes. Therefore the ORIGEN2.2 output was subsequently manipulated. The procedures performed in this data reduction process are also described herein. A listing of the ORIGEN2.2 input deck for two-cycle MOX is provided in the appendix. The final output from this data reduction process was three tables containing the radionuclide inventories for LEU/MOX in elemental form. Masses, thermal powers, and activities were reported for each category.
The Phebus and VERCORS data have played an important role in contemporary understanding and modeling of fission product release and transport from damaged LWR fuel. The data from these test programs have allowed improvement of MELCOR modeling of release and transport processes for both low enrichment uranium fuel as well as high burnup and MOX fuels. The following paper describes the derivation, testing and incorporation of improved radionuclide release models into the MELCOR severe accident code.
The Phebus and VERCORS data have played an important role in contemporary understanding and modeling of fission product release and transport from damaged light water reactor fuel. The data from these test programs have allowed improvement of MELCOR modeling of release and transport processes for both low enrichment uranium fuel as well as high burnup and mixed oxide (MOX) fuels. This paper discusses the synthesis of these findings in the MELCOR severe accident code. Based on recent assessments of MELCOR 1.8.5 fission product release modeling against the Phebus FPT-1 test and on observations from the ISP-46 exercise, modifications to the default MELCOR 1.8.5 release models are recommended. The assessments identified an alternative set of Booth diffusion parameters recommended by ORNL (ORNL-Booth), which produced significantly improved release predictions for cesium and other fission product groups. Some adjustments to the scaling factors in the ORNL-Booth model were made for selected fission product groups, including UO{sub 2}, Mo and Ru in order to obtain better comparisons with the FPT-1 data. The adjusted model, referred to as 'Modified ORNL-Booth,' was subsequently compared to original ORNL VI fission product release experiments and to more recently performed French VERCORS tests, and the comparisons was as favorable or better than the original CORSOR-M MELCOR default release model. These modified ORNL-Booth parameters, input to MELCOR 1.8.5 as 'sensitivity coefficients' (i.e. user input that over-rides the code defaults) are recommended for the interim period until improved release models can be implemented into MELCOR. For the case of ruthenium release in air-oxidizing conditions, some additional modifications to the Ru class vapor pressure are recommended based on estimates of the RuO{sub 2} vapor pressure over mildly hyperstoichiometric UO{sub 2}. The increased vapor pressure for this class significantly increases the net transport of Ru from the fuel to the gas stream. A formal model is needed. Deposition patterns in the Phebus FPT-1 circuit were also significantly improved by using the modified ORNL-Booth parameters, where retention of lower volatile Cs{sub 2}MoO{sub 4} is now predicted in the heated exit regions of the FPT-1 test, bringing down depositions in the FPT-1 steam generator tube to be in closer alignment with the experimental data. This improvement in 'RCS' deposition behavior preserves the overall correct release of cesium to the containment that was observed even with the default CORSOR-M model. Not correctly treated however is the release and transport of Ag to the FPT-1 containment. A model for Ag release from control rods is presently not available in MELCOR. Lack of this model is thought to be responsible for the underprediction by a factor of two of the total aerosol mass to the FPT-1 containment. It is suggested that this underprediction of airborne mass led to an underprediction of the aerosol agglomeration rate. Underprediction of the agglomeration rate leads to low predictions of the aerosol particle size in comparison to experimentally measured ones. Small particle size leads low predictions of the gravitational settling rate relative to the experimental data. This error, however, is a conservative one in that too-low settling rate would result in a larger source term to the environment. Implementation of an interim Ag release model is currently under study. In the course of this assessment, a review of MELCOR release models was performed and led to the identification of several areas for future improvements to MELCOR. These include upgrading the Booth release model to account for changes in local oxidizing/reducing conditions and including a fuel oxidation model to accommodate effects of fuel stoichiometry. Models such as implemented in the French ELSA code and described by Lewis are considered appropriate for MELCOR. A model for ruthenium release under air oxidizing conditions is also needed and should be included as part of a fuel oxidation model since fuel stoichiometry is a fundamental parameter in determining the vapor pressure of ruthenium oxides over the fuel. There is also a need to expand the MELCOR architecture for tracking fission product classes to allow for more speciation of fission products. An example is the formation of CsI and Cs{sub 2}MoO{sub 4} and potentially CsOH if all Mo is combined with Cs such that excess Cs exists in the fuel. Presently, MELCOR can track only one class combination (CsI) accurately, where excess Cs is assumed to be CsOH. Our recommended interim modifications map the CsOH (MELCOR Radionuclide Class 2) and Mo (Class 7) vapor pressure properties to Cs{sub 2}MoO{sub 4}, which approximates the desired formal class combination of Cs and Mo. Other extensions to handle properly iodine speciation from pool/gas chemistry are also needed.
In this study, risk-significant pressurized-water reactor severe accident sequences are examined using MELCOR 1.8.5 to explore the range of fission product releases to the reactor containment building. Advances in the understanding of fission product release and transport behavior and severe accident progression are used to render best estimate analyses of selected accident sequences. Particular emphasis is placed on estimating the effects of high fuel burnup in contrast with low burnup on fission product releases to the containment. Supporting this emphasis, recent data available on fission product release from high-burnup (HBU) fuel from the French VERCOR project are used in this study. The results of these analyses are treated as samples from a population of accident sequences in order to employ approximate order statistics characterization of the results. These trends and tendencies are then compared to the NUREG-1465 alternative source term prescription used today for regulatory applications. In general, greater differences are observed between the state-of-the-art calculations for either HBU or low-burnup (LBU) fuel and the NUREG-1465 containment release fractions than exist between HBU and LBU release fractions. Current analyses suggest that retention of fission products within the vessel and the reactor coolant system (RCS) are greater than contemplated in the NUREG-1465 prescription, and that, overall, release fractions to the containment are therefore lower across the board in the present analyses than suggested in NUREG-1465. The decreased volatility of Cs2MoO4 compared to CsI or CsOH increases the predicted RCS retention of cesium, and as a result, cesium and iodine do not follow identical behaviors with respect to distribution among vessel, RCS, and containment. With respect to the regulatory alternative source term, greater differences are observed between the NUREG-1465 prescription and both HBU and LBU predictions than exist between HBU and LBU analyses. Additionally, current analyses suggest that the NUREG-1465 release fractions are conservative by about a factor of 2 in terms of release fractions and that release durations for in-vessel and late in-vessel release periods are in fact longer than the NUREG-1465 durations. It is currently planned that a subsequent report will further characterize these results using more refined statistical methods, permitting a more precise reformulation of the NUREG-1465 alternative source term for both LBU and HBU fuels, with the most important finding being that the NUREG-1465 formula appears to embody significant conservatism compared to current best-estimate analyses.
As part of a Nuclear Regulatory Commission (NRC) research program to evaluate the impact of using mixed-oxide (MOX) fuel in commercial nuclear power plants, a study was undertaken to evaluate the impact of the usage of MOX fuel on the consequences of postulated severe accidents. A series of 23 severe accident calculations was performed using MELCOR 1.8.5 for a four-loop Westinghouse reactor with an ice condenser containment. The calculations covered five basic accident classes that were identified as the risk- and consequence-dominant accident sequences in plant-specific probabilistic risk assessments for the McGuire and Catawba nuclear plants, including station blackouts and loss-of-coolant accidents of various sizes, with both early and late containment failures. Ultimately, the results of these MELCOR simulations will be used to provide a supplement to the NRC's alternative source term described in NUREG-1465. Source term magnitude and timing results are presented consistent with the NUREG-1465 format. For each of the severe accident release phases (coolant release, gap release, in-vessel release, ex-vessel release, and late in-vessel release), source term timing information (onset of release and duration) is presented. For all release phases except for the coolant release phase, magnitudes are presented for each of the NUREG-1465 radionuclide groups. MELCOR results showed variation of noble metal releases between those typical of ruthenium (Ru) and those typical of molybdenum (Mo); therefore, results for the noble metals were presented for Ru and Mo separately. The collection of the source term results can be used as the basis to develop a representative source term (across all accident types) that will be the MOX supplement to NUREG-1465.
A simplified ESBWR MELCOR model was developed to track the transport of iodine released from damaged reactor fuel in a hypothesized core damage accident. To account for the effects of iodine pool chemistry, radiolysis of air and cable insulation, and surface coatings (i.e., paint) the iodine pool model in MELCOR was activated. Modifications were made to MELCOR to add sodium pentaborate as a buffer in the iodine pool chemistry model. An issue of specific interest was whether iodine vapor removed from the drywell vapor space by the PCCS heat exchangers would be sequestered in water pools or if it would be rereleased as vapor back into the drywell. As iodine vapor is not included in the deposition models for diffusiophoresis or thermophoresis in current version of MELCOR, a parametric study was conducted to evaluate the impact of a range of iodine removal coefficients in the PCCS heat exchangers. The study found that higher removal coefficients resulted in a lower mass of iodine vapor in the drywell vapor space.
Analyses were performed using MELCOR and RADTRAD to investigate main steam isolation valve (MSIV) leakage behavior under design basis accident (DBA) loss-of-coolant (LOCA) conditions that are presumed to have led to a significant core melt accident. Dose to the control room, site boundary and LPZ are examined using both approaches described in current regulatory guidelines as well as analyses based on best estimate source term and system response. At issue is the current practice of using containment airborne aerosol concentrations as a surrogate for the in-vessel aerosol concentration that exists in the near vicinity of the MSIVs. This study finds current practice using the AST-based containment aerosol concentrations for assessing MSIV leakage is non-conservative and conceptually in error. A methodology is proposed that scales the containment aerosol concentration to the expected vessel concentration in order to preserve the simplified use of the AST in assessing containment performance under assumed DBA conditions. This correction is required during the first two hours of the accident while the gap and early in-vessel source terms are present. It is general practice to assume that at {approx}2hrs, recovery actions to reflood the core will have been successful and that further core damage can be avoided. The analyses performed in this study determine that, after two hours, assuming vessel reflooding has taken place, the containment aerosol concentration can then conservatively be used as the effective source to the leaking MSIV's. Recommendations are provided concerning typical aerosol removal coefficients that can be used in the RADTRAD code to predict source attenuation in the steam lines, and on robust methods of predicting MSIV leakage flows based on measured MSIV leakage performance.
The primary objective of this report is to provide the technical basis for development of recommendations for updates to the NUREG-1465 Source Term for BWRs that will extend its applicability to accidents involving high burnup (HBU) cores. However, a secondary objective is to re-examine the fundamental characteristics of the prescription for fission product release to containment described by NUREG-1465. This secondary objective is motivated by an interest to understand the extent to which research into the release and behaviors of radionuclides under accident conditions has altered best-estimate calculations of the integral response of BWRs to severe core damage sequences and the resulting radiological source terms to containment. This report, therefore, documents specific results of fission product source term analyses that will form the basis for the HBU supplement to NUREG-1465. However, commentary is also provided on observed differences between the composite results of the source term calculations performed here and those reflected NUREG-1465 itself.
Before this LDRD research, no single tool could simulate a very high temperature reactor (VHTR) that is coupled to a secondary system and the sulfur iodine (SI) thermochemistry. Furthermore, the SI chemistry could only be modeled in steady state, typically via flow sheets. Additionally, the MELCOR nuclear reactor analysis code was suitable only for the modeling of light water reactors, not gas-cooled reactors. We extended MELCOR in order to address the above deficiencies. In particular, we developed three VHTR input models, added generalized, modular secondary system components, developed reactor point kinetics, included transient thermochemistry for the most important cycles [SI and the Westinghouse hybrid sulfur], and developed an interactive graphical user interface for full plant visualization. The new tool is called MELCOR-H2, and it allows users to maximize hydrogen and electrical production, as well as enhance overall plant safety. We conducted validation and verification studies on the key models, and showed that the MELCOR-H2 results typically compared to within less than 5% from experimental data, code-to-code comparisons, and/or analytical solutions.
A two-year effort focused on applying ASCI technology developed for the analysis of weapons systems to the state-of-the-art accident analysis of a nuclear reactor system was proposed. The Sandia SIERRA parallel computing platform for ASCI codes includes high-fidelity thermal, fluids, and structural codes whose coupling through SIERRA can be specifically tailored to the particular problem at hand to analyze complex multiphysics problems. Presently, however, the suite lacks several physics modules unique to the analysis of nuclear reactors. The NRC MELCOR code, not presently part of SIERRA, was developed to analyze severe accidents in present-technology reactor systems. We attempted to: (1) evaluate the SIERRA code suite for its current applicability to the analysis of next generation nuclear reactors, and the feasibility of implementing MELCOR models into the SIERRA suite, (2) examine the possibility of augmenting ASCI codes or alternatives by coupling to the MELCOR code, or portions thereof, to address physics particular to nuclear reactor issues, especially those facing next generation reactor designs, and (3) apply the coupled code set to a demonstration problem involving a nuclear reactor system. We were successful in completing the first two in sufficient detail to determine that an extensive demonstration problem was not feasible at this time. In the future, completion of this research would demonstrate the feasibility of performing high fidelity and rapid analyses of safety and design issues needed to support the development of next generation power reactor systems.