Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “microstructure simulation”

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 433 records · Page 24

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗

Modeling Ni redistribution in the hydrogen electrode of solid oxide cells through Ni(OH) 2 diffusion and Ni-YSZ wettability change

Ni redistribution in the hydrogen electrodes of solid oxide cells is an important degradation mechanism. Its driving force is still under debate. This work focuses on the interplay between Ni(OH) 2 diffusion and Ni-YSZ wettability change. With a Ni(OH) 2 diffusion model developed in our previous work, we further employ three models for the Ni-YSZ contact angle to quantify the Ni-YSZ wettability change. The microstructure evolutions in a reconstructed Ni-YSZ electrode are simulated under selected experimental operating conditions. It is shown that the phenomenological model can capture the driving force of Ni spreading/detachment on YSZ surface and the Ni migration reported in experiments can be qualitatively reproduced through the competition between Ni(OH) 2 diffusion and Ni-YSZ wettability change. It is also found that the initial microstructure of the active layer and the microstructure characteristics in the support layer can strongly affect the distribution of steam partial pressure and overpotential, hence the Ni redistribution, in the active layer, which may explain the inconsistencies in experiments. Further, our results show that there are two missing pieces in current theory of Ni redistribution: a mechanism of fast Ni diffusion under humid atmosphere and the physical reason behind the Ni spreading/detachment.

25 ENERGY STORAGE↗

Link statistics of dislocation network during strain hardening

Dislocations are line defects in crystals that multiply and self-organize into a complex network during strain hardening. The length of dislocation links, connecting neighboring nodes within this network, contains crucial information about the evolving dislocation microstructure. By analyzing data from Discrete Dislocation Dynamics (DDD) simulations in face-centered cubic (fcc) Cu, we characterize the statistical distribution of link lengths of dislocation networks during strain hardening on individual slip systems. Here, our analysis reveals that link lengths on active slip systems follow a double-exponential distribution, while those on inactive slip systems conform to a single-exponential distribution. The distinctive long tail observed in the double-exponential distribution is attributed to the stress-induced bowing out of long links on active slip systems, a feature that disappears upon removal of the applied stress. We further demonstrate that both observed link length distributions can be explained by extending a one-dimensional Poisson process to include different growth functions. Specifically, the double-exponential distribution emerges when the growth rate for links exceeding a critical length becomes super-linear, which aligns with the physical phenomenon of long links bowing out under stress. This work advances our understanding of dislocation microstructure evolution during strain hardening and elucidates the underlying physical mechanisms governing its formation.

Crystal plasticity↗

Canister Centerline Cooling Experiments for DPF5-336 Reference Material Made with ERV3b Salt Simulant

This report provides experimental details and results on the elemental distributions (compositions), crystalline/amorphous phase distribution, and microstructures for DPF5-336 reference materials made with ERV3b salt that were subjected to the canister centerline cooling (CCC) heat treatment process. To fabricate these materials, ammonium dihydrogen phosphate was mixed with ERV3b salt simulant and a dechlorination procedure was run in the generation-2 dechlorinator at the Pacific Northwest National Laboratory up to 600°C in an alumina crucible. After the dechlorination process, the product was removed from the crucible, lightly crushed, loaded into a new alumina crucible, Fe 2 O 3 was added, and this was vitrified at 1100°C for 1 hour and quenched. The quenched material was mostly amorphous but showed some small diffraction peaks attributed to Li 3 Fe 2 (PO 4 ) 3 . This quenched material was ground to a fine particle size in a milling chamber, and this was subjected to the CCC cooling profile in duplicate within alumina and fused quartz crucibles. The CCC experiment is designed to simulate the slowest cooling profile (the vertical centerline) of a 0.61-m diameter canister where the slower cooling profile often results in the production of a variety of different crystalline phases upon formation of which result in a residual glass of unknown composition and properties. The information gleaned from these types of experiments will lead to improved composition-property predictions for materials with similar compositions. The results of the alumina and fused quartz experiments show different microstructures but similar diffraction patterns providing evidence of similar crystalline phases being present in each sample. Both samples had monazite [i.e., (Ce,Nd)PO 4 ] and mixtures of phases containing P, Fe, and alkali metal elements. This report completes the milestone M4FT-23PN030104042 with details provided in Appendix B.

