Search NASA⌕ Search

SEARCH · Search NASA

Results for “Defect kinetics”

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 127 records · Page 7

Effect of void positioning on the detonation sensitivity of a heterogeneous energetic material

We show although it is well-established that voids profoundly influence the initiation and reaction behaviors of heterogeneous energetic materials such as polymer-bonded explosives (PBX) and propellants, there has been little study of how void location in different constituents in the microstructures of such materials affect the macroscale behavior. Here, we use three-dimensional (3D) mesoscale simulations to study how void placement within the reactive grains versus the polymer binder influences the shock-to-detonation transition (SDT) in a polymer-bonded explosive. The material studied here has a microstructure comprised of 75% PETN (pentaerythritol tetranitrate) grains and 25% HTPB (hydroxyl-terminated polybutadiene) polymer binder by volume. Porosities up to 10% in the form of spherical voids distributed in both the grains and polymer are considered. An Arrhenius reactive burn relation is used to model the chemical kinetics of the PETN grains under shock loading, thereby resolving the heterogeneous detonation behavior of the PBX. The influence of void location on the shock initiation sensitivity of the material is quantitatively ranked by comparing the predicted run distance to detonation (RDD) for each sample. The analysis includes inherent quantification of uncertainties arising from the stochastic variations in the microstructure morphologies and void distributions by using statistically equivalent microstructure sample sets (SEMSS), leading to probabilistic formulations for the RDD as a function of shock pressure. The calculations reveal that the location of voids in the composite microstructure significantly affects the RDD. Specifically, voids exclusively within the grains cause the PBX to be more sensitive (having shorter RDD) than voids in the polymer binder. Unique probabilistic relationships are derived to map the probability of observing RDD for each void location material case, allowing for prediction of initiation behavior anywhere in the shock pressure – RDD space. These findings agree with trends reported in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evidence of non-isentropic release from high residual temperatures in shocked metals measured with ultrafast x-ray diffraction

Shock experiments are widely used to understand the mechanical and electronic properties of matter under extreme conditions. However, after shock loading to a Hugoniot state, a clear description of the post-shock thermal state and its impacts on materials is still lacking. We used diffraction patterns from 100-fs x-ray pulses to investigate the temperature evolution of laser-shocked Al–Zr metal film composites at time delays ranging from 5 to 75 ns driven by a 120-ps short-pulse laser. We found significant heating of both Al and Zr after shock release, which can be attributed to heat generated by inelastic deformation. A conventional hydrodynamic model that employs (i) typical descriptions of Al and Zr mechanical strength and (ii) elevated strength responses (which might be attributed to an unknown strain rate dependence) did not fully account for the measured temperature increase, which suggests that other strength-related mechanisms (such as fine-scale void growth) could play an important role in thermal responses under shock wave loading/unloading cycles. Our results suggest that a significant portion of the total shock energy delivered by lasers becomes heat due to defect-facilitated plastic work, leaving less converted to kinetic energy. This heating effect may be common in laser-shocked experiments but has not been well acknowledged. High post-shock temperatures may induce phase transformation of materials during shock release. Another implication for the study is the preservability of magnetic records from planetary surfaces that have a shock history from frequent impact events.

42 ENGINEERING↗

On the prediction of equilibrium states in homogeneous turbulence

A comparison of several commonly used turbulence models (including the Kappa-epsilon and two second-order closures) is made for the test problem of homogeneous turbulent shear flow in a rotating frame. The time evolution of the turbulent kinetic energy and dissipation rate is calculated for a variety of models and comparisons are made with previously published experiments and numerical simulations. Particular emphasis is placed on examining the ability of each model to accurately predict equilibruim states for a range of the parameter Omega/S (the ratio of the rotation rate to the shear rate). It is found that none of the commonly used second-order closure models yield substantially improved predictions for the time evolution of the turbulent kinetic energy and dissipation rate over the somewhat defective results obtained from the simpler Kappa-epsilon model for the turbulent flow regime. There is also a problem with the equilibrium states predicted by the various models. For example, the Kappa-epsilon model erroneously yields equilibrium states that are independent of Omega/S while the Launder, Reece, and Rodi model predicts a flow relaminarization when Omega/S is greater than 0.39 - a result which is contrary to numerical simulations and linear spectral analysis which indicate flow instability for at least the range 0 less than or = Omega/S less than or = 0.5. The physical implications of the results obtained from the various turbulence models considered here are discussed in detail along with proposals to remedy the deficiencies based on a dynamical systems approach.

