Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat Conduction”

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 1,027 records · Page 57

Numerical studies of the thermal design sensitivity calculation for a reaction-diffusion system with discontinuous derivatives

The aim of this study is to find a reliable numerical algorithm to calculate thermal design sensitivities of a transient problem with discontinuous derivatives. The thermal system of interest is a transient heat conduction problem related to the curing process of a composite laminate. A logical function which can smoothly approximate the discontinuity is introduced to modify the system equation. Two commonly used methods, the adjoint variable method and the direct differentiation method, are then applied to find the design derivatives of the modified system. The comparisons of numerical results obtained by these two methods demonstrate that the direct differentiation method is a better choice to be used in calculating thermal design sensitivity.

Hou, Jean W.↗

The loss and depth of CO2 ice in comet nuclei

An analytical model was developed to simulate the material differentiation of a cometary nucleus composed of water ice, putative unclathrated CO2 ice and silicate dust in specified proportions. Selective sublimation of any free CO2 ice present in a comet would produce a surface layer of water ice and dust overlying the original CO2 rich material. This surface layer reduces the temperature of buried CO2 ice and restricts the outflow of gaseous CO2. On each orbit, water sublimation at smaller heliocentric distances temporarily reduces the thickness of the water ice and dust layer and liberates dust. The model includes the effects of nucleus rotation, arbitrary orientation of the rotation axis, latitude, heat conduction into the deep interior of the nucleus and restriction of CO2 gas outflow by the water ice and dust layer. The effects of the permeability of the surface water ice layer, the nucleus rotation rate, and the latitude were investigated. Comparison of these and similar results with observations could yield information regarding the permeability and chemical composition of cometary material and suggest sampling strategies to minimize fractionation effects.

Fanale, F. P.↗

Numerical solutions of Navier-Stokes equations for a Butler wing

The flow field is simulated on the surface of a Butler wing in a uniform stream. Results are presented for the Mach number 3.5 and Reynolds number of 2,000,000. The simulation is done by integrating the viscous Navier-Stokes equations. These equations govern the unsteady, viscous, compressible and heat conducting flow of an ideal gas. The equations are written in curvilinear coordinates so that the wing surface is represented accurately. O-type and H-type grids have been used for this study, and results are compared. The governing equations are solved by the McCormack time-split method, and the results are compared with other theoretical and experimental results. The codes are written in FORTRAN, vectorized and currently run on the CDC Vector Processing System (VPS-32) computer.

Abolhassani, Jamshid S.↗

Brightness curve and porosity of cometary nuclei

The brightness curve of a comet is analyzed in terms of porosity and nature of volatiles in its nucleus. Mobility of pores, densification, phase changes, differential thermal expansion and heat flow into and out of the nucleus are discussed leading to a method for estimating the effective heat conductivity and porosity of the nucleus. The method is applied to comets P/Giacobini-Zinner, P/Halley, P/Kopff, P/Encke, P/Tempel 2 and P/d'Arrest.

Smoluchowski, R.↗

Thermal shaft effects on load-carrying capacity of a fully coupled, variable-properties cryogenic journal bearing

The purpose of this work was to perform a rather complete analysis for a cryogenic (oxygen) journal bearing. The Reynolds equation required coupling and simultaneous solution with the fluid energy equation. To correctly account for the changes in the fluid viscosity, the fluid energy equation was coupled with the shaft and bearing heat conduction energy equations. The effects of pressure and temperature on the density, viscosity, and load-carrying capacity were further discussed as analysis parameters, with respect to relative eccentricity and the angular velocity. The isothermal fluid case and the adiabatic fluid case represented the limiting boundaries. The discussion was further extrapolated to study the Sommerfeld number dependency on the fluid Nusselt number and its consequence on possible total loss of load-carrying capacity and/or seizure (catastrophic failure).

Braun, M. J.↗

Numerical solutions of Navier-Stokes equations for a Butler wing

The flow field is simulated on the surface of a Butler wing in a uniform stream. Results are presented for Mach number 3.5 and Reynolds number of 2,000,000. The simulation is done by integrating the viscous Navier-Stokes equations. These equations govern the unsteady, viscous, compressible and heat conducting flow of an ideal gas. The equations are written in curvilinear coordinates so that the wing surface is represented accurately. The O-type and H-type grids have been used for this study, and results are compared. The governing equations are solved by the MacCormack time-split method, and the results are compared with other theoretical and experimental results. The codes are written in FORTRAN, vectorized and currently run on the CDC Vector Processing System (VPS-32) computer.

Abolhassani, Jamshid S.↗

Multiply scaled constrained nonlinear equation solvers

