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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Packaging Technology for SiC High Temperature Circuits Operable up to 500 Degrees Centigrade

New high temperature low power 8-pin packages have been fabricated using commercial fabrication service. These packages are made of aluminum nitride and 96 percent alumina with Au metallization. The new design of these packages provides the chips inside with EM shielding. Wirebond geometry control has been achieved for precise mechanical tests. Au wirebond samples with 45 degree heel-angle have been tested using wireloop test module. The geometry control improves the consistency of measurement of the wireloop breaking point.Also reported on is a parametric study of the thermomechanical reliability of a Au thick-film based SiC die-attach assembly using nonlinear finite element analysis (FEA) was conducted to optimize the die-attach thermo-mechanical performance for operation at temperatures from room temperature to 500 degrees Centigrade. This parametric study centered on material selection, structure design and process control.

Chen, Lian-Yu↗
NASA NTRS1981

Thermal iron ions in high speed solar wind streams. II - Temperatures and bulk velocities

Mitchel and Roelof (1980) reported the detection of iron in high speed solar wind flows using the small, but finite sensitivity of solid state detectors to Fe ions in the low energy (50-200 keV protons) L1 channel of the NOAA/JHU energetic particle experiment (EPE). In the current investigation, the EPE response is modeled to a convected Maxwellian to obtain the thermal velocity, flow angle, and bulk velocity of the iron distribution. It is assumed that the iron bulk flow velocity can be represented as a vector sum of the hydrogen bulk velocity and an interplanetary magnetic field (IMF) aligned velocity increment. It is found that the velocity increment is smaller than the local Alfven speed in magnitude, and that the iron thermal velocity is comparable with or greater than the proton thermal velocity, with the 'thermal' velocity defined as the square root of 2kT/m.

Mitchell, D. G.↗
NASA NTRS1988

Analytical study of He II flow characteristics in the SHOOT transfer line

Pressure and temperature profiles for He II flow in the transfer line for the Superfluid Helium On-Orbit Transfer (SHOOT) experiment are determined using a finite difference method to solve the one-dimensional two-fluid conservation equations. For practical purposes, forced flow of He II with Re greater than 10 to the 6th can be treated as a classical turbulent flow For a maximum SHOOT transfer rate of 800 cu dm/h, the temperature rise between the pump outlet and the receiver tank is about 80 mk with a 7-W heat load imposed on the 4-m-long transfer line. A pump outlet pressure of 160 torr is required to sustain the 800 cu dm/h flow without causing vapor formation within the transfer line.

Lee, J. H.↗
NASA NTRS1978

Evaluation of a hybrid, anisotropic, multilayered, quadrilateral finite element

A multilayered finite element with bending-extensional coupling is evaluated for: (1) buckling of general laminated plates; (2) thermal stresses of laminated plates cured at elevated temperatures; (3) displacements of a bimetallic beam; and (4) displacement and stresses of a single-cell box beam with warped cover panels. Also, displacements and stresses for flat and spherical orthotropic and anisotropic segments are compared with results from higher order plate and shell finite-element analyses.

Robinson, J. C.↗
NASA NTRS1995

An evaluation of a coupled microstructural approach for the analysis of functionally graded composites via the finite-element method

A comparison is presented between the predictions of the finite-element analysis and a recently developed higher-order theory for functionally graded materials subjected to a thorough-thickness temperature gradient. In contrast to existing micromechanical theories that utilize classical (i.e., uncoupled) homogenization schemes to calculate micro-level and macro-level stress and displacement fields in materials with uniform or nonuniform fiber spacing (i.e., functionally graded materials), the new theory explicitly couples the microstructural details with the macrostructure of the composite. Previous thermo-elastic analysis has demonstrated that such coupling is necessary when: the temperature gradient is large with respect to the dimension of the reinforcement; the characteristic dimension of the reinforcement is large relative to the global dimensions of the composite and the number of reinforcing fibers or inclusions is small. In these circumstances, the standard micromechanical analyses based on the concept of the representative volume element used to determine average composite properties produce questionable results. The comparison between the predictions of the finite-element method and the higher-order theory presented herein establish the theory's accuracy in predicting thermal and stress fields within composites with a finite number of fibers in the thickness direction subjected to a thorough-thickness thermal gradient.