Speziale, Charles G.↗

On the prediction of equilibrium states in homogeneous turbulence

A comparison of several commonly used turbulence models (including the Kappa-epsilon and two second-order closures) is made for the test problem of homogeneous turbulent shear flow in a rotating frame. The time evolution of the turbulent kinetic energy and dissipation rate is calculated for a variety of models and comparisons are made with previously published experiments and numerical simulations. Particular emphasis is placed on examining the ability of each model to accurately predict equilibrium states for a range of the parameter Omega/S (the ratio of the rotation rate to the shear rate). It is found that none of the commonly used second-order closure models yield substantially improved predictions for the time evolution of the turbulent kinetic energy and dissipation rate over the somewhat defective results obtained from the simpler Kappa-epsilon model for the turbulent flow regime. There is also a problem with the equilibrium states predicted by the various models. For example, the Kappa-epsilon model erroneously yields equilibrium states that are independent of Omega/S while the Launder, Reece, and Rodi model predicts a flow relaminarization when Omega/S is greater than 0.39 - a result which is contrary to numerical simulations and linear spectral analysis which indicate flow instability for at least the range 0 less than or = Omega/S less than or = 0.5. The physical implications of the results obtained from the various turbulence models considered here are discussed in detail along with proposals to remedy the deficiencies based on a dynamical systems approach.

Speziale, C. G.↗

The Structure of Ice Nanoclusters and Thin-films of Water Ice: Implications for Icy Grains in Cold Molecular Clouds

The cubic to hexagonal phase transformation in water ice (I(sub c) yields I(sub h)) is used to measure the extent to which surface structure and impurities control bulk properties. In pure crystalline (I(sub c)) water ice nanoclusters and in thin-films of impure water ice, I(sub c) yields I(sub h) occurs at lower temperatures than in thin-films of pure water ice. The disordered surface of the 20 nm diameter nanoclusters promotes transformations or reactions which would otherwise be kinetically hindered. Likewise, impurities such as methanol introduce defects into the ice network, thereby allowing sluggish structural transitions to proceed. Such surface-related phenomena play an important role in promoting chemical reactions on interstellar ice grains within cold molecular clouds, where the first organic compounds are formed.

Delzeit, Lance↗

Assessment of Shape Memory Alloys - From Atoms To Actuators - Via In Situ Neutron Diffraction

As shape memory alloys (SMAs) become an established actuator technology, it is important to identify the fundamental mechanisms responsible for their performance by understanding microstructure performance relationships from processing to final form. Yet, microstructural examination of SMAs at stress and temperature is often a challenge since structural changes occur with stress and temperature and microstructures cannot be preserved through quenching or after stress removal, as would be the case for conventional materials. One solution to this dilemma is in situ neutron diffraction, which has been applied to the investigation of SMAs and has offered a unique approach to reveal the fundamental micromechanics and microstructural aspects of bulk SMAs in a non-destructive setting. Through this technique, it is possible to directly correlate the micromechanical responses (e.g., internal residual stresses, lattice strains), microstructural evolutions (e.g., texture, defects) and phase transformation properties (e.g., phase fractions, kinetics) to the macroscopic actuator behavior. In this work, in situ neutron diffraction was systematically employed to evaluate the deformation and transformation behavior of SMAs under typical actuator conditions. Austenite and martensite phases, yield behavior, variant selection and transformation temperatures were characterized for a polycrystalline NiTi (49.9 at. Ni). As the alloy transforms under thermomechanical loading, the measured textures and lattice plane-level variations were directly related to the cyclic actuation-strain characteristics and the dimensional instability (strain ratcheting) commonly observed in this alloy. The effect of training on the shape memory characteristics of the alloy and the development of two-way shape memory effect (TWSME) were also assessed. The final conversion from a material to a useful actuator, typically termed shape setting, was also investigated in situ during constrained heatingcooling and subsequent shape recovery experiments. Neutron diffraction techniques are also being applied to the investigation of novel high temperature SMAs with the objective of designing alloys with better stability, higher transition temperatures and ultimately superior durability.

