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 469 records · Page 26

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↗

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↗

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↗

Oxidation behaviors of matrix-grade graphite during water vapor ingress accidents for high temperature gas-cooled reactors

Water ingress into cores of high temperature gas-cooled reactors (HTGRs) can cause serious safety problems in graphitic components, especially for matrix graphite, which is more vulnerable than nuclear graphite. Here in this study, oxidation behaviors of a matrix-grade graphite, ARB-B1, were investigated under simulated water ingress accidental conditions. The oxidation tests were conducted in water vapor/helium mixed atmospheres at elevated temperatures up to 1200 °C. Oxidation rates and activation energy were evaluated by mass loss measurements. The activation energy was 121.5 kJ/mol in the considered temperature range, indicating only one oxidation regime-in-pore diffusion assisted by cracks. The cracks originated from the fabrication process, acting as penetration paths for the oxidant, played an important role in the oxidation behaviors. Ungraphitized binder oxidation occurred at all temperatures, while apparent oxidation of filler particles started from 1100 °C. Moreover, oxidation effects on compressive properties and hardness of the matrix graphite were discussed in terms of oxidized microstructures. This first detailed study on matrix graphite oxidation in water vapor provides important kinetics knowledge and microstructural and mechanical data for accidental scenarios of HTGRs as well as for development of new matrix graphite materials.

36 MATERIALS SCIENCE↗

Contaminant source identification using semi-supervised machine learning

Identification of the original groundwater types present in geochemical mixtures observed in an aquifer is a challenging but very important task. Frequently, some of the groundwater types are related to different infiltration and/or contamination sources associated with various geochemical signatures and origins. The characterization of groundwater mixing processes typically requires solving complex inverse models representing groundwater flow and geochemical transport in the aquifer, where the inverse analysis accounts for available site data. Usually, the model is calibrated against the available data characterizing the spatial and temporal distribution of the observed geochemical types. Numerous different geochemical constituents and processes may need to be simulated in these models which further complicates the analyses. In this paper, we propose a new contaminant source identification approach that performs decomposition of the observation mixtures based on Non-negative Matrix Factorization (NMF) method for Blind Source Separation (BSS), coupled with a custom semi-supervised clustering algorithm. Our methodology, called NMFk, is capable of identifying (a) the unknown number of groundwater types and (b) the original geochemical concentration of the contaminant sources from measured geochemical mixtures with unknown mixing ratios without any additional site information. NMFk is tested on synthetic and real-world site data. Finally, the NMFk algorithm works with geochemical data represented in the form of concentrations, ratios (of two constituents; for example, isotope ratios), and delta notations (standard normalized stable isotope ratios).

54 ENVIRONMENTAL SCIENCES↗

Water vapor sorption isotherms, pore structure, and moisture transport characteristics of alkali-activated and Portland cement-based binders

Moisture transport plays a key role in determining the different durability-related features of cementitious materials. In this paper, moisture sorption in a range of low-calcium (geopolymer-type) and calcium-rich alkali-activated binders are studied and compared with that of Portland cement-based binders. Through the analysis of water vapor sorption isotherms (WVSI) and mercury intrusion porosimetry (MIP) test results, two vastly different pore structures were observed. Fly ash-based geopolymer-type binders showed a very porous structure where a large volume of mesopores coexisted with a significant volume of macropores. Alkali-activated slag binders, however, had a very fine pore structure, with a relative lack of large macropores. The different pore structure of fly ash-based and slag-based binders led to amplification of pore blocking and cavitation in blended systems by the addition of slag. Analysis of the sorption kinetics showed the prominent effect of the presence of calcium in the matrix to reduce the permeability.

36 MATERIALS SCIENCE↗

Calculation of Transfer Functions of Multilayer Biotissues in the Problems of Correction of Their Fluorescence Spectra

