Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix multiplication”

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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Developing Platinum-Group-Metal-Free Catalysts for Oxygen Reduction Reaction in Acid: Beyond the Single Metal Site

This project is to develop M (x) -N-C catalysts with dense multiple metal center (MMC) sites to meet the DOE 2025 activity target of 0.044 mA/cm 2 at 0.9 V IR-free ., as well as other goals such as durability. We have made important contributions to both catalyst development and fundamental understandings of the active sites in M-N-C catalysts in this project. We successfully made M (x) -N-C catalysts with some multiple metal center (MMC) sites by combining ionothermal carbonization with chemical vapor deposition (CVD). These catalysts, however, are not as active as the most active single-atom Fe-N-C catalysts made by the similar CVD process, owing likely to the low site density and the presence of inorganic Fe species such as iron carbides and nanoparticles. The most significant accomplishments we achieved in this project are: (1) we unraveled the formation pathway of Fe-N 4 sites during pyrolysis step-by-step and identified the trans-metalation mechanism, in collaboration with Deborah Myers and her colleagues at Argonne National Laboratory (ANL); (2) inspired by this finding, we pioneered the CVD synthesis of M-N-C catalysts (M = Mn, Fe, and Co), in which the Fe-N-C catalyst by CVD demonstrated an ORR activity of 33 mA/cm 2 at 0.9 V in H 2 -O 2 proton exchange membrane fuel cells (PEMFCs), very close to the ultimate goal of 35 mA/cm 2 of our project. This catalyst is the first Fe-N-C catalyst that contains only D1 sites without the D2 sites; whereas D1 and D2 sites have been always identified by Mossbauer in previous Fe-N-C catalysts. This finding helps to understand what the D1 and D2 sites are and their roles in catalyzing the ORR. (3) by improving the mass transport of the carbon matrix prior to the CVD process, the revised Fe-N-C catalyst made by CVD delivered a maximum power density of 0.53 W/cm 2 in H 2 -air PEMFCs. The improvement strategy was partly inspired by the computational modeling work by Adam Weber from LBNL, the Co-PI of this project, by developing, coding, and exercising a continuum level model of transport phenomena within a PGM-free catalyst layer. The model demonstrated that local resistances combined with limited site density of the PGM-free catalyst can result in limiting currents and poor polarization performance. The model also gave design guidance for impact of ECSA and overall catalyst-layer thickness. However, both Fe-N-C and Co-N-C catalysts developed by CVD showed poor durability in PEMFCS, in comparison with the traditional M-N-C catalysts synthesized via regular pyrolysis process. Consequently, we did not achieve the proposed durability targets. Despite so, we believe the FeNC-CVD catalysts with the poor durability and D1 sites only provides an excellent platform to understand the degradation mechanism of Fe-N-C in PEMFCs, the most important challenge in the development of M-N-C catalysts.

08 HYDROGEN↗
NASA NTRS1997

Determination of Residual Stress in Composite Materials Using Ultrasonic Waves

The performance of high temperature composites can be significantly affected by the presence of residual stresses. These stresses arise during cooling processes from fabrication to room temperature due to mismatch of thermal expansion coefficients between matrix and fiber materials. This effect is especially pronounced in metal matrix and intermetallic composites. It can lead to plastic deformations, matrix cracking and fiber/matrix interface debonding. In this work the feasibility of ultrasonic techniques for residual stress assessment in composites is addressed. A novel technique for absolute stress determination in orthotropic materials from angular dependencies of ultrasonic velocities is described. The technique is applicable for determination of both applied and residual stresses and does not require calibration measurements on a reference sample. The important advantage of this method is that stress is determined simultaneously with stress-dependent elastic constants and is thus decoupled from the material texture. It is demonstrated that when the principal plane stress directions coincide with acoustical axes, the angular velocity data in the plane perpendicular to the stress plane may be used to determine both stress components. When the stress is off the acoustical axes, the shear and the difference of the normal stress components may be determined from the angular dependence of group velocities in the plane of stresses. Synthetic sets of experimental data corresponding to materials with different anisotropy and stress levels are used to check the applicability of the technique. The method is also verified experimentally. A high precision ultrasonic wave transmission technique is developed to measure angular dependence of ultrasonic velocities. Examples of stress determination from experimental velocity data are given. A method is presented for determination of velocities of ultrasonic waves propagating through the composite material with residual stresses. It is based on the generalized self-consistent multiple scattering model. Calculation results for longitudinal and shear ultrasonic wave velocities propagating perpendicular to the fibers direction in SCS-6/Ti composite with and without residual stresses are presented. They show that velocity changes due to presence of stresses are of order 1%.

Rokhlin, S. I.↗
DOE OSTI2024

FY24 Progress Report: Disposition of Cracks & Features Observed in the Inner Can Closure Weld Region (ICCWR) of the 3013 Package

