Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer codes”

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 937 records · Page 52

Finite-difference solution of the compressible stability eigenvalue problem

A compressible stability analysis computer code is developed. The code uses a matrix finite difference method for local eigenvalue solution when a good guess for the eigenvalue is available and is significantly more computationally efficient than the commonly used initial value approach. The local eigenvalue search procedure also results in eigenfunctions and, at little extra work, group velocities. A globally convergent eigenvalue procedure is also developed which may be used when no guess for the eigenvalue is available. The global problem is formulated in such a way that no unstable spurious modes appear so that the method is suitable for use in a black box stability code. Sample stability calculations are presented for the boundary layer profiles of a Laminar Flow Control (LFC) swept wing.

Malik, M. R.↗

Automatic Processing of Reactive Polymers

A series of process modeling computer codes were examined. The codes use finite element techniques to determine the time-dependent process parameters operative during nonisothermal reactive flows such as can occur in reaction injection molding or composites fabrication. The use of these analytical codes to perform experimental control functions is examined; since the models can determine the state of all variables everywhere in the system, they can be used in a manner similar to currently available experimental probes. A small but well instrumented reaction vessel in which fiber-reinforced plaques are cured using computer control and data acquisition was used. The finite element codes were also extended to treat this particular process.

Roylance, D.↗

Fluid management technology: Liquid slosh dynamics and control

Flight experiments were defined for the Cryogenic On-Orbit Liquid Depot Storage, Acquisition and Transfer Satellite (COLD-SAT) test bed satellite and the Shuttle middeck to help establish the influence of the gravitational environment on liquid slosh dynamics and control. Several analytical and experimental studies were also conducted to support the experiments and to help understand the anticipated results. Both FLOW-3D and NASA-VOF3D computer codes were utilized to simulate low Bond number, small amplitude sloshing, for which the motions are dominated by surface forces; it was found that neither code provided a satisfactory simulation. Thus, a new analysis of low Bond number sloshing was formulated, using an integral minimization technique that will allow the assumptions made about surface physics phenomena to be modified easily when better knowledge becomes available from flight experiments. Several examples were computed by the innovative use of a finite-element structural code. An existing spherical-pendulum analogy of nonlinear, rotary sloshing was also modified for easier use and extended to low-gravity conditions. Laboratory experiments were conducted to determine the requirements for liquid-vapor interface sensors as a method of resolving liquid surface motions in flight experiments. The feasibility of measuring the small slosh forces anticipated in flight experiments was also investigated.

Dodge, Franklin T.↗

Creating a Simple Single Computational Approach to Modeling Rarefied and Continuum Flow About Aerospace Vehicles

We proposed to create a single computational code incorporating methods that can model both rarefied and continuum flow to enable the efficient simulation of flow about space craft and high altitude hypersonic aerospace vehicles. The code was to use a single grid structure that permits a smooth transition between the continuum and rarefied portions of the flow. Developing an appropriate computational boundary between the two regions represented a major challenge. The primary approach chosen involves coupling a four-speed Lattice Boltzmann model for the continuum flow with the DSMC method in the rarefied regime. We also explored the possibility of using a standard finite difference Navier Stokes solver for the continuum flow. With the resulting code we will ultimately investigate three-dimensional plume impingement effects, a subject of critical importance to NASA and related to the work of Drs. Forrest Lumpkin, Steve Fitzgerald and Jay Le Beau at Johnson Space Center. Below is a brief background on the project and a summary of the results as of the end of the grant.

Goldstein, David B.↗

Grid generation for turbomachinery problems

The development of a computer code to generate numerical grids for complex internal flow systems such as the fluid flow inside the space shuttle main engine is outlined. The blending technique for generating a grid for stator-rotor combination at a particular radial section is examined. The computer programs which generate these grids are listed in the Appendices. These codes are capable of generating grids at different cross sections and thus providng three dimensional stator-rotor grids for the turbomachinery of the space shuttle main engine.

Steinhoff, J.↗

Simulation of Hypervelocity Impact on Aluminum-Nextel-Kevlar Orbital Debris Shields

An improved hybrid particle-finite element method has been developed for hypervelocity impact simulation. The method combines the general contact-impact capabilities of particle codes with the true Lagrangian kinematics of large strain finite element formulations. Unlike some alternative schemes which couple Lagrangian finite element models with smooth particle hydrodynamics, the present formulation makes no use of slidelines or penalty forces. The method has been implemented in a parallel, three dimensional computer code. Simulations of three dimensional orbital debris impact problems using this parallel hybrid particle-finite element code, show good agreement with experiment and good speedup in parallel computation. The simulations included single and multi-plate shields as well as aluminum and composite shielding materials. at an impact velocity of eleven kilometers per second.

