Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differentiability”

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 739 records · Page 41

Differentiable modeling and optimization of non-aqueous Li-based battery electrolyte solutions using geometric deep learning

Electrolytes play a critical role in designing next-generation battery systems, by allowing efficient ion transfer, preventing charge transfer, and stabilizing electrode-electrolyte interfaces. In this work, we develop a differentiable geometric deep learning (GDL) model for chemical mixtures, DiffMix, which is applied in guiding robotic experimentation and optimization towards fast charging battery electrolytes. In particular, we extend mixture thermodynamic and transport laws by creating GDL-learnable physical coefficients. We evaluate our model with mixture thermodynamics and ion transport properties, where we show improved prediction accuracy and model robustness of Diff-Mix than its purely data-driven variants. Furthermore, with a robotic experimentation setup, Clio, we improve ionic conductivity of electrolytes by over 18.8% within 10 experimental steps, via differentiable optimization built on DiffMix gradients. By combining GDL, mixture physics laws, and robotic experimentation, DiffMix expands the predictive modeling methods for chemical mixtures and enables efficient optimization in large chemical spaces.

25 - ENERGY STORAGE↗

Covariance Matrix of a Double-Differential Doppler-broadened Elastic Scattering Cross Section

Legendre moments of a double-differential Doppler-broadened elastic neutron scattering cross section on {sup 238}U are computed near the 6.67 eV resonance at temperature T = 10{sup 3} K up to angular order 14. A covariance matrix of these Legendre moments is computed as a functional of the covariance matrix of the elastic scattering cross section. A variance of double-differential Doppler-broadened elastic scattering cross section is computed from the covariance of Legendre moments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The two modes extension to the Berk-Breizman equation: Delayed differential equations and asymptotic solutions

The integro-differential Berk-Breizman (BB) equation, describing the evolution of particle-driven wave mode is transformed into a simple delayed differential equation form {nu}{partial_derivative}a({tau})/{partial_derivative}{tau}=a({tau}) -a{sup 2}({tau}- 1) a({tau}- 2). This transformation is also applied to the two modes extension of the BB theory. The obtained solutions are presented together with the derived asymptotic analytical solutions and the numerical results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High Bandwidth Differential Amplifier for Shock Experiments

We developed a high bandwidth differential amplifier for gas gun shock experiments. The circuit has a bandwidth up to 1 GHz, and is capable of measuring signals of ≤1.5 V with a common mode rejection of 250 V. Conductivity measurements of gas gun targets are measured by flowing high currents through the targets. The voltage is measured across the target using a technique similar to a four-point probe. Because of the design of the current source and load, the target voltage is approximately 250 V relative to ground. Since the expected voltage change in the target is < 1 V, the differential amplifier must have a large common mode rejection. Various amplifying designs are shown, although the increased amplification decreases bandwidth. Bench tests show that the amplifier can withstand significant common mode DC voltage and measure 10 ns, and 50 mV signals.

42 ENGINEERING↗

Fully quantum state-resolved inelastic scattering of NO(X) + Kr: Differential cross sections and product rotational alignment

Fully quantum state selected and resolved inelastic scattering of NO(X) by krypton has been investigated. Initial Λ-doublet state selection is achieved using an inhomogeneous hexapole electric field. Differential cross sections and even-moment polarization dependent differential cross sections have been obtained at a collision energy of 514 cm{sup −1} for both spin-orbit and parity conserving and changing collisions. Experimental results are compared with those obtained from quantum scattering calculations and are shown to be in very good agreement. Hard shell quantum scattering calculations are also performed to determine the effects of the different parts of the potential on the scattering dynamics. Comparisons are also made with the NO(X) + Ar system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electron collisions with phenol: Total, integral, differential, and momentum transfer cross sections and the role of multichannel coupling effects on the elastic channel

