Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coding”

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,153 records · Page 64

Further Validation of a CFD Code for Calculating the Performance of Two-Stage Light Gas Guns

Earlier validations of a higher-order Godunov code for modeling the performance of two-stage light gas guns are reviewed. These validation comparisons were made between code predictions and experimental data from the NASA Ames 1.5" and 0.28" guns and covered muzzle velocities of 6.5 to 7.2 km/s. In the present report, five more series of code validation comparisons involving experimental data from the Ames 0.22" (1.28" pump tube diameter), 0.28", 0.50", 1.00" and 1.50" guns are presented. The total muzzle velocity range of the validation data presented herein is 3 to 11.3 km/s. The agreement between the experimental data and CFD results is judged to be very good. Muzzle velocities were predicted within 0.35 km/s for 74% of the cases studied with maximum differences being 0.5 km/s and for 4 out of 50 cases, 0.5 - 0.7 km/s.

Ballistic Range↗

Performance of the Dot Product Function in Radiative Transfer Code SORD

The successive orders of scattering radiative transfer (RT) codes frequently call the scalar (dot) product function. In this paper, we study performance of some implementations of the dot product in the RT code SORD using 50 scenarios for light scattering in the atmosphere-surface system. In the dot product function, we use the unrolled loops technique with different unrolling factor. We also considered the intrinsic Fortran functions. We show results for two machines: ifort compiler under Windows, and pgf90 under Linux. Intrinsic DOT_PRODUCT function showed best performance for the ifort. For the pgf90, the dot product implemented with unrolling factor 4 was the fastest. The RT code SORD together with the interface that runs all the mentioned tests are publicly available from ftp:maiac.gsfc.nasa.govpubskorkinSORD_IP_16B (current release) or by email request from the corresponding (first) author.

polarized radiative transfer↗

Development and Use of the SPACE Computer Code for Analyzing the Space Station Electrical Power System

This special publication tells the story of the dedicated efforts of very talented individuals to create a preeminent space electrical power system modeling and simulation tool called SPACE, short for System Power Analysis for Capability Evaluation. This computer model has evolved for 30 years, withstanding the test of time and obsolescence, and garnering international recognition for its ability to simulate complex space electrical power systems. The analytical results from this model have saved millions of dollars in hardware redesign, testing, and verification for NASA's International Space Station (ISS) and its European and Russian partners. SPACE has played a pivotal role in the station's design and development and continues to support its ongoing operation. It has also extended its reach beyond the ISS to other key NASA programs, where it guides the design and planned operation of NASA's Multi-Purpose Crew Vehicle Orion and simulates electric power system operation in a dusty atmosphere on Mars' surface. The SPACE lineage was created by a core civil servant staff, supplemented by a cadre of interns and other temporary helpers. They created a tightly integrated tool that includes all phenomena that impact a solar array and battery space power system performance. SPACE is self-contained, requiring no other software modules and associated license fees. SPACE "rings true" in that is has been extensively validated with ISS on-orbit telemetry data. This report is being released as the generation of engineers who created it are nearing retirement, passing the baton to a new generation. This next generation will carry the code into the future, no doubt further evolving it to be able to assure mission planners that newly conceived systems will successfully power NASA's next endeavors. As a previous SPACE code developer and analyst, I have worked alongside many of the people mentioned in this report. The engineers who created the code, along with those just now learning it, are among the best and brightest at NASA. It is an honor to write this foreword as the present Branch Chief under which the legacy of SPACE continues to thrive.

Electrical Power System↗

Scaling Evaluation of Ice-Crystal Icing on a Modern Turbofan Engine in PSL Using the COMDES-MELT Code

