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,495 records · Page 83

Radiation from an array of gray circular fins of trapezoidal profile

A numerical method is developed to calculate the temperature distribution and radiation heat transfer for an annular fin and tube radiator, with fins having trapezoidal profiles. All surfaces are assumed gray and to emit and reflect diffusely. Radiative interactions between adjacent fins and between the fins and tube are included. The thermal conductivity of the fin material may vary linearly with temperature. Results of a parametric study of the special case of circular fins of triangular profile having constant thermal conductivity are presented and used to optimize a fin array with respect to minimum weight.

Schnurr, N. M.↗

A frost formation model and its validation under various experimental conditions

A numerical model that was used to calculate the frost properties for all regimes of frost growth is described. In the first regime of frost growth, the initial frost density and thickness was modeled from the theories of crystal growth. The 'frost point' temperature was modeled as a linear interpolation between the dew point temperature and the fog point temperature, based upon the nucleating capability of the particular condensing surfaces. For a second regime of frost growth, the diffusion model was adopted with the following enhancements: the generalized correlation of the water frost thermal conductivity was applied to practically all water frost layers being careful to ensure that the calculated heat and mass transfer coefficients agreed with experimental measurements of the same coefficients.

Dietenberger, M. A.↗

Radiation Heat Transfer Between Diffuse-Gray Surfaces Using Higher Order Finite Elements

This paper presents recent work on developing methods for analyzing radiation heat transfer between diffuse-gray surfaces using p-version finite elements. The work was motivated by a thermal analysis of a High Speed Civil Transport (HSCT) wing structure which showed the importance of radiation heat transfer throughout the structure. The analysis also showed that refining the finite element mesh to accurately capture the temperature distribution on the internal structure led to very large meshes with unacceptably long execution times. Traditional methods for calculating surface-to-surface radiation are based on assumptions that are not appropriate for p-version finite elements. Two methods for determining internal radiation heat transfer are developed for one and two-dimensional p-version finite elements. In the first method, higher-order elements are divided into a number of sub-elements. Traditional methods are used to determine radiation heat flux along each sub-element and then mapped back to the parent element. In the second method, the radiation heat transfer equations are numerically integrated over the higher-order element. Comparisons with analytical solutions show that the integration scheme is generally more accurate than the sub-element method. Comparison to results from traditional finite elements shows that significant reduction in the number of elements in the mesh is possible using higher-order (p-version) finite elements.

Gould, Dana C.↗

Engineering Analysis Studies for Preliminary Design of Lightweight Cryogenic Hydrogen Tanks in UAV Applications

A series of engineering analysis studies were conducted to investigate the potential application of nanoclay-enhanced graphite/epoxy composites and polymer cross-linked silica aerogels in cryogenic hydrogen storage tank designs. This assessment focused on the application of these materials in spherical tank designs for unmanned aeronautic vehicles with mission durations of 14 days. Two cryogenic hydrogen tank design concepts were considered: a vacuum-jacketed design and a sandwiched construction with an aerogel insulating core. Analyses included thermal and structural analyses of the tank designs as well as an analysis of hydrogen diffusion to specify the material permeability requirements. The analyses also provided material property targets for the continued development of cross-linked aerogels and nanoclay-enhanced graphite/epoxy composites for cryogenic storage tank applications. The results reveal that a sandwiched construction with an aerogel core is not a viable design solution for a 14-day mission. A vacuum-jacketed design approach was shown to be far superior to an aerogel. Aerogel insulation may be feasible for shorter duration missions. The results also reveal that the application of nanoclay-enhanced graphite/epoxy should be limited to the construction of outer tanks in a vacuum-jacketed design, since a graphite/epoxy inner tank does not provide a significant weight savings over aluminum and since the ability of nanoclay-enhanced graphite/epoxy to limit hydrogen permeation is still in question.

Sullivan, Roy M.↗

The Performance of Current Atmospheric Radiation Codes in Phase I of CIRC

