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.

SourceResearch connectionFrom yearThrough year

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

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Water droplet impingement on airfoils and aircraft engine inlets for icing analysis

This paper includes the results of a significant research program for verification of computer trajectory codes used in aircraft icing analysis. Experimental water droplet impingement data have been obtained in the NASA Lewis Research Center Icing Research Tunnel for a wide range of aircraft geometries and test conditions. The body whose impingement characteristics are required is covered at strategic locations by thin strips of moisture absorbing (blotter) paper and then exposed to an airstream containing a dyed-water spray cloud. Water droplet impingement data are extracted from the dyed blotter strips by measuring the optical reflectance of the dye deposit on the strips with an automated reflectometer. Impingement characteristics for all test geometries have also been calculated using two recently developed trajectory computer codes. Good agreement is obtained with experimental data. The experimental and analytical data show that maximum impingement efficiency and impingement limits increase with mean volumetric diameter for all geometries tested. For all inlet geometries tested, as the inlet mass flow is reduced, the maximum impingement efficiency is reduced and the location of the maximum impingement shifts toward the inlet inner cowl.

Papadakis, Michael↗

NASA NTRS2013

Prediction of Ablation Rates from Solid Surfaces Exposed to High Temperature Gas Flow

A mathematical model and a solution algorithm is developed to study the physics of high temperature heat transfer and material ablation and identify the problems associated with the flow of hydrogen gas at very high temperatures and velocities through pipes and various components of Nuclear Thermal Rocket (NTR) motors. Ablation and melting can be experienced when the inner solid surface of the cooling channels and the diverging-converging nozzle of a Nuclear Thermal Rocket (NTR) motor is exposed to hydrogen gas flow at temperatures around 2500 degrees Kelvin and pressures around 3.4 MPa. In the experiments conducted on typical NTR motors developed in 1960s, degradation of the cooling channel material (cracking in the nuclear fuel element cladding) and in some instances melting of the core was observed. This paper presents the results of a preliminary study based on two types of physics based mathematical models that were developed to simulate the thermal-hydrodynamic conditions that lead to ablation of the solid surface of a stainless steel pipe exposed to high temperature hydrogen gas near sonic velocities. One of the proposed models is one-dimensional and assumes the gas flow to be unsteady, compressible and viscous. An in-house computer code was developed to solve the conservations equations of this model using a second-order accurate finite-difference technique. The second model assumes the flow to be three-dimensional, unsteady, compressible and viscous. A commercial CFD code (Fluent) was used to solve the later model equations. Both models assume the thermodynamic and transport properties of the hydrogen gas to be temperature dependent. In the solution algorithm developed for this study, the unsteady temperature of the pipe is determined from the heat equation for the solid. The solid-gas interface temperature is determined from an energy balance at the interface which includes heat transfer from or to the interface by conduction, convection, radiation, and ablation. Two different ablation models are proposed to determine the heat loss from the solid surface due to the ablation of the solid material. Both of them are physics based. Various numerical simulations were carried out using both models to predict the temperature distribution in the solid and in the gas flow, and then predict the ablation rates at a typical NTR motor hydrogen gas temperature and pressure. Solid mass loss rate per foot of a pipe was also calculated from these predictions. The results are presented for fully developed turbulent flow conditions in a sample SS pipe with a 6 inch diameter.

Akyuzlu, Kazim M.↗

NASA NTRS1982

Flow aerodynamics modeling of an MHD swirl combustor - Calculations and experimental verification

This paper describes a computer code for calculating the flow dynamics of constant density flow in the second stage trumpet shaped nozzle section of a two stage MHD swirl combustor for application to a disk generator. The primitive pressure-velocity variable, finite difference computer code has been developed to allow the computation of inert nonreacting turbulent swirling flows in an axisymmetric MHD model swirl combustor. The method and program involve a staggered grid system for axial and radial velocities, and a line relaxation technique for efficient solution of the equations. Turbulence simulation is by way of a two-equation Kappa-epsilon model. The code produces as output the flowfield map of the nondimensional stream function, axial, and swirl velocity. Good agreement was obtained between the theoretical predictions and the qualitative experimental results. The best seed injector location for uniform seed distribution at combustor exit is with injector located centrally on the combustor axis at entrance to the second stage combustor.

