Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thermal diffusion”

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,549 records · Page 86

Finite-element solution for a combined radiative-conductive analysis with mixed diffuse-specular surface characteristics

The finite-element method has been applied to solve a combined radiative-conductive heat transfer problem for a large space telescope similar to those used in orbiting satellites. The derivation of the underlying matrices and associated solution algorithm for a 2-D triangular element is presented in detail. The resulting expressions for this triangular element typify such an analysis, which yields constitutive matrices when the heat equation is cast in the matrix form. The relevant matrices include those pertaining to thermal conductance, internal heat generation, radiative exchanges, and all possible external thermal loadings. Emphasis is placed on the treatment of non-linear radiative interchange between surfaces in an enclosure having mixed diffuse-specular surface characteristics. Essential differences in governing equations describing these distinctive surface characteristics are identified. Concluding remarks are drawn from an example simulating a Cassegrainian space telescope.

Lee, H.-P.↗

Three-dimensional microstructural characterization of FBR MOX fuel and the contribution of microstructural features to the thermal conductivity of the fuel

Combination of microstructural characterization, property measurements, and phase field modeling is used to investigate fast breeder reactor (FBR) mixed oxide (MOX) fuel irradiated to burnup of 13.7% fissions per initial metal atom (FIMA). Here, the fuel was characterized at different radial locations, which revealed that grey phase can be present in the central region if it nucleates on five metal precipitates (FMPs). In addition, in the mid-radial region FMPs do not diffuse out of the region once formed and the size of Pd–Te precipitates is dictated by the porosity present in the region. Thermal conductivity measurements were conducted as a function of radial location and the microstructure of the fuel was correlated with the observed trend. Reconstructions of the 3D solid and gaseous fission product structures in different regions of the fuel were used to simulate the effective thermal conductivity (ETC) of the respective regions and determine which microstructural feature has the strongest impact on thermal conductivity. Based on conducted assessment, FMPs and Pd-Te precipitates improve local conductivity of central and mid-radial regions even in the presence of grey phase, but defects are primary contributor to the degradation of thermal conductivity on the periphery of the fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Dynamical Evolution of A Tubular Leonid Persistent Train

The dynamical evolution of the persistent train of a bright Leonid meteor was examined for evidence of the source of the luminosity and the physical conditions in the meteor path. The train consisted of two parallel somewhat diffuse luminous tracks, interpreted as the walls of a tube. A general lack of wind shear along the trail allowed these structures to remain intact for nearly 200 s, from which it was possible to determine that the tubular structure expanded at a near constant 10.5 m/s, independent of altitude between 86 and 97 km. An initial fast decrease of train intensity below 90 km was followed by an increase in intensity and then a gradual decrease at longer times, whereas at high attitudes the integrated intensity was nearly constant with time. These results are compared to a model that describes the dynamical evolution of the train by diffusion, following an initial rapid expansion of the hot gaseous trail behind the meteoroid. The train luminosity is produced by O ((sup 1)S) emission at 557 nm, driven by elevated atomic O levels produced by the meteor impact, as well as chemiluminescent reactions of the ablated metals Na and Fe with O3. Ozone is rapidly removed within the train, both by thermal decomposition and catalytic destruction by the metallic species. Hence, the brightest emission occurs at the edge of the train between outwardly diffusing metallic species and inwardly diffusing O3. Although the model is able to account plausibly for a number of characteristic features of the train evolution, significant discrepancies remain that cannot easily be resolved.

Jenniskens, Peter↗

Vacuum plasma coatings for turbine blades

Turbine blades, vacuum plasma spray coated with NiCrAlY, CoCrAlY or NiCrAlY/Cr2O3, were evaluated and rated superior to standard space shuttle main engine (SSME) coated blades. Ratings were based primarily on 25 thermal cycles in the MSFC Burner Rig Tester, cycling between 1700 F (gaseous H2) and -423 F (liquid H2). These tests showed no spalling on blades with improved vacuum plasma coatings, while standard blades spalled. Thermal barrier coatings of ZrO2, while superior to standard coatings, lacked the overall performance desired. Fatigue and tensile specimens, machined from MAR-M-246(Hf) test bars identical to the blades were vacuum plasma spray coated, diffusion bond treated, and tested to qualify the vacuum plasma spray process for flight hardware testing and application. While NiCrAlY/Cr2O3 offers significant improvement over standard coatings in durability and thermal protection, studies continue with an objective to develop coatings offering even greater improvements.

Holmes, R. R.↗

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Growth of Solid Solution Single Crystals

