Search NASA⌕ Search

SEARCH · Search NASA

Results for “matrix diffusion”

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 451 records · Page 25

Conductance and resonant tunneling in multi-channel double barrier structures under transverse and longitudinal electric fields

Transport properties of electrons through biased double barrier semiconductor structures with finite transverse width w{sub y}, in the presence of a channel-mixing transverse electric field E{sub T} (along the y-axis), were studied. We solve the multichannel Schrödinger equation using the transfer matrix method and transport properties, like the conductance G and the transmission coefficients T{sub ij} have been evaluated as functions of the electrons' energy E and the transverse and longitudinal (bias) electric forces, f{sub T} and f{sub b}. We show that peak-suppression effects appear, due to the applied bias. Similarly, coherent interference of wave-guide states induced by the transverse field is obtained. We show also that the coherent interference of resonant wave-guide states gives rise to resonant conductance, which can be tuned to produce broad resonant peaks, implying operation frequencies of the order of 10 THz or larger.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Extension of modified power method to two-dimensional problems

In this study, the generalized modified power method was extended to two-dimensional problems. A direct application of the method to two-dimensional problems was shown to be unstable when the number of requested eigenmodes is larger than a certain problem dependent number. The root cause of this instability has been identified as the degeneracy of the transfer matrix. In order to resolve this instability, the number of sub-regions for the transfer matrix was increased to be larger than the number of requested eigenmodes; and a new transfer matrix was introduced accordingly which can be calculated by the least square method. The stability of the new method has been successfully demonstrated with a neutron diffusion eigenvalue problem and the 2D C5G7 benchmark problem. - Graphical abstract:.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nitrogen effects on radiation response in 12Cr ferritic/martensitic alloys

Ferritic/martensitic steels are one of the best candidates for structural materials for high dose applications in next generation nuclear reactors. The composition of structural materials must be optimized for reliable service during irradiation. This study reports the effect of interstitial nitrogen on radiation response in 12Cr ferritic/martensitic HT9 steels having a controlled nitrogen concentration. Results show that a high amount of ‘free’ nitrogen in the matrix stabilizes the interstitial clusters which leads to (i) larger loop sizes (ii) lower loop density and (iii) slightly reduced radiation induced hardening. It also affects diffusion mechanism of Ni and formation of Ni/Si-rich precipitates.

36 MATERIALS SCIENCE↗

Water sub-diffusion in membranes for fuel cells

We investigate the dynamics of water confined in soft ionic nano-assemblies, an issue critical for a general understanding of the multi-scale structure-function interplay in advanced materials. We focus in particular on hydrated perfluoro-sulfonic acid compounds employed as electrolytes in fuel cells. These materials form phase-separated morphologies that show outstanding proton-conducting properties, directly related to the state and dynamics of the absorbed water. We have quantified water motion and ion transport by combining Quasi Elastic Neutron Scattering, Pulsed Field Gradient Nuclear Magnetic Resonance, and Molecular Dynamics computer simulation. Effective water and ion diffusion coefficients have been determined together with their variation upon hydration at the relevant atomic, nanoscopic and macroscopic scales, providing a complete picture of transport. We demonstrate that confinement at the nanoscale and direct interaction with the charged interfaces produce anomalous sub-diffusion, due to a heterogeneous space-dependent dynamics within the ionic nanochannels. This is irrespective of the details of the chemistry of the hydrophobic confining matrix, confirming the statistical significance of our conclusions. Our findings turn out to indicate interesting connections and possibilities of crossfertilization with other domains, including biophysics. They also establish fruitful correspondences with advanced topics in statistical mechanics, resulting in new possibilities for the analysis of Neutron scattering data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Computational modeling of charge hopping dynamics along a disordered one-dimensional wire with energy gradients in quantum environments