A method for rapid calculation of a flux of stimulated fluorescence of a multilayer optically dense medium with inhomogeneous distribution of the fluorophore has been developed. The light field in the medium at the excitation wavelength of fluorescence is represented by a superposition of incident collimated, incident diffuse, and reflected diffuse fluxes. A two-stream approximation is used to describe the light field in the medium at the wavelength of emission of the fluorescence. Fluxes in adjacent elementary layers of the medium and on its surface are connected by simple matrix operators that are obtained using a combination of engineering approaches of radiation-transfer theory and single-scattering approximation. The calculations of fluorescence fluxes of a four-layer biotissue that are excited and recorded at 400–800 nm are compared with their Monte Carlo simulation with a discrepancy of 1%. The effect of the propagation medium on the fluorescence spectra of 5-ALA-induced protoporphyrin IX that are recorded from human skin was studied, and a technique for their correction that is based on measurements and quantitative analysis of the diffuse reflectance spectrum of skin was proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Molecular dynamics insight to phase transition in n-alkanes with carbon nanofillers

The present work aims to investigate the phase transition, dispersion and diffusion behavior of nanocomposites of carbon nanotube (CNT) and straight chain alkanes. These materials are potential candidates for organic phase change materials(PCMs) and have attracted flurry of research recently. Accurate experimental evaluation of the mass, thermal and transport properties of such composites is both difficult as well as economically taxing. Additionally it is crucial to understand the factors that results in modification or enhancement of their characteristic at atomic or molecular level. Classical molecular dynamics approach has been extended to elucidate the same. Bulk atomistic models have been generated and subjected to rigorous multistage equilibration. To reaffirm the approach, both canonical and constant-temperature, constant- pressure ensembles were employed to simulate the models under consideration. Explicit determination of kinetic, potential, non-bond and total energy assisted in understanding the enhanced thermal and transport property of the nanocomposites from molecular point of view. Crucial parameters including mean square displacement and simulated self diffusion coefficient precisely define the balance of the thermodynamic and hydrodynamic interactions. Radial distribution function also reflected the density variation, strength and mobility of the nanocomposites. It is expected that CNT functionalization could improve the dispersion within n-alkane matrix. This would further ameliorate the mass and thermal properties of the composite. Additionally, the determined density was in good agreement with experimental data. Thus, molecular dynamics can be utilized as a high throughput technique for theoretical investigation of nanocomposites PCMs.

36 MATERIALS SCIENCE↗

Platinum-Coated Hollow Graphene Nanocages as Cathode Used in Lithium-Oxygen Batteries

One of the formidable challenges facing aprotic lithium-oxygen (Li-O-2) batteries is the high charge overpotential, which induces the formation of byproducts, loss in efficiency, and poor cycling performance. Herein, the synthesis of the ultrasmall Pt-coated hollow graphene nano cages as cathode in Li-O-2 batteries is reported. The charge voltage plateau can reduce to 3.2 V at the current density of 100 mA g(-1), even maintain below 3.5 V when the current density increased to 500 mA g(-1). The unique hollow graphene nanocages matrix can not only provide numerous nanoscale tri-phase regions as active sites for efficient oxygen reduction, but also offer sufficient amount of mesoscale pores for rapid oxygen diffusion. Furthermore, with strong atomic-level oxygen absorption into its subsurface, ultrasmall Pt catalytically serves as the nucleation site for Li2O2 growth. The Li2O2 is subsequently induced into a favorable form with small size and amorphous state, decomposed more easily during recharge. Meanwhile, the conductive hollow graphene substrate can enhance the catalytic activity of noble metal Pt catalysts due to the graphene-metal interfacial interaction. Benefiting from the above synergistic effects between the hollow graphene nanocages and the nanosized Pt catalysts, the ultrasmall Pt-decorated graphene nanocage cathode exhibits enhanced electrochemical performances.

Wu, Feng↗

The Effect of Multistage Aging on Mechanical Properties and Microstructure of Forged 7050 Aluminum Alloys