This paper presents preliminary ice-crystal icing (ICI) altitude scaling evaluation results of a Honeywell Uncertified Research Engine (HURE) that was tested in the NASA Glenn Research Center Propulsion Systems Laboratory (PSL) during January of 2018. This engine geometry features a hidden core design to keep the core less exposed. The engine was fitted with internal video cameras to observe various ice buildup processes at multiple selected locations within the engine core flow path covering the fan stator, the splitter-lip/shroud/strut, and the high pressure compressor (HPC) variable inlet guide vane (IGV) regions. The potential ice accretion risk was pre-determined to occur by using NASA’s in-house 1D Engine Icing Risk assessment code, COMDES-MELT. The code was successful in predicting the risk of ice accretion in adiabatic regions like the fan-stator of the HURE at specific engine operating points. However at several operating points during the test, liquid water was observed running along the shroud toward the variable IGV of the HPC regions with an air temperature well below freezing, thus no particle melting could have occurred due to heating from the air alone. It was reasoned that other sources of heat were present in that region. To account for these heat sources the inlet total temperature was adjusted to give a wet bulb temperature of 24 F below the standard minimum wet bulb temperature of 492 R to allow ice to accrete in the splitter-lip/shroud/strut region, which was determined from a reference case where hard ice was observed in that region. With that adjustment the COMDES-MELT code was successful in providing operating points where there was a risk of ice accretion during the test campaign. In addition to calculating possible conditions at different selected lower altitudes, simulations were run to determine potential inlet conditions that could lead to ice-crystal accretion along the prescribed stations where the cameras were available. From there, scaled test conditions were determined by best matching the following three icing related parameters of the reference condition: (1) the local air total wet bulb temperature, (2) the local ice crystal cloud melt ratio and (3) the engine fan face ice/water to air mass flux ratio of the ice crystal cloud. Instantaneous images taken from the time-lapsed movies of ice buildup were used along with the relevant thermodynamic data of air, water vapor and local icing condition to help evaluate how closely the proposed altitude scaling method could be used in ground based test facility to duplicate selected reference ICI features observed at specific location inside this engine at different scale altitudes. Discussions on observed limitation for engine icing scaling application from this test campaign and needed improvement are provided. A scaling test procedure to help identify potential ICI risk conditions and possible ice accretion locations of a new turbofan engine is evaluated in PSL.

Engine Testing in PSL↗

CoCoSim, a Code Generation Framework for Control/command Applications: An Overview of CoCoSim for Multi-Periodic Discrete Simulink Models

We present CoCoSim, a framework to support the design, code generation and analysis of discrete dataflow model expressed in Simulink. In this work, we specifically focus on the analysis and code generation of multi-periodic systems. For that CoCoSim provides two complementary approaches: the first amounts to encode the multiperiodic semantics in a pure-synchronous one – à la Lustre–, enabling the use of model checker for verifying properties. The second provides a faithful code generation into multiple communicating (mono)synchronous components – à la Prelude– that can be then simulated or embedded in the final platform with any real-time scheduler. These approaches have been experimented in various settings.

Bourbouh, Hamza↗

Light Ion Double-Differential Cross Section Parameterization and Results from the SHIELD Transport Code

Light ions and neutrons have been shown to make large contributions to space radiation dose equivalent for realistic shielding scenarios. Efficient and accurate calculations of light ion double-differential cross sections are required for input into space radiation transport codes. A thermal proton cross section model is developed which includes proton production from the three sources of projectile, central fireball, and target. It is shown that this three-source model is able to explain the low momentum shoulder seen in proton spectra. Using the coalescence model, the thermal proton model is used to calculate light ion double-differential cross sections employed in space radiation transport codes. The three-source model is also seen to be essential to explain the light ion shoulders, which are even more pronounced than the proton shoulders. Comparisons are also made to the cross section models used in the SHIELD transport code.

Space radiation↗

Introduction and Overview of the Object-Oriented Turbomachinery Analysis Code (OTAC)

This is an overview of the Object-Oriented Turbomachinery Analysis Code (OTAC). OTAC is an object-oriented code developed to perform 2-D meridional flowthrough design and analysis of turbomachines. OTAC uses an implicit formulation of the governing equations relevant to meanline and streamline simulations to solve for the thermodynamic and velocity state of the fluid at the exit of each blade row. OTAC is written using the interpreted language capability of the Numerical Propulsion System Simulation (NPSS) code; as such, OTAC requires a version of NPSS to run. A series of models and test cases have been developed to demonstrate the capability and features of OTAC.

Turbomachinery↗

New Features of the NEQAIR Radiation Code

