Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled Simulation”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Size- and temperature-dependent magnetization of iron nanoclusters

In this paper, the magnetic behavior of bcc iron nanoclusters, with diameters between 2 and 8 nm, is investigated by means of spin dynamics simulations coupled to molecular dynamics, using a distance-dependent exchange interaction. Finite-size effects in the total magnetization as well as the influence of the free surface and the surface/core proportion of the nanoclusters are analyzed in detail for a wide temperature range, going beyond the cluster and bulk Curie temperatures. Comparison is made with experimental data and with theoretical models based on the mean-field Ising model adapted to small clusters, and taking into account the influence of low coordinated spins at free surfaces. Our results for the temperature dependence of the average magnetization per atom M (T), including the thermalization of the transnational lattice degrees of freedom, are in very good agreement with available experimental measurements on small Fe nanoclusters. In contrast, significant discrepancies with experiment are observed if the translational degrees of freedom are artificially frozen. The finite-size effects on M (T) are found to be particularly important near the cluster Curie temperature. Simulated magnetization above the Curie temperature scales with cluster size as predicted by models assuming short-range magnetic ordering. Analytical approximations to the magnetization as a function of temperature and size are proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Development of Multiresolution Capabilities for the Holistic Energy Resource Optimization Network (HERON) tool A progress update

INL researchers work on technoeconomic analyses for integrated energy systems (IES) using the Framework for Optimization of ResourCes and Economics (FORCE). Within FORCE, researchers use the Holistic Energy Resource Optimization Network (HERON) tool to conduct optimization of grid portfolios under uncertain market conditions. These optimizations determine optimal capacities for all IES components and strategies for resource dispatch which maximize some economic metric (e.g., net present value). Resource dispatch occurs on finer timescales (typically hours) and thus are asked to respond to a given time series (e.g. hourly load demand profiles for a grid, or pre-determined electricity prices). Volatile and complex bidding dynamics as well as poorly forecasted weather events within deregulated markets add uncertainty to the time series; FORCE can address this uncertainty by training a reduced order model on historical time series and generate unique synthetic time series which represent individual scenarios or realizations of the market. The IES configuration can be simulated under these different sampled realizations and a stochastic optimization is conducted which optimizes the expected value of the desired economic metric. The training of a synthetic time series generator is limited by the chosen time resolution; dynamics can occur on different time scales. Seasonal demand trends can dominate faster dynamical events (such as power outages from certain sectors or severe weather events) which might not get captured correctly by the trained model. In this report, we investigate different ways of addressing the training and generation of time series on multiple time scales using three main algorithms: wavelet decomposition, dynamic mode decomposition, and generative adversarial networks for time series. We demonstrate a time series analysis that yields information on not just the frequency space but also temporal space: where a fast Fourier transform can provide what frequencies dominate, the new algorithms can provide when the frequencies dominate as well. These analyses can help improve IES optimization by allowing researchers to couple simulations at different timescales when it is most needed - seasonal, day-ahead, and real time optimization - with greater computational efficiency. Future work will include implementation of a subset of the proposed algorithms into the FORCE toolset and application of these analyses into multiple timescale optimization.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Embedded symmetric positive semi-definite machine-learned elements for reduced-order modeling in finite-element simulations with application to threaded fasteners

Here, we present a machine-learning strategy for finite element analysis of solid mechanics wherein we replace complex portions of a computational domain with a data-driven surrogate. In the proposed strategy, we decompose a computational domain into an “outer” coarse-scale domain that we resolve using a finite element method (FEM) and an “inner” fine-scale domain. We then develop a machine-learned (ML) model for the impact of the inner domain on the outer domain. In essence, for solid mechanics, our machine-learned surrogate performs static condensation of the inner domain degrees of freedom. This is achieved by learning the map from displacements on the inner-outer domain interface boundary to forces contributed by the inner domain to the outer domain on the same interface boundary. We consider two such mappings, one that directly maps from displacements to forces without constraints, and one that maps from displacements to forces by virtue of learning a symmetric positive semi-definite (SPSD) stiffness matrix. We demonstrate, in a simplified setting, that learning an SPSD stiffness matrix results in a coarse-scale problem that is well-posed with a unique solution. We present numerical experiments on several exemplars, ranging from finite deformations of a cube to finite deformations with contact of a fastener-bushing geometry. We demonstrate that enforcing an SPSD stiffness matrix drastically improves the robustness and accuracy of FEM–ML coupled simulations, and that the resulting methods can accurately characterize out-of-sample loading configurations with significant speedups over the standard FEM simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Towards performing high-resolution inelastic X-ray scattering measurements at hard X-ray free-electron lasers coupled with energetic laser drivers

