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.
Our immediate objective is to understand the limitations of guided-wave and grating coupler devices in their application to optical data storage. Our long-range goal is to develop and validate design codes for integrated optic devices. The principal research activity was in the development of numerical models for the design of a blue wavelength integrated optical source for data storage applications.
LinAir is a vortex lattice aerodynamic prediction method similar to Weissinger's extended lifting-line theory, except that the circulation around a wing is represented by discrete horseshoe vortices, not a continuous distribution of vorticity. The program calculates subsonic longitudinal and lateral/directional aerodynamic forces and moments for arbitrary aircraft geometries. It was originally written by Dr. Ilan Kroo of Stanford University, and subsequently modified by the author to simplify modeling of complex configurations. The Polhamus leading-edge suction analogy was added by the author to extend the range of applicability of LinAir to low aspect ratio (i.e., fighter-type) configurations. A brief discussion of the theory of LinAir is presented, and details on how to run the program are given along with some comparisons with experimental data to validate the code. Example input and output files are given in the appendices to aid in understanding the program and its use. This version of LinAir runs in the VAX/VMS, Cray UNICOS, and Silicon Graphics Iris workstation environments at the time of this writing.
The current research is aimed at developing and extending numerical methods to accurately predict the high Reynolds number flow about the NASA F-18 HARV at large angles of attack. The resulting codes are validated by comparison of the numerical results with in-flight aerodynamic measurements and flow visualization obtained on the HARV. Further, computations have been used to provide an analysis and numerical optimization of a pneumatic slot blowing concept, and a mechanical strake concept, for use as potential forebody flow control devices in improving high-alpha maneuverability.
This test was initiated to provide validation data on low aspect ratio wings at transonic speeds. The test was conducted so that the data obtained would be useful in the validation of codes, and all boundary condition data required would be measured as part of the test. During the conduct of the test, the measured quantities were checked for repeatability, and when the data would not repeat, the cause was tracked down and either eliminated or included in the measurement uncertainty. The accuracy of the data was in the end limited by wall imperfections of the wind tunnel in which the test was run.
Film/dump cooling a rocket nozzle with fuel rich gas, as in the National Launch System (NLS) Space Transportation Main Engine (STME), adds potential complexities for integrating the engine with the vehicle. The chief concern is that once the film coolant is exhausted from the nozzle, conditions may exist during flight for the fuel-rich film gases to be recirculated to the vehicle base region. The result could be significantly higher base temperatures than would be expected from a regeneratively cooled nozzle. CFD analyses were conduced to augment classical scaling techniques for vehicle base environments. The FDNS code with finite rate chemistry was used to simulate a single, axisymmetric STME plume and the NLS base area. Parallel calculations were made of the Saturn V S-1 C/F1 plume base area flows. The objective was to characterize the plume/freestream shear layer for both vehicles as inputs for scaling the S-C/F1 flight data to NLS/STME conditions. The code was validated on high speed flows with relevant physics. This paper contains the calculations for the NLS/STME plume for the baseline nozzle and a modified nozzle. The modified nozzle was intended to reduce the fuel available for recirculation to the vehicle base region. Plumes for both nozzles were calculated at 10kFT and 50kFT.
The Regenerator, Cooler and Heater for the NASA Space Power Research Engine (SPRE) have been analyzed in detail for laminar, incompressible and oscillatory flow conditions. Each component has been analyzed independently and in detail with the regenerator being modeled as two-parallel-plates channel with a solid wall. The ends of the channel are exposed to two reservoir maintained at different temperature thus providing an axial temperature gradient along the channel. The cooler and heater components have been modeled as circular pipes with isothermal walls. Two different types of thermal boundary conditions have been investigated for the cooler and heater, namely, symmetric and asymmetric temperature inflow. In symmetric temperature inflow the flow enters the channel with the same temperature in throughout the velocity cycle whereas, in asymmetric temperature inflow the flow enters with a different temperature in each half cycle. The study was conducted over a wide range of Maximum Reynolds number (RE(max) varying from 75 to 60000, Valensi number (Va) from 2.5 to 800, and relative amplitude of fluid displacement (A(sub r) from 0.357 to 1.34. A two dimensional Finite volume method based on the SIMPLE algorithm was used to solve the governing partial differential equations. Post processing programs were developed to effectively describe the heat transfer mechanism under oscillatory flows. The computer code was validated by comparing with existing analytical solutions for oscillating flows. The thermal field have been studied with the help of temperature contour and three dimensional plots. The instantaneous friction factor, wall heat flux and heat transfer coefficient have been examined. It has been concluded that in general, the frictional factor and heat transfer coefficient are higher under oscillatory flow conditions when the Valensi number is high. Also, the thermal efficiency decreases for lower A(r) values. Further, the usual steady state definition for the heat transfer coefficient does not seem to be valid.
Spectral calibration of the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) as data are acquired in flight is essential to quantitative analysis of the measured upwelling spectral radiance. In each spectrum measured by AVIRIS in flight, there are numerous atmospheric gas absorption bands that drive this requirement for accurate spectral calibration. If the surface and atmospheric properties are measured independently, these atmospheric absorption bands may be used to deduce the in-flight spectral calibration of an imaging spectrometer. Both the surface and atmospheric characteristics were measured for a calibration target during an in-flight calibration experiment held at Lunar Lake, Nevada on April 5, 1994. This paper uses upwelling spectral radiance predicted for the calibration target with the MODTRAN radiative transfer code to validate the spectral calibration of AVIRIS in flight.
This paper presents an extension of the nonlinear least squares fitting technique of Vinas and Scudder (1986) (VS), which finds the physical and geometrical properties of nondissipational magnetohydrodynamic (MHD) shocks. The new method incorporates plasma temperature observations in the form of normal momentum flux and energy density flux conservation as well as plasma density, velocity, and magnetic field data. The new technique is capable of using known standard deviations in the individual measurement points to properly weight the fitting procedure. The new fitting code is validated through the analysis of synthetic shocks with known physical and geometrical properties. Finally, it is compared to the original VS method and the preaveraged velocity coplanarity technique.
A study has been conducted to develop and to analyze a FORTRAN computer code for performing agility analysis on fighter aircraft configurations. This program is one of the modules of the NASA Ames ACSYNT (AirCraft SYNThesis) design code. The background of the agility research in the aircraft industry and a survey of a few agility metrics are discussed. The methodology, techniques, and models developed for the code are presented. The validity of the existing code was evaluated by comparing with existing flight test data. A FORTRAN program was developed for a specific metric, PM (Pointing Margin), as part of the agility module. Example trade studies using the agility module along with ACSYNT were conducted using a McDonnell Douglas F/A-18 Hornet aircraft model. Tile sensitivity of thrust loading, wing loading, and thrust vectoring on agility criteria were investigated. The module can compare the agility potential between different configurations and has capability to optimize agility performance in the preliminary design process. This research provides a new and useful design tool for analyzing fighter performance during air combat engagements in the preliminary design.
A gradient-based shape optimization methodology, that is intended for practical three-dimensional aerodynamic applications, has been developed. It is based on the quasi-analytical sensitivities. The flow analysis is rendered by a fully implicit, finite volume formulation of the Euler equations.The aerodynamic sensitivity equation is solved using the alternating-direction-implicit (ADI) algorithm for memory efficiency. A flexible wing geometry model, that is based on surface parameterization and platform schedules, is utilized. The present methodology and its components have been tested via several comparisons. Initially, the flow analysis for for a wing is compared with those obtained using an unfactored, preconditioned conjugate gradient approach (PCG), and an extensively validated CFD code. Then, the sensitivities computed with the present method have been compared with those obtained using the finite-difference and the PCG approaches. Effects of grid refinement and convergence tolerance on the analysis and shape optimization have been explored. Finally the new procedure has been demonstrated in the design of a cranked arrow wing at Mach 2.4. Despite the expected increase in the computational time, the results indicate that shape optimization, which require large numbers of grid points can be resolved with a gradient-based approach.
There is a need for more faithful simulation of space debris impacts on various space vehicles. Space debris impact velocities can range up to 14 km/sec and conventional two-stage light gas guns with moderately heavy saboted projectiles are limited to launch velocities of 7-8 km/sec. Any increases obtained in the launch velocities will result in more faithful simulations of debris impacts. It would also be valuable to reduce the maximum gun and projectile base pressures and the gun barrel erosion rate. In this paper, the results of a computational fluid dynamics (CFD) study designed to optimize the performance of the NASA Ames 0.5' gun by systematically varying seven gun operating parameters are reported. Particularly beneficial effects were predicted to occur if (1) the piston mass was decreased together with the powder mass and the hydrogen fill pressure and (2) the pump tube length was decreased. The optimum set of changes in gun operating conditions were predicted to produce an increase in muzzle velocity of 0.7-1.0 km/sec, simultaneously with a substantial decrease in gun erosion. Preliminary experimental data have validated the code predictions. Velocities of up to 8.2 km/sec with a 0.475 cm diameter saboted aluminum sphere have been obtained, along with large reductions in gun erosion rates.
Probabilistic method is not a universally accepted approach for the design and analysis of aerospace structures. The validity of this approach must be demonstrated to encourage its acceptance as it viable design and analysis tool to estimate structural reliability. The objective of this Study is to develop a well characterized finite population of similar aerospace structures that can be used to (1) validate probabilistic codes, (2) demonstrate the basic principles behind probabilistic methods, (3) formulate general guidelines for characterization of material drivers (such as elastic modulus) when limited data is available, and (4) investigate how the drivers affect the results of sensitivity analysis at the component/failure mode level.
The objective of the project is to determine the relative contributions of gravity-driven fluid flows to the compositional distribution, incorporation of impurities and defects, and deviation from stoichiometry observed in the crystals grown by vapor transport as results of buoyance-driven convection and growth interface fluctuations caused by irregular fluid-flows. ZnSe and related ternary compounds, such as ZnSeS and ZnSeTe, were grown by vapor transport technique with real time in-situ non-invasive monitoring techniques. The grown crystals were characterized extensively to correlate the grown crystal properties with the growth conditions. The following are the research progress in the past two years. In-situ monitoring of partial pressure by optical absorption technique and visual observation of the growing crystal were performed during vapor growth of ZnSe. Low-temperature photoluminescence (PL) spectra and glow discharge mass spectroscopy (GDMS) were measured on ZnSe starting materials provided by various vendors and on bulk crystals grown from these starting materials by physical vapor transport (PVT) to study the effects of purification and contamination during crystal growth process. Optical characterization was performed on wafers sliced from the grown crystals of ZnSe, ZnTe and ZnSe(1-x),Te(x), (0<x<0.4). Energy band gaps at room temperature were determined from optical transmission measurements and a best fit curve to the band gap vs. composition, x, data gives a bowing parameter of 1.45. Low-temperature photoluminescence (PL) spectra of ZnSe and ZnTe were dominated by near band edge emissions and no deep donor-acceptor pairs were observed. The PL spectrum exhibited a broad emission for the ZnSe(1-x),Te(x), samples, 0.09<x<0.39. The single broad PL emission spectra and the spectra measured as a function of temperature were interpreted as being associated with the exciton bound to Te clusters because of the high Te content in these samples. To validate numerical codes, in-situ monitoring during the PVT of HgI2 was designed. Optical absorption spectra of the vapor phase over HgI2 were measured for wavelengths between 200 and 600nm at sample temperatures between 349 and 610K. The Beer's Law constants for 15 wavelengths between 200 and 440 nm were determined. From these constants the vapor pressure of HgI2 was established as a function of temperature for the liquid and the solid Beta-phases To characterize the growth conditions during the PVT growth of In-doped ZnSe the optical absorbance of the vapor phase over the In-Se system were measured and were used to obtain the partial pressures of Se2(g) and In2Se(g).
A self-consistent ring current (RC) model has been developed that couples electron and ion magnetospheric dynamics with the calculation of the electric field. Two new features were taken into account in order to close the self-consistent magnetosphere- ionosphere coupling loop. First, in addition to the RC ions, we solve an electron kinetic equation in our model. Second, using the relation of Galand and Richmond [2001], we calculate the height integrated ionospheric conductances as a function of the precipitated high energy magnetospheric electrons and ions that are produced by our code. To validate the results of our model we simulate the magnetic storm of May 2, 1986, a storm that has has been comprehensively studied, and compare our results with different theoretical approaches. The self-consistent inclusion of the hot electrons and, their effect on the conductance results in deeper penetration of the magnetospheric electric field. In addition, a slight westward rotation of the potential pattern (compared to previous self-consistent results) is evident in the inner magnetosphere. This changes the hot plasma distribution, especially by allowing increased access of plasma sheet ions and electrons to low L shells.
Macrosegregation in metal casting can be caused by thermal and solutal melt convection, and the transport of unattached solid crystals. These free grains can be a result of, for example, nucleation in the bulk liquid or dendrite fragmentation. In an effort to develop a comprehensive numerical model for the casting of alloys, an experimental study has been conducted to generate benchmark data with which such a solidification model could be tested. The specific goal of the experiments was to examine equiaxed solidification in situations where sinking of grains is (and is not) expected. The objectives were: 1) experimentally study the effects of solid transport and thermosolutal convection on macrosegregation and grain size distribution patterns; and 2) provide a complete set of controlled thermal boundary conditions, temperature data, segregation data, and grain size data, to validate numerical codes. The alloys used were Al-1 wt. pct. Cu, and Al-10 wt. pct. Cu with various amounts of the grain refiner TiB2 added. Cylindrical samples were either cooled from the top, or the bottom. Several trends in the data stand out. In attempting to model these experiments, concentrating on experiments that show clear trends or differences is recommended.
A thermomechanical model for shape memory alloy (SMA) actuators and SMA hybrid composite (SMAHC) structures has been recently implemented in the commercial finite element codes MSC.Nastran and ABAQUS. The model may be easily implemented in any code that has the capability for analysis of laminated composite structures with temperature dependent material properties. The model is also relatively easy to use and requires input of only fundamental engineering properties. A brief description of the model is presented, followed by discussion of implementation and usage in the commercial codes. Results are presented from static and dynamic analysis of SMAHC beams of two types; a beam clamped at each end and a cantilevered beam. Nonlinear static (post-buckling) and random response analyses are demonstrated for the first specimen. Static deflection (shape) control is demonstrated for the cantilevered beam. Approaches for modeling SMAHC material systems with embedded SMA in ribbon and small round wire product forms are demonstrated and compared. The results from the commercial codes are compared to those from a research code as validation of the commercial implementations; excellent correlation is achieved in all cases.
The Higher Order Theory - Structural/Micro Analysis Code (HOT-SMAC) software package is applied to analyze the linearly elastic and elasto-plastic response of adhesively bonded tee joints. Joints of this type are finding an increasing number of applications with the increased use of composite materials within advanced aerospace vehicles, and improved tools for the design and analysis of these joints are needed. The linearly elastic results of the code are validated vs. finite element analysis results from the literature under different loading and boundary conditions, and new results are generated to investigate the inelastic behavior of the tee joint. The comparison with the finite element results indicates that HOT-SMAC is an efficient and accurate alternative to the finite element method and has a great deal of potential as an analysis tool for a wide range of bonded joints.
One of the primary uses of the in-flight icing research performed aboard NASA Glenn s DHC-6 Twin Otter is for Icing Research Tunnel (IRT) and icing prediction code (Lewice) validation. Using the in-flight data to establish the IRT and Lewice as accurate simulators of actual icing conditions is crucial for supporting the research done in the Icing Branch. During test flights during the 2003 and 2004 flight season, a Natural Ice Shape Database was collected. For flights where conditions were appropriate, the aircraft is flown in an icing cloud with all ice protection systems deactivated. The duration of this period is usually determined by the pilot s ability to safely control the aircraft. When safe flight is no longer possible, the aircraft is maneuvered into clear air above the cloud layer. At this point several photographs are taken of the ice shape that was accreted on the wing test section during this icing encounter using a stereo photograph system (Figure 1). The stereo photograph system utilizes two cameras located at different locations on the fuselage that are both pointed at the same location on the wing. When both cameras take photographs of the same location at the same time, the negatives can be combined digitally to generate a two dimensional plot describing the cross-section of the ice shape. After these photographs are taken, the wing de-icing boots are activated and the ice shape is removed.