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 523 records · Page 29

A Method for Aircraft Concept Selection Using Multicriteria Interactive Genetic Algorithms

The problem of aircraft concept selection has become increasingly difficult in recent years as a result of a change from performance as the primary evaluation criteria of aircraft concepts to the current situation in which environmental effects, economics, and aesthetics must also be evaluated and considered in the earliest stages of the decision-making process. This has prompted a shift from design using historical data regression techniques for metric prediction to the use of physics-based analysis tools that are capable of analyzing designs outside of the historical database. The use of optimization methods with these physics-based tools, however, has proven difficult because of the tendency of optimizers to exploit assumptions present in the models and drive the design towards a solution which, while promising to the computer, may be infeasible due to factors not considered by the computer codes. In addition to this difficulty, the number of discrete options available at this stage may be unmanageable due to the combinatorial nature of the concept selection problem, leading the analyst to arbitrarily choose a sub-optimum baseline vehicle. These concept decisions such as the type of control surface scheme to use, though extremely important, are frequently made without sufficient understanding of their impact on the important system metrics because of a lack of computational resources or analysis tools. This paper describes a hybrid subjective/quantitative optimization method and its application to the concept selection of a Small Supersonic Transport. The method uses Genetic Algorithms to operate on a population of designs and promote improvement by varying more than sixty parameters governing the vehicle geometry, mission, and requirements. In addition to using computer codes for evaluation of quantitative criteria such as gross weight, expert input is also considered to account for criteria such as aeroelasticity or manufacturability which may be impossible or too computationally expensive to consider explicitly in the analysis. Results indicate that concepts resulting from the use of this method represent designs which are promising to both the computer and the analyst, and that a mapping between concepts and requirements that would not otherwise be apparent is revealed.

Buonanno, Michael↗

Calculation of two-dimensional inlet flow fields in a supersonic free stream: Program documentation and test cases

The use of a computer code for the calculation of two dimensional inlet flow fields in a supersonic free stream and a nonorthogonal mesh-generation code are illustrated by specific examples. Input, output, and program operation and use are given and explained for the case of supercritical inlet operation at a subdesign Mach number (M Mach free stream = 2.09) for an isentropic-compression, drooped-cowl inlet. Source listings of the computer codes are also provided.

Biringen, S. H.↗

Automated Development of Accurate Algorithms and Efficient Codes for Computational Aeroacoustics

The simulation of sound generation and propagation in three space dimensions with realistic aircraft components is a very large time dependent computation with fine details. Simulations in open domains with embedded objects require accurate and robust algorithms for propagation, for artificial inflow and outflow boundaries, and for the definition of geometrically complex objects. The development, implementation, and validation of methods for solving these demanding problems is being done to support the NASA pillar goals for reducing aircraft noise levels. Our goal is to provide algorithms which are sufficiently accurate and efficient to produce usable results rapidly enough to allow design engineers to study the effects on sound levels of design changes in propulsion systems, and in the integration of propulsion systems with airframes. There is a lack of design tools for these purposes at this time. Our technical approach to this problem combines the development of new, algorithms with the use of Mathematica and Unix utilities to automate the algorithm development, code implementation, and validation. We use explicit methods to ensure effective implementation by domain decomposition for SPMD parallel computing. There are several orders of magnitude difference in the computational efficiencies of the algorithms which we have considered. We currently have new artificial inflow and outflow boundary conditions that are stable, accurate, and unobtrusive, with implementations that match the accuracy and efficiency of the propagation methods. The artificial numerical boundary treatments have been proven to have solutions which converge to the full open domain problems, so that the error from the boundary treatments can be driven as low as is required. The purpose of this paper is to briefly present a method for developing highly accurate algorithms for computational aeroacoustics, the use of computer automation in this process, and a brief survey of the algorithms that have resulted from this work. A review of computational aeroacoustics has recently been given by Lele.

Goodrich, John W.↗

