Search NASA⌕ Search

SEARCH · Search NASA

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 199 records · Page 11

Multiscale Multiphysics Developments for Accident Tolerant Fuel Concepts

U 3 Si 2 and iron-chromium-aluminum (Fe-Cr-Al) alloys are two of many proposed accident-tolerant fuel concepts for the fuel and cladding, respectively. The behavior of these materials under normal operating and accident reactor conditions is not well known. As part of the Department of Energy’s Accident Tolerant Fuel High Impact Problem program significant work has been conducted to investigate the U 3 Si 2 and FeCrAl behavior under reactor conditions. This report presents the multiscale and multiphysics effort completed in fiscal year 2015. The report is split into four major categories including Density Functional Theory Developments, Molecular Dynamics Developments, Mesoscale Developments, and Engineering Scale Developments. The work shown here is a compilation of a collaborative effort between Idaho National Laboratory, Los Alamos National Laboratory, Argonne National Laboratory and Anatech Corp.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

COMSOL Multiphysics Model for HLW Canister Filling

The U.S. Department of Energy (DOE) is building a Tank Waste Treatment and Immobilization Plant (WTP) at the Hanford Site in Washington to remediate 55 million gallons of radioactive waste that is being temporarily stored in 177 underground tanks. Efforts are being made to increase the loading of Hanford tank wastes in glass while meeting melter lifetime expectancies and process, regulatory, and product quality requirements. Wastes containing high concentrations of Al 2 O 3 and Na 2 O can contribute to nepheline (generally NaAlSiO 4 ) crystallization, which can sharply reduce the chemical durability of high level waste (HLW) glass. Nepheline crystallization can occur during slow cooling of the glass within the stainless steel canister. The purpose of this work was to develop a model that can be used to predict temperatures of the glass in a WTP HLW canister during filling and cooling. The intent of the model is to support scoping work in the laboratory. It is not intended to provide precise predictions of temperature profiles, but rather to provide a simplified representation of glass cooling profiles within a full scale, WTP HLW canister under various glass pouring rates. These data will be used to support laboratory studies for an improved understanding of the mechanisms of nepheline crystallization. The model was created using COMSOL Multiphysics, a commercially available software. The model results were compared to available experimental data, TRR-PLT-080, and were found to yield sufficient results for the scoping nature of the study. The simulated temperatures were within 60 ºC for the centerline, 0.0762m (3 inch) from centerline, and 0.2286m (9 inch) from centerline thermocouples once the thermocouples were covered with glass. The temperature difference between the experimental and simulated values reduced to 40 ºC, 4 hours after the thermocouple was covered, and down to 20 ºC, 6 hours after the thermocouple was covered. This level of precision is considered acceptable for the scoping nature of the model and the subsequent laboratory glass studies Using the model, two additional glass pouring cycles were conducted. Representative thermocouple data were plotted to show the variations between the two cycles. This provides preliminary data that will be used in laboratory experiments to determine the potential for controlling nepheline crystallization in glass by varying the glass pouring conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Research and Development of Multiphysics Models in Support of the Conversion of the High Flux Isotope Reactor to Low Enriched Uranium Fuel

The findings presented in this report are results of a five year effort led by the RRD Division of the ORNL, which is focused on research and development toward the conversion of the High Flux Isotope Reactor (HFIR) fuel from high-enriched uranium (HEU) to low-enriched uranium (LEU). This report focuses on the tasks accomplished by the University of Tennessee Knoxville (UTK) team from the Department of Mechanical, Aerospace, and Biomedical Engineering (MABE) that provided expert support in multiphysics modeling of complex problems associated with the LEU conversion of the HFIR reactor. The COMSOL software was used as the main computational modeling tool, whereas Solidworks was also used in support of computer-aided-design (CAD) modeling of the proposed LEU fuel design. The UTK research has been governed by a statement of work (SOW), which was updated annually to clearly define the specific tasks reported herein. Ph.D. student Isaac T. Bodey has focused on heat transfer and fluid flow modeling issues and has been aided by his major professor Dr. Rao V. Arimilli. Ph.D. student Franklin G. Curtis has been focusing on modeling the fluid-structure interaction (FSI) phenomena caused by the mechanical forces acting on the fuel plates, which in turn affect the fluid flow in between the fuel plates, and ultimately the heat transfer, is also affected by the FSI changes. Franklin Curtis has been aided by his major professor Dr. Kivanc Ekici. M.Sc. student Adam R. Travis has focused two major areas of research: (1) on accurate CAD modeling of the proposed LEU plate design, and (2) reduction of the model complexity and dimensionality through interdimensional coupling of the fluid flow and heat transfer for the HFIR plate geometry. Adam Travis is also aided by his major professor, Dr. Kivanc Ekici. We must note that the UTK team, and particularly the graduate students, have been in very close collaboration with Dr. James D. Freels (ORNL technical monitor and mentor) and have benefited greatly from his leadership and expertise in COMSOL modeling of complex physical phenomena. Both UTK and ORNL teams have used COMSOL releases 3.4 through 5.0 inclusive (with particular emphasis on 3.5a, 4.3a, 4.3b, and 4.4) for most of the work described in this report, except where stated otherwise. Just as in the performance of the research, each of the respective sections has been originally authored by respective authors. Therefore, the reader will observe a contrast in writing style throughout this document.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Request for Information FEARCE: An extensible finite element multiphysics software suite for complex and reactive modeling

