Search NASASearch

SEARCH · Search NASA

Results for “energy efficient computing”

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 37 records · Page 2

Computational Investigation of Novel Tip Leakage Mitigation Methods for High Pressure Turbine Blades

This paper presents preliminary findings on a possible approach to reducing tip leakage losses. In this paper a computational study was conducted on the Energy Efficient Engine (EEE) High Pressure Turbine (HPT) rotor tip geometry using the commercial numerical solver ANSYS FLUENT. The flow solver was validated against aerodynamic data acquired in the NASA Transonic Turbine Blade Cascade facility. The scope of the ongoing study is to computationally investigate how the tip leakage and overall blade losses are affected by (1) injection from the tip near the pressure side, (2) injection from the tip surface at the camber line, and (3) injection from the tip surface into the tip separation bubble. The objective is to identify the locations on the tip surface at which to place appropriately configured blowing keeping in mind the film cooling application of tip blowing holes. The validation was conducted at Reynolds numbers of 85,000, 343,000, and 685,000 and at engine realistic flow conditions. The coolant injection simulations were conducted at a Reynolds number of 343,000 based on blade chord and inlet velocity and utilized the SST turbulence model in FLUENT. The key parameters examined are the number of jets, jet angle and jet location. A coolant to inlet pressure ratio of 1.0 was studied for angles of +30 deg, -30 deg, and 90 deg to the local free stream on the tip. For the 3 hole configuration, 3 holes spaced 3 hole diameters apart with length to diameter ratio of 1.5 were used. A simulation including 11 holes along the entire mean camber line is also presented (30 deg toward suction side). In addition, the effect of a single hole is also compared to a flat tip with no injection. The results provide insight into tip flow control methods and can be used to guide further investigation into tip flow control. As noted in past research it is concluded that reducing leakage flow is not necessarily synonymous with reducing losses due to leakage.

Ibrahim, Mounir

Computational Investigation of Novel Tip Leakage Mitigation Methods for High Pressure Turbine Blades

This paper presents preliminary findings on a possible approach to reducing tip leakage losses. In this paper a computational study was conducted on the EEE (Energy Efficient Engine) HPT (High Pressure Turbine) rotor tip geometry using the commercial numerical solver ANSYS FLUENT. The flow solver was validated against aerodynamic data acquired in the NASA Transonic Turbine Blade Cascade facility. The scope of the ongoing study is to computationally investigate how the tip leakage and overall blade losses are affected by 1. injection from the tip near the pressure side, 2. injection from the tip surface at the camber line, and 3. injection from the tip surface into the tip separation bubble. The objective is to identify the locations on the tip surface at which to place appropriately configured blowing keeping in mind the film cooling application of tip blowing holes. The validation was conducted at Reynolds numbers of 85,000, 343,000 and 685,000 and at engine realistic flow conditions. The coolant injection simulations were conducted at a Reynolds number of 343,000 based on blade chord and inlet velocity and utilized the SST turbulence model in FLUENT. The key parameters examined are the number of jets, jet angle and jet location. A coolant to inlet pressure ratio of 1.0 was studied for angles of +30 deg., -30 deg. and 90 deg. to the local free stream on the tip. For the 3 hole configuration, 3 holes spaced 3 hole diameters apart with length to diameter ratio of 1.5 were used. A simulation including 11 holes along the entire mean camber line is also presented (30 degrees toward suction side). In addition, the effect of a single hole is also compared to a flat tip with no injection. The results provide insight into tip flow control methods and can be used to guide further investigation into tip flow control. As noted in past research it is concluded that reducing leakage flow is not necessarily synonymous with reducing losses due to leakage.

Ibrahim, Mounir

Rapid calculation of radiative heating rates and photodissociation rates in inhomogeneous multiple scattering atmospheres

The solution of the generalized two-stream approximation for radiative transfer in homogeneous multiple scattering atmospheres is extended to vertically inhomogeneous atmospheres in a manner which is numerically stable and computationally efficient. It is shown that solar energy deposition rates, photolysis rates, and infrared cooling rates all may be calculated with the simple modifications of a single algorithm. The accuracy of the algorithm is generally better than 10 percent, so that other uncertainties, such as in absorption coefficients, may often dominate the error in calculation of the quantities of interest to atmospheric studies.