High-resolution inelastic X-ray scattering is an established technique in the synchrotron community, used to investigate collective low-frequency responses of materials. When fielded at hard X-ray free-electron lasers (XFELs) and combined with high-intensity laser drivers, it becomes a promising technique for investigating matter at high temperatures and high pressures. This technique gives access to important thermodynamic properties of matter at extreme conditions, such as temperature, material sound speed, and viscosity. The successful realization of this method requires the acquisition of many identical laser-pump/X-ray-probe shots, allowing the collection of a sufficient number of photons necessary to perform quantitative analyses. Here, a 2.5-fold improvement in the energy resolution of the instrument relative to previous works at the Matter in Extreme Conditions (MEC) endstation, Linac Coherent Light Source (LCLS), and the High Energy Density (HED) instrument, European XFEL, is presented. Some aspects of the experimental design that are essential for improving the number of photons detected in each X-ray shot, making such measurements feasible, are discussed. A careful choice of the energy resolution, the X-ray beam mode provided by the XFEL, and the position of the analysers used in such experiments can provide a more than ten-fold improvement in the photometrics. The discussion is supported by experimental data on 10 µm-thick iron and 50 nm-thick gold samples collected at the MEC endstation at the LCLS, and by complementary ray-tracing simulations coupled with thermal diffuse scattering calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Investigating PVC polymer–plasticizer interactions with atomistic MD simulations and potential of mean force calculations

For this work, atomistic molecular dynamics (MD) simulations coupled with potential of mean force (PMF) calculations were employed to investigate the interactions between PVC polymer chains containing 6–20 repeating units and various plasticizers, with the goal of identifying potential replacements for the toxic plasticizer di(2-ethylhexyl)phthalate (DEHP) used in blood bags. The selected plasticizers belong to various chemical families, such as orthophthalates, citrates, adipates, and the terephthalate DEHT. Both the polymer and the plasticizers lack ionizable groups, and their interactions are primarily governed by van der Waals and electrostatic forces. A model correlating PMF profiles with interaction forces was developed and validated across polyvinyl chloride (PVC) polymers of different lengths and all investigated plasticizers. This model provides insight into how structural variations in plasticizers influence their respective PMF values. The study ranks the investigated plasticizers based on binding affinity, identifying TOTM (tris(2-ethylhexyl) trimellitate, TEHTM), BTHC (butyryl trihexyl citrate, Citroflex B-6), ATHC (acetyl trihexyl citrate, Citroflex A-6, CA-6), and DEHT (bis(2-ethylhexyl)terephthalate, DOTP/DEHTP) as promising alternatives for further investigation. This comprehensive study encompasses PVC polymers of four different lengths and 14 plasticizers, with all data averaged over 30 independent simulations. The approach provides a deeper understanding of molecular interactions, enabling the tailoring of polymer–plasticizer systems for diverse applications, including extractables and leachables, with relevance spanning materials science to biomedical engineering, and also serves as a basis for developing coarse-grained simulation protocols. Overall, this study provides valuable insights for designing safer and more efficient plasticizer substitutes and is well supported by other studies.

