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 703 records · Page 39

Three-Dimensional Ignition and Flame Propagation Above Liquid Fuel Pools: Computational Analysis

A three-dimensional unsteady reactive Navier-Stokes code is developed to study the ignition and flame spread above liquid fuels initially below the flashpoint temperature. Opposed air flow to the flame spread due to forced and/or natural convection is considered. Pools of finite width and length are studied in air channels of prescribed height and width. Three-dimensional effects of the flame front near the edge of the pool are captured in the computation. The formation of a recirculation zone in the gas phase similar to that found in two-dimensional calculations is also present in the three-dimensional calculations. Both uniform spread and pulsating spread modes are found in the calculated results.

Cai, Jinsheng↗

Nonlinear Thermal Analyses of a Liquid Hydrogen Tank Wall

A thermal evaluation of a composite tank wall design for a liquid hydrogen tank was performed in the present study. The primary focus of the current effort was to perform one-dimensional, temperature nonlinear, transient thermal analyses to determine the through-the-thickness temperature profiles. These profiles were used to identify critical points within the flight envelope that could have detrimental effects on the adhesive bondlines used in the construction of the tank wall. Additionally, this paper presents the finite element models, analysis strategies, and thermal analysis results that were determined for several vehicle flight conditions. The basic tank wall configuration used to perform the thermal analyses consisted of carbon-epoxy facesheets and a Korex honeycomb core sandwich that was insulated with an Airex cryogenic foam and an Alumina Enhanced Thermal Barrier (AETB-12). Nonlinear, transient thermal analyses were conducted using the ABAQUS finite element code. Tank wall models at a windward side location on the fuel tank were analyzed for three basic flight conditions: cold-soak (ground-hold), ascent, and re-entry. Additionally, three ambient temperature boundary conditions were applied to the tank wall for the cold-soak condition, which simulated the launch pad cooldown process. Time-dependent heating rates were used in the analyses of the ascent and reentry segments of the flight history along with temperature dependent material properties. The steady-state through-the-thickness temperature profile from the cold-soak condition was used as the initial condition for the ascent analyses. Results from the nonlinear thermal analyses demonstrated very good correlation with results from similar models evaluated by Northrop- Grumman using a different analysis tool. Wall through-the-thickness temperature gradients as a function of flight time were obtained for future incorporation into a full-scale thermostructural analysis to evaluate the adhesive bondlines. As a result of the thermal analyses conducted, a sufficient level of confidence was demonstrated in the thermal modeling and analysis capabilities of ABAQUS to warrant future use as a thermo-structural analysis tool to evaluate cryogenic tank wall designs.

Smeltzer, Stanley S., III↗

Why Is the Overheating Problem Difficult: the Role of Entropy

The development of computational fluid dynamics over the last few decades has yielded enormous successes and capabilities being routinely employed today; however there remain some open problems to be properly resolved-some are fundamental in nature and some resolvable by operational changes. These two categories are distinguished and broadly explored previously. One, that belongs to the former, is the so-called overheating problem, especially in rarefying flow. This problem up to date still dogs every method known to the author; a solution to it remains elusive. The study in this paper concludes that: (1) the entropy increase is quantitatively linked to the increase in the temperature increase, (2) it is argued that the overheating is inevitable in the current shock capturing or traditional finite difference framework, and (3) a simple hybrid method is proposed that removes the overheating problem in the rarefying problems, but also retains the property of accurate shock capturing. This remedy (enhancement of current numerical methods) can be included easily in the present Eulerian codes.

Liou, Meng-Sing↗

Application of the Quadrupole Method for Simulation of Passive Thermography

Passive thermography has been shown to be an effective method for in-situ and real time nondestructive evaluation (NDE) to measure damage growth in a composite structure during cyclic loading. The heat generation by subsurface flaw results in a measurable thermal profile at the surface. This paper models the heat generation as a planar subsurface source and calculates the resultant temperature profile at the surface using a three dimensional quadrupole. The results of the model are compared to finite element simulations of the same planar sources and experimental data acquired during cyclic loading of composite specimens.

Winfree, William P.↗

A comprehensive model to determine the effects of temperature and species fluctuations on reaction rates in turbulent reacting flows

A computationally-viable model describing the interaction between fluid-mechanical turbulence and finite-rate combustion reactions, principally in high-speed flows was developed. Chemical kinetic mechanisms, complete and global, were developed describing the finite rate reaction of fuels of interest to NASA. These fuels included principally hydrogen and silane, although a limited amount of work involved hydrocarbon fuels as well.

Chinitz, W.↗

High Temperature Propulsion System Structural Seals for Future Space Launch Vehicles

