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 397 records · Page 22

Influence of Tabula Rasa on Process- and Light-Induced Degradation of Solar Cells Fabricated From Czochralski Silicon

Monocrystalline Si solar cells are fabricated from Czochralski (Cz) Si, which contains 10 17 –10 18 cm -3 oxygen atoms. Cz Si undergoes degradation during high-temperature thermal processing steps, such as dopant diffusion to form the p - n junction. This degradation in the bulk minority carrier lifetime can be related to the formation of oxygen precipitates. We found that a high-temperature annealing process known as tabula rasa (TR) not only mitigates process-induced degradation via oxygen precipitate nuclei dissolution, but also modifies subsequent light-induced degradation. We report on the bulk carrier lifetime of n - and p -type Cz Si after TR, which homogenizes the interstitial oxygen in the bulk Si to its monoatomic form in either an N 2 or O 2 environment. A control sample, which was not subjected to a TR processing step, experienced severe process-induced degradation during a boron emitter thermal budget as oxygen precipitates were formed in the Si bulk. These precipitates could be imaged using photoluminescence. Additionally, samples that underwent a TR processing step prior to the boron emitter thermal budget show efficient gettering of metallic impurities compared to the control sample, which showed a decline in the implied open-circuit voltage after the gettering step. Furthermore, modification of the interstitial oxygen bonding by TR had a strong effect on the light-induced degradation kinetics. Instead of a typically observed monotonic decrease, minority carrier lifetime increases first, followed by a nonmonotonic decrease over a ~20 h period. Finally, we conclude that by modifying the interstitial oxygen bonding via TR pretreatment, p -type Cz Si wafers become substantially resistant to harsh solar cell processes and strongly modified light-induced degradation, which would open alternative ways to mitigate this degradation mechanism.

14 SOLAR ENERGY↗

Application of BISON to UO 2 MiniFuel fission gas release analysis

There has been a recent push to accelerate fuel qualification by developing revolutionary capabilities to reduce irradiation periods, and thereby, reduce the time required to qualify a new fuel system. One such capability is the MiniFuel irradiation capsule designed to miniaturize fuel samples and irradiate “mini” fuel samples under isothermal temperature conditions. MiniFuel allows steady-state irradiations to decouple the traditionally coupled fission rate (i.e., power) and temperature parameters to understand and generate microstructures observed in fuel operated in a commercial reactor. Furthermore, this process offers the possibility to gather in situ data as well as postirradiation or transient data such as thermal conductivity, specific heat, fission gas diffusion and release, etc. However, accelerating fuel qualification is not solely reliant on generating large amounts of data but also on developing an informed test matrix designed to rapidly generate impactful data. Additionally, this process is reliant on fuel performance codes, such as BISON, to evaluate MiniFuel irradiations using existing material models. This process pinpoints model/data gaps, identifies desired irradiation conditions, and subsequently supports model validation and development. This work describes the use of BISON to perform a number of sensitivity studies designed to understand conditions that lead to fission gas release (FGR) under steady-state isothermal irradiation conditions and temperature transient conditions. The model is applied to a UO 2 MiniFuel example and shows an overall good qualitative agreement with experimental FGR annealing tests under different temperature conditions. It also accounts well for microstructural effects on FGR. When quantitatively compared with FGR data from previously irradiated 103 MWd/kgU UO 2 discs under thermal annealing, the model shows a less satisfactory agreement with the experimental data. Finally, a UO 2 MiniFuel test matrix is proposed to help to extend the model's operational range and validate the new FGR model capabilities to higher burnups and transient conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Thermal Conductivity of Graphene Coated Copper Under Uniaxial Tensile Mechanical Strain

Graphene continues to demonstrate promise as a highly effective barrier coating, even at only one atom thick. The thermal properties of this coating are also promising to allow diffusion of heat across the surface, as the isolated graphene is an intrinsically good thermal conductor, yet this and its behavior under mechanical deformation has been less extensively studied. This report demonstrates that the in-plane thermal conductivity and interfacial thermal conductance of graphene coatings on copper are affected by mechanical strain. By inducing a strain in the copper substrate, the Raman-active 2D peak exhibits a change in position as well as a change in laser power dependence as the copper substrate is uniaxially elongated to a maximum of 0.5%. Non-linear trends in thermal conductivity are observed with tensile strain in samples with differing strain transfer rates from the substrate, indicating the close correlation between intrinsic thermal conduction and interfacial properties in atomically thin coatings transferred onto metals.

