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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Sulfur vacancies in photorefractive Sn{sub 2}P{sub 2}S{sub 6} crystals

A photoinduced electron paramagnetic resonance (EPR) spectrum in single crystals of Sn{sub 2}P{sub 2}S{sub 6} (SPS) is assigned to an electron trapped at a sulfur vacancy. These vacancies are unintentionally present in undoped SPS crystals and are expected to play an important role in the photorefractive behavior of the material. Nonparamagnetic sulfur vacancies are formed during the initial growth of the crystal. Subsequent illumination below 100 K with 442 nm laser light easily converts these vacancies to EPR-active defects. The resulting S = 1/2 spectrum shows well-resolved and nearly isotropic hyperfine interactions with two P ions and two Sn ions. Partially resolved interactions with four additional neighboring Sn ions are also observed. Principal values of the g matrix are 1.9700, 1.8946, and 1.9006, with the corresponding principal axes along the a, b, and c directions in the crystal. The isotropic parts of the two primary {sup 31}P hyperfine interactions are 19.5 and 32.6 MHz and the isotropic parts of the two primary Sn hyperfine interactions are 860 and 1320 MHz (the latter values are each an average for {sup 117}Sn and {sup 119}Sn). These hyperfine results suggest that singly ionized sulfur vacancies have a diffuse wave function in SPS crystals, and thus are shallow donors. Before illumination, sulfur vacancies are in the doubly ionized charge state because of compensation by unidentified acceptors. They then trap an electron during illumination. The EPR spectrum from the sulfur vacancy is destroyed when a crystal is heated above 120 K in the dark and reappears when the crystal is illuminated again at low temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Oxidation and Ignition Behaviors of Molten Mg-Nd-Zr Alloy in Resin-Sand Mold

High-temperature oxidation behavior of molten Mg-Nd-Zr in the pyrolysis atmosphere of resin-bind sand mold was investigated. The results showed that pouring temperature higher than 1023 K would sharply increase the ignition risk of melt in mold. The oxide formed on Mg-Nd-Zr melt below 973 K consists of two layers: mainly MgO for outer layer; mixture of MgO and Nd{sub 2}O{sub 3} for inner layer. The inner oxide mixture layer played an essential role in the protection of Mg-Nd-Zr matrix. Based on the experimental results and thermodynamic analysis, the oxidation process containing three stages was proposed. Nd- and Zn-enriched zone was found at the subsurface of oxide scale, without the enrichment of Zr. DICTRA calculation results indicated that this phenomenon was caused by the depletion of Mg and up-hill diffusion of Zr. Thus, Zr showed limited protection to the melt in long-term oxidation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Constraining Hamiltonians from chiral effective field theory with neutron-star data

Multi-messenger observations of neutron stars (NSs) and their mergers have placed strong constraints on the dense-matter equation of state (EOS). The EOS, in turn, depends on microscopic nuclear interactions that are described by nuclear Hamiltonians. These Hamiltonians are commonly derived within chiral effective field theory (EFT). Ideally, multi-messenger observations of NSs could be used to directly inform our understanding of EFT interactions, but such a direct inference necessitates millions of model evaluations. This is computationally prohibitive because each evaluation requires us to calculate the EOS from a Hamiltonian by solving the quantum many-body problem with methods such as auxiliary-field diffusion Monte Carlo (AFDMC), which provides very accurate and precise solutions but at a significant computational cost. Additionally, we need to solve the stellar structure equations for each EOS which further slows down each model evaluation by a few seconds. In this work, we combine emulators for AFDMC calculations of neutron matter, built using parametric matrix models, and for the stellar structure equations, built using multilayer perceptron neural networks, with the PyCBC data-analysis framework to enable a direct inference of coupling constants in an EFT Hamiltonian using multi-messenger observations of NSs. We find that astrophysical data can provide informative constraints on two-nucleon couplings despite the high densities probed in NS interiors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Biopolymer-supramolecular polymer hybrids for photocatalytic hydrogen production