The solidification of a solid solution semiconductor, having a wide separation between liquidus and serious has been extensively studied in ground based, high magnetic field and Spacelab experiments. Two alloys of mercury cadmium telluride have been studied; mercury cadmium telluride with 80.0 mole percent of HgTe and 84.8 mole percent respectively. These alloys are extremely difficult to grow by directional solidification on earth due to high solutal and thermal density differences that give rise to fluid flow and consequent loss of interface shape and composition. Diffusion controlled growth is therefore impossible to achieve in conventional directional solidification. The ground based experiments consisted of growing crystals in several different configurations of heat pipe furnaces, NASA's Advanced Automated Directional Solidification Furnace (AADSF), and a similar furnace incorporated in a superconducting magnet capable of operating at up to 5T. The first microgravity experiment took place during the flight of STS-62 in March 1994, with the AADSF installed on the second United States Microgravity Payload (USMP-2). The alloy was solidified at 3/4 inch per day over a 9 day period, and for the first time a detailed evaluation was performed of residual acceleration effects. The second flight experiment took place in the fourth United States Microgravity Payload Mission (USMP-4) in November 1997. Due to contamination of the furnace system by a previously processed sample, the sample was not received until May 1998, and the preliminary analysis shows that the conditions prevailing during the experiment were quite different from the requirements requested prior to the mission. Early results are indicating that the sample may not accomplish the desired objectives. As with the USMP-2 mission, the results of the ground based experiments were compared with the crystal grown in orbit under microgravity conditions. On the earth, it has been demonstrated that the application of the magnetic field leads to a significant reduction in fluid flow, with improved homogeneity of composition. The field strength required to suppress flow increases with diameter of the material. The 8 mm diameter sample used here was less than the upper diameter limit for a ST magnet. The configuration for USMP-4 was changed so that the material was seeded and other processing techniques were also modified. It was decided to examine the effects of a strong magnetic field under the modified configuration and parameters. A further change from USMP-2 was that a different composition of material was grown, namely with 0.152 mole fraction of cadmium telluride rather than the 0.200 of the USMP-2 experiment. The objective was to grow highly homogeneous, low defect density material of a composition at which the conduction band and the valence band of the material impinge against each other. As indicated, the furnace was contaminated during the mission. As a result of solid debris remaining in the furnace bore, the cartridge in this experiment, denoted as SL1-417, was significantly bent during the insertion phase. During translation the cartridge scraped against the plate which isolates the hot and cold zones of the furnace. Thermocouples indicated that a thermal assymetry resulted. The scraping in the slow translation or crystal growth part of the processing was not smooth and it is probable that the jitter was sufficient to give rise to convection in the melt. Early measurements of composition from the surface of the sample have shown that the composition varies in an oscillatory manner.

Lehoczky, Sandor L.↗

The role of Si in determining the stability of the θ; precipitate in Al-Cu-Mn-Zr alloys

Precipitate-strengthened Al-Cu-Mn-Zr (ACMZ) alloys demonstrate improved microstructural stability compared to conventional Al-Cu alloys. Si is typically considered an impurity in these alloys, but a systematic analysis of ACMZ alloys with a range of Si levels revealed that there is a Si content range (~0.05 wt. % < Si < ~0.10 wt. %) for achieving optimal high temperature microstructural stability. A combination of hardness testing, scanning electron microscopy, and atom probe tomography has been used to understand the relationship between Si content, aging response, and thermal stability in ACMZ alloys. It is shown that the optimum Si content range corresponds to a reduced as-aged hardness, but a greater hardness value retained by the alloy after thermal exposure to 350 °C for 200 hours. A mechanism is proposed in which optimum Si levels reduce the number density of nucleation sites for θ' precipitates, resulting in larger precipitates that, on average, provide a reduced as-aged strength but are more coarsening resistant. This inherent θ' coarsening resistance allows more time for slow-diffusing Mn and Zr to reach the θ' interface and further stabilize the precipitates. Designing the optimal as-aged microstructure for improved thermal stability of the alloy by controlling impurity levels is a concept that has implications beyond the investigated Al-Cu alloy system.

36 MATERIALS SCIENCE↗

Multi-Dimensional, Non-Pyrolyzing Ablation Test Problems

