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 559 records · Page 31

Validation of NASA Thermal Ice Protection Computer Codes: The Validation of Antice - Part 3

An experimental program was generated by the Icing Technology Branch at NASA Glenn Research Center to validate two ice protection simulation codes: (1) LEWICE/Thermal for transient electrothermal de-icing and anti-icing simulations, and (2) ANTICE for steady state hot gas and electrothermal anti-icing simulations. An electrothermal ice protection system was designed and constructed integral to a 36 inch chord NACA0012 airfoil. The model was fully instrumented with thermo-couples, RTD'S, and heat flux gages. Tests were conducted at several icing environmental conditions during a two week period at the NASA Glenn Icing Research Tunnel. Experimental results of running-wet and evaporative cases were compared to the ANTICE computer code predictions and are presented in this paper.

Al-Khalil, Kamel M.↗

Prediction of High-Latitude Ionospheric Electrodynamics Using the Machine Learning Based Auroral Ionospheric Electrodynamics Model

We introduce a new framework for Machine-Learning (ML) based Auroral Ionosphere Model (ML-AIM). ML-AIM solves a current continuity equation by utilizing the ML model of Field Aligned Currents (FACs) of Kunduri et al., 2020 (https://doi.org/10.1029/2020JA027908), the FAC-derived aurora conductance model of Robinson et al., 2020 (https://doi.org/10.1029/2020JA028008), and the solar irradiance conductance model of Moen & Brekke (1993). The ML-AIM inputs are 60min time histories of solar wind plasma, interplanetary magnetic fields (IMF), and geomagnetic indices, and its outputs are ionospheric electric potential, electric fields, Pederson/Hall currents, and Joule Heating. We conduct two ML-AIM simulations for a weak geomagnetic activity on 14 May 2013 and a geomagnetic storm on 7-8 September 2017. ML-AIM produces reasonable ionospheric potential patterns such as two cell convection patterns and the enhancement of electric potentials during active times. The cross polar cap potential drop from ML-AIM is also comparable to the ones from the Weimer 2005 model, Super Dual Auroral Radar Network (SuperDARN), and Defense Meteorological Satellite Program (DMSP) F17 observations. ML-AIM is unique in a sense that it predicts ionospheric responses to the time-varying solar wind and geomagnetic conditions, while other traditional empirical model like Weimer 2005 is designed to provide static ionospheric conditions under steady solar wind/IMF conditions. In future, ML-AIM will include ML-based models of aurora precipitation and ionospheric conductance, improving its performance during active times.

H. K. Connor↗

Microsample Characterization of Coatings for Grcop-84 for High Heat Flux Applications

A multidisciplinary Johns Hopkins University-NASA Glenn team is undertaking a collaborative research program to elucidate and model the thermal stability and mechanical integrity of candidate coatings for GRCop-84. GRCop-84 is a high conductivity, high strength copper alloy that was recently developed at NASA Glenn for use in high temperature, high heat flux applications. With potential applications in rocket motor combustion chamber liners, nozzle ramps and other actively cooled structures, this new material offers great potential for decreasing weight and increasing reliability of third generation reusable launch vehicles. Current emphasis has turned toward the development of environmentally resistant and thermal barrier coatings for this alloy. Metallic coatings such as NiCrAlY and Cu-8-30%Cr have shown promise in: prohibiting blanching, reducing dog-house failures, increasing operating temperatures and decreasing cooling requirements. The focus of this research program is to develop a fundamental understanding of the substrate-coating interactions that occur during thermal cycling (inter-diffusion, viscoplasticity, morphological evolution, crack formation, etc.) and to derive a science-based protocol for future coating selection, optimization and reliability assurance. The microsample tensile testing approach adopted for this study allows us to characterize small-scale and highly scale-specific coatings and properties in a way not possible by conventional means. In addition to providing much needed design data, the integration of microsample testing with detailed microstructural observations provides a mechanistic foundation for coating optimization and life prediction modeling.

Hemker, Kevin↗

The effects of turbulence and stator/rotor interactions on turbine heat transfer. I - Design operating conditions

A combined experimental and analytical program was conducted to examine the effects of inlet turbulence, stator-rotor axial spacing, and relative circumferential spacing of first and second stators on turbine airfoil heat transfer. The experimental portion of the study was conducted in a large-scale (approximately 5X engine), ambient temperature, stage-and-a-half rotating turbine model. The data indicate that while turbine inlet turbulence can have a very strong impact on the first stator heat transfer, its impact in downstream rows is minimal. The effects on heat transfer produced by relatively large changes in stator/rotor spacing or by changing the relative row-to-row circumferential positions of stators were very small. Analytical results consist of airfoil heat transfer distributions computed with a finite-difference boundary layer code.

Blair, M. F.↗

Heat Transfer Analysis in Wire Bundles for Aerospace Vehicles

Design of wiring for aerospace vehicles relies on an understanding of "ampacity" which refers to the current carrying capacity of wires, either, individually or in wire bundles. Designers rely on standards to derate allowable current flow to prevent exceedance of wire temperature limits due to resistive heat dissipation within the wires or wire bundles. These standards often add considerable margin and are based on empirical data. Commercial providers are taking an aggressive approach to wire sizing which challenges the conventional wisdom of the established standards. Thermal modelling of wire bundles may offer significant mass reduction in a system if the technique can be generalized to produce reliable temperature predictions for arbitrary bundle configurations. Thermal analysis has been applied to the problem of wire bundles wherein any or all of the wires within the bundle may carry current. Wire bundles present analytical challenges because the heat transfer path from conductors internal to the bundle is tortuous, relying on internal radiation and thermal interface conductance to move the heat from within the bundle to the external jacket where it can be carried away by convective and radiative heat transfer. The problem is further complicated by the dependence of wire electrical resistivity on temperature. Reduced heat transfer out of the bundle leads to higher conductor temperatures and, hence, increased resistive heat dissipation. Development of a generalized wire bundle thermal model is presented and compared with test data. The steady state heat balance for a single wire is derived and extended to the bundle configuration. The generalized model includes the effects of temperature varying resistance, internal radiation and thermal interface conductance, external radiation and temperature varying convective relief from the free surface. The sensitivity of the response to uncertainties in key model parameters is explored using Monte Carlo analysis.

S L Rickman↗

The three-dimensional steady radial expansion of a viscous gas from a sonic source into a vacuum.

The three-dimensional steady radial expansion of a viscous, heat-conducting, compressible fluid from a spherical sonic source into a vacuum is analyzed using the Navier-Stokes equations as a basis. It is assumed that the model fluid is a perfect gas having constant specific heats, a constant Prandtl number of order unity, and viscosity coefficients varying as a power of the absolute temperature. Limiting forms for the flow variable solutions are studied for the Reynolds number based on the sonic source conditions, going to infinity and the Newtonian parameter both fixed and going to zero. For the case of the viscosity-temperature exponent between .5 and 1, it is shown that the velocity as well as the pressure approach zero as the radial distance goes to infinity. The formulations of the distinct regions that span the domain extending from the sonic source to the vacuum are presented.

Bush, W. B.↗

Structural and spectral studies of sunspots

Observations of umbral cores, both by multicolor photometry and by narrow band photometry in the vicinity of the sodium D lines, are described, and evidence is given which supports the validity of many umbral models, each of which describes different aspects of the observed umbral cores. Theoretical studies carried on at the observatory include the following: (1) Zeeman profiles of the sodium D sub 2 line and other lines; (2) turbulent heat conduction, sound waves, and the missing flux in sunspots; (3) chromospheric heating above spots by Alfven waves; (4) magnetic convection in the sun and solar neutrinos; (5) models of starspots on flare stars; (5) starspots on the primaries of contact binary systems; and (6) implications of starspots on red dwarfs.

Wyller, A. A.↗

Theoretical Limits of Vapor Cooling Large Cylindrical Structures Using Boil-Off

The Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) is planning to demonstrate the performance benefits of using boil-off to reduce the heat load on various upper stage structural members by using the boil-off vapor to intercept some of the heat load on the structures. A first order simplified analytical model was constructed in order to understand sensitivity of various parameters to performance as well as to understand the theoretical maximum performance of the vapor cooling of a thin-walled cylindrical structure. Constant fluid properties and a temperature dependent thermal conductivity of the structure material were assumed. The model showed that vapor cooling of an aluminum cylindrical structure with hydrogen can reduce the heat load along the cylinder by as much as 70%. Also, cooling a portion of the cylinder is quite effective in comparison with cooling it entirely.