Solution of 3-dimensional time-dependent viscous flows. Part 3: Application to turbulent and unsteady flows

A numerical scheme is developed for solving the time dependent, three dimensional compressible viscous flow equations to be used as an aid in the design of helicopter rotors. In order to further investigate the numerical procedure, the computer code developed to solve an approximate form of the three dimensional unsteady Navier-Stokes equations employing a linearized block implicit technique in conjunction with a QR operator scheme is tested. Results of calculations are presented for several two dimensional boundary layer flows including steady turbulent and unsteady laminar cases. A comparison of fourth order and second order solutions indicate that increased accuracy can be obtained without any significant increases in cost (run time). The results of the computations also indicate that the computer code can be applied to more complex flows such as those encountered on rotating airfoils. The geometry of a symmetric NACA four digit airfoil is considered and the appropriate geometrical properties are computed.

Weinberg, B. C.↗

Particle trajectory computation on a 3-dimensional engine inlet

A 3-dimensional particle trajectory computer code was developed to compute the distribution of water droplet impingement efficiency on a 3-dimensional engine inlet. The computed results provide the essential droplet impingement data required for the engine inlet anti-icing system design and analysis. The droplet trajectories are obtained by solving the trajectory equation using the fourth order Runge-Kutta and Adams predictor-corrector schemes. A compressible 3-D full potential flow code is employed to obtain a cylindrical grid definition of the flowfield on and about the engine inlet. The inlet surface is defined mathematically through a system of bi-cubic parametric patches in order to compute the droplet impingement points accurately. Analysis results of the 3-D trajectory code obtained for an axisymmetric droplet impingement problem are in good agreement with NACA experimental data. Experimental data are not yet available for the engine inlet impingement problem analyzed. Applicability of the method to solid particle impingement problems, such as engine sand ingestion, is also demonstrated.

Kim, J. J.↗

Application of an optimization method to high performance propeller designs

The application of an optimization method to determine the propeller blade twist distribution which maximizes propeller efficiency is presented. The optimization employs a previously developed method which has been improved to include the effects of blade drag, camber and thickness. Before the optimization portion of the computer code is used, comparisons of calculated propeller efficiencies and power coefficients are made with experimental data for one NACA propeller at Mach numbers in the range of 0.24 to 0.50 and another NACA propeller at a Mach number of 0.71 to validate the propeller aerodynamic analysis portion of the computer code. Then comparisons of calculated propeller efficiencies for the optimized and the original propellers show the benefits of the optimization method in improving propeller performance. This method can be applied to the aerodynamic design of propellers having straight, swept, or nonplanar propeller blades.

Li, K. C.↗

Application of an optimization method to high performance propeller designs

The application of an optimization method to determine the propeller blade twist distribution which maximizes propeller efficiency is presented. The optimization employs a previously developed method which has been improved to include the effects of blade drag, camber and thickness. Before the optimization portion of the computer code is used, comparisons of calculated propeller efficiencies and power coefficients are made with experimental data for one NACA propeller at Mach numbers in the range of 0.24 to 0.50 and another NACA propeller at a Mach number of 0.71 to validate the propeller aerodynamic analysis portion of the computer code. Then comparisons of calculated propeller efficiencies for the optimized and the original propellers show the benefits of the optimization method in improving propeller performance. This method can be applied to the aerodynamic design of propellers having straight, swept, or nonplanar propeller blades.

Li, K. C.↗

HOTCFGM-1D: A Coupled Higher-Order Theory for Cylindrical Structural Components with Through-Thickness Functionally Graded Microstructures

