Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphysics”

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 1,189 records · Page 66

SIERRA Code Coupling Module: Arpeggio User Manual (Version 4.46)

The SNL Sierra Mechanics code suite is designed to enable simulation of complex multiphysics scenarios. The code suite is composed of several specialized applications which can operate either in standalone mode or coupled with each other. Arpeggio is a supported utility that enables loose coupling of the various Sierra Mechanics applications by providing access to Framework services that facilitate the coupling. More importantly Arpeggio orchestrates the execution of applications that participate in the coupling. This document describes the various components of Arpeggio and their operability. The intent of the document is to provide a fast path for analysts interested in coupled applications via simple examples of its usage.

97 MATHEMATICS AND COMPUTING↗

Bison Robustness and Performance Improvements

In last year’s Bison robustness report, the authors demonstrated that improving the Jacobian matrix can drastically increase solution convergence rates. The development of an automatic differentiation (AD) capability in the MOOSE framework has created the possibility of having exact Jacobian matrices for all of our material models. A large focus in this year’s report is the application of AD to MOOSE physics modules, such as the tensor mechanics module, and the application to Bison fuel and cladding material models. A second main focus in this year’s report is the investigation of a set of algorithmic solver changes to improve the robustness and performance of Bison LWR simulations. Beyond providing increased solver robustness, the ability to solve strongly coupled multiphysics problems, and solving large deformation problems with ease, these automatically generated perfect Jacobians also drastically reduce development time for new models. In models converted to AD, Jacobian computation made up the majority of the code base. With this code removed, both readability and maintainability of the MOOSE/Bison code base are improved. The work documented here significantly improved the robustness of the Bison LWR assessment case suite. Through algorithmic changes in solver infrastructure, such as the application of predictors and improvements in timestepping, tangible improvements in performance were observed. These improvements are documented using a new standard set of metrics that are collected in all assessment cases, and a suite of statistical analysis tools developed for this report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SIERRA Code Coupling Module: Arpeggio User Manual (V. 4.40)

The SNL Sierra Mechanics code suite is designed to enable simulation of complex multiphysics scenarios. The code suite is composed of several specialized applications which can operate either in standalone mode or coupled with each other. Arpeggio is a supported utility that enables loose coupling of the various Sierra Mechanics applications by providing access to Framework services that facilitate the coupling. More importantly Arpeggio orchestrates the execution of applications that participate in the coupling. This document describes the various components of Arpeggio and their operability. The intent of the document is to provide a fast path for analysts interested in coupled applications via simple examples of its usage.

42 ENGINEERING↗

Modeling Continuous-Flow Reactor Improvements over Batch Reactions for Enzyme Catalyzed Microsphere Surface Reactions

Our research objective was to study enzymatic reactions occurring on microsphere surfaces. Specifically, our work involved enzymatic polymerization reactions (chain lengthening) taking place both in batch and in a continuous flow reactor. My research involved creating simulations using COMSOL Multiphysics 5.4a to further study variable interactions within the continuous flow reactor. Our goal was to use suitable software to create a simulation platform suite that could be used for advanced research studies of continuous flow enzymatic surfaces reactions. In addition, we studied various characteristics of our present enzymatic polymerization reactor assembly, modeled inlet species flow (Species M) and estimated outlet species production (Species P). Our simulations spanned three levels of detail: Single-Bead, Packed-Bead Cell and Full Reactor Column. Furthermore, the Packed-Bead Cells were studied using three possible regular arrangement of solid spheres: Face-Centered Cubic (FCC), Hexagonal Close Packed (HCP) and Cubic Close Packed (CCP).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Diagnostics and Models for Birthing Stress

Understanding the stress development in fluids as they transition to solids is not well-understood. Computational models are needed to represent "birthing stress" for multiphysics applications such as polymer encapsulation around sensitive electronics and additive manufacturing where these stresses can lead to defects such as cracking and voids. The local stress state is also critical to understand and predict the net shape of parts formed in the liquid phase. In this one-year exploratory LDRD, we have worked towards a novel experimental diagnostic to measure the fluid rheology, degree of solidification, and the solid stress development simultaneously. We debugged and made viable a "first-generation" Rheo-Raman system and used it to characterize two types of solidifying systems: paraffin wax, which crystalizes as it solidifies, and thermoset polymers, which form a network of covalent bonds. We used the paraffin wax as a model system to perform flow visualization studies and did some preliminary modeling of the experiment, to demonstrate the inadequacy of the current modeling approaches. This work will inform an advanced fluid constitutive equation that includes a yield stress, temperature dependence, and an evolving viscosity when we pursue the full proposal, which was funded for FY20.

