Engineering Papers⌕ Search

Engineering topics

Lindsay, Alexander D

Publications and source records attributed to Lindsay, Alexander D.

At least 19 records

New systems in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) serves as a common library of classes between applications developed for advanced reactor analysis, fusion device engineering, spent fuel cask analysis, geochemistry studies, among other fields. These applications drive the development of the framework to meet their needs. Systems in MOOSE group capabilities that share a common purpose and generally common code. They can be leveraged by all downstream applications, providing extensive code re-use and shared maintenance. They facilitate the discovery by new users of the classes meeting at least partially their needs, and offer the same opportunities for customization as other systems. The addition of a new system to MOOSE opens new ways of solving or discretizing nonlinear problems, of performing distributed postprocessing, and a plethora of other needs. While new systems can be introduced in downstream applications rather than at the framework level, the framework team monitors common needs across the community and often triggers their addition. Documentation, training material, development needs can be centralized, limiting duplicated work across the community. The last three years have seen a large expansion in the capabilities of MOOSE. The supporting role of the framework in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has created numerous feature requests to support neutronics, thermal hydraulics, computational fluid dynamics and thermo-mechanics simulations in the Griffin, SAM, Pronghorn and Bison applications respectively. Similarly, laboratory-directed research and development (LDRD) projects in additive manufacturing, high-Reynolds flow simulations, structure optimization also necessitate an expansion of the framework capabilities. This summary reports on the new systems created in MOOSE, their design, their capabilities and some of the relevant interfaces.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Developing Multiphysics, Integrated, High-Fidelity, Massively Parallel Computational Capabilities for Fusion Applications Using MOOSE

As the need for fusion as a clean, sustainable, and abundant energy source grows internationally, so does the need for multiphysics, computational tools to model, study, and predict the complex interactions between plasma, materials, and engineering processes. These tools have a crucial role to play in solving scientific and engineering challenges and accelerating fusion energy deployment. To address these needs, modeling capabilities should enable massively parallel, multiphysics, fully integrated high-fidelity simulations of fusion systems. Additional attributes, such as being open source and modular while maintaining high software quality assurance standards will maximize impact by ensuring accessibility for all and wide acceptance, rapid expansion and development, as well as reliability, efficiency, and robustness. In this paper, we describe how the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, which has a track record of success in the fission space thanks to the attributes listed above, can be leveraged in the fusion energy field. We highlight key successes of the MOOSE application in the fission space and describe how MOOSE has been and is being applied to fusion applications in the United States---e.g., Tritium Migration Analysis Program, version 8 (TMAP8), MOOSE Fusion Module, Fusion ENergy Integrated multiphys-X (FENIX)---and the United Kingdom---e.g., AURORA, Achlys, Apollo. These efforts aim to establish a suite of tools that can be further extended to accelerate fusion energy deployment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hybridized Discontinuous Galerkin Methods for Computational Fluid Dynamics

Hybridizable Discontinuous Galerkin (HDG) methods hold promise for any applications with significant advection character, including thermal hydraulics in light water reactors and advanced reactor concepts and fluid models of plasmas in magnetic confinement fusion. Its features include natural upwinding, local element conservation, and extensibility to arbitrarily high order accuracy. In the last fiscal year we have implemented HDG in the Multiphysics Object-Oriented Simulation Environment (MOOSE). We developed a first-of-its-kind automatic static condensation system in MOOSE’s underlying finite element library libMesh which can condense out arbitrarily many internal variables. Finally, we developed the first preconditioner for HDG discretizations of the Navier-Stokes equations which shows robust performance across a wide range of problem sizes and Reynolds numbers. This preconditioner yields solution times that are equivalent to the fastest developed for industry standard finite volume methods. Moreover, the arbitrarily high-order nature of HDG makes it a prime candidate for acceleration via graphical processing units (GPUs). We believe these developments will hold significant importance for future DOE Nuclear Energy (NE) and Fusion Energy Science (FES) programs.

97 MATHEMATICS AND COMPUTING↗

Development of Segregated Thermal-Hydraulics Solvers in MOOSE