Solar generation of H 2 is a promising strategy for dense energy storage. Supramolecular polymers composed of chromophore amphiphile monomers containing perylene monoimide (PMI) have been reported as crystalline light-harvesting assemblies for aqueous H 2 -evolving catalysts. Gelation of these supramolecular polymers with multivalent ions creates hydrogels with high diffusivity but insufficient mechanical stability and catalyst retention for reusability. We report here on using sodium alginate (SA) biopolymer to both induce supramolecular polymerization of PMI and co-immobilize them with catalysts in a robust hydrogel with high diffusivity that can also be 3D-printed. Faster mass transfer was achieved by controlling the material macrostructure by reducing gel diameter and microstructure by reducing biopolymer loading. Optimized gels produce H 2 at rates rivaling solution-based PMI and generate H 2 for up to 6 days. The PMI assemblies in the SA matrix create a percolation network capable of bulk-electron transfer under illumination. These PMI–SA materials were then 3D-printed on conductive substrates to create 3D hydrogel photoelectrodes with optimized porosity. The design of these versatile hybrid materials was bioinspired by the soft matter environment of natural photosynthetic systems and opens the opportunity to carry out light-to-fuel conversion within soft matter with arbitrary shapes and particular local environments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

MFiX: Fractional-Step Method Implementation

A comprehensive, multiphase computational fluid dynamics (CFD) simulation solves several coupled transport equations including continuity, momentum, species, and energy. Chemical reactions further couple these equations through heats of reaction and rates of formation of products and rates of destruction of reactants. A fractional-step method separates changes attributed to chemical reactions from transport phenomena like convection and diffusion. When the governing equations are split into the transport and reacting components, efficient and independent methodologies can be exploited to solve the different systems. Specifically, discretization of field variable transport equations results in large, sparse matrices which are loosely coupled. These systems are solved in succession using iterative techniques that take advantage of the matrix structure. In contrast, chemical reactions tightly couple field variables locally within the domain (e.g., within a single computational cell) resulting in low dimensional but dense, nonlinear systems that are better solved using direct integration techniques.

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

Facilitating Screening of MOFs for Mixed Matrix Membranes Using Machine Learning and the Maxwell Model

Metal organic framework (MOF)-based mixedmatrix membranes (MMMs), which embed MOF particles in polymer matrices, combine the advantages of polymeric and inorganic membranes. Multiple previous studies have used the Maxwell model together with molecular simulations and machine learning (ML) to predict the performance of MOF/polymer MMMs. However, the assumption of rigid MOF frameworks in molecular simulations limited the accuracy of the data used in the predictions, particularly in predicting molecular diffusivities. We developed a novel workflow integrating ML models with consideration of MOF flexibility to predict the permeability and selectivity of 131,722 MMMs for CO 2 /CH 4 , O 2 /N 2 and He/H 2 separations. The full range of achievable MMM performance within the Maxwell model was analyzed, and several promising MOFs were identified using this workflow. This approach offers an efficient tool for screening any polymer and MOF combination in gas separation applications.

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

Design-order, non-conformal low-Mach fluid algorithms using a hybrid CVFEM/DG approach

Here, a hybrid, design-order sliding mesh algorithm, which uses a control volume finite element method (CVFEM), in conjunction with a discontinuous Galerkin (DG) approach at non-conformal interfaces, is outlined in the context of a low-Mach fluid dynamics equation set. This novel hybrid DG approach is also demonstrated to be compatible with a classic edge-based vertex centered (EBVC) scheme. For the CVFEM, element polynomial, P , promotion is used to extend the low-order P = 1 CVFEM method to higher-order, i.e., P = 2. An equal-order low-Mach pressure-stabilized methodology, with emphasis on the non-conformal interface boundary condition, is presented. A fully implicit matrix solver approach that accounts for the full stencil connectivity across the non-conformal interface is employed. A complete suite of formal verification studies using the method of manufactured solutions (MMS) is performed to verify the order of accuracy of the underlying methodology. The chosen suite of analytical verification cases range from a simple steady diffusion system to a traveling viscous vortex across mixed-order non-conformal interfaces. Results from all verification studies demonstrate either second- or third-order spatial accuracy and, for transient solutions, second-order temporal accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Enhanced information content for three-dimensional localization and tracking using the double-helix point spread function with variable-angle illumination epifluorescence microscopy