36 MATERIALS SCIENCE↗

Derivation of heterogeneous material distributions and their sensitivity to HM-coupled two-phase flow models exemplified with the LASGIT experiment

Abstract Advective gas transport in bentonite, a possible buffer material in repositories for radioactive materials, is difficult to simulate in numerical continuum models, partly due to the complicated microstructure of bentonite. To generate reliable models of repositories nevertheless, spatially distributed heterogeneous material properties can be used to allow localization of gas flow. In this study, a pore-size-dependent stochastic approach of the gas entry pressure is derived from Mercury Intrusion Porosimetry, which is used to replicate measurements from the LASGIT experiment. In addition, three benchmark tests are simulated to investigate the dependence of heterogeneous distributions of material properties on the mesh discretization, the temporal dependence, and the coupling between the processes influenced by the heterogeneous parameters. The numerical modeling results of the LASGIT experiment show that the onset of gas flow into the system and the subsequent increase in pressure and stress can be well reproduced using heterogeneous distributions. Compared to a model with homogeneous material properties, heterogeneous distributions may allow the generation of dilatancy-controlled microfractures—an important feature with regard to the advective gas flow in bentonites. However, it can be observed that the heterogeneous distributions in LASGIT are less significant, as technical gaps or differences in material types could have a greater impact.

Environmental Sciences & Ecology↗

Effect of heat treatment on microstructural evolution and hardness homogeneity in laser powder bed fusion of alloy 718

The thermal history developed in laser powder bed fusion (LPBF) processes has been shown to be complex resulting in equally complex microstructures and mechanical properties. Microscopic observations and Vickers micro-hardness mapping measurements were carried out on diff ;erent section planes of LPBF alloy 718 cuboids. Three-dimensional finite element analysis was used to simulate thermal history and to predict the residual stress distribution in the as-built material. Computational thermodynamics was used to predict the micro-segregation and nucleation driving force of various phases in the bulk and in segregated regions. Varied heat-treatments such as simulated hot isostatic pressing, and double aging were applied. Their influence on the microstructure, microsegregation, precipitate formation, and micro-hardness variations of LPBF alloy 718 were investigated. Hardness map results showed heterogeneous micro-hardness on the xy- and xz-planes of the as-built parts where the bottom plane and center regions had larger hardness of ~315 HV 0.5 while the top plane and contours showed hardness of ~300 HV 0.5 . It was found that the aging treatment increased the overall hardness of the as-built condition from ~310 HV 0.5 to 470 HV 0.5 but also increased the hardness gradient throughout the coupon. After simulated hot isostatic pressing process (i.e., without applied pressure) at 1020 °C for 4 h followed by water quench (HIPWQ), the hardness gradient and hardness was minimized (~210 HV 0.5 ) as the microstructure transitioned from heterogeneous columnar grains in the as-built condition to more uniform recrystallized grains. A double aging treatment was applied to enhance hardness from ~210 HV 0.5 to ~440 HV 0.5 . Finally, HIPWQ followed by double aging produced a homogeneous microstructure and more uniform hardness map with enhanced mechanical properties in LPBF alloy 718 coupons.

36 MATERIALS SCIENCE↗

Design of a separate effects MiniFuel irradiation experiment investigating microstructure evolution in high burnup UO 2