The long-term integrity of the 3013 containers is of interest for the safe storage of Pu materials. Although the 3013 standard embodies multiple barrier concept, the integrity of the inner container is considered to be crucial since it protects the Safety Class outer container from its contents. The container is designed to withstand high pressures that could result from the complete radiolysis of the maximum water content permissible by the 3013 standard. Shelf-life studies and destructive examination have not found high gas pressures but have found that corrosive gases are generated. Pitting and stress corrosion cracking (SCC) are considered to be critical corrosion modes for the performance of the inner container and have been observed during destructive examinations. Considerable research has been conducted on the possibility of aqueous electrolytes condensing on the inner container that can promote SCC. These studies indicate the uncertainties surrounding the formation of corrosive environments and the corrosion behavior of container materials. In this initial report, a Bayesian network (BN) model is described that can consider the uncertainties and the causal connections between various factors influencing the corrosion modes of the inner container. The BN model is preliminary and provides an initial framework to identify the necessary information. The report also provides initial experimental results on the electrochemical behavior of stainless steels in anticipated condensed environments from gas phase migration of acidic gases. The experimental results are consistent with the corrosion model. Recommendation for further work on the BN model include assembling an expert group to provide input to the BN structure and quantification of the conditional probability matrix, experimental studies to characterize the microstructure of the container, electrochemical studies to identify critical potentials for localized corrosion and SCC, and crack growth rate studies

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Nuclear data uncertainty propagation applied to the versatile test reactor conceptual design

We report the Versatile Test Reactor (VTR) currently under development is a 300 MWth sodium-cooled fast reactor (SFR) fueled with ternary metal alloy fuel, which aims to accelerate the testing of advanced nuclear fuels, materials, instrumentation, and sensors in high flux environments that are necessary to license the next generation of advanced reactor concepts. To support the VTR design process, uncertainties associated with the nuclear data has been propagated through the reactor core neutronics calculation to global parameters of interest, such as the core multiplication factor, kinetic parameters, and various reactivity feedback coefficients, following the sensitivity based uncertainty propagation approach. By folding the sensitivity coefficients, separately computed by the generalized perturbation theory code PERSENT and Monte Carlo code Serpent 2, with the variance-covariance matrices from COMMARA-2.0, we obtain the reaction-wise, isotope-wise, and overall uncertainties for each response of interest due to nuclear data uncertainty. With Serpent 2, the statistical error of the uncertainty is obtained by propagating the statistical error of the sensitivity coefficients through the same process using a newly developed uncertainty propagation method. From both codes, the overall top uncertainty contributors are found to be the cross section of Fe-56 elastic scattering, Na-23 elastic scattering, and U 238 inelastic scattering. The large contributions of the Fe-56 elastic scattering cross sections to global parameters are due to its relatively large relative uncertainty of 5–10% in nuclear data and the large volume of Fe-containing reflector assemblies in the fairly compact VTR core design. Both codes agreed well for the overall uncertainty estimates of all responses of interest, except the delayed neutron fraction, prompt neutron generation time, and the coolant density feedback coefficient, where Serpent 2 yielded a much larger value than PERSENT due to the large statistical error of sensitivity coefficients. The calculated uncertainties are also compared to those associated with other SFR cores. Another outcome of this study is a variance-covariance matrix of reactivity coefficients, which can be used in the subsequent uncertainty propagation to the system level to investigate the impact of identified uncertainties on system responses in the safety analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2004

Using Redundancy To Reduce Errors in Magnetometer Readings

A method of reducing errors in noisy magnetic-field measurements involves exploitation of redundancy in the readings of multiple magnetometers in a cluster. By "redundancy"is meant that the readings are not entirely independent of each other because the relationships among the magnetic-field components that one seeks to measure are governed by the fundamental laws of electromagnetism as expressed by Maxwell's equations. Assuming that the magnetometers are located outside a magnetic material, that the magnetic field is steady or quasi-steady, and that there are no electric currents flowing in or near the magnetometers, the applicable Maxwell 's equations are delta x B = 0 and delta(raised dot) B = 0, where B is the magnetic-flux-density vector. By suitable algebraic manipulation, these equations can be shown to impose three independent constraints on the values of the components of B at the various magnetometer positions. In general, the problem of reducing the errors in noisy measurements is one of finding a set of corrected values that minimize an error function. In the present method, the error function is formulated as (1) the sum of squares of the differences between the corrected and noisy measurement values plus (2) a sum of three terms, each comprising the product of a Lagrange multiplier and one of the three constraints. The partial derivatives of the error function with respect to the corrected magnetic-field component values and the Lagrange multipliers are set equal to zero, leading to a set of equations that can be put into matrix.vector form. The matrix can be inverted to solve for a vector that comprises the corrected magnetic-field component values and the Lagrange multipliers.

Kulikov, Igor↗
NASA NTRS2011

Post-Fragmentation Whole Genome Amplification-Based Method