36 MATERIALS SCIENCE↗

Progress Towards A New CTH Code Verification & Validation Test Suite

The CTH multiphysics hydrocode is used in a wide variety of important calculations. An essential part of ensuring hydrocode accuracy and credibility is thorough code verification and validation (V&V). In the past, CTH V&V work (particularly verification) has not been consistently well documented. In FY19, we have made substantial progress towards addressing this need. In this report, we present a new CTH V&V test suite composed of traditional hydrocode verification problems used by similar ASC codes as well as validation problems for some of the most frequently used materials models and capabilities in CTH. For the verification problems, we present not only results and computed errors, but also convergence rates. Validation problems include mesh refinement studies, providing evidence that results are converging.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SIERRA Code Coupling Module: Arpeggio User Manual (V.4.54)

The SNL Sierra Mechanics code suite is designed to enable simulation of complex multiphysics scenarios. The code suite is composed of several specialized applications which can operate either in standalone mode or coupled with each other. Arpeggio is a supported utility that enables loose coupling of the various Sierra Mechanics applications by providing access to Framework services that facilitate the coupling. More importantly Arpeggio orchestrates the execution of applications that participate in the coupling. This document describes the various components of Arpeggio and their operability. The intent of the document is to provide a fast path for analysts interested in coupled applications via simple examples of its usage.

42 ENGINEERING↗

Micro Reactor Simulation Using the PROTEUS Suite in FY19

Two micro reactor cores were simulated using the MOC solver of PROTEUS: MegaPower for a heat pipe micro reactor and Holos-Quad for a high-temperature gas-cooled micro reactor. Verification of PROTEUS standalone was performed by simulating 2D and 3D core benchmark problems developed based on MegaPower and Holos-Quad. Simulation results indicated that PROTEUS successfully modelled 2D and 3D cores of the MegaPower and Holos-Quad micro reactors and PROTEUS eigenvalues agreed well with the SERPENT Monte Carlo solutions. For a heat pipe cooled micro reactor, multiphysics simulations were performed using PROTEUS / FLUENT / ANLHTP. The coupling system of the three codes were developed using Python-based drivers, which exchange the data (power, temperature, and heat rate) required for the coupling and control the boundary conditions of FLUENT and ANLHTP. The coupled system was qualitatively verified using a 3D unit assembly problem composed of six fuel rods and seven heat pipes, which was developed based on the MegaPower compositions. The steady-state problem was solved using FLUENT / ANLHTP with given asymmetric fuel powers, demonstrating the reasonable convergence of temperatures of 7 heat pipes. Transient simulations were also performed by making one out of 7 heat pipes fail. For the transient, the coupled system of PROTEUS / FLUENT / ANLHTP demonstrated reasonable changes of total power and heat pipe temperatures with time, accounting for temperature feedback effects. In this simulation, no effect of structural change was accounted for, which will be performed using ANSYS-mechanical in future along with further verification with larger or whole-core 3D problems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DiaMonD: An Integrated Multifaceted Approach to Mathematics at the Interfaces of Data, Models, and Decisions

DiaMonD is a multi-institutional DOE MMICC effort that aims to address the challenges of end-to-end, data-to-decisions modeling and simulation for complex computational science and engineering problems in a unified, systematic, and integrated fashion. Institutions involved are Colorado State, Florida State, Los Alamos, MIT, Oak Ridge, Stanford, and UT-Austin. The goals of DiaMonD are (1) to develop advanced mathematical methods and analysis for multimodel, multiphysics, multiscale model problems driven by frontier DOE applications, including those in subsurface energy and environmental flows, materials for energy storage and conversion, and ice sheet dynamics; (2) to create theory and algorithms for integrated inversion, optimization, and uncertainty quantification for these complex problems; and (3) to disseminate the philosophy of an integrated end-to-end, data-to-decisions approach to modeling and simulation of complex problems to the broader applied math and computational science communities through workshops and other outreach.

42 ENGINEERING↗

Updates and Verifications of the PROTEUS Suite in FY19