Durable, flexible sliding seals are required in advanced hypersonic engines to seal the perimeters of movable engine ramps for efficient, safe operation in high heat flux environments at temperatures of 2000 to 2500 F. Current seal designs do not meet the demanding requirements for future engines, so NASA s Glenn Research Center is developing advanced seals and preloading devices to overcome these shortfalls. An advanced ceramic wafer seal design and two types of seal preloading devices were evaluated in a series of compression, scrub, and flow tests. Silicon nitride wafer seals survived 2000 in. 1000 cycles) of scrubbing at 1600 F against an Inconel 625 rub surface with no chips or signs of damage. Flow rates measured for the wafers before and after scrubbing were almost identical and were up to 32 times lower than those recorded for the best braided rope seal flow blockers. Canted coil springs and silicon nitride compression springs showed promise conceptually as potential seal preloading devices to help maintain seal resiliency. A finite element model of the canted coil spring revealed that it should be possible to produce a spring out of high temperature materials for applications at 2000+ F.

Dunlap, Patrick H., Jr.↗

Numerical investigation of the thermal stratification in cryogenic tanks subjected to wall heat flux

The flow pattern and thermal stratification of a cryogenic cylindrical tank are numerically studied. The tank sidewall is subjected to either a uniform heat-flux or two discrete levels of uniform heat-flux at the upper and lower halves of the tank wall. The tank bottom is kept at a constant temperature controlled by the heat exchanger of a thermodynamic vent system. The tank pressure is also assumed constant resulting in a constant saturation temperature at the interface which is higher than the tank bottom temperature. The effects of vapor motion and vapor superheat on the mass and heat transfer processes at the interface are assumed negligible such that the calculations of liquid region can be decoupled from the vapor region. Dimensionless steady-state conservation equations are solved by a finite-difference method. The effects of modified Rayleigh number, Prandtl number, tank aspect ratio, wall heat-flux parameter, and wall heat-flux distribution on the liquid velocity and temperature fields are investigated. Also, their effects on the rate of heat transfer through the interface and the tank bottom are examined.

Lin, Chin-Shun↗

Numerical investigation of the thermal stratification in cryogenic tanks subjected to wall heat flux

The flow pattern and thermal stratification of a cryogenic cylindrical tank are numerically studied. The tank sidewall is subjected to either a uniform heat-flux or two discrete levels of uniform heat-flux at the upper and lower halves of the tank wall. The tank bottom is kept at a constant temperature controlled by the heat exchanger of a thermodynamic vent system. The tank pressure is also assumed constant resulting in a constant saturation temperature at the interface which is higher than the tank bottom temperature. The effects of vapor motion and vapor superheat on the mass and heat transfer processes at the interface are assumed negligible such that the calculations of liquid region can be decoupled from the vapor region. Dimensionless steady-state conservation equations are solved by a finite-difference method. The effects of modified Rayleigh number, Prandtl number, tank aspect ratio, wall heat-flux parameter, and wall heat-flux distribution on the liquid velocity and temperature fields are investigated. Also, their effects on the rate of heat transfer through the interface and the tank bottom are examined.

Lin, Chin-Shun↗

Computation of the wind by forced adjustment to the height field.

Experimental study of the problem of the forced adjustment of the wind field to the height field with the Mintz-Arakawa two-level atmospheric general circulation model. In all but one of the experiments, the height field was assumed to be perfectly observed at 6-hr intervals, over a time period of one day or less, and from these height data the vector wind field was computed by forced dynamical adjustment. In one experiment, the temperature alone was prescribed. The winds computed in these experiments were compared with the 'control' winds of the general circulation simulation. The best agreement between the computed and the control winds was obtained when the time-differencing scheme in the governing finite-difference equations of motion had a large rate of damping of high-frequency motions. This damping rate also determined the optimum fraction and frequency of restoration of the height (or temperature) fields. With strong damping, total restoration at every time step gave the most rapid rate of wind error reduction and the smallest asymptotic limit of the wind error.

Mesinger, F.↗

EZDESIT: A Computer Program for Structural Element Sizing and Vehicle Weight Prediction

A method is presented for formulating stiffness terms and thermal coefficients of fiber reinforced composite, stiffened panels for input to finite element analysis (FEA). The method is robust enough to handle panels with general cross sectional shapes, including those which are unsymmetric and/or unbalanced. New thermal coefficients are introduced to quantify panel response from through-the-thickness temperature gradients. Equations are defined for stiffness, thermal expansion, and thermal bending that consider the full complement of membrane, bending, membrane-bending coupling, and transverse shear behavior. A technique of implementing this capability with a single plane of shell finite elements using the MSC/NASTRAN™ FEA program is revealed. Finally, an example of a composite, hat-stiffened panel is included to demonstrate errors that occur when an unsymmetric panel is symmetrically formulated as traditionally done. These erroneous values and the correct ones produced from the presented method are listed for the panel's stiffness matrices, thermal expansion and bending vectors, and the forces and moments that are thermally induced from in-plane and through-the-thickness temperature gradients.

