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 595 records · Page 33

Numerical Simulation of Supersonic Compression Corners and Hypersonic Inlet Flows Using the RPLUS2D Code

A two-dimensional computational code, PRLUS2D, which was developed for the reactive propulsive flows of ramjets and scramjets, was validated for two-dimensional shock-wave/turbulent-boundary-layer interactions. The problem of compression corners at supersonic speeds was solved using the RPLUS2D code. To validate the RPLUS2D code for hypersonic speeds, it was applied to a realistic hypersonic inlet geometry. Both the Baldwin-Lomax and the Chien two-equation turbulence models were used. Computational results showed that the RPLUS2D code compared very well with experimentally obtained data for supersonic compression corner flows, except in the case of large separated flows resulting from the interactions between the shock wave and turbulent boundary layer. The computational results compared well with the experiment results in a hypersonic NASA P8 inlet case, with the Chien two-equation turbulence model performing better than the Baldwin-Lomax model.

Kapoor, Kamlesh↗

Numerical studies of the deposition of material released from fixed and rotary wing aircraft

The computer code AGDISP (AGricultural DISPersal) 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. In this report, the equations governing the motion of aerially released particles are developed, including a description of the evaporation model used. A series of case studies, using AGDISP, are included.

Bilanin, A. J.↗

Heat Transfer Computations of Internal Duct Flows With Combined Hydraulic and Thermal Developing Length

This study investigated the Navier-Stokes computations of the surface heat transfer coefficients of a transition duct flow. A transition duct from an axisymmetric cross section to a non-axisymmetric cross section, is usually used to connect the turbine exit to the nozzle. As the gas turbine inlet temperature increases, the transition duct is subjected to the high temperature at the gas turbine exit. The transition duct flow has combined development of hydraulic and thermal entry length. The design of the transition duct required accurate surface heat transfer coefficients. The Navier-Stokes computational method could be used to predict the surface heat transfer coefficients of a transition duct flow. The Proteus three-dimensional Navier-Stokes numerical computational code was used in this study. The code was first studied for the computations of the turbulent developing flow properties within a circular duct and a square duct. The code was then used to compute the turbulent flow properties of a transition duct flow. The computational results of the surface pressure, the skin friction factor, and the surface heat transfer coefficient were described and compared with their values obtained from theoretical analyses or experiments. The comparison showed that the Navier-Stokes computation could predict approximately the surface heat transfer coefficients of a transition duct flow.

Wang, C. R.↗

Development, Verification and Use of Gust Modeling in the NASA Computational Fluid Dynamics Code FUN3D

This paper presents the implementation of gust modeling capability in the CFD code FUN3D. The gust capability is verified by computing the response of an airfoil to a sharp edged gust. This result is compared with the theoretical result. The present simulations will be compared with other CFD gust simulations. This paper also serves as a users manual for FUN3D gust analyses using a variety of gust profiles. Finally, the development of an Auto-Regressive Moving-Average (ARMA) reduced order gust model using a gust with a Gaussian profile in the FUN3D code is presented. ARMA simulated results of a sequence of one-minus-cosine gusts is shown to compare well with the same gust profile computed with FUN3D. Proper Orthogonal Decomposition (POD) is combined with the ARMA modeling technique to predict the time varying pressure coefficient increment distribution due to a novel gust profile. The aeroelastic response of a pitch/plunge airfoil to a gust environment is computed with a reduced order model, and compared with a direct simulation of the system in the FUN3D code. The two results are found to agree very well.

Bartels, Robert E.↗

Numerical, Analytical, Experimental Study of Fluid Dynamic Forces in Seals: Executive Summary and Description of Knowledge-Based System - Volume 1

