Search NASA⌕ Search

SEARCH · Search NASA

Results for “Software dependability”

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 1,405 records · Page 78

High-Fidelity Numerical Investigation on Elucidating Sodium Heat Transfer Characteristics for 37-Pin Wire-Wrapped Fuel Bundle in the PLANDTL Facility

This study involved a Reynolds-averaged Navier-Stokes- (RANS-) based computational fluid dynamics (CFD) analysis of the 37-pin wire-wrapped fuel bundle of the PNC Plant dynamics test loop (PLANDTL) facility. Previously, mainly the hydrodynamic phenomena of the wire-wrapped fuel bundle were analyzed, but the present study additionally included heat transfer analysis through conjugate heat transfer. The main purpose of the study was to benchmark the experimental data of the PLANDTL 37-pin wire-wrapped fuel bundle to investigate the heat transfer phenomena. In addition, the aim was to verify the accuracy of the RANS-based CFD analysis method using the STAR-CCM+ simulation software in comparison with the experimental data. The grid used for verification was an innovative grid system consisting of hexahedra using Fortran-based code. The development of the RANS-based CFD methodology included grid sensitivity analysis, turbulence model sensitivity analysis, and turbulent Prandtl number sensitivity analysis. Information on the temperature, mass flow rate, and area of the CFD results for each subchannel was provided for the top of the heated section and is expected to serve as a reference for future studies aiming to perform the validation and verification of a PLANDTL facility. In addition, the dependence of the peak temperature on the azimuth angle of each pin was analyzed.

97 MATHEMATICS AND COMPUTING↗

Design for Validation

The use of computer hardware and software in life-critical applications, such as for civil air transports, demands the use of rigorous formal mathematical validation procedures. However, formal specification and verification will only be tractable if the system is designed in a manner that lends itself to formal methods. Likewise, accurate reliability analysis will only be tractable if the number of interacting components that must be individually included in a single reliability model is kept to a low number and if their failure behavior interactions can be modeled simply. Also, the system must be designed such that the system reliability does not directly depend on system parameters that cannot be accurately determined. This paper presents a design methodology based on the concept of designing a system in such a manner that it can be rigorously validated, or "design for validation."

Johnson, Sally C.↗

Lessons Learned from OSIRIS-Rex Autonomous Navigation Using Natural Feature Tracking

The Origins, Spectral Interpretation, Resource Identification, Security-Regolith Explorer (Osiris-REx) spacecraft is scheduled to launch in September, 2016 to embark on an asteroid sample return mission. It is expected to rendezvous with the asteroid, Bennu, navigate to the surface, collect a sample (July 20), and return the sample to Earth (September 23). The original mission design called for using one of two Flash Lidar units to provide autonomous navigation to the surface. Following Preliminary design and initial development of the Lidars, reliability issues with the hardware and test program prompted the project to begin development of an alternative navigation technique to be used as a backup to the Lidar. At the critical design review, Natural Feature Tracking (NFT) was added to the mission. NFT is an onboard optical navigation system that compares observed images to a set of asteroid terrain models which are rendered in real-time from a catalog stored in memory on the flight computer. Onboard knowledge of the spacecraft state is then updated by a Kalman filter using the measured residuals between the rendered reference images and the actual observed images. The asteroid terrain models used by NFT are built from a shape model generated from observations collected during earlier phases of the mission and include both terrain shape and albedo information about the asteroid surface. As a result, the success of NFT is highly dependent on selecting a set of topographic features that can be both identified during descent as well as reliably rendered using the shape model data available. During development, the OSIRIS-REx team faced significant challenges in developing a process conducive to robust operation. This was especially true for terrain models to be used as the spacecraft gets close to the asteroid and higher fidelity models are required for reliable image correlation. This paper will present some of the challenges and lessons learned from the development of the NFT system which includes not just the flight hardware and software but the development of the terrain models used to generate the onboard rendered images.

Navigation↗

Generalized and efficient algorithm for computing multipole energies and gradients based on Cartesian tensors

