Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion Couple Experiments”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Localizing Clinical Patterns of Blast Traumatic Brain Injury Through Computational Modeling and Simulation

Blast traumatic brain injury is ubiquitous in modern military conflict with significant morbidity and mortality. Yet the mechanism by which blast overpressure waves cause specific intracranial injury in humans remains unclear. Reviewing of both the clinical experience of neurointensivists and neurosurgeons who treated service members exposed to blast have revealed a pattern of injury to cerebral blood vessels, manifested as subarachnoid hemorrhage, pseudoaneurysm, and early diffuse cerebral edema. Additionally, a seminal neuropathologic case series of victims of blast traumatic brain injury (TBI) showed unique astroglial scarring patterns at the following tissue interfaces: subpial glial plate, perivascular, periventricular, and cerebral gray-white interface. The uniting feature of both the clinical and neuropathologic findings in blast TBI is the co-location of injury to material interfaces, be it solid-fluid or solid-solid interface. This motivates the hypothesis that blast TBI is an injury at the intracranial mechanical interfaces. In order to investigate the intracranial interface dynamics, we performed a novel set of computational simulations using a model human head simplified but containing models of gyri, sulci, cerebrospinal fluid (CSF), ventricles, and vasculature with high spatial resolution of the mechanical interfaces. Simulations were performed within a hybrid Eulerian—Lagrangian simulation suite (CTH coupled via Zapotec to Sierra Mechanics). Because of the large computational meshes, simulations required high performance computing resources. Twenty simulations were performed across multiple exposure scenarios—overpressures of 150, 250, and 500 kPa with 1 ms overpressure durations—for multiple blast exposures (front blast, side blast, and wall blast) across large variations in material model parameters (brain shear properties, skull elastic moduli). All simulations predict fluid cavitation within CSF (where intracerebral vasculature reside) with cavitation occurring deep and diffusely into cerebral sulci. These cavitation events are adjacent to high interface strain rates at the subpial glial plate. Larger overpressure simulations (250 and 500kPa) demonstrated intraventricular cavitation—also associated with adjacent high periventricular strain rates. Additionally, models of embedded intraparenchymal vascular structures—with diameters as small as 0.6 mm—predicted intravascular cavitation with adjacent high perivascular strain rates. The co-location of local maxima of strain rates near several of the regions that appear to be preferentially damaged in blast TBI (vascular structures, subpial glial plate, perivascular regions, and periventricular regions) suggest that intracranial interface dynamics may be important in understanding how blast overpressures leads to intracranial injury.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Effective kinetics driven by dynamic concentration gradients under coupled transport and reaction

Biogeochemical reaction kinetics are generally established from batch reactors where concentrations are uniform. In natural systems, many biogeochemical processes are characterized by spatially and temporally variable concentration gradients that often occur at scales which are not resolved by field measurements or biogeochemical and reactive transport models. Yet, it is not clear how these sub-scale chemical gradients affect reaction kinetics compared to batch kinetics. Here we investigate this question by studying the paradigmatic case of localized pulses of solute reacting with a solid or a dissolved species in excess. Additionally, we consider non-linear biogeochemical reactions, representative of mineral dissolution, adsorption and redox reactions, which we quantify using simplified power-law kinetics. The combined effect of diffusion and reaction leads to effective kinetics that differ quantitatively and qualitatively from the batch kinetics. Depending on the nonlinearity (reaction order) of the local kinetics, these effects lead to either enhancement or decrease of the overall reaction rate, and result in a rich variety of reaction dynamics. We derive analytical results for the effective kinetics, which are validated by comparison to direct numerical simulations for a broad range of Damköhler numbers and reaction order. Our findings provide new insights into the interpretation of imperfectly mixed lab experiments, the effective kinetics of field systems characterized by intermittent reactant release and the integration of sub-scale concentration gradients in reactive transport models.

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

Investigation of Drift Effects in UEDGE Simulations of NSTX-U Edge Plasma With Lithium Divertors

