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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Experimental Aerothermodynamics In Support Of The Columbia Accident Investigation

The technical foundation for the most probable damage scenario reported in the Columbia Accident Investigation Board's final report was largely derived from synergistic aerodynamic/aerothermodynamic wind tunnel measurements and inviscid predictions made at NASA Langley Research Center and later corroborated with engineering analysis, high fidelity numerical viscous simulations, and foam impact testing near the close of the investigation. This report provides an overview of the hypersonic aerothermodynamic wind tunnel program conducted at NASA Langley and illustrates how the ground-based heating measurements provided early insight that guided the direction and utilization of agency resources in support of the investigation. Global surface heat transfer mappings, surface streamline patterns, and shock shapes were measured on 0.0075 scale models of the Orbiter configuration with and without postulated damage to the thermal protection system. Test parametrics include angle of attack from 38 to 42 degs, sideslip angles of 38 to 42 degs, sideslip angles of plus or minus 1 deg, Reynolds numbers based upon model length from 0.05 x 10(exp 6) to 6.5 x 10(exp 6), and normal shock density ratios of 5 (Mach 6 Air) and 12 (Mach 6 CF4). The primary objective of the testing was to provide surface heating characteristics on scaled Orbiter models with outer mold line perturbations to simulate various forms of localized surface damage to the thermal protection system. Initial experimental testing conducted within two weeks of the accident simulated a broad spectrum of thermal protection system damage to the Orbiter windward surface and was used to refute several hypothesized forms of thermal protection system damage, which included gouges in the windward thermal protection system tiles, breaches through the wing new the main landing gear door, and protuberances along the wing leading edge that produced asymmetric boundary layer transition. As the forensic phase of the investigation developed and the condition of recovered debris was examined, increasing emphasis was placed on identifying wing leading edge damage (partially and fully missing reinforced carbon-carbon panels, and eventually holes in the wing leading edge with venting to the wing upper surface) that produced off-nominal heating trends consistent with extracted Orbiter flight recorder temperature data.

Horvath, Thomas J.↗
NASA NTRS2005

Energy Transport in the Thermosphere During the Solar Storms of April 2002

The dramatic solar storm events of April 2002 deposited a large amount of energy into the Earth's upper atmosphere, substantially altering the thermal structure, the chemical composition, the dynamics, and the radiative environment. We examine the flow of energy within the thermosphere during this storm period from the perspective of infrared radiation transport and heat conduction. Observations from the SABER instrument on the TIMED satellite are coupled with computations based on the ASPEN thermospheric general circulation model to assess the energy flow. The dominant radiative response is associated with dramatically enhanced infrared emission from nitric oxide at 5.3 microns from which a total of approx. 7.7 x 10(exp 23) ergs of energy are radiated during the storm. Energy loss rates due to NO emission exceed 2200 Kelvin per day. In contrast, energy loss from carbon dioxide emission at 15 microns is only approx. 2.3% that of nitric oxide. Atomic oxygen emission at 63 microns is essentially constant during the storm. Energy loss from molecular heat conduction may be as large as 3.8% of the NO emission. These results confirm the "natural thermostat" effect of nitric oxide emission as the primary mechanism by which storm energy is lost from the thermosphere below 210 km.

SOLAR STORMS↗
NASA NTRS1995

Solar wind theory and modelling

