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 649 records · Page 36

NASA Multidimensional Stirling Convertor Code Developed

A high-efficiency Stirling Radioisotope Generator (SRG) for use on potential NASA Space Science missions is being developed by the Department of Energy, Lockheed Martin, Stirling Technology Company, and the NASA Glenn Research Center. These missions may include providing spacecraft onboard electric power for deep space missions or power for unmanned Mars rovers. Glenn is also developing advanced technology for Stirling convertors, aimed at substantially improving the specific power and efficiency of the convertor and the overall power system. Performance and mass improvement goals have been established for second- and third-generation Stirling radioisotope power systems. Multiple efforts are underway to achieve these goals, both in house at Glenn and under various grants and contracts. These efforts include the development of a multidimensional Stirling computational fluid dynamics code, high-temperature materials, advanced controllers, an end-to-end system dynamics model, low-vibration techniques, advanced regenerators, and a lightweight convertor. Under a NASA grant, Cleveland State University (CSU) and its subcontractors, the University of Minnesota (UMN) and Gedeon Associates, have developed a twodimensional computer simulation of a CSUmod Stirling convertor. The CFD-ACE commercial software developed by CFD Research Corp. of Huntsville, Alabama, is being used. The CSUmod is a scaled version of the Stirling Technology Demonstrator Convertor (TDC), which was designed and fabricated by the Stirling Technology Company and is being tested by NASA. The schematic illustrates the structure of this model. Modeled are the fluid-flow and heat-transfer phenomena that occur in the expansion space, the heater, the regenerator, the cooler, the compression space, the surrounding walls, and the moving piston and displacer. In addition, the overall heat transfer, the indicated power, and the efficiency can be calculated. The CSUmod model is being converted to a two-dimensional model of the TDC at NASA Glenn. Validation of the multidimensional Stirling code is an important part of the grant effort. UMN has been generating data in an oscillating-flow test facility using two different test sections: a 90 turn and a cooler/regenerator/heater test section. CSU has created computational fluid dynamics models of both these test sections and has been making comparisons with the data, then improving their models to improve the agreement with the test data. CSU has also been using data available in the literature for code validation. UMN is now preparing to begin fabrication of a new 180 turn test section that will be more representative of certain portions of the Stirling engine geometry. Simulations to almost periodic steady state with the two-dimensional CSUmod model indicate that, to reach periodic steady state on a single 2-GHz desktop computer, 75 to 100 complete simulation cycles would be required and between 1 and 2 months of computer time. Therefore, Glenn has purchased the first 8 computers, of a 64-computer cluster, to be run in parallel to accelerate the simulation. On the basis of CFD Research Corp.'s experience with running the parallelized version of CFD-ACE on their clusters, we estimate that the complete 64-computer cluster will reduce simulation computing time by a factor of about 40. Plans are to continue development of these multidimensional Stirling codes and to use them to study the fluid-flow and heat-transfer phenomena that occur inside Stirling convertors. This is expected to lead to improved thermodynamic loss understanding, onedimensional design and performance codes, and engine performance.

Tew, Roy C.↗

Component-specific modeling

A series of interdisciplinary modeling and analysis techniques that were specialized to address three specific hot section components are presented. These techniques will incorporate data as well as theoretical methods from many diverse areas including cycle and performance analysis, heat transfer analysis, linear and nonlinear stress analysis, and mission analysis. Building on the proven techniques already available in these fields, the new methods developed will be integrated into computer codes to provide an accurate, efficient and unified approach to analyzing combustor burner liners, hollow air-cooled turbine blades and air-cooled turbine vanes. For these components, the methods developed will predict temperature, deformation, stress and strain histories throughout a complete flight mission.

Mcknight, R. L.↗

Intelligent Patching of Conceptual Geometry for CFD Analysis