Non-pyrolyzingcarbonaceous materials represent a class of candidate material for hypersonic vehicle components providing both structural and thermal protection system capabilities. Two problems relevant to this technology are presented. The first considers the one-dimensional ablation of a carbon material subject to convective heating. The second considers two-dimensional conduction in a rectangular block subject to radiative heating. Surface thermochemistry for both problems includes finite-rate surface kinetics at low temperatures, diffusion limited ablation at intermediate temperatures, and vaporization at high temperatures. The first problem requires the solution of both the steady-state thermal profile with respect to the ablating surface and the transient thermal history for a one-dimensional ablating planar slab with temperature-dependent material properties. The slab front face is convectively heated and also reradiates to a room temperature environment. The back face is adiabatic. The steady-state temperature profile and steady-state mass loss rate should be predicted. Time-dependent front and back face temperature, surface recession and recession rate along with the final temperature profile should be predicted for the time-dependent solution. The second problem requires the solution for the transient temperature history for an ablating, two-dimensional rectangular solid with anisotropic, temperature-dependent thermal properties. The front face is radiatively heated, convectively cooled, and also reradiates to a room temperature environment. The back face and sidewalls are adiabatic. The solution should include the following 9 items: final surface recession profile, time-dependent temperature history of both the front face and back face at both the centerline and sidewall, as well as the time-dependent surface recession and recession rate on the front face at both the centerline and sidewall. The results of the problems from all submitters will be collected, summarized, and presented at a later conference.

Ablation↗

M558 radioactive tracer diffusion

This experiment was performed in Skylab 3 with two objectives in mind. First, the experimental self-diffusion coefficients for liquid zinc were to be determined in a convection-free environment. Secondly the reduction in convective mixing in earth gravity by going into the zero-gravity environment of space was to be estimated. The experiment was designed to utilize high temperatures and linear thermal gradients provided by the M518 Multipurpose Electric Furnace, and the radioactivity of zinc-65 of 245-day half-life to investigate self-diffusion in liquid zinc. The distribution of zinc-65 tracer, after melting, maintaining at soak temperature for 1 hour of soak time and then resolidifying, was obtained by sample sectioning. The concentration of activity of each section (microcurie-gram) was plotted against positions along the sample axial and radial position. Experimental data and theoretical results from solution of Fick's law of diffusion in one dimensional were compared. Samples tested on earth showed very rapid diffusion. Diffusion coefficient in unit gravity was 50 times the zero-gravity diffusion coefficient of Skylab.

Ukanwa, A. O.↗

Mass and energy balance of the cold Io torus

A new model of the cold Io torus is described. Ions and energy are injected into the system by independent processes so that the mass balance is isolated from the energy balance. The primary source of energy is local ionization and acceleration of hot pickup ions resulting from charge exchange between thermal ions and an extended cloud of Iogenic sulfur and oxygen atoms. The primary energy loss mechanism of the plasma is collisionally excited line emission at optical wavelengths. The primary ion source is radial diffusion inward from the hot torus on a time scale of 140-710 days. The primary ion loss mechanism is a novel two-step enhanced recombination mechanism involving charge exchange between thermal ions and an extended cloud of neutral SO2 molecules, followed by rapid dissociative recombination of the resultant molecular ion. The model provides a self-consistent solution which reconciles a number of diverse observations with known physical processes.

Moreno, M. A.↗

Observations of the diffuse interstellar EUV radiation field

Results are reported for observations and a preliminary theoretical analysis of the diffuse interstellar EUV radiation field exhibiting positive measurable fluxes in the wavelength range from 912 to 1150 A. The observations were made with the UV spectrometers on the two Voyager spacecraft, and a spectral half-width of 30 A was obtained. The procedures used in the analysis are described in detail, particularly those relating to the elimination of thermal and cosmic-ray noise, internal instrumental scattering from bright lines, and possible contributions from direct starlight. Measurements in several directions reveal diffuse intensities at 975 A ranging from a maximum of about 10 to the -7th erg/sq cm per (sec sr A) near the galactic plane to an upper limit of 10 to the -8th erg/sq cm per (sec sr A). A comparison of the present measurements with various theoretical predictions and previous measurements indicates reasonable agreement with some theoretical predictions of the interstellar EUV intensity.

Sandel, B. R.↗

An Investigation of Thermal Properties of 2D Materials [Dissertation]