The longest-lived code for predicting shock layer radiation, NEQAIR, is now in its 5th decade of service. Substantial changes to the code have been made over the previous decade, the most recent report of which was at the 5th Workshop on Radiation in High Temperature Gases in 2014, for the version referred to as NEQAIR14. This paper will review some of the improvements made to the NEQAIR code since then, which is now at v15.2. Some of these features are discussed briefly below. NEQAIR15 and subsequent versions have enabled parallel evaluation of multiple lines of sight. This is accomplished by utilizing the HDF5 file format and placing multiple lines into a single file, LOS.h5, which is used for both input and output. This approach enables straightforward parallel execution both over the number of lines of sight and the number of points per line. For large problems, runtime reduces linearly with the number of nodes deployed since each line is processed independently by a subset of MPI ranks. Three applications of the multi-line solver are discussed. The first has to do with performing loosely coupled radiation-flowfield solutions. In this case the computed absorption and emission coefficients are used to evaluate the total energy absorbed or emitted at each point, allowing evaluation of the volumetric source term in the flowfield. The second computation is for obtaining heat flux from nonuniform flows, which require integration over spherical co-ordinates. These are of particular interest for evaluating radiation on the vehicle backshell. This 3D option improves the angular integration scheme and allows adaptive line selection that together reduce the number of lines required by about an order of magnitude. The final application is for remote observation, which is essentially the 3D integration problem over a small solid angle. For all three of these computations, data can be stored in the HDF5 file which allows a NEQAIR run to be restarted when it times out, or to add atmospheric absorption or instrument scan functions. An additional level of parallelism is enabled in NEQAIR15.2 using GPU routines. The GPU parallelism has realized up to 8x speed-up when running on a single core but diminishes as CPU parallelism is increased. For running multi-line simulations, it may be easier to reserve a large number of CPU nodes than to obtain the number of GPU nodes required for similar performance. A GUI, known as NEQTPY, allows for reading and creating input files, running NEQAIR, and displaying results. A significant feature of NEQTPY is the ability to perform spectral fits to data. The fits can operate on a single line spectrum (radiance vs. wavelength) or a 3D input file with multiple columns of data. Other new features include improved constants, additional species, more detailed non-Boltzmann modelling, advanced user controls, the ability to read and calculate spectra from HITRAN datafiles, photodissociation and photoionization cross-sections. A “fast” automatic grid option may reduce the size and time of spectral calculations while still maintaining good accuracy for total heat flux.

Brett A Cruden↗

Mercury-T: A New Code to Study Tidally Evolving Multi-Planet Systems. Applications to Kepler-62

A large proportion of observed planetary systems contain several planets in a compact orbital configuration, and often harbor at least one close-in object. These systems are then most likely tidally evolving. We investigate how the effects of planet-planet interactions influence the tidal evolution of planets. We introduce for that purpose a new open-source addition to the Mercury N-body code, Mercury-T, which takes into account tides, general relativity and the effect of rotation-induced flattening in order to simulate the dynamical and tidal evolution of multi-planet systems. It uses a standard equilibrium tidal model, the constant time lag model. Besides, the evolution of the radius of several host bodies has been implemented (brown dwarfs, M-dwarfs of mass 0:1 M(sub ⨀), Sun-like stars, Jupiter). We validate the new code by comparing its output for one-planet systems to the secular equations results. We find that this code does respect the conservation of total angular momentum. We applied this new tool to the planetary system Kepler-62. We find that tides influence the stability of the system in some cases. We also show that while the four inner planets of the systems are likely to have slow rotation rates and small obliquities, the fifth planet could have a fast rotation rate and a high obliquity. This means that the two habitable zone planets of this system, Kepler-62e ad f are likely to have very different climate features, and this of course would influence their potential at hosting surface liquid water..

Dynamical evolution and stability↗

GPU Implementation of the OVERFLOW CFD Code

The high-performance computing (HPC) landscape is quickly changing to systems where most of the performance comes from specialized chips, specifically graphics processing units (GPUs). Such GPU systems are throughput machines, where efficient use of the GPU often requires code refactoring to expose a few orders of magnitude more fine grain parallelism than was previously used on the CPU. Recent modifications to OVERFLOW, an overset, structured grid, computational fluid dynamics flow solver, written in Fortran will be presented. These modifications include both code modernization efforts and algorithmic changes to enable OVERFLOW to efficiently utilize GPUs. Many of these algorithmic changes would likely also be applicable for other structured grid, stencil-based codes wanting to utilize GPUs. The capabilities that have been ported to run on the GPUs are presented, along with the performance gains of the GPU version relative the CPU version of OVERFLOW.

