Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code Case”

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 109 records · Page 6

MRP: Design Rules for Refractory Metals

Presentation on Design rules for refractory metals to be given at the Joint ART Materials/AMMT Program Review at DOE Headquarters, Germantown, MD, June 5-8, 2023. Includes microreactor program work for structural materials, failure points of interest, technology, technology maturation and de-risk designs, refractory metals, plan for developing code case as prototype by advanced reactor developers which meets design requirements, and also provides ASME Section III rules.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Irradiation Testing of Additively Manufactured Materials for LWR Applications: Alloy 718 and 316L Stainless Steel

The objective of this NSUF Project is to assess the changes in irradiated additively manufactured (AM) material properties as compared to non-irradiated material. Type 316L stainless steel and Alloy 718 samples were produced using Direct Metal Laser Melting (DMLM) fabrication. Materials produced from this fabrication method have several potential applications within the nuclear industry as reactor internal repair parts, fuel debris resistant filters, or fuel spacers within existing light water reactors (LWRs). AM materials have been shown to achieve equivalent mechanical behavior in simulated reactor environments as compared to wrought materials, but have significantly more flexibility when it comes to unique design features. The increased component design flexibility makes these AM materials an attractive choice for both current LWR applications as well as for small modular reactor (SMR) designs. Prior to use of these materials in reactor fleet operation, the industry as a whole must evaluate the effects of irradiation on their material properties. Standard 0.4 inch thick Compact Tension specimens and SSJ3 type tensile bars were neutron irradiated at the Advanced Test Reactor to ~1 dpa for the purpose of performing a variety of mechanical tests in a range of simulated environments applicable to LWRs. For the ductile austenitic Type 316L stainless steel, the irradiated data will be used to confirm that the AM process produces materials with properties that are equivalent to wrought materials under testing conditions applicable to LWR operation. Transmission electron microscopy analysis was also performed in order to understand microstructural and microchemical changes induced in each material in response to neutron irradiation. If possible, data collected from these AM 316L samples will be used to remove fluence limits from specifications of ASME code cases for this alloy, which will give vendors much more flexibility in building future components.

36 - MATERIALS SCIENCE↗

ECAR-6580 Rev 0 ASME Section III, Division 5 Analysis of the MARVEL PCS and GVS Top Corner

The purpose of this Engineering Calculations and Analysis Report (ECAR) is to document the structural evaluation for part of the Microreactor Applications Research Validation and Evaluation (MARVEL) Primary Coolant System (PCS) and Guard Vessel System (GVS). For the PCS this relates specifically to the Distribution Plenum (DP), The Intermediate Heat Exchanger (IHX), and the Upper Downcomer (UD). For the GVS, this will include only the top corner that is machined into the DP Top Plate. These components will be evaluated using the 2021 version of ASME Section III, Division 5 [1] which is a design code that governs the construction of vessels, storage tanks, piping, pumps, valves, supports, core support structures and nonmetallic core components for use in high temperature reactor systems and their supporting systems. Materials at high temperature are subject to creep and fatigue mechanisms that require additional analyses that aren’t covered in ASME Section III, Division 1 rules. Division 5 contains two approaches: elastic or inelastic, however, additional Code Cases specific to Division 5, allow for an Elastic-Perfectly Plastic (EPP) approach. The components analyzed in this ECAR will use a combination of the elastic and EPP approach. Only Design and Service Levels A and B are evaluated in this ECAR. Service Level D evaluations for the entire PCS are documented in ECAR-6564, “MARVEL Project Primary Coolant System Pressure Vessel Stress Documentation” [2] and for the entire GVS are documented in ECAR-6574, “MARVEL Guard Vessel System FEA and ASME Analysis” [3]. The Design and Service Level A and B analyses for the PCS Downcomer piping and Core Barrel are documented in ANL-24/36, "Engineering Calculations and Analysis of the Core Barrel and Downcomer Piping in the MARVEL PCS” [4], and the GVS (except for the top corner) is in ECAR-6574, “MARVEL Guard Vessel System FEA and ASME Analysis.” See Section 2.0 for more detail on the analysis boundaries.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

HPC-driven computational reproducibility in numerical relativity codes: a use case study with IllinoisGRMHD