Shape Memory Alloys↗

Defect‐Induced, Ferroelectric‐Like Switching and Adjustable Dielectric Tunability in Antiferroelectrics

Abstract Antiferroelectrics, which undergo a field‐induced phase transition to ferroelectric order that manifests as double‐hysteresis polarization switching, exhibit great potential for dielectric, electromechanical, and electrothermal applications. Compared to their ferroelectric cousins, however, considerably fewer efforts have been made to understand and control antiferroelectrics. Here, it is demonstrated that the polarization switching behavior of an antiferroelectric can be strongly influenced and effectively regulated by point defects. In films of the canonical antiferroelectric PbZrO 3 , decreasing oxygen pressure during deposition (and thus increasing adatom kinetic energy) causes unexpected “ferroelectric‐like” polarization switching although the films remain in the expected antiferroelectric orthorhombic phase. This “ferroelectric‐like” switching is correlated with the creation of bombardment‐induced point‐defect complexes which pin the antiferroelectric–ferroelectric phase boundaries, and thus effectively delay the phase transition under changing field. The effective pinning energy is extracted via temperature‐dependent switching‐kinetics studies. In turn, by controlling the concentration of defect complexes, the dielectric tunability of the PbZrO 3 can be adjusted, including being able to convert between “positive” and “negative” tunability near zero field. This work reveals the important role and strong capability of defects to engineer antiferroelectrics for new performance and functionalities.

Pan, Hao↗

Point Defects Control Guest Molecule Diffusion in the 1D Pores of Zn(tbip)

Molecular diffusion plays a critical role in metal-organic frameworks (MOFs) within the application of kinetic chemical separations. We carefully study in this work the unexpected role of point defects for short-chain alkanes diffusing in Zn(tbip), an MOF with rigid one-dimensional (1D) channels. Inside a defect-free Zn(tbip), guest molecules are expected to follow single-file diffusion along 1D channels. It has been found previously by Heinke et al. that these parallel 1D channels are connected for molecular diffusion. Our density functional theory (DFT) calculations suggest that linker vacancy defects could arise under experimentally relevant conditions by removing a pair of linkers. Further climbing-image nudged elastic band (cNEB) DFT calculations indicate that hopping of short-chain alkanes between adjacent 1D channels over defect windows can occur at moderate temperatures. In addition to providing insights into connected adjacent 1D channels in Zn(tbip), Heinke et al. also inferred that most 1D pores are blocked from a microkinetic model. Additionally, we explored the influence of hydrolyzed linker created by the formation of linker removal inside 1D pores. Our DFT calculations show that the linkers can effectively block the pores and the linker diffusion in 1D channels is slow. Our results, for the first time, offer a mechanistic explanation of the unexpected molecular diffusion behavior in this MOF with 1D channels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Complex Intergrowths of Non-Stoichiometric Defect-Structured Hibonite and Al-Rich Spinel in an Allende Ca-Al-Rich Inclusion

Hibonite is a primary ultra-refractory mineral occurring in many Ca-Al-rich inclusions (CAIs) and is predicted to condense as the second major phase from a cooling gas of solar composition. Our previous microstructural studies of hibonite in carbonaceous chondrites revealed its unique microstructures consisting of numerous defects that contain Mg-enriched, wider spinel blocks in stoichiometric hibonite. Han, J. et al. [2 items] demonstrated experimentally that defect-structured hibonites can grow easily in the presence of minor Mg and are kinetically more stable than equilibrium assemblages predicted by thermodynamic calculations. However, a thermo-dynamic vs. structural stability of defect-structured hibonite relative to other early-condensed Al-rich phases such as corundum, grossite, and spinel remains poorly constrained. Here, we present the results of atomic resolution TEM (Transmission Electron Microscopy) imaging of hibonite in a compact Type A CAI in the Allende CV3 chondrite in order to better understand the crystal structure and chemistry of defect-structured hibonite and its associated Al-rich phases, especially non-stoichiometric, Al-rich spinel, in the context of the formation of first refractory solids in the early solar nebula.

