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 433 records · Page 24

Comparison of Numerical Model Results with Diffusion Flames in Microgravity

The effects of gravity on methane diffusion flames were studied using a finite difference numerical model. Variables under consideration included velocity, pressure, temperature, enthalpy, reactive species, and inert carrier gas. Results obtained showed that flames in zero gravity were wider and taller than those in normal gravity. The use of a six reaction scheme produced similar predictions of flame width and height compared to those of the single reaction case. The numerical models predicted a considerable difference in the flame base region between normal gravity and low gravity. This region is considered to be beyond the boundaries of current analytical models.

Hyer, Paul V.↗

Burning of a spherical fuel droplet in a uniform flowfield with exact property variation

An analytical/numerical model is developed for single droplet evaporation and burning in a convective flowfield. The model is based on the boundary-layer approach, and chemical reaction kinetics are represented by a one-step, finite-rate reaction mechanism, while variation of gas properties with temperature and gas composition is based on the kinetic theory of gases. Four droplet models differing in the degree of complexity concerning property variation and chemistry are compared. Comparisons are also provided with existing empirical correlations for convective droplet evaporation and burning.

Madooglu, K.↗

Interlaminar shear test method development for long term durability testing of composites

The high speed civil transport is a commercial aircraft that is expected to carry 300 passengers at Mach 2.4 over a range of more than 6000 nautical miles. With the existing commercial structural material technology (i.e., aluminum) the performance characteristics of the high speed civil transport would not be realized. Therefore there has been a concerted effort in the development of light weight materials capable of withstanding elevated temperatures for long duration. Thermoplastic composite materials are such candidate materials and the understanding of how these materials perform over the long term under harsh environments is essential to safe and effective design. The matrix dominated properties of thermoplastic composites are most affected by both time and temperature. There is currently an effort to perform short term testing to predict long term behavior of in-plane mechanical properties E22 (transverse modulus of elasticity) and G12 (shear modulus). Out-of-plane properties such as E33, G13, and G23 are inherently more difficult to characterize. This is especially true for the out-of-plane shear modulus G23 and hence there is no existing acceptable standard test method. Since G23 is the most matrix dominated property, it is essential that a test method be developed. A shear test methodology is developed to do just that. The test method, called the double notched specimen, along with the previously developed shear gage was tested at room temperature. Mechanical testing confirmed the attributes of the methodology. A finite element parametric study was conducted for specimen optimization. Moire interferometry, a high sensitivity laser optical method, was used for full-field analysis of the specimen. From this work, material parameters will be determined and thus enable the prediction of long term material behavior of laminates subjected to general loading states.

Ifju, Peter G.↗

A PV view of the zonal mean distribution of temperature and wind in the extratropical troposphere

The dependence of the temperature and wind distribution of the zonal mean flow in the extratropical troposphere on the gradient of pontential vorticity along isentropes is examined. The extratropics here refer to the region outside the Hadley circulation. Of particular interest is whether the distribution of temperature and wind corresponding to a constant potential vorticity (PV) along isentropes resembles the observed, and the implications of PV homogenization along isentropes for the role of the tropics. With the assumption that PV is homogenized along isentropes, it is found that the temperature distribution in the extratropical troposphere may be determined by a linear, first-order partial differential equation. When the observed surface temperature distribution and tropical lapse rate are used as the boundary conditions, the solution of the equation is close to the observed temperature distribution except in the upper troposphere adjacent to the Hadley circulation, where the troposphere with no PV gradient is considerably colder. Consequently, the jet is also stronger. It is also found that the meridional distribution of the balanced zonal wind is very sensitive to the meridional distribution of the tropopause temperature. The result may suggest that the requirement of the global momentum balance has no practical role in determining the extratropical temperature distribution. The authors further investigated the sensitivity of the extratropical troposphere with constant PV along isentropes to changes in conditions at the tropical boundary (the edge of the Hadley circulation). It is found that the temperature and wind distributions in the extratropical troposphere are sensitive to the vertical distribution of PV at the tropical boundary. With a surface distribution of temperature that decreases linearly with latitude, the jet maximum occurs at the tropical boundary and moves with it. The overall pattern of wind distribution is not sensitive to the change of the position of the tropical boundary. Finally, the temperature and wind distributions of an extratropical troposphere with a finite PV gradient are calculated. It is found that the larger the isentropic PV gradient, the warmer the troposphere and the weaker the jet.