The simulation of fluid flows is an essential part of the design and analysis of nuclear systems. Algorithms able to simulate flows at different fidelity levels are available in the Multiphysics Object-Oriented Simulation Environment (MOOSE) and MOOSE-based applications such as Pronghorn \cite{novak2018pronghorn}, Pronghorn-Subchannel, RELAP-7, and SAM. Currently, significant effort is being invested in the development of coarse-mesh Computational Fluid Dynamics (CFD) capabilities within MOOSE and Pronghorn for the simulation of Generation IV nuclear reactors. Traditionally, the solution algorithms in MOOSE have relied on Newton or quasi-Newton methods (such as the preconditioned Jacobian-free Newton-Krylov method) where residuals and Jacobians (or approximations thereof) are constructed. Both Newton and quasi-Newton methods require the solution of a linear system at each nonlinear Newton iteration with the Jacobian as the system matrix. The Jacobian contains blocks originating from all variables in the problem (i.e., for thermal-hydraulics at least pressure, velocities, and temperature). Due to the formulation of the problem in a general multiphysics setting on unstructured mesh, creating a good preconditioner for the linear system can be challenging, thus many fluid applications have utilized direct solver-based methods such as LU factorization. However, with increasing system size and complexity in multi-dimensional problems, the direct solution of linear systems becomes computationally expensive both in execution time and and memory. For this reason, recent effort has focused on adapting segregated solution algorithms for CFD problems in MOOSE. These algorithms use fixed-point iteration between segregated systems whose assembly and preconditioning are easier those of the monolithic system. Initial results show that the segregated solution algorithm outperforms the monolithic approach in terms of memory usage and for large 3D problems in terms of CPU time as well.

42 ENGINEERING↗

Development of Segregated Thermal-Hydraulics Solvers in MOOSE

The simulation of fluid flows is an essential part of the design and analysis of nuclear systems. Algorithms able to simulate flows at different fidelity levels are available in the Multiphysics Object-Oriented Simulation Environment (MOOSE) and MOOSE-based applications such as Pronghorn \cite{novak2018pronghorn}, Pronghorn-Subchannel, RELAP-7, and SAM. Currently, significant effort is being invested in the development of coarse-mesh Computational Fluid Dynamics (CFD) capabilities within MOOSE and Pronghorn for the simulation of Generation IV nuclear reactors. Traditionally, the solution algorithms in MOOSE have relied on Newton or quasi-Newton methods (such as the preconditioned Jacobian-free Newton-Krylov method) where residuals and Jacobians (or approximations thereof) are constructed. Both Newton and quasi-Newton methods require the solution of a linear system at each nonlinear Newton iteration with the Jacobian as the system matrix. The Jacobian contains blocks originating from all variables in the problem (i.e., for thermal-hydraulics at least pressure, velocities, and temperature). Due to the formulation of the problem in a general multiphysics setting on unstructured mesh, creating a good preconditioner for the linear system can be challenging, thus many fluid applications have utilized direct solver-based methods such as LU factorization. However, with increasing system size and complexity in multi-dimensional problems, the direct solution of linear systems becomes computationally expensive both in execution time and and memory. For this reason, recent effort has focused on adapting segregated solution algorithms for CFD problems in MOOSE. These algorithms use fixed-point iteration between segregated systems whose assembly and preconditioning are easier those of the monolithic system. Initial results show that the segregated solution algorithm outperforms the monolithic approach in terms of memory usage and for large 3D problems in terms of CPU time as well.

42 ENGINEERING↗

Modeling and Simulation of Advanced Manufacturing Techniques using MOOSE and MALAMUTE

Advanced manufacturing techniques offer increased geometry complexity, energy and material usage efficiency improvements, and an expanded palette of materials as compared to conventional manufacturing approaches. Advanced-manufacturing-produced parts can experience wide variations in the final microstructure, and these microstructure variations significantly impact the parts’ performance. In this chapter, we present recent code developments within Multiphysics Object-Oriented Simulation Environment (MOOSE) and in the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE). Here we demonstrate applying these modeling and simulation codes to two advanced manufacturing process types: advanced sintering techniques and laser-based additive manufacturing techniques. The multiphysics and multiscale capabilities of these codes enable prediction of the microstructure evolution resulting from variations in the Advanced manufacturing process parameters.

36 MATERIALS SCIENCE↗

