SEARCH · Engineering Papers

Results for “Statistical mechanics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Effect of Groundwater Age and Recharge Source on Nitrate Concentrations in Domestic Wells in the San Joaquin Valley

Nitrate is one of the most abundant contaminants in groundwater globally, in the United States, and in California (CA). We studied well construction information, water chemistry, stable isotopes, and noble gases to understand how groundwater travel time and recharge source and mechanism control nitrate concentrations in domestic wells in the San Joaquin Valley (SJV), CA, a large semiarid, irrigated agricultural region. Using nonparametric statistics, we find a decreasing trend in nitrates with groundwater travel time and well depth. Samples collected from wells that are closer to rivers and that show indications of river water recharge, either low recharge temperature or low δ 18 O signature, have lower concentrations of nitrates than samples with isotopic signatures indicating mixed source or local precipitation recharge. The curbing effect of river water recharge on nitrate concentrations in domestic wells is similar for direct river recharge and water applied as irrigation. This suggests that irrigation with river water also has a diluting effect that reduces the concentration of nitrate found in groundwater. This conclusion supports the idea that flood-managed aquifer recharge may be considered for remediation of groundwater nitrate when designing replenishment of aquifers.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Small-Magnitude Seismic Swarms in Central Utah (US): Interactions of Regional Tectonics, Local Structures and Hydrothermal Systems

Swarms in Central Utah are situated in the complex transition between the Basin and Range (BR) province and the Colorado Plateau. Transecting transverse structures, volcanic deposits, and hydrothermal systems complicate the extensional BR horst and graben structures and provide a multitude of plausible triggering mechanisms. Revisiting the catalog of the University of Utah Seismograph Stations (1981–2022), we analyze spatio-temporal patterns and characteristic features of seismic sequences. Swarms with alternating seismicity rates, bursts, and longer swarms with persistent moment release exhibit a remarkable diversity in temporal evolution. Swarm durations do not scale with cumulative seismic moment: swarms lasting less than 1 day can have similar cumulative seismic moments as month-long swarms. Here, we observe stationary swarms re-occurring for years (e.g., Mineral Mountains), as well as singular swarms in low-seismicity areas (e.g., activating a local structure). The swarms show a pronounced heterogeneity in triggering and driving mechanisms, observed in the detailed analysis of exemplary sequences (detections, relocations, moment tensors, waveform-based clustering, and repeater analysis). The 2022 Sevier Valley sequence activated a BR-related normal fault, the first resolved fault plane in the valley since 1983. The 2011 Circleville sequence is interpreted as a swarm triggered by mainshock-aftershock activity characterized by increasing magnitudes, changing rupture mechanisms, and a concentration of highly similar events in the second part of the sequence. By jointly discussing exemplary sequences and catalog statistics, we draw a comprehensive picture of swarm activity and its relation to geothermal and tectonic activity.

58 GEOSCIENCES↗

DOE OSTI2022

Calculating adsorption isotherms using the two-phase thermodynamic method and molecular dynamics simulations

We describe the calculation of adsorption isotherms from molecular dynamics simulations based on the two-phase thermodynamic (2PT) model. The 2PT model developed for bulk fluid phases treats the gas-like components as hard spheres (HSs), which correctly recovers the limiting behaviors of unconfined fluids. We showed that this treatment, however, does not always lead to the correct zero-loading behavior in strongly confining systems. For methane adsorption into zeolite MFI, the HS reference state underestimates entropy by up to 20% at low loadings and leads to an order-of-magnitude increase in the adsorption onset pressure. To fix these issues, we propose the use of ideal adsorbed gas (IAG) as the gas reference model, the properties of which can be computed using the Widom insertion method on an empty adsorbent. We further describe three routes to compute adsorption isotherms from the Helmholtz free energy at different loadings. Comparing against established Monte Carlo (MC) methods, we found that the adsorption isotherms obtained using the IAG reference state agrees to within 40%, which corresponds to deviations of <5% in adsorption free energy. The isotherms calculated using the HS reference state underestimate the adsorption uptake at low to medium loadings in strongly confining systems, but its accuracy improves at higher loadings and as the pore size increases relative to the sorbate diameter. The methods described here provide an alternative approach for computing adsorption isotherms when MC simulations in an open ensemble are undesirable and enable a direct comparison of computed adsorption thermodynamics with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Robust estimation of charge carrier diffusivity using transient photoluminescence microscopy