Gupta, A. K.↗

NASA NTRS1993

CNSFV code development, virtual zone Navier-Stokes computations of oscillating control surfaces and computational support of the laminar flow supersonic wind tunnel

The work performed during the past year on this cooperative agreement covered two major areas and two lesser ones. The two major items included further development and validation of the Compressible Navier-Stokes Finite Volume (CNSFV) code and providing computational support for the Laminar Flow Supersonic Wind Tunnel (LFSWT). The two lesser items involve a Navier-Stokes simulation of an oscillating control surface at transonic speeds and improving the basic algorithm used in the CNSFV code for faster convergence rates and more robustness. The work done in all four areas is in support of the High Speed Research Program at NASA Ames Research Center.

Klopfer, Goetz H.↗

NASA NTRS1988

Optimization of large matrix calculations for execution on the Cray X-MP vector supercomputer

A considerable volume of large computational computer codes were developed for NASA over the past twenty-five years. This code represents algorithms developed for machines of earlier generation. With the emergence of the vector supercomputer as a viable, commercially available machine, an opportunity exists to evaluate optimization strategies to improve the efficiency of existing software. This result is primarily due to architectural differences in the latest generation of large-scale machines and the earlier, mostly uniprocessor, machines. A sofware package being used by NASA to perform computations on large matrices is described, and a strategy for conversion to the Cray X-MP vector supercomputer is also described.

Hornfeck, William A.↗

NASA NTRS2011

Modeling the Effects of Ice Accretion on the Low Pressure Compressor and the Overall Turbofan Engine System Performance

The focus of this study is on utilizing a mean line compressor flow analysis code coupled to an engine system thermodynamic code, to estimate the effects of ice accretion on the low pressure compressor, and quantifying its effects on the engine system throughout a notional flight trajectory. In this paper a temperature range in which engine icing would occur was assumed. This provided a mechanism to locate potential component icing sites and allow the computational tools to add blockages due to ice accretion in a parametric fashion. Ultimately the location and level of blockage due to icing would be provided by an ice accretion code. To proceed, an engine system modeling code and a mean line compressor flow analysis code were utilized to calculate the flow conditions in the fan-core and low pressure compressor and to identify potential locations within the compressor where ice may accrete. In this study, an "additional blockage" due to the accretion of ice on the metal surfaces, has been added to the baseline aerodynamic blockage due to boundary layer, as well as the blade metal blockage. Once the potential locations of ice accretion are identified, the levels of additional blockage due to accretion were parametrically varied to estimate the effects on the low pressure compressor blade row performance operating within the engine system environment. This study includes detailed analysis of compressor and engine performance during cruise and descent operating conditions at several altitudes within the notional flight trajectory. The purpose of this effort is to develop the computer codes to provide a predictive capability to forecast the onset of engine icing events, such that they could ultimately help in the avoidance of these events.

Veres, Joseph P.↗

NASA NTRS1986

Radiation transport in dust in disk geometry

The main objective of the research program is twofold: (1) to develop a computer code to solve the problem of scattering, absorption and emission of photons by dust grains in a dusty medium with 2 dimensional disk geometry, and (2) to study the various physical and geometrical effects of 2 dimensional radiation transport on the thermal structure and radiation field. These tasks were accomplished and are briefly summarized. The method for solving the radiation transport problem in disk geometry is a generalization of the quasi-diffusion method (QDM) previously developed by the author.

Chun, Ming Leung↗

NASA NTRS1993

Computations of spray, fuel-air mixing, and combustion in a lean-premixed-prevaporized combustor

A code was developed for computing the multidimensional flow, spray, combustion, and pollutant formation inside gas turbine combustors. The code developed is based on a Lagrangian-Eulerian formulation and utilizes an implicit finite-volume method. The focus of this paper is on the spray part of the code (both formulation and algorithm), and a number of issues related to the computation of sprays and fuel-air mixing in a lean-premixed-prevaporized combustor. The issues addressed include: (1) how grid spacings affect the diffusion of evaporated fuel, and (2) how spurious modes can arise through modelling of the spray in the Lagrangian computations. An upwind interpolation scheme is proposed to account for some effects of grid spacing on the artificial diffusion of the evaporated fuel. Also, some guidelines are presented to minimize errors associated with the spurious modes.

