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,189 records · Page 66

The gabbro-eclogite phase transition and the elevation of mountain belts on Venus

The hypothesis is explored that the crust-mantle boundary of Venus is not in phase equilibrium but rather is rate-limited by the temperature-dependent volume diffusion of the slowest ionic species. The 1D thermal evolution problem is solved assuming that the mountains formed by uniform horizontal shortening of the crust and the lithospheric mantle at a constant rate. The time-dependent density structure and surface elevation are calculated by assuming a temperature-dependent reaction rate and local Airy isostatic compensation. For a horizontal strain rate of 10 exp -15/s or greater, the temperature increase at the base of the crust during mountain formation is modest to negligible, the deepening lower crust is metastable, and the surface elevation increases as the crust thickens. For strain rates less than 10 exp -16/s, crustal temperature increases with time because of internal heat production and the lower crust is more readily transformed to the dense eclogite assemblage. For such models, a maximum elevation is reached during crustal shortening.

Namiki, Noriyuki↗

An Analysis of Neptune's Stratospheric Haze Using High-Phase-Angle Voyager Images

We have inverted high-phase-angle Voyager images of Neptune to determine the atmospheric extinction coefficient as a function of altitude and the scattering phase function at a reference altitude. Comparisons between theoretical model and observations help separate the contributions from molecular Rayleigh and aerosol scattering and help determine the variation of the aerosol size, concentration, and scattering properties with altitude. Further comparisons between models and data allow us to place constraints on the location and composition of the hazes, the concentration and downward flux of certain condensible hydrocarbon gases, the eddy diffusion coefficient in the lower stratosphere, and the thermal profile in parts of Neptune's stratosphere. We find that a distinct stratospheric haze layer exists near 12(sub -1, sup +1) mbar in Neptune's lower stratosphere, most probably due to condensed ethane. The derived stratospheric haze production rate of 1.0(sub -0.3, sup +0.2) x 10(exp -15) g cm(exp -2) sec(exp -1) is substantially lower than photochemical model predictions. Evidence for hazes at higher altitudes also exists. Unlike the situation on Uranus, large particles (0.08-0.11 microns) may be present at high altitudes on Neptune (e.g., near 0.5 mbar), well above the region in which we expect the major hydrocarbon species to condense. Near 28 mbar, the mean particle size is about 0.13(sub -0.02, sup +0.02) microns with a concentration of 5(sub -3, sup +3) particles cm(exp -3). The cumulative haze extinction optical depth above 15 mbar in the clear filter is approx. 3 x 10(exp -3), and much of this extinction is due to scattering rather than absorption; thus, if our limb-scan sites are typical, the hazes cannot account for the stratospheric temperature inversion on Neptune and may not contribute significantly to atmospheric heating. We compare the imaging results with the results from other observations, including those of the Voyager Photopolarimeter Subsystem, and discuss differences between Neptune and Uranus.

Moses, Julianne I.↗

Multi-Absorber Transition-Edge Sensors for X-Ray Astronomy Applications

We are developing multi-absorber Transition-Edge Sensors (TESs) for applications in x-ray astronomy. These position-sensitive devices consist of multiple x-ray absorbers each with a different thermal coupling to a single readout TES. Heat diffusion between the absorbers and the TES gives rise to a characteristic pulse shape corresponding to each absorber element and enables position discrimination. The development of these detectors is motivated by a desire to maximize focal plane arrays with the fewest number of readout channels. In this contribution we report on the first results from devices consisting of nine) 65 X 65 sq. microns Au x-ray absorbers) 5 microns thick. These are coupled to a single 35 X 35 sq. microns Mo/Au bilayer TES. These devices have demonstrated full-width-half-maximum (FWHM) energy resolution of 2.1 eV at 1.5 keV) 2.5 eV at 5.9 keV and 3.3 eV at 8 keV. This is coupled with position discrimination from pulse shape over the same energy range. We use a finite-element model to reproduce the measured pulse shapes and investigate the detector non-linearity with energy) which impacts on the devices position sensitivity and energy resolution.