To improve the numerical stability of nonlinear equation solvers, a partitioned multiply scaled constraint scheme is developed. This scheme enables hierarchical levels of control for nonlinear equation solvers. To complement the procedure, partitioned convergence checks are established along with self-adaptive partitioning schemes. Overall, such procedures greatly enhance the numerical stability of the original solvers. To demonstrate and motivate the development of the scheme, the problem of nonlinear heat conduction is considered. In this context the main emphasis is given to successive substitution-type schemes. To verify the improved numerical characteristics associated with partitioned multiply scaled solvers, results are presented for several benchmark examples.

Padovan, Joe↗

Electrically conductive, black thermal control coatings for space craft application. II - Silicone matrix formulation

Five black electrically conductive thermal-control coatings have been formulated and tested for application on the Galileo spacecraft. The coatings consisted of organic and inorganic systems applied on titanium and aluminum surfaces. The coatings were tested under simulated space environment conditions. Coated specimens were subjected to thermal radiation and convective and conductive heating from -196 to 538 C. Mechanical, physical, thermal, electrical, and optical characteristics, formulation, mixing, application, surface preparation of substrates, and a method of determining electrical resistance are presented for the silicone matrix formulation designated as GF-580.

Hribar, V. F.↗

The normal modes of the thermosphere

The linearized momentum, energy, and continuity equations for the thermosphere can be reduced to a form that gives the vertical structure for each horizontal wave mode. The vertical structure equation can be described in terms of the normal modes, or eigenmodes, of the thermosphere. The latter are obtained by using a 27-layer model that includes a realistic temperature profile and the effects of the Lorentz force, viscosity, and heat conduction. The normal modes have one real eigenfrequency for every two complex conjugate eigenfrequency values. The real modes have a dominant rotational wind component and are nonpropagating. The complex modes have comparable divergent and rotational wind components. The complex eigenvalues give vertically propagating modes, primarily associated with the transient response to forcing, and are significantly affected by the dissipation in the upper E region and F region. Results show that the rotational wind component dominates in the steady state when the forcing is due to the two-cell convection pattern at high latitudes and that the normal modes explain the large shears and large winds speeds that are typically observed in the high-latitude E region. The vertical energy flux for the normal modes is also calculated. The results show that the flux is upward above 130 km but downward in the lower E region for the total solution. The downward energy flux is a contribution from the real eigenmode structure.

Larsen, M. F.↗

The evolution of cooling flows. I - Self-similar cluster flows

The evolution of a cooling flow from an initial state of hydrostatic equilibrium in a cluster of galaxies is investigated. After gas mass and energy are injected into the cluster at an early phase, the gas approaches hydrostatic equilibrium over most of the cluster and cooling becomes important in the dense central regions. As time passes, cooling strongly affects an increasing amount of gas. The effects of mass removal from the flow, the inclusion of magnetic or cosmic-ray pressure, and heat conduction are considered individually.

Chevalier, Roger A.↗

Large-scale volcano-ground ice interactions on Mars

In the present investigation of Martian volcano-ground ice interaction processes, a numerical model is developed that encompasses conductive heat transport, surface radiation, heat transfer to the atmosphere, and H2O phase-changes in an ice-rich permafrost over which lave flows erupt, followed by the intrusion of sills. An examination is made of the two large scale interaction regions formed near Aeolis Mensae and near the volcano Hadriaca Patera, northeast of Hellas. The inferred channel discharges are compared to discharge rates calculated for lava-ground ice interactions, showing that meltwater (probably accumulated under the surface) was rapidly released at a discharge rate that was limited by soil permeability. The volcano-ground ice interactions have been an important Martian geologic process, and could account for the palagonites constituting the Martian dust.

Squyres, Steven W.↗

An approximation theory for the identification of nonlinear distributed parameter systems

An abstract approximation framework for the identification of nonlinear distributed parameter systems is developed. Inverse problems for nonlinear systems governed by strongly maximal monotone operators (satisfying a mild continuous dependence condition with respect to the unknown parameters to be identified) are treated. Convergence of Galerkin approximations and the corresponding solutions of finite dimensional approximating identification problems to a solution of the original finite dimensional identification problem is demonstrated using the theory of nonlinear evolution systems and a nonlinear analog of the Trotter-Kato approximation result for semigroups of bounded linear operators. The nonlinear theory developed here is shown to subsume an existing linear theory as a special case. It is also shown to be applicable to a broad class of nonlinear elliptic operators and the corresponding nonlinear parabolic partial differential equations to which they lead. An application of the theory to a quasilinear model for heat conduction or mass transfer is discussed.

Banks, H. T.↗

Hermetically sealable package for hybrid solid-state electronic devices and the like

