Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model compounds”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Hybrid molecular beam epitaxy of germanium-based oxides

Abstract Germanium-based oxides such as rutile GeO 2 are garnering attention owing to their wide band gaps and the prospects of ambipolar doping for application in high-power devices. Here, we present the use of germanium tetraisopropoxide (GTIP), a metal-organic chemical precursor, as a source of germanium for the demonstration of hybrid molecular beam epitaxy for germanium-containing compounds. We use Sn 1- x Ge x O 2 and SrSn 1- x Ge x O 3 as model systems to demonstrate our synthesis method. A combination of high-resolution X-ray diffraction, scanning transmission electron microscopy, and X-ray photoelectron spectroscopy confirms the successful growth of epitaxial rutile Sn 1- x Ge x O 2 on TiO 2 (001) substrates up to x = 0.54 and coherent perovskite SrSn 1- x Ge x O 3 on GdScO 3 (110) substrates up to x = 0.16. Characterization and first-principles calculations corroborate that germanium occupies the tin site, as opposed to the strontium site. These findings confirm the viability of the GTIP precursor for the growth of germanium-containing oxides by hybrid molecular beam epitaxy, thus providing a promising route to high-quality perovskite germanate films.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond

Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within “chemical accuracy” (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob’s ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 − N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.

Chemistry↗
DOE OSTI2023

Characterizing random-singlet state in two-dimensional frustrated quantum magnets and implications for the double perovskite Sr 2 CuTe 1 – x W x O 6

Motivated by the experimental observation of a nonmagnetic phase in compounds with frustration and disorder, we study the ground state of a spin-1/2 square-lattice Heisenberg model with randomly distributed nearest-neighbor J 1 and next-nearest-neighbor J 2 couplings. By using the density matrix renormalization group (DMRG) calculation on a cylinder system with a circumference of up to ten lattice sites, we identify a disordered phase between the Néel and stripe magnetic phase with growing J 2 /J 1 in the presence of strong bond randomness. The vanished spin-freezing parameter indicates the absence of spin-glass order. The large-scale DMRG results unveil the size-scaling behaviors of the spin-freezing parameter, the power-law decay of the average spin correlation, and the exponential decay of the typical spin correlation, which all agree with the corresponding behavior in the one-dimensional random-singlet (RS) state and characterize the RS nature of this disordered phase. The DMRG simulation also provides insights and opportunities for characterizing a class of nonmagnetic states in two-dimensional frustrated magnets with disorder. Here, we also compare with existing experiments and suggest more measurements for understanding the spin-liquid-like behaviors in the double perovskite Sr 2 CuTe 1–x W x O 6 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Partitioning of Complex Fluids at Mineral Surfaces

This report summarizes the results obtained during the LDRD project entitled "Partitioning of Complex Fluids at Mineral Interfaces." This research addressed fundamental aspects of such interfaces, which are relevant to energy-water applications in the subsurface, including fossil energy extraction and carbon sequestration. This project directly addresses the problem of selectivity of complex fluid components at mineral-fluid interfaces, where complex fluids are defined as a mixture of hydrophobic and hydrophilic components: e.g., water, aqueous ions, polar/nonpolar organic compounds. Specifically, this project investigates how adsorption selectivity varies with surface properties and fluid composition. Both experimental and molecular modeling techniques were used to better understand trends in surface wettability on mineral surfaces. The experimental techniques spanned the macroscale (contact angle measurements) to the nanoscale (cryogenic electronic microscopy and vibrational spectroscopy). We focused on an anionic surfactant and a well-characterized mineral phase representative of clay phases present in oil- and gas-producing shale deposits. Collectively, the results consistently demonstrate that the presence of surfactant in the aqueous fluid significantly affects the mineral-fluid interfacial structure. Experimental and molecular modeling results reveal details of the surfactant structure at the interface, and how this structure varies with surfactant coverage and fluid composition.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2022

Less Common Topological Phenomena in Bulk Materials

"Prediction and subsequent discovery of topological insulators is considered to be one of the main results in condensed matter physics in the last decades. Not surprisingly, it has received major attention of both researchers and funding agencies. This attention is well-deserved; yet, one cannot but note that nearly all this research revolves about essentially the same concept: electronic excitations with linear dispersion, covering, of course, such diverse and intriguing phenomena as topological insulator, bulk Dirac states (or Weyl, if not spin degenerate), Mayorana fermions. In this project, we will address, mainly, other topological phenomena, such as topologically nontrivial magnetic patterns (as, for instance, topological Hall and related phenomena). Specifically, we propose three interrelated trusts: (1) Time-reversal symmetry breaking nonrelativistic antiferromagnets, called altermagnets. These are materials that break Kramers degeneracy of electronic bands, despite having zero net magnetization by symmetry and being fully collinear, and not necessarily non-centrosymmetric. The corresponding band structure is very similar to the band structure in non-centrosymmetric spin-orbital materials, but materials that we propose to study are distinctly different, first and foremost in the sense that despite sharing many aspects of their electronic properties with the latter, they break the time-reversal symmetry without either spin-orbit coupling or lack of inversion symmetry. (2) Topologically nontrivial magnetic spirals. The PI has been engaged with the experimental group of Dr. Ghimire at GMU investigating Dirac materials with helical magnetism, based on stacked magnetic Kagome layers, with a generic formula of RMn6Sn6. In particular, Y Mn6Sn6 demonstrates a component of the Hall effect that is naturally interpreted in terms of a topological spin texture, as well as linear magnetoresistance. Our calculation identify Dirac states that are robust with respect to the spiral formation, and let us derive an advance mean-field model explaining the observed phase diagram. This model predicted four distinct phases, with very distinct properties, which have now been seen in neutron experiments. The same compound is known to demonstrate topological Hall effect in a particular magnetic phase, and only at elevated temperature. Based on our understanding of the phase diagram, we have worked out a phenomenological theory of a chiral (skyrmionic) response to an external magnetic field, similar to the nematic response to external strain in Fe-based superconductors, which is possible in a centrosymmetric lattice and without interplanar Dzyaloshinskii-Moriya interaction. This phenomenological theory agrees quantitatively with the experiment. It is in our plans to research other similar materials for this effect. (3) Search for 3D analogues of Fe-based superconductors. We want to investigate materials that can be viewed as 3D analogues of FeSe. Specifically, we want materials that are good metals and host antiferromagnetism, which can be suppressed by pressure and generate an s-wave superconductivity, as in Fe-based superconductors. We have in mind some candidates already. This work will proceed in close collaboration with the experimental group of Prof. Nirmal Ghimire in the same department, whose expertise lies in sample making, magnetometry and transport measurements of materials with complex magnetic structures."

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Housing-Performance Atlas of Baltimore Row Homes: Archetype-Based Multi-Hazard Baseline of Energy, Heat, Survivability, and Durability

Baltimore’s historic row-home neighborhoods face escalating risks to energy, heat, and durability under intensifying climate stress. This study develops a Housing-Performance Atlas that quantifies multi-hazard performance for eight representative archetypes using DesignBuilder/EnergyPlus Version 7.3.1.003, under Baltimore TMY3 boundary conditions. Performance is evaluated across the following four adaptation domains: energy use intensity, passive survivability during 72 h outage events, roof overheating exposure (>150 °F exceedance hours), and material service life derived from ISO 15686 and synthesized into Lean and Full Deficit Indices for comparative resilience ranking. Results show that EUI ranged from 46.7 to 67.6 kBtu ft −2 ·yr −1 , survivability from 0 to 23 h, and roof temperatures exceeded 150 °F for 150–210 h, shortening roof service life by up to 10 years. Composite Lean and Full Deficit Indices ranged 7.8–92.4, ranking Model 5 (end-unit, flat roof, two-story with basement) as the most resilient configuration and Model 8 (end-unit, pitched roof, three-story above-grade) as the least resilient due to compounded overheating and energy losses. Heat-related domains accounted for nearly 70% of overall resilience deficits, confirming thermal safety and roof reflectivity as retrofit priorities. The Housing-Performance Atlas establishes a reproducible diagnostic framework linking simulation, service life, and resilience metrics to guide cost-effective, climate-responsive retrofits in Baltimore’s aging urban housing stock.

Housing-Performance atlas↗
DOE OSTI2025

Large language model-driven database for thermoelectric materials

Thermoelectric materials have the ability to convert waste heat into electricity, offering a valuable solution for energy harvesting. However, their widespread use is hindered by low conversion efficiency, the reliance on expensive rare earth elements, and the environmental and regulatory concerns associated with lead-based materials. A fast and cost-effective way to identify highly efficient thermoelectric materials is through data-driven methods. These approaches rely on robust and comprehensive datasets to train models. Although there are several databases on thermoelectric materials, there is still a need to collect and integrate experimental data from peer-reviewed research articles to capture diverse compositions and properties of materials. Here, in this work, we developed a comprehensive database of 7,123 thermoelectric compounds, containing key information such as chemical composition, structural detail, seebeck coefficient, electrical and thermal conductivity, power factor, and figure of merit (ZT). We used the GPTArticleExtractor workflow, powered by large language models (LLM), to extract and curate data automatically from the scientific literature published in Elsevier journals. This process enabled the creation of a structured database that addresses the challenges of manual data collection. The open access database could stimulate data-driven research and advance thermoelectric material analysis and discovery.

Database↗
DOE OSTI2021

Modeling of Peptides with Classical and Novel Machine Learning Force Fields: A Comparison

The replacement of classical force fields (FFs) with novel neural-network-based frameworks is an emergent topic in molecular dynamics (MD) simulations. In contrast to classical FFs, which have proven their capability to provide insights into complex soft matter systems at an atomistic resolution, the machine learning (ML) potentials have yet to demonstrate their applicability for soft materials. However, the underlying philosophy, which is learning the energy of an atom in its surrounding chemical environment, makes this approach a promising tool. In particular for the exploration of novel chemical compounds, which have not been considered in the original parametrization of classical FFs. In this article, we study the performance of the ANI-2x ML model and compare the results with those of two classical FFs, namely, CHARMM27 and the GROMOS96 43a1 FF. We explore the performance of these FFs for bulk water and two model peptides, trialanine and a 9-mer of the α-aminoisobutyric acid, in vacuum and water. The results for water describe a highly ordered water structure, with a structure similar to those using ab initio molecular dynamics simulations. The energy landscape of the peptides described by Ramachandran maps show secondary structure basins similar to those of the classical FFs but differ in the position and relative stability of the basins. Details of the sampled structures show a divergent performance of the different models, which can be related either to the short-ranged nature of the ML potentials or to shortcomings of the underlying data set used for training. These findings highlight the current state of the applicability of ANI-2x ML potential for MD simulations of soft matter systems. Simultaneously, they provide insights for future improvements of current ML potentials.

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

First-principles equation of state of CHON resin for inertial confinement fusion applications

A wide-range (0 to 1044.0 g/ cm 3 and 0 to 10 9 K) equation-of-state (EOS) table for a CH 1.72 O 0.37 N 0.086 quaternary compound has been constructed based on density-functional theory (DFT) molecular-dynamics (MD) calculations using a combination of Kohn-Sham DFT MD, orbital-free DFT MD, and numerical extrapolation. The first-principles EOS data are compared with predictions of simple models, including the fully ionized ideal gas and the Fermi-degenerate electron gas models, to chart their temperature-density conditions of applicability. The shock Hugoniot, thermodynamic properties, and bulk sound velocities are predicted based on the EOS table and compared to those of C-H compounds. The Hugoniot results show the maximum compression ratio of the C-H-O-N resin is larger than that of CH polystyrene due to the existence of oxygen and nitrogen; while the other properties are similar between CHON and CH. Radiation hydrodynamic simulations have been performed using the table for inertial confinement fusion targets with a CHON ablator and compared with a similar design with CH. Here, the simulations show CHON outperforms CH as the ablator for laser-direct-drive target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Impact of trapping on tritium self-sufficiency and tritium inventories in fusion power plant fuel cycles

The dynamic analysis of fusion power plant (FPP) fuel cycles highlights the challenge of achieving tritium self-sufficiency in future FPPs. While state-of-the-art fuel cycle models offer valuable insights into the necessary design parameters for attaining tritium self-sufficiency, none of these models currently consider the impact of tritium trapping within fuel cycle components. However, detailed analysis of individual components reveals that substantial amounts of tritium can be trapped within the first wall, divertors, and breeding blanket systems, suggesting that tritium trapping may significantly influence the FPP ability to achieve self-sufficiency. The compounded effects of additional tritium traps generated by irradiation effects and component replacements further exacerbate this challenge. The novelty of this work is the integration of an explicit, physics-based model for tritium trapping, evolution of damage-induced traps, and component replacements into a dynamic, system-level model of a fuel cycle. The results show an increase of a factor 10 3 – 10 4 of tritium inventory in the first wall and vacuum vessel of an ARC-class FPP when accounting for the aforementioned phenomena. This, coupled with the replacement of components subject to significant tritium trapping, slows down fuel cycle dynamics, resulting in an extended tritium doubling time (50% increase), higher start-up inventory (30% increase), and higher required tritium breeding ratio (2%–5%) compared to a scenario without tritium trapping.

fuel cycle↗
DOE OSTI2023

A standardized quantitative analysis strategy for stable isotope probing metagenomics

ABSTRACT Stable isotope probing (SIP) facilitates culture-independent identification of active microbial populations within complex ecosystems through isotopic enrichment of nucleic acids. Many DNA-SIP studies rely on 16S rRNA gene sequences to identify active taxa, but connecting these sequences to specific bacterial genomes is often challenging. Here, we describe a standardized laboratory and analysis framework to quantify isotopic enrichment on a per-genome basis using shotgun metagenomics instead of 16S rRNA gene sequencing. To develop this framework, we explored various sample processing and analysis approaches using a designed microbiome where the identity of labeled genomes and their level of isotopic enrichment were experimentally controlled. With this ground truth dataset, we empirically assessed the accuracy of different analytical models for identifying active taxa and examined how sequencing depth impacts the detection of isotopically labeled genomes. We also demonstrate that using synthetic DNA internal standards to measure absolute genome abundances in SIP density fractions improves estimates of isotopic enrichment. In addition, our study illustrates the utility of internal standards to reveal anomalies in sample handling that could negatively impact SIP metagenomic analyses if left undetected. Finally, we present SIPmg , an R package to facilitate the estimation of absolute abundances and perform statistical analyses for identifying labeled genomes within SIP metagenomic data. This experimentally validated analysis framework strengthens the foundation of DNA-SIP metagenomics as a tool for accurately measuring the in situ activity of environmental microbial populations and assessing their genomic potential. IMPORTANCE Answering the questions, “who is eating what?” and “who is active?” within complex microbial communities is paramount for our ability to model, predict, and modulate microbiomes for improved human and planetary health. These questions can be pursued using stable isotope probing to track the incorporation of labeled compounds into cellular DNA during microbial growth. However, with traditional stable isotope methods, it is challenging to establish links between an active microorganism’s taxonomic identity and genome composition while providing quantitative estimates of the microorganism’s isotope incorporation rate. Here, we report an experimental and analytical workflow that lays the foundation for improved detection of metabolically active microorganisms and better quantitative estimates of genome-resolved isotope incorporation, which can be used to further refine ecosystem-scale models for carbon and nutrient fluxes within microbiomes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Uncertainties in Simulating Flooding During Hurricane Harvey Using 2D Shallow Water Equations

Abstract Flooding is one of the most impactful weather‐related natural hazards. Numerical models that solve the two dimensional (2D) shallow water equations (SWE) represent the first‐principles approach to simulate all types of spatial flooding, such as pluvial, fluvial, and coastal flooding, and their compound dynamics. High spatial resolution (e.g., () m) is needed in 2D SWE simulations to capture flood dynamics accurately, resulting in formidable computational challenges. Thus, relatively coarser spatial resolutions are used for large‐scale simulations of flooding, which introduce uncertainties in the results. It is unclear how the uncertainty associated with the model resolution compares to the uncertainties in precipitation data sets and assumptions regarding boundary conditions when channelized flows interact with other water bodies. In this study, we compare these three sources of uncertainties in 2D SWE simulations for the 2017 Houston flooding event. Our results show that precipitation uncertainty and mesh resolution have more significant impacts on the simulated streamflow and inundation dynamics than the choice of the downstream boundary condition at the watershed outlet. We point out the viability to confine the uncertainty of coarsening mesh resolution by using the variable resolution mesh (VRM) which refines critical topographic features with far fewer grid cells. Specifically, in simulations with VRM, the simulated inundation depths over the refined region are comparable to that use the finest uniform mesh. This study contributes to understanding the challenges and pathways for applying 2D SWE models to improve the realism of flood simulations over large scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

DFT Calculation of Thermal Expansion within Debye-Grüneisen Framework Made Simple and Accurate

The density functional theory (DFT) based Debye-Grüneisen approach is a computational method used to calculate the thermodynamic properties of materials by relating the volume dependence of the Debye frequency to the Grüneisen parameter. In this study, we propose a new model to describe the volume dependence of the Grüneisen parameter. This model offers a simple yet accurate approach to predict the thermal expansion of materials within the framework of the Debye-Grüneisen theory. Specifically, for pure elements or species used as a reference system, we calibrated the parameters to reproduce the experimental coefficient of thermal expansion, heat capacity, and entropy. For compounds with arbitrary compositions, the parameters are determined by taking the arithmetic or geometric average of those obtained from the pure elements or species based on their composition. We then compared the model prediction with experiments. This allowed for a comprehensive characterization of the thermal expansion behavior of various materials.

Wang, Yi↗
DOE OSTI2023

Systematic study of short- and long-range correlations in RE 3 TaO 7 weberite-type compounds by neutron total scattering and X-ray diffraction

The atomic structures of the lanthanide tantalates, Ln 3 TaO 7 , series (Ln = Pr, Tb, Dy, Ho, Tm, Yb) were systematically investigated using total scattering techniques. High-energy X-ray and neutron diffraction analysis revealed that the long-range structures can be grouped into three distinct families: (1) ordered Cmcm (Ln = Pr), (2) ordered Ccmm (Ln = Tb, Dy, Ho), and (3) disordered, defect-fluorite Fm$\bar{3}$m (Ln = Ho, Tm, Yb). These findings help to clarify the symmetry discrepancy for the already reported long-range structures in the literature. The short-range analysis of neutron total scattering data via pair distribution functions reveals a high degree of structural heterogeneity across length scales for all compounds, with distinct local atomic arrangements that are not fully captured by the average, long-range structure. The short-range structures at the level of coordination polyhedra are better captured by a set of alternative non-centrosymmetric structural models: (1) C2cm, (2) C222 1 , and (3) C2mm. This establishes a short-range multiferroic character for weberite-type tantalates because ferroelectric interactions compete with magnetic correlations. These ferroelectric interactions are particularly pronounced for the disordered compounds Tm 3 TaO 7 and Yb 3 TaO 7 . In conclusion, the structural differences among the three families are the result of changes in TaO 6 polyhedral tilt (transition between families 1 and 2) and dipolar interactions of off-centered Ta cations (transition between families 2 and 3).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling Novel Aqueous Particle and Cloud Chemistry Processes of Biomass Burning Phenols and Their Potential to Form Secondary Organic Aerosols

Phenols emitted from biomass burning contribute significantly to secondary organic aerosol (SOA) formation through the partitioning of semivolatile products formed from gas-phase chemistry and multiphase chemistry in aerosol liquid water and clouds. The aqueous-phase SOA (aqSOA) formed via hydroxyl radical (˙OH), singlet molecular oxygen ( 1 O 2 *), and triplet excited states of organic compounds ( 3 C*), which oxidize dissolved phenols in the aqueous phase, might play a significant role in the evolution of organic aerosol (OA). However, a quantitative and predictive understanding of aqSOA has been challenging. Here, we develop a stand-alone box model to investigate the formation of SOA from gas-phase ˙OH chemistry and aqSOA formed by the dissolution of phenols followed by their aqueous-phase reactions with ˙OH, 1 O 2 *, and 3 C* in cloud droplets and aerosol liquid water. Here, we investigate four phenolic compounds, i.e., phenol, guaiacol, syringol, and guaiacyl acetone (GA), which represent some of the key potential sources of aqSOA from biomass burning in clouds. For the same initial precursor organic gas that dissolves in aerosol/cloud liquid water and subsequently reacts with aqueous phase oxidants, we predict that the aqSOA formation potential (defined as aqSOA formed per unit dissolved organic gas concentration) of these phenols is higher than that of isoprene-epoxydiol (IEPOX), a well-known aqSOA precursor. Cloud droplets can dissolve a broader range of soluble phenols compared to aqueous aerosols, since the liquid water contents of aerosols are orders of magnitude smaller than cloud droplets. Further, our simulations suggest that highly soluble and reactive multifunctional phenols like GA would predominantly undergo cloud chemistry within cloud layers, while gas-phase chemistry is likely to be more important for less soluble phenols. But in the absence of clouds, the condensation of low-volatility products from gas-phase oxidation followed by their reversible partitioning to organic aerosols dominates SOA formation, while the SOA formed through aqueous aerosol chemistry increases with relative humidity (RH), approaching 40% of the sum of gas and aqueous aerosol chemistry at 95% RH for GA. Our model developments of biomass-burning phenols and their aqueous chemistry can be readily implemented in regional and global atmospheric chemistry models to investigate the aqueous aerosol and cloud chemistry of biomass-burning organic gases in the atmosphere.

09 BIOMASS FUELS↗
DOE OSTI2023

Preliminary Nuclear Containment Vessel Modeling for Multi-Hazard Probabilistic Risk Assessment under Seismic Hazards and Concrete Degradation

The current practice for natural phenomena hazards (NPH) risk assessment of nuclear facilities is to compute the risk for each hazard independently and then compound the total risk as a combination of single hazard risks. This state of practice does not consider correlations between hazards and the cascading impacts to structures, systems, and components (SSCs), and could thus underestimate the NPH risk or overestimate the nuclear facility safety. Events such as the Fukushima Daiichi accident have highlighted the importance of multi-hazard risk considerations to nuclear power plants (NPPs) that quantify the cascading damage effects to SSCs in the risk models. Moreover, the current fleet of NPPs in the United States is aging; these NPPs are now expected to operate well beyond their initially planned design life. Aging-related deterioration can potentially decrease the capacity of critical structures such as containment vessels to withstand NPH. Such aging considerations may not be adequately accounted for by the current NPH risk assessment guidelines. This paper presents a preliminary modeling and simulation of a representative reinforced concrete containment vessel subjected to seismic mainshock and aftershock considering concrete degradation due to alkali silica reaction. The broader aim is to develop multi-hazard time-dependent fragility functions that could be subsequently used in the probabilistic risk assessment (PRA) model. The multi-hazard component comes into play due to the consideration of damage to the containment vessel under seismic loads and concrete degradation. Consideration of concrete degradation also brings into play the time-dependent nature of the containment vessel response. The response of the containment vessel under varying degrees of concrete degradation to seismic loads is investigated. The results presented are simulated using the Multi-hazard Analysis for STOchastic time-DOmaiN phenomena (MASTODON) software for seismic analysis and the Blackbear software for concrete degradation and damage modeling. Both software are open source and developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE).