Han, J.↗

Correlations Among Microstructure, Morphology, Chemistry, and Isotopic Systematics of Hibonite in CM Chondrites

Introduction: Hibonite is a primary refractory phase occurring in many CAIs, typically with spinel and perovskite. Our microstructural studies of CAIs from carbonaceous chondrites reveal a range of stacking defect densities and correlated non-stoichiometry in hibonite. We also conducted a series of annealing experiments, demonstrating that the Mg-Al substitution stabilized the formation of defect-structured hibonite. Here, we continue a detailed TEM analysis of hibonite-bearing inclusions from CM chondrites that have been well-characterized isotopically. We examine possible correlations of microstructure, morphology, mineralogy, and chemical and isotopic systematics of CM hibonites in order to better understand the formation history of hibonite in the early solar nebula. Methods: Fifteen hibonite-bearing inclusions from the Paris CM chondrite were analyzed using a JEOL 7600F SEM and a JEOL 8530F electron microprobe. In addition to three hibonite-bearing inclusions from the Murchison CM chondrite previously reported, we selected three inclusions from Paris, Pmt1-6, 1-9, and 1-10, representing a range of 26Al/27Al ratios and minor element concentrations for a detailed TEM study. We extracted TEM sections from hibonite grains using a FEI Quanta 3D field emission gun SEM/FIB. The sections were then examined using a JEOL 2500SE field-emission scanning TEM equipped with a Thermo-Noran thin window EDX spectrometer. Results and Discussion: A total of six hibonite-bearing inclusions, including two platy hibonite crystals (PLACs) and four spinel-hibonite inclusions (SHIBs), were studied. There are notable differences in chemical and isotopic compositions between the inclusions (Table 1), indicative of their different formation environment or timing. Our TEM observations show perfectly-ordered, stoichiometric hibonite crystals without stacking defects in two PLACs, 2-7-1 and 2-8-2, and in three SHIBs, Pmt1-6, 1-9, and 1-10. In contrast, SHIB 1-9-5 hibonite grains contain a low density of stacking defects linked to an increase in MgO contents, indicating complex, disordered intergrowths of stoichiometric and MgO-enriched hibonites. From the data collected to date, we find no clear correlation between the microstructures of hibonite and its morphological and mineralogical types that reflect distinct chemical and isotopic systematics [6-8,10]. Interestingly, the presence of no or few stacking defects in hibonite from the PLACs and SHIBs are in contrast to our experimental studies that produced very high densities of stacking defects in hibonite [3-5]. Unlike our experi-ments, electron microprobe data from the PLACs and SHIBs hibonite grains show a strong correlation between (Ti4++Si4+) and Mg2+ cations, suggesting that coupled substitutions of (Ti4++Mg2+) and (Si4++Mg2+) for 2Al3+ inhibit the formation of defect-structured hibonite. However, our experimental studies suggest that kinetics (e.g., cooling rate) or other thermal effects also exert a strong control on the microstructures and chemical compositions of hibonite. In Pmt1-6, elongated perovskite grains present at the hibonite grain boundaries display (121) twinning, indicative of a fast cooling (>50degC/min) after high-temperature events. Therefore, the nebular microstructural characteristics of hibonite, at least in this inclusion, would not have destroyed by subsequent high-temperature annealing. Conclusions: Our TEM observations thus far show no clear correlation in microstructures, morphological and mineralogical characteristics, and chemical and isotopic systematics of hibonites from CM chondrites. The observed variation in stacking defect densities in the hibonites may be controlled by thermal processes in the early solar nebula. A detailed TEM analysis of additional CM hibonite samples is underway to evaluate this hypothesis.

