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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

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↗
DOE OSTI2015

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↗
DOE OSTI2023

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↗
DOE OSTI2017

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↗
NASA NTRS2002

Techniques Optimized for Reducing Instabilities in Advanced Nickel-Base Superalloys for Turbine Blades

The High-Speed Research (HSR) Airfoil Alloy program developed fourth-generation single-crystal superalloys with up to an 85 F increase in creep rupture capability over current production airfoil alloys. Recent results have been generated at the NASA Glenn Research Center on these fourth-generation alloys, but in coated form, for subsonic turbine blade applications under NASA's Ultra-Efficient Engine Technology (UEET) Program. One goal for UEET is to optimize the airfoil alloy/thermal barrier coating system for 3100 F turbine inlet temperatures. The state-of-the art turbine blade airfoil system consists of a superalloy single crystal that provides the basic mechanical performance of the airfoil. A thermal barrier coating is used to reduce the temperature of the base superalloy, and a bondcoat is deposited between the base material and the thermal barrier coating. The bondcoat improves the oxidation and corrosion resistance of the base superalloy and improves the spallation resistance of the thermal barrier coating. A commercial platinum aluminide bondcoat was applied to the HSR-developed alloys, and a diffusion zone developed as a result of interaction between the bondcoat and the superalloy. Optimized strength is obtained for superalloys when the refractory element content is high and the limits of microstructural stability are approached or exceeded slightly. For fourthgeneration alloys, instability leads to the formation of topologically close packed (TCP) phases, which form internally in the superalloy, and a secondary reaction zone (SRZ), which forms under the diffusion zone. There was a concern that excessive quantities of either TCP or SRZ might decrease the mechanical properties of the superalloy, with SRZ thought to be particularly detrimental and its formation unpredictable. Thus, an SRZreduction effort was initiated in the NASA UEET Program so that methods developed during the HSR project could be optimized further to reduce or eliminate the SRZ. An SRZ is a three-phase constituent composed of TCP and stringers of gamma phase in a matrix of gamma prime. An incoherent grain boundary separates the SRZ from the gammagamma prime microstructure of the superalloy. The SRZ is believed to form as a result of local chemistry changes in the superalloy due to the application of the diffusion aluminide bondcoat. Locally high surface stresses also appear to promote the formation of the SRZ. Thus, techniques that change the local alloy chemistry or reduce surface stresses have been examined for their effectiveness in reducing SRZ. These SRZ-reduction steps are performed on the test specimen or the turbine blade before the bondcoat is applied. Stressrelief heat treatments developed at NASA Glenn have been demonstrated to reduce significantly the amount of SRZ that develops during subsequent high-temperature exposures. Stress-relief heat treatments reduce surface stresses by recrystallizing a thin surface layer of the superalloy. However, in alloys with very high propensities to form SRZ, stress relief heat treatments alone do not eliminate SRZ entirely. Thus, techniques that modify the local chemistry under the bondcoat have been emphasized and optimized successfully at Glenn. One such technique is carburization, which changes the local chemistry by forming submicron carbides near the surface of the superalloy. Detailed characterizations have demonstrated that the depth and uniform distribution of these carbides are enhanced when a stress relief treatment and an appropriate surface preparation are employed in advance of the carburization treatment. Even in alloys that have the propensity to develop a continuous SRZ layer beneath the diffusion zone, the SRZ has been completely eliminated or reduced to low, manageable levels when this combination of techniques is utilized. Now that the techniques to mitigate SRZ have been established at Glenn, TCP phase formation is being emphasized in ongoing work under the UEET Program. The limitsf stability of the fourth-generation alloys with respect to TCP phase formation are currently being defined along with high-temperature creep rupture properties. In addition, a regression model is being developed at Glenn for the prediction of the presence of TCP phase in the microstructure and SRZ under the diffusion zone. The model is based on a design-of-experiments methodology with emphasis on the potential synergistic effects of alloying elements.

MacKay, Rebecca A.↗
NASA NTRS1986

New experimental and analytical results for diffusion and swelling of resins used in graphite/epoxy composite materials

The epoxy resins currently in use can slowly absorb moisture from the atmosphere over a long period. This reduces those mechanical properties of composites which depend strongly on the matrix, such as compressive strength and buckling instabilities. The effect becomes greater at elevated temperatures. The paper will discuss new phenomena which occur under simultaneous temperature and moisture variations. An analytical model will also be discussed and documented.

Hiel, C. C.↗
NASA NTRS1991

Thermal histories of CO3 chondrites - Application of olivine diffusion modelling to parent body metamorphism

The petrologic sequence observed in the CO3 chondrite group has been suggested to be the result of thermal metamorphism on a parent body. A model developed to examine the possibility that chondrule and matrix olivines equilibrated in situ, during parent body metamorphism is presented. The model considers Fe-Mg interdiffusion between chondrule and matrix olivines. Zoning profiles comparable to those observed in chondrule olivines from partially equilibrated members of the series are reproduced successfully. Metamorphism of CO3 chondrites on a parent body is therefore a viable model for the observed equilibration. Results indicate that peak metamorphic temperatures experienced by the CO3 chondrites were around 500 C, and that the range of peak temperatures between unequilibrated and equilibrated subtypes was relatively narrow, around 100 C.

