Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity coefficient”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Effective Fragment Potential-Based Molecular Dynamics Studies of Diffusion in Acetone and Hexane

To facilitate more reliable descriptions of transport properties in liquids, molecular dynamics (MD) simulations are performed based on the effective fragment potential (EFP) method derived from first-principles quantum mechanics (in contrast to MD based upon empirically fitted potentials). The EFP method describes molecular interactions in terms of Coulomb, polarization/induction, dispersion, exchange-repulsion, and charge-transfer interactions. The EFP MD simulations described in this paper, performed on hexane and acetone, are able to track the mean-square displacement of molecules for sufficient time to reliably extract translational diffusion coefficients. Finally, the results reported here are in reasonable agreement with experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Kinetics and transport of hydrogen in graphite at high temperature and the effects of oxidation, irradiation and isotopics

The kinetics of uptake and desorption impact the performance of graphite as a vector for tritium in high-temperature fission reactors and in the blanket of fusion reactors. Graphite components in these reactors are exposed to temperatures > 500 °C and H 2 partial pressures of few Pa and desorption temperatures are limited to < 1600 °C; limited data is available at these conditions. Here we review the mechanisms for uptake in, transport and desorption of hydrogen from graphite at high temperature, compiling data on uptake rates, diffusion coefficients and activation energies and providing a discussion of the impact of irradiation, pre-oxidation and isotope. At FHR conditions, trapping impacts uptake rates, leading to a reduction in apparent diffusivity by 35 to 80% compared to higher partial-pressure uptake. Timelines for desorption are not clearly defined; extrapolating from available data, at 1150 °C desorbing 80% of tritium uptaken at FHR conditions may take from 100 to 10,000 h.

36 MATERIALS SCIENCE↗
DOE OSTI2020

STRUCTURE AND PROPERTIES STUDY ON ENERGY MATERIALS: THERMOELECTRIC MATERIAL TETRAHEDRITE AND LITHIUM ION CONDUCTOR

Development of efficient energy materials is critical in order to ease the energy demand and reduce our dependence of fossil fuel. Thermoelectric materials are promising due to their capability of generating electrical power by recovering waste heat. The performance of thermoelectric materials is quantified by a dimensionless figure of merit zT, which depends on their properties such as electrical conductivity, Seebeck coefficient and thermal conductivity. Tetrahedrites, a copper antimony sulfosalt mineral, typified by Cu 12-x MxSb 4 S 13 , where M is a transition metal element such as Ni, Zn, Fe or Mn, have great potential for thermoelectric application due to their relatively high zT (close to 1 at 700 K), earth-abundance, environmental friendliness, favorable electrical properties, and most importantly intrinsic low lattice thermal conductivity (less than 1 W m -1 K -1 ) in wide temperature. In addition to energy recovery, reliable energy storage devices are also emerging to relieve the energy demand and improve the efficiency of consuming energy resources. Lithium-ion batteries are known to be reliable and successful electrochemical energy storage devices and appliable in various aspects, including laptops, smartphones and electrical vehicles. Lithium phosphorous oxynitride (LiPON) are widely used as thin-film solid-state electrolytes in Li-ion battery, which is the only demonstrated solid-state electrolyte that is quite stable in direct contact with Li metal at potentials from 0-5 V. However, the structure of LiPON, the effects of N doping, and the origin of its good electrochemical stability remains inconclusive. In this thesis, reliable modeling techniques accompanied with experimental tools, are applied to study the thermoelectric material tetrahedrite and the ionic conductor LiPON, in order to study their structural and dynamical properties. Accurate and efficient density-functional theory (DFT) and density-functional tight-binding (DFTB) methods, combined with molecular dynamics (MD) simulations are utilized in order to investigate the structures and properties of these energy materials. The incoherent and coherent atomic dynamics study of tetrahedrite Cu 10.5 NiZn 0.5 Sb 4 S 13 provides the origin of softening upon cooling by investigate the motion of Cu12e at different temperatures. The dynamic structure factors in the longitudinal and transverse direction will also be discussed. The Cu movement of Cu-rich tetrahedrite Cu 14 Sb 4 S 13 is revealed by Cu self-diffusivity, nuclear density map and “nudged elastic band” (NEB). Moreover, we investigate the effect of simulation cell size and basis sets on the DFT-based MD simulation results using tetrahedrite Cu1 0 Zn 2 Sb 4 S 13 thermoelectric as a model material, showing the advantage of larger cell by accessing smaller Q range. In addition, the low-temperature structural properties of Cu 12 Sb 4 S 13 is measured by neutron diffraction, which indicates that no cubic to tetragonal transition occurs at metal-semiconductor transition (MST) temperature. Thermoelectric properties such as Seebeck coefficient, electrical resistivity and electrical thermal conductivity will also be investigated. DFTB method is implemented to study the structure and transport properties of Li 3 PO 4 and LiPON, while the exploration of N doping effect is included. Lastly, the LiPON/Li interphase will be revealed in order to study the origin of electrochemical stability.