Toon, Owen B.

The role of energetic particle precipitation in Jovian magnetospherics. I - Secondary electrons from the ionosphere of Jupiter

Escape of photoelectrons from the Jovian ionosphere produces only a meager source of thermal plasma for the inflated centrifugally unstable magnetosphere and is unable to account for the 100 eV thermal plasma temperatures of the magnetosphere. Since the Jovian magnetosphere is well populated with highly energetic electrons, the creation of secondaries and the energy degradation of the primaries precipitating into the lower ionosphere provide additional sources of magnetospheric thermal plasma as well as of the energy for further elevating the plasma temperatures in the top side ionosphere. The efficiencies with which escaping electrons are created by precipitating electrons with energies up to several MeV are computed using energetic electron transport and thermalization codes. The more energetic incident fluxes are far less efficient in creating escaping electrons than the lower energy fluxes with only 0.001% of the secondaries escaping for a 1 MeV source versus 3% for a 1 keV source. Incident fluxes of the order of 10 to 100 per sq cm per sec per eV between 100 eV and 100 keV are required to produce 50 eV escape fluxes comparable to those generated by solar EUV.

Swartz, W. E.

Photocell heat engine solar power systems

A combined photocell heat engine concept is proposed for high efficiency solar energy conversion in space. In this concept the short wavelength portion of the solar spectrum is split by a dichroic filter and sent to a bank of photocells. The long wave-length remainder of the spectrum is used by the heat engine. This technique allows the photocells to operate with the minimum amount of waste heat, increasing their efficiency and reducing the amount of cooling required. The heat engine operates by direct absorption in a working fluid containing broadband absorber molecules or particulates. A window in the heat engine admits the long wave-lengths from the solar spectrum. The window may also reflect a portion of the internal gaseous reradiation spectrum (e.g., a heat mirror) to help reduce radiation losses. Flow-induced thermal gradients may also reduce reradiation losses in the case of optically thick working fluids. The efficiencies computed for the photocell heat engine solar energy converter can be as high as 42 percent.

Taussig, R. T.

Phase change energy storage for solar dynamic power systems

This paper presents the results of a transient computer simulation that was developed to study phase change energy storage techniques for Space Station Freedom (SSF) solar dynamic (SD) power systems. Such SD systems may be used in future growth SSF configurations. Two solar dynamic options are considered in this paper: Brayton and Rankine. Model elements consist of a single node receiver and concentrator, and takes into account overall heat engine efficiency and power distribution characteristics. The simulation not only computes the energy stored in the receiver phase change material (PCM), but also the amount of the PCM required for various combinations of load demands and power system mission constraints. For a solar dynamic power system in low earth orbit, the amount of stored PCM energy is calculated by balancing the solar energy input and the energy consumed by the loads corrected by an overall system efficiency. The model assumes an average 75 kW SD power system load profile which is connected to user loads via dedicated power distribution channels. The model then calculates the stored energy in the receiver and subsequently estimates the quantity of PCM necessary to meet peaking and contingency requirements. The model can also be used to conduct trade studies on the performance of SD power systems using different storage materials.

Chiaramonte, F. P.

Constraining Cometary Crystal Shapes from IR Spectral Features

