Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential transformation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Hutchinson Trace Estimation for high-dimensional and high-order Physics-Informed Neural Networks

Physics-Informed Neural Networks (PINNs) have proven effective in solving partial differential equations (PDEs), especially when some data are available by seamlessly blending data and physics. However, extending PINNs to high-dimensional and even high-order PDEs encounters significant challenges due to the computational cost associated with automatic differentiation in the residual loss function calculation. Herein, we address the limitations of PINNs in handling high-dimensional and high-order PDEs by introducing the Hutchinson Trace Estimation (HTE) method. Starting with the second-order high-dimensional PDEs, which are ubiquitous in scientific computing, HTE is applied to transform the calculation of the entire Hessian matrix into a Hessian vector product (HVP). This approach not only alleviates the computational bottleneck via Taylor-mode automatic differentiation but also significantly reduces memory consumption from the Hessian matrix to an HVP’s scalar output. We further showcase HTE’s convergence to the original PINN loss and its unbiased behavior under specific conditions. Comparisons with the Stochastic Dimension Gradient Descent (SDGD) highlight the distinct advantages of HTE, particularly in scenarios with significant variability and variance among dimensions. We further extend the application of HTE to higher-order and higher-dimensional PDEs, specifically addressing the biharmonic equation. By employing tensor-vector products (TVP), HTE efficiently computes the colossal tensor associated with the fourth-order high-dimensional biharmonic equation, saving memory and enabling rapid computation. The effectiveness of HTE is illustrated through experimental setups, demonstrating comparable convergence rates with SDGD under memory and speed constraints. Additionally, HTE proves valuable in accelerating the Gradient-Enhanced PINN (gPINN) version as well as the Biharmonic equation. Overall, HTE opens up a new capability in scientific machine learning for tackling high-order and high-dimensional PDEs.

Curse of dimensionality↗
DOE OSTI2024

A Comprehensive Review of Latent Space Dynamics Identification Algorithms for Intrusive and Non-Intrusive Reduced-Order-Modeling

Numerical solvers of partial differential equations (PDEs) have been widely employed for simulating physical systems. However, the computational cost remains a major bottleneck in various scientific and engineering applications, which has motivated the development of reduced-order models (ROMs). Recently, machine-learning-based ROMs have gained significant popularity and are promising for addressing some limitations of traditional ROM methods, especially for advection dominated systems. In this chapter, we focus on a particular framework known as Latent Space Dynamics Identification (LaSDI), which transforms the high-fidelity data, governed by a PDE, to simpler and low-dimensional latent-space data, governed by ordinary differential equations (ODEs). These ODEs can be learned and subsequently interpolated to make ROM predictions. Each building block of LaSDI can be easily modulated depending on the application, which makes the LaSDI framework highly flexible. In particular, we present strategies to enforce the laws of thermodynamics into LaSDI models (tLaSDI), enhance robustness in the presence of noise through the weak form (WLaSDI), select high-fidelity training data efficiently through active learning (gLaSDI, GPLaSDI), and quantify the ROM prediction uncertainty through Gaussian processes (GPLaSDI). We demonstrate the performance of different LaSDI approaches on Burgers equation, a non-linear heat conduction problem, and a plasma physics problem, showing that LaSDI algorithms can achieve relative errors of less than a few percent and up to thousands of times speed-ups.

