Modeling Internal Material Melting
The purpose of the SEED LDRD was to create new functionality to model internal material melting and solidification inside the state-of-the-art nuclear fuel performance simulation code BISON.
SEARCH · Engineering Papers
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.
The purpose of the SEED LDRD was to create new functionality to model internal material melting and solidification inside the state-of-the-art nuclear fuel performance simulation code BISON.
Pronghorn is an engineering-scale, coarse-mesh, thermal-hydraulics tool for supporting reactor-core simulations of advanced nuclear reactors. Most of the current efforts in Pronghorn have been devoted in developing porous finite-volume capabilities and adapting closure correlations for coarse-mesh thermal-hydraulics modeling. However, for liquid-metal reactors (LMRs) with wire-wrapped fuel pin assemblies, a pin-level thermal-hydraulic resolution is required for most safety case studies (pin rupture, channel blockage, etc.). For this purpose, a new Subchannel application is developed in MOOSE, which affords the required flow field resolution, while still preserving an engineering-scale approach. This new solver can be natively coupled to Pronghorn and other MOOSE objects to enable full-core, multi-physics, multi-scale engineering studies. This presentation presents the main features of the SCM code and demonstrates a validation case based on the EBR-II SHRT tests.
The DOE National Nuclear Security Administration (NNSA) Stockpile Stewardship Program (SSP) is an integrated technical program for maintaining the safety, surety, and reliability of the U.S. nuclear stockpile. The SSP incorporates nuclear test data, computational modeling and simulation, and experimental facilities to advance understanding of nuclear weapons. The suite of data analyzed comes from activities including stockpile surveillance, experimental research, and development and engineering programs. This integrated national program requires the continued use of experimental facilities and the computational capabilities to support the SSP missions. These component parts, in addition to an appropriately scaled production capability, enable NNSA to support stockpile requirements. The ultimate goal of the SSP, and thus of the Advanced Simulation and Computing (ASC) Program, is to ensure that the U.S. maintains a safe, secure, and effective strategic deterrent. Specific work activities and scope contained in this Implementation Plan (IP) represent the full-year annual operating plan for FY22. The Initial IP, effective , should be consistent with the Department’s Base Table when operating under a Continuing Resolution (CR). The final IP, effective date TBD, is consistent with the final, enacted appropriation.
To assist both the advanced reactor industry and U.S. Nuclear Regulatory Commission (NRC) in the pursuit of reactor design and licensing, the U.S. Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has established a mechanistic source term (MST) research project under the auspices of the Multiphysics Applications technical area. MST analysis is central to the reactor licensing process and a focus of the NRC regarding their mission to provide reasonable assurance of adequate protection of public health and safety and environment. The MST research project is a collaboration between Argonne National Laboratory (Argonne) and Sandia National Laboratories (SNL) with a high-level objective to coordinate the development of comprehensive advanced reactor MST mod/sim capabilities to support risk-informed design and licensing decisions.
To assist both the advanced reactor industry and U.S. Nuclear Regulatory Commission (NRC) in the pursuit of reactor design and licensing, the U.S. Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has established a mechanistic source term (MST) research project under the auspices of the Application Drivers technical area. The MST research project is a collaboration between Argonne National Laboratory (Argonne) and Sandia National Laboratories (SNL) with a high-level objective to coordinate the development of comprehensive advanced reactor MST mod/sim capabilities to support risk-informed design and licensing decisions. As part of this effort, an MST mod/sim development pathway was developed in FY21, which outlines the high-level objectives and near-term tasks necessary to achieve the project objectives. Since publication of the development pathway, Argonne and SNL have been addressing the “near-term tasks” outlined in the report. The current provides a status update of the progress achieved in the fourth quarter of FY21. The report structure follows that of ref and is divided by advanced reactor type: High temperature gas reactor (HTGR), molten salt reactor (MSR) including fluoride salt-cooled high temperature reactor (FHR), and sodium fast reactor (SFR). Each section reviews the near-term tasks associated with the reactor type, current progress, and future plans. The near-term tasks associated with consequence modeling are reviewed in an upcoming work by SNL.
To assist both the advanced reactor industry and U.S. Nuclear Regulatory Commission (NRC) in the pursuit of reactor design and licensing, the U.S. Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has established a mechanistic source term (MST) research project under the auspices of the Multiphysics Applications technical area. MST analysis is central to the reactor licensing process and a focus of the NRC regarding their mission to provide reasonable assurance of adequate protection of public health and safety and environment. The MST research project is a collaboration between Argonne National Laboratory (Argonne) and Sandia National Laboratories (SNL) with a high-level objective to coordinate the development of comprehensive advanced reactor MST mod/sim capabilities to support risk-informed design and licensing decisions. As part of this effort, an MST mod/sim development pathway was developed in FY21, which outlines the high-level objectives and near-term tasks necessary to achieve the project objectives. Since publication of the development pathway, Argonne and SNL have focused initial efforts on addressing the “near-term tasks” outlined in the report, and the current work provides a status update of the progress achieved in FY22. This report is a progress update on 6 tasks supporting molten salt reactor (MSR) and sodium fast reactor (SFR) MST analyses which were initiated in FY22 as recommended high priority near-term tasks. These tasks were continued into FY23 as they also support previously recommended long-term development objectives to advance MST analyses for these advanced reactors.
Nuclear data are a vital component of predictive simulations used in applications like experiment design, stockpile stewardship, nuclear nonproliferation/safeguards, health physics, and criticality safety. A singular simulation requires the coalescence of different areas of nuclear data such as cross sections, angular distributions, and energy distributions of emitted neutrons for different materials and energy ranges. Improving nuclear data and thus reducing the uncertainty in simulated parameters could enable smaller, better-informed safety factors and ultimately reduce operational and procedural costs. There is a constant effort to garner a better understanding of the physical quantities represented by nuclear data through experiments. Integral experiment benchmarks use simulated and measured results to validate current nuclear data values. In the past, benchmarks primarily focused on the effective multiplication factor (k eff ); however, this limited scope has caused compensating errors and areas of nuclear data that lack validation. Compensating errors are inaccuracies in nuclear data that are obfuscated by cancellation when observing integrated values such as k eff . Diverse integral benchmark experiments that look for quantities of interest other than k eff and include multiple responses minimize the possibility of compensating errors and provides validation to areas of nuclear data previously lacking experimental validation. Benchmark experiments can be optimized during the design process to be highly dependent on specific areas of nuclear data. The dependence of a response in an experiment to a specific area/type of nuclear data is defined as sensitivity. A larger sensitivity means that nuclear data uncertainties will play a larger role in the response(s) resulting in larger bias. Currently, the sensitivity capabilities of the Monte Carlo N-Particle (MCNP ®1 ) transport code are limited to responses of k eff and tallied values (e.g., flux, surface current). As a part of the EUCLID project, this work explores estimating list-mode nuclear data sensitivities that can be used to design experiments aimed to constrain and reduce compensating errors in nuclear data by focusing on responses other than k eff . Tallied values are ideal quantities that are estimated with detectors during experiments. List-mode data (a list of neutron collection times) are the direct output of detector systems in subcritical neutron noise experiments. Expanding MCNP sensitivity capabilities to include the sensitivity of responses estimated from list-mode data, such as the prompt neutron decay constant (α) and multiplicity estimates (S and D), enables more direct comparison of simulated and measured experimental quantities. Additionally, deterministic tools such as SENSMG are capable of obtaining sensitivities to a wide variety of responses; however, these tools cannot handle complex geometries due to the assumptions made in discretizing the phase-space variables of the Boltzman transport equation.
Nuclear power plant components operate in environments that promote multiple degradation mecha- nisms, several of which involve crack initiation and growth. An ongoing effort in the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is developing a general capability within the Multiphysics Object Oriented Simulation Environment (MOOSE) framework for simulating three-dimensional crack growth under a range of driving conditions, including fatigue, stress corrosion cracking (SCC), brittle fracture, and stress-relaxation cracking. This report demonstrates an end-to-end workflow that uses this capability to model weld-residual-stress-driven SCC in the J-groove weld of a pressurized-water reactor control rod drive mechanism penetration in the vessel head. The workflow consists of a thermomechanical welding simulation with temperature-dependent plasticity, followed by cooldown to ambient conditions, and a restart of the simulation using the MOOSE extended finite element method (XFEM) module to propagate a three-dimensional crack through the residual stress field. New welding capabilities were developed to properly initialize newly activated elements in the weld region, and robustness improvements were made to the mesh-based algorithm for defining cutting planes in the 3D XFEM algorithm, allowing it to handle complex crack fronts and stress fields. Together these advances allowed the simulated SCC crack to grow from an initial elliptical flaw in the weld, across the weld, through the tube wall, and almost to the triple point (where the weld, tube, and reactor pressure vessel head intersect) over roughly 36 years of simulated service. These results demonstrate a workflow that can be extended to fully three-dimensional welding simulations and more complex crack interaction problems.
The paper, drawing on the example of simulation codes used in nuclear physics and high-energy physics, seeks to call attention to the ethical implications of discontinuing support for simulation codes and the loss of knowledge embodied in them. Predicated on the concept of trading zones and actor network models, the paper addresses the problem of extinction of simulation codes and attempts to understand their evolution and development within those frameworks. We show that simulation codes of closed type develop to the level of creoles, becoming local languages and standards of scientific centers and disappearing as their few main developers leave, whereas codes of open types become universal languages, imposing problemsolving patterns on the entire community and crowding out other codes. The paper suggests that because of simulations’ reliance on tacit knowledge, practices entrenched in codes cannot be exhaustively explicated or transmitted through writing alone; on the contrary, the life cycle of a simulation code is determined by the life cycle of its trading zone. We examine the extent to which both of these phenomena pose a risk to the preservation of knowledge. Bearing upon intergenerational ethics, we draw analogies between the pure intergenerational problem (PIP) and the problem of preserving the knowledge implemented in simulation codes and transmitting it to future generations. We argue that for the complete transfer of knowledge, it is necessary to develop and maintain inhabitability and sustainability of simulation trading zones in a controllable way, at least until the demand for these codes is warranted to cease in the future.
Experiments were conducted to develop test configurations and technology necessary to simulate the thermal environment and fuel region expected to exist in in-reactor tests of small models of nuclear light bulb configurations. Particular emphasis was directed at rf plasma tests of approximately full-scale models of an in-reactor cell suitable for tests in Los Alamos Scientific Laboratory's Nuclear Furnace. The in-reactor tests will involve vortex-stabilized fissioning uranium plasmas of approximately 200-kW power, 500-atm pressure and equivalent black-body radiating temperatures between 3220 and 3510 K.
The Pacific Northwest National Laboratory (PNNL) is conducting confirmatory research for the U.S. Nuclear Regulatory Commission (NRC) to evaluate commercially available nondestructive examination (NDE) modeling and simulation software used in the nuclear industry. Simulation results from ultrasonic testing (UT) models can inform the design and qualification of inspection techniques and help interpret inspection results. CIVA is a modeling and simulation package developed by the French Alternative Energies and Atomic Energy Commission (CEA). CIVA was selected for this study because it is readily available and has been used for NDE in the US commercial nuclear power industry. This report is focused on completing the efforts initiated in the previous PNNL report to evaluate UT modeling and simulation performance, reliability, and accuracy in relation to common inservice inspection (ISI) scenarios in nuclear power plants (NPP). This work will be used to provide guidance when establishing methods to perform and evaluate simulations for more standardized model implementation, simulation analysis, and interpretation of results.
The Short-Baseline Near Detector (SBND) is the first Liquid Argon Time Projection Chamber (LArTPC) with high enough resolution and large enough neutrino flux to measure the cross-section for production of heavier-than-proton fragments in neutrino interactions. SBND is an 112-ton LArTPC, and lies 110 m downstream from the Booster Neutrino Beam target, where it is collecting more than two million neutrino interactions per year. As a result of both the intranuclear cascade and nuclear de-excitation, neutrino interactions with argon can produce nuclear clusters such as deuterons, tritons, helions and alpha particles, which have so far only been seen as vertex activity or mis-reconstructed as protons. With LArTPC technology's excellent reconstruction capabilities and SBND's unprecedented neutrino statistics, a measurement of these nuclear clusters could distinguish between nuclear models and improve neutrino energy reconstruction. This poster presents the status of simulation and reconstruction of nuclear clusters in SBND.
To improve the economics of commercial nuclear energy generation, U.S. utilities are currently seeking licensing approval to operate UO2 fuel to higher burnups. One significant safety issue that must be addressed to obtain approval is the potential for fine fragmentation/pulverization of the fuel during a loss-of-coolant accident (LOCA). The cause of pulverization has been hypothesized to be the rapid increase of pressure in fission gas bubbles in the high burnup region of the fuel during a LOCA. To better understand this phenomenon, a novel phase-field model of the fission gas bubble microstructure in UO2 has been developed and implemented in Idaho National Laboratory (INL)'s Marmot application for phase-field simulation of nuclear materials. Simulations of the bubble response to steady-state and transient conditions were conducted. Simulation results were used to inform a mechanistic model of pulverization in BISON, INL’s fuel performance simulation code.
The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) framework underpins most of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) physics applications and coupling methods. Users interact with MOOSE in multiple ways to enable the construction of complex multiphysics models, including but not limited to compilation, input creation and syntax validation, meshing, solving the physics problem via MOOSE-based solvers and coupling, and output inspection. In FY23, numerous enhancements have been made to MOOSE to enhance the user experience. Reactor-oriented meshing capabilities in MOOSE have been expanded, an online tutorial was created and hosted on the MOOSE site, and a hands-on workshop was held for more than 100 users. The Language Server Protocol (LSP) capability has been implemented in MOOSE to better communicate correct syntax to the NEAMS Workbench user interface which is commonly used to validate input and submit jobs by users. Finally, several other user-facing improvements were made including the implementation of improved screen output and logging control, enhancements to the initial condition system, and enhancements of the MultiApp system commonly used for coupling. Collectively, these enhancements were driven by user needs and directly improve user ability to construct complex multiphysics models for nuclear reactors.
Microgravity induces inflammatory responses and modulates immune functions that may increase oxidative stress. Exposure to a microgravity environment induces adverse neurological effects; however, there is little research exploring the etiology of these effects resulting from exposure to such an environment. It is also known that spaceflight is associated with increase in oxidative stress; however, this phenomenon has not been reproduced in land-based simulated microgravity models. In this study, an attempt has been made to show the induction of reactive oxygen species (ROS) in mice brain, using ground-based microgravity simulator. Increased ROS was observed in brain stem and frontal cortex with concomitant decrease in glutathione, on exposing mice to simulated microgravity for 7 d. Oxidative stress-induced activation of nuclear factor-kappaB was observed in all the regions of the brain. Moreover, mitogen-activated protein kinase kinase was phosphorylated equally in all regions of the brain exposed to simulated microgravity. These results suggest that exposure of brain to simulated microgravity can induce expression of certain transcription factors, and these have been earlier argued to be oxidative stress dependent.
The multi-physics STAR-CCM+ code has been used to simulate the temperature and flow fields in canister gas leakage experiments conducted by using a 1/4.5-scale model storage cask. The simulations were conducted at Argonne National Laboratory’s Laboratory Computing Resource Center, utilizing high-performance computing resources. Development of STAR-CCM+ simulation models for the 1/4.5-scale model cask is described herein, followed by validation of the simulation results against experimental data. Canister depressurization and thermal response during gas leakage are discussed, along with analyses of the leakage path and allowable leakage rate. The insights gained from the STAR-CCM+ simulations and leakage analyses will help guide future experiments and actual industry applications of Argonne’s Remote Area Modular Monitoring system for canister surface temperature measurement (RAMM-TM) to enable effective aging management of spent fuel during extended dry storage, as well as help reduce risks to public safety and health and protect the environment.
Recirculation zone in fluid mechanics experiments on cylindrical cavity simulating gas-core nuclear rocket as heavy gas reservoir, discussing flow patterns
We present 3D general-relativistic neutrino-radiation hydrodynamics simulations of two asymmetric binary neutron star mergers producing long-lived neutron stars remnants and spanning a fraction of their cooling time scale. The mergers are characterized by significant tidal disruption with neutron rich material forming a massive disc around the remnant. The latter develops one-armed dynamics that is imprinted in the emitted kilo-Hertz gravitational waves. Angular momentum transport to the disc is initially driven by spiral-density waves and enhanced by turbulent viscosity and neutrino heating on longer timescales. The mass outflows are composed by neutron-rich dynamical ejecta of mass ∼10 −3 –10 −2 M ⊙ followed by a persistent spiral-wave/neutrino-driven wind of ≳ 10 −2 M ⊙ with material spanning a wide range of electron fractions, ∼0.1–0.55. Dynamical ejecta (winds) have fast velocity tails up to ∼0.8 (∼0.4) c. The outflows are further evolved to days timescale using 2D ray-by-ray radiation-hydrodynamics simulations that include an online nuclear network. We find complete r-process yields and identify the production of 56 Ni and the subsequent decay chain to 56 Co and 56 Fe. Synthetic kilonova light curves predict an extended (near-) infrared peak a few days postmerger originating from r-process in the neutron-rich/high-opacity ejecta and UV/optical peaks at a few hours (ten minutes) postmerger originating from weak r-process (free-neutron decay) in the faster ejecta components. Additionally, the fast tail of tidal origin generates kilonova afterglows potentially detectable in radio and X band on a few to ten years time scale. Quantitative effects originating from the tidal disruption merger dynamics are reflected in the multimessenger emissions.