Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Chain rule”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Development of a Framework and Methodology for an Advanced Reactor Materials Environmental Effects Design Guide

Advanced non-light-water reactor components may operate at elevated temperature while experiencing cyclic loading, significant neutron irradiation, and exposure to reactor coolant. ASME Boiler and Pressure Vessel Code, Section III, Division 5, provides design rules for elevated-temperature service but does not include specific procedures to account for environmental effects on material properties. This report develops an initial framework and methodology for an Environmental Effects Design Guide (EEDG) focused on neutron irradiation; coolant-environment effects are reserved for future work. The proposed approach treats irradiation as a property-based overlay on the existing Division 5 design process, with two routes: a sparse-data route applying two reduction factors — FCR on creep-rupture strength and FF on fatigue life — for the creep-fatigue evaluations that typically control the design of advanced high-temperature reactor components, and a fuller framework developing the property-to-rule chain across the four Division 5 checks (primary load, strain limits and ratcheting, creep-fatigue, and buckling), together with swelling and weldments as scope items. Both routes are scoped by an in-pile qualification that restricts the use of post-irradiation-examination-derived properties in regimes where an in-pile mechanism could control the design outcome. Illustrative outputs derived on a compiled annealed Type 316 database — FCR ≈ 0.78–0.86 and FF ≈ 0.4 — demonstrate the calculation method within that specific dataset. The framework is an initial, testable design-rule concept; it identifies a practical path for preliminary design evaluations under sparse data and the material data and testing needed to develop the framework further.

Barua, Bipul (ORCID:0000000247184113)↗
DOE OSTI2021

Deterministic symbolic regression with derivative information: General methodology and application to equations of state

Symbolic regression methods simultaneously determine the model functional form and the regression parameter values by generating expression trees. Symbolic regression can capture the complexity of real–world phenomena but the use of deterministic optimization for symbolic regression has been limited due to the complexity of the search space of existing formulations. Herein we present a novel deterministic mixed–integer nonlinear programming formulation for symbolic regression that incorporates derivative constraints through auxiliary expression trees. By applying the chain rule to mathematical operations, binary expression trees are capable of representing the calculation of first and second derivatives. We apply this formulation to illustrative examples using derivative information to show increased model discrimination capability. In addition, we perform a case study of a thermodynamic equation of state to gain insight on valid functional forms with thermodynamics–based constraints on the first and second derivatives.

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

Automatic Differentiation in MetaPhysicL and Its Applications in MOOSE

Efficient solution via Newton’s method of nonlinear systems of equations requires an accurate representation of the Jacobian, corresponding to the derivatives of the component residual equations with respect to the degrees of freedom. In practice these systems of equations often arise from spatial discretization of partial differential equations used to model physical phenomena. These equations may involve domain motion or material equations that are complex functions of the systems’ degrees of freedom. Computing the Jacobian by hand in these situations is arduous and prone to error. Finite difference approximations of the Jacobian or its action are prone to truncation error, especially in multiphysics settings. Symbolic differentiation packages may be used, but often result in an excessive number of terms in realistic model scenarios. An alternative to symbolic and numerical differentiation is automatic differentiation (AD), which propagates derivatives with every elementary operation of a computer program, corresponding to continual application of the chain rule. Automatic differentiation offers the guarantee of an exact Jacobian at a relatively small overhead cost. In this work, we outline the adoption of AD in the Multiphysics Object Oriented Simulation Environment (MOOSE) via the MetaPhysicL package. We describe the application of MOOSE’s AD capability to several sets of physics that were previously infeasible to model via hand-coded or Jacobian-free simulation techniques, including arbitrary Lagrangian-Eulerian and level-set simulations of laser melt pools, phase-field simulations with free energies provided through neural networks, and metallic nuclear fuel simulations that require inner Newton loop calculation of nonlinear material properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A Trajectory-Driven Algorithm for Differentiating SRB Measures on Unstable Manifolds