A major challenge in deriving the silicate mineralogy of comets is ascertaining how the anisotropic nature of forsterite crystals affects the spectral features' wavelength, relative intensity, and asymmetry. Forsterite features are identified in cometary comae near 10, 11.05-11.2, 16, 19, 23.5, 27.5 and 33 microns [1-10], so accurate models for forsterite's absorption efficiency (Qabs) are a primary requirement to compute IR spectral energy distributions (SEDs, lambdaF lambda vs. lambda) and constrain the silicate mineralogy of comets. Forsterite is an anisotropic crystal, with three crystallographic axes with distinct indices of refraction for the a-, b-, and c-axis. The shape of a forsterite crystal significantly affects its spectral features [13-16]. We need models that account for crystal shape. The IR absorption efficiencies of forsterite are computed using the discrete dipole approximation (DDA) code DDSCAT [11,12]. Starting from a fiducial crystal shape of a cube, we systematically elongate/reduce one of the crystallographic axes. Also, we elongate/reduce one axis while the lengths of the other two axes are slightly asymmetric (0.8:1.2). The most significant grain shape characteristic that affects the crystalline spectral features is the relative lengths of the crystallographic axes. The second significant grain shape characteristic is breaking the symmetry of all three axes [17]. Synthetic spectral energy distributions using seven crystal shape classes [17] are fit to the observed SED of comet C/1995 O1 (Hale-Bopp). The Hale-Bopp crystalline residual better matches equant, b-platelets, c-platelets, and b-columns spectral shape classes, while a-platelets, a-columns and c-columns worsen the spectral fits. Forsterite condensation and partial evaporation experiments demonstrate that environmental temperature and grain shape are connected [18-20]. Thus, grain shape is a potential probe for protoplanetary disk temperatures where the cometary crystalline forsterite formed. The forsterite crystal shapes (equant, b-platelets, c-platelets, b-columns - excluding a- and c-columns) derived from our modeling [17] of comet Hale- Bopp, compared to laboratory synthesis experiments [18], suggests that these crystals are high temperature condensates. By observing and modeling the crystalline features in comet ISON, we may constrain forsterite crystal shape(s) and link to their formation temperature(s) and environment(s).

anisotropic crystal

Pragmatic Stress Prediction on Additively Manufactured Coupons

Prediction of residual stresses from process parameters for additively manufactured large metal parts is computationally expensive. NASA is currently developing meter-scale parts with direct energy deposition. Practically, the predictive computational methods need to efficiently scale-up to meter-scale parts. Coupled thermal-mechanical multi-physics simulations have been developed with the pragmatic method using ABAQUS, COMSOL Multiphysics, ALE3D software. The residual stresses are a result of the manufacturing process which creates thermal cycling of the build layers. The pragmatic method uses lumped thermal layers for stress predictions to reduce computational costs. The stress predictions as well as deformations of the different codes are compared with each other and with ANSYS Additive using identical material models, boundary and initial conditions. The codes were used to simulate three different geometries: a thin wall, hollow cylinder and twin-cantilever part. The coupon parts were then manufactured with Inconel-625. The residual stresses in these parts were measured using X-ray diffraction as well as neutron beam diffraction at NIST. The stress measurements for the two technologies are compared. The pragmatic stress prediction method enabled predictions of the multi-centimeter scale parts using desktop computer workstations in only a few hours for each coupon. The results of the simulated stress predictions compared favorably with the measured stresses even though thermally lumped layers were employed. Finally, a two-meter scale nozzle was simulated using ANSYS Additive. The simulations were used to examine the build orientation trade-space with respect to resulting geometric deformation. The predicted deformations were compared to measurements of an actual subscale part manufactured with direct energy deposition.

Additive Manufacturing

Pragmatic Stress Prediction on Additively Manufactured Coupons

Prediction of residual stresses from process parameters for additively manufactured large metal parts is computationally expensive. NASA is currently developing meter-scale parts with direct energy deposition. Practically, the predictive computational methods need to efficiently scale-up to meter-scale parts. Coupled thermal-mechanical multi-physics simulations have been developed with the pragmatic method using ABAQUS, COMSOL Multiphysics, ALE3D software. The residual stresses are a result of the manufacturing process which creates thermal cycling of the build layers. The pragmatic method uses lumped thermal layers for stress predictions to reduce computational costs. The stress predictions as well as deformations of the different codes are compared with each other and with ANSYS Additive using identical material models, boundary and initial conditions. The codes were used to simulate three different geometries: a thin wall, hollow cylinder and twin-cantilever part. The coupon parts were then manufactured with Inconel-625. The residual stresses in these parts were measured using X-ray diffraction as well as neutron beam diffraction at NIST. The stress measurements for the two technologies are compared. The pragmatic stress prediction method enabled predictions of the multi-centimeter scale parts using desktop computer workstations in only a few hours for each coupon. The results of the simulated stress predictions compared favorably with the measured stresses even though thermally lumped layers were employed. Finally, a two-meter scale nozzle was simulated using ANSYS Additive. The simulations were used to examine the build orientation trade-space with respect to resulting geometric deformation. The predicted deformations were compared to measurements of an actual subscale part manufactured with direct energy deposition.