A new four-stage aging treatment, which is composed of pre-aging, retrogression, natural aging (NA) and re-aging, is designed to increase the tensile properties and the fatigue life of forged 7050 aluminum alloy in this work. As compared to traditional T74 and RRA treatments, the four-stage aging-treated sample exhibits the highest tensile properties and ductility. The tensile strength and yield strength are 627 and 608 MPa, respectively, and the elongation is 11.3%. The results of TEM and DSC show that an increased G.P zones and fine η′ phase distributed in the matrix of the alloy are subjected by this new aging processing as compared with T74 and RRA treatments. The natural aging in four-stage aging treatment accelerates the diffusion of Mg vacancies, and promotes the precipitation of G.P zones in grain interior during re-aging stage. On the basis of dislocation shearing mechanism, more dense and fine G.P zones of four-stage aging temper result in the superior tensile strength as compared to T74 and RRA temper. In addition, a combination of shearable precipitates, fine GBPs and narrow PFZs is beneficial to the slip reversibility in the plastic zone at the crack tip, thereby reducing fatigue damage accumulation and increasing the fatigue life.

36 MATERIALS SCIENCE↗

A cellular automaton method to simulate the microstructure and evolution of low-enriched uranium (LEU) U–Mo/Al dispersion type fuel plates

Low-enriched uranium (LEU) fuel plates for high power materials test reactors (MTR) are composed of nominally spherical uranium-molybdenum (U-Mo) particles within an aluminum matrix. Fresh U-Mo particles typically range between 10 and 100 mu m in diameter, with particle volume fractions up to 50%. As the fuel ages, reaction-diffusion processes cause the formation and growth of interaction layers that surround the fuel particles. The growth rate depends upon the temperature and radiation environment. The cellular automaton algorithm described in this paper can synthesize realistic random fuel-particle structures and simulate the growth of the intermetallic interaction layers. Examples in the present paper pack approximately 1000 particles into three-dimensional rectangular fuel structures that are approximately 1 mm on each side. The computational approach is designed to yield synthetic microstructures consistent with images from actual fuel plates and is validated by comparison with empirical data on actual fuel plates. (C) 2014 Elsevier B.V. All rights reserved.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Processing and microstructure of ZrB 2 –SiC composite prepared by reactive spark plasma sintering

In-situ formation of ZrB 2 –SiC composites was investigated by reactive spark plasma sintering of precursor powders according to the reaction B4C + 2ZrC + 3Si → 3SiC + 2ZrB 2 . The reaction and process presented here involves a diffusion reaction between B4C and ZrC which facilitates the formation of ZrB 2 , while liquid phase sintering of silicon facilitates atomic diffusion and combines with free C from the B4C and ZrC reaction to form SiC within minutes of heating and there were some residual unreacted precursor materials. An interpenetrating matrix of ZrB 2 –SiC was formed that shows increased fracture toughness (6.03 ± 0.45 MPa m 1/2 ) despite relatively low density (95 %).

36 MATERIALS SCIENCE↗

Morphology and Flow Behavior of Cellulose Nanofibers Dispersed in Glycols

Understanding the morphology and flow behavior of cellulose nanofibers (CNFs) dispersed in organic solvents can improve the process of fabricating new cellulose-based nanocomposites. In this study, jute-based 2,2,6,6-tetramethylpiperidinyl-1-oxyl (TEMPO)-oxidized CNFs with two different charge densities (0.64 and 1.03 mmol/g) were dispersed in ethylene glycol (EG) and propylene glycol (PG) using the solvent exchange method. The morphology and dimensions of CNFs in dry and suspension states were characterized using transmission electron microscopy, atomic force microscopy, and small-angle X-ray scattering techniques. The results showed that the cross-sectional dimensions remained the same in different solvents. Rheological measurements showed that CNF suspensions in water or glycol (EG and PG) behaved similar to typical polymer solutions with a solvent-independent overlap concentration corresponding to the crowding factor of about 14. Moreover, a thixotropic behavior was found in the concentrated CNF/glycol systems as determined in typical CNF aqueous suspensions. The fact that TEMPO-oxidized CNFs can be well dispersed in organic solvents opens up new possibilities to improve the CNF–polymer matrix blending, where the use of a viscous solvent can delay the transition to turbulence in processing and improve the control of fiber orientation because of a slower Brownian diffusive motion.