Fahrenthold, Eric P.↗

A theoretical study of a laminar diffusion flame

Theoretical models of an axisymmetric laminar diffusion flame are discussed, with an emphasis on the behavior of such flames at increasing pressures. The flame-sheet or Burke-Schumann model (in terms of Bessel functions) and various boundary layer numerical solutions are presented and their results compared with experimental data. The most promising theoretical model combines the numerical flow field solution of the Patankar-Spalding computer code with the Pratt-Wormeck chemical reaction subroutine. The flame shapes for pressures of 1, 5, 10, 20, and 50 atmospheres were computed and agree remarkably well with experimental data. There is a noticeable shape change with pressure, believed to be a result of buoyancy effects. The chemical concentration profiles do not exhibit much dependence on pressure, a reflection of the fact that only one chemical mechanism was utilized at all pressures.

Frair, K. L.↗

Recent advances in the modelling of crack growth under fatigue loading conditions

Fatigue crack growth associated with cyclic (secondary) plastic flow near a crack front is modelled using an incremental formulation. A new description of threshold behaviour under small load cycles is included. Quasi-static crack extension under high load excursions is described using an incremental formulation of the R-(crack growth resistance)- curve concept. The integration of the equations is discussed. For constant amplitude load cycles the results will be compared with existing crack growth laws. It will be shown that the model also properly describes interaction effects of fatigue crack growth and quasi-static crack extension. To evaluate the more general applicability the model is included in the NASGRO computer code for damage tolerance analysis. For this purpose the NASGRO program was provided with the CORPUS and the STRIP-YIELD models for computation of the crack opening load levels. The implementation is discussed and recent results of the verification are presented.

Dekoning, A. U.↗

Program Helps To Determine Chemical-Reaction Mechanisms

General Chemical Kinetics and Sensitivity Analysis (LSENS) computer code developed for use in solving complex, homogeneous, gas-phase, chemical-kinetics problems. Provides for efficient and accurate chemical-kinetics computations and provides for sensitivity analysis for variety of problems, including problems involving honisothermal conditions. Incorporates mathematical models for static system, steady one-dimensional inviscid flow, reaction behind incident shock wave (with boundary-layer correction), and perfectly stirred reactor. Computations of equilibrium properties performed for following assigned states: enthalpy and pressure, temperature and pressure, internal energy and volume, and temperature and volume. Written in FORTRAN 77 with exception of NAMELIST extensions used for input.

Bittker, D. A.↗

The Use of Steady and Pulsed Detonations for Propulsion Systems

Objectives of the ODWE concept studies are: demonstrate the feasibility of the oblique detonation wave engine (ODWE) for hypersonic propulsion; demonstrate the existance and stability of an oblique detonation wave in hypersonic wind tunnels; develop engineering codes which predict the performance characteristics of the ODWE including specific impulse and thrust coefficients for various operating conditions; develop multi-dimensional computer codes which can model all aspects of the ODWE including fuel injection, mixing, ignition, combustion and expansion with fully detailed chemical kinetics and turbulence models; and validate the codes with experimental data use the simulations to predict the ODWE performance for conditions not easily obtained in wind tunnels.

Adelman, Henry G.↗

Anomalous Results from PO Applied to Reflector Antennas: The importance of Near Field Computations

There is no doubt among the designers of reflector antennas that the physical optics (PO) analysis technique is the most popular numerical technique. Powerful computer codes are available for the analysis of single or multi reflector antenna systems. Additionally, ever increasingly demand on the antenna performance necessitates the computation of antenna far field patterns under various situations. For example, in using multi reflector antennas such as, Gregorian or Cassegrain, it may become necessary to determine the total fields including the feed radiation pattern, subreflector scattered pattern and the main reflector scattered pattern. In these situations, the common practice is to sum up various scattered fields and the incident field contributions to obtain the desired total field, it is the purpose of this paper to demonstrate that the typical approach based on the far field pattern of the feed would result into erroneous result and special care must be exercised to obtain the correct result. This will be demonstrated through a detailed investigation of a representative test case.

Rahmat-Samii, Yahya↗

Efficient computation of the stability of three-dimensional compressible boundary layers

