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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Novel NASICON‐Type Na 4 MnCr(PO 4 ) 3 Demonstrating the Energy Density Record of Phosphate Cathodes for Sodium‐Ion Batteries

Abstract Sodium‐ion batteries (SIBs) have attracted incremental attention as a promising candidate for grid‐scale energy‐storage applications. To meet practical requirements, searching for new cathode materials with high energy density is of great importance. Herein, a novel Na superionic conductor (NASICON)‐type Na 4 MnCr(PO 4 ) 3 is developed as a high‐energy cathode for SIBs. The Na 4 MnCr(PO 4 ) 3 nanoparticles homogeneously embedded in a carbon matrix can present an extraordinary reversible capacity of 160.5 mA h g −1 with three‐electron reaction at ≈3.53 V during the Na + extraction/insertion process, realizing an unprecedentedly high energy density of 566.5 Wh kg −1 in the phosphate cathodes for SIBs. It is intriguing to reveal the underlying mechanism of the unique Mn 2+ /Mn 3+ , Mn 3+ /Mn 4+ , and Cr 3+ /Cr 4+ redox couples via X‐ray absorption near‐edge structure spectroscopy. The whole electrochemical reaction undergoes highly reversible single‐phase and biphasic transitions with a moderate volume change of 7.7% through in situ X‐ray diffraction and ex situ high‐energy synchrotron X‐ray diffraction. Combining density functional theory (DFT) calculations with the galvanostatic intermittent titration technique, the superior performance is ascribed to the low ionic‐migration energy barrier and desirable Na‐ion diffusion kinetics. The present work can offer a new insight into the design of multielectron‐reaction cathode materials for SIBs.

Zhang, Jian↗
DOE OSTI2014

Feature-space assessment of electrical impedance tomography coregistered with computed tomography in detecting multiple contrast targets

Purpose: Fusion of electrical impedance tomography (EIT) with computed tomography (CT) can be useful as a clinical tool for providing additional physiological information about tissues, but requires suitable fusion algorithms and validation procedures. This work explores the feasibility of fusing EIT and CT images using an algorithm for coregistration. The imaging performance is validated through feature space assessment on phantom contrast targets. Methods: EIT data were acquired by scanning a phantom using a circuit, configured for injecting current through 16 electrodes, placed around the phantom. A conductivity image of the phantom was obtained from the data using electrical impedance and diffuse optical tomography reconstruction software (EIDORS). A CT image of the phantom was also acquired. The EIT and CT images were fused using a region of interest (ROI) coregistration fusion algorithm. Phantom imaging experiments were carried out on objects of different contrasts, sizes, and positions. The conductive medium of the phantoms was made of a tissue-mimicking bolus material that is routinely used in clinical radiation therapy settings. To validate the imaging performance in detecting different contrasts, the ROI of the phantom was filled with distilled water and normal saline. Spatially separated cylindrical objects of different sizes were used for validating the imaging performance in multiple target detection. Analyses of the CT, EIT and the EIT/CT phantom images were carried out based on the variations of contrast, correlation, energy, and homogeneity, using a gray level co-occurrence matrix (GLCM). A reference image of the phantom was simulated using EIDORS, and the performances of the CT and EIT imaging systems were evaluated and compared against the performance of the EIT/CT system using various feature metrics, detectability, and structural similarity index measures. Results: In detecting distilled and normal saline water in bolus medium, EIT as a stand-alone imaging system showed contrast discrimination of 47%, while the CT imaging system showed a discrimination of only 1.5%. The structural similarity index measure showed a drop of 24% with EIT imaging compared to CT imaging. The average detectability measure for CT imaging was found to be 2.375 ± 0.19 before fusion. After complementing with EIT information, the detectability measure increased to 11.06 ± 2.04. Based on the feature metrics, the functional imaging quality of CT and EIT were found to be 2.29% and 86%, respectively, before fusion. Structural imaging quality was found to be 66% for CT and 16% for EIT. After fusion, functional imaging quality improved in CT imaging from 2.29% to 42% and the structural imaging quality of EIT imaging changed from 16% to 66%. The improvement in image quality was also observed in detecting objects of different sizes. Conclusions: The authors found a significant improvement in the contrast detectability performance of CT imaging when complemented with functional imaging information from EIT. Along with the feature assessment metrics, the concept of complementing CT with EIT imaging can lead to an EIT/CT imaging modality which might fully utilize the functional imaging abilities of EIT imaging, thereby enhancing the quality of care in the areas of cancer diagnosis and radiotherapy treatment planning.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2017

