Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat conduction modeling”

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,009 records · Page 56

Conductive Sphere in a Radio Frequency Field: Theory and Applications to Positioners, Heating, and Noncontact Measurements

An electrically conductive spherical sample located in an electromagnetic field excited by rf (radio frequency) current in a system of coaxial coils is treated theoretically. Maxwell's equations are solved exactly and all integrals in the formulas for the fields are evaluated analytically for the case where the sphere is on the axis and the coil system is modeled by a stack of filamentary circular loops. Formulas are also derived for electromagnetic force exerted on the sphere, excess impedance in the coil system due to the presence of the sphere, and power absorbed by the sphere. All integrals in those formulas have been evaluated analytically. Force measurements are presented and they are in excellent agreement with the new theory. A low-power electromagnetic levitator that is accurately described by the theory has been demonstrated and is discussed. Experimental measurements of excess impedance are presented and compared with theory, and those results are used to demonstrate an accurate noncontact method for determining electrical conductivity. Theoretical formulas for power absorption are evaluated numerically and their usefulness in both rf heating and in making noncontact measurements of a number of thermophysical properties of materials is discussed.

Jackson, H. W.↗

The energetics of Titan's ionosphere

We have developed a comprehensive model to study the dynamics and energetics of the ionosphere of Titan. We solved the one-dimensional, time-dependent, coupled continuity and momentum equations for several ion species, together with single ion and electron energy equations, in order to calculate density, velocity, and temperature profiles. Calculations were carried out for several cases corresponding to different local times and configurations of the Titan-Saturn system. In our model the effects of horizontal magnetic fields were assumed to be negligible, except for their effect on reducing the electron and ion thermal conductivities and inhibiting vertical transport in the subram region. The ionospheric density peak was found to be at an altitude of about 1100 km, in accordance with earlier model calculations. The ionosphere is chemically controlled below an altitude of about 1500 km. Above this level, ion densities differ significantly from their chemical equilibrium values due to strong upward ion velocities. Heat is deposited in a narrow region around the ionospheric peak, resulting in temperature profiles increasing sharply and reaching nearly constant values of 800-1000 deg K for electrons and 300 deg K for ions in the topside, assuming conditions appropriate for the wake region. In the subram region magnetic correction factors make the electron heat conductivities negligible, resulting in electron temperatures increasing strongly with altitude and reaching values in the order of 5000 deg K at our upper boundary located at 2200 km. Ion chemical heating is found to play an important role in shaping the ion energy balance in Titan's ionosphere.

Roboz, A.↗

Heat transfer tests of the McDonnell-Douglas delta wing orbiter mated with -17A booster at Mach number 8

A wind tunnel test program to determine aerodynamic heat transfer distributions on the McDonnell-Douglas configurations is reported. The tests were conducted at the Arnold Engineering Development Center (AEDC) in Tunnel B of the von Karman Gas Dynamics Facility (VKF). Heat-transfer rates were determined by the phase-change paint technique on 0.011-scale Stycast models using Tempilaq as the surface temperature indicator. The nominal test conditions were; Mach 8, freestream unit Reynolds numbers of 0.8 x one million, 2.5 x one million, and 3.7 x one million, and angles of attack of -5 deg, 0 deg, +5deg. Model details, test conditions, phase-change paint photographs and reduced heat-transfer coefficients are presented.

Matthews, R. K.↗

Results of tests of a Rockwell International space shuttle orbiter (-139 configuration) 0.0175-scale model (no. 29-0) in AEDC tunnel F to determine hypersonic heating effects (OH11)

Results from wind tunnel tests to determine hypersonic aerodynamic heating rates on a NASA/Rockwell Space Shuttle Orbiter are reported. The tests were to determine Mach number effects, if any, and to obtain overall heating rate data at high Mach numbers from 10.5 to 16. The model used was a 0.0175-scale model built to Rockwell Orbiter lines VL70-000139. The model identity number is 29-0. These tests, designated OH11, were conducted in the AEDC Tunnel F.

Quan, M.↗

Effect of propellant deformation on ignition and combustion processes in solid propellant cracks