The objectives of the program were to develop computational fluid dynamics (CFD) codes and simpler industrial codes for analyzing and designing advanced seals for air-breathing and space propulsion engines. The CFD code SCISEAL is capable of producing full three-dimensional flow field information for a variety of cylindrical configurations. An implicit multidomain capability allows the division of complex flow domains to allow optimum use of computational cells. SCISEAL also has the unique capability to produce cross-coupled stiffness and damping coefficients for rotordynamic computations. The industrial codes consist of a series of separate stand-alone modules designed for expeditious parametric analyses and optimization of a wide variety of cylindrical and face seals. Coupled through a Knowledge-Based System (KBS) that provides a user-friendly Graphical User Interface (GUI), the industrial codes are PC based using an OS/2 operating system. These codes were designed to treat film seals where a clearance exists between the rotating and stationary components. Leakage is inhibited by surface roughness, small but stiff clearance films, and viscous pumping devices. The codes have demonstrated to be a valuable resource for seal development of future air-breathing and space propulsion engines.

Liang, Anita D.↗

XSECT: A computer code for generating fuselage cross sections - user's manual

A computer code, XSECT, has been developed to generate fuselage cross sections from a given area distribution and wing definition. The cross sections are generated to match the wing definition while conforming to the area requirement. An iterative procedure is used to generate each cross section. Fuselage area balancing may be included in this procedure if desired. The code is intended as an aid for engineers who must first design a wing under certain aerodynamic constraints and then design a fuselage for the wing such that the contraints remain satisfied. This report contains the information necessary for accessing and executing the code, which is written in FORTRAN to execute on the Cyber 170 series computers (NOS operating system) and produces graphical output for a Tektronix 4014 CRT. The LRC graphics software is used in combination with the interface between this software and the PLOT 10 software.

Ames, K. R.↗

Onboard magnetic field modeling for Solar Maximum Mission /SMM/

Analysis and simulation results are presented for magnetic field models for use in attitude acquisition onboard Solar Maximum Mission (SMM). A study was made of the degree of the spherical harmonic expansion of the magnetic field required, considering mission requirements, modeling errors, and magnetometer quantization and biases. It is shown that a fifth-degree field is sufficient to provide two-degree roll angle determination accuracy with a residual magnetic bias of 10 milligauss. Also, a spherical harmonic expansion for the McIlwain L-parameter is included for the first time. This parameter will be telemetered to ground with experimental data. The fifth-degree expansion will provide the L-parameter to within two percent of accepted values. The additional onboard computational burden is the storage of 36 coefficients and an increase of about 15% in computation time. Prototype flight code was developed which is anticipated to require about 2000 bytes of core storage and 30 milliseconds of computation time per orbit point on the NSSC-1 computer.

Headrick, R. D.↗

Converging-Diverging Nozzle Thruster Code for Nuclear or Chemical Rocket Performance Computations

This report describes development of a variable area ratio C-D nozzle numerical computational code "THRST", applicable to the evaluation of nuclear or chemical rocket performance. The propellant for the nuclear rockets is para-hydrogen with gas properties supplied by a sub-routine, while properties for chemical thrusters need to be user input. Detailed documentation of input and output variables and sample output listings enable the report to function as a "User Manual".

Juhasz, Albert J.↗

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 Modelling Program. This will be accomplished by reviewing existing methodologies for thermal data transfer and selecting three heat transfer codes for application in this program, 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, developing the automated thermal load transfer module, and verifying and documenting the module. The overall objectives of this thermal transfer module are that it handle independent mesh configurations, perform the transfer in an accurate and efficient fashion and that the total system be flexible for future improvements.

Maffeo, R. J.↗

Experimental Verification of Computational Models for Laminated Composites