Sun, De-Zheng↗

Thermal/Fluid Analysis of a Composite Heat Exchanger for Use on the RLV Rocket Engine

As part of efforts to design a regeneratively cooled composite nozzle ramp for use on the reusable vehicle (RLV) rocket engine, a C-SiC composite heat exchanger concept was proposed for thermal performance evaluation. To test the feasibility of the concept, sample heat exchanger panels were made to fit the Glenn Research Center's cell 22 for testing. Operation of the heat exchanger was demonstrated in a combustion environment with high heat fluxes similar to the RLV Aerospike Ramp. Test measurements were reviewed and found to be valuable for the on-going fluid and thermal analysis of the actual RLV composite ramp. Since the cooling fluid for the heat exchanger is water while the RLV Ramp cooling fluid is LH2, fluid and therma models were constructed to correlate to the specific test set-up. The knowledge gained from this work will be helpful for analyzing the thermal response of the actual RLV Composite Ramp. The coolant thermal properties for the models are taken from test data. The heat exchanger's cooling performance was analyzed using the Generalized Fluid System Simulation Program (GFSSP). Temperatures of the heat exchanger's structure were predicted in finite element models using Patran and Sinda. Results from the analytical models and the tests show that RSC's heat exchanger satisfied the combustion environments in a series of 16 tests.

Nguyen, Dalton↗

Thermal/Fluid Analysis of a Composite Heat Exchanger for Use on the RLV Rocket Engine

As part of efforts to design a regeneratively cooled composite nozzle ramp for use on the reusable vehicle (RLV) rocket engine, an C-SiC composites heat exchanger concept was proposed for thermal performance evaluation. To test the feasibility of the concept, sample heat exchanger panels were made to fit the Glenn Research Center's cell 22 for testing. Operation of the heat exchanger was demonstrated in a combustion environment with high heat fluxes similar to the RLV Aerospike Ramp. Test measurements were reviewed and found to be valuable for the on going fluid and thermal analysis of the actual RLV composite ramp. Since the cooling fluid for the heat exchanger is water while the RLV Ramp cooling fluid is LH2, fluid and thermal models were constructed to correlate to the specific test set-up. The knowledge gained from this work will be helpful for analyzing the thermal response of the actual RLV Composite Ramp. The coolant thermal properties for the models are taken from test data. The heat exchanger's cooling performance was analyzed using the Generalized Fluid System Simulation Program (GFSSP). Temperatures of the heat exchanger's structure were predicted in finite element models using Patran and Sinda. Results from the analytical models and the tests show that RSC's heat exchanger satisfied the combustion environments in a series of 16 tests.

Nguyen, Dalton↗

Simulation of Impact on a Ductile Polymer Plate

Explicit finite element calculations were used to visualize the deformation and temperature rise in an elastic-plastic plate impacted by a rigid projectile. Results were compared to results of experiments involving ballistic penetration of a "self-healing" thermoplastic. The calculated temperature rise agreed well with the experimental observation, but the total energy absorbed in the penetration event was underestimated in the calculation, which neglected friction.

Cremona, Rebecca L.↗

Thermal-Structural Analysis of PICA Tiles for Solar Tower Test

Thermal protection materials used in spacecraft heatshields are subjected to severe thermal and mechanical loading environments during re-entry into earth atmosphere. In order to investigate the reliability of PICA tiles in the presence of high thermal gradients as well as mechanical loads, the authors designed and conducted solar-tower tests. This paper presents the design and analysis work for this tests series. Coupled non-linear thermal-mechanical finite element analyses was conducted to estimate in-depth temperature distribution and stress contours for various cases. The first set of analyses performed on isolated PICA tile showed that stresses generated during the tests were below the PICA allowable limit and should not lead to any catastrophic failure during the test. The tests results were consistent with analytical predictions. The temperature distribution and magnitude of the measured strains were also consistent with predicted values. The second test series is designed to test the arrayed PICA tiles with various gap-filler materials. A nonlinear contact method is used to model the complex geometry with various tiles. The analyses for these coupons predict the stress contours in PICA and inside gap fillers. Suitable mechanical loads for this architecture will be predicted, which can be applied during the test to exceed the allowable limits and demonstrate failure modes. Thermocouple and strain-gauge data obtained from the solar tower tests will be used for subsequent analyses and validation of FEM models.