Computational Engineering, Finance, and Science (c↗
DOE OSTI2025

Phylogenomic discovery and engineering of nitrogen fixation into the bioenergy woody crop poplar

Biological nitrogen fixation (BNF) is a key process enabling plants in specific lineages to convert atmospheric dinitrogen (N₂) into bioavailable ammonia through symbioses with diazotrophic microbes. Expanding this capability beyond native nitrogen-fixing clades into non-nodulating crops would reduce synthetic fertilizer use, lowering energy inputs and environmental impacts in agriculture. Supported by DOE Funding Award DE-SC0018247, the NitFix project advanced foundational knowledge required to engineer root-nodule symbioses in new host species. The team generated the most comprehensive phylogenomic analysis to date of all known nodulating lineages, resolving the evolutionary history of nitrogen-fixing symbiosis and identifying core gene suites retained across nodulating taxa. Through multimodal genomics, transcriptomics, and functional analyses in Medicago truncatula and related species, the project mapped regulatory networks underlying nodule organogenesis, bacterial infection, and nitrogen-fixation efficiency. Key discoveries include the identification of conserved signaling modules for rhizobial recognition, transcription factors controlling nodule differentiation, and metabolic pathways integrating fixed nitrogen into plant growth. The project also developed enabling tools—including optimized transformation pipelines, gene-editing workflows, and imaging-based phenotyping—to accelerate engineering efforts in emerging models. Together, these results refine the mechanistic framework of symbiotic nitrogen fixation and highlight transferable components essential for rewiring these traits into non-nodulating crops.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A molecular mechanism for the generation of ligand-dependent differential outputs by the epidermal growth factor receptor

The epidermal growth factor receptor (EGFR) is a receptor tyrosine kinase that couples the binding of extracellular ligands, such as EGF and transforming growth factor-α (TGF-α), to the initiation of intracellular signaling pathways. EGFR binds to EGF and TGF-α with similar affinity, but generates different signals from these ligands. To address the mechanistic basis of this phenomenon, we have carried out cryo-EM analyses of human EGFR bound to EGF and TGF-α. We show that the extracellular module adopts an ensemble of dimeric conformations when bound to either EGF or TGF-α. The two extreme states of this ensemble represent distinct ligand-bound quaternary structures in which the membrane-proximal tips of the extracellular module are either juxtaposed or separated. EGF and TGF-α differ in their ability to maintain the conformation with the membrane-proximal tips of the extracellular module separated, and this conformation is stabilized preferentially by an oncogenic EGFR mutation. Close proximity of the transmembrane helices at the junction with the extracellular module has been associated previously with increased EGFR activity. Our results show how EGFR can couple the binding of different ligands to differential modulation of this proximity, thereby suggesting a molecular mechanism for the generation of ligand-sensitive differential outputs in this receptor family.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Determining the hydration energetics on carbon-supported Ru catalysts: An adsorption calorimetry and density functional theory study

Fundamental knowledge on the energetics at the interface between a water layer and a metal catalyst is essential so as to understand the roles that water can play in the synthesis, activation and regeneration of noble metal-based catalysts. Here, we report the direct measurement of the enthalpy of water adsorption (Δh ads ) on activated carbon (C) and activated C-supported Ru nanoparticles, which are promising catalyst as applied to the hydrogenation/hydrodeoxygenation (HDO) of oxygenates (phenolics, aldehydes, etc.). Specifically, the near-zero coverage enthalpy of water adsorption on a C-supported Ru catalyst is -75.3 ± 0.4 kJ/(mol water), suggesting favorable water–metal binding. This is much more exothermic than that on C, which has an enthalpy of adsorption of -50.3 ± 1.3 kJ/(mol water). Despite the favorable initial binding, the magnitudes of enthalpies of water condensation on C and Ru-C indicate that overall, their surfaces are both hydrophobic. Moreover, the experimentally-measured near-zero coverage water adsorption enthalpy at the Ru sites is in very good agreement with our density functional theory based calculations. At low coverages, we obtain a water binding energy of -61.7 kJ/(mol water), which increases to -78.1 kJ/(mol water) at saturation. Complementary results are also obtained from a thermal analysis, which employed a thermogravimetric analysis–differential scanning calorimetry–mass spectrometry (TG-DSC-MS), a spectroscopic investigation using ex situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) and a morphological evaluation with transmission electron microscopy (TEM). We point out that in carbon-supported metal catalysts, such as Ru-C, a strong hydration at near-zero coverage and relative weak water-surface interactions occurs upon saturation. Such heterogeneity is essential and crucial for catalytic hydrogenation/HDO reactions that involve balanced interactions among the water-rich reactant mixture and nonpolar organic products.