Lithium is a low-Z material, and lithium-based plasma-facing components (PFCs) are planned for the National Spherical Torus Experiment Upgrade (NSTX-U) to explore potential benefits for divertor power exhaust and core plasma management. NSTX-U is a medium-sized spherical tokamak with up to 12 MW of auxiliary heating, capable of generating reactor-relevant plasma conditions. This work presents boundary plasma simulations for NSTX-U with lithium PFCs using the UEDGE code, incorporating full magnetic and 𝐄 ×𝐁 drift physics. The simulations show that drifts strongly influence heat and particle transport: they enhance convective transport, broaden the scrape-off layer heat-flux width 𝜆 𝑞 , and reduce the anomalous heat diffusivity 𝜒 required to reproduce predicted SOL heat-flux width. 𝐄 ×𝐁 drifts provide poloidal transport, while ∇𝐵 (which includes both gradB and curvature) drifts provide radial heat and particle transport. Lithium transport is also affected by drifts, with lithium ions migrating from the outer divertor to the inner divertor through the private flux region (PFR) following the 𝐄 ×𝐁 drifts flow, lowering upstream impurity lithium densities. UEDGE is self-consistently coupled with the Wall-Li model to study plasma lithium PFC interactions depending on the local lithium sourcing based on local plasma conditions and lithium surface temperature. In these simulations, lithium evaporation shows a vapor-shielding effect that reduces divertor heat flux and increases radiative losses once surface temperatures exceed 450°C. This research work provides a first step toward self-consistent modeling of lithium PFCs in NSTX-U, demonstrating the impact of drift-driven plasma transport in SOL and divertor regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2002

Pulsed Electromagnetic Acceleration of Plasmas

A major shift in paradigm in driving pulsed plasma thruster is necessary if the original goal of accelerating a plasma sheet efficiently to high velocities as a plasma "slug" is to be realized. Firstly, the plasma interior needs to be highly collisional so that it can be dammed by the plasma edge layer not (upstream) adjacent to the driving 'vacuum' magnetic field. Secondly, the plasma edge layer needs to be strongly magnetized so that its Hall parameter is of the order of unity in this region to ensure excellent coupling of the Lorentz force to the plasma. Thirdly, to prevent and/or suppress the occurrence of secondary arcs or restrike behind the plasma, the region behind the plasma needs to be collisionless and extremely magnetized with sufficiently large Hall parameter. This places a vacuum requirement on the bore conditions prior to the shot. These requirements are quantified in the paper and lead to the introduction of three new design parameters corresponding to these three plasma requirements. The first parameter, labeled in the paper as gamma (sub 1), pertains to the permissible ratio of the diffusive excursion of the plasma during the course of the acceleration to the plasma longitudinal dimension. The second parameter is the required Hall parameter of the edge plasma region, and the third parameter the required Hall parameter of the region behind the plasma. Experimental research is required to quantify the values of these design parameters. Based upon fundamental theory of the transport processes in plasma, some theoretical guidance on the choice of these parameters are provided to help designing the necessary experiments to acquire these data.

Thio, Y. C. Francis↗
DOE OSTI2024

Phase field-volumetric lattice Boltzmann model of ion uptake in porous nuclear waste form materials under continuous flow

The flow field within the mesopores of sorbent particles plays a crucial role in radionuclide diffusion and ion uptake kinetics, thus, impacting the overall performance of porous nuclear waste form materials. To fundamentally understand the influence of microstructures and material properties on the radionuclide absorption and retention processes requires a coupled multi-physics model that considers the advection and diffusion within the flow field, the reaction at liquid-solid interfaces, and finally, the solid-state diffusion within a complex nanoporous medium. Here, this study employs the volumetric lattice Boltzmann method (VLBM) to accurately and efficiently calculate the steady state velocity field inside the mesopores of sorbent particles. The obtained velocity field is then utilized to calculate the advection of ions in the steady flow. A phase field (PF) model of ion uptake is used to describe the reaction occurring at the solid-liquid interface and diffusion inside the porous medium. The integrated PF-VLBM model is verified in terms of the mass conservation and numerical efficiency and validated qualitatively with experimental observation data. Then, it is applied to study the influence of thermodynamic and kinetic properties, as well as flow field conditions on the ion uptake kinetics. The numerical results demonstrate that the ion uptake kinetics in porous particles has three distinct stages, which is in agreement with the observations in continuous flow experiments. In the first stage, the kinetics is predominantly controlled by the flow field and ion diffusivity in the liquid phase. The kinetics in the second stage is primarily governed by ion diffusivity in the solid phase. In the third stage the system reaches a dynamic equilibrium with a net zero uptake flux at the interface. It is also found that porous structures significantly affect the efficiency and capacity of ion uptake. The simulation results can help to understand the physics behind the observed ion uptake kinetics in experiments and to facilitate the development of constitutive equations that can account for heterogeneous microstructures in engineering performance codes.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Seismic anisotropy, dominant slip systems and phase transitions in the lowermost mantle

