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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Paleomagnetic Evidence for a Partially Differentiated Ordinary Chondrite Parent Asteroid

Abstract The textures and accretion ages of chondrites have been used to argue that their parent asteroids never differentiated. Without a core, undifferentiated planetesimals could not have generated magnetic fields through dynamo activity, so chondrites are not expected to have experienced such fields. However, the magnetic remanence carried by the CV chondrites is consistent with dynamo‐generated fields, hinting that partially differentiated asteroids consisting of an unmelted crust atop a differentiated interior may exist. Here, we test this hypothesis by applying synchrotron X‐ray microscopy to metallic veins in the slowly cooled H6 chondrite Portales Valley. The magnetic remanence carried by nanostructures in these veins indicates that this meteorite recorded a magnetic field over a period of tens to hundreds of years at ∼100 Myr after solar system formation. These properties are inconsistent with external field sources such as the nebula, solar wind, or impacts, but are consistent with dynamo‐generated fields, indicating that the H chondrite parent body contained an advecting metallic core and was therefore partially differentiated. We calculate the thermal evolution of the chondritic portions of partially differentiated asteroids that form through incremental accretion across 10 5 to 10 6 years, finding this can agree with the measured ages and cooling rates of multiple H chondrites. We also predict that the cores of these bodies could have been partially liquid and feasibly generating a dynamo at 100 Myr after solar system formation. These observations contribute to a growing body of evidence supporting a spectrum of internal differentiation within some asteroids with primitive surfaces.

Bryson, J. F. J.↗
DOE OSTI2023

Differentiable modelling to unify machine learning and physical models for geosciences

Process-based modelling offers interpretability and physical consistency in many domains of geosciences but struggles to leverage large datasets efficiently. Machine-learning methods, especially deep networks, have strong predictive skills yet are unable to answer specific scientific questions. Here, in this Perspective, we explore differentiable modelling as a pathway to dissolve the perceived barrier between process-based modelling and machine learning in the geosciences and demonstrate its potential with examples from hydrological modelling. ‘Differentiable’ refers to accurately and efficiently calculating gradients with respect to model variables or parameters, enabling the discovery of high-dimensional unknown relationships. Differentiable modelling involves connecting (flexible amounts of) prior physical knowledge to neural networks, pushing the boundary of physics-informed machine learning. It offers better interpretability, generalizability, and extrapolation capabilities than purely data-driven machine learning, achieving a similar level of accuracy while requiring less training data. Additionally, the performance and efficiency of differentiable models scale well with increasing data volumes. Under data-scarce scenarios, differentiable models have outperformed machine-learning models in producing short-term dynamics and decadal-scale trends owing to the imposed physical constraints. Differentiable modelling approaches are primed to enable geoscientists to ask questions, test hypotheses, and discover unrecognized physical relationships. Future work should address computational challenges, reduce uncertainty, and verify the physical significance of outputs.

58 GEOSCIENCES↗
DOE OSTI2026

Study of fully coupled three-dimensional envelope instability using automatic differentiation

Automatic differentiation is a powerful tool for computing derivatives of simulation results with respect to given parameters. In this Letter, we have applied this tool to investigate the instability of a dynamical system governed by 21 ordinary differential equations. This second-order instability (named envelope instability) is driven by space-charge effects and has a significant impact on the operational regimes of particle accelerators. Our study delves into the three-dimensional envelope instability, incorporating both transverse and longitudinal coupling. Conventionally, analyzing this complex system would necessitate solving 441 ordinary differential equations, which is computationally intractable. However, by employing automatic differentiation, we were able to track only 21 equations. This approach allowed us to uncover an additional instability stopband, which arises from space-charge-induced coupling and has not been reported in previous studies. This research highlights the significant advantages of automatic differentiation in analyzing complicated dynamical systems involving a large number of ordinary differential equations.