36 MATERIALS SCIENCE↗

Computational Modeling of Molten Salt Infiltration and Oxidation in Nuclear Graphite

Graphite is utilized as a moderator and reflector in advanced nuclear reactor designs due to its high thermal conductivity, neutron moderation properties, and resistance to radiation damage. However, its longterm performance and reliability are challenged by degradation mechanisms such as molten salt infiltration in molten salt reactors (MSRs) and oxidation in gas-cooled reactors (GCRs). These mechanisms can compromise the structural integrity and operational lifetime of graphite components, necessitating a more detailed assessment of their physical behavior. This report focuses on the development of computational models for molten salt infiltration and oxidation of graphite to aid the design and performance analysis of graphite components. For molten salt infiltration, a computational framework is developed that couples incompressible Navier-Stokes and phase-field model to simulate the penetration of molten salt into graphite?s interconnected pore structure. Initial model verification is performed using two-phase flows in two dimensions, demonstrating the models ability to capture fundamental physical behavior and agree with analytical solution. This framework is then applied to a realistic IG110 nuclear graphite , where a computed tomography extracted pore geometry is used to analyse the infiltration behavior of FLiNaK molten salt. This model provides insights into how the microstructure and other relevant parameters influence the transport pathways of molten salt into graphite, potentially offering a means to rapidly evaluate a graphite grade?s resistance to infiltration. For oxidation, the report details pore-scale mass and heat transport models, describing the diffusion of gases, reaction kinetics, and thermal effects. Additionally, this report highlights inconsistencies in the existing volume-averaged macroscopic model, particularly in upscaling of reaction kinetics and flux terms, and surface to volume transformations. These inconsistencies suggest that current formulations may not accurately capture the experimentally observed graphite oxidation process, highlighting the need for improved model development. This work advances the development of physics-based computational models for graphite degradation, contributing to improved predictive models for next-generation nuclear reactor designs. Future efforts will focus on refining the infiltration model to address non-physical behaviors and enhance its robustness. Additionally, for oxidation, further studies will employ the principles of volume averaging to rigorously derive the upscaled equations, potentially in collaboration with subject matter experts.

Computational Modeling of Molten Salt Infiltration↗

Thermal tomography 3D imaging of additively manufactured metallic structures

Thermal tomography is a computational method for heat diffusion-based imaging of solids, which provides 3D visualization of data from flash thermography measurements. We investigate thermal tomography imaging and nondestructive evaluation of stainless steel and nickel super alloy metallic structures produced with the laser powder bed fusion (LPBF) additive manufacturing (AM) process. Metallic structures produced with LPBF contain defects, and there are limited capabilities to evaluate these structures non-destructively. Thermal tomography reconstruction of 3D apparent spatial effusivity provides information about AM structure geometry and internal material flaws. We study performance of thermal tomography in imaging of metallic structures through COMSOL computer simulations of transient heat transfer and through reconstruction of data obtained from experimental measurements. Thermal tomography reconstructions of structure shape and dimensions are shown for the Inconel 718 AM structure which has variations in the horizontal plane but is uniform along the depth dimension. Reconstruction of internal defects is investigated using a stainless steel 316L specimen with flat bottom hole (FBH) indentations, and the Inconel 718 plate is produced with the LPBF method, which contains imprinted hemispherical shape low density regions containing non-sintered metallic powder. The FBHs have the same sizes as the imprinted defects in the LPBF specimens but offer better imaging contrast. Thermal tomography reconstructions provide visualizations of internal defects and allow for estimation of their sizes and locations. Results of this study demonstrate that thermal tomography can be used for visualization and quality control in AM.

36 MATERIALS SCIENCE↗

Pulsed Thermal Tomography Nondestructive Examination of Additively Manufactured Reactor Materials (Second Annual Progress Report)