Li, Junchao↗
DOE OSTI2021

A non‐intrusive domain‐decomposition model reduction method for linear steady‐state partial differential equations with random coefficients

Abstract Domain decomposition methods have been proved to be an effective strategy to reduce the dimension of parametric partial differential equations (PDEs). However, existing domain decomposition methods for parametric PDEs are usually intrusive, which means domain decomposition based solvers need to be implemented from scratch for each target parametric PDE. To address this issue, we develop a new non‐intrusive domain‐decomposition model reduction method for linear steady‐state PDEs with random‐field coefficients. As a variant of our previous work by Mu and Zhang, the new method only needs access to the final linear system, that is, the global stiffness matrix and the right hand side, of a deterministic PDE solver, in order to build a domain‐decomposition‐based reduced model without intrusive implementation from scratch. The key idea is to remove the interface condition between sub‐domains and rely on the correlation between columns of the linear system to couple the sub‐domains. The non‐intrusive feature enables the applicability of the proposed method to a broader class of uncertainty quantification problems, where many legacy codes/solvers can be fully reused by our method. Two numerical examples including diffusion equations with random diffusivity and convection‐dominated transport with random velocity, are provided to demonstrate the effectiveness and efficiency of our method.

Zhang, Guannan↗
DOE OSTI2023

Irradiation Effects on Stability of δ-UZr2 phase in U-50 wt% Zr Alloy

U-50wt%Zr is a candidate metallic nuclear fuel with potential application in light water reactors due to its excellent thermal properties and high radiation tolerance. The Zr-rich UZr fuels possess greater swelling resistance and fission gas release characteristics compared with U-rich UZr fuels. In this current study, the δ-phase U-50wt%Zr is proton irradiated at room temperature to 1 displacement per atom (dpa) to provide insights on phase stability under irradiation conditions. High resolution characterization of Transmission Electron Microscopy (TEM) and Atom Probe Tomography (APT) characterization techniques are used to elucidate microstructural changes due to irradiation. TEM and APT results show highly oriented bcc β-Zr-rich platelet precipitates nucleating adjacent to α-U phases inside the UZr 2 matrix. Formation of this platelet morphology is characteristic of Widmänstatten structure which can be attributed to a variety of factors such as differences in thermal expansion coefficient between the phases, grain size, alloy composition, and cooling rate. The phases present are distinctly different than those observed through in situ annealing, but irradiation accelerates diffusion and phase separation kinetics. These microstructural changes in U-50wt%Zr are different from those achieved by pure thermodynamic or high temperature heavy ion irradiation experiments. Finally, our work, together with previous ones, highlight the necessity to study the U-Zr phase diagram under non-equilibrium thermodynamics conditions to support this material's deployment as a viable nuclear fuel form.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Impact of Ion–Ion Correlated Motion on Salt Transport in Solvated Ion Exchange Membranes

The influence of dynamical ion-ion correlations and ion pairing on salt transport in ion exchange membranes remain poorly understood. In this study, we use the framework of Onsager transport coefficients within atomistic molecular dynamics sim- ulations to study the impact of ion-ion correlated motion on salt transport in hydrated polystyrene sulfonate membranes and compare with the results from aqueous salt so- lutions. At sufficiently high salt concentrations, cation-anion dynamical correlations exert a significant influence on both salt diffusivities and conductivities. Anion-anion distinct correlations, arising from the imbalance between the concentration of free (mo- bile) cations and anions, and the retarding effect of the fixed charge groups on cations, proves to be an additional important feature for polymer membranes. Furthermore, our results demonstrate that dynamical correlations should become an important consideration in experimental measurements of salt diffusivities and conductivities for non-dilute salt solutions in polymer membranes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Interfacial Effects on Transport Coefficient Measurements in Li-ion Battery Electrolytes