R. Balasubramaniam↗

The Multi-Spectral Solar Telescope Array (MSSTA)

In 1987, our consortium pioneered the application of normal incidence multilayer X-ray optics to solar physics by obtaining the first high resolution narrow band, "thermally differentiated" images of the corona', using the emissions of the Fe IX/Fe X complex at ((lambda)lambda) approx. 171 A to 175 A, and He II Lyman (beta) at 256 A. Subsequently, we developed a rocket borne solar observatory, the Multi Spectral Solar Telescope Array (MSSTA) that pioneered multi-thermal imaging of the solar atmosphere, using high resolution narrow band X-ray, EUV and FUV optical systems. Analysis of MSSTA observations has resulted in four significant insights into the structure of the solar atmosphere: (1) the diameter of coronal loops is essentially constant along their length; (2) models of the thermal and density structure of polar plumes based on MSSTA observations have been shown to be consistent with the thesis that they are the source of high speed solar wind streams; (3) the magnetic structure of the footpoints of polar plumes is monopolar, and their thermal structure is consistent with the thesis that the chromosphere at their footpoints is heated by conduction from above; (4) coronal bright points are small loops, typically 3,500 - 20,000 km long (5 sec - 30 sec); their footpoints are located at the poles of bipolar magnetic structures that are are distinguished from other network elements by having a brighter Lyman a signature. Loop models derived for 26 bright points are consistent with the thesis that the chromosphere at their footpoints is heated by conduction from the corona.