Jeffrey A Cerro↗

Composite Stiffened Panels: Part I - Stiffness, Thermal Expansion, and Thermal Bending Formulation for Finite Element Analysis

A method is presented for formulating stiffness terms and thermal coefficients of fiber-reinforced composite, stiffened panels for input to finite element analysis (FEA). The method is robust enough to handle panels with general cross sectional shapes, including those which are unsymmetric and/or unbalanced. New thermal coefficients are introduced to quantify panel response from through-the-thickness temperature gradients. Equations are defined for stiffness, thermal expansion, and thermal bending that consider the full complement of membrane, bending, membrane-bending coupling, and transverse shear behavior. A technique of implementing this capability with a single plane of shell finite elements using the MSC/NASTRAN(TradeMark) FEA program is revealed. Finally, an example of a composite, hat-stiffened panel is included to demonstrate errors that occur when an unsymmetric panel is symmetrically formulated as traditionally done. These erroneous values and the correct ones produced from the presented method are listed for the panel's stiffness matrices, thermal expansion and bending vectors, and the forces and moments that are thermally induced from in-plane and through-the-thickness temperature gradients.

Craig S Collier↗

Finite Element Thermal Model for Ultrasonic Welding of Thermoplastic Composites

Ultrasonic welding, UW, is a fast and energy-efficient technique for joining thermoplastic composites. It involves the use of high-frequency mechanical vibrations and a static welding force to melt and join adherends. Ultrasonic welding is an enabling technology to reduce the cost and complexity of in-space construction because lightweight thermoplastic composite components can be packaged compactly for launch and then efficiently assembled using supervised autonomous robotic technologies on site. However, the temperatures in space present challenges to UW, and it is critical that the efficacy of process parameters selected for manufacturing in space is understood prior to launch. To this end, a three-dimensional finite element model is presented in this technical presentation. The model incorporates equations for effects of viscoelastic heating and heat transfer on the welding process. The proposed model is applied to predict the temperature distribution in single lap shear, SLS, samples composed of AS4/PEEK (TC1200) composite that were welded using a terrestrial machine as part of a comprehensive weldability study. Thermocouple and infrared spot sensor data from the SLS samples provide empirical temperature measurements for calibration and validation of the UW thermal model. Calibration and validation of the model is an important step given the significant uncertainties in material properties such as the loss modulus and necessary assumptions in the physics implementations which allow the model to converge in an acceptable amount of time. The validated thermal model can be used to simulate the process for the space environment. Future validation sample testing is planned in a vacuum chamber. The result will be a model capable of guiding process parameter selection to ensure acceptable weld bonds when manufacturing in space.

Josh Fody↗

Modeling Orifice Pulse Tube Coolers

We have developed a calculational model that treats all the components of an orifice pulse tube cooler. We base our analysis on 1-dimensional thermodynamic equations for the regenerator and we assume that all mass flows, pressure oscillations and temperature oscillations are small and sinusoidal. Non-linear pressure drop effects are included in the regenerator to account for finite pressure amplitude effects. The resulting mass flows and pressures are matched at the boundaries with the other components of the cooler: compressor, aftercooler, cold heat exchanger, pulse tube, hot heat exchanger, orifice and reservoir. The results of the calculation are oscillating pressures, mass flows and enthalpy flows in the main components of the cooler. By comparing with the calculations of other available models, we show that our model is very similar to REGEN 3 from NIST and DeltaE from Los Alamos National Lab. Our model is much easier to use than other available models because of its simple graphical interface and the fact that no guesses are required for the operating pressures or mass flows. In addition, the model only requires a few minutes of running time allowing many parameters to be optimized in a reasonable time. A version of the model is available for use over the World Wide Web at http://irtek.arc.nasa.gov. Future enhancements include adding a bypass orifice and including second order terms in steady mass streaming and steady heat transfer. A two-dimensional anelastic approximation of the fluid equations will be used as the basis for the latter analysis. Preliminary results are given in dimensionless numbers appropriate for oscillating compressible flows. The model shows how transverse heat transfer reduces enthalpy flow, particularly for small pulse tubes. The model also clearly shows mass recirculation in the open tube on the order of the tube length. They result from the higher order Reynolds stresses. An interesting result of the linearized approach is that the steady mass streaming does not affect the enthalpy flow at second order. The major effect of recirculating mass streaming is to increase transverse temperature gradients, which leads to higher entropy production and reduced efficiency.