The iPatch computer code for intelligently patching surface grids was developed to convert conceptual geometry to computational fluid dynamics (CFD) geometry (see figure). It automatically uses bicubic B-splines to extrapolate (if necessary) each surface in a conceptual geometry so that all the independently defined geometric components (such as wing and fuselage) can be intersected to form a watertight CFD geometry. The software also computes the intersection curves of surface patches at any resolution (up to 10.4 accuracy) specified by the user, and it writes the B-spline surface patches, and the corresponding boundary points, for the watertight CFD geometry in the format that can be directly used by the grid generation tool VGRID. iPatch requires that input geometry be in PLOT3D format where each component surface is defined by a rectangular grid {(x(i,j), y(i,j), z(i,j)):1less than or equal to i less than or equal to m, 1 less than or equal to j less than or equal to n} that represents a smooth B-spline surface. All surfaces in the PLOT3D file conceptually represent a watertight geometry of components of an aircraft on the half-space y greater than or equal to 0. Overlapping surfaces are not allowed, but could be fixed by a utility code "fixp3d". The fixp3d utility code first finds the two grid lines on the two surface grids that are closest to each other in Hausdorff distance (a metric to measure the discrepancies of two sets); then uses one of the grid lines as the transition line, extending grid lines on one grid to the other grid to form a merged grid. Any two connecting surfaces shall have a "visually" common boundary curve, or can be described by an intersection relationship defined in a geometry specification file. The intersection of two surfaces can be at a conceptual level. However, the intersection is directional (along either i or j index direction), and each intersecting grid line (or its spine extrapolation) on the first surface should intersect the second surface. No two intersection relationships will result in a common intersection point of three surfaces. The output files of iPatch are IGES, d3m, and mapbc files that define the CFD geometry in VGRID format. The IGES file gives the NURBS definition of the outer mold line in the geometry. The d3m file defines how the outer mold line is broken into surface patches whose boundary curves are defined by points. The mapbc file specifies what the boundary condition is on each patch and the corresponding NURBS surface definition of each non-planar patch in the IGES file.

Li, Wu↗

SEDS Tether M/OD Damage Analyses

The Small Expendable Deployer System (SEDS) was designed to deploy an endmass at the end of a 20-km-long tether which acts as an upper stage rocket, and the threats from the meteoroid and orbital debris (M/OD) particle environments on SEDS components are important issues for the safety and success of any SEDS mission. However, the possibility of severing the tether due to M/OD particle impacts is an even more serious concern, since the SEDS tether has a relatively large exposed area to the M/OD environments although its diameter is quite small. The threats from the M/OD environments became a very important issue for the third SEDS mission, since the project office proposed using the shuttle orbiter as a launch platform instead of the second stage of a Delta II expendable rocket, which was used for the first two SEDS missions. A series of hyper-velocity impact tests were performed at the Johnson Space Center and Arnold Engineering Development Center to help determine the critical particle sizes required to sever the tether. The computer hydrodynamic code or hydrocode called CTH, developed by the Sandia National Laboratories, was also used to simulate the damage on the SEDS tether caused by both the orbital debris and test particle impacts. The CTH hydrocode simulation results provided the much needed information to help determine the critical particle sizes required to sever the tether. The M/OD particle sizes required to sever the tether were estimated to be less than 0.1 cm in diameter from these studies, and these size particles are more abundant in low-Earth orbit than larger size particles. Finally, the authors performed the M/OD damage analyses for the three SEDS missions; i.e., SEDS-1, -2, and -3 missions, by using the information obtained from the hypervelocity impact test and hydrocode simulations results.

Hayashida, K. B.↗

An improved supersonic, three-dimensional, external, inviscid flow field code

A numerical procedure was developed to compute the inviscid super/hypersonic flow fields about complex vehicle geometries accurately and efficiently. A second-order accurate finite difference scheme is used to integrate the three-dimensional Euler equations in regions of continuous flow, while all shock waves are computed as discontinuities via the Rankine-Hugoniot jump conditions. Conformal mappings are used to develop a computational grid. The effects for equilibrium air are included using curve fits of Mollier charts. This report deals only with modifications to these procedures in four specific areas: inlet mass ingestion, subsonic axial Mach number, improved conformal mappings, and vehicles flying at yaw. In each area both the modifications to the computational procedures and computer code are discussed.

Marconi, F.↗

First principles numerical model of avalanche-induced arc discharges in electron-irradiated dielectrics

The model consists of four phases: single electron dynamics, single electron avalanche, negative streamer development, and tree formation. Numerical algorithms and computer code implementations are presented for the first three phases. An approach to developing a code description of fourth phase is discussed. Numerical results are presented for a crude material model of Teflon.

Beers, B. L.↗

RWF rotor-wake-fuselage code software reference guide

The RWF (Rotor-Wake-Fuselage) code was developed from first principles to compute the aerodynamics associated with the complex flow field of helicopter configurations. The code is sized for a single, multi-bladed main rotor and any configuration of non-lifting fuselage. The mathematical model for the RWF code is based on the integration of the momentum equations and Green's theorem. The unknowns in the problem are the strengths of prescribed singularity distributions on the boundaries of the flow. For the body (fuselage) a surface of constant strength source panels is used. For the rotor blades and rotor wake a surface of constant strength doublet panels is used. The mean camber line of the rotor airfoil is partitioned into surface panels. The no-flow boundary condition at the panel centroids is modified at each azimuthal step to account for rotor blade cyclic pitch variation. The geometry of the rotor wake is computers at each time step of the solution. The code produces rotor and fuselage surface pressures, as well as the complex geometry of the evolving rotor wake.

Berry, John D.↗

The computation of supersonic combustor flows using multi-computers