Walker, A. B. C., Jr.↗

Energy balance in the solar transition region. I - Hydrostatic thermal models with ambipolar diffusion

The energy balance in the lower transition region is analyzed by constructing theoretical models which satisfy the energy balance constraint. The energy balance is achieved by balancing the radiative losses and the energy flowing downward from the corona. This energy flow is mainly in two forms: conductive heat flow and hydrogen ionization energy flow due to ambipolar diffusion. Hydrostatic equilibrium is assumed, and, in a first calculation, local mechanical heating and Joule heating are ignored. In a second model, some mechanical heating compatible with chromospheric energy-balance calculations is introduced. The models are computed for a partial non-LTE approach in which radiation departs strongly from LTE but particles depart from Maxwellian distributions only to first order. The results, which apply to cases where the magnetic field is either absent, or uniform and vertical, are compared with the observed Lyman lines and continuum from the average quiet sun. The approximate agreement suggests that this type of model can roughly explain the observed intensities in a physically meaningful way, assuming only a few free parameters specified as chromospheric boundary conditions.

Fontenla, J. M.↗

Validation of a Two-Phase CFD Model for Predicting Tank Self-Pressurization in the Ground-Based K-Site Experiment

A two-phase CFD model for self-pressurization of a cryogenic storage tank partially filled with liquid hydrogen is presented using the Volume-Of-Fluid approach for representing the phase boundary and the associated interfacial heat, mass and momentum transfer between the liquid and the vapor regions. The CFD model is validated against self-pressurization experiment performed in the K-site flightweight hydrogen storage tank at NASA Glenn Research Center. Laminar and turbulent simulations together with conjugated heat transfer analysis are performed. Effects of turbulence, as well as tank wall conduction are presented and discussed. Predicted tank pressures and fluid temperatures are compared with the experimental data at two different heat loads for validating the CFD model.

phase change↗

Numerical Modeling of Anti-icing Systems and Comparison to Test Results on a NACA 0012 Airfoil

A series of experimental tests were conducted in the NASA Lewis IRT on an electro-thermally heated NACA 0012 airfoil. Quantitative comparisons between the experimental results and those predicted by a computer simulation code were made to assess the validity of a recently developed anti-icing model. An infrared camera was utilized to scan the instantaneous temperature contours of the skin surface. Despite some experimental difficulties, good agreement between the numerical predictions and the experiment results were generally obtained for the surface temperature and the possibility for each runback to freeze. Some recommendations were given for an efficient operation of a thermal anti-icing system.

Al-Khalil, Kamel M.↗

Numerical modeling of anti-icing systems and comparison to test results on a NACA 0012 airfoil

A series of experimental tests were conducted in the NASA Lewis IRT on an electro-thermally heated NACA 0012 airfoil. Quantitative comparisons between the experimental results and those predicted by a computer simulation code were made to assess the validity of a recently developed anti-icing model. An infrared camera was utilized to scan the instantaneous temperature contours of the skin surface. Despite some experimental difficulties, good agreement between the numerical predictions and the experimental results were generally obtained for the surface temperature and the possibility for the runback to freeze. Some recommendations were given for an efficient operation of a thermal anti-icing system.

Al-Khalil, Kamel M.↗

Compatibility of Segments of Thermoelectric Generators