Smith, S. J.↗

Metallized Nanotube Polymer Composite (MNPC) and Methods for Making Same

A novel method to develop highly conductive functional materials which can effectively shield various electromagnetic effects (EMEs) and harmful radiations. Metallized nanotube polymer composites (MNPC) are composed of a lightweight polymer matrix, superstrong nanotubes (NT), and functional nanoparticle inclusions. MNPC is prepared by supercritical fluid infusion of various metal precursors (Au, Pt, Fe, and Ni salts), incorporated simultaneously or sequentially, into a solid NT-polymer composite followed by thermal reduction. The infused metal precursor tends to diffuse toward the nanotube surface preferentially as well as the surfaces of the NT-polymer matrix, and is reduced to form nanometer-scale metal particles or metal coatings. The conductivity of the MNPC increases with the metallization, which provides better shielding capabilities against various EMEs and radiations by reflecting and absorbing EM waves more efficiently. Furthermore, the supercritical fluid infusion process aids to improve the toughness of the composite films significantly regardless of the existence of metal.

Park, Cheol↗

Recent Developments to the Porous Microstructure Analysis (PuMA) Software

The Porous Microstructure Analysis (PuMA) software is a suite of tools for the analysis of porous materials and generation of material microstructures. From microstructural data, often obtained through X-ray microtomography, PuMA can determine a number of effective material properties and perform material response simulations. Version 2.2 includes capabilities for computing volume fractions, porosity, specific surface area, effective thermal and electrical conductivities, and continuum and rarefied diffusive tortuosity. PuMA can also simulate competitive diffusion/reaction processes at the micro-scale, such as surface oxidation. In this poster, recent advancements to the PuMA software are detailed, including the full refactoring of PuMA into v3.0, a new module to compute heat conduction in anisotropic materials, a particle method for simulating molecular beam experiments, a new finite-volume Laplace solver, complex fibrous material generation, woven material generation, and a coupling of PuMA with the DAKOTA software for advanced statistics.

PuMA↗

Ice Growth Model and Lunar Tanker for Polar Water Extraction

In-Situ Resource Utilization (ISRU) is a major component of sustainable and cost-effective missions to the Moon and beyond. Ever since it has been identified in the Permanently Shadowed Regions (PSRs) of the Moon, water has become a resource of interest for it’s wide range of uses such as life support and hydrogen and oxygen rocket propellant. It is ice-bound in the regolith, so separation typically involves heating to sublimation temperatures and capturing at very low pressure. One option is to then store it in ice form and transport it out of the PSR and into a sunlit ridge. From here, processes such as electrolysis can be performed. An understanding of frost growth dynamics under rarefied conditions is necessary for successful water production, though. An analytical model was developed using a Diffusion-Limited Aggregation (DLA) approach to evaluate bulk thermal conductivity changes of frost under varying conditions. This model will be evaluated using a one-dimensional frost growth experiment on a horizontal, flat plate under vacuum and cryogenic conditions. The tests will measure frost growth rate and bulk density, and the thermal conductivity of the ice layer as a function of thickness for different plate temperatures. The test data will then be used to validate and improve the analytical model.

Lunar water↗

Loci/STREAM Sharp Interface Model Validation of Cryogenic Propellant Tank Self-Pressurization

