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.
In this summary we present extremal average mutual information values carried by the extrinsic loglikelihood under the constraint of a given mean and variance while accounting for the consistency feature of loglikelihoods.
A simple analytical method has been introduced for aerohydrodynamic load analysis of an airborne configuration during water ditching. The method employs an aerodynamic panel code, based on linear potential flow theory, to simulate the flow of air and water around an aircraft configuration. The free surface separating the air and water region is represented by doublet" singularities. Although all the theoretical load distributions are computed for air, provisions are made to correct the pressure coefficients obtained on the configuration-wetted surfaces to account for the water density. As an analytical tool, a vortex aerodynamic code is chosen to carry out the present investigation. After assessing· the validity of the method, an application is presented for the water ditching of the Space Shuttle configuration at a 12-deg attitude.
An account is given of the numerical methods employed in a code for the simulation of supersonic combustion, which is then applied to the simulation of attached detonations and flames associated with the oblique-detonation wave supersonic combustor concept. The addition of heat by a detonation wave results in a shorter combustor than can be obtained in more conventional scramjet designs. Pure oblique detonations have been produced in a stoichiometric, uniformly mixed hydrogen/air stream; the wave rotates upstream with energy release, according to simple analytical arguments. Flow visualization maps for Mach number and temperature are presented.
The application of the computer code IPACS (Integrated Probabilistic Assessment of Composite Structures) to aircraft wing type structures is described. The code performs a complete probabilistic analysis for composites taking into account the uncertainties in geometry, boundary conditions, material properties, laminate lay-ups, and loads. Results of the analysis are presented in terms of cumulative distribution functions (CDF) and probability density function (PDF) of the fatigue life of a wing type composite structure under different hygrothermal environments subjected to the random pressure. The sensitivity of the fatigue life to a number of critical structural/material variables is also computed from the analysis.
Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape,spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars. Ongoing field programs such as Department of Energy's Atmospheric Radiation Measurement (DoE/ARM) are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS-AM1 instruments which will provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", G,which depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple IPA estimates of cloud radiative properties. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of adiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international "Intercomparison of 3D Radiation Codes" or I3RC, program is beginning to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.qsfc.nasa.crov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large- scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate.
Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape, spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars and lidars. Ongoing field programs such as DoE/ARM are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS (Earth Observing System) instruments which are beginning to provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", (G). G depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple estimates of cloud radiative properties that are based on the "Independent Pixel Approximation" or IPA. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international "Intercomparison of 3D Radiation Codes" or I3RC, program is underway to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.gsfc.nasa.gov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large-scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate the effects of cloud inhomogeneity and 3D radiation.
The satellite selection method as utilized by the spaceborne Global Positioning System receiver provides navigational solutions and is designed for use in low Earth orbit. The satellite selection method is a robust algorithm that can be used a GPS receiver to select appropriate GPS satellites for use in calculating point solutions or attitude solutions. The method is takes into account the difficulty of finding a particular GPS satellite phase code, especially when the search range in greatly increased due to Doppler shifts introduced into the carrier frequency. The method starts with an update of the antenna pointing and spacecraft vectors to determine the antenna backplane direction. Next, the GPS satellites that will potentially be in view of the antenna are ranked on a list, whereby the list is generated based on the estimated attitude and position of each GPS satellite. Satellites blocked by the Earth are not entered on this list. A second list is created, whereby the GPS satellites are ranked according to their desirability for use in attitude determination. GPS satellites are ranked according to their orthogonality to the antenna backplane, and according to geometric dilution of precision considerations. After the lists are created, the channels of the spaceborne GPS receiver are assigned to various GPS satellites for acquisition and lock. Preliminary Doppler frequencies for searching are assigned to the various channels.
Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape, spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars and lidars. Ongoing field programs such as DoE/ARM (atmospheric radiation measurement) are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS (Earth observing system) instruments which are beginning to provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of 'radiative smoothing' of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. 'Offbeam' Lidar returns are now being used to estimate the cloud 'radiative Green's function' (G). G depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple estimates of cloud radiative properties that are based on the 'Independent Pixel Approximation' or IPA. This and other measurements of 3D (three dimensional) transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international 'Intercomparison of 3D Radiation Codes' (I3RC) program is underway to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.gsfc.nasa.gov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large-scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate the effects of cloud inhomogeneity and 3D radiation.
In the development of advanced aircraft-engine fan cases and containment systems, composite materials are beginning to be used due to their low weight and high strength. The design of these structures must include the capability of withstanding impact loads from a released fan blade. Relatively complex triaxially braided fiber architectures have been found to yield the best performance for the fan cases. To properly work with and design these structures, robust analytical tools are required that can be used in the design process. A new analytical approach models triaxially braided carbon fiber composite materials within the environment of a transient dynamic finite-element code, specifically the commercially available transient dynamic finite-element code LS-DYNA. The geometry of the braided composites is approximated by a series of parallel laminated composites. The composite is modeled by using shell finite elements. The material property data are computed by examining test data from static tests on braided composites, where optical strain measurement techniques are used to examine the local strain variations within the material. These local strain data from the braided composite tests are used along with a judicious application of composite micromechanics- based methods to compute the stiffness properties of an equivalent unidirectional laminated composite required for the shell elements. The local strain data from the braided composite tests are also applied to back out strength and failure properties of the equivalent unidirectional composite. The properties utilized are geared towards the application of a continuum damage mechanics-based composite constitutive model available within LS-DYNA. The developed model can be applied to conduct impact simulations of structures composed of triaxially braided composites. The advantage of this technology is that it facilitates the analysis of the deformation and damage response of a triaxially braided polymer matrix composite within the environment of a transient dynamic finite-element code such as LS-DYNA in a manner which accounts for the local physical mechanisms but is still computationally efficient. This methodology is tightly coupled to experimental tests on the braided composite, which ensures that the material properties have physical significance. Aerospace or automotive companies interested in using triaxially braided composites in their structures, particularly for impact or crash applications, would find the technology useful. By the development of improved design tools, the amount of very expensive impact testing that will need to be performed can be significantly reduced.
The ALESEP program for the analysis of the inviscid/viscous interaction which occurs due to the presence of a closed laminar transitional separation bubble on an airflow is presented. The ALESEP code provides a iterative solution of the boundary layer equations expressed in an inverse formulation coupled to a Cauchy integral representation of the inviscid flow. This interaction analysis is treated as a local perturbation to a known solution obtained from a global airfoil analysis. Part of the required input to the ALESEP code are the reference displacement thickness and tangential velocity distributions. Special windward differencing may be used in the reversed flow regions of the separation bubble to accurately account for the flow direction in the discretization of the streamwise convection of momentum. The ALESEP code contains a forced transition model based on a streamwise intermittency function and a natural transition model based on a solution of the integral form of the turbulent kinetic energy equation. Instructions for the input/output, and program usage are presented.
A program called ALESEP is presented for the analysis of the inviscid-viscous interaction which occurs due to the presence of a closed laminar-transitional separation bubble on an airfoil or infinite swept wing. The ALESEP code provides an iterative solution of the boundary layer equations expressed in an inverse formulation coupled to a Cauchy integral representation of the inviscid flow. This interaction analysis is treated as a local perturbation to a known solution obtained from a global airfoil analysis; hence, part of the required input to the ALESEP code are the reference displacement thickness and tangential velocity distributions. Special windward differencing may be used in the reversed flow regions of the separation bubble to accurately account for the flow direction in the discretization of the streamwise convection of momentum. The ALESEP code contains a forced transition model based on a streamwise intermittency function, a natural transition model based on a solution of the integral form of the turbulent kinetic energy equation, and an empirical natural transition model.
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.
Program NSEG is a rapid mission analysis code based on the use of approximate flight path equations of motion. Equation form varies with the segment type, for example, accelerations, climbs, cruises, descents, and decelerations. Realistic and detailed vehicle characteristics are specified in tabular form. In addition to its mission performance calculation capabilities, the code also contains extensive flight envelope performance mapping capabilities. For example, rate-of-climb, turn rates, and energy maneuverability parameter values may be mapped in the Mach-altitude plane. Approximate take off and landing analyses are also performed. At high speeds, centrifugal lift effects are accounted for. Extensive turbojet and ramjet engine scaling procedures are incorporated in the code.
A numerical code is developed for evolutionary calculations of the thermal structure and composition of a porous comet nucleus made of water ice, in amorphous or crystalline form, other volatiles, dust, and gases trapped in amorphous ice. Bulk evaporation, crystallization, gas release, and free (Knudsen) flow of gases through the pores are taken into account. The numerical scheme yields exact conservation laws for mass and energy. The code is used to study the effect of bulk evaporation of ice in the interior of a comet nucleus during crystallization. It is found that evaporation controls the temperature distribution; the vapor prevents cooling of the crystallized layer of ice, by recondensation and release of latent heat. Thus high temperatures are maintained below the surface of the nucleus and down to depths of tens or hundreds of meters, even at large heliocentric distances, as long as crystallization goes on. Gas trapped in the ice and released during the phase transition flows both toward the interior and toward the surface and out of the nucleus. The progress of crystallization is largely determined by the contribution of gas fluxes to heat transfer.
The utility of multidisciplinary analysis tools for aeropropulsion applications is being investigated at the NASA Lewis Research Center. The goal of this project is to apply Spectrum, a multiphysics code developed by Centric Engineering Systems, Inc., to simulate multidisciplinary effects in turbomachinery components. Many engineering problems today involve detailed computer analyses to predict the thermal, aerodynamic, and structural response of a mechanical system as it undergoes service loading. Analysis of aerospace structures generally requires attention in all three disciplinary areas to adequately predict component service behavior, and in many cases, the results from one discipline substantially affect the outcome of the other two. There are numerous computer codes currently available in the engineering community to perform such analyses in each of these disciplines. Many of these codes are developed and used in-house by a given organization, and many are commercially available. However, few, if any, of these codes are designed specifically for multidisciplinary analyses. The Spectrum code has been developed for performing fully coupled fluid, thermal, and structural analyses on a mechanical system with a single simulation that accounts for all simultaneous interactions, thus eliminating the requirement for running a large number of sequential, separate, disciplinary analyses. The Spectrum code has a true multiphysics analysis capability, which improves analysis efficiency as well as accuracy. Centric Engineering, Inc., working with a team of Lewis and AlliedSignal Engines engineers, has been evaluating Spectrum for a variety of propulsion applications including disk quenching, drum cavity flow, aeromechanical simulations, and a centrifugal compressor flow simulation.
A multidimensional numerical simulation framework aimed at investigating the process of plasma detachment from a magnetic nozzle is introduced. An existing numerical code based on a magnetohydrodynamic formulation of the plasma flow equations that accounts for various dispersive and dissipative processes in plasmas was significantly enhanced to allow for the modeling of axisymmetric domains containing three.dimensiunai momentum and magnetic flux vectors. A separate magnetostatic solver was used to simulate the applied magnetic field topologies found in various nozzle experiments. Numerical results from a magnetic diffusion test problem in which all three components of the magnetic field were present exhibit excellent quantitative agreement with the analytical solution, and the lack of numerical instabilities due to fluctuations in the value of del(raised dot)B indicate that the conservative MHD framework with dissipative effects is well-suited for multi-dimensional analysis of magnetic nozzles. Further studies will focus on modeling literature experiments both for the purpose of code validation and to extract physical insight regarding the mechanisms driving detachment.
Work was conducted over a ten-year period to address the central issue of damage in primary load-bearing aircraft composite structure, specifically fuselage structure. This included the three facets of damage resistance, damage tolerance, and damage arrest. Experimental, analytical, and numerical work was conducted in order to identify and better understand the mechanisms that control the structural behavior of fuselage structures in their response to the three aspects of damage. Furthermore, work was done to develop straightforward design methodologies that can be employed by structural designers in preliminary design stages to make intelligent choices concerning the material, layup, and structural configurations so that a more efficient structure with structural integrity can be designed and built. Considerable progress was made towards achieving these goals via this work. In regard to damage tolerance considerations, the following were identified as important effects: composite layup and associated orthotropy/structural anisotropy, specifics of initial local damage mechanisms, role of longitudinal versus hoop stress, and large deformation and associated geometric nonlinearity. Means were established to account for effects of radius and for the nonlinear response. In particular, nondimensional parameters were identified to characterize the importance of nonlinearity in the response of pressurized cylinders. This led to the establishment of a iso-nonlinear-error plot for reference in structural design. Finally, in the case of damage tolerance, the general approach of the original methodology to predict the failure pressure involving extending basic plate failure data by accounting for the local stress intensification was accomplished for the general case by accounting for the mechanisms noted by utilizing the capability of the STAGS finite element code and numerically calculating the local stress intensification for the particular configuration to be considered. For the issue of damage arrest, placement of and configuration of stiffeners (including stiffener curvature), and magnitude and orientation of principal strains due to local bending were found to be key considerations. Means were established to account for stiffener effectiveness quantitatively based on radius, slit size, stiffener curvature' and relative bending stifffiesses involved. Geometric nonlinearity was also found to play an - 24 - important role here. Furthermore, it was determined that damage propagation is controlled by different mechanisms (hoop stress versus flapping stress and the associated factors involved in each) depending upon the direction of damage propagation. This latter item results in an inability to scale these phenomena in one test due to the different factors involved. Finally, the importance of shell curvature and associated instability in response to transverse loading including impact were found to be important considerations in damage resistance. A technique, involving asymmetric meshing of a finite element mesh, was developed to predict this behavior and showed excellent correlation with experimental results. Further details of these ten years of work are presented herein with references made to the fourteen documents produced during this work where full details can be found. Implications of this work are discussed and recommendations made. Although it is clear that there is more work to be done to fully understand composite fuselage technology and specifically the overall issue of damage in primary load-bearing composite structures, important understanding and capability has been extended via this work.
Prior studies implying significant benefits from using Navstar/GPS in civil aviation are seen as meriting further technical evaluation. While an uncorrupted C/A code and a P-code are seen as highly desirable in civil aviation, the differential GPS using C/A code only that is considered here and that also takes into account rather significant intentional degradation is found to hold promise. It is noted that such a technique may even approach the requirements associated with the conventional ILS window. It is stressed that since the receiver sequencing period/update rate appears to be an important parameter in differential GPS, sophisticated algorithms should be developed and evaluated in order to determine the full potential of the differential technique.