Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermomechanical analysis”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Correlating processing variables to material properties in recycled polypropylene: A data‐driven approach

Abstract Polypropylene (PP) is one of the most widely used plastics, yet its recycling remains limited, with less than 1% of solid waste PP being reprocessed. Mechanical recycling through extrusion is the most practical method, but inconsistent reprocessing conditions introduce variability in material properties. While temperature, screw speed, and residence time influence the thermomechanical stress applied during reprocessing, there are no standardized guidelines for optimizing these parameters. This study examines how these factors shape the properties of recycled PP, using conditions designed to mimic post‐industrial recycled (PIR) scrap. Residence time was measured using colorimetric tracking and correlated with molecular weight, viscosity, and mechanical properties over multiple extrusion cycles. Data‐driven modeling, including response surface methodology, support vector machines, and artificial neural networks, identified processing temperature as the dominant factor in material degradation, followed by residence time. Mechanical properties remained stable, while viscosity decreased predictably with increasing residence time. By linking reprocessing conditions to property evolution, this study provides a method to optimize processing parameters and reduce variability in recycled PP. These findings help manufacturers improve process control, making recycled PP more predictable for reuse in manufacturing. Highlights Study of PIR‐quality PP without additives or compatibilizers. Residence time analysis shows processing temperature drives PP property changes. Mark‐Houwink enables quick molecular weight checks for quality control. Models predict mechanical and rheological shifts in reprocessing. Optimized processing parameters minimize property degradation in recycling.