A comprehensive theoretical model was formulated to study the development of convective burning in a solid propellant crack which continually deforms due to burning and pressure loading. In the theoretical model, the effect of interrelated structural deformation and combustion processes was taken into account by considering (1) transient, one dimensional mass, momentum, and energy conservation equations in the gas phase; (2) a transient, one dimensional heat conduction equation in the solid phase; and (3) quasi-static deformation of the two dimensional, linear viscoelastic propellant crack caused by pressure loading. Partial closures may generate substantial local pressure peaks along the crack, implying a strong coupling between chamber pressurization, crack combustion, and propellant deformation, especially when the cracks are narrow and the chamber pressurization rates high. The maximum pressure in the crack cavity is generally higher than that in the chamber. The initial flame-spreading process is not affected by propellant deformation.

Kumar, M.↗

Comparison of Heat Transfer from Airfoil in Natural and Simulated Icing Conditions

An investigation of the heat transfer from an airfoil in clear air and in simulated icing conditions was conducted in the NACA Lewis 6- by 9-foot icing-research tunnel in order to determine the validity of heat-transfer data as obtained in the tunnel. This investiation was made on the same model NACA 65,2-016 airfoil section used in a previous flight study, under similar heating, icing, and operating conditions. The effect of tunnel turbulence, in clear air and in icingwas indicated by the forward movement of transition from laminar to turbulent heat transfer. An analysis of the flight results showed the convective heat transfer in icing to be considerably different from that measured in clear air and. only slightly different from that obtained in the icing-research tunnel during simulated icing.

Gelder, Thomas F.↗

Experimental feasibility study of a thermoelectric heat flux gage

An experiment was conducted to determine the feasibility of using a commercially available thermoelectric device as a heat flux gage at near ambient conditions. In certain research applications, the thermoelectric heat flux gage can provide a relatively simple means to model a warm fluid cold wall convection environment. The experiment showed that heat flux through the gage could be correlated within 2.5 percent with a simple algebraic equation which considered the thermoelectric current through the device and the hot and cold side temperatures.

Vanfossen, G. J.↗

Effect of an oscillating flow direction on leading edge heat transfer

An experimental investigation was conducted to examine the effect of a periodic variation in the angle of attack on heat transfer at the leading edge of a gas turbine blade. A circular cylinder was used as a large-scale model of the leading edge region. The cylinder was placed in a wind tunnel and was oscillated rotationally about its axis. The incident flow Reynolds number and the Strouhal number of oscillation were chosen to model an actual turbine condition. Incident turbulence levels up to 4.9 percent were produced by grids placed upstream of the cylinder. The transfer rate was measured using a mass transfer technique and heat transfer rates inferred from the results. A direct comparison of the unsteady and steady results indicate that the effect is dependent on the Strouhal number, turbulence level, and the turbulence length scale, but that the largest observed effect was only a 10 percent augmentation at the nominal stagnation position.

Marziale, M. L.↗

A spacecraft cooling system for a charged coupled device

This paper describes the thermal analysis, design, and testing of a dedicated cooling system for a Spartan spacecraft payload. A simple reliable design that requires minimum power consumption and minimum weight was developed. The payload has a CCD detector that must be maintained at a temperature of approximately -40 C or colder. The cooling system consists of a fin radiator, dual redundant heat pipes, and a thermal electric device (TED). The system was analytically modeled through the use of the Simplified Shuttle Payload Thermal Analyzer (SSPTA) computer program. A thermal test of the system simulating flight conditions was conducted to correlate the computer model and verify performance specifications.

Walker, Mary S.↗

Natural melting within a spherical shell

Fundamental heat transfer experiments were performed on the melting of a phase change medium in a spherical shell. Free expansion of the medium into a void space within the sphere was permitted. A step function temperature jump on the outer shell wall was imposed and the timewise evolution of the melting process and the position of the solid-liquid interface was photographically recorded. Numerical integration of the interface position data yielded information about the melted mass and the energy of melting. It was found that the rate of melting and the heat transfer were significantly affected by the movement of the solid medium to the base of the sphere due to gravity. The energy transfer associated with melting was substantially higher than that predicted by the conduction model. Furthermore, the radio of the measured values of sensible energy in the liquid melt to the energy of melting were nearly proportional to the Stefan number. The experimental results are in agreement with a theory set forth in an earlier paper.

Bahrami, Parviz A.↗

Thermographic inspection of fatigue cracks in riveted plates