We report theoretical and experimental total cross sections for electron scattering by phenol (C{sub 6}H{sub 5}OH). The experimental data were obtained with an apparatus based in Madrid and the calculated cross sections with two different methodologies, the independent atom method with screening corrected additivity rule (IAM-SCAR), and the Schwinger multichannel method with pseudopotentials (SMCPP). The SMCPP method in the N{sub open}-channel coupling scheme, at the static-exchange-plus-polarization approximation, is employed to calculate the scattering amplitudes at impact energies ranging from 5.0 eV to 50 eV. We discuss the multichannel coupling effects in the calculated cross sections, in particular how the number of excited states included in the open-channel space impacts upon the convergence of the elastic cross sections at higher collision energies. The IAM-SCAR approach was also used to obtain the elastic differential cross sections (DCSs) and for correcting the experimental total cross sections for the so-called forward angle scattering effect. We found a very good agreement between our SMCPP theoretical differential, integral, and momentum transfer cross sections and experimental data for benzene (a molecule differing from phenol by replacing a hydrogen atom in benzene with a hydroxyl group). Although some discrepancies were found for lower energies, the agreement between the SMCPP data and the DCSs obtained with the IAM-SCAR method improves, as expected, as the impact energy increases. We also have a good agreement among the present SMCPP calculated total cross section (which includes elastic, 32 inelastic electronic excitation processes and ionization contributions, the latter estimated with the binary-encounter-Bethe model), the IAM-SCAR total cross section, and the experimental data when the latter is corrected for the forward angle scattering effect [Fuss et al., Phys. Rev. A 88, 042702 (2013)].

74 ATOMIC AND MOLECULAR PHYSICS↗

Existence of a coupled system of fractional differential equations

We manage the existence and uniqueness of a fractional coupled system containing Schrödinger equations. Such a system appears in quantum mechanics. We confirm that the fractional system under consideration admits a global solution in appropriate functional spaces. The solution is shown to be unique. The method is based on analytic technique of the fixed point theory. The fractional differential operator is considered from the virtue of the Riemann-Liouville differential operator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Differential geometric invariants for time-reversal symmetric Bloch-bundles: The “Real” case

Topological quantum systems subjected to an even (resp. odd) time-reversal symmetry can be classified by looking at the related “Real” (resp. “Quaternionic”) Bloch-bundles. If from one side the topological classification of these time-reversal vector bundle theories has been completely described in De Nittis and Gomi [J. Geom. Phys. 86, 303–338 (2014)] for the “Real” case and in De Nittis and Gomi [Commun. Math. Phys. 339, 1–55 (2015)] for the “Quaternionic” case, from the other side it seems that a classification in terms of differential geometric invariants is still missing in the literature. With this article and its companion [G. De Nittis and K. Gomi (unpublished)] we want to cover this gap. More precisely, we extend in an equivariant way the theory of connections on principal bundles and vector bundles endowed with a time-reversal symmetry. In the “Real” case we generalize the Chern-Weil theory and we show that the assignment of a “Real” connection, along with the related differential Chern class and its holonomy, suffices for the classification of “Real” vector bundles in low dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Symmetries of stochastic differential equations: A geometric approach

A new notion of stochastic transformation is proposed and applied to the study of both weak and strong symmetries of stochastic differential equations (SDEs). The correspondence between an algebra of weak symmetries for a given SDE and an algebra of strong symmetries for a modified SDE is proved under suitable regularity assumptions. This general approach is applied to a stochastic version of a two dimensional symmetric ordinary differential equation and to the case of two dimensional Brownian motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Symmetries of the Gas Dynamics Equations using the Differential Form Method

Here, a brief review of the theory of exterior differential systems and isovector symmetry analysis methods is presented in the context of the one-dimensional inviscid compressible flow equations. These equations are formulated as an exterior differential system with equation of state (EOS) closure provided in terms of an adiabatic bulk modulus. The scaling symmetry generators—and corresponding EOS constraints—otherwise appearing in the existing literature are recovered through the application and invariance under Lie derivative dragging operations.

97 MATHEMATICS AND COMPUTING↗

Learning effective stochastic differential equations from microscopic simulations: Linking stochastic numerics to deep learning