Accurate representation of intermolecular forces has been the central task of classical atomic simulations, known as molecular mechanics. Recent advancements in molecular mechanics models have put forward the explicit representation of permanent and/or induced electric multipole (EMP) moments. The formulas developed so far to calculate EMP interactions tend to have complicated expressions, especially in Cartesian coordinates, which can only be applied to a specific kernel potential function. For example, one needs to develop a new formula each time a new kernel function is encountered. The complication of these formalisms arises from an intriguing and yet obscured mathematical relation between the kernel functions and the gradient operators. Here, I uncover this relation via rigorous derivation and find that the formula to calculate EMP interactions is basically invariant to the potential kernel functions as long as they are of the form f(r), i.e., any Green’s function that depends on inter-particle distance. I provide an algorithm for efficient evaluation of EMP interaction energies, forces, and torques for any kernel f(r) up to any arbitrary rank of EMP moments in Cartesian coordinates. The working equations of this algorithm are essentially the same for any kernel f(r). Recently, a few recursive algorithms were proposed to calculate EMP interactions. Depending on the kernel functions, the algorithm here is about 4–16 times faster than these algorithms in terms of the required number of floating point operations and is much more memory efficient. I show that it is even faster than a theoretically ideal recursion scheme, i.e., one that requires 1 floating point multiplication and 1 addition per recursion step. This algorithm has a compact vector-based expression that is optimal for computer programming. The Cartesian nature of this algorithm makes it fit easily into modern molecular simulation packages as compared with spherical coordinate-based algorithms. A software library based on this algorithm has been implemented in C++11 and has been released.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CRPropa 3.1—a low energy extension based on stochastic differential equations

The propagation of charged cosmic rays through the Galactic environment influences all aspects of the observation at Earth. Energy spectrum, composition and arrival directions are changed due to deflections in magnetic fields and interactions with the interstellar medium. Today the transport is simulated with different simulation methods either based on the solution of a transport equation (multi-particle picture) or a solution of an equation of motion (single-particle picture). We developed a new module for the publicly available propagation software CRPropa 3.1, where we implemented an algorithm to solve the transport equation using stochastic differential equations. This technique allows us to use a diffusion tensor which is anisotropic with respect to an arbitrary magnetic background field. The source code of CRPropa is written in C++ with python steering via SWIG which makes it easy to use and computationally fast. In this paper, we present the new low-energy propagation code together with validation procedures that are developed to proof the accuracy of the new implementation. Furthermore, we show first examples of the cosmic ray density evolution, which depends strongly on the ratio of the parallel κ{sub ∥} and perpendicular κ{sub ⊥} diffusion coefficients. This dependency is systematically examined as well the influence of the particle rigidity on the diffusion process.

79 ASTRONOMY AND ASTROPHYSICS↗

Computational Modeling of a 3D Printed Recuperator and Subsequent Experimental Loop for Supercritical Carbon Dioxide Cycles

Oak Ridge National Laboratory (ORNL), in collaboration with mechanical-thermal energy storage (mTES) provider EarthEn, a US Department of Energy (DOE) Lab-Embedded Entrepreneurship Program (LEEP) recipient at ORNL’s Innovation Crossroads 2023, is utilizing a state-of-the-art patented 3D printing technique to design an additively manufactured (AM) supercritical CO2 (sCO2) recuperator (REC) for EarthEn’s charge/discharge cycle. The AM REC will be printed at ORNL’s Manufacturing Demonstration Facility using Inconel Alloy 718 and tested on a closed-loop, ∼100 kW scale experimental facility that is under construction. The testing will compare the printed design against a commercial-off-the-shelf Printed Circuit Heat Exchanger (PCHE) REC. The design of the sCO2 facility is guided by a Modelica-based system model which is primarily dependent on the open-source TRANSFORM library developed at ORNL and uses the open-source CoolProp library for thermophysical properties of sCO2 via the External Media library. It is envisioned that an iterative process will be followed between the physical loop and the system model wherein the initial experimental data will be used to tune the model, which in turn will be used to guide future loop operation. Simultaneously, the AM REC is being designed using computer-aided design models, and it is also being analyzed for hydraulic and thermomechanical response using commercial computational fluid dynamics software, Simcenter STAR-CCM+, on highperformance computing resources.1

See, Nate [ORNL] (ORCID:0000000178581202)↗

SU-E-T-108: An Investigation of Cerenkov Light Production in the Exradin W1 Scintillator Under Various Measurement Conditions

