SEARCH · Search NASA

Results for “Code Robustness”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Liquid propellant rocket engine combustion simulation with a time-accurate CFD method

Time-accurate computational fluid dynamics (CFD) algorithms are among the basic requirements as an engineering or research tool for realistic simulations of transient combustion phenomena, such as combustion instability, transient start-up, etc., inside the rocket engine combustion chamber. A time-accurate pressure based method is employed in the FDNS code for combustion model development. This is in connection with other program development activities such as spray combustion model development and efficient finite-rate chemistry solution method implementation. In the present study, a second-order time-accurate time-marching scheme is employed. For better spatial resolutions near discontinuities (e.g., shocks, contact discontinuities), a 3rd-order accurate TVD scheme for modeling the convection terms is implemented in the FDNS code. Necessary modification to the predictor/multi-corrector solution algorithm in order to maintain time-accurate wave propagation is also investigated. Benchmark 1-D and multidimensional test cases, which include the classical shock tube wave propagation problems, resonant pipe test case, unsteady flow development of a blast tube test case, and H2/O2 rocket engine chamber combustion start-up transient simulation, etc., are investigated to validate and demonstrate the accuracy and robustness of the present numerical scheme and solution algorithm.

Chen, Y. S.↗

NASA NTRS2018

Initial Results of the Software-Driven Navigation for Station Experiment

To enable the next generation of robotic and human exploration of the solar system, improvements are needed to enable robust and accurate autonomous navigation. The purpose of this work is to take advantage of the growth in and use of software-defined platforms to incorporate additional navigation capability on existing assets, while also incorporating with new vehicle designs. The Software-driven Navigation for Station Experiment focuses on implementing two soft solutions to this: transmitting pseudolite signals to perform ranging and Doppler measurements as part of the signal coding (similar to underlying Global Navigation Satellite System approaches), and the Multi-spacecraft Autonomous Positioning System, which uses existing communication protocols to embed navigation and timing information to be shared among all assets in a peer-to-peer network. These technologies were implemented on the SCaN Testbed onboard the International Space Station and exercised over the course of mid-June and late-July 2018. This paper will discuss the operational architecture, experiment plan, and initial results from the data collected. One of the key conclusions of this work is the strong need for stable accurate clock synchronization across the dispersed space network.

Anzalone, Evan J.↗

NASA NTRSDate not supplied

Performance Assessment of LunaNet’s Augmented Forward Signal

LunaNet provides a common set of interoperable specifications for communication and position, navigation and time (PNT) services and interfaces soon to be implemented in lunar vicinity. The LunaNet Interoperability Specification (LNIS) provides the design for the GNSS-like Augmented Forward Signal (AFS), which enables orbiting and surface users in lunar space, such as Artemis, to estimate their position, velocity, and time. The specification of AFS defines two orthogonal signal components on a single carrier: the in-phase component (AFS-I), a lower-chip-rate data channel tailored for applications where low SWaP (Size, Weight, and Power) is critical (e.g., IoT devices or search and rescue), and the quadrature component (AFS-Q), a high-chip-rate data-less pilot signal for high-precision, robust lunar navigation and positioning applications. An initial description of AFS was provided in LNIS 2023, with initial analysis results shown in Dafesh 2024 and Dafesh 2025, and the current signal in space description provided in LNIS 2025. As part of NASA's Lunar Communication Relay and Navigation Systems (LCRNS) project, this work expands upon the initial analysis results and proposes a new expanded set of AFS-Q spreading codes that exceed the cross-correlation and autocorrelation sidelobe performance of L1C and other GNSS signals, while providing additional expansion capabilities for future service satellites. A set of 420 codes was selected from a Weil-based code derived from the prime number 10247, which is larger than the 10243 prime number used to derive BeiDou’s B1C Weil sequences. Both the initial set of 210 codes and the expanded set of 420 codes are shown to provide the best cross-correlation of any 10230-chip satellite navigation codes. The performance is demonstrated for hierarchical sets of spreading codes optimized and organized in sets of 30 codes. The work also compares LunaNet’s AFS to terrestrial GNSS signals in terms of acquisition, tracking, and data demodulation performance. Performance is evaluated for receivers that only track the 1.023 MCPS data channel spreading code for low SWaP IoT use cases, as well as for receivers that track both the 1.023 MCPS data channel and the 5.115 MCPS pilot channel spreading code for high-performance use cases. Performance is assessed in the presence of interference and thermal noise. The analysis is performed in terms of expected operating conditions on the lunar surface. Several unique flexibility aspects of the augmented forward signal are described, including the use of the Q channel’s secondary and tertiary codes to enable variable coherent integrations during acquisition. This is compared to GNSS signals such as L5/E5 and MBOC in terms of achievable processing gain for interference mitigation versus acquisition complexity. The work details acquisition and tracking techniques used to optimally acquire and track the primary, secondary, and tertiary codes on the Q channel, as well as acquisition of the I channel spreading code. Acquisition of the 8 ms, Q channel spreading code is also compared to joint acquisition of the I and Q channel primary codes in noise and interference environments