Abstract Reproducibility of results is a cornerstone of the scientific method. Scientific computing encounters two challenges when aiming for this goal. Firstly, reproducibility should not depend on details of the runtime environment, such as the compiler version or computing environment, so results are verifiable by third-parties. Secondly, different versions of software code executed in the same runtime environment should produceconsistent numerical results for physical quantities. In this manuscript, we test the feasibility of reproducing scientific results obtained using theIllinoisGRMHDcode that is part of an open-source community software for simulation in relativistic astrophysics, theEinstein Toolkit. We verify that numerical results of simulating a single isolated neutron star withIllinoisGRMHDcan be reproduced, and compare them to results reported by the code authors in 2015. We use two different supercomputers: Expanse at SDSC, and Stampede2 at TACC. By compiling the source code archived along with the paper on both Expanse and Stampede2, we find thatIllinoisGRMHDreproduces results published in its announcement paper up to errors comparable to round-off level changes in initial data parameters. We also verify that a current version ofIllinoisGRMHDreproduces these results once we account for bug fixes which have occurred since the original publication.

Astronomy & Astrophysics↗

The Continuous Intercomparison of Radiation Codes (CIRC): Phase I Cases

CIRC aspires to be the successor to ICRCCM (Intercomparison of Radiation Codes in Climate Models). It is envisioned as an evolving and regularly updated reference source for GCM-type radiative transfer (RT) code evaluation with the principle goal to contribute in the improvement of RT parameterizations. CIRC is jointly endorsed by DOE's Atmospheric Radiation Measurement (ARM) program and the GEWEX Radiation Panel (GRP). CIRC's goal is to provide test cases for which GCM RT algorithms should be performing at their best, i.e, well characterized clear-sky and homogeneous, overcast cloudy cases. What distinguishes CIRC from previous intercomparisons is that its pool of cases is based on observed datasets. The bulk of atmospheric and surface input as well as radiative fluxes come from ARM observations as documented in the Broadband Heating Rate Profile (BBHRP) product. BBHRP also provides reference calculations from AER's RRTM RT algorithms that can be used to select the most optimal set of cases and to provide a first-order estimate of our ability to achieve radiative flux closure given the limitations in our knowledge of the atmospheric state.

Oreopoulos, Lazaros↗

HINCOF-1: a Code for Hail Ingestion in Engine Inlets

One of the major concerns during hail ingestion into an engine is the resulting amount and space- and time-wise distribution of hail at the engine face for a given geometry of inlet and set of atmospheric and flight conditions. The appearance of hail in the capture streamtube is invariably random in space and time, with respect to size and momentum. During the motion of a hailstone through an inlet, a hailstone undergoes several processes, namely impact with other hailstones and material surfaces of the inlet and spinner, rolling and rebound following impact; heat and mass transfer; phase change; and shattering, the latter three due to friction and impact. Taking all of these factors into account, a numerical code, designated HINCOF-I, has been developed for determining the motion hailstones from the atmosphere, through an inlet, and up to the engine face. The numerical procedure is based on the Monte-Carlo method. The report presents a description of the code, along with several illustrative cases. The code can be utilized to relate the spinner geometry - conical or, more effective, elliptical - to the possible diversion of hail at the engine face into the bypass stream. The code is also useful for assessing the influence of various hail characteristics on the ingestion and distribution of hailstones over the engine face.

Gopalaswamy, N.↗

High speed corner and gap-seal computations using an LU-SGS scheme

The hybrid Lower-Upper Symmetric Gauss-Seidel (LU-SGS) algorithm was added to a widely used series of 2D/3D Euler/Navier-Stokes solvers and was demonstrated for a particular class of high-speed flows. A limited study was conducted to compare the hybrid LU-SGS for approximate Newton iteration and diagonalized Beam-Warming (DBW) schemes on a work and convergence history basis. The hybrid LU-SGS algorithm is more efficient and easier to implement than the DBW scheme originally present in the code for the cases considered. The code was validated for the hypersonic flow through two mutually perpendicular flat plates and then used to investigate the flow field in and around a simplified scramjet module gap seal configuration. Due to the similarities, the gap seal flow was compared to hypersonic corner flow at the same freestream conditions and Reynolds number.

Coirier, William J.↗

High speed corner and gap-seal computations using an LU-SGS scheme

The hybrid Lower-Upper Symmetric Gauss-Seidel (LU-SGS) algorithm was added to a widely used series of 2D/3D Euler/Navier-Stokes solvers and was demonstrated for a particular class of high-speed flows. A limited study was conducted to compare the hybrid LU-SGS for approximate Newton iteration and diagonalized Beam-Warming (DBW) schemes on a work and convergence history basis. The hybrid LU-SGS algorithm is more efficient and easier to implement than the DBW scheme originally present in the code for the cases considered. The code was validated for the hypersonic flow through two mutually perpendicular flat plates and then used to investigate the flow field in and around a simplified scramjet module gap seal configuration. Due to the similarities, the gap seal flow was compared to hypersonic corner flow at the same freestream conditions and Reynolds number.

Coirier, William J.↗

