Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analysis codes”

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 127 records · Page 7

Preliminary Analysis of a Light Water Ingress Accident for the French High Flux Reactor

The consequences of a light water ingress accident scenario have been analyzed for the French High Flux Reactor, an heavy-water cooled, moderated and reflected reactor. It is found that the use of safety analysis code such as RELAP or CATHARE are required and neutronic tool alone cannot provide an accurate description of all phenomena involved. It is also found that even a loose coupling between neutronic and safety analysis code can be implemented and used to analyze the accident scenario considered. As a result, despite the substantial reactivity increase induced by the ingress of light water in the core, it is concluded that the RHF reactor would remain safe (peak cladding temperature well below melting temperature) in any conditions considered and at any time.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

U.S. Efforts in Support of Examinations at Fukushima Daiichi: 2019 Evaluations

Much is still not known about the end-state of core materials in each unit that was operating on March 11, 2011 at the Fukushima Daiichi Nuclear Power Station (Daiichi). Information obtained from Daiichi is required to inform Decontamination and Decommissioning (D&D) activities, improving the ability of the Tokyo Electric Power Company Holdings, Incorporated (TEPCO Holdings) to characterize potential hazards and to ensure the safety of workers involved with cleanup activities. This document, which has been updated to include Fiscal Year 2019 (FY2019) information, summarizes results from a U.S. effort to review Daiichi information and extract insights to enhance the safety of existing and future nuclear power plant designs. This U.S. effort, which was initiated in 2014 by the Department of Energy Office of Nuclear Energy (DOE-NE), is completed by a group of experts in reactor safety and plant operations that identify examination needs and evaluate recent Daiichi examination data to address these needs. Each year, annual reports are issued that document significant safety insights being obtained in areas of special emphasis. For FY2019, the areas emphasized by the expert panel include component performance, radionuclide surveys and sampling, debris end-state location, combustible gas effects, and plant operations and maintenance. In addition to reducing uncertainties related to severe accident modeling progression, these insights have and continue to be used to update guidance for severe accident prevention, mitigation, and emergency planning. Reduced uncertainties in modeling the events at Daiichi improve the realism of reactor safety evaluations that inform future D&D activities. In addition, the process established by this U.S. effort provides an example that Japan may wish to implement in future international efforts to support D&D. As emphasized in this FY2019 report, important insights are being gained from examinations of the affected units at Daiichi. Although many results from prior years have been moved to appendices, key areas of interest are highlighted in the main document. For example, larger intact components and fuel containing materials observed on the floor of the primary containment and retained on ex-vessel structures continue to be emphasized because of insights that may be gained with respect to vessel failure and ex-vessel interactions with structures and concrete. In addition, the main report presents insights gained by focused evaluations of plant data because of potential use of systems and components during severe accident conditions and alternate perspectives regarding accident progression that may be inferred from such data. The report also illustrates how systems analysis codes can be applied to reduce uncertainties about the progression of events at Daiichi and how practical insights obtained from Daiichi can be used to enhance radiation safety for the existing fleet. Recent and planned updates to systems analysis codes are also discussed. As in prior years, this report contains an appendix with a list of examination needs that was updated after U.S. experts reviewed recently obtained information from Daiichi. Additional details for higher priority, near-term, examination activities are also provided. This report also includes an appendix with a description of a website that was developed to assist U.S. experts by providing information in an archived retrievable location.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

High-Temperature Gas-Cooled Pebble-Bed Reactors Running In And Transient Modeling Capabilities Demonstration

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor (HTGR) applications. The first part of this effort focuses on the running-in verification of Griffin's multiphysics capabilities, specifically for simulating the evolution of Pebble Bed reactor cores from startup to equilibrium. In the absence of validation data, code-to-code comparisons are conducted with Kugelpy, showing good agreement for key quantities like maximum power density and fresh core k-eigenvalue predictions. However, discrepancies in equilibrium core predictions suggest potential issues with cross sections, underscoring the need for further refinement and evaluation. The HTTF system analysis code benchmark involves RELAP5-3D, SAM, and GAMMA+ to assess their predictive capabilities for HTTF behavior under both normal operation and pressurized conduction cooldown (PCC) transient conditions. While there is good agreement in predicting major parameters such as coolant temperature, solid temperature, and flow distribution, discrepancies in transient behavior highlight differences in modeling approaches, nodalizations, and heat transfer models. The HTTF lower plenum CFD benchmark employs nekRS to simulate flow mixing phenomena, successfully capturing relevant flow physics and demonstrating mesh independence in complex geometries. Preliminary results suggest a relatively uniform temperature field but significant unsteadiness in the flow, requiring time-averaging analyses. The GPBR200 system analysis code benchmark uses SAM's core channel and porous media models, incorporating an RCCS loop for decay heat removal. During steady-state and transient conditions, including protected de-pressurized and pressurized loss of forced cooling (DLOFC and PLOFC), both models show good agreement in predicting temperature profiles and key parameters. Notably, while the core channel model underpredicts convective heat transfer effects, both models maintain temperatures well below the TRISO fuel safety limit. These benchmarking efforts collectively enhance the predictive capabilities of the tools used in HTGR design and safety analysis, guiding developments to improve their accuracy and applicability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Process metallurgy simulation for metal drawing process optimization by using two-scale finite element method