Qiang, Ji [Lawrence Berkeley National Laboratory (↗
DOE OSTI2021

Informing nuclear physics via machine learning methods with differential and integral experiments

Information from differential nuclear-physics experiments and theory is often too uncertain to accurately define nuclear-physics observables such as cross sections or energy spectra. Integral experimental data, representing the applications of these observables, are often more precise but depend simultaneously on too many of them to unambiguously identify issues in the observable with human expert analysis alone. Here, we explore how we can leverage physics knowledge gained from differential experimental data, nuclear theory, integral experiments, and neutron-transport calculations to better understand nuclear-physics observables in the context of the application area represented by integral experiments. We support this task with machine-learning methods to discern trends in a large amount of convoluted data. Differential and integral information was used in an analysis augmented by the random forest and the Shapley additive explanations metric. We chose as an application area one that is represented by criticality measurements and pulsed-sphere neutron-leakage spectra. We show one representative example ( 241 Pu fission observables) where the combination of differential and integral information allowed to resolve issues in data representing these observables. As a starting point, the machine learning (ML) algorithms highlighted several observables as leading potentially to bias in simulating integral experiments. Differential information, paired with sensitivity to integral quantities, allowed us then to pinpoint one specific observable ( 241 Pu fission cross section) as the main driver of bias. The comparison to integral experiments, on the other hand, allowed us to indicate a likely reliable experiment among several discrepant ones for this observables. In other cases (e.g., 239 Pu observables), we were not able to resolve the confounding introduced by integral experiments but instead highlighted the need for targeted new experiments and theory developments to better constrain the nuclear-physics space for the application area represented by integral experiments. We were able to combine information from differential experimental data, nuclear-physics theory, integral experiments, and neutron-transport simulations of the latter experiments with the help of the random forest algorithm and expert judgment. This combination of knowledge allows to improve our description of nuclear-physics observables as applied to a particular application area.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

The Memory Scaling of Reverse-Mode Differentiation in Particle Accelerator Simulations with Space Charge

The recent development of differentiable simulation codes for particle accelerators has enabled gradient-based workflows that promise finer control and more realistic modeling of accelerator facilities. However, when using reverse-mode automatic differentiation, the memory usage continuously increases during the simulation, and can potentially exceed the available hardware memory - especially when costly space charge computation is included. To study the memory requirements for differentiable simulations, we have implemented space charge in Cheetah, a PyTorch-based beam tracking code that supports reverse-mode differentiation. We find that the memory usage for reverse-mode differentiation grows linearly with the number of macroparticles and cells, and that it is proportional to the number of space charge kicks involved in the simulation. This general scaling can be used to evaluate whether a given differentiable simulation is feasible given hardware memory constraints.

Dhamrait, Arjun↗
DOE OSTI2019

Application of symmetries to differential equations: symmetry reduction and solution transformation examples

In the following lecture notes, we provide an outline of what we mean by symmetries of differential equations and illustrate their utility in simplification of differential equations which may lead to solutions. Loosely, a symmetry of a differential equation is a local group of transformations acting on the space of independent and dependent variables that transform solutions of a differential equation to new solutions, i.e., if Δ(x,u,p)=0 denotes a differential equation with a solution prescribed by [x,u,p], then knowlege of a symmetry facilitates the attainment of a new solution, [$\tilde{x}$, $\tilde{u}$, $\tilde{p}$], such that Δ(x,u,p)=Δ($\tilde{x}$, $\tilde{u}$, $\tilde{p}$)=0. In Equations (1) and (2), x, u and p denote the independent variables, dependent variables and partial derivatives, respectively and the tildes denote the transformations. Equation (2) will subsequently be referred to as the symmetry criterion. We can apply the knowlede of the existence of a symmetry in two ways: 1) If not solution is known apriori, the symmetry can be applied to simplify the differential equation through the method of symmetry reduction. Often, the resultant simplification may enable a solution to be determined; 2) If we already possess a solution to a different equation, knowledge of an admissible symmetry enables us to generate new solutions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1978

Absolute elastic differential electron scattering cross sections in the intermediate energy region. IV - CO