FEARCE is a modular multiphysics software package configured as a suite of fully-integrated tools that leverages custom templates and models alongside finite element analysis and computational fluid dynamics. FEARCE produces fast, high-fidelity solutions for any aspiring designer needing to model and predict turbulent, reactive, and otherwise complex fluid flows.

97 MATHEMATICS AND COMPUTING↗

SUNDIALS Multiphysics+MPIManyVector Performance Testing

In this report we document performance test results on a SUNDIALS-based multiphysics demonstration application. We aim to assess the large-scale parallel performance of new capabilities that have been added to the SUNDIALS suite of time integrators and nonlinear solvers in recent years under funding from both the Exascale Computing Project (ECP) and the Scientific Discovery through Advanced Scientific (SciDAC) program, specifically: 1. SUNDIALS’ new MPIManyVector module, that allows extreme flexibility in how a solution "vector" is staged on computational resources. 2. ARKode’s new multirate integration module, MRIStep, allowing high-order accurate calculations that subcycle "fast" processes within "slow" ones. 3. SUNDIALS’ new flexible linear solver interfaces, that allow streamlined specification of problem-specific linear solvers. 4. SUNDIALS’ new N_Vector additions of "fused" vector operations (to increase arithmetic intensity) and separation of reduction operations into "local" and "global" versions (to reduce latency by combining multiple reductions into a single MPI_Allreduce call). We anticipate that subsequent reports will extend this work to investigate a variety of other new features, including SUNDIALS’ generic SUNNonlinearSolver interface and accelerator-enabled N Vector modules, and upcoming MRIStep extensions to support custom "fast" integrators (that leverage problem structure) and IMEX integration of the "slow" time scale (to add diffusion).

97 MATHEMATICS AND COMPUTING↗

Benchmark problems for the Mesoscale Multiphysics Phase Field Simulator (MEMPHIS)

This report details the current benchmark results to verify, validate and demonstrate the capabilities of the in-house multi-physics phase-field modeling framework Mesoscale Multiphysics Phase Field Simulator (MEMPHIS) developed at the Center for Integrated Nanotechnologies (CINT). MEMPHIS is a general phase-field capability to model various nanoscience and materials science phenomena related to microstructure evolution. MEMPHIS has been benchmarked against a suite of reported classical phase-field benchmark problems to verify and validate the correctness, accuracy and precision of the models and numerical methods currently implemented into the code.

74 ATOMIC AND MOLECULAR PHYSICS↗

Towards an IMEX Monolithic ALE Method with Integrated UQ for Multiphysics Shock-hydro

A number of critical science, engineering and advanced technology applications require predictive analysis of complex shock-hydrodynamics of fluid/solid materials with possible electromagnetic interaction. The physical mechanisms include wave-phenomena, material-transport, diffusion, chemical reactions, and electromagnetics. The highly nonlinear multiple-time and length-scale response of these mechanisms include discontinuities formed from shocks, contact surfaces, and complex tabular equations-of-state (EOS). Current dominant computational solution strategies use ad-hoc combinations of operator- splitting, semi-implicit, and explicit time-integration methods and decoupled nonlinear-solvers. While these approaches have enabled progress in forward simulation, the inherited mathematical structure has not provided stability, accuracy and efficiency to resolve all the dynamical time-scales of interest, nor has it enabled integrated fast sensitivity analysis and uncertainty quantification (UQ). This draft report describes initial progress towards developing a new multiphysics shock-hydro capability that is intended to be more robust, mathematically well-structured and can readily be combined with advanced higher-order implicit/explicit (IMEX) time integration and efficient adjoint-enhanced uncertainty quantification (UQ) techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Coupling of MCNP and ABAQUS for Reactor Multiphysics (MARM) [Slides]