We developed two-scale FE analysis procedure based on the crystallographic homogenization method by considering the hierarchical structure of poly-crystal aluminium alloy metal. It can be characterized as the combination of two-scale structure, such as the microscopic polycrystal structure and the macroscopic elastic plastic continuum. Micro polycrystal structure can be modeled as a three dimensional representative volume element (RVE). RVE is featured as by 3×3×3 eight-nodes solid finite elements, which has 216 crystal orientations. This FE analysis code can predict the deformation, strain and stress evolutions in the wire drawing processes in the macro- scales, and further the crystal texture and hardening evolutions in the micro-scale. In this study, we analyzed the texture evolution in the wire drawing processes by our two-scale FE analysis code under conditions of various drawing angles of dice. We evaluates the texture evolution in the surface and center regions of the wire cross section, and to clarify the effects of processing conditions on the texture evolution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Covariance Applications in Criticality Safety, Light Water Reactor Analysis, and Spent Fuel Characterization

A new covariance data library based on ENDF/B-VII.1 was recently processed for the SCALE nuclear analysis code system. The multigroup covariance data are discussed here, along with testing and application results for critical benchmark experiments. The cross section covariance library, along with covariances for fission product yields and decay data, is used to compute uncertainties in the decay heat produced by a burned reactor fuel assembly.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Statistical Study of Primary Knock-on Atoms under Neutron Irradiation in RMC Code

Materials in high flux neutron condition will suffer from an evolution of damage which ultimately restricts the lifetime of the fuel clad and reactor pressure vessel. Atoms in irradiated material transform to primary knock-on atoms (PKAs) by absorbing the kinetic energy transferred from neutrons, and the energy released from nuclear reactions. The secondary collision cascades which are induced by PKAs create point defects within the material, with residual vacancies and interstitials eventually causing swelling and embrittlement that degrades mechanical behavior. For a long time, Displacements per Atom (DPA) index has been identified as a measure of radiation damage in material. Nowadays, several advanced methods, such as Molecular Dynamics (MD) and Kinetic Monte Carlo (KMC) simulations, have been developed for more accurately simulating the evolution of the damage. However, all of these techniques need more comprehensive information of PKA, such as the type, energy, direction and spatial distribution. The limitations of DPA emerge and a method should be provided to simulate the generative process of PKAs. Recently, a deterministic analysis code SPECTRA-PKA has been developed by Gilbert to calculate the energy spectra of almost all types PKAs by solving the group-to-group cross-section matrices equation which includes neutron scattering, light charged particles, and recoils of the residual nuclei. But this deterministic analysis code could not give the direction, spatial and time distribution of PKAs. Based on the Geant4 program, Cui Zhenguo has simulated the irradiation by neutrons with 1 MeV energy in zirconium, while only the two adjacent elastic collisions between injected neutron and target atoms are taken into consideration. In this paper, Monte Carlo method to calculate PKAs information during neutron transport in RMC code is introduced. RMC is a continuous-energy Reactor Monte Carlo neutron and photon transport code being developed by Department of Engineering Physics at Tsinghua University, Beijing. PKAs in zirconium are calculated in fixed neutron source mode. Statistical information of PKAs in PWR fuel cell and SFR fuel cell is shown and compared. The results show that Monte Carlo method can provide more comprehensive information of PKAs. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Experimental study on the guillotine break accident of the main steam line with Atlas