In FY19 efforts were made on updating and verifying the PROTEUS suite in terms of depletion and fuel cycle capabilities, updated MSR capabilities, improving code usability, and supporting simulation of micro reactors. The fuel reloading and shuffling schemes as well as the built-in depletion solver were implemented to add the fuel cycle capability to PROTEUS-NODAL. PERSENT can run with PROTEUS-NODAL which now produce the CCCC interface files. The MSR simulation capability was further updated by making PROTEUS-NODAL coupled with SAM via the MOOSE wrapper. The high fidelity solvers of PROTEUS were updated to support the multiphysics simulation of SFRs and micro reactors. The coupled system of PROTEUS / FLUENT / ANLHTP was developed using Python-based external drivers, demonstrating the coupled simulation for small 3D test problems with heat pipe failure transients. In addition, the gamma transport capability was implemented to compute the detailed gamma distributions that determines the primary heat sources for the fuel assembly structure and non-fueled assemblies. Verification tests using fuel assembly and 2D core problems indicated very good agreement in heating distributions between PROTEUS and VARIANT / GAMSOR solutions. For user support, the NEAMS Workbench was updated to include PROTEUS-NODAL and further extended to handle PROTEUS-MOC for which the support of Workbench is limited, instead requiring the off-line generation of mesh files and corresponding cross section sets. The utility tools of PROTEUS were also updated to support users to model and simulate complicated and non-conventional geometry cores.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Use of Neutron Flux Calculated by Shift in a Grizzly Reactor Pressure Vessel Fracture Simulation

The reactor pressure vessel (RPV) plays the critical role of containing the reactor and coolant in a light water reactor nuclear power plant, and must be able to safely perform this function under a variety of normal and off-normal transient loading conditions. The RPV is subjected to harsh environmental conditions (elevated temperature and high neutron flux) whenever the plant is operating, which leads to embrittlement of the steel over time. RPVs contain populations of flaws introduced during the manufacturing process, and the primary safety concern is that during a transient event, a fracture could initiate at the site of one of these flaws and lead to rupture of the RPV. RPVs are very rugged structures, and the probability of this occurring is very low, but as the RPV's steel becomes embrittled over time, its resistance to fracture decreases and this becomes more likely. Assessing the likelihood of RPV fracture during a transient event involves performing a probabilistic fracture mechanics analysis, in which random sampling procedures are used to assess the probability of fracture of a population of flaws. Grizzly offers unique capabilities for modeling the RPV using 1D, 2D, or 3D representations, taking advantage of parallel computing for this type of analysis. One of the important inputs to Grizzly is the spatial distribution of the neutron fluence, which has a large impact on the material embrittlement, and should thus be characterized as accurately as possible. The VERA core simulator developed by the CASL program takes advantage of high performance computing resources to perform high-fidelity multiphysics simulations of reactor physics, including the effects of neutron transport, fuel performance, and thermal hydraulics. Recently a capability was developed to also evaluate neutron fluxes away from the reactor core through the Shift code. This allows for a detailed 3D map of the spatial variation of the fluence in components such as the RPV to be computed, which is ideal for use with Grizzly. This report documents a procedure developed to use the fluences computed in Shift in a Grizzly simulation, which allows for improved realisim in the prediction of the effects of environmental exposure on propensity for reactor pressure vessel fracture. This capability demonstrated on RPVs with multiple types of configurations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Induction Heating for Tritium Storage Beds

Metal hydride beds used for storing hydrogen isotopes accumulate helium due to beta-decay of tritium. Helium alters the hydrogen capacity and plateau pressures of these storage beds and ultimately results in the beds being disqualified for use in tritium facilities. Removal of helium and hydrogen isotopes for disposal or storage bed regeneration can be achieved with sufficient heating; however, existing Gen 1 and Gen 2 storage beds are not designed to reach helium removal temperatures and face disposal challenges. In this report, we describe experiments and Multiphysics models that were used to develop a method to achieve induction heating of existing tritium storage beds using induction coils placed on the outside of the beds. The heating method involves low-frequency (< 100 Hz) induction heating through flexible solenoid coils that wrap tightly around the outer jacket of a tritium storage bed. This method of heating offers a pragmatic path to heat existing storage beds for metal hydride regeneration, helium removal for disposal, and possibly routine gas cycling operations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DRACO and Nonequilibrium Statistical Mechanics of Aging [Slides]

Topics covered include: Types of problems which DRACO (Multiphysics simulation) was designed to address; Tools from nonequilibrium statistical mechanics; Rare events (Extreme value statistics); Transients (Diverging time scales); and, Thresholds (Absorbing phase transitions).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Transient and Steady-State Inverse Problems in Sierra/Aria