Additive manufacturing (AM) is an emerging method for cost-efficient fabrication of nuclear reactor parts. AM of metallic structures for nuclear energy applications is currently based on laser powder bed fusion (LPBF) process, which can introduce internal material flaws, such as pores and anisotropy. Integrity of AM structures needs to be evaluated nondestructively because material flaws could lead to premature failures due to exposure to high temperature, radiation and corrosive environment in a nuclear reactor. Thermal tomography (TT) provides a capability for non-destructive evaluation of sub-surface defects in arbitrary size structures. Thermal tomography is a computational method for heat diffusion-based imaging of solids, which provides 3D visualization of data from flash thermography measurements. We investigate thermal tomography imaging and nondestructive evaluation of stainless steel and nickel super alloy metallic structures produced with laser powder bed fusion (LPBF) additive manufacturing (AM) process. Metallic structures produced with LPBF contain defects, and there are limited capabilities to evaluate these structures non-destructively. Thermal tomography reconstruction of 3D apparent spatial effusivity provides information about AM structure geometry and internal material flaws. We study performance of thermal tomography in imaging of metallic structures through COMSOL computer simulations of transient heat transfer, and through reconstruction of data obtained from experimental measurements. Reconstruction of internal defects is investigated using a stainless steel 316L specimen with flat bottom hole (FBH) indentations, and Inconel 718 plate produced with laser powder bed fusion (LPBF) method, which contains imprinted hemispherical shape low density regions containing non-sintered metallic powder. The FBH’s have the same sizes as the imprinted defects in the LPBF specimens, but offer better imaging contrast. Thermal tomography reconstructions provide visualizations of internal defects, and allow for estimation of their sizes and locations. Detection sensitivity of TT is limited by noises. We investigate separation of signal from noise in thermography images using several machine learning (ML) methods, including new spatio-temporal blind source separation (STBSS) and spatio-temporal sparse dictionary learning (STSDL) methods. Performance of the ML methods is benchmarked using thermography data obtained from imaging stainless steel 316L and Inconel 718 specimens produced LPBF method with imprinted calibrated porosity defects. The ML methods are ranked by F-score and execution runtime. Finally, we investigate TT of AM stainless steel 316L specimen with imprinted internal porosity defects using relatively low-cost, small form factor infrared (IR) camera based on uncooled micro bolometer detector. Sparse coding related K-means singular value decomposition (SVD) machine learning, image processing algorithms are developed to improve quality of TT images through removal of Additive white Gaussian noise without blurring the images. Following initial qualification of an AM component for deployment in a nuclear reactor, a compact TT can also be used for in-service nondestructive evaluation (NDE). With capability to perform in-service NDE of the AM component lifecycle, TT data can be used for development of a component digital twin.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Molecular tail chemistry controls thermal transport in fullerene films

Herein we report on the thermal conductivities of alkyl- and indene-group functionalized fullerene derivative thin films as measured via time domain thermoreflectance and steady state thermoreflectance. The thermal conductivities vary from 0.064 ± 0.007 W m – 1 K – 1 for [6,6]-phenyl C 61 -butyric acid methyl ester (PCBM) to 0.15 ± 0.017 W m – 1 K – 1 for indene- C 60 bisadduct at room temperature and do not exhibit significant temperature dependence from 300 to 375 K. In comparison to the thermal conductivity of PCBM, increasing the length of the alkyl chain, as in the case of [6,6]-phenyl C 61 -butyric acid butyl ester, and [6,6]-phenyl C 61 -butyric acid octyl ester leads to higher thermal conductivities. Likewise, increasing the number of alkyl chains attached to the fullerenes as in the case of bisadduct PCBM leads to a higher thermal conductivity compared to that of PCBM. We present atomistic insights into the role of chemical functionalization on the overall heat transfer in these fullerene derivatives by conducting molecular dynamics simulations and lattice dynamics calculations. The thermal conductivities predicted via our atomistic simulations qualitatively agree with the experimental trends for our fullerene derivatives. Lattice dynamics calculations reveal that one of the main factors dictating the ultralow thermal conductivities in fullerene derivatives is the large reduction in modal diffusivities in the molecular crystals as calculated from the Allen-Feldman model, thus providing an explanation for their largely reduced thermal conductivities as compared to that of C 60 crystals. The low diffusivities result from high degrees of localization of Einstein-like vibrations in fullerene derivatives due to the molecular side chains, providing the ability to dial-in the properties of these low thermal conductivity solids via molecular engineering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

3D–Integrated, Multi–Functional Carbon Fibers for Stable, High–Areal–Capacity Batteries