Using a crossed electron beam-molecular beam scattering geometry and a relative-flow technique, ratios of elastic differential cross sections of CO to those of He have been measured at electron impact energies of 3, 5, 7.5, 9.9, 15, 20, 30, 50, 75, and 100 eV. At each energy, an angular range of 15 to 130 deg has been covered. These ratios have been multiplied by previously known He elastic differential cross sections to obtain elastic differential cross sections for CO. Since pure rotational excitations were not resolved, the elastic differential cross sections are a sum of elastic and pure rotational excitations at room temperature. From a knowledge of differential cross sections (DCS), integral and momentum transfer cross sections have been calculated. Both the DCS and integral cross sections are compared at 50, 75, and 100 eV to a recent two-potential theory of e-molecule scattering. Present results show that the isoelectronic molecules CO and N2 have very similar magnitudes and shapes of their differential cross sections.

Tanaka, H.↗
NASA NTRS1985

Differentiating matrices for arbitrarily spaced grid points

Differentiating matrices allow the numerical differentiation of functions defined at points of a discrete grid. A type of differentiating matrix based on local approximation on a sequence of sliding subgrids is considered. Previous derivations of this type of matrix have been restricted to grids with uniformly spaced points, and the resulting derivative approximations have lacked precision, especially at endpoints. The new formulation allows grids which have arbitrarily space points. It is shown that high accuracy can be achieved through use of differentiating matrices on non-uniform grids which include near-boundary points. Use of the differentiating matrix as an operator to solve eigenvalue problems involving ordinary differential equations is also considered.

Lakin, W. D.↗
NASA NTRS1986

Differentiating matrices for arbitrarily spaced grid points

Differentiating matrices allow the numerical differentiation of functions defined at points of a discrete grid. Previous derivations of these matrices have been restricted to grids with uniformly spaced points, and the resulting derivative approximations have lacked precision, especially at endpoints. The present work derives differentiating matrices on grids with arbitrarily spaced points. It is shown that high accuracy can be achieved through use of differentiating matrices on non-uniform grids through the expedient of including 'near boundary' points. Use of the differentiating matrix as an operator in the solution of problems involving ordinary differential equations is also considered.

Lakin, W. D.↗
NASA NTRS1991

A unified approach to the helioseismic forward and inverse problems of differential rotation

A general, degenerate perturbation theoretic treatment of the helioseismic forward and inverse problem for solar differential rotation is presented. For the forward problem, differential rotation is represented as the axisymmetric component of a general toroidal flow field using velocity spherical harmonics. This approach allows each degree of differential rotation to be estimated independently from all other degrees. In the inverse problem, the splitting caused by differential rotation is expressed as an expansion in a set of orthonormal polynomials that are intimately related to the solution of the forward problem. The combined use of vector spherical harmonics as basis functions for differential ratio and the Clebsch-Gordon coefficients to represent splitting provides a unified approach to the forward and inverse problems of differential rotation which greatly simplify inversion.

Ritzwoller, Michael H.↗
NASA NTRS1992

Differential phase acoustic microscope for micro-NDE

A differential phase scanning acoustic microscope (DP-SAM) was developed, fabricated, and tested in this project. This includes the acoustic lens and transducers, driving and receiving electronics, scanning stage, scanning software, and display software. This DP-SAM can produce mechanically raster-scanned acoustic microscopic images of differential phase, differential amplitude, or amplitude of the time gated returned echoes of the samples. The differential phase and differential amplitude images provide better image contrast over the conventional amplitude images. A specially designed miniature dual beam lens was used to form two foci to obtain the differential phase and amplitude information of the echoes. High image resolution (1 micron) was achieved by applying high frequency (around 1 GHz) acoustic signals to the samples and placing two foci close to each other (1 micron). Tone burst was used in this system to obtain a good estimation of the phase differences between echoes from the two adjacent foci. The system can also be used to extract the V(z) acoustic signature. Since two acoustic beams and four receiving modes are available, there are 12 possible combinations to produce an image or a V(z) scan. This provides a unique feature of this system that none of the existing acoustic microscopic systems can provide for the micro-nondestructive evaluation applications. The entire system, including the lens, electronics, and scanning control software, has made a competitive industrial product for nondestructive material inspection and evaluation and has attracted interest from existing acoustic microscope manufacturers.

Waters, David D.↗
NASA NTRS1996

Differential Canard deflection for generation of yawing moment on the X-31 with and without the vertical tail