Development of a CFD code for internal flows in liquid fueled engines

To support the design efforts of engines in which liquid propellants are used, one is often required to analyze incompressible, two-dimensional or axisymmetric flows within ducts and cavities with rotating walls and complicated geometries. The steady-state solution is of interest in most cases. This code is intended to provide a tool for efficient computational fluid dynamics (CFD) analysis of such flow problems, taking advantage of the artificial compressibility concept and the Beam-Warming numerical scheme modified for second order accurate, implicit boundary conditions. These concepts ensure a stable, robust, and accurate algorithm due to the reduced speed of sound and the accuracy of the boundary conditions. The code is dedicated only to two-dimensional or axisymmetric flows with or without swirl. Three-dimensional computation is excluded to increase efficiency and speed. This paper briefly presents the theory of the code, as well as several benchmark applications with comparison to well known analytical solutions. In all these test cases, the code produced accurate results.

Dakhoul, Youssef↗

Perturbation of a normal shock - A test case for unsteady Euler codes

Shock tracking capabilities of modern nonoscillatory schemes for solving the Euler equations are studied by comparing their solutions to solutions obtained using the method of characteristics on two test problems. Both test problems involve accelerating and decelerating shock motions. The schemes tested include first order upwind, second order SONIC-A, third order ENO, and third order TONIC-A. Results show that the second and higher order schemes give nearly identical shock trajectories that are closer to the characteristic solution than the first order upwind scheme. All schemes approach the characteristic solution upon mesh refinement.

Suresh, Ambady↗

Deriving Safety Cases for the Formal Safety Certification of Automatically Generated Code

We present an approach to systematically derive safety cases for automatically generated code from information collected during a formal, Hoare-style safety certification of the code. This safety case makes explicit the formal and informal reasoning principles, and reveals the top-level assumptions and external dependencies that must be taken into account; however, the evidence still comes from the formal safety proofs. It uses a generic goal-based argument that is instantiated with respect to the certified safety property (i.e., safety claims) and the program. This will be combined with a complementary safety case that argues the safety of the framework itself, in particular the correctness of the Hoare rules with respect to the safety property and the trustworthiness of the certification system and its individual components. Keywords: Automated code generation, Hoare logic, formal code certification, safety case, Goal Structuring Notation.

Basir, Nurlida↗

Porting a fast-running wildfire plume model for use in national security and climate applications

In this project, the code for the fast-running Freitas wildfire plume model was successfully ported into the NARAC system, creating an operations-ready capability for predicting the transport of hazardous smoke from large wildfires. “Porting” refers to copying a section of code from one software package to another while accounting for unique dependencies in each software package. In this case, the code was ported from the Weather Research and Forecasting (WRF) Model, which is used extensively by the project team. To verify correct implementation, model results from the updated NARAC system were compared to those from WRF and the High-Resolution Rapid Refresh (HRRR) forecast model, which also employs the Freitas scheme. The NARAC User Interface (UI) and code documentation were also updated to include the new wildfire plume option, and a wildfire case was added to the NARAC test suite. The code was used in a real-time, semi-operational setup to predict smoke transport during the Cerro Pelado fire near Los Alamos National Laboratory in May 2022. The Freitas plume model was also tested for LLNL climate applications. Although the code was not ported into the Energy Exascale Earth System Model (E3SM), it was tested offline using E3SM atmospheric inputs, and results were compared to existing high-resolution plume predictions from large-eddy simulations in WRF. The code implementation and testing completed during this project will enable future work predicting smoke transport from large wildfires in LLNL national security and climate mission spaces.

54 ENVIRONMENTAL SCIENCES↗

A supersonic three-dimensional code for flow over blunt bodies: Program documentation and test cases

The use of a computer code for the calculation of steady, supersonic, three dimensional, inviscid flow over blunt bodies is illustrated. Input and output are given and explained for two cases: a pointed code of 20 deg half angle at 15 deg angle of attack in a free stream with M sub infinite = 7, and a cone-ogive-cylinder at 10 deg angle of attack with M sub infinite = 2.86. A source listing of the computer code is provided.

Chaussee, D. S.↗

Summary of LANL Critical Benchmark Comparison Study and Revisions for Cases Involving HEU, LEU, MIX, and Pu