Development of Li + -containing electrolytes with improved transport properties requires reliable, reproducible, and ideally low volume techniques to rigorously understand ion-transport with varying composition. Precisely measuring the complete set of transport coefficients in liquid electrolytes under battery-relevant operating conditions is difficult and the reliability of these methods are sparsely described in electrolyte transport literature. In this work, we apply a potentiostatic polarization-based transport characterization approach typically used for polymer electrolytes to liquid electrolyte systems in an attempt to fully measure all transport coefficients (conductivity, total salt diffusion coefficient, thermodynamic factor and transference number) for the model system of LiPF 6 in an ethylene carbonate—ethyl methyl carbonate (EC:EMC) mixture. Using systematic timescale and statistical analyses, we find that transport coefficients measured using potentiostatic polarization of Li-Li symmetric cells exhibit strong correlation to Li electrode interfacial resistance, indicating that such methods are probing both bulk and interfacial phenomena. This reveals a major roadblock in characterizing electrolyte systems where the interfacial resistance is significantly larger than ohmic electrolyte resistance. As a result, we find that methods that rely on potentiostatic Li metal stripping/plating do not readily result in reliable liquid electrolyte transport coefficients, unlike similar methods for solid polymer electrolytes, where interfacial resistances are typically smaller than electrolyte resistances at the elevated temperatures typically of interest for such electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2024

Compatibility of LaFe13−x−yMnxSiyH1.6 and Eutectic Liquid GaInSn Alloy

The heat transfer rate of magnetocaloric regenerators is a topic of extensive research and the cyclability of these regenerators is critical to the operation of systems with a high coefficient of performance (e.g., potentially >22, significantly higher than typical vapor compression cooling technologies). To enable a high operating frequency that will result in a high specific cooling power, the heat transfer fluid should have high thermal conductivity and lower specific heat, i.e., higher thermal diffusivity. Eutectic metal alloys possess these qualities, such as gallium–indium–tin (Galinstan), whose thermal diffusivity has been found to be approximately an order of magnitude higher than water. For this study, the effects of eutectic liquid Galinstan exposure on the phase stability of LaFe13−x−yMnxSiyH1.6 magnetocaloric powders in an active magnetic regenerator device were investigated. The powders were characterized before and after exposure to Galinstan using X-ray diffraction, in which the phases were determined using the Rietveld refinement technique and X-ray fluorescence. It was found that after Galinstan exposure, hydrogen containing phases were present in the powder, suggesting that the hydrogen was lost from the magnetocaloric phase. The magnetocaloric phase degradation indicates that the powder was incompatible with the Galinstan metal in an environment with moisture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Oxygen diffusion in garnet: Experimental calibration and implications for timescales of metamorphic processes and retention of primary O isotopic signatures