Han, J.↗

Reaction–drift–diffusion models from master equations: application to material defects

We present a general method to produce well-conditioned continuum reaction–drift–diffusion equations directly from master equations on a discrete, periodic state space. We assume the underlying data to be kinetic Monte Carlo models (i.e. continuous-time Markov chains) produced from atomic sampling of point defects in locally periodic environments, such as perfect lattices, ordered surface structures or dislocation cores, possibly under the influence of a slowly varying external field. Our approach also applies to any discrete, periodic Markov chain. Here, the analysis identifies a previously omitted non-equilibrium drift term, present even in the absence of external forces, which can compete in magnitude with the reaction rates, thus being essential to correctly capture the kinetics. To remove fast modes which hinder time integration, we use a generalized Bloch relation to efficiently calculate the eigenspectrum of the master equation. A well conditioned continuum equation then emerges by searching for spectral gaps in the long wavelength limit, using an established kinetic clustering algorithm to define a proper reduced, Markovian state space.

36 MATERIALS SCIENCE↗

A cryogenic white light absorption spectroscopy setup with in situ gamma irradiation and thermo-optical annealing for optical fiber radiation-induced attenuation characterization

Silica optical fiber sensors offer a fast, distributed measurement solution in various cryogenic and radiation environments, such as magnets for fusion power. Under these conditions, light-absorbing point defects limit lifetime via radiation-induced attenuation (RIA). To support RIA kinetics prediction, we present an in situ, broadband absorption spectroscopy setup combining gamma irradiation with liquid nitrogen cooling. The apparatus enables continuous monitoring of narrowband RIA levels and the use of secondary optical annealing light sources, alongside broadband spectrum measurements to characterize RIA defects through spectrum decomposition. Initial results confirm the inevitable photobleaching effect of the probe light source, which must be accounted for. In addition, we report RIA kinetics during cycles of gamma irradiation at 77 K and isochronal thermal annealing steps from liquid nitrogen to room temperature, showcasing the setup’s ability to replicate real fiber operation scenarios. These instruments form an ideal platform to further study the kinetics of RIA coupled with thermal and optical annealing.

Absorption spectroscopy↗

Single-Crystallization of O3-Type Layered Oxide Cathode for Na-Ion Battery

The development of high-energy-density Na-ion batteries places significant demands on single-crystal layered oxide cathodes, especially for further high-voltage, solid-state battery scenarios. In the O3-type structure, due to the original sluggish Na ion diffusion kinetics (approximately 1 order of magnitude lower than that of Li-ion), and further hindrance against diffusion kinetics caused by single-crystal architecture, these inherent defects lead to the decline in the electrochemical performance. Herein, we demonstrated that the single crystallization of O3-type NaNi 1/3 Fe 1/3 Mn 1/3 O 2 cathode (d 50 = 5.04 μm) aggravates surface-to-bulk phase inhomogeneity distribution, which is attributed to the uneven Na ions extraction. Moreover, the Na-depletion of the surface/shell region not only aggravates Na ion diffusion resistance but also leads to a higher valence state of transition-metal elements (e.g., Ni/Fe) near the surface of the single-crystal particle, which further compromises the cathode-electrolyte interface stability. Furthermore, not limited to revealing the challenges, tuning the particle size and moderating quasi-single-crystal strategies have been proven to effectively mitigate the negative uneven distributions of Na ions, phases, and valence/oxidative states, resulting in efficient modification for single crystallization of Na-layered oxide cathodes.

36 MATERIALS SCIENCE↗

Directed self-assembly of block copolymers for high-precision patterning in the era of extreme ultraviolet lithography

