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 865 records · Page 48

NASA Engineering and Safety Center Technical Bulletin No. 21-05: Industry Recommendations from Arecibo Observatory Zinc Spelter Socket Joint Failure Analysis

A structural analysis and forensic investigation concluded that the Arecibo Observatory M4N socket joint failure in August 2020 was primarily due to cumulative damage caused by initially low structural design margins and a high percentage of sustained load, resulting in zinc creep deformation, progressive internal socket wire damage, and eventual loss of joint capacity. Open spelter sockets of this type are used throughout industry in stay cables. Recommendations are proposed to prevent failures of similar socket joints, including verification of positive stress margins in socket joint wires for all failure modes, periodic visual inspections with pass/fail criteria for zinc extrusion that are tied to structural qualification, and revisiting codes/ industry standards to capture lessons learned.

Arecibo Observatory↗

SAM Two-Phase Flow Model Development and Applications for Operational Transients in Advanced Reactors

As advanced nuclear technologies continue to develop, the need for the flexible operation and generation of these advanced reactors becomes necessary to maximize economic potential. As large-scale experiments are not always feasible, modeling and simulations of advanced reactors play a crucial role in design optimization and analysis. The SAM (System Analysis Module) code developed at Argonne National Laboratory is a state-of-the-art system-level thermal-hydraulic code aimed at simulating advanced reactor systems. Recent code developments have implemented two-phase flow modeling using the homogeneous equilibrium model, and a new steam generator component has been developed to utilize the two-phase flow implementation. In addition to verification tests, a load-following simulation was performed to model a realistic load-following transient in a proposed integrated system consisting of a conceptual advanced reactor known as the Advanced Burner Test Reactor (ABTR) and thermal energy storage (TES) tanks. The integrated system model uses two large TES tanks designed for sodium and a model helical coil steam generator to simulate the operational load-following transient. The flow rates of the feedwater and secondary loops are regulated to meet a prescribed steam generator load consistent with the electricity demand over a 24-h period. In conclusion, the results found the ABTR system was able to maintain stable reactor conditions and primary- and secondary-side characteristics over the course of the load-following transient.

Advanced Burner Test Reactor (ABTR)↗

HARD: A performance portable radiation hydrodynamics code based on FleCSI framework

Hydrodynamics And Radiation Diffusion (HARD) is an open-source application for high-performance simulations of compressible hydrodynamics with radiation-diffusion coupling. Built on the FleCSI (Bergen et al., 2021 [1]) (Flexible Computational Science Infrastructure) framework, HARD expresses its computational units as tasks whose execution can be orchestrated by multiple back-end runtimes, including Legion (Bauer et al., 2012 [2]), MPI (Forum, 1994 [3]), and HPX (Kaiser et al., 2020 [4]). Node-level parallelism is handled through Kokkos (Edwards et al., 2014 [5]), providing a single-source, portable code base that runs efficiently on laptops, small homogeneous clusters, and the largest heterogeneous supercomputers currently available. To ensure scientific reliability, HARD includes a regression test suite that automatically reproduces canonical verification problems such as the Sod and LeBlanc shock tubes, and the Sedov blast wave, comparing numerical solutions against known analytical results. The project is distributed under an OSI-approved license, hosted on GitHub, and accompanied by reproducible build scripts and continuous integration workflows. This combination of performance portability, verification infrastructure, and community-focused development makes HARD a sustainable platform for advancing radiation hydrodynamics research across multiple domains.

97 MATHEMATICS AND COMPUTING↗

An Integral-based Technique to Accelerate the Monte Carlo Radiative Transfer Computation for Supernovae