Shet, Sai Athmeeya G. [Sri Sathya Sai Institute of↗
DOE OSTI2025

Acidities of MgO surface sites: implications for the formation mechanism of Mg(OH) 2

The hydroxylation of periclase (MgO) to brucite (Mg(OH) 2 ) is thought to be an important intermediate step when using MgO to capture CO 2 from the atmosphere. However, the mechanism of hydroxylation of MgO to form Mg(OH) 2 is poorly understood. In this work, we used atomic-scale density functional tight binding simulations coupled with the metadynamics rare event method to analyze the surface chemistry of MgO and the acid dissociation equilibrium constants (pK a ) of its surface sites. The method and parameters were validated by calculating the pK a for hydroxylation of the first shell water bound to aqueous Mg 2+ ion. The pK a value derived using a probabilistic method was 12.3, which is in fair agreement with the accepted value of 11.4, with the difference between them equal to a ∼5 kJ mol −1 error in the calculations. We then extended these pK a calculations to probe the hydroxylation reactions of the surface sites of the MgO(100)–water interface, arriving at pK a s of 5.4 to deprotonate terminal water molecules bound to the surface magnesium sites (η-OH 2 or 〉MgOH 2 ), and 13.9 to deprotonate hydroxylated bridging oxygen sites (μ 5 -oxo or 〉O). Hydroxide (OH − ) adsorption on the surface was also probed and found to be less thermodynamically favorable than deprotonation of the terminal water molecule. The plausibility of the computed pK a s was verified using an activity-based speciation model and compared to pH measurements of water equilibrated with MgO nanoparticles and single crystals. The model predicted a solution pH of 7.1 when surface sites buffered and the pH of 12.0 when MgO dissolution dominated. These are close to the experimental initial solution pHs of 7–7.5 and the long term pHs of ∼10.5. The similarity suggests that the calculated pK a values from the DFTB+/metadynamics simulations are plausible and that these methods can be a useful tool to probe reaction mechanisms involving covalent bonds.

Adapa, Sai Krishna Reddy [Oak Ridge National Labor↗
DOE OSTI2025

Sensitivities of Large Eddy Simulations of Aerosol Plume Transport and Cloud Response

Abstract Cloud responses to surface‐based sources of aerosol perturbation partially depend on how turbulent transport of the aerosol to cloud base affects the spatial and temporal distribution of aerosol. Here, scenarios of plume injection below a marine stratocumulus cloud are modeled using large eddy simulations coupled to a prognostic bulk aerosol and cloud microphysics scheme. Both passive plumes, consisting of an inert tracer, and active plumes are investigated, where the latter are representative of saltwater droplet plumes such as have been proposed for marine cloud brightening. Passive plume scenarios show higher in‐plume cloud brightness (relative to out‐of‐plume) due to the predominant transport of the passive plume tracer from the near‐surface to the cloud layer within updrafts. These updrafts rise into brighter areas within the cloud deck, even in the absence of an aerosol perturbation associated with an active plume. Comparing albedo at in‐plume to out‐of‐plume locations associates the inert plume with the brightest cloud locations, without any causal effect of the plume on the cloud. Numerical sensitivities are first assessed to establish a suitable model configuration. Then sensitivity to particle injection rate is investigated. Trade‐offs are identified between the number of injected particles and the suppressive effect of droplet evaporation on plume loft and spread. Furthermore, as the near‐field in‐plume brightening effect does not depend significantly on injection rate given a suitable definition of perturbed versus unperturbed regions of the flow, plume area is a key controlling factor on the overall cloud brightening effect of an aerosol perturbation.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Impact of Calibrated Land Surface Model Parameters on the Accuracy and Uncertainty of Land-Atmosphere Coupling in WRF Simulations

Land-atmosphere (L-A) interactions play a critical role in determining the diurnal evolution of both planetary boundary layer (PBL) and land surface temperature and moisture budgets, as well as controlling feedbacks with clouds and precipitation that lead to the persistence of dry and wet regimes. Recent efforts to quantify the strength of L-A coupling in prediction models have produced diagnostics that integrate across both the land and PBL components of the system. In this study, we examine the impact of improved specification of land surface states, anomalies, and fluxes on coupled WRF forecasts during the summers of extreme dry (2006) and wet (2007) land surface conditions in the U.S. Southern Great Plains. The improved land initialization and surface flux parameterizations are obtained through the use of a new optimization and uncertainty estimation module in NASA's Land Information System (LIS-OPT/UE), whereby parameter sets are calibrated in the Noah land surface model and classified according to a land cover and soil type mapping of the observation sites to the full model domain. The impact of calibrated parameters on the a) spinup of the land surface used as initial conditions, and b) heat and moisture states and fluxes of the coupled WRF simulations are then assessed in terms of ambient weather and land-atmosphere coupling along with measures of uncertainty propagation into the forecasts. In addition, the sensitivity of this approach to the period of calibration (dry, wet, average) is investigated. Finally, tradeoffs of computational tractability and scientific validity, and the potential for combining this approach with satellite remote sensing data are also discussed.

Santanello, Joseph A., Jr.↗
DOE OSTI2025

Multiscale Thermal-hydraulic analysis of the MARVEL microreactor using a coupled SAM and SCM simulation. (PPTX)

This is .pptx document presenting a summary of the paper of the same name that has been submitted and accepted to NURETH-21: MARVEL is a natural-convection-cooled sodium-potassium microreactor that is anticipated to generate 85 kilowatts of thermal energy. It will operate within Idaho National Laboratory Transient Reactor Test Facility and is being developed by the DOE Microreactor Program. MARVEL will be used to test microreactor applications, evaluate systems for remote monitoring, and develop autonomous control technologies. A thermal-hydraulic computational model of this facility is a valuable tool to study important transients and calculate the safety limits of the micro-reactor design. For this purpose, the authors have chosen to use a multiscale coupled simulation: SCM for modeling the reactor core and SAM for the reactor's primary cooling system. SCM is MOOSE physics module for subchannel analysis, which was designed to model single-phase flows through liquid-metal cooled, wire-wrapped fuel pin sub-assemblies, ordered in a triangular lattice. The SCM code was modified to be able to model MARVEL?s unique geometry. SAM is a systems analysis module based on the MOOSE framework. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for various advanced reactor types. The coupling between the two SCM and SAM for MARVEL modeling is done implementing a domain over-lapping approach. The resulting coupled simulation can model transients such as reactor startup/shutdown and provide an intermediate fidelity picture of the temperature field and other variables, in the core. Results for the steady-state simulations are presented in the article as well as flow blockage transient.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2023

Iridescence from Total Internal Reflection at 3D Microscale Interfaces: Mechanistic Insights and Spectral Analysis

An experimental investigation and the optical modeling of the structural coloration produced from total internal reflection interference within 3D microstructures are described. Ray-tracing simulations coupled with color visualization and spectral analysis techniques are used to model, examine, and rationalize the iridescence generated for a range of microgeometries, including hemicylinders and truncated hemispheres, under varying illumination conditions. An approach to deconstruct the observed iridescence and complex far-field spectral features into its elementary components and systematically link them to ray trajectories that emanate from the illuminated microstructures is demonstrated. The results are compared with experiments, wherein microstructures are fabricated with methods such as chemical etching, multiphoton lithography, and grayscale lithography. Microstructure arrays patterned on surfaces with varying orientation and size lead to unique color-traveling optical effects and highlight opportunities for how total internal reflection interference can be used to create customizable reflective iridescence. The findings herein provide a robust conceptual framework for rationalizing this multibounce interference mechanism and establish approaches for characterizing and tailoring the optical and iridescent properties of microstructured surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Downwind dispersion of CO 2 from a major subsea blowout in shallow offshore waters

Growing interest in offshore geologic carbon sequestration (GCS) motivates risk assessment of large-scale subsea CO 2 well blowouts or pipeline ruptures. For major leaks of CO 2 from wells or pipelines, significant fluxes of CO 2 may occur from the sea surface depending on water depth. In the context of risk assessment of human health and safety, we have used previously simulated coupled well-reservoir and water column model results as a source term for dense gas dispersion of CO 2 above the sea surface. The models are linked together by one-way coupling, that is, output of one model is used as input to the next model. These first-of-their-kind coupled flow results are applicable to assessing the hazard of CO 2 to people at and downwind of the sea surface location of emission. Hazard is quantified by plotting the downwind dispersion length (DDL), which we define in the study as the distances from the emission source to the point at which the emitted CO 2 has been diluted to 5% and 1.5% in air by volume. Here results suggest that large-scale blowouts in shallow water (10 m) may cause hazardous CO 2 plumes extending on the order of several hundred meters downwind. Details of the modeling show DDL has a maximum for windspeed (at an elevation of 10 m) of approximately 5 m/s, with smaller DDL for both weaker and stronger winds. This is explained by the fact that wind favors transport but also causes dispersion; therefore there is a certain wind speed that maximizes DDL.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Mobility assessment of the BCC and carbide phases in the C-Nb, C-U and Nb-U systems

Uranium carbides with refractory metal additions are considered for Gen IV nuclear reactors and nuclear thermal propulsion as fuels for their high-temperature and corrosion resistant properties. Understanding kinetic effects that dictate microstructural evolution during fabrication and operating conditions is essential to advance technological development of these fuels. This work presents the development of an atomic mobility database for C-Nb-U systems based off available experimental data supported with ab-initio methods. The mobility assessments and uncertainty quantification (using Markov chain Monte Carlo) were conducted in the Kawin software. Carbon diffusion is considered dominant, as metal diffusion is much slower, with niobium diffusion being even slower and rate limiting than uranium metal. We provide a comprehensive and self-consistent thermo-kinetic database that is validated by diffusion couple simulations through Kawin. In conclusion, this enables prediction of microstructural and phase evolution critical for the development and lifetime assessment of next generation nuclear fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

SA508 low alloy steel to 316L stainless steel dissimilar metal joint made by powder metallurgy hot isostatic pressing

Joining ferritic SA508 low alloy steel (LAS) and austenitic 316 L stainless steel (SS) via powder metallurgy hot isostatic pressing (PM-HIP) was evaluated as an alternative method to welding. This study investigated the mechanical and microstructural evolutions of the bimetallic interface under different joint designs and heat treatments. Here, the direct joining of dissimilar metal alloys by PM-HIP method resulted in two designs: 1) powder SA508 to solid bar 316 L (P508–B316L) and 2) powder 316 L solid bar SA508 (P316L-B508). In both cases, P508–B316L and P316L-demonstrated satisfactory tensile strength, however, high hardness and severe depreciation in toughness were located on the bimetallic interface. The mechanisms responsible for the detrimental mechanical properties were verified. Large oxides were observed in P508–B316L due to the prior powder boundary (PPB) oxides present in SA508 powder. The intense sensitization occurred from the formation of M 23 C 6 carbides, consequently from the slow cooling after PM-HIP in P316L-B508. Post-HIP heat treatments were explored to reduce the distance of carbide formation; however, the heat treatment could not eliminate the carbides. The experimental results were compared to the diffusion couple simulation as a function of carbide formation with distance. The analysis also showed the high hardness at the bimetallic interface was primarily contributed by solid solution strengthening. In conclusion, the direct joining of P316L-B508 and P508–B316L via PM-HIP was deemed to be unfeasible, and a transitional material is necessary to impede the diffusion of carbon.