The outflow of coronal plasma into interplanetary space is a consequence of the coronal heating process. Therefore the formation of the corona and the acceleration of the solar wind should be treated as a single problem. Traditionally the mass or particle flux emanating from the extended corona has been thought of as being determined by the coronal temperature or scale height and the coronal (base) density. This argument follows from considerations of the momentum balance of the corona-wind system from which one obtains models of a close to hydrostatic corona out to the critical point where the flow becomes supersonic. With this approach to the acceleration of the wind is has been difficult to reconcile the relatively small variation observed in the proton flux at 1 AU with the predicted exponential dependence of the proton flux on the coronal temperature. In this talk we would like to emphasize another approach in which coronal energetics play the primary role. The deposition of energy into the corona through some 'mechanical' energy flux is balanced by the various energy sinks available to the corona and the sum of these processes determine the coronal structure, i.e. its temperature and density. The corona loses energy through heat conduction into the transition region, through radiative losses, and through the gravitational potential energy and kinetic energy put into the solar wind itself. We will show from a series of models of the chromosphere transition region-corona-solar wind system that most of the energy deposited in a magnetically open region will go into the solar wind, with roughly half going into kinetic energy and half into lifting the plasma out of the solar gravity field. The coronal base density will adjust itself in such a way that the heat conductive flux flowing into the transition region is radiated away in the upper chromosphere. The coronal temperature is set by the requirements that most of the deposited energy goes into accelerating the solar wind; the coronal scale height will adjust itself so that the solar wind energy losses conform to the amplitude of the input energy. These processes are modified by the 'mode' of energy deposition, and we will show the effects on coronal structure of changing the parameters describing coronal heating as well as the effects of including a helium fluid in the models. However, the location, scale height and/or form of the energy deposition (i.e. heating or direct acceleration) are not too important for the solar wind, the coronal density and temperature structure will vary with the 'mode' of energy deposition, but the solar wind mass flux depends mainly on the amplitude of the energy flux.

Hansteen, Viggo H.↗
DOE OSTI2023

COMSOL Models and Simulation Results for Closed-Loop Geothermal Systems with Thermal Enhancements

The Closed-Loop Geothermal Working Group and the National Renewable Energy Laboratory (NREL) have studied the performance of disc-shaped and linear thermal enhancements (conductive material introduced around a wellbore) in closed-loop geothermal systems. COMSOL models were developed to investigate the increase in heat production in systems with thermal enhancements compared to those without, across a range of scenarios. This dataset includes two example COMSOL models, a summary of the simulation results, and a link to an accompanying paper. Studied are idealized versions of two thermal enhancement designs: (1) existing or newly created fractures surrounding the wellbore, filled with highly conductive material, and (2) a 'fishbone structure,' where short wellbores are drilled radially outward from the main wellbore and filled with highly conductive material.

15 GEOTHERMAL ENERGY↗
NASA NTRS1988

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↗
DOE OSTI2023

Self-consistent integrated modeling of combined hybrid discharge-laser produced plasma devices for extreme ultraviolet metrology

Discharge- and laser-produced plasma (DLPP) devices are being used as light sources for extreme ultraviolet (EUV) generation. A key challenge for both, DPP and LPP, is achieving sufficient brightness to support the throughput requirements of nanometrology tools. To simulate the environment of a hybrid DLPP device and optimize the EUV output, we have developed an integrated HEIGHTS-DLPP computer simulation package. The package integrates simulation of two evolving plasmas (DPP and LPP) and includes modeling of a set of integrated self-consistent processes: external power source and plasma energy balance, plasma resistive magnetohydrodynamics, plasma heat conduction, detailed radiation transport (RT), and laser absorption and refraction. We simulated and optimized DLPP devices using Xe gas as a target material. We synchronized the external circuit parameters, chamber gas parameters, and laser beam temporal and spatial profiles to achieve maximum EUV output. The full 3D Monte Carlo scheme was integrated for detailed RT and EUV output calculations in Xe using more than 3600 spectral groups. The modeling results are in good agreement with Julich Forschungszentrum experimental data. Theoretical models, developed and integrated into the HEIGHTS package, showed wide capabilities and flexibility. In conclusion, the models and package can be used for optimization of the experimental parameters and settings, investigation of DLPP devices with complex design, analyzing the impact of integrated spatial effects and working timeline arrangement on the final EUV output, and EUV source size, shape, and angular distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1989

Thermal modelling of various thermal barrier coatings in a high heat flux rocket engine

