Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RELAP”

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 109 records · Page 6

DireWolf

DireWolf is an integrated micro-reactor capability being developed using an advanced multiphysics modeling and simulation capability for heat pipe-cooled micro-reactor concepts. The DireWolf application is built on INL’s Multiphysics Object-Oriented Simulation Environment (MOOSE) development and runtime framework and is composed of the following MOOSE- based applications as submodules, coupled together using the MOOSE MultiApps and Transfers capability and executed as a single application: Griffin (reactor physics and radiation transport), BISON (nuclear fuel performance), Grizzly (structural mechanics), Sockeye (heat pipe technology), and RELAP-7 (open-air and super-critical CO2 Brayton cycle).

Permann, Cody↗

Dynamic Event Tree advancements and control logic improvements

The RAVEN code has been under development at the Idaho National Laboratory since 2012. Its main goal is to create a multi-purpose platform for the deploying of all the capabilities needed for Probabilistic Risk Assessment, uncertainty quantification, data mining analysis and optimization studies. RAVEN is currently equipped with three different sampling categories: Forward samplers (Monte Carlo, Latin Hyper Cube, Stratified, Grid Sampler, Factorials, etc.), Adaptive Samplers (Limit Surface search, Adaptive Polynomial Chaos, etc.) and Dynamic Event Tree (DET) samplers (Deterministic and Adaptive Dynamic Event Trees). The main subject of this document is to report the activities that have been done in order to: start the migration of the RAVEN/RELAP-7 control logic system into MOOSE, and develop advanced dynamic sampling capabilities based on the Dynamic Event Tree approach. In order to provide to all MOOSE-based applications a control logic capability, in this Fiscal Year an initial migration activity has been initiated, moving the control logic system, designed for RELAP-7 by the RAVEN team, into the MOOSE framework. In this document, a brief explanation of what has been done is going to be reported. The second and most important subject of this report is about the development of a Dynamic Event Tree (DET) sampler named “Hybrid Dynamic Event Tree” (HDET) and its Adaptive variant “Adaptive Hybrid Dynamic Event Tree” (AHDET). As other authors have already reported, among the different types of uncertainties, it is possible to discern two principle types: aleatory and epistemic uncertainties. The classical Dynamic Event Tree is in charge of treating the first class (aleatory) uncertainties; the dependence of the probabilistic risk assessment and analysis on the epistemic uncertainties are treated by an initial Monte Carlo sampling (MCDET). From each Monte Carlo sample, a DET analysis is run (in total, N trees). The Monte Carlo employs a pre-sampling of the input space characterized by epistemic uncertainties. The consequent Dynamic Event Tree performs the exploration of the aleatory space. In the RAVEN code, a more general approach has been developed, not limiting the exploration of the epistemic space through a Monte Carlo method but using all the forward sampling strategies RAVEN currently employs. The user can combine a Latin Hyper Cube, Grid, Stratified and Monte Carlo sampling in order to explore the epistemic space, without any limitation. From this pre-sampling, the Dynamic Event Tree sampler starts its aleatory space exploration. As reported by the authors, the Dynamic Event Tree is a good fit to develop a goal-oriented sampling strategy. The DET is used to drive a Limit Surface search. The methodology that has been developed by the authors last year, performs a Limit Surface search in the aleatory space only. This report documents how this approach has been extended in order to consider the epistemic space interacting with the Hybrid Dynamic Event Tree methodology.

97 MATHEMATICS AND COMPUTING↗

Year End Progress Report on Rattlesnake Improvements