Transient microscopy has emerged as a powerful tool for imaging the diffusion of excitons and free charge carriers in optoelectronic materials. In many excitonic materials, extraction of diffusion coefficients can be simplified because of the linear relationship between signal intensity and local excited state population. However, in materials where transport is dominated by free charge carriers, extracting diffusivities accurately from multidimensional data is complicated by the nonlinear dependence of the measured signal on the local charge carrier density. To obtain accurate estimates of charge carrier diffusivity from transient microscopy data, statistically robust fitting algorithms coupled to efficient 3D numerical solvers that faithfully relate local carrier dynamics to raw experimental measurables are sometimes needed. Here, we provide a detailed numerical framework for modeling the spatiotemporal dynamics of free charge carriers in bulk semiconductors with significant solving speed reduction and for simulating the corresponding transient photoluminescence microscopy data. To demonstrate the utility of this approach, we apply a fitting algorithm using a Markov chain Monte Carlo sampler to experimental data on bulk CdS and methylammonium lead bromide (MAPbBr 3 ) crystals. Parameter analyses reveal that transient photoluminescence microscopy can be used to obtain robust estimates of charge carrier diffusivities in optoelectronic materials of interest, but that other experimental approaches should be used for obtaining carrier recombination constants. Additionally, simplifications can be made to the fitting model depending on the experimental conditions and material systems studied. Our open-source simulation code and fitting algorithm are made freely available to the scientific community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Filtering, averaging, and scale dependency in homogeneous variable density turbulence

We investigate relationships between statistics obtained from filtering and from ensemble or Reynolds-averaging turbulence flow fields as a function of length scale. Generalized central moments in the filtering approach are expressed as inner products of generalized fluctuating quantities, q ' ( ξ , x ) = q ( ξ ) - q ¯ ( x ) , representing fluctuations of a field q ( ξ ) , at any point ξ, with respect to its filtered value at x. For positive-definite filter kernels, these expressions provide a scale-resolving framework, with statistics and realizability conditions at any length scale. In the small-scale limit, scale-resolving statistics become zero. In the large-scale limit, scale-resolving statistics and realizability conditions are the same as in the Reynolds-averaged description. Using direct numerical simulations (DNS) of homogeneous variable density turbulence, we diagnose Reynolds stresses, T i j , resolved kinetic energy, kr, turbulent mass-flux velocity, a i , and density-specific volume covariance, b, defined in the scale-resolving framework. These variables, and terms in their governing equations, vary smoothly between zero and their Reynolds-averaged definitions at the small and large scale limits, respectively. At intermediate scales, the governing equations exhibit interactions between terms that are not active in the Reynolds-averaged limit. For example, in the Reynolds-averaged limit, b follows a decaying process driven by a destruction term; at intermediate length scales, it is a balance between production, redistribution, destruction, and transport, where b grows as the density spectrum develops and then decays when mixing becomes strong enough. This work supports the notion of a generalized, length-scale adaptive model that converges to DNS at high resolutions and to Reynolds-averaged statistics at coarse resolutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Exploring Si Heterojunction Solar Cell Degradation: Bulk and Interface Processes Analyzed by Simulations and Experiments in Order to Develop Mitigation Strategies

The Si technology with the highest conversion efficiency is a-Si/c-Si heterojunction (HJ) PV. Its market penetration, however, is slowed by reports that fielded HJ modules degraded at twice the rate of regular c-Si modules. Our very recent work confirmed this degradation and attributed it to enhanced recombination at the a-Si/c-Si interface, caused by the slow, order-of-magnitude increase of the defect density. We propose to comprehensively explore degradation mechanisms in Si HJ cells by combining simulations and experiments. Theoretically, we will: (1) simulate structure of a-Si and a-Si/c-Si interfaces, identify defects; (2) determine statistics of energy barriers that control defect formation; (3) compute growth of defect density from the barrier distribution, and the resulting degradation of Voc. Experimentally, we will: (1) create a series of HJ-cell-representative stacks with varying layer thicknesses and deposition conditions by using PECVD tools; (2) use temperature and injection-dependent lifetime spectroscopy to determine effective lifetimes; (3) deconvolve the data to separate bulk and interface effects, and the effects of charge density and interface defects to analyze long time degradation. The goal is to identify material and device degradation mechanisms, to develop mitigation strategies for improved stability, such as the introduction of capping layers and hydrogen diffusion control.

