Engineering topics

Martin Recuero, Antonio

Publications and source records attributed to Martin Recuero, Antonio.

7 records

DOE OSTI2023

Interface Problem Formulation Improvements with Application to Nuclear Fuel Performance Analysis

The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation Program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. This report documents new developments and robustness improvements in mechanical and thermal (gap heat transfer) contact formulations. The improvements range from the migration of industrial level (“assessment”) nuclear fuel model setups to the usage of mortar formulations, the addition of frictional contact to one-dimensional layered representations of fuel and cladding components, and the addition of the Petrov-Galerkin approach to dual mortar, which improves performance on curved, relatively coarse meshes. In addition, the Lagrange-multiplier enforcement of mechanical mortar contact constraints has been extended to two additional types of enforcement: penalty and augmented Lagrange-Uzawa. We show that the latter approach yields the same interface results as dual mortar in the Multiphysics Object-Oriented Simulation Environment, with the advantage of not worsening the condition number of the system matrix—thereby enabling the use of some general implementations of iterative preconditioners, at the expense of additional system evaluations (i.e., augmentations).

42 ENGINEERING↗

DOE OSTI2023

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↗

DOE OSTI2023

Engineering-scale Modeling of High-Temperature Creep and Creep Crack Growth in Alloy 316H

This document demonstrates completion of the goals described in the technical narrative of the Department of Energy’s Industry Funding Opportunity Announcement (iFOA) project entitled “Modeling and Simulation Development Pathways to Accelerating KP-FHR Licensing,” which relates to the development and demonstration of capabilities for conducting engineering-scale simulations of Alloy 316H components under high-temperature conditions. This work encompassed two major aspects: integrating and testing constitutive models for the creep response of 316H at high temperatures, and developing tools for modeling creep crack growth (CCG) in 316H. Two classes of constitutive models were used in this effort: phenomenological models based on behavior observed at the engineering scale, and reduced-order models (ROMs) that represent the nonlinear response of mesoscale models that capture the sensitivity to material microstructure and processing. Likewise, the approaches employed for CCG modeling considered both simplified engineering approaches and detailed simulations of creep and damage ahead of the crack tip. These developments, which were performed by utilizing the Grizzly code as well as the open-source libraries it depends on, strengthen Grizzly’s ability to support licensing and safety analyses of high-temperature reactor components.

316H↗

DOE OSTI2022

Implement and Test 3D Mortar Contact in BISON

We leverage the extension of the generation of mortar segment meshes to three dimensions in MOOSE’s framework to extend thermomechanical modeling capabilities to problems with three dimensions. A modular approach to gap heat transfer physics using the mortar finite element method was created and documented, mechanical contact was extended to three dimensions—including frictional behavior, performance and ease of use were improved, and steps towards scalability of solid mechanics problems involving contact were taken. Many of these new developments are demonstrated in the simulation of 3D light-water reactor (LWR) problems, where the thermomechanical interface problem is solved using the mortar finite element method. Usage of the mortar framework has improved convergence in 2D problems and has enabled employing friction in 3D problems, of which we show results of a short, local stack of 3D pellets. Consequently, the benefits of mortar in terms of solution convergence and quality are extended to three dimensions. Section 2 discusses fundamental developments that enabled the simulation of practical mortar problems in three dimensions and other general improvements, including the reduction of the derivative container size, the modification of dual basis computations when edge dropping (lack of secondary element projection) takes place, the improvement of conditioning when employing the VCP in-edge dropping conditions, and code usability and quality improvements. These latter code enhancements include the migration of tests using “old” mortar contact constraints to using dual mortar with a semi-smooth Newton solution strategy and the reuse of lower dimensional domains for straightforwardly setting up a mortar thermomechanical LWR problem, i.e. the MOOSE action is employed for mechanical contact and the thermal LWR action is employed to capture the gas conductance, contact, and radiation components of gap heat transfer physics. Independently of the mortar LWR thermal action, we developed a modular approach to gap heat transfer that resides in MOOSE and can be leveraged, e.g., in metallic fuel problems. This approach, whose code design based on MOOSE’s user objects to model specific physics was proposed by the maintenance activity, is detailed in Section 3. Based on the dual mortar finite element method, the frictional contact constraints were extended to three dimensions. A block sheared in two directions in and out of contact with a rigid plane is employed in Section 4 to show the way the approach handles changes in frictional states (e.g. stick to slip) within a competitive number of Newton iterations. Equations and numerical results on the use of the VCP with Cartesian Lagrange multipliers, whose combination enables their direct condensation, are described in Section 5.3. Two-dimensional and three-dimensional BISON LWR simulations are discussed in Section 6. Particularly, a stack of five eccentric pellets with a surface defect is simulated and the effect of pellet-cladding friction is assessed. Finally, conclusions are outlined in Section 7.

42 ENGINEERING↗

DOE OSTI2021

Anisotropic inelastic materials milestone 2021

The focus on new fuel and cladding models are drivers for several new capabilities at the MOOSE and BISON levels. Since these materials are anisotropic, they require an extension of our inelastic mechanics capabilities, which currently heavily rely on radial return mapping and the assumption of isotropic stiffness. We plan to add additional code paths to support fully anisotropic materials, without degrading the performance for modeling isotropic materials.

42 ENGINEERING↗

DOE OSTI2021

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↗

DOE OSTI2020

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↗