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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Development of SAM Code Capabilities for Safety Analysis of GCR Air-ingress Events

Air-ingress following a depressurized loss-of-forced-cooling (DLOFC) event is a challenging, multiphysics safety scenario for High-Temperature Gas-Cooled Reactors (HTGRs), involving coupled gas composition transport, buoyancy-driven flow redistribution, graphite oxidation, and structural heat-up. Despite its importance — air ingress is a key scenario identified in the PIRT process for the HTGRs — existing system-level safety codes have lacked the integrated capability to simulate the complete event sequence with high confidence. This report documents the development, validation, and demonstration of three new capabilities in the SAM code to address this gap: (1) a multi-component gas mixture flow model with binary diffusion to track the helium-air composition and its effect on system density and flow; (2) a 0-D graphite oxidation model based on the Roes correlation, including oxygen consumption and exothermic heat release; and (3) an isentropic critical flow model for accurate representation of primary system depressurization through a break. These capabilities are validated against two benchmark experiments. The NSTF heavy-gas ingress experiment validates the multi-component flow model: SAM correctly reproduces the rapid buoyancydriven flow stagnation and subsequent natural circulation recovery driven by composition-dependent density changes. The NACOK graphite oxidation experiment validates the oxidation model: SAM predicts a bottom-level graphite weight loss of 25%, in close agreement with the measured 24%, and reproduces the strong axial nonuniformity and block-geometry dependence of oxidation, at a level comparable to the SPECTRA and TINTE codes. The validated capabilities are then exercised together in an integrated, reactor-scale simulation of a DLOFC air-ingress transient in a simplified HTR-PM pebble-bed reactor. In a single calculation spanning approximately 8 days, SAM reproduces the complete accident sequence: rapid depressurization, densityand diffusion-driven air ingress over ˜15 hours, onset of buoyancy-driven natural circulation, exothermic graphite oxidation with a peak fuel temperature at ˜62 hours, and eventual passive cooldown. These results demonstrate that SAM now provides the nuclear community with a preliminarily validated, modern systemlevel tool for HTGR air-ingress safety analysis, filling a recognized capability gap. Future extensions to broaden species tracking, improve oxidation chemistry, and refine the reactor model are discussed.

Yang, Gang↗
DOE OSTI2022

A self-consistent hybrid model of kinetic striations in low-current argon discharges

A self-consistent hybrid model of standing and moving striations was developed for low-current DC discharges in noble gases. We introduced the concept of surface diffusion in phase space ( r , u ) (where u denotes the electron kinetic energy) described by a tensor diffusion in the nonlocal Fokker–Planck kinetic equation for electrons in the collisional plasma. Electrons diffuse along surfaces of constant total energy ε = u - eφ ( r ) between energy jumps in inelastic collisions with atoms. Numerical solutions of the 1d1u kinetic equation for electrons were obtained by two methods and coupled to ion transport and Poisson solver. We studied the dynamics of striation formation in Townsend and glow discharges in argon gas at low discharge currents using a two-level excitation-ionization model and a ‘full-chemistry’ model, which includes stepwise and Penning ionization. Standing striations appeared in Townsend and glow discharges at low currents, and moving striations were obtained for the discharge currents exceeding a critical value. These waves originate at the anode and propagate towards the cathode. We have seen two types of moving striations with the two-level and full-chemistry models, which resemble the s and p striations previously observed in the experiments. Simulations indicate that processes in the anode region could control moving striations in the positive column plasma. The developed model helps clarify the nature of standing and moving striations in DC discharges of noble gases at low discharge currents and low gas pressures.

Physics↗
DOE OSTI2026

Collapse of Jahn-Teller phonons in La 1−x Sr x MnO 3 with weak magnetoresistance

Perovskite manganites are quantum materials exhibiting competing interactions inducing colossal magnetoresistance (CMR). The prevailing theory of CMR highlights the essential role of electron-phonon coupling (EPC), but mounting evidence suggests the underlying mechanism is more complicated. Here, we investigate phonons and spin-phonon coupling in ferromagnetic CMR manganites La 1−x Sr x MnO 3 (x=0.2,0.3) with relatively small CMR associated with melting of the magnetic order above room temperature. High-resolution neutron scattering experiments combined with density functional theory (DFT) show that the low-temperature ferromagnetic phase is conventional: neutron scattering from phonons agrees with DFT predictions and magnons follow sinusoidal dispersions. Fluctuating magnetic moments and low-energy phonons remain conventional in the high-temperature paramagnetic phase, indicating the Mn and La/Sr sublattices are not strongly perturbed by melting of ferromagnetism. In contrast, the Jahn–Teller-active optical oxygen vibrations collapse entirely above the Curie temperature, despite low CMR in these compositions, with some of the lost spectral weight reappearing as quasielastic scattering. We attribute this highly anomalous behavior to giant EPC in the charge and/or orbital channel. It drives cooperative diffusive motion of quasistatic carrier-trapping oxygen sublattice distortions once ferromagnetism disappears. We hypothesize the magnitude of magnetoresistance correlates with the rate of diffusion rather than with the strength of Jahn–Teller EPC.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Particle charge distributions in the effluent of a flow-through atmospheric pressure low temperature plasma