A theoretical and experimental analysis of the inspection of reiveted plates using the thermoelastic effect is presented. Theoretically, a model is employed to estimate the thermoelastic response expected from a cracked and uncracked rivet hole. The stress invariant, which causes the thermoelastic response, is calculated for the cracked case using an FEM model and for the uncracked case using a theoretical approximation. The differences in the two thermal responses are reduced when the effects of optical spatial averaging, thermal conduction, and a concentric blind area representing the rivet heat are considered. It is suggested that discrimination between the cracked and uncracked cases for short cracks will improve with increasing frequency of stress oscillation. Measurements are shown for aerospace aluminum samples with and without fatigue cracks propagating from holes in a plate. The effect of mean stress on the thermoelastic response is illustrated, and the results are qualitatively compared with the finite element analysis.

Welch, Christopher S.↗

The microscopic state of the solar wind: Links between composition, velocity distributions, waves and turbulence

An overview is given of the microscopic state of the solar wind with emphasis on recent Ulysses high-latitude observations and previous Helios in-ecliptic observations. The possible links between composition, ionization state. velocity distribution functions of electrons, protons and heavy ions. kinetic plasma waves and MHD-scale turbulence are elaborated. Emphasis is placed on a connection of interplanetary kinetic-scale phenomena with their generating microscopic processes in the corona. The fast streams seem to consist of mesoscale pressure-balanced plasma filaments and magnetic flux tubes, reminiscent of the supergranular-size structures building the open corona, from which copious Alfven waves emanate. The wind from the magnetically structured and active corona shows developed compressive turbulence and considerable abundance and ionization state variations. Some modelling attempts to explain the observed element fractionation are briefly reviewed. The causes of the nonthermal particle features. such as proton-ion differential streaming, ion beams. temperature anisotropies, and skewed distributions associated with collisionless heat conduction, are ultimately to be searched in the fact, that the corona is never quiescent but fundamentally variable in space and time. Consequently, the radial evolution of the internal state of the wind resembles at all latitudes a complicated relaxation process, in the course of which the free (in comparison with LTE conditions) particle kinetic energy is converted into plasma waves and MHD turbulence on a wide range of scales. This leads to intermittent wave-particle interactions and unsteady anomalous transport, mixed with the weak effects of the rare Coulomb collisions. Spherical expansion and large-scale inhomogeneity forces the wind to attain microscopically a complex internal state of dynamic statistical equilibrium.

Marsch, E.↗

Method of Forming a Hot Film Sensor System on a Model

A method of forming a hot film sensor directly on a model is provided. A polyimide solution is sprayed onto the model. The model so sprayed is then heated in air. The steps of spraying and heating are repeated until a polyimide film of desired thickness is achieved on the model. The model with the polyimide film thereon is then thoroughly dried in air. One or more hot film sensors and corresponding electrical conducting leads are then applied directly onto the polyimide film.

Tran, Sang Q.↗

Using In-Situ Deposition of Metallic Thin Films on Mars to Monitor Atmospheric H20

We describe the use of a novel sensor, designed to characterize the reactive nature of the Martian atmosphere, to also characterize the instantaneous abundance of atmospheric H,O. The sensor deposits, in situ, a thin silver film onto a sapphire substrate and monitors oxidation by measuring resistance of the Ag film during both deposition and subsequent oxidation [1,2]. . The evaporation source is placed in the center of a hollow, open-ended tube. On flash heating the source, evaporated metallic silver rapidly reacts with oxidizing gases in the tube, depositing on the tube walls as a resisitive Ag oxide film. Unoxidized silver deposits over the Ag oxides, once the oxidizing gases in the vicinity of the source have been consumed. Since the ends of the tube are open, atmospheric gases in the tube cause a time delay as the evaporated silver initially reacts with the available O2, and H2O; once oxidizing gases in the tube are significantly depleted, a silver film closes the chemiresistor circuit. In principle this operation mode provides a repeatable measure of the variable H2O abundance, since O2, levels are constant in the atmosphere. To explore the utility of the instrument as an H2O sensor, predict its behavior, and identify possible failure modes, we developed a numerical model of the instrument, and exercised it for Mars conditions. The model predicts the behavior of the electrical circuit, the temperature dependence of resistivity for each component, and the resultant Joule heating. The model balances Joule heating against radiation, sensible heat loss to ambient CO2, and, latent heat loss from sublimating Ag. The flux of Ag atoms from the source is tracked continuously. In the gas phase, the number density of Ag, O2, and H2O is calculated, along with their reaction rates on collision. The model suggests that Ag substantially depletes H2O only near the center of the tube. The model supports the hypothesis that the onset of high conductivity between the A u electrodes is a function of the H2O abundance in the ambient gas, but only to a limit; when H2O abundances are buffered at or above 240K (an H2O abundance that would not have occurred on Mars), metallic Ag never deposits in excess of Ago, and the circuit does not close. A significant difference in the Mars simulation is that the Ag source sublimates more quickly, primarily because of differences in sensible heat fluxes from the source. This allows subsequent heating of the W filament used to support the Ag source. Had the originally-designed experimental sequence been carried out on Mars, the filament might have melted, making subsequent investigations impossible.