Rattlesnake is a MOOSE-based radiation transport application developed at INL to support modern multi-physics simulations. At the beginning of the last year, Rattlesnake was able to perform steady-state, transient and eigenvalue calculations for the multigroup radiation transport equations. Various discretization schemes, including continuous finite element method (FEM) with discrete ordinates method (SN) and spherical harmonics expansion method (PN) for the self-adjoint angular flux (SAAF) formulation, continuous FEM (CFEM) with SN for the least square (LS) formulation, diffusion approximation with CFEM and discontinuous FEM (DFEM), have been implemented. A separate toolkit, YAKXS, for multigroup cross section management was developed to support Rattlesnake calculations with feedback both from changes in the field variables, such as fuel temperature, coolant density, and etc., and in isotope inventory. The framework for doing nonlinear diffusion acceleration (NDA) within Rattlesnake has been set up, and both NDA calculations with SAAF-SN-CFEM scheme and Monte Carlo with OpenMC have been performed. It was also used for coupling BISON and RELAP-7 for the full-core multiphysics simulations. Within the last fiscal year, significant improvements have been made in Rattlesnake. Rattlesnake development was migrated into our internal GITLAB development environment at the end of year 2014. Since then total 369 merge requests has been accepted into Rattlesnake. It is noted that the MOOSE framework that Rattlesnake is based on is under continuous developments. Improvements made in MOOSE can improve the Rattlesnake. It is acknowledged that MOOSE developers spent efforts on patching Rattlesnake for the improvements made on the framework side. This report will not cover the code restructuring for better readability and modularity and documentation improvements, which we have spent tremendous effort on. It only details some of improvements in the following sections.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Improved Sampling Algorithms in the Risk-Informed Safety Margin Characterization Toolkit

The RISMC approach is developing advanced set of methodologies and algorithms in order to perform Probabilistic Risk Analyses (PRAs). In contrast to classical PRA methods, which are based on Event-Tree and Fault-Tree methods, the RISMC approach largely employs system simulator codes applied to stochastic analysis tools. The basic idea is to randomly perturb (by employing sampling algorithms) timing and sequencing of events and internal parameters of the system codes (i.e., uncertain parameters) in order to estimate stochastic parameters such as core damage probability. This approach applied to complex systems such as nuclear power plants requires to perform a series of computationally expensive simulation runs given a large set of uncertain parameters. These types of analysis are affected by two issues. Firstly, the space of the possible solutions (a.k.a., the issue space or the response surface) can be sampled only very sparsely, and this precludes the ability to fully analyze the impact of uncertainties on the system dynamics. Secondly, large amounts of data are generated and tools to generate knowledge from such data sets are not yet available. This report focuses on the first issue and in particular employs novel methods that optimize the information generated by the sampling process by sampling unexplored and risk-significant regions of the issue space: adaptive (smart) sampling algorithms. They infer system response from surrogate models constructed from existing samples and predict the most relevant location of the next sample. It is therefore possible to understand features of the issue space with a small number of carefully selected samples. In this report, we will present how it is possible to perform adaptive sampling using the RISMC toolkit and highlight the advantages compared to more classical sampling approaches such Monte-Carlo. We will employ RAVEN to perform such statistical analyses using both analytical cases but also another RISMC code: RELAP-7.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Reduced Order Model Implementation in the Risk-Informed Safety Margin Characterization Toolkit

The RISMC project aims to develop new advanced simulation-based tools to perform Probabilistic Risk Analysis (PRA) for the existing fleet of U.S. nuclear power plants (NPPs). These tools numerically model not only the thermo-hydraulic behavior of the reactor primary and secondary systems but also external events temporal evolution and components/system ageing. Thus, this is not only a multi-physics problem but also a multi-scale problem (both spatial, µm-mm-m, and temporal, ms-s-minutes-years). As part of the RISMC PRA approach, a large amount of computationally expensive simulation runs are required. An important aspect is that even though computational power is regularly growing, the overall computational cost of a RISMC analysis may be not viable for certain cases. A solution that is being evaluated is the use of reduce order modeling techniques. During the FY2015, we investigated and applied reduced order modeling techniques to decrease the RICM analysis computational cost by decreasing the number of simulations runs to perform and employ surrogate models instead of the actual simulation codes. This report focuses on the use of reduced order modeling techniques that can be applied to any RISMC analysis to generate, analyze and visualize data. In particular, we focus on surrogate models that approximate the simulation results but in a much faster time (µs instead of hours/days). We apply reduced order and surrogate modeling techniques to several RISMC types of analyses using RAVEN and RELAP-7 and show the advantages that can be gained.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