This computational study investigates the effects of energy gradients on charge hopping dynamics along a one-dimensional chain of discrete sites coupled to quantum bath, which is modeled at the level of Pauli master equation (PME). This study also assesses the performance of different approximations for the hopping rates. Three different methods for solving the PME, a fourth order Runge–Kutta method, numerical diagonalization of the rate matrix followed by analytic propagation, and kinetic Monte Carlo simulation method, are tested and confirmed to produce virtually identical values of time dependent mean square displacement, diffusion constant, and mobility. Five different rate expressions, exact numerical evaluation of Fermi’s Golden Rule (FGR) rate, stationary phase interpolation (SPI) approximation, semiclassical approximation, classical Marcus rate, and Miller–Abrahams rate, are tested to help understand the effects of approximations in representing quantum environments in the presence of energy gradients. The results based on direct numerical evaluation of FGR rate exhibit transition from diffusive to non-diffusive behavior with the increase in the gradient and show that the charge transport in the quantum bath is more sensitive to the magnitude of the gradient and the disorder than in the classical bath. Among all the four approximations for the hopping rates, the SPI approximation is confirmed to work best overall. A comparison of two different methods to calculate the mobility identifies drift motion of the population distribution as the major source of non-diffusive behavior and provides more reliable information on the contribution of quantum bath.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assembly of Rabbit Capsules for Irradiation of Pyrolytic Carbon / Silicon Carbide Diffusion Couples in the High Flux Isotope Reactor

Tristructural isotropic (TRISO)–coated particle fuel is a promising advanced fuel concept being considered for several advanced reactor applications and for accident-tolerant fuel for light water reactors. One of the aspects studied in the development of this advanced fuel concept is the release of specific fission products (Ag, Eu, and Sr). The silicon carbide (SiC) layer of TRISO fuel serves as the primary barrier to metallic fission products and actinides not retained in the fuel kernel. The goal of this project is to evaluate the effect of irradiation on the diffusion of these fission products in the SiC layer of the fuel. For this purpose, rabbit capsules containing small slab diffusion couple specimens have been assembled to be irradiated in the High Flux Isotope Reactor (HFIR). The diffusion couple specimens have been fabricated using similar processes and equipment as those used to make TRISO particles; the desired fission products have been implanted in the specimens using an ion accelerator. Moreover, the effect of temperature on the fission products diffusion will be studied separately by performing thermal experiments in the absence of irradiation. This report describes the irradiation experiment design concept, summarizes the irradiation test matrix, and reports on the successful assembly of two rabbit capsules that will be irradiated in the HFIR.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Measurement of Fission Product Concentration Profiles in AGR-3/4 TRISO Fuel Graphitic Matrix and Nuclear Graphites