Dasgupta, A.↗

NASA NTRS1986

Variable-Conductance Heat Pipes

In response to need to accurately and efficiently predict performance of variable-conductance heat pipes (VCHP's) incorporated in spacecraft thermalcontrol systems, computer code VCHPDA developed to interact with thermal analyzer programs such as SINDA (Systems Improved Numerical Differencing Analyzer). Calculates length of gas-blocked region and vapor temperature in active portion. Advantages of VCHPDA over prior programs improved accuracy, unconditional stability, and increased efficiency of solution resulting from novel approach and use of state-of-the-art numerical techniques for solving VCHP mathematical model. Code valuable tool in design and evaluation of advanced thermal-control systems using variable-conductance heat pipes. Written in FORTRAN IV for use on CDC 600 computers.

Antoniuk, D.↗

NASA NTRS1991

Heat pipe radiation cooling of advanced hypersonic propulsion system components

Heat transfer, heat pipe, and system studies were performed to assess the newly proposed heat pipe radiation cooling (HPRC) concept. With an HPRC system, heat is removed from the ramburner and nozzle of a hypersonic aircraft engine by a surrounding, high-temperature, heat pipe nacelle structure, transported to nearby external surfaces, and rejected to the environment by thermal radiation. With HPRC, the Mach number range available for using hydrocarbon fuels for aircraft operation extends into the Mach 4 to Mach 6 range, up from the current limit of about Mach 4. Heat transfer studies using a newly developed HPRC computer code determine cooling system and ramburner and nozzle temperatures, heat loads, and weights for a representative combined-cycle engine cruising at Mach 5 at 80,000 ft altitude. Heat pipe heat transport calculations, using the Los Alamos code HTPIPE, reveal that adequate heat trasport capability is available using molybdenum-lithium heat pipe technology. Results show that the HPRC system radiator area is limited in size to the ramburner-nozzle region of the engine nacelle; reasonable system weights are expected; hot section temperatures are consistent with advanced structural materials development goals; and system impact on engine performance is minimal.

Martin, R. A.↗

NASA NTRS1988

Computation of the tip vortex flowfield for advanced aircraft propellers

The tip vortex flowfield plays a significant role in the performance of advanced aircraft propellers. The flowfield in the tip region is complex, three-dimensional and viscous with large secondary velocities. An analysis is presented using an approximate set of equations which contains the physics required by the tip vortex flowfield, but which does not require the resources of the full Navier-Stokes equations. A computer code was developed to predict the tip vortex flowfield of advanced aircraft propellers. A grid generation package was developed to allow specification of a variety of advanced aircraft propeller shapes. Calculations of the tip vortex generation on an SR3 type blade at high Reynolds numbers were made using this code and a parametric study was performed to show the effect of tip thickness on tip vortex intensity. In addition, calculations of the tip vortex generation on a NACA 0012 type blade were made, including the flowfield downstream of the blade trailing edge. Comparison of flowfield calculations with experimental data from an F4 blade was made. A user's manual was also prepared for the computer code (NASA CR-182178).

Tsai, Tommy M.↗

NASA NTRS1990

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 1 is the Analysis Description, and describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Towne, Charles E.↗

NASA NTRS1990

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 2: User's guide

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 2 is the User's Guide, and describes the program's general features, the input and output, the procedure for setting up initial conditions, the computer resource requirements, the diagnostic messages that may be generated, the job control language used to run the program, and several test cases.

Towne, Charles E.↗

NASA NTRS1990

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

A new computer code was developed to solve the 2-D or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating-direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 3 is the Programmer's Reference, and describes the program structure, the FORTRAN variables stored in common blocks, and the details of each subprogram.

Towne, Charles E.↗

NASA NTRS1984

STAEBL: Structural tailoring of engine blades, phase 2

The Structural Tailoring of Engine Blades (STAEBL) program was initiated at NASA Lewis Research Center in 1980 to introduce optimal structural tailoring into the design process for aircraft gas turbine engine blades. The standard procedure for blade design is highly iterative with the engineer directly providing most of the decisions that control the design process. The goal of the STAEBL program has been to develop an automated approach to generate structurally optimal blade designs. The program has evolved as a three-phase effort with the developmental work being performed contractually by Pratt & Whitney Aircraft. Phase 1 was intended as a proof of concept in which two fan blades were structurally tailored to meet a full set of structural design constraints while minimizing DOC+I (direct operating cost plus interest) for a representative aircraft. This phase was successfully completed and was reported in reference 1 and 2. Phase 2 has recently been completed and is the basis for this discussion. During this phase, three tasks were accomplished: (1) a nonproprietary structural tailoring computer code was developed; (2) a dedicated approximate finite-element analysis was developed; and (3) an approximate large-deflection analysis was developed to assess local foreign object damage. Phase 3 is just beginning and is designed to incorporated aerodynamic analyses directly into the structural tailoring system in order to relax current geometric constraints.

Hirschbein, M. S.↗

NASA NTRS1986

Numerical studies of motion and decay of vortex filaments

A computational code is developed for the integro-differential equations governing the motion of the centerlines of vortex filaments submerged in a background potential flow. These equations, which are derived from the method of matched asymptotic analysis, include the effect of decaying large-magnitude circumferential and axial velocity components in the vortical cores. Numerical examples are presented to assess the effect of large axial velocity and of nonsimilar initial profiles in vortical cores. The initial configurations of the filaments are chosen so as to fulfill the basic assumption of asymptotic analysis, which is the effective vortical core size is much smaller than all other length scales in the flowfield, e.g., the radius of curvature and interfilament distance. The computations are continued until the basic assumption is no longer valid, that is, when the merging or intersection of filaments have begun. Various types of local or global merging or intersection of filaments are classified and demonstrated by numerical examples.

Liu, C. H.↗

NASA NTRS1982

A simple method to compute the change in earth-atmosphere radiative balance due to a stratospheric aerosol layer

A computer code was developed in terms of a three-layer model for the earth-atmosphere system, using a two-stream approximation for the troposphere and stratosphere. The analysis was limited to variable atmosphere loading by solar radiation over an unperturbed section of the atmosphere. The scattering atmosphere above a Lambertian ground layer was considered in order to derive the planar albedo and the spherical albedo. Attention was given to the influence of the aerosol optical thickness in the stratosphere, the single scattering albedo and asymmetry factor, and the sublayer albedo. Calculations were performed of the zonal albedo and the planetary radiation balance, taking into account a stratospheric aerosol layer containing H2SO4 droplets and volcanic ash. The resulting ground temperature disturbance was computed using a Budyko (1969) climate model. Local decreases in the albedo in the summer were observed in high latitudes, implying a heating effect of the aerosol. An accompanying energy loss of 23-27 W/sq m was projected, which translates to surface temperature decreases of either 1.1 and 0.45 C, respectively, for background and volcanic aerosols.

Lenoble, J.↗

NASA NTRS1984

Numerical studies of motion of vortex filaments - Implementing the asymptotic analysis

A computational code is developed for the integro-differential equations governing the motion of the centerlines of vortex filaments submerged in a background potential flow. These equations, which are derived from the method of matched asymptotic analysis, include the effect of the decaying large-magnitude circumferential and axial velocity components in the vortical cores. Numerical examples are presented to assess the effect of a large axial velocity and that of nonsimilar initial profiles in the vortical cores. The initial configurations of the filaments are chosen so as to fulfill the basic assumption of the asymptotic analysis, which is that the effective vortical core size is much smaller than all the other length scales in the flowfield, e.g., the radius of curvature and the interfilament distance. The computations are continued until the basic assumption is no longer valid, that is when the merging or intersection of filaments has begun. A classification of the various types of local or global merging or intersection of filaments is made and demonstrated by numerical examples. It is then shown that the asymptotic solution not only provides the initial data but also can be used to formulate the appropriate boundary conditions for the numerical solution of a merged region.

Liu, C. H.↗