Pindera, Marek-Jerzy↗
NASA NTRS1993

An analytical/numerical correlation study of the multiple concentric cylinder model for the thermoplastic response of metal matrix composites

The utility of a recently developed analytical micromechanics model for the response of metal matrix composites under thermal loading is illustrated by comparison with the results generated using the finite-element approach. The model is based on the concentric cylinder assemblage consisting of an arbitrary number of elastic or elastoplastic sublayers with isotropic or orthotropic, temperature-dependent properties. The elastoplastic boundary-value problem of an arbitrarily layered concentric cylinder is solved using the local/global stiffness matrix formulation (originally developed for elastic layered media) and Mendelson's iterative technique of successive elastic solutions. These features of the model facilitate efficient investigation of the effects of various microstructural details, such as functionally graded architectures of interfacial layers, on the evolution of residual stresses during cool down. The available closed-form expressions for the field variables can readily be incorporated into an optimization algorithm in order to efficiently identify optimal configurations of graded interfaces for given applications. Comparison of residual stress distributions after cool down generated using finite-element analysis and the present micromechanics model for four composite systems with substantially different temperature-dependent elastic, plastic, and thermal properties illustrates the efficacy of the developed analytical scheme.

Pindera, Marek-Jerzy↗
NASA NTRS1989

Finite element flow-thermal-structural analysis of aerodynamically heated leading edges

Hypersonic vehicles operate in a hostile aerothermal environment which has a significant impact on their aerothermostructural performance. Significant coupling occurs between the aerodynamic flow field, structural heat transfer, and structural response creating a multidisciplinary interaction. A long term goal of the Aerothermal Loads Branch at the NASA Langley Research Center is to develop a computational capability for integrated fluid, thermal and structural analysis of aerodynamically heated structures. The integrated analysis capability includes the coupling between the fluid and the structure which occurs primarily through the thermal response of the structure, because: (1) the surface temperature affects the external flow by changing the amount of energy absorbed by the structure, and (2) the temperature gradients in the structure result in structural deformations which alter the flow field and attendant surface pressures and heating rates. In the integrated analysis, a finite element method is used to solve: (1) the Navier-Stokes equations for the flow solution, (2) the energy equation of the structure for the temperature response, and (3) the equilibrium equations of the structure for the structural deformation and stresses.

Dechaumphai, Pramote↗
NASA NTRS2003

Integrated Propulsion/Vehicle System Structurally Optimized

Ongoing research and testing are essential in the development of air-breathing hypersonic propulsion technology, and this year some positive advancement was made at the NASA Glenn Research Center. Recent work performed for GTX, a rocket-based combined-cycle, single-stage-to-orbit concept, included structural assessments of both the engine and flight vehicle. In the development of air-breathing engine technology, it is impractical to design and optimize components apart from the fully integrated system because tradeoffs must be made between performance and structural capability. Efforts were made to control the flight trajectory, for example, to minimize the aerodynamic heating effects. Structural optimization was applied to evaluate concept feasibility and was instrumental in the determination of the gross liftoff weight of the integrated system. Achieving low Earth orbit with even a small payload requires an aggressive approach to weight minimization through the use of lightweight, oxidation-resistant composite materials. Assessing the integrated system involved investigating the flight trajectory to determine where the critical load events occur in flight and then generating the corresponding environment at each of these events. Structural evaluation requires the mapping of the critical flight loads to finite element models, including the combined effects of aerodynamic, inertial, combustion, and other loads. NASA s APAS code was used to generate aerodynamic pressure and temperature profiles at each critical event. The radiation equilibrium surface temperatures from APAS were used to predict temperatures through the thickness. Heat transfer solutions using NASA's MINIVER code and the SINDA code (Cullimore & Ring Technologies, Littleton, CO) were calculated at selective points external to the integrated vehicle system and then extrapolated over the entire exposed surface. FORTRAN codes were written to expedite the finite element mapping of the aerodynamic heating effects for the internal structure.