We present an integral-based technique (IBT) algorithm to accelerate supernova (SN) radiative transfer calculations. The algorithm utilizes “integral packets,” which are calculated by the path integral of the Monte Carlo (MC) energy packets, to synthesize the observed spectropolarimetric signal at a given viewing direction in a 3D time-dependent radiative transfer program. Compared to the event-based technique (EBT) proposed by M. Bulla et al., our algorithm significantly reduces the computation time and increases the MC signal-to-noise ratio (S/N). Using a 1D spherical symmetric Type Ia SN ejecta model DDC10 and its derived 3D model, the IBT algorithm has successfully passed the verification of spherical symmetry and cross comparison on a 3D SN model with the direct-counting technique and EBT. Notably, with our algorithm implemented in the 3D MC radiative transfer code SEDONA, the computation time is faster than EBT by a factor of 10−30, and the S/N is better by a factor of 1.5−3, with the same number of MC quanta.

79 ASTRONOMY AND ASTROPHYSICS↗

Dispersed Fringe Sensing Analysis - DFSA

Dispersed Fringe Sensing (DFS) is a technique for measuring and phasing segmented telescope mirrors using a dispersed broadband light image. DFS is capable of breaking the monochromatic light ambiguity, measuring absolute piston errors between segments of large segmented primary mirrors to tens of nanometers accuracy over a range of 100 micrometers or more. The DFSA software tool analyzes DFS images to extract DFS encoded segment piston errors, which can be used to measure piston distances between primary mirror segments of ground and space telescopes. This information is necessary to control mirror segments to establish a smooth, continuous primary figure needed to achieve high optical quality. The DFSA tool is versatile, allowing precise piston measurements from a variety of different optical configurations. DFSA technology may be used for measuring wavefront pistons from sub-apertures defined by adjacent segments (such as Keck Telescope), or from separated sub-apertures used for testing large optical systems (such as sub-aperture wavefront testing for large primary mirrors using auto-collimating flats). An experimental demonstration of the coarse-phasing technology with verification of DFSA was performed at the Keck Telescope. DFSA includes image processing, wavelength and source spectral calibration, fringe extraction line determination, dispersed fringe analysis, and wavefront piston sign determination. The code is robust against internal optical system aberrations and against spectral variations of the source. In addition to the DFSA tool, the software package contains a simple but sophisticated MATLAB model to generate dispersed fringe images of optical system configurations in order to quickly estimate the coarse phasing performance given the optical and operational design requirements. Combining MATLAB (a high-level language and interactive environment developed by MathWorks), MACOS (JPL s software package for Modeling and Analysis for Controlled Optical Systems), and DFSA provides a unique optical development, modeling and analysis package to study current and future approaches to coarse phasing controlled segmented optical systems.

Sigrist, Norbert↗

Dynamic Modelling of HiveWind Floating Wind Substructure in OpenFAST

Current aero-elastic modes tend to consider the wind turbine floating substructures as rigid. However, the new generation of offshore wind platforms are leaner, less rigid. The OpenFAST code has been updated to consider the flexibility of the floating substructure. This presentation shows how a numerical model of the HiveWind design was created within OpenFAST as well as some verifications performed. When comparing the flexible model to a simplified rigid model (previous state-of-the-art), there are some significant differences in key performance indicators. The new flexibility capability within OpenFAST allows to calculate ultimate and fatigue loads within the floating substructure, allowing for the development of more economical designs.

ENGINEERING,WIND ENERGY↗

Application Program Interface for the Orion Aerodynamics Database