Agrawal, Parul↗

JWST ISIM Distortion Analysis Challenge

Very tight distortion requirements are imposed on the JWST's ISM structure due to the sensitivity of the telescope's mirror segment and science instrument positioning. The ISIM structure is a three dimensional truss with asymmetric gusseting and metal fittings. One of the primary challenges for ISIM's analysis team is predicting the thermal distortion of the structure both from the bulk cooldown from ambient to cryo, and the smaller temperature changes within the cryogenic operating environment. As a first cut to estimate thermal distortions, a finite element model of bar elements was created. Elements representing joint areas and metal fittings use effective properties that match the behavior of the stack-up of the composite tube, gusset and adhesive under mechanical and thermal loads. These properties were derived by matching tip deflections of a solid model simplified T-joint. Because of the structure s asymmetric gusseting, this effective property model is a first attempt at predicting rotations that cannot be captured with a smeared CTE approach. In addition to the finite element analysis, several first order calculations have been performed to gauge the feasibility of the material design. Because of the stringent thermal distortion requirements at cryogenic temperatures, a composite tube material with near zero or negative CTE is required. A preliminary hand analysis of the contribution of the various components along the distortion path between FGS and the other instruments, neglecting second order effects were examined. A plot of bounding tube longitudinal and transverse CTEs for thermal stability requirements was generated to help determine the feasibility of meeting these requirements. This analysis is a work in progress en route to a large degree of freedom hi-fidelity FEA model for distortion analysis. Methods of model reduction, such as superelements, are currently being investigated.

Cifie, Emmanuel↗

The computation of thermo-chemical nonequilibrium hypersonic flows

Several conceptual designs for vehicles that would fly in the atmosphere at hypersonic speeds have been developed recently. For the proposed flight conditions the air in the shock layer that envelops the body is at a sufficiently high temperature to cause chemical reaction, vibrational excitation, and ionization. However, these processes occur at finite rates which, when coupled with large convection speeds, cause the gas to be removed from thermo-chemical equilibrium. This non-ideal behavior affects the aerothermal loading on the vehicle and has ramifications in its design. A numerical method to solve the equations that describe these types of flows in 2-D was developed. The state of the gas is represented with seven chemical species, a separate vibrational temperature for each diatomic species, an electron translational temperature, and a mass-average translational-rotational temperature for the heavy particles. The equations for this gas model are solved numerically in a fully coupled fashion using an implicit finite volume time-marching technique. Gauss-Seidel line-relaxation is used to reduce the cost of the solution and flux-dependent differencing is employed to maintain stability. The numerical method was tested against several experiments. The calculated bow shock wave detachment on a sphere and two cones was compared to those measured in ground testing facilities. The computed peak electron number density on a sphere-cone was compared to that measured in a flight test. In each case the results from the numerical method were in excellent agreement with experiment. The technique was used to predict the aerothermal loads on an Aeroassisted Orbital Transfer Vehicle including radiative heating. These results indicate that the current physical model of high temperature air is appropriate and that the numerical algorithm is capable of treating this class of flows.

Candler, Graham↗

Thermal distortion analysis of the space station solar dynamic concentrator

A method was developed to evaluate the thermal distortion of the Space Station Solar Dynamic Concentrator and the effects of thermal distortion on concentrator optical performance. The analytical method includes generating temperature distributions with TRASYS and SINDA models, interfacing the SINDA results with the SINDA-NASTRAN Interface Program (SNIP), calculating thermal distortion with a NASTRAN/PATRAN finite element model, and providing flux distribution maps within the receiver with the ray tracing OFFSET program. Temperature distributions, thermally induced slope errors, and flux distribution maps within the receiver are discussed. Results during a typical orbit indicate that temperatures of the hexagonal panels and triangular facets range between -18 and 99 C (-1 to 210 F), facet rotations are less than 0.2 mrad, and a change in facet radius due to thermal flattening is less than 5 percent. The predicted power loss with thermal distortion effects was less than 0.3 percent. The thermal distortion of the Solar Dynamic concentrator has negligible effect on the flux distribution within the receiver cavity.