R&D Plan for RISMC Industry Application #1: ECCS/LOCA Cladding Acceptance Criteria

The Nuclear Regulatory Commission (NRC) is finalizing a rulemaking change that would revise the requirements in 10 CFR 50.46. In the proposed new rulemaking, designated as 10 CFR 50.46c, the NRC proposes a fuel performance-based equivalent cladding reacted (ECR) criterion as a function of cladding hydrogen content before the accident (pre-transient) in order to include the effects of higher burnup on cladding performance as well as to address other technical issues. A loss of operational margin may result due to the more restrictive cladding embrittlement criteria. Initial and future compliance with the rule may significantly increase vendor workload and licensee costs as a spectrum of fuel rod initial burnup states may need to be analyzed to demonstrate compliance. The Idaho National Laboratory (INL) has initiated a project, as part of the DOE Light Water Reactor Sustainability Program (LWRS), to develop analytical capabilities to support the industry in the transition to the new rule. This project is called the Industry Application 1 (IA1) within the Risk-Informed Safety Margin Characterization (RISMC) Pathway of LWRS. The general idea behind the initiative is the development of an Integrated Evaluation Model (IEM). The motivation is to develop a multiphysics framework to analyze how uncertainties are propagated across the stream of physical disciplines and data involved, as well as how risks are evaluated in a LOCA safety analysis as regulated under 10 CFR 50.46c. This IEM is called LOTUS which stands for LOCA Toolkit for US, and it represents the LWRS Program’s response to the proposed new rule making. The focus of this report is to complete an R&D plan to describe the demonstration of the LOCA/ECCS RISMC Industry Application # 1 using the advanced RISMC Toolkit and methodologies. This report includes the description and development plan for a RISMC LOCA tool that fully couples advanced MOOSE tools already in development in order to characterize and optimize plant safety and operational margins. Advanced MOOSE tools that are needed to complete this integrated evaluation model are: RAVEN, RELAP-7, BISON, and MAMMOTH.

10CFR50.46C↗

Data Analysis Approaches for the Risk-Informed Safety Margins Characterization Toolkit

In the past decades, several numerical simulation codes have been employed to simulate accident dynamics (e.g., RELAP5-3D, RELAP-7, MELCOR, MAAP). In order to evaluate the impact of uncertainties into accident dynamics, several stochastic methodologies have been coupled with these codes. These stochastic methods range from classical Monte-Carlo and Latin Hypercube sampling to stochastic polynomial methods. Similar approaches have been introduced into the risk and safety community where stochastic methods (such as RAVEN, ADAPT, MCDET, ADS) have been coupled with safety analysis codes in order to evaluate the safety impact of timing and sequencing of events. These approaches are usually called Dynamic PRA or simulation-based PRA methods. These uncertainties and safety methods usually generate a large number of simulation runs (database storage may be on the order of gigabytes or higher). The scope of this paper is to present a broad overview of methods and algorithms that can be used to analyze and extract information from large data sets containing time dependent data. In this context, “extracting information” means constructing input-output correlations, finding commonalities, and identifying outliers. Some of the algorithms presented here have been developed or are under development within the RAVEN statistical framework.

97 MATHEMATICS AND COMPUTING↗

IAEA CRP on HTGR Uncertainties in Modeling: Assessment of Phase I Lattice to Core Model Uncertainties