The objective of this three-year project was to develop and deliver to NASA Lewis one-dimensional and two-dimensional higher-order theories, and related computer codes, for the analysis, optimization and design of cylindrical functionally graded materials/structural components for use in advanced aircraft engines (e.g., combustor linings, rotor disks, heat shields, blisk blades). To satisfy this objective, a quasi one-dimensional version of the higher-order theory, HOTCFGM-1D, and four computer codes based on this theory, for the analysis, design and optimization of cylindrical structural components functionally graded in the radial direction were developed. The theory is applicable to thin multi-phased composite shell/cylinders subjected to macroscopically axisymmetric thermomechanical and inertial loading applied uniformly along the axial direction such that the overall deformation is characterized by a constant average axial strain. The reinforcement phases are uniformly distributed in the axial and circumferential directions, and arbitrarily distributed in the radial direction, thereby allowing functional grading of the internal reinforcement in this direction.

Pindera, Marek-Jerzy↗

Probabilistic Thermomechanical Fatigue of Polymer Matrix Composites

Traditional computational approaches for predicting the life and long-term behavior of materials rely on empirical data and are neither generic nor unique in nature. Also, those approaches are not easy to implement in a design procedure in an effective, integrated manner. The focus of ongoing research at the NASA Lewis Research Center has been to develop advanced integrated computational methods and related computer codes for a complete reliability-based assessment of composite structures. These methods - which account for uncertainties in all the constituent properties, fabrication process variables, and loads to predict probabilistic micromechanics, ply, laminate, and structural responses - have already been implemented in the Integrated Probabilistic Assessment of Composite Structures (IPACS) computer code. The main objective of this evaluation is to illustrate the effectiveness of the methodology to predict the long-term behavior of composites under combined mechanical and thermal cyclic loading conditions.

Chamis, Christos C.↗

Progressive Fracture of Composite Structures

A new approach is described for evaluating fracture in composite structures. This approach is independent of classical fracture mechanics parameters like fracture toughness. It relies on computational simulation and is programmed in a stand-alone integrated computer code. It is multiscale, multifunctional because it includes composite mechanics for the composite behavior and finite element analysis for predicting the structural response. It contains seven modules; layered composite mechanics (micro, macro, laminate), finite element, updating scheme, local fracture, global fracture, stress based failure modes, and fracture progression. The computer code is called CODSTRAN (Composite Durability Structural ANalysis). It is used in the present paper to evaluate the global fracture of four composite shell problems and one composite built-up structure. Results show that the composite shells and the built-up composite structure global fracture are enhanced when internal pressure is combined with shear loads.

Chamis, Christos C.↗

Multiscale Multifunctional Progressive Fracture of Composite Structures

A new approach is described for evaluating fracture in composite structures. This approach is independent of classical fracture mechanics parameters like fracture toughness. It relies on computational simulation and is programmed in a stand-alone integrated computer code. It is multiscale, multifunctional because it includes composite mechanics for the composite behavior and finite element analysis for predicting the structural response. It contains seven modules; layered composite mechanics (micro, macro, laminate), finite element, updating scheme, local fracture, global fracture, stress based failure modes, and fracture progression. The computer code is called CODSTRAN (Composite Durability Structural ANalysis). It is used in the present paper to evaluate the global fracture of four composite shell problems and one composite built-up structure. Results show that the composite shells. Global fracture is enhanced when internal pressure is combined with shear loads. The old reference denotes that nothing has been added to this comprehensive report since then.

Chamis, C. C.↗

Impact induced stresses, strains, and delaminations in composite plates

A method is presented for calculating the locations and sizes of delaminations which occur in a rectangular, fiber reinforced composite plate subjected to nonpenetrating (low velocity) impact of a solid object. The plate may be simply supported or clamped along its edges. In-plane loads or in-plane strains may be imposed on the plate during the impact. The method includes two steps. First, the stresses and strains in the plate are calculated by a three-dimensional, transient finite element method using 8-node brick elements with incompatible modes. Second, the locations, lengths, and widths of delaminations inside the plate are predicted by means of a proposed failure criterion, which is based on the concept of dimensional analysis. The finite element method and the failure criterion were implemented by a computer code which can be used to calculate the impactor position and velocity, the displacements of the plate, the stresses and strains inside the plate during the impact, and the locations and dimensions of the delaminations after the impact. Parametric studies were performed to illustrate the information which can be generated by the computer code.