The microstructural evolution of UO 2 fuel pellets during commercial operation in light water reactors (LWRs) is known to vary significantly across the pellet radius due to spatial variations in local temperature and burnup. The primary obstacle to extending LWR refueling cycles to 24-month intervals is the susceptibility of certain high burnup fuel microstructures to fuel fragmentation, relocation, and dispersal (FFRD) during a loss of coolant accident (LOCA). Although FFRD of the high burnup structure in the rim region of a pellet is well studied, the fine fragmentation that has been observed in a second region, near the midradius of the pellet (termed the “dark zone”) following mock LOCA testing of high burnup commercial fuel rods is less understood. This paper describes the design, analysis, and execution of a separate effects MiniFuel irradiation experiment that aims to identify the specific temperature and burnup regimes under which FFRD-susceptible dark zone microstructures form. The small disc specimens (3 mm diameter by ∼0.3 mm thick) enable more precise control of the relatively uniform temperature and burnup conditions. A total of 42 specimens were fabricated with typical LWR fuel densities (∼96%–98% of theoretical density) and grain sizes (∼12 μm) and are being irradiated over a range of temperatures (600°C–1000°C) and discharge burnups (50–72 MWd/kg-U) that bound the midradius region of high burnup LWR fuel. Fuel specimens with identical 235 U enrichments were inserted in two irradiation locations in the High Flux Isotope Reactor and are currently undergoing irradiation to further evaluate the impact of rate effects (fission rate, time at temperature) on the microstructural evolution. The fuel fabrication and the thermal and neutronic simulations used for designing the experiment are detailed in this paper. A secondary objective of the experiment is to observe fission gas release (FGR) under the various irradiation conditions, and this work provides first-order predictions of FGR from all fuel specimens. The insights gained from these experiments will inform future high burnup core designs that could minimize the formation of susceptible microstructures and ultimately enable 24-month refueling cycles while minimizing the fraction of the fuel susceptible to FFRD.

FFRD↗

Physics-Informed Machine Learning Model for Ceramic Matrix Composite Creep

A physics-informed recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear, time-dependent constitutive behavior of ceramic matrix composites (CMCs) driven by matrix damage and constituent creep at the microscale. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the high-fidelity generalized method of cells (HFGMC) approach which calls appropriate creep and damage models for each of the constituents. This coupling permits simulating the nonlinear behavior of CMCs based on constituent response at the microscale along with microstructural features such as fiber and porosity volume fraction and fiber radius. The microscale repeating unit cell is loaded under creep fatigue conditions to replicate the material loading experienced in a turbine engine. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input stress sequence, temperature, and microstructural features, the resulting strain history response while satisfying physical constraints related to creep rate, isochoric inelastic deformation, and strain energy density. The trained surrogate model is shown to effectively match the strain history over quantified distributions of microstructural features and relevant loading regimes and temperatures. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore, the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex microstructures, nonlinear time-dependent material response, and under non-monotonic loading conditions.

ceramic matrix composites↗

Microstructure Generation via Generative Adversarial Network for Heterogeneous, Topologically Complex 3D Materials

Using a large-scale, experimentally captured 3D microstructure dataset, we implement the generative adversarial network (GAN) framework to learn and generate 3D microstructures of solid oxide fuel cell electrodes. The generated microstructures are visually, statistically, and topologically realistic, with distributions of microstructural parameters, including volume fraction, particle size, surface area, tortuosity, and triple phase boundary density, being highly similar to those of the original microstructure. These results are compared and contrasted with those from an established, grain-based generation algorithm (DREAM.3D). Importantly, simulations of electrochemical performance, using a locally resolved finite element model, demonstrate that the GAN generated microstructures closely match the performance distribution of the original, while DREAM.3D leads to significant differences. Finally, the ability of the generative machine learning model to recreate microstructures with high fidelity suggests that the essence of complex microstructures may be captured and represented in a compact and manipulatable form.

36 MATERIALS SCIENCE↗

On residual stress analysis and microstructural evolution for stainless steel type 304 spent nuclear fuel canisters weld joint: Numerical and experimental studies