Hunter, James E.↗
NASA NTRSDate not supplied

Analysis of Heat Transfer from Local Heating and Cooling Sources at Cryogenic Temperatures

Understanding the dispersion of heat around a cryogenic fluid tank, specifically the interaction between the cryogenic fluid and the tank wall is critical in the analysis of long duration cryogen storage in microgravity. The heat transfer interaction between a cryogenic storage tank and heat sources from external spacecraft structures is also one of the many factors that determine how much heat enters a tank. Recent flight experiments with two-phase fluids have indicated that local concentrations of heat input (also known as “hot spots”) can cause unwanted affects including local boiling. Computational fluid dynamic (CFD) models can provide a detailed assessment of the heat transfer occurring across a cryogenic storage system. However, CFD modeling takes time to construct and run. A simpler approach that can act as initial guidance for later CFD modeling analyzes external “hot spots” as point or finite heat sources. A radial, finite element network or a local direct solution can effectively estimate the heat spread across a cryogenic storage tank by calculating the temperature and heat load as a function of distance from the heat source. This calculation accounts for the convective heat transfer between the cryogenic fluid and storage tank surface. Similar approaches can be used to determine the effectiveness of cooling from a cryocooler as a finite, local heat sink. This approach allows for quick approximations of the thermal map across a cryogenic tank as well as sensitivity analysis under a wide range of design parameters including gravitational fields as implied through natural convection coefficients.

Cryogenic Fluid Management↗
NASA NTRS1972

Neutral-particle wake method for measuring the atmospheric temperature from a satellite.

Description of a method that would permit a satellite-borne neutral mass spectrometer to measure the atmospheric temperature. The spectrometer examines the partial pressure variations that occur as the wake of a small rectangular baffle is swept across the entrance orifice of the spectrometer. For a given baffle size and for a mounting distance from the orifice, the depth of the resulting pressure minimum depends only on the thermal velocity or temperature of the observed species. The validity of the method can be checked by measuring the wake characteristics of more than one species and/or by employing each of several baffle sizes. The theory includes the effect of a finite orifice size, finite baffle length, and the backscattering of particles from the baffle into the orifice. It is found that a suitable baffle arrangement can be achieved that will permit the temperature to be measured over at least the range normally encountered in the thermosphere (200 to 2000 K) and, depending on the sensitivity and background pressure of the spectrometer, over an altitude range of about 140 to 600 km.

Brace, L. H.↗
NASA NTRS2023

Polar Firn Properties in Greenland and Antarctica and Related Effects on Microwave Brightness Temperatures

In studying the mass balance of polar ice sheets, fluctuations of firn density near the surface is a major uncertainty. In this paper, we explore these variations at locations in the Greenland Ice Sheet and at the Dome C location in Antarctica. Borehole in situ measurements, Snow Radar echoes, microwave brightness temperatures, and modelling results from the Community firn model (CFM) are used. It is shown that firn density profiles can be represented using 3 processes: “long” and “short” length scale density variations and “refrozen layers”. Consistency with this description is observed in the dynamic range of airborne 0.5-2 GHz brightness temperatures and Snow Radar echo peaks in measurements performed in Greenland in 2017. Based on these insights, a new analytical partially coherent model is implemented to explain the microwave brightness temperatures using the three scale description of the firn. Short and long scale firn processes are modelled as a 3D continuous random medium with finite vertical and horizontal correlation lengths as opposed to past 1D random layered medium descriptions. Refrozen layers are described as deterministic sheets with planar interfaces, with the number of refrozen layer interfaces determined by radar observations. Firn density and correlation length parameters used in forward modelling to match measured 0.5-2 GHz brightness temperatures in Greenland show consistency with similar parameters in CFM predictions. Model predictions also are in good agreement with multi-angle 1.4 GHz vertically and horizontally polarized brightness temperature measured by the SMOS satellite at DOME C, Antarctica. This work shows that co-located active and passive microwave measurements can be used to infer polar firn properties that can be compared with predictions of the CFM. In particular, 0.5-2 GHz brightness temperature measurements are shown to be sensitive to long scale firn density fluctuations with density standard deviations in the range 0.01-0.06 g/cm^3 and vertical correlation lengths of 6-20 cm.