14 SOLAR ENERGY↗

DOE OSTI2019

Prediction of Probabilistic Detonation Threshold via Millimeter-Scale Microstructure-Explicit and Void-Explicit Simulations

We present an approach and relevant models for predicting the probabilistic shoc-to-detonation transition (SDT) behavior and Pop plot (PP) of heterogeneous energetic materials (HEM) via mesoscopic microstructure-explicit (ME) and void explicit (VE) simulations at the millimeter (mm) sample size scale. Although the framework here is general, the particular material considered in this paper is pressed Octahydro-1,3,5,7-tetranitro-1,2,3,5-tetrazocine (HMX). To systematically delineate the effects of material heterogeneities, four material cases are considered. These cases are homogeneous material, material with granular microstructure but no voids, homogeneous material with voids, and material with both granular microstructure and voids. Statistically equivalent microstructure sample sets (SEMSS) are generated and used. Eulerian hydrocode simulations explicitly resolve the material heterogeneities, voids, and the coupled mechanical-thermal-chemical processes. In particular, it is found that both microstructure and voids strongly influence the SDT behavior and PP. The effects of different combinations of microstructure heterogeneity and voids on the SDT process and PP are quantified and rank-ordered. The overall framework uses the Mie–Grüneisen equation of state and a history variable reactive burn model (HVRB). A novel probabilistic representation for quantifying the PP is developed, allowing the calculation of (1) the probability of observing SDT at a given combination of shock pressure and run distance, (2) the run-distance to detonation under a given combination of shock pressure and prescribed probability, and (3) the shock pressure required for achieving SDT at a given run distance with a prescribed probability. The results are in agreement with general trends in experimental data in the literature.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2021

Core–Multishell-Structured Digital-Gradient Cathode Materials with Enhanced Mechanical and Electrochemical Durability

Ni-rich cathode materials provide high energy density, but their structural and surface instability limits their cyclability and thermal stability. As one of the approaches to mitigate this problem, cathode materials comprising Ni-rich high-capacity core wrapped in Mn-rich multiple shells are produced successfully. In contrast to the conventional batch-type process for concentration-gradient materials, a digital-gradient cascade coprecipitation process described here achieves the improvements in productivity and quality consistency needed to move toward large-scale manufacturing. The core-multishell cathode materials produced in this manner not only have longer cycle life and improved rate performance compared to homogeneous Ni-rich cathode materials having the same overall composition, but also show remarkably enhanced thermal stability and low impedance growth characteristics. In a novel attempt to determine the correlation between the mechanical properties of the core-multishell cathode particles and their electrochemical cyclabilities, their breaking force and elasticity were successfully measured using a statistical approach, which indicates that a cathode particle with stable surface composition as well as high breaking force has improved capacity retention and durability. These results guide the realization of long life and high thermal stability in Ni-rich cathode materials through heterogeneous particle engineering.

25 ENERGY STORAGE↗

DOE OSTI2022

Imaging Phase Segregation in Nanoscale Li x CoO 2 Single Particles

Li x CoO 2 (LCO) is a common battery cathode material that has recently emerged as a promising material for other applications including electrocatalysis and as electrochemical random access memory (ECRAM). During charge– discharge cycling LCO exhibits phase transformations that are significantly complicated by electron correlation. While the bulk phase diagram for an ensemble of battery particles has been studied extensively, it remains unclear how these phases scale to nanometer dimensions and the effects of strain and diffusional anisotropy at the single-particle scale. Understanding these effects is critical to modeling battery performance and for predicting the scalability and performance of electrocatalysts and ECRAM. Here we investigate isolated, epitaxial LiCoO 2 islands grown by pulsed laser deposition. After electrochemical cycling of the islands, conductive atomic force microscopy (c-AFM) is used to image the spatial distribution of conductive and insulating phases. Above 20 nm island thicknesses, we observe a kinetically arrested state in which the phase boundary is perpendicular to the Li-planes; we propose a model and present image analysis results that show smaller LCO islands have a higher conductive fraction than larger area islands, and the overall conductive fraction is consistent with the lithiation state. Thinner islands (14 nm), with a larger surface to volume ratio, are found to exhibit a striping pattern, which suggests surface energy can dominate below a critical dimension. When increasing force is applied through the AFM tip to strain the LCO islands, significant shifts in current flow are observed, and underlying mechanisms for this behavior are discussed. The c-AFM images are compared with photoemission electron microscopy images, which are used to acquire statistics across hundreds of particles. Finally, the results indicate that strain and morphology become more critical to electrochemical performance as particles approach nanometer dimensions.