Purpose: To investigate the effects of depth, fiber-optic cable bends, and incident radiation angle on Cerenkov production in the Standard Imaging Exradin W1. Methods: Measurements were completed using a Varian Clinac 21EX linear accelerator with an Exradin W1 scintillator as well as a cable-only scintillator (no scintillation material) to isolate the Cerenkov signal. The effects of cable bend radius and location were investigated by bending the fiber-optic cable into a circle with radii ranging from 1.0 to 10.8 cm and positioning the center of the coil at distances ranging from 10.0 to 175.0 cm from the photodiode. The effects of depth and incident radiation angle were investigated by performing measurements in water at depths ranging from 1.0 cm to 25.0 cm and angles ranging from 0° to 80°. Eclipse treatment-planning software was utilized to ensure a consistent dose was delivered to the W1 regardless of depth or angle. Results: Measured signal in both channels of the cable-only scintillator decreased as the bend radius decreased and as the distance between the bend and photodiode increased. A fiber bend of 1.0 cm radius produced a 17.1% decrease in the green channel response in the cable-only scintillator. The effect of depth was less severe; a maximum increase of 6.6% in the green channel response was observed at a depth of 25.0 cm in the W1. In the angular dependence investigation, the signal in both channels of the W1 peaked at an angle of 40°; which is in agreement with the nominal Cerenkov emission angle of 45°. Conclusion: The green channel response in the W1 (mainly Cerenkov signal) varied with depth, fiber-optic cable bends, and incident radiation angle. Fully characterizing Cerenkov production is essential to ensure it is properly accounted for in scintillator measurements. Research funding and materials received by Standard Imaging, Inc. (Middleton WI)

07 ISOTOPE AND RADIATION SOURCES↗

Current-Voltage Analysis Tool for Solar Fuel Production (CATS) v0.1

The program is expected to be a useful tool during operation of a range of solar-driven electrochemical devices, such as PV-electrolyzers and photoelectrochemical (PEC) devices. It makes use of current-voltage data collected during operation, operation parameters that are easily accessible for many devices of interest. The program captures the time-dependence of loss mechanisms at play during PEC device operation and opens the door for real-time optimizations of operational parameters like the feed humidity to minimize performance losses. Currently, the program is useful for applications where the current-voltage polarization response of the photoabsorber and electrolyzer are pre-recorded, and changes over time can be estimated based on the pre-recorded current-voltage scans. More complicated scenarios that involve permanent degradation of the PV component or complex calculations of the electrolyzer current-voltage characteristics may be included in future versions of this software. Aside from the possibility of deconvoluting loss mechanisms, this software can also reveal possible performance benefits by convective PV cooling when the photoabsorber is integrated into the electrolysis cell. The latter can highlight why certain device architectures may be beneficial over others. No similar technologies are known to the developer.

Kistler, Tobias↗

Structural Reliability Using Probability Density Estimation Methods Within NESSUS

A reliability analysis studies a mathematical model of a physical system taking into account uncertainties of design variables and common results are estimations of a response density, which also implies estimations of its parameters. Some common density parameters include the mean value, the standard deviation, and specific percentile(s) of the response, which are measures of central tendency, variation, and probability regions, respectively. Reliability analyses are important since the results can lead to different designs by calculating the probability of observing safe responses in each of the proposed designs. All of this is done at the expense of added computational time as compared to a single deterministic analysis which will result in one value of the response out of many that make up the density of the response. Sampling methods, such as monte carlo (MC) and latin hypercube sampling (LHS), can be used to perform reliability analyses and can compute nonlinear response density parameters even if the response is dependent on many random variables. Hence, both methods are very robust; however, they are computationally expensive to use in the estimation of the response density parameters. Both methods are 2 of 13 stochastic methods that are contained within the Numerical Evaluation of Stochastic Structures Under Stress (NESSUS) program. NESSUS is a probabilistic finite element analysis (FEA) program that was developed through funding from NASA Glenn Research Center (GRC). It has the additional capability of being linked to other analysis programs; therefore, probabilistic fluid dynamics, fracture mechanics, and heat transfer are only a few of what is possible with this software. The LHS method is the newest addition to the stochastic methods within NESSUS. Part of this work was to enhance NESSUS with the LHS method. The new LHS module is complete, has been successfully integrated with NESSUS, and been used to study four different test cases that have been proposed by the Society of Automotive Engineers (SAE). The test cases compare different probabilistic methods within NESSUS because it is important that a user can have confidence that estimates of stochastic parameters of a response will be within an acceptable error limit. For each response, the mean, standard deviation, and 0.99 percentile, are repeatedly estimated which allows confidence statements to be made for each parameter estimated, and for each method. Thus, the ability of several stochastic methods to efficiently and accurately estimate density parameters is compared using four valid test cases. While all of the reliability methods used performed quite well, for the new LHS module within NESSUS it was found that it had a lower estimation error than MC when they were used to estimate the mean, standard deviation, and 0.99 percentile of the four different stochastic responses. Also, LHS required a smaller amount of calculations to obtain low error answers with a high amount of confidence than MC. It can therefore be stated that NESSUS is an important reliability tool that has a variety of sound probabilistic methods a user can employ and the newest LHS module is a valuable new enhancement of the program.