We identify effective stochastic differential equations (SDEs) for coarse observables of fine-grained particle- or agent-based simulations; these SDEs then provide useful coarse surrogate models of the fine scale dynamics. We approximate the drift and diffusivity functions in these effective SDEs through neural networks, which can be thought of as effective stochastic ResNets. The loss function is inspired by, and embodies, the structure of established stochastic numerical integrators (here, Euler–Maruyama and Milstein); our approximations can thus benefit from backward error analysis of these underlying numerical schemes. They also lend themselves naturally to “physics-informed” gray-box identification when approximate coarse models, such as mean field equations, are available. Existing numerical integration schemes for Langevin-type equations and for stochastic partial differential equations can also be used for training; we demonstrate this on a stochastically forced oscillator and the stochastic wave equation. Our approach does not require long trajectories, works on scattered snapshot data, and is designed to naturally handle different time steps per snapshot. We consider both the case where the coarse collective observables are known in advance, as well as the case where they must be found in a data-driven manner.

97 MATHEMATICS AND COMPUTING↗

Parallelizable adjoint stencil computations using transposed forward-mode algorithmic differentiation

Algorithmic differentiation (AD) is a tool for generating discrete adjoint solvers, which efficiently compute gradients of functions with many inputs, for example for use in gradient-based optimization. AD is often applied to large computations such as stencil operators, which are an important part of most structured-mesh PDE solvers. Stencil computations are often parallelized, for example by using OpenMP, and optimized by using techniques such as cache-blocking and tiling to fully utilize multicore CPUs and many-core accelerators and GPUs. Differentiating these codes with conventional reverse-mode AD results in adjoint codes that cannot be expressed as stencil operations and may not be easily parallelizable. They thus leave most of the compute power of modern architectures unused. We present a method that combines forward-mode AD and loop transformation to generate adjoint solvers that use the same memory access pattern as the original computation that they are derived from and can benefit from the same optimization techniques. The effectiveness of this method is demonstrated by generating a scalable adjoint CFD solver for multicore CPUs and Xeon Phi accelerators.

algorithmic differentiation↗

Differential cross sections for excitation of H 2 by low-energy electron impact

Experimental and theoretical differential cross sections (DCS) for the electron-impact excitation of molecular hydrogen to the $B{}^{1}{{\rm{\Sigma }}}_{u}^{+}$, $c{}^{3}{{\rm{\Pi }}}_{u}$, $a{}^{3}{{\rm{\Sigma }}}_{g}^{+}$, $C{}^{1}{{\rm{\Pi }}}_{u}$, and the $E(F){}^{1}{{\rm{\Sigma }}}_{g}^{+}$ states are presented at incident energies near to threshold. The experimental DCSs were taken at incident energies of 14, 15, 16 and 17.5 eV and for scattering angles from 10° to 130°. In conclusion, the theoretical DCSs are from the convergent close-coupling method which has recently successfully modeled differential electron scattering from H 2 when compared with available experiment at energies of 17.5 eV and above.

74 ATOMIC AND MOLECULAR PHYSICS↗

Differentiable physics-enabled closure modeling for Burgers’ turbulence

Abstract Data-driven turbulence modeling is experiencing a surge in interest following algorithmic and hardware developments in the data sciences. We discuss an approach using the differentiable physics paradigm that combines known physics with machine learning to develop closure models for Burgers’ turbulence. We consider the one-dimensional Burgers system as a prototypical test problem for modeling the unresolved terms in advection-dominated turbulence problems. We train a series of models that incorporate varying degrees of physical assumptions on an a posteriori loss function to test the efficacy of models across a range of system parameters, including viscosity, time, and grid resolution. We find that constraining models with inductive biases in the form of partial differential equations that contain known physics or existing closure approaches produces highly data-efficient, accurate, and generalizable models, outperforming state-of-the-art baselines. Addition of structure in the form of physics information also brings a level of interpretability to the models, potentially offering a stepping stone to the future of closure modeling.

97 MATHEMATICS AND COMPUTING↗

Alignment of low-dose X-ray fluorescence tomography images using differential phase contrast