Sinai-Ruelle-Bowen (SRB) measures are limiting stationary distributions describing the statistical behavior of chaotic dynamical systems. Directional derivatives of SRB measure densities conditioned on unstable manifolds are critical in the sensitivity analysis of hyperbolic chaos. These derivatives, known as the SRB density gradients, are by-products of the regularization of Lebesgue integrals appearing in the original linear response expression. In this paper, we propose a novel trajectory- driven algorithm for computing the SRB density gradient defined for systems with high-dimensional unstable manifolds. We apply the concept of measure preservation together with the chain rule on smooth manifolds. Due to the recursive one-step nature of our derivations, the proposed procedure is memory-efficient and can be naturally integrated with existing Monte Carlo schemes widely used in computational chaotic dynamics. Here, we numerically show the exponential convergence of our scheme, analyze the computational cost, and present its use in the context of Monte Carlo integration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Efficient Active Learning for Gaussian Process Classification by Error Reduction

Active learning sequentially selects the best instance for labeling by optimizing an acquisition function to enhance data/label efficiency. The selection can be either from a discrete instance set (pool-based scenario) or a continuous instance space (query synthesis scenario). In this work, we study both active learning scenarios for Gaussian Process Classification (GPC). The existing active learning strategies that maximize the Estimated Error Reduction (EER) aim at reducing the classification error after training with the new acquired instance in a onestep-look-ahead manner. The computation of EER-based acquisition functions is typically prohibitive as it requires retraining the GPC with every new query. Moreover, as the EER is not smooth, it can not be combined with gradient-based optimization techniques to efficiently explore the continuous instance space for query synthesis. To overcome these critical limitations, we develop computationally efficient algorithms for EER-based active learning with GPC. Further, we derive the joint predictive distribution of label pairs as a one-dimensional integral, as a result of which the computation of the acquisition function avoids retraining the GPC for each query, remarkably reducing the computational overhead. We also derive the gradient chain rule to efficiently calculate the gradient of the acquisition function, which leads to the first query synthesis active learning algorithm implementing EER-based strategies. Our experiments clearly demonstrate the computational efficiency of the proposed algorithms. We also benchmark our algorithms on both synthetic and real-world datasets, which show superior performance in terms of sampling efficiency compared to the existing state-of-the-art algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Decay rates of almost strong modes in Floquet spin chains beyond Fermi's Golden Rule

The stability and dynamics of almost strong zero and π modes in weakly nonintegrable Floquet spin chains are investigated. Such modes can also be viewed as localized Majorana modes at the edge of a topological superconductor. Perturbation theory in the strength of integrability breaking interaction J z is employed to estimate the decay rates of these modes, and compared to decay rates obtained from exact diagonalization. Here, the structure of the perturbation theory and thus the lifetime of the modes is governed by the conservation of quasienergy modulo 2π/T, where T is the period of the Floquet system. If there are 4n – 1 bulk excitations whose quasienergies add up to zero (or π/T for a π mode), one obtains a decay channel of the Majorana mode with a decay rate proportional to $J$$^{2n}_{z}$. Thus the lifetime is sensitively controlled by the width of the single-particle Floquet bands. For regimes where the decay rates are quadratic in J z , an analytic expression for the decay rate in terms of an infinite temperature autocorrelation function of the integrable model is derived, and shown to agree well with exact diagonalization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Side chain engineering in indacenodithiophene- co -benzothiadiazole and its impact on mixed ionic–electronic transport properties