316L stainless Steel↗
DOE OSTI2023

First measurement of prompt and non-prompt D* + vector meson spin alignment in pp collisions at $\sqrt{s}$ = 13 TeV

This letter reports the first measurement of spin alignment, with respect to the helicity axis, for D* + vector mesons and their charge conjugates from charm-quark hadronisation (prompt) and from beauty meson decays (non-prompt) in hadron collisions. The measurements were performed at midrapidity (|y| < 0.8) as a function of transverse momentum (p T ) in proton–proton (pp) collisions collected by ALICE at the centre-of-mass energy $\sqrt{s}$ = 13 TeV. The diagonal spin density matrix element $ρ$00 of D* + mesons was measured from the angular distribution of the D* + → D 0 (→ K - π + )π + decay products, in the D* + rest frame, with respect to the D* + momentum direction in the pp centre of mass frame. The $ρ$00 value for prompt D* + mesons is consistent with 1/3, which implies no spin alignment. However, for non prompt D* + mesons an evidence of $ρ$00 larger than 1/3 is found. The measured value of the spin density element is $ρ$00 = 0.455±0.022(stat.)±0.035(syst.) in the 5 < p T < 20 GeV/c interval, which is consistent with a Pythia 8 Monte Carlo simulation coupled with the EvtGen package, which implements the helicity conservation in the decay of D* + meson from beauty mesons. In non-central heavy-ion collisions, the spin of the D* + mesons may be globally aligned with the direction of the initial angular momentum and magnetic field. Based on the results for pp collisions reported in this letter it is shown that alignment of non-prompt D* + mesons due to the helicity conservation coupled to the collective anisotropic expansion may mimic the signal of global spin alignment in heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Molecular Mechanisms of Sorbed Ion Effects during Boehmite Particle Aggregation