A main steam line break (MSLB) accident is defined as a pipe break in the main steam system of a nuclear power plant. The MSLB accident is characterized as cool-down event due to the increased steam flow rate, which causes excessive energy removal from the steam generators and eventually reactor coolant system. This results in a decrease in reactor coolant temperature. The cool-down causes an increase in core reactivity due to the negative moderator and Doppler reactivity coefficients. The thermalhydraulic integral effect test was performed to simulate a double-ended guillotine break accident of the main steam line of a pressurized water reactor. With an aim of simulating the MSLB accident as realistically as possible, a pertinent scaling approach was taken, especially focusing on a break flow. The main objectives of this test were not only to provide physical insight into the system response of a nuclear power plant during the MSLB accident but also to produce an integral effect test data to validate the safety analysis codes. In the present test, asymmetric thermal-hydraulic behavior was clearly observed between the loops. In the cold leg-1A and -1B where the affected steam generator is connected, the flow rates were higher and the fluid temperatures were lower than those in the cold leg-2A and -2B. This could be attributed to the excessive heat removal in the cold leg-1A and -1B through the affected steam generator. This integral effect test data will be used to evaluate the prediction capability of existing safety analysis codes. Furthermore, this data can be utilized to identify any code deficiency for a MSLB simulation, especially for DVI-adapted plants. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Covariance Structure of Pu Isotopic Compositions Obtained Using FRAM 5.2

In the nondestructive assay of special nuclear materials, isotopic determination based on the analysis of high-purity Ge high-resolution gamma-ray spectra plays a pivotal role. The isotopic results, which are based on a relative measurement of peak intensities, are used to interpret data acquired using quantitative (mass) nondestructive assay methods. In the case of plutonium measurements, this includes both passive and active neutron techniques, and also calorimetry (heat). Historically, precision estimates on the Pu- isotopic composition have been generated, on a case-by-case basis, by least squares fitting algorithms embedded in the analysis codes, but all other sources of uncertainty have generally been estimated using a top-down approach by item category. A largely neglected consideration is the covariance between the reported abundances and derived quantities – notably Pu-239-effective mass, Pu-240-effective mass and specific thermal power. In this work we study the covariance structure by statistically evaluating repeat spectrum measurements of two known low burn-up heavily filtered Pu items reduced using the spectrum analysis code FRAM 5.2 and alternate standard parameter sets. For routine measurements, however, only a single spectrum from an item will typically be acquired. We therefore discuss how to generate the sought after covariance information from an individual spectrum by spectrum re-sampling, an approach which could be added to provide item specific reports, and compare the two approaches to test its reliability.

Croft, Stephen↗

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↗

Application of the new IAPWS Guideline on the fast and accurate calculation of steam and water properties with the Spline-Based Table Look-Up Method (SBTL) in RELAP-7

The numerical simulation of thermalhydraulic processes in nuclear power plants requires very accurate and extremely fast algorithms for calculating the thermophysical properties of water and steam. In order to provide such algorithms, the International Association for the Properties of Water and Steam (IAPWS) has adopted the new “IAPWS Guideline on the Fast Calculation of Steam and Water Properties with the Spline-Based Table Look-Up Method (SBTL)”. In this article, the SBTL method is applied to property functions of specific volume and specific internal energy (v,e) based on the scientific formulation IAPWS-95 and the latest IAPWS formulations for transport properties. From the newly generated SBTL functions, thermodynamic and transport properties as well as their derivatives and inverse functions are calculable in the fluid range of state for pressures up to 100 MPa and for temperatures up to 1273 K, including the metastable liquid and the metastable vapor regions. The SBTL functions reproduce the underlying formulations with an accuracy of 10–100 ppm and significantly reduced computing times. The SBTL method has been implemented into the nuclear reactor system safety analysis code RELAP-7 [2] to consider the real fluid behavior of water and steam in a novel 7-equation two-phase flow model.

42 ENGINEERING↗

Magnet Design Considerations for Fusion Nuclear Science Facility

The Fusion Nuclear Science Facility (FNSF) is a nuclear confinement facility that provides a fusion environment with components of the reactor integrated together to bridge the technical gaps of burning plasma and nuclear science between the International Thermonuclear Experimental Reactor (ITER) and the demonstration power plant (DEMO). Compared with ITER, the FNSF is smaller in size but generates much higher magnetic field, i.e., 30 times higher neutron fluence with three orders of magnitude longer plasma operation at higher operating temperatures for structures surrounding the plasma. Input parameters to the magnet design from system code analysis include magnetic field of 7.5 T at the plasma center with a plasma major radius of 4.8 m and a minor radius of 1.2 m and a peak field of 15.5 T on the toroidal field (TF) coils for the FNSF. Both low-temperature superconductors (LTS) and high-temperature superconductors (HTS) are considered for the FNSF magnet design based on the state-of-the-art fusion magnet technology. The higher magnetic field can be achieved by using the high-performance ternary restacked-rod process Nb3Sn strands for TF magnets. The circular cable-in-conduit conductor (CICC) design similar to ITER magnets and a high-aspect-ratio rectangular CICC design are evaluated for FNSF magnets, but low-activation-jacket materials may need to be selected. The conductor design concept and TF coil winding pack composition and dimension based on the horizontal maintenance schemes are discussed. Neutron radiation limits for the LTS and HTS superconductors and electrical insulation materials are also reviewed based on the available materials previously tested. The material radiation limits for FNSF magnets are defined as part of the conceptual design studies for FNSF magnets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dataset_for_Conserved_macromolecular_architecture_of_Poplar_secondary_cell_walls_revealed_by_ssNMR_and_atomistic_modeling