Chamis, Chrisos C.↗

Mars Express Forward Link Capabilities for the Mars Relay Operations Service (MaROS)

This software provides a new capability for landed Mars assets to perform forward link relay through the Mars Express (MEX) European Union orbital spacecraft. It solves the problem of standardizing the relay interface between lander missions and MEX. The Mars Operations Relay Service (MaROS) is intended as a central point for relay planning and post-pass analysis for all Mars landed and orbital assets. Through the first two phases of implementation, MaROS supports relay coordination through the Odyssey orbiter and the Mars Reconnaissance Orbiter (MRO). With this new software, MaROS now fully integrates the Mars Express spacecraft into the relay picture. This new software generates and manages a new set of file formats that allows for relay request to MEX for forward and return link relay, including the parameters specific to MEX. Existing MEX relay planning interactions were performed via email exchanges and point-to-point file transfers. By integrating MEX into MaROS, all transactions are managed by a centralized service for tracking and analysis. Additionally, all lander missions have a single, shared interface with MEX and do not have to integrate on a mission-by mission basis. Relay is a critical element of Mars lander data management. Landed assets depend largely upon orbital relay for data delivery, which can be impacted by the availability and health of each orbiter in the network. At any time, an issue may occur to prevent relay. For this reason, it is imperative that all possible orbital assets be integrated into the overall relay picture.

Allard, Daniel A.↗

Towards elastic in situ analysis for high-performance computing simulations

In situ analysis and visualization have grown increasingly popular for enabling direct access to data from high-performance computing (HPC) simulations. As a simulation progresses and interesting physical phenomena emerge, however, the data produced may become increasingly complex, and users may need to dynamically change the type and scale of in situ analysis tasks being carried out and consequently adapt the amount of resources allocated to such tasks. To date, none of the production in situ analysis frameworks offer such an elasticity feature, and for good reason: the assumption that the number of processes could vary during run time would force developers to rethink software and algorithms at every level of the in situ analysis stack. In this paper we present Colza, a data staging service with elastic in situ visualization capabilities. We demonstrate the use of Colza with the Deep Water Impact and the AMR-Wind simulations, coupling them with the ParaView Catalyst and Ascent in situ libraries, and show that Colza enables dynamic rescaling of these widely-used frameworks with no interruption to the simulation or staging service. Here, we highlight the challenges of enabling such elasticity, which requires overcoming these frameworks' reliance on MPI, using distinct engineering approaches, namely dependency injection and dependency overload. To the best of our knowledge, this work is the first to enable elastic in situ visualization capabilities for HPC applications on top of existing production analysis tools.

97 MATHEMATICS AND COMPUTING↗

Rapid and semi-analytical design and simulation of a toroidal magnet made with YBCO and MgB 2 superconductors