The objective of the research reported herein is to develop a progressive damage methodology capable of predicting the residual strength of continuous fiber-reinforced, laminated, polymer matrix composites with through-penetration damage. The fracture behavior of center-notch tension panels with thin crack-like slits was studied. Since fibers are the major load-carrying constituent in polymer matrix composites, predicting the residual strength of a laminate requires a criterion for fiber fracture. The effects on fiber strain due to other damage mechanisms such as matrix cracking and delaminations must also be modeled. Therefore, the research herein examines the damage mechanisms involved in translaminate fracture and identifies the toughening mechanisms responsible for damage growth resistance in brittle epoxy matrix systems. The mechanics of matrix cracking and fiber fracture are discussed as is the mathematical framework for the progressive damage model developed by the authors. The progressive damage analysis algorithms have been implemented into a general purpose finite element code developed by NASA, the Computational Structural Mechanics Testbed (COMET). Damage growth is numerically simulated and the analytical residual strength predictions are compared to experimental results for a variety of notched panel configurations and materials systems.

Harris, Charles E.↗

Prediction of sound radiation from different practical jet engine inlets

The computer codes necessary for this study were developed and checked against exact solutions generated by the point source method using the NASA Lewis QCSEE inlet geometry. These computer codes were used to predict the acoustic properties of the following five inlet configurations: the NASA Langley Bellmouth, the NASA Lewis JT15D-1 Ground Test Nacelle, and three finite hyperbolic inlets of 50, 70 and 90 degrees. Thirty-five computer runs were done for the NASA Langley Bellmouth. For each of these computer runs, the reflection coefficient at the duct exit plane was calculated as was the far field radiation pattern. These results are presented in both graphical and tabular form with many of the results cross plotted so that trends in the results verses cut-off ratio (wave number) and tangential mode number may be easily identified.

Zinn, B. T.↗

Combustor-Turbine Interactions By Using the Open National Combustion Code (OpenNCC) and the Glenn-HT Code

We develop the methodology of coupling two different computational fluid dynamic (CFD) codes, OpenNCC and Glenn-HT, in order to investigate the unsteady flow fields inside the combustor and around the first-stage stator of a high-pressure turbine (HPT) from the energy efficient engine (E 3 ) program. In our coupling strategy, OpenNCC, a time accurate unstructured mesh multi-phase combustion CFD code, is applied to the combustor region and Glenn-HT, a structured mesh CFD code designed for turbine applications, is utilized inside the HPT region. As a proof of concept, we first simulate the three-dimensional airflow over the backward facing step and compared the predicted pressure coefficient against experimental data. We then considered a Jet-A/Air combustor with the NASA single learn direct injector (LDI) as a reacting flow test case. Here, Glenn-HT is used in the downstream of the combustor and assumes that the working fluid is a single gas (i.e., air) where the properties of the gas are taken from time- and area-averaged values at the combustor exit. A fully coupled combustor-turbine simulation using both OpenNCC and Glenn-HT was conducted on the E 3 combustor and a detailed analysis done pertaining to how a realistic combustor outflow affects migration of hot-streaks and its impact on the aerodynamic behavior inside the HPT. It is found that the combustion dynamics (i.e., the switching main and pilot flame strength) significantly alters the aerodynamic behavior inside the HPT including hot-streak migration. Variations of the temperature and the velocity magnitude in the passage could be up to ±75 and ±15. These large variations are an important part of the combustor-turbine interactions, which are overlooked by a single component simulation of the HPT while imposing the steady inflow boundary condition at the HPT inlet. Finally, it is observed the strong density gradient associated with the hot-streaks and the pressure gradient at the passage significantly enhances the baroclinic torque, affecting the vorticity field.

CFD hot-streaks↗

Computational simulation of transition to turbulence through inverse modeling