Organic semiconductors are increasingly being decorated with hydrophilic solubilising chains to create materials that can function as mixed ionic–electronic conductors, which are promising candidates for interfacing biological systems with organic electronics. While numerous organic semiconductors, including p- and n-type materials, small molecules and polymers, have been successfully tailored to encompass mixed conduction properties, common to all these systems is that they have been semicrystalline materials. Here, we explore how side chain engineering in the nano-crystalline indacenodithiophene-co-benzothiadiazole (IDTBT) polymer can be used to instil ionic transport properties and how this in turn influences the electronic transport properties. This allows us to ultimately assess the mixed ionic–electronic transport properties of these new IDTBT polymers using the organic electrochemical transistor as the testing platform. Using a complementary experimental and computational approach, we find that polar IDTBT derivatives can be infiltrated by water and solvated ions, they can be electrochemically doped efficiently in aqueous electrolyte with fast doping kinetics, and upon aqueous swelling there is no deterioration of the close interchain contacts that are vital for efficient charge transport in the IDTBT system. Despite these promising attributes, mixed ionic–electronic charge transport properties are surprisingly poor in all the polar IDTBT derivatives. Albeit a “negative” result, this finding clearly contradicts established side chain engineering rules for mixed ionic–electronic conductors, which motivated our continued investigation of this system. We eventually find this anomalous behaviour to be caused by increasing energetic disorder in the polymers with increasing polar side chain content. We have investigated computationally how the polar side chain motifs contribute to this detrimental energetic inhomogeneity and ultimately use the learnings to propose new molecular design criteria for side chains that can facilitate ion transport without impeding electronic transport.

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

Closing the loop between microstructure and charge transport in conjugated polymers by combining microscopy and simulation

A grand challenge in materials science is to identify the impact of molecular composition and structure across a range of length scales on macroscopic properties. We demonstrate a unified experimental–theoretical framework that coordinates experimental measurements of mesoscale structure with molecular-level physical modeling to bridge multiple scales of physical behavior. Here we apply this framework to understand charge transport in a semiconducting polymer. Spatially-resolved nanodiffraction in a transmission electron microscope is combined with a self-consistent framework of the polymer chain statistics to yield a detailed picture of the polymer microstructure ranging from the molecular to device relevant scale. Using these data as inputs for charge transport calculations, the combined multiscale approach highlights the underrepresented role of defects in existing transport models. Short-range transport is shown to be more chaotic than is often pictured, with the drift velocity accounting for a small portion of overall charge motion. Local transport is sensitive to the alignment and geometry of polymer chains. At longer length scales, large domains and gradual grain boundaries funnel charges preferentially to certain regions, creating inhomogeneous charge distributions. While alignment generally improves mobility, these funneling effects negatively impact mobility. The microstructure is modified in silico to explore possible design rules, showing chain stiffness and alignment to be beneficial while local homogeneity has no positive effect. This combined approach creates a flexible and extensible pipeline for analyzing multiscale functional properties and a general strategy for extending the accesible length scales of experimental and theoretical probes by harnessing their combined strengths.

42 ENGINEERING↗
DOE OSTI2022

Small-world complex network generation on a digital quantum processor

Abstract Quantum cellular automata (QCA) evolve qubits in a quantum circuit depending only on the states of their neighborhoods and model how rich physical complexity can emerge from a simple set of underlying dynamical rules. The inability of classical computers to simulate large quantum systems hinders the elucidation of quantum cellular automata, but quantum computers offer an ideal simulation platform. Here, we experimentally realize QCA on a digital quantum processor, simulating a one-dimensional Goldilocks rule on chains of up to 23 superconducting qubits. We calculate calibrated and error-mitigated population dynamics and complex network measures, which indicate the formation of small-world mutual information networks. These networks decohere at fixed circuit depth independent of system size, the largest of which corresponding to 1,056 two-qubit gates. Such computations may enable the employment of QCA in applications like the simulation of strongly-correlated matter or beyond-classical computational demonstrations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Enhancing ACPF Analysis: Integrating Newton-Raphson Method with Gradient Descent and Computational Graphs

This paper presents a new method for enhancing Alternating Current Power Flow (ACPF) analysis. The method integrates the Newton-Raphson (NR) method with Enhanced-Gradient Descent (GD) and computational graphs. The integration of renewable energy sources in power systems introduces variability and unpredictability, and this method addresses these challenges. It leverages the robustness of NR for accurate approximations and the flexibility of GD for handling variable conditions, all without requiring Jacobian matrix inversion. Furthermore, computational graphs provide a structured and visual framework that simplifies and systematizes the application of these methods. The goal of this fusion is to overcome the limitations of traditional ACPF methods and improve the resilience, adaptability, and efficiency of modern power grid analyses. We validate the effectiveness of our advanced algorithm through comprehensive testing on established IEEE benchmark systems. Furthermore, our findings demonstrate that our approach not only speeds up the convergence process but also ensures consistent performance across diverse system states, representing a significant advancement in power flow computation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Kensor: Coordinated Intelligence from Co-Located Sensors