Estela‐García, John E. [Polymer Engineering Center↗
DOE OSTI2025

Microstructure and bonding between calcium aluminate cement‐containing gahnite–alumina matrix and refractory aggregates

Calcium aluminate cement enhances the thermomechanical properties of refractory castables through the formation of acicular calcium hexaluminate (CaO·6Al 2 O 3 ), Ca 2 Mg 2 Al 28 O 46 (CAM-I), and CaMg 2 Al 16 O 27 (CAM-II) phases in MgO- or MgAl 2 O 4 -containing castables. The compatibility of CA 6 with gahnite (ZnAl 2 O 4 ), and acicular Ca 2 Zn 2 Al 28 O 46 (CAZ-I) and CaZn 2 Al 16 O 27 (CAZ-II) phases formation have been previously reported. Here, in this work, the interaction between a CAC binder containing ZnAl 2 O 4 -Al 2 O 3 matrix with commonly used refractory aggregates such as tabular alumina (TA), alumina-rich (AR90, AR78) and stoichiometric (SM72) MgAl 2 O 4 spinels, and fused and dead-burned magnesia (FM, DBM, respectively) were investigated at 1650°C for 5 h. Microstructural analysis, using digital microscopy, scanning electron microscopy, and energy dispersive spectroscopy, revealed the formation of acicular CaZn 0.18 Al 11.82 O 18.91 , CAZ-I and CAZ-II grains, and strong interfacial bonding between the matrix and TA and spinel aggregates. FM and DBM were found to debond from the matrix. Thick interface layers were observed between the matrix and all the aggregates but TA. Null hypothesis significance testing (NHST) shows that the difference in the number of acicular grains between the interface zone and the bulk matrix (Z) is statistically significant for AR90/Z, SM72/Z, FM/Z, and DBM/Z interfaces, but not significant for TA/Z and AR78/Z. The role of the aggregates’ chemistry on the interfacial bonding and microstructure evolution is discussed.

Ramteke, Rajat Durgesh [Univ. of Alabama, Birmingh↗
DOE OSTI2023

Sensitivity of UO 2 fuel performance to microstructural evolutions driven by dilute additives

Use of dilute additives to nuclear fuel is being considered to increase the security of commercial fuel management through traceability of fabricated fuel elements. Taggants, as additives are denoted when included for traceability purposes, may also improve fuel performance, as demonstrated in Cr-containing uranium dioxide as described in the literature, and they may also improve fuel safety. In fact, studies have shown that some additives affect fuel material properties such as grain size and density after sintering. Given the possible range of elements that could be used as additives, the impact of such fuel property variations on the fuel’s thermomechanical behavior becomes relevant. These effects can be evaluated through a sensitivity study of standard fuel models to analyze changes in these properties using a fuel performance code. In this work, the BISON code is being used to investigate these effects through a 2D axisymmetric model of smeared UO 2 fuel pellets and ZIRLO® cladding under realistic pressurized water reactor core irradiation conditions. Here, randomly sampled densities and grain sizes within specified ranges are used as input parameters in the simulations, and several fuel model-related outputs are evaluated. The thermomechanical response of the cladding is also addressed in this study. The simultaneous variation of both input parameters offers a more comprehensive path to identify key sensitivities. Outputs explored include temperature, fission gas release, creep, and radial stress. Results show that although most of these outputs are sensitive to grain size to a certain extent, density mainly affects fuel temperature and elastic strain. Furthermore, sensitivities can vary depending on the radial position within the fuel pellet.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Towards Commercialization of Low-Cost, Crack-Tolerant, Screen-Printable Metallization by Full-Size Module Testing and Field Characterization

This project is motivated by the need to develop a materials engineering solution to reduce solar module degradation caused by cell cracks. The cell-crack-induced power loss is a long-term degradation mechanism and one of the main causes of solar panel field failures. Cell cracks can occur during module fabrication, transportation, installation, and long-term operation due to thermomechanical stressors. Our team’s internal estimation – based on the national weather pattern, frequency of severe weather, and 39- GW asset survey by Heliolytics – reveals that the cell-crack-induced module degradation translates to lost revenues >$17B over an average 5-year period for solar farm owners and to reduced reliability (<25 years panel lifetime) for solar energy consumers. In response to this challenge, the prime recipient offers a metal matrix composite (MMC) silver paste that is tailor-engineered for screen-printing gridlines and busbars, which serve as the electrical contacts on solar cells. We formulate our MMC paste by adding low-cost (∼0.02¢/WDC for research grade), surface-engineered carbon nanotubes to commercially available silver paste. The MMC paste offers a drop-in, cost-effective solution to cell cracks for solar cell and module manufacturers. The main goal of this project was to conduct field-relevant, module-level analysis and qualification of MMC metallization.

14 SOLAR ENERGY↗
DOE OSTI2020

Drivers for the cracking of multilayer polyamide‐based backsheets in field photovoltaic modules: In‐depth degradation mapping analysis

Abstract There is a lack of understanding on the root cause of cracking of photovoltaic (PV) backsheets due to the challenge of multilayer characterization and the complicated failure modes at the submodule level. In this work, the in‐depth degradation mapping of field‐exposed polyamide‐based (PA‐based) PV module backsheets was studied, with the major focus on the identification of underlying drivers for the through cracking (in between the solar cells). PV modules were retrieved from five different locations, comprising a variety of climates, including humid subtropical, hot‐summer Mediterranean, tropical savanna climate, and hot arid. A suite of microscale cross‐sectional characterizations, including chemical changes, fluorescence intensity, and modulus as a function of distance from the air surface of the backsheet, was performed. Results showed more advanced signs of degradation of the inner layer than the outer layer in the cracking region. Increases in the modulus were identified as the major indicator for the cracking. Moreover, a rudimentary test by immersion in acetic acid, which forms during photodegradation of the ethylene‐vinyl acetate copolymer (EVA) encapsulant, showed the first‐time direct evidence that acetic acid can largely accelerate the chemical degradation and facilitate the cracking of PA inner layer. This study suggests that the field cracking of PA‐based backsheet can be attributed to the combined effects of chemical degradation and physical reorganization (chemi‐crystallization) under cyclic thermomechanical stresses. Therefore, it is important to consider effects of local microclimates and interlayer infection to understand the heterogeneous nature of backsheet failures.

Lyu, Yadong↗
DOE OSTI2025

Ultra-high temperature testing and performance of L-PBF C103

Additive Manufacturing (AM) of refractory alloys is gaining traction as a materials processing route for components subject to extreme temperature environments. Due to the low oxidation resistance of refractory alloys, novel methods for evaluating their elevated temperature performance must be developed. In this work, a Gleeble® 3800 thermomechanical load frame was modified to evaluate the mechanical properties of laser powder bed fusion (L-PBF) consolidated niobium alloy C103 ranging from room temperature (RT) to 1400 °C. The fixturing and sample geometry were designed to accommodate Joule heating and prevent damage to the test chamber. Oxidation of the samples was minimized via testing in vacuum level of 1E-5 Torr. Ultimate tensile strength (UTS), yield strength (YS), elongation, and strain-hardening behavior were determined as a function of temperature. L-PBF C103 presented an average UTS of ∼650 MPa and over 25 % elongation at RT. Above RT, the UTS and YS dropped then leveled off from 500 °C to 1000 °C with values ranging from ∼400 MPa to ∼460 MPa, which is consistent with dynamic strain aging observed in this class of alloys. The strength rapidly declined after 1200 °C to ∼150 MPa at 1400 °C. Fractography indicated ductile fractures for the C103 at all test temperatures, and Electron Backscatter Diffraction (EBSD) analysis revealed a textured microstructure and the presence of dynamic recrystallization within the necked region of the sample tested at 1400 °C.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Wastewater Recycling Using a Hygroscopic Cooling System

This project by the Energy & Environmental Research Center (EERC), Baltimore Aircoil Company (BAC), and Great River Energy (GRE) evaluated the concept of recycling wastewater at a coal-fired power plant using a hygroscopic cooling system, which is an evaporative cooling technology analogous to conventional cooling towers, except that sparingly soluble dissolved solids are precipitated and removed as waste solids instead of purging them with a liquid blowdown stream. This technology can maximize the use of plant makeup water by obtaining useful evaporative cooling from wastewater while minimizing the volume of wastes needing disposal. Experimental activities were conducted in two phases, a laboratory-based evaluation of the properties of wastewater from the host site, GRE’s Coal Creek Station near Underwood, North Dakota, and a field test of a small pilot hygroscopic cooling system at the host site power plant. Findings from the laboratory study informed the design of the pilot system and the system’s field test performance served as the basis for a techno-economic analysis (TEA) of the hygroscopic recycling concept. At the preferred operating conditions identified during the TEA, the wet-bulb approach temperature of the tower was 7.3°C (13°F) and the volume of blowdown produced by the plant was reduced to 5.4% of its incoming volume. Waste solids produced during field testing were classified as nonhazardous waste based on the measured hazardous element content and evaluation of their leaching potential. However, to qualify as a solid for landfill disposal i.e., as determined by the U.S. Environmental Protection Agency’s paint filter test, it appears that a dewatering step beyond hydrocyclone separation is needed. The baseline levelized cost of wastewater disposal (LCWD) for hygroscopic wastewater recycling was estimated to be $\$ $3.69–$\$ $3.72 per m3 of plant blowdown. Capital cost was estimated to contribute over 54% to the LCWD, and parameters that impact capital cost such as the heat exchange coil material of construction and the tower’s wet-bulb approach temperature were identified as having the greatest impact on overall LCWD. A LCWD estimate prepared for the same application but using thermomechanical brine evaporation was almost 40% higher than that calculated for hygroscopic cooling, despite recovering distilled-quality water for reuse, while the LCWD for disposal-only, deep well injection was estimated to be 30% lower compared to hygroscopic wastewater recycling.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Adhesion in thermomechanically processed seaweed-lignocellulosic composite materials

The increasing concerns associated with petroleum-derived polymers motivate the development of sustainable, renewably sourced alternatives. In ubiquitous applications such as structural materials for infrastructure, the built environment as well as packaging, where natural materials such as wood are used, we rely on nonrenewable and nondegradable polymers to serve as adhesives. In wood panels, such as medium density fiberboards (MDFs), formaldehyde-based resins are predominantly used to bond wood fibers and to provide strength to the materials. To further mitigate the environmental impact of construction materials, more sustainable adhesives need to be investigated. Here, in this article, we introduce Ulva seaweed as an adhesive to enable cohesion and strength in hot-pressed wood panels. Upon hot-pressing, powdered Ulva flows in between the wood particles, generating a matrix, which provides strong binding. We show that the flexural strength of Ulva-bonded wood biocomposites increases with increasing Ulva concentrations. At an Ulva concentration of 40 wt%, our composites reach an average elastic modulus of 6.1 GPa, and flexural strength of 38.2 MPa (compared to 4.7 GPa and 22.6 MPa, respectively, for pure wood compressed at the same pressing conditions). To highlight the bonding mechanisms, we performed infrared and x-ray photoelectron spectroscopy and identified indications of fatty acid mobility during hot-pressing. In addition, we demonstrate that the presence of Ulva improves other properties of the composites such as water resistance and fame retardancy. Ulva is also shown to behave as an excellent adhesive agent between two prepressed beams. Finally, we perform an in-depth analysis of the environmental impact of wood-Ulva biocomposites.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First considerations on the supporting structures of FCC-ee booster and collider in the arc regions

In 2022, the FCC Feasibility Study management mandated a working group to analyse the best configuration of the FCC-ee tunnel in the arc regions, in view of the construction of a mock-up of the arc half-cell. One of the main and most challenging goals of the study, named FCC-ee Arc Half-Cell Mock-up Project, was to perform a preliminary investigation on the principles of supporting the Short-Straight Sections and dipoles of the half-cells, both for the booster and for the collider machines. This is an important input needed for the choice of the best configuration of the relative placement of the booster with respect to the collider. The structural stiffness, mass and stability of the supporting structures must be optimized to minimize the vibrations transmitted/transferred to the magnetic system of the accelerators by elements such as pumps, water cooling system, beam thermomechanical stresses, powering elements, etc. To perform the study, tools such as CAD software, FEM and analytical techniques were employed. This paper summarizes the preliminary design concepts and the results of the simulations performed.

Accelerator Subsystems and Technologies↗
DOE OSTI2024

Stress Relaxation Cracking in 347H Stainless Steel Arc Welds: Susceptibility Evaluation of Heat-Affected Zone

Stress relaxation cracking (SRC) is considered one of the major failure mechanisms for 347H stainless steel welds at elevated service temperatures or during post weld heat treatment (PWHT), especially within the heat-affected zone (HAZ). This work focuses on the characterization of SRC susceptibility within 347H physically simulated arc welded HAZ at elevated temperatures. A four-step SRC thermomechanical test in combination with finite element modeling (FEM) of the welding and testing processes is developed to establish a susceptibility map for HAZ. The test first runs a thermal cycle with three different peak temperatures (1335, 1275, and 1150 °C) to duplicate representative HAZ subzone microstructures, followed by time-to-failure examination under a variety of pre-stress (260–600 MPa) and pre-strain conditions (0.03–0.19) as a function of reheat temperatures between 750 and 1050 °C. With the aid of FEM, SRC susceptibility maps are generated to identify the threshold stress, plastic strain, and creep strain as a function of test temperature. It was found out that HAZ subzone with a lower peak temperature (1150 °C) appears to be slightly less susceptible to SRC than the other two subzones that experienced higher peak temperatures. Generally, time-to-fracture reduces with increasing initially applied stress and strain for all test temperatures. The pre-stress thresholds decrease from about 500 to 330 MPa as the testing temperature increases from 800 to 1050 °C, while the corresponding initial plastic strain thresholds reduces from 0.15 to 0.06. The SRC susceptibility was also evaluated through the Larson–Miller Parameter (LMP) analysis as a function of plastic strain, initial stress and starting stress upon reaching the testing temperature, respectively. The 1050 °C test with a high pre-applied strain (0.1) exhibits an extremely short time to failure (t = 3 s) that lies outside the general trend in LMP analysis. Additionally, it was identified that a plastic strain above 0.07 is identified to significantly reduce the bulk creep strain tolerance to fracture and therefore increases SRC susceptibility. Hardness measurement and fractography analysis indicated that the strain aging of niobium carbonitrides and other potential phases in conjunction with intergranular precipitates contributes to an increase in microhardness and increased intergranular cracking susceptibility.

14 SOLAR ENERGY↗
DOE OSTI2025

Step-loaded creep testing of Zircaloy-4 cladding at higher temperatures in the α-phase

A refined understanding of zirconium-based cladding thermomechanical performance during rapid transients is essential for enhancing the safety and operation of light-water reactors. Traditional models for zirconium alloys under accident conditions generally assume that creep dominates fuel cladding performance. Here, these historic models have largely remained unchanged and serve as the basis for safety criteria development. As the U.S. nuclear industry pursues higher burnup levels, the increased release of fission gases during transients raises the risk of cladding failure in the low-temperature hcp α-phase, making the fidelity of these models of greater importance. Creep testing was conducted from 550–700°C with 25–120 MPa applied hoop stresses to investigate Zircaloy-4 deformation at accident-relevant temperatures in the α-phase. Step-loading was employed to capture creep behavior across a wide stress range from a single sample. The stress-strain rate data at higher temperatures (650 and 700°C) were well-described by isotropic versions of the Erbacher and Kaddour models, while the lower temperature data (550 and 600°C) were underpredicted by both anisotropic and isotropic model variants. Greater strain rates during the initial loading step at 650 and 700°C were attributed to recrystallization and grain growth of sub-micron crystallites. Yet, texture analysis revealed the basal split texture remained after testing. These observations produced results suggesting Zircaloy-4 claddings experience higher creep rates across the α-phase than previously thought, possibly related to dynamic anisotropy due to temperature dependent activation of deformation mechanisms, effects of biaxial loading, and variation in material condition between the current testing used in previous model development.

36 MATERIALS SCIENCE↗
DOE OSTI2021

NRC Multiphysics Analysis Capability Deployment FY21: Part 3

This report details the progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Support for Advanced Non-Light Water Reactors.” The deliverables completed for this report are: Deliverable 1c: the capability to model gas mixtures was added to Pronghorn. A test problem mimicking the conditions achieved in a depressurized loss of forced cooling (DLOFC) event was solved with both RELAP-5 and Pronghorn. Pronghorn employed a finite vol ume method with the Kurganov-Tadmor discretization. The comparison between the mass fraction spatial profiles computed with RELAP-5 and Pronghorn clearly shows the presence of numerical artifacts (i.e., overly diffusive behavior at low Mach numbers). We confirmed that the problem disappears at higher Mach numbers. We recommend future work on the implementation of a low Mach finite volume formulation to better treat low Mach number problems. Deliverable 2a: we demonstrated two approaches to model the radiation/conduction/natural convection heat transfer across a stagnant gas for the PBMR-400 design using Pronghorn. The first approach is based on the net radiation method, which relies on the computation of view factors with the Multiphysics Object-Oriented Simulation Environment (MOOSE) ray tracing capability. The second method is a traditional thermal resistance approach. The test problems include both 2D and 3D geometries. In all cases, the results show very good agreement during a DLOFC transient. This confirms that the faster thermal resistance method produces solutions that are equivalent to the net radiation method for this geometry. Deliverable 3d: we demonstrated the use of the advection kernel for the delayed neutron precursor equation in Griffin with a 2D MSFR model. The results appear physical but further verification is recommended. We also recommend the addition of conjugate heat transfer to compute the temperatures and model the thermomechanic behavior of the reflectors and other structures. Significant memory and performance issues were encountered in the 3D axisymmetric model. Future work is recommended in this area. Task 8g: this task allows multidimensional MOOSE applications to be coupled to system codes (RELAP-7 and SAM). We implemented a faster multiphysics iteration coupling algorithm, which provides an overall 6× acceleration of the 3D-1D coupling of the core multidi- mensional fluid flow solver and the 1D primary and secondary loop model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Topology-Informed Design Rules for Deconstructable Thermoset Copolymer Networks

Existing models of thermoset deconstruction facilitated by incorporating cleavable comonomers rely on a mean-field reverse gel point paradigm, which predicts network dissolution once cleavable bonds reach a critical stoichiometric threshold, but does not account for where those bonds reside within the network architecture. Using reactive coarse-grained molecular dynamics simulations coupled with graph-theoretic analysis, we extend this stoichiometric picture to show that deconstructability is governed by the curing-imprinted network topology rather than stoichiometry alone. This topological organization is hierarchical: at the local scale, the elastic effectiveness of cross-link junctions determines which cross-links constitute the load-bearing scaffold; at the mesoscale, the cross-linking rate kinetically templates that scaffold into topologically modular communities─densely cross-linked clusters connected by sparse bridging strands that sustain network connectivity. Using betweenness centrality to identify nodes that disproportionately lie on intercommunity shortest paths, we demonstrate that effective deconstruction of the network into macromolecular fragments requires cleavable comonomers to intercept these high-centrality bridging strands. We further find that under uniform, disassortative comonomer incorporation, this topological requirement provides a mechanistic basis for extending the reverse gel point to incorporate network topology. We also show that modularity imposes a fundamental limit on fragment uniformity that persists even when the centrality requirement is met. Finally, we demonstrate that chain stiffness provides a nearly independent lever to suppress mechanically redundant cross-links and raise the glass transition temperature without significantly altering the deconstruction outcome. Together, these findings reframe the thermoset design space around network topology and provide actionable guidelines for engineering thermoset copolymers with predictable deconstructability and targeted thermomechanical performance.

coarse-grained molecular dynamics↗
DOE OSTI2023

Discovery of multi-functional polyimides through high-throughput screening using explainable machine learning

Polyimides have been widely used in modern industries because of their excellent mechanical and thermal properties, e.g., high-temperature fuel cells, displays, and aerospace composites. However, it usually takes decades of experimental efforts to develop a successful product. Aiming to expedite the discovery of high-performance polyimides, we utilize computational methods of machine learning (ML) and molecular dynamics (MD) simulations. Our study provides compelling evidence for the effectiveness of a data-driven approach in discovering novel polyimides. We first build a comprehensive library of more than 8 million hypothetical polyimides based on the polycondensation of existing dianhydride and diamine/diisocyanate molecules. Then we establish multiple ML models for the thermal and mechanical properties of polyimides based on their experimentally reported values, including glass transition temperature, Young’s modulus, and tensile yield strength. The obtained ML models demonstrate excellent predictive performance in identifying the key chemical substructures influencing the thermal and mechanical properties of polyimides. The use of explainable machine learning describes the effect of chemical substructures on individual properties, from which human experts can understand the cause of the ML model decision. Applying the well-trained ML models, we obtain property predictions of the 8 million hypothetical polyimides. Then, we screen the whole hypothetical dataset and identify three (3) best-performing novel polyimides that have better-combined properties than existing ones through Pareto frontier analysis. For an easy query of the discovered high-performing polyimides, we also create an online platform https://polyimide-explorer.herokuapp.com/ that embeds the developed ML model with interactive visualization. Furthermore, we validate the ML predictions through all-atom MD simulations and examine their synthesizability. The MD simulations are in good agreement with the ML predictions and the three novel polyimides are predicted to be easy to synthesize via Schuffenhauer’s synthetic accessibility score. Following the proposed ML guidance, we successfully synthesized a novel polyimide and the experimentally obtained high glass transition/thermal decomposition temperature demonstrated its excellent thermal stability. Here our study demonstrates an efficient way to expedite the discovery of novel polymers using ML prediction and MD validation. The high-throughput screening of a large computational dataset can serve as a general approach for new material discovery in other polymeric material exploration problems, such as organic photovoltaics, polymer membranes, and dielectrics.

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

Bayesian Calibration of Nuclear Graphite Property Models Accounting for Model Inadequacy and Impacts on Component Performance

Nuclear-grade structural graphite is extensively utilized in the core designs of various advanced nuclear reactors. In the reactor environment, graphite is subjected to prolonged exposure to extreme conditions, including high temperatures, radiation, and potentially molten salt and oxygen. Such exposure can induce several degradation mechanisms in graphite, including nonuniform volumetric strains caused by irradiation and thermal expansion, leading to stresses that may compromise the performance of graphite components. Assessing component integrity requires accurate models of graphite's thermomechanical response. This report documents the Bayesian calibration of thermomechanical properties for nuclear-grade graphite and their application to graphite component modeling and simulation using the Grizzly code. As part of this work, uncertainty-quantified models were developed for the elastic modulus, coefficient of thermal expansion, irradiation-induced dimensional change, and irradiation-induced creep for graphite grades IG-110, NBG-18, NBG-17, PCEA, and 2114. Using a hierarchical Bayesian approach, multiple experimental data sources were combined to develop Gaussian process models for the properties. Using the Kennedy O'Hagan framework, the uncertainties due to inadequacies in the model and the inherent spread in the experimental data were quantified for three different models. These uncertainty-quantified models, with a model-form correction, were subsequently applied to a coupled-physics simulation of representative graphite components, revealing that the uncertainties have a large impact on the components' deformation.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Three-Dimensional Heat Flux and Thermal Analysis of Angled Tungsten Samples on DIII-D

ITER-grade tungsten and dispersoid-strengthened tungsten samples with the top surface angled at ~15° towards the incident plasma flux were exposed to 9 H-mode discharges with edge-localized modes (ELMs) in the lower divertor of DIII-D tokamak using the Divertor Material Evaluation System (DiMES). Surface damage included cracking and flaking of material on the two samples farthest away from the plasma strike point, and significant melting of the two samples closest to the strike point. Heat flux and thermal analysis tools new to DIII-D have been applied to better understand this material response and to help optimize the exposure conditions for future experiments. SMITER field-line tracing simulations based on IRTV data and EFIT equilibria estimate an average inter-ELM perpendicular heat flux, 𝑞⊥,𝑖nter−𝐸LM , on the angled surfaces of 10.1 – 19.6 MW/m² for a majority of the 9 discharges, increasing to 15.6 – 24.5 MW/m² for the single, higher-power shot where samples melted. Fast camera data showed shallow intra-ELM melting and re-solidification, which transitioned to bulk inter-ELM melting with melt motion in the 𝐽⃗ 𝑥 𝐵⃗ direction. About 50% of the protruding volume of the most affected sample was displaced via melt-motion. SIERRA thermal modeling software was able to reproduce an onset time of melting consistent with fast camera data and final sample conditions, within < 200 ms. Maximum surface temperatures of 3122 K and 2787 K are estimated for the samples farthest away from the strike point, while the closest samples achieve melting at 4067 ms and 4750 ms into the ~5000 ms plasma exposure. A +10% increase in both the SMITER 𝑞⊥,𝑖nter−𝐸LM calculations and the estimated ELM heat loads 𝑞⊥, 𝐸LM was required to achieve this result, which is within the uncertainty of the diagnostic data but likely accounts for non-ideal geometry effects plus other physics uncertainties not included in this first iteration of modeling. This work provided valuable estimates of the 3D temperature evolution to help better understand the observed surface morphology and internal recrystallization of samples, which are discussed in detail in a complementary manuscript [1]. Benchmarking efforts with more diagnosed DIII-D experiments are underway to further refine the SMITER and SIERRA models for DiMES. Future use of these tools will enable researchers to precisely target heat flux exposure conditions in DIII-D to test, but not exceed, the thermomechanical limitations of novel plasma-facing materials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