Classical theories of particle aggregation, such as Derjaguin–Landau–Verwey–Overbeek (DLVO), do not explain recent observations of ion-specific effects or the complex concentration dependence for aggregation. Thus, here, we probe the molecular mechanisms by which selected alkali nitrate ions (Na + , K + , and NO 3 –) influence aggregation of the mineral boehmite (γ-AlOOH) nanoparticles. Nanoparticle aggregation was analyzed using classical molecular dynamics (CMD) simulations coupled with the metadynamics rare event approach for stoichiometric surface terminations of two boehmite crystal faces. Calculated free energy landscapes reveal how electrolyte ions alter aggregation on different crystal faces relative to pure water. Consistent with experimental observations, we find that adding an electrolyte significantly reduces the energy barrier for particle aggregation (~3–4×). However, in this work, we show this is due to the ions disrupting interstitial water networks, and that aggregation between stoichiometric (010) basal–basal surfaces is more favorable than between (001) edge–edge surfaces (~5–6×) due to the higher interfacial water densities on edge surfaces. The interfacial distances in the interlayer between aggregated particles with electrolytes (~5–10 Å) are larger than those in pure water (a few Ångströms). Together, aggregation/disaggregation in salt solutions is predicted to be more reversible due to these lower energy barriers, but there is uncertainty on the magnitudes of the energies that lead to aggregation at the molecular scale. By analyzing the peak water densities of the first monolayer of interstitial water as a proxy for solvent ordering, we find that the extent of solvent ordering likely determines the structures of aggregated states as well as the energy barriers to move between them. Further, the results suggest a path for developing a molecular-level basis to predict the synergies between ions and crystal faces that facilitate aggregation under given solution conditions. Such fundamental understanding could be applied extensively to the aggregation and precipitation utilization in the biological, pharmaceutical, materials design, environmental remediation, and geological regimes.

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