Roach, Kittel P.↗

Analytical simulation of weld effects in creep range

The inelastic analysis procedure used to investigate the effect of welding on the creep rupture strength of a typical Liquid Metal Fast Breeder Reactor (LMFBR) nozzle is discussed. The current study is part of an overall experimental and analytical investigation to verify the inelastic analysis procedure now being used to design LMFBR structural components operating at elevated temperatures. Two important weld effects included in the numerical analysis are: (1) the residual stress introduced in the fabrication process; and (2) the time-independent and the time-dependent material property variations. Finite element inelastic analysis was performed on a CRAY-1S computer using the ABAQUS program with the constitutive equations developed for the design of LMFBR structural components. The predicted peak weld residual stresses relax by as much as 40% during elevated temperature operation, and their effect on creep-rupture cracking of the nozzle is considered of secondary importance.

Dhalla, A. K.↗

Contact thermal shock test of ceramics

A novel quantitative thermal shock test of ceramics is described. The technique employs contact between a metal-cooling rod and hot disk-shaped specimen. In contrast with traditional techniques, the well-defined thermal boundary condition allows for accurate analyses of heat transfer, stress, and fracture. Uniform equibiaxial tensile stresses are induced in the center of the test specimen. Transient specimen temperature and acoustic emission are monitored continuously during the thermal stress cycle. The technique is demonstrated with soda-lime glass specimens. Experimental results are compared with theoretical predictions based on a finite-element method thermal stress analysis combined with a statistical model of fracture. Material strength parameters are determined using concentric ring flexure tests. Good agreement is found between experimental results and theoretical predictions of failure probability as a function of time and initial specimen temperature.

Rogers, W. P.↗

Near-surface temperature gradients and their effects on thermal-infrared emission spectra of particulate planetary surfaces

The infrared energy emitted from a planetary surface is generated within a finite depth determined by the material's absorption skin depth. This parameter varies significantly with wavelength in the infrared but has an average value of around 50 microns for most geologic materials. In solid rock, heat transfer is efficient enough so that this 50 micron zone of the near surface from which the radiation emanates will be more or less isothermal. In particulate materials, however, heat transfer is more complicated and occurs via a combination of mechanisms, including solid conduction within grains and across grain contacts, conduction through the interstitial gas, and thermal radiation within individual particles and across the void spaces in between grains. On planets with substantial atmospheres, the gas component dominates the heat transfer and tends to mitigate near-surface thermal gradients. However, on airless bodies, the gas component is absent and heat transfer occurs via solid conductions and radiation. If the particles are small relative to the average absorption skin depth, then the top 50-100 microns or so of the surface will be cooled by radiation to space allowing the creation of significant near-surface thermal gradients. In those regions of the spectrum where the absorption coefficient is low, the emission will come from the deeper, warmer parts of the medium, whereas in regions of high absorption, the emission will emanate from shallower, cooler parts of the medium. The resulting emission spectrum will show non-compositional features as a result of the thermal structure in the material. We have modeled the heat transfer in a particulate medium in order to determine the magnitude of near-surface thermal gradients for surfaces on airless bodies and on Mars. We use the calculated thermal structure to determine the effects it has on the infrared emission spectrum of the surface.

Henderson, B. G.↗

Comparison of three-dimensional nonequilibrium PNS codes

A comparison study has been conducted using four recently developed parabolized Navier-Stokes (PNS) codes which have the capability of predicting finite-rate, chemically reacting flows over three-dimensional bodies. These are the (1) UPS code, (2) the STUFF code, (3) the TONIC code, and (4) the VRA-PNS code. All of the codes use the same seven-species, single-temperature air chemistry model, but otherwise they are unique, with different capabilities and characteristics. The differences include upwinding vs central differencing, strongly-coupled vs weakly-coupled chemistry, shock capturing vs shock fitting, finite volume vs finite difference, and full PNS vs thin-layer PNS equations. Three test cases were utilized to compare the codes. The comparisons presented indicate a good agreement among the codes tested.

Buelow, Philip E.↗

Finite Element Method for Thermal Analysis

A two- and three-dimensional, finite-element thermal-analysis program which handles conduction with internal heat generation, convection, radiation, specified flux, and specified temperature boundary conditions is presented. Elements used in the program are the triangle and tetrahedron for two- and three-dimensional analysis, respectively. The theory used in the program is developed, and several sample problems demonstrating the capability and reliability of the program are presented. A guide to using the program, description of the input cards, and program listing are included.

Heuser, J.↗