Inverse problems arise in a wide range of applications, whenever unknown model parameters cannot be measured directly. Instead, the unknown parameters are estimated using experimental data and forward simulations. Thermal inverse problems, such as material characterization problems, are often large-scale and transient. Therefore, they require intrusive adjoint-based gradient implementations in order to be solved efficiently. The capability to solve large-scale transient thermal inverse problems using an adjoint-based approach was recently implemented in SNL Sierra Mechanics, a massively parallel capable multiphysics code suite. This report outlines the theory, optimization formulation, and path taken to implement thermal inverse capabilities in Sierra within a unit test framework. The capability utilizes Sierra/Aria and Sierra/Fuego data structures, the Rapid Optimization Library, and an interface to the Sierra/InverseOpt library. The existing Sierra/Aria time integrator is leveraged to implement a time-dependent adjoint solver.

97 MATHEMATICS AND COMPUTING↗

ECP ST Project 2.3.1.06-STPM08-RAJA (Final Report)

Modern high-performance computing architectures are diverse, with hardware designs that feature multicore CPUs and CPUs plus accelerators, such as GPUs or FPGAs. Achieving good application performance often requires using a programming model and managing application data in ways that are best suited to a particular platform. For large codes that are used daily in production, and are under continual development, architecture-specific ports are untenable. Maintainability requires single-source code that is performance portable across a range of architectures and programming model back-ends. This report describes RAJA, Umpire, and CHAI, three software libraries that provide complementary capabilities to enable application portability for high performance computing (HPC) platforms. In particular, these libraries provide software abstractions that insulate application source code from hardware architecture and programming model-specific implementation details. These development of these libraries was motivated by the needs of multiphysics applications in the Weapons Simulation and Computing Program (WSC) at Lawrence Livermore National Laboratory (LLNL). Funding though the Exascale Computing Project (ECP) Software Technology (ST) project 2.3.1.06 RAJA has fostered outreach and collaboration with the ECP application and software technology development communities. In this report, we describe ECP adoption and integration of these tools to this point. The larger goal of the report is to describe to the ECP community what these projects offer, how they are used, and how to approach adoption and integration of them.

97 MATHEMATICS AND COMPUTING↗

Select Transient Metallic Fuel Bison Benchmarks

To explore the current capabilities of Bison with respect to transient metallic fuel modeling, four transient experiments are compared with their respective Bison model. The four experiments are the OPT-1 experiment of the Experimental Breeder Reactor II~(EBR-II) and the M5, M6, and M7 experiments of the Transient Reactor Test facility (TREAT). Each benchmark is a 2D-RZ geometry of a single pin employing coupled multiphysics. Bison models of the benchmark agree well and illustrate areas needing development. The areas of suggested development deal with better material models, better coolant channel module for transients, proper treatment of material melting, improved feedback with failure models, and an ability to set initial conditions from post irradiation examination (PIE) data. Two Fast Flux Test Facility (FFTF) pins are also modeled with M7 transient conditions to illustrate Bison being applied to investigate potential experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Advancing Scientific Productivity through Better Scientific Software: Developer Productivity and Software Sustainability Report

The Exascale Computing Project (ECP) provides a unique opportunity to advance computational science and engineering (CSE) through an accelerated growth phase in extreme-scale computing. Central to the project is the development of next-generation applications and software technologies that can exploit emerging architectures for optimal performance and provide high-fidelity, multiphysics, multiscale capabilities. However, disruptive changes in computer architectures and the complexities of tackling new frontiers in extreme-scale modeling, simulation, and analysis present daunting challenges to the productivity of software developers and the sustainability of software artifacts. Members of the CSE community - especially at extreme scales but more broadly at all scales of computing - face an urgent need to improve developer productivity, positively impacting product quality, development time, and staffing resources, and software sustainability, reducing the cost of maintaining, sustaining, and evolving software capabilities.

97 MATHEMATICS AND COMPUTING↗

SIERRA Code Coupling Module: Arpeggio User Manual - Version 4.56

The SNL Sierra Mechanics code suite is designed to enable simulation of complex multiphysics scenarios. The code suite is composed of several specialized applications which can operate either in standalone mode or coupled with each other. Arpeggio is a supported utility that enables loose coupling of the various Sierra Mechanics applications by providing access to Framework services that facilitate the coupling. More importantly Arpeggio orchestrates the execution of applications that participate in the coupling. This document describes the various components of Arpeggio and their operability. The intent of the document is to provide a fast path for analysts interested in coupled applications via simple examples of its usage.

97 MATHEMATICS AND COMPUTING↗