Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Modelling Capabilities”

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,117 records · Page 62

Modeling spatial patterns of longleaf pine needle dispersal using long-term data

Abstract Background Predicting patterns of fire behavior and effects in frequent fire forests relies on an understanding of fine-scale spatial patterns of available fuels. Leaf litter is a significant canopy-derived fine fuel in fire-maintained forests. Litter dispersal is dependent on foliage production, stand structure, and wind direction, but the relative importance of these factors is unknown. Results Using a 10-year litterfall dataset collected within eighteen 4-ha longleaf pine ( Pinus palustris Mill.) plots varying in canopy spatial pattern, we compared four spatially explicit models of annual needle litter dispersal: a model based only on basal area, an overstory abundance index (OAI) model, both isotropic and anisotropic litter kernel models, and a null model that assumed no spatial relationship. The best model was the anisotropic model (R 2 = 0.656) that incorporated tree size, location, and prevailing wind direction, followed by the isotropic model (R 2 = 0.612), basal area model (R 2 = 0.488), OAI model (R 2 = 0.416), and the null model (R 2 = 0.08). Conclusions As with previous studies, the predictive capability of the litter models was robust when internally validated with a subset of the original dataset (R 2 = 0.196–0.549); however, the models were less robust when challenged with an independent dataset (R 2 = 0.122–0.319) from novel forest stands. Our model validation underscores the need for rigorous tests with independent, external datasets to confirm the validity of litter dispersal models. These models can be used in the application of prescribed fire to estimate fuel distribution and loading, as well as aid in the fine tuning of fire behavior models to better understand fire outcomes across a range of forest canopy structures.

54 ENVIRONMENTAL SCIENCES↗

Diversifying Isoprenoid Platforms via Atypical Carbon Substrates and Non-model Microorganisms

Isoprenoid compounds are biologically ubiquitous, and their characteristic modularity has afforded products ranging from pharmaceuticals to biofuels. Isoprenoid production has been largely successful in Escherichia coli and Saccharomyces cerevisiae with metabolic engineering of the mevalonate (MVA) and methylerythritol phosphate (MEP) pathways coupled with the expression of heterologous terpene synthases. Yet conventional microbial chassis pose several major obstacles to successful commercialization including the affordability of sugar substrates at scale, precursor flux limitations, and intermediate feedback-inhibition. Now, recent studies have challenged typical isoprenoid paradigms by expanding the boundaries of terpene biosynthesis and using non-model organisms including those capable of metabolizing atypical C1 substrates. Conversely, investigations of non-model organisms have historically informed optimization in conventional microbes by tuning heterologous gene expression. Here, we review advances in isoprenoid biosynthesis with specific focus on the synergy between model and non-model organisms that may elevate the commercial viability of isoprenoid platforms by addressing the dichotomy between high titer production and inexpensive substrates.

59 BASIC BIOLOGICAL SCIENCES↗

A data-driven method to estimate the antiproton background in the Mu2e experiment

The Mu2e experiment at Fermilab will search for the Charged Lepton Flavour Violating (CLFV) process of coherent, neutrinoless µ− → e − conversion in the field of an aluminum nucleus. The expected signal is a monochromatic electron with the energy of 104.97 MeV, slightly below the muon rest mass. Observation of a CLFV process would provide unambiguous evidence for Beyond the Standard Model (BSM) physics. Mu2e is sensitive to a wide range of BSM models and has the capability to distinguish between them, guiding us towards the most accurate models. The key features of the Mu2e experiment are: (1) a high intensity pulsed negative muon beam with about 1010 stopped µ −/s, and (2) a sophisticated superconducting solenoid system with a gradient magnetic field to form and guide the intense muon beam to the target. The Mu2e physics data taking is expected to begin in 2027. For Run I, the expected 5σ discovery sensitivity is Rµe = 1.2 × 10−15, with a total expected background of 0.11 ± 0.03 events. In the absence of a signal, the expected upper limit is Rµe < 6.2 × 10−16 at 90% CL. The success of this experiment hinges on the accurate estimation of the background from various SM processes that could provide signal-like electrons. One of the background processes is antiprotons annihilating in the stopping target to produce signal like electrons through π0 → γγ decays followed by γ conversions, and π− → µ−ν¯ decays followed by µ− decay. It is a relatively small background with large uncertainty (100%) due to the lack of antiproton production cross section information for the Mu2e proton beam energy of 8 GeV. We have developed a novel methodology to estimate the antiproton background in-situ. This forms the main theme of the thesis. We observed that at Mu2e energies, antiproton annihilation in the stopping target is the only source of events with multiple, simultaneous particle trajectories. From Geant4 simulations, only about 0.2% of the simulated antiproton annihilation events have a signal-like electron. Meanwhile, ∼ 5% of events have multiple reconstructible particle tracks per event. Therefore, we have devised a methodology to reconstruct the multi-track events and estimate the antiproton background by exploiting the large ratio of the production rates of the two final states.