36 MATERIALS SCIENCE↗

Neutrino physics with JUNO

The Jiangmen Underground Neutrino Observatory (JUNO), a 20 kton multi-purpose underground liquid scintillator detector, was proposed with the determination of the neutrino mass hierarchy (MH) as a primary physics goal. The excellent energy resolution and the large fiducial volume anticipated for the JUNO detector offer exciting opportunities for addressing many important topics in neutrino and astro-particle physics. In this document, we present the physics motivations and the anticipated performance of the JUNO detector for various proposed measurements. Following an introduction summarizing the current status and open issues in neutrino physics, we discuss how the detection of antineutrinos generated by a cluster of nuclear power plants allows the determination of the neutrino MH at a 3–4σ significance with six years of running of JUNO. The measurement of antineutrino spectrum with excellent energy resolution will also lead to the precise determination of the neutrino oscillation parameters ${\mathrm{sin}}^{2}{\theta }_{12}$, ${\rm{\Delta }}{m}_{21}^{2}$, and $| {\rm{\Delta }}{m}_{{ee}}^{2}| $ to an accuracy of better than 1%, which will play a crucial role in the future unitarity test of the MNSP matrix. The JUNO detector is capable of observing not only antineutrinos from the power plants, but also neutrinos/antineutrinos from terrestrial and extra-terrestrial sources, including supernova burst neutrinos, diffuse supernova neutrino background, geoneutrinos, atmospheric neutrinos, and solar neutrinos. As a result of JUNO's large size, excellent energy resolution, and vertex reconstruction capability, interesting new data on these topics can be collected. For example, a neutrino burst from a typical core-collapse supernova at a distance of 10 kpc would lead to ~5000 inverse-beta-decay events and ~2000 all-flavor neutrino–proton ES events in JUNO, which are of crucial importance for understanding the mechanism of supernova explosion and for exploring novel phenomena such as collective neutrino oscillations. Detection of neutrinos from all past core-collapse supernova explosions in the visible universe with JUNO would further provide valuable information on the cosmic star-formation rate and the average core-collapse neutrino energy spectrum. Antineutrinos originating from the radioactive decay of uranium and thorium in the Earth can be detected in JUNO with a rate of ~400 events per year, significantly improving the statistics of existing geoneutrino event samples. Atmospheric neutrino events collected in JUNO can provide independent inputs for determining the MH and the octant of the ${\theta }_{23}$ mixing angle. Detection of the 7Be and 8B solar neutrino events at JUNO would shed new light on the solar metallicity problem and examine the transition region between the vacuum and matter dominated neutrino oscillations. Regarding light sterile neutrino topics, sterile neutrinos with ${10}^{-5}\;{{\rm{eV}}}^{2}\lt {\rm{\Delta }}{m}_{41}^{2}\lt {10}^{-2}\;{{\rm{eV}}}^{2}$ and a sufficiently large mixing angle ${\theta }_{14}$ could be identified through a precise measurement of the reactor antineutrino energy spectrum. Meanwhile, JUNO can also provide us excellent opportunities to test the eV-scale sterile neutrino hypothesis, using either the radioactive neutrino sources or a cyclotron-produced neutrino beam. The JUNO detector is also sensitive to several other beyondthe-standard-model physics. Examples include the search for proton decay via the $p\to {K}^{+}+\bar{\nu }$ decay channel, search for neutrinos resulting from dark-matter annihilation in the Sun, search for violation of Lorentz invariance via the sidereal modulation of the reactor neutrino event rate, and search for the effects of non-standard interactions. The proposed construction of the JUNO detector will provide a unique facility to address many outstanding crucial questions in particle and astrophysics in a timely and cost-effective fashion. Furthermore, it holds the great potential for further advancing our quest to understanding the fundamental properties of neutrinos, one of the building blocks of our Universe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Space-Time Block Preconditioning for Incompressible Flow