Trudell, Jeffrey J.↗

Thermal distortion analysis of the Space Station solar dynamic concentrator

A method was developed to evaluate the thermal distortion of the Space Station Solar Dynamic Concentrator and the effects of thermal distortion on concentrator optical performance. The analytical method includes generating temperature distributions with TRASYS and SINDA models, interfacing the SINDA results with the SINDA-NASTRAN Interface Program (SNIP), calculating thermal distortion with a NASTRAN/PATRAN finite element model, and providing flux distribution maps within the receiver with the ray tracing OFFSET program. Temperature distributions, thermally induced slope errors, and flux distribution maps within the receiver are discussed. Results during a typical orbit indicate that temperatures of the hexagonal panels and triangular facets range between -18 and 99 C (-1 to 210 F), facet rotations are less than 0.2 mrad, and a change in facet radius due to thermal flattening is less than 5 percent. The predicted power loss with thermal distortion effects was less than 0.3 percent. The thermal distortion of the Solar Dynamic concentrator has negligible effect on the flux distribution within the receiver cavity.

Trudell, Jeffery J.↗

The Oxidation Kinetics of Continuous Carbon Fibers in a Cracked Ceramic Matrix Composite

Experimental observations and results suggest two primary regimes as a function of temperature, i.e., diffusion and reaction controlled kinetics. Thermogravimetric analysis of carbon fiber in flowing oxygen gave an activation energy of 64.1 kJ/mol in the temperature range of 500 to 600 C and an apparent activation energy of 7.6 kJ/mol for temperatures from 600 to 1400 C. When C/SiC composite material was unstressed, matrix effects at temperatures from 900 to 1400 C protected the internal fibers. When under stress, self-protection was not observed. Increasing the stress from 10 to 25 ksi caused a 67 to 82 percent reduction in times to failure at temperatures from 750 to 1500 C. Based on experimental results, observation, and theory, a finite difference model was developed, which simulates the diffusion of oxygen into a matrix crack that is bridged by carbon fibers. The model allows the influence of important variables on oxidation kinetics to be studied systematically, i.e., temperature, reaction rate constant, diffusion coefficient, environment, and sample geometry.

Halbig, Michael C.↗

Ascent Heating Thermal Analysis on Spacecraft Adaptor Fairings

When the Crew Exploration Vehicle (CEV) is launched, the spacecraft adaptor (SA) fairings that cover the CEV service module (SM) are exposed to aero heating. Thermal analysis is performed to compute the fairing temperatures and to investigate whether the temperatures are within the material limits for nominal ascent aeroheating case. The ascent heating is analyzed by using computational fluid dynamics (CFD) and engineering codes at Marshall Space Flight Center. The aeroheating environment data used for this work is known as Thermal Environment 3 (TE3) heating data. One of the major concerns is with the SA fairings covering the CEV SM and the SM/crew launch vehicle (CLV) flange interface. The TE3 heating rate is a function of time, wall temperature, and the spatial locations. The implementation of the TE3 heating rate as boundary conditions in the thermal analysis becomes challenging. The ascent heating thermal analysis on SA fairings and SM/CLV flange interface are performed using two commercial software packages: Cullimore & Ring (C&R) Thermal Desktop (TD) 5.1 and MSC Patran 2007r1 b. TD is the pre-and post-processor for SINDA, which is a finite-difference-based solver. In TD, the geometry is built and meshed, the boundary conditions are defined, and then SINDA is used to compute temperatures. MSC Pthermal is a finite-element- based thermal solver. MSC Patran is the pre- and post-processor for Pthermal. Regarding the boundary conditions, the convection, contact resistance, and heat load can be imposed in different ways in both programs. These two software packages are used to build the thermal model for the same analysis to validate each other and show the differences in the modeling details.

Wang, Xiao Yen↗

Numerical Prediction of Transient Axial Thrust and Internal Flows in a Rocket Engine Turbopump