Studying the thermal conductivity of 2D materials is important due to the applications of 2D materials in fields such as thermal management, thermoelectricity, renewable energy, and sensors. As such, measurements of the thermal conductivity of these 2D materials become important to measure. Thermal conductivity is often difficult to measure for 2D materials due to their atomically thin nature and many experimental methods for doing so requiring contact with the sample, which can alter the thermal properties. A non-contact method for calculating the thermal conductivity of 2D materials supported on substrates in order to model the thermal conductivity of 2D materials for devices, is proposed and experimentally performed in this dissertation. The optothermal Raman technique is a useful non-contact diagnostic technique useful in determining the thermal conductivity of 2D materials. The optothermal Raman typically does not account for heat losses due to convection or radiation or substrate resistance, which are shown to be important factors to consider when developing an optothermal Raman model. Additionally, the calculation of the interfacial thermal conductance between the bottom surface of the sample and the top surface of the substrate, plays an important role in determining the final value of the thermal conductivity of a supported sample, and will yield differing results based on whether or not the conductance is calculated using an approach such as the Diffuse Mismatch Model (DMM) or calculated directly by varying the laser heating profile (usually done by changing the laser objective). This is shown to be the case for both graphene on Ni, graphene on Cu, and SnSe 2 on Cu. In addition to experimentally calculating the thermal conductivity of a 2D material with the optothermal Raman technique, the thermal conductivity of 2D materials can also be calculated using computational methods. The three-phonon method is a method which can be used to simulate phonon scattering processes and determine the thermal conductivity of semiconductors, wherein phonon scattering is the dominant mechanism which determines the thermal conductivity. The three-phonon method uses relaxation times for phonon scattering with other phonons, electrons, and other material system elements, such as isotopes or material defects, in order to create a single-mode relaxation time approximation (SMRTA), which is used to calculate the final value of the thermal conductivity. An important consideration when determining the thermal conductivity of a 2D material using this method is the device geometry, which is reflected in this work as the phonon-boundary scattering relaxation time. This inclusion is important along with the inclusion of phonon-electron scattering in accurately determining the thermal conductivity of a 2D material. In both the optothermal Raman experiments and the three-phonon method computations, strain is shown to have a demonstrable effect on the thermal conductivity of 2D materials. When a 1.1% strain was applied to the mechanical properties of SnSe, the three-phonon processes yielded a lower thermal conductivity than the no-strain case. For the optothermal Raman experiments, the strain induced in the Cu substrate and transferred to a single-layer graphene (SLG) sample yields a trend where the thermal conductivity of the SLG decreases with respect to strain applied. In the case where the interfacial thermal conductance was calculated directly, the conductance increased with respect to strain applied. This presents strain as a reliable and viable method for tuning the thermal properties of 2D materials for device applications.

36 MATERIALS SCIENCE↗

Numerical thermalization in 2D PIC simulations: Practical estimates for low-temperature plasma simulations

The process of numerical thermalization in particle-in-cell (PIC) simulations has been studied extensively. It is analogous to Coulomb collisions in real plasmas, causing particle velocity distributions (VDFs) to evolve toward a Maxwellian as macroparticles experience polarization drag and resonantly interact with the fluctuation spectrum. This paper presents a practical tutorial on the effects of numerical thermalization in 2D PIC applications. Scenarios of interest include simulations, which must be run for many thousands of plasma periods and contain a population of cold electrons that leave the simulation space very slowly. This is particularly relevant to many low-temperature plasma discharges and materials processing applications. We present numerical drag and diffusion coefficients and their associated timescales for a variety of grid resolutions, discussing the circumstances under which the electron VDF is modified by numerical thermalization. Though the effects described here have been known for many decades, direct comparison of analytically derived, velocity-dependent numerical relaxation timescales to those of other relevant processes has not often been applied in practice due to complications that arise in calculating thermalization rates in 1D simulations. Using these comparisons, we estimate the impact of numerical thermalization in several examples of low-temperature plasma applications including capacitively coupled plasma discharges, inductively coupled plasma discharges, beam plasmas, and hollow cathode discharges. Finally, we discuss possible strategies for mitigating numerical relaxation effects in 2D PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radio continuum emission from T Tauri stars and associated structures

The current status of radio continuum observations of low mass, pre-main sequence stars (T Tauri stars) is discussed. The majority of T Tauri radio sources have been detected by the VLA telescope and a guide to both these and earlier observations is presented. Both nonthermal and thermal radio continuum emission has been detected from T Tauri stars. The likely emission processes and specific examples of their occurrence are outlined. The detected radio sources are often associated with the surroundings of the star, such as bipolar flows or diffuse extended optical nebulosity, rather than regions very close to the stellar surface. Only recently have radio continuum observations reached a sensitivity appropriate to the expected level of thermal free-free emission from the winds of ordinary T Tauri stars and it is possible that many more radio sources will now be found associated with these stars.

Brown, A.↗