GPU Programming↗

GPU Implementation of the OVERFLOW CFD Code

The high-performance computing (HPC) landscape is quickly changing to systems where most of the performance comes from specialized chips, specifically graphics processing units (GPUs). Such GPU systems are throughput machines, where efficient use of the GPU often requires code refactoring to expose a few orders of magnitude more fine grain parallelism than was previously used on the CPU. Recent modifications to OVERFLOW, an overset, structured grid, computational fluid dynamics flow solver, written in Fortran will be presented. These modifications include both code modernization efforts and algorithmic changes to enable OVERFLOW to efficiently utilize GPUs. Many of these algorithmic changes would likely also be applicable for other structured grid, stencil-based codes wanting to utilize GPUs. The capabilities that have been ported to run on the GPUs are presented, along with the performance gains of the GPU version relative the CPU version of OVERFLOW.

GPU Programming↗

Code Migration and Optimization for Calculating Atmospheric Attenuation in RF Link Budgets Within Glenn Research Center Communication Analysis Suite

This technical report encompasses development completed on the Glenn Research Center Communication Analysis Suite (GCAS) as it relates to converting code from the original MATLAB® (The MathWorks, Inc.) GCAS code base to C++. It focuses on a majority of the functions, structures, and intermediate code developed for radiofrequency (RF) links calculations. The report serves as a comprehensive record of the current progress and outlines the technical challenges overcome during this research. Throughout this project, numerous functions and structures were migrated end to end, with validation test cases conducted in accordance with NASA’s verification and validation efforts for standard models and simulations. This report also provides high-level documentation of the completed work; detailed low-level documentation has been compiled to provide a line-by-line explanation in alignment with NASA’s goals. This development has largely been focused on atmospheric effects on calculating RF link budgets for both terrestrial and space links. Some developments outside of atmospheric functions are also covered throughout this report; these are largely related to mathematical operations and MATLAB® structure conversions.

Brady M. Phelps↗

Updates and Modernization of NASA’s Chemical Equilibrium with Applications (CEA) Code

NASA’s Chemical Equilibrium with Applications (CEA) code is a foundational tool for propulsion system analysis. It provides equilibrium chemistry, rocket performance, shock, and detonation calculations used across NASA and the broader aerospace community. NASA Engineering and Safety Center (NESC) Activity TI-22-01730 modernized the legacy CEA2 Fortran code into CEA v3, a Fortran 2008, object-oriented software package with expanded interface support, updated thermochemical data, improved maintainability, and substantially improved workflow integration. The modernized code preserves backward compatibility with legacy CEA input workflows while enabling direct use from modern analysis environments, including Python, C, MATLAB, and automated design studies.

Mark K Leader↗

Updates and Modernization of the Chemical Equilibrium with Applications (CEA) Code

NASA’s Chemical Equilibrium with Applications (CEA) code is a foundational tool for propulsion system analysis. It provides equilibrium chemistry, rocket performance, shock, and detonation calculations used across NASA and the broader aerospace community. NASA Engineering and Safety Center (NESC) Activity TI-22-01730 modernized the legacy CEA2 Fortran code into CEA v3, a Fortran 2008, object-oriented software package with expanded interface support, updated thermochemical data, improved maintainability, and substantially improved workflow integration. The modernized code preserves backward compatibility with legacy CEA input workflows while enabling direct use from modern analysis environments, including Python, C, MATLAB, and automated design studies.

Combustion↗

New codes for the optical channel

We present new codes and modulation formats for the deep space optical channel. By taking maximum avdvantage of inherent physical constraints on the communications link, these codes and modulation formats substantially close the large gap between state of the art systems and the Shannon limit.

channel codes modulation optical channel turbo cod↗

Optimization of Particle-in-Cell Codes on RISC Processors

General strategies are developed to optimize particle-cell-codes written in Fortran for RISC processors which are commonly used on massively parallel computers. These strategies include data reorganization to improve cache utilization and code reorganization to improve efficiency of arithmetic pipelines.

particle-cell-codes cache utilization code reorgan↗