This innovation is derived from a proprietary amplification scheme that is based upon random fragmentation of the genome into a series of short, overlapping templates. The resulting shorter DNA strands (<400 bp) constitute a library of DNA fragments with defined 3 and 5 termini. Specific primers to these termini are then used to isothermally amplify this library into potentially unlimited quantities that can be used immediately for multiple downstream applications including gel eletrophoresis, quantitative polymerase chain reaction (QPCR), comparative genomic hybridization microarray, SNP analysis, and sequencing. The standard reaction can be performed with minimal hands-on time, and can produce amplified DNA in as little as three hours. Post-fragmentation whole genome amplification-based technology provides a robust and accurate method of amplifying femtogram levels of starting material into microgram yields with no detectable allele bias. The amplified DNA also facilitates the preservation of samples (spacecraft samples) by amplifying scarce amounts of template DNA into microgram concentrations in just a few hours. Based on further optimization of this technology, this could be a feasible technology to use in sample preservation for potential future sample return missions. The research and technology development described here can be pivotal in dealing with backward/forward biological contamination from planetary missions. Such efforts rely heavily on an increasing understanding of the burden and diversity of microorganisms present on spacecraft surfaces throughout assembly and testing. The development and implementation of these technologies could significantly improve the comprehensiveness and resolving power of spacecraft-associated microbial population censuses, and are important to the continued evolution and advancement of planetary protection capabilities. Current molecular procedures for assaying spacecraft-associated microbial burden and diversity have inherent sample loss issues at practically every step, particularly nucleic acid extraction. In engineering a molecular means of amplifying nucleic acids directly from single cells in their native state within the sample matrix, this innovation has circumvented entirely the need for DNA extraction regimes in the sample processing scheme.

Benardini, James↗
DOE OSTI2023

DEEPEN: Newberry Volcano MT and Gravity Data - 2022 and 2023 Acquisition and Processing

As part of DEEPEN (DE-risking Exploration of geothermal Plays in magmatic ENvironments), a 3D play fairway analysis (PFA) was conducted at Newberry Volcano in Central Oregon for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). For use in this PFA, combined full tensor broadband magnetotelluric (MT) and gravity data were acquired, processed, and inverted by Enthalpion Energy LLC (Enthalpion) with support from NREL staff. The data collection efforts took place from June 19th to July 24. Data were collected with the goal of gaining an improved understanding of the South Flank and the extent of the magma chamber. This GDR submission includes the raw data, single inversions, joint inversions, resolution matrices, and a report on these processes. More detailed information about the folder structure of the datasets included in this submission may be found on pages 94-97 of Magnetotelluric and Gravity Survey Report.pdf below.

15 GEOTHERMAL ENERGY↗
NASA NTRS2008

Update on Development of SiC Multi-Chip Power Modules

Progress has been made in a continuing effort to develop multi-chip power modules (SiC MCPMs). This effort at an earlier stage was reported in 'SiC Multi-Chip Power Modules as Power-System Building Blocks' (LEW-18008-1), NASA Tech Briefs, Vol. 31, No. 2 (February 2007), page 28. The following recapitulation of information from the cited prior article is prerequisite to a meaningful summary of the progress made since then: 1) SiC MCPMs are, more specifically, electronic power-supply modules containing multiple silicon carbide power integrated-circuit chips and silicon-on-insulator (SOI) control integrated-circuit chips. SiC MCPMs are being developed as building blocks of advanced expandable, reconfigurable, fault-tolerant power-supply systems. Exploiting the ability of SiC semiconductor devices to operate at temperatures, breakdown voltages, and current densities significantly greater than those of conventional Si devices, the designs of SiC MCPMs and of systems comprising multiple SiC MCPMs are expected to afford a greater degree of miniaturization through stacking of modules with reduced requirements for heat sinking; 2) The stacked SiC MCPMs in a given system can be electrically connected in series, parallel, or a series/parallel combination to increase the overall power-handling capability of the system. In addition to power connections, the modules have communication connections. The SOI controllers in the modules communicate with each other as nodes of a decentralized control network, in which no single controller exerts overall command of the system. Control functions effected via the network include synchronization of switching of power devices and rapid reconfiguration of power connections to enable the power system to continue to supply power to a load in the event of failure of one of the modules; and, 3) In addition to serving as building blocks of reliable power-supply systems, SiC MCPMs could be augmented with external control circuitry to make them perform additional power-handling functions as needed for specific applications. Because identical SiC MCPM building blocks could be utilized in such a variety of ways, the cost and difficulty of designing new, highly reliable power systems would be reduced considerably. This concludes the information from the cited prior article. The main activity since the previously reported stage of development was the design, fabrication, and testing a 120- VDC-to-28-VDC modular power-converter system composed of eight SiC MCPMs in a 4 (parallel)-by-2 (series) matrix configuration, with normally-off controllable power switches. The SiC MCPM power modules include closed-loop control subsystems and are capable of operating at high power density or high temperature. The system was tested under various configurations, load conditions, load-transient conditions, and failure-recovery conditions. Planned future work includes refinement of the demonstrated modular system concept and development of a new converter hardware topology that would enable sharing of currents without the need for communication among modules. Toward these ends, it is also planned to develop a new converter control algorithm that would provide for improved sharing of current and power under all conditions, and to implement advanced packaging concepts that would enable operation at higher power density.

Lostetter, Alexander↗
NASA NTRS2012

Analysis of Direct Samples of Early Solar System Aqueous Fluids