Radial diffusion of low-energy plasma ions in Saturn's magnetosphere

Radial diffusion of low-energy plasma ions in Saturn's magnetosphere is investigated using a comprehensive set of equations for radial diffusion that incorporate distributed sources and sinks of ions. The results of calculations indicate that the radial-diffusion transport of low-energy O(+) ions with a source in the neutral H2O cloud of the satellites Dione and Tethys can account for Voyager observations of thermal heavy ions in Saturn's magnetosphere. The source rate was calculated to be about 10 to the 26th O(+) ions/sec, in good agreement with the sputtering calculations of Johnson et al. (1989). It is estimated that, due to fast radial diffusion, the residence time of O(+) ions in the Dione-Tethys torus is about 30 days, sufficiently short to account for the plasma density observed there. The densities of hot H(+) and N(+) resulting from the ionization and pickup of Titan's neutral clouds in the outer magnetosphere can also be accounted for within the framework of diffusive ion transport.

Barbosa, D. D.↗

Molecular simulation of the separation of toluene and p-xylene with the thermally-robust ionic liquid triphenyl-p-phenyl sulfonyl phenyl phosphonium

In this work, we use molecular dynamics simulations to evaluate the behavior of toluene (TOL) and p-xylene (XYL) within three different solvents: n -hexane (HEX), n -heptane (HEP), and a thermally-robust ionic liquid (triphenyl-p-phenyl sulfonyl phenyl phosphonium + bis(trifluoromethylsulfonyl)imide (TPSP + Tf 2 N)). Several different temperatures are explored, while the solvation structures, energetics, and dynamics of TOL and XYL are evaluated. The solvation free energy (ΔG solv ) and transfer free energy (ΔG transfer ) of TOL and XYL are calculated using thermodynamic integration (TI). The transfer free energy of XYL from HEP/HEX to TPSP + Tf 2 N is more favorable than that of TOL over the entire temperature range, while the transfer of both TOL and XYL becomes less favorable as the temperature increases. The diffusion rates of TOL and XYL are lower in the ionic liquid when compared to HEP and HEX, but the thermal stability of TPSP + Tf 2 N may allow for higher operating temperatures and accelerated mass transport rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cosmic-ray heating of cooling flows - A critical analysis

It is shown that a combination of MHD wave-mediated cosmic ray heating and thermal conduction could balance cooling in intracluster media and substantially reduce the rate of inflow. The appropriate system of steady state equations is solved, including a new self-consistent formulation for the cosmic-ray diffusivity. Models which can produce substantial positive temperature gradients in static configurations are found when conduction is reduced by a factor of 10 or more. These models have too-flat thermal pressure profiles compared with observations. It is found that cosmic-ray heating is unlikely either to stabilize positive density perturbations against condensation or to contribute appreciably to the powering of the optical filaments.

Loewenstein, Michael↗

Higher-Order Finite Elements for Computing Thermal Radiation

Two variants of the finite-element method have been developed for use in computational simulations of radiative transfers of heat among diffuse gray surfaces. Both variants involve the use of higher-order finite elements, across which temperatures and radiative quantities are assumed to vary according to certain approximations. In this and other applications, higher-order finite elements are used to increase (relative to classical finite elements, which are assumed to be isothermal) the accuracies of final numerical results without having to refine computational meshes excessively and thereby incur excessive computation times. One of the variants is termed the radiation sub-element (RSE) method, which, itself, is subject to a number of variations. This is the simplest and most straightforward approach to representation of spatially variable surface radiation. Any computer code that, heretofore, could model surface-to-surface radiation can incorporate the RSE method without major modifications. In the basic form of the RSE method, each finite element selected for use in computing radiative heat transfer is considered to be a parent element and is divided into sub-elements for the purpose of solving the surface-to-surface radiation-exchange problem. The sub-elements are then treated as classical finite elements; that is, they are assumed to be isothermal, and their view factors and absorbed heat fluxes are calculated accordingly. The heat fluxes absorbed by the sub-elements are then transferred back to the parent element to obtain a radiative heat flux that varies spatially across the parent element. Variants of the RSE method involve the use of polynomials to interpolate and/or extrapolate to approximate spatial variations of physical quantities. The other variant of the finite-element method is termed the integration method (IM). Unlike in the RSE methods, the parent finite elements are not subdivided into smaller elements, and neither isothermality nor other unrealistic physical conditions are assumed. Instead, the equations of radiative heat transfer are integrated numerically over the parent finite elements by use of a computationally efficient Gaussian integration scheme.

Gould, Dana C.↗