The present investigation has focused on a computational methodology for the fundamental case of transition in channel flow, in which recently published experimental data are utilized both as a stimulus and as a measure of merit of the method. The research has proceeded along three avenues in parallel. The first task has consisted of the development and verification of a computer code which calculates the mean evolution of flow in a channel similar to the one employed experimentally by Blair and Anderson. An analytical test case was created for the dual purposes of code verification and of highlighting the interactions between the Reynolds stress and the mean velocity profile. This test case generated a Reynolds stress by the residue in the momentum equation which is produced by a typical analytical velocity profile. By a substitution of this Reynolds stress into the appropriate code module, the correctness of the code may be verified, along with the accuracy of the computational method. The second task pursued has involved the development of a triple layer model for the Reynolds stress profile, which was suggested and derived from experimental velocity profiles. It is demonstrated that the innermost length scale is based on the local friction velocity, the intermediate layer corresponds to the usual logarithmic law of the wall region in which the normalized Reynolds stress is approximately unity, and the outermost layer is represented by a closed mathematical form depending explicitly on the velocity profile in the wake region. The third task was comprised of scrutiny of the excellent databases developed by Blair and others, and the planning of its incorporation into the transition analysis. These extensive measurements indicate that turbulent statistics in the transition regime may be considered to alternate between laminar and fully turbulent types, the proportions of which are quantified by a measured intermittency function.

Sepri, Paavo↗

JAE: A Jupiter Atmospheric Entry Probe Heating Code

The strong gravitational attraction of Jupiter on probes approaching the planet results in very high atmospheric entry velocities. The values relative to the rotating atmosphere can vary from about 47 to 60 km/sec, depending on the latitude of the entry. Therefore, the peak heating rates and heat shield mass fractions exceed those for any other atmospheric entries. For example, the Galileo probe's heat shield mass fraction was 50%, of which 45% was devoted to the forebody. Although the Galileo probe's mission was very successful, many more scientific questions about the Jovian atmosphere remain to be answered and additional probe missions are being planned. Recent developments in microelectronics have raised the possibility of building smaller and less expensive probes than Galileo. Therefore, it was desirable to develop a code that could quickly compute the forebody entry heating environments when performing parametric probe sizing studies. The Jupiter Atmospheric Entry (JAE) code was developed to meet this requirement. The body geometry consists of a blunt-nosed conical shape of arbitrary nose and base radius and cone angles up to about 65 deg at zero angle of attack.

Wercinski, Paul F.↗

Propulsion stability codes for liquid propellant propulsion systems developed for use on a PC computer

Research into component modeling and system synthesis leading to the analysis of the major types of propulsion system instabilities and the characterization of various components characteristics are presented. Last year, several programs designed to run on a PC were developed for Marshall Space Flight Center. These codes covered the low, intermediate, and high frequency modes of oscillation of a liquid rocket propulsion system. No graphics were built into these programs and only simple piping layouts were supported. This year's effort was to add run time graphics to the low and intermediate frequency codes, allow new types of piping elements (accumulators, pumps, and split pipes) in the low frequency code, and develop a new code for the PC to generate Nyquist plots.

Doane, George B., III↗

A benchmark for galactic cosmic ray transport codes

A nontrivial analytic benchmark solution for galactic cosmic ray transport is presented for use in transport code validation. Computational accuracy for a previously-developed cosmic ray transport code is established to within one percent by comparison with this exact benchmark. Hence, solution accuracy for the transport problem is mainly limited by inaccuracies in the input spectra, input interaction databases, and the use of a straight ahead/velocity-conserving approximation.

Wilson, John W.↗

A benchmark for galactic cosmic-ray transport codes

A nontrivial analytic benchmark solution for galactic cosmic ray transport is presented for use in transport code validation. Computational accuracy for a previously-developed cosmic ray transport code is established to within one percent by comparison with this exact benchmark. Hence, solution accuracy for the transport problem is mainly limited by inaccuracies in the input spectra, input interaction databases, and the use of a straight ahead/velocity-conserving approximation.

Wilson, John W.↗

Guidelines for development structured FORTRAN programs

Computer programming and coding standards were compiled to serve as guidelines for the uniform writing of FORTRAN 77 programs at NASA Langley. Software development philosophy, documentation, general coding conventions, and specific FORTRAN coding constraints are discussed.

Earnest, B. M.↗