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.
Physical and Thermomechanical Properties of Yttrium Hydride from Large Scale Bulk Metal Hydriding Furnace
Given the superior thermal stability and highly attainable hydrogen density, yttrium hydride is an excellent high-temperature moderator material in advanced thermal neutron spectrum reactors that require small core volumes. Yttrium hydride has been selected as the moderator material for the Transformational Challenge Reactor, which was launched at Oak Ridge National Laboratory (ORNL) in 2019. However, fabrication of large-scale crack-free yttrium hydride is challenging and very limited efforts have been committed to the characterization of bulk yttrium hydride in response to the need to establish a complete database of the thermomechanical properties of YHx. In this report, the challenges associated with fabricating large-scale crack-free yttrium hydride are discussed herein. In response to those challenges, a hydriding system was designed and constructed at ORNL and was used to successfully fabricate crack-free yttrium hydride in complex geometries at large scales. This was accomplished by precisely controlling the hydrogen’s partial pressure and the retort temperature, which was informed by the well-established thermodynamic properties of the binary H-Y system. Hydrogen content in as-fabricated hydride was determined by the weight change method and vacuum hot extraction technique, complemented by the X-ray diffraction (XRD). In addition, significant efforts are being dedicated to establishing a complete database of the thermomechanical properties of as-fabricated yttrium hydride. In FY2020, we investigated the thermophysical properties of yttrium hydrides as a function of temperature (room temperature to 700°C) and hydrogen concentration (H/Y ratio ranges from 1.52 to 1.93). The results indicate that at the temperatures below 300 °C, the hydrogen content did not have a significant influence on the thermal expansion, while the specific heat capacity, the thermal diffusivity, and the calculated thermal conductivity were slightly higher for the higher H/Y ratio. Between 300°C and 700 °C, a reversible second-order endothermic transition in all measured thermal properties was observed. It was also found that the onset temperatures of the observed transition varied, with the composition having inverse dependence on the hydrogen content. An attempt was made to explain the behavior of the thermophysical properties at higher temperatures by considering the order– disorder transition as a result of hydrogen redistribution. In addition, nanoindentation was employed to determine the elastic modulus and hardness and to capture the crystal orientation dependence of these parameters. Vickers hardness was also reported. The final section of the report introduces ongoing neutron irradiation campaign of yttrium hydride.
Intra-Particle Diffusion-Limited Thermal Runaway Predictions in Lithium-Ion Systems .
Abstract not provided.
Towards inverse microstructure-centered materials design using generative phase-field modeling and deep variational autoencoders
The field of Integrated Computational Materials Engineering (ICME) combines a broad range of methods to study materials’ responses over a spectrum of length scales. A relatively unexplored aspect of microstructure-sensitive materials design is uncertainty propagation and quantification (UP/UQ) of materials’ microstructure, as well as establishing process-structure–property (PSP) relationships for inverse material design. In this study, an efficient UP technique built on the idea of changing probability measures and a deep generative unsupervised representative machine learning method for microstructure-based design of thermal conductivity of materials is proposed. Probability measures are used to represent microstructure space, and Wasserstein metrics are used to test the efficiency of the UP method. By using deep Variational AutoEncoder (VAE), we identify the correlations between the material/process parameters and the thermal conductivity of heterogeneous dual-phase microstructures. Through high-throughput screening, UP, and the deep-generative VAE method, PSP relationships that are too complex can be revealed by exploiting the materials’ design space with an emphasis on microstructures. As a last point, we demonstrate generative machine learning serves as a useful tool for inverse microstructure-centered materials design, and we demonstrate this by examining the inverse design of thermal conductivity in nano-structured materials. Here, the results reveal the effects of morphology, volume fraction, characteristic length scale, and the individual thermal diffusivity of phases on the thermal conductivity of dual-phase alloys. Our findings emphasize the advantages of high-throughput phase-field modeling and generative deep learning for linking PSP and inverse microstructure-centered materials design.
Thermoacoustic Soret separation
In a fluid mixture in a channel with an axial time-averaged temperature gradient, high-amplitude oscillating flow can greatly increase the axial flux of thermal diffusion (Soret) separation of the components of the mixture. The enhancement occurs when the oscillating lateral temperature gradient greatly exceeds the axial gradient, causing a large oscillating concentration that can be favorably time-phased with the oscillating flow. This process can occur even with a negligible pressure oscillation or with a negligible temperature response to pressure, as is the case in most liquid solutions. The thermal boundary condition imposed by realistic solids on thermoacoustic liquids is imperfect, adding mathematical complications that are absent for typical gases, for which the solid surface is temporally isothermal. Compared with gas mixtures, the high Lewis number in typical liquid solutions reduces the separation flux associated with the time-averaged temperature gradient, but it also reduces the remixing associated with the time-averaged mole-fraction gradient. For large enough channels, the second-law separation efficiency is only slightly reduced from that of steady liquid Soret separation.
Investigation of species-mass diffusion in binary-species boundary layers at high pressure using direct numerical simulations
Here, direct numerical simulations of single-species and binary-species temporal boundary layers at high pressure are performed with special attention to species-mass diffusion. The working fluids are nitrogen or a mixture of nitrogen and methane. Mean profiles and turbulent fluctuations of mass fraction show that their qualitative characteristics are different from those of streamwise velocity and temperature, due to the different boundary conditions. In a wall-parallel plane near the wall, the streamwise velocity and temperature have streaky patterns and the fields are similar. However, the mass fraction field at the same location is different from the streamwise velocity and temperature fields indicating that species-mass diffusion is not similar to the momentum and thermal diffusion. In contrast, at the centre and near the edge of the boundary layer, the mass fraction and temperature fields have almost the same pattern, indicating that the similarity between thermal and species-mass diffusion holds away from the wall. The lack of similarity near the wall is traced to the Soret effect that induces a temperature-gradient-dependent species-mass flux. As a result, a new phenomenon has been identified for a non-isothermal binary-species system – uphill diffusion, which in its classical isothermal definition can only occur for three or more species. A quadrant analysis for the turbulent mass flux reveals that near the wall the Soret effect enhances the negative contributions of the quadrants. Due to the enhancement of the negative contributions, small species-concentration fluid tends to be trapped near the wall.
Large eddy simulation of a supersonic lifted hydrogen flame: Impacts of Lewis, turbulent Schmidt and Prandtl numbers
Parametric large eddy simulations (LES) of a supersonic lifted hydrogen flame are reported. The emphases are on two aspects: impacts of (1) Lewis number (Lei of the ith species) and (2) turbulent Schmidt and Prandtl numbers (Sc t and Pr t ) on supersonic turbulent flame and flow structures. Five cases are considered: species-specific Le i , Sc t = Sr t = 1.0 (C 0 ); unity Le i , Sc t = Pr t = 1.0 (C 1 ); species-specific Sc t = 0.5, Pr t = 1.0 (C 2 ); species-specific Le i , Sc t = 1.0, Pr t = 0.5 (C 3 ); and species-specific Le i , Sc t = Pr t = 0.5 (C 4 ). Numerical results of instantaneous and/or time-averaged species mole fractions, mixture fraction, heat release rate, flame base location, and mixed modes of premixed and diffusion combustion are compared between cases C 0 and C 1 . Differences in auto-ignition locations and strengths and flame structures and stabilization specify the impacts of Lewis number. They are triggered by different predictions of species mass and thermal diffusions at fuel-coflow and/or coflow-ambient air mixing layers. These differences are rationalized by a scale analysis of mass/thermal diffusion and convection for case C 0 , which suggests the relatively low but non-negligible former against the latter. Cases C 0 and C 2 –C 4 barely see differences in terms of instantaneous and/or time-averaged temperature, velocity, and mixed combustion modes except for further downstream areas where combustion occurs. Both Sc t and Pr t impose less significant influences than Lewis number, as sub-grid scale mass/thermal diffusion is subordinate to its resolved counterpart according to their scale analysis for case C 4 .
Finite temperature properties of uranium mononitride
Uranium mononitride (UN) is a promising nuclear fuel that combines the advantageous properties of readily used UO 2 and uranium alloys, such as high melting temperature and high uranium density, and thermal conductivity, respectively. A better understanding of UN behavior at operating temperatures can be obtained from finite temperature data, such as elastic properties. To get this information, ab initio molecular dynamics (AIMD) simulations were performed at five different temperatures using constant volume (NVT) and constant pressure (NPT) ensembles. Initially, the performance of PBE functional in reproducing experimental crystallographic properties and magnetic ordering is assessed. The finite temperature phonon dispersions are calculated using NVT simulation results, which show a softening of the phonon modes with increasing temperature. The NPT results are used to obtain the thermal expansion of UN and finite temperature electronic properties. The calculated thermal expansion is compared with our measurements using neutron diffraction. Additionally, the temperature dependent elastic properties of UN are evaluated using the strain-stress method in AIMD simulations, indicating that UN becomes softer and more compressible with increasing temperature. Also, the calculated Young’s modulus slope is in very good agreement with the experiment. The finite temperature heat capacity and electronic thermal conductivity are calculated from AIMD simulations, which are in better agreement with the experiment than the heat capacity and thermal conductivity calculated using the structures relaxed at 0K. Finally, the thermal diffusivity from AIMD has opposite temperature dependence compared to experimental results, which we argued comes from the underestimated electronic thermal conductivity.
Plasma thermal-chemical instability of low-temperature dimethyl ether oxidation in a nanosecond-pulsed dielectric barrier discharge
Plasma stability in reactive mixtures is critical for various applications from plasma-assisted combustion to gas conversion. To generate stable and uniform plasmas and control the transition towards filamentation, the underlying physics and chemistry need a further look. Here, this work investigates the plasma thermal-chemical instability triggered by dimethyl-ether (DME) low-temperature oxidation in a repetitive nanosecond pulsed dielectric barrier discharge. First, a plasma-combustion kinetic mechanism of DME/air is developed and validated using temperature and ignition delay time measurements in quasi-uniform plasmas. Then the multi-stage dynamics of thermal-chemical instability is experimentally explored: the DME/air discharge was initially uniform, then contracted to filaments, and finally became uniform again before ignition. By performing chemistry modeling and analyzing the local thermal balance, it is found that such nonlinear development of the thermal-chemical instability is controlled by the competition between plasma-enhanced low-temperature heat release and the increasing thermal diffusion at higher temperature. Further thermal-chemical mode analysis identifies the chemical origin of this instability as DME low-temperature chemistry. This work connects experiment measurements with theoretical analysis of plasma thermal-chemical instability and sheds light on future chemical control of the plasma uniformity.
Influence of pre-existing defects on thermal transport in nuclear graphite
Nuclear graphite is a critical material in high-temperature nuclear reactors due to its superior thermal and mechanical properties. The manufacturing process leaves multi-scale ‘pre-existing’ defects that can affect thermal transport characteristics. Because these defects are remnant of graphitization temperature, they cannot be thermally annealed. This study employs a non-thermal electron wind force (EWF) annealing technique to avoid this obstacle. 2 min of EWF treatment of the as-received graphite IG-110 at temperatures below 100 °C led up to 67% increase in thermal diffusivity and ~ 35% decrease in electrical resistivity in average. Differential scanning calorimetry also showed similar outcome for specific heat. X-ray diffraction characterization was performed by fitting a bi-modal crystallite size distribution model to reveal the enhancement in crystallinity after the EWF treatment. The findings emphasize the potential of EWF annealing for optimizing thermal performance in nuclear graphite and its implications for reactor efficiency and safety.
Thermophysical properties of high-density, sintered monoliths of yttrium dihydride in the range 373–773 K
Yttrium dihydride is a promising nuclear reactor moderator for high-temperature, thermal micro-reactors due to its ability to retain hydrogen to high temperatures while having a relatively low impact on neutron economy. However, it is difficult to fabricate yttrium dihydride in high-density, near net-shape monoliths for moderator applications due to challenges associated with the hydrogen absorption process. Here, high-density monoliths of yttrium dihydride were prepared by powder metallurgical methods. The lattice strain and thermal expansion coefficient were measured using neutron diffraction of powders of directly-hydrided yttrium, while the molar heat capacity and thermal diffusivity of sintered yttrium dihydride were measured using differential scanning calorimetry (DSC) and laser flash analysis (LFA), respectively. The molar heat capacity of yttrium dihydride was also calculated using density functional theory (DFT) for comparison. These parameters were used to then calculate the thermal conductivity and resistivity of yttrium dihydride as a function of temperature. The thermophysical properties of materials produced by both methods were observed to be consistent with the values from literature for yttrium dihydride. Due to the novelty of producing yttrium dihydride by powder metallurgical methods, this result indicated that the sintered monoliths were of high quality and that powder metallurgy is a viable method for large-scale production of yttrium dihydride monoliths for nuclear reactor moderator applications.
Observation of reduced thermal conductivity in a metal-organic framework due to the presence of adsorbates
Whether the presence of adsorbates increases or decreases thermal conductivity in metalorganic frameworks (MOFs) has been an open question. Here we report observations of thermal transport in the metal-organic framework HKUST-1 in the presence of various liquid adsorbates: water, methanol, and ethanol. Experimental thermo-reflectance measurements were performed on single crystals and thin films, and theoretical predictions were made using molecular dynamics simulations. We find that the thermal conductivity of HKUST-1 decreases by 40 – 80% depending on the adsorbate, a result that cannot be explained by effective medium approximations. Our findings demonstrate that adsorbates introduce additional phonon scattering in HKUST-1, which particularly shortens the lifetimes of low frequency phonon modes. As a result, the system thermal conductivity is lowered to a greater extent than the increase expected by the creation of additional heat transfer channels. Finally, we show that thermal diffusivity is even more greatly reduced than thermal conductivity by adsorption.
Diffuse scattering from dynamically compressed single-crystal zirconium following the pressure-induced α → ω phase transition
The prototypical α → ω phase transition in zirconium is an ideal test bed for our understanding of polymorphism under extreme loading conditions. After half a century of study, a consensus had emerged that the transition is realized via one of two distinct displacive mechanisms, depending on the nature of the compression path. However, recent dynamic-compression experiments equipped with diffraction diagnostics performed in the past few years have revealed new transition mechanisms, demonstrating that our understanding of the underlying atomistic dynamics and transition kinetics is in fact far from complete. We present classical molecular dynamics simulations of the α → ω phase transition in single-crystal zirconium shock compressed along the [0001] axis using a machine-learning-class potential. The transition is predicted to proceed primarily via a modified version of the two-stage Usikov-Zilberstein mechanism, whereby the high-pressure ω phase heterogeneously nucleates at boundaries between grains of an intermediate β phase. We further observe the fomentation of atomistic disorder at the junctions between β grains, leading to the formation of highly defective interstitial material between the ω grains. We directly compare synthetic x-ray diffraction patterns generated from our simulations with those obtained using femtosecond diffraction in recent dynamic-compression experiments, and show that the simulations produce the same unique, anisotropic diffuse scattering signal unlike any previously seen from an elemental metal. Our simulations suggest that the diffuse signal arises from a combination of thermal diffuse scattering, nanoparticlelike scattering from residual kinetically stabilized α and β grains, and scattering from interstitial defective structures. Published by the American Physical Society 2024
Direct numerical simulation of single-species and binary-species boundary layers at high pressure
Direct numerical simulations of single-species and binary-species temporal boundary layers at high pressure are performed. The main objective is to investigate flow physics in a binary-species boundary layer at high pressure, with special attention to mass diffusion. The working fluids are nitrogen in the single-species cases, and a mixture of nitrogen and methane in the binary-species cases. An investigation of mean profiles shows that velocity and temperature profiles have steep gradients near the wall, whereas the mass fraction profiles do not have a steep gradient. This result indicates that the similarity law between velocity, temperature and mass fraction does not hold true in the binary-species boundary layer. The comparison of turbulent fluctuations shows that the qualitative characteristics of velocity and temperature are similar to each other. In contrast, profiles of mass fraction fluctuations are largely different from those of velocity and temperature. This result indicates that mass diffusion is not similar to momentum and thermal diffusions in binary-species boundary layers.
A Computational Study of the Thermodynamic Conditions Leading to Autoignition in Nanosecond Pulsed Discharges
Nanosecond pulsed discharges have attracted the attention of engine manufacturers due to the possibility of attaining distributed ignition sites that accelerate burn rates while resulting in very little electrode erosion. Multidimensional modeling tools currently capture the electrical structure of such discharges accurately, but resolving the chemical structure remains a challenging problem owing to the disparity of time-scales in streamer propagation (nanoseconds) and ignition phenomena (microseconds). The purpose of this study is to extend multidimensional results toward resolving the chemical structure in the wake of streamers (or the afterglow) by using a batch reactor model (BRM). This can afford the use of very detailed chemical kinetic information. The full nonequilibrium nature of the electrons is taken into account, along with fast gas heating, shock wave propagation, and thermal diffusion. The results shed light on ignition phenomena brought about by such discharges.
Towards performing high-resolution inelastic X-ray scattering measurements at hard X-ray free-electron lasers coupled with energetic laser drivers
High-resolution inelastic X-ray scattering is an established technique in the synchrotron community, used to investigate collective low-frequency responses of materials. When fielded at hard X-ray free-electron lasers (XFELs) and combined with high-intensity laser drivers, it becomes a promising technique for investigating matter at high temperatures and high pressures. This technique gives access to important thermodynamic properties of matter at extreme conditions, such as temperature, material sound speed, and viscosity. The successful realization of this method requires the acquisition of many identical laser-pump/X-ray-probe shots, allowing the collection of a sufficient number of photons necessary to perform quantitative analyses. Here, a 2.5-fold improvement in the energy resolution of the instrument relative to previous works at the Matter in Extreme Conditions (MEC) endstation, Linac Coherent Light Source (LCLS), and the High Energy Density (HED) instrument, European XFEL, is presented. Some aspects of the experimental design that are essential for improving the number of photons detected in each X-ray shot, making such measurements feasible, are discussed. A careful choice of the energy resolution, the X-ray beam mode provided by the XFEL, and the position of the analysers used in such experiments can provide a more than ten-fold improvement in the photometrics. The discussion is supported by experimental data on 10 µm-thick iron and 50 nm-thick gold samples collected at the MEC endstation at the LCLS, and by complementary ray-tracing simulations coupled with thermal diffuse scattering calculations.
Faster ablative Kelvin–Helmholtz instability growth in a magnetic field
Shear flows along a plasma interface will quickly grow unstable due to the Kelvin–Helmholtz instability. If there is a concurrent temperature gradient across the interface, higher modes are stabilized by the thermal diffusion. These ablative effects must be considered in, for example, jet features in inertial confinement fusion hot-spots, or plasma plumes in young supernovae. We show that magnetization of the plasma can greatly affect the instability, even if magnetic pressure is small. This is because electrons are localized by their gyromotion, reducing the heat flux and material ablation. We use a two-dimensional numerical extended-magnetohydrodynamics approach to assess this effect for dense fusion conditions. In comparison with the unmagnetized case, self-generated Biermann fields make only a minor difference to growth rates. However, simulations with a large 50 kT external field found that the growth rate of the least stable mode increased by 40%. In conclusion, this has implications for mix processes in Z-pinches and magnetized inertial confinement fusion concepts.
Preliminary Reporting of Thermophysical Property Measurements for the Ghareb Formation
Accurate knowledge of thermophysical properties of rock is vital to develop meaningful models of high level nuclear waste emplacement scenarios. The Israel Atomic Energy Commission is considering storing high level nuclear waste in the Ghareb formation, a porous kerogen bearing chalk. Sandia is supporting this effort with an evolving lab- based geomechanics testing program. We have completed measurements of thermal properties up to 275C and room temperature hydrostatic compaction measurements. We report thermal conductivity, thermal diffusivity, specific heat, and mass loss from our thermal measurements, and we report bulk moduli and porosity loss from our compaction measurements. These values are crucial for the numerical models to simulate heat transfer and formation compressibility around a heat generating repository.