The Continual Intercomparison of Radiation Codes (CIRC) is intended as an evolving and regularly updated reference source for evaluation of radiative transfer (RT) codes used in Global Climate Models and other atmospheric applications. In our presentation we will discuss our evaluation of the performance of 13 shortwave and 11 longwave RT codes that participated in Phase I of CIRC. CIRC differs from previous intercomparisons in that it relies on an observationally validated catalogue of cases. The seven CIRC Phase I baseline cases, five cloud-free, and two with overcast liquid clouds, are built around observations by the Atmospheric Radiation Measurements (ARM) program that satisfy the goals .of Phase I, namely to examine RT model performance in realistic, yet not overly complex, atmospheric conditions. Besides the seven baseline cases, additional idealized "subcases" are also examined to facilitate interpretation of model errors. We will quantify individual model performance with respect to reference line-by-line calculations, and will also highlight RT code behavior for conditions of doubled CO2 , aspects of utilizing a spectral specification of surface albedo, and the impact of the inclusion of scattering in the thermal infrared. Our analysis suggests that RT codes should work towards improving their calculation of diffuse shortwave flux, shortwave absorption, treatment of spectral surface albedo, and shortwave CO2 forcing. Despite practical difficulties in comparing our results to previous results by the Intercomparison of Radiation Codes in Climate Models (ICRCCM) conducted about 20 years ago, it appears that the current generation of RT codes do indeed perform better than the codes of the ICRCCM era. By enhancing the range of conditions under which participating codes are tested, future CIRC phases will hopefully allow even more rigorous examination of RT code performance.

Oreopoulos, L.↗

Hydroxylation of Apollo 17 Soil Sample 78421 by Solar Wind Protons

Hydroxylation by solar wind protons has been simulated in our laboratory on Apollo 17 lunar sample 78421, a very mature regolith sample that is rich with agglutinates (68%). The goal of this study was to determine the rate of hydroxyl formation and their thermal stability by monitoring changes in the SiOH (hydroxyl) stretching band near 3 μm using diffuse reflectance FTIR spectroscopy (DRIFTS). A 2 keV H2+ ion beam was used to simulate proton implantation on 78421 and on a crushed fused silica sample. We find that the OH band does not change unless the samples have been annealed in vacuum prior to irradiation. Qualitatively, the OH bands for the fused silica and 78421 are very different. The OH band for fused silica is centered at 2.74 μm and is relatively sharp ranging from 2.67 - 3.1 μm at full-width-at-half-maximum (FWHM), while the OH band for 78421 is centered at 3.0 μm and ranges from 2.74 - 3.37 μm at FWHM. The increase in wavelength and broadened nature of the OH band in 78421 may be associated with the OH’s proximity to surface defects and/or lattice vacancies. The lack of the H2O bending mode at 6.1 μm indicates that any adsorbed terrestrial H2O is below our detection limit, and therefore the H2O stretching mode at 2.9 μm is not significantly contributing to the broad 3 μm OH band and implies that proton implantation by itself does not lead to water formation. To simulate the maximum dayside temperature on the lunar surface, the lunar sample was heated after proton irradiation. The proton induced OH concentration was reduced by as much as 25% after heating to 400 K (127 °C).

Apollo 17↗

J-OCTA Applications to in Space Manufacturing and Nanofluidic Technologies

We consider examples of J-OCTA software application to physics-based analysis of in-space manufacturing technologies. Firstly, we consider an atomistic model of thermal welding at the polymer-polymer interface for applications to 3D manufacturing by following diffusion of polymer chains. We show that each component of the polymer blend has its own characteristic time of diffusion and that both strain–stress and shear viscosity curves agree with experimental data. Next, nonlinear shrinkage of the metal part during manufacturing by bound metal deposition under microgravity, is considered using a multi-scale physics-based approach that spans from atomistic dynamics to full-part shrinkage. At the smallest scale we use atomistic dynamics to simulated sintering of a few Ti6Al4V nanoparticles, study sintering mechanisms and estimate important parameters including grain boundary width and diffusion rates of the atoms at the surface that are used at the other scales of modelling. Finally, we analyse selective ionic conduction through an artificial nanopore in a graphene sheet motivated by applications to desalination and energy harvesting. We build distributions of the number and orientation of water molecules in the ion’s hydration shells providing novel insight into the ion-water-carbon interaction at the nanopore.

molecular dynamics↗

Fundamental Studies of Tritium Diffusivity in Pure and Defective Zircaloy-4 Getters (Tritium Science Program FY2020 Report)