25 ENERGY STORAGE↗

DOE OSTI2021

Phase-field model of precipitation processes with coherency loss

A phase-field model is proposed to simulate coherency loss coupled with microstructure evolution. A special field variable is employed to describe the degree of coherency loss of each particle and its evolution is governed by a Ginzburg-Landau type kinetic equation. For the sake of computational efficiency, a flood-fill algorithm is introduced that can drastically reduce the required number of field variables, which allows the model to efficiently simulate a large number of particles sufficient for characterizing their statistical features during Ostwald ripening. The model can incorporate size dependence of coherency loss, metastability of coherent particles, and effectively incorporate the underlying mechanisms of coherency loss by introducing a so-called differential energy criterion. The model is applied to simulate coarsening of Al 3 Sc precipitates in aluminum alloy and comprehensively compared with experiments. Our results clearly show how the particle size distribution is changed during coherency loss and affects the coarsening rate.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Raman spectrum and polarizability of liquid water from deep neural networks

We introduce a scheme based on machine learning and deep neural networks to model the environmental dependence of the electronic polarizability in insulating materials. Application to liquid water shows that training the network with a relatively small number of molecular configurations is sufficient to predict the polarizability of arbitrary liquid configurations in close agreement with ab initio density functional theory calculations. In combination with a neural network representation of the interatomic potential energy surface, the scheme allows us to calculate the Raman spectra along 2-nanosecond classical trajectories at different temperatures for H 2 O and D 2 O. The vast gains in efficiency provided by the machine learning approach enable longer trajectories and larger system sizes relative to ab initio methods, reducing the statistical error and improving the resolution of the low-frequency Raman spectra. Decomposing the spectra into intramolecular and intermolecular contributions elucidates the mechanisms behind the temperature dependence of the low-frequency and stretch modes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Surface deformation coupled with self-organized pattern on a liquid anode of 1 atm DC glow discharge

Intricate, self-organized plasma structures observed above the surface of a liquid anode of atmospheric DC glow discharge were found to give rise to coherent, organized surface deformation and mechanical wave formation at the plasma–liquid interface. This new phenomenon indicates that the liquid is closely coupled to the plasma by the anode sheath’s electrohydrodynamic (EHD) force. A scientific question then arises: Do surface perturbations, coupled with the nonuniform surface charge distribution, enhance the electric field and induce self-organization? Using the reflective background-oriented schlieren technique, the liquid surface profile under the plasma pattern was measured for the first time. The results show that surface distortions are driven by the repulsive Coulomb force of nonuniform net-negative surface charge acted by the anode sheath field. The impacts of various operating parameters on the patterns and surface waves were examined, revealing the significance of gas heating and liquid charge relaxation time in the pattern formation mechanism. Time-resolved dynamics of a pulsed DC discharge indicated that the surface deformation only appeared after the establishment of plasma patterns. Statistically, the surface wave under the plasma has high wave numbers (8000–16000 m −1 ) and small amplitudes ($<$ 10 µm), generally found in the capillary wave regime. Yet the motion of surface deformations is in tandem with the plasma pattern and exhibits a nondispersive nature of constant phase velocity (0.1–0.4 m s −1 ), suggesting the dominant role of EHD force over the surface tension in the observed surface wave. These results indicate that the deformed liquid surface is driven by the EHD force of the plasma sheath. Although they share a similar geometry, the deformation and wave dynamics of the liquid surface do not stimulate a plasma pattern. Importantly, the complex EHD coupling in the plasma–liquid system raises awareness and new challenges for plasma control engineering, and the quantitative characteristics of the nondispersive surface wave are informative for advancing relevant theory and modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

(Anti-)deuteron production in pp collisions at $\sqrt{s}=13 \ \text {TeV}$