Best-estimate plus uncertainty analysis of reactors is replacing the traditional conservative (stacked uncertainty) method for safety and licensing analysis. To facilitate uncertainty analysis applications, a comprehensive approach and methodology must be developed and applied. High temperature gas cooled reactors (HTGRs) have several features that require techniques not used in light-water reactor analysis (e.g., coated-particle design and large graphite quantities at high temperatures). The International Atomic Energy Agency has therefore launched the Coordinated Research Project on HTGR Uncertainty Analysis in Modeling to study uncertainty propagation in the HTGR analysis chain. The benchmark problem defined for the prismatic design is represented by the General Atomics Modular HTGR 350. The main focus of this report is the compilation and discussion of the results obtained for various permutations of Exercise I 2c and the use of the cross section data in Exercise II 1a of the prismatic benchmark, which is defined as the last and first steps of the lattice and core simulation phases, respectively. The report summarizes the Idaho National Laboratory (INL) best estimate results obtained for Exercise I 2a (fresh single-fuel block), Exercise I 2b (depleted single-fuel block), and Exercise I 2c (super cell) in addition to the first results of an investigation into the cross section generation effects for the super-cell problem. The two dimensional deterministic code known as the New ESC based Weighting Transport (NEWT) included in the Standardized Computer Analyses for Licensing Evaluation (SCALE) 6.1.2 package was used for the cross section evaluation, and the results obtained were compared to the three dimensional stochastic SCALE module KENO VI. The NEWT cross section libraries were generated for several permutations of the current benchmark super-cell geometry and were then provided as input to the Phase II core calculation of the stand alone neutronics Exercise II 1a. The steady state core calculations were simulated with the INL coupled-code system known as the Parallel and Highly Innovative Simulation for INL Code System (PHISICS) and the system thermal-hydraulics code known as the Reactor Excursion and Leak Analysis Program (RELAP) 5 3D using the nuclear data libraries previously generated with NEWT. It was observed that significant differences in terms of multiplication factor and neutron flux exist between the various permutations of the Phase I super-cell lattice calculations. The use of these cross section libraries only leads to minor changes in the Phase II core simulation results for fresh fuel but shows significantly larger discrepancies for spent fuel cores. Furthermore, large incongruities were found between the SCALE NEWT and KENO VI results for the super cells, and while some trends could be identified, a final conclusion on this issue could not yet be reached. This report will be revised in mid 2016 with more detailed analyses of the super-cell problems and their effects on the core models, using the latest version of SCALE (6.2). The super-cell models seem to show substantial improvements in terms of neutron flux as compared to single-block models, particularly at thermal energies.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Initial Testing of the Microscopic Depletion Implementation in the MAMMOTH Reactor Physics Application

Present and new nuclear fuels that will be tested at the Transient Reactor Test (TREAT) facility will be analyzed with the MAMMOTH reactor physics application, currently under development, at Idaho National Laboratory. MAMMOTH natively couples the BISON, RELAP-7, and Rattlesnake applications within the MOOSE framework. This system allows the irradiation of fuel from beginning of life in a nuclear reactor until it is placed in TREAT for fuel testing within the same analysis mesh and, thus, retaining a very high level of resolution and fidelity. The calculation of the isotopic distribution in fuel requires the solution to the decay and transmutation equations coupled to the neutron transport equation. The Chebyshev Rational Approximation Method (CRAM) is the current state-of-the-art in the field, as was chosen to be the solver for the decay and transmutation equations. This report shows that the implementation of the CRAM solver within MAMMOTH is correct with various analytic benchmarks for decay and transmutation of nuclides. The results indicate that the solutions with CRAM order 16 achieve the level of precision of the benchmark. The CRAM solutions show little sensitivity to the time step size and consistently produce a high level of accuracy for isotopic decay for time steps of 1x10^11 years. Comparisons to DRAGON5 with 297 isotopes yield comparable results, but some differences need to be further analyzed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Long-term Prediction of Emissivity of Structural materials for High Temperature Reactor Systems

The purpose of this research was to extend the recent short-term studies of emissivity of reactor materials performed independently at the University of Missouri and University of Wisconsin to tests aimed at longerterm predictions of emissivity. The University of Missouri measured hemispherical emissivity, which provides an integrated measurement, while the University of Wisconsin measured spectral emissivity (2 to 10 μm wavelength) from which the integrated quantity can be obtained. Hemispherical emissivity is needed for MELCOR, RELAP and FLUENT type heat transfer code calculations, while spectral emissivity is needed for a detailed Monte Carlo Transport theory computation of heat transfer and fundamental modeling verification. The objectives of the study were: Objective 1: Defining materials for long-term emissivity in cross-cutting reactor applications - procurement, preparation, and analyses of materials. Objective 2: High temperature materials testing in various environments and post-test surface analysis Objective 3: Independent hemispherical and spectral emissivity testing. Objective 4: Development of models to gain a predictive capability for long-term emissivity and application and validation in current RCCS research programs. The research has been performed at two universities over a four-year term with the participation of students, scientific staff, and faculty. The report describes progress in realizing the objectives, and the accomplishments at the two universities. There was considerable progress with respect to Objectives 1-3. The model development was redirected towards modeling of interrelationships between the spectra and hemispherical emissivities and uncertainty analysis of total hemispherical emissivity measurements.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

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↗

