Search NASA⌕ Search

SEARCH · Search NASA

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Analysis of differential scanning calorimetry data for aged plutonium

Differential scanning calorimetry data for samples of a 52 year old plutonium alloy with 3.3 at. % Ga that were heated beyond the melting point is analyzed using transition state theory to find activation energies for the δ to ε and ε to liquid phase transitions. A Bayesian statistical method involving a Gaussian process model is used to find mean values and confidence intervals for the activation energies. The activation energy for the δ to ε phase transition increases by 3.3 ± 3.8% per decade, relative to the case when all age related plutonium lattice point defects have been removed through annealing. The corresponding increase in activation energy for the ε to liquid transition is shown to be 7.1 ± 1.8% per decade. It is postulated that the change in activation energy with age for both phase transitions is caused, in part, by the accumulation of the same type of lattice point defects associated with the observed increase in elastic bulk modulus over time.

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

Stable perturbations of linear differential equations generating a uniformly bounded group

Stability problems for solutions of the differential equation u{sup ′}(t)=Au+εB(t,u) in a Banach space are considered. It is assumed that for ε=0 this equation generates a uniformly bounded group of class C{sub 0}. Sufficient conditions on B and A are found under which the solutions of this equation are bounded for small ε. A linearization principle is proved for this equation under certain conditions on the operator B. Bibliography: 9 titles. (paper)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Direct Numerical Simulation of Turbulent Flow Inside a Differentially Heated Composite Cavity

In this work, turbulent flow inside a modified differentially heated cavity at high Rayleigh number (Ra ~ 10 9 ) has been studied through fully resolved direct numerical simulation (DNS) using the high-order spectral element method code Nek5000. The flow configuration includes two separate physical phenomena: the natural recirculation itself, and the flow inside a curved channel. Simulations have been carried out using both the Boussinesq approximation and the low-Mach compressible formulation. Significant discrepancies between the two methods inform of the extreme caution that should be exercised when using the Boussinesq approximation in the limits of its applicability. The DNS solutions are analyzed in terms of polynomial-order convergence and Reynolds stress budgets, and the turbulence quantities and velocity profiles are presented as a reference for the validation of turbulence models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

A physics-constrained neural ordinary differential equations approach for robust learning of stiff chemical kinetics

The high computational cost associated with solving for detailed chemistry poses a significant challenge for predictive computational fluid dynamics (CFD) simulations of turbulent reacting flows. While deep learning techniques have been explored to develop faster surrogate models, they often fail to integrate reliably with CFD solvers. This instability arises because traditional deep learning approaches optimize for training error without ensuring compatibility with ordinary differential equation (ODE) solvers, resulting in accumulation of errors over time. Recently, neuralODE (NODE) based approaches have been shown to be a promising technique to emulate and accelerate detailed chemistry computations. Here, in the present work, we extend this NODE framework for stiff chemical kinetics by incorporating mass conservation constraints directly into the loss function during training. This ensures that the total mass as well as the individual elemental species masses are conserved in an a-posteriori manner. Proof-of-concept studies are performed with the novel physics-constrained NODE (PC-NODE) approach for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions. It is demonstrated that the PC-NODE framework not only improves the physical consistency of the resulting data-driven model with respect to mass conservation criteria, but also improves training efficiency. PC-NODE is shown to achieve 2–100× speedup relative to the hydrogen-air detailed chemical mechanism depending on the type of the ODE solver (implicit or explicit) used during autoregressive inference tests. Lastly, a-posteriori studies are performed wherein the trained PC-NODE model is coupled with a CFD solver. It is shown that higher accuracy is achieved with PC-NODE relative to the purely data-driven NODE approach. Moreover, PC-NODE also exhibits robustness and generalizability to unseen initial conditions from within (interpolative capability) as well as outside (extrapolative capability) the training regime.

computational combustion↗
DOE OSTI2016

Design and fabrication of a differential scanning nanocalorimeter