Knowledge of oxygen diffusion in garnet is crucial for a correct interpretation of oxygen isotope signatures in natural samples. A series of experiments was undertaken to determine the diffusivity of oxygen in garnet, which remains poorly constrained. The first suite included high-pressure (HP), nominally dry experiments performed in piston-cylinder apparatus at: (1) T = 1050–1600 °C and P = 1.5 GPa and (2) T = 1500 °C and P = 2.5 GPa using yttrium aluminum garnet (YAG; Y3Al5O12) cubes. Second, HP H2O-saturated experiments were conducted at T = 900 °C and P = 1.0–1.5 GPa, wherein YAG crystals were packed into a YAG + Corundum powder, along with 18O-enriched H2O. Third, 1 atm experiments with YAG cubes were performed in a gas-mixing furnace at T = 1500–1600 °C under Ar flux. Finally, an experiment at T = 900 °C and P = 1.0 GPa was done using a pyrope cube embedded into pyrope powder and 18O-enriched H2O. Experiments using grossular were not successful. Profiles of 18O/(18O+16O) in the experimental charges were analyzed with three different secondary ion mass spectrometers (SIMS): sensitive high-resolution ion microprobe (SHRIMP II and SI), CAMECA IMS-1280, and NanoSIMS. Considering only the measured length of 18O diffusion profiles, similar results were obtained for YAG and pyrope annealed at 900 °C, suggesting limited effects of chemical composition on oxygen diffusivity. However, in both garnet types, several profiles deviate from the error function geometry, suggesting that the behavior of O in garnet cannot be fully described as simple concentration-independent diffusion, certainly in YAG and likely in natural pyrope as well. The experimental results are better described by invoking O diffusion via two distinct pathways with an inter-site reaction allowing O to move between these pathways. Modeling this process yields two diffusion coefficients (D values) for O, one of which is approximately two orders of magnitude higher than the other. Taken together, Arrhenius relationships are:logDm2s-1=-7.2(±1.3)+(-321(±32)kJmol-12.303RT) for the slow pathway, andlogDm2s-1=-5.4(±0.7)+(-321(±20)kJmol-12.303RT) for the fast pathway. We interpret the two pathways as representing diffusion following vacancy and inter-stitial mechanisms, respectively. Regardless, our new data suggest that the slow mechanism is prevalent in garnet with natural compositions, and thus is likely to control the retentivity of oxygen isotopic signatures in natural samples. The diffusivity of oxygen is similar to Fe-Mn diffusivity in garnet at 1000–1100 °C and Ca diffusivity at 850 °C. However, the activation energy for O diffusion is larger, leading to lower diffusivities at P-T conditions characterizing crustal metamorphism. Therefore, original O isotopic signatures can be retained in garnets showing major element zoning partially re-equilibrated by diffusion, with the uncertainty caveat of extrapolating the experimental data to lower temperature conditions.

Geochemistry & Geophysics↗
DOE OSTI2023

Intrinsic Berry curvature driven anomalous Nernst thermopower in the semimetallic Heusler alloy CoFeVSb

Understanding of spin-heat coupling mechanisms and magnetothermoelectric phenomena, including the anomalous Nernst effect (ANE), in emergent quaternary Heusler alloys is of practical importance for applications in thermal management and energy harvesting. Here, we demonstrate an intrinsic Berry curvature mediated anomalous Nernst thermopower in CoFeVSb, which orders magnetically at high temperature ( T C ≈ 850 K ) with a large saturation magnetization of ≈ 2.2 μ B / f . u . at room temperature. We show that the electron-electron elastic and electron-magnon inelastic scattering dominate longitudinal electrical transport at low temperatures ( T ≤ 50 K ), whereas the electron-phonon and electron-magnon scatterings govern it at higher T . The longitudinal thermopower is resulted mainly from the diffusive contribution with a very large longitudinal Seebeck coefficient ( 42 μ V K - 1 at 395 K). The value of the anomalous Nernst coefficient ( S ANE ) for CoFeVSb at room temperature is 0.039 μ V K - 1 which is higher than the compressively strained SrRu O 3 film ( 0.03 μ V K - 1 ) as well as the spin gapless semiconductor CoFeCrGa ( 0.018 μ V K - 1 ). On lowering T , both the ordinary Nernst coefficient and carrier mobility increase but an opposite trend is found for S ANE . Our ab initio simulations reveal the topological semimetallic nature of CoFeVSb with a pair of Weyl points. These Weyl crossings result in a significant contribution to the Berry curvature, leading to an intrinsic anomalous Hall conductivity ( σ x y AHE ) of ≈ 85 S/cm, which matches well with experiment (77 S/cm at 2 K). Our experimental findings and ab initio calculations support the dominance of the intrinsic Berry curvature in the observed ANE. The ratio of σ x y AHE to the transverse anomalous thermoelectric conductivity ( α x y ANE ) shows an increasing trend with T attaining a sizable fraction of k B e ( ≈ 0.35 k B e ) at room temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of particle shape on stratification in drying films of binary colloidal mixtures