Zirconium (Zr) alloys are extensively used as tritium (T) getters in nuclear reactors due to their low absorption cross section to thermal neutron, good mechanical properties at high temperature and resistant to corrosion in different environmental conditions. Zr alloys have better thermal properties than many other refractory alloys including stainless steel. The nuclear characteristic is simulated, and T is produced in tritium-producing burnable absorber rods (TPBAR) by the irradiation of neutron flux in pressurized water reactor (PWR). T thus produced diffuses through the pellets and is captured by the getter where it chemically reacts with Zr to form metal hydride (ZrTx). These hydrides are brittle and affect adversely the mechanical properties of alloys. In addition to that, mismatch in the lattice structures after hydride formation creates a stress which needs to be considered in component design and their life evaluation. Therefore, understanding the behavior of T and its species becomes significant as fuel burnup is increased, which leads to increase in hydrogen pickup and oxide formation. The transport of T and its species in pure Zr alloys, and in alloys with certain alloying elements and defects are important to understand in order to enhance the performance of the materials.

36 MATERIALS SCIENCE↗

Nature of fluid flows in differentially heated cylindrical container filled with a stratified solution

Semiconductor crystals such as Hg(1-x)Cd(x)Te grown by unidirectional solidification Bridgmann method have shown compositional segregations in both the axial and radial directions. Due to the wide separation between the liquidus and the solidus of its pseudobinary phase diagram, there is a diffusion layer of higher HgTe content built up in the melt near the melt-solid interface which gives a solute concentration gradient in the axial direction. Because of the higher thermal conductivity in the melt than that in the crystal there is a thermal leakage through the fused silica crucible wall near the melt-solid interface. This gives a thermal gradient in the radial direction. Hart (1971), Thorpe, Hutt and Soulsby (1969) have shown that under such condition a fluid will become convectively unstable as a result of different diffusivities of temperature and solute. It is quite important to understand the effects of this thermosolute convection on the compositional segregation in the unidirectionally solidified crystals. To reach this goal, we start with a simplified problem. We study the nature of fluid flows of a stratified solution in a cylindrical container with a radial temperature gradient. The cylindrical container wall is considered to be maintained at a higher temperature than that at the center of the solution and the solution in the lower gravitational direction has higher solute concentration which decrease linearly to a lower concentration and then remain constant to the top of the solution. The sample solution is taken to be salt water.

Wang, Jai-Ching↗

Intrinsic thermal expansion and tunability of thermal expansion coefficient in Ni-substituted Co 2 V 2 O 7

Abstract Framework oxide materials are well-known for exhibiting not only negative thermal expansion (NTE), but also demonstrating thermal expansion that can be controlled using composition as a tuning parameter. In this work, we study the intrinsic thermal expansion properties of Co 2 V 2 O 7 , which has shown bulk linear NTE, and attempt to understand how substituting Ni 2+ for Co 2+ will affect the thermal expansion. The isomorphic solid solution is synthesized through solid-state methods and characterized using x-ray diffraction (XRD), diffuse reflectance spectroscopy, and neutron diffraction. The size difference between Ni 2+ and Co 2+ as well as the polyhedral volume of each Co 2+ metal coordination environment in the crystal structure allows Ni 2+ to partially be directed toward one crystallographic site over the other. Variable temperature synchrotron XRD data are employed to understand intrinsic thermal expansion. Across the solid solution, no intrinsic NTE is observed at the microscopic level, yet a degree of tunability in the thermal expansion coefficient with Ni substitution is demonstrated. The disparities between the intrinsic and bulk thermal expansion properties suggest that a morphological mechanism may have resulted in NTE in the bulk.

Materials Science↗

Diffusion in liquid metal systems

Physical properties of twenty liquid metals are reported; some of the data on such liquid metal properties as density, electrical resistivity, thermal conductivity, and heat capacity are summarized in graphical form. Data on laboratory handling and safety procedure are summarized for each metal; heat-transfer-correlations for liquid metals under various conditions of laminar and turbulent flow are included. Where sufficient data were available, temperature equations of properties were obtained by the method of least-squares fit. All values of properties given are valid in the given liquid phase ranges only. Additional tabular data on some 40 metals are reported in the appendix. Included is a brief description of experiments that were performed to investigate diffusion in liquid indium-gallium systems.

Ukanwa, A. O.↗

Monte-Carlo modeling of phonon thermal transport using DFT-based anisotropic dispersion relations over the full Brillouin zone

