Search NASA⌕ Search

SEARCH · Search NASA

Results for “computer code development”

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 505 records · Page 28

On the Coupling Between a Supersonic Turbulent Boundary Layer and a Flexible Structure

A mathematical model and a computer code have been developed to fully couple the vibration of an aircraft fuselage panel to the surrounding flow field, turbulent boundary layer and acoustic fluid. The turbulent boundary layer model is derived using a triple decomposition of the flow variables and applying a conditional averaging to the resulting equations. Linearized panel and acoustic equations are used. Results from this model are in good agreement with existing experimental and numerical data. It is shown that in the supersonic regime, full coupling of the flexible panel leads to lower response and radiation from the panel. This is believed to be due to an increase in acoustic damping on the panel in this regime. Increasing the Mach number increases the acoustic damping, which is in agreement with earlier work.

Frendi, Abdelkader↗

An Investigation of High-Order Shock-Capturing Methods for Computational Aeroacoustics

This project is motivated by the desire to develop numerical methods that will be useful in the study of compressible flows that exhibit aeroacoustic phenomena. Solutions to linear problems have been investigated through the development of a computer code based on the recent dispersion-relation-preserving (DRP) methodology. In regard to nonlinear problems, the class of essentially nonoscillatory (ENO) schemes have been considered as the primary candidates for solving aeroacoustic problems in which discontinuities are involved. Discontinuities in the solution itself (e.g. shocks) as well as in the geometry on which the problem is defined have been studied. Two-dimensional nonlinear problems were considered in order to determine if the one-dimensional results obtained in the first phase of this project were extendable to more realistic problems. Conclusions have been drawn in regard to the ability to numerically predict solutions of nonlinear problems with shocks to high-order accuracy.

Casper, Jay↗

Computational Modeling Program

An Integrated Product Team (IPT) has been formed at NASA Ames Research Center which has set objectives to investigate devices and processes suitable for meeting NASA requirements on ultrahigh performance computers, fast and low power devices, and high temperature wide bandgap materials. These devices may ultimately be sub-100nm feature-size. Processes and equipment must meet the stringent demands posed by the fabrication of such small devices. Until now, the reactors for Chemical Vapor Deposition (CVD) and plasma processes have been designed by trial and error procedures. Further, once the reactor is in place, optimum processing parameters are found through expensive and time-consuming experimentation. If reliable models are available that describe processes and the operation of the reactors, that chore would be reduced to a routine task while being a cost-effective option. The goal is to develop such a design tool, validate that tool using available data from current generation processes and reactors, and then use that tool to explore avenues for meeting NASA needs for ultrasmall device fabrication. Under the present grant, ARL/Penn State along with other IPT members has been developing models and computer code to meet IPT goals. Some of the accomplishments achieved during the first year of the grant are described in this report

Govindan, T. R.↗

Calculation of Nonequilibrium Radiation in the Flight Regimes of Aeroassisted Orbital Transfer Vehicles

A computer code has been developed that calculates radiative properties of nonequilibrium air in the low-density regimes expected during the flight of aeroassisted, orbital transfer vehicles. From the given nonequilibrium thermodynamic state variables, the code calculates number densities of internal states and the accompanying emission and absorption characteristics. In addition, the code calculates the number density of the hypothetical gas in radiative equilibrium that produces the same radiation emission as the given nonequilibrium gas. Sample results are shown to demonstrate how the code predicts the non-linear variation of radiation with density at low densities.

Park, Chul↗

Significance of Strain in Formulation in Theory of Solid Mechanics

The basic theory of solid mechanics was deemed complete circa 1860 when St. Venant provided the strain formulation or the field compatibility condition. The strain formulation was incomplete. The missing portion has been formulated and identified as the boundary compatibility condition (BCC). The BCC, derived through a variational formulation, has been verified through integral theorem and solution of problems. The BCC, unlike the field counterpart, do not trivialize when expressed in displacements. Navier s method and the stiffness formulation have to account for the extra conditions especially at the inter-element boundaries in a finite element model. Completion of the strain formulation has led to the revival of the direct force calculation methods: the Integrated Force Method (IFM) and its dual (IFMD) for finite element analysis, and the completed Beltrami-Michell formulation (CBMF) in elasticity. The benefits from the new methods in elasticity, in finite element analysis, and in design optimization are discussed. Existing solutions and computer codes may have to be adjusted for the compliance of the new conditions. Complacency because the discipline is over a century old and computer codes have been developed for half a century can lead to stagnation of the discipline.

Patnaik, Surya N.↗