The role of particle shape in evaporation-induced auto-stratification in polydisperse colloidal suspensions is explored with molecular dynamics simulations of mixtures of spheres and aspherical particles. A unified framework based on the competition between diffusion and diffusiophoresis is proposed to understand the effects of shape and size dispersity. In general, particles diffusing more slowly (e.g., larger particles) tend to accumulate more strongly at the evaporation front. However, larger particles have larger surface areas and therefore greater diffusiophoretic mobility. Hence, they are more likely to be driven away from the evaporation front via diffusiophoresis. For a rapidly dried bidisperse suspension containing small and large spheres, the competition leads to “small-on-top” stratification. Here, we employ a computational model in which the diffusion coefficient is inversely proportional to particle mass. For a mixture of spheres and aspherical particles with similar mass, the diffusion contrast is reduced, and the spheres are always enriched at the evaporation front as they have the smallest surface area for a given mass and, therefore, the lowest diffusiophoretic mobility. Furthermore, for a mixture of solid and hollow spheres that have the same outer radius and thus the same surface area, the diffusiophoretic contrast is suppressed, and the system is dominated by diffusion. Consequently, the solid spheres, which have a larger mass and diffuse more slowly, accumulate on top of the hollow spheres. Finally, for a mixture of thin disks and long rods that differ significantly in shape but have similar mass and surface area, both diffusion and diffusiophoresis contrasts are suppressed, and the mixture does not stratify.

Classical molecular dynamic simulations↗
DOE OSTI2025

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A quantitative model of charge injection by ruthenium chromophores connecting femtosecond to continuous irradiance conditions

A kinetic framework for the ultrafast photophysics of tris(2,2-bipyridine)ruthenium(II) phosphonated and methyl-phosphonated derivatives is used as a basis for modeling charge injection by ruthenium dyes into a semiconductor substrate. By including the effects of light scattering, dye diffusion, and adsorption kinetics during sample preparation and the optical response of oxidized dyes, quantitative agreement with multiple transient absorption datasets is achieved on timescales spanning femtoseconds to nanoseconds. In particular, quantitative agreement with important spectroscopic handles—the decay of an excited state absorption signal component associated with charge injection in the UV region of the spectrum and the dynamical redshift of a ~500 nm isosbestic point—validates our kinetic model. Pseudo-first-order rate coefficients for charge injection are estimated in this work, with an order of magnitude ranging from 10 11 to 10 12 s –1 . The model makes the minimalist assumption that all excited states of a particular dye have the same charge injection coefficient, an assumption that would benefit from additional theoretical and experimental exploration. Here, we have adapted this kinetic model to predict charge injection under continuous solar irradiation and find that as many as 68 electron transfer events per dye per second take place, significantly more than prior estimates in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Thermal properties of HfN-MoW surrogate cermet fuel for nuclear thermal propulsion

Here, the thermal properties of MoW-HfN, a surrogate cermet for MoW-UN nuclear thermal propulsion fuel, were characterized over a wide range of elevated temperatures. Thermal diffusivity, coefficient of thermal expansion (CTE), and heat capacity were measured. Optical and scanning electron microscopy were performed to characterize the microstructure and draw structure–property correlations. The thermal diffusivity was obtained using the laser flash method. Diffusivity values ranged from about 0.18 cm 2 /s at 200°C to 0.15 cm 2 /s at 1800°C. The CTE was measured using push-rod dilatometry up to 1600°C, giving values between 6.0 and 9.0 μm/m. A scientific rationalization of the effective material properties is made using the rule-of-mixtures and other effective properties models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Sensitivity of Global Ocean Deoxygenation to Vertical and Isopycnal Mixing in an Ocean Biogeochemistry Model

Large-scale loss of oxygen under global warming is termed “ocean deoxygenation” and is caused by the imbalance between physical supply and biological consumption of oxygen in the ocean interior. Significant progress has been made in the theoretical understanding of ocean deoxygenation; however, many questions remain unresolved. The oxygen change in the tropical thermocline is poorly understood, with diverging projections among different models. Physical oxygen supply is controlled by a suite of processes that transport oxygen-rich surface waters into the interior ocean, which is expected to weaken due to increasing stratification under global warming. Using a numerical model and a series of sensitivity experiments, the role of ocean mixing is examined in terms of effects on the mean state and the response to a transient warming. Both vertical and horizontal (isopycnal) mixing coefficients are systematically varied over a wide range, and the resulting oxygen distributions in equilibrated and transient simulations are examined. The spatial patterns of oxygen loss are sensitive to both vertical and isopycnal mixing, and the sign of tropical oxygen trend under climate warming can reverse depending on the choice of mixing parameters. An elevated level of isopycnal mixing disrupts the vertical advective-diffusive balance of the tropical thermocline, increasing the mean state oxygen as well as the magnitude of the transient oxygen decline. These results provide first-order explanations for the diverging behaviors of simulated tropical oxygen with respect to ocean mixing parameters.