The MOOSE fluid properties module

The fluid properties module in MOOSE~\cite{lindsay2022moose} serves a variety of fluid simulation applications based on MOOSE, including the MOOSE Navier Stokes module~\cite{moose_ns}, Pronghorn~\cite{pgh}, the MOOSE thermal hydraulics module, SAM~\cite{sam}, RELAP-7~\cite{relap7}, Sockeye~\cite{sockeye}, Pronghorn-subchannel~\cite{subchannel} and the MOOSE porous flow module~\cite{porous}. These applications are used to solve coarse mesh multi-dimensional thermal-hydraulics~\cite{pgh}, 1D systems analysis~\cite{sam,relap7} in nuclear reactor analysis, heat pipe modeling~\cite{sockeye} and porous flow simulations~\cite{porous} for underground gas storage and water seepage. The use of consistent fluid properties across fluid flow applications facilitates coupled flow simulations~\cite{anl_sam_pgh,osti_1889653}. Each application has historically driven the implementation of several fluid properties, which were later extended to be compatible with other applications. The unique diversity of applications of the module, due to its presence in MOOSE, has driven its expansion to new fluids, such as advanced nuclear reactor coolants and, more recently, arbitrary functions or tables-based property definitions, as detailed in section~\ref{content}, as well as numerous thermophysical properties and variable sets, as detailed in subsection~\ref{sec:prop}. The need for different discretizations of flow equations based on the compressibility of the fluid has motivated support for both primitive (pressure- and temperature-based) and conservative (internal energy- and specific volume-based) flow variables in the module. Thermodynamic relations are used to convert between these two formulations, as needed. The module serves a dual purpose of providing fluid properties for direct use in flow simulations and facilitating the implementation of user-specific fluid properties. Contributions of new properties for existing fluids or new fluids are strongly encouraged.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Latest developments in the MOOSE fluid properties module

The fluid properties module in MOOSE serves a variety of fluid simulation applications based on MOOSE, including the MOOSE Navier Stokes module~\cite{moose_ns}, Pronghorn~\cite{pgh}, SAM~\cite{sam}, the MOOSE thermal hydraulics module, RELAP-7~\cite{relap7} and subchannel~\cite{subchannel}. It is used for coarse mesh multi-dimensional thermal-hydraulics~\cite{pgh}, 1D systems analysis~\cite{sam,relap7} in nuclear reactor analysis, and porous flow simulations~\cite{porous} for underground gas storage and water seepage. The use of consistent fluid properties across fluid flow applications facilitates coupled simulations~\cite{anl_sam_pgh}. The module offers a consistent set of interfaces to implement to create a new fluid property. There are numerous fluid properties of interest in the entirety of all fields of fluid flow simulations, and this is exacerbated by the use of different variable sets depending on the compressibility of the fluid. For single-phase fluids, the following variable sets may be used to compute fluid properties: (pressure, temperature) and (specific volume, specific internal energy). Some properties may also be computed using the (pressure, density) or the (specific volume, specific enthalpy) variable sets. In order to reduce the challenge of adding a new fluid property, properties may be implemented partially.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

On Practical Aspects of Variational Consistency in Contact Dynamics

Usage of contact mechanics methodologies is a pervasive modeling requirement in dynamic simulations. While for some trivial problems, solutions taken from analytical geometry are available, use of a finite element framework is common to achieve formulation generality. This work explores two dynamic contact formulations: one based on the traditional node-to-segment (NTS) approach, and a variationally consistent segment-to-segment (STS) mortar formulation. The NTS formulation employed here enforces the constraints kinematically (i.e., the interpenetration is enforced to the solver tolerance), whereas the mortar approach uses Lagrange multipliers to enforce the contact constraints. Both approaches are implemented in the open-source finite element framework Multiphysics Object-Oriented Simulation Environment (MOOSE). The results highlight two relevant contact-interface-related dynamic phenomena in finite element simulations. First, stabilization of contact constraints is discussed, taking into account the evolution of the total energy in a benchmark problem. Second, the influence of finite element discretization on both of the aforementioned contact formulations is analyzed by exercising a large-deformation example with continuous relative sliding. Variationally consistent contact approaches such as the mortar formulation lead to improved energy preservation and avoid spurious excitation of the system's frequencies. This is especially relevant in settings where inertia and vibrations are of importance.