Haokui Xu↗
NASA NTRS1980

Application of finite-element-techniques to the interaction of conduction and radiation in an absorbing, scattering and emitting medium

In this paper, the authors demonstrate that a Galerkin finite element method of analysis, utilizing isoparametric elements, offers a viable means of solving continuum thermal radiation problems with conduction in a participating medium. The participating medium was considered to be a gray radiation medium exhibiting isotropic absorption, emission, and scattering characteristics and optical properties that are independent of temperature. The medium was considered to be bounded by infinite parallel opaque, gray surfaces with diffuse emission and reflection characteristics. In solving this problem, a finite element formulation was developed to describe a system in radiative equilibrium. Then the results of this first analysis were linked with a second finite element model which incorporated conduction into the analysis. The results of this study were found to be in good agreement with existing published data. The model offers the following advantageous features: geometric generality, a computational algorithm which is 'convenient' and 'computable', and a functional basis for extension of the radiation model to higher order approximation.

Wu, S. T.↗
NASA NTRS1994

Rotational Relaxation in Nonequilibrium Freejet Expansions of Heated Nitrogen

Rotational temperatures have been measured in rarefied, nonequilibrium, heated freejet expansions of nitrogen using the electron beam fluorescence technique at the University of California at Berkeley Low Density Wind Tunnel facility. Spectroscopic measurements of the (0,0) band of the first negative system of nitrogen reveal the nonequilibrium behavior in the flowfield upstream of, and through the Mach disk, which forms as the freejet expands into a region of finite back pressure. Results compare well with previous freejet expansion data and computations regarding location of the Mach disk and terminal rotational temperature in the expansion. Measurements are also presented for shock thickness based on the rotational temperature changes in the flow. Thickening shock layers, departures of rotational temperature from equilibrium in the expansion region, and downstream rotational temperature recovery much below that of an isentropic normal shock provide indications of the rarefied, nonequilibrium flow behavior. The data are analyzed to infer constant values of the rotational-relaxation collision number from 2.2 to 6.5 for the various flow conditions. Collision numbers are also calculated in a consistent manner for data from other investigations for which is seen a qualitative increase with increasing temperature. Rotational-relaxation collision numbers are seen as not fully descriptive of the rarefied freejet flows. This may be due to the high degree of nonequilibrium in the flowfields, and/or to the use of a temperature-insensitive rotational-relaxation collision number model in the data analyses.

Gochberg, Lawrence A.↗
NASA NTRS1995

Nozzle flow with vibrational nonequilibrium

This research concerns the modeling and numerical solutions of the coupled system of compressible Navier-Stokes equations in cylindrical coordinates under conditions of equilibrium and nonequilibrium thermodynamics. The problem considered was the modeling of a high temperature diatomic gas N2 flowing through a converging-diverging high expansion nozzle. The problem was modeled in two ways. The first model uses a single temperature with variable specific heats as functions of this temperature. For the second model we assume that the various degrees of freedom all have a Boltzmann distribution and that there is a continuous redistribution of energy among the various degrees of freedom as the gas passes through the nozzle. Each degree of freedom is assumed to have its own temperature and, consequently, each system state can be characterized by these temperatures. This suggests that formulation of a second model with a vibrational degree of freedom along with a rotational-translation degree of freedom, each degree of freedom having its own temperature. Initially the vibrational degree of freedom is excited by heating the gas to a high temperature. As the high temperature gas passes through the nozzle throat there is a sudden drop in temperature along with a relaxation time for the vibrational degree of freedom to achieve equilibrium with the rotational-translation degree of freedom. That is, we assume that the temperature change upon passing through the throat is so great that the changes in the vibrational degree of freedom occur at a much slower pace and consequently lags behind the rotational-translational energy changes. This lag results in a finite relaxation time. In this context the term nonequilibrium is used to denote the fact that the energy content of the various degrees of freedom are characterized by two temperatures. We neglect any chemical reactions which could also add nonequilibrium effects. We develop the energy equations for the nonequilibrium model from first principles. The resulting equations, which model the nozzle flow, can be expressed in various forms. In most forms the resulting equations are coupled systems of nonlinear partial differential equations subject to certain boundary conditions. To solve the resulting coupled system of nonlinear partial differential equations, several numerical techniques were investigated: (1) the explicit MacCormack method, (2) the explicit-implicit MacCormack method, (3) the method of operator splitting, (4) factorization schemes, and (5) the Steger-Warming scheme.