Summer 2021 Internship Objectives: Learn MCNP (Monte Carlo N-Particle) transport code - Unstructured Mesh & Burnup Features, and Implement Python Code for MARM framework - MARM = MCNP and Abaqus based Reactor Multiphysics. Slides provide background and details for achieving learning objectives, as well as ideas for future work in the area.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Assessment of the Griffin Reactor Multiphysics Application Using the Empire Micro Reactor Design Concept

In late 2019, INL and ANL agreed to jointly develop the reactor physics code named Griffin based on the integration of the two code suites, MAMMOTH/Rattlesnake (INL) and MC2 - 3/PROTEUS (ANL). Griffin is being developed based on the MOOSE framework and MOOSE quality assurance procedures. This decision was made to be able to allow DOE-NE to efficiently invest funding to this area and to provide effective and timely support for existing and potential users; the latter includes industry and government organizations who are developing various types of advanced reactors in the near and long term. Since MAMMOTH/Rattlesnake has been developed based on the MOOSE framework, the INL/ANL Griffin development team agreed to build Griffin beginning with a merger of MAMMOTH and Rattlesnake into a single code and moving forward by implementing capabilities from the PROTEUS suite into Griffin. Moving forward, both ANL and INL efforts are equally invested in the Griffin project, with management support, to provide an advanced reactor multiphysics tool to assist in reactor design, optimization, and safety analysis. Much work remains in moving Griffin forward to migrate PROTEUS capabilities and to optimize performance to meet user needs. The main objective of this work is to assess the current status of Griffin capabilities in terms of performance and accuracy, to determine priorities for PROTEUS migration, and to identify capabilities and features to improve for supporting the code integration effort. For this assessment, the Empire micro reactor problem that was developed in the ARPA-E MEITNER program was selected as an advance reactor concept of interest to the technical community. The Empire reactor problem was expanded from its original incomplete specification to be a small heat-pipe-cooled micro reactor core with ~113 cm radius and 70 cm in height, composed of 18 fuel assemblies, 12 control drums, and beryllium radial and axial reflectors. In the current model, using 5 cm axial reflectors specified in the original Empire assembly model, more than 10% of neutrons leak axially and through the empty center safety hole, as well as through heat pipe channels in fuel assembly elements that extend through the top reflector region. Several calculation models of the core were defined for systematic assessment, including 2-D and 3-D fuel assemblies and whole cores with cylindrical boundaries. Cross sections were generated using Serpent 2, and meshes were produced using the Argonne mesh tool or the INL neutronics meshing tools combined with CUBIT. Cross sections and meshes were converted to the ISOXML and Exodus formats, respectively, so that Griffin and PROTEUS could use consistent data for solving the reactor problems. With the prepared cross sections and meshes, PROTEUS was run first to ensure that all input data were correctly generated and input options in terms of angle, mesh, and energy group were accurately determined. Comparisons against Serpent 2 solutions were made in terms of eigenvalue and pin power. The same calculations and comparisons were then conducted using Griffin. For the fuel assembly and whole core problems, the PROTEUS eigenvalues agreed well with reference Serpent 2 solutions within 100 and 30 pcm, respectively, and pin power differences relative to Serpent 2 were overall less than 2.2% and RMS 0.8% for the whole core models. This indicated that all input data were properly prepared. Using the same data, Griffin was run selecting the SAAF-CFEM SN solver with Legendre-Gaussian quadrature and NDA and DSA for acceleration. It was found that the SAAF-CFEM solver of Griffin required finer meshes to achieve eigenvalue and pin power solutions in good agreement with Serpent 2, consequently requiring more memory requirement and longer computation time. On the other hand, the SPH-Diffusion 2-D core calculations performed using Griffin were able to recover the exact eigenvalue from the reference Serpent 2 solutions, resulting in a pin-power distribution with an RMS of 0.6% and maximum absolute difference of less than 1.4%. The runtimes for SPH-Diffusion for the 2-D core were less than 3 minutes on 40 cores. During this evolution of this evaluation, many updates were made in Griffin by the Griffin development team of INL (focusing on software updates) and ANL (reviewing and supporting software updates) to complete this assessment. Observations from the code assessment are presented in the conclusion section of this report, followed by a discussion of recommendations for future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Artificial Intelligence for Multiphysics Nuclear Design Optimization with Additive Manufacturing