Studies of Pressure-Broadening of Alkali Atom Resonance Lines for Modeling Atmospheres of Extrasolar Giant Planets and Brown Dwarfs

In L-dwarfs and T-dwarfs the resonance lines of sodium and potassium are so profoundly pressure-broadened that their wings extend several hundred nanometers from line center. With accurate knowledge of the line profiles as a function of temperature and pressure: such lines can prove to be valuable diagnostics of the atmospheres of such objects. We have initiated a joint program of theoretical and experimental research to study the line-broadening of alkali atom resonance lines due to collisions with species such as helium and molecular hydrogen. Although potassium and sodium are the alkali species of most interest in the atmospheres of cool brown dwarfs and extrasolar giant planets, some of our theoretical focus this year has involved the calculation of pressure-broadening of lithium resonance lines by He, as a test of a newly developed suite of computer codes. In addition, theoretical calculations have been carried out to determine the leading long range van der Waals coefficients for the interactions of ground and excited alkali metal atoms with helium atoms, to within a probable error of 2%. Such data is important in determining the behavior of the resonance line profiles in the far wings. Important progress has been made on the experimental aspects of the program since the arrival of a postdoctoral fellow in September. A new absorption cell has been designed, which incorporates a number of technical improvements over the previous cell, including a larger cell diameter to enhance the signal, and fittings which allow for easier cleaning, thereby significantly reducing the instrument down-time.

Kirby, Kate↗

Carbon Dioxide Dispersion in the Combustion Integrated Rack Simulated Numerically

When discharged into an International Space Station (ISS) payload rack, a carbon dioxide (CO2) portable fire extinguisher (PFE) must extinguish a fire by decreasing the oxygen in the rack by 50 percent within 60 sec. The length of time needed for this oxygen reduction throughout the rack and the length of time that the CO2 concentration remains high enough to prevent the fire from reigniting is important when determining the effectiveness of the response and postfire procedures. Furthermore, in the absence of gravity, the local flow velocity can make the difference between a fire that spreads rapidly and one that self-extinguishes after ignition. A numerical simulation of the discharge of CO2 from PFE into the Combustion Integrated Rack (CIR) in microgravity was performed to obtain the local velocity and CO2 concentration. The complicated flow field around the PFE nozzle exits was modeled by sources of equivalent mass and momentum flux at a location downstream of the nozzle. The time for the concentration of CO2 to reach a level that would extinguish a fire anywhere in the rack was determined using the Fire Dynamics Simulator (FDS), a computational fluid dynamics code developed by the National Institute of Standards and Technology specifically to evaluate the development of a fire and smoke transport. The simulation shows that CO2, as well as any smoke and combustion gases produced by a fire, would be discharged into the ISS cabin through the resource utility panel at the bottom of the rack. These simulations will be validated by comparing the results with velocity and CO2 concentration measurements obtained during the fire suppression system verification tests conducted on the CIR in March 2003. Once these numerical simulations are validated, portions of the ISS labs and living areas will be modeled to determine the local flow conditions before, during, and after a fire event. These simulations can yield specific information about how long it takes for smoke and combustion gases produced by a fire to reach a detector location, how large the fire would be when the detector alarms, and the behavior of the fire until it has been extinguished. This new capability could then be used to optimize the location of fire detectors and fire-suppression ports as well as to evaluate the effectiveness of fire suppressants and response strategies. Numerical data collected from these simulations could also be used to develop a virtual reality fire event for crew training and fire safety awareness. This work is funded by NASA's Bioastronautics Initiative, which has the objective of ensuring and enhancing the health, safety, and performance of humans in space. As part of this initiative, the Microgravity Combustion Science Branch at the NASA Glenn Research Center is conducting spacecraft fire safety research to significantly improve fire safety on inhabited spacecraft.

Wu, Ming-Shin↗

Performance Bounds on Two Concatenated, Interleaved Codes