SUMMARY The presence of seismic anisotropy at the base of the Earth's mantle is well established, but there is no consensus on the deformation mechanisms in lower mantle minerals that could explain it. Strong anisotropy in magnesium post-perovskite (pPv) has been invoked, but different studies disagree on the dominant slip systems at play. Here, we aim to further constrain this by implementing the most recent results from atomistic models and high-pressure deformation experiments, coupled with a realistic composition and a 3-D geodynamic model, to compare the resulting deformation-induced anisotropy with seismic observations of the lowermost mantle. We account for forward and reverse phase transitions from bridgmanite (Pv) to pPv. We find that pPv with either dominant (001) or (010) slip can both explain the seismically observed anisotropy in colder regions where downwellings turn to horizontal flow, but only a model with dominant (001) slip matches seismic observations at the root of hotter large-scale upwellings. Allowing for partial melt does not change these conclusions, while it significantly increases the strength of anisotropy and reduces shear and compressional velocities at the base of upwellings.

58 GEOSCIENCES↗
DOE OSTI2021

Efficient First-Principles Methodology for the Calculation of the All-Phonon Inelastic Scattering in Solids

Inelastic scattering experiments are key methods for mapping the full dispersion of fundamental excitations of solids in the ground as well as nonequilibrium states. A quantitative analysis of inelastic scattering in terms of phonon excitations requires identifying the role of multiphonon processes. Here, we develop an efficient first-principles methodology for calculating the all-phonon quantum mechanical structure factor of solids. We demonstrate our method by obtaining excellent agreement between measurements and calculations of the diffuse scattering patterns of black phosphorus, showing that multiphonon processes play a substantial role. The present approach constitutes a step towards the interpretation of static and time-resolved electron, x-ray, and neutron inelastic scattering data.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Coupling between interfacial charge and mechanical deformation at high temperatures in ceramics. Final report

Basic research, under this project, has led to new fundamental discoveries in the effect of electric field (and temperature) on diffusion kinetics, electroluminescence and electronic conductivity in ceramic materials. For example, yttria stabilized zirconia sinters in just seconds well below 1000°C, while, at the same time, exhibiting a highly non-linear rise in conductivity, and intense electroluminescence. The phenomenon has been shown to occur in nearly all ceramics whether they are ionic or electronic conductors, or insulators. The rise in conductivity leads to Joule heating but in-situ measurements with a platinum standard at APS and Brookhaven (NSLS-II) have shown that the temperatures are much below what would have been needed to sinter in a few seconds. We have explored the hypothesis that huge concentrations of Frenkel pairs of vacancies and interstitials are produced. While these point defects accelerate sintering (by the interstitials being transported to the pores and the vacancies to the grain boundaries), they may also produce electron-hole pairs which can impart high electrical conductivity or recombine to emit photons. Many measurements, including in-situ measurements of lattice expansion are consistent with this hypothesis. However, the energy barrier for generating Frenkels in zirconia for example is much larger than can be expected at temperature and electrical fields in these experiments. Most recent work suggests that short wavelength phonons play a critical role in the genesis of the “flash” phenomenon. For example, the Debye temperature has been shown to be the lower bound for initiating flash. Molecular dynamics simulations have shown that large concentrations of Frenkel pairs can form under the proliferation of “optical” phonons at the edge of the Brillouin zone, that is with a wavelength at or below the lattice parameter. This project has led to the development of new experimental techniques as well as a software hardware interface for process and control and data analysis in real time, with a time scale of less than one second. This approach spells a new paradigm in materials science research. It potential impact on high rate manufacturing is self evident.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A model of thermal creep and annealing in finite domains based on coupled dislocation climb and vacancy diffusion