Methods for the computer analysis of the stability of three-dimensional compressible boundary layers are discussed and the user-oriented Compressible Stability Analysis (COSAL) computer code is described. The COSAL code uses a matrix finite-difference method for local eigenvalue solution when a good guess for the eigenvalue is available and is significantly more computationally efficient than the commonly used initial-value approach. The local eigenvalue search procedure also results in eigenfunctions and, at little extra work, group velocities. A globally convergent eigenvalue procedure is also developed which may be used when no guess for the eigenvalue is available. The global problem is formulated in such a way that no unstable spurious modes appear so that the method is suitable for use in a black-box stability code. Sample stability calculations are presented for the boundary layer profiles of an LFC swept wing.

Malik, M. R.↗

Linear stability analysis of three-dimensional compressible boundary layers

A compressible stability analysis computer code is developed. The code uses a matrix finite-difference method for local eigenvale solution when a good guess for the eigenvalue is available and is significantly more computationally efficient than the commonly used inital-value approach. The local eigenvalue search procedure also results in eigenfunctions and, at little extra work, group velocities. A globally convergent eigenvalue procedure is also developed that may be used when no guess for the eigenvalue is available. The global problem is formulated in such a way that no unstable spurious modes appear so that the method is suitable for use in a black-box stability code. Sample stability calculations are presented for the boundary layer profiles of an LFC swept wing.

Malik, Mujeeb R.↗

Flow field analysis for a class of waverider configurations

A package of computer codes for analysis of flow fields for waverider configurations is described. The package consists of a surface/volume grid generator and a finite-volume flow solver. The grid generator defines body geometries and computational grids by an algebraic homotopy procedure. The algebraic procedure is versatile in its application and can readily generate configurations in the class of blended wing-body geometries. This code has the ability to produce a wide variety of geometries in the given class with varying geometrical attributes. The flow solver employs a finite-volume formation and solves the explicit, Runge-Kutta integration technique. The method or flow simulation incorporates several techniques for acceleration of the convergence of the interaction process and an entropy corrected enthalpy damping procedure for efficient computation of high Mach number flows.

Moitra, Anutosh↗

The free-wake prediction of rotor hover performance using a vortex embedding method

A method is developed to predict the rotor hover performance. This method solves the compressible mass conservation equation much like current full potential codes and can therefore predict the transonic flows on a rotor. However, the newly developed approach also allows for the free convection of shed vorticity and permits the computation of the entire hover wake system. The method uses a vortex embedding scheme in potential flow and has been implemented in a computer code, HELIX -I. To predict power we implement a simple boundary layer and two different induced-drag integration schemes. The induced-drag is obtained from surface pressure integration and an energy flux integral. Comparisons between computations and experiment show good agreement for the prediction of power polars, surface pressure distribution, and tip vortex geometry.

Ramachandran, K.↗

Parallelized real-time physics codes for plasma control on DIII-D

A real-time safe multi-threading library was developed on the DIII-D plasma control system to optimize the real-time TORBEAM and real-time STRIDE physics codes. These physics codes are crucial for future fusion power plant operation as they provide information about electron cyclotron wave propagation and heating as well as inform about ideal plasma stability limits. The real-time TORBEAM code executed consistently in under 20 ms while the real-time STRIDE code computes in 100 ms. The multi-threading library developed in this work can be applied to other real-time physics-based codes that will be crucial for the next generation of fusion devices.

DIII-D↗

Computational fluid dynamics - The coming revolution

The development of aerodynamic theory is traced from the days of Aristotle to the present, with the next stage in computational fluid dynamics dependent on superspeed computers for flow calculations. Additional attention is given to the history of numerical methods inherent in writing computer codes applicable to viscous and inviscid analyses for complex configurations. The advent of the superconducting Josephson junction is noted to place configurational demands on computer design to avoid limitations imposed by the speed of light, and a Japanese projection of a computer capable of several hundred billion operations/sec is mentioned. The NASA Numerical Aerodynamic Simulator is described, showing capabilities of a billion operations/sec with a memory of 240 million words using existing technology. Near-term advances in fluid dynamics are discussed.

Graves, R. A., Jr.↗

Ray tracing optical analysis of offset solar collector for space station solar dynamic system

OFFSET, a detailed ray tracing computer code, was developed at NASA Lewis Research Center to model the offset solar collector for the Space Station solar dynamic electric power system. This model traces rays from 50 points on the face of the Sun to 10 points on each of the 456 collector facets. The triangular facets are modeled with spherical, parabolic, or toroidal reflective surface contour and surface slope errors. The rays are then traced through the receiver aperture to the walls of the receiver. Images of the collector and of the Sun within the receiver produced by this code provide insight into the collector receiver interface. Flux distribution on the receiver walls, plotted by this code, is improved by a combination of changes to aperture location and receiver tilt angle. Power loss by spillage at the receiver aperture is computed and is considerably reduced by using toroidal facets.

Jefferies, Kent S.↗