42 ENGINEERING↗
DOE OSTI2022

A New Perspective of Post-Weld Baking Effect on Al-Steel Resistance Spot Weld Properties through Machine Learning and Finite Element Modeling

The root cause of post-weld baking on the mechanical performance of Al-steel dissimilar resistance spot welds (RSWs) has been determined by machine learning (ML) and finite element modeling (FEM) in this study. A deep neural network (DNN) model was constructed to associate the spot weld performance with the joint attributes, stacking materials, and other conditions, using a comprehensive experimental dataset. The DNN model positively identified that the post-weld baking reduces the joint performance, and the extent of degradation depends on the thickness of stacking materials. A three-dimensional finite element (FE) model was then used to investigate the root cause and the mechanism of the baking effect. It revealed that the formation of high thermal stresses during baking, from the mismatch of thermal expansion between steel and Al alloy, causes damage and cracking of the brittle intermetallic compound (IMC) formed at the interface of the weld nugget during welding. This in turn reduces the joint performance by promoting undesirable interfacial fracture when the welds were subjected to externally applied loads. The FEM model further revealed that increase in structural stiffness, because of increase in steel sheet thickness, reduces the thermal stresses at the interface caused by the thermal expansion mismatch and consequently lessens the detrimental effect of post-weld baking on the joint performance.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Enter Gaussian Mixture Modeling Extensions for Improved False Discovery Rate Estimation in GC-MS Metabolomics