This dataset contains solid-state 13C NMR data and atomistic molecular dynamics simulation files supporting the study of nanoscale secondary cell wall architecture across 13 genetically diverse Populus trichocarpa genotypes grown under uniform greenhouse conditions in 13C-enriched CO2 atmospheres (~89% 13C enrichment).The dataset contains two collections of solid-state 13C NMR data. (1) 200 MHz data (Bruker Avance III HD, 4 mm HX probe, 10 kHz MAS): raw Bruker TopSpin experiment folders and DMFIT-exported ascii spectra for selective and non-selective 1D 13C-13C spin diffusion experiments (3000 ms mixing) used to quantify inter-polymer spatial proximities, and short-mixing (1 ms) reference spectra used for polymeric abundance quantification by spectral deconvolution. (2) 600 MHz data (Bruker Avance III, 1.6 mm PhoenixNMR HXY probe, 30 kHz MAS): raw Bruker TopSpin experiment folders containing 2D CORD, 2D CP-INADEQUATE, and 13C/1H relaxation (T1, T1rho) experiments for all 13 genotypes, with processed Excel workbooks per experiment type. Molecular dynamics simulation code, coordinate files, and analysis scripts (NAMD/CHARMM/Python) for six atomistic cell wall models are included. Summarized ssNMR data are compiled into a single excel file and subjected to statistical analysis. Multivariate analysis code (PCA, Pearson correlation) and summary data are provided as excel worksheets and Jupyter notebooks (Python 3).

09 BIOMASS FUELS↗

Undergraduate education at the University of Kentucky accelerator laboratory

The University of Kentucky accelerator laboratory provides undergraduate students an opportunity to experience scientific research in nuclear physics. They are introduced to our research program, which focuses on the examination of basic nuclear properties using the 7 MV CN Van de Graaff accelerator and the time-of-flight facility to perform inelastic neutron scattering measurements with g-ray and neutron detection. Furthermore, students are given a series of projects to support the measurements including experimental equipment setup, kinematic calculations, data collection, experimental diagnostics, MCNP simulations, analysis code development, and preliminary analysis of the data. Undergraduate physics majors from the University of Dallas took part in (n,n'g) measurements for LiF. Neutron scattering cross sections for these nuclei, which are of significant interest in coolant designs for molten salt reactors, were determined for incident neutron energies between 0.8 and 4.5 MeV. By measuring the g-ray production cross sections and feeding contributions from higher-lying levels, students extracted the neutron level cross sections for Li and F. Additional (n,n'g) measurements for Li metal will be performed in the near future to supplement the g-ray data.

43 PARTICLE ACCELERATORS↗

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↗

Simulated hydrological dynamics and coupled iron redox cycling impact methane production in an Arctic soil: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic publication "Simulated hydrological dynamics and coupled iron redox cycling impact methane production in an Arctic soil" in the Journal of Geophysical Research-Biogeosciences. We simulated biogeochemical cycling in arctic soils using the PFLOTRAN geochemical model combined with measurements from previous NGEE Arctic incubations of polygonal permafrost soils in northern Alaska (Zheng et al., 2018). Simulated iron cycling, carbon dioxide production, and methane production were compared with incubation measurements and the parameterized model was then used to simulate coupled iron and carbon cycling over repeated oxic-anoxic cycles at different levels of carbon substrate availability and pH. The most recent data version (2.0) in the archive incorporates changes to the model and simulations as suggested by reviewers during the manuscript review process. These changes include an updated parameterization of the model; a new set of simulations omitting the iron cycle for direct evaluation of how iron cycle processes affect modeled outcomes; and a set of simulations testing different scenarios of carbon substrate availability in addition to scenarios of initial soil pH. This archive contains simulation code, model output, and analysis code for PFLOTRAN simulations. All scripts are python except the batch script for submitting multiprocessor jobs. Note that the model also requires compiled versions of the Alquimia interface and the NGEE Arctic fork of the PFLOTRAN geochemical simulator (see the README_INSTALL document for basic instructions). The Output directory contains eight data files in netCDF format generated by the model. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 10-year research effort (2012-2022) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