Increasing lithium-ion batteries' (LIBs) electrode areal capacity can boost energy density and lower manufacturing costs, but faces challenges in manufacturing, rate performance, and cycling stability. A conductive framework made of commercial micro-sized carbon fibers (Cfs) is presented that serves as a host for both the LiNi 0.5 Mn 0.3 Co 0.2 O 2 (NMC 532) cathode and Cfs anode. The Cf framework has multiple functions that offer high electronic conductivity (270 mS cm –1 ), low tortuosity (1.7), low Li + diffusion resistance (22 Ω), and high thermal conductivity (200 W mK –1 ). Additionally, the Cf-integrated electrodes can have an extremely high mass loading of NMC 532 (70 mg cm –2 ) with a theoretical capacity of 14 mAh cm –2 . Thus, the practical full cells assembled with the Cfs-enabled electrodes exhibit an initial areal capacity of 4.1 mAh cm –2 and capacity retention of 90.4% at 500 cycles at a cycling rate of C/3, 1.5 mA cm –2 . Data collected from the operando isothermal microcalorimetry suggest that full cells utilizing the Cf anode experience less heat release from side reactions compared to cells utilizing a conventional graphite anode. Finally, this present approach is scalable and cost-effective and can fabricate practical LIBs that boast high areal capacity, rate performance, and a lengthy cycling lifetime.

25 ENERGY STORAGE↗

Quantitative Morphological Characterization of Carbide Inclusions in Uranium Metal

Uranium carbides (UCs) are prevalent inclusions in U metal that form during melting operations from interactions with crucible walls and the casting chamber atmosphere. Although UCs have been studied extensively since the beginning of U metal foundry operations, there are still unknowns regarding the effects of thermomechanical processing on their sizes and morphology. Here, we present the results of a series of controlled cooling experiments with molten uranium to elucidate the effect of cooling rate on inclusion morphology in a-U. Samples were melted using a vacuum induction melter and manually cooled at rates of 2.5, 1.7, 1.1, 0.8, 0.3 (± 1%) K/s from ~1600 K to <700 K in under 1 hour. Subsequent scanning electron microscopy (SEM) was performed on cross-sections of the samples, revealing a complex mixture of UC morphologies that are indicative of diffusion and growth influenced by the thermal processing of the U matrix. Image analysis using the Morphological Analysis of Materials (MAMA) software showed that UC sizes generally grew larger with slower cooling rates, and the two slowest cooling rates noticeably impacted the inclusion circularity and ellipse aspect ratio. These results indicate that UC morphology is sensitive to short cooling rates (<1 hour) and could therefore be controlled in the production of metallic nuclear fuels. Additionally, inclusion speciation and morphologies could potentially provide forensic clues about processing history of unknown metal samples. Understanding the driving forces involved in UC morphology evolution is beneficial for evaluating metal fuels for next generation nuclear reactors and for identifying signatures for nuclear forensics.

Athon, Matthew T.↗

Effects of Polymer Morphology on Solvent and Catalyst Accessibility during Polyethylene and Polystyrene Autoxidation

Efficient catalytic deconstruction of plastics requires facile solvent and catalyst access to polymer substrates to minimize mass transfer effects. Autoxidation using Co(II) acetate, Mn(II) acetate, and a radical carrier in acetic acid is a promising strategy to deconstruct mixed plastic waste, yet the role of polymer morphology in governing solvent and catalyst accessibility remains poorly understood. Here, in situ simultaneous small- and wide-angle X-ray scattering (SAXS/WAXS), complemented by X-ray fluorescence (XRF) imaging and high-pressure differential scanning calorimetry (DSC), were used to elucidate interactions between acetic acid, a Co/Mn catalyst solution, and semicrystalline polyethylene (PE) and amorphous polystyrene (PS) from room temperature to 160 °C. In PE, acetic acid and catalyst access were confined to amorphous regions and cryomilled particle interfaces at room temperature, while crystalline lamellae remained intact after soaking for up to 34 h. Increasing temperature enabled solvent uptake into PE, followed by solvent-assisted softening above 100 °C, and a modest melting-point depression that removed lamellar transport barriers upon melting. Conversely for PS, acetic acid penetrated the glassy polymer without inducing chain mobility until the glass transition was reached, above which the observed structural changes were consistent with enhanced segmental mobility which enabled bulk penetration. These results suggest that polymer morphology and thermally activated physical transitions arising from diffusion and polymer–solvent interactions can influence whether autoxidation of plastics is transport-limited or kinetically controlled, providing a framework for aligning reaction conditions with reaction outcomes.

09 BIOMASS FUELS↗

Comparison of fixed charge and polarizable models for predicting the structural, thermodynamic, and transport properties of molten alkali chlorides