Identifying small molecules (e.g., metabolites) is key towards driving scientific advancement in metabolomics, and gas chromatography–mass spectrometry (GC-MS) is an analytic method that may be applied to facilitate this process. The typical GC-MS identification workflow involves quantifying the similarity of an observed sample spectrum and other features (e.g. retention index) to that of several references, noting the compound of the best-matching reference spectrum as the identified metabolite. While a deluge of similarity metrics exists, none characterize the error rate of generated identifications, thereby presenting an unknown risk of false identification or discovery. To quantify this unknown risk, we propose a model-based framework for estimating the false discovery rate (FDR) among a set of identifications. Extending the traditional mixture modeling framework, our method incorporates both similarity score and experimental information in estimating the FDR. We apply these models to identification lists derived from across 548 samples of varying complexity and sample type (e.g., fungal species, standard mixtures, etc.), comparing their performance to that of the traditional Gaussian mixture model (GMM). Through simulation, we additionally assess the impact of reference library size on the accuracy of FDR estimates. In comparing the best performing model extensions to the GMM, our results indicate relative decreases in median absolute estimation error (MAE) ranging from 12% to 70%, based on comparisons of the median MAEs across all hit-lists. Results indicate that these relative performance improvements generally hold despite library size, however FDR estimation error typically worsens as the set of reference compounds diminishes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