The signal-to-noise ratio (SNR) and three-dimensional localization precision of a double helix point spread function (DH-PSF) can be significantly improved by applying variable-angle illumination epifluorescence microscopy (VAI, also commonly known as “pseudo-TIRF” or “quasi-TIRF”). In this paper, we performed a quantitative analysis of the dependence of SNR and localization precision on the number of measured photons and the incident angle for static particles under both low (at a planar index-matched interface) and high (within a porous silica matrix) fluorescent background conditions. We found that under noisier imaging conditions, the SNR and localization precision obtained using VAI are up to fivefold and threefold greater, respectively, than those obtained using epi-illumination. Finally and moreover, we demonstrate that the combination of DH-PSF and VAI can significantly improve the accuracy of the measured diffusion coefficient for mobile particles, even at a relatively large distance (50 μm) from the boundary of the optical cell.

42 ENGINEERING↗
DOE OSTI2025

Code Coverage Status of the ARC Code PERSENT

The Argonne Reactor Code (ARC) software system supports users in their fast reactor design goals by providing neutronic, thermal-hydraulic, and structural analysis capabilities. PERSENT fulfills the role of generating reactivity coefficients for a given time point of a REBUS calculation usable in a point kinetics based safety analysis capability. PERSENT also provides a sensitivity coefficient capability on eigenvalue, reactivity worth, and several other key coefficients that are used in the follow-on safety analysis. Given a co-variance matrix, PERSENT can carry out the uncertainty quantification to indicate the amount of error in the reactivity coefficients derived from the errors in the cross section measurements. With continued improvement of computational resources, many of the geometry modeling capabilities in DIF3D that were primarily used in low order schemes are not really needed anymore. Today, the diffusion and transport capabilities of DIF3D-VARIANT are primarily used in the reactor design process with some scattered usage of DIF3D-FD and DIF3D-Nodal. PERSENT is part of the ARC code system and is built around DIF3D-VARIANT and the flux solution it provides. The purpose of the present work is to identify a set of test problems for PERSENT and assess the code coverage of PERSENT for those test problems. PERSENT treats the DIF3D executable as an external executable and thus the code coverage considerations only need to focus on the PERSENT source code and only a fraction of the connected modules in the existing ARC software library. The goal is to document what parts of the existing PERSENT code are touched by the set of test problems and which are not. Because the verification work done on PERSENT was focused on the most common uses of PERSENT for fast reactor analysis, the code coverage assessment of those capabilities is the highest priority. This will ensure that nothing is being missed by the existing verification test problems that users of PERSENT rely upon. The code coverage analysis of PERSENT was performed with the Code Coverage Tool of the Intel Fortran compiler which requires modifications to the compilation of PERSENT. The detailed coverage tables are given for each submodule of PERSENT. Most of the uncovered parts/files could be easily ignored because they are either for error message and debugging output or not needed by PERSENT today. Only a few uncovered parts of PERSENT deserve extending the verification test suite.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Effects of Fe atoms on hardening of a nickel matrix: Nanoindentation experiments and atom-scale numerical modeling

Significant hardening effect due to Fe concentrations in Ni-based alloys with face-centered-cubic structure has been studied by using a combined experimental and atomistic-based computational approach via nanoindentation tests. The obtained experimental load–displacement data for the [0 0 1] crystal orientation reached a qualitative good agreement with molecular dynamics simulations results, leading to strong evidence that the main strengthening factors are associated to sluggish dislocation diffusion, reduced defect sizes and the nucleation of tetrahedral stacking faults. Here, interstitial type prismatic dislocation loops mainly formed by $\frac{1}{6}$$\langle$1 1 2$\rangle$ Shockley dislocations are nucleated during the loading process, where their interaction leads to the formation of pyramidal shaped stacking fault which are mainly created by $\frac{1}{3}$$\langle$1 0 0$\rangle$ Hirth dislocations lines. Observing both types of defects coexisting in the same plastic deformation zone by both approaches. Reported mechanical data, measured experimentally and interpreted numerically, are also in accordance with microstructural SEM and TEM investigations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Operator-level quantum acceleration of non-logconcave sampling

