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 415 records · Page 23

Design, Generation and Tooth Contact Analysis (TCA) of Asymmetric Face Gear Drive With Modified Geometry

A new type of face gear drive for application in transmissions, particularly in helicopters, has been developed. The new geometry differs from the existing geometry by application of asymmetric profiles and double-crowned pinion of the face gear mesh. The paper describes the computerized design, simulation of meshing and contact, and stress analysis by finite element method. Special purpose computer codes have been developed to conduct the analysis. The analysis of this new type of face gear is illustrated with a numerical example.

Litvin, Faydor L.↗

A 3D-PNS computer code for the calculation of supersonic combusting flows

A computer code has been developed based on the three-dimensional parabolized Navier-Stokes (PNS) equations which govern the supersonic combusting flow of the hydrogen-air system. The finite difference algorithm employed was a hybrid of the Schiff-Steger algorithm and the Vigneron, et al., algorithm which is fully implicit and fully coupled. The combustion of hydrogen and air was modeled by the finite-rate two-step combustion model of Rogers-Chinitz. A new dependent variable vector was introduced to simplify the numerical algorithm. Robustness of the algorithm was considerably enhanced by introducing an adjustable parameter. The computer code was used to solve a premixed shock-induced combustion problem and the results were compared with those of a full Navier-Stokes code. Reasonably good agreement was obtained at a fraction of the cost of the full Navier-Stokes procedure.

Chitsomboon, Tawit↗

Wake curvature and trailing edge interaction effects in viscous flow over airfoils

A theory developed for analyzing viscous flows over airfoils at high Reynolds numbers is described. The theory includes a complete treatment of viscous interaction effects induced by the curved wake behind the airfoil and accounts for normal pressure gradients across the boundary layer in the trailing edge region. A brief description of a computer code that was developed to solve the extended viscous interaction equations is given. Comparisons of the theoretical results with wind tunnel data for two rear loaded airfoils at supercritical conditions are presented.

Melnik, R. E.↗

A comparison of computed Space Shuttle Orbiter surface pressures with flight measurements

Vehicles such as the Shuttle Orbiter operate over a wide range of angle of attack when entering the earth's atmosphere. This angle-of-attack range can be divided into two classes, a moderate regime, and a high regime for angles of attack greater than 30 deg. The HALIS (High Alpha Inviscid Solution) code was developed for computing flowfields in the high angle-of-attack regime. The present investigation is concerned with conducting a comparison of HALIS results with wind-tunnel data and the limited flight data from STS-1. The investigation shows that at high angle-of-attack, the axial flow over the entire windward surface can be subsonic, and, for those conditions, where the subsonic axial flow is confined to the nose of the vehicle, embedded regions of subsonic axial flow may be found along the vehicle strake tips.

Weilmuenster, K. J.↗

Influence of bulk turbulence and entrance boundary layer thickness on the curved duct flow field

The objective of this investigation was the experimental evaluation of bulk turbulence and boundary thickness influence on the secondary flow development in a square, 90 deg turning duct. A three dimensional laser velocimetry system was utilized to measure the mean and fluctuating components of velocity in the large curved duct facility. The three dimensional development of the viscous shear layers in the curved duct has a strong influence on the complete flow field. Since ducted three dimensional flows are found in many engineering applications, including gas turbine engines, and contain high turbulence levels and high wall heat transfer rates, they present a difficult challenge to computational fluid mechanics codes. Turbulence modeling remains one of constraints to CFD advance due to inadequate physical understanding and experimental definition of turbulent shear flows. The results of this investigation expand the curved duct data base to higher turbulence levels and thicker entrance boundary layers. The experimental results provide a challenging benchmark data base for computational fluid dynamics code development and validation. The variation of inlet bulk turbulence intensity provides additional information to aid in turbulence model evaluation.

Crawford, Roger A.↗

Influence of bulk turbulence and entrance boundary layer thickness on the curved duct flow field