The Application Programming Interface (API) for the Crew Exploration Vehicle (CEV) Aerodynamic Database has been developed to provide the developers of software an easily implemented, fully self-contained method of accessing the CEV Aerodynamic Database for use in their analysis and simulation tools. The API is programmed in C and provides a series of functions to interact with the database, such as initialization, selecting various options, and calculating the aerodynamic data. No special functions (file read/write, table lookup) are required on the host system other than those included with a standard ANSI C installation. It reads one or more files of aero data tables. Previous releases of aerodynamic databases for space vehicles have only included data tables and a document of the algorithm and equations to combine them for the total aerodynamic forces and moments. This process required each software tool to have a unique implementation of the database code. Errors or omissions in the documentation, or errors in the implementation, led to a lengthy and burdensome process of having to debug each instance of the code. Additionally, input file formats differ for each space vehicle simulation tool, requiring the aero database tables to be reformatted to meet the tool s input file structure requirements. Finally, the capabilities for built-in table lookup routines vary for each simulation tool. Implementation of a new database may require an update to and verification of the table lookup routines. This may be required if the number of dimensions of a data table exceeds the capability of the simulation tools built-in lookup routines. A single software solution was created to provide an aerodynamics software model that could be integrated into other simulation and analysis tools. The highly complex Orion aerodynamics model can then be quickly included in a wide variety of tools. The API code is written in ANSI C for ease of portability to a wide variety of systems. The input data files are in standard formatted ASCII, also for improved portability. The API contains its own implementation of multidimensional table reading and lookup routines. The same aerodynamics input file can be used without modification on all implementations. The turnaround time from aerodynamics model release to a working implementation is significantly reduced

Robinson, Philip E.↗

Verification Test Suite for Spalart-Allmaras QCR2000 Turbulence Model

The paper presents three benchmark cases for verification of Reynolds-averaged Navier-Stokes solvers. The verification studies focus on a one-equation Spalart-Allmaras model, SA-[neg]-QCR2000, that uses a version of quadratic constitutive relations.The benchmark cases are a two-dimensional subsonic flow around a Joukowski airfoil, a three-dimensional supersonic flow through a square duct, and a three-dimensional flow over a wing-fuselage configuration. The turbulence-model formulation, geometry, flow conditions, grids, and the expected output are described in detail. Reference solutions computed by several established codes are shown

Computational Aerodynamics↗

Calculation of collisionless pitch-angle scattering of runaway electrons with synchrotron radiation via high-order guiding-centre equation

Recently, the collisionless pitch-angle scattering for relativistic runaway electrons (REs) in toroidal geometries such as tokamaks was discovered through a full orbit simulation approach (Liu et al., Nucl. Fusion, vol. 56, 2016, p. 064002), and it was then theoretically investigated that a new expression for the magnetic moment, including the second-order corrections, could essentially reproduce the so-called collisionless pitch-angle scattering process (Liu et al., Nucl. Fusion, vol. 58, 2018, p. 106018). In this paper, with synchrotron radiation, extensive numerical verification of the validity of the high-order guiding-centre theory is given for simulations involving REs by incorporating such an expression for the magnetic moment into our particle tracing code. A high-order guiding-centre simulation approach with synchrotron radiation (HGSA) is applied. Synchrotron radiation plays an essential role in the life cycle of REs. The energy of REs first increases and then becomes saturated until the electric field acceleration is balanced by the radiation dissipation. Unfortunately, the process cannot be simulated accurately with the standard guiding-centre model, i.e. the first-order guiding-centre model. Remarkably, it is found that the HGSA can effectively produce the fundamental process of REs. Since the time scale of the energy saturation of REs is close to seconds, the computational cost becomes significant. In order to save costs, it is necessary to estimate the time of energy saturation. An analytical estimate is derived for the time it takes for synchrotron drag to balance an accelerating electric field and the provided formula has been numerically verified. Finally, test calculations reveal that HGSA is favourable for exploiting the dynamics of REs in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Navigating United States Standards and Regulation for Digital Energy Systems

This report provides an analysis of the U.S. standards and regulatory landscape for digital energy systems, focusing on cybersecurity, safety, and reliability requirements. It examines the interplay between federal mandates, state regulations, voluntary industry standards, and utility-specific policies, highlighting critical gaps between compliance and real-world risk mitigation. While NERC CIP standards enforce cybersecurity for Bulk Electric System assets, distribution-level infrastructure and emerging technologies often fall outside mandatory oversight, creating vulnerabilities. The report identifies systemic challenges such as reliance on self-attestation, uneven state adoption of safety codes, and lagging standards for advanced technologies like battery energy storage and inverter-based resources. Through a detailed gap analysis, it underscores the urgency of proactive risk-based approaches, independent verification, and strategic engagement with state and federal entities. Recommendations include adopting tiered security frameworks, strengthening procurement practices, and addressing emerging technology risks to ensure resilient and secure digital energy infrastructure. This guidance is intended for utilities, regulators, and stakeholders navigating compliance obligations and seeking to enhance cybersecurity and safety beyond minimum standards.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