A method has been developed of computing bounds on the performance of a code comprised of two linear binary codes generated by two encoders serially concatenated through an interleaver. Originally intended for use in evaluating the performances of some codes proposed for deep-space communication links, the method can also be used in evaluating the performances of short-block-length codes in other applications. The method applies, more specifically, to a communication system in which following processes take place: At the transmitter, the original binary information that one seeks to transmit is first processed by an encoder into an outer code (Co) characterized by, among other things, a pair of numbers (n,k), where n (n > k)is the total number of code bits associated with k information bits and n k bits are used for correcting or at least detecting errors. Next, the outer code is processed through either a block or a convolutional interleaver. In the block interleaver, the words of the outer code are processed in blocks of I words. In the convolutional interleaver, the interleaving operation is performed bit-wise in N rows with delays that are multiples of B bits. The output of the interleaver is processed through a second encoder to obtain an inner code (Ci) characterized by (ni,ki). The output of the inner code is transmitted over an additive-white-Gaussian- noise channel characterized by a symbol signal-to-noise ratio (SNR) Es/No and a bit SNR Eb/No. At the receiver, an inner decoder generates estimates of bits. Depending on whether a block or a convolutional interleaver is used at the transmitter, the sequence of estimated bits is processed through a block or a convolutional de-interleaver, respectively, to obtain estimates of code words. Then the estimates of the code words are processed through an outer decoder, which generates estimates of the original information along with flags indicating which estimates are presumed to be correct and which are found to be erroneous. From the perspective of the present method, the topic of major interest is the performance of the communication system as quantified in the word-error rate and the undetected-error rate as functions of the SNRs and the total latency of the interleaver and inner code. The method is embodied in equations that describe bounds on these functions. Throughout the derivation of the equations that embody the method, it is assumed that the decoder for the outer code corrects any error pattern of t or fewer errors, detects any error pattern of s or fewer errors, may detect some error patterns of more than s errors, and does not correct any patterns of more than t errors. Because a mathematically complete description of the equations that embody the method and of the derivation of the equations would greatly exceed the space available for this article, it must suffice to summarize by reporting that the derivation includes consideration of several complex issues, including relationships between latency and memory requirements for block and convolutional codes, burst error statistics, enumeration of error-event intersections, and effects of different interleaving depths. In a demonstration, the method was used to calculate bounds on the performances of several communication systems, each based on serial concatenation of a (63,56) expurgated Hamming code with a convolutional inner code through a convolutional interleaver. The bounds calculated by use of the method were compared with results of numerical simulations of performances of the systems to show the regions where the bounds are tight (see figure).

Moision, Bruce↗

Hall-Effect Thruster Simulations with 2-D Electron Transport and Hydrodynamic Ions

A computational approach that has been used extensively in the last two decades for Hall thruster simulations is to solve a diffusion equation and energy conservation law for the electrons in a direction that is perpendicular to the magnetic field, and use discrete-particle methods for the heavy species. This "hybrid" approach has allowed for the capture of bulk plasma phenomena inside these thrusters within reasonable computational times. Regions of the thruster with complex magnetic field arrangements (such as those near eroded walls and magnets) and/or reduced Hall parameter (such as those near the anode and the cathode plume) challenge the validity of the quasi-one-dimensional assumption for the electrons. This paper reports on the development of a computer code that solves numerically the 2-D axisymmetric vector form of Ohm's law, with no assumptions regarding the rate of electron transport in the parallel and perpendicular directions. The numerical challenges related to the large disparity of the transport coefficients in the two directions are met by solving the equations in a computational mesh that is aligned with the magnetic field. The fully-2D approach allows for a large physical domain that extends more than five times the thruster channel length in the axial direction, and encompasses the cathode boundary. Ions are treated as an isothermal, cold (relative to the electrons) fluid, accounting for charge-exchange and multiple-ionization collisions in the momentum equations. A first series of simulations of two Hall thrusters, namely the BPT-4000 and a 6-kW laboratory thruster, quantifies the significance of ion diffusion in the anode region and the importance of the extended physical domain on studies related to the impact of the transport coefficients on the electron flow field.

plasma↗

Optical Testing Using Portable Laser Coordinate Measuring Instruments

High precision, portable coordinate measuring instruments (CMI) such as laser radars (LR) and laser trackers (LT) have been used for optical system alignment and integration. The LRs ability to perform a non-contact scan of surfaces was previously utilized to characterize large spherical and aspheric mirrors. In this paper, we explore the use of a CMI as an accurate, fast, robust, and non-contact tool for prescription characterization of powered optical surfaces. Using Nikons MV-224350 LR and Leicas Absolute Tracker AT401402 instruments, proof of concept measurements were performed to characterize a variety of optical components by measuring the actual and apparent, or equivalently the direct and through (DT), coordinates of calibrated metrology targets. Custom macros in metrology software and other data reduction code were developed to compute surface-ray intercepts and surface slopes from the DT shots. The calculated data is fit to an aspheric surface formula to obtain the optimum prescription. The results were compared to the nominal parameters and were crosschecked using LR scans or other approaches. We discuss potential applications across the fields of optical component fabrication and system alignment and testing.

