u ltra-High Temperature SiC-matrix Solar Thermal Air Receiver (HOTSSTAR), Enabled by Additive Manufacturing
Progress report on the design and initial on-sun test results of HOTSSTAR
SEARCH · Engineering Papers
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.
Progress report on the design and initial on-sun test results of HOTSSTAR
Poster presentation on HOTSSTAR concept, design and initial on-sun test and evaluation of SiC air receiver design.
The U.S. Department of Energy Office of Nuclear Energy (DOE NE) and the Idaho National Laboratory (INL) Advanced Reactor Technologies (ART) Advanced Gas Reactor (AGR) Fuel Development and Qualification program (referred to as AGR Fuel program hereafter) are pursuing qualification of tristructural isotropic (TRISO) coated particle fuel for use in high temperature gas cooled reactors (HTGRs). The AGR Fuel program was established to provide a fuel qualification data set in support of the licensing and operation of an HTGR. BWX Technologies Nuclear Operations Group (BWXT-NOG) was subcontracted to fabricate the fuel for the AGR program. Several investments and innovations were realized in preparation to fabricate fuel for the AGR-5/6/7 irradiation experiments that brought fuel fabrication fully out of the laboratory and into engineering-scale operations. These included: • Increased the kernel fabrication line capacity and uniformity • Upgraded ancillary support equipment and processes for the tristructural isotropic (TRISO) coating furnace • Demonstrated efficient production of the matrix precursor powder by dry jet milling of co mingled components • Demonstrated an engineering-scale method for quick and efficient overcoating TRISO particles with the matrix precursor • Demonstrated an automated, multi cavity compacting system with a volumetric feed system • Demonstrated a combined-cycle thermal treatment furnace These changes increased production rates of some of these processes by an order of magnitude or more while eliminating the use of flammable solvents, multiple grinding and sorting operations, and the weighing out of individual die charges. Three fuel kernel lots were fabricated for production of the fuel for AGR-5/6/7. The initial lot (J52R-16-39316) was certified to fuel specifications but was not used because the kernels had a high fraction of internal fissures that caused an unacceptable fraction of the kernels to fragment when charged to the coating furnace where the TRISO coating would be deposited. Fragmented kernels increased the dispersed uranium in the particles and produced a worrisome fraction of dimpled particles with an elevated probability of in-pile failure. After some efforts to identify the cause of the fissure formation, two additional lots were produced with much lower fissure fractions, J52R-16-69317 and 69318. The latter kernel lot was a backup to the first and was not needed. Multiple kernel batches were composited to form each of the lots so as to simulate a commercial-scale operation where kernel batches would also be composited. Multiple TRISO coating runs were performed and the product characterized so that several could be composited into a TRISO lot. TRISO lot J52R-16-98005 conformed to all fuel specifications except the mean outer pyrocarbon (OPyC) layer thickness was thinner than specified. Furthermore, it was determined that the TRISO lot had a dispersed uranium fraction (DUF) that might result in the compacts not meeting the DUF specification. A review of the role of the OPyC layer and consequences of the DUF by the Technical Coordination Team and INL concluded that the fuel was acceptable for use in the AGR-5/6/7 irradiation experiment. The TRISO particles were overcoated with the matrix precursor that had been produced in a jet-milling operation. The overcoating was performed in equipment originally designed to coat pharmaceuticals. The overcoating process performed well; producing highly spherical and uniform overcoats requiring little upgrading and no recycle or rework. TRISO particles were overcoated with the matrix precursor to achieve nominal volumetric packing fractions (PFs) of TRISO particles of 25% and 40% for the irradiation experiments. The 40% PF compacts occupy the first and fifth test capsule in the test train while the inner three capsules are loaded with 25% PF compacts. The resinated graphite matrix precursor powder was a derivative of the German A3-27 matrix formulation, which differs from previous AGR irradiation campaigns that used an A3-3 formulation. Jet milling of the matrix powder precursor produced a finer mean graphite particle size than the milling operations used for the A3-3 matrix powder precursor. Changes made in the matrix formula and equipment yielded compacts with significantly higher matrix density than was attained in previous AGR irradiation campaigns. The changes in the matrix formulation and the means of milling the powders also complicated resolution of the three fuel compact defect fractions, DUF, exposed kernel fraction (EKF), and the silicon carbide defect fraction (SDF). Characterization data from BWXT-NOG had some anomalous results, so samples of the fuel compacts and overcoated TRISO particles were also analyzed by Oak Ridge National Laboratory (ORNL) to ensure that the defect fractions were accurately characterized.
A natural uranium-based, unirradiated tristructural isotropic (TRISO) fuel compact was nondestructively imaged using both X-ray (XCT) and neutron computed tomography (nCT). While XCT of compacts can provide information on fuel kernels, imaging artifacts preclude examination of the graphite matrix. In this work, nCT was used for the first time on a TRISO compact to examine the graphite matrix. A crack was clearly resolved within the graphite matrix, proving that nCT is a viable tool for nondestructive volumetric examination of the matrix material in TRISO fuel compacts. The XCT and nCT data were then fused together to create a more comprehensive dataset containing both matrix and fuel kernels.
We generalize the self-consistent cooperative hopping theory for a dilute spherical penetrant or tracer activated dynamics in dense metastable hard sphere fluids and glasses to address the effect of external stress, the consequences of which are systematically established as a function of matrix packing fraction and penetrant-to-matrix size ratio. All relaxation processes speed up under stress, but the difference between the penetrant and matrix hopping (alpha relaxation) times decreases significantly with stress corresponding to less time scale decoupling. A dynamic crossover occurs at a critical “slaving onset” stress beyond which the matrix activated hopping relaxation time controls the penetrant hopping time. This characteristic stress increases (decreases) exponentially with packing fraction (size ratio) and can be well below the absolute yield stress of the matrix. Below the slaving onset, the penetrant hopping time is predicted to vary exponentially with stress, differing from the power law dependence of the pure matrix alpha time due to system-specificity of the stress-induced changes in the penetrant local cage and elastic barriers. An exponential growth of the penetrant alpha relaxation time with size ratio under stress is predicted, and at a fixed matrix packing fraction, the exponential relation between penetrant hopping time and stress for different size ratios can be collapsed onto a master curve. Direct connections between the short- and long-time activated penetrant dynamics and between the penetrant (or matrix) alpha relaxation time and matrix thermodynamic dimensionless compressibility are also predicted. The presented results should be testable in future experiments and simulations.
High strength and ductility are highly desired in fiber-reinforced composites, yet achieving both simultaneously remains elusive. A hierarchical architecture is developed utilizing high aspect ratio chemically transformable thermoplastic nanofibers that form covalent bonding with the matrix to toughen the fiber-matrix interphase. The nanoscale fibers are electrospun on the micrometer-scale reinforcing carbon fiber, creating a physically intertwined, randomly oriented scaffold. Unlike conventional covalent bonding of matrix molecules with reinforcing fibers, here, the nanofiber scaffold is utilized - interacting non-covalently with core fiber but bridging covalently with polymer matrix - to create a high volume fraction of immobilized matrix or interphase around core reinforcing elements. This mechanism enables efficient fiber-matrix stress transfer and enhances composite toughness. Molecular dynamics simulation reveals enhancement of the fiber-matrix adhesion facilitated by nanofiber-aided hierarchical bonding with the matrix. The elastic modulus contours of interphase regions obtained from atomic force microscopy clearly indicate the formation of stiffer interphase. These nanoengineered composites exhibit a ≈60% and ≈100% improved in-plane shear strength and toughness, respectively. This approach opens a new avenue for manufacturing toughened high-performance composites.
This work presents a new formulation of the axial expansion transport method explicitly using Legendre polynomials for arbitrarily high-order expansions. This new formulation also features an alternative method of axial leakage calculation to allow for nonextruded flat source region meshes. This alternative axial leakage is introduced alongside a balance equation requirement to ensure that neutron balance is preserved in the coarse mesh for a given axial leakage formulation, which allows for effective coarse mesh finite difference acceleration. A matrix exponential table method is derived to allow for fast computations of arbitrarily high-order matrix exponentials for this work and precludes the need for further research into matrix exponential calculations for this method. Numerical results are presented that demonstrate the stability of the axial expansion method in systems with voidlike regions, showcase the speedup from matrix exponential tables, and investigate the axial convergence of the method in terms of both expansion order and mesh size.
The position operator 𝑟̂ appears as 𝑖∂ 𝑝 in wave mechanics, while its matrix form (e.g., under a Bloch basis) is well known diverging in diagonals, causing difficulties in basis transformation, observable yielding, etc. We aim to find a convergent r-matrix (CRM) to improve the existing divergent r-matrix (DRM), and investigate its influence at both the conceptual and the application levels. A key modification is increasing the familiar substitution of 𝑟̂ by 𝑖∂ 𝑝 to 𝑖∑ 𝑗 ∂ 𝑘 𝑗 , namely the N-th Weyl algebra. Resolving the divergence makes r-matrix rigorously defined, and we are able to show r-matrix is distinct from a spin matrix in terms of its defining principles, transformation behavior, and the observable it yields. Conceptually, the CRM fills the logical gap between the r-matrix and the Berry connection (this unremarked vagueness has caused the diagonal divergence). In application, we focus on transport, and discover that the Hermitian matrix is not identical with the associative Hermitian operator, i.e., 𝑟 𝑚,𝑛 = 𝑟$^∗_{𝑛,𝑚}$ ⇎ 𝑟̂ = 𝑟̂ † , which subtly affects the celebrated Berry curvature formula for adiabatic current. We also discuss how such a non-representation CRM can contribute to building a unified transport theory.
Partitioned Global Address Space (PGAS) models exhibit tremendous promise in developing efficient and productive distributed-memory parallel applications. They have been used extensively in scientific computations due to conveniently offering a ``shared-memory''-like model and convenient interfaces that separate communication with synchronization. Traditionally, PGAS communication models have been applied to dense/contiguously distributed data, but most modern applications depict varied levels of sparsity. Existing PGAS models require certain adaptations to support distributed sparse computations, since associated computations often require matrix arithmetic, in addition to data movement. The Global Arrays toolkit from Pacific Northwest National Laboratory (PNNL) is one of the earliest PGAS models to combine one-sided data communication and distributed matrix operations and is still used in the popular NWChem quantum chemistry suite. Recently, we have expanded the Global Arrays toolkit to support common sparse operations, like sparse matrix-dense matrix multiplies (SpMM), sparse matrix-sparse matrix multiplication (SpGEMM) and Sampled Dense-Dense Matrix Multiplication (SDDMM). As it turns out, these operations are the bedrock of sparse Deep Learning (DL); sparse deep neural networks and Graph Neural Networks (GNNs) have gained increasing attention recently in achieving speedups on training and inference with reduced memory footprints. Unlike scientific applications in High Performance Computing (HPC), modern (distributed-memory capable) DL toolkits often rely on non-standardized and closed-source vendor software optimizations, creating challenges in software-hardware co-design at scale. Our goal is to support a variety of distributed-memory sparse matrix operations and helper functions in the newly created Sparse Global Arrays (SGA), such that it is possible to build portable and productive Machine Learning scenarios for algorithm/software and hardware codesign purposes. Contemporary data-parallel schemes for training/inference are undergoing a major overhaul since model replication limits scalability and causes resource inefficiencies. As such, we have adopted tensor parallelism in decomposing the model and inputs, to mitigate memory issues. Current implementation is built on top of MPI and uses CPUs to maximize the portability across the platforms.
Lipid membranes, the primary matrix of cell membranes, exhibit a hierarchy of dynamics from molecular motions to collective undulations, which control spatiotemporal membrane phenomena and regulate various cellular functions. Thus, understanding how various dynamic modes influence membrane behavior is crucial to uncovering the design rules of cell membranes and implementing them in future technologies. Previous membrane dynamics studies using neutron spin-echo spectroscopy (NSE) have shed light on the dependence of collective membrane fluctuations on viscoelastic membrane properties, demonstrating that thickness fluctuation modes (100ns) are dictated by membrane viscosity, which directly governs fast (ps) molecular diffusive motions. Here, we focus on phase-separating lipid membranes, composed of DMPC:DSPC mixtures, to understand the effect of domain formation/growth on hierarchical matrix dynamics using selective lipid deuteration and different neutron spectroscopy techniques. NSE measurements on DMPC:DSPC-d83 (70:30 mol%) show that bending rigidity of the DMPC-rich matrix graduallyincreases with decreasing temperature - as the fully-fluid membrane approaches the upper phase transition (i.e. with the formation of transient DSPC clusters) - and increases further in the phase coexistence region with the growth of DSPC domains. Comparatively, NSE measurements on DMPC-d54:DSPC-d83 vesicles show pronounced matrix thickness fluctuations in the fully-fluid phase and drastic suppression of this fluctuation mode in the gel-fluid coexistence phase, even when the DMPC-rich matrix is still in its fluid phase. This suppression in fluctuations suggests an increase in the matrix viscosity with the growth of gel DSPC domains, which is further investigated by selectively measuring lipid diffusion in the DMPC-rich matrix using neutron backscattering spectroscopy. Comparisons of the experimental results with coarse-grained MD simulationsillustrate the underlying coupling mechanisms. Overall, these observations indicate a clear mechanical and viscoelastic coupling betweenlipid domains and their host matrix, suggesting that domain formation influences membrane properties beyond local effects.
Hybrid materials of metal chalcogenide colloidal quantum dots (QDs) embedded in metal halide perovskites (MHPs) have led to composites with synergistic properties. Here, we investigate how QD size, surface chemistry, and MHP film formation methods affect the resulting optoelectronic properties of QD/MHP “dot-in-matrix” systems. We monitor the QD absorption and photoluminescence throughout synthesis, ligand exchange, and transfer into the MHP ink, and we characterize the final QD/MHP films via electron microscopy and transient absorption. In addition, we are the first to globally map how PbS QDs are distributed on the micrometer scale within these dot-in-matrix systems, using three-dimensional (3D) tomography time-of-flight secondary ion mass spectrometry. The surface chemistry imparted during synthesis directly affects the optical properties of the dot-in-matrix composites. Pb-halide passivation leads to QD/MHP dot-in-matrix samples with optical properties that are well-described by a theoretical model, based on a Type I finite-barrier heterostructure between the PbS QD and the MHP matrix. Samples without Pb-halide passivation show complicated size-dependent behavior, indicating a transition from a Type I heterostructure between the PbS QD wells and MHP barriers for small-sized QDs to PbS QDs that are electronically decoupled from the MHP matrix for larger QDs. Furthermore, the choice in perovskite antisolvent crystallization method leads to a difference in the spatial QD distribution within the perovskite matrix, differences in carrier lifetime, and photoluminescence shifts of up to 180 meV for PbS in methylammonium lead iodide. This work establishes an understanding of such emerging synergistic systems relevant for technologies such as photovoltaics, infrared emitters and detectors, and other unexplored technological applications.
This study investigates the role of discontinuous reinforcement on the microcracking of a preceramic polymer matrix during polymer to ceramic conversion. The material system comprised a carbosilane-based, preceramic polymer reinforced with mullite particles. The preceramic resin slurry was formulated for photopolymerization on a digital light projection 3D printer. Micro X-ray computed tomography, using a synchrotron light source, monitored a printed composite part during pyrolysis. The conversion of the carbosilane matrix into Si(O)C generated microcracks in the matrix that radially extended from particles. The likelihood of microcrack formation positively correlated with particle size. The largest particles (>10 4 µm 3 volume or >60 µm side length) always abutted microcracks, while the matrix surrounding smaller particles (<5 µm) was free of microcracks. Numerical calculations showed the particles resist the shrinking matrix during its conversion. This created tensile stresses within the matrix. The driving force for microcrack growth was also analyzed with regard to the influence of matrix material parameters, particle form-factor and crack length. These results reveal the need to consider the form factor of discontinuous reinforcements, with regard to the material properties of the polymer-derived ceramic, in order to reduce or eliminate defects in the final part.
We illustrate a technique to calculate the gradient of scalar functions that are defined on any arbitrary matrix subspace. It generalizes our earlier work titled “What is the gradient of a scalar function of a symmetric matrix ?”(Indian Journal of Pure and Applied Mathematics (2022), https://doi.org/10.1007/s13226-022-00313-x), in which we considered the special case of the subspace of symmetric matrices. Extant methods to calculate the gradient in such cases have an inherent flaw which leads to spurious results that populate several publications, as well as respected textbooks and handbooks on matrix calculus. Here, we examine these sources and results in a rigorous and concrete mathematical setting of a finite-dimensional inner-product space and discover the inherent flaw and also a remedy. We demonstrate two ways to calculate the derivative/gradient and second derivative for scalar functions of matrices defined over an arbitrary matrix subspace; the first method is by considering any (differentiable) extension to the space of square matrices and projection of its gradient onto the given subspace. The second method utilizes an ordered basis and computes each component of the gradient through evaluation of the directional derivative. All the ideas presented are illustrated by non-trivial examples, namely, considering the subspace of 3 x 3 circulant and Toeplitz matrices and presenting the results of gradient-descent with both the spurious and correct gradients. Moreover, our bibliography makes it clear that a rigorous approach to matrix calculus is not common in practice, and our presentation of matrix calculus in the language of inner-product spaces will be significant and meaningful for applied mathematicians, engineers and researchers working in inter-disciplinary fields to avoid the conceptual pitfalls that exist.
We illustrate a technique to calculate the gradient of scalar functions that are defined on any arbitrary matrix subspace. It generalizes our earlier work titled “What is the gradient of a scalar function of a symmetric matrix ?”(Indian Journal of Pure and Applied Mathematics (2022), doi:10.1007/s13226-022-00313-x), in which we considered the special case of the subspace of symmetric matrices. Extant methods to calculate the gradient in such cases have an inherent flaw which leads to spurious results that populate several publications, as well as respected textbooks and handbooks on matrix calculus. We examine these sources and results in a rigorous and concrete mathematical setting of a finite-dimensional inner-product space and discover the inherent flaw and also a remedy. We demonstrate two ways to calculate the derivative/gradient and second derivative for scalar functions of matrices defined over an arbitrary matrix subspace; the first method is by considering any (differentiable) extension to the space of square matrices and projection of its gradient onto the given subspace. The second method utilizes an ordered basis and computes each component of the gradient through evaluation of the directional derivative. All the ideas presented are illustrated by non-trivial examples, namely, considering the subspace of 3 × 3 circulant and Toeplitz matrices and presenting the results of gradient-descent with both the spurious and correct gradients. Moreover, our bibliography makes it clear that a rigorous approach to matrix calculus is not common in practice, and our presentation of matrix calculus in the language of inner-product spaces will be significant and meaningful for applied mathematicians, engineers and researchers working in inter-disciplinary fields to avoid the conceptual pitfalls that exist.
Numerical simulations using the staggered-grid hydrodynamics (SGH) discretization suffer from hourglass instabilities. In this work, we develop a stabilization method to suppress the hourglass instabilities using techniques from the virtual element method (VEM). The stiffness matrix of the VEM consists of two terms: the consistency matrix which is rank deficient and the stability matrix. Here, we first show that in two dimensions and on general polygons, the stiffness matrix of the SGH is identical to the consistency matrix of the linear VEM for both the diffusion equation and the linear elasticity equation. These analyses explain the origin of the hourglass instabilities of the SGH discretization method, and establish a theoretical foundation for our proposed stabilization method by augmenting the stiffness matrix of the SGH discretization using the VEM stability matrix. Then, we present numerical examples using Lagrangian SGH simulations. The numerical experiments demonstrate that the proposed VEM stabilization method is effective at eliminating hourglass modes in the SGH discretization.
Geopolymers, as a potentially environmentally friendly alternative to Portland cement, are increasingly attracting attention in the construction industry. Various methods have been applied for customizing the properties of geopolymers and improving their commercial viability. One of the promising methods for refining the properties of geopolymers such as their toughness is the use of short fibers. The effectiveness of a high-strength short fiber in the geopolymer matrix is largely dependent on the interfacial bonding between the fiber and its surrounding matrix. While the importance of this interfacial chemistry is highlighted in the literature, the characteristics of this bonding structure have not been fully understood. In this paper, we aim to investigate the bonding mechanism between the carbon fiber and metakaolin-based geopolymer matrix. For the first time, the existence and nature of the chemical bonding at the interfacial region (interphase) between carbon fiber and geopolymer matrix has been revealed. X-ray pair distribution function computed tomography (PDF-CT), field emission-scanning electron microscopy imaging, and nanoindentation techniques are employed to discern the chemo-mechanical properties of the interphase. PDF-CT results show the emergence of a new atom-atom correlation at the interfacial region (around 1.82 angstrom). This correlation is a characteristic of interfacial bonding between the fiber and its surrounding matrix, where the existence of chemical linkages (potentially Al-V-O-C) between fibers and the matrix contributes to the adhesion between the two constituents making up the composite. Due to such chemical bonding, the nanomechanical properties of the interfacial region fall between that of the carbon fiber and geopolymer. The combination of advanced techniques is proved useful for enhancing our understanding of the interfacial chemistry between fibers and the binding matrix. This level of knowledge facilitates the engineering of composite systems through the manipulation of their nanostructure.
This report describes the development of a comprehensive catalogue of generic features, events, and processes (FEPs) that are potentially important for the post-closure performance of a repository for high-level radioactive waste (HLW) and spent nuclear fuel (SNF) in salt (halite) host rock. The FEPs and other supporting information have been entered into a “SaltFEP” Database. The generic salt repository FEPs include consideration of relevant FEPs from a number of U.S., Dutch, German, and international FEP lists and should be a suitable starting point for any repository program in salt host rock. The salt FEP catalogue and database employ a FEP classification matrix approach that is based on the concept that a FEP is typically a process or event acting upon or within a feature. The FEP matrix provides a two-dimensional structure consisting of a Features/Components axis that defines the “rows” and a Processes/Events axis that defines the “columns” of the matrix. The design of the FEP classification matrix is consistent with repository performance assessment – the Features/Components axis is organized vertically to generally correspond to the direction of potential radionuclide migration (from the waste to the biosphere) and the Processes/Events axis is designed to represent the common two-way couplings between thermal processes and other processes (such as thermal-mechanical or thermal-hydrologic processes). Related FEPs can be easily identified – related FEPs will typically be grouped in a single matrix cell or aligned along a common row (Feature/Component) or column (Process/Event). The online SaltFEP database can be downloaded from www.saltfep.org. It contains the FEP matrix, the FEPs, and the associated processes for each FEP. It provides a starting point to create and document site-specific individual FEPs. Furthermore, the FEP matrix is connected to the Salt Knowledge Archive (SKA), a database of about 20,000 references and documents representing the historical knowledge on radioactive disposal in salt. This work is the result of an ongoing collaboration between researchers in the U.S., the Netherlands, and Germany, and supports the NEA Salt Club Mandate. It builds upon prior work which is documented.
A thermal neutron absorber material composed of Al 3 Hf particles in an aluminum matrix is under development for the Advanced Test Reactor. This metal matrix composite was fabricated via hot pressing of high-purity aluminum and micrometer-size Al 3 Hf powders at volume fractions of 20.0, 28.4, and 36.5%. Room temperature tensile and hardness testing of unirradiated specimens revealed a linear relationship between volume fraction and strength, while the tensile data showed a strong decrease in elongation between the 20 and 36.5% volume fraction materials. Tensile tests conducted at 200 °C on unirradiated material revealed similar trends. Evaluations were then conducted on specimens irradiated at 66 to 75 °C to four dose levels ranging from approximately 1 to 4 dpa. Tensile properties exhibited the typical increase in strength and decrease in ductility with dose that are common for metallic materials irradiated at ≤0.4T m . Hardness also increased with neutron dose. The difference in strength between the three different volume fraction materials was roughly constant as the dose increased. Nanoindentation measurements of Al 3 Hf particles in the 28.4 vol% material showed the expected trend of increased hardness with irradiation dose. Transmission electron microscopy revealed oxygen at the interface between the Al 3 Hf particles and aluminum matrix in the irradiated material. Scanning electron microscopy of the exterior surface of tensile tested specimens revealed that deformation of the material occurs via plastic deformation of the Al matrix, cracking of the Al 3 Hf particles, and to a lesser extent, tearing of the matrix away from the particles. The fracture surface of an irradiated 28.4 vol% specimen showed failure by brittle fracture in the particles and ductile tearing of the aluminum matrix with no loss of cohesion between the particles and matrix. The coefficient of thermal expansion decreased upon irradiation, with a maximum change of –6.3% for the annealed irradiated 36.5 vol% specimen.