pragmatic method

Use of Transition Modeling to Enable the Computation of Losses for Variable-Speed Power Turbine

To investigate the penalties associated with using a variable speed power turbine (VSPT) in a rotorcraft capable of vertical takeoff and landing, various analysis tools are required. Such analysis tools must be able to model the flow accurately within the operating envelope of VSPT. For power turbines low Reynolds numbers and a wide range of the incidence angles, positive and negative, due to the variation in the shaft speed at relatively fixed corrected flows, characterize this envelope. The flow in the turbine passage is expected to be transitional and separated at high incidence. The turbulence model of Walters and Leylek was implemented in the NASA Glenn-HT code to enable a more accurate analysis of such flows. Two-dimensional heat transfer predictions of flat plate flow and two-dimensional and three-dimensional heat transfer predictions on a turbine blade were performed and reported herein. Heat transfer computations were performed because it is a good marker for transition. The final goal is to be able to compute the aerodynamic losses. Armed with the new transition model, total pressure losses for three-dimensional flow of an Energy Efficient Engine (E3) tip section cascade for a range of incidence angles were computed in anticipation of the experimental data. The results obtained form a loss bucket for the chosen blade.

Ameri, Ali A.

Digital Technologies at NASA for Science and Engineering

While scientific and engineering advancements used to rely primarily on theoretical studies and physical experiments, today digital technology enabled by petaflops-scale supercomputers is an equal, if not a greater, contributor to such achievements. In addition, computational modeling and simulation serves as a predictive tool that is not otherwise available. As a result, the use of high performance computing is integral to NASA's work in all mission areas such as space exploration, aeronautics, and scientific discovery. But traditional supercomputing alone is not sufficient for all of the space agency's needs. The success of many NASA missions depends on solving complex computing challenges, some of which are NP-hard (decision theory) if using classical solution methods. Quantum computing promises an unprecedented ability to solve such intractable problems by harnessing quantum mechanical effects such as tunneling, superposition, and entanglement. Another disruptive digital technology is neuromorphic computing that uses brain-inspired lessons to generate new architectures that are much more energy efficient, and capable of massive parallel processing and learning in-situ. Finally, with large amounts of observational and computational data sets, the opportunities of big data and data analytics can be leveraged to enable deep learning and knowledge discovery - it's all a massive digital transformation. This talk will be an overview how NASA utilizes digital technologies for its science and engineering efforts.

Biswas, Rupak

Effects of aircraft and flight parameters on energy-efficient profile descents in time-based metered traffic

Attention is given to a computer algorithm yielding the data required for a flight crew to navigate from an entry fix, about 100 nm from an airport, to a metering fix, and arrive there at a predetermined time, altitude, and airspeed. The flight path is divided into several descent and deceleration segments. Results for the case of a B-737 airliner indicate that wind and nonstandard atmospheric properties have a significant effect on the flight path and must be taken into account. While a range of combinations of Mach number and calibrated airspeed is possible for the descent segments leading to the metering fix, only small changes in the fuel consumed were observed for this range of combinations. A combination that is based on scheduling flexibility therefore seems preferable.

Dejarnette, F. R.

Electron energy deposition in carbon monoxide gas

A comprehensive set of electron impact cross sections for carbon monoxide molecules is presented on the basis of the most recent experimental measurements and theoretical calculations. The processes by which energetic electrons lose energy in CO gas are analyzed with these input cross sections. The efficiencies are computed of vibrational and electronic excitation, dissociation, ionization, and heating for CO gas with fractional ionization ranging from 0% to 10%. The calculated mean energy per ion pair for neutral CO gas is 32.3 eV, which is in excellent agreement with the experimental value of 32.2 eV. It increases to 35.6 eV at a fractional ionization of 1%, typical of supernovae ejecta.

Liu, Weihong

Bidirectional Reflectance of a Macroscopically Flat, High-Albedo Particulate Surface: An Efficient Radiative Transfer Solution and Applications to Regoliths

