Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite temperature”

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 325 records · Page 18

Shear flow beneath oceanic plates - Local nonsimilarity boundary layers for olivine rheology

The principle of local similarity, which has been used to model the two-dimensional boundary layers in the oceanic upper mantle, permits calculation of the temperature, velocity, and stress fields with essentially analytic techniques. Finite difference numerical methods are hard pressed to resolve the detail required by the large variation of viscosity between the lithosphere and the asthenosphere. In this paper the local similarity approximation has been justified by quantitatively evaluating the effect of nonsimilarity due to viscous heating, nonlinear temperature- and pressure-dependent rheology, buoyancy, adiabatic cooling, etc. Nonsimilar effects produce only small modifications of the locally similar boundary layers; important geophysical observables such as surface heat flux and ocean floor topography are given to better than 10 percent by the locally similar solution. A posteriori evaluations of the terms neglected in the boundary layer simplification of the complete equations have been conducted on the locally similar temperature and velocity profiles close to the spreading ridge. The boundary layer models are valid to depths of 100 km at 3 m.y. and 10 km at 0.3 m.y.

Yuen, D. A.↗

The effects of uneven fiber spacing on thermal residual stresses in a unidirectional SCS-6Ti-15-3 laminate

High residual stresses develop in SCS-6/Ti-15-3 composites during cooldown from the fabrication temperature; these residual stresses can effect the mechanical and physical properties of the composite. Discrete fiber-matrix finite element models were used to study the residual stresses due to the temperature change during the fabrication process, including the effects of uneven fiber spacing, the free surface, and increased fiber volume fractions. To accurately model the effects of the free surface, it is only necessary to model one fiber through the thickness. Below the first ply, the analysis predicts stress distributions that are identical to the infinite array predictions. For uneven fiber space less then 0.042 mm in an interior ply, the maximum hoop stress was predicted to occur between fibers within a ply and to increase as the fiber spacing decreased. The maximum hoop stress correlated well with the observed radial cracking between fibers. For the case of touching fibers, the analysis predicted tensile radial stresses at the fiber-matrix debonding during the fabrication cooldown. Identical trends were predicted for uneven fiber spacing in surface plies with slightly greater values of maximum stresses. The analysis predicted matrix yielding to occur upon cooldown when the edge-to-edge fiber spacing was less than or equal to 0.022 mm. The stress distributions predicted for increasing fiber volume fractions were similar to those predicted for decreasing the fiber spacing for two adjacent fibers within a ply.

Bigelow, C. A.↗

The effects of uneven fiber spacing on thermal residual stresses in a unidirectional SCS-6/Ti-15-3 laminate

High residual stresses develop in SCS-6/Ti-15-3 composites during cooldown from the fabrication temperature; these residual stresses can effect the mechanical and physical properties of the composite. Discrete fiber-matrix finite element models were used to study the residual stresses due to the temperature change during the fabrication process, including the effects of uneven fiber spacing, the free surface, and increased fiber volume fractions. To accurately model the effects of the free surface, it is only necessary to model one fiber through the thickness. Below the first ply, the analysis predicts stress distributions that are identical to the infinite array predictions. For uneven fiber space less than 0.042 mm in an interior ply, the maximum hoop stress was predicted to occur between fibers within a ply and to increase as the fiber spacing decreased. The maximum hoop stress correlated well with the observed radial cracking between fibers. For the case of touching fibers, the analysis predicted tensile radial stresses at the fiber-matrix debonding during the fabrication cooldown. Identical trends were predicted for uneven fiber spacing in surface plies with slightly greater values of maximum stresses. The analysis predicted matrix yielding to occur upon cooldown when the edge-to-edge fiber spacing was less than or equal to 0.022 mm. The stress distributions preidcted for increasing fiber volume fractions were similar to those predicted for decreasing the fiber spacing for two adjacent fibers within a ply.

Bigelow, Catherine A.↗

Icing Formation Strains on Airfoil Skin