Sockeye: A 1-D Heat Pipe Modeling Tool

Sockeye is a 1-D MOOSE-based application for modeling heat pipes. Sockeye is aimed at modeling heat pipes with annular wick and contains two-phase fluid property packages for sodium and potassium working fluids. Sockeye adapts the two-phase flow model from RELAP-7 to the simulation of heat pipes, by accounting for the capillary pressure in pressure relaxation terms. Additionally, a simple, aggregate heat transfer model exists for transitioning from temperatures below the melting temperature of the working fluid. Results are presented for an example test problem, and some future work is discussed.

42 ENGINEERING↗

Architecture Improvements and Retesting of the Microscopic Depletion Implementation in the MAMMOTH Reactor Physics Application

Nuclear fuel experiments conducted at the Transient Reactor Test Facility (TREAT) will be analyzed with the MAMMOTH reactor physics application, currently under development at Idaho National Laboratory. MAMMOTH natively couples with the BISON, RELAP-7, and Rattlesnake applications within the MOOSE framework. The implemented methods allow for simulating the irradiation of fuel in a nuclear reactor from the beginning of life in a nuclear reactor until it is placed in TREAT for fuel testing using the same finite element mesh throughout the analysis retaining a very high level of resolution and fidelity. The calculation of the isotopic distribution in fuel requires the solution to the decay and transmutation equations coupled to the neutron transport equation. This report focuses on three significant improvements to the MAMMOTH micro-depletion capability. First, the depletion architecture within MAMMOTH is modified to improve the flexibility and speed of burnup calculations. The results of various benchmark tests demonstrate that the refactoring of the previously implemented capability does not introduce any regression. The second improvement involves the expansion of the data functionality. MAMMOTH can now read the decay data files used by the Oak Ridge Isotope Generation (ORIGEN) code, thus enabling access to data of more than 1,600 isotopes. The third improvement is the capability to perform depletion at a constant power level which was verified via a comparison to ORIGEN.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

FY21 Status report on the CMVB and CNWG International Collaborations

The ART-GCR Methods area includes an international collaboration work package that covers the tasks defined for the Computational Methods Validation and Benchmark (CMVB)and Civil Nuclear Energy Research and Development Working Group (CNWG) projects. This report summarizes the status of the FY21 tasks and planned FY22 DOE contributions. The CMVB Project Arrangement (PA) is not yet formally approved by all signatories, and no work has therefore been performed in FY21 at INL related to this activity. The CNWG activities inFY21 consisted of the simulation of the High Temperature Test Reactor (HTTR) Loss Of Forced Cooling (LOFC) experimental with the Idaho National Laboratory (INL) codes Griffin, BISON and RELAP-7 based on the Multiphysics Object-Oriented Simulation Environment (MOOSE).It was found that the multiphysics coupled suite is capable of simulating all the important phenomena occurring during the Depressurized Loss of Forced Cooling (DLOFC) experiments showing promising agreement with the measurements given the number of uncertainties and approximations introduced into the model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

FY22 Status Report on the ART-GCR CMVB and CNWG International Collaborations