The feasibility of augmenting the available yaw control power on the X-31 through differential deflection of the canard surfaces was studied as well as the possibility of using differential canard control to stabilize the X-31 with its vertical tail removed. Wind-tunnel tests and the results of departure criteria and linear analysis showed the destabilizing effect of the reduction of the vertical tail on the X-31. Wind-tunnel testing also showed that differential canard deflection was capable of generating yawing moments of roughly the same magnitude as the thrust vectoring vanes currently in place on the X-31 in the post-stall regime. Analysis showed that the X-31 has sufficient aileron roll control power that with the addition of differential canard as a yaw controller, the wind-axis roll accelerations will remain limited by yaw control authority. It was demonstrated, however, that pitch authority may actually limit the maximum roll rate which can be sustained. A drop model flight test demonstrated that coordinated, wind axis rolls could be performed with roll rates as high as 50 deg/sec (full scale equivalent) at 50 deg angle of attack. Another drop model test was conducted to assess the effect of vertical tail reduction, and an analysis of using differential canard deflection to stabilize the tailless X-31 was performed. The results of six-degree-of-freedom, non-linear simulation tests were correlated with the drop model flights. Simulation studies then showed that the tailless X-31 could be controlled at angles of attack at or above 20 deg using differential canard as the only yaw controller.

Whiting, Matthew Robert↗
NASA NTRS1988

Kinetics of determination in the differentiation of isolated mesophyll cells of Zinnia elegans to tracheary elements

Mechanically isolated mesophyll cells of Zinnia elegans L. cv Envy differentiate to tracheary elements when cultured in inductive medium containing 0.5 micromolar alpha-naphthaleneacetic acid and 0.5 micromolar benzyladenine. The cells do not differentiate when cultured in medium in which the concentration of auxin and/or cytokinin has been reduced to 0.005 micromolar. Cells require an initial 24-hour exposure to inductive cytokinin and 56-hour exposure to inductive auxin for differentiation at 72 hours of culture. Freshly isolated Zinnia cells can be maintained in medium having low concentrations of both auxin and cytokinin for only 1 day without significant loss of potential to differentiate upon transfer to inductive medium. Initial culture for up to 2 days in medium having high auxin and low cytokinin, or low auxin and high cytokinin, allows full differentiation on the third day after transfer to inductive medium and potentiates the early differentiation of some cells.

Non-NASA Center↗
NASA NTRS1998

Fibronectin is a survival factor for differentiated osteoblasts

The skeletal extracellular matrix produced by osteoblasts contains the glycoprotein fibronectin, which regulates the adhesion, differentiation and function of various adherent cells. Interactions with fibronectin are required for osteoblast differentiation in vitro, since fibronectin antagonists added to cultures of immature fetal calvarial osteoblasts inhibit their progressive differentiation. To determine if fibronectin plays a unique role in fully differentiated osteoblasts, cultures that had already formed mineralized nodules in vitro were treated with fibronectin antagonists. Fibronectin antibodies caused >95% of the cells in the mature cultures to display characteristic features of apoptosis (nuclear condensation, apoptotic body formation, DNA laddering) within 24 hours. Cells appeared to acquire sensitivity to fibronectin antibody-induced apoptosis as a consequence of differentiation, since antibodies failed to kill immature cells and the first cells killed were those associated with mature nodules. Intact plasma fibronectin, as well as fragments corresponding to the amino-terminal, cell-binding, and carboxy-terminal domains of fibronectin, independently induced apoptosis of mature (day-13), but not immature (day-4), osteoblasts. Finally, transforming growth factor-beta1 partially protected cells from the apoptotic effects of fibronectin antagonists. Thus, in the course of maturation cultured osteoblasts switch from depending on fibronectin for differentiation to depending on fibronectin for survival. These data suggest that fibronectin, together with transforming growth factor-beta1, may affect bone formation, in part by regulating the survival of osteoblasts.

NASA Center ARC↗
DOE OSTI2017

Calculating Control Rod Differential Worth Using Adjoint-Weighted Method