Sampling from probability distributions of the form 𝝈 ∝ e −𝜷V , where V is a continuous potential, is a fundamental task across physics, chemistry, biology, computer science, and statistics. However, when V is nonconvex, the resulting distribution becomes non-logconcave, and classical methods such as Langevin dynamics often exhibit poor performance. We introduce a quantum algorithm that provably accelerates a broad class of continuous-time sampling dynamics. For Langevin dynamics, our method encodes the target Gibbs measure into the amplitudes of aquantum state, identified as the kernel of a block matrix derived from a factorization of the Witten Laplacian operator. This connection enables Gibbs sampling via singular value thresholding and yields up to a quartic quantum speedup over best-knownclassical Langevin-based methods in the non-logconcave setting. Building on this framework, we further develop the first quantum algorithm that accelerates replica exchange Langevin diffusion, a widely used method for sampling from complex, rugged energy landscapes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Modified Dirac Hamiltonian for efficient quantum mechanical simulations of micron sized devices

Representing massless Dirac fermions on a spatial lattice poses a potential challenge known as the Fermion Doubling problem. Addition of a quadratic term to the Dirac Hamiltonian provides a possible way to circumvent this problem. We show that the modified Hamiltonian with the additional term results in a very small Hamiltonian matrix when discretized on a real space square lattice. The resulting Hamiltonian matrix is considerably more efficient for numerical simulations without sacrificing on accuracy and is several orders of magnitude faster than the atomistic tight binding model. Using this Hamiltonian and the non-equilibrium Green's function formalism, we show several transport phenomena in graphene, such as magnetic focusing, chiral tunneling in the ballistic limit, and conductivity in the diffusive limit in micron sized graphene devices. The modified Hamiltonian can be used for any system with massless Dirac fermions such as Topological Insulators, opening up a simulation domain that is not readily accessible otherwise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Advanced Methodology and Preliminary Measurements of Molecular and Mechanical Properties of Heart Valves under Dynamic Strain

Mammalian heart valves are soft tissue assemblies with multi-scale material properties. This is because they are constructs comprising both muscle and non-contractile extracellular matrix proteins (such as collagens and proteoglycans) and transition regions where one form of tissue structure becomes another, significantly different form. The leaflets of the mitral and tricuspid valves are connected to chordae tendinae which, in turn, bind through papillary muscles to the cardiac wall of the ventricle. The transition regions between these tissue subsets are complex and diffuse. Their material composition and mechanical properties have not been previously described with both micro and nanoscopic data recorded simultaneously, as reported here. Annotating the mechanical characteristics of these tissue transitions will be of great value in developing novel implants, improving the state of the surgical simulators and advancing robot-assisted surgery. We present here developments in multi-scale methodology that produce data that can relate mechanical properties to molecular structure using scanning X-ray diffraction. We correlate these data to corresponding tissue level (macro and microscopic) stress and strain, with particular emphasis on the transition regions and present analyses to indicate points of possible failure in these tissues.

valve failure↗
DOE OSTI2016

Fast Multilevel Solvers for a Class of Discrete Fourth Order Parabolic Problems

In this paper, we study fast iterative solvers for the solution of fourth order parabolic equations discretized by mixed finite element methods. We propose to use consistent mass matrix in the discretization and use lumped mass matrix to construct efficient preconditioners. We provide eigenvalue analysis for the preconditioned system and estimate the convergence rate of the preconditioned GMRes method. Furthermore, we show that these preconditioners only need to be solved inexactly by optimal multigrid algorithms. Our numerical examples indicate that the proposed preconditioners are very efficient and robust with respect to both discretization parameters and diffusion coefficients. We also investigate the performance of multigrid algorithms with either collective smoothers or distributive smoothers when solving the preconditioner systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The j and k dependencies of the line coupling and line mixing effects: Theoretical studies of the relaxation matrices of N 2 -broadened CH 3 D