The study of (anti-)deuteron production in pp collisions has proven to be a powerful tool to investigate the formation mechanism of loosely bound states in high-energy hadronic collisions. In this paper the production of $\text {(anti-)deuterons}$ is studied as a function of the charged particle multiplicity in inelastic pp collisions at $\sqrt{s}=13$ TeV using the ALICE experiment. Thanks to the large number of accumulated minimum bias events, it has been possible to measure (anti-)deuteron production in pp collisions up to the same charged particle multiplicity ( ${\mathrm {d} N_{ch}/\mathrm {d} \eta } \sim 26$ ) as measured in p–Pb collisions at similar centre-of-mass energies. Within the uncertainties, the deuteron yield in pp collisions resembles the one in p–Pb interactions, suggesting a common formation mechanism behind the production of light nuclei in hadronic interactions. In this context the measurements are compared with the expectations of coalescence and statistical hadronisation models (SHM).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Synthesis of ARM User Facility Surface Rainfall Datasets to Construct a Best Estimate Value Added Product (PrecipBE)

Surface precipitation measurements are essential for Earth system model (ESM) evaluation and understanding cloud processes. An ever-growing need for robust, temporally evolving, and easy-to-use statistical datasets provides motivation for a baseline ground-based precipitation properties data product. The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility operates an extensive suite of precipitation instruments with various sensitivities and operating mechanisms, which render the decision of which instrument to use based on one or more fixed thresholds challenging and prone to errors and bias. Using a long-term instrument inter-comparison from a unique per-precipitation event perspective, rather than instantaneous sample comparison, we demonstrate that ARM rainfall-measuring instruments are generally consistent with each other at the statistical level. Inter-instrument deviations at the single event level can be large, especially for specific rainfall event properties such as maximum precipitation rates. A machine-learning (ML) analysis using a random forest regressor indicates that in some cases, depending on instrument, local site climatology, and/or specific deployment configuration, certain atmospheric state variables influence the measured quantities in an unpredictable manner. Thus, a-priori weighting of different instruments does not necessarily lead to more accurate and less biased synthesis of instrument data. These results motivate the design of the ARM precipitation best-estimate (PrecipBE) value-added product, which incorporates all valid precipitation data while considering data quality and other instrument limitations. PrecipBE consists of time series and tabular statistics datasets in an easy-to-use and insightful per-precipitation event format. It provides a large set of precipitation event properties supplemented with ancillary data from ARM datasets that correspond to the detected precipitation events. We describe the PrecipBE algorithm and demonstrate its use via the examination of a single-day output as well as a long-term trend analysis of precipitation events at the ARM Southern Great Plains (SGP) site, covering more than 30 years of data. The trend analysis tentatively suggests a long-term temporal tendency for mainly shorter and less intense precipitation events at the SGP site, but a long-term increase in annual rainfall by more than 36 mm (5 %) per decade. This rainfall trend is catalyzed primarily by more extreme event properties of relatively rare, intense precipitation events, with event total and 1 min maximum precipitation rate at a 1 year timeframe increasing up to 5 mm and 9 mm h −1 (several percent) per decade, respectively. While the currently available PrecipBE datasets (at https://adc.arm.gov/discovery/, last access: 8 December 2025) cover rainfall from multiple ARM deployments up to March 2025, PrecipBE is planned to be expanded to include solid-phase precipitation and will soon become an operational product with a several-day lag from real-time. We invite the ARM user community to leverage this new product and welcome user feedback to further enhance the dataset.

Silber, Israel [Pacific Northwest National Laborat↗

DOE OSTI2021

Unraveling the discontinuous plastic flow of a Co-Cr-Fe-Ni-Mo multiprincipal-element alloy at deep cryogenic temperatures

We report an analysis of the discontinuous plastic flow of a multiprincipal-element alloy, Co 17.5 Cr 12.5 Fe 55 Ni 10 Mo 5 (atomic percent, at. %), in the temperature range of 0.5–4.2 K showing serrated deformation curves. Using the analytical techniques, we studied the statistics of the stress drops associated with the unstable plastic flow. The analysis showed that the complexity and heterogeneity of a discontinuous plastic flow were reduced when the temperature was lowered. This behavior was associated with the effects of dynamic recovery and adiabatic heating on the dislocation-density evolution.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Systems Analysis of the Physiological and Molecular Mechanisms of Sorghum Nitrogen Use Efficiency, Water Use Efficiency and Interactions with the Soil Microbiome (Final Report for DE-SC0014395)

The specific project objectives were to: 1) Conduct deep census surveys of root microbiomes concurrent with phenotypic characterizations of a diverse panel of sorghum genotypes across multiple years to define the microbes associated with the most productive lines under drought and low nitrogen conditions. 2) Associate systems-level genotypic, microbial, and environmental factors with improved sorghum performance using robust statistical approaches. 3) Develop culture collections of sorghum root/leaf associated microbes that recapitulate root-enriched sequences defined in the census. 4) Perform controlled environment experiments for in-depth characterization and hypothesis testing of G sorghum x G microbe x E interactions . Validate physiological mechanisms, map genetic loci for stress tolerance, and determine the persistence of optimal microbial strains under greenhouse and field conditions.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Quantitative analysis of Raman spectral parameters for carbon fibers: practical considerations and connection to mechanical properties