Traditional Air Plasma Sprayed (APS) ZrO2-Y2O3 Thermal Barrier Coatings (TBC's) and Low Pressure Plasma Sprayed (LPPS) ZrO2-Y2O3/Ni-Cr-Al-Y cermet coatings were tested in a H2/O2 rocked engine. The traditional ZrO2-Y2O3 (TBC's) showed considerable metal temperature reductions during testing in the hydrogen-rich environment. A thermal model was developed to predict the thermal response of the tubes with the various coatings. Good agreement was observed between predicted temperatures and measured temperatures at the inner wall of the tube and in the metal near the coating/metal interface. The thermal model was also used to examine the effect of the differences in the reported values of the thermal conductivity of plasma sprayed ZrO2-Y2O3 ceramic coatings, the effect of 100 micron (0.004 in.) thick metallic bond coat, the effect of tangential heat transfer around the tube, and the effect or radiation from the surface of the ceramic coating. It was shown that for the short duration testing in the rocket engine, the most important of these considerations was the effect of the uncertainty in the thermal conductivity of temperatures (greater than 100 C) predicted in the tube. The thermal model was also used to predict the thermal response of the coated rod in order to quantify the difference in the metal temperatures between the two substrate geometries and to explain the previously-observed increased life of coatings on rods over that on tubes. A thermal model was also developed to predict heat transfer to the leading edge of High Pressure Fuel Turbopump (HPFTP) blades during start-up of the space shuttle main engines. The ability of various TBC's to reduce metal temperatures during the two thermal excursions occurring on start-up was predicted. Temperature reductions of 150 to 470 C were predicted for 165 micron (0.0065 in.) coatings for the greater of the two thermal excursions.

Nesbitt, James A.↗
DOE OSTI2022

High-Temperature Latent Heat Storage System using Transportable Heat Pipes for Versatile Integration with Emerging Microreactors

With the growing interest in the development of micro nuclear reactors and their application to microgrid, high-temperature thermal energy storage (TES) has also been receiving great attention as a technology to reinforce the energy utilization and profitability of this class of reactors. Given the importance and research needs of the TES technology, a research project is underway at Idaho National Laboratory (INL) to develop high-temperature TES that can be easily interfaced with emerging micro nuclear reactors. Currently under development is the Heat pipe-Integrated Thermal Battery (HITB), a novel latent heat storage design that enables flexible and robust integration with any design of currently considered micro nuclear reactors. The development of detailed design based on the proposed TES concept, HITB, requires not only development and validation of modeling and simulation capabilities, but also research that spans thermal performance analysis, component- and system-scale design optimization, and experimental demonstrations including high-temperature feasibility testing. The current research efforts (and plans) are divided into the following four areas: (i) Thermal and structural modeling and analysis at both system- and component-scales for evaluating performance and conducting design optimization, (ii) Small-scale experimental demonstrations for the key components of the proposed HITB design, (iii) Modelica-based dynamic model development for developing a control strategy under dynamic operation within integrated energy systems configuration, (iv) Integrated experimental testing and demonstration. This paper introduces the overall research efforts that have been and are being conducted to develop and demonstrate the novel design of the high-temperature latent heat storage device, HITB. The key design features, ongoing research efforts including the development of modeling and simulation capabilities, small-scale demonstrations of the HITB’s key functionalities prior to scaled-up experimental demonstration, are discussed. Then, this paper concludes with a remark on the future research plans for this project.

25 ENERGY STORAGE↗
DOE OSTI2022

Analytical benchmark solution for 1-D neutron transport coupled with thermal conduction and material expansion