Results from extensive molecular dynamics simulations of molten LiCl, NaCl, KCl, and RbCl over a wide range of temperatures are reported. Comparison is made between the “Polarizable Ion Model” (PIM) and the non-polarizable “Rigid Ion Model” (RIM). Densities, self-diffusivities, shear viscosities, ionic conductivities, and thermal conductivities are computed and compared with experimental data. In addition, radial distribution functions are computed from ab initio molecular dynamics simulations and compared with the two sets of classical simulations as well as experimental data. Here, the two classical models perform reasonably well at capturing structural and dynamic properties of the four molten alkali chlorides, both qualitatively and often quantitatively. With the singular exception of liquid density, for which the PIM is more accurate than the RIM, there are few clear trends to suggest that one model is more accurate than the other for the four alkali halide systems studied here.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transport Properties of Magnetized High-Energy Density Plasma (Final Project Report)

This report summarizes results of project DE-SC0016159 "Transport Properties of Magnetized High-Energy-Density Plasma,'' which ran from 7/15/2016 - 7/14/2021. The primary objectives of the work, as stated in the original proposal, were to develop a theory that describes transport coefficients in magnetized high energy density plasmas and to test the theory with molecular dynamics (MD) simulations. The science challenge is that strong ion coupling, strong magnetization of electrons, and partial degeneracy of electrons, cause the system to be in a regime that is not well described by current theory. The particular processes that were to be investigated include: ion stopping power, diffusion, electron-ion temperature relaxation, thermal conduction and viscosity. All of the primary research objectives were accomplished during the course of this work. In addition, some unexpected results were discovered along the way that led to new and productive research directions. This work resulted in 17 publications in well-respected peer-reviewed journals (primarily Physics of Plasmas and Physical Review E), and 2 more are being prepared for publication. A few of these were selected as Editor's Picks and one was published as a Rapid Communication. Highlights of the research results include: Transport phase space: The first identification of the parameter space that defines fundamental transport regimes in terms of the Coulomb coupling and magnetization parameters. Mean force kinetic theory: Systematic derivation of a kinetic theory for strongly coupled plasmas based on a new expansion parameter of the BBGKY hierarchy. This provided the derivation of a theory that we had previously posed phenomenologically, and also revealed a new term that captures the equation of state properties at all coupling strengths. Transverse friction force: Discovery that the friction force on a test charge in a strongly magnetized plasma includes a component that is perpendicular to the motion of the test charge in the plane defined by the velocity and magnetic field vectors. Gyrofriction force: Discovery that the friction force on a test charge in a plasma that is both strongly magnetized and strongly coupled includes a component of the force in the direction of the Lorentz force. Kinetic theory for warm dense matter: Extension of the mean force kinetic theory to include partial electron degeneracy, so that it applies to warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparison of Zirconium Redistribution in BISON EBR-II Models Using FIPD and IMIS Databases with Experimental Post Irradiation Examination

Metallic fuels have seen increased interest for future sodium fast reactors due to their material properties: high thermal conductivities and advantageous neutronic properties allow for greater fission densities. One drawback to typical metallic fuels is zirconium redistribution, which impacts this advantageous material and its neutronic properties. Unfortunately, the processes behind zirconium migration behavior are understood using first principles, so before these fuels are implemented in future fast reactors, characterization and fuel qualification regimes must be completed. These activities can be supported through the use of robust modeling using the most accurate empirical models currently available to fuel researchers around the world. The tool that allows researchers to model this complex coupled thermo-mechanical behavior and nuclear properties is BISON. Additionally, BISON model parameters need to be compared against PIE measurements. The current work utilizes two fuel pins from EBR-II experiment X441 to optimize various model parameters, including porosity correction factor, thermal conductivity, phase transition temperature, and diffusion coefficient multipliers, before implementing the final model for seven fuel pins with differing characteristics. To properly evaluate the BISON simulations, the results are compared to PIE metallography data for each fuel pin, to ensure the zirconium redistribution is properly reflected in the simulation results. Six out of seven analyzed fuel pins demonstrate good agreement between the metallography images and BISON results, showing alignment of the Zr-rich, Zr-depleted, and moderately Zr-enriched zones at various axial heights along the fuel pins. Further work is needed to refine the model parameters for general pin use.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Superionic-like diffusion in yttrium dihydride