This paper describes the design, fabrication, and characterization of a differential scanning nanocalorimeter that significantly reduces the sample volume to microliters and can potentially improve the temperature sensitivity to 10 µK. The nanocalorimeter consists of a polymeric freestanding membrane, four high-sensitive low-noise thermistors based on silicon carbide (SiC), and a platinum heater and temperature sensor. With the integrated heater and sensors, temperature scanning and power compensation can be achieved for calorimetric measurement. Temperature sensing SiC film was prepared by using sintered SiC target and DC magnetron sputtering under different gas pressures and sputtering power. The SiC sensing material is characterized through the measurement of current–voltage curves and noise levels. The thermal performance of a fabricated nanocalorimeter is studied in simulation and experiment. The experiment results show the device has excellent thermal isolation to hold thermal energy. As a result, the noise test together with the simulation show the device is promising for micro 10 µK temperature sensitivity and nanowatt resolution which will lead to low-volume ultra-sensitive nanocalorimetry for biological processes, such as protein folding and ligand binding.

42 ENGINEERING↗
DOE OSTI2019

Negative differential resistance and quantum oscillations in FeSb 2 with embedded antimony

We present a systematical study on single crystalline FeSb 2 using electrical transport and magnetic torque measurements at low temperatures. Nonlinear magnetic field dependence of Hall resistivity demonstrates a multi-carrier transport instinct of the electronic transport. Current-controlled negative differential resistance (CC-NDR) observed in current–voltage characteristics below ~ 7 K is closely associated with the intrinsic transition ~ 5 K of FeSb 2 , which is, however, mediated by extrinsic current-induced Joule heating effect. The antimony crystallized in a preferred orientation within the FeSb 2 lattice in the high-temperature synthesis process leaves its fingerprint in the de Haas-Van Alphen (dHvA) oscillations, and results in the regular angular dependence of the oscillating frequencies. Nevertheless, possible existence of intrinsic non-trivial states cannot be completely ruled out. Our findings call for further theoretical and experimental studies to explore novel physics on flux-free grown FeSb 2 crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Optimized finite-build stellarator coils using automatic differentiation

A new stellarator coil design code is introduced that optimizes the position and winding pack orientation of finite-build coils. The new code, called flexible optimized curves in space using automatic differentiation (AD) and finite build (FOCUSADD), performs gradient-based optimization in a high-dimensional, non-convex space. The derivatives with respect to parameters of finite-build coils are easily and efficiently computed using AD. FOCUSADD parametrizes coil positions in free space using a Fourier series and uses a multi-filament approximation to the coil winding pack. The orientation of the winding pack is parametrized with a Fourier series and can be optimized as well. Optimized finite-build coils for a Wendelstein 7-X (W7-X)-like stellarator are found, and compared with filamentary coil results. The final positions of optimized finite-build W7-X-like coils are shifted, on average, by approximately 2.5 mm relative to optimized filamentary coils. These results suggest that finite-build effects should be accounted for in the optimization of stellarators with low coil tolerances.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Flow annealed importance sampling bootstrap meets differentiable particle physics

High-energy physics requires the generation of large numbers of simulated data samples from complex but analytically tractable distributions called matrix elements. Surrogate models, such as normalizing flows, are gaining popularity for this task due to their computational efficiency. We adopt an approach based on flow annealed importance sampling bootstrap (FAB) that evaluates the differentiable target density during training and helps avoid the costly generation of training data in advance. We show that FAB reaches higher sampling efficiency with fewer target evaluations in high dimensions in comparison to other methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Spatiotemporal control of laser intensity using differentiable programming

Optical techniques for spatiotemporal control can produce laser pulses with custom amplitude, phase, or polarization structure. In nonlinear optics and plasma physics, the use of structured pulses typically follows a forward design approach, in which the efficacy of a known structure is analyzed for a particular application. Inverse approaches, in contrast, enable the discovery of new structures with the potential for superior performance. Here, an implementation of the unidirectional pulse propagation equation that supports automatic differentiation is combined with gradient-based optimization to design structured pulses with features that are advantageous for a range of nonlinear optical and plasma-based applications: (1) a longitudinally uniform intensity over an extended region, (2) a superluminal intensity peak that travels many Rayleigh ranges with constant duration, spot size, and amplitude, and (3) a laser pulse that ionizes a gas to form a uniform column of plasma. In the final case, optimizing the full spatiotemporal structure improves the performance by a factor of 15 compared to optimizing only spatial or only temporal structure, highlighting the advantage of spatiotemporal control.

automatic differentiation↗
DOE OSTI2025

From disorganized data to emergent dynamic models: Questionnaires to partial differential equations