This paper documents a model validation effort for the prediction of cryogenic propellant tank selfpressurization physics using the Loci/STREAM sharp interface CFD model developed at NASA MSFC. Cryogenic fluid management (CFM) applications for long duration spaceflight missions require accurate modeling of heat transfer and thermodynamics. Therefore, detailed validation of those physics is important to ensure simulation accuracy. In the past, the MSFC Propulsion Fluid Dynamics Branch has validated and used the Loci/STREAM Volume of Fluid (VOF) module extensively for propellant slosh applications, which are dominated by fluid motion. For cases where the interface gas/liquid interface is essentially static, the Loci/STREAM sharp interface was developed for more rapid and efficient CFM analysis. In the present study, Loci/STREAM sharp interface simulations were compared with the NASA K-site self-pressurization experiment of a flight weight, partially full, LH2 tank in normal gravity. This study investigated several parametric effects and made a number of important observations on the performance of the Loci/STREAM sharp interface method: - First, a homogeneous model was derived to verify the CFD model. The homogeneous model assumes both gas and liquid phases are uniform and are at saturation temperature. The predicted pressure rise rates at 3.5 W/m2 and 2.0 W/m2 with a 49% fill level are all lower than the experimental measurement. This is expected due to the observed temperature stratification near the gas-liquid interface during testing. - The k-w and k-e turbulence models predict comparable pressurization rates, which are much lower than the experiment measurement. On the other hand, the pressure values and pressurization rates from both models are almost the same as those of the homogeneous model. The current turbulence model is too diffusive and is not capable of capturing the thermal stratification in the liquid and ullage. - When using the laminar flow model, the pressurization rate is found to be very close to the experimental value when a quasi-steady flow field is allowed to develop by continuously venting the tank prior to pressurization as was done during testing. On the other hand, without initial venting, the rate is not constant. Additionally, the development time to reach a quasi-steady flow field differed from the experiment due to the lack of turbulent mixing away from the gas-liquid interface that is expected in reality. - Using the laminar flow model and appropriate initial conditioning, the pressurization rate is within 15% of the experimental value for 49% fill level at 3.5 W/m 2 and is within 18% for 2.0 W/m 2 heat loading. - Using the laminar flow model and appropriate initial conditioning, the pressurization rate is within 18% of the experimental value for a 29% fill level and within 15% for an 83% fill level.

H. Q. Yang↗

Dynamic System Scaling Application to Accelerated Nuclear Fuel Testing

The development of nuclear fuel and materials requires a continuous effort to investigate the acute and prolonged effects of irradiation, thermal-material stress, chemical change, or other conceivable damage mechanics acquired during normal operation or accident scenarios throughout its lifetime. As in-core fuel property measurement techniques advance to support in the real-time, non-invasive, and enhanced accuracy realm, it is the future of fuel development to pursue to higher degrees of control, predictability of integral test behavior via separate effects test (SET), and shorter test time intervals. The fuel development life cycle from initial concept to commercial licensing is approximated to be 20 years and current literature suggests by optimizing fuel performance codes with SETs, the process could possibly be compressed to 5 to 10 years. Recently, a research group in the Idaho National Laboratory (INL) is testing reduced scale fuel rods and increased power density to accelerate evolution of fuel phenomena in metallic fuels. In support of nuclear fuel rod development, compressing fuel test process, and accelerating fuel phenomena, it was the purpose of this study to investigate nuclear fuel performance phenomena via literature review and effectively scale the initial conditions, boundary conditions, and geometric properties to describe the time-dependent response including to fuel burnup, thermal and mechanical stress, transmutations and inter-diffusion, and other relevant observed phenomena. The study was based on BISON simulations of historic EBR-II metallic fuel experiments and Dynamical System Scaling (DSS) method are utilized to assess effects of scaling activity including fuel phenomena acceleration and calculations of time-dependent distortions. The research successfully scaled metallic fuel phenomena, accelerated fuel testing, and assessed the distortions for each scaled case derived.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Monte Carlo analysis of lobular gas-surface scattering in tubes applied to thermal transpiration

A model of rarefied gas flow in tubes was developed which combines a lobular distribution with diffuse reflection at the wall. The model with Monte Carlo techniques was used to explain previously observed deviations in the free molecular thermal transpiration ratio which suggest molecules can have a greater tube transmission probability in a hot-to-cold direction than in a cold-to-hot direction. The model yields correct magnitudes of transmission probability ratios for helium in Pyrex tubing (1.09 to 1.14), and some effects of wall-temperature distribution, tube surface roughness, tube dimensions, gas temperature, and gas molecular mass.

Smith, J. D.↗

Carbon doping in molecular beam epitaxy of GaAs from a heated graphite filament