Extreme ultraviolet (EUV) lithography enables unprecedented resolution in semiconductor patterning but faces critical challenges in developing resist materials that achieve high-precision at economically viable throughput. Directed self-assembly (DSA) of block copolymers (BCPs) offers a promising solution for pattern rectification by leveraging thermodynamically determined domain structures to decouple BCP pattern quality from the imperfect original lithographic pattern. This prospective presents an overview of recent progress on the EUV + DSA strategy, covering advances in BCP material design, processing, metrology, and pattern transfer. We highlight recent advances in high-χ BCPs with perpendicular orientation and domain spacings compatible with EUV dimensions, leveraging A-b-(B-r-C) architectures. We also discuss progress in chemical pre-pattern fabrication using both positive- and negative tone resists, along with processing strategies to minimize defects and roughness based on BCP thermodynamics and assembly kinetics. We further examine metrology platforms for characterizing the thermodynamics of BCP materials and quantifying the size and shape of BCP domains. Lastly, we review pattern transfer strategies for generating functional inorganic masks suitable for semiconductor manufacturing. Together, these advances highlight the potential of DSA to complement EUV lithography, offering a pathway to address critical challenges in achieving high-precision patterning for the semiconductor industry.

Lee, Kyunghyeon [Univ. of Chicago, IL (United Stat↗

Optimizing ensemble NV − spin properties of fluorescent diamond microparticles by systematic low pressure high temperature annealing

Low pressure high temperature annealing is a means for driving nitrogen and defect diffusion in diamond to reduce internal lattice damage without the need for technically complicated high-pressure cells. Herein, we perform a systematic time (5, 15, and 30 min) and temperature (1200 °C–1800 °C) study of effects of low-pressure high temperature annealing on photoluminescence, spin concentrations, and spin relaxation properties of NV centers in ca. 3 μm synthetic type 1b diamond particles. Annealing in the temperature range of ca. 1400 °C–1700 °C for even 5 min leads to a higher optically detected magnetic resonance contrast as compared to standard annealing at 900 °C for 2 h. Particles annealed at 1700 °C for 5 min exhibit a contrast close to about 13% as compared to about 9% for those annealed at 900 °C for 2 h. A reduction in the zero-field splitting strain parameter from E ≈ 4.5 MHz to ≈ 2.5 MHz and spectral linewidth from Δν ≈ 7 MHz to ≈ 4 MHz are observed even after 5 min annealing at 1700 °C. Improvements in these spectral parameters resulted in a roughly 2-fold reduction in the noise level of temperature monitoring experiment utilizing an ensemble of NV centers in the particles. Annealing in the temperature range of 1600 °C for 15 or 30 min or 1700 °C for 5 min resulted in NV T 1 relaxation times approaching ca. 5 ms typically observed for bulk diamond. Quantitative electron paramagnetic resonance (EPR) allowed for estimations of thermal activation energies of paramagnetic center annihilation. Monitoring the primary defect concentration (P1 and other defects with half integer spins) and utilizing second order kinetic modeling, an activation energy of 3.63 ± 0.28 eV was estimated. Alternatively, using the NV half field EPR signal and first order kinetic modeling, a similar activation energy 3.89 ± 0.29 eV was estimated.

NV centers↗

Thermophysical Property Measurements of Molten Semiconductors in 1-g and Reduced-g Condition

Understanding and controlling the formation kinetics of varieties of crystal imperfections such as point defects, non uniform distribution of doping atoms, and impurity atoms in growing crystals are very important. Theoretical (numerical) modeling of the crystal growth process is an essential step to achieving these objectives. In order to obtain reliable modeling results, input parameters, i.e. various thermophysical parameters, must be accurate. The importance of accurate thermophysical properties of semiconductors in crystal growth cannot be overly emphasized. The total hemispherical emissivity, for instance, has a dramatic impact on the thermal environment. It determines the radiative emission from the surface of the melt which determines to a large extent the profile of the solidified crystal. In order to understand the convection and the turbulence in a melt, viscosity becomes an important parameter. The liquid surface tension determines the shape of the liquid-atmosphere interface near the solid-liquid-atmosphere triple point. Currently used values for these parameters are rather inaccurate, and this program intends to provide more reliable measurements of these thermophysical properties. Thus, the objective of this program is in the accurate measurements of various thermophysical properties which can be reliably used in the modeling of various crystal growth processes. In this program, thermophysical properties of molten semiconductors, such as Si, Ge, Si-Ge, and InSb will be measured as a function of temperature using the High Temperature Electrostatic Levitator at JPL. Each material will be doped by different kinds of impurities at various doping levels. Thermophysical properties which will be measured include: density, thermal expansion coefficient, surface tension, viscosity, specific heat, hemispherical total emissivity, and perhaps electrical and thermal conductivities. Many molten semiconductors are chemically reactive with crucibles. As a result, these dispersed impurities in the melts tend to substantially modify the properties of pure semiconductors. Sample levitation done in a vacuum clearly helps maintain the sample purity. However, in the 1-g environment, all gravity caused effects such as convection, sedimentation and buoyancy are still present in the sample. In addition, large forces needed to levitate a sample in the presence of the gravity can cause additional flows in the melt. The use of the High Temperature Electrostatic Levitator (HTESL) for the present research is a recent development and little is known about the flows induced by the electrostatic forces. In this ground base program, we will define the limits of HTESL technology as various thermophysical properties of molten semiconductors are measured.

Rhim, Won-Kyu↗

Cluster Dynamics Modeling Needs for the Advanced Materials and Manufacturing Technologies Program

This milestone report aims to identify and assess the cluster dynamics (CD) modeling requirements within the Department of Energy's Office of Nuclear Energy (DOE-NE) Advanced Materials and Manufacturing Technologies (AMMT) program and to communicate these needs to the DOE-NE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The goal is to ensure NEAMS is well-informed about the CD modeling requirements to support AMMT's mission of accelerating the development, qualification, demonstration, and deployment of advanced structural materials and manufacturing for nuclear energy applications. CD modeling is an essential tool for predicting the degradation of structural materials under irradiation, which is a key component of AMMT's accelerated qualification process. The AMMT program focuses on both additively manufactured and wrought structural alloys, such as laser powder-bed fusion 316H austenitic stainless steel, alloy 709, Haynes 244, and alloy 617. These materials require a generalized CD modeling framework to facilitate rapid model development and computational simulation. A flexible, generalized CD software, similar to the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework, would enable modeling of various cluster types, including defect clusters, defect-solute clusters, and multicomponent clusters, incorporating thermodynamics and kinetics parameters. Radiation effects, microstructural feature evolution, and multi-dimensional modeling are critical considerations for the CD model. The usability of the CD code should allow for easy modification and coupling with MOOSE-based simulations. Additionally, the software should adhere to Nuclear Quality Assurance-1 standards, include a testing suite for verification and validation, and be version-controlled within a national laboratory-managed Git repository. Benchmark problems are needed to assess code predictions and performance.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Nonlinear Response of Layer Growth Dynamics in the Mixed Kinetics-Bulk-Transport Regime

In situ high-resolution interferometry on horizontal facets of the protein lysozyme reveal that the local growth rate R, vicinal slope p, and tangential (step) velocity v fluctuate by up to 80% of their average values. The time scale of these fluctuations, which occur under steady bulk transport conditions through the formation and decay of step bunches (macrosteps), is of the order of 10 min. The fluctuation amplitude of R increases with growth rate (supersaturation) and crystal size, while the amplitude of the v and p fluctuations changes relatively little. Based on a stability analysis for equidistant step trains in the mixed transport-interface-kinetics regime, we argue that the fluctuations originate from the coupling of bulk transport with nonlinear interface kinetics. Furthermore, step bunches moving across the interface in the direction of or opposite to the buoyancy-driven convective flow increase or decrease in height, respectively. This is in agreement with analytical treatments of the interaction of moving steps with solution flow. Major excursions in growth rate are associated with the formation of lattice defects (striations). We show that, in general, the system-dependent kinetic Peclet number, Pe(sub k) , i.e., the relative weight of bulk transport and interface kinetics in the control of the growth process, governs the step bunching dynamics. Since Pe(sub k) can be modified by either forced solution flow or suppression of buoyancy-driven convection under reduced gravity, this model provides a rationale for the choice of specific transport conditions to minimize the formation of compositional inhomogeneities under steady bulk nutrient crystallization conditions.

Vekilov, Peter G.↗