Over the past three decades we have become increasingly aware of the fundamental importance of water, and aqueous alteration, on primitive solar-system bodies. Some carbonaceous and ordinary chondrites have been altered by interactions with liquid water within the first 10 million years after formation of their parent asteroids. Millimeter to centimeter-sized aggregates of purple halite containing aqueous fluid inclusions were found in the matrix of two freshly-fallen brecciated H chondrite falls, Monahans (1998, hereafter simply "Monahans") (H5) and Zag (H3-6) (Zolensky et al., 1999; Whitby et al., 2000; Bogard et al., 2001) In order to understand origin and evolution of the aqueous fluids inside these inclusions we much measure the actual fluid composition, and also learn the O and H isotopic composition of the water. It has taken a decade for laboratory analytical techniques to catch up to these particular nanomole-sized aqueous samples. We have recently been successful in (1) measuring the isotopic composition of H and O in the water in a few fluid inclusions from the Zag and Monahans halite, (2) mineralogical characterization of the solid mineral phases associated with the aqueous fluids within the halite, and (3) the first minor element analyses of the fluid itself. A Cameca ims-1270 equipped with a cryo-sample-stage of Hokkaido University was specially prepared for the O and H isotopic measurements. The cryo-sample-stage (Techno. I. S. Corp.) was cooled down to c.a. -190 C using liquid nitrogen at which the aqueous fluid in inclusions was frozen. We excavated the salt crystal surfaces to expose the frozen fluids using a 15 keV Cs+ beam and measured negative secondary ions. The secondary ions from deep craters of approximately 10 m in depth emitted stably but the intensities changed gradually during measurement cycles because of shifting states of charge compensation, resulting in rather poor reproducibility of multiple measurements of standard fluid inclusions of +/- 90 0/00(2 sigma) for delta D, and +/- 29 0/00 (2 sigma) for delta O-18. On the other hand, the reproducibility of Delta O-17 is plus or minus 8 /00 (2 sigma ) because the observed variations of isotope ratios follow a mass dependent fractionation law. Variations of delta D of the aqueous fluids range over sog,a 330(90; 2 sigma ) to +1200(90) 0/00 for Monahans and delta 300(96) 0/00 to +90(98)0/00 for Zag. Delta O-17 of aqueous fluids range over delta 16(22) 0/00 to +18(10) 0/00 for Monahans and +3(10) 0/00 to +27(11) 0/00 for Zag. These variations are larger than the reproducibility of standard analyses and suggest that isotope equilibria were under way in the fluids before trapping into halite. The mean values of delta D and Delta O-17 are +290 0/00 and +9 0/00, respectively. The mean values and the variations of these fluids are different from the representative values of ordinary chondrites, verifying our working hypothesis that the fluid inclusion-bearing halites were not indigenous to the H chondrite parent-asteroid but rather represent exogenous material delivered onto the asteroid from a separate cryovolcanically-active body. This initial isotopic work has demonstrated the feasibility of the measurements, but also revealed sample processing and analytical shortcomings that are now being addressed. Examination of solid mineral inclusions within Monahans and Zag halite grains by confocal Raman spectroscopy at the Carnegie Geophysical Laboratory has revealed them to be metal, magnetite, forsteritic olivine (Fo.98), macromolecular carbon (MMC), pyroxenes, feldspar with Raman spectral affinity to anorthoclase and, probably, fine-grained lepidocrocite (FeO(OH)). In addition, one inclusion features aliphatic material with Raman spectral features consistent with a mixture of short-chain aliphatic compounds. We have initiated analyses of the bulk composition of the fluids within the inclusions in Zag and Monahans halites at Virginia Tech by LA ICPMS using angilent 7500ce quadrupole ICPMS and a Lambda Physik GeoLas 193 nm excimer laser ablation system. Preliminary results reveal that the inclusion aqueous fluids contain highly charged cations of Ca, Mg and Fe. The minerals and compounds discovered thus far within Monahans/Zag halites are indicative of an originating body at least partly composed of unequilibrated anhydrous materials (high Fo olivine, pyroxenes, feldspars, possibly the metal) which were subjected to aqueous alteration (the halite parent brine) and containing a light organic component (the short-chain aliphatic compounds). This material was ejected from the originating body with little or no disruption, as evidenced with the presence of fluid inclusions. An actively geysering body similar to modern Enceladus (Postberg et al., 2011) may be a reasonable analogue in this respect. Also, the originating body should have been within close proximity to the H chondrite parent in order to generate the number of halite grains seen in Monahans and Zag. Other candidates for Monahans/Zag halite parent bodie(s) may include a young Ceres with its possible liquid ocean, or Main Belt comets.

Zolensky, Michael E.↗
DOE OSTI2021

Noise and error analysis and optimization in particle-based kinetic plasma simulations