Chithirasreemadam, Namitha [Pisa U.] (ORCID:000000↗

Can a shock-induced phonon up-pumping model relate to impact sensitivity of molecular crystals, polymorphs and cocrystals?

Impact sensitivity engineering of high-energy molecular crystals requires accurate predictive models. For this purpose, the promising multi-phonon based approach is selected, assessing a bit more its strengths and weaknesses. Presently used with high-quality phonon calculations of 22 molecular crystals, using a physics-based criterion to determine the phonon bath extent, the resulting intrinsic shock sensitivity index (SSI) is compared to the most common marker of impact sensitivity, h 50 , as determined from drop-weight impact tests. Selecting a data subset from experiments performed under very similar conditions (2.5 kg hammer with grit and 30–40 mg samples), the model can predict h 50 values for mono-molecular crystals with very good accuracy, including the ability to discriminate the polymorphs of HMX and CL20. This very good agreement validates an initial indirect up-pumping mechanism occurring under these conditions, where the doorway modes also interact with the phonon bath. However, the phonon bath criterion for mono-molecular crystals does not transfer well to cocrystals. Owing to the vibrational coupling of the co-molecules, it seems a broader phonon bath should be considered. Additionally recalling experimental uncertainty and various experimental factors affecting h 50 values for a given compounds, we recommend that the density of the sample, granularity and morphology be systematically considered and reported along with measurements, which will in turn allow for more systematic data and predictive capabilities for sensitivity models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Declining carbohydrate solubilization with increasing solids loading during fermentation of cellulosic feedstocks by Clostridium thermocellum: documentation and diagnostic tests

Abstract Background For economically viable 2nd generation biofuels, processing of high solid lignocellulosic substrate concentrations is a necessity. The cellulolytic thermophilic anaerobe Clostridium thermocellum is one of the most effective biocatalysts for solubilization of carbohydrate harbored in lignocellulose. This study aims to document the solubilization performance of Clostridium thermocellum at increasing solids concentrations for two lignocellulosic feedstocks, corn stover and switchgrass, and explore potential effectors of solubilization performance. Results Monocultures of Clostridium thermocellum demonstrated high levels of carbohydrate solubilization for both unpretreated corn stover and switchgrass. However, fractional carbohydrate solubilization decreases with increasing solid loadings. Fermentation of model insoluble substrate (cellulose) in the presence of high solids lignocellulosic spent broth is temporarily affected but not model soluble substrate (cellobiose) fermentations. Mid-fermentation addition of cells ( C. thermocellum ) or model substrates did not significantly enhance overall corn stover solubilization loaded at 80 g/L, however cultures utilized the model substrates in the presence of high concentrations of corn stover. An increase in corn stover solubilization was observed when water was added, effectively diluting the solids concentration mid-fermentation. Introduction of a hemicellulose-utilizing coculture partner, Thermoanaerobacterium thermosaccharolyticum, increased the fractional carbohydrate solubilization at both high and low solid loadings. Residual solubilized carbohydrates diminished significantly in the presence of T. thermosaccharolyticum compared to monocultures of C. thermocellum, yet a small fraction of solubilized oligosaccharides of both C 5 and C 6 sugars remained unutilized. Conclusion Diminishing fractional carbohydrate solubilization with increasing substrate loading was observed for C. thermocellum -mediated solubilization and fermentation of unpretreated lignocellulose feedstocks. Results of experiments involving spent broth addition do not support a major role for inhibitors present in the liquid phase. Mid-fermentation addition experiments confirm that C. thermocellum and its enzymes remain capable of converting model substrates during the middle of high solids lignocellulose fermentation. An increase in fractional carbohydrate solubilization was made possible by (1) mid-fermentation solid loading dilutions and (2) coculturing C. thermocellum with T. thermosaccharolyticum , which ferments solubilized hemicellulose. Incomplete utilization of solubilized carbohydrates suggests that a small fraction of the carbohydrates is unaffected by the extracellular carbohydrate-active enzymes present in the culture.

09 BIOMASS FUELS↗

Novel Chalcopyrites for Advanced Photoelectrochemical Water Splitting

With the support of DoE’s EERE office, our team has established a unique tool-chest of capabilities, including theoretical modeling (Lawrence Livermore National Laboratory: LLNL), state-of-the-art synthesis (Hawaii Natural Energy Institute: HNEI, Stanford, and the National Renewable Energy Laboratory: NREL) and advanced materials and interfaces characterization (University of Nevada, Las Vegas: UNLV, and Lawrence Berkeley National Laboratory: LBNL), to accelerate the development of high efficiency and durable chalcopyrite materials for advanced photoelectrochemical (PEC) water splitting. Using this synergistic approach, we have successfully created new wide bandgap chalcopyrite photocathodes generating over 10 mA/cm 2 , developed innovative strategies to protect them from corrosion, and engineered novel integration methods to circumvent thin film materials mechanical, chemical and thermal incompatibility. In Task 1 “Modeling and synthesis of chalcopyrite photocathodes”, we expanded our library of wide bandgap chalcopyrites for PEC water splitting. With support from LLNL’s “Computational Materials Diagnostics and Optimization of PEC Devices”, LBNL’s “photophysical” and NREL’s “I-III-VI Compound Semiconductors for Water-Splitting” nodes, we investigated two new chalcopyrite candidates for PEC water splitting: Cu(In,Al)Se 2 and Cu(In,B)Se 2 . We also further developed ordered vacancy compounds, such as CuGa 3 Se 5 , with unprecedented durability during PEC waters splitting in acidic solutions. In Task 2 “Interfaces engineering for enhanced efficiency and durability”, we addressed both the non-ideal band-edge positions of chalcopyrites with respect to water redox potentials, as well as their chemical instability under PEC water splitting, with a buried-junctions approach. With help from NREL’s “High-Throughput Experimental Thin Film Combinatorial Capabilities” and “Corrosion Analysis of Materials” nodes, we engineered environmentally friendly n-type buffers, including Mn x Zn 1-x O, to adjust the chalcopyrite band-edge positions and achieved photovoltages as high as 925 mV. Also, we integrated non-precious catalytic-protecting layers, such as WO 3 , to enhance the water splitting long-term stability of chalcopyrite absorbers. Finally, in Task 3 “Hybrid photoelectrode device integration”, we proposed an innovative method to bond wide bandgap photocathodes onto narrow bandgap PV drivers at room temperature using conductive polymers. Our semi-monolithic approach addressed fundamental processing incompatibility issues, as both the photocathode and the PV driver are processed separately. Proof-of-concept whole-chalcopyrite tandems were obtained by consecutive exfoliation and transfer of fully integrated 1.85 eV CuGa 3 Se 5 and 1.13 eV CuInGaSe 2 stacks from their Mo/SLG substrates onto a new single FTO host substrate.

08 HYDROGEN↗

A general predictive methodology for fuel-mixture properties up to supercritical conditions

A predictive thermodynamic model is utilized for the calculation of fuel properties of oxymethylene dimethyl ethers (OME 3–4 ), surrogates for gasoline, diesel and aviation fuel, as well as alcohol blends with gasoline and diesel. The alcohols used for these blends are methanol, ethanol, propanol, butanol and pentanol; their mixing ratio ranges from 10 to 50% by volume. The model is based on the Perturbed-Chain Statistical Association Fluid Theory (PC-SAFT) equation of state (EoS) and Vapor Liquid Equilibrium (VLE) calculations at constant temperature, density and composition. The model includes the association term, with the assumption of two association sites (2B scheme), to enable the modeling of alcohols. The pure-component parameters are estimated based on the Group Contribution (GC) method of various sources, as well as a parametrization model specifically designed for the case of OME 3–4 . The results of the computational model for the density, vapor pressure and distillation curves at various conditions, including high-pressure, high-temperature (HPHT), are compared to experimental and computational data available in the literature. In the cases where no measurements are available for the surrogates, experimental data for the corresponding target fuel are used, taking into consideration the inherent deviation in properties between real and surrogate fuel. Overall, the results are in good agreement with the data from the literature, with the average deviation not exceeding 12% for temperature (Kelvin) on the distillation curves, 10% for density and 46% for vapor pressure and the general trend being captured successfully. The use of different pure component parameter estimation techniques can further improve the prediction quality in the cases of OME3–4 and the aviation fuel surrogate, especially for the vapor pressure, leading to an average deviation lower than 18%. These results demonstrate the predictive capabilities of the model, which extend to a wide range of fuel types and pressure/temperature conditions. Through this investigation, the present work aims to establish the limits of applicability of this thermodynamic property prediction methodology.

33 ADVANCED PROPULSION SYSTEMS↗

Requirements Description of GAMSOR and GAMSRC

This report presents the modeling and simulation capabilities of Argonne National Laboratory’s GAMSOR and GAMSRC code that is used in present reactor analysis activities. These capabilities will then be used to establish the set of verification tasks necessary to verify GAMSOR for use on commercial projects. A similar approach was taken for the PERSENT, REBUS and DIF3D software packages. The GAMSOR program is a modified version of DIF3D used to obtain the coupled neutron gamma flux distributions for a given reactor system. GAMSOR was first developed in the 1980s and its work has progressed in parallel with DIF3D throughout its various development phases. While GAMSOR is compatible with DIF3D-FD and DIF3D-Nodal, DIF3D-VARIANT is most commonly used today and will be the focus of the present verification effort. As covered in [1], GAMSOR was built to provide more accurate heating distributions by explicitly solving for both the gamma and neutron heating distributions. Without this consideration, gamma heat is normally assumed to be deposited locally at the site of gamma production. This is acceptable for modeling some phenomena because the total energy production is preserved, but it can lead to inaccurate approximations of fluence-dependent behaviors such as fuel performance, swelling, cladding embrittlement, radiation induced creep, and power peaking. The outputs from GAMSOR are contained in multiple intermediate and final output files: FIXSRC, VARSRC, GTFLUX, GHFLUX, PWDINT, NPDINT, and GPDINT. DIF3D normally produces a PWDINT file containing power density results, with the assumption that gamma rays are absorbed where they are emitted (i.e., energy deposition is local with no modeled transport). With GAMSOR, neutron power density (NPDINT) and gamma power density (GPDINT) are both produced, with their summed power densities replacing the results in PWDINT for total heating (i.e., recoverable) power density.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Enhancing occupant behavior representation for interoperability between building information modeling and building energy modeling

Building Performance Simulation (BPS) has been adopted as an essential tool for designing, operating, and retrofitting buildings to optimize energy efficiency throughout the building life cycle. The Green Building XML (gbXML) schema facilitates seamless data exchange between Building Information Modeling (BIM) and Building Energy Modeling (BEM) software tools. However, limited occupant behavior (OB) representation in BIM often leads to inconsistent and inaccurate energy simulation in BEM software. This paper presents 154 systematic enhancements to the existing occupant behavior XML (obXML) schema v1.3.4, initially developed for standardizing OB representation for BEM, to address existing limitations and improve interoperability with BIM models. The enhancements encompass improved integration with BIM models through extended building representations and system operations, expanded support for advanced OB models with additional environmental parameters and mathematical capabilities, and implementation of a standardized model documentation framework. To facilitate seamless data transformation between gbXML and obXML schemas, we developed a publicly available gb-obXML Schema Converter. Three case studies demonstrate the enhanced schema’s capabilities: representation of building information using a two-story office building model, documentation of a window operation behavior model, and validation of the schema converter’s functionality. The enhanced obXML schema v1.4 enables sophisticated modeling of occupant-building interactions while maintaining consistency with industry-standard BIM schemas. The standardized documentation framework facilitates reproducibility and knowledge sharing in the OB research community, while the schema converter automates the integration of building information into OB simulation workflows. These enhancements establish a foundation for more accurate building performance simulation by supporting sophisticated representation of occupant behavior within the BIM-to-BEM simulation workflows.

Chung, Jihoon↗

Advancing quantum simulations of the nuclear shell model with Gray-code–based resource-efficient protocols

Background: Some of the computational limitations in solving the nuclear many-body problem could be overcome by utilizing quantum computers. The nuclear shell-model calculations providing deeper insights into the properties of atomic nuclei are one such case with high demand for resources, as the size of the Hilbert space grows exponentially with the number of particles involved. Quantum algorithms are being developed to overcome these challenges and advance such calculations. Purpose: To develop quantum circuits for the nuclear shell-model, leveraging the capabilities of noisy intermediate-scale quantum (NISQ) devices. Here, we aim to minimize resource requirements (specifically in terms of qubits and gates) and strive to reduce the impact of noise by employing relevant mitigation techniques. Methods: We achieve noise resilience by designing an optimized Ansatz for the variational quantum eigensolver (VQE) based on Givens rotations and incorporating qubit-ADAPT-VQE in combination with variational quantum deflation (VQD) to compute ground and excited states, incorporating the zero-noise extrapolation mitigation technique. Furthermore, the qubit requirements are significantly reduced by mapping the basis states to qubits using Gray-code encoding and generalizing transformations of fermionic operators to efficiently represent many-body states. Results: By employing the resource-efficient protocols, we achieve the ground and excited state energy levels of 38 Ar and 6 Li with better accuracy. These energy levels are presented for noiseless simulations, noisy conditions, and after applying noise mitigation techniques. Results are compared for Jordan-Wigner and Gray-code encoding using VQE, qubit-ADAPT-VQE, and VQD. Conclusions: Our work highlights the potential of resource-efficient protocols to leverage the full potential of NISQ devices in scaling the nuclear shell model calculations, offering a pathway toward more complex quantum simulations in nuclear physics. This approach establishes a framework for studying other nuclear systems with improved quantum resource efficiency, marking a significant advancement in applying quantum computing to realistic nuclear physics applications.

Physics - Nuclear physics and radiation physics↗

Requirements Description of DASSH-F

This report reviews the modeling and simulation capabilities of Argonne National Laboratory’s DASSH code that is used in present reactor analysis activities. These capabilities will be used to establish the set of verification tasks necessary to verify DASSH for use on commercial projects. A similar approach was taken for the PERSENT, REBUS and DIF3D software packages. The DASSH program is a thermal analysis code designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. DASSH is a follow-on development to the SE2-ANL software and SUPERENERGY-2 software that it is based upon. DASSH was designed to account for both neutron and gamma heating and is inherently connected to the GAMSOR part of the ARC suite of fast reactor analysis software. SE2-ANL is a developed piece of software from the 1980s while DASSH is a modern implementation with notable improvements in geometry handling. The most important upgrade of DASSH relative to SE2-ANL is that it can analyze multiple time points in a single run where SE2-ANL can only treat a single time point. This allows the user to understand the impact of and search the flow distribution for the entire operational period of a reactor design considering pressure drop, peak coolant and fuel temperatures, and thermal striping. DASSH has three input paths that have to be verified. The first input path builds the geometry and power distribution based upon the DIF3D model but ignores the gamma heating aspects of the problem. The second input path also builds the geometry from the DIF3D model but it takes the neutron and gamma heating distributions from GAMSOR. The third input path is to take the geometry and power distribution directly from user input (i.e. not coupled to DIF3D or GAMSOR). DASSH also has many built in correlations for material properties along with a user defined specification of the fuel, structure, and coolant properties. There are correlations for flow split, mixing, pressure drop, and heat transfer coefficients (subchannel rather than a direct methodology). In total, verification of DASSH will require an extensive testing to cover all possible user features of the software.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Dry in-pile fracture test (DRIFT) for separate-effects validation of ceramic fuel fracture models

Fracture is an important component of nuclear fuel behavior, and significant efforts have been invested into developing fuel performance models that are capable of accurately representing fracture. Usable data on the process of fracture propagation in nuclear fuel under realistic service conditions are very limited. To address this need, a series of separate-effects experiments were developed and performed at Idaho National Laboratory's Transient Reactor Test (TREAT) facility. These experiments employ a heat sink to radially remove heat from the fuel in a manner that approximates the effect of coolant in an operating light-water reactor (LWR). The test holder for these experiments is known as the Dry In-pile Fracture Test (DRIFT). A series of experiments employing DRIFT and TREAT were performed to provide data on the extent and nature of fracture in fresh fuel at various points during a ramp to full power. Novel aspects of these experiments include the way they employ a heat sink to replicate steady-state LWR conditions, as well as the use of fiber optic sensors for in-reactor thermal instrumentation. Details on the development of this experiment, experimental conditions, and resulting data (including in situ thermal measurements and post-irradiation imaging of fracture) are provided in this work. LWR-equivalent powers ranging from 10 to 25 kW/m were tested using this apparatus. Cracking was visible at all power levels, with increasing cracking extent as the power level increased, although there was little difference in the cracking between the two highest-power tests, which had LWR-equivalent powers of 20 and 25 kW/m.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

pyTCR: A tropical cyclone rainfall model for python

pyTCR is a climatology software package developed in the Python programming language. It integrates the capabilities of several legacy physical models and increases computational efficiency to allow rapid estimation of tropical cyclone (TC) rainfall consistent with the large-scale environment. Specifically, pyTCR implements a horizontally distributed and vertically integrated model [Zhu et al., 2013] for simulating rainfall driven by TCs. Along storm tracks, rainfall is estimated by computing the cross-boundary-layer, upward water vapor transport caused by different mechanisms including frictional convergence, vortex stretching, large-scale baroclinic effect (i.e., wind shear), topographic forcing, and radiative cooling [Lu et al., 2018]. The package provides essential functionalities for modeling and interpreting spatio-temporal TC rainfall data. pyTCR requires a limited number of model input parameters, making it a convenient and useful tool for analyzing rainfall mechanisms driven by TCs. To sample rare (most intense) rainfall events that are often of great societal interest, pyTCR adapts and leverages outputs from a statistical-dynamical TC downscaling model [Lin et al., 2023] capable of rapidly generating a large number of synthetic TCs given a certain climate. As a result, pyTCR significantly reduces computational effort and improves the efficiency in capturing extreme TC rainfall events at the tail of the distributions from limited datasets. Furthermore, the TC downscaling model is forced entirely by large-scale environmental conditions from reanalysis data or coupled General Circulation Models (GCMs), simplifying the projection of TC-induced rainfall and wind speed under future climate using pyTCR. Finally, pyTCR can be coupled with hydrological and wind models to assess risks associated with independent and compound events (e.g., storm surges and freshwater flooding).

54 ENVIRONMENTAL SCIENCES↗

Cost Analysis of Annual and Monthly Net Zero Energy Performance for Multifamily Buildings in the United States

As the United States continues its progress toward sustainable construction, net zero building design is becoming an increasingly important and popular topic. Two definitions of net zero energy performance, monthly and annual, provide different levels of energy autonomy and efficiency within a building. This article analyzes the viability and incremental cost for two-, three-, and four-story multifamily apartment buildings to reach both annual and monthly net zero energy performance throughout four climate zones in the United States using baseline reference buildings that represent current construction practices. Building size plays a large role in determining the capability for a building model to reach annual or monthly net zero. Two-story buildings are capable of reaching annual net zero with an increase in construction cost of about 4.4–5.6%. Three-story buildings in warmer climates can achieve annual net zero with an increase of about 5.1%, but models in colder climates cannot reach annual net zero performance. Four-story buildings cannot achieve net zero, owing to insufficient building area for solar arrays to produce electricity. Monthly net zero energy performance was significantly harder to achieve, with the only case reaching this goal being the two-story Houston based model, at an added cost of 8.2%. Generally speaking, only small buildings in warm climates will be able to achieve monthly net zero without vastly oversizing photovoltaic systems and increasing costs without adequate payback.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Examining Graphite Degradation in Molten Salt Environments: A Chemical, Physical, and Material Analysis

Molten-salt reactors (MSRs) are Generation IV nuclear reactors that use liquid salt as a coolant and/or fuel. In several MSR designs, graphite serves as a moderator and/or reflector. However, due to limited experimental data and operational experience, our understanding of graphite behavior in molten salt environments remains incomplete. This report aims to identify the degradation mechanisms of nuclear graphite in MSRs, detail the mechanisms of each factor, and provide an initial assessment of their impact on the structural integrity of graphite components. This assessment is based on an extensive literature review and insights from subject matter experts. Furthermore, given the limited data, a modeling strategy using existing Grizzly software is proposed for a more thorough analysis where appropriate. Additionally, it presents mitigation strategies where applicable. The report covers physical degradation mechanisms such as infiltration, erosion, and abrasion, as well as chemical degradation mechanisms including fluorination, intercalation, corrosion, and oxidation. Molten salt can infiltrate the porous structure of graphite, leading to several detrimental effects. Entrapment of fissile products within the graphite pores can cause radiation damage and could pose challenges in the handling and disposal of contaminated components. The differential thermal expansion between the infiltrated salt and graphite, along with internal stress from pressurized molten salt and volumetric heating, can compromise the structural integrity of graphite. To mitigate these effects, employing ultra-fine graphite grades and applying sealants and coatings are effective strategies. A computational model based on coupled solid mechanics and heat transfer phenomena could be used to predict the internal stresses using Grizzly software. In pebble-bed MSRs, graphite fuel pebbles can cause abrasion against reactor components due to friction and wear. The severity of wear is influenced by various factors such as temperature, environment, and the presence of lubricants. Tribological studies reveal that higher temperatures and molten salt environments, such as FLiBe, significantly reduce wear rates compared to dry conditions. Additionally, the chemical composition of the salt can further optimize graphite's tribological performance. Long-term wear effects can be modeled by incorporating surface defects into the geometry and predict stresses under thermal and radiation effects using Grizzly software. Chemical degradation of graphite in a molten salt environment can occur through fluorination and intercalation. Fluorination can occur via replacement of hydrogen or oxygen atoms, or at the active sites, but does not cause structural degradation. Intercalation, on the other hand, can lead to exfoliation, where layers of graphite separate and peel away, damaging the graphite. Protective coatings can enhance graphite's resistance to intercalation. Graphite generally exhibits good chemical stability in molten salt environments, though it can corrode under specific conditions, particularly in the presence of impurities or oxidants. Studies have shown that protective coatings, such as plasma-sprayed partially stabilized zirconia (PSZ), can effectively prevent such degradation. Corrosion behavior varies significantly with different graphite grades and coating applications, underscoring the need for detailed studies on uncoated and coated graphite to understand and mitigate corrosion mechanisms in MSRs. Research indicates that the presence of oxidants and impurities can accelerate graphite degradation in molten salts, making it essential to explore acceptable impurity limits. Oxidation is another critical degradation mechanism, leading to weight loss and structural damage due to the formation of CO and CO 2 from the reaction of carbon atoms with oxygen. This process creates new porosity and compromises graphite's integrity. While extensive research on graphite oxidation has been conducted for gas-cooled reactors, studies specific to MSRs are limited. Findings from the coal industry suggest that molten alkali metal salts can significantly accelerate graphite oxidation, a hypothesis worth exploring for fluoride salts in MSRs. Understanding oxidation behavior in MSRs is vital for developing protective measures. The analysis of post-irradiated graphite from the MSRE experiment demonstrated exceptional chemical compatibility with molten fluoride salt, suggesting that the extent of chemical attack on graphite largely depends on the salt's infiltration capability. Therefore, the use of ultra-fine grade graphite could help mitigate chemical degradation effects. Existing oxidation modeling capabilities in Grizzly, which use reaction-diffusion equations to model graphite-air interactions, could be adapted to simulate the chemical degradation effects of graphite in molten salt environments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Variability of near-fault seismic risk to reinforced concrete buildings based on high-resolution physics-based ground motion simulations

Broadband physics-based simulated earthquake ground motions are utilized to characterize the regional-scale seismic risk to modern reinforced concrete (RC) structures. A highly dense dataset of ground motions covering a 100-km urn:x-wiley:00988847:media:eqe3413:eqe3413-math-0001 40-km domain was generated using kinematic fault rupture models with varying rupture characteristics to represent shallow crustal earthquakes and resolved up to frequencies of 5 Hz. Over 40,000 nonlinear response history simulations of short- and mid-rise RC special moment frame buildings were conducted using simulation models that are capable of representing nonlinear behavior and component deterioration effects. The spatial variability of structural risk within a single earthquake scenario and between different rupture scenarios is examined, and the regions of strongest directivity effects and highest structural demands are identified. The structural demands may vary by factor of up to 8.0 at very short distances from the fault, and the large dispersion in the demands decreases significantly beyond a distance of 15 km. The interstory drift and member rotation demands are substantially impacted by important features of the geological structure and the characteristics of the rupture scenarios, particularly the presence of localized high-slip regions. The frequency characteristics of the structures are found to play an important role in determining the effects of near-fault ground motions on the structural response and expected damage. The results of this study suggest that the simulated ground motions, particularly those generated using the hybrid rupture approach, may offer reasonable structural risk estimates for low-frequency structures and conservative estimates for high-frequency structures.

58 GEOSCIENCES↗

Uncovering acoustic signatures of pore formation in laser powder bed fusion

Abstract We present a machine learning workflow to discover signatures in acoustic measurements that can be utilized to create a low-dimensional model to accurately predict the location of keyhole pores formed during additive manufacturing processes. Acoustic measurements were sampled at 100 kHz during single-layer laser powder bed fusion (LPBF) experiments, and spatio-temporal registration of pore locations was obtained from post-build radiography. Power spectral density (PSD) estimates of the acoustic data were then decomposed using non-negative matrix factorization with custom $$\varvec{k}$$ k -means clustering (NMF $$\varvec{k}$$ k ) to learn the underlying spectral patterns associated with pore formation. NMF $$\varvec{k}$$ k returned a library of basis signals and matching coefficients to blindly construct a feature space based on the PSD estimates in an optimized fashion. Moreover, the NMF $$\varvec{k}$$ k decomposition led to the development of computationally inexpensive machine learning models which are capable of quickly and accurately identifying pore formation with classification accuracy of supervised and unsupervised label learning greater than 95% and 90%, respectively. The intrinsic data compression of NMF k , the relatively light computational cost of the machine learning workflow, and the high classification accuracy makes the proposed workflow an attractive candidate for edge computing toward in-situ keyhole pore prediction in LPBF.

36 MATERIALS SCIENCE↗

Spatiotemporal 4D Whole-cell Modeling of a Minimal Autotroph Reveals Central Carbon Metabolism Regulated Locally by Protein Megacomplexes via Post-translational Modifications under Light Disturbance

Photosynthetic microorganisms rely on multiple pathways in central carbon metabolism to adapt to fluctuating light and energy availability across diel cycles. Mechanistic insight into the regulatory dynamics of this adaptation requires integrating processes spanning disparate timescales, from rapid redox-dependent post-translational modifications (PTMs) to slower changes in protein expression and metabolic pathway usage. To address this complexity beyond genome-based inference and traditional modeling, we develop a whole-cell four-dimensional (3D + time) model of the marine cyanobacterium Prochlorococcus marinus MED4 that explicitly represents the spatial organization of enzymatic and molecular processes in central carbon metabolism under light perturbation. We employ a perturbation-based research design to experimentally generate time-series, multi-omics measurements that provide molecular descriptors and cryo-ET derived 3D segmented volumes as constraints for this dynamic 4D framework. The integration of experiments and modeling across defined light regimes enables quantitative validation of system-level responses and forecasting under distinct light disturbances. We test the hypothesis that light-dependent redox PTMs regulating the structural assembly of a protein megacomplex, the “dark complex,” modulate metabolic flux at a conserved regulatory node of the Calvin–Benson cycle (CBC) in cyanobacteria. Our model shows that subcellular spatial organization buffers rapid light-induced changes in thylakoid reaction rates, which are followed by redox-PTM-mediated sequestration or release of CBC enzymes in the dark complex, ultimately impacting carbon fixation dynamics within carboxysomes. Comparison with an equivalently parameterized well-mixed stochastic model demonstrates that post-translational regulation not only buffers transcriptional noise and diffusion-driven fluctuations but also stabilizes phenotypic outcomes, underscoring the importance of spatial heterogeneity in phenotypic robustness. This ability to probe adaptive, spatiotemporally resolved mechanisms in photosynthetic machinery and central carbon metabolism addresses a critical gap in genotype-to-phenotype inference and expands modeling and design capabilities for understudied or genetically intractable autotrophs such as P. marinus MED4.

Johnson, Connah G.↗