Although the literature on the Raman spectra of carbon fibers is vast, no consistent, robust predictive relationship between mechanical properties of carbon fibers and spectral parameters exists. This shortcoming is due to the use of numerous fitting functions to evaluate Raman spectra of carbon fibers and the inconsistencies in establishing the best fitting models in a statistically robust fashion. To address this gap, we present a comprehensive work on the Raman spectra of carbon fibers that combines a vast library of experimental data with a robust numerical analysis and a statistical evaluation of a wide range of suggested fitting models. This manuscript begins with a brief review of the commonly applied fitting models. Then, the Raman spectra of 32 commercially available polyacrylonitrile-based carbon fibers collected at excitation wavelengths 532, 633, and 785 nm are presented and the best fit for all fibers is evaluated based on several statistical criteria in conjunction with numerical calculations and physical arguments. The results suggest that high-performance fibers must be fit with at least five peaks, whereas high-tensile modulus fibers are best fit with at least six distinct peaks. Finally, we employ simultaneous fitting of the Raman spectra of specific fibers and wavelengths and demonstrate that strong correlations exist between mechanical properties and the D1 peak position and shape across the range of evaluated mechanical properties. Finally, we suggest straightforward improvements in fitting analysis procedures that can be implemented to increase coherency in the understanding of the underlying carbon fiber microstructure intuited from Raman spectroscopy.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Scaling microstructural processes in the sintering of ionic ceramics

A multi-scale framework, combining a multiphase field formulation and large deformation mechanics, was developed as a stepping stone to perform the data analytics of the microstructural level kinetics of a sintering solid. Relevant microstructural information from this framework, such as grain, stress, and porosity statistics, was scaled up to describe the macroscopic level sintering kinetics. Here, the developed formulation was applied to describe the electric field assisted sintering of Y 2 O 3 . Microstructural inhomogeneities in a multi-granular solid result in the formation of a field of compressive stress networks, which interleave with low compression and weakly tensile regions, defining a scaffolding for sintering concentration regions to develop. A Poisson effect-induced lateral stress network is also naturally self-induced as a result of the mechanical constraints imposed by the sintering apparatus. For long sintering times, localized shear stresses enhancing mass flow along grain boundaries and internal surfaces develop. Three-sided pores are removed by either vacancy transport to the surrounding pores, or move towards the external surfaces through grain boundary diffusion. Four- and higher order-sided pores stabilize because an equal amount of vacancies are gained and lost through the connecting grain boundaries. Grain dewetting contributes to pore coalescence, suggesting that pore kinetics and grain growth are coupled and should be analyzed in concert. The combined sintering and grain growth kinetics define six regimes of sintering behavior: (1) T, the transient regime; (2) E$_Υ$, the surface energy dominated, early sintering regime, where the grain growth exponent, p = 1, and the stress concentration factor, $f$ ~ $1/\hat{ρ}^{4.6}$; (3) E S , the stress dominated, early sintering regime, where p = 1 and $f$ ~ $1/\hat{ρ}^{4}$; (4) I$_Υ$, the surface energy dominated, intermediate sintering regime, where p = 2 and $f$ ~ $1/\hat{ρ}^{4.6}$; (5) I S , the stress dominated, intermediate sintering regime, where p = 2 and $f$ ~ $1/\hat{ρ}^{4}$; and (6) L, the late sintering regime, where p = 3 and $f$ ~ 1. At the macroscopic level, the rapid densification and suppression of grain growth observed in the electric field assisted sintering process is a consequence of the compounding effects of the underlying stress-, transport-, and interfacial-energy-induced energy minimization kinetics, as predicted by the multi-scale framework.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