In this paper we analyze the noise in macro-particle methods used in plasma physics and fluid dynamics, leading to approaches for minimizing the total error, focusing on electrostatic models in one dimension. We begin by describing kernel density estimation for continuous values of the spatial variable x, expressing the kernel in a form in which its shape and width are represented separately. The covariance matrix of the noise in the density is computed, first for uniform true density. The bandwidth of the covariance matrix C(x,y) is related to the width of the kernel. A feature that stands out is the presence of constant negative terms in the elements of the covariance matrix both on and off-diagonal. These negative correlations are related to the fact that the total number of particles is fixed at each time step; they also lead to the property ∫C(x,y)dy = 0. We investigate the effect of these negative correlations on the electric field computed by Gauss's law, finding that the noise in the electric field is related to a process called the Ornstein-Uhlenbeck bridge, leading to a covariance matrix of the electric field with variance significantly reduced relative to that of a Brownian process. For non-constant density, p(x), still with continuous x, we analyze the total error in the density estimation and discuss it in terms of bias-variance optimization (BVO). For some characteristic length l, determined by the density and its second derivative, and kernel width h, having too few particles within h leads to too much variance; for h that is large relative to l, there is too much smoothing of the density. The optimum between these two limits is found by BVO. For kernels of the same width, it is shown that this optimum (minimum) is weakly sensitive to the kernel shape. Next, we repeat the analysis for x discretized on a grid. In this case the charge deposition rule is determined by a particle shape. An important property to be respected in the discrete system is the exact preservation of total charge on the grid; this property is necessary to ensure that the electric field is equal at both ends, consistent with periodic boundary conditions. We find that if the particle shapes satisfy a partition of unity property, the particle charge deposited on the grid is conserved exactly. Further, if the particle shape is expressed as the convolution of a kernel with another kernel that satisfies the partition of unity, then the particle shape obeys the partition of unity. This property holds for kernels of arbitrary width, including widths that are not integer multiples of the grid spacing. Furthermore, we show results relaxing the approximations used to do BVO optimization analytically, by doing numerical computations of the total error as a function of the kernel width, on a grid in x. The comparison between numerical and analytical results shows good agreement over a range of particle shapes. We discuss the practical implications of our results, including the criteria for design and implementation of computationally efficient particle shapes that take advantage of the developed theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Sparse Control Synthesis for Uncertain Responsive Loads With Stochastic Stability Guarantees

In this report, recent studies have demonstrated the potential of flexible loads in providing frequency response services, predominantly due to their availability and cost-effectiveness. However, uncertainty and variability in various weather-related and end-use behavioral factors often impact the reliability of demand-side control performance. This work addresses this problem with the design of a demand-side control to achieve frequency response under load uncertainties. Our approach involves modeling the load uncertainties via stochastic processes that appear as both multiplicative and additive in the power system dynamics. Recently developed mean square exponential stability (MSES) results for continuous-time linear stochastic systems are applied to pose the control synthesis problem which results in an LMI-based optimization problem. Additional costs and constraints are added to the LMI-based controller synthesis to ensure MSES, improve closed-loop transient performance, maximize tolerable uncertainties, and promote sparsity in the controller. Additionally, the fundamental limitations between the tolerable uncertainties and control efforts while ensuring MSES are discussed. Further, the control synthesis problem for the case of the full-state measurement is generalized to the case of partial-state measurements. The proposed control synthesis is illustrated on an IEEE 39 bus system with rigorous studies to demonstrate the role of sparsity, closed-loop transient performance, tolerable uncertainties, and control efforts while ensuring MSES and achieving frequency response.

42 ENGINEERING↗
NASA NTRS2019

High Velocity Impact Performance of a Dual Layer Thermal Protection System for the Mars Sample Return Earth Entry Vehicle

The Mars Sample Return (MSR) Earth Entry Vehicle (EEV) is currently planned on being released from its Micro-Meteorite/Orbital Debri (MM/OD) shielding housing about two days before the Earth entry phase. This leaves the EEV exposed to incoming MM/OD impacts, potentially damaging the heat shield and compromising its Entry, Decent, and Landing (EDL) integrity. Currently, two materials are proposed to comprise the MSR-EEV heat shield, a dual layer material Heat-shield for Extreme Entry Environment Technology (HEEET) and Phenolic Impregnated Carbon Ablator (PICA). PICA has been well characterized for OD class impacts, ~7km/s high mass impacts, from previous testing done in the Orion program, but hasn't under-gone extensive MM impact testing. HEEET is a rela-tively new material with minimal prior testing in re-gards to High Velocity Impacts (HVI). In order to in-form selectability of a material, it is crucial to under-stand the material performance when faced with an HVI, directly affecting mission success probability.The current measure of a Thermal Protection Sys-tem's (TPS) performance against an HVI is evaluating a thermally sized material against its derived Balistic Limit Equation (BLE). A BLE is generated empirically from multiple shots of HVI testing, and is used as a first order method in evaluating TPS's performance against the expected MM/OD environment. This method has proved useful for previous uniform densi-ty TPS materials, but has never been validated against a dual layer recession material such as HEEET. Testing in the MSR program for FY19 has a re-quirement to assess the effects of a dual layer TPS by testing various thicknesses of HEEET's Recession Layer, seen in Figure 1. From this data, BLE's will be derived for the individual thickness ratio samples, as well as the material as a whole to evaluate if a heritage form of the BLE can capture the complex physics associated with a dual layer system. Crater morpholo-gy will also be assessed with post-test Non-Destructive Evaluation (NDE) methods such as CT scanning to visualize if a BLE can well predict the associated pene-tration depths, since a BLE assumes full disinigration of the impacting particle and is generally only used to size a spherical crater ? disregarding any shrapenel effects from a high density impactor.To inform this analysis, expected MM environ-ments from the Meteoroid Engineering Model (MEM) and analytical equations for mass flux of incoming MM were evaluated against the notional MSR-EEV trajectory [1]. Using those dispersions, a monte-carlo was run to determine the most probable particle pa-rameters, as well as the riskiest in terms of full bondline penetration. From these probabilities, a test matrix was designed to test against bounding cases for the various parameters of the BLE: projectile density, projectile mass, projectile velocity, and the impact angle.This presentation will discuss the performance of the dual layer TPS material HEEET against a wide range of impact kinetic energies and densitites, as well as the comparison of HEEET to PICA in terms of MM/OD performance and selectability criteria.