Integration of CeO 2 -Based Memristor with Vertically Aligned Nanocomposite Thin Film: Enabling Selective Conductive Filament Formation for High-Performance Electronic Synapses

The CeO 2 -based memristor has attracted significant attention due to its intrinsic resistive switching (RS) properties, large on/off ratio, and great plasticity, making it a promising candidate for artificial synapses. However, significant challenges such as high power consumption and poor device reliability hinder its broad application in neuromorphic microchips. To tackle these issues, in this work, we design a novel bilayer (BL) memristor by integrating a CeO 2 -based memristor with a Co-CeO 2 vertically aligned nanocomposite (VAN) layer and compare it with the single layer (SL) memristor. Preliminary electrical testing reveals that the BL memristor offers a reduced set/reset voltage (~67% lower), a higher on/off ratio (~5 × 10 2 ), enhanced device reliability, and improved device-to-device variation compared to the SL memristor. Insight from COMSOL simulation, coupled with microstructural analysis, provides a comprehensive elucidation on how the VAN layer facilitates the selective conductive filament (CF) formation. Subsequently, the plasticity of the BL memristor is evaluated through long-term potentiation/depression (LTP/LTD), paired-pulse facilitation (PPF), and spike-time-dependent plasticity (STDP). The spiking neural network (SNN) built upon the BL memristor achieves remarkable accuracy (~94%) after only 12 iterations, underscoring its potential for high-performance neural networks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