For the next-generation high temperature microreactors, yttrium dihydride (YH 2 ) is an attractive solid state neutron moderator. Despite a number of recent investigations, the mechanism of hydrogen transport remains poorly understood. Experimental evaluations of diffusivity are inconclusive with large variations in diffusivities and activation energies. In this work, we perform ab initio molecular dynamics (AIMD) simulations on YH 2 for temperatures spanning 300 K to 1200 K. Our main finding is that YH 2 shows a superionic-like behavior with hydrogen atoms hopping from one native site to another above a characteristic temperature of 800 K. This correlated motion results in quasi-one-dimensional string-like displacements that enable the hydrogen atoms to diffuse rapidly. We confirm that the octahedral sites are mostly unoccupied, although channeling through them is the most favored pathway between lattice hops above 800 K. At the highest temperature of 1200 K, the string relaxation time is merely of the order of a few picoseconds, which indicates a liquid-like diffusive behavior. Based on the formation of spontaneous thermal vacancies, an order-disorder crossover temperature T α ~ 800 K is established for YH 2 with an activation energy of 0.83 eV for hydrogen diffusion in the superionic-like state.

Superionic-like Diffusion↗

Onset of anharmonicity and thermal conductivity in SnSe

The anharmonicity in SnSe is investigated through the analysis of moments of 119Sn nuclear resonant inelastic x-ray scattering and ab initio molecular dynamics calculations. Experimental evidences show that the anharmonic behavior started around 300 K, substantially lower than the usually suggested structural transition at 800 K. Both theory and experiments reveal substantial lifetime broadening and frequency renormalization of the optical phonons. Additionally, thermal conductivities calculated from the temporal energy moment using the Einstein diffusion equation are in good agreement with previous experiments. The abrupt increase of the thermal power near 800 K is driven by an electronic factor and not by the enhanced anharmonicity due to structural change.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A New Route to Quantum-Scale Structures through a Novel Enhanced Germanium Diffusion Mechanism

This project sought to develop a fundamental understanding of the mechanisms underlying a newly observed enhanced germanium (Ge) diffusion process in silicon germanium (SiGe) semiconductor nanostructures during thermal oxidation. Using a combination of oxidationdiffusion experiments, high resolution imaging, and theoretical modeling, a model for the enhanced Ge diffusion mechanism was proposed. Additionally, a nanofabrication approach utilizing this enhanced Ge diffusion mechanism was shown to be applicable to arbitrary 3D shapes, leading to the fabrication of stacked silicon quantum dots embedded in SiGe nanopillars. A new wet etch-based method for preparing 3D nanostructures for highresolution imaging free of obscuring material or damage was also developed. These results enable a new method for the controlled and scalable fabrication of on-chip silicon nanostructures with sub-10 nm dimensions needed for next generation microelectronics, including low energy electronics, quantum computing, sensors, and integrated photonics.

36 MATERIALS SCIENCE↗

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗

Investigating the Linear Thermal Expansion of Additively Manufactured Multi-Material Joining between Invar and Steel

This work investigated the linear thermal expansion properties of a multi-material specimen fabricated with Invar M93 and A36 steel. A sequence of tests was performed to investigate the viability of additively manufactured Invar M93 for lowering the coefficient of thermal expansion (CTE) in multi-material part tooling. Invar beads were additively manufactured on a steel base plate using a fiber laser system, and samples were taken from the steel, Invar, and the interface between the two materials. The CTE of the samples was measured between 40 °C and 150 °C using a thermomechanical analyzer, and the elemental composition was studied with energy dispersive X-ray spectroscopy. The CTE of samples taken from the steel and the interface remained comparable to that of A36 steel; however, deviations between the thermal expansion values were prevalent due to element diffusion in and around the heat-affected zone. The CTEs measured from the Invar bead were lower than those from the other sections with the largest and smallest thermal expansion values being 10.40 μm/m-K and 2.09 μm/m-K. In each of the sections, the largest CTE was measured from samples taken from the end of the weld beads. An additional test was performed to measure the aggregate expansion of multi-material tools. Invar beads were welded on an A36 steel plate. The invar was machined, and the sample was heated in an oven from 40 °C and 160 °C. Strain gauges were placed on the surface of the part and were used to analyze how the combined thermal expansions of the invar and steel would affect the thermal expansion on the surface of a tool. There were small deviations between the expansion values measured by gauges placed in different orientations, and the elongation of the sample was greatest along the dimension containing a larger percentage of steel. On average, the expansion of the machined Invar surface was 42% less than the expansion of the steel surface. The results of this work demonstrate that additively manufactured Invar can be utilized to decrease the CTE for multi-material part tooling.

36 MATERIALS SCIENCE↗