The geometric flexibility of additively manufactured metals and ceramics generates a very large and open design space that requires advanced modeling and simulation tools for physics simulations and the rigorous definition of design problems. This effort deploys artificial intelligence (AI) and machine learning (ML) algorithms to understand the design space, evaluate potential designs, and more efficiently generate optimized results. The Transformational Challenge Reactor (TCR) program is leveraging advances in several scientific areas—including materials, manufacturing, sensors and control systems, data analytics, and high-fidelity modeling and simulation—to accelerate the design, manufacturing, qualification, and deployment of advanced nuclear energy systems. Through a manufacturing-informed design approach, the TCR program seeks to integrate digital data for rapid nuclear innovation; accelerate the adoption of advances in manufacturing, materials, and computational sciences for nuclear applications; and dramatically reduce deployment costs and timelines for new nuclear reactor technologies. This report documents efforts under the TCR program to leverage advanced modeling and simulation techniques driven by AI/ML algorithms on high-performance computing (HPC) systems to yield more optimized TCR core designs. A multiphysics ML surrogate model was developed to run on the HPC architectures. The surrogate model is trained on high-fidelity simulation data of coupled neutronics and thermofluidics and is used to quickly evaluate thousands of candidate core designs in parallel, which drives the evolution of the cooling channel shapes to minimize temperature peaking and material stress. Outcomes from these activities provide design information and feedback into the core design efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Final Review of FY23 ASC ATDM L2 Milestone, MRT #8541, ATDM Multiphysics Scaling on EAS-3

As displayed in the FY23 ASC Implementation Plan (IP) and the Milestone Reporting Tool (MRT), the milestone description and completion criteria state: Demonstrate the capability to execute mission-relevant calculations with a multiphysics code from the Ristra project at scale on the pre-El Capitan Early Access System 3 (EAS-3). These calculations must use the GPU hardware available on EAS-3, the AMD MI-250s, to achieve any reasonable performance.

97 MATHEMATICS AND COMPUTING↗

On-the-Fly Energy Condensation for Whole-Core Multiphysics Simulations

As the nuclear energy industry considers ways to achieve improved economics in the current fleet of light-water reactors (LWRs), one possible approach is to operate each cycle longer. This causes a greater portion of the fuel to be burned and reduces the frequency of outages, which ultimately reduces the cost to operate the reactor. However, this also leads to higher burnup fuels than have traditionally been allowed in these reactors. There are concerns about integrity of high-burnup (HBu) fuel, especially during accident conditions such as loss-of-coolant accidents (LOCAs), as shown by Capps et al. To investigate these concerns, advanced modeling and simulation capabilities are under development to determine the susceptibility of HBu fuel to fuel fragmentation, relocation, and dispersion (FFRD). Improvements have previously been made to fuel performance capabilities to model these phenomena more accurately; multiphysics simulations have also been conducted to determine the power and burnup histories of the HBu fuel, which are needed as inputs to the fuel performance calculations. Most recently, new statistical approaches have been developed to identify a subset of fuel rods that have greater FFRD susceptibility, reducing the total number of fuel performance simulations required.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Green's function methods for multiphysics simulations (Final Report)

Green's functions are important tools for analyzing mathematical properties of partial differential equations (PDEs), and for numerically solving PDEs, especially when equations for the same operator but with multiple right hand sides need to be solved simultaneously. This proposal aims at developing efficient and accurate numerical methods for computing Green's functions, which can be used to tackle a challenging question in multiphysics simulation of DOE-mission science problems: how to couple quantum physics with classical physics. The key mathematical difficulty is properly formulate a "boundary condition" for the region described by quantum physics, and conventional approaches often model such boundary conditions in an empirical way. The proposed Green's function methods use the Dirichlet-to-Neumann map to formulate a boundary condition that is non-empirical and can couple the quantum and classical regions in an in principle exact way. The key ingredient of the new methods is to construct the Dirichlet-to-Neumann map in an efficient, accurate and versatile manner. The new methods have provably low complexity and are ideally suited for massively parallel and emerging many-core computational systems.

97 MATHEMATICS AND COMPUTING↗

PhILMs: Collaboratory on Mathematics and Physics-Informed Learning Machines for Multiscale and Multiphysics Problems