Jones, Rhian H.↗
DOE OSTI2016

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↗
NASA NTRS2016

Scattering of Gaussian Beams by Disordered Particulate Media

A frequently observed characteristic of electromagnetic scattering by a disordered particulate medium is the absence of pronounced speckles in angular patterns of the scattered light. It is known that such diffuse speckle-free scattering patterns can be caused by averaging over randomly changing particle positions and/or over a finite spectral range. To get further insight into the possible physical causes of the absence of speckles, we use the numerically exact superposition T-matrix solver of the Maxwell equations and analyze the scattering of plane-wave and Gaussian beams by representative multi-sphere groups. We show that phase and amplitude variations across an incident Gaussian beam do not serve to extinguish the pronounced speckle pattern typical of plane-wave illumination of a fixed multi-particle group. Averaging over random particle positions and/or over a finite spectral range is still required to generate the classical diffuse speckle-free regime.

radiative transfer theory↗
DOE OSTI2023

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↗
DOE OSTI2017

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↗
DOE OSTI2020

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↗
DOE OSTI2017

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↗
DOE OSTI2018

Final Report - Subcontract B629114

During my visits to LLNL August 20-24, 2018, and September 10-14, 2018, I worked on linear system solvers. The two level hierarchical solver that initiated our study was developed to solve linear systems arising from $hp$ adaptive finite element calculations, and is implemented in the PLTMG software package, version 13. 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 indefinite Helmholtz equations, among others). The convergence rate is not independent of the polynomial degree $p$ as $p → ∞$, but but remains strong for $p ≤$ 9, which is the highest polynomial degree allowed in PLTMG , due to limitations of the numerical quadrature rules implemented in the software package. A more complete description of the method and some numerical experiments illustrating its effectiveness appear. Like traditional geometric multilevel methods, this scheme relies on knowledge of the underlying finite element space in order to construct the smoother and the coarse grid correction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Subcontract B635525 (Final Report)

During my visits to LLNL July 8–12 and July 22–26, 2019, I worked on linear system solvers. The two level hierarchical solver that initiated our study was developed to solve linear systems arising from $\textit{hp}$ adaptive finite element calculations, and is implemented in the $\textit{PLTMG}$ software package, version 13. 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 indefinite Helmholtz equations, among others). The convergence rate is not independent of the polynomial degree $\textit{p}$ as $\textit{p} → ∞$, but remains strong for $\textit{p}$ ≤ 9, which is the highest polynomial degree allowed in $\textit{PLTMG}$, due to limitations of the numerical quadrature rules implemented in the soft- ware package. A more complete description of the method and some numerical experiments illustrating its effectiveness appear in. Like traditional geometric multilevel methods, this scheme relies on knowledge of the underlying finite element space in order to construct the smoother and the coarse grid correction.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

A finite element method for the thermochemical decomposition of polymeric materials. I - Theory

The governing differential equations are developed to model the thermomechanical behavior of chemically decomposing, polymeric materials. These equations account for thermal and gaseous diffusion through a poroelastic, transversely isotropic solid. The Bubnov-Galerkin finite element method is applied to the governing equations to cast the coupled set into a single matrix equation. A method for solving these equations simultaneously at each time step is discussed.

Sullivan, R. M.↗
NASA NTRS2011

Evaluation of Lightning Induced Effects in a Graphite Composite Fairing Structure

Defining the electromagnetic environment inside a graphite composite fairing due to lightning is of interest to spacecraft developers. This paper is the first in a two part series and studies the shielding effectiveness of a graphite composite model fairing using derived equivalent properties. A frequency domain Method of Moments (MoM) model is developed and comparisons are made with shielding test results obtained using a vehicle-like composite fairing. The comparison results show that the analytical models can adequately predict the test results. Both measured and model data indicate that graphite composite fairings provide significant attenuation to magnetic fields as frequency increases. Diffusion effects are also discussed. Part 2 examines the time domain based effects through the development of a loop based induced field testing and a Transmission-Line-Matrix (TLM) model is developed in the time domain to study how the composite fairing affects lightning induced magnetic fields. Comparisons are made with shielding test results obtained using a vehicle-like composite fairing in the time domain. The comparison results show that the analytical models can adequately predict the test and industry results.

Trout, Dawn H.↗
DOE OSTI2018

A meshless Galerkin method for non-local diffusion using localized kernel bases

Here, we introduce a meshless method for solving both continuous and discrete variational formulations of a volume constrained, nonlocal diffusion problem. We use the discrete solution to approximate the continuous solution. Our method is nonconforming and uses a localized Lagrange basis that is constructed out of radial basis functions. By verifying that certain inf-sup conditions hold, we demonstrate that both the continuous and discrete problems are well-posed, and also present numerical and theoretical results for the convergence behavior of the method. The stiffness matrix is assembled by a special quadrature routine unique to the localized basis. Combining the quadrature method with the localized basis produces a well-conditioned, symmetric matrix. This then is used to find the discretized solution.

97 MATHEMATICS AND COMPUTING↗