09 BIOMASS FUELS↗
DOE OSTI2023

Revisiting two thiophosphate compounds constituting d 0 transition metal HfP 2 S 6 and d 10 transition metal α-Ag 4 P 2 S 6 as multifunctional materials for combining second harmonic generation response and photocurrent response

Two acentric thiophosphate compounds, HfP 2 S 6 and α-Ag 4 P 2 S 6 , are revisited and studied as infrared nonlinear optical materials. HfP 2 S 6 and α-Ag 4 P 2 S 6 were structurally characterized without any property measurements. Here, in this work, HfP 2 S 6 and α-Ag 4 P 2 S 6 were synthesized via high temperature salt flux reactions. Low-temperature polymorph acentric α-Ag 4 P 2 S 6 was purified and grown as mm-sized crystals with the aid of AgBr flux. The AgBr flux was revealed to play an important role in stabilizing the acentric α-Ag 4 P 2 S 6 . The acentric α-Ag 4 P 2 S 6 transforms to centrosymmetric β-Ag 4 P 2 S 6 at 850(5) K, which is revealed by differential scanning calorimetry (DSC) analysis and powder X-ray diffraction experiments. HfP 2 S 6 and α-Ag 4 P 2 S 6 are discovered by UV-vis spectrum measurements as indirect bandgap semiconductors with bandgaps of 2.2(1) eV and 2.5(1) eV, respectively, which is supported by DFT calculations and TB-LMTO-ASA calculations. The bonding pictures of α-Ag 4 P 2 S 6 were studied by crystal orbital Hamilton population calculations (COHP) coupled with electron localization function (ELF) analysis. DFT calculations predict that HfP 2 S 6 and α-Ag 4 P 2 S 6 would exhibit different optical performances regardless of being constructed from identical [P 2 S 6 ] motifs. HfP 2 S 6 exhibits a low second harmonic generation (SHG) response, ~0.21 × AGS (for the sample of particle size of 25 μm). α-Ag 4 P 2 S 6 possesses moderate SHG response, ~0.61 × AGS (for the sample of particle size of 225 μm) coupled with a high laser damage threshold (LDT) of ~3.2 × AGS. Characteristics of high ambient stability, moderate bandgap and SHG response, type-I phase-matching capability, and high LDT together with the easy growth of large crystals make α-Ag 4 P 2 S 6 attractive for future infrared nonlinear optical applications. Photocurrent measurements found that α-Ag 4 P 2 S 6 and β-Ag 4 P 2 S 6 have high photocurrent response, 165 nA cm -2 and 135 nA cm -2 , respectively. α-Ag 4 P 2 S 6 is a new multifunctional material of the ternary Ag–P–S system, which combines nonlinear optical (NLO) properties and photocurrent response.

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

Gauges, loops, and polynomials for partition functions of graphical models