Starting with sets of disorganized observations of spatially varying and temporally evolving systems, obtained at different (also disorganized) sets of parameters, we demonstrate the data-driven derivation of parameter dependent, evolutionary partial differential equation (PDE) models capable of generating the data. This tensor type of data is reminiscent of shuffled (multidimensional) puzzle tiles. The independent variables for the evolution equations (their “space” and “time”) as well as their effective parameters are all emergent , i.e. determined in a data-driven way from our disorganized observations of behavior in them. We use a diffusion map based questionnaire approach to build a smooth parametrization of our emergent space/time/parameter space for the data. This approach iteratively processes the data by successively observing them on the “space,” the “time” and the “parameter” axes of a tensor. Once the data become organized, we use machine learning (here, neural networks) to approximate the operators governing the evolution equations in this emergent space. Our illustrative examples are based (i) on a simple advection–diffusion model; (ii) on a previously developed vertex-plus-signaling model of Drosophila embryonic development; and (iii) on two complex dynamic network models (one neuronal and one coupled oscillator model) for which no obvious smooth embedding geometry is known a priori. This allows us to discuss features of the process like symmetry breaking, translational invariance, and autonomousness of the emergent PDE model, as well as its interpretability.

generative models↗
DOE OSTI2024

A differentiable simulation package for performing inference of synchrotron-radiation-based diagnostics

The direction of particle accelerator development is ever-increasing beam quality, currents and repetition rates. This poses a challenge to traditional diagnostics that directly intercept the beam due to the mutual destruction of both the beam and the diagnostic. An alternative approach is to infer beam parameters non-invasively from the synchrotron radiation emitted in bending magnets. However, inferring the beam distribution from a measured radiation pattern is a complex and computationally expensive task. To address this challenge we present SYRIPY ( SYnchrotron Radiation In PYthon ), a software package intended as a tool for performing inference of synchrotron-radiation-based diagnostics. SYRIPY has been developed using PyTorch , which makes it both differentiable and able to leverage the high performance of GPUs, two vital characteristics for performing statistical inference. The package consists of three modules: a particle tracker, Lienard–Wiechert solver and Fourier optics propagator, allowing start-to-end simulation of synchrotron radiation detection to be carried out. SYRIPY has been benchmarked against SRW , the prevalent numerical package in the field, showing good agreement and up to a 50× speed improvement. Finally, we have demonstrated how SYRIPY can be used to perform Bayesian inference of beam parameters using stochastic variational inference.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

A Semi-Analytical Solution Approach for Solving Constant-Coefficient First-Order Partial Differential Equations

Simulation and control of many dynamic systems involve solving partial differential equations (PDE). This letter proposes a semi-analytical solution (SAS) approach for fast and high-quality solution of first-order PDEs. The region of interest of the studied PDE is divided into a grid, and an SAS is derived for each grid cell in the form of the multivariate polynomials, of which the coefficients are identified using initial value and boundary value conditions. The solutions are solved in a “time-stepping” manner, i.e. within one time step, the coefficients of the SAS are identified and the initial value of the next time step is evaluated. This approach achieves a significantly larger grid cell than the widely used finite difference method, and thus enhances the computational efficiency significantly. Furthermore, the simulation result on the natural gas pipeline model demonstrates the advantages of SAS in accuracy and computational efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Dynamic Parameter Estimation with Physics-based Neural Ordinary Differential Equations

Accurate estimation of dynamic parameters of gen-erators is crucial to building a reliable model for dynamical studies and reliable operation of the power system. This paper develops a physics-based neural ordinary differential equations (ODE) approach to learn the parameters of generator dynamic model using phasor measurement units (PMU) data. We design a physics-based neural network to represent the swing equations of the power system dynamics. A loss function is defined as the difference between dynamic simulation results from the physics-based neural networks and pseudo PMU measurements. The parameters of generator dynamic model are iteratively updated using the neural ODEs and the adjoint method. By exploiting the mini-batch scheme in neural ODE training, the parameter estimation performance is significantly improved. Numerical study results on a 3-machine 9-bus system show that the proposed algorithm outperforms state-of-the-art baseline method in both computation time and dynamic parameter estimation accuracy.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Optimization of Distribution Feeder Topology: A Differential Programming Learning Approach