Many remote sensing applications rely on accurate knowledge of the bidirectional reflection function (BRF) of surfaces composed of discrete, randomly positioned scattering particles. Theoretical computations of BRFs for plane-parallel particulate layers are usually reduced to solving the radiative transfer equation (RTE) using one of existing exact or approximate techniques. Since semi-empirical approximate approaches are notorious for their low accuracy, violation of the energy conservation law, and ability to produce unphysical results, the use of numerically exact solutions of RTE has gained justified popularity. For example, the computation of BRFs for macroscopically flat particulate surfaces in many geophysical publications is based on the adding-doubling (AD) and discrete ordinate (DO) methods. A further saving of computer resources can be achieved by using a more efficient technique to solve the plane-parallel RTE than the AD and DO methods. Since many natural particulate surfaces can be well represented by the model of an optically semi-infinite, homogeneous scattering layer, one can find the BRF directly by solving the Ambartsumian's nonlinear integral equation using a simple iterative technique. In this way, the computation of the internal radiation field is avoided and the computer code becomes highly efficient and very accurate and compact. Furthermore, the BRF thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. In this paper, we discuss numerical aspects and the computer implementation of this technique, examine the applicability of the Henyey-Greenstein phase function and the sigma-Eddington approximation in BRF and flux calculations, and describe sample applications demonstrating the potential effect of particle shape on the bidirectional reflectance of flat regolith surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web, and can be applied to a wide range of remote sensing problems. BRF computations for undulated (macroscopically rough) surfaces are more complicated and often rely on time consuming Monte Carlo procedures. This approach is especially inefficient for optically thick, weakly absorbing media (e.g., snow and desert surfaces at visible wavelengths since a photon may undergo many internal scattering events before it exists the medium or is absorbed. However, undulated surfaces can often be represented as collections of locally flat tilted facets characterized by the BRF found from the traditional plane parallel RTE. In this way the MOnte Carlo procedure could be used only to evaluate the effects of surface shadowing and multiple surface reflections, thereby bypassing the time-consuming ray tracing inside the medium and providing a great savings of CPU time.

Mishchenko, Michael I.

How to quickly predict the overall TWT and the multistage depressed collector efficiency

The study deals with an empirical, simple formula extracted from a three-dimensional helical-TWT computer program that expresses the lowest energy in a spent beam in terms of beam perveance and electronic efficiency. The formula has a general validity down to 4 - 5 dB below saturation and gives 1 - delta V/V with less than 20% error down to 10 dB below saturation.

Kosmahl, H. G.

Computer Model Of Fragmentation Of Atomic Nuclei

High Charge and Energy Semiempirical Nuclear Fragmentation Model (HZEFRG1) computer program developed to be computationally efficient, user-friendly, physics-based program for generating data bases on fragmentation of atomic nuclei. Data bases generated used in calculations pertaining to such radiation-transport applications as shielding against radiation in outer space, radiation dosimetry in outer space, cancer therapy in laboratories with beams of heavy ions, and simulation studies for designing detectors for experiments in nuclear physics. Provides cross sections for production of individual elements and isotopes in breakups of high-energy heavy ions by combined nuclear and Coulomb fields of interacting nuclei. Written in ANSI FORTRAN 77.

Wilson, John W.

Computation of strain energy release rates for skin-stiffener debonds modeled with plate elements

An efficient method for predicting the strength of debonded composite skin-stiffener configurations is presented. This method, which is based on fracture mechanics, models the skin and the stiffener with two-dimensional (2D) plate elements instead of three-dimensional (3D) solid elements. The skin and stiffener flange nodes are tied together by two modeling techniques. In one technique, the corresponding flange and skin nodes are required to have identical translational and rotational degrees-of-freedom. In the other technique, the corresponding flange and skin nodes are only required to have identical translational degrees-of-freedom. Strain energy release rate formulas are proposed for both modeling techniques. These formulas are used for skin-stiffener debond cases with and without cylindrical bending deformations. The cylindrical bending results are compared with plane-strain finite element results. Excellent agreement between the two sets of results is obtained when the second technique is used. Thus, from these limited studies, a preferable modeling technique for skin-stiffener debond analysis using plate elements is established.

Wang, J. T.