Chloride-induced stress corrosion cracking (CISCC) in the weldments of austenitic stainless steel canisters is one of the primary safety concerns during the dry storage of used nuclear fuel at Independent Spent Fuel Storage Installations in coastal areas. In order to evaluate the CISCC behavior in the canister, a 3D sequentially coupled thermo-mechanical finite element model was built to simulate the residual stresses induced by two intersecting longitudinal and circumferential multi-pass welds in austenitic stainless steel mockup canister. Weld-induced residual stresses from simulation were in good agreement with the experimental measurement results by deep-hole drilling and contour methods. Through-wall high tensile axial and hoop stresses were observed in both the fusion zone and heat affected zone (HAZ) in the longitudinal and circumferential weld, respectively. On the other hand, the microstructure within the weldments of 304 stainless steel mockup was characterized using optical microscopy, scanning electron microscopy, electron back-scattered diffraction and energy dispersive X-ray spectroscopy methods. A few amount of MnS inclusions in both the base metal and HAZ, and a lot of (Mn,Si) x O y inclusions in the fusion zone were found. Meanwhile, there is some α'-martensite in the base metal, while the martensite was reversed in the HAZ. In the fusion zone, there is a lot of δ-ferrite. The martensite reversion phenomenon was also verified by the combination of Feritscope test and thermal simulation, which could predict a martensite reversion temperature as 400C°. The welding residual stresses and microstructure analysis results could offer the instructive information for the prediction of CISCC behavior of the spent nuclear fuel canisters.

36 MATERIALS SCIENCE↗

Understanding the Space Weathering of Mercury Through Laboratory Experiments

