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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Predicting Thermal Response in a Li-ion Cell on a UAV Fight Profile

As the energy storage devices continue to "pack" more energy in a small space, any damage, battery component failure, manufacturing defect, or electrically abusing the battery can lead to catastrophic thermal runaway events. A catastrophic thermal event in a cell leads to high temperature, in some instances spewing of battery materials due to gas development from side reactions initiated due to high internal temperatures. Also, a thermal runaway event can propagate from a single "failed" cell to the pack in a battery pack, leading to a more significant event. Mitigating a thermal runaway event is important in the commercial and automotive sectors. However, preventing such events in an electric aircraft (or air taxis) is paramount due to the lack of alternatives in the event of a failure. Battery prognostics algorithms allow the prediction of state-of-charge (SOC) and end-of-life (EOL) of a Li-ion battery in a UAV (unmanned air vehicle) [1]. For this presentation, we will extend this two-level battery predictive algorithm to predict SOC, EOL, and estimated maximum temperature during a simulated flight. The model is extended by integrating a lumped physics-driven thermal model for high current densities [2]. The parameters used to control SOC and EOL are maximum storable charge, time constant for Li-ion diffusivity in the carbon particles, and internal cell resistance. Cycling leads to an increase in the heat generated by an aged Li-ion cell with a LiyCoO2 (LCO) cathode and a LixC6 (MCMB) anode. The aging of a cell leads to increase in SEI layer thickness, the diffusion time for the lithium ions inside the electrodes, and the local reaction rates, in addition to the thermodynamic abuse caused by fixed cycling voltages controlled by a Battery Management System. As the battery ages, the cell resistance increases, while the onset temperature of the thermal runaway decreases (depends on the cell chemistry and cell abuse history). Any large deviation of the cell temperature from the estimated (expected) value can identify a faulty cell. Since SEI decomposition has the lowest onset temperature in the series of reactions leading to thermal runaway, the model considers the self-heating rate of the SEI decomposition as onset temperature (similar to Ref. [3]). The parameters in the Arrhenius equation for the SEI heating rate depend on the number of cycles, the cell's operating temperature, and the cell's abuse history [4,5]. Coupling the electrochemical, thermal, and aging model allow the prognostic algorithm to estimate a typical cell voltage and temperature as a function of age (cycling and calendar), whose departure from measured values from the BMS is used to identify a safety event. In addition, we will present the results from two simulated flight scenarios for a UAV: typical and extreme, since the power requirements vary significantly during take-off, landing, and changing altitudes, while the power requirements remain low during the cruise. For this presentation, the power requirement for a battery pack in a UAV is scaled to a single cell. This cell is cycled through a simulated profile, and the data is collected and used to predict a safety event.

Li-ion↗
NASA NTRS1981

Equatorial superrotation in a thermally driven zonally symmetric circulation

Near the equator where the Coriolis force vanishes, the momentum balance for the axially symmetric circulation is established between horizontal and vertical diffusion, which, a priori, does not impose constraints on the direction or magnitude of the zonal winds. Solar radiation absorbed at low latitudes is a major force in driving large scale motions with air rising near the equator and falling at higher latitudes. In the upper leg of the meridional cell, angular momentum is redistributed so that the atmosphere tends to subrotate (or corotate) at low latitudes and superrotate at high latitudes. In the lower leg, however, the process is reversed and produces a tendency for the equatorial region to superrotate. The outcome depends on the energy budget which is closely coupled to the momentum budget through the thermal wind equation; a pressure (temperature) maximum is required to sustain equatorial superrotation. Such a condition arises in regions which are convectively unstable and the temperature lapse rate is superadiabatic. It should arise in the tropospheres of Jupiter and Saturn; planetary energy from the interior is carried to higher altitudes where radiation to space becomes important. Upward equatorial motions in the direct and indirect circulations (Ferrel-Thomson type) imposed by insolation can then trap dynamic energy for equatorial heating which can sustain the superrotation of the equatorial region.