An expert system development methodology which supports verification and validation

Expert systems have demonstrated commercial viability in a wide range of applications, but still face some obstacles to widespread use. A major stumbling block is the lack of well defined verification and validation (V and V) techniques. The primary difficulty with expert system V and V is the use of development methodologies which do not support V and V. As with conventional code, the key to effective V and V is the development methodology. An expert system development methodology is described which is based upon a panel review approach, that allows input from all parties concerned with the expert system.

Culbert, Chris↗

Overview of the instrumentation program

This program is aimed at developing sensors and measurement systems capable of obtaining the data necessary for the verification of computational models of the structural behavior, the fatigue life, and the environmental conditions pertinent to advanced reusable space propulsion systems. One of the characteristics of measurement systems needed to verify codes is that the sensors must be nonintrusive or at least minimally intrusive so as not to significantly perturb the conditions being measured. This leads to a heavy emphasis on laser optical techniques and on thin-film sensors. Another characteristic of such instruments is that they must be highly accurate and produce very high spatial and temporal resolution of the parameter being measured. The measurement systems needed generally fall into a number of broad categories. First there are the measurements needed on the surfaces of components such as turbine blades and vanes. Some of the desired parameters are temperature, strain, and heat flux. Another broad category encompasses those measurements needed in the flow environment around these components. Here, the desired results are high resolution maps of such parameters as flow velocity, temperature, density, pressure, and species concentration. The remaining category deals with measurements necessary for monitoring the health of the engine. This category has loomed ever more important since the Challenger disaster. An optical method for determining the characteristics of the plume is presented. Holographic measurement of structural damage is also presented.

Nieberding, William C.↗

Helium 2 slosh in low gravity

This paper describes the status and plans for the work being performed under NASA NRA contract NASW-4803 so that members of the Microgravity Fluid Dynamics Discipline Working Group are aware of this program. The contract is a cross-disciplinary research program and is administered under the Low Temperature Microgravity Research Program at the Jet Propulsion Laboratory. The purpose of the project is to perform low-gravity verification experiments on the slosh behavior of He II to use in the development of a CFD model that incorporates the two-fluid physics of He II. The two-fluid code predicts a different fluid motion response in low-gravity environment from that predicted by a single-fluid model, while the 1g response is identical for the both types of model.

Ross, Graham O.↗

Precise and Scalable Static Program Analysis of NASA Flight Software

Recent NASA mission failures (e.g., Mars Polar Lander and Mars Orbiter) illustrate the importance of having an efficient verification and validation process for such systems. One software error, as simple as it may be, can cause the loss of an expensive mission, or lead to budget overruns and crunched schedules. Unfortunately, traditional verification methods cannot guarantee the absence of errors in software systems. Therefore, we have developed the CGS static program analysis tool, which can exhaustively analyze large C programs. CGS analyzes the source code and identifies statements in which arrays are accessed out of bounds, or, pointers are used outside the memory region they should address. This paper gives a high-level description of CGS and its theoretical foundations. It also reports on the use of CGS on real NASA software systems used in Mars missions (from Mars PathFinder to Mars Exploration Rover) and on the International Space Station.

Brat, G.↗

Modeling of resistive plasma response in toroidal geometry using an asymptotic matching approach