In this work we present fully analytical solutions for a class of finite, homogeneous, 1-D slab benchmark problems with nonlinear temperature feedback effects. The proposed class of benchmarks include multiplicative 1-D neutron transport (limited to quasistatic S{sub 2} with μ = ±1) coupled with thermal conduction, convection, Doppler broadening, and expansion effects along the length of the slab. This class of benchmark models, along with the corresponding analytical solutions, are valuable for validating multiphysics analysis frameworks that support coupled neutronics/thermal/structural calculations. Analytical solutions for the benchmarks are obtained by introducing an ansatz that the equilibrium flux and temperature distributions in the slab have the same shape. Specific values for the total microscopic cross section, σ{sub t,0}, and conductive heat transfer coefficient, h, that satisfy the assumed ansatz are then determined. A discussion of the procedure for generating benchmark models and analytical solutions is provided, along with numerical results for an example set of model parameters and thoughts on practical applications of the benchmarks for multiphysics code validation. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Non-Contact Mass Density and Thermal Conductivity Measurements of Organic Thin Films Using Frequency–Domain Thermoreflectance

Thin-film organic materials are broadly used to study amorphous stabilization of active pharmaceuticals, control explosive detonation phenomena, and introduce insulation in novel thermal barriers. Their synthesis, however, introduces defects and thickness variations that warrant careful characterization of local thermophysical properties such as thermal conductivity and mass density. In this study, wide bandwidth (200 Hz to 20 MHz) frequency–domain thermoreflectance (FDTR) is demonstrated to simultaneously extract the thermal conductivity and mass density of 1 μm physical vapor-deposited indomethacin films on Si and SiO 2 substrates, as well as 10 and 100 μm films on Si. By assuming a bulk specific heat capacity, mass densities are determined with FDTR measurements of volumetric heat capacity and are in good agreement with the literature, as well as models based upon a dependence on porosity and the kinetic theory for phonons. Lastly, it is found that for broad-band FDTR measurements, insulating substrates provide improved fidelity for the extraction of thermal conductivity and volumetric heat capacity in organic thin films. Overall, this work demonstrates the potential for FDTR as a non-contact method to determine microscale mass density variations across the surface and thickness of organic thin films.

36 MATERIALS SCIENCE↗
NASA NTRS1987

The chromospheres and coronae of five G-K main-sequence stars

Previous data and high and low dispersion IUE observations of Chi 1 Ori, Alpha Cen A, Xi Boo A, Alpha Cen B, and Epsilon Eri are used to model the atmospheres of these stars from the high chromosphere to the corona. The energy lost by radiation and transferred by thermal conduction is investigated to establish the heating requirements, and results similar to those found for the solar atmosphere are found. The observed emission measurement distribution below 100,000 K can be described by a model in which Alfven wave energy input, observed through nonthermal line broadening, is balanced by radiation losses. Comparison of the coronal, transition region, and chromospheric pressures shows scaling relations which are compatible with previous flux correlations.

Jordan, C.↗
DOE OSTI2021

An efficient numerical model for predicting residual stress and strain in parts manufactured by laser powder bed fusion

Abstract Computational modeling of additively manufactured structures plays an increasingly important role in product design and optimization. For laser powder bed fusion processes, the accurate modeling of stress and distortion requires large amount of computational cost due to very localized heat input and evolving complex geometries. The current study takes advantage of a graphics processing unit accelerated explicit finite element analysis code and approximated heat conduction analysis to predict the macroscopic thermo-mechanical behavior in laser selective melting. Adjacent layers and tracks were lumped to reduce the number of time steps and elements in the finite element model. The effects of track and layer grouping on prediction accuracy and solution efficiency are investigated to provide a guidance for a cost-effective simulation. Thin-wall builds from Inconel alloy 625 (IN625) powders were simulated by applying the developed modeling approach to get the detailed residual stress and distortion at a computational speed 50 times higher than conventional approach. Under repeated heating and cooling cycles, a high tensile stress was produced near surfaces of a build due to a larger shrinkage on surface than that in central area. It is also shown that horizontal stresses concentrate near the root and top layers of the IN625 build. The predicted residual elastic strain distribution was validated by the experimental measurement using x-ray synchrotron diffraction.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2012

NOR-USA Scientific Traverse of East Antarctica: Science and Logistics on a Three-Month Expedition Across Antarctica's Farthest Frontier