Zent, A. P.↗

Measurement of Effective Canopy Temperature: The Missing Link to Modeling Transpiration in Controlled Environments

The objectives were to apply energy balance principles to plant canopies, and to determine which parameters are essential for predicting plant canopy transpiration (E) in controlled environments. Transpiration was accurately measured in a gas-exchange system. Absorbed radiation (R(sub abs)) by the canopy was measured with a net radiometer and calculated from short and long-wave radiation components. Average canopy foliar temperature T(sub L) can be measured with an infrared radiometer, but since T(sub L) is seldom uniform, a weighed average measurement of T(sub L) must be made. The effective canopy temperature T(sub C) is that temperature that balances the energy flux between absorbed radiation and latent heat L(sub E) and sensible heat (H) fluxes. TC should exactly equal air temperature T(sub A) when L(sub E) equals R(sub abs). When unnecessary thermal radiation from the lighting system is removed by a water filter, the magnitude of L(sub E) from transpiration approaches Rabs and T(sub C) is close to T(sub A). Unlike field models, we included the energy used in photosynthesis and found that up to 10% of Rabs was used in photosynthesis. We calculated aerodynamic conductance for H from measurements of wind speed and canopy height using the wind profile equation. Canopy aerodynamic conductance ranged from.03 to.04 m/s for wind speeds from.6 to 1 m/s; thus a 0.1 C canopy to air temperature difference results in a sensible heat flux of about 4 W/sq m, which is only 1% of R(sub abs). We examined the ability of wide angle infrared transducers to accurately integrate T(sub L) from the top to the bottom of the canopy. We measured evaporation from the hydroponic media to be approximately 1 micro mol/sq m s or 10% of R(sub abs). This result indicates that separating evaporation from transpiration is more important than exact measurement of canopy temperature.

Monje, O. A.↗

Modeling of Aerobrake Ballute Stagnation Point Temperature and Heat Transfer to Inflation Gas

A trailing Ballute drag device concept for spacecraft aerocapture is considered. A thermal model for calculation of the Ballute membrane temperature and the inflation gas temperature is developed. An algorithm capturing the most salient features of the concept is implemented. In conjunction with the thermal model, trajectory calculations for two candidate missions, Titan Explorer and Neptune Orbiter missions, are used to estimate the stagnation point temperature and the inflation gas temperature. Radiation from both sides of the membrane at the stagnation point and conduction to the inflating gas is included. The results showed that the radiation from the membrane and to a much lesser extent conduction to the inflating gas, are likely to be the controlling heat transfer mechanisms and that the increase in gas temperature due to aerodynamic heating is of secondary importance.

Bahrami, Parviz A.↗

The importance of properties in modeling