Libben, B. J.↗
DOE OSTI2023

Impact of phase state and non-ideal mixing on equilibration timescales of secondary organic aerosol partitioning

Abstract. Evidence has accumulated that secondary organic aerosols (SOAs) exhibit complex morphologies with multiple phases that can adopt amorphous semisolid or glassy phase states. However, experimental analysis and numerical modeling on the formation and evolution of SOA still often employ equilibrium partitioning with an ideal mixing assumption in the particle phase. Here we apply the kinetic multilayer model of gas–particle partitioning (KM-GAP) to simulate condensation of semi-volatile species into a core–shell phase-separated particle to evaluate equilibration timescales of SOA partitioning. By varying bulk diffusivity and the activity coefficient of the condensing species in the shell, we probe the complex interplay of mass transfer kinetics and the thermodynamics of partitioning. We found that the interplay of non-ideality and phase state can impact SOA partitioning kinetics significantly. The effect of non-ideality on SOA partitioning is slight for liquid particles but becomes prominent in semisolid or solid particles. If the condensing species is miscible with a low activity coefficient in the viscous shell phase, the particle can reach equilibrium with the gas phase long before the dissolution of concentration gradients in the particle bulk. For the condensation of immiscible species with a high activity coefficient in the semisolid shell, the mass concentration in the shell may become higher or overshoot its equilibrium concentration due to slow bulk diffusion through the viscous shell for excess mass to be transferred to the core phase. Equilibration timescales are shorter for the condensation of lower-volatility species into semisolid shell; as the volatility increases, re-evaporation becomes significant as desorption is faster for volatile species than bulk diffusion in a semisolid matrix, leading to an increase in equilibration timescale. We also show that the equilibration timescale is longer in an open system relative to a closed system especially for partitioning of miscible species; hence, caution should be exercised when interpreting and extrapolating closed-system chamber experimental results to atmosphere conditions. Our results provide a possible explanation for discrepancies between experimental observations of fast particle–particle mixing and predictions of long mixing timescales in viscous particles and provide useful insights into description and treatment of SOA in aerosol models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

SR19021 Tritium Aging of Regenerated LANA.75 (Final Project Report)

The Savannah River Tritium Enterprise (SRTE) has used the metal hydride LaNi 4.25 Al 0.75 (LANA.75) in the Tritium Facilities for over two decades. LANA.75 beds store significant quantities of tritium but have a limited service life due to the radiolytic decay of tritium to He-3 within the metal matrix. It has been shown that heating tritium-aged LANA.75 under vacuum can reverse tritium aging effects, eliminating the heel of trapped hydrogen, and restoring the reversible capacity. Additional investigation is needed to ensure there are no unexpected changes to the hydride before this restoration technique is employed in full scale beds in the Tritium Facilities. This project was to be comprised of three distinct scopes: obtain tritium aging data on the sample regenerated in 2018, regenerate a tritium-aged sample at 600 °C, and perform thermal stability testing on a non-tritiated sample. Isotherms were collected on the previously regenerated sample after approximately two years of tritium aging. Isotherms were collected at 80, 100, and 120 °C. As expected, there was a decrease in the plateau pressure, an increase in plateau slope, and a portion of the “heel” of tritium trapped in the metal had been reestablished. Unexpectedly, it appeared that the plateau had shortened at the higher tritium to metal ratios as well. This is typically seen in older samples. The second scope, to regenerate a second tritium-aged LANA.75 sample, was not completed. A Task Technical and Quality Assurance Plan was written and approved, a high temperature test cell was fabricated, several pre-job briefs were held, and the hydride sample was passivated with air. Despite these successes, the hydride sample was not recovered from the legacy test cell. The third scope was to perform thermal stability testing on a “cold” LANA.75 sample. A non-tritiated sample of LANA.75 was held at 750 °C under vacuum for 200 hours to simulate exposure to multiple regeneration evolutions. Hydride isotherm performance, chemical composition, crystallinity, particle size, and morphology are compared between the pre- and post-regeneration samples. No significant changes were observed in composition, crystallinity, or particle size. Comparison of before and after isotherms showed that performance improved rather than deteriorated during the evolution. Scanning Electron Microscopy (SEM) analysis showed small growths on the particle surface after exposure to regeneration conditions. Additional testing will be required to determine the cause of these growths.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2026

Mechanical Behavior of Neutron Irradiated Refractory Multi-Principal Element Alloys Processed via Spark Plasma Sintering