Introduction: Airless surfaces across the solar system are continually modified by energetic particles from solar wind and micrometeoroid bombardment [1,2]. This process is known as space weathering, and it alters the chemical, microstructural, and optical properties of surface regoliths on airless bodies, including Mercury. On the Moon and S-type asteroids, the reflectance spectral signatures of space weathering include reddening (increasing reflectance with increasing wavelength), darkening (lowering of reflectance), and the attenuation of characteristic absorption bands [2]. Such spectral changes are driven by the production of Fe-bearing nanoparticles(npFe) through both solar wind irradiation and micrometeoroid bombardment. While our understanding of space weathering for the Moon and near-Earth S-types asteroids is advanced, insight into how these processes operate on other planetary bodies is limited. In particular, Mercury experiences a uniquely intense space weathering environment than planetary counterparts at 1 AU, including a moreintense solar wind flux and higher velocity micrometeoroid impacts [4]. Additionally, Mercury has a surface composition unique in the inner solar system, including regions of the surface with very low albedo known as the low reflectance material (LRM), which is enriched in carbon, likely graphite, up to 4wt.% [5].In addition, the concentration of Fe across Mercury’s surface islow (<2 wt.%) compared to the Moonor S-type asteroids asteroids[6]. Our understanding of the effects of space weathering on C-rich and Fe-poor phases is limited. Since Fe plays a critical role inthe development of space weathering characteristicson other airless surfaces(e.g., npFe), its limited availability may significantly affect the development of space weathering features in Mercury surface materials. We can simulate space weathering processes in the laboratory to explore their effects on the microstructural, chemical, and spectral characteristics of Mercury surface materials[7]. Here we used pulsed laser irradiation to simulate the short duration, high-temperature events associated with micrometeoroid impacts. We performed coordinated analyses including reflectance spectroscopy and electron microscopy to investigate the spectral, chemical, and microstructural changes in these mercurian analog samples. Methods: For these experiments, we usedforsteritic olivine with varying FeOcontents, a mineral phase proposed to be abundant on the surface of Mercury. We mixed each sample with graphite to simulate LRM regions of the surface. Wesynthesized the olivinesamplesat 1-bar at NASA’s Johnson Space Centerand prepared pressed powder pellets for laser irradiation [8].We prepared three samples, each with a base layer of olivineto maintain structural integrity and topped witha surface layer containing the graphite-olivine mixture: 1) Sample SC-001 San Carlos olivine(Fo90.91),2)Sample F-S-002 with0.05 wt.% FeO olivine, and 3) F-T-004 with 0.53 wt.% FeO olivine. Each sample was mixed with 5 wt.% powdered graphite and had grain sizes ranging from 45to 125μm. We irradiated each sample usinga pulsed Nd-YAG laser, (l=1064 nm, ~6 ns pulse duration, energy of 48 mJ/pulse) while undervacuumat Northern Arizona University. The laser was rastered1x and then 5x over the surfaceof each sample to simulate progressive space weathering. We collected in situreflectance spectra from the samples after each laser pulse witha Nicolet IS50 Fourier-Transform Infrared spectrometer (lfrom 0.65-2.5 μm). We used an FEI Nova NanoSEM200scanning electron microscope (SEM) and a Hitachi TM4000 Plus benchtop SEM at Purdue University to image the surface morphology and topography of the samples. We extracted thin sections for analysis in the transmission electron microscope(TEM)using the FEI Helios NanoLab 660 focused ion beam (FIB) SEM at the University of Arizona.We performed analysis of the microstructural and chemical characteristics of the samples using the 200 keV JEOL 2500 scanning TEM at Johnson Space Center. Reflectance Spectroscopy Results:Reflectance spectra for each sample are shown in Fig. 1.SC-001:The spectrum of the unirradiated sample exhibits a weak 1.0 μm absorption feature, associated with Fe2+in the olivine,and low overall reflectance (Fig. 1a). Thereflectance and the depth of the absorption band increases after 1x laser raster but are at their lowest after 5x laser rasters.F-T-004:The unirradiated sample has a blue-sloped spectrum with low reflectance without identifiable absorption features (Fig. 1b). With progressive laser irradiation, the sample reflectance increases and becomes strongly red-sloped.F-S-002:The unirradiated sample exhibits a dark, blue-slopedspectrum. The brightness of the sample increases significantly from <0.2average reflectance over >0.8 reflectance in the most irradiated sampleand thespectral slope also becomesslightly reddened (Fig. 1c). Microstructural and Chemical Analysis: Two primary alteration textures were observed in the samples exposed to simulated space weathering: 1) fluffy C-rich,and 2) vesiculated melt. The fluffy C-rich texture is composed oflow-densitydeposits distributed across the surface of the sample(Fig. 2A). Analysis of a FIB section extracted from a low-density C-rich region in sample SC-001 reveals multiple globule-type deposits, discrete from stacked graphite, likely produced via melting from the laser irradiation [9].The vesiculated melt textureis smooth and uniformly distributed across isolated regions of the sample surface. The vesicles measure up to 100s of nmin diameter. Analysis of a FIB section from this texture was extracted from sample F-T-004 reveals a layer of amorphous melt material, close to 100 nm thick and uniform across the FIB section (Fig. 2B). Isolated regions of this melt layer contain small nanoparticles, <5 nm in diameter. Chemical analysis through energy dispersive X-ray spectroscopy reveals the composition of this layer is enriched in Si and depleted in Mg and O compared to the underlying sample. Implications for Space Weathering on Mercury: Previous experiments simulating space weathering of Mercury have showndarkening and reddening of spectra[7,10].However, our use of low-Fe materials and graphite to create a sample set more analogous to the mercurian surface. Our results indicate that sample composition plays a significant and important role in the space weathering of Mercury. In particular, our spectral data demonstrates a strong correlation between spectral slope, Fe content, and simulated space weathering. While the variation in FeO content between samples F-S-002 and F-T-004 is <0.6 wt.%, the spectra deviate from flat to strongly red-sloped(F-T-004). This reddeningmay be linked to the presence of very small nanoparticles observed in the melt textures extracted from sample F-T-004. For the SC-001 sample, the fluffy C-rich textures may be developed by the amalgamation of small graphite particles into these unique morphologies. Such observations indicate that space weathering on Mercury may result in both familiar and new microstructural and chemical characteristics. References: [1]Hapke B. (2001) J. Geophys. Res.-Planet.,106,10039–10073. [2]Pieters C.M. and Noble S.K. (2016) J. Geophys. Res-Planet., 121, 1865–1884. [3] Lucey P.G., and Riner, M.A. (2011) Icarus,212, 451-462.[4]CintalaM.J.(1992)J. Geophys. Res.-Planet.,97,947–973.[5]Klima R.L.et al.(2018)Geophys.Res.Letters, 45, 2945–2953. [6]Nittler L.R., et al. (2011) Science 333, 1847-1850.[7]Sasaki S. and Kurahashi E. (2004) Space weathering on Mercury, Adv.Space Res., 33, 2152-2155.[8] Vander KaadenK.E., et al. (2018) LPSCXLIX, Abstract 1230. [9] McGlaun M.L. et al. (2019) LPSCL, Abstract 2019. [10] TrangD.et al. (2018)LPSCXLIX,Abstract2083