Carbon doping of GaAs grown by molecular beam epitaxy has been obtained for the first time by use of a heated graphite filament. Controlled carbon acceptor concentrations over the range of 10 to the 17th-10 to the 20th/cu cm were achieved by resistively heating a graphite filament with a direct current power supply. Capacitance-voltage, p/n junction and secondary-ion mass spectrometry measurements indicate that there is negligible diffusion of carbon during growth and with postgrowth rapid thermal annealing. Carbon was used for p-type doping in the base of Npn AlGaAs/GaAs heterojunction bipolar transistors. Current gains greater than 100 and near-ideal emitter heterojunctions were obtained in transistors with a carbon base doping of 1 x 10 to the 19th/cu cm. These preliminary results indicate that carbon doping from a solid graphite source may be an attractive substitute for beryllium, which is known to have a relatively high diffusion coefficient in GaAs.

Malik, R. J.↗

Ab-initio informed cluster dynamics simulation of self- and Xe diffusivity in uranium mononitride under irradiation

Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuels considered for light water reactor and advanced reactor designs, as it presents significant advantages such as high uranium density (better economics) and high thermal conductivity and melting point (increased safety). Self- and fission gas diffusivities need to be better understood, given that they influence key fuel performance phenomena like swelling and fission gas release. Recently, radiation enhanced diffusivity was investigated in UN by means of cluster dynamics simulations relying on empirical potential-based parameterizations, the reliability of which highly depends on the interatomic potential accuracy. Here, in this work, we refine this approach by determining, using ab-initio calculations, the properties of defect clusters containing vacancies, self-interstitials and Xe impurities. We also consider larger clusters than previous studies. The obtained dataset (formation enthalpies, entropies, and migration barriers) is used to parameterize a cluster dynamics model of mobile clusters, and to calculate the defect cluster concentrations under irradiation. This gives us access to the radiation enhanced self- and fission gas diffusivities. Although the resulting diffusivities are close to the values reported in the literature, we find important qualitative differences in the diffusion mechanisms. Capturing the correct mechanisms is crucial to properly describe the chemistry and fission rate dependence of the model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling the Thermostructural Stability of Melt-infiltrated Sic/sic Composites

SiC/SiC composites developed by NASA with Sylramic-iBN fibers and melt-infiltrated (MI) SiC-Si matrices have demonstrated 1000-hour rupture life in air at 100 MPa and 1315OC. Recently it has been determined that a major factor controlling the long-term rupture life of these composites is not environment or stress, but an intrinsic microstructural and strength instability caused by a thermally-induced silicon attack of the Sic fibers. The objective of this paper is to present a simple diffusion-based analytical model which predicts well the observed effects of stress-free thermal exposure on the residual tensile strength of Sylramic-iBN/SiC-Si composites. The practical implications of the model for SiC/SiC composites with MI matrices are discussed.

DiCarlo, James A.↗

Shape fluctuations and optical transition of He 2 * excimer tracers in superfluid He 4

Metastable He$_{2}^{*}$ excimer molecules have been utilized as tracer particles of the normal component in superfluid $^4$He (He II) which can be imaged via laser-induced fluorescence. These excimer molecules form tiny bubbles in He II and can bind to quantized vortices at sufficiently low temperatures, thereby allowing for direct visualization of vortex dynamics in an inviscid superfluid. However, the $a^{3}\Sigma^+_{u}$${\rightarrow}$$c^{3}\Sigma^+_{g}$ optical absorption line, which is responsible for the fluorescence imaging of the He$_{2}^{*}$ molecules, is controlled by fluctuations on the bubble shape, and its exact line profile is not known at low temperatures. In this paper, we present a bubble model for evaluating the surface fluctuation eigenmodes of the excimers in He II. The line profile of the $a^{3}\Sigma^+_{u}{\rightarrow}c^{3}\Sigma^+_{g}$ transition is calculated at different temperatures by considering both the zero-point and thermal fluctuations on the bubble shape. We show that, as the temperature drops from 2~K to 20 mK, the peak absorption strength is enhanced by a factor of about five, accompanying a blueshift of the peak location by about 2 nm. A double-peak line profile due to the rotational levels of the molecular core can be resolved. This bubble model also allows us to evaluate the stiffness of the He$_{2}^{*}$ bubbles and hence their diffusion constant in He II due to scattering off thermal phonons. Our results will aid the design of future experiments on imaging quantized vortices in He II using He$_{2}^{*}$ tracers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Comparison of the Radiative Two-Flux and Diffusion Approximations