42 ENGINEERING↗

HI-STORM Overpack and MPC-32 Thermal-Hydraulic Model with MOOSE Framework

Nuclear power is a significant source of electricity in the United States, but the average age of nuclear power plants is around 40 years old. The safe management of the spent nuclear fuel (SNF) is a key aspect of the back-end of the nuclear fuel cycle. Spent fuel dry storage systems are becoming a popular and effective solution in this area, given the absence of a final disposal system. The spent fuel cask system (dry cask method) provides a feasible solution to maintain spent fuel for 60 years before final disposal. Dry cask storage has many characteristics that make it attractive. It fulfills the safety requirements of the Nuclear Regulatory Commission (NRC) while providing modularity and flexibility to contractors. The HI-STORM overpack and MPC-32 canister are the main parts of the HI-STORM 100 dry cask storage system. These components remove heat from the system using natural circulation, requiring no human intervention. This is the characteristic that provides passive heat removal and low maintenance features in dry cask storage systems. To develop a thermal model for a dry cask storage system, the physics behind the system should be defined clearly. There are two natural circulation loops in the system; circulation of helium cools down the nuclear assemblies in the MPC, while circulation of air cools down the walls of the MPC. This work aims to develop a thermal model of the MPC-32 canister and HI-STORM overpack using the Multiphysics Object-Oriented Simulation Environment (MOOSE). MOOSE is an open-source framework developed by Idaho National Laboratory (INL) for multiscale, multiphysics simulations. In this study, we will investigate and demonstrate the thermal-hydraulics modeling capabilities of the MOOSE framework, including natural circulation, heat transfer, and porous flow.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Electromagnetics and Advanced Manufacturing Simulation Development within the MOOSE Ecosystem

The MOOSE electromagnetics module (EMM) was developed to provide a general simulation tool for modeling a wide range of problems encountered in computational electromagnetics and enable MOOSE and MOOSE-based application coupling to Maxwell's equations to solve a wide class of new problems. Development of the electromagnetics module has also enabled development of new advanced manufacturing models in the MOOSE ecosystem, including that of the electric field-assisted sintering (EFAS) process. Recent results from this effort show good agreement with experimental results, and these two capabilities will be used in future projects to develop models for novel materials development and manufacturing for extreme environments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Tritium Transport Modeling at the Pore Scale in Ceramic Breeder Materials Using TMAP8

Fusion reactors depend on the blanket material to breed and release tritium at the same rate or faster than it is consumed by the fusion reaction. Cellular ceramic breeders (CCBs) are dense materials that can maintain a high tritium breeding ratio while promoting tritium release because of highly connected pores. Assessing the tritium breeding capabilities of these materials requires a combination of extensive experimental and modeling efforts. In this work, we develop and calibrate a multiphysics model of tritium transport. This novel model accounts for ceramic and pore diffusion, trapping and detrapping, and several surface reactions at the pore surface. We perform a sensitivity analysis and calibrate the model by comparing its predictions against experimental measurements of deuterium absorption. The calibrated model is then used to model tritium absorption in samples with different pore microstructures to investigate the effect of pore interconnectivity on tritium absorption. The model is part of the development of the multiscale, multiphysics framework for tritium transport [i.e., the Tritium Migration Analysis Program (TMAP8)], which is itself built on top of the finite-element multiphysics framework multiphysics object-oriented simulation environment (MOOSE). This study demonstrates some of TMAP8’s capabilities and is the first step toward assessing the tritium breeding capabilities of ceramic breeder material designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

TRITIUM TRANSPORT MODELING AT THE PORE SCALE IN CERAMIC BEEDER MATERIALS USING TMAP8