Recent advances in second-generation (YBCO) high-temperature superconducting wire could potentially enable the design of super high performance energy storage devices that combine the high energy density of chemical storage with the high power of superconducting magnetic storage. However, the high aspect ratio and the considerable filament size of these wires require the concomitant development of dedicated optimization methods that account for the critical current density in type-II superconductors. In this study, we report on the novel application and results of a CPU-efficient semianalytical computer code based on the Radia 3-D magnetostatics software package. Our algorithm is used to simulate and optimize the energy density of a superconducting magnetic energy storage device model, based on design constraints, such as overall size and number of coils. The rapid performance of the code is pivoted on analytical calculations of the magnetic field based on an efficient implementation of the Biot-Savart law for a large variety of 3-D “base” geometries in the Radia package. The significantly reduced CPU time and simple data input in conjunction with the consideration of realistic input variables, such as material-specific, temperature, and magnetic-field-dependent critical current densities, have enabled the Radia -based algorithm to outperform finite-element approaches in CPU time at the same accuracy levels. Comparative simulations of MgB 2 and YBCO-based devices are performed at 4.2 K, in order to ascertain the realistic efficiency of the design configurations.

2G↗

DEMographic MicrOSimulation (DEMOS) v2.0

DEMOS is an agent-based simulation framework used to evolve population demographic characteristics or lifecycle events such as education, marital status etc. DEMOS modules are designed to capture the interdependencies of short-term and long-term lifecycle events often influential in downstream transportation and land use modeling. An important facet of DEMOS is that it is dynamic, meaning it captures the impact of an agent's demographic characteristics in year 't' on their demographic status in year 't+1'. This has important consequences on medium- and long-term downstream modeling in different domains, such as transportation where decisions such as household vehicle transactions (i.e., buying, selling, or replacing a vehicle) depend current or past household vehicle holdings as well as household transitions. Other downstream modeling that can be enabled by DEMOS include building technology adoption or retrofit decisions, and more. DEMOS modules include: individual age progression year by year, individual mortality events, household restructuring events including marriage match-making; household birth events; individual education participation; individual labor force participation; child leaving household; household in-migration and out-migration for the region; household residence location choice and household mandatory destination choice (work and school). This version of DEMOS is tailored to be integrated with the land use simulation software UrbanSim.

Caicedo, Juan [UrbanSim, Inc.]↗

Mathematical model of the metal mould surface temperature optimization

The article is focused on the problem of generating a uniform temperature field on the inner surface of shell metal moulds. Such moulds are used e.g. in the automotive industry for artificial leather production. To produce artificial leather with uniform surface structure and colour shade the temperature on the inner surface of the mould has to be as homogeneous as possible. The heating of the mould is realized by infrared heaters located above the outer mould surface. The conceived mathematical model allows us to optimize the locations of infrared heaters over the mould, so that approximately uniform heat radiation intensity is generated. A version of differential evolution algorithm programmed in Matlab development environment was created by the authors for the optimization process. For temperate calculations software system ANSYS was used. A practical example of optimization of heaters locations and calculation of the temperature of the mould is included at the end of the article.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measuring Thread Timing to Assess the Feasibility of Early-Bird Message Delivery Across Systems and Scales

Early-bird communication is a communication/computation overlap technique that leverages fine-grained communication to improve application run-time. Communication is divided such that each individual thread can initiate transmission of its portion of the data upon completion rather than waiting for a dedicated communication phase. The benefit of early-bird communication depends on the completion timing of the individual threads: On the one hand, if all threads are complete at nearly the same time, the overheads of sending multiple messages will accumulate, leading to performance that is worse than if a single message had been sent. On the other hand, if thread completions are spread out in time, those that complete earlier can send data while others continue working, leading to performance that is better than if a single message had been sent. The challenge is that the completion times are currently unknown and can vary based on application, problem size, system software, and underlying hardware. In this paper, we address this lacuna by measuring and evaluating the potential overlap afforded by early-bird communication for a selection of proxy applications. These measurements help us understand whether a given application could benefit from early-bird communication. Here, we present our technique for gathering this data and evaluate data collected from three proxy applications: MiniFE, MiniMD, and MiniQMC. Each application is run on three systems with distinct CPU architectures and strong scales across three run sizes. To characterize the behavior of these workloads, we study the trends of thread timings at both a macro level, across all threads across all runs of an application, and a micro level, that is, within a single process of a single run. We observe that our tested applications exhibit significantly different thread arrival distributions. The machine used had a significant impact, with the window of potential overlap varying by as much as an order of magnitude.