Approximate solutions are sometimes used to determine the heat transfer and temperatures in a semitransparent material in which conduction and thermal radiation are acting. A comparison of the Milne-Eddington two-flux approximation and the diffusion approximation for combined conduction and radiation heat transfer in a ceramic material was preformed to determine the accuracy of the diffusion solution. A plane gray semitransparent layer without a substrate and a non-gray semitransparent plane layer on an opaque substrate were considered. For the plane gray layer the material is semitransparent for all wavelengths and the scattering and absorption coefficients do not vary with wavelength. For the non-gray plane layer the material is semitransparent with constant absorption and scattering coefficients up to a specified wavelength. At higher wavelengths the non-gray plane layer is assumed to be opaque. The layers are heated on one side and cooled on the other by diffuse radiation and convection. The scattering and absorption coefficients were varied. The error in the diffusion approximation compared to the Milne-Eddington two flux approximation was obtained as a function of scattering coefficient and absorption coefficient. The percent difference in interface temperatures and heat flux through the layer obtained using the Milne-Eddington two-flux and diffusion approximations are presented as a function of scattering coefficient and absorption coefficient. The largest errors occur for high scattering and low absorption except for the back surface temperature of the plane gray layer where the error is also larger at low scattering and low absorption. It is shown that the accuracy of the diffusion approximation can be improved for some scattering and absorption conditions if a reflectance obtained from a Kubelka-Munk type two flux theory is used instead of a reflection obtained from the Fresnel equation. The Kubelka-Munk reflectance accounts for surface reflection and radiation scattered back by internal scattering sites while the Fresnel reflection only accounts for surface reflections.

Spuckler, Charles M.↗

The Chandra M10l Megasecond: Diffuse Emission

Because MIOl is nearly face-on, it provides an excellent laboratory in which to study the distribution of X-ray emitting gas in a typical late-type spiral galaxy. We obtained a Chandra observation with a cumulative exposure of roughly 1 Ms to study the diffuse X-ray emission in MlOl. The bulk of the X-ray emission is correlated with the star formation traced by the FUV emission. The global FUV/Xray correlation is non-linear (the X-ray surface brightness is roughly proportional to the square root of the FUV surface brightness) and the small-scale correlation is poor, probably due to the delay between the FUV emission and the X-ray production ill star-forming regions. The X-ray emission contains only minor contributions from unresolved stars (approximates less than 3%), unresolved X-ray point sources (approximates less than 4%), and individual supernova remnants (approximates 3%). The global spectrum of the diffuse emission can be reasonably well fitted with a three component thermal model, but the fitted temperatures are not unique; many distributions of emission measure can produce the same temperatures when observed with the current CCD energy resolution. The spectrum of the diffuse emission depends on the environment; regions with higher X-ray surface brightnesses have relatively stronger hard components, but there is no significant evidence that the temperatures of the emitting components increase with surface brightness.

Kuntz, K. D.↗

Pressure effects on reactivity and extinction of n -dodecane diffusion cool flame