The objective of this investigation was the experimental evaluation of bulk turbulence and boundary thickness influence on the secondary flow development in a square, 90 deg turning duct. A three dimensional laser velocimetry system was utilized to measure the mean and fluctuating components of velocity in the large curved duct facility. The three dimensional development of the viscous shear layers in the curved duct has a strong influence on the complete flow field. Since ducted three dimensional flows are found in many engineering applications, including gas turbine engines, and contain high turbulence levels and high wall heat transfer rates, they present a difficult challenge to computational fluid mechanics codes. Turbulence modeling remains one of constraints to CFD advance due to inadequate physical understanding and experimental definition of turbulent shear flows. The results of this investigation expand the curved duct data base to higher turbulence levels and thicker entrance boundary layers. The experimental results provide a challenging benchmark data base for computational fluid dynamics code development and validation. The variation of inlet bulk turbulence intensity provides additional information to aid in turbulence model evaluation.

Shohadaee, Ahmad A.↗

Transonic Airfoil Development

This lecture consists of three parts, in which discussions are presented of the current state of development of transonic or supercritical airfoils designed for fully turbulent boundary layers on the surfaces, previous research on subcritical airfoils designed to achieve laminar boundary layers on all or parts of the surfaces, and current research on supercritical airfoils designed to achieve laminar boundary layers. In the first part the use of available two dimensional computer codes in the development of supercritical airfoils and the general trends in the design of such airfoils with turbulent boundary layers are discussed. The second part provides the necessary background on laminar boundary layer phenomena. The last part, which constitutes the major portion of the lecture, covers research by NASA on supercritical airfoils utilizing both decreasing pressure gradients and surface suction for stabilizing the laminar boundary layer. An investigation of the former has been recently conducted in fight using gloves on the wing panels of the U.S. Air Force F111 TACT airplane, research on the later is currently being conducted in a transonic wind tunnel which has been modified to greatly reduce the stream turbulence and noise levels in the tests section.

Whitcomb, R. T.↗

Seals Flow Code Development

In recognition of a deficiency in the current modeling capability for seals, an effort was established by NASA to develop verified computational fluid dynamic concepts, codes, and analyses for seals. The objectives were to develop advanced concepts for the design and analysis of seals, to effectively disseminate the information to potential users by way of annual workshops, and to provide experimental verification for the models and codes under a wide range of operating conditions.

Source record↗

On the generation of secondary motion in circular to rectangular transition ducts

A computer code has been developed for the calculation of three-dimensional turbulent flow in circular to rectangular transition ducts. The governing equations are written in a general curvilinear coordinate system so that all singly-connected geometrical configurations can be accommodated. Turbulence closure is through the standard k-epsilon model or a Reynolds stress closure model. Computed results show qualitative agreement with experimental data. Pressure-induced secondary flow is the dominant mode through the transition, but turbulence-induced secondary flow becomes prominent some distance downstream of the transition.

Demuren, A. O.↗

Component aerodynamic characteristics of banked cruciform missiles with arbitrary control deflections

A comprehensive computer code has been developed for the prediction of the static aerodynamic characteristics of cruciform wing-body combinations with one or two sets of fins. The code may be used for subsonic, transonic, and supersonic speeds. Arbitrary bank angles and control deflections are allowed. The present version of the code may be used for angles of attack up to 20 deg. The code computes the forces and moments acting on each component of the missile including the effects of vortices shed by the body and fins. Comparisons between experiment and theory indicate good agreement for longitudinal characteristics and fair to good agreement for lateral characteristics.

Hemsch, M. J.↗

Shielding materials for highly penetrating space radiations