Wu, Hsi-Yung T.↗

Analyzing rocket plume spectral data with neural networks

The Optical Plume Anomaly Detection (OPAD) system is under development to provide early-warning failure detection in support of ground-level testing of the Space Shuttle Main Engine (SSME). Failure detection is to be achieved through the acquisition of spectrally resolved plume emissions and subsequent identification of abnormal levels indicative of engine corrosion or component failure. Two computer codes (one linear and the other non-linear) are used by the OPAD system to iteratively determine specific element concentrations in the SSME plume, given emission intensity and wavelength information. Since this analysis is extremely labor intensive, a study was initiated to develop neural networks that would model the 'inverse' of these computer codes. Optimally connected feed-forward networks with imperceptible prediction error have been developed for each element modeled by the linear code, SPECTRA4. Radial basis function networks were developed for the non-linear code, SPECTRA5, and predict combustion temperature in addition to element concentrations.

Whitaker, Kevin W.↗

Prediction of sound radiation from different practical jet engine inlets

Computer codes which were based upon a special integral representation of the external solutions of the Helmholtz equation were upgraded so that they would yield accurate results for the acoustic radiation patterns in the field surrounding an axisymmetric body for nondimensional wave numbers, based on duct radius, of up to twenty. The accuracy of these computer programs was checked by the use of the point source method for the generation of exact solutions and then by comparison with the results of other experimental and theoretical studies. These computer codes were used in a parametric study of the dependence of the radiated sound field on input modal distribution, wave number, and inlet lip shape. The results of this study show that: (1) as the wave number (i.e., cut off ratio) is increased for a given input modal distribution that the acoustic radiation peak moves towards the inlet centerline and becomes more compact (i.e., narrows); (2) as the input mode number is increased for a given cut off ratio the acoustic radiation peak moves away from the inlet centerline and becomes more compact.

Meyer, W. L.↗

On 3-D inelastic analysis methods for hot section components (base program)

A 3-D Inelastic Analysis Method program is described. This program consists of a series of new computer codes embodying a progression of mathematical models (mechanics of materials, special finite element, boundary element) for streamlined analysis of: (1) combustor liners, (2) turbine blades, and (3) turbine vanes. These models address the effects of high temperatures and thermal/mechanical loadings on the local (stress/strain)and global (dynamics, buckling) structural behavior of the three selected components. Three computer codes, referred to as MOMM (Mechanics of Materials Model), MHOST (Marc-Hot Section Technology), and BEST (Boundary Element Stress Technology), have been developed and are briefly described in this report.

Wilson, R. B.↗

Advances in 3-D Inelastic Analysis Methods for hot section components

3-D Inelastic Analysis Methods are described. These methods consist of a series of new computer codes embodying a progression of mathematical models (mechanics of materials, specialty finite element, boundary element) for streamlined analysis of: (1) combustor liners, (2) turbine blades, and (3) turbine vanes. These models address the effects of high temperatures and thermal/mechanical loadings on the local (stress/strain) and global (dynamics, buckling) structural behavior of the three selected components. Three computer codes, referred to as MOMM (Mechanics of Materials Model), MHOST (MARC-Hot Section Technology), and BEST (Boundary Element Stress Technology), have been developed and are briefly described in this paper.

Chamis, Christos C.↗

1025:1 Area Ratio Nozzle Evaluated at High Combustion Chamber Pressure