Two-Level Convergence Theory for Multigrid Reduction in Time (MGRIT)

Here we develop a two-grid convergence theory for the parallel-in-time scheme known as multigrid reduction in time (MGRIT), as it is implemented in the open-source package [XBraid: Parallel Multigrid in Time, http://llnl.gov/casc/xbraid]. MGRIT is a scalable and multilevel approach to parallel-in-time simulations that nonintrusively uses existing time-stepping schemes, and in a specific two-level setting it is equivalent to the widely known parareal algorithm. The goal of this paper is twofold. First, we present a two-level MGRIT convergence analysis for linear problems where the spatial discretization matrix can be diagonalized, and then apply this analysis to our two basic model problems, the heat equation and the advection equation. One important assumption is that the coarse and fine time-grid propagators can be diagaonalized by the same set of eigenvectors, which is often the case when the same spatial discretization operator is used on the coarse and fine time grids. In many cases, the MGRIT algorithm is guaranteed to converge, and we demonstrate numerically that the theoretically predicted convergence rates are sharp in practice for our model problems. Second, we explore how the convergence of MGRIT compares to the stability of the chosen time-stepping scheme. In particular, we demonstrate that a stable time-stepping scheme does not necessarily imply convergence of MGRIT, although MGRIT with FCF-relaxation always converges for the diffusion dominated problems considered here.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Regenerative Carbonate-Based Thermochemical Energy Storage System for Concentrating Solar Power

Southern Research has developed a thermochemical energy storage (TCES) technology that utilizes the endothermic-exothermic reversible carbonation of calcium oxide (lime) to store thermal energy at high-temperatures, such as those achieved by next generation concentrating solar power (CSP) facilities. The major challenges addressed in the development of this system include refining a high capacity, yet durable sorbent material and designing a low thermal resistance low-cost heat exchanger reactor system to move heat between the sorbent and a heat transfer fluid under conditions relevant for CSP operation (e.g., energy density, reaction kinetics, heat flow). The proprietary stabilized sorbent was developed by Precision Combustion, Inc. (PCI). A factorial matrix of sorbent compositions covering the design space was tested using accelerated high throughput screening in a thermo-gravimetric analyzer. Several promising formulations were selected for more thorough evaluation and one formulation with high capacity (0.38 g CO 2 /g sorbent) and durability (>99.7% capacity retention over 100 cycles) was chosen as a basis for further development of the energy storage reactor system. In parallel with this effort, a full range of currently available commercial and developmental heat exchange reactor systems and sorbent loading methods were examined through literature research and contacts with commercial vendors. Process models were developed to examine if a heat exchange reactor system and balance of plant can meet required TCES performance and cost targets, optimizing tradeoffs between thermal performance, exergetic efficiency, and cost. Reactor types evaluated included many forms, from microchannel reactor, to diffusion bonded heat exchanger, to shell and tube heat exchangers. The most viable design for application to a supercritical CO 2 power cycle operating at 200-300 bar pressure and >700°C was determined to be a combination of a diffusion bonded heat exchanger with a shell and tube reactor. A bench scale reactor system was then designed and constructed to test sorbent performance under more commercially relevant conditions. This system utilizes a tube-in tube reactor design containing approximately 250 grams sorbent and is able to operate under a wide range of temperature, pressure and flow conditions as needed to explore system performance under a variety of operating conditions. A variety of sorbent loading methods may be tested using the reactor design. Initial bench test results over 25 cycles showed very high sorbent stability (>99%) and sufficient capacity (>0.28 g CO 2 /g sorbent) for an economical commercial-scale system. Initial technoeconomic evaluation of the proposed storage system show that the sorbent cost should not have a significant impact on overall system cost, and that the largest cost impacts come from the heat exchanger reactor and balance of plant equipment, including compressors and gas storage, due to the high temperatures for sCO 2 cycles. Current estimated system costs are $47/kWhth based on current material and equipment cost estimates.

14 SOLAR ENERGY↗
DOE OSTI2015

Bimetallic oxamato complexes synthesized into mesoporous matrix as precursor to tunable nanosized oxide

Highlights: • The bimetallic oxamato complexes as single-source precursor. • We prepared into a porous silica glass tunable nanosized oxide powders. • X-ray diffraction shows the formation of CeO{sub 2}/CuO and spinel cobaltite. • The different number of IDC allows control of the nanoparticle size. - Abstract: The bimetallic complexes were employed to prepare into a porous silica glass tunable nanosized oxide powders through the single source precursor (SSP) method. These materials were prepared by first anchoring of [Cu(opba)]{sup 2−} [opba = ortho-phenylenebis(oxamato)], second by reaction in situ with second metal [Co(II) or Ce(III)] and followed by a thermal treatment. The different number of impregnation–decomposition cycles (IDC) allows control of the nanoparticle size. X-ray diffraction shows the formation of mixture CeO{sub 2}–CuO and spinel copper cobaltite. Raman spectroscopy confirmed the formation of such phases. Transmission electron microscopy images revealed that spinel cobaltite particles (8 IDC) present a mean size of about 9 nm, whereas for the CeO{sub 2}–CuO phase the particle diameters are 4 nm (2 IDC) and 8 nm (6 IDC). For CeO{sub 2}–CuO the diffuse reflectance spectroscopy indicates a consistent red shift in band gap from 3.41 to 2.87 eV with increasing of particle size due to quantum confinement effect.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Microstructural and infiltration properties of woven preforms during chemical vapor infiltration

Interface-resolved direct numerical simulations (DNSs) of chemical vapor infiltration (CVI) have been performed over a range of furnace-operating conditions (Thiele moduli) and for practical woven preform geometries. A level-set method is used to resolve the geometry of the initial preform at tow scale. The interface between the vapor and solid phase is then evolved in time through the entire CVI densification cycle, fully resolving the time-varying topology between the two phases. In contrast to previous level-set methods for CVI simulation, the physical reaction and diffusion processes govern the level-set movement in the current approach. The surface deposition kinetics is described by the usual one-step model. In this paper, the DNS data are used to study the evolving porosity, surface-to-volume ratio, and flow infiltration properties (permeability and effective diffusivities). Additionally, comparisons are made to popularly-assumed structure functions and the standard, Kozeny–Carmen porous media model commonly employed in modeled CFD simulations of CVI. The virtual DNS experiments reveal a Thiele modulus and preform geometry (fabric layup) dependence which the existing microstructural and infiltration models are not able to describe throughout the entire densification process. The DNS-based, woven geometry-specific correlations can be applied directly to mean-field, furnace-scale CFD simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Assessment the leaching characteristics and long-term leaching behavior of some radionuclides from synthesized zeolite cement matrix

Highlights: • Synthesis and characterization of zeolite cement synthesized from fly ash • Leaching tests were used to assess the release rate of {sup 85}Sr and {sup 60}Co radionuclides. • Simple mathematical models were used to detect leaching mechanisms. • Long-term leaching behavior of {sup 85}Sr and {sup 60}Co from cementitious matrix was assessed. The development of cementitious materials remains a vital goal to produce valuable products with good mechanical, physical and chemical properties suitable for the safe disposal of concentrated radionuclides resulted from the treatment of contaminated solutions. In this paper, zeolite cement with good properties was laboratory synthesized from industrial by-products fly ash and characterized using various characterization techniques. The leaching characteristics of strontium and cobalt radionuclides from the synthesized zeolite cement were scrutinized conforming with the standard leaching methodology of International Atomic Energy Authority (IAEA). A mechanical strength assessment was executed to characterize the extent of immobilization process of the solidified matrices. The cumulative leaching fraction (φ) of the two studied radionuclides was found to be less than 5% in all examined conditions, which implies the applicability of the IAEA-recommended methodology for estimating the diffusion coefficient. The experimental leaching data were regressed nonlinearly to various mathematical kinetic models to assess the controlling leaching mechanism and to determine the leaching parameters. The regression results indicated that strontium ({sup 85}Sr) and cobalt ({sup 60}Co) leaching resulted from two succeeding mechanisms: as first order kinetic reaction, and then diffusion. The calculated values of leachability indices signify that the performance of each of the studied matrices is within an acceptable range. A simplified mathematical model, rooted in the first order reaction and diffusion mechanisms, was simulated to predict the radionuclides leaching rates from zeolite cement matrices. By comparing the synthesized zeolite cement with other cementitious materials, it can be concluded that the synthesized material can be classified as an efficient material suitable to immobilize {sup 85}Sr and {sup 60}Co from radioactive wastes. The acquired findings demonstrated that the studied immobilized waste matrices have acceptable mechanical effectiveness.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Monte Carlo based investigation of Berry phase for depth resolved characterization of biomedical scattering samples

The propagation of light in turbid media is an active area of research with relevance to numerous investigational fields, e.g., biomedical diagnostics and therapeutics. The statistical random-walk nature of photon propagation through turbid media is ideal for computational based modeling and simulation. Ready access to super computing resources provide a means for attaining brute force solutions to stochastic light-matter interactions entailing scattering by facilitating timely propagation of sufficient (>10million) photons while tracking characteristic parameters based on the incorporated physics of the problem. One such model that works well for isotropic but fails for anisotropic scatter, which is the case for many biomedical sample scattering problems, is the diffusion approximation. In this report, we address this by utilizing Berry phase (BP) evolution as a means for capturing anisotropic scattering characteristics of samples in the preceding depth where the diffusion approximation fails. We extend the polarization sensitive Monte Carlo method of Ramella-Roman, et al.,1 to include the computationally intensive tracking of photon trajectory in addition to polarization state at every scattering event. To speed-up the computations, which entail the appropriate rotations of reference frames, the code was parallelized using OpenMP. The results presented reveal that BP is strongly correlated to the photon penetration depth, thus potentiating the possibility of polarimetric depth resolved characterization of highly scattering samples, e.g., biological tissues.

back scattered Mueller matrix↗
DOE OSTI2017

Influence of Thermal Annealing on Free Carrier Concentration in (GaN) 1– x (ZnO) x Semiconductors

It has been previously demonstrated that the efficiency of (GaN) 1–x (ZnO) x semiconductors for solar water splitting can be improved by thermal annealing, though the origin of this improvement was not resolved. In the present work, it is shown that annealing reduces the free carrier (electron) concentration of (GaN) 1–x (ZnO) x . The time-, temperature-, and atmosphere-dependent changes were followed through two simple techniques: indirect diffuse reflectance measurements from 0.5 to 3.0 eV which show very high sensitivity to the free carrier response at the lowest energies and EPR measurements which directly probe the number of unpaired electrons. For the thermal annealing of investigated compositions, it is found that temperatures of 250 °C and below do not measurably change the free carrier concentration, a gradual reduction of the free carrier concentration occurs over a time period of many hours at 350 °C, and the complete elimination of free carriers happens within an hour at 550 °C. These changes are driven by an oxidative process which is effectively suppressed under actively reducing atmospheres (H 2 , NH 3 ) but which can still occur under nominally inert gases (N 2 , Ar). Surprisingly, it is found that the N 2 gas released during thermal oxidation of (GaN) 1–x (ZnO) x samples remains trapped within the solid matrix and is not expelled until temperatures of about 900 °C, a result directly confirmed through neutron pair-distribution fuction (PDF) measurements which show a new peak at the 1.1 Å bond length of molecular nitrogen after annealing. Preliminary comparative photoelectrochemical (PEC) measurements of the influence of free carrier concentration on photoactivity for water oxidation were carried out for a sample with x = 0.64. Samples annealed to eliminate free carriers exhibited no photoactivity for water oxidation, while a complex dependence on carrier concentration was observed for samples with intermediate free carrier concentrations. The methods demonstrated here provide an important approach for quantifying (and controlling) the carrier concentrations of semiconductors which are only available in the form of loose powders, as is commonly the case for oxynitride compounds.

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

Engineering defect energy landscape of CoCrFeNi high-entropy alloys by the introduction of additional dopants

The concept of high-entropy alloys (HEAs) focusing on tuning the overall chemical complexity represents a novel alloy design strategy. In contrast, alloying of a metallic matrix with minor doping elements with limited and localized tunability has been a common practice to improve material performance. Combining the idea of globally engineering defect energy landscape in HEAs and the localized doping strategy in dilute alloys, in this work, we explore doping effects of minor elements in a HEA matrix to further enhance the overall and localized chemical tunability, aiming to improve its irradiation resistance. Specifically, we study the influence of minor Al, Cu, Ti, and Pd substitutional doping elements on defect energetics in a CoCrFeNi model HEA based on density-functional theory (DFT) calculations. The DFT results indicate that the formation and migration energies of vacancies can be strongly influenced when a dopant is introduced at the first nearest neighbor shells around a vacancy. On the other hand, interstitial energetics are only slightly affected. Among the four elements, Ti and Pd generally decrease vacancy formation energies and increase vacancy migration energies more significantly than Al and Cu. The doping effects become more pronounced when the concentration of the substitutional dopants increases. Based on the energy distributions obtained from DFT, we build a kinetic Monte Carlo (kMC) model to assess the impact of dopants on vacancy-mediated diffusivity in the doped HEAs. Our results suggest that Ti and Pd can lower the tracer diffusivity in the considered HEAs and act as trapping sites, whereas Cu may enhance the atomic transport. This work indicates that substitutional doping in HEAs is an effective strategy in metallurgy to further tune the defect and transport properties of complex alloys.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Learning in Modal Space: Solving Time-Dependent Stochastic PDEs Using Physics-Informed Neural Networks

One of the open problems in scientific computing is the long-time integration of nonlinear stochastic partial differential equations (SPDEs), especially with arbitrary initial data. We address this problem by taking advantage of recent advances in scientific machine learning and the spectral dynamically orthogonal (DO) and borthogonal (BO) methods for representing stochastic processes. The recently introduced DO/BO methods reduce the SPDE to solving a system of deterministic PDEs and a system of stochastic ordinary differential equations. Specifically, we propose two new physics-informed neural networks (PINNs) for solving time-dependent SPDEs, namely the neural network (NN)-DO/BO methods. The proposed methods incorporate the DO/BO constraints into the loss function (along with the modal decomposition of the SPDE) with an implicit form instead of generating explicit expressions for the temporal derivatives of the DO/BO modes. Hence, the NN-DO/BO methods can overcome some of the drawbacks of the original DO/BO methods. For example, we do not need the assumption that the covariance matrix of the random coefficients is invertible as in the original DO method, and we can remove the assumption of no eigenvalue crossing as in the original BO method. Moreover, the NN-DO/BO methods can be used to solve time-dependent stochastic inverse problems with the same formulation and same computational complexity as for forward problems. Furthermore, we demonstrate the capability of the proposed methods via several numerical examples, namely: (1) A linear stochastic advection equation with deterministic initial condition: we obtain good results with the proposed methods, while the original DO/BO methods cannot be applied directly in this case. (2) Long-time integration of the stochastic Burgers' equation: we show the good performance of NN-DO/BO methods, especially the effectiveness of the NN-BO approach for such problems with many eigenvalue crossings during the whole time evolution, while the original BO method fails. (3) Nonlinear reaction diffusion equation: we consider both the forward problem and the inverse problems, including very noisy initial point values, to investigate the flexibility of the NN-DO/BO methods in handling inverse and mixed type problems. Taken together, these simulation results demonstrate that the NN-DO/BO methods can be employed to effectively quantify uncertainty propagation in a wide range of physical problems, but future work should address the efficiency issue of PINNs for forward problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Chitosan-shelled oxygen-loaded nanodroplets abrogate hypoxia dysregulation of human keratinocyte gelatinases and inhibitors: New insights for chronic wound healing

Background: : In chronic wounds, efficient epithelial tissue repair is hampered by hypoxia, and balances between the molecules involved in matrix turn-over such as matrix metalloproteinases (MMPs) and tissue inhibitors of metalloproteinases (TIMPs) are seriously impaired. Intriguingly, new oxygenating nanocarriers such as 2H,3H-decafluoropentane-based oxygen-loaded nanodroplets (OLNs) might effectively target chronic wounds. Objective: : To investigate hypoxia and chitosan-shelled OLN effects on MMP/TIMP production by human keratinocytes. Methods: : HaCaT cells were treated for 24 h with 10% v/v OLNs both in normoxia or hypoxia. Cytotoxicity and cell viability were measured through biochemical assays; cellular uptake by confocal microscopy; and MMP and TIMP production by enzyme-linked immunosorbent assay or gelatin zymography. Results: : Normoxic HaCaT cells constitutively released MMP-2, MMP-9, TIMP-1 and TIMP-2. Hypoxia strongly impaired MMP/TIMP balances by reducing MMP-2, MMP-9, and TIMP-2, without affecting TIMP-1 release. After cellular uptake by keratinocytes, nontoxic OLNs abrogated all hypoxia effects on MMP/TIMP secretion, restoring physiological balances. OLN abilities were specifically dependent on time-sustained oxygen diffusion from OLN core. Conclusion: : Chitosan-shelled OLNs effectively counteract hypoxia-dependent dysregulation of MMP/TIMP balances in human keratinocytes. Therefore, topical administration of exogenous oxygen, properly encapsulated in nanodroplet formulations, might be a promising adjuvant approach to promote healing processes in hypoxic wounds. - Highlights: • Hypoxia impairs MMP9/TIMP1 and MMP2/TIMP2 balances in HaCaT human keratinocytes. • Chitosan-shelled oxygen-loaded nanodroplets (OLNs) are internalised by HaCaT cells. • OLNs are not toxic to HaCaT cells. • OLNs effectively counteract hypoxia effects on MMP/TIMP balances in HaCaT cells. • OLNs appear as promising and cost-effective therapeutic tools for hypoxic wounds.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Simulation of Coupled Processes of Flow, Transport, and Storage of CO 2 in Saline Aquifers

This report is the final scientific one for the award DE- FE0000988 entitled “Simulation of Coupled Processes of Flow, Transport, and Storage of CO 2 in Saline Aquifers.” The work has been divided into six tasks. In task, “Development of a Three-Phase Non-Isothermal CO 2 Flow Module,” we developed a fluid property module for brine-CO 2 mixtures designed to handle all possible phase combinations of aqueous phase, sub-critical liquid and gaseous CO 2 , supercritical CO 2 , and solid salt. The thermodynamic and thermophysical properties of brine-CO 2 mixtures (density, viscosity, and specific enthalpy of fluid phases; partitioning of mass components among the different phases) use the same correlations as an earlier fluid property module that does not distinguish between gaseous and liquid CO 2 -rich phases. We verified the fluid property module using two leakage scenarios, one that involves CO 2 migration up a blind fault and subsequent accumulation in a secondary “parasitic” reservoir at shallower depth, and another investigating leakage of CO 2 from a deep storage reservoir along a vertical fault zone. In task, “Development of a Rock Mechanical Module,” we developed a massively parallel reservoir simulator for modeling THM processes in porous media brine aquifers. We derived, from the fundamental equations describing deformation of porous elastic media, a momentum conservation equation relating mean stress, pressure, and temperature, and incorporated it alongside the mass and energy conservation equations from the TOUGH2 formulation, the starting point for the simulator. In addition, rock properties, namely permeability and porosity, are functions of effective stress and other variables that are obtained from the literature. We verified the simulator formulation and numerical implementation using analytical solutions and example problems from the literature. For the former, we matched a one-dimensional consolidation problem and a two-dimensional simulation of the Mandel-Cryer effect. For the latter, we obtained a good match of temperature and gas saturation profiles, and surface uplift, after injection of hot fluid into a model of a caldera structure. In task, “Incorporation of Geochemical Reactions of Selected Important Species,” we developed a novel mathematical model of THMC processes in porous and fractured saline aquifers, simulating geo-chemical reactions associated with CO 2 sequestration in saline aquifers. Two computational frameworks, sequentially coupled and fully coupled, were used to simulate the reactions and transport. We verified capabilities of the THMC model to treat complex THMC processes during CO 2 sequestration by analytical solutions and we constructed reactive transport models to analyze the THMC process quantitatively. Three of these are 1D reactive transport under chemical equilibrium, a batch reaction model with equilibrium chemical reactions, and a THMC model with CO 2 dissolution. In task “Study of Instability in CO 2 Dissolution-Diffusion-Convection Processes,” We reviewed literature related to the study of density driven convective flows and on the instability of CO 2 dissolution-diffusion-convection processes. We ran simulations that model the density-driven flow instability that would occur during CO 2 sequestration. CO 2 diffused through the top of the system and dissolved in the aqueous phase there, increasing its density. Density fingers formed along the top boundary, and coalesced into a few prominent ones, causing convective flow that forced the fluid to the system bottom. These simulations were in two and three dimensions. We ran additional simulations of convective mixing with density contrast caused by variable dissolved CO 2 concentration in saline water, modeled after laboratory experiments in which supercritical CO2 was circulated in the headspace above a brine saturated packed sand in a pressure vessel. As CO 2 dissolved into the upper part of the saturated sand, liquid phase density increases causing instability and setting off convective mixing. We obtained good agreement with the laboratory experiments, which were characterized by finger development and associated mixing of dissolved CO 2 into the system. We then varied a wide range of parameters and conceptual models in order to analyze the possibility of convective mixing under different conditions, such as various boundary conditions, and chemical reaction conditions. The CO 2 fingers from different simulations showed great differences as time progressed, caused by permeability heterogeneity. The early time diffusive phenomenon was captured by fine grid resolution, and the permeability heterogeneity affected the pattern of the CO 2 fingers. In addition, the fingers from three-dimensional simulations tended to be larger and flatter than the two-dimensional ones. In task “Implementation of Efficient Parallel Computing Technologies,” we made enhancements and modifications to our code in order to substantially increase the grid size that could be run. We installed and ran it on various platforms, including a multi-core PC and a cluster, and verified the numerical implementation and parallel code using an example problem from the literature. This problem, with a grid size of sixty million, utilized the cluster’s entire memory and all of its processors. In task “Implementation of General Fracture Conceptual Models,” we used the MINC approach, a generalization of the double-porosity concept, to model flow through porous and fractured media. In this approach, flow within the matrix is described by subdividing the matrix into nested volumes, with flow occurring between adjacent nested matrix volumes as well as between the fractures and the outer matrix volume. We generalized Hooke’s law to a thermo-multi- poroelastic medium, and derived from the fundamental equations describing deformation of porous and fractured elastic media a momentum conservation equation for thermo-multi- poroelastic media. This equation is a generalization to multi-poroelastic media of the one derived in Task 3.0 for single porosity media. We describe two simulations to provide model verification and application examples. The first, one-dimensional consolidation of a double-porosity medium, is compared to an analytical solution. The second is a match of published results from the literature, a simulation of CO 2 injection into hypothetical aquifer-caprock systems.

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

Nonnegative methods for bilinear discontinuous differencing of the S N equations on quadrilaterals

Historically, matrix lumping and ad hoc flux fixups have been the only methods used to eliminate or suppress negative angular flux solutions associated with the unlumped bilinear discontinuous (UBLD) finite element spatial discretization of the two-dimensional S N equations. Though matrix lumping inhibits negative angular flux solutions of the S N equations, it does not guarantee strictly positive solutions. In this paper, we develop and define a strictly nonnegative, nonlinear, Petrov-Galerkin finite element method that fully preserves the bilinear discontinuous spatial moments of the transport equation. Additionally, we define two ad hoc fixups that maintain particle balance and explicitly set negative nodes of the UBLD finite element solution to zero but use different auxiliary equations to fully define their respective solutions. We assess the ability to inhibit negative angular flux solutions and the accuracy of every spatial discretization that we consider using a glancing void test problem with a discontinuous solution known to stress numerical methods. Though significantly more computationally intense, the nonlinear Petrov-Galerkin scheme results in a strictly nonnegative solution and is a more accurate solution than all the other methods considered. One fixup, based on shape preserving, results in a strictly nonnegative final solution but has increased numerical diffusion relative to the Petrov-Galerkin scheme and is less accurate than the UBLD solution. The second fixup, which preserves as many spatial moments as possible while setting negative values of the unlumped solution to zero, is less accurate than the Petrov-Galerkin scheme but is more accurate than the other fixup. However, it fails to guarantee a strictly nonnegative final solution. As a result, the fully lumped bilinear discontinuous finite element solution is the least accurate method, with significantly more numerical diffusion than the Petrov-Galerkin scheme and both fixups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2014

A Novel Algorithm for Solving the Multidimensional Neutron Transport Equation on Massively Parallel Architectures

We employ the Integral Transport Matrix Method (ITMM) as the kernel of new parallel solution methods for the discrete ordinates approximation of the within-group neutron transport equation. The ITMM abandons the repetitive mesh sweeps of the traditional source iterations (SI) scheme in favor of constructing stored operators that account for the direct coupling factors among all the cells' fluxes and between the cells' and boundary surfaces' fluxes. The main goals of this work are to develop the algorithms that construct these operators and employ them in the solution process, determine the most suitable way to parallelize the entire procedure, and evaluate the behavior and parallel performance of the developed methods with increasing number of processes, P. The fastest observed parallel solution method, Parallel Gauss-Seidel (PGS), was used in a weak scaling comparison with the PARTISN transport code, which uses the source iteration (SI) scheme parallelized with the Koch-baker-Alcouffe (KBA) method. Compared to the state-of-the-art SI-KBA with diffusion synthetic acceleration (DSA), this new method- even without acceleration/preconditioning-is completitive for optically thick problems as P is increased to the tens of thousands range. For the most optically thick cells tested, PGS reduced execution time by an approximate factor of three for problems with more than 130 million computational cells on P = 32,768. Moreover, the SI-DSA execution times's trend rises generally more steeply with increasing P than the PGS trend. Furthermore, the PGS method outperforms SI for the periodic heterogeneous layers (PHL) configuration problems. The PGS method outperforms SI and SI-DSA on as few as P = 16 for PHL problems and reduces execution time by a factor of ten or more for all problems considered with more than 2 million computational cells on P = 4.096.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Bowl breakout: Escaping the positive region when searching for saddle points

We present a scheme improving the minimum-mode following method for finding first order saddle points by confining the displacements of atoms to the subset of those subject to the largest force. By doing so it is ensured that the displacement remains of a local character within regions where all eigenvalues of the Hessian matrix are positive. However, as soon as a region is entered where an eigenvalue turns negative all atoms are released to maintain the ability of determining concerted moves. Applying the proposed scheme reduces the required number of force calls for the determination of connected saddle points by a factor two or more compared to a free search. Furthermore, a wider distribution of the relevant low barrier saddle points is obtained. Finally, the dependency on the initial distortion and the applied maximal step size is reduced making minimum-mode guided searches both more robust and applicable.

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

Advanced surface plasma nitriding for development of corrosion resistant and accident tolerant fuel cladding: Final Report

Although various surface coating techniques have been proposed to increase oxidation and corrosion resistance of fuel cladding materials, the de-bonding of the coating layer with the original cladding matrix under exposure to coolants makes such approaches unsuitable for reactor applications. Furthermore, the feasibility of techniques for large scale processing on cladding tubes remains another technological bottleneck. This project aims to develop a hollow cathode plasma nitriding technique to solve the above issues. Starting with alloys of DOE interest, the team will apply an advanced surface plasma nitriding technique to convert alloy surface layers into nitride layers for better structural integrity and compatibility with both coolants and nuclear fuels. The technique can form uniform nitride layers on both inner and outer surfaces of cylindrical hollow tubes. Both treated and untreated samples will be irradiated by using Fe self ions or dual beams of Fe ions and He ions. Various microstructural characterizations and mechanical property measurements will be performed. For outside surfaces, samples will be tested in high temperature water and liquid molten salt to evaluate their corrosion resistance and compatibility with coolants for LWR and fast reactor applications. For inside surfaces, solid-to-solid diffusion couple experiments will be used to study interface reactions of nitrided samples with uranium, to test whether a nitrided layer can slow down fuel-cladding interactions. The project will impact both the development of advanced methods for manufacturing and the development of advanced reactor in-core structural materials. The proposed technique can form nitride layers on alloys of irregular shapes without an edging effect - an issue encountered in traditional plasma nitriding processes, and is suitable for surface processing of other in-core components.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Behavior of helium cavities in ion-irradiated W-Ni-Fe ductile-phase toughened tungsten

This study reports on the distribution of helium (He) cavities in a hot-rolled W-Ni-Fe ductile-phase toughened tungsten (DPT W) composite irradiated to a dose and a helium concentration that are comparable to those in the material after 5-year irradiation in a conceptual fusion power plant. Here, the DPT W sample consists of W particles embedded in a ductile-phase NiFeW matrix with a nominal composition of 90W-7Ni-3Fe by weight. It was hot-rolled to a thickness reduction by 87% (87R DPT W). Sequential irradiation of the material with 1.2 MeV Ni + ions to a fluence of 2.15×10 16 Ni + /cm 2 and 90 keV He + ions to 6.5×10 15 He + /cm 2 was performed at 973 K. Larger He cavities with a lower number density are observed in NiFeW than W. Helium cavities are aggregated preferentially along the NiFeW/W interphase boundary. This behavior is not observed along the W/W grain boundary under the same irradiation conditions. A similar distribution of He cavities with comparable sizes is observed in 87R DPT W irradiated with He + ions only, suggesting that the He cavity distribution is not significantly affected by pre-irradiation of Ni + ions at 973 K. The data from this study does not provide any evidence for formation of visible voids or Ni precipitates in W phase irradiated with Ni + ions only. Diffusion and clustering of He atoms and mobile vacancies during ion irradiation at 973 K are responsible for the formation and growth of the He cavities.

36 MATERIALS SCIENCE↗