Parallel-in-time methods have become increasingly popular in the simulation of time-dependent numerical PDEs, allowing for the efficient use of additional message passing interface processes when spatial parallelism saturates. Most methods treat the solution and parallelism in space and time separately. In contrast, all-at-once methods solve the full space-time system directly, largely treating time as simply another spatial dimension. All-at-once methods offer a number of benefits over separate treatment of space and time, most notably significantly increased parallelism and faster time to solution (when applicable). However, the development of fast, scalable all-at-once methods has largely been limited to time-dependent (advection-)diffusion problems. This paper introduces the concept of space-time block preconditioning for the all-at-once solution of incompressible flow. By extending well-known concepts of spatial block preconditioning to the space-time setting, we develop a block preconditioner whose application requires the solution of a space-time (advection-)diffusion equation in the velocity block, coupled with a pressure Schur complement approximation consisting of independent spatial solves at each time-step, and a space-time matrix-vector multiplication. The new method is tested on four classical models in incompressible flow. Finally, the results indicate perfect scalability in refinement of spatial and temporal mesh spacing, perfect scalability in nonlinear Picard iteration count when applied to a nonlinear Navier--Stokes problem, and minimal overhead in terms of number of preconditioner applications compared with sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗

Final Report - Subcontract B623760

During my visit to LLNL during July 17{27, 2017, 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 12. 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 de nite 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 ! 1, 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 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↗

LEU-Mo Fuel Out-of-Pile Characterization for TUM (Final Report)

The Technische Universität München (TUM) located in Germany requested that Pacific Northwest National Laboratory (PNNL) perform thermal-physical property measurement services for determining the thermal conductivities of irradiated U-Mo dispersion fuel samples. Two dispersion fuel segments containing U-7Mo dispersion fuel in an Al-2wt% Si matrix harvested from the AFIP-1 irradiation experiment were provided to PNNL for this effort. The segments were sectioned into samples for pycnometry (density) and laser flash analysis (thermal diffusivity), differential scanning calorimetry (specific heat capacity), and optical metallography (thickness, image analysis). This report summarizes the findings of the analysis performed at TUM’s request.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

TU-CD-BRA-01: A Novel 3D Registration Method for Multiparametric Radiological Images

Purpose: Multiparametric and multimodality radiological imaging methods, such as, magnetic resonance imaging(MRI), computed tomography(CT), and positron emission tomography(PET), provide multiple types of tissue contrast and anatomical information for clinical diagnosis. However, these radiological modalities are acquired using very different technical parameters, e.g.,field of view(FOV), matrix size, and scan planes, which, can lead to challenges in registering the different data sets. Therefore, we developed a hybrid registration method based on 3D wavelet transformation and 3D interpolations that performs 3D resampling and rotation of the target radiological images without loss of information Methods: T1-weighted, T2-weighted, diffusion-weighted-imaging(DWI), dynamic-contrast-enhanced(DCE) MRI and PET/CT were used in the registration algorithm from breast and prostate data at 3T MRI and multimodality(PET/CT) cases. The hybrid registration scheme consists of several steps to reslice and match each modality using a combination of 3D wavelets, interpolations, and affine registration steps. First, orthogonal reslicing is performed to equalize FOV, matrix sizes and the number of slices using wavelet transformation. Second, angular resampling of the target data is performed to match the reference data. Finally, using optimized angles from resampling, 3D registration is performed using similarity transformation(scaling and translation) between the reference and resliced target volume is performed. After registration, the mean-square-error(MSE) and Dice Similarity(DS) between the reference and registered target volumes were calculated. Results: The 3D registration method registered synthetic and clinical data with significant improvement(p<0.05) of overlap between anatomical structures. After transforming and deforming the synthetic data, the MSE and Dice similarity were 0.12 and 0.99. The average improvement of the MSE in breast was 62%(0.27 to 0.10) and prostate was 63%(0.13 to 0.04;p<0.05). The Dice similarity was in breast 8%(0.91 to 0.99) and for prostate was 89%(0.01 to 0.90;p<0.05) Conclusion: Our 3D wavelet hybrid registration approach registered diverse breast and prostate data of different radiological images(MR/PET/CT) with a high accuracy.

60 APPLIED LIFE SCIENCES↗