Internet of Things (IoT) is becoming more pervasive in many installations, including homes, manufacturing plants, and industrial facilities of all kinds. The data that IoT produces is a reflection of usual behavior such as daily routines and scheduled tasks, but also from unexpected behavior due to unintentional or undesirable abnormalities. Here, we focus on achieving coordinated intelligence about normal and abnormal phenomena from multiple sensors that are geographically colocated in close proximity, monitoring and controlling a set of co-located devices. Given a set of co-located sensors, we seek an intelligent approach that would automatically determine the “normal” patterns of behaviors among the correlated sensors. After normal behavior is extracted, later monitoring should detect any deviant variations over time. An example application is an entry monitoring and alert system for facilities such as nuclear reactors, where badge readers, door locks, lights, weight trackers and other co-located sensors at the entry point are collectively tracked. To address this problem, we identify the possible solution approach that can be used to solve its different variants. The implemented model is developed as a combination of rules and Markov Chain methods.

Kotevska, Olivera↗
DOE OSTI2023

Linear Viscoelastic Properties of Putative Cyclic Polymers Synthesized by Reversible Radical Recombination Polymerization (R3P)

Linear viscoelastic properties in both melt and solution states are reported for a series of poly(3,6-dioxa-1,8-octanedithiol) (polyDODT) made by reversible radical recombination polymerization (R 3 P) under conditions designed to produce linear (LDODT), cyclic (RDODT) and linear-cyclic mixtures (LRDODT). PolyDODT is amorphous (T g < -50 °C) and highly flexible (entanglement molecular weight M e,lin ≈ 1850 g/mol for LDODT). PolyDODT’s low T g and low M e,lin enable characterization over a wide dynamic range and a wide range of dimensionless weight-average molecular weight Z w = M w /M e,lin . Measurements at temperatures from -57 °C to 100 °C provide up to 18 deades of reduced frequency, which is necessary to characterize RDODT melts with Z w from 23 to 300. The two highest molecular weight polymers in the present RDODT series have such high M w (406k and 556k g/mol) that mass spectrometry, NMR spectroscopy, and even chemical assays for chain ends are unable to rule out up to 2 mol% linear contaminant. By studying the samples in solution (using dilution to reduce Z w ) we could compare their dynamics with those of previously established high-purity polystyrene (PS) rings (limited to Z w ≤ 13.6). RDODT solutions with Z w < 15 (concentrations <5wt% for RDODT-406k and 556k) have dynamic moduli G^* that accord with LCCC-purified PS rings in terms of the frequency dependence (including the absence of a plateau), the progression of shapes of G* as a function of Z w , and the linear scaling of their zero-shear viscosity η0 with M w . The shape of G* as a function of Z w for solutions of RDODT-406k and -556k also accord with lower M w RDODT melts (which have ≤ 1.3mol% linear contaminant). Thus, measurement of the linear viscoelastic properties of appropriate concentrations of high M w (>200k g/mol) putative cyclic polymers, in which linear chains evade spectroscopic detection, may provide an alternative means (though not fully proven) of validation of sample purity. When Z w > 15 (including all seven RDODT melts and eight of their solutions), G* has a rubbery plateau. This suggests that the onset of entanglement-like behavior in rings requires 4-5-fold greater Z w than is required for linear chains. Further, the plateau moduli of RDODT samples are indistinguishable from G N o of the corresponding LDODT (melt or matched-concentration solutions). In entangled linear polymers, the observation that G N o is independent of Z w follows from limitations on lateral fluctuations due to neighboring chains becoming independent of position along a given chain. The present results for RDODT suggest that this holds for sufficiently long endless chains, too. While the RDODT have the same G N o as entangled LDODT, when Z w > 60 the terminal relaxation, if reached at all, of RDODT extends to orders of magnitude lower frequency than an entangled linear polymer of the same Z w . Consequently, the viscosity of RDODT with Z w > 60 increases with Z w much more strongly than the 3.4-power observed for entangled linear polymers. Lastly, these novel polymers, with disulfide-linked backbone and broad relaxation time distribution may prove important in relation to biodegradable elastomers and materials with exceptional low-frequency dissipation, extending at least 12 decades below the onset of the rubbery plateau.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Topological electronic properties of silicon