Interplanetary travel involves the transfer from an Earth orbit to a solar orbit. Once outside the Earth's magnetosphere, the major sources of particulate radiation are solar cosmic rays (SCR's) and galactic cosmic rays (GCR's). Intense fluxes of SCR's come from solar flares and consist primarily of protons with energies up to 1 GeV. The GCR consists of a low flux of nuclei with energies up to 10(exp 10) GeV. About 70 percent of the GCR are protons, but a small amount (0.6 percent) are nuclei with atomic numbers greater than 10. High energy charged particles (HZE) interact with matter by transferring energy to atomic electrons in a Coulomb process and by reacting with an atomic nucleus. Energy transferred in the first process increases with the square of the atomic number, so particles with high atomic numbers would be expected to lose large amounts of energy by this process. Nuclear reactions produced by (HZE) particles produce high-energy secondary particles which in turn lose energy to the material. The HZE nuclei are a major concern for radiation protection of humans during interplanetary missions because of the very high specific ionization of both primary and secondary particles. Computer codes have been developed to calculate the deposition of energy by very energetic charged particles in various materials. Calculations show that there is a significant buildup of secondary particles from nuclear fragmentation and Coulomb dissociation processes. A large portion of these particles are neutrons. Since neutrons carry no charge, they only lose energy by collision or reaction with a nucleus. Neutrons with high energies transfer large amounts of energy by inelastic collisions with nuclei. However, as the neutron energy decreases, elastic collisions become much more effective for energy loss. The lighter the nucleus, the greater the fraction of the neutron's kinetic energy that can be lost in an elastic collision. Thus, hydrogen-containing materials such as polymers are most effective in reducing the energy of neutrons. Once neutrons are reduced to very low energies, the probability for undergoing a reaction with a nucleus (the cross section) becomes very high. The product of such a reaction is often radioactive and can involve the release of a significant amount of energy. Thus, it is important to provide protection from low energy neutrons during a long duration space flight. Among the light elements, lithium and boron each have an isotope with a large thermal neutron capture cross section, Li-6 and B-10. However, B-10 is more abundant in the naturally-occurring element than Li-6, has a thermal neutron capture cross section four times that of Li-6, and produces the stable products, He-4 and Li-7 in the interaction while Li-6 produces radioactive tritium (H-3). Thus, boron is the best light-weight material for thermal neutron absorption in spacecraft. The work on this project was focused in two areas: computer design where existing computer codes were used, and in some cases modified, to calculate the propagation and interactions of high energy charged particles through various media, and materials development where boron was incorporated into high performance materials.

Kiefer, Richard L.↗

Analysis and Characterization of Damage Utilizing an Orthotropic Generalized Composite Material Model Suitable for Use in Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive communities. In order to address a series of issues identified by the aerospace community as being desirable to include in a next generation composite impact model, an orthotropic, macroscopic constitutive model incorporating both plasticity and damage suitable for implementation within the commercial LS-DYNA computer code is being developed. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain hardening-based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used in which a load in one direction results in a stiffness reduction in multiple material coordinate directions. A detailed analysis is carried out to ensure that the strain equivalence assumption is appropriate for the derived plasticity and damage formulations that are employed in the current model. Procedures to develop the appropriate input curves for the damage model are presented and the process required to develop an appropriate characterization test matrix is discussed.

Impact↗

Analysis and Characterization of Damage Utilizing an Orthotropic Generalized Composite Material Model Suitable for Use in Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive communities. In order to address a series of issues identified by the aerospace community as being desirable to include in a next generation composite impact model, an orthotropic, macroscopic constitutive model incorporating both plasticity and damage suitable for implementation within the commercial LS-DYNA computer code is being developed. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain hardening-based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used in which a load in one direction results in a stiffness reduction in multiple material coordinate directions. A detailed analysis is carried out to ensure that the strain equivalence assumption is appropriate for the derived plasticity and damage formulations that are employed in the current model. Procedures to develop the appropriate input curves for the damage model are presented and the process required to develop an appropriate characterization test matrix is discussed

Impact↗

Three-dimensional inviscid analysis of the scramjet inlet flow field

A computer code has been developed to analyze the inviscid flow field in a supersonic combustion ramjet (scramjet) inlet. The code uses the three-dimensional Euler equations in full conservation form to describe the inlet flow. An algebraic numerical coordinate transformation is used to generate a set of boundary-fitted curvilinear coordinates. The governing equations are solved by a time-asymptotic, unsplit, two-step, finite-difference method. This method is highly efficient on the vector processing computers for which the current code is written. Detailed results are presented for two scramjet inlet configurations over a range of Mach numbers. The calculated results are compared with the available experimental and theoretical results.

Kumar, A.↗

Numerical study of boundary layer interaction with shocks: Method and code validation

A major problem in modeling of turbulent supersonic flows is the correct assessment of viscous-inviscid interaction problems. Of particular interest is the interaction of boundary layers with shocks. Present turbulence models give in most cases unsatisfactory results in the region of rapid distortion and in the separation region (if one is present) in particular with regard to mean flow profiles and turbulence quantities. The objective of the present work is the direct numerical simulation of shock boundary layer interaction. This report summarizes the first phase during which a numerical method suitable for this problem has been developed and a computer code has been written and tested.

Adams, Nikolaus A.↗

Computing interference flows past cylinder-fin-sting assembly in and near a cavity

The internal carriage of stores by the military aircraft is an option, particularly at high speeds, for a possible reduction in the aerodynamic drag and in the radar observability. Trade studies of this option, which include the safe separation of stores from their bays and aerodynamic complications of a cavity flow, require studying the aircraft and the stores together. With this motivation, the computer code, VISCC, is developed to solve the three-dimensional, Reynolds-averaged, unsteady, compressible, and complete Navier-Stokes equations. The store considered here for simulation is a missile configuration with four tail fins and an L-shaped offset sting. The missile is placed first directly above the opening of a rectangular box cavity, then submerged completely inside the cavity. The freestream Mach and the unit Reynolds numbers are 2.75 and 2.97 million per foot, respectively. The results of the time-accurate computations depict these complex flows and help understanding the interference effects. The time-averaged surface pressures compare favorably with the available wind tunnel data.

Baysal, O.↗

Unsteady three-dimensional thermal field prediction in turbine blades using nonlinear BEM

A time-and-space accurate and computationally efficient fully three dimensional unsteady temperature field analysis computer code has been developed for truly arbitrary configurations. It uses boundary element method (BEM) formulation based on an unsteady Green's function approach, multi-point Gaussian quadrature spatial integration on each panel, and a highly clustered time-step integration. The code accepts either temperatures or heat fluxes as boundary conditions that can vary in time on a point-by-point basis. Comparisons of the BEM numerical results and known analytical unsteady results for simple shapes demonstrate very high accuracy and reliability of the algorithm. An example of computed three dimensional temperature and heat flux fields in a realistically shaped internally cooled turbine blade is also discussed.

Martin, Thomas J.↗

Algorithm and code development for unsteady three-dimensional Navier-Stokes equations

In the last two decades, there have been extensive developments in computational aerodynamics, which constitutes a major part of the general area of computational fluid dynamics. Such developments are essential to advance the understanding of the physics of complex flows, to complement expensive wind-tunnel tests, and to reduce the overall design cost of an aircraft, particularly in the area of aeroelasticity. Aeroelasticity plays an important role in the design and development of aircraft, particularly modern aircraft, which tend to be more flexible. Several phenomena that can be dangerous and limit the performance of an aircraft occur because of the interaction of the flow with flexible components. For example, an aircraft with highly swept wings may experience vortex-induced aeroelastic oscillations. Also, undesirable aeroelastic phenomena due to the presence and movement of shock waves occur in the transonic range. Aeroelastically critical phenomena, such as a low transonic flutter speed, have been known to occur through limited wind-tunnel tests and flight tests. Aeroelastic tests require extensive cost and risk. An aeroelastic wind-tunnel experiment is an order of magnitude more expensive than a parallel experiment involving only aerodynamics. By complementing the wind-tunnel experiments with numerical simulations the overall cost of the development of aircraft can be considerably reduced. In order to accurately compute aeroelastic phenomenon it is necessary to solve the unsteady Euler/Navier-Stokes equations simultaneously with the structural equations of motion. These equations accurately describe the flow phenomena for aeroelastic applications. At Ames a code, ENSAERO, is being developed for computing the unsteady aerodynamics and aeroelasticity of aircraft and it solves the Euler/Navier-Stokes equations. The purpose of this contract is to continue the algorithm enhancements of ENSAERO and to apply the code to complicated geometries. During the last year, the geometric capability of the code was extended to simulate transonic flows, a wing with oscillating control surface. Single-grid and zonal approaches were tested. For the zonal approach, a new interpolation technique was introduced. The key development of the algorithm was an interface treatment between moving zones for a control surface using the virtual-zone concept. The work performed during the period, 1 Apr. 1992 through 31 Mar. 1993 is summarized. Additional details on the various aspects of the study are given in the Appendices.

Obayashi, Shigeru↗