A method of calculating (usually for the purpose of maximizing) the power-conversion efficiency of a segmented thermoelectric generator is based on equations derived from the fundamental equations of thermoelectricity. Because it is directly traceable to first principles, the method provides physical explanations in addition to predictions of phenomena involved in segmentation. In comparison with the finite-element method used heretofore to predict (without being able to explain) the behavior of a segmented thermoelectric generator, this method is much simpler to implement in practice: in particular, the efficiency of a segmented thermoelectric generator can be estimated by evaluating equations using only hand-held calculator with this method. In addition, the method provides for determination of cascading ratios. The concept of cascading is illustrated in the figure and the definition of the cascading ratio is defined in the figure caption. An important aspect of the method is its approach to the issue of compatibility among segments, in combination with introduction of the concept of compatibility within a segment. Prior approaches involved the use of only averaged material properties. Two materials in direct contact could be examined for compatibility with each other, but there was no general framework for analysis of compatibility. The present method establishes such a framework. The mathematical derivation of the method begins with the definition of reduced efficiency of a thermoelectric generator as the ratio between (1) its thermal-to-electric power-conversion efficiency and (2) its Carnot efficiency (the maximum efficiency theoretically attainable, given its hot- and cold-side temperatures). The derivation involves calculation of the reduced efficiency of a model thermoelectric generator for which the hot-side temperature is only infinitesimally greater than the cold-side temperature. The derivation includes consideration of the ratio (u) between the electric current and heat-conduction power and leads to the concept of compatibility factor (s) for a given thermoelectric material, defined as the value of u that maximizes the reduced efficiency of the aforementioned model thermoelectric generator.

Snyder, G. Jeffrey↗

Skylab extravehicular mobility unit thermal simulator

The analytical methods, thermal model, and user's instructions for the Skylab Extravehicular Mobility Unit (SEMU) routine are presented. This digital computer program was developed for detailed thermal performance predictions of the SEMU on the NASA-JSC Univac 1108 computer system. It accounts for conductive, convective, and radiant heat transfer as well as fluid flow and special component characterization. The program provides thermal performance predictions for a 967 node thermal model in one thirty-sixth (1/36) of mission time when operated at a calculating interval of three minutes (mission time). The program has the operational flexibility to: (1) accept card or magnetic tape data input for the thermal model describing the SEMU structure, fluid systems, crewman and component performance, (2) accept card and/or magnetic tape input of internally generated heat and heat influx from the space environment, and (3) output tabular or plotted histories of temperature, flow rates, and other parameters describing system operating modes.

Hixon, C. W.↗

Candle Flames in Non-Buoyant Atmospheres

This paper addresses the behavior of a candle flame in a long-duration, quiescent microgravity environment both on the space Shuttle and the Mir Orbiting Station. On the Shuttle, the flames became dim blue after an initial transient where there was significant yellow (presumably soot) in the flame. The flame lifetimes were typically less than 60 seconds. The safety-mandated candlebox that contained the candle flame inhibited oxygen transport to the flame and thus limited the flame lifetime. The flames on the Mir were similar, except that the yellow luminosity persisted longer into the flame lifetime because of a higher initial oxygen concentration, The Mir flames burned for as long as 45 minutes. The difference in the flame lifetime between the Shuttle and Mir flames was primarily the redesigned candlebox that did not inhibit oxygen transport to the flame. In both environments, the flame intensity and the height-to-width ratio gradually decreased as the ambient oxygen content in the sealed chamber slowly decreased. Both sets of experiments showed spontaneous, axisymmetric flame oscillations just prior to extinction. The paper also presents a numerical model of a candle flame. The formulation is two-dimensional and time-dependent in the gas phase with constant specific heats, thermal conductivity and Lewis number (although different species can have different Lewis numbers), one-step finite-rate kinetics, and gas-phase radiative losses from CO2 and H2O. The treatment of the liquid/wick phase assumes that the, fuel evaporates from a constant diameter sphere connected to an inert cone. The model predicts a steady flame with a shape and size quantitatively similar to the Shuttle and Mir flames. The computation predicts that the flame size will increase slightly with increasing ambient oxygen mole fraction. The model also predicts pre-extinction flame oscillations if the rate of decrease in ambient oxygen is small enough, such as that which would occur for a flame burning in a sealed ambient.

Dietrich, D. L.↗

Calculations of the moon's heat history at different concentrations of radioactive elements taking account of the material differentiation with melting

A mathematical procedure for analyzing the heat conductivity of the lunar surface is discussed. The solution is based on homogeneous and laminated moon models and considers the effects of radioactive elements conveyed to the lunar surface by melting. The various parameters which introduce uncertainties into the numerical analysis are identified. The application of data obtained from radio astronomy and from analyses of lunar samples returned by the Apollo flights is explained. Tables of data are included to show the types and amounts of radioactive materials which have been identified.

Arnatskaya, O. I.↗