SEARCH · Engineering Papers
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.
Program developed for computing aerodynamic design of multistage axialflow compressor and associated blading geometry input for internal flow analysis. Aerodynamic solution gives velocity diagrams on selected streamlines of revolution at blade row edges. Program written in FORTRAN IV.
The development of a three-dimensional thermal analysis model of the Lidar Atmospheric Sensing Experiment (LASE) is the subject of this paper. The use of an interactive computer graphics and finite element generation program to define the geometry information for the thermal model is discussed. The methods used in calculating the heat transfer parameters are explained. The results of the thermal analysis are given, and these results are compared with actual flight data.
The necessary information for an aerodynamic investigation requiring load cell force measurements at the National Full-Scale Aerodynamics Complex (NFAC) is provided. Included are details of the Ames 40x80 three component load cells; typical model/load cell installation geometries; transducer signal conditioning; a description of the Ames Standard Computations Wind Tunnel Data Reduction Program for Load Cells Forces and Moments (SCELLS), and the inputs required for SCELLS. The Outdoor Aerodynamic Facilities Complex (OARF), a facility within the NFAC where three axes load cells serve as the primary balance system, is used as an example for many of the techniques, but the information applies equally well to other static and wind tunnel facilities that make use of load cell balances.
This paper describes progress in validating incompressible Navier-Stokes codes using image-processed flow-field data. Solutions of the Navier-Stokes equations using a predictor/corrector scheme are compared directly with experimental data. A uniform asymptotic expansion is used as an initial condition in order to correctly capture the initial growth of the vortex sheet. With the aid of an experimental study in a tow tank using the particle streak visualization method, measured and computed velocity vectors are compared directly. Comparison of the velocity data with computations shows little effect of initial conditions and varying grid geometry on the predictive capabilities.
Rotorcraft aerodynamic interaction investigations conducted experimentally and computationally by the U.S. Army at the NASA Langley subsonic tunnel are described. Consideration is given to fuselage velocity computation, an interactional aerodynamics method, a rotor/wake/fuselage method, predicted wake geometry, views of tip vortex trajectory, and an inflow velocity study. Calculations were found to predict the periodic geometry of a rotor wake. Details of its interactions with the fuselage surface are not well modeled with inviscid panel methods. For the cases investigated in the present study, the overall effects of the fuselage are relatively small. However, increasing the size of the fuselage or decreasing the fuselage-rotor space will amplify these effects.
A new numerical procedure for determining load proportioning in multiple-hole mechanical joints in composite plates is presented. The joints are loaded in double lap fashion in tension with pins through the holes. The commercial finite-element program ABAQUS is used to predict the load proportioning among holes using two independent plane stress finite-element models, one representing the composite inner lap and one representing the two steel outer laps, interacting through rigid circular surfaces. The circular surfaces effectively represent rigid pins. Load proportioning is predicted for a number of geometries. Excellent correlation with experimental data is obtained. Experimental and computed surface strains are also found to compare well. The assumption of a radial cosine distribution of contact stress between pin and hole is shown to be substantially in error for some holes.
A large-field laser holographic focusing schlieren system for high-speed flow visualization has been built and evaluated. This system is based on a recently improved large-field focusing schlieren technique and is combined with laser holography methods to record three-dimensional flows. A coordinated experimental and computational study of supersonic flows over wedge, cone, and sphere geometries was conducted to evaluate the capabilities and limitations of the system. The ability of the system to focus on planes normal to its optical axis is demonstrated. The sharpness of focus of the present system was found to be limited. Issues regarding obtaining quantitative measurements of the density gradient are discussed.
Solution-adaptive CFD codes based on unstructured methods for 3-D complex geometries in subsonic to supersonic regimes were investigated, and the computed solution data were analyzed in conjunction with experimental data obtained from wind tunnel measurements in order to assess and validate the predictability of the code. Specifically, the FELISA code was assessed and improved in cooperation with NASA Langley and Imperial College, Swansea, U.K.
This paper provides a description of procedures typically used to predict convective heating rates to hypersonic reentry vehicles using the two-layer method. These procedures were used to compute the pitch-plane heating distributions to the Apollo geometry for a wind tunnel test case and for three flight cases. Both simple engineering methods and coupled inviscid/boundary layer solutions were used to predict the heating rates. The sensitivity of the heating results in the choice of metrics, pressure distributions, boundary layer edge conditions, and wall catalycity used in the heating analysis were evaluated. Streamline metrics, pressure distributions, and boundary layer edge properties were defined from perfect gas (wind tunnel case) and chemical equilibrium and nonequilibrium (flight cases) inviscid flow-field solutions. The results of this study indicated that the use of CFD-derived metrics and pressures provided better predictions of heating when compared to wind tunnel test data. The study also showed that modeling entropy layer swallowing and ionization had little effect on the heating predictions.
A multigrid procedure that makes use of coarse grids generated by the agglomeration of control volumes is advocated as a practical approach for solving the three dimensional Euler equations on unstructured grids about complex configurations. It is shown that the agglomeration procedure can be tailored to achieve certain coarse grid properties such as the sizes of the coarse grids and aspect ratios of the coarse grid cells. The agglomeration is done as a preprocessing step and runs in linear time. The implications for multigrid of using arbitrary polyhedral coarse grids are discussed. The agglomeration multigrid technique compares very favorably with existing multigrid procedures both in terms of convergence rates and elapsed times. The main advantage of the present approach is the ease with which coarse grids of any desired degree of coarseness may be generated in three dimensions, without being constrained by considerations of geometry. Inviscid flows over a variety of complex configurations are computed using the agglomeration multigrid strategy.
This report summarizes research done over the past two years as part of NASA Grant NCC 2-729. This research has been aimed at validating numerical methods for computing the flow about the complete F-18 HARV at alpha = 30 deg and alpha = 45 deg. At 30 deg angle of attack, the flow about the F-18 is dominated by the formation, and subsequent breakdown, of strong vortices over the wing leading-edge extensions (LEX). As the angle of attack is increased to alpha = 45 deg, the fuselage forebody of the F-18 contains significant laminar and transitional regions which are not present at alpha = 30 deg. Further, the flow over the LEX at alpha = 45 deg is dominated by an unsteady shedding in time, rather than strong coherent vortices. This complex physics, combined with the complex geometry of a full aircraft configuration, provides a challenge for current computational fluid dynamics (CFD) techniques. The following sections present the numerical method and grid generation scheme that was used, a review of prior research done to numerically model the F-18 HARV, and a discussion of the current research. The current research is broken into two main topics: the effect of engine-inlet mass-flow rate on the F-18 vortex breakdown position, and the results using a refined F-18 computational model to compute the flow at alpha = 30 deg and alpha = 45 deg.
The reliability with bounded distribution parameters (mean, standard deviation) was maximized and the reliability-based cost was minimized for adaptive intra-ply hybrid fiber composites by using a probabilistic method. The probabilistic method accounts for all naturally occurring uncertainties including those in constituent material properties, fabrication variables, structure geometry, and control-related parameters. Probabilistic sensitivity factors were computed and used in the optimization procedures. For actuated change in the angle of attack of an airfoil-like composite shell structure with an adaptive torque plate, the reliability was maximized to 0.9999 probability, with constraints on the mean and standard deviation of the actuation material volume ratio (percentage of actuation composite material in a ply) and the actuation strain coefficient. The reliability-based cost was minimized for an airfoil-like composite shell structure with an adaptive skin and a mean actuation material volume ratio as the design parameter. At a O.9-mean actuation material volume ratio, the minimum cost was obtained.
Boundary-layer bleed in supersonic inlets is typically used to avoid separation from adverse shock-wave/boundary-layer interactions and subsequent total pressure losses in the subsonic diffuser and to improve normal shock stability. Methodologies used to determine bleed requirements are reviewed. Empirical sonic flow coefficients are currently used to determine the bleed hole pattern. These coefficients depend on local Mach number, pressure ratio, hole geometry, etc. A new analytical bleed method is presented to compute sonic flow coefficients for holes and narrow slots and predictions are compared with published data to illustrate the accuracy of the model. The model can be used by inlet designers and as a bleed boundary condition for computational fluid dynamic studies.
The goal of our research is to produce a flexible, general grid generator for automated use by other programs, such as numerical optimizers. The current trend in the gridding field is toward interactive gridding. Interactive gridding more readily taps into the spatial reasoning abilities of the human user through the use of a graphical interface with a mouse. However, a sometimes fruitful approach to generating new designs is to apply an optimizer with shape modification operators to improve an initial design. In order for this approach to be useful, the optimizer must be able to automatically grid and evaluate the candidate designs. This paper describes and intelligent gridder that is capable of analyzing the topology of the spatial domain and predicting approximate physical behaviors based on the geometry of the spatial domain to automatically generate grids for computational fluid dynamics simulators. Typically gridding programs are given a partitioning of the spatial domain to assist the gridder. Our gridder is capable of performing this partitioning. This enables the gridder to automatically grid spatial domains of wide range of configurations.
A probabilistic design method is described and demonstrated using a smart composite wing. Probabilistic structural design incorporates naturally occurring uncertainties including those in constituent (fiber/matrix) material properties, fabrication variables, structure geometry and control-related parameters. Probabilistic sensitivity factors are computed to identify those parameters that have a great influence on a specific structural reliability. Two performance criteria are used to demonstrate this design methodology. The first criterion requires that the actuated angle at the wing tip be bounded by upper and lower limits at a specified reliability. The second criterion requires that the probability of ply damage due to random impact load be smaller than an assigned value. When the relationship between reliability improvement and the sensitivity factors is assessed, the results show that a reduction in the scatter of the random variable with the largest sensitivity factor (absolute value) provides the lowest failure probability. An increase in the mean of the random variable with a negative sensitivity factor will reduce the failure probability. Therefore, the design can be improved by controlling or selecting distribution parameters associated with random variables. This can be implemented during the manufacturing process to obtain maximum benefit with minimum alterations.
The production rates of cosmogenic nuclides depend on the primary cosmic-ray particles, the irradiated-body's bulk composition, size, and shape, and the sample's composition and shielding depth. Although much work has been done on some of these dependencies, more detailed studies still need to be done on others. This work describes the influence of irradiation geometry on nuclide production in lunar rocks. In most cases, computer simulations of cosmogenic nuclide production were restricted to spherical objects irradiated with a 4 pi isotropic flux (meteoroids) or in lunar core samples irradiated by a 2 pi flux incident on semi-infinite layers or cylinders of huge sizes. Many lunar samples are rocks found on top of the lunar surface. For these rocks, neither of the above-mentioned models correspond to the real conditions. We present results of our simulations of cosmogenic nuclide production in models simulating the irradiation of rocks sitting on top of the lunar surface. The Galactic Cosmic Rays (GCR) production profiles in lunar rocks were calculated using the Los Alamos 3-D Monte Carlo LAHET Code System (LCS). The irradiated object was modeled as the union of a sphere with the radius of the Moon and a small hemisphere with radii varying from 10 to 100 g/sq cm simulating the lunar rock. These calculations for the production of cosmogenic nuclides in lunar rocks by GCR particle show that there are important differences between the results obtained by commonly used geometric irradiation models and the lunar-rock models presented. The steeper GCR production profiles for a rock could help to explain the poor agreement for Be-10 in rock 68815, where slab models give GCR profiles flatter than the observed profiles.
Large shifts of liquid arising from small changes in certain container shapes in zero gravity can be used as a basis for accurately determining contact angle. Canonical geometries for this purpose, recently developed mathematically, are investigated here computationally. It is found that the desired nearly-discontinuous behavior can be obtained and that the shifts of liquid have sufficient volume to be readily observed.
An aerodynamic shape optimization procedure based on discrete sensitivity analysis is extended to treat three-dimensional geometries. The function of sensitivity analysis is to directly couple computational fluid dynamics (CFD) with numerical optimization techniques, which facilitates the construction of efficient direct-design methods. The development of a practical three-dimensional design procedures entails many challenges, such as: (1) the demand for significant efficiency improvements over current design methods; (2) a general and flexible three-dimensional surface representation; and (3) the efficient solution of very large systems of linear algebraic equations. It is demonstrated that each of these challenges is overcome by: (1) employing fully implicit (Newton) methods for the CFD analyses; (2) adopting a Bezier-Bernstein polynomial parameterization of two- and three-dimensional surfaces; and (3) using preconditioned conjugate gradient-like linear system solvers. Whereas each of these extensions independently yields an improvement in computational efficiency, the combined effect of implementing all the extensions simultaneously results in a significant factor of 50 decrease in computational time and a factor of eight reduction in memory over the most efficient design strategies in current use. The new aerodynamic shape optimization procedure is demonstrated in the design of both two- and three-dimensional inviscid aerodynamic problems including a two-dimensional supersonic internal/external nozzle, two-dimensional transonic airfoils (resulting in supercritical shapes), three-dimensional transport wings, and three-dimensional supersonic delta wings. Each design application results in realistic and useful optimized shapes.