This paper presents the application of the Generalized Fluid System Simulation Program (GFSSP) to model the time-dependent flow in a complex secondary flow circuit of the turbopump of the Fastrac engine currently under development at Marshall Space Flight Center. GFSSP is a general purpose computer program for analyzing steady-state and time-dependant flowrates, pressures, temperatures, and concentrations in a complex flow network. The program employs a finite volume formulation of mass, momentum and energy conservation equations in conjunction with the thermodynamic equation of state of real fluids. GFSSP was used to calculate the axial thrust and internal flow distribution of the Fastrac engine turbopump during the start and shut down transients. The models discussed in this paper use boundary conditions that were extracted from turbopump test data. The GFSSP predicted turbopump secondary flow passage pressures and temperatures were compared with actual measured values.

VanHooser, Katherine↗

Fatigue Characterization of Alloy 10: a 1300F Disk Alloy for Small Gas Turbine Engines

A detailed fatigue characterization of Alloy 10, a high strength nickel-based disk alloy, was conducted on test coupons machined from a 'pancake' forging. Smooth bar, strain controlled fatigue testing at various R-ratios was run at representative bore, 750 F, and rim, 1300 F, temperatures. This was followed by notch fatigue testing (Kt=2) run under load control. Analysis of the fatigue data using a Smith-Watson-Topper approach and finite element analysis of the notch root was employed to understand material behavior in these tests. Smooth bar fatigue data showed a significant R-ratio dependence at either test temperature which could be accounted for using a Smith-Watson-Topper parameter (SWT). In general, fatigue life was longer at 750 F than 1300 F for a given SWT. For notch fatigue tests, life was longer at 750 F than 1300 F but only at higher stresses. This was attributed to differences in alloy strength. At lower stresses, finite element analysis suggested that convergence of fatigue life at both temperatures resulted from relaxation of stresses at the notch root in the 1300 F tests.

Gayda, John↗

Thermal Effects on the Compressive Behavior of IM7/PET15 Laminates

The effect of changing operating temperature on the compressive response of IM7/PETI5 composite laminates is investigated within this paper. The three temperatures evaluated for this study were 129 C, 21 C, and 177 C, a spectrum from cryogenic to an elevated operating temperature. Laminate compressive strength property testing was conducted using the Wyoming Combined Load Compression fixture to generate strength data at the three operating temperatures of interest for several lay-ups. A three-dimensional finite element analysis model of a [90/0]8s composite laminate subject to compressive loading is developed. The model is used to study the key attributes of the laminate that significantly influence the state of stress in the laminate. Both the resin rich layer located between lamina and the thermal residual stresses present in the laminate due to curing are included in the analysis model. For the laminate modeled, the effect of modeling temperature dependent material properties was determined to be insignificant for the operating temperatures studied. Simply using the material properties measured at the operating temperature of interest was sufficient for predicting stresses accurately in a linear analysis for the current problem. The three-dimensional analysis results revealed that the application of an applied compressive axial load in the 0-degree direction decreased the interlaminar stresses present in the laminate initially due to curing. Therefore, failure was concluded not be attributable to the interlaminar stresses in the composite laminate being studied when a compressive load is applied. The magnitude of the measured laminate compressive strength change with a change in temperature is concluded to be dominated by the change in the lamina compressive axial strength with a change in temperature.

Walker, Sandra Polesky↗

3-D Modeling of Directional Solidification of a Non-Dilute Alloy with Temperature and Concentration Fields Coupling via Materials Properties Dependence and via Double Diffusive Convection

Numerical simulation of the HgCdTe growth by the vertical Bridgman method was performed using FIDAP finite element code. Double-diffusive melt convection is analyzed, as the primary factor at controls inhomogeneity of the solidified material. Temperature and concentration fields in the model are also coupled via material properties, such as thermal and solutal expansion coefficients with the dependence on both temperature and concentration, and melting temperature evaluation from pseudobinary CdTe-HgTe phase diagram. Experimental measurements were used to obtain temperature boundary conditions. Parametric study of the melt convection dependence on the gravity conditions was undertaken. It was found, that the maximum convection velocity in the melt can be reduced under certain conditions. Optimal conditions to obtain a near flat solidified interface are discussed. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D calculations are compared with previous 2- D findings. A video film featuring 3-D melt convection will be presented.

Bune, Andris V.↗