M S Thompson↗

Microstructural characterization and assessment of mechanical properties of concrete based on combined elemental analysis techniques and Fast-Fourier transform-based simulations

This paper presents a pixel-based modeling approach of concrete which combines an experimental characterization of concrete and the Fast-Fourier transform simulations. High-resolution phase maps created from experimental characterization by micro X-ray fluorescence, energy dispersive X-ray analysis, and X-ray diffraction contain 9 different phases, including 22.17 vol.% of hydrated cement paste, 54.21 vol.% of minerals, and 23.62 vol.% of interfaces. These phases provide the input for determining the effective elastic properties, coupled with Fast-Fourier transform-based simulation. The simulation results show that the effective range of Young’s modulus of concrete is comparable with the range of experimental values ~ 37 ± 4 GPa with the assumption of realistic properties of interfaces.

36 MATERIALS SCIENCE↗

A coupled fluid-mechanical workflow to simulate the directed energy deposition additive manufacturing process

Simulation of additive manufacturing processes can provide essential insight into material behavior, residual stress, and ultimately, the performance of additively manufactured parts. In this work, we describe a new simulation based workflow utilizing both solid mechanics and fluid mechanics based formulations within the finite element software package SIERRA (Sierra Solid Mechanics Team in Sierra/Solid Mechanics 4.52 User’s Guide SAND2019-2715. Technical report, Sandia National Laboratories, 2011) to enable integrated simulations of directed energy deposition (DED) additive manufacturing processes. In this methodology, a high-fidelity fluid mechanics based model of additive manufacturing is employed as the first step in a simulation workflow. This fluid model uses a level set field to track the location of the boundary between the solid material and background gas and precisely predicts temperatures and material deposition shapes from additive manufacturing process parameters. Next, the resulting deposition shape and temperature field from the fluid model are then mapped into a solid mechanics formulation to provide a more accurate surface topology for radiation and convection boundary conditions and a prescribed temperature field. Solid mechanics simulations are then conducted to predict the evolution of material stresses and microstructure within a part. By combining thermal history and deposition shape from fluid mechanics with residual stress and material property evolutions from solid mechanics, additional fidelity and precision are incorporated into additive manufacturing process simulations providing new insight into complex DED builds.

36 MATERIALS SCIENCE↗

Predicting Melt Properties Using Atomistic Simulations With A Highly Accurate Physically Informed Neural Network Interatomic Potential