X-ray fluorescence nanotomography provides unprecedented sensitivity for studies of trace metal distributions in whole biological cells. Dose fractionation, in which one acquires very low dose individual projections and then obtains high statistics reconstructions as signal from a voxel is brought together (Hegerl & Hoppe, 1976), requires accurate alignment of these individual projections so as to correct for rotation stage runout. It is shown here that differential phase contrast at 10.2 keV beam energy offers the potential for accurate cross-correlation alignment of successive projections, by demonstrating that successive low dose, 3 ms per pixel, images acquired at the same specimen position and rotation angle have a narrower and smoother cross-correlation function (1.5 pixels FWHM at 300 nm pixel size) than that obtained from zinc fluorescence images (25 pixels FWHM). The differential phase contrast alignment resolution is thus well below the 700 nm × 500 nm beam spot size used in this demonstration, so that dose fractionation should be possible for reduced-dose, more rapidly acquired, fluorescence nanotomography experiments.

differential phase contrast↗

MapsTorch : automatic differentiation for X-ray fluorescence data analysis

X-ray fluorescence (XRF) is a popular spectroscopy technique for elemental analysis. Spectrum fitting and parameter tuning are at the core of XRF analysis and are conventionally manually intensive, especially for synchrotron experiments involving large amounts of diverse samples. This work introduces the automatic differentiation (AD) technique to XRF and an open-source package called MapsTorch. By transforming an analytical model of the XRF spectrum into a differentiable computation graph with AD, MapsTorch enables robust optimization of parameters and elemental intensities. We evaluate MapsTorch by conducting computational experiments on a large number of historical synchrotron XRF datasets and compare its performance with the currently practiced fitting tool NLopt. The results show that MapsTorch consistently achieves high-quality fits and often leads to better fitting quality than NLopt, particularly in tasks such as initial spectrum fitting and elemental intensity refinement. The robust performance of MapsTorch paves the way for developing automated and high-throughput XRF data analysis workflows to handle the increasing data volumes expected from next-generation synchrotron facilities.

X-ray fluorescence↗

Differentially Private Generation of Social Networks via Exponential Random Graph Models

Many social networks contain sensitive relational information. One approach to protect the sensitive relational information while offering flexibility for social network research and analysis is to release synthetic social networks at a pre-specified privacy risk level, given the original observed network. In this work, we propose the DP-ERGM procedure that synthesizes networks that satisfy the differential privacy (DP) via the exponential random graph model (EGRM). We apply DP-ERGM to a college student friendship network and compare its original network information preservation in the generated private networks with two other approaches: differentially private DyadWise Randomized Response (DWRR) and Sanitization of the Conditional probability of Edge given Attribute classes (SCEA). The results suggest that DP-EGRM preserves the original information significantly better than DWRR and SCEA in both network statistics and inferences from ERGMs and latent space models. In addition, DP-ERGM satisfies the node DP, a stronger notion of privacy than the edge DP that DWRR and SCEA satisfy.

97 MATHEMATICS AND COMPUTING↗

Efficient GPU Implementation of Automatic Differentiation for Computational Fluid Dynamics

Many scientific and engineering applications require repeated calculations of derivatives of output functions with respect to input parameters. Automatic Differentiation (AD) is a method that automates derivative calculations and can significantly speed up code development. In Computational Fluid Dynamics (CFD), derivatives of flux functions with respect to state variables (Jacobian) are needed for efficient solutions of the nonlinear governing equations. AD of flux functions on graphics processing units (GPUs) is challenging as flux computations involve many intermediate variables that create high register pressure and require significant memory traffic because of the need to store the derivatives. This paper presents a forward-mode AD method based on multivariate dual numbers that addresses these challenges and simultaneously reduces the floating-point operation count. The dimension of the multivariate dual numbers is optimized for performance. The flux computations are restructured to minimize the number of temporary variables and reduce register pressure. For effective utilization of memory bandwidth, shared memory is used to store the local flux Jacobian. This AD implementation is compared with several other Jacobian implementations on an NVIDIA V100 GPU (V100). For three-dimensional perfect-gas compressible-flow equations implemented in a practical CFD code, the AD implementation of a flux Jacobian based on multivariate dual numbers of dimension 5 outperforms all other GPU AD implementations on V100. Its performance is comparable with the optimized hand-differentiated version. Finally, the implementation achieves 75% of the peak floating-point throughput and 61 % of the peak global device memory bandwidth usage.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