The landscape of computational science and engineering is continually evolving, with the challenge of high-dimensional regression problems standing as a significant hurdle in numerous scientific endeavors. Addressing this challenge, our research, funded by this award, has led to the development of an innovative computational framework known as Probabilistic Partition of Unity Networks (PPOU-Nets). This initiative represents a collaborative effort to harness the potential of mathematics and physics-informed machine learning in tackling multiscale and multiphysics problems prevalent in high-dimensional spaces. Through this work, we have proposed a novel methodology that seamlessly integrates adaptive dimensionality reduction and a mixture of experts model, thereby facilitating a more efficient and accurate approximation of complex functions. This research effort has not only advanced the state of computational science but also opened new avenues for exploration in quantum computing and beyond. This report outlines the motivation, methodology, key findings, and implications of our work, underscoring our contributions to the broader scientific community and the potential pathways for future research.

97 MATHEMATICS AND COMPUTING↗

In-House Developed Multiphysics Simulation for the Performance of Solid Oxide Cells (SOCs)

Reversible solid oxide cells(rSOCs) are an enabling energy storage/production technology for a dynamic grid environment. Reversible operation requires strong working knowledge of fuel cell (SOFC) and electrolyzer (SOEC). Performance degradation of SOECs has been observed; however, the details of physical processes related to the performance degradation remain unknown. Multiphysics simulations were performed to investigate the performance degradation of solid oxide electrolysis cells under various working conditions.

Yang, Tao↗

Multiscale thermal properties prediction in the Multiphysics Object Oriented Simulation Environment (MOOSE) via a general Boltzmann solver [Poster]

The project objectives were: (1) Engineering-scale fuel performance modeling relies on accurate thermal properties; (2) Thermal properties are inherently multi-scale, arising from atomistic processes and interactions with a material’s microstructure; (3) Heat transport in solids via conduction occurs through transport and scattering of electrons and phonons; (4) Use the Boltzmann transport equation (BTE) to predict the macroscopic behavior of a materials system in terms of the microscopic dynamics of its heat carriers; and (5) This project establishes a new MOOSE (Multiphysics Object Oriented Simulation Environment) module, Boltzmann, dedicated to phonon and thermal electron transport.

36 MATERIALS SCIENCE↗

Ristra Project FY23 L2 Milestone Report, Rev.1: MRT #8541: Multiphysics Scaling on EAS-3

The findings of this report were used to close out the ATDM milestone MRT# 8541, which was designed to demonstrate readiness of ATDM multiphysics codes for mission-relevant work on ATS-4, El Capitan. To this end, the closure criteria were to run a 3D shaped charge problem at scale up to 50% of the El Capitan early-access system, RZVernal (AMD Trento CPUs and AMD MI-250X GPUs), demonstrate scalability, and document challenges with the software stack and environment. LANL’s approach to this milestone was to test our modular software capability by developing an entirely new code, Moya, built upon our FleCSI framework. The physics capability and the GPU infrastructure needed for the shaped charge problem on GPUs was added to Moya, and the required calculations were performed at scale. Moya showed good scaling without any fine-tuning of GPU kernels; there is still significant room for performance enhancements, especially for the Legion backend. Tied up in this L2 milestone was a closeout of KPP-3s for the ECP ST Projects at LANL; this material will be covered in a separate document.

97 MATHEMATICS AND COMPUTING↗

Multiphysics Simulations of Thermal Shock Testing of Nanofibrous High Power Targets

Increase of primary beam power for neutrino beam-lines leads to a reduced lifespan for production targets. The field of High Power Targetry (HPT) is generating new concepts to meet the need for robust targets. One idea being investigated by the HPT Research and Development Group at Fermilab is an electrospun nanofiber target. As part of their evaluation, samples with different densities were sent to the HiRadMat facility at CERN for thermal shock tests. The samples with the higher density, irradiated under a high intensity beam pulse, exhibit major damages at the impact site whereas those with the lower density show no apparent damages. The exact cause of this failure was unclear at the time. In this paper, we present the results of multiphysics simulations of the thermal shock experienced by the nanofiber targets that suggest the failure originates from the reduced permeability of the high density sample to airflow. The air present in the porous target expands due to heating from the beam, but is unable to flow freely in the high density sample, resulting in a larger back pressure that blows apart the mat. We close with a discussion on how to further validate this hypothesis.

43 PARTICLE ACCELERATORS↗