The third Advanced Gas Reactor (AGR) irradiation experiment, AGR-3/4, was designed to investigate the migration of fission products in fuel compact graphitic matrix and reactor graphite components. Using destructive methods, radial fission product concentration profiles were measured for gamma-emitting fission products (e.g., Ag-110m, Cs-134, and Eu-154) and beta-emitting Sr-90 in irradiated graphitic and graphite components from six different AGR-3/4 irradiation capsules. These new measured concentration profiles can now be compared to non-destructive measurements and fission product transport simulations and will be used to derive new diffusivities and sorptivities to support refinement of fission product transport models and high-temperature gas-cooled reactor (HTGR) source-term analyses. Each capsule in the AGR-3/4 experiment had four fuel compacts in the middle of two concentric rings of graphitic matrix material, PCEA graphite, or IG-110 graphite. In addition to the approximately 1898 tristructural isotropic (TRISO) coated particles in each compact, there were 20 designed to fail (DTF) particles coated only in pyrocarbon so that they released fission products into the surrounding cylindrical rings of carbonaceous materials. Destructive sampling of the rings involved machining/milling material from around the circumference of the rings, collecting that material, and performing radiochemical analyses on it. Milling operations were performed in multiple steps or segments, and each segment was generally 0.508 mm (0.020 in) thick. Knowing the radial position at which each segment was milled, the volume of the milled material at each segment, and the fission product content in each segment, the radial fission product concentrations were constructed for select isotopes in each ring. Ag-110m profiles had the most variation. Some profiles were peaked at an inner or outer surface. Some were peaked at the middle of the ring wall thickness. Some increased radially outward, and some decreased radially outward. These types of variations and the fact that the measured profiles do not generally compare favorably with the transport model employed for AGR-3/4 may adversely impact the ability to extract reasonable transport parameters for this isotope. In many cases, the Cs-134 profiles decreased somewhat linearly in the outward radial direction, and in cursory comparisons, the shapes of these profiles appeared similar to those from model predictions. The step changes in concentration across the inner-outer ring gap were generally consistent with the model predictions as well. In some cases, there were local maxima in concentration at the outer surface of the rings. This suggests that fission products could have transported in the small gaps between the inner ring and the outer ring and between the outer ring and the sink ring such that some portion of a given fission product can bypass diffusion through the ring itself. The analysis of the small nubs on the outer surfaces of some of the outer rings revealed fission product concentrations in the nubs that were often higher than in the outermost segments of the rings. This further supports the hypothesis that short-circuit, gap transport occurred, causing relatively high surface concentrations on the outer surfaces of the rings. Eu-154 and Sr-90 profiles tend to have very similar shapes, suggesting that they transport via the same mechanisms. The observed profiles were indicative of a transport process where the isotopes are sorbed on the inner surface of the ring, but diffusion into the ring from that surface is quite slow. Some elevated concentrations of Sr-90 (relative to Eu-154) on the outer surface of a ring suggested that rapid, gas-gap transport of gaseous precursor Kr-90 and volatile Rb-90 could have occurred prior to their decaying to Sr-90. Overall, the Eu 154 and Sr-90 profiles were still very similar, which indicates that the transport of short-lived Sr-90 precursors is not a major effect. In some capsules, the qualitative Sr-90 behavior across the ring gaps was consistent with the model (using the available legacy Sr-90 transport parameters), but in other capsules the model was inconsistent with the measurements and seems to underestimate the amount of Sr-90 in the outer rings. The total ring Sr-90 inventories were estimated for all the rings that were subject to physical sampling. These results will be used to adjust the predicted particle and/or compact releases used in the AGR-3/4 fission product transport model. Given the different irradiation temperatures among the capsules and the rings, it was not possible to discern fundamental differences in the transport of isotopes within the different carbon materials, i.e., graphitic matrix, IG-110, or PCEA. It may be possible to do this in the course of determining transport from the concentration profiles in future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Elucidating the role of Cr migration in Ni-Cr exposed to molten FLiNaK via multiscale characterization

Structural materials used in nuclear reactor environments are subjected to coupled extremes such as high temperature, irradiation, and corrosion which act in concert to degrade their functional performance. Connecting alloy microstructure such as grain boundaries, and accumulating point defects with corrosion attack and pore morphology is critical to understanding underlying mechanisms. We uncover the compositional variations and morphology at multiple length scales in corrosion-damaged Ni-Cr alloy after exposure to oxidants in molten fluoride salts. A complex network of dense corrosion pores is detected by surface-level SEM observations. The corrosion pores take on a 1-dimensional morphology and are enriched with Ni and depleted of Cr 1–5 µm from the pore surface. STEM-EDS and 4D-STEM strain maps acquired simultaneously highlight the local variations in composition and structure of a ≤ 200 nm Cr-rich layer identified from a cross-section taken at the bottom of an isolated corrosion pore between the Ni-Cr alloy matrix and the embedded salt. However, the absence of an observed interface between the Ni-Cr alloy matrix and the FCC-structured Cr-rich layer suggests that Cr plating from the salt did not transpire. These findings support a proposed Cr lattice diffusion mechanism rather than Cr-precipitation from the salt to accommodate temperature transient conditions during sample cooling. Through the development of a 1D phase field model, these results are rationalized by formation energies for the Ni- and Cr-oxidation into the molten salt. This study reveals the locally altered microstructure caused by high temperature corrosion in non-steady-state molten salt nuclear reactor environments.