In this work, we present a Monte Carlo (MC) approach to solve the phonon Boltzmann transport equation (BTE) in which the anisotropic phonon dispersion relations over the full Brillouin zone (BZ) are used. In this approach, the discretization of the BZ used to compute the phonon relaxation time places constraints on the direction of scattered phonons in the real-space simulation domain. The phonon dispersion and phonon relaxation times are calculated using the density functional theory (DFT) approach. The modified MC approach is validated by a close examination of its ability to simulate phonon transport in both the ballistic and diffusive regimes for multiple materials including GaAs, InAs, ThO 2 , and α-U. In doing so, the phonon thermal conductivities from 100 K to 1000 K are calculated and compared with traditional non-transport solution of the phonon BTE. It is found that the phonon thermal conductivities of α-U and ThO 2 obtained from MC simulations using isotropic dispersion are larger than the values obtained using anisotropic phonon dispersion relations over the full BZ. The effect of phonon-defect scattering on the thermal conductivity of ThO 2 is also studied as an application of the current MC approach and found to agree with previously computed values in the literature. The MC solver developed here has been parallelized as a step to demonstrate its potential to solving computationally intensive phonon thermal transport problems at the mesoscale.

36 MATERIALS SCIENCE↗

Analysis of radiative sources with both specular and diffuse characteristics in the visible and infrared spectrum

The application of analytical concepts developed to facilitate the design and analysis of the thermal simulator for the Apollo Soyuz docking module is discussed. The functions used to mathematically formulate the irradiance from a quartzline lamp and reflector are described with examples of applications to other geometries. General information is denoted for utilization of the functions in computer programs with approximate computer storage and statement requirements.

Williams, F. U.↗

Permeation of halide anions through phospholipid bilayers occurs by the solubility-diffusion mechanism

Two alternative mechanisms are frequently used to describe ionic permeation of lipid bilayers. In the first, ions partition into the hydrophobic phase and then diffuse across (the solubility-diffusion mechanism). The second mechanism assumes that ions traverse the bilayer through transient hydrophilic defects caused by thermal fluctuations (the pore mechanism). The theoretical predictions made by both models were tested for halide anions by measuring the permeability coefficients for chloride, bromide, and iodide as a function of bilayer thickness, ionic radius, and sign of charge. To vary the bilayer thickness systematically, liposomes were prepared from monounsaturated phosphatidylcholines (PC) with chain lengths between 16 and 24 carbon atoms. The fluorescent dye MQAE (N-(ethoxycarbonylmethyl)-6-methoxyquinolinium bromide) served as an indicator for halide concentration inside the liposomes and was used to follow the kinetics of halide flux across the bilayer membranes. The observed permeability coefficients ranged from 10(-9) to 10(-7) cm/s and increased as the bilayer thickness was reduced. Bromide was found to permeate approximately six times faster than chloride through bilayers of identical thickness, and iodide permeated three to four times faster than bromide. The dependence of the halide permeability coefficients on bilayer thickness and on ionic size were consistent with permeation of hydrated ions by a solubility-diffusion mechanism rather than through transient pores. Halide permeation therefore differs from that of a monovalent cation such as potassium, which has been accounted for by a combination of the two mechanisms depending on bilayer thickness.

NASA Discipline Exobiology↗

Tidal influences on vertical diffusion and diurnal variability of ozone in the mesosphere

Possible dynamical influences on the diurnal behavior of ozone are investigated. A time dependent one-dimensional photochemical model is developed for this purpose; all model calculations are made at 70 deg N during summer. It is shown that the vertical diffusion can vary as much as 1 order of magnitude within a day as a result of large changes in the zonal wind induced by atmospheric thermal tides. It is found that by introducing a dissipation time scale for turbulence produced by breaking gravity waves, the agreement with Poker Flat echo data is improved. Comparisons of results from photochemical model calculations, where the vertical diffusion is a function of height only, with those in which the vertical diffusion coefficient is changing in time show large differences in the diurnal behavior of ozone between 70 and 90 km. By including the dynamical effect, much better agreement with the Solar Mesosphere Explorers data is obtained. The results are, however, sensitive to the background zonally averaged wind. The influence of including time-varying vertical diffusion coefficient on the OH densities is also large, especially between 80 and 90 km. This suggests that dynamical effects are important in determining the diurnal behavior of the airglow emission from the Meinel bands.

Bjarnason, Gudmundur G.↗

Dynamics and chemical mode analysis of plasma thermal-chemical instability