Aircraft icing has long presented a safety hazard to the aerospace industry. One problem inherent to aircraft icing is predicting the adhesive properties of the ice-aircraft interface; values commonly reported span multiple orders of magnitude, making shedding prediction and surface evaluation difficult. An effort has been undertaken at the NASA Glenn Research Center to understand the cause of the scatter in the data. Geometry dependencies, variation in testing procedures, and residual stresses in ice samples have previously been identified as significant contributors to the problem. Residual stresses were singled out for this study, which investigated the strains at the leading edge of a Common Research Model (CRM) wind tunnel model in the NASA Glenn Icing Research Tunnel. Sprays were conducted at a particle size median volumetric diameter (MVD) of 20 µm, with liquid water contents (LWC) spanning from 0.4 to 3.0 g/m3. Total temperatures were varied between 0 and -25 °C, and test section velocity between 51 and 118 m/s. The leading edge was instrumented with custom strain gages, resistance temperature detectors (RTDs), optical fibers, and a pressure mat. Results indicate Δ strains, which were the difference between recorded strains and those predicted from isotropic thermal expansion of aluminum, below -400 µε with consistent patterns of strain across a wide variety of tunnel conditions. Isotropic, linear elastic finite element analysis (FEA) was conducted by varying the temperature of simplified airfoil and ice geometries, producing similar strains to those observed in testing. This allowed induced stresses to be predicted. Predicted stresses were above the expected tensile strength for ice, and suggest a source for previously observed cracking in tunnel ice.

Ice Adhesion↗

A Computational Fluid Dynamic and Heat Transfer Model for Gaseous Core and Gas Cooled Space Power and Propulsion Reactors

A computational model based on the axisymmetric, thin-layer Navier-Stokes equations is developed to predict the convective, radiation and conductive heat transfer in high temperature space nuclear reactors. An implicit-explicit, finite volume, MacCormack method in conjunction with the Gauss-Seidel line iteration procedure is utilized to solve the thermal and fluid governing equations. Simulation of coolant and propellant flows in these reactors involves the subsonic and supersonic flows of hydrogen, helium and uranium tetrafluoride under variable boundary conditions. An enthalpy-rebalancing scheme is developed and implemented to enhance and accelerate the rate of convergence when a wall heat flux boundary condition is used. The model also incorporated the Baldwin and Lomax two-layer algebraic turbulence scheme for the calculation of the turbulent kinetic energy and eddy diffusivity of energy. The Rosseland diffusion approximation is used to simulate the radiative energy transfer in the optically thick environment of gas core reactors. The computational model is benchmarked with experimental data on flow separation angle and drag force acting on a suspended sphere in a cylindrical tube. The heat transfer is validated by comparing the computed results with the standard heat transfer correlations predictions. The model is used to simulate flow and heat transfer under a variety of design conditions. The effect of internal heat generation on the heat transfer in the gas core reactors is examined for a variety of power densities, 100 W/cc, 500 W/cc and 1000 W/cc. The maximum temperature, corresponding with the heat generation rates, are 2150 K, 2750 K and 3550 K, respectively. This analysis shows that the maximum temperature is strongly dependent on the value of heat generation rate. It also indicates that a heat generation rate higher than 1000 W/cc is necessary to maintain the gas temperature at about 3500 K, which is typical design temperature required to achieve high efficiency in the gas core reactors. The model is also used to predict the convective and radiation heat fluxes for the gas core reactors. The maximum value of heat flux occurs at the exit of the reactor core. Radiation heat flux increases with higher wall temperature. This behavior is due to the fact that the radiative heat flux is strongly dependent on wall temperature. This study also found that at temperature close to 3500 K the radiative heat flux is comparable with the convective heat flux in a uranium fluoride failed gas core reactor.

Anghaie, S.↗

The NATA code: Theory and analysis, volume 1

A computer program for calculating quasi-one-dimensional gas flow in axisymmetric and two-dimensional nozzles and rectangular channels is presented. Flow is assumed to start from a state of thermochemical equilibrium at a high temperature in an upstream reservoir. The program provides solutions based on frozen chemistry, chemical equilibrium, and nonequilibrium flow with finite reaction rates. Electronic nonequilibrium effects can be included using a two-temperature model. An approximate laminar boundary layer calculation is given for the shear and heat flux on the nozzle wall. Boundary layer displacement effects on the inviscid flow are considered also. Chemical equilibrium and transport property calculations are provided by subroutines. The code contains precoded thermochemical, chemical kinetic, and transport cross section data for high-temperature air, CO2-N2-Ar mixtures, helium, and argon. It provides calculations of the stagnation conditions on axisymmetric or two-dimensional models, and of the conditions on the flat surface of a blunt wedge. The primary purpose of the code is to describe the flow conditions and test conditions in electric arc heated wind tunnels.