Heinbockel, J. H.↗
NASA NTRS1992

Finite difference solution for transient cooling of a radiating-conducting semitransparent layer

Transient solutions were obtained for cooling a semitransparent material by radiation and conduction. The layer is in a vacuum environment so the only means for heat dissipation is by radiation from within the medium leaving through the boundaries. Heat conduction serves only to partially equalize temperatures across the layer. As the optical thickness is increased, steep temperature gradients exist near the boundaries when conduction is relatively small. A solution procedure is required that will provide accurate temperature distributions adjacent to the boundaries, or radiative heat losses will be in error. The approach utilized numerical Gaussian integration to obtain the local radiative source term, and a finite difference procedure with variable space and time increments to solve the transient energy equation.

Siegel, Robert↗
NASA NTRS1991

Optimal design of solidification processes

An optimal design algorithm is presented for the analysis of general solidification processes, and is demonstrated for the growth of GaAs crystals in a Bridgman furnace. The system is optimal in the sense that the prespecified temperature distribution in the solidifying materials is obtained to maximize product quality. The optimization uses traditional numerical programming techniques which require the evaluation of cost and constraint functions and their sensitivities. The finite element method is incorporated to analyze the crystal solidification problem, evaluate the cost and constraint functions, and compute the sensitivities. These techniques are demonstrated in the crystal growth application by determining an optimal furnace wall temperature distribution to obtain the desired temperature profile in the crystal, and hence to maximize the crystal's quality. Several numerical optimization algorithms are studied to determine the proper convergence criteria, effective 1-D search strategies, appropriate forms of the cost and constraint functions, etc. In particular, we incorporate the conjugate gradient and quasi-Newton methods for unconstrained problems. The efficiency and effectiveness of each algorithm is presented in the example problem.

Dantzig, Jonathan A.↗
NASA NTRS1992

A numerical model including PID control of a multizone crystal growth furnace

This paper presents a 2D axisymmetric combined conduction and radiation model of a multizone crystal growth furnace. The model is based on a programmable multizone furnace (PMZF) designed and built at NASA Lewis Research Center for growing high quality semiconductor crystals. A novel feature of this model is a control algorithm which automatically adjusts the power in any number of independently controlled heaters to establish the desired crystal temperatures in the furnace model. The control algorithm eliminates the need for numerous trial and error runs previously required to obtain the same results. The finite element code, FIDAP, used to develop the furnace model, was modified to directly incorporate the control algorithm. This algorithm, which presently uses PID control, and the associated heat transfer model are briefly discussed. Together, they have been used to predict the heater power distributions for a variety of furnace configurations and desired temperature profiles. Examples are included to demonstrate the effectiveness of the PID controlled model in establishing isothermal, Bridgman, and other complicated temperature profies in the sample. Finally, an example is given to show how the algorithm can be used to change the desired profile with time according to a prescribed temperature-time evolution.

Panzarella, Charles H.↗
NASA NTRS1986

Thermal analysis of an orthotropic engineering body

This work is an effort in determining the temperature distribution of an orthotropic engineering body. First the governing partial differential equation is derived, then, a finite element procedure using variational approach is developed for the solution of the governing equation. Finally, a computer program provides a set of satisfactory results.

Fu, K.-C.↗