97 MATHEMATICS AND COMPUTING↗

Spectral Effects in Albedo and Rearside Irradiance Measurement for Bifacial Performance Estimation

We investigate the impact of spectral dependence of ground surface reflectivity on albedo and rearside irradiance measurements necessary for bifacial photovoltaic (PV) module performance estimation and monitoring. Because PV modules are spectrally selective, albedo and irradiance measurements performed with common irradiance sensors may require spectral mismatch corrections when used for performance prediction. We investigate via simulation the differences in spectrally responsive albedo measured with thermopile pyranometers and crystalline silicon PV reference cells in comparison to a typical crystalline-silicon bifacial PV module. Simulations are performed for nine different representative ground surface materials using simulated solar spectra together with spectral reflectivity data distributed with the SMARTS simulation software. For the materials considered, the results show that albedo spectral mismatch relative to the bifacial module is distributed over a range of ±9.2% for thermopile pyranometers versus only ±3.7% for a typical PV reference cell. We consider the impact of this spectrally-responsive albedo mismatch on bifacial PV module rearside irradiance measurements. Using synthesized rearside spectral irradiance distributions, we find that for the nine different ground surface materials the predicted rearside irradiance measurement deviates from the effective irradiance observed by the PV module by on the order of 16.5 W/m 2 for the pyranometer and 3.6 W/m 2 for the PV reference cell. We discuss the implications for bifacial albedo and irradiance measurement.

41 EE - Solar Energy Technologies Office (EE-4S)↗

First evaluation of the feasibility of MLC tracking using ultrasound motion estimation

Purpose: To quantify the performance of the Clarity ultrasound (US) imaging system (Elekta AB, Stockholm, Sweden) for real-time dynamic multileaf collimator (MLC) tracking. Methods: The Clarity calibration and quality assurance phantom was mounted on a motion platform moving with a periodic sine wave trajectory. The detected position of a 30 mm hypoechogenic sphere within the phantom was continuously reported via Clarity’s real-time streaming interface to an in-house tracking and delivery software and subsequently used to adapt the MLC aperture. A portal imager measured MV treatment field/MLC apertures and motion platform positions throughout each experiment to independently quantify system latency and geometric error. Based on the measured range of latency values, a prostate stereotactic body radiation therapy (SBRT) delivery was performed with three realistic motion trajectories. The dosimetric impact of system latency on MLC tracking was directly measured using a 3D dosimeter mounted on the motion platform. Results: For 2D US imaging, the overall system latency, including all delay times from the imaging and delivery chain, ranged from 392 to 424 ms depending on the lateral sector size. For 3D US imaging, the latency ranged from 566 to 1031 ms depending on the elevational sweep. The latency-corrected geometric root-mean squared error was below 0.75 mm (2D US) and below 1.75 mm (3D US). For the prostate SBRT delivery, the impact of a range of system latencies (400–1000 ms) on the MLC tracking performance was minimal in terms of gamma failure rate. Conclusions: Real-time MLC tracking based on a noninvasive US input is technologically feasible. Current system latencies are higher than those for x-ray imaging systems, but US can provide full volumetric image data and the impact of system latency was measured to be small for a prostate SBRT case when using a US-like motion input.

60 APPLIED LIFE SCIENCES↗

Peregrine Software Development: Report on the Code Conversion From Python to C++

This work package seeks to convert the Peregrine software tool from its original Python implementation to a production version based on the C++ language. Peregrine is a powerful research platform with a multitude of advanced data analytics and data visualization functionalities. Developed by scientists to explore multimodal and multidimensional data related to the production of components using powder bed additive manufacturing processes, the tool implements state-of-the-art algorithms to assist machine users in making build or part quality determinations. Given that Peregrine is data-intensive, the goal of this conversion is to enhance the tool’s flexibility and interactivity and reduce the number of code dependencies to facilitate its deployment as part of the ongoing technology transfer campaign. This brief document provides an overview of Peregrine’s functionalities and capabilities, along with a detailed description of the core functionalities that have been implemented to date in the new C++ version. This document serves as a development update at the end of the first year of the ongoing conversion and will be regularly updated as progress continues.

97 MATHEMATICS AND COMPUTING↗