Atmospheric pressure low-temperature plasmas are often utilized to perform particle synthesis, treatment, and removal. It is well-known that dust particles are highly negatively charged in these plasmas; however, little is known about dust particle charging behavior as particles leave the plasma volume and pass through the spatial afterglow region. In this work, monodisperse particles of various sizes and work functions were introduced into an atmospheric pressure radiofrequency capacitively coupled flow-through plasma. Dust particle electrical mobility distributions downstream of the flow-through plasma were measured utilizing a differential mobility analyzer in conjunction with a condensation particle counter at various gas flow velocities. Charge distributions were determined from the measured electrical mobility distributions. Experiments confirm that particles become less negatively charged, and even net-positively charged after leaving the plasma volume, with a distribution that follows a shifted Boltzmann charge distribution. Additionally, particle charge in the effluent of the flow-through plasma is negligibly dependent on work function but highly size and flow velocity dependent. Larger particles were shown to have a higher magnitude of charge under all studied conditions; however, particle polarity was switchable by varying gas flow velocity. The charging dynamics were simulated utilizing a constant number Monte Carlo model that accounts for electron temperature decay and the transition from ambipolar to free diffusion of electrons and ions in the spatial afterglow. Simulation results also suggest that, at the same flow velocity, larger particles obtain a greater magnitude of charge, negative or positive. Furthermore, the decrease in electron mobility and the difference between ion and electron convective loss rates create an ion-rich region in the plasma effluent that promotes ion–particle collisions and drives particle charge removal and even reversal of polarity. Larger particles more favorably collide with energetic species in these environments, which results in higher charge states.

atmospheric pressure low temperature plasma↗
DOE OSTI2024

Structure and Dynamics of Aqueous Electrolytes at Quartz (001) and (101) Surfaces

Here, understanding and describing reactivity at mineral-water interfaces such as ion adsorption, the kinetics of dissolution, or surface charge development depends on our ability to improve the accuracy of electrical double layer (EDL) models. While molecular dynamics (MD) simulations are routinely used to investigate the structure and energetics of adsorbed ions comprising the EDL, less attention is paid to their self-diffusion dynamics, which can uniquely inform on coupling to interfacial reactions. Here we use MD to investigate both the organization and diffusion dynamics of water and electrolyte ions (NaCl, KCl, CaCl 2 ) at hydroxylated quartz (001) and (101) surfaces, a comparison which allowed us to assess surface structural effects of corrugation and silanol density. We found that inner- versus outer-sphere complex formation depends on cation size and charge but not necessarily hydration energies. Participation of surface silanols in the hydration spheres of Na + and K + generally indicated their preference for inner-sphere complexation, but this depends strongly on the orientation of the surface considered through its influence over the organization and dynamics of adsorbed water layers. In particular, surface orientation substantially affects the diffusive behavior of the near-surface water. Na + was found to decrease the mobility of water in the first layer, consistent with an increasing frequency of hydrolysis implied by faster quartz dissolution rates observed in experiments via the well known salt effect. Our results are also in good agreement with the observed dissolution rate of quartz vs. surface adsorption strength measure by Dove and Nix. This study sets the stage for a forthcoming paper examining how the dynamics at quartz/electrolyte interfaces are influenced by externally applied electric fields.

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

Pervasive symmetry-lowering nanoscale structural fluctuations in the cuprate La 2−𝑥 ⁢Sr 𝑥 ⁢Cu⁢O 4