Khreishi, M.↗

Optical Testing Using Portable Laser Coordinate Measuring Instruments

High precision, portable coordinate measuring instruments (CMI) such as laser radars (LR) and laser trackers (LT) have been used for optical system alignment and integration. The LR's ability to perform a non-contact scan of surfaces was previously utilized to characterize large spherical and aspheric mirrors. In this paper, we explore the use of a CMI as an accurate, fast, robust, and non-contact tool for prescription characterization of powered optical surfaces. Using Nikon's MV-224/350 LR and Leica's Absolute Tracker AT401/402 instruments, proof of concept measurements were performed to characterize a variety of optical components by measuring the actual and apparent, or equivalently the "direct and through" (D&T), coordinates of calibrated metrology targets. Custom macros in metrology software and other data reduction code were developed to compute surface-ray intercepts and surface slopes from the D&T shots. The calculated data is fit to an aspheric surface formula to obtain the optimum prescription. The results were compared to the nominal parameters and were crosschecked using LR scans or other approaches. We discuss potential applications across the fields of optical component fabrication and system alignment and testing.

Khreishi, Manal↗

The SPACE Computer Code for Analyzing the International Space Station Electrical Power System: Past, Present, and Future

The System Power Analysis for Capability Evaluation (SPACE) computer code was initially developed by NASA in 1988 to assess the Space Station Freedom electric power system and later adapted to support contractor electrical power system capability analyses for the International Space Station (ISS). Over time, the code has supported many efforts such as ISS redesign activities in the early 1990s, assessment of time-phased loads against power system operating limits for future ISS assembly flights (including Certification of Flight Readiness reviews by the ISS program office), and determining the optimum solar array gimbal positions while respecting keep-out zones which minimize both solar array contamination and structural loads. The code has been validated by comparisons with ISS on-orbit data in multiple validation episodes. Recent updates to the code include the incorporation of a Lithium-Ion battery model in addition to the Nickel Hydrogen battery model and modifications to the solar array degradation model to better match on-orbit test results. SPACE has also been extended beyond the ISS to include modeling of the Orion Multi-Purpose Crew Vehicle electrical power system (SPACE-MPCV) and Mars Surface Electrical Power Systems (MSEPS). Portions of SPACE were integrated with a trajectory code to form a Solar Electric Propulsion Simulation (SEPSim), which can be used for analyzing solar electric propulsion missions. In addition, SPACE methods and subroutines have been adapted to a multitude of other projects. This paper summarizes the initial code development and subsequent code utilization in the context of the overall ISS program development and on-orbit operations. Recent updates and results from the code are discussed, including preliminary analyses for the Orion power system.

Electrical Power System↗

Microscale Analysis of Spacecraft Heat Shields

Imagine entering Earth’s atmosphere after returning from the outer solar system. A heat shield less than 2 inches thick protects you from temperatures up to 2,900° Celsius (5,252° Fahrenheit). Such conditions were experienced by NASA’s Stardust capsule during reentry in 2006. The only materials capable of providing the necessary protection are composites with complex microstructures. Evaluating these materials is difficult, requiring precise knowledge of their properties. To this end, NASA scientists are developing research codes to compute material properties and simulate ablation at the microscale using agency supercomputers. Utilizing these tools, along with experiments, researchers are working to push the limits of spaceflight, allowing for greater flexibility in future space missions.

Thornton, John M.↗

Bumper Implementation of LMEEM

The Hypervelocity Impact Technology (HVIT) group at the NASA Johnson Space Center (JSC) maintains the Bumper 3 computer program (referred to as “Bumper”) to perform Micrometeoroid and Orbital Debris (MMOD) risk assessments for spacecraft. To perform its calculations, this program requires detailed mathematical models of the space environments and can include those models representing orbital debris, micrometeoroids, and lunar ejecta. The lunar ejecta environment is a projectile environment specific to the lunar surface. Due to the Moon’s lack of an atmosphere, meteoroids frequently strike the lunar surface. When this happens, the impact ejects material radially away from the impact point, which can pose a potential impact hazard to surface-based hardware. This lunar ejecta environment was previously described in 1969, in NASA Special Publication SP-8013, which was used as the basis for a model that will be referred to in this document as the “SP-8013 environment model,” or simply, “SP-8013.” Since around 2020, this model has been used in MMOD risk assessments for the Human Landing System (HLS) program and other lunar surface programs. A new lunar meteoroid ejecta environment was developed via a computer code named Lunar Meteoroid Ejecta Engineering Model (LMEEM), produced by the Natural Environments Branch (NEB) at Marshall Space Flight Center (MSFC). A global lunar surface environment definition using LMEEM is to be provided in the Cross-Program Design Specification for Natural Environments (DSNE) Revision J in 2024. Contact the NASA MSFC Natural Environments Branch for an in-depth discussion of the LMEEM environment model and its derivation. This environment model will supersede SP-8013 for MMOD risk assessments. To perform such assessments, Bumper must read data from text files output by LMEEM and use that data in risk calculations. The purpose of this document is to: •Demonstrate verification that the Bumper code has properly implemented the LMEEM environment model (Section 3) •Discuss the LMEEM environment as understood through an initial environment file (Section 4) •Compare the calculated risk due to the SP-8013 and LMEEM models (Section 5)

