Mesoscale Multiphysics Phase Field Simulator (MEMPHIS): Recent progresses in phase-field modeling at Sandia National Laboratories.
Abstract not provided.
SEARCH · Search NASA
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.
Abstract not provided.
ENRICO is a coupled application developed under the U.S. Department of Energy's Exascale Computing Project (ECP) targeting the modeling of advanced nuclear reactors. It couples radiation transport with heat and fluid simulation, including the high-fidelity, highresolution Monte-Carlo code Shift and the Computational fluid dynamics code NekRS. NekRS is a highly-performant open-source code for simulation of incompressible and low-Mach fluid flow, heat transfer, and combustion with a particular focus on turbulent flows in complex domains. It is based on rapidly convergent high-order spectral element discretizations that feature minimal numerical dissipation and dispersion. State-of-the-art multilevel preconditioners, efficient high-order time-splitting methods, and runtime-adaptive communication strategies are built on a fast OCCA-based kernel library, libParanumal, to provide scalability and portability across the spectrum of current and future high-performance computing platforms. On Frontier, Nek5000/RS has recently achieved an unprecedented milestone in breaching over 1 billion spectral elements and 350 billion degrees of freedom. Shift has demonstrated the capability to transport upwards of 1 billion particles per second in full core nuclear reactor simulations featuring complete temperature-dependent, continuous-energy physics on Frontier. Shift achieved a weak-scaling efficiency of 97.8% on 8192 nodes of Frontier and calculated 6 reactions in 214,896 fuel pin regions below 1% statistical error yielding first-of-a-kind resolution for a Monte Carlo transport application.
An overview of Idaho National Laboratory's MOOSE: Multiphysics Object-Oriented Simulation Environment.
Multiphysics modeling capabilities have a crucial role to play in the accelerated deployment of fusion energy. To that end, we developed new multiphysics fusion blanket simulation capabilities in the Multiphysics Object-Oriented Simulation Environment (MOOSE). Firstly, we expanded on the existing capabilities of the previously published work, by coupling 3D tritium transport modeling capabilities using the Tritium Migration Analysis Program, version 8 (TMAP8) to an existing tool including thermal hydraulics, fully three-dimensional (3D) heat transfer, and loosely coupled neutronics analysis. Secondly, we performed a thorough verification of the new capabilities and increased testing code coverage to meet MOOSE’s software quality standards. The MOOSE framework follows a strict software quality assurance plan to be Nuclear Quality Assurance, Level 1 compliant. The new multiphysics fusion blanket simulation capabilities are now held to the same standard. Thirdly, to demonstrate MOOSE’s new fusion blanket modeling capabilities, we performed a fully integrated, multiphysics simulation of a 3D solid ceramic breeder blanket design. This proof-of-concept simulation provides the temperature and tritium distribution across the blanket. In conclusion, the combined efforts towards software quality and the development of multiphysics coupling capabilities provide an effective and reliable framework for modeling solid ceramic fusion blankets using MOOSE.
Griffin is a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). This work discloses the methods, algorithms, and implementation for simulating heterogeneous reactor dynamics models. Griffin utilizes a discontinuous finite-element method with discrete ordinates (DFEM-S ) to discretize the field variable of the multigroup neutron transport equation. Multiphysics feedback is handled using two-step tabulated cross-section methodology. Feedback quantities are evaluated using the MOOSE-MultiApp system to couple various engineering phenomena, such as heat conduction and thermal fluids. The multiphysics DFEM-S system is solved using fixed-point iteration with a fully asynchronous parallel sweeper, unstructured coarse-mesh finite difference acceleration, and a multi-timescale improved quasi-static method scheme. The implementation is applied to a multiphysics microreactor model, with two transients: one initiated by a single heat-pipe failure and another by control drum rotation. Importantly, these examples demonstrate the ability of Griffin to tractably solve the neutron transport equation considering seven independent variables and feedback.
As the goal of achieving fusion power on the grid comes closer to fruition, fully coupled multiphysics models of fusion devices will be crucial. These models must incorporate the interconnected phenomena of these devices, including plasma physics, neutronics, first wall interactions, and tritium transport. Currently, there are two main approaches to developing these platforms: (1) loosely coupled, where one couples existing codes and solvers together through input and output parameters and data, and (2) tightly coupled, where one develops the necessary models within a singular, integrated framework. This work focuses on the latter approach for magnetically confined fusion devices by developing a fluid-based plasma-edge model within the Multiphysics Object Oriented Simulation Environment (MOOSE) Framework. This effort is coordinated with other efforts to develop, test, demonstrate, and deploy fusion relevant multiphysics capabilities including electromagnetics, particle-in-cell plasma, tritium transport, and fusion blanket design. This new model is an expansion of the MOOSE-based plasma application, Zapdos, which was originally formulated to model low-temperature, non-magnetized plasma processes. Verification studies have been conducted using newly developed magnetic plasma capabilities. These involved convergence analyses utilizing the method of manufactured solutions to verify new operators and case studies. A modular approach was taken here to demonstrate increasingly complicated simulation scenarios, which included a singular fluid with uniform magnetic field case, a singular fluid with spatially varying magnetic field case, and a coupled multifluid case.
In this work we present a new approach to solving a multiphysics steady-state simulation for the High Temperature Test Reactor (HTTR) using one Super Homogenization (SPH) corrected homogenized 3-D neutronics model in MAMMOTH, multiple heterogeneous 2-D axisymmetric fuel pin models in BISON and several 1-D annular flow thermohydraulics models in RELAP-7. The coupling of the different physics is complicated by the fundamentally disparate models, potentially affecting the stability of the global system. A pseudo-transient to converge the temperature fields is used to circumvent this difficulty. The benefit of doing so as opposed to having fully 3-D heterogenous models for all physics is that the computational burden for reactors as complex as HTTR may quickly become unbearable while coarse models may not capture all the essential phenomena of the problem. The driving purpose of this work is to establish a self-sustaining steady-state solution that can be used to start safety transients.
In this work we present a new approach to solving a multiphysics steady-state simulation for the High Temperature Test Reactor (HTTR) using one Super Homogenization (SPH) corrected homogenized 3-D neutronics model in MAMMOTH, multiple heterogeneous 2-D axisymmetric fuel pin models in BISON and several 1-D annular flow thermohydraulics models in RELAP-7. The coupling of the different physics is complicated by the fundamentally disparate models, potentially affecting the stability of the global system. A pseudo-transient to converge the temperature fields is used to circumvent this difficulty. The benefit of doing so as opposed to having fully 3-D heterogenous models for all physics is that the computational burden for reactors as complex as HTTR may quickly become unbearable while coarse models may not capture all the essential phenomena of the problem. The driving purpose of this work is to establish a self-sustaining steady-state solution that can be used to start safety transients.
In this work we present a new approach to solving a multiphysics steady-state simulation for the High Temperature Test Reactor (HTTR) using one Super Homogenization (SPH) corrected homogenized 3-D neutronics model in MAMMOTH, multiple heterogeneous 2-D axisymmetric fuel pin models in BISON and several 1-D annular flow thermohydraulics models in RELAP-7. The coupling of the different physics is complicated by the fundamentally disparate models, potentially affecting the stability of the global system. A pseudo-transient to converge the temperature fields is used to circumvent this difficulty. The benefit of doing so as opposed to having fully 3-D heterogenous models for all physics is that the computational burden for reactors as complex as HTTR may quickly become unbearable while coarse models may not capture all the essential phenomena of the problem. The driving purpose of this work is to establish a self-sustaining steady-state solution that can be used to start safety transients.
Core bowing due to thermal gradients and irradiation induced swelling and creep introduces significant reactivity feedback effects in liquid metal-cooled fast spectrum reactors. During startup, normal operations, and transient events, differential thermal and flux gradients cause expansion, creep, and swelling, which results in bowing in preferential directions depending on local material properties and load-pad and restraint-ring design. The bowing phenomenon produces negative reactivity during accident conditions provided the restraint system has been properly designed to optimally guide the deformation outwards in the active core region. Under the Department of Energy Nuclear Energy Advanced Modeling and Simulation program, a Multiphysics Object Oriented Simulation Environment (MOOSE)-based multiphysics approach to model core bowing is being developed. The present work expands on previous modeling of simpler International Atomic Energy Agency (IAEA) verification problems with these tools by modeling IAEA Verification Problem 4, which involves a symmetric sector of a reactor core with ducted assemblies undergoing differential thermal expansion due to thermal gradients bowing outward with duct-to-duct contact. This model is available on the National Reactor Innovation Center Virtual Test Bed repository. This example was verified against benchmark participant results, which includes bowing deformation evaluation and duct-to-duct interactions at load pads with mechanical contact.
Here, this study identifies the phenomena that impact the accuracy of the multiphysics analysis of fluoride salt-cooled high temperature reactors (FHR) and categorizes the level of coupling needed to model these phenomena with sufficient accuracy in the context of three scenarios, namely, normal operation, station blackout, and simultaneous withdrawal of all control rods. Such a study is a necessary step prior to developing or extending existing multiphysics modeling and simulation tools. The identified multiphysics phenomena are categorized into the groups requiring ‘‘tight” and ‘‘loose” coupling. ‘‘Tight” coupling refers to phenomena that require iterative feedback between multiple physics modules/methods whereas ‘‘loose” coupling refers to phenomena that require only data from other physics that can be pre-calculated and shared in tabular or functional form. For the station blackout scenario, it was determined that no identified phenomena required tight coupling. For the rod withdrawal and normal operation scenarios, it is found that there is a large number of phenomena that require tight coupling between neutronics and thermal hydraulics.
With the recent development of advanced numerical algorithms, software design, and low-cost high-performance computer hardware, reliance on coupled multiphysics to predict the behavior of complex physical systems is beginning to become standard practice. This is especially true in nuclear energy applications where strong nonlinear interdependencies exist between reactor physics, radiation transport, multi-scale nuclear fuels performance, thermal fluids, etc. Resolving these nonlinear dependencies requires choices in multiphysics software approaches. Two main multiphysics modeling and simulation approaches have emerged. The first is based upon "code coupling" where disparate physics codes of different software design, code languages, and spatial and temporal integration schemes are coupled together with relatively complex data passing interfaces. The second multiphysics software approach is to employ a "cohesive" framework where all physics applications are developed with a common software design, i.e., data structures, syntax, input format, integrated spatial and temporal discretization schemes, etc. In this paper we present the Multiphysics Object-Oriented Simulation Environment (MOOSE) development and runtime framework and describe the framework's cohesive modeling and simulation multiphysics approach. Then, a "cohesive-like" extension of the MOOSE framework is presented where MOOSE-based physics software applications are efficiently coupled to non-MOOSE (external) physics codes to form multiphysics applications using MOOSE's unique interface capabilities. Finally, several examples of MOOSE's cohesive and cohesive-like multiphysics applications will be demonstrated. These multiphysics demonstrations will incorporate both MOOSE-based applications and external codes, including Nek5000, RELAP-7, TRACE, BISON, and Pronghorn.
The Multiphysics Object-Oriented Simulation Environment (MOOSE) Framework, as well as MOOSE-based simulation tools, have accelerated the development of fission energy and advanced reactor technologies through the United States Department of Energy, Office of Nuclear Science, Nuclear Energy Advanced Modeling & Simulation (NEAMS) Program. MOOSE contains a complete platform of multiphysics simulation capabilities, capable of running on massively parallel systems, and is developed in an open-source manner with great attention paid to high-quality software quality assurance practices. This overall approach could greatly benefit the fusion energy community, which requires rapid design iteration and improvement in order to facilitate the successful development of fusion as an alternative energy source to fossil fuels. In the first half of this talk, applications of MOOSE and MOOSE-based tools for advanced reactor designs will be showcased, as well as MOOSE ecosystem infrastructure (such as the NEAMS Virtual Test Bed) that enables and accelerates fission reactor design. In the second half, a discussion of how the MOOSE approach to modeling and simulation is currently being applied internationally in fusion energy research and development at the United Kingdom Atomic Energy Authority will be discussed, and ongoing/future domestic research efforts will be highlighted.
The coupling effect of stray AC (alternating current) + DC (direct current) and tensile stress-induced X70 pipeline corrosion is complicated with its mechanism uncovered. In this work, the corrosion behavior of X70 pipeline steel under three different corrosion (including an individual AC corrosion, mixed AC + 150 MPa tensile stress corrosion, coupling effect of AC + a 30 A/m{sup 2} DC interference and 150 MPa tensile stress corrosion) conditions in Dagang soil simulated solution was investigated by COMSOL Multiphysics software simulation, electrochemical tests, scanning electron microscope and other techniques. The results indicated that with the increase in AC density, the open-circuit potential of the X70 steel samples was negatively shifted and corrosion rate increased in three different corrosion conditions. The corrosion rates of X70 steel samples under individual AC interference were almost unchanged when 150.9 MPa tensile stress was applied. When applying 30 A/m{sup 2} DC interference, the corrosion rate of X70 steel increases significantly. Therefore, the factors affecting corrosion performance of X70 pipeline steel were DC > AC > stress. Under individual AC interference, uniform corrosion occurred on the X70 steel at an AC current density of less than 30 A/m{sup 2}, while the pitting corrosion occurred at an AC current density of greater than 100 A/m{sup 2}. The application of tensile stress and DC had a certain inhibitory effect on pitting corrosion.
Application packaging and distribution are the final steps for delivering software to end-users; both are frequently neglected when creating scientific software. Commercial businesses rely on electronic distribution systems that have rendered disk drives obsolete. Still, national laboratories continue to rely heavily on removable media to distribute and limit access to controlled applications. With increasing concerns of unauthorized copying of sensitive applications, a modern distribution system that utilizes cryptographically secure communication and authentication protocols has been developed. This new distribution system will secure the chain of custody for nuclear software while simultaneously simplifying access to these tools. This report summarizes four primary advancements made toward the secure distribution of Nuclear Energy Advanced Modeling and Simulation (NEAMS)-developed, Multiphysics Object Oriented Simulation Environment (MOOSE)-based applications: application installation, package distribution, automated package building, and distribution of documentation. NEAMS is currently developing more than ten separate applications based on the open-source MOOSE Framework. Distribution of these applications has primarily been accomplished by distributing source code, with end-users compiling the applications themselves. This work created a mechanism where MOOSE applications can be installed in a similar way to any other software. This allows both administrators and end-users simplified access to runnable executables. With this new installation capability, it was then possible to rethink distribution. A new, secure capability for delivering MOOSE-based applications over the internet has been created. This system requires unique cryptographic tokens for authentication, greatly securing the custody chain for software. Once granted access, installation of any NEAMS code can be accomplished with these terminal commands: "conda install ncrc" "ncrc install ncrc-bison." After these two commands (and authenticating) the BISON application will be securely down- loaded from Idaho National Laboratory (INL)’s servers, installed, and ready to use. To enable this new distribution capability to be successful, the open-source Continuous Integration, Verification, Enhancement, and Testing (CIVET) Continuous Integration (CI) capability was augmented to add Continuous Delivery (CD). CD enables the automated building and packaging of MOOSE-based applications as they are modified by development teams, ensuring that our customers can obtain up-to-date versions of the software at any time. The need for instruction on how to use these applications was addressed through modifications to the MOOSE documentation system. The MooseDocs capability, which enables robust documentation of MOOSE-based applications, has been extended to allow both for the installation of documentation and the packaging of documentation with installed applications. Together, these enhancements form the core of a new, secure distribution mechanism for nuclear simulation tools. In concert with the Nuclear Computational Resource Center (NCRC), NEAMS- developed applications will now be straightforward to obtain securely.
VERA-CS (the Virtual Environment for Reactor Applications, Core Simulator) is a high-fidelity multiphysics core simulator under development by the Consortium for Advanced Simulation of Light Water Reactors (CASL). VERA-CS includes three main components: MPACT for reactor physics and neutron transport, COBRA-TF for thermal-hydraulics, and ORIGEN for isotopic depletion. VERA-CS has been applied to an array of problems including core physics analysis, and the modeling for all of the fuel cycles of Watts Bar Nuclear Unit 1 (WBN1). MPACT is a 3D pin-resolved reactor transport code developed by Oak Ridge National Laboratory (ORNL) and the University of Michigan. COBRA-TF (CTF) is a three-field two-phase flow subchannel analysis code developed by ORNL and North Carolina State University. CTF solves the mass, momentum, and energy equations for the liquid, vapor, and droplet phases for each subchannel in the core and captures the axial flow within each channel as well as cross-flow between channels. Direct coupling between neutronics (MPACT) and thermal-hydraulics (CTF) in VERA-CS provides temperature and density feedback. CTF includes a number of models important for LWR safety and design analysis including flow regime dependent two-phase wall heat transfer, inter-phase heat transfer and drag, droplet breakup, and quench-front tracking. ORIGEN is developed by ORNL and released through the SCALE package. ORIGEN is directly integrated in MPACT through its API to calculate the changes in nuclides for every depletable cross-section region. Uncertainties exist in safety calculations with computer codes for nuclear reactor modeling and simulation, and they come from code limitations, scaling inaccuracies embedded in the experimental data, and uncertainties associated with the state of the reactor at the initiation of the transient. Consequently, for applications of high fidelity calculations using VERA-CS, proper uncertainty quantification (UQ) and sensitivity analysis (SA) of the computing models of reactor neutronics and thermal-hydraulics are important to comprehensively understand the results. New methods in UQ and SA have to be developed in multi-physics calculations using VERA-CS. Toward this end, a software package called VUSAT which stands for VERA-CS Uncertainty and Sensitivity Analysis Toolkit has been developed to facilitate the applications of VERA-CS within the context of best estimate plus uncertainty (BEPU) analysis framework. In this work, the VUSAT application to analyze departure from nucleate boiling (DNB) is presented.
The Multiphysics Object Oriented Simulation Environment (MOOSE) [1] is a massively parallel finite-element/volume package for multiphysics simulation in science and engineering. The package focuses on providing rapid-development capabilities for engineering applications by leveraging well-built features from libMesh [2] and the Portable Extensible Toolkit for Scientific Computation (PETSc) [3]. Fiscal year 2021 (FY-21) was the first year with funding dedicated to supporting MOOSE-derived applications relevant to the Nuclear Engineering Advanced Modeling and Simulation (NEAMS) program. In this report we outline the work done to support NEAMS applications such as BISON, Griffin, Pronghorn, and System Analysis Module (SAM).
Cardinal is an open-source application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics (CFD) to the Multiphysics Object-Oriented Simulation Environment (MOOSE), closing neutronics and thermal-fluid gaps in conducting high-resolution multiscale and multiphysics analyses of nuclear systems. We first provide a brief introduction to Cardinal's software design, data mapping, and coupling strategy to highlight our approach to overcoming common challenges in high-fidelity multiphysics simulations. Here we then present two Cardinal simulations for hexagonal pin bundles. The first is a validation of Cardinal's conjugate heat transfer coupling of NekRS's Reynolds-Averaged Navier Stokes model with MOOSE's heat conduction physics for a bare seven-pin Freon-12 bundle flow experiment. Predictions for pin surface temperatures under three different heating modes agree reasonably well with experimental data and similar CFD modeling from the literature. The second simulation is a multiphysics coupling of OpenMC, NekRS, and BISON for a reduced-scale, seven-pin wire-wrapped version of an Advanced Burner Reactor bundle. Wire wraps are approximated using a momentum source model, and coupled predictions are provided for velocity, temperature, and power distribution.