LCRNS↗

DOE OSTI2021

TRITON: A Multi-GPU open source 2D hydrodynamic flood model

A new open source multi-GPU 2D flood model called TRITON is presented in this work. The model solves the 2D shallow water equations with source terms using a time-explicit first order upwind scheme based on an Augmented Roe's solver that incorporates a careful estimation of bed strengths and a local implicit formulation of friction terms. Here, the scheme is demonstrated to be first order accurate, robust and able to solve for flows under various conditions. TRITON is implemented such that the model effectively utilizes heterogeneous architectures, from single to multiple CPUs and GPUs. Different test cases are shown to illustrate the capabilities and performance of the model, showing promising runtimes for large spatial and temporal scales when leveraging the computer power of GPUs. Under this hardware configuration, communication and input/output subroutines may impact the scalability. The code is developed under an open source license and can be freely downloaded in https://code.ornl.gov/hydro/triton.

2D flood model↗

DOE OSTI2017

User Guidelines and Best Practices for CASL VUQ Analysis Using Dakota

Sandia’s Dakota software (available at http://dakota.sandia.gov) supports science and engineering transformation through advanced exploration of simulations. Specifically it manages and analyzes ensembles of simulations to provide broader and deeper perspective for analysts and decision makers. This enables them to enhance understanding of risk, improve products, and assess simulation credibility. In its simplest mode, Dakota can automate typical parameter variation studies through a generic interface to a physics-based computational model. This can lend efficiency and rigor to manual parameter perturbation studies already being conducted by analysts. However, Dakota also delivers advanced parametric analysis techniques enabling design exploration, optimization, model calibration, risk analysis, and quantification of margins and uncertainty with such models. It directly supports verification and validation activities. Dakota algorithms enrich complex science and engineering models, enabling an analyst to answer crucial questions of - Sensitivity: Which are the most important input factors or parameters entering the simulation, and how do they influence key outputs?; Uncertainty: What is the uncertainty or variability in simulation output, given uncertainties in input parameters? How safe, reliable, robust, or variable is my system? (Quantification of margins and uncertainty, QMU); Optimization: What parameter values yield the best performing design or operating condition, given constraints? Calibration: What models and/or parameters best match experimental data? In general, Dakota is the Consortium for Advanced Simulation of Light Water Reactors (CASL) delivery vehicle for verification, validation, and uncertainty quantification (VUQ) algorithms. It permits ready application of the VUQ methods described above to simulation codes by CASL researchers, code developers, and application engineers.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1988

Aerodynamic applications of an efficient incompressible Navier-Stokes solver

A robust and accurate upwind scheme based on flux-difference splitting for three-dimensional, incompressible flows in general coordinates is described. Viscous effects are included through a thin-layer approximation to the Navier-Stokes equations. Turbulent flows are simulated with extensions to the Baldwin-Lomax (1978) turbulence model. The effects of accuracy and gridding on the computed flow field results are assessed. The results obtained from the code for vortical flows over a delta wing, a double delta wing, and a tangent-ogive forebody compare well with experimental data.

Hartwich, Peter M.↗

NASA NTRS2019

Weighted Least-Squares Cell-Average Gradient Construction Methods for the VULCAN-CFD Second-Order Accurate Unstructured-Grid Cell-Centered Finite-Volume Solver

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. The construction of cell-average gradients using a weighted linear least-squares method and the use of these gradients in the construction of the inviscid fluxes is the focus of this paper. A comparison of least-squares stencil construction methodologies is presented and approaches to augment the number of cells participating in the stencil while preserving accuracy are explored. Due to our interest in hypersonic flow, a robust multidimensional cell-average gradient limiter procedure that is consistent with the stencil used to construct the cell-average gradients is described and investigated. Canonical problems are computed to illustrate the challenges and investigate the accuracy, robustness and convergence behavior of the cell-average gradient methods on unstructured cell-centered finite-volume grids. Finally, thermally perfect, chemically frozen, Mach 8 turbulent flow of air around a blunt wedge is computed to demonstrate the robustness and convergence behavior of the new method for constructing stencils of use in a weighted linear least-squares gradient method for a hypersonic flow.

White, Jeffery A.↗

DOE OSTI2018

TENDL-2017: the making of multi-faceted technological nuclear data library - 25340

Multi-purpose methodical nuclear data libraries generated by TALYS/T6 as Evaluated Nuclear Data Library TENDL have been released regularly since 2008. Considerable experience has been acquired during those early production of such unique, truly multi-faceted, multi-particle nuclear data libraries based on the feedback from the developers, evaluators, processing experts and, most importantly, users. The backbone of this achievement is simple and robust: completeness, quality, upgradability and, most of all, reproducibility and methodology. Since TENDL has been seamlessly adopted by many different applications (accelerator, astrophysics, fusion, fission, medical) that require multi-faceted nuclear reaction data forms for shielding, radio-protection, transmutation, materials or earth sciences, it is felt necessary to understand its strengths and remaining weaknesses. The essential knowledge is not the TENDL library itself, but rather the necessary methods, processes, codes, tools and know-how that go into the making of every evaluations of such a library. Current efforts are focused on a proper evaluation of the underlying physics and incorporation of information and metrics into the scientific TALYS/T6 system. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2018

Wide ranging equation of state with Tartarus: A hybrid Green’s function/orbital based average atom code

Average atom models are widely used to make equation of state tables and for calculating other properties of materials over a wide range of conditions, from zero temperature isolated atom to fully ionized free electron gases. The numerical challenge of making these density functional theory based models work for any temperature, density or nuclear species is formidable. Here we present in detail a hybrid Green’s function/orbital based approach that has proved to be stable and accurate for wide ranging conditions. Algorithmic strategies are discussed. In particular the decomposition of the electron density into numerically advantageous parts is presented and a robust and rapid self consistent field method based on a quasi-Newton algorithm is given. Example application to the equation of state of lutetium (Z = 71) is explored in detail, including the effect of relativity, finite temperature exchange and correlation, and a comparison to a less approximate method. Furthermore, the hybrid scheme is found to be numerically stable and accurate for lutetium over at least 6 orders of magnitude in density and 5 orders of magnitude in temperature.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2026

Documentation, assessment, and improvements of fission gas modeling capabilities in BISON

Fission gas behavior critically affects nuclear fuel performance, influencing thermal conductivity, microstructure, swelling, and fuel-cladding interactions. The Simple Integrated Fission Gas Release and Swelling (Sifgrs) model, implemented within the BISON fuel performance code, provides advanced capabilities for modeling these phenomena. This manuscript reviews and assesses Sifgrs’ modular structure and submodels describing intragranular, intergranular, and release mechanisms. The modular implementation enables rapid improvements as our understanding advances and new fuel forms emerge. After an overview of capabilities, we present a comprehensive validation suite, comparing predictions to experimental data from tests under normal and transient conditions, including reproduction of the Vitanza curve. Results demonstrate Sifgrs’ accuracy, robustness, and applicability. Current limitations and ongoing efforts are discussed, including mechanistic, multiscale modeling, uncertainty quantification, and validation within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, supporting the safe and efficient use of nuclear fuels.

BISON↗

DOE OSTI2020

MHD stability and disruptions in the SPARC tokamak

SPARC is being designed to operate with a normalized beta of $β_N$ = 1.0, a normalized density of $n_G$ = 0.37 and a safety factor of $q_{95}$ ≈ 3.4, providing a comfortable margin to their respective disruption limits. Further, a low beta poloidal $β_p$ = 0.19 at the safety factor $\textit{q}$ = 2 surface reduces the drive for neoclassical tearing modes, which together with a frozen-in classically stable current profile might allow access to a robustly tearing-free operating space. Although the inherent stability is expected to reduce the frequency of disruptions, the disruption loading is comparable to and in some cases higher than that of ITER. The machine is being designed to withstand the predicted unmitigated axisymmetric halo current forces up to 50 MN and similarly large loads from eddy currents forced to flow poloidally in the vacuum vessel. Runaway electron (RE) simulations using GO+CODE show high flattop-to-RE current conversions in the absence of seed losses, although NIMROD modelling predicts losses of ~80 %; self-consistent modelling is ongoing. A passive RE mitigation coil designed to drive stochastic RE losses is being considered and COMSOL modelling predicts peak normalized fields at the plasma of order 10 –2 that rises linearly with a change in the plasma current. Massive material injection is planned to reduce the disruption loading. A data-driven approach to predict an oncoming disruption and trigger mitigation is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Unveiling the chemical fingerprint of phosphorus-rich stars: I. In the infrared region of APOGEE-2

The origin of phosphorus, one of the essential elements for life on Earth, is currently unknown. Prevalent models of Galactic chemical evolution (GCE) systematically underestimate the amount of P compared to observations, especially at low metallicities. The recently discovered P-rich ([P/Fe] ≳ 1.2 dex) and metal-poor ([Fe/H] ≃ –1.0 dex) giants further challenge the GCE models, calling current theories on stellar nucleosynthesis into question. Since the observed low-mass giants are not expected to produce their high P contents themselves, our primary goal is to find clues on their progenitor or polluter. By increasing the number of known P-rich stars, we aim to narrow down a statistically reliable chemical abundance pattern that defines these peculiar stars. In this way, we place more robust constraints on the nucleosynthetic mechanism that causes the unusually high P abundances. In the long term, identifying the progenitor of the P-rich stars may contribute to the search for the source of P in our Galaxy. We performed a detailed chemical abundance analysis based on the high-resolution near-infrared (H band) spectra from the latest data release (DR17) of the APOGEE-2 survey. Employing the BACCHUS code, we measured the abundances of 13 elements in the inspected sample, which is mainly composed of a recent collection of Si-enhanced giants. We also analyzed the orbital motions and compared the abundance results to possible nucleosynthetic formation scenarios, and also to detailed GCE models. These models were produced with the OMEGA+ chemical evolution code, using four different massive star yield sets to investigate different scenarios for massive star evolution. We enlarged the sample of confirmed P-rich stars from 16 to a group of 78 giants, which represents the largest sample of P-rich stars to date. The sample includes the first detection of a P-rich star in a Galactic globular cluster. Significant enhancements in O, Al, Si, and Ce, as well as systematic correlations among the studied elements, unveil the unique chemical fingerprint of the P-rich stars. In contrast, the high [Mg/Fe] and [(C+N)/Fe] found in some of the P-rich stars with respect to P-normal stars is not confirmed over the full sample because of the current uncertainties. Strikingly, the strong overabundance in the α-element Si is accompanied by normal Ca and S abundances. This is at odds with current stellar nucleosynthesis models of massive stars. Our analysis of the orbital motion showed that the P-rich stars do not belong to a locally specific population in the Galaxy. In addition, we confirm that the majority of the sample stars are not part of binary systems.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

MAESTROeX: A Massively Parallel Low Mach Number Astrophysical Solver

@e present MAESTROeX, a massively parallel solver for low Mach number astrophysical flows. The underlying low Mach number equation set allows for efficient, long-time integration for highly subsonic flows compared to compressible approaches. MAESTROeX is suitable for modeling full spherical stars as well as well as planar simulations of dynamics within localized regions of a star, and can robustly handle several orders of magnitude of density and pressure stratification. Previously, we have described the development of the predecessor of MAESTROeX, called MAESTRO, in a series of papers. Here, we present a new, greatly simplified temporal integration scheme that retains the same order of accuracy as our previous approaches. We also explore the use of alternative spatial mapping of the one-dimensional base state onto the full Cartesian grid. The code leverages the new AMReX software framework for block-structured adaptive mesh refinement (AMR) applications, allowing for scalability to large fractions of leadership-class machines. Finally, using our previous studies on the convective phase of single-degenerate progenitor models of SNe Ia as a guide, we characterize the performance of the code and validate the new algorithmic features. Like MAESTRO, MAESTROeX is fully open source.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1995

The coupling of fluids, dynamics, and controls on advanced architecture computers

This grant provided for the demonstration of coupled controls, body dynamics, and fluids computations in a workstation cluster environment; and an investigation of the impact of peer-peer communication on flow solver performance and robustness. The findings of these investigations were documented in the conference articles.The attached publication, 'Towards Distributed Fluids/Controls Simulations', documents the solution and scaling of the coupled Navier-Stokes, Euler rigid-body dynamics, and state feedback control equations for a two-dimensional canard-wing. The poor scaling shown was due to serialized grid connectivity computation and Ethernet bandwidth limits. The scaling of a peer-to-peer communication flow code on an IBM SP-2 was also shown. The scaling of the code on the switched fabric-linked nodes was good, with a 2.4 percent loss due to communication of intergrid boundary point information. The code performance on 30 worker nodes was 1.7 (mu)s/point/iteration, or a factor of three over a Cray C-90 head. The attached paper, 'Nonlinear Fluid Computations in a Distributed Environment', documents the effect of several computational rate enhancing methods on convergence. For the cases shown, the highest throughput was achieved using boundary updates at each step, with the manager process performing communication tasks only. Constrained domain decomposition of the implicit fluid equations did not degrade the convergence rate or final solution. The scaling of a coupled body/fluid dynamics problem on an Ethernet-linked cluster was also shown.

Atwood, Christopher↗

NASA NTRS2008

A Numerical Investigation of Turbine Noise Source Hierarchy and Its Acoustic Transmission Characteristics

Understanding the relative importance of the various turbine noise generation mechanisms and the characteristics of the turbine acoustic transmission loss are essential ingredients in developing robust reduced-order models for predicting the turbine noise signature. A computationally based investigation has been undertaken to help guide the development of a turbine noise prediction capability that does not rely on empiricism. The investigation relies on highly detailed numerical simulations of the unsteady flowfield inside a modern high-pressure turbine (HPT). The simulations are developed using TURBO, which is an unsteady Reynolds-averaged Navier-Stokes (URANS) code capable of multi-stage simulations. The purpose of this study is twofold. First, to determine an estimate of the relative importance of the contributions to the coherent part of the acoustic signature of a turbine from the three potential sources of turbine noise generation, namely, blade-row viscous interaction, potential field interaction, and entropic source associated with the interaction of the blade rows with the temperature nonuniformities caused by the incomplete mixing of the hot fluid and the cooling flow. Second, to develop an understanding of the turbine acoustic transmission characteristics and to assess the applicability of existing empirical and analytical transmission loss models to realistic geometries and flow conditions for modern turbine designs. The investigation so far has concentrated on two simulations: (1) a single-stage HPT and (2) a two-stage HPT and the associated inter-turbine duct/strut segment. The simulations are designed to resolve up to the second harmonic of the blade passing frequency tone in accordance with accepted rules for second order solvers like TURBO. The calculations include blade and vane cooling flows and a radial profile of pressure and temperature at the turbine inlet. The calculation can be modified later to include the combustor pattern factor at the turbine inlet to include that contribution to turbine noise. We shall present preliminary analysis of the results obtained so far in order to assess the validity of such an approach and to seek feedback on improving the approach. This work addresses both Area 1 (Turbine Tone Noise) and Area 5 (Influence of the Turbine on Combustor Noise) topics.

VanZante, Dale↗

NASA NTRS2012

Open Rotor Aeroacoustic Modelling

Owing to their inherent fuel efficiency, there is renewed interest in developing open rotor propulsion systems that are both efficient and quiet. The major contributor to the overall noise of an open rotor system is the propulsor noise, which is produced as a result of the interaction of the airstream with the counter-rotating blades. As such, robust aeroacoustic prediction methods are an essential ingredient in any approach to designing low-noise open rotor systems. To that end, an effort has been underway at NASA to assess current open rotor noise prediction tools and develop new capabilities. Under this effort, high-fidelity aerodynamic simulations of a benchmark open rotor blade set were carried out and used to make noise predictions via existing NASA open rotor noise prediction codes. The results have been compared with the aerodynamic and acoustic data that were acquired for this benchmark open rotor blade set. The emphasis of this paper is on providing a summary of recent results from a NASA Glenn effort to validate an in-house open noise prediction code called LINPROP which is based on a high-blade-count asymptotic approximation to the Ffowcs-Williams Hawkings Equation. The results suggest that while predicting the absolute levels may be difficult, the noise trends are reasonably well predicted by this approach.

Envia, Edmane↗

NASA NTRS2012

Open Rotor Aeroacoustic Modeling

Envia, Edmane↗

DOE OSTI2016

Characterization of a Polyethylene Moderated Highly Enriched Uranium System

The polyethylene interstitial materials for the Class Foils experiment are one of the teaching tools used to teach criticality safety at the National Criticality Experiments Research Center (NCERC). This set of polyethylene plates is designed to provide the optimal H/{sup 235}U ratio on the critical mass curve allowing the experiment to have the minimum amount of {sup 235}U. This makes a great teaching tool to show that either additional {sup 235}U or additional moderator would increase the critical mass. The optimal H/{sup 235}U ratio and increased hydrogen density compared to the previously used Lucite plates has reduced the critical mass by 34%, to 997.9 ± 0.65 g. The robustness of this experiment makes it ideal for hands on training. The fuel foils consist of thin highly enriched uranium (93 wt% {sup 235}U) metal foils laminated in plastic; each foil weighs approximately 68 grams. The fuel is interleaved between interstitial plates to attain a critical configuration. The critical configuration requires 14.5 fuel foils. Rossi-α experiments have been proposed as a method of benchmarking kinetics properties given by neutronics codes such as MCNP. (authors)

07 ISOTOPE AND RADIATION SOURCES↗