The use of a recently developed machine learning (ML) interatomic potential for molecular dynamics simulations of aluminum melt properties will be presented. Such properties are critical for process modeling in additive manufacturing, including the melt pool size, solidification, and formation of solidification microstructures. Direct first-principles modeling of these processes is computationally prohibitive whereas simulations employing ML potentials combine the high accuracy of quantum-mechanical methods with high computational speeds. The physically-informed neural network (PINN) method used herein, integrates a high-dimensional regression implemented by an artificial neural network with a physics-based bond-order interatomic potential. PINN potentials can accurately reproduce many properties of aluminum in both crystalline-solid and liquid phases. We examine the accuracy of a PINN Al potential in predicting the density, self-diffusivity, viscosity, and the tension of the liquid surface and liquid-solid interfaces. Comparison with experimental data and ab initio molecular dynamics calculations shows very good agreement for all properties tested.

molecular dynamics↗

Melt-processing of lunar ceramics

The goal of this project is to produce useful ceramics materials from lunar resources using the by products of lunar oxygen production processes. Emphasis is being placed on both fabrication of a variety of melt-processed ceramics, and on understanding the mechanical properties of these materials. Previously, glass-ceramics were formed by casting large glass monoliths and heating these to grow small crystallites. The strengths of the resulting glass-ceramics were found to vary with the inverse square root of the crystal grain size. The highest strengths (greater than 300 MPa) were obtained with the smallest crystal sizes (less than 10 microns). During the past year, the kinetics of crystallization in simulated lunar regolith were examined in an effort to optimize the microstructure and, hence, mechanical properties of glass ceramics. The use of solar energy for melt-processing of regolith was examined, and strong (greater than 630 MPa) glass fibers were successfully produced by melt-spinning in a solar furnace. A study of the mechanical properties of simulated lunar glasses was completed during the past year. As on Earth, the presence of moisture was found to weaken simulated lunar glasses, although the effects of surface flaws was shown to outweigh the effect of atmospheric moisture on the strength of lunar glasses. The effect of atmospheric moisture on the toughness was also studied. As expected, toughness was found to increase only marginally in an anhydrous atmosphere. Finally, our efforts to involve undergraduates in the research lab fluorished this past year. Four undergraduates worked on various aspects of these projects; and two of them were co-authors on papers which we published.

Fabes, B. D.↗

A molecular dynamics survey study of impact-induced amorphization in yttria nanoparticles

The impact of a single yttria (Y 2 O 3 ) nanoparticle onto a yttria substrate was studied as a function of particle velocity (300–1200 m/s) and diameter (5–12 nm) for three particle orientations and two substrate orientations using molecular dynamics simulations. Impacts resulted in minimal dislocation motion but extensive disordering of the crystal lattice without approaching the melting temperature. The resulting amorphous microstructure was both precisely coincident with plastic strain and retained until the end of the simulation, suggesting a deformation mechanism of amorphization-facilitated viscous flow. Further analysis revealed that the degree of amorphization depended primarily on the impact velocity and secondarily on the orientations of the particle and substrate, with significant second-order effects from variable interactions.

74 ATOMIC AND MOLECULAR PHYSICS↗

Three-dimensional microstructural characterization of FBR MOX fuel and the contribution of microstructural features to the thermal conductivity of the fuel

Combination of microstructural characterization, property measurements, and phase field modeling is used to investigate fast breeder reactor (FBR) mixed oxide (MOX) fuel irradiated to burnup of 13.7% fissions per initial metal atom (FIMA). Here, the fuel was characterized at different radial locations, which revealed that grey phase can be present in the central region if it nucleates on five metal precipitates (FMPs). In addition, in the mid-radial region FMPs do not diffuse out of the region once formed and the size of Pd–Te precipitates is dictated by the porosity present in the region. Thermal conductivity measurements were conducted as a function of radial location and the microstructure of the fuel was correlated with the observed trend. Reconstructions of the 3D solid and gaseous fission product structures in different regions of the fuel were used to simulate the effective thermal conductivity (ETC) of the respective regions and determine which microstructural feature has the strongest impact on thermal conductivity. Based on conducted assessment, FMPs and Pd-Te precipitates improve local conductivity of central and mid-radial regions even in the presence of grey phase, but defects are primary contributor to the degradation of thermal conductivity on the periphery of the fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