Here, we develop a framework to investigate thermal creep and annealing in finite domains, where the climb motion of discrete dislocations is coupled to the diffusion of a continuum vacancy field. The model is first formulated in a continuum finite-deformation setting. All governing equations and boundary conditions are obtained from a unified irreversible thermodynamics principle. The resulting model couples a mechanical boundary value problem (BVP), a vacancy diffusion BVP, and the climb and glide motion of the discrete dislocation network within the crystal. The framework is then linearized for implementation in three-dimensional (3D) discrete dislocation dynamics (DDD) simulations for arbitrary anisotropic crystals. A solution scheme is developed based on the superposition principle, which is imposed weakly on the dislocation network to obtain a Galerkin solution for the nodal climb velocities. The framework includes diffusional (Nabarro–Herring) creep deformation as well as dislocation creep by climb-assisted-glide. The method is applied to simulate the annealing of vacancy loops in Al, with good agreement to experimental measurements by Silcox and Hirsch. We further consider the effects of annealing under stress, and of the proximity of the vacancy loops to loaded and free boundaries Simulations in polycrystalline materials are carried out to highlight the effects of the grain size on dislocation climb and vacancy loop annealing. The method is also applied to estimate the creep rate due to climb-assisted glide of jogged-screw dislocations in γ-TiAl, and results are compared to experiments by Viswanathan et al. Finally, we discuss the effects of uniaxial and hydrostatic stresses on the two diffusive deformation pathways of the material, namely Nabarro–Herring creep and dislocation climb.

36 MATERIALS SCIENCE↗
NASA NTRS2023

Microstructural Changes and Pb Mobility During the Zircon to Reidite Transformation: Implications for Planetary Impact Chronology

Impact events modify and leave behind a complex history of rock metamorphism on terrestrial planets. Evidence for an impact event may be recorded in physical changes to minerals, such as mineral deformation and formation of high P-T polymorphs, but also in the form of chemical fingerprints, such as enhanced elemental diffusion and isotopic mixing. Here we explore laboratory shock-induced physical and chemical changes to zircon and feldspar, the former of which is of interest because its trace elements abundances and isotope ratios are used extensively in geochemistry and geochronology. To this end, a granular mixture of Bishop Tuff sanidine and Kuehl Lake zircon, both with well characterized Pb isotope compositions, was prepared and then shocked via a flat plate accelerator. The peak pressure of the experiment, as calculated by the impedance matching method, was ~24 GPa although a broader range of P-T conditions is anticipated due to starting sample porosity. Unshocked and shocked materials were characterized via scanning electron microscopy (SEM), electron backscatter diffraction (EBSD), and Raman spectroscopy. These methods show that the starting zircon material had abundant metamict regions, and the conversion of the feldspar to glass in the post-shock material. Analyses of the shocked product also yielded multiple occurrences of the high-pressure ZrSiO 4 polymorph reidite, with some domains up to 300 μm across. The possibility of U-Pb system disturbance was evaluated via laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) and secondary ion mass spectrometry (SIMS). The isotopic data reveal that disturbance of the U-Pb geochronometer in the reidite was minimal (<2% for the main U-Pb geochronometers). To better constrain the P-T conditions during the shock experiment, we complement impedance matching pressure calculations with iSALE2D impact simulations. The simulated results yield a range of P-T conditions experienced during the experiment and show that much of the sample may have reached >30 GPa, which is consistent with formation of reidite. In the recovered shocked material, we identified lamellae of reidite, some of which interlock with zircon lamellae. Reidite {112} twins were identified, which we interpret to have formed to reduce stress between the crystal structure of the host zircon and reidite. These two findings support the interpretation that shear transformation enabled the transition of zircon to reidite. The size and presence of reidite found here indicate that this phase is probably common in impact-shocked crustal rocks that experienced ~25 to ~35 GPa, especially when the target material has porosity. Additionally, shock loading of the zircon and transformation to reidite at these pressures in porous materials is unlikely to significantly disturb the U-Pb system in zircon and that the reidite inherits the primary U and Pb elemental and isotopic ratios from the zircon.

Bishop Tuff↗
DOE OSTI2025

The role of sample–window coupling in shock temperature experiments