The stability of the weakly ionized plasma and the transition from a stable homogeneous discharge to unstable filaments play an important role in gas laser physics, plasma-assisted combustion, chemical reforming, and material synthesis. Here, theoretical stability analysis and thermal-chemical mode analysis were performed to understand the mechanism of plasma thermal-chemical instability by using a zero-dimensional plasma system with both simplified and detailed chemical kinetics of H 2 /O 2 /N 2 mixtures. The plasma dynamic and kinetic models accounted for multiple physical mechanisms in the chemically-reactive weakly ionized plasma, including ionization, attachment/detachment, recombination, vibrational and electronic energy relaxation, convective and diffusive species/heat removal, Joule heating, and detailed chemical kinetics. An analytical criterion and the explosive mode species/temperature pointers were formulated while the representative active species were identified for different thermal-chemical modes. The results showed that in addition to the classical thermal-ionization mechanism, various chemical modes from chemical heat imbalance and elementary kinetics significantly modified the time dynamics and the stability of the weakly ionized plasma. The present analysis provides insights and guidance to control plasma instability using chemical kinetics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of surface diffusion in electrocatalytic CO 2 reduction on Cu revealed by kinetic Monte Carlo simulations

We report Kinetic Monte Carlo (KMC) methods are frequently used for mechanistic studies of thermally driven heterogeneous catalysis systems but are underused for electrocatalysis. Here, we develop a lattice KMC approach for electrocatalytic CO 2 reduction. The work is motivated by a prior experimental report that performed electroreduction of a mixed feed of 12 CO 2 and 13 CO on Cu; differences in the 13 C content of C2 products ethylene and ethanol (Δ 13 C) were interpreted as evidence of site selectivity. The lattice KMC model considers the effect of surface diffusion on this system. In the limit of infinitely fast diffusion (mean-field approximation), the key intermediates 12 CO* and 13 CO* would be well mixed on the surface and no evidence of site selectivity could have been observed. Using a simple two-site model and adapting a previously reported microkinetic model, we assess the effects of diffusion on the relative isotope fractions in the products using the estimated surface diffusion rate of CO* from literature reports. We find that the size of the active sites and the total surface adsorbate coverage can have a large influence on the values of Δ 13 C that can be observed. Δ 13 C is less sensitive to the CO* diffusion rate as long as it is within the estimated range. We further offer possible methods to estimate surface distribution of intermediates and to predict intrinsic selectivity of active sites based on experimental observations. This work illustrates the importance of considering surface diffusion in the study of electrochemical CO 2 reduction to multi-carbon products. Our approach is entirely based on a freely available open-source code, so will be readily adaptable to other electrocatalytic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Study of the Non-Thermal X-Ray Emission of Shell-Type Supernova Remnants

We present an analysis of the X-ray spectrum of the Galactic shell-type SNR G347.3-0.5 (RX 51713.7-3946). This SNR is a member of a growing class of SNRs which are dynamically young, shell-type sources that emit non-thermal X-rays from specific regions on their outer shells. By performing a joint spectral analysis of data from observations made of G347.3-0.5 using the ROSAT PSPC, the ASCA GIS and the RXTE PCA, we have fit the spectra of particular regions of this SNR (including the bright northwestern and southwestern rims, the northeast rim and the interior diffuse emission) over the approximate energy range of 0.5 through 30 keV. We find that fits to the spectra of this SNR over this energy range using the SRCUT model were superior to a simple power law model or the SRESC model. We find that the inclusion of a thermal model with the SRCUT model helps to improve the fit to the observed X-ray spectrum: this represents the first detection of thermal X-ray emission from G347.3-0.5. Thermal emission appears to be more clearly associated with the diffuse emission in the interior of the SNR than with the bright X-ray emitting rims. A weak emission feature seen near 6.4 keV in the RXTE PCA spectrum most likely originates from diffuse X-ray emission from the surrounding Galactic Ridge rather than from G347.3-0.5 itself. We have analyzed our RXTE PCA data to search for pulsations from a recently discovered radio pulsar (PSR 51713-3949) which may be associated with G347.3-0.5, and we do not detect any X-ray pulsations at the measured radio period of 392 ms. Using the best-fit parameters obtained from the SRCUT model, we estimate the maximum energy of cosmic-ray electrons accelerated by the rims of G347.3-0.5 to be 19-25 TeV (assuming a magnetic field strength of B = 10muG), consistent with the results of Ellison et al. We present a broadband (radio to gamma ray) photon energy-flux spectrum for the northwestern rim of G347.3-0.5, where we have fit the spectrum using a more sophisticated synchrotron-inverse Compton model with a variable magnetic field strength. Our fit derived from this model yields a maximum energy of only 8.8 (+4.1) (-3.4) TeV for the accelerated cosmic-ray electrons and a much greater magnetic field strength of 150 (+250) (-80). Finally, we compare the gross properties of G347.3-0.5 with other SNRs known to possess X-ray spectra dominated by non-thermal emission.

Allen, Glenn E.↗