Low-temperature combustion at high pressure has gained increasing interest owing to the growing demand for clean and efficient propulsion and energy systems. In particular, cool flame structure, extinction limits, and fuel reactivities at high pressure are critical in affecting the engine performance at near-limit conditions. Here, this paper investigates the effects of pressure on the cool flame extinction limit, structure, radical index, reactivity, and oxygen concentration dependence by experiments, analysis, and modeling. A large n-alkane, namely n-dodecane (n C 12 H 26 ), is selected to study its diffusion cool flame dynamics and reactivity in a high-pressure counterflow burner up to 10 atm. The experimental results show that higher pressure increases the cool flame extinction strain rates, and that the pressure-weighted extinction strain rate (aP) is proportional to the square of pressure, P 2 . A scaling analysis explains the relationship between the dependence of flame structure, heat release rate, and pressure-weighted cool flame strain rate on pressure. Furthermore, radical indexes at different pressures are measured by isolating the thermal and transport effects from the chemical contribution to diffusion cool flame extinction. The radical index clearly shows that the low-temperature reactivity increases with pressure. In addition, due to the critical role of multiple oxygen addition reactions in low-temperature chemistry, the relationship between the cool flame extinction limit and the oxygen concentration is explored. It is found that the cool flame extinction limits are proportional to the n th power of the oxygen concentration, [O 2 ] n , and increasing pressure leads to stronger extinction limit dependence (larger n) on the oxygen concentration. The present experiment and detailed kinetic analysis show clearly that increasing pressure promotes the low-temperature chemistry including the oxygen addition reactions, while the scaling analysis explains well the experimental results.

33 ADVANCED PROPULSION SYSTEMS↗

Simulations and Experiments Reveal Effect of Nanopores on Helium Diffusion in Quartz

The diffusion properties of noble gases in minerals are widely used to reconstruct the thermal histories of rocks. In this work, we combine density functional theory (DFT) calculations with laboratory experiments to investigate controls on helium diffusion in quartz. DFT calculations for perfect α-quartz predict substantially lower activation energies and frequency factors for helium diffusion than observed in laboratory experiments, especially in the [001] direction. These results imply that no helium could be retained in quartz at Earth surface temperatures, which conflicts with observations of partial cosmogenic 3 He retention over geologic time scales. Here, we implement a model of helium diffusion in α-quartz modulated by nanopore defects that disrupt energetically favorable diffusion pathways. In this model, we find that laboratory-determined diffusivities can be most closely reproduced when a helium atom encounters ~70 nanopore sites per million interstitial sites. The results of our model indicate that diffusion of helium in natural quartz, like other noble gases in other minerals, can be significantly modulated by extended defects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparison of DeePMD, MTP, GAP, ACE and MACE Machine‐Learned Potentials for Radiation‐Damage Simulations: A User Perspective

Accurate and efficient interatomic potentials are essential for molecular dynamics (MD) simulations of radiation damage, gas diffusion, and phase stability in complex ceramics such as LiAlO 2 , especially under extreme conditions relevant to tritium production. Here, we evaluate the performance of six machine-learned interatomic potentials (MLIPs), moment tensor potential (MTP), Gaussian approximation potential, deep potential (DeePMD), atomic cluster expansion (ACE), message-passing ACE (multilayer atomic cluster expansion (MACE) pretrained) and MACE (trained from-scratch), all trained on the same density functional theory dataset with inclusion of tritium. The MLIPs are benchmarked against traditional Buckingham and ReaxFF potentials in terms of energy accuracy, density predictions, thermal equilibration behavior, threshold displacement energy (E d ), tritium diffusivity, and computational cost. Among the models, MTP shows the best overall balance between efficiency and accuracy, with low force and energy errors and realistic E d values for Li and Al. The ACE and MACE (pretrained and trained from scratch) models exhibit high E d (>200 eV) and unphysical pair interactions. DeePMD underestimates Ed due to overly repulsive behavior even at equilibrium distances. All models over-estimate tritium diffusion but the pretrained MACE model behaves well during tritium-diffusion simulations up to 500 K, maintaining diffusivities in the physically consistent 10 −11 m 2 /s range. Finally, we quantify the computational cost of each potential in large-scale atomic/molecular massively parallel simulator, finding that only MTP is more efficient than traditional empirical potentials, while others are significantly more expensive. These findings explain the trade-offs between accuracy and computational cost in MLIP development and provide essential guidance for use in high-throughput radiation damage and gas diffusion simulations in nuclear ceramics.

74 ATOMIC AND MOLECULAR PHYSICS↗