This report documents results obtained for revisions made to cases involving Highly Enriched Uranium (HEU), Intermediate Enriched Uranium (IEU), a mixture of Pu and Uranium (MIX), as well as Pu cases. A previous summary of revisions for HEU an Pu cases was reported and additional investigations into four cases originally presented therein uncovered further revisions which led to better agreement with other transport codes, those cases are updated in this report. The summary of all cases reported in Reference 2 is updated in this report. In addition, a previous summary of revisions for LEU and MIX was reported, a summary of those revisions in reproduced in this report for a comprehensive summary of changes to benchmarks beginning in fiscal year 2020 to current date. The report focuses on the changes made to LANL benchmarks modeled with MCNP6 using ENDF/B-VII.1 nuclear data that appeared to have discrepant results when compared with results of other codes. Feedback was used to pinpoint review of benchmark input files and to revise them when necessary. This report documents the results of review and revision of specific benchmarks highlighted as possibly discrepant in the comparison study. In addition, there is an effort tied to this work involving collaboration between LANL XCP and NCS Divisions in the development of a shared review/revision procedure and use of a new benchmark repository. LANL has a benchmark library of critical experiments from the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook modeled for use with MCNP. This collection is now over 1100 benchmarks, referred to as the Whisper-1.1 library because it is used with the sensitivity/uncertainty package, Whisper, which supports nuclear criticality safety validation and is released with MCNP6.2. The collection, originally created several decades ago, is a combination of smaller collections, which has been revised and expanded, by various groups at LANL over the years. The original authors are no longer at the laboratory and little formal documentation of review and revision of these benchmarks exists today. A branch of the benchmark collection was already the subject of a formal review undertaken by the LANL NCS Division and expanded to include XCP Division.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Takeoff predictions for powered-lift aircraft

Takeoff predictions for powered-lift short takeoff (STO) and conventional takeoff (CTO) aircraft have been added to NASA Ames Research Center's Aircraft Synthesis (ACSYNT) code. The new computer code predicts the aircraft engine and nozzle settings required to achieve the minimum takeoff roll. As a test case, the code predicted takeoff ground rolls and nozzle settings for the YAV-8B Harrier that compared well with measured values. Brief analysis of takeoff performance for an Ejector, Remote Augmented Lift, Hybrid-Tandem Fan, and Vectored Thrust STO aircraft using the new routine will be presented.

Wardwell, Douglas A.↗

Verification of the Shafranov shift in free-boundary VMEC, DESC, and SPEC calculations with elliptical geometry

Here, given a set of current-carrying filaments that create a magnetic field with rotating elliptical flux surfaces, we present verification calculations using the VMEC, DESC, and SPEC free-boundary magnetohydrostatic equilibrium codes for both the “vacuum” equilibrium, for which there is no plasma pressure and no plasma currents, and for a non-zero pressure, zero toroidal current equilibrium. For the vacuum case, these codes are quantitatively compared to the magnetic field produced by the coils. For the non-zero pressure case, they are compared to each other. As the stepped pressure profile used in SPEC approaches the continuous profile used in VMEC and DESC, the SPEC magnetic axis approaches the VMEC and DESC axes. To within an error given by the difference between the pressure profiles, the VMEC, DESC, and SPEC calculations give the same pressure-induced geometric shift of the magnetic axis, known as the Shafranov shift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the proper Mach number and ratio of specific heats for modeling the Venus bow shock

Observational data from the Pioneer Venus Orbiter are used to investigate the physical characteristics of the Venus bow shock, and to explore some general issues in the numerical simulation of collisionless shocks. It is found that since equations from gas-dynamic (GD) models of the Venus shock cannot in general replace MHD equations, it is not immediately obvious what the optimum way is to describe the desired MHD situation with a GD code. Test case analysis shows that for quasi-perpendicular shocks it is safest to use the magnetospheric Mach number as an input to the GD code. It is also shown that when comparing GD predicted temperatures with MHD predicted temperatures total energy should be compared since the magnetic energy density provides a significant fraction of the internal energy of the MHD fluid for typical solar wind parameters. Some conclusions are also offered on the properties of the terrestrial shock.

Tatrallyay, M.↗

Toward a CFD nose-to-tail capability - Hypersonic unsteady Navier-Stokes code validation

Computational fluid dynamics (CFD) research for hypersonic flows presents new problems in code validation because of the added complexity of the physical models. This paper surveys code validation procedures applicable to hypersonic flow models that include real gas effects. The current status of hypersonic CFD flow analysis is assessed with the Compressible Navier-Stokes (CNS) code as a case study. The methods of code validation discussed to beyond comparison with experimental data to include comparisons with other codes and formulations, component analyses, and estimation of numerical errors. Current results indicate that predicting hypersonic flows of perfect gases and equilibrium air are well in hand. Pressure, shock location, and integrated quantities are relatively easy to predict accurately, while surface quantities such as heat transfer are more sensitive to the solution procedure. Modeling transition to turbulence needs refinement, though preliminary results are promising.

Edwards, Thomas A.↗