Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “code verification”

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

Demonstration of Model-Based Design for Digital Controller Using Formal Methods

This report describes work originally performed in FY19 that assembled a workflow enabling formal verification of high-consequence digital controllers. The approach builds on an engineering analysis strategy using multiple abstraction levels (Model-Based Design) and performs exhaustive formal analysis of appropriate levels – here, state machines and C code – to assure always/never properties of digital logic that cannot be verified by testing alone. The operation of the workflow is illustrated using example models and code, including expected failures of verification when properties are violated.

97 MATHEMATICS AND COMPUTING↗

Viscoelastic Modeling of Polymers in ALEGRA with the GAP Model

The Glassy Amorphous Polymer (GAP) model is a viscoelastic/plastic model developed at Los Alamos National Laboratory to accurately model a variety of polymers across a wide range of conditions and loading rates, including shock loading. In the present report we introduce and assess this model, newly implemented in the ALEGRA shock and multiphysics code, using a series of verification and application-related validation problems. We describe the mathematical and theoretical formulation of the model, as well as its implementation in ALEGRA, in detail. We provide verification results that assess the model implementation against published computational results, as well as validation results which we compare to existing experimental results when possible. These comparisons instill confidence in the implementation and indicate that the addition of the GAP model to the ALEGRA code provides users with a high-fidelity polymer modeling capability that is capable of recreating complex polymer phenomena.

36 MATERIALS SCIENCE↗

A process to verify numerical models for seismic fluid-structure interaction in advanced reactor vessels

Seismic design and qualification of a liquid-filled advanced nuclear reactor will have to account for fluid-structure interaction (FSI). Interaction between the tank, internal components, and contained liquid will rely on analysis of numerical models that must be verified and validated. Here this study demonstrates a verification process for models of a base-supported cylindrical tank by comparing numerical predictions and analytical solutions. The numerical models are consistent with the assumptions made to derive analytical solutions, namely, either a rigid or a linear elastic tank, ideal fluid, and small-amplitude, unidirectional, horizontal inputs. One software platform is used to illustrate the process. Seismic FSI analysis is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Reported responses are those used for design, including hydrodynamic pressures on the tank wall, shear forces and moments at the tank base, and wave heights of the contained liquid. The accuracy of the numerical results is discussed. The numerical models are verified for calculating the pressures on the tank wall and reactions at its base. Accurate simulation of wave action is challenging for both solvers. Recommendations for modeling, code development, and steps for verification are provided. Although focused on reactor vessels and one software platform, the verification process described herein is broadly applicable to liquid-filled vessels and other finite element codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Towards Ultra-high-resolution E3SM Land Modeling on Exascale Computers

Here we present an ultra-high-resolution E3SM land model (uELM) for high-fidelity land simulations targeting new Exascale computers. After considering modeling infrastructure compatibility and ELM software features, we designed a parallel model for the uELM development targeting hybrid architectures of new US Exascale computers. We also described a function unit test framework to expedite the piece-wise code porting (with compiler directives), verification, and global variable management. Furthermore, in this study, we report an early uELM model development using OpenACC within a function unit test framework on a pre-Exascale computer, demonstrate the performance of a uLEM submodel with a 3.0-time speedup, and summarize the code porting experience regarding global variable handling, deepcopy, memory reduction, and parallel loop reconstruction.

97 MATHEMATICS AND COMPUTING↗

Black Box Equations of State: Creating Semi-analytic Solutions to the Noh Problem and Verifying Equation of State Interfaces

The objective of this report is threefold. First, it details a method for deriving a semi-analytic solution to the Noh Problem when using a “black-box” equation of state. Such capability allows us to perform verification on complicated, more realistic equations of state. Examples include Steinberg equations of state for materials and tabulated equations of state. The second objective is to apply the methodology to verify the singularity-eos equation of state library. We do so by solving the Rankine-Hugoinot jump conditions for the Noh Problem, ensuring singularity derives the correct solution and comparing the error to an exact implementation of the equation of state. The third objective is to perform verification of the xRAGE Eulerian hydrodynamics code when interfaced with singularity. We provide the theory, analysis, documentation for a python implementation of the proposed solver, and verification results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Review of Experimental Data for Validating Computer Codes Used in Shielding Calculations for Spent Fuel Storage and Transportation Systems