58 GEOSCIENCES↗
DOE OSTI2026

Measurements of the electron energy distribution function in partially magnetized low temperature plasmas

While Langmuir probes (LPs) are relatively simple and inexpensive plasma diagnostics for the electron density, temperature, and the electron energy distribution function (EEDF), the interpretation of the measured current–voltage (I–V) characteristic is complicated considerably by the presence of a magnetic field. In regimes where the electron gyroradius is comparable to the probe radius, the electron flux to the probe surface is retarded by reduced mobility across field lines and can no longer be described by a thermal model. Predicting the current collected by the probe in these regimes requires accurate estimates of the plasma diffusion coefficients, which are usually difficult to obtain. In this work, we measure electron energy distribution functions in E×B plasmas with magnetized electrons and non-magnetized ions in argon and krypton gases, using both a LP and laser Thomson scattering (LTS) at various magnetic fields. Using the LTS measurements to provide a robust benchmark for comparison, we compare existing theories describing the flux to the probe under magnetized conditions. We find that even when the electron gyroradius associated with the effective electron temperature is small compared to the probe radius, the EEDF computed using classical probe theory is still robust at energies higher than the energy at which the gyroradius becomes larger than the probe size. For plasmas that are not strongly non-Maxwellian, we formulate a method to extract robust density and temperature measurements using physics-informed fitting techniques to analyze EEDFs computed using classical theory.