Graphical models represent multivariate and generally not normalized probability distributions. Computing the normalization factor, called the partition function, is the main inference challenge relevant to multiple statistical and optimization applications. The problem is #P-hard that is of an exponential complexity with respect to the number of variables. Here, aimed at approximating the partition function, we consider multi-graph models where binary variables and multivariable factors are associated with edges and nodes, respectively, of an undirected multi-graph. We suggest a new methodology for analysis and computations that combines the Gauge function technique from Chertkov and Chernyak with the technique developed in Anari and Oveis Gharan 2017 arXiv:1702.02937; Gurvits 2011 arXiv:1106.2844; Straszak and Vishnoi 2017 55th Annual Allerton Conf. on Communication, Control, and Computing, based on the recent progress in the field of real stable polynomials. We show that the Gauge function, representing a single-out term in a finite sum expression for the partition function which achieves extremum at the so-called belief-propagation gauge, has a natural polynomial representation in terms of gauges/variables associated with edges of the multi-graph. Moreover, Gauge function can be used to recover the partition function through a sequence of transformations allowing appealing algebraic and graphical interpretations. Algebraically, one step in the sequence consists of the application of a differential operator over gauges associated with an edge. Graphically, the sequence is interpreted as a repetitive elimination/contraction of edges resulting in multi-graph models on decreasing in size (number of edges) graphs with the same partition function as in the original multi-graph model. Even though the complexity of computing factors in the sequence of the derived multi-graph models and respective Gauge functions grow exponentially with the number of eliminated edges, polynomials associated with the new factors remain bi-stable if the original factors have this property. Moreover, we show that BP estimations in the sequence do not decrease, each low-bounding the partition function.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Linear Solvers for Collector Systems of Generalized Large-scale Inverter-Based Resources

Collector systems for inverter-based resources (IBRs) are typically represented by equivalent circuits for electromagnetic transient (EMT) simulations. Recent studies have revealed that modeling a detailed collector system is essential to accurately represent the behavior of IBRs, especially when dealing with partial tripping during external disturbances. However, there are several challenges in simulating a detailed EMT model of a collector system due to the time required to simulate such systems. Thus, this paper investigates the modeling of a detailed collector system, taking into account its configuration and components as defined in IEEE standard 2800. The configurations include the collector systems of generalized large-scale IBR plants. The components include the main IBR transformer, collector bus, and feeders with lines and/or cables. The EMT model of the collector system is represented by differential algebraic equations (DAEs) that are discretized to form linear equations that are solved using linear solvers. In this paper, linear solvers are proposed based on the Schur complement method, which are utilized for simulation of the EMT model of collector systems of generalized large-scale IBRs to accelerate simulation speed while maintaining the accuracy of the results. The proposed solvers are verified by comparing the performance to that of linear solvers provided in MATLAB.

Choi, Jongchan↗
DOE OSTI2024

Swift

Swift is a fast Fourier transform based spectral solver based on the MOOSE framework. It supports GPU accelerated semi-implicit solves of partial differential equations, such as those used for phase field mesoscale microstructure evolution.