Casting and welding of superalloys, stainless steel and titanium alloys are processes which can be improved through modeling of heat flow, fluid flow, residual stress development, and microstructural evolution. These simulations require inputs of thermophysical data, some of which involves the partially or totally liquid state. In particular, these processes involve melting, flow in the liquid, and solidification. Modeling of such processes can lead to an improved understanding of defects such as shrinkage, inclusions, cracks, incomplete filling (or penetration), macrosegregation, improper grain structure, and deviations from dimensional specifications. Effective modeling can shorten process development time and improve quality. An approach to these problems is to develop efficient models; validate through correlations with thermal, distortion, and microstructural data; run parametric studies; extract knowledge based rules; and apply to adaptive closed loop control systems. With the appropriate pre- and post-processing, such analyses can be made 'user friendly'. This would include graphical user interfaces as well as realistic images and color maps. In such form, these models can be used for sensitivity analyses, which are useful in defining appropriate sensors and in the development of control strategies. Such modeling can be done at several levels, e.g., the MARO level, modeling large scale phenomena such as heat and fluid flow or material deformation; the MICRO level, modeling the development of dendrites, grains or precipitates; or at the NANO level, modeling point defects, dislocations, stacking faults, etc. There are many computational issues associated with these simulations, e.g. computational efficiency and accuracy. In addition, there are many materials issues, not the least of which is the availability of accurate high temperature thermophysical data for complex alloys. This would include latent heat of fusion, temperature dependent heat capacity and thermal conductivity (for liquid and solid), viscosity, surface tension, thermal expansion, mechanical properties, etc. Preliminary data is frequently gathered from the literature; however, this is often not available for modern alloys. If additional data are required, measurements can be used; however, these are costly, time consuming and can be erroneous due to a lack of testing standards or impure materials. Microstructural predictors can be extracted from thermal information, e.g. cooling rate and thermal gradient; the prediction of microstructure is dependent on solidus and liquidus temperature, mushy zone permeability, the solidification curve, volume changes, phase transformations, alloying effects (such as surface tension or viscosity), mold/metal reactions, metal/environment reactions, etc. Defect maps may be needed to predict the onset of shrinkage, hot cracking or 'freckling'. Constants may be needed for stress relaxation, dendrite coarsening, vaporization, etc. Visualization was used as a tool to better comprehend complex data sets associated with the analysis of directional solidification (including crystal growth) and welding. Examples include not only isotherms, but also cooling rate, growth rate and thermal gradient. The latter two are not single valued scalars, but rather time and space dependent vector fields. Efficient models were developed for both casting and welding to predict heat flow and the relationship to dendrite and grain growth. These codes include many of the non-linear effects, e.g. radiation, which dominate these processes. The home-built FDM code(s) were designed to be useful not only to the scientist, but also to the process engineer. Special output can be requested to compare directly to experimental data. Visualization procedures were developed to visualize critical results, e.g. fusion zone width at the surface opposite that where the arc is applied ('penetration'). Both elaborate and simplified distortion analyses were carried out. It is clear that extensive mechanical property data are critical in order to accurately predict residual stress patterns. A scheme is currently being developed to integrate these modeling tools into a set of control algorithms; however, the success of this approach is critically dependent on the availability of accurate high temperature thermophysical data.

Giamei, A. F.↗

Uncertainty Determination for Aeroheating in Uranus and Saturn Probe Entries by the Monte Carlo Method

The 2013-2022 Decaedal survey for planetary exploration has identified probe missions to Uranus and Saturn as high priorities. This work endeavors to examine the uncertainty for determining aeroheating in such entry environments. Representative entry trajectories are constructed using the TRAJ software. Flowfields at selected points on the trajectories are then computed using the Data Parallel Line Relaxation (DPLR) Computational Fluid Dynamics Code. A Monte Carlo study is performed on the DPLR input parameters to determine the uncertainty in the predicted aeroheating, and correlation coefficients are examined to identify which input parameters show the most influence on the uncertainty. A review of the present best practices for input parameters (e.g. transport coefficient and vibrational relaxation time) is also conducted. It is found that the 2(sigma) - uncertainty for heating on Uranus entry is no more than 2.1%, assuming an equilibrium catalytic wall, with the uncertainty being determined primarily by diffusion and H(sub 2) recombination rate within the boundary layer. However, if the wall is assumed to be partially or non-catalytic, this uncertainty may increase to as large as 18%. The catalytic wall model can contribute over 3x change in heat flux and a 20% variation in film coefficient. Therefore, coupled material response/fluid dynamic models are recommended for this problem. It was also found that much of this variability is artificially suppressed when a constant Schmidt number approach is implemented. Because the boundary layer is reacting, it is necessary to employ self-consistent effective binary diffusion to obtain a correct thermal transport solution. For Saturn entries, the 2(sigma) - uncertainty for convective heating was less than 3.7%. The major uncertainty driver was dependent on shock temperature/velocity, changing from boundary layer thermal conductivity to diffusivity and then to shock layer ionization rate as velocity increases. While radiative heating for Uranus entry was negligible, the nominal solution for Saturn computed up to 20% radiative heating at the highest velocity examined. The radiative heating followed a non-normal distribution, with up to a 3x variation in magnitude. This uncertainty is driven by the H(sub 2) dissociation rate, as H(sub 2) that persists in the hot non-equilibrium zone contributes significantly to radiation.

Palmer, Grant↗