A light-weight, inexpensively fabricated, hermetically sealable, repairable package for small electronic or electromechanical units, having multiple connections, is described. A molded ring frame of polyamide-imide plastic (Torlon) is attached along one edge to a base plate formed of a highly heat conducting material, such as aluminum or copper. Bores are placed through a base plate within the area of the edge surface of ring frame which result in an attachment of the ring frame to the base plate during molding. Electrical leads are molded into the ring frame. The leads are L-shaped gold-plated copper wires imbedded within widened portions of the side wall of the ring frame. Within the plastic ring frame wall the leads are bent (typically, though not necessarily at 90 deg) so that they project into the interior volume of the ring frame for connection to the solid state devices.

Miller, Wilson N.↗

Vibrational Power Flow Analysis of Rods and Beams

A new method to model vibrational power flow and predict the resulting energy density levels in uniform rods and beams is investigated. This method models the flow of vibrational power in a manner analogous to the flow of thermal power in a heat conduction problem. The classical displacement solutions for harmonically excited, hysteretically damped rods and beams are used to derive expressions for the vibrational power flow and energy density in the rod and beam. Under certain conditions, the power flow in these two structural elements will be shown to be proportional to the energy density gradient. Using the relationship between power flow and energy density, an energy balance on differential control volumes in the rod and beam leads to a Poisson's equation which models the energy density distribution in the rod and beam. Coupling the energy density and power flow solutions for rods and beams is also discussed. It is shown that the resonant behavior of finite structures complicates the coupling of solutions, especially when the excitations are single frequency inputs. Two coupling formulations are discussed, the first based on the receptance method, and the second on the travelling wave approach used in Statistical Energy Analysis. The receptance method is the more computationally intensive but is capable of analyzing single frequency excitation cases. The traveling wave approach gives a good approximation of the frequency average of energy density and power flow in coupled systems, and thus, is an efficient technique for use with broadband frequency excitation.

Wohlever, James Christopher↗

Steady and transient least square solvers for thermal problems

This paper develops a hierarchical least square solution algorithm for highly nonlinear heat transfer problems. The methodology's capability is such that both steady and transient implicit formulations can be handled. This includes problems arising from highly nonlinear heat transfer systems modeled by either finite-element or finite-difference schemes. The overall procedure developed enables localized updating, iteration, and convergence checking as well as constraint application. The localized updating can be performed at a variety of hierarchical levels, i.e., degree of freedom, substructural, material-nonlinear groups, and/or boundary groups. The choice of such partitions can be made via energy partitioning or nonlinearity levels as well as by user selection. Overall, this leads to extremely robust computational characteristics. To demonstrate the methodology, problems are drawn from nonlinear heat conduction. These are used to quantify the robust capabilities of the hierarchical least square scheme.

Padovan, Joe↗

Thermal stresses of a wind turbine blade made of orthotropic material

This study is to investigate the thermal stress of a wind turbine blade made of wood composite material. First, the governing partial differential equation on heat conduction is stated, then, a finite element procedure using a variational approach is employed for the solution of the governing equation. Thus, the temperature field throughout the blade is determined. Next, based on the temperature field, a finite element procedure using potential energy approach is applied to determine the thermal stress field. A set of results is obtained through the use of a computer, which is considered to be satisfactory.

Fu, Kuan-Chen↗

Overview of numerical codes developed for predicted electrothermal deicing of aircraft blades

An overview of the deicing computer codes that have been developed at the University of Toledo under sponsorship of the NASA-Lewis Research Center is presented. These codes simulate the transient heat conduction and phase change occurring in an electrothermal deicier pad that has an arbitrary accreted ice shape on its surface. The codes are one-dimensional rectangular, two-dimensional rectangular, and two-dimensional with a coordinate transformation to model the true blade geometry. All modifications relating to the thermal physics of the deicing problem that have been incorporated into the codes will be discussed. Recent results of reformulating the codes using different numerical methods to increase program efficiency are described. In particular, this reformulation has enabled a more comprehensive two-dimensional code to run in much less CPU time than the original version. The code predictions are compared with experimental data obtained in the NASA-Lewis Icing Research Tunnel with a UH1H blade fitted with a B. F. Goodrich electrothermal deicer pad. Both continuous and cyclic heater firing cases are considered. The major objective in this comparison is to illustrate which codes give acceptable results in different regions of the airfoil for different heater firing sequences.

Keith, Theo G.↗

The loss and depth of CO2 ice in comet nuclei

The chemical differentiation of a homogeneous, initially unmantled cometary nucleus whose composition encompasses water ice, unclathrated CO2 ice, and silicate dust, is presently simulated by an analytical model to test whether the selective sublimation of any free CO2 ice in the new comet will produce a surface layer of water ice and dust covering the undifferentiated core. The model encompasses the effects of nucleus rotation, arbitrary rotatiuon axis orientation, latitude, nucleus heat conduction, and CO2 gas outflow restriction by the water ice and dust layer.

Fanale, Fraser P.↗