Line coupling and line mixing effects in parallel and perpendicular bands of CH 3 D perturbed by N 2 have been studied. The work focuses on exhibiting the j and k dependencies of these two processes. The calculations were based on a previously reported anisotropic intermolecular potential including both the long-range multipole, induction, and dispersion forces and a short-range atom-atom model. Throughout this work, it is shown that components with L 1 = 3 of the atom-atom model are dominant contributions to the diffusion operator. As a consequence, in comparison with other molecular systems such as the CH 3 Cl-N 2 and CH 3 I-N 2 , theoretically predicted line coupling and line mixing effects exhibit completely new j and k dependencies. In general, the theoretically calculated halfwidths and intra-doublets’ off-diagonal elements of the relaxation matrix are in reasonable agreement with measurements.

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

Elements distribution and deformation features of a 4.5% Re nickel-based single crystal superalloy during creep at high temperature

Highlights: • Atom probe is used to investigate the element concentration changes in γ and γ' phase during creep. • Mo, W and Re diffuse from γ' to γ phase and enrich in γ/γ' interface during high temperature creep. • Super-dislocations in γ' phase may cross-slip from {111} to {100} plane due to the elements diffusions. • Dislocations shearing into γ' phase may decompose into K-W locks plus APB due to the refractory elements enrichment. - Abstract: By means of creep property measurement, microstructure observation and atom probe tomography, the elements distribution and deformation features of a 4.5% Re nickel-based single crystal superalloy during creep are investigated. Results show that after fully heat treated, the elements Al and Ta are mainly distributed in γ′ phase of alloy. After the alloy is crept up to fracture at 1100 °C/137 MPa, the change of the elements distributing in γ and γ′ phases occurs, wherein, the concentration decrement of the elements Al and Cr in γ′ and γ phases is attributed to their oxidation which may consume the ones. Moreover, some of the atoms Re, W and Mo dissolving in γ′ phase are repelled into γ matrix to reduce its internal energy. And the peak concentration of the atoms Re and W appears in the γ phase near the transition region, the atoms Re, W in the peak region cause the lattice distortion to increase the resistance of dislocations motion, which may improve the creep resistance of alloy. In the later stage of creep, dislocations shearing into the rafted γ′ phase may cross-slip from {111} to {100} plane to form the configuration of K-W locks, which may restrain the slipping of dislocations to improve the creep resistance of alloy. Wherein, the atom Re with lower diffusion coefficient may delay the dislocations in K-W locks for slipping and cross-slipping, which is thought to be the main reason of the K-W locks keeping in 4.5%Re alloy during creep at 1100 °C.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Annealing temperature and barrier thickness effect on the structural and optical properties of silicon nanocrystals/SiO₂ superlattices

The effect of the annealing temperature and the SiO₂ barrier thickness of silicon nanocrystal (NC)/SiO₂ superlattices (SLs) on their structural and optical properties is investigated. Energy-filtered transmission electron microscopy (TEM) revealed that the SL structure is maintained for annealing temperatures up to 1150 °C, with no variation on the nanostructure morphology for different SiO₂ barrier thicknesses. Nevertheless, annealing temperatures as high as 1250 °C promote diffusion of Si atoms into the SiO₂ barrier layers, which produces larger Si NCs and the loss of the NC size control expected from the SL approach. Complementary Raman scattering measurements corroborated these results for all the SiO₂ and Si-rich oxynitride layer thicknesses. In addition, we observed an increasing crystalline fraction up to 1250 °C, which is related to a decreasing contribution of the suboxide transition layer between Si NCs and the SiO₂ matrix due to the formation of larger NCs. Finally, photoluminescence measurements revealed that the emission of the superlattices exhibits a Gaussian-like lineshape with a maximum intensity after annealing at 1150 °C, indicating a high crystalline degree in good agreement with Raman results. Samples submitted to higher annealing temperatures display a progressive emission broadening, together with an increase in the central emission wavelength. Both effects are related to a progressive broadening of the size distribution with a larger mean size, in agreement with TEM observations. On the other hand, whereas the morphology of the Si NCs is unaffected by the SiO₂ barrier thickness, the emission properties are slightly modified. These observed modifications in the emission lineshape allow monitoring the precipitation process of Si NCs in a direct non-destructive way. All these experimental results evidence that an annealing temperature of 1150 °C and 1-nm SiO₂ barrier can be reached whilst preserving the SL structure, being thus the optimal structural SL parameters for their use in optoelectronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