This paper presents a gradient based method for optimizing distribution feeder network topology under load un- certainty. We recast the optimal network reconfiguration problem as a learning problem where edge weights of a graph are learned to produce an optimized spanning tree for a distribution network. Using recent methods published on differentiable programming, we provide a data driven method for learning these weights. We test our method on 100 variations of an IEEE 15-bus test system. Our results show that our method outperforms more traditional mathematical programming-based approaches.

differentiable programming↗
DOE OSTI2021

The UPR regulator IRE1 promotes balanced organ development by restricting TOR-dependent control of cellular differentiation in Arabidopsis

Proteostasis of the endoplasmic reticulum (ER) is controlled by sophisticated signaling pathways that are collectively called the unfolded protein response (UPR) and are initiated by specialized ER membrane-associated sensors. The evidence that complete loss-of-function mutations of the most conserved of the UPR sensors, inositol-requiring enzyme 1 (IRE1), dysregulates tissue growth and development in metazoans and plants raises the fundamental question as to how IRE1 is connected to organismal growth. To address this question, we interrogated the Arabidopsis primary root, an established model for organ development, using the tractable Arabidopsis IRE1 mutant ire1a ire1b, which has marked root development defects in the absence of exogenous stress. We demonstrate that IRE1 is required to reach maximum rates of cell elongation and root growth. We also established that in the actively growing ire1a ire1b mutant root tips the Target of Rapamycin (TOR) kinase, a widely conserved pro-growth regulator, is hyperactive, and that, unlike cell proliferation, the rate of cell differentiation is enhanced in ire1a ire1b in a TOR-dependent manner. By functionally connecting two essential growth regulators, these results underpin a novel and critical role of IRE1 in organ development and indicate that, as cells exit an undifferentiated state, IRE1 is required to monitor TOR activity to balance cell expansion and maturation during organ biogenesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Fully quantum algorithm for mesoscale fluid simulations with application to partial differential equations

Fluid flow simulations marshal our most powerful computational resources. In many cases, even this is not enough. Quantum computers provide an opportunity to speed up traditional algorithms for flow simulations. We show that lattice-based mesoscale numerical methods can be executed as efficient quantum algorithms due to their statistical features. This approach revises a quantum algorithm for lattice gas automata to reduce classical computations and state preparation at every time step. For this, the algorithm approximates the qubit relative phases and subtracts them at the end of each time step. Phases are evaluated using the iterative phase estimation algorithm and subtracted using single-qubit rotation phase gates. Further, this method optimizes the quantum resource required and makes it more appropriate for near-term quantum hardware. We also demonstrate how the checkerboard deficiency that the D1Q2 scheme presents can be resolved using the D1Q3 scheme. The algorithm is validated by simulating two canonical partial differential equations: the diffusion and Burgers' equations on different quantum simulators. We find good agreement between quantum simulations and classical solutions for the presented algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Differential Cross Section for Elastic Deuteron–Proton Scattering at the Energy of 700 MeV per Nucleon

The results of measurements of the differential cross section for elastic deuteron–proton scattering at the energy of 700 MeV per nucleon that were performed at the internal target station of the nuclotron at the Joint Institute for Nuclear Research (JINR, Dubna) are presented. These data were obtained at angles in the range between 70◦ and 120◦ in the c.m. frame. The angular dependence obtained in this way is compared with world-averaged data at similar energies and with the results of theoretical calculations performed within the relativistic theory of multiple scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Description of Differential Cross Sections for {sup 63}Cu + p Nuclear Reactions Induced by High-Energy Cosmic-Ray Protons

The results of calculation of {sup 63}Cu + p differential cross sections at incident-proton energies between 10 and 200 MeV and a comparative analysis of these results are presented as a continuation of the earlier work of our group on developing methods for calculating the contribution of nuclear reactions to radiative effects arising in the onboard spacecraft electronics under the action of high-energy cosmic-ray protons on {sup 63}Cu nuclei (generation of single-event upsets) and as a supplement to the earlier calculations performed on the basis of the TАLYS code in order to determine elastic- and inelastic-scattering cross sections and charge, mass, and energy distributions of recoil nuclei (heavy products of the {sup 63}Cu + p nuclear reaction). The influence of various mechanisms of the angular distributions of particles emitted in the {sup 63}Cu + p nuclear reaction is also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