Dr. Albert's current research is centered on transfer processes in porous media, including air-snow exchange in the Polar Regions and in soils in temperate areas. Her research includes field measurements, laboratory experiments, and theoretical modeling. Mary conducts field and laboratory measurements of the physical properties of natural terrain surfaces, including permeability, microstructure, and thermal conductivity. Mary uses the measurements to examine the processes of diffusion and advection of heat, mass, and chemical transport through snow and other porous media. She has developed numerical models for investigation of a variety of problems, from interstitial transport to freezing of flowing liquids. These models include a two-dimensional finite element code for air flow with heat, water vapor, and chemical transport in porous media, several multidimensional codes for diffusive transfer, as well as a computational fluid dynamics code for analysis of turbulent water flow in moving-boundary phase change problems.

Albert, Mary R.↗
NASA NTRS1981

Liquid spray cooling of a heated surface

The lowest surface temperature possible for the existance of spray evaporative cooling is determined experimentally to be a linear function of the impinging spray mass flux. A conduction-controlled analytical model of droplet evaporation gives fairly good agreement with experimental measurements at atmospheric pressure. At reduced pressures droplet evaporation rates are decreased significantly such that an optimum operating pressure exists for each desired surface heat flux. The initiation of the 'Leidenfrost state' provides the upper surface temperature bound for spray evaporative cooling.

Grissom, W. M.↗
DOE OSTI2022

Low-temperature nanoscale heat transport in a gadolinium iron garnet heterostructure probed by ultrafast x-ray diffraction

Time-resolved x-ray diffraction has been used to measure the low-temperature thermal transport properties of a Pt/Gd3Fe5O12//Gd3Ga5O12 metal/oxide heterostructure relevant to applications in spin caloritronics. A pulsed femtosecond optical signal produces a rapid temperature rise in the Pt layer, followed by heat transport into the Gd3Fe5O12 (GdIG) thin film and the Gd3Ga5O12 (GGG) substrate. The time dependence of x-ray diffraction from the GdIG layer was tracked using an accelerator-based femtosecond x-ray source. The ultrafast diffraction measurements probed the intensity of the GdIG (1 −1 2) x-ray reflection in a grazing-incidence x-ray diffraction geometry. The comparison of the variation of the diffracted x-ray intensity with a model including heat transport and the temperature dependence of the GdIG lattice parameter allows the thermal conductance of the Pt/GdIG and GdIG//GGG interfaces to be determined. Complementary synchrotron x-ray diffraction studies of the low-temperature thermal expansion properties of the GdIG layer provide a precise calibration of the temperature dependence of the GdIG lattice parameter. The interfacial thermal conductance of the Pt/GdIG and GdIG//GGG interfaces determined from the time-resolved diffraction study is of the same order of magnitude as previous reports for metal/oxide and epitaxial dielectric interfaces. The thermal parameters of the Pt/GdIG//GGG heterostructure will aid in the design and implementation of thermal transport devices and nanostructures.

36 MATERIALS SCIENCE↗
NASA NTRS1979

Correlation of predicted and measured thermal stresses on an advanced aircraft structure with similar materials

A laboratory heating test simulating hypersonic heating was conducted on a heat-sink type structure to provide basic thermal stress measurements. Six NASTRAN models utilizing various combinations of bar, shear panel, membrane, and plate elements were used to develop calculated thermal stresses. Thermal stresses were also calculated using a beam model. For a given temperature distribution there was very little variation in NASTRAN calculated thermal stresses when element types were interchanged for a given grid system. Thermal stresses calculated for the beam model compared similarly to the values obtained for the NASTRAN models. Calculated thermal stresses compared generally well to laboratory measured thermal stresses. A discrepancy of signifiance occurred between the measured and predicted thermal stresses in the skin areas. A minor anomaly in the laboratory skin heating uniformity resulted in inadequate temperature input data for the structural models.

Jenkins, J. M.↗