Bade, W. L.↗

Analysis of film cooling in rocket nozzles

Progress during the reporting period is summarized. Analysis of film cooling in rocket nozzles by computational fluid dynamics (CFD) computer codes is desirable for two reasons. First, it allows prediction of resulting flow fields within the rocket nozzle, in particular the interaction of the coolant boundary layer with the main flow. This facilitates evaluation of potential cooling configurations with regard to total thrust, etc., before construction and testing of any prototype. Secondly, CFD simulation of film cooling allows for assessment of the effectiveness of the proposed cooling in limiting nozzle wall temperature rises. This latter objective is the focus of the current work. The desired objective is to use the Finite Difference Navier Stokes (FDNS) code to predict wall heat fluxes or wall temperatures in rocket nozzles. As prior work has revealed that the FDNS code is deficient in the thermal modeling of boundary conditions, the first step is to correct these deficiencies in the FDNS code. Next, these changes must be tested against available data. Finally, the code will be used to model film cooling of a particular rocket nozzle. The third task of this research, using the modified code to compute the flow of hot gases through a nozzle, is described.

Woodbury, Keith A.↗

Numerical modeling and passive thermal control of external lighting systems for Space Station Freedom

Consideration is given to three generic families of luminaries with lamp power ranging from 11 to 150 watts. A concept of an equivalent radiation node boundary temperature was used to impose worst hot and cold environments, and transient finite difference models were developed to study the effects of geometry and optical properties of thermal control coatings. Minimum and maximum transient temperatures were computed at the critical location during 90 minute orbit and were compared with allowable limits. Results show that with the proper choice of optical properties, the luminaries can be passively controlled to within acceptable limits.

Buck, Gregory A.↗

Solar Sail Topology Variations Due to On-Orbit Thermal Effects

The objective of this research was to predict the influence of non-uniform temperature distribution on solar sail topology and the effect of such topology variations on sail performance (thrust, torque). Specifically considered were the thermal effects due to on orbit attitude control maneuvers. Such maneuvers are expected to advance the sail to a position off-normal to the sun by as much as 35 degrees; a solar sail initially deformed by typical pre-tension and solar pressure loads may suffer significant thermally induced strains due to the non-uniform heating caused by these maneuvers. This on-orbit scenario was investigated through development of an automated analytical shape model that iterates many times between sail shape and sail temperature distribution before converging on a final coupled thermal structural affected sail topology. This model utilizes a validated geometrically non-linear finite element model and a thermal radiation subroutine. It was discovered that temperature gradients were deterministic for the off-normal solar angle cases as were thermally induced strains. Performance effects were found to be moderately significant but not as large as initially suspected. A roll torque was detected, and the sail center of pressure shifted by a distance that may influence on-orbit sail control stability.

Banik, Jeremy A.↗

Effects of forced and free convections on structural temperatures of Space Shuttle Orbiter during reentry flight

Structural performance and resizing (SPAR) finite-element thermal analysis computer program was used in the heat transfer analysis of the Space Shuttle Orbiter wing subjected to reentry aerodynamic heating. With sufficient external forced convective cooling near the end of the heating cycle, the calculated surface temperatures of the thermal protection system (TPS) agree favorably with the flight data for the entire flight profile. However, the effects of this external forced convective cooling on the structural temperatures were found to be negligible. Both free convection and forced convection elements were introduced to model the internal convection effect of the cool air entering the Shuttle interior. The introduction of the internal free convection effect decreased the calculated wing lower skin temperatures by 20 F at most, 1200 sec after touchdown. If the internal convection is treated as forced convection, the calculated wing lower skin temperatures after touchdown can be reduced to match the flight-measured data very closely. By reducing the TPS thicknesses to certain effective thickness to account for the TPs gap heating, the calculated wing lower skin temperatures prior to touchdown can be raised to agree with the flight data perfectly.