36 MATERIALS SCIENCE↗

Estimation and modeling of coal pore accessibility using small angle neutron scattering

Gas diffusion in coal is controlled by nano-structure of the pores. The interconnectivity of pores not only determines the dynamics of gas transport in the coal matrix but also influences the mechanical strength. In this study, small angle neutron scattering (SANS) was employed to quantify pore accessibility for two coal samples, one of sub-bituminous rank and the other of anthracite rank. Moreover, a theoretical pore accessibility model was proposed based on scattering intensities under both vacuum and zero average contrast (ZAC) conditions. Our results show that scattering intensity decreases with increasing gas pressure using deuterated methane (CD 4 ) at low Q values for both coals. Pores smaller than 40 nm in radius are less accessible for anthracite than sub-bituminous coal. On the contrary, when the pore radius is larger than 40 nm, the pore accessibility of anthracite becomes larger than that of sub-bituminous coal. Only 20% of pores are accessible to CD 4 for anthracite and 37% for sub-bituminous coal, where the pore radius is 16 nm. For these two coals, pore accessibility and pore radius follows a power-law relationship.

01 COAL, LIGNITE, AND PEAT↗

Multiscale Reactive Model for 1,3,5-Triamino-2,4,6-trinitrobenzene Inferred by Reactive MD Simulations and Unsupervised Learning

When high-energy-density materials are subjected to thermal or mechanical insults at extreme conditions (shock loading), a coupled response between the thermo-mechanical and chemical behaviors is systematically induced. Herein we develop a reaction model for the fast chemistry of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) at the mesoscopic scale, where the chemical behavior is determined by underlying microscopic reactive simulations. The slow carbon cluster formation is not discussed in the present work. All-atom reactive molecular dynamics (MD) simulations are performed with the ReaxFF potential, and a reduced-order chemical kinetics model for TATB is fitted to isothermal and adiabatic simulations of single crystal chemical decomposition. Unsupervised machine learning techniques based on non-negative matrix factorization are applied to MD trajectories to model the decomposition kinetics of TATB in terms of a four-component model. The associated heats of reaction are fit to the temperature evolution from adiabatic decomposition trajectories. Using a chemical species analysis, we show that non-negative matrix factorization captures the main chemical decomposition steps of TATB and provides an accurate estimation of their evolution with temperature. The final analytical formulation, coupled to a diffusion term, is incorporated into a continuum formalism, and simulation results are compared one-to-one against MD simulations of 1D reaction propagation along different crystallographic directions and with different initial temperatures. A good agreement is found for both the temporal and spatial evolution of the temperature field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Delayed Slater determinant update algorithms for high efficiency quantum Monte Carlo

Within ab initio Quantum Monte Carlo simulations, the leading numerical cost for large systems is the computation of the values of the Slater determinants in the trial wavefunction. Each Monte Carlo step requires finding the determinant of a dense matrix. This is most commonly iteratively evaluated using a rank-1 Sherman-Morrison updating scheme to avoid repeated explicit calculation of the inverse. The overall computational cost is therefore formally cubic in the number of electrons or matrix size. To improve the numerical efficiency of this procedure, we propose a novel multiple rank delayed update scheme. This strategy enables probability evaluation with application of accepted moves to the matrices delayed until after a predetermined number of moves, K . The accepted events are then applied to the matrices en bloc with enhanced arithmetic intensity and computational efficiency via matrix-matrix operations instead of matrix-vector operations. Here this procedure does not change the underlying Monte Carlo sampling or its statistical efficiency. For calculations on large systems and algorithms such as diffusion Monte Carlo where the acceptance ratio is high, order of magnitude improvements in the update time can be obtained on both multi- core CPUs and GPUs.

97 MATHEMATICS AND COMPUTING↗

A Monte Carlo synthetic-acceleration method for solving the thermal radiation diffusion equation