This work presents the numerical model for the high temperature test reactor (HTTR) loss of forced cooling (LOFC) experiment with the INL codes Griffin, BISON, and RELAP-7, based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Promising results were obtained, with the overall behavior of the reactor successfully captured. Changes in the heat transfer coupling, as compared to the fiscal year (FY)-21 model, enabled drastic improvement of the steady-state solution, both in terms of computational time and global energy discrepancies. The former was reduced by a factor of roughly 60, whereas the latter decreased from 11% to less than 2%. Furthermore, the discrepancy in the steady-state multiplication factor was improved from +2,300 and +2,900 pcm (for the 30 and 9 MW cases) to -700 and +1,200 pcm, respectively, and now falls well within the large measurement uncertainties stemming from graphite impurities. Validation of the Monte Carlo model used to generate cross sections was also performed against available measurements. Though significant, the discrepancies remain acceptable overall in light of the large uncertainty. Numerous improvements are still needed to better compare with the experiment involving the 9 MW case and to instill greater confidence in the model’s ability to accurately predict the 30 MW behavior. Specifically, the power levels predicted by the 9 MW transient simulation following re-criticality remain low, pointing to an underestimation of the passive cooling of the core. A key aspect of future work will be to better understand the flow pattern during the LOFC event, particularly to determine if natural or forced convection is occurring inside the reactor pressure vessel (RPV). More generally, additional validation data would be immensely useful for further enhancing the numerical model and better matching the experiments. In addition, a more sophisticated thermal-hydraulics model that simulates all the channels as a single system model should be considered to take into account the rest of the primary loop. Finally, even if the results are in better agreement with the experiments, sensitivity analysis and uncertainty quantification will be necessary to evaluate the model uncertainty.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Accuracy Enhancement of Nuclear Power Plant Simulators Utilizing High Accuracy Simulation Predictions

More recently, reactor core simulators for core designs associated with commercial nuclear power plants that utilize what is believed to be higher fidelity models have been developed. Features such as neutronics models that utilize transport equation solvers with fine spatial meshes and many energy-groups, thermal-hydraulic models that utilize sub-channel solvers with fine spatial mesh and capable of treating a wide range of fluid conditions, and fuel-coolant chemistry interaction models capable of treating CRUD deposition are to be found in these higher fidelity core simulators. These reactor core simulators require access to higher performance computers, characterized by many processors, cores and large memory. So associated with utilization of these simulators is access to high performance computers and ability to accommodate in one’s workflow longer execution times. By contrast, currently used core simulators by the nuclear industry can execute on engineering workstations and have execution times of seconds to minutes. The desirability for having short execution times is not only desired for support of time critical tasks but supports the mental process of decision making by engineers. The goal of the work reported upon here has the objective of retaining the fidelity of higher fidelity models while retaining the ability to utilize engineering workstations. Beyond the core simulator goal, additional goals of this work include incorporating the just described core simulator capability into a Nuclear Steam Supply System (NSSS) simulator, and to incorporate the resulting capability into an environment supportive of design and operational decision making associated with nuclear power stations. The model selected for the core neutronics model is the NESTLE code, for the core thermal-hydraulic model is the CTF code utilizing coarse mesh, and for the NSSS model is the RELAP5-3D code. WSC’s proprietary 3KEYMASTERTM platform is being used to provide software coupling, user interface, visualization, and reporting. The NESTLE core neutronics simulator was first integrated with the CTF core thermal-hydraulic simulator using CTF developed communication commands which are also used for CTF to communicate with RELAP5-3D under WSC’s proprietary 3KEYMASTERTM platform. To assure NESTLE prediction consistency with higher fidelity core neutronic simulators, buffer codes have been created to automatically generate from output files written by the VERA core simulator the NESTLE nodal neutronic parameter’ library, geometry, and pin-power reconstruction input files, thereby avoiding a number of challenges associated with utilizing lattice physics codes and providing consistency with VERA predictions. To treat absorber rod effects a multi-set library is utilized, where a set refers to a specific absorber rod fully inserted pattern. A coarse spatial mesh CTF model was developed with features added that support using CTF as envisioned in the engineering quality simulator. A hybrid meshing approach was implemented to allow for automated construction of models with mixed levels of refinement. Specifically, a core model could resolve some assemblies at a nodal level (4 subchannels per assembly) and others at a pin-resolution (one subchannel per coolant subchannel in the assembly). The intention is that this will allow for better resolution of limiting conditions such as DNBR and PCT, which are based on local rod and subchannel conditions. Further development was done of features that enhance the capabilities for the envisioned engineering quality simulator that has been developed, but now for RELAP-3D. The RELAP5-3D code development includes ability to model more than 999 components and the addition of the cross-channels turbulence mixing model and the void drift model that are implemented in CTF, aiming to achieve closer prediction agreement of the two codes for transient simulations, specifically, more accurate matches of the overall mass, momentum, and energy exchanges of both the liquid and gas phases between the neighboring core assemblies. Graphics were also developed for the Instructor Station for this project under WSC’s proprietary 3KEYMASTERTM platform to facilitate design and operational decision making.