A recently completed experimental test program obtained performance data on an optimally contoured nozzle with an exit-to-throat area ratio of 1025:1 and on a truncated version of this nozzle with an area ratio of 440:1. The nozzles were tested with gaseous hydrogen and liquid oxygen propellants at combustion chamber pressures of 12.4 to 16.5 mPa (1800 to 2400 psia). Testing was conducted in the altitude test capsule at the NASA Lewis Research Center's Rocket Engine Test Facility (RETF), and results were compared with analytical performance predictions. This testing builds on previous work with this nozzle at Lewis, where testing was completed at nominal chamber pressure of 350 psia. High-area-ratio nozzles have long been sought as a means to increase the performance of spacebased rocket engines. However, as the area ratio increases, the physical size and weight of the nozzle also increase. As a result, engine and vehicle designers must make tradeoffs between nozzle size and performance enhancement. Until this test program, very little experimental data existed on the performance of the high-area-ratio nozzles used in rocket engine designs. The computer codes being used by rocket engine designers rely on data extrapolated from tests of low-area-ratio nozzles, and these extrapolations do not always provide the accuracy needed for a reliable design assessment. Therefore, we conducted this high-area-ratio nozzle testing program to provide performance data for use in rocket engine design and analysis computer codes. The nozzle had a nominal 2.54-cm- (1-in.-) diameter throat, an exit diameter of 81.3-cm (32.0-in.) at an exit-to-throat area ratio of 1025, and a length of 128.6 cm (50.6 in.). Testing was conducted in an altitude test capsule to simulate the static pressure at altitude by vacuum pumping. Data such as propellant mass flow, oxidizer-to-fuel mixture, and thrust were measured. These measurements were then used to calculate performance factors such as the thrust coefficient, the characteristic exhaust velocity efficiency, and the vacuum specific impulse. In addition, the nozzle temperature was measured to calculate the amount of heat transferred from the combustion gases to the nozzle.

Source record↗

Development of Computational Aeroacoustics Code for Jet Noise and Flow Prediction

Accurate prediction of jet fan and exhaust plume flow and noise generation and propagation is very important in developing advanced aircraft engines that will pass current and future noise regulations. In jet fan flows as well as exhaust plumes, two major sources of noise are present: large-scale, coherent instabilities and small-scale turbulent eddies. In previous work for the NASA Glenn Research Center, three strategies have been explored in an effort to computationally predict the noise radiation from supersonic jet exhaust plumes. In order from the least expensive computationally to the most expensive computationally, these are: 1) Linearized Euler equations (LEE). 2) Very Large Eddy Simulations (VLES). 3) Large Eddy Simulations (LES). The first method solves the linearized Euler equations (LEE). These equations are obtained by linearizing about a given mean flow and the neglecting viscous effects. In this way, the noise from large-scale instabilities can be found for a given mean flow. The linearized Euler equations are computationally inexpensive, and have produced good noise results for supersonic jets where the large-scale instability noise dominates, as well as for the tone noise from a jet engine blade row. However, these linear equations do not predict the absolute magnitude of the noise; instead, only the relative magnitude is predicted. Also, the predicted disturbances do not modify the mean flow, removing a physical mechanism by which the amplitude of the disturbance may be controlled. Recent research for isolated airfoils' indicates that this may not affect the solution greatly at low frequencies. The second method addresses some of the concerns raised by the LEE method. In this approach, called Very Large Eddy Simulation (VLES), the unsteady Reynolds averaged Navier-Stokes equations are solved directly using a high-accuracy computational aeroacoustics numerical scheme. With the addition of a two-equation turbulence model and the use of a relatively coarse grid, the numerical solution is effectively filtered into a directly calculated mean flow with the small-scale turbulence being modeled, and an unsteady large-scale component that is also being directly calculated. In this way, the unsteady disturbances are calculated in a nonlinear way, with a direct effect on the mean flow. This method is not as fast as the LEE approach, but does have many advantages to recommend it; however, like the LEE approach, only the effect of the largest unsteady structures will be captured. An initial calculation was performed on a supersonic jet exhaust plume, with promising results, but the calculation was hampered by the explicit time marching scheme that was employed. This explicit scheme required a very small time step to resolve the nozzle boundary layer, which caused a long run time. Current work is focused on testing a lower-order implicit time marching method to combat this problem.

Keith, Theo G., Jr.↗