The central role that materials play in human history is exemplified by the three-age division of prehistory into the stone, bronze, and iron ages. References to our present time as the information age or silicon age epitomizes the important role that this semiconducting material came to play in the development of computers and devices that permeate our daily lives. Here we show that the electronic states in silicon have nontrivial topological structures that are captured by a network of Berry flux lines that link at points of high symmetry in the Brillouin zone. Further, This complex network has ice-nodal points where fluxes satisfy ice rules, making silicon a nodal-chain insulator. Fixing the longitudinal momentum parallel to such flux lines yields a two-dimensional Dirac Hamiltonian for the transverse degrees of freedom. Similar to nodal-chain semimetals, we find drumheadlike states in the regions that are delimited by the projections of the bulk Berry flux network on the surface Brillouin zone.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Tuning the Mechanical and Electric Properties of Conjugated Polymer Semiconductors: Side–Chain Design Based on Asymmetric Benzodithiophene Building Blocks

Understanding molecular design rules for stretchable polymer semiconductors is important for enabling next generation stretchable electronic circuits. To simultaneously improve both electrical properties and mechanical stretchability, a design strategy is reported in introducing conjugated rigid fused-rings with bulky side groups in semiconducting polymers. In this work, the understanding of this design concept is improved by systematically investigating the effect of different types of bulky side groups asymmetrically substituted on conjugated polymer semiconductor backbones. Specifically, four types of side groups are investigated, including naphthalene (NaPh), biphenyl (PhPh), thienylphenyl (ThPh), and alkylphenyl (C4Ph), asymmetrically substituted on benzodithiophene units, namely asy-BDT. With the four types of side groups installed on BDT-containing conjugated polymers in an asymmetrical fashion, it is observed that they reduced the polymer chain aggregation and film crystallinity, hence improving the film stretchability. Furthermore, the fully conjugated polymer back-bone allows maintenance of good charge carrier mobilities. Specifically, polymer PDPP-C4Ph (with C4Ph side groups) shows the highest mobility in the fully stretchable transistor and maintained its mobility even after being subjected to hundreds of stretching-releasing cycles at 25% strain. Altogether, the results provide anunderstanding of the use of asymmetrically substituted fused-ring conjugated polymer structures to tune mechanical and charge transport properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-Pressure Rate Rules for Ether Alkylperoxy Radical Isomerization

The first isomerization reaction of an alkylperoxy (RO 2 ) radical holds significant importance in low-temperature oxidation, as it governs the branching ratios of the hydroperoxyalkyl (QOOH) radicals, which influence the competition between the chain-propagation and chain-branching reactions. In this study, we systematically calculated high-pressure rate rules for the RO 2 isomerization reaction of monoethers, exploring 5-, 6-, 7-, and 8-membered ring transition states. Primary, secondary, and tertiary carbon sites, where both the abstracting peroxy group and the abstracted hydrogen are located, were considered, with particular emphasis on distinguishing between secondary carbons adjacent (alpha) and nonadjacent to the ether functional group. Using the G4//B3LYP/6-311++G(2df,2pd) level of theory and the transition state theory, we estimated the rate constants and the Arrhenius coefficient for over 120 possible isomerization reactions. We examined the effect of ring size and ring atoms, revealing that 6- and 7-membered ring isomerizations were generally the fastest. The impact of the ether functional group on transition states was investigated by comparing reactions with identical ring size, peroxy, and radical positions, but with the ether functional group positioned either outside (i.e., out) or inside (i.e., in) the transition state ring, leading to differences in the rate constants. When comparing to analogous alkane rate constants, differences of up to an order of magnitude were observed, underscoring the need for caution when assigning rate rules by analogy. We applied our rate constants in the di-iso-butyl ether kinetic model and evaluated their influence on low-temperature chemistry finding that they altered the branching ratios by up to a factor of 9, highlighting the significance of site-specific rate constants for more accurate low-temperature modeling.

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