Recently, there are more and more research on using Monte Carlo method to calculate reactor important parameters. Control rod worth is one of the important parameters as it can determine the critical control rod position. There are two presentations of the control rod worth for utilization purpose, integral and differential. The former is the total reactivity change, while the latter is the reactivity change caused by per unit movement. In Monte Carlo calculation, differential control rod worth of a particular position is estimated by finite difference of reactivity between a small gap divided by gap distance. Though easy to implemented, the problem of this approach lies in that too large gap will enlarge error between 'differential' and 'difference', while too small gap will be disturbed by error caused by random fluctuation. Monte Carlo perturbation theory has been used to calculate the k{sub eff} sensitivity to certain parameters for a long time. By calculating the differential coefficient along with the k{sub eff}, this method is quite time-consuming and relatively accurate compared to finite difference method. Nevertheless, most of perturbation research focus on the perturbation of nuclide density or nuclear data. For geometrical perturbation, limited research was done before. Takahashi estimated the k{sub eff} sensitivity to reflector movement of pulsed fast reactor. MCNP6 used adjoint-weighted method for Godiva, Zeus and uranium oxyfluoride solution container benchmarks. This paper tries to apply the perturbation method to the determination of control rod differential worth. In this paper, based on RMC, a continuous-energy Reactor Monte Carlo transport code, Iterated Fission Probability method (IFP) is researched to calculate control rod differential worth on assembly level. IFP was originally developed by Hurwitz, and has been adopted by most of Monte Carlo code to compute adjoint-weighted parameters. This paper is organized with four parts. The first is introduction. Adjoint-weighted method is introduced in the second part. Problem 4 of benchmark from CASL is built to verify the method. Conclusion is given in final part.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Regularized Differentiation for Bioburden Density Estimation in Planetary Protection

In this paper, we propose and investigate the performance of two novel shrinkage estimators for bioburden density estimation in planetary protection. The estimators are based on the regularized differentiation of a cumulative count of colony forming units collected throughout the data collecting session or the life cycle of the entire mission. The regularized differentiation recasts the problem of bioburden density estimation as a linear least squares problem. The least squares problem is then solved through regularization techniques, such as truncated singular value decomposition and penalized least squares. The regularization is necessary to avoid noise amplification during the differentiation of noisy data. The two regularization estimators are compared with four other commonly used estimators to simultaneously evaluate the means of multivariable independent Poisson distributions: the maximum likelihood, noninformative Bayes estimator with Jeffreys prior, Empirical Bayes using conjugate gamma-Poisson model with gamma parameters selected by method of moments, and the Clevenson-Zidek estimator. It is shown through computer-simulated data that the regularized differentiation based on ridge regression has the smallest mean-squared error among all estimators. The analysis of shrinkage mechanism implemented by regularized differentiation is performed, and it is shown that the regularized differentiation amounts to performing a weighted averaging of all the samples. The weights are determined by the regularization parameter automatically selected by the L-curve technique. Since the method of least squares makes no distributional assumptions about the data, it presents an attractive technique for bioburden density estimation when there are concerns about the misspecification of the distributional model. The paper concludes with the analysis of the bioburden data collected during InSight mission and directions for future work.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2018

On the periodic solutions of the Michelson continuous and discontinuous piecewise linear differential system

Applying new results from the averaging theory for continuous and discontinuous differential systems, we study the periodic solutions of two distinct versions of the Michelson differential system: a Michelson continuous piecewise linear differential system and a Michelson discontinuous piecewise linear differential system. The tools here used can be applied to general nonsmooth differential systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Chapter 4 - Recent Advances in Identification of Differential Equations from Noisy Data: IDENT Review

Differential equations and numerical methods are extensively used to model various real-world phenomena in science and engineering. With modern developments, we aim to find the underlying differential equation from a single observation of time-dependent data. If we assume that the differential equation is a linear combination of various linear and nonlinear differential terms, then the identification problem can be formulated as solving a linear system. The goal then reduces to finding the optimal coefficient vector that best represents the time derivative of the given data. We review some recent works on the identification of differential equations. We find some common themes for the improved accuracy: (i) The formulation of linear system with proper denoising is important, (ii) how to utilize sparsity and model selection to find the correct coefficient support needs careful attention, and (iii) there are ways to improve the coefficient recovery. We present an overview and analysis of recent developments on the topic.

97 MATHEMATICS AND COMPUTING↗