Ko, William L.↗

Creep Life Prediction of Ceramic Components Using the Finite Element Based Integrated Design Program (CARES/Creep)

The desirable properties of ceramics at high temperatures have generated interest in their use for structural applications such as in advanced turbine systems. Design lives for such systems can exceed 10,000 hours. Such long life requirements necessitate subjecting the components to relatively low stresses. The combination of high temperatures and low stresses typically places failure for monolithic ceramics in the creep regime. The objective of this work is to present a design methodology for predicting the lifetimes of structural components subjected to multiaxial creep loading. This methodology utilizes commercially available finite element packages and takes into account the time varying creep stress distributions (stress relaxation). In this methodology, the creep life of a component is divided into short time steps, during which, the stress and strain distributions are assumed constant. The damage, D, is calculated for each time step based on a modified Monkman-Grant creep rupture criterion. For components subjected to predominantly tensile loading, failure is assumed to occur when the normalized accumulated damage at any point in the component is greater than or equal to unity.

Jadaan, Osama M.↗

Creep Life of Ceramic Components Using a Finite-Element-Based Integrated Design Program (CARES/CREEP)

The desirable properties of ceramics at high temperatures have generated interest in their use for structural application such as in advanced turbine engine systems. Design lives for such systems can exceed 10,000 hours. The long life requirement necessitates subjecting the components to relatively low stresses. The combination of high temperatures and low stresses typically places failure for monolithic ceramics in the creep regime. The objective of this paper is to present a design methodology for predicting the lifetimes of structural components subjected to creep rupture conditions. This methodology utilizes commercially available finite element packages and takes into account the time-varying creep strain distributions (stress relaxation). The creep life, of a component is discretized into short time steps, during which the stress and strain distributions are assumed constant. The damage is calculated for each time step based on a modified Monkman-Grant creep rupture criterion. Failure is assumed to occur when the normalized accumulated damage at any point in the component is greater than or equal to unity. The corresponding time will be the creep rupture life for that component. Examples are chosen to demonstrate the Ceramics Analysis and Reliability Evaluation of Structures/CREEP (CARES/CREEP) integrated design program, which is written for the ANSYS finite element package. Depending on the component size and loading conditions, it was found that in real structures one of two competing failure modes (creep or slow crack growth) will dominate. Applications to benchmark problems and engine components are included.

Powers, L. M.↗

Creep Life of Ceramic Components Using a Finite-Element-Based Integrated Design Program (CARES/CREEP)

The desirable properties of ceramics at high temperatures have generated interest in their use for structural applications such as in advanced turbine systems. Design lives for such systems can exceed 10,000 hours. The long life requirement necessitates subjecting the components to relatively low stresses. The combination of high temperatures and low stresses typically places failure for monolithic ceramics in the creep regime. The objective of this paper is to present a design methodology for predicting the lifetimes of structural components subjected to creep rupture conditions. This methodology utilized commercially available finite element packages and takes into account the time-varying creep strain distributions (stress relaxation). The creep life of a component is discretized into short time steps, during which the stress and strain distributions are assumed constant. The damage is calculated for each time step based on a modified Monkman-Grant creep rupture criterion. Failure is assumed to occur when the normalized accumulated damage at any point in the component is greater than or equal to unity. The corresponding time will be the creep rupture life for that component. Examples are chosen to demonstrate the CARES/CREEP (Ceramics Analysis and Reliability Evaluation of Structures/CREEP) integrated design programs, which is written for the ANSYS finite element package. Depending on the component size and loading conditions, it was found that in real structures one of two competing failure modes (creep or slow crack growth) will dominate. Applications to benechmark problems and engine components are included.

Gyekenyesi, J. P.↗

Evaporating Spray in Supersonic Streams Including Turbulence Effects