Bumper↗

Near-Horizontal, Two-Phase Flow Patterns of Nitrogen and Hydrogen at Low Mass Heat and Flux (on CD-ROM)

One reason for NASA's interest in cryogenic two-phase flow with low mass and heat flux is the need to design spacecraft heat exchangers used for vaporizing cryogenic propellants. The CD-ROM provides digitized movies of particular flow patterns observed in experimental work. The movies have been provided in (QuickTime9Trademark) format, encoded at 320w x 240h pixels, 15 fps, using the Sorenson(Trademark) Video Codec for compression. Experiments were conducted to obtain data on the two-phase (liquid and vapor) flow behavior of cryogenic nitrogen and hydrogen under low mass and heat flux conditions. Tests were performed in normal gravity with a 1.5 degree up flow configuration. View ports in the apparatus permitted visual observation of the two-phase flow patterns. Computer codes to predict flow patterns were developed from theoretical/empirical models reported in the literature. Predictions from the computer codes were compared with experimental flow pattern observations. Results are presented employing the traditional two-dimensional flow pattern map format using the liquid and gas superficial velocities as coordinates. In general, the agreement between the experimental results and the analytical predictive methods is reasonably good. Small regions of the flow pattern maps are identified where the models are deficient as a result of neglecting phase change phenomena. Certain regions of the maps were beyond the range of the experiments and could not be completely validated. Areas that could benefit from further work include modeling of the transition from separated flow, collection of additional data in the bubble and annular flow regimes, and collection of experimental data at other inclination angles, tube diameters and high heat flux.

TWO PHASE FLOW↗

Burner liner thermal/structural load modelling

The objective of this program is to develop a thermal data transfer computer program module for the burner liner thermal structural load modeling program. This will be accomplished by (1) reviewing existing methodologies for thermal data transfer and selecting three heat transfer codes for application in this program; (2) evaluating the selected codes to establish criteria for developing a computer program module to transfer thermal data from the heat transfer codes to selected stress analysis codes; (3) developing the automated thermal load transfer module; and (4) verifying and documenting the module. In aircraft turbine engine hot section components, cyclic thermal stresses are the most important damage mechanism. Consequently, accurate and reliable prediction of thermal loads is essential to improving durability. To achieve this goal, a considerable effort over the past 20 years has been devoted to the acquisition of engine temperature test data, as well as the development of accurate, reliable, and efficient computer codes for the prediction of steady state and transient temperatures and for the calculation of elastic and inelastic cyclic stresses and strains in hot section components. There is a need for continued development of these codes, because the availability of more accurate analysis techniques for complex configurations has enabled engine designers to use more sophisticated designs to achieve higher cycle efficiency and reduce weight.

Maffeo, R. J.↗

A three-dimensional upwind parabolized Navier-Stokes code for chemically reacting flows

A new upwind, parabolized Navier-Stokes (PNS) code has been developed to compute the three-dimensional flow of chemically reacting air around hypersonic vehicles. The code is a modification of the perfect gas, three-dimensional UPS code of Lawrence et al. (1986) which has been extended in the present study to permit the calculation of hypersonic, viscous flows in chemical nonequilibrium. The algorithm solves the PNS equations using a finite-volume, upwind TVD method based on Roe's approximate Riemann solver that has been modified to account for real gas effects. The present code solves the fluid dynamic and species continuity equations in a loosely-coupled manner. The fluid medium is assumed to be a chemically reacting mixture of thermally perfect (but calorically imperfect) gases in thermal equilibrium. Results are presented for the hypersonic laminar flow over a cone at 0- and 10-deg angles of attack and for a generic hypersonic vehicle. Calculations are performed assuming either perfect gas, equilibrium air, or finite-rate chemistry.

Buelow, Philip E.↗