The search for advanced materials capable of withstanding the extreme conditions of Generation IV reactors is a critical area in materials science research. These reactors operate under severe environments, including high temperatures, corrosion, stress, and irradiation damage. Consequently, there is a need for innovative alloy systems to ensure the reliability and longevity of proposed Generation IV reactor components. Refractory multi-principal-element alloys (RMPEA) have emerged as a promising candidate due to their exceptional properties. These alloys, characterized by their composition of multiple principal elements in near-equiatomic ratios, exhibit superior resistance to irradiation damage, reduced void swelling, enhanced microstructural stability, and minimal irradiation-induced hardening. While initial studies on RMPEAs have shown promising results, most research has been limited to thin films, nanocrystalline microstructures, and ion irradiation, which do not accurately represent the behavior of bulk materials. To address this gap, our research focused on the neutron irradiation of bulk RMPEAs. We aim to conduct comprehensive post-irradiation examinations (PIE) of RMPEAs irradiated at the Advanced Test Reactor at Idaho National Laboratory. The RMPEAs were synthesized using spark plasma sintering (SPS) with mechanically alloyed metallurgical powder. The RMPEA specimens are a MoNbTi alloy system with additions of -Zr, and -ZrV. Furthermore, PIE consisted of mechanical testing and advanced materials characterization. The mechanical testing consisted of sub-sized tensile testing, micro- and nano- indentation. Microstructural characterization included scanning electron microscopy and transmission electron microscopy. Mechanical testing coupled with advanced microscopy techniques provides insight into phase morphology and its effects on the mechanical properties of the RMPEA specimens. The results indicate that both pristine and irradiated RMPEA specimens exhibited brittle behavior during tensile testing, which can be attributed to their heterogeneous microstructure. The SPS manufacturing process did not include any post treatment, which resulted in a heterogeneous microstructure. Energy-dispersive X-ray spectroscopy revealed the presence of intermetallic such as laves phases within the microstructure. Specifically, Ti-rich precipitates were observed in the MoNbTi specimen, while Mo-rich precipitates were found in the MoNbTiZrV specimen. Nano-hardness testing of pristine samples showed that the laves phases exhibited higher hardness values compared to the matrix phase, suggesting that precipitate hardening is likely the dominant hardening mechanism in these specimens. The results from this work will be used to build a finite element model to predict mechanical behavior of future MPEA compositions. Thus, enabling for a streamlined approach to developing novel MPEAs for the nuclear industry.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

Carbon Transport and Storage Planning and Viability Support Tools

The EDX disCO2ver Carbon Transport and Storage Planning and Viability Support Tools are made up of the Carbon Storage Planning Inquiry Tool (CS PlanIT, Justman et al. 2024) and the Carbon Storage Technical Viability Approach Support Tool (CS TVA). Together, these tools support data access to support understanding data availability to support planning efforts for carbon transport and storage. The Carbon Storage Planning Inquiry Tool (CS PlanIT) is an online web mapping application designed to help users explore, query, and evaluate multiple data layers to support and accelerate carbon storage resource and feasibility assessments and planning efforts. CS PlanIT currently contains a range of datasets associated with geologic, technical, and infrastructure factors. The data sets can be filtered geographically for an area of interest to update statistics and charts within the dashboard. The dashboard is divided into different sections called widgets, relating to different steps in the carbon storage planning process. The resources in this submission include a link to PlanIT, as well as a data catalog and link to user documentation. The original citation for the CS PlanIT tool, which has now been integrated into the toolset here, was: - Devin Justman, Scott Pantaleone, Maneesh Sharma, Lucy Romeo, Paige Morkner, CS PlanIT (Carbon Storage Planning Inquiry Tool) , 6/28/2024, https://edx.netl.doe.gov/dataset/cs-planit-carbon-storage-planning-inquiry-tool, DOI: 10.18141/2377953 The Carbon Storage Technical Viability Approach Support (CS TVA) Tool displays spatial data availability for the many components of Geologic Carbon Storage (GCS) technical viability assessment (Creason al 2025). Identifying sites suitable for GCS requires evaluating the intersection of myriad factors, including reservoir conditions, subsurface and surface hazards, infrastructure requirements, and energy community metrics. The technical viability of a site can only be confirmed for instances where all these factors have data available, and where those data support viability. Additional Resources related to the Technical Viability Assessment Tool: - Julia Mulhern, Casey White, Araceli Lara, Neyda Cordero Rodriguez, Zachary Jackson, Jacob Shay, Gabriel Creason, MacKenzie Mark-Moser, Paige Morkner, Kelly Rose, Carbon Storage Technical Viability Approach (CS TVA) Database, 3/26/2025, https://edx.netl.doe.gov/dataset/edx4ccs-carbon-storage-technical-viability-approach-database , DOI:10.18141/1984655 - Julia Mulhern, MacKenzie Mark-Moser, Gabriel Creason, Casey White, Araceli Lara, Neyda Cordero Rodriguez, Zach Jackson, Paige Morkner, Kelly Rose, Carbon Storage Technical Viability Approach (CS TVA) Matrix, 3/27/2025, https://edx.netl.doe.gov/dataset/carbon-storage-technical-viability-approach-cs-tva-matrix , DOI: 10.18141/2539979 - Gabriel Creason, Zach Jackson, Neyda Cordero Rodriguez, Julia Mulhern, Casey White, Araceli Lara, MacKenzie Mark-Moser, Paige Morkner, Kelly Rose, Carbon Storage Technical Viability Approach (CS TVA) Data Availability Results Database, 3/27/2025, https://edx.netl.doe.gov/dataset/carbon-storage-technical-viability-approach-cs-tva-data-availability-results-database, DOI:10.18141/2538557