Devin, E. G. [Princeton Plasma Physics Laboratory ↗
DOE OSTI2023

Description of the multinucleon transfer mechanism for Ca 48 + Pu 244 and Kr 86 + Pt 198 reactions in a quantal transport approach

Multinucleon transfer (MNT) reactions involving heavy projectile and target combinations stand as a promising method for synthesizing new neutron-rich exotic nuclei, which may not be possible using hot or cold fusion reactions or fragmentation. Exploring the mechanisms behind MNT reactions is essential and it requires a comprehensive theoretical framework that can explain the physical observables in these reactions. This work aims to show that the quantal diffusion approach based on the stochastic mean-field (SMF) theory is capable of explaining the reaction dynamics observed in MNT reactions. Primary product mass distributions in 48 Ca + 244 Pu reaction at E c.m. = 203.2 MeV and 86 Kr + 198 Pt reaction at E c.m. = 324.2 MeV are calculated and compared with the available experimental data. In this work, we utilize the time-dependent Hartree-Fock (TDHF) calculations to analyze the mean-field reaction dynamics computationally in the reactions 48 Ca + 244 Pu and 86 Kr + 198 Pt for a broad range of initial angular momenta. Quantal transport description based on the SMF approach is used to calculate quantal diffusion coefficients and mass variances in 48 Ca + 244 Pu and 86 Kr + 198 Pt systems. The primary products arising from quasifission reactions are described by joint probability distribution in the SMF approach and those arising from fusion-fission are estimated by using the statistical deexcitation code gemini + +. Mean values of charge and mass numbers, scattering angles of the primary reaction products, and the total kinetic energies after the collision are calculated within the TDHF framework for a broad range of initial angular momenta. Throughout all the collisions, drift toward the mass symmetry and large mass dispersion associated with this drift are observed. Here, the calculated primary fragment and mass distributions using the SMF approach successfully explain experimental observations for the 48 Ca + 244 Pu and 86 Kr + 198 Pt systems. The primary mass distributions, mean values of binary products, and mass dispersions are determined and results are compared with the available experimental data. The observed agreement between the experimental data and SMF results highlights the effectiveness of the quantal diffusion mechanism based on the SMF approach, which does not include any adjustable parameters other than standard parameters of Skyrme energy density functional.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Rapid Bayesian High Entropy Alloy Designs Fabricated via Wire Arc Additive Manufacturing

Purpose: This project seeks to demonstrate a new high-throughput (rapid) alloy design technique applied to creating new high entropy alloys (HEAs) for extreme environments. High entropy alloys shift the design paradigm from being focused on a single principal element (e.g. nickel-based alloys) to target alloys that include high atomic fractions (X >10%) of multiple elements. These HEA materials can exhibit sluggish diffusion and enhanced corrosion resistance, ideal for potential applications in advanced ultra supercritical (A-USC) steam cycles for power generation. Scope: The addition of multiple elements in high atomic fractions creates an enormous design space that cannot easily be investigated by traditional material design strategies such as designed of experiments (DOE). This project utilizes a Bayesian machine learning algorithm that has been modified to work with calculation of phase diagrams (CALPHAD) software. This Bayesian algorithm reduces manual inputs and increase the likelihood of achieving an optimal solution. Compositional inputs to this algorithm will be assessed using existing material property models for high temperature strength and corrosion resistance. The target for alloy performance will be a 15% (~100 ⁰C) increase in allowable service temperature beyond heat-resistant stainless steels while maintaining or improving alloy cost and corrosion resistance. Haynes 230 was selected as a baseline, which is 57 wt% Ni with 22 wt% Cr 14 wt% W, and 2 wt% Mo as solid solution strengtheners. In addition to rapid design via Bayesian machine learning, the alloys were rapidly fabricated using a multi-wire arc additive manufacturing (mWAAM) technique which allows for precise control of alloy composition and assessing of alloy design “windows” to study composition effects. Build speeds for wire-arc additive processes are among the highest for additive technologies enabling rapid and reliable sample fabrication when compared to conventional methods such as arc button melting. The mWAAM samples will be rapidly characterized via instrumented indentation for room temperature modulus and strength and for elevated temperature strength via hot hardness tests. After being screened with hardness testing, potential alloys will be further evaluated with conventional microscopy techniques including scanning electron microscopy (SEM) and transmission electron microscopy (TEM) to assess agreement with modeling results. The most promising compositions will also be evaluated by printing full sized tensile specimens for mechanical behavior tests at elevated temperatures. Results: Bayesian machine learning of a single performance function was initially used to optimize five performance metrics: 1) single phase stability, 2) yield strength, 3) creep resistance (low diffusion coefficient), 4) freezing range (weldability), and 5) material cost. The single performance function was suboptimal as assumptions had to be made about the results while formulating the optimization. A goal-oriented Bayesian optimization strategy (Hanaoka, 2021) was implemented with CALPHAD for use with the five metrics above. This multi-objective Bayesian optimization (MOBO) enabled the design of NiCrCoFe alloys with V and W additions. A base composition of NiCoCr was selected as Ni provides a stable FCC matrix, Cr aids corrosion/oxidation resistance, and Co is a solid-solutions strengthener that also improves creep by increasing the activation energy. Fe helps reduce diffusion coefficients and cost. Finally, V and W were selected for their reasonable solubility and high atomic misfit to aid in solid solution strengthening. Cracking of the mWAAM specimens was an early issue, and the Easton solidification cracking model (Easton et al., 2014a) was selected for addition to the MOBO function. High performing alloys fabricated by mWAAM included Ni 28 Cr 25 Co 26 Fe 15 V 8 and Ni 62 Cr 18 Co 1 Fe 3 W 15 . It was observed that even after adapting the mWAAM process for W, the W did not fully dissolve. To fully evaluate the Ni 62 Cr 18 Co 1 Fe 3 W 15 composition, a cored wire (80-20 NiCr sheath/powder core) was manufactured and printed via WAAM, and HIP’ing was utilized to homogenize and densify the printed alloy. The V and W alloys produced met metrics 1 (solid solution), 4 (solidification cracking), and 5 (cost). However, an unmodeled mechanism of thermal stress cracking was identified in the WAAM produced materials, perhaps exacerbated by the lack of grain boundary strengthening elements (B, C). Conclusions & Recommendations: A high-throughput (rapid) alloy design technique was applied to designing and manufacturing new high entropy alloys (HEAs) for extreme environments utilizing MOBO and mWAAM. The developed process was rapid and effective in addressing the mechanisms included in the model. The lack of grain boundary strengthening element additions (e.g., B, C) was a simplification that likely produced thermal stress cracking that turned into a large part of the investigation. Additions on the order of 0.005 wt% B and 0.05 wt% C likely would have minimized thermal stress grain boundary cracking. Overall, the high throughput design strategy is promising for rapid design of metrics-driven alloys for advanced ultra supercritical (A-USC) steam cycles for power generation. The MOBO and mWAAM process could be commercialized to accelerate metrics-driven alloy design. In addition, the cored-wire process utilized for scale-up is a promising high-volume process for WAAM alloy development and scale-up.

36 MATERIALS SCIENCE↗