Mayr, H. G.↗
NASA NTRS1975

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.↗
NASA NTRS2000

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↗
NASA NTRS1985

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.↗
NASA NTRS1999

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.↗
NASA NTRS2016

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↗
NASA NTRS1974

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.↗
NASA NTRS1986

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.↗
NASA NTRS1979

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.↗
NASA NTRS1985

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.↗
NASA NTRS1990

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.↗
NASA NTRS1991

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↗
NASA NTRS2004

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.↗
NASA NTRSDate not supplied

Updates on the Predictive Materials Modeling Software Tools

Updates on NASA‘s efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. The PMM effort is part of the Entry Systems Modeling (ESM) project under NASA’s Game Changing Development (GCD) program. To reduce the need for extensive testing and accelerate the design cycle process, ESM is developing simulation and modeling tools that enable the characterization of the properties of thermal protection materials and their response to extremely hot plasma. The Porous Microstructure Analysis (PuMA) software has been developed to compute effective material properties and perform material response simulations on digitized microstructures of porous media. PuMA is able to import three-dimensional digital images obtained from X-ray microtomography or to generate artificial microstructures that mimic real materials. PuMA also provides a module for interactive 3D visualizations. Version 3, which was recently released as open-source, includes modules to compute simple morphological properties such as porosity, volume fractions, pore diameter, and specific surface area. Additional capabilities include the determination of effective thermal and electrical conductivity (both radiative and solid conduction - including the ability to simulate local anisotropy for the latter); effective diffusivity and tortuosity from the continuum to the rarefied regime; techniques to determine the local material orientation, as well as mechanical properties (elasticity coefficients), and permeability. Computed properties are then used to inform a macro-scale material response model, such as those implemented in the Porous material Analysis Toolbox based on OpenFOAM (PATO) software developed within ESM. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale, assuming Local Thermal Equilibrium. Recent efforts include the development of a mechanical erosion model as well as a unified model allowing an intrinsic coupling between fluid and material. Comparison to flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI] and Mars 2020 MEDLI2) is critical in order to validate these computational tools. Examples of ablative material response using the code will be presented, including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrated the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains through the use of massively parallel computations.

material modeling↗
NASA NTRS1991

Titanium honeycomb panel testing

The paper describes the procedures of thermal mechanical tests carried out at the NASA Dryden Flight Research Facility on two tianium honeycomb wing panels bonded using liquid interface diffusion (LID) technique, and presents the results of these tests. The 58.4 cm square panels consisted of two 0.152-cm-thick Ti 6-2-4-2 face sheets LID-bonded to a 1.9-cm-thick honeycomb core, with bearing plates fastened to the perimeter of the upper and the lower panel surfaces. The panels were instrumented with sensors for measuring surface temperature, strain, and deflections to 315 C and 482 C. Thermal stress levels representative of those encountered during aerodynamic heating were produced by heating the upper panel surface and restraining all four edges. After more than 100 thermal cycles from room temperature to 315 C and 50 cycles from room temperature to 482 C, no significant structural degradation was detected in the panels.

Richards, W. L.↗
NASA NTRS1979

Ar-40/Ar-39 age of the Shergotty achondrite and implications for its post-shock thermal history

Ar-40/Ar-39 measurements are used to determine the age of the Shergotty achondrite and the chronology of the shock event responsible for the complete conversion of its plagioclase to maskelynite is discussed. Apparent ages are found to vary between 240 and 640 million years for the whole rock sample, with a plateau age of 254 million years for a maskelynite separate. The Rb-Sr age of 165 million years determined by Nyquist at al (1978) suggests that the maskelynite as well as the whole rock was incompletely degassed. Argon diffusion characteristics indicate a post-shock cooling time greater than 1000 years and a burial depth greater than 300 m for a thermal model of a cooling ejecta blanket of variable thickness. It is concluded that the shock event which degassed the argon and reset the Rb-Sr systematics occurred between 165 and 250 million years ago when the parent body experienced a collision in the asteroid belt.

Bogard, D. D.↗