Evaporating spray plays an important role in spray combustion processes. This paper describes the development of a new finite-conductivity evaporation model, based on the two-temperature film theory, for two-phase numerical simulation using Eulerian-Lagrangian method. The model is a natural extension of the T-blob/T-TAB atomization/spray model which supplies the turbulence characteristics for estimating effective thermal diffusivity within the droplet phase. Both one-way and two-way coupled calculations were performed to investigate the performance of this model. Validation results indicate the superiority of the finite-conductivity model in low speed parallel flow evaporating sprays. High speed cross flow spray results indicate the effectiveness of the T-blob/T-TAB model and point to the needed improvements in high speed evaporating spray modeling.

Balasubramanyam, M. S.↗

Combined Loads Test Fixture for Thermal-Structural Testing Aerospace Vehicle Panel Concepts

A structural test requirement of the National Aero-Space Plane (NASP) program has resulted in the design, fabrication, and implementation of a combined loads test fixture. Principal requirements for the fixture are testing a 4- by 4-ft hat-stiffened panel with combined axial (either tension or compression) and shear load at temperatures ranging from room temperature to 915 F, keeping the test panel stresses caused by the mechanical loads uniform, and thermal stresses caused by non-uniform panel temperatures minimized. The panel represents the side fuselage skin of an experimental aerospace vehicle, and was produced for the NASP program. A comprehensive mechanical loads test program using the new test fixture has been conducted on this panel from room temperature to 500 F. Measured data have been compared with finite-element analyses predictions, verifying that uniform load distributions were achieved by the fixture. The overall correlation of test data with analysis is excellent. The panel stress distributions and temperature distributions are very uniform and fulfill program requirements. This report provides details of an analytical and experimental validation of the combined loads test fixture. Because of its simple design, this unique test fixture can accommodate panels from a variety of aerospace vehicle designs.

Fields, Roger A.↗

NEXT Ion Thruster Thermal Model

As the NEXT ion thruster progresses towards higher technology readiness, it is necessary to develop the tools that will support its implementation into flight programs. An ion thruster thermal model has been developed for the latest prototype model design to aid in predicting thruster temperatures for various missions. This model is comprised of two parts. The first part predicts the heating from the discharge plasma for various throttling points based on a discharge chamber plasma model. This model shows, as expected, that the internal heating is strongly correlated with the discharge power. Typically, the internal plasma heating increases with beam current and decreases slightly with beam voltage. The second is a model based on a finite difference thermal code used to predict the thruster temperatures. Both parts of the model will be described in this paper. This model has been correlated with a thermal development test on the NEXT Prototype Model 1 thruster with most predicted component temperatures within 5 to 10 C of test temperatures. The model indicates that heating, and hence current collection, is not based purely on the footprint of the magnet rings, but follows a 0.1:1:2:1 ratio for the cathode-to-conical-to-cylindrical-to-front magnet rings. This thermal model has also been used to predict the temperatures during the worst case mission profile that is anticipated for the thruster. The model predicts ample thermal margin for all of its components except the external cable harness under the hottest anticipated mission scenario. The external cable harness will be re-rated or replaced to meet the predicted environment.

VanNoord, Jonathan L.↗

Three-dimensional vector modeling and restoration of flat finite wave tank radiometric measurements

In this paper, a three-dimensional Fourier transform inversion method describing the interaction between water surface emitted radiation from a flat finite wave tank and antenna radiation characteristics is reported. The transform technique represents the scanning of the antenna mathematically as a correlation. Computation time is reduced by using the efficient and economical fast Fourier transform algorithm. To verify the inversion method, computations have been made and compared with known data and other available results. The technique has been used to restore data of the finite wave tank system and other available antenna temperature measurements made at the Cape Cod Canal. The restored brightness temperatures serve as better representations of the emitted radiation than the measured antenna temperatures.

Truman, W. M.↗

An exact solution for the solidification of a liquid slab of binary mixture

The time dependent temperature and concentration profiles of a one dimensional finite slab of a binary liquid alloy is investigated during solidification. The governing equations are reduced to a set of coupled, nonlinear initial value problems using the method outlined by Meyer. Two methods will be used to solve these equations. The first method uses a Runge-Kutta-Fehlberg integrator to solve the equations numerically. The second method comprises of finding closed form solutions of the equations.

Antar, B. N.↗