42 ENGINEERING↗

Towards a NEAMS-based high-fidelity model of the MARVEL reactor

This report outlines the progress of Idaho National Laboratory in developing a high-fidelity and high-resolution model of the Microreactor Applications Research Validation and Evaluation reactor. The model was developed under the Nuclear Energy Advanced Modeling and Simulation microreactor application driver at Idaho National Laboratory. The overarching objective of this activity is the development of a high-fidelity multiphysics MARVEL model using NEAMS tools, and to verify and validate NEAMS tools against MARVEL reference simulation and experimental data, respectively. This is a unique opportunity to conduct multiphysics analysis on a soon-to-be-deployed microreactor. This multiphysics model developed under the NEAMS-funded INL microreactor application driver leverages three single-physics models coupled via the MOOSE’s MultiApp and Transfer systems. The latter systems enable in-memory data transfer between MOOSE-based and MOOSE-wrapped applications. The first single-physics model, that functions as main application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). Several optimization flags that were developed by the Griffin developer team were beta-tested to enhance the solver’s performance. These include the combined use of using_average_xs and update_averaged_xs_on that enable to avoid expensive on-the-fly cross sections evaluations at each linear iterations in favor of evaluations of the macroscopic cross sections at each Picard iteration. The second single-physics model uses BISON to handle solid heat transfer and asymptotic hydrogen redistribution analysis in the fuel. While the model returns consistent results for the temperature and hydrogen distribution in the fuel, a mismatch was noticed in the calculated temperature in the reflector due to the value of the gap conductance used in our model. Ongoing investigations are being performed to assess the origin of this discrepancy. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. A first verification was also performed showing good agreement in terms of mass flow rate and inlet temperature. All mesh files were generated using the MOOSE Reactor module, removing the need for external meshing tools. Notably, this workscope represents one of the initial applications of the MOOSE Reactor module for modeling highly irregular geometries. The use of the reactor module significantly streamlined the mesh generation process. The full multiphysics mode, that combines all the single physics models, was leveraged to conduct initial steady-state multiphysics simulations to compute power, and temperature distribution in the reactor. Initial testing was performed for transient simulations as well. In this case, the new checkpoint restart capability for eigenvalue calculations was tested showing the capability for streamlined restart of transient calculations. Future work will focus on improving the fidelity of the model by performing comprehensive code-to-code comparisons. For instance, the full-core Griffin neutronics model will be benchmarked against MCNP reference results, that were provided by the MARVEL design team. Additionally, the SAM T/H model will be verified against reference RELAP-5 results for selected accident scenarios. Besides code-to-code verification exercises, the model fidelity will be improved by replacing the single-channel SAM model with a more complex SAM-Pronghorn coupled model, in which the sub-channel capability is deployed to obtain radial temperature resolution in the coolant. This model will be developed in synergy with the NEAMS thermal hydraulics team.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Apparatus and method for lapping an edge surface of an object

An apparatus for lapping an edge surface of an object comprises a block having a side adapted to engage a wide surface of an object, adjustable spacers disposed on the block and adapted to engage a lap plate, and a weighted spring disposed on the block for urging the spacers and the object edge surface towards the lap plate. A method for lapping comprises setting surfaces of adjustable spacers disposed on a block to be substantially the same distance from the block, affixing a wide surface of an object to the block, urging an edge surface of the object and the spacers towards a lap plate, lapping the edge of the object, inspecting the edge for parallelism to a reference line, resetting the spacers and relapping the edge surface.

Rossi, Vito N.↗