Schwen, Daniel [Idaho National Laboratory (INL), I↗
DOE OSTI2023

Transformative Technology for FLASH Radiation Therapy

The general concept of radiation therapy used in conventional cancer treatment is to increase the therapeutic index by creating a physical dose differential between tumors and normal tissues through precision dose targeting, image guidance, and radiation beams that deliver a radiation dose with high conformality, e.g., protons and ions. However, the treatment and cure are still limited by normal tissue radiation toxicity, with the corresponding side effects. A fundamentally different paradigm for increasing the therapeutic index of radiation therapy has emerged recently, supported by preclinical research, and based on the FLASH radiation effect. FLASH radiation therapy (FLASH-RT) is an ultra-high-dose-rate delivery of a therapeutic radiation dose within a fraction of a second. Experimental studies have shown that normal tissues seem to be universally spared at these high dose rates, whereas tumors are not. While dose delivery conditions to achieve a FLASH effect are not yet fully characterized, it is currently estimated that doses delivered in less than 200 ms produce normal-tissue-sparing effects, yet effectively kill tumor cells. Despite a great opportunity, there are many technical challenges for the accelerator community to create the required dose rates with novel compact accelerators to ensure the safe delivery of FLASH radiation beams.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Spectroscopic Investigation of Catalyst Inks and Thin Films Toward the Development of Ionomer Quality Control

As the production of polymer electrolyte fuel cells expands, novel quality control methods must be invented or adapted in order to support expected rates of production. Ensuring the quality of deposited catalyst layers is an essential step in the fuel cell manufacturing process, as the efficiency of a fuel cell is reliant on the catalyst layer being uniform at both the target platinum loading and the target ionomer content. Implementing a quality control method that is sensitive to these aspects is imperative, as wasting precious metals and other catalyst materials is expensive, and represents a potential barrier to entry into the field for manufacturers experimenting with novel deposition processes. In this work, we analyzed catalyst inks to determine if their ionomer content could be quantized spectroscopically. Attenuated total reflection (ATR) Fourier transform infrared spectroscopic technique was investigated producing a signal proportional to the ionomer content. ATR spectroscopy was able to quantitatively differentiate samples in which the ionomer to carbon mass ratio (I/C) varied between 0.9 and 3.0. The I/C ratio was correlated to the measured ATR signal near the CF 2 vibrational bands located between 1100 cm −1 and 1400 cm −1 . The experimental results obtained constitute a step toward the development of novel quality control methodologies for catalyst inks utilized by the fuel cell industry.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Dynamic State Estimation for Radial Microgrid Protection

Dynamic state estimation (DSE) is a generalization of differential protection that offers a reduced likelihood of misoperation, particularly in the case of devices with nonlinear characteristics such as transformers which are being energized. It is also useful in cases where distance protection performs poorly, such as transmission lines with series compensation or mutually coupled transmission lines. DSE has been previously applied to microgrid branch protection. This study investigates the use of DSE for protection of radial portions of a microgrid or distribution system. This can be a challenge in microgrids or distribution systems with distributed generation on account of lack of fault current from inverter-interfaced generation, varying fault current between grid- connected and islanded modes, the potential for normally-meshed operation and unbalanced operation due to single-phase loads. Admittance relaying has been investigated as a solution for protection of microgrids, though it has been observed to have issues with grounded-wye connected loads, and additional relaying is necessary to prevent misoperation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Temporally quasiperiodic data, propagating in the laboratory frame, can be rendered periodic by Galilean transformation

For a broad class of distributions of temperature, concentration, or another quantity propagating rectilinearly, we show that temporally quasiperiodic behavior in the laboratory frame can be rendered periodic by Galilean transformation. The approach is illustrated analytically and numerically using as an example a closed-form model distribution generated from a one-dimensional partial differential equation, and a detailed process is developed to determine frame speed from more general quasiperiodic, one-dimensional, temporally- and spatially-discretized data. Furthermore, the approach is extended to two- and three-dimensional rectilinear propagation, and its application to nonrectilinear propagation, along with implications for interpreting noise-corrupted data, are also discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Spectral scheme for atomic structure calculations in density functional theory

In this study, we present a spectral scheme for atomic structure calculations in pseudopotential Kohn-Sham density functional theory. In particular, after applying an exponential transformation of the radial coordinates, we employ global polynomial interpolation on a Chebyshev grid, with derivative operators approximated using the Chebyshev differentiation matrix, and integrations using Clenshaw-Curtis quadrature. We demonstrate the accuracy and efficiency of the scheme through spin-polarized and unpolarized calculations for representative atoms, while considering local, semilocal, and hybrid exchange-correlation functionals. In particular, we find that $\mathcal{O}$(200) grid points are sufficient to achieve an accuracy of 1 microhartree in the eigenvalues for optimized norm conserving Vanderbilt pseudopotentials spanning the periodic table from atomic number Ζ = 1 to 83.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Latent Twins

Over the past decade, scientific machine learning has transformed the development of mathematical and computational frameworks for analyzing, modeling, and predicting complex systems. From inverse problems to numerical partial differential equations (PDEs), dynamical systems, and model reduction, these advances have pushed the boundaries of what can be simulated. Yet they have often progressed in parallel, with representation learning and algorithmic solution methods evolving largely as separate pipelines. With Latent Twins, we propose a unifying mathematical framework that creates a hidden surrogate in latent space for the underlying equations. Whereas digital twins mirror physical systems in the digital world, Latent Twins mirror mathematical systems in a learned latent space governed by operators. Through this lens, classical modeling, inversion, model reduction, and operator approximation all emerge as special cases of a single principle. We establish the fundamental approximation properties of Latent Twins for both ordinary differential equations (ODEs) and PDEs and demonstrate the framework across three representative settings: (i) canonical ODEs, capturing diverse dynamical regimes; (ii) a PDE benchmark using the shallow-water equations, contrasting Latent Twin simulations with deep operator network and forecasts with a four-dimensional variational method baseline; and (iii) a challenging real-data geopotential reanalysis dataset, reconstructing and forecasting from sparse, noisy observations. Latent Twins provide a compact, interpretable surrogate for solution operators that evaluate across arbitrary time gaps in a single-shot, while remaining compatible with scientific pipelines such as assimilation, control, and uncertainty quantification. Looking forward, this framework offers scalable, theory-grounded surrogates that bridge data-driven representation learning and classical scientific modeling across disciplines.

Latent Twins↗
DOE OSTI2026

3D reconstruction and neural rendering for adversarial machine learning

While evasion attacks on computer vision systems have been widely studied, creating attacks that remain effective under significant changes in viewpoint continues to be challenging. Traditional approaches often rely on affine transformations of images, but these approaches degrade at larger perspective shifts and often produce unrealistic or ineffective perturbations. Recent methods use differentiable renderers to improve viewpoint robustness, but they typically depend on manually constructed 3D models. We introduce a semi-automated pipeline that generates physically printable and perspective-invariant adversarial patches using only a small set of 2D images. Our method integrates 3D reconstruction, neural rendering, adversarial patch optimization, and an object detection victim model into a unified workflow. We use 2D Gaussian Splatting for high fidelity mesh reconstruction and FlexPara for surface parameterization that produces texture maps suitable for patch editing. Together, these components form a fully differentiable pipeline in PyTorch3D that links texture modification to model outputs, enabling efficient optimization of patches that remain effective across many viewpoints. The complete process, from image capture to patch printing and physical evaluation, can be completed within a few hours. We demonstrate the effectiveness of the resulting patches through attacks on the YOLOv8 object detection model and discuss remaining challenges and opportunities for improving robustness and scalability.

Singhvi, Vivaan [ORNL] (ORCID:0009000586288221)↗
DOE OSTI2026

Time at Temperature Experiments on Neutron Irradiated Zircaloy-2 using Conventional & Flash DSC

Dryout events in Boiling Water Reactors (BWRs) are currently treated by NRC regulations as automatic disqualification for continued fuel rod operation, even though this criterion does not account for the rate or duration of power increases, associated changes in material behavior, or the possibility of rewetting. Operational history from Anticipated Operational Occurrences (AOOs) shows that short, transient power excursions often demand only modest heat removal, and industry experience suggests that fuel can briefly enter dryout yet return to safe, stable operation. The lack of detailed understanding of the material response during such events motivates the present series of experiments. This study uses unique and innovative methods to investigate microstructural evolution in irradiated Zircaloy-2 exposed to high temperatures in inert environments. Differential Scanning Calorimetry (DSC) and FlashDSC are combined to build a comprehensive experimental framework capable of identifying the a–ß phase transformation in zirconium and examining defect annealing under steep thermal gradients. FlashDSC enables rapid heating and cooling of Focused Ion Beam (FIB)–prepared large-area lift-outs (LALOs) of irradiated Zircaloy-2 at rates of 1,000 K/s to peak temperatures of 600°C, 750°C, and 900°C. The goal is to determine whether these conditions produce measurable microstructural changes that could influence cladding performance in typical BWR environments. Microstructural characterization includes quantifying dislocation density and assessing secondary phase particle (SPP) size and distribution using Transmission Electron Microscopy (TEM). Ongoing analysis, such as diffraction pattern indexing and 4D STEM processing, will further refine these observations.

36 - MATERIALS SCIENCE↗