Carbon storage↗
DOE OSTI2021

Pore-Scale Dynamics of Liquid CO2–Water Displacement in 2D Axisymmetric Porous Micromodels Under Strong Drainage and Weak Imbibition Conditions: High-Speed μPIV Measurements

Resolving pore-scale transient flow dynamics is crucial to understanding the physics underlying multiphase flow in porous media and informing large-scale predictive models. Surface properties of the porous matrix play an important role in controlling such physics, yet interfacial mechanisms remain poorly understood, in part due to a lack of direct observations. This study reports on an experimental investigation of the pore-scale flow dynamics of liquid CO 2 and water in two-dimensional (2D) circular porous micromodels with different surface characteristics employing high-speed microscopic particle image velocimetry (μPIV). The design of the micromodel minimized side boundary effects due to the limited size of the domain. The high-speed μPIV technique resolved the spatial and temporal dynamics of multiphase flow of CO 2 and water under reservoir-relevant conditions, for both drainage and imbibition scenarios. When CO 2 displaced water in a hydrophilic micromodel (i.e., drainage), unstable capillary fingering occurred and the pore flow was dominated by successive pore-scale burst events (i.e., Haines jumps). When the same experiment was repeated in a nearly neutral wetting micromodel (i.e., weak imbibition), flow instability and fluctuations were virtually eliminated, leading to a more compact displacement pattern. Energy balance analysis indicates that the conversion efficiency between surface energy and external work is less than 30%, and that kinetic energy is a disproportionately smaller contributor to the energy budget. This is true even during a Haines jump event, which induces velocities typically two orders of magnitude higher than the bulk velocity. These novel measurements further enabled direct observations of the meniscus displacement, revealing a significant alteration of the pore filling mechanisms during drainage and imbibition. While the former typically featured burst events, which often occur only at one of the several throats connecting a pore, the latter is typically dominated by a cooperative filling mechanism involving simultaneous invasion of a pore from multiple throats. This cooperative filling mechanism leads to merging of two interfaces and releases surface energy, causing instantaneous high-speed events that are similar, yet fundamentally different from, burst events. Finally, pore-scale velocity fields were statistically analyzed to provide a quantitative measure of the role of capillary effects in these pore flows.

Li, Yaofa↗
NASA NTRS1993

Relationships between pigment composition variation and reflectance for plant species from a coastal savannah in California

Advances in imaging spectroscopy have indicated that remotely sensed reflectance measurements of the plant canopy may be used to identify and qualify some classes of canopy biochemicals; however, the manner in which differences in biochemical compositions translate into differences is not well understood. Most frequently, multiple linear regression routines have been used to correlate narrow band reflectance values with measured biochemical concentrations. Although some success has been achieved with such methods for given data sets, the bands selected by multiple regression are not consistent between data sets, nor is it always clear what physical or biological basis underlies the correlation. To examine the relationship between biochemical concentration and leaf reflectance signal we chose to focus on the visible spectrum where the primary biochemical absorbances are due to photosynthetic pigments. Pigments provide a range of absorbance features, occur over a range of concentrations in natural samples, and are ecophysiologically important. Concentrations of chlorophyll, for example, have been strongly correlated to foliar nitrogen levels within a species and to photosynthetic capacity across many species. In addition pigments effectively absorb most of the photosynthetically active radiation between 400-700 nm, a spectral region for which silicon detectors have good signal/noise characteristics. Our strategy has been to sample a variety of naturally occurring species to measure leaf reflectance and pigment compositions. We hope to extend our understanding of pigment reflectance effects to interpret small overlapping absorbances of other biochemicals in the infrared region. For this reason, selected samples were also tested to determine total nitrogen, crude protein, cellulose, and lignin levels. Leaf reflectance spectra measured with AVIRIS bandwidths and wavelengths were compared between species and within species and for differences between seasons, for changes in the the shape of the spectra. We attempt to statistically correlate these shape changes with differences in pigment compositions. In parallel with our comparisons of pigment composition and leaf reflectance, we have modified the PROSPECT leaf reflectance model to test the contributions of pigments or pigment group concentrations. PROSPECT considers a leaf as a multi-layer dielectric plane with an uneven surface. Jacquemoud adapted the basic analysis of Allen for surface effects, a leaf thickness factor, and the absorption of water and chlorophyll (actually all pigments) and the plant matrix. Our modifications to PROSPECT in the forward direction include breaking out the pigment concentration parameter into separate components for chlorophyll a and b and a number of xanthophylls and carotenes, and introducing a shift and convolution function to model the spread and shift from their in vitro measurements to their in vivo state. Further, we have considered how the matrix elements (i.e., all biochemicals and structural effects not modeled explicity) vary with species.

Ustin, Susan L.↗