We present a novel synthetic-acceleration-based Monte Carlo method for solving the equilibrium thermal radiation diffusion equation in three spatial dimensions. The algorithm performance is compared against traditional solution techniques using a Marshak benchmark problem and a more complex multiple material problem. Our results show that our Monte Carlo method is an effective solver for sparse matrix systems. For solutions converged to the same tolerance, it performs competitively with deterministic methods including preconditioned conjugate gradient and GMRES. We also discuss various aspects of preconditioning the method and its general applicability to broader classes of problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

High Temperature Fuel Cladding Chemical Interactions Between TRIGA Fuels and 304 Stainless Steel

High-temperature fuel-cladding chemical interactions (FCCI) between TRIGA (Training, Research, Isotopes, General Atomics) fuel elements and the 304 stainless steel (304SS) are of interest to develop an understanding of the fuel behavior during transient reactor scenarios. TRIGA fuels are composed of uranium (U) particles dispersed in a zirconium-hydride (Zr-H) matrix. In reactor, the fuel is encased in 304-stainless-steel (304SS) or Incoloy 800 clad tubes. At high temperatures, the fuel can readily interact with the cladding, resulting in FCCI. A number of FCCI can take place in this system. Interactions can be expected between the cladding and the Zr-H matrix, and/or between the cladding and the U-particles. Other interactions may be expected between the Zr-H matrix and the U-particles. Furthermore, the fuel contains erbium-oxide (Er-O) additions. Interactions can also be expected between the Er-O, the cladding, the Zr-H and the U-particles. The overall result is that very complex interactions may take place as a result of fuel and cladding exposures to high temperatures. This report discusses the characterization of the baseline fuel microstructure in the as-received state (prior to exposure to high temperature), characterization of the fuel after annealing at 950C for 24 hours and the results from diffusion couple experiments carries out at 1000C for 5 and 24 hours. Characterization was carried out via scanning electron microscopy (SEM) and transmission electron microscopy (TEM) with sample preparation via focused ion beam in situ-liftout-technique.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Diffusion study of uranium mononitride/zirconium carbide composite for space nuclear propulsion

Here, the next generation of space exploration will require extensive developments in rocket technology. Space nuclear propulsion is of interest due to its high fuel density and power, however, it also has high temperature and stability demands. This study examines a uranium nitride/zirconium carbide ceramic-ceramic particle composite as a potential fuel for these missions. Based on the constituent properties, this fuel composite is expected to be highly thermally efficient and resistant to the hot hydrogen propellant. One of the main concerns with this composite is unknown diffusion behavior between UN and ZrC over time. Diffusion couples and composites were used to determine how these constituents will behave when left in contact at high temperatures. UN and ZrC were found to have limited but observable diffusion at the phase boundary. The resulting phase is a UZr(CN) quaternary phase. UZr(CN) has been examined for other high temperature and gas nuclear reactors with preliminary success; however, there lacks sufficient data to fully qualify this phase. Based on the limited information available, the resultant quaternary phase could be deemed acceptable and even provide better adhesion for the fuel particles (UN) to the matrix (ZrC).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

New theory of diffusive and coherent nature of optical wave via a quantum walk

We propose a new theory on a relation between diffusive and coherent nature in one dimensional wave mechanics based on a quantum walk. It is known that the quantum walk in homogeneous matrices provides the coherent property of wave mechanics. Using the recent result of a localization phenomenon in a one-dimensional quantum walk (Konno, 2010), we numerically show that the randomized localized matrices suppress the coherence and give diffusive nature. - Highlights: • A relation between diffusive and coherent nature of wave is obtained. • A rigorous decoherent process is observed in quantum walk. • Quantum walk for a point-process is numerically considered. • Transfer matrix theory with S-matrix in optics is connected with quantum walk.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

UO 2 + 5 vol% ZrB 2 nano composite nuclear fuels with full boron retention and enhanced oxidation resistance