The method of solving linear resistive plasma response, based on the asymptotic matching approach, is developed for full toroidal tokamaks by upgrading the Resistive DCON code [A.H. Glasser, Z.R. Wang and J.-K. Park, Physics of Plasmas, \textbf{23}, 112506 (2016)]. The derived matching matrix, asymptotically matching the outer and inner regions, indicates that the applied three dimension (3-D) magnetic perturbations contribute additional small solutions at each resonant surface due to the toroidal coupling of poloidal modes. In contrast, the resonant harmonic only affects the corresponding resonant surface in the cylindrical plasma. Since the solution of ideal outer region is critical to the asymptotic matching and is challenging to be solved in the toroidal geometry due to the singular power series solution at the resonant surfaces, systematic verification of the outer region $\Delta^\prime$ matrix is made by reproducing the well known analytical $\Delta^{\prime}$ result in [H.P. Furth, P.H. Rutherford and H. Selberg, The Physics of Fluids, \textbf{16}, 1054-1063 (1073)] as well as by making a quantitative benchmark with the PEST3 code [A. Pletzer and R.L. Dewar, J. Plasma Physics, \textbf{45}, 427-451 (1991)]. Finally, the reconstructed numerical solution of resistive plasma response from the toroidal matching matrix is presented. Comparing with the ideal plasma response, the global structure of the response can be affected by the small finite island at the resonant surfaces.

MHD simulation↗

Partnership Center for High-Fidelity Boundary Plasma Simulation (Final Report)

Within the Partnership Center for High-Fidelity Boundary Plasma Simulation (HBPS), work at UT-Austin was aimed at improved verification, validation, and uncertainty quantification (VVUQ) for edge plasma simulations and on performing gyrokinetics simulations of pedestal instabilities and turbulence in order to expand foundational understanding of pedestal transport. Regarding VVUQ, the accomplishments can be summarized as follows. First, it was shown that the Moment Preserving Constrained Resampling technique, when applied periodically in particle-in-cell simulations in the XGC code, can dramatically improve the accuracy of the simulation at essentially equivalent computational cost. Second, a technique for estimating model correlations, which are required to solve the model selection and sample allocation problem in multifidelity UQ techniques, without sampling the highest fidelity, most computationally expensive model, was developed and demonstrated. Third, previously developed methods for estimating statistical and discretization errors were applied to numerical methods relevant to edge plasma simulations, namely in particle-in-cell-based approaches, and shown to work. Finally, benchmark studies for comparing gyrokinetic codes were developed and performed, leading to reasonable agreement between four commonly used codes. Regarding physics studies, gyrokinetic simulations to investigate microtearing modes in the DIII-D pedestal were performed using the GENE code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Derivation and verification of the direct-sampling method for simulating Monte Carlo flight paths in tetrahedral meshes with linear finite-element cross sections

This paper provides a derivation of a direct-sampling approach for modeling continuously varying cross sections in tetrahedral-mesh-based Monte Carlo codes. Specifically, cross sections are spatially approximated using linear nodal finite elements. A linearization strategy is provided for non-linearly varying cross sections. The method is verified against seven analytical pure-absorber test problems. These test problems also highlight the benefit of using linear finite elements over element-wise-constant cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluating wind speed and power forecasts for wind energy applications using an open-source and systematic validation framework

Building on the verification and validation work developed under the Second Wind Forecast Improvement Project, this work exhibits the value of a consistent procedure to evaluate wind power forecasts. We established an open-source Python code base tailored for wind speed and wind power forecast validation, WE-Validate. The code base can evaluate model forecasts with observations in a coherent manner. To demonstrate the systematic validation framework of WE-Validate, we designed and hosted a forecast evaluation benchmark exercise. We invited forecast providers in industry and academia to participate and submit forecasts for two case studies. We then evaluated the submissions with WE-Validate. Our findings suggest that ensemble means have reasonable skills in time series forecasting, whereas they are often inferior to single ensemble members in wind ramp forecasting. Adopting a voting scheme in ramp forecasting that allows ensemble members to detect ramps independently leads to satisfactory skill scores. Throughout this document, we also emphasize the importance of using statistically robust and resistant metrics as well as equitable skill scores in forecast evaluation.

17 WIND ENERGY↗