An explicit CFD computer code with parallel processing is developed for the purpose of simulating internal high-speed reacting flows. The code solves the 3D Navier-Stokes equations for compressible flows. Preliminary results indicate that parallel processing should solve the resource and speed shortages presently encountered with sequential computers for problems of interest.

Underwood, Mark↗

Procedures for the computation of unsteady transonic flows including viscous effects

Modifications of the code LTRAN2, developed by Ballhaus and Goorjian, which account for viscous effects in the computation of planar unsteady transonic flows are presented. Two models are considered and their theoretical development and numerical implementation is discussed. Computational examples employing both models are compared with inviscid solutions and with experimental data. Use of the modified code is described.

Rizzetta, D. P.↗

Computer program for prediction of the deposition of material released from fixed and rotary wing aircraft

This is a user manual for the computer code ""AGDISP'' (AGricultural DISPersal) which has been developed to predict the deposition of material released from fixed and rotary wing aircraft in a single-pass, computationally efficient manner. The formulation of the code is novel in that the mean particle trajectory and the variance about the mean resulting from turbulent fluid fluctuations are simultaneously predicted. The code presently includes the capability of assessing the influence of neutral atmospheric conditions, inviscid wake vortices, particle evaporation, plant canopy and terrain on the deposition pattern.

Teske, M. E.↗

Alfven wave resonances and flow induced by non-linear Alfven waves in a stratified atmosphere

A nonlinear, time-dependent, ideal MHD code has been developed and used to compute the flow induced by nonlinear Alfven waves propagating in an isothermal, stratified, plane-parallel atmosphere. The code is based on characteristic equations solved in a Lagrangian frame and is highly accurate. Results show that resonance behavior of Alfven waves exists in the presence of a continuous density gradient and that the waves with periods corresponding to resonant peaks exert considerably more force on the medium than off-resonance periods; this leads to enhanced flow. If only off-peak periods are considered, the relationship between the wave period and induced longitudinal velocity shows that short period WKB waves push more on the background medium than longer period, non-WKB, waves. The results also show the development of the longitudinal waves produced by the finite amplitude of the Alfven waves. The longitudinal wave becomes strong as the Alfven wave relative amplitude grows above 10 percent and will lead to strong damping of the Alfven waves.

Stark, B. A.↗

Developing a gas rocket performance prediction technique

A simple, semi-empirical performance correlation/prediction technique applicable to gaseous and liquid propellant rocket engines is presented. Excellent correlations were attained for over 100 test firings by adjusting the computation of the gaseous mixing of an unreactive, coaxial jet using a correlation factor, F, which resulted in prediction of the experimental combustion efficiency for each firing. Static pressure, mean velocity and turbulence intensity in the developing region of non-reactive coaxial jets, typical of those of coaxial injector elements were determined. Detailed profiles were obtained at twelve axial locations (extending from the nozzle exit for a distance of five diameters) downstream from a single element of the Bell Aerospace H2/O2 19-element coaxial injector. These data are compared with analytical predictions made using both eddy viscosity and turbulence kinetic energy mixing models and available computer codes. Comparisons were disappointing, demonstrating the necessity of developing improved turbulence models and computational techniques before detailed predictions of practical coaxial free jet flows are attempted.

Morgenthaler, J. H.↗

Effectiveness evaluation of STOL transport operations

A short-takeoff and landing (STOL) systems simulation model has been developed and implemented in a computer code (known as STOL OPS) which permits evaluation of the operation of a STOL aircraft and its avionics in a commercial airline operating environment. STOL OPS concentrated on the avionics functions of navigation, guidance, control, communication, hazard aviodance, and systems management. External world factors influencing the operation of the STOL aircraft include each airport and its geometry, air traffic at each airport, air traffic control equipment and procedures, weather (including winds and visibility), and the flight path between each airport served by the route. The development of the STOL OPS program provides NASA a set of computer programs which can be used for detailed analysis of a STOL aircraft and its avionics and permit establishment of system requirements as a function of airline mission performance goals.

Hitt, E. F.↗

Navier-Stokes Computations on Commodity Computers

In this paper we discuss and demonstrate the feasibility of solving high-fidelity, nonlinear computational fluid dynamics (CFD) problems of practical interest on commodity machines, namely Pentium Pro PC's. Such calculations have now become possible due to the progress in computational power and memory of the off-the-shelf commodity computers, along with the growth in bandwidth and communication speeds of networks. A widely used CFD code known as TLNS3D, which was developed originally on large shared memory computers was selected for this effort. This code has recently been ported to massively parallel processor (MPP) type machines, where natural partitioning along grid blocks is adopted in which one or more blocks are distributed to each of the available processors. In this paper, a similar approach is adapted to port this code to a cluster of Pentium Pro computers. The message passing among the processors is accomplished through the use of standard message passing interface (MPI) libraries. Scaling studies indicate fairly high level of parallelism on such clusters of commodity machines, thus making solutions to Navier-Stokes equations for practical problems more affordable.