Fine‐Tuning Alkyl Chains on Quinoxaline Nonfullerene Acceptors Enables High‐Efficiency Ternary Organic Solar Cells with Optimizing Molecular Stacking and Reducing Energy Loss

Abstract Material design of guest acceptor is always a big challenge for improving the efficiency of ternary organic solar cells (OSCs). Here, a pair of isomeric nonfullerene acceptors based on quinoxaline core, Qx–p‐C 7 H 8 O and Qx–m‐C 7 H 8 O, is designed and synthesized. By moving the alkoxy chain attached on side phenyl frommeta‐position topara‐position, both π–π stacking distance and crystallinity are enhanced simultaneously. They obtain the uplifted lowest unoccupied molecular orbital level. Compared to Qx–m‐C 7 H 8 O, Qx–p‐C 7 H 8 O exhibits wider absorption spectrum and higher extinction coefficient. Using D18‐Cl:N3 as host materials, the addition of guest acceptor Qx–p‐C 7 H 8 O significantly improves the power conversion efficiency (PCE) from 17.61% to 18.49% because of higher open‐circuit voltage (0.875 V) and short‐circuit current density (27.85 mA cm −2 ). This can be attributed to the faster exciton dissociation, more balanced carrier mobility, fine fiber morphology, and lower energy loss in the ternary devices. However, Qx–m‐C 7 H 8 O‐based ternary device achieves relatively low PCE of 17.17% because this device shows extremely low electron mobility. The results indicate that molecular stacking, film morphology, etc., can be effectively modulated by fine‐tuning the side chains of guest materials, which may be an effective design rule for further improving the PCE of OSCs.

Chemistry↗
DOE OSTI2020

Effect of Molecular Weight and Layer Thickness on Dielectric Breakdown Strength of Neat and Homopolymer Swollen Lamellar Block Copolymer Films

Designing next-generation lightweight pulsed power devices hinges on understanding the factors influencing the energy storage performance of dielectric materials. Polymer dielectric films have a quadratic dependence of energy storage on the voltage breakdown strength and strategies to enhance the breakdown strength are expected to yield a path toward high energy storage densities. Highly stratified lamellar block copolymer (L-BCP) films of model polystyrene-b-polymethylmethacrylate (PS-b-PMMA) exhibited as much as ~50% enhancement in breakdown voltage (E BD ), (225 % increase in stored energy density, U~ E BD 2 ), compared to unordered as-cast L-BCP films. Such an energy density using amorphous polymer is on par with industry-standard semi-crystalline biaxially oriented polypropylene (BOPP), and as such a notable development in the field. This work develops a deeper understanding of the molecular mechanisms of E BD enhancement in L-BCP films, due to combination of interface formation and spatial isolation of the chain ends into segregated zones. We further confirm the role of chain ends in the breakdown process blending a low M n L-BCP with matched M n homopolymers to attain same layer spacing as neat L-BCP of higher M n . E BD shows a significant decrease at low homopolymer fractions due to increased net chain end density within swollen ordered L-BCP domains in wet-brush regime, followed by increased E BD due to layer thickness increase via segregated “interphase layer” formation by excess homopolymers. Notably, E BD of homopolymer swollen L-BCPs is always lower than neat L-BCP of the same domain spacing, due to overall adverse chain ends contribution from homopolymers. These findings provide important selection rules for L-BCPs for designing next-generation flexible electronics with high energy density solid-state BCP film capacitors.

36 MATERIALS SCIENCE↗