Radiometric studies of opaque materials under shock compression probe a sample–window interface that may be very different than the bulk state. Subtle experimental details, particularly an initial sample–window gap vs an adhesively bonded stack, have tremendous influence on measured temperature. This work experimentally tests gapped vs bonded targets and theoretically investigates observed differences in the measured temperature. The results indicate that perturbations created by a thin adhesive are smaller than those suggested by a previous work, whereas sample–window gaps alter the measured state more significantly than that previously realized.

Adhesives↗
DOE OSTI2022

Fast redox switches lead to rapid transformation of goethite in humid tropical soils: A Mössbauer spectroscopy study

Abstract Humid tropical forest soils experience frequent rainfall, which limits oxygen diffusion and creates redox heterogeneity in upland soils. In this study we gauged the effect of short‐term anoxic conditions on Fe mineralogy of relatively Fe‐ and C‐rich surface soils (C/Fe mole ratio ∼5) from a humid tropical forest in Puerto Rico. Soils subjected to 4‐d oxic/anoxic oscillation were characterized by selective chemical extractions, Mössbauer spectroscopy (MBS), and X‐ray diffraction. Chemical extraction data suggested that rapidly switching redox conditions had subtle effects on bulk Fe mineralogy. Mössbauer, on the other hand, indicated that (a) the soil Fe is a mixture of goethites of varying characteristics with minor contributions from ferrihydrite (<5%) and Fe(III)‐organic matter (OM) phases (∼10%), and (b) anoxic conditions rapidly transformed all forms of goethite to relatively stable Fe oxides. Such fast changes in goethite features could be due to rapid depletion and sorption of bio‐released structural Al or sorbed OM onto residual soil components. The rapid temporal changes in MBS parameters and corresponding pore water nominal oxidation state of C values suggest that Fe‐C transformations in these upland tropical soils are complex and intricately coupled. A comprehensive understanding of the fate of Fe, Al, and OM (Fe‐organic moieties) during redox switches and concurrent changes in pore water chemistry is critical for development of robust transport models in humid tropical soils, which are subject to episodic low‐redox events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Low-Order Preconditioning for the High-Order Finite Element de Rham Complex

Here, we present a unified framework for constructing spectrally equivalent low-order-refined discretizations for the high-order finite element de Rham complex. This theory covers diffusion problems in H 1 , H(curl), and H(div) and is based on combining a low-order discretization posed on a refined mesh with a high-order basis for Nédélec and Raviart–Thomas elements that makes use of the concept of polynomial histopolation (polynomial fitting using prescribed mean values over certain regions). This spectral equivalence, coupled with algebraic multigrid methods constructed using the low-order discretization, results in highly scalable matrix-free preconditioners for high-order finite element problems in the full de Rham complex. Additionally, a new lowest-order (piecewise constant) preconditioner is developed for high-order interior penalty discontinuous Galerkin (DG) discretizations, for which spectral equivalence results and convergence proofs for algebraic multigrid methods are provided. In all cases, the spectral equivalence results are independent of polynomial degree and mesh size; for DG methods, they are also independent of the penalty parameter. These new solvers are flexible and easy to use; any “black-box” preconditioner for low-order problems can be used to create an effective and efficient preconditioner for the corresponding high-order problem. A number of numerical experiments are presented, based on an implementation in the finite element library MFEM. A range of challenging three-dimensional problems are used to corroborate the theoretical properties and demonstrate the flexibility and scalability of the method.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1996

Extinguishment of a Diffusion Flame Over a PMMA Cylinder by Depressurization in Reduced-Gravity

Extinction of a diffusion flame burning over horizontal PMMA (Polymethyl methacrylate) cylinders in low-gravity was examined experimentally and via numerical simulations. Low-gravity conditions were obtained using the NASA Lewis Research Center's reduced-gravity aircraft. The effects of velocity and pressure on the visible flame were examined. The flammability of the burning solid was examined as a function of pressure and the solid-phase centerline temperature. As the solid temperature increased, the extinction pressure decreased, and with a centerline temperature of 525 K, the flame was sustained to 0.1 atmospheres before extinguishing. The numerical simulation iteratively coupled a two-dimensional quasi-steady, gas-phase model with a transient solid-phase model which included conductive heat transfer and surface regression. This model employed an energy balance at the gas/solid interface that included the energy conducted by the gas-phase to the gas/solid interface, Arrhenius pyrolysis kinetics, surface radiation, and the energy conducted into the solid. The ratio of the solid and gas-phase conductive fluxes Phi was a boundary condition for the gas-phase model at the solid-surface. Initial simulations modeled conditions similar to the low-gravity experiments and predicted low-pressure extinction limits consistent with the experimental limits. Other simulations examined the effects of velocity, depressurization rate and Phi on extinction.

Goldmeer, Jeffrey Scott↗
DOE OSTI2026

In situ catalyst activation and regeneration enable energy-efficient high-current CO 2 reduction to ethanol-rich C 2+ mixtures

Electrochemical conversion of dissolved CO 2 in bicarbonate electrolytes, i.e., bicarbonate electrolysis, offers distinct advantages over gas diffusion electrode systems by enabling direct utilization of the CO 2 capture electrolyte while bypassing the energy-intensive CO 2 release step. However, bicarbonate electrolysis faces challenges such as CO 2 mass-transfer limitation, local pH-driven CO 2 depletion, and high cathodic potentials. The higher potential often causes catalyst surface reorganization, leading to a gradual loss of active sites and variations in selectivity during CO 2 reduction. Here, we report a directed, in situ activation and regeneration method that allows precatalysts to equilibrate under dynamic (pulsed) electrolysis conditions. We demonstrate in situ activation of a scalable Cu 2 O/Cu mesh that, under short-width (t = 4 s) pulsed electrolysis, provides stable mixed oxidation states of Cu, favoring the formation of an ethanol-rich crude mixture. The pulsed electrolysis waveform, consisting of six distinct segments, is tuned to form Cu + oxides, which are then reduced to generate local alkaline conditions favoring C–C coupling. This synergistic effect results in FEs of 73% for C2+ products and 39% for ethanol at an applied current density of −150 mA cm −2 and a cathodic potential of −1.45 V (vs. RHE). The overall half-cell energy efficiency is ∼30% for C 2+ products. The in situ Raman experiments confirm the role of pCO 2 R in dynamically regenerating Cu+-containing surface species during pulsed operation, thereby steering selectivity towards C 2+ products. A comprehensive multiscale, multiphysics model is developed to investigate the dynamic behavior of copper surface species (Cu, Cu + , and Cu 2+ ) and local microenvironmental conditions during the pCO 2 R. The results reveal that the coexistence of different copper oxidation states, especially the Cu+ intermediate, is critical in steering selectivity towards multicarbon (C 2+ ) products. The dynamic modulation of surface redox states via tailored pulsing strategies favors C–C coupling pathways by inducing localized alkaline conditions and stabilizing reactive intermediates. This work establishes a predictive modeling platform that links pulse waveform design with mechanistic insights into catalyst state evolution and product selectivity. Overall, this study provides valuable insights into the synergistic effect of in situ activation of pre-catalysts and pulsed electrolysis for higher selectivity towards C 2+ products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1996

Soot and Radiation Measurements in Microgravity Jet Diffusion Flames

The subject of soot formation and radiation heat transfer in microgravity jet diffusion flames is important not only for the understanding of fundamental transport processes involved but also for providing findings relevant to spacecraft fire safety and soot emissions and radiant heat loads of combustors used in air-breathing propulsion systems. Our objectives are to measure and model soot volume fraction, temperature, and radiative heat fluxes in microgravity jet diffusion flames. For this four-year project, we have successfully completed three tasks, which have resulted in new research methodologies and original results. First is the implementation of a thermophoretic soot sampling technique for measuring particle size and aggregate morphology in drop-tower and other reduced gravity experiments. In those laminar flames studied, we found that microgravity soot aggregates typically consist of more primary particles and primary particles are larger in size than those under normal gravity. Comparisons based on data obtained from limited samples show that the soot aggregate's fractal dimension varies within +/- 20% of its typical value of 1.75, with no clear trends between normal and reduced gravity conditions. Second is the development and implementation of a new imaging absorption technique. By properly expanding and spatially-filtering the laser beam to image the flame absorption on a CCD camera and applying numerical smoothing procedures, this technique is capable of measuring instantaneous full-field soot volume fractions. Results from this technique have shown the significant differences in local soot volume fraction, smoking point, and flame shape between normal and reduced gravity flames. We observed that some laminar flames become open-tipped and smoking under microgravity. The third task we completed is the development of a computer program which integrates and couples flame structure, soot formation, and flame radiation analyses together. We found good agreements between model predictions and experimental data for laminar and turbulent flames under both normal and reduced gravity. We have also tested in the laboratory the techniques of rapid-insertion fine-wire thermocouples and emission pyrometry for temperature measurements. These techniques as well as laser Doppler velocimetry and spectral radiative intensity measurement have been proposed to provide valuable data and improve the modeling analyses.

Ku, Jerry C.↗
DOE OSTI2026

Effects of temperature and dose rate on ion-irradiated γ-LiAlO 2 pellets

Defect accumulation and microstructural evolution during ion irradiation at elevated temperatures are governed by competing processes of defect production, driven by the dose rate, and defect recovery, controlled by diffusion, interaction, and annihilation. Here, this study investigates the effects of irradiation temperature and the dose rate on microstructural evolution, deuterium retention, and lithium volatilization in γ-LiAlO 2 pellets subjected to sequential He + and D + ion irradiation. Experiments were performed to a total fluence of 3 × 10 17 (He + + D + )/cm 2 at 623, 673, 723, and 773 K with an average He + dose rate of 7.7 × 10 −4 dpa/s, and to 2 × 10 17 (He + + D + )/cm 2 at 773 K with dose rates of 6.8 × 10 −5 , 2.9 × 10 −4 , and 7.3 × 10 −4 dpa/s. At 623 K, the microstructure was dominated by cavities and fractures with no observable precipitate formation, while small precipitates emerged at 673 K. Increasing the irradiation temperature to 723–773 K promoted the formation of larger, faceted LiAl 5 O 8 precipitates, and surface amorphization, accompanied by pronounced lithium depletion and H–D isotopic exchange. At 773 K, medium and high dose rates produced an amorphized surface layer over a crystalline subsurface containing LiAl 5 O 8 precipitates and blisters at the crystalline–amorphous interface, whereas low-dose-rate irradiation preserved surface crystallinity with cavities distributed in the matrix, around precipitates, and along grain boundaries. Precipitate morphology was anisotropic with limited size dependence on the dose rate. These results elucidate the coupled effects of temperature and the dose rate and demonstrate that sequential He + and D 2 + irradiation at 773 K reproduces key microstructural features and H isotope behavior observed in neutron-irradiated γ-LiAlO 2 at 573 K.

dose rate effects↗
DOE OSTI2023

Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies

This report represents the milestone deliverable M2SF-23SN010301072 “Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies” The report provides a status update of FY23 activities for the work package Argillite Disposal work packages for the DOE-NE Spent Fuel Waste Form Science and Technology (SFWST) Program. Clay-rich geological media (often referred as shale or argillite) are among the most abundant type of sedimentary rock near the Earth’s surface. Argillaceous rock formations have the following advantageous attributes for deep geological nuclear waste disposal: widespread geologic occurrence, found in stable geologic settings, low permeability, self-sealing properties, low effective diffusion coefficient, high sorption capacity, and have the appropriate depth and thickness to host nuclear waste repository concepts. The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress (through experiment, modeling, and testing) in the study of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in clay/shale/argillaceous rock. International collaboration activities comprising field-scale heater tests, field data monitoring, and laboratory-scale experiments provide key information on changes to the engineered barrier system (EBS) material exposed high thermal loads. Moreover, consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the engineered barrier system (EBS) design concepts. Chemical and structural analyses of sampled bentonite material from laboratory tests at elevated temperatures are key to the characterization of thermal effects affecting bentonite clay barrier performance. The knowledge provided by these experiments is crucial to constrain the extent of sacrificial zones in the EBS design during the thermal period. Thermal, hydrologic, mechanical, and chemical (THMC) data collected from heater tests and laboratory experiments have been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches to assess issues on coupled processes involving porous media flow, transport, geomechanical phenomena, chemical interactions with barrier/geologic materials, and the development of EBS concepts. These lines of knowledge are central to the design of deep geological backfilled repository concepts where temperature plays a key role in the EBS behavior, potential interactions with host rock, and long-term performance in the safety assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