The boron isotope ( 10 B) can be used as a neutron absorber in UO 2 to control the reactivity of nuclear fuel pellets, however, the boron source can react with oxygen source in UO 2 to form B 2 O 3 that vaporize readily at temperatures above 1200 °C. Unfortunately, the sintering of UO 2 fuel requires hours holding at high temperature (>1700 °C), resulting in an inevitable B loss during sintering and unpredictable B concentration in final product. It is challenging to incorporate boron through a conventional sintering method. In this work, we demonstrated that spark plasma sintering (SPS), a field assisted sintering technology, can effectively densify UO 2 + 5 vol% ZrB 2 composite fuel pellets by rapid consolidation at 1600 °C for a short duration of 5 min under an applied pressure of 40 MPa. Further, thermogravimetric analysis (TGA) measurements confirm that ZrB 2 is fully retained inside the composite fuel pellets. Inside the composite fuel pellets, nano sized ZrB 2 particles are uniformly distributed along the grain boundaries of the UO 2 matrix. The ZrB 2 particle transforms to a glassy B 2 O 3 phase covering the sample surface and grain boundaries of UO 2 matrix after a simple post-sintering annealing at 1000 °C in flowing Argon gas for 4 h. The formed glassy B 2 O 3 slows down the diffusion of oxygen ions and postpones the onset temperature for oxidation of UO 2 from 400 °C to 550 °C. This study demonstrates the capability of SPS, an advanced fuel manufacturing technique, to achieve a full retention of ZrB 2 in UO 2 oxide fuel and increase oxidation resistance through a simple post-sintering annealing. The reported work holds great engineering potential for development of advanced oxide fuel for nuclear application.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A transmission electron microscopy and X-ray photoelectron spectroscopy study of annealing induced γ-phase nucleation, clustering, and interfacial dynamics in reactively sputtered amorphous alumina thin films

Pure α-Al{sub 2}O{sub 3} exhibits a very high degree of thermodynamical stability among all metal oxides and forms an inert oxide scale in a range of structural alloys at high temperatures. We report that amorphous Al{sub 2}O{sub 3} thin films sputter deposited over crystalline Si instead show a surprisingly active interface. On annealing, crystallization begins with nuclei of a phase closely resembling γ-Alumina forming almost randomly in an amorphous matrix, and with increasing frequency near the substrate/film interface. This nucleation is marked by the signature appearance of sharp (400) and (440) reflections and the formation of a diffuse diffraction halo with an outer maximal radius of ≈0.23 nm enveloping the direct beam. The microstructure then evolves by a cluster-coalescence growth mechanism suggestive of swift nucleation and sluggish diffusional kinetics, while locally the Al ions redistribute slowly from chemisorbed and tetrahedral sites to higher anion coordinated sites. Chemical state plots constructed from XPS data and simple calculations of the diffraction patterns from hypothetically distorted lattices suggest that the true origins of the diffuse diffraction halo are probably related to a complex change in the electronic structure spurred by the a-γ transformation rather than pure structural disorder. Concurrent to crystallization within the film, a substantially thick interfacial reaction zone also builds up at the film/substrate interface with the excess Al acting as a cationic source.

36 MATERIALS SCIENCE↗

Final Report - Subcontract B6183769

During my visit to LLNL during July 5{15, 2016, we worked on linear system solvers. The two level hierarchical solver that initiated our study was developed to solve linear system arising from hp adaptive finite element calculations, and is implemented in the PLTMG software package, version 12 [1]. This preconditioner typically requires 3-20% of the space used by the stiffness matrix for higher order elements. It has multigrid like convergence rates for a wide variety of PDEs (self-adjoint positive definite elliptic equations, convection dominated convection-diffusion equations, and highly inde nite Helmholtz equations, among others). The convergence rate is not independent of the polynomial degree p as p ! 1, but but remains strong for p 9, which is the highest polynomial degree allowed in PLTMG, due to limitations of numerical quadrature formulae. A more complete description of the method and some numerical experiments illustrating its effectiveness appear in [2]. Like traditional geometric multilevel methods, this scheme relies on extensive knowledge of the underlying finite element space in order to construct both the smoother and the coarse grid correction components.

97 MATHEMATICS AND COMPUTING↗