This report presents a review of available radiochemical assay data and shielding benchmarks applicable to spent nuclear fuel (SNF) shielding calculations. The relevant information reviewed herein includes the Spent Fuel Composition (SFCOMPO) database, the Shielding Integral Benchmark Archive and Database (SINBAD), the International Handbook of Evaluated Criticality Safety Benchmark Experiments, and published measurements of external dose rates of casks loaded with SNF. The relevant experimental data identified in this report may be used to support verification and validation of computer codes used in SNF cask/transport shielding applications, as well as development of calculation uncertainties. It should be noted that a relatively small subset of the identified experimental data (e.g., criticality alarm experiments) is available in a standard format established by the international community participating in experimental isotopic and shielding data evaluations. An effort of the SFCOMPO Technical Review Group (TRG) is underway to publish first isotopic evaluations of individual assay data using a standard data evaluation format. The SINBAD TRG has recently initiated benchmark evaluations and modernization of the database. Therefore, more relevant information is expected in the future that will enable users to select quality experimental data in depletion code and shielding code validations for SNF applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling of Macroscopic Dynamics in Three-Dimensional Plasma Configurations (Final Technical Report)

A numerical simulation-based study to investigate the stability properties of inherently three-dimensional plasma magnetic configurations is presented. These configurations, called stellarators and heliotrons, use external helical coils to provide the twisting magnetic field that is required to confine ionized plasma particles. In contrast to the tokamak configuration, which nominally has the property of axisymmetry, electrical current does not need to be induced in stellarator and heliotron plasma. Experiments have shown that these three-dimensional configurations are remarkably robust to macroscopic instability, even when driven to relatively high levels of plasma pressure. Our understanding of how the plasma is able to withstand macroscopic instability will be enhanced by the work proposed here. The effort will generalize existing numerical modeling capability, and it will apply the new capability to model macroscopic dynamics in these three-dimensional configurations. The development work has been applied to the NIMROD simulation code (https://nimrodteam.org) that was originally based on axisymmetric problem domains. It has been generalized to accommodate the three-dimensional shaping of the plasma volume in stellarators and heliotrons. Its use of equilibrium data has also been generalized, and the implementation allows different numerical representation of the magnetic field. Analysis has been performed to inform the selection of numerical methods, and a method to accelerate the iterative solution of the necessary algebraic systems has been developed and implemented. Verification of the new NIMSTELL code on benchmark problems has also been an important part of the effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Verification and Validation Approach for COMSOL Multiphysics to Support High Flux Isotope Reactor (HFIR)

Over the last several decades, many reactors have successfully been converted from high enriched uranium (HEU) to low enriched uranium (LEU) fuels in United States in support of its global non-proliferation objectives. Of the reactors slated for conversion, five high-performance research reactors (HPRRs) remain. The high-flux isotope reactor (HFIR) at Oak Ridge National Laboratory (ORNL) is one of the five. Conversion of HFIR requires the qualification of a new fuel. To aid in the conversion process as well as improve safety margins, COMSOL Multiphysics was chosen to support and eventually supplement the steady-state heat-transfer code (SSHTC). In order to use COMSOL as a supplementary tool, i.e. one used for safety basis calculations, the code needs to undergo a verification and validation process. Verification and validation are important processes for the use of any software; however, physics-based solvers have a unique challenge, in that the validation is much more rigorous and its domain is only applicable to the general scope of the problem. This document details the verification and validation of COMOSL for analysis of HFIR. This report is intended for two different users: (1) current and future users of COMSOL for HFIR and (2) for individuals interested in assessing the scope of the validation of COMSOL for HFIR. This report is divided into six sections. The first of which details the process of verification and validation. Section one provides an overview of key concepts from ASME V&V 10, 20, and 40 and how these concepts are used to provide a validation of any multiphysics software. Section two covers the approach needed in order to validate multiphysics software for a given HFIR analysis. Parts three through five split up the key concepts highlighted in section two. The final section provides concluding remarks which demonstrate the scope and limitations of the validation study and how these may be expanded and improved.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

BlueCRAB Domain Overlapping Coupling: Theory & Verification

Coupling low- and high-fidelity codes is a useful way to model complex engineering systems. Thanks to the high-fidelity code, complex phenomena can be resolved in areas of the system where this is required, yet the efficiency of the lower-fidelity code is still retained in modeling the rest of the system. This document details the theory and implementation of coupling two different thermal-hydraulic codes: the system thermal-hydraulics (STH) code System Analysis Module (SAM) and the coarse-mesh computational fluid dynamics (CFD) code Pronghorn. Both applications are included in the Comprehensive Reactor Analysis Bundle (BlueCRAB) code suite and are based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. The domain overlapping (DO) coupling approach was adopted, as it offers proven advantages over more conventional domain decomposition methods. In the DO coupling, SAM provides Pronghorn with boundary conditions that depend on the system-level simulation of the entire plant. In return, the overlapping coupled SAM components–termed “surrogate components”–receive friction factors and source terms computed online based on the Pronghorn simulation. The framework is designed to be generic and enable coupling regardless of geometry and the number of inlet/outlet boundaries in the DO coupled domain. The developed method leads to consistent pressure drops, enthalpies, and scalar concentrations when comparing the coupled SAM and Pronghorn simulations. This document presents the DO coupling approach, along with two verification and two demonstration cases. The proposed problems explore different physical aspects relevant to nuclear reactor analysis, including buoyancy-driven flows, complex flow patterns, and multiple inlets and outlets. Periodically, new versions of this “BlueCRAB Domain Overlapping Coupling: Theory & Verification” report will be issued to reflect future developments and verification tests.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

An Educational Guide for 2D Stellar Structure Calculations of Rapidly Rotating Stars using the ESTER code

The Evolution STEllaire en Rotation (ESTER) code is the first 2D stellar structure code to be made open-source and freely available to the astronomy and astrophysics community. An important and novel advancement of this code is that it can reproduce the distorted shape and observable signatures (e.g., gravity darkening) of rapidly rotating stars. ESTER also calculates the steady-state large-scale flows within the star, namely their differential rotation and associated meridional circulation. In this report, we explore and document the physics implemented within version 1.1.0rc2 of the ESTER code, in a way that complements published descriptions. We illustrate this physics by plotting how stellar structure parameters vary through stellar interiors at a range of latitudes and at different angular velocities. We investigate how the thin convective envelopes of intermediate mass stars vary with latitude when rapidly rotating, becoming deeper and thicker near the equator. Simple comparisons of ESTER model predictions (e.g., central temperature and density, luminosity) with the output from the Modules for Experiments in Stellar Astrophysics (MESA) code [Paxton et al., 2010] shows generally good agreement. Additional comparisons provide important benchmarking and verification for ESTER as a comparatively young code. Finally, we provide a guide for installing and running the code on our local university cluster, aimed at helping students to begin work.

79 ASTRONOMY AND ASTROPHYSICS↗

Verification of the PERSENT Software

Ongoing commercial design activities require a thorough verification of the Argonne Reactor Computation codes to be performed. DIF3D is central to this system and substantial work has been done to verify its accuracy on several identified commercial needs. This manuscript details the verification work done on PERSENT which relies upon the DIF3D for its forward and adjoint flux solution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Effect of net direct current on the properties of radio frequency sheaths: simulation and cross-code comparison

In order to understand, predict and control ion cyclotron range of frequency (ICRF) interactions with tokamak scrape-off layer plasmas, computational tools which can model radio frequency (RF) sheaths are needed. In particular, models for the effective surface impedance and DC rectified sheath potentials may be coupled with full wave RF simulation codes to predict self-consistent wave fields near surfaces and the resulting power dissipation and plasma-material interactions from ion sputtering. In this study, previous work assuming zero net DC current flow through the sheath is generalized to allow the surface to collect net positive or negative current, as is often observed in experiments. The waveforms, DC potential and RF admittance are investigated by means of analytical theory, nonlinear fluid and particle-in-cell (PIC) codes. Cross-code comparisons provide detailed model verification and elucidate the roles of ion and electron kinetics. When the sheath draws negative (positive) DC current, the voltage rectification is reduced (increased) compared with the zero-current case, and the magnitude of both the real and imaginary parts of the admittance are increased (reduced). In conclusion, a previous four-input parametrization of the sheath rectification and admittance properties is generalized to include a fifth parameter describing the DC sheath current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CFDverify

Estimating the discretization error of computational fluid dynamics (CFD) or other scientific codes as part of solution verification is often a non-trivial part of the analysis process. Methods can be complicated and may include assumptions/qualifications that need to be checked during analysis. CFD analysts, therefore, are likely to make errors when trying to conduct this necessary analysis on their own in not knowing about the best method for their problem, not correctly implementing a method, or in not have diagnostic tools to determine if the method was correctly applied.

Weinmeister, Justin [Oak Ridge National Laboratory↗

FY2021 Progress Report on BISON Metallic Fuel Model Development and V&V Using EBR-II Legacy Data

In this report, the activities and achievements made by Argonne National Laboratory for the Nuclear Energy Advanced Modeling and Simulation (NEAMS) BISON code metallic fuel validation and verification project in FY2021 are summarized. The cladding degradation model based on the FCCI/CCCI wastage calculations has been developed and implemented into BISON. A comprehensive evaluation of the cladding degradation model was performed based on FIPD data of the IFR experiment X447. BISON objects were also developed to enable direct use of time-varying cladding outer surface temperature profile as temperature boundary conditions, which proved to provide more accurate temperature predictions for the metallic fuel pins irradiated in EBR-II. Additionally, a new BISON object was implemented to enable direct comparison between BISON predicted data and FIPD-based post-irradiation examination (PIE) results, which would significantly facilitate BISON metallic fuel verification and validation (V&V) activities. These new BISON-FIPD integration features were used in the establishment of a low-burnup fuel swelling evaluation framework as demonstration. The framework was successfully used to evaluate current fuel swelling models based on the IFR experiment X423.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Core Design of the Holos-Quad Microreactor

The Holos-Quad micro-reactor concept, developed by HolosGen LLC, is equipped with a 22 MWt (Mega-Watt thermal) core and an integral power conversion system converting the core thermal energy into approximately 10 MWe (Mega-Watt electric). This design can be configured to support a wide range of applications. It is a very innovative high-temperature gas-cooled reactor concept using TRI-structural ISOtropic particle fuel (TRISO) distributed in graphite hexagonal blocks, cooled with helium in a direct Brayton cycle independently executed by four Subcritical Power Modules (SPMs) fitted into a hardened 40-foot container whose dimensions are in compliance with ISO shipping containers requirements. In FY2019 HolosGen LLC was awarded by the Department of Energy Advanced Research Project Agency-Energy (DOE ARPA-E) under the MEITNER funding program. As part of the MEITNER award, the Argonne National Laboratory (ANL) contributed expertise through two specialized teams: The “Design Team” and the “Resource Team”. The Design Team was dedicated to validate feasibility of the Holos-Quad core and to optimize its core design through neutronics analyses. The Resource Team was dedicated to feasibility verification via high-fidelity codes of Holos-Quad thermal-hydraulic, heat transfer, shielding, and structural aspects. This report summarizes the activities conducted by ANL Design Team. A rigorous design approach based on multi-criteria optimization and code-to-code comparison involving stochastic and high-fidelity deterministic solutions was developed and employed at several evolutionary stages of the Holos-Quad design. Several generations of the Holos-Quad core were designed within this project before converging to the current full-scale Gen 2+ design that is detailed in this report. Figure EA-1 illustrates a cross-sectional view of Gen 2+ Holos-Quad core configuration, and Figure EA-2 provides a simplified perspective view of 1-of-4 SPMs. The Holos-Quad uses four thermal-hydraulically independent SPMs locked into stationary positions during power operation, surrounded by BeO reflector and structural component fully comprised within the dimensional constraints represented by traditional ISO containers. One of the benefits of this approach is to enable transportation of each SPM promptly after irradiation in shielded containers. The core is designed to operate for approximately 8 full-power years while the reactivity controls and power conversion system enable load-following operations. The reactivity controls are represented by independent, diversified, and redundant reactivity control systems based on control drums and redundant sets of shutdown rods. The high-fidelity simulation tools were used to assess detailed power and flux distributions of the three-dimensional full-core or quarter-core of the Gen 2+ configuration. Single-physics and multi-physics simulations of the neutronics code PROTEUS and the thermal-hydraulic code System Analysis Module (SAM) were performed to analyze the Holos design configurations with detailed high-fidelity solutions. The design work performed confirmed feasibility of the Holos-Quad concept, provided realistic design description for detailed design of the operational system, and identified several core design improvements to be further considered for future reactor development activities.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Update of CEA DES Criticality-Safety Activities and Perspectives [Slides]

Neutronics staff at CEA DES are in charge of a whole set of modeling and simulation tools for neutronics: nuclear data evaluation and processing, transport codes development, calculation sequences edition, verification, validation and uncertainty analysis. Advanced user groups are in charge of specific nuclear analysis in reactor physics, fuel cycle, criticality-safety, radiation shielding and nuclear instrumentation for CEA and its partners, at all stages of nuclear facilities life cycle. CEA is willing to participate/collaborate on the following issues: detector design, experiments, analysis and evaluation, DH validation by combining TAGS measurements with Fission Yield variance-covariance matrices.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

FY22 Progress Report on BISON Metallic Fuel Model Development and V&V Using EBR-II Legacy Data

In this report, the activities and achievements made by Argonne National Laboratory for the Nuclear Energy Advanced Modeling and Simulation (NEAMS) BISON code metallic fuel validation and verification project in FY2022 are summarized. The FIPD-BISON integration powered metallic fuel low-burnup swelling framework has been enhanced to cover both radial and axial swelling analyses. Axial-dependent as-fabricated fuel radius profiles have been used to improve the accuracy in evaluating radial swelling strain. The simulations of the IFR experiment X447, which was focused on FCCI/CCCI induced cladding degradation and failure, have been converted into a BISON metallic fuel assessment case with detailed documentation, which is also the first FIPD-BISON integration powered assessment case in BISON. A FIPD-BISON integration repository has been established to be an optional submodule of BISON to support FIPD-BISON integration powered assessment cases. Preliminary trials to simulate the in-cell out-of-pile transient experiment have been made after implementing some essential models such as cladding plasticity models and liquid phase penetration models. The initial results are promising and help identify gaps that need to be filled in FY2023 and beyond, to eventually achieve OPTD-BISON integration powered assessment cases.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Enhance BISON Metal Fuel Transient Models and Enable Assessment Based on Out-of-Pile Transient Experiments

In this report, the activities and achievements made by Argonne National Laboratory for the Nuclear Energy Advanced Modeling and Simulation (NEAMS) BISON code metallic fuel validation and verification (V&V) project in FY2023, with a focus on the capability development and assessment case development for the out-of-pile transient experiments, are summarized. The liquid cladding penetration model was improved based on a comprehensive set of fuel behavior test apparatus (FBTA) test results. The new model was implemented into BISON along with a simple model to account for fuel melting effects. The metallic fuel meshing capability was also upgraded to allow local refinement using unstructured triangular mesh. Based on the previous success in developing the BISON-FIPD integration powered BISON assessment case, this database enabled BISON evaluation and V&V capability was expanded to the out-of-pile transient experiments by introducing BISON-OPTD integration. A series of initial out-of-pile transient assessment cases were implemented into BISON to provide powerful platforms to evaluate current BISON’s capabilities in simulating metallic fuel performance during out-of-pile transients. The established capabilities will also benefit the calibration and V&V of future advanced mechanistic models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