Vatsa, Veer N.↗

Unsteady Full Annulus Simulations of a Transonic Axial Compressor Stage

Two recent research endeavors in turbomachinery at NASA Glenn Research Center have focused on compression system stall inception and compression system aerothermodynamic performance. Physical experiment and computational research are ongoing in support of these research objectives. TURBO, an unsteady, three-dimensional, Navier-Stokes computational fluid dynamics code commissioned and developed by NASA, has been utilized, enhanced, and validated in support of these endeavors. In the research which follows, TURBO is shown to accurately capture compression system flow range-from choke to stall inception-and also to accurately calculate fundamental aerothermodynamic performance parameters. Rigorous full-annulus calculations are performed to validate TURBO s ability to simulate the unstable, unsteady, chaotic stall inception process; as part of these efforts, full-annulus calculations are also performed at a condition approaching choke to further document TURBO s capabilities to compute aerothermodynamic performance data and support a NASA code assessment effort.

Herrick, Gregory P.↗

Mathematical Modeling of Loop Heat Pipes with Multiple Capillary Pumps and Multiple Condensers: Stead State Stimulations - Part 1

Loop Heat Pipes (LHPs) have proven themselves as reliable and robust heat transport devices for spacecraft thermal control systems. So far, the LHPs in earth-orbit satellites perform very well as expected. Conventional LHPs usually consist of a single capillary pump for heat acquisition and a single condenser for heat rejection. Multiple pump/multiple condenser LHPs have shown to function very well in ground testing. Nevertheless, the test results of a dual pump/condenser LHP also revealed that the dual LHP behaved in a complicated manner due to the interaction between the pumps and condensers. Thus it is redundant to say that more research is needed before they are ready for 0-g deployment. One research area that perhaps compels immediate attention is the analytical modeling of LHPs, particularly the transient phenomena. Modeling a single pump/single condenser LHP is difficult enough. Only a handful of computer codes are available for both steady state and transient simulations of conventional LHPs. No previous effort was made to develop an analytical model (or even a complete theory) to predict the operational behavior of the multiple pump/multiple condenser LHP systems. The current research project offered a basic theory of the multiple pump/multiple condenser LHP operation. From it, a computer code was developed to predict the LHP saturation temperature in accordance with the system operating and environmental conditions.

Hoang, Triem T.↗

A rapid perturbation procedure for determining nonlinear flow solutions: Application to transonic turbomachinery flows

Perturbation procedures and associated computational codes for determining nonlinear flow solutions were developed to establish a method for minimizing computational requirements associated with parametric studies of transonic flows in turbomachines. The procedure that was developed and evaluated was found to be capable of determining highly accurate approximations to families of strongly nonlinear solutions which are either continuous or discontinuous, and which represent variations in some arbitrary parameter. Coordinate straining is employed to account for the movement of discontinuities and maxima of high gradient regions due to the perturbation. The development and results reported are for the single parameter perturbation problem. Flows past both isolated airfoils and compressor cascades involving a wide variety of flow and geometry parameter changes are reported. Attention is focused in particular on transonic flows which are strongly supercritical and exhibit large surface shock movement over the parametric range studied; and on subsonic flows which display large pressure variations in the stagnation and peak suction pressure regions. Comparisons with the corresponding 'exact' nonlinear solutions indicate a remarkable accuracy and range of validity of such a procedure.

Stahara, S. S.↗

Optimization of thermal protection systems for the space shuttle vehicle. Volume 1: Final report

A study performed to continue development of computational techniques for the Space Shuttle Thermal Protection System is reported. The resulting computer code was used to perform some additional optimization studies on several TPS configurations. The program was developed in Fortran 4 for the CDC 6400, and it was converted to Fortran 5 to be used for the Univac 1108. The computational methodology is developed in modular fashion to facilitate changes and updating of the techniques and to allow overlaying the computer code to fit into approximately 131,000 octal words of core storage. The program logic involves subroutines which handle input and output of information between computer and user, thermodynamic stress, dynamic, and weight/estimate analyses of a variety of panel configurations. These include metallic, ablative, RSI (with and without an underlying phase change material), and a thermodynamic analysis only of carbon-carbon systems applied to the leading edge and flat cover panels. Two different thermodynamic analyses are used. The first is a two-dimensional, explicit precedure with variable time steps which is used to describe the behavior of metallic and carbon-carbon leading edges. The second is a one-dimensional implicity technique used to predict temperature in the charring ablator and the noncharring RSI. The latter analysis is performed simply by suppressing the chemical reactions and pyrolysis of the TPS material.

Source record↗