Fusion reactors depend on the blanket material to breed and release tritium at the same rate or faster than it is consumed by the fusion reaction. Cellular ceramic breeders are dense materials that can maintain a high tritium breeding ratio while promoting tritium release thanks to highly connected pores. Assessing the tritium breeding capabilities of these materials requires a combination of extensive experimental and modeling efforts. In this work, we develop and calibrate a multiphysics model of tritium transport. This novel model accounts for ceramic and pore diffusion, trapping and detrapping, and several surface reactions at the pore surface. We perform a sensitivity analysis and calibrate the model by comparing its predictions against experimental measurements of deuterium absorption. The calibrated model is then used to model tritium absorption in samples with different pore microstructures to investigate the effect of pore interconnectivity on tritium absorption. The model is part of the development of the multiscale, multiphysics framework for tritium transport (i.e., the Tritium Migration Analysis Program [TMAP8]), which is itself built on top of the finite element multiphysics framework MOOSE. This study demonstrates some of TMAP8’s capabilities, and is the first step toward assessing the tritium breeding capabilities of ceramic breeder material designs.

08 HYDROGEN↗

Calibration of a mesoscale tritium transport model for ceramic breeder materials in TMAP8 using experimental data

Due to the scarcity of the long-term external tritium supply, lithium-containing breeder materials are used in the blanket of fusion reactors to produce tritium faster than the deuterium-tritium fusion reaction consumes it. A cellular breeder is a promising material with higher lithium density and thermal conductivity than conventional ceramic pebbles and enhanced tritium extraction with interconnected pores. Although several experimental and modeling efforts have improved our understanding of tritium transport, the exact mechanisms governing tritium release from breeder materials are still unknown. As a result, current models cannot reliably assess the tritium breeding capabilities of cellular breeder materials. This work presents a tritium transport multiphysics model and calibrates it using experimental data of deuterium absorption and desorption from cellular breeder materials at different temperatures and pressures. This model accounts for ceramic and pore diffusion, trapping and detrapping, and several surface reactions at the pore surface. After calibrating the model and performing sensitivity analysis, we discuss pre-dominant mechanisms governing tritium release from cellular breeder materials. The model is part of the Tritium Migration Analysis Program [TMAP8], a multiscale, multiphysics framework for tritium transport based on the finite element multiphysics framework MOOSE. This study demonstrates some of TMAP8’s capabilities and provides insight into the mechanisms governing tritium transport in cellular ceramic breeder materials.

08 HYDROGEN↗

User-Oriented Improvements in the MOOSE framework in support of Multiphysics Simulation

The MOOSE Framework is a foundational capability used by NEAMS to create over 15 different simulation tools for advanced nuclear reactors. Due to this ubiquity, improvements to the framework in support of modeling and simulation goals are critical to the program. These improvements can take many forms including optimization, improved user experience, streamlined APIs, parallelism, and new capability. The work transcribed in the report was in direct support of the simulation tools and is already deployed or will be deployed in the coming months. The capabilities implemented were, in the same order as this report, chainable execution objects or executors, support for transfers between applications at the same level in a coupling scheme, support for boundary/subdomain restricted transfers, support for transfers between applications with different coordinate or unit systems, support for MOOSE applications in the NEAMS workbench, deployment of MOOSE application of the INL HPC OnDemand platform, addition of a triangular meshing library in libMesh and increased support of face variables.

97 MATHEMATICS AND COMPUTING↗

M3 milestone: Advanced contact 2021

This report describes advanced contact approaches and implementations in FY21 for use in emerging reactors with geometries that differ from LWR/EBR-II cylindrical fuel. Contact remains a challenging and crucially important component of our modeling toolkit. While significant improvements have been made through mortar and automatic differentiation, we explored advanced approaches (i.e., dual mortar) that yields better conditioned systems and will improve the robustness and performance of BISON contact problems. This effort is relevant to emerging reactor types, such as micro reactors and the corresponding complex and numerous contact surfaces. In this report, improvements of the mortar-based contact formulation are described. Benchmark problems and BISON assessment cases are included to demonstrate the improvements and advantages that are brought by the recent capability. Remaining issues are discussed and corresponding development plans are outlined.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Summary of Bison documentation and UX milestones - NEAMS FY20 Report

This summary report contains an overview of work performed under the work package entitled “MS- 20IN020107 - BISON advanced numerical model development and usability improvements - INL”, which is focused on the development and support of the fuel performance code BISON [1]. The second chapter lists FY20 milestones titles, completion schedule, and milestone level. Subsequent chapters summarize and demonstrate completion of milestones and activities. The last chapter outlines FY21 proposed future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