The cuprate superconductors are among the most widely studied quantum materials, yet there remain fundamental open questions regarding their electronic properties and the role of structural degrees of freedom. Recent neutron and x-ray scattering measurements uncovered exponential scaling with temperature of the strength of orthorhombic fluctuations in the tetragonal phase of La 2−𝑥⁢ Sr 𝑥 ⁢Cu⁢O 4 and Tl 2 ⁢Ba 2 ⁢Cu⁢O 6+𝛿 , unusual behavior that closely resembles prior results for the emergence of superconducting fluctuations, and that points to a common origin rooted in inherent correlated structural inhomogeneity. Here, in this study, we use neutron and x-ray diffuse scattering to provide further insight into structural fluctuations in La 2−𝑥 ⁢Sr 𝑥 ⁢Cu⁢O 4 . We perform measurements up to temperatures that approach the melting point of the undoped parent compound, and we investigate the effects of in situ in-plane uniaxial elastic stress and plastic deformation. Our neutron scattering results for the parent compound La 2 ⁢Cu⁢O 4 reveal that short-range orthorhombic fluctuations persist to the maximum experimental temperature of nearly 1000 K, i.e., to a significant fraction of the crystallization temperature. At this temperature, the spatial characteristic length extracted from the momentum-space data is still about three lattice constants. The neutron scattering experiment enables quasistatic discrimination and reveals that the response is increasingly dynamic at higher temperatures. We also find that neither compressive elastic stress along the tetragonal [110] direction nor irreversible plastic deformation activating the ⟨100⟩⁢[001] slip system significantly alter the robust orthorhombic fluctuations. Overall, these results support the notion that these fluctuations couple to subtle, underlying inhomogeneity that underpins the cuprate phase diagram. Finally, in La 1.8⁢ Sr 0.2⁢ Cu⁢O 4 , we uncover low-energy structural fluctuations at a nominally forbidden reflection that are distinct from the orthorhombic distortions. While the origin of these fluctuations is not clear, they might be related to the presence of extended defects such as dislocations or stacking faults.

Spieker, R. J. [Univ. of Minnesota, Minneapolis, M↗
DOE OSTI2022

Nuclear Theory from First Principles to Forefront Experiments (Final Report for DE-SC0018638)

The Lee research group is a part of the Nuclear Lattice Effective Field Theory Collaboration, which has developed and performed ab initio lattice simulations of nuclear structure, scattering, reactions, and many-body systems. The other senior members of the collaboration include Ulf-G. Meißner at Bonn/Julich, Evgeny Epelbaum and Hermann Krebs at Bochum, Timo Lahde and Thomas Luu at Julich, and Gautam Rupak at Mississippi State. Our letter “Ab initio alpha-alpha scattering” was featured in a Nature News and Views article. Another letter “Nuclear binding near a quantum phase transition” was highlighted in a Viewpoint article in the online APS journal Physics as well as a news article in the IOP publication Physics World (September 21, 2016). Our letter “Eigenvector continuation with subspace learning” was also highlighted a Synopsis article in Physics. There are many promising ab initio approaches being used to calculate the properties of few-and many-nucleon systems. This includes the no-core shell model, symmetry-adapted no-core shell model quantum Monte Carlo, auxiliary-field diffusion Monte Carlo, self-consistent Green’s functions, many-body perturbation theory, in-medium similarity renormalization group, and coupled cluster methods. Lattice effective field theory is another ab initio approach which combines the framework of effective field theory with lattice Monte Carlo methods to allow favorable scaling from few- to many-body systems. Perhaps the most important aspect of lattice effective field theory is that its strengths and weaknesses are orthogonal to that of other ab initio methods. For example, lattice effective field theory has little difficulty in probing cluster structures and collectivity, while such features are much more difficult using other methods. Furthermore it can be used to compute superfluid condensate fractions as well as the phase diagram of strongly matter and the density and temperature dependence of clustering. Lattice effective field theory was first used in simulations of infinite nuclear matter and infinite neutron matter with pions and without pions. In addition to the efforts by our collaboration, there have been recent efforts by other groups as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Numerical thermalization in 2D PIC simulations: Practical estimates for low-temperature plasma simulations

The process of numerical thermalization in particle-in-cell (PIC) simulations has been studied extensively. It is analogous to Coulomb collisions in real plasmas, causing particle velocity distributions (VDFs) to evolve toward a Maxwellian as macroparticles experience polarization drag and resonantly interact with the fluctuation spectrum. This paper presents a practical tutorial on the effects of numerical thermalization in 2D PIC applications. Scenarios of interest include simulations, which must be run for many thousands of plasma periods and contain a population of cold electrons that leave the simulation space very slowly. This is particularly relevant to many low-temperature plasma discharges and materials processing applications. We present numerical drag and diffusion coefficients and their associated timescales for a variety of grid resolutions, discussing the circumstances under which the electron VDF is modified by numerical thermalization. Though the effects described here have been known for many decades, direct comparison of analytically derived, velocity-dependent numerical relaxation timescales to those of other relevant processes has not often been applied in practice due to complications that arise in calculating thermalization rates in 1D simulations. Using these comparisons, we estimate the impact of numerical thermalization in several examples of low-temperature plasma applications including capacitively coupled plasma discharges, inductively coupled plasma discharges, beam plasmas, and hollow cathode discharges. Finally, we discuss possible strategies for mitigating numerical relaxation effects in 2D PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
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↗
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↗