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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A Quantum Approach for Implementing Fixed-Point Arithmetic in Solving Ordinary Differential Equations

Differential equations (DEs) serve as fundamental tools in mathematical modeling across scientific disciplines, yet classical numerical solvers face limitations with large-scale or computationally intensive problems. This study explores a quantum-inspired approach to solving DEs, combining quantum- inspired techniques with classical methods. It focuses on fixed- point arithmetic on quantum circuits, utilizing basic quantum gates to manipulate DE solutions. We expand upon the techniques introduced by Zanger et al. [Quantum, 5, 502 (2021)] by offering a precise computation for a fixed-point signed multiplication scheme, while also presenting a quantum circuit capable of executing the fixed-point division algorithm. We demonstrate the feasibility of our approach through the simulation of a linear Ordinary Differential Equation (ODE), where initial conditions and parameters are encoded into quantum circuits using fixed- point representation. By executing sequences of quantum gates mimicking numerical integration steps, we obtain approximate solutions to the ODE with specified fixed-point precision.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Hypercomplex Automatic Differentiation in the Eulerian Hydrocode PAGOSA

Enabling the computation of partial derivatives or sensitivities in production hydrocodes is beneficial for design, optimization, sensitivity analysis, and uncertainty quantification. Traditional finite difference approximations of these sensitivities are inefficient since convergence studies of the step size is required for each parameter of interest. For these reasons, HYPercomplex Automatic Differentiation (HYPAD) was implemented in the Eulerian hydrocode PAGOSA. HYPAD is analogous to forward-mode automatic differentiation except hypercomplex numbers (numbers with multiple imaginary parts) are used instead of dual numbers. Accurate partial derivatives can be computed of all state variables with respect to multiple input variables in a single run. The method was implemented using operator overloading to handle hypercomplex algebra. HYPAD was demonstrated and verified on Sod’s shock tube problem to compute derivatives of the state variables with respect to a material parameter, initial conditions, and geometry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

The Poisson tensor completion non-parametric differential entropy estimator

We introduce the Poisson tensor completion (PTC) estimator, a non-parametric differential entropy estimator. The PTC estimator leverages inter-sample relationships to compute a low-rank Poisson tensor decomposition of the frequency histogram. Our crucial observation is that the histogram bins are an instance of a space partitioning of counts and thus can be identified with a spatial Poisson process. The Poisson tensor decomposition leads to a completion of the intensity measure over all bins—including those containing few to no samples—and leads to our proposed PTC differential entropy estimator. A Poisson tensor decomposition models the underlying distribution of the count data and guarantees non-negative estimated values and so can be safely used directly in entropy estimation. Our estimator is the first tensor-based estimator that exploits the underlying spatial Poisson process related to the histogram explicitly when estimating the probability density with low-rank tensor decompositions for the purpose of tensor completion. Furthermore, we demonstrate that our PTC estimator is a substantial improvement over standard histogram-based estimators for sub-Gaussian probability distributions because of the concentration of norm phenomenon.

42 ENGINEERING↗
DOE OSTI2025

Explainable and Differentiable Reinforcement Learning for Multi-objective Optimization in Particle Accelerators

Operating particle accelerators involves optimizing multiple goals simultaneously, which can be challenging due to trade-offs among objectives. While evolutionary algorithms like the genetic algorithm (GA) have been used for various Multi-Objective Optimization (MOO) tasks, they are not inherently suited for complex control problems. This talk highlights two variations of Reinforcement Learning (RL) for concurrently optimizing heat load and trip rates at the Continuous Electron Beam Accelerator Facility (CEBAF). The problem involves strict constraints on individual states, actions, and overall energy requirements of the beam. First, this talk highlights how differentiability can be harnessed through a Deep Differentiable Reinforcement Learning (DDRL) approach to address MOO issues within particle accelerators. We examine the DDRL method alongside Model Free Reinforcement Learning (MFRL), GA, and Bayesian Optimization (BO). The performance of these methods is assessed by generating a Pareto-front for two objectives. Our findings indicate that DDRL excels in handling high-dimensional problems more effectively than MFRL, BO, and GA. Next, we will show integration of explainable physics-based constraints into RL algorithms to enhance trans- parency and trust in decision-making processes by enabling users to verify that agents adhere to established physical principles. This surrogate function can be modeled using neural networks or sparse dictionary mod- els. By examining the mathematical form of the learned constraint function, we are able to confirm the agent has learned to use the established physics of each environment provided but the surrogate model. In addi- tion, we find that the introduction of a mathematical functional dictionary based surrogate model enables our reinforcement learning algorithms to reliably converge for difficult high-dimensional accelerator controls environments.

Rajput, Kishansingh [Thomas Jefferson National Acc↗
DOE OSTI2026

Progress toward double-differential cross-section measurements of single charged-pion production in charged current muon neutrino interactions on argon with SBND

Resonant neutrino–nucleus processes constitute a significant portion of neutrino interactions in the few-GeV energy region. Charged-current interactions with a muon and one single charged pion in the final state (CC1$\pi$) are primarily sensitive to resonant processes, while also receiving contributions from non-resonant processes and being strongly affected by nuclear effects and final-state interactions. A precise understanding of this channel is essential for improving neutrino interaction models and reducing systematic uncertainties in oscillation measurements. The Short-Baseline Near Detector (SBND) at Fermilab has collected the largest neutrino–argon dataset to date, providing an excellent opportunity to study CC1$\pi$ interactions on argon. This analysis aims to improve existing single-differential cross-section measurements by reporting results in a variety of kinematic variables and to perform the first double-differential CC1$\pi$ cross-section measurement on argon, using novel techniques for pion energy reconstruction.

Pelegrina-Gutiérrez, Luis [Granada U.] (ORCID:0000↗
DOE OSTI2023

Improved quantum algorithms for linear and nonlinear differential equations

We present substantially generalized and improved quantum algorithms over prior work for inhomogeneous linear and nonlinear ordinary differential equations (ODE). Specifically, we show how the norm of the matrix exponential characterizes the run time of quantum algorithms for linear ODEs opening the door to an application to a wider class of linear and nonlinear ODEs. In [1], a quantum algorithm for a certain class of linear ODEs is given, where the matrix involved needs to be diagonalizable. The quantum algorithm for linear ODEs presented here extends to many classes of non-diagonalizable matrices including singular matrices. The algorithm here is also exponentially faster than the bounds derived in [1] for certain classes of diagonalizable matrices. Our linear ODE algorithm is then applied to nonlinear differential equations using Carleman linearization (an approach taken recently by us in [2]). The improvement over that result is two-fold. First, we obtain an exponentially better dependence on error. This kind of logarithmic dependence on error has also been achieved by [3], but only for homogeneous nonlinear equations. Second, the present algorithm can handle any sparse matrix (that models dissipation) if it has a negative log-norm (including non-diagonalizable matrices), whereas [2] and [3] additionally require normality.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Secondary Wall Regulating NACs Differentially Bind at the Promoter at a CELLULOSE SYNTHASE A4 Cis -eQTL

Arabidopsis thaliana CELLULOSE SYNTHASE A4/7/8 (CESA4/7/8) are three non-redundant subunits of the secondary cell wall cellulose synthase complex. Transcript abundance of these genes can vary among genotypes and expression quantitative trait loci (eQTL) were identified in a recombinant population of the accessions Bay-o and Shahdara. Genetic mapping and analysis of the transcript levels of CESAs between two distinct near isogenic lines (NILs) confirmed a change in CESA4 expression that segregates within that interval. We sequenced the promoters and identified 16 polymorphisms differentiating CESA4 Sha and CESA4 Bay . In order to determine which of these SNPs could be responsible for this eQTL, we screened for transcription factor protein affinity with promoter fragments of CESA4 Bay , CESA4 Sha , and the reference genome CESA4 Col . The wall thickening activator proteins NAC SECONDARY WALL THICKENING PROMOTING FACTOR2 (NST2) and NST3 exhibited a decrease in binding with the CESA4 Sha promoter with a tracheary element-regulating cis-element (TERE) polymorphism. While NILs harboring the TERE polymorphisms exhibited significantly different CESA4 expression, cellulose crystallinity and cell wall thickness were indistinguishable. These results suggest that the TERE polymorphism resulted in differential transcription factor binding and CESA4 expression; yet A. thaliana is able to tolerate this transcriptional variability without compromising the structural elements of the plant, providing insight into the elasticity of gene regulation as it pertains to cell wall biosynthesis and regulation. We also explored available DNA affinity purification sequencing data to resolve a core binding site, C(G/T)TNNNNNNNA(A/C)G, for secondary wall NACs referred to as the VNS element.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Differential Regulation of Maize and Sorghum Orthologs in Response to the Fungal Pathogen Exserohilum turcicum

Pathogens that infect more than one host offer an opportunity to study how resistance mechanisms have evolved across different species. Exserohilum turcicum infects both maize and sorghum and the isolates are host-specific, offering a unique system to examine both compatible and incompatible interactions. We conducted transcriptional analysis of maize and sorghum in response to maize-specific and sorghum-specific E. turcicum isolates and identified functionally related co-expressed modules. Maize had a more robust transcriptional response than sorghum. E. turcicum responsive genes were enriched in core orthologs in both crops, but only up to 16% of core orthologs showed conserved expression patterns. Most changes in gene expression for the core orthologs, including hub genes, were lineage specific, suggesting a role for regulatory divergent evolution. We identified several defense-related shared differentially expressed (DE) orthologs with conserved expression patterns between the two crops, suggesting a role for parallel evolution of those genes in both crops. Many of the differentially expressed genes (DEGs) during the incompatible interaction were related to quantitative disease resistance (QDR). This work offers insights into how different hosts with relatively recent divergence interact with a common pathogen. Our results are important for developing resistance to this critical pathogen and understanding the evolution of host-pathogen interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Perennial grass root system specializes for multiple resource acquisitions with differential elongation and branching patterns

Roots optimize the acquisition of limited soil resources, but relationships between root forms and functions have often been assumed rather than demonstrated. Furthermore, how root systems co-specialize for multiple resource acquisitions is unclear. Theory suggests that trade-offs exist for the acquisition of different resource types, such as water and certain nutrients. Measurements used to describe the acquisition of different resources should then account for differential root responses within a single system. To demonstrate this, we grew Panicum virgatum in split-root systems that vertically partitioned high water availability from nutrient availability so that root systems must absorb the resources separately to fully meet plant demands. We evaluated root elongation, surface area, and branching, and we characterized traits using an order-based classification scheme. Plants allocated approximately 3/4th of primary root length towards water acquisition, whereas lateral branches were progressively allocated towards nutrients. However, root elongation rates, specific root length, and mass fraction were similar. Our results support the existence of differential root functioning within perennial grasses. Similar responses have been recorded in many plant functional types suggesting a fundamental relationship. Root responses to resource availability can be incorporated into root growth models via maximum root length and branching interval parameters.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Deletion of the EP402R Gene from the Genome of African Swine Fever Vaccine Strain ASFV-G-ΔI177L Provides the Potential Capability of Differentiating between Infected and Vaccinated Animals

The African swine fever virus (ASFV) mutant ASFV-G-ΔI177L is a safe and efficacious vaccine which induces protection against the challenge of its parental virus, the Georgia 2010 isolate. Although a genetic DIVA (differentiation between infected and vaccinated animals) assay has been developed for this vaccine, still there is not a serological DIVA test for differentiating between animals vaccinated with ASFV-G-ΔI177L and those infected with wild-type viruses. In this report, we describe the development of the ASFV-G-ΔI177L mutant having deleted the EP402R gene, which encodes for the viral protein responsible for mediating the hemadsorption of swine erythrocytes. The resulting virus, ASFV-G-ΔI177L/ΔEP402R, does not have a decreased ability to replicates in swine macrophages when compared with the parental ASFV-G-ΔI177L. Domestic pigs intramuscularly (IM) inoculated with either 10 2 or 10 6 HAD 50 of ASFV-G-ΔI177L/ΔEP402R remained clinically normal, when compared with a group of mock-vaccinated animals, indicating the absence of residual virulence. Interestingly, an infectious virus could not be detected in the blood samples of the ASFV-G-ΔI177L/ΔEP402R-inoculated animals in either group at any of the time points tested. Furthermore, while all of the mock-inoculated animals presented a quick and lethal clinical form of ASF after the intramuscular inoculation challenge with 10 2 HAD 50 of highly virulent parental field isolate Georgia 2010 (ASFV-G), all of the ASFV-G-ΔI177L/ΔEP402R-inoculated animals were protected, remaining clinically normal until the end of the observational period. Most of the ASFV-G-ΔI177L/ΔEP402R-inoculated pigs developed strong virus-specific antibody responses against viral antigens, reaching maximum levels at 28 days post inoculation. Importantly, all of the sera collected at that time point in the ASFV-G-ΔI177L/ΔEP402R-inoculated pigs did not react in a direct ELISA coated with the recombinant EP402R protein. Conversely, the EP402R protein was readily recognized by the pool of sera from the animals immunized with recombinant live attenuated vaccine candidates ASFV-G-ΔI177L, ASFV-G-ΔMGF, or ASFV-G-Δ9GL/ΔUK. Therefore, ASFV-G-ΔI177L/ΔEP402R is a novel, safe and efficacious candidate with potential to be used as an antigenically DIVA vaccine.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

SPIRAL STRUCTURE AND DIFFERENTIAL DUST SIZE DISTRIBUTION IN THE LkHα 330 DISK

Dust trapping accelerates the coagulation of dust particles, and, thus, it represents an initial step toward the formation of planetesimals. We report H-band (1.6 μm) linear polarimetric observations and 0.87 mm interferometric continuum observations toward a transitional disk around LkHα 330. As a result, a pair of spiral arms were detected in the H-band emission, and an asymmetric (potentially arm-like) structure was detected in the 0.87 mm continuum emission. We discuss the origin of the spiral arm and the asymmetric structure and suggest that a massive unseen planet is the most plausible explanation. The possibility of dust trapping and grain growth causing the asymmetric structure was also investigated through the opacity index (β) by plotting the observed spectral energy distribution slope between 0.87 mm from our Submillimeter Array observation and 1.3 mm from literature. The results imply that grains are indistinguishable from interstellar medium-like dust in the east side (β=2.0±0.5) but are much smaller in the west side β=0.7 -0.4 +0.5 , indicating differential dust size distribution between the two sides of the disk. Combining the results of near-infrared and submillimeter observations, we conjecture that the spiral arms exist at the upper surface and an asymmetric structure resides in the disk interior. Future observations at centimeter wavelengths and differential polarization imaging in other bands (Y–K) with extreme AO imagers are required to understand how large dust grains form and to further explore the dust distribution in the disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

DOES A DIFFERENTIATED, CARBONATE-RICH, ROCKY OBJECT POLLUTE THE WHITE DWARF SDSS J104341.53+085558.2?

We present spectroscopic observations of the dust- and gas-enshrouded, polluted, single white dwarf star SDSS J104341.53+085558.2 (hereafter SDSS J1043+0855). Hubble Space Telescope Cosmic Origins Spectrograph far-ultraviolet spectra combined with deep Keck HIRES optical spectroscopy reveal the elements C, O, Mg, Al, Si, P, S, Ca, Fe, and Ni and enable useful limits for Sc, Ti, V, Cr, and Mn in the photosphere of SDSS J1043+0855. From this suite of elements we determine that the parent body being accreted by SDSS J1043+0855 is similar to the silicate Moon or the outer layers of Earth in that it is rocky and iron-poor. Combining this with comparison to other heavily polluted white dwarf stars, we are able to identify the material being accreted by SDSS J1043+0855 as likely to have come from the outermost layers of a differentiated object. Furthermore, we present evidence that some polluted white dwarfs (including SDSS J1043+0855) allow us to examine the structure of differentiated extrasolar rocky bodies. Enhanced levels of carbon in the body polluting SDSS J1043+0855 relative to the Earth–Moon system can be explained with a model where a significant amount of the accreted rocky minerals took the form of carbonates; specifically, through this model the accreted material could be up to 9% calcium-carbonate by mass.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

On Differential Rotation and Overshooting in Solar-like Stars

We seek to characterize how the change of global rotation rate influences the overall dynamics and large-scale flows arising in the convective envelopes of stars covering stellar spectral types from early G to late K. We do so through numerical simulations with the ASH code, where we consider stellar convective envelopes coupled to a radiative interior with various global properties. As solar-like stars spin down over the course of their main sequence evolution, such a change must have a direct impact on their dynamics and rotation state. We indeed find that three main states of rotation may exist for a given star: anti-solar-like (fast poles, slow equator), solar-like (fast equator, slow poles), or a cylindrical rotation profile. Under increasingly strict rotational constraints, the last profile can further evolve into a Jupiter-like profile, with alternating prograde and retrograde zonal jets. We have further assessed how far the convection and meridional flows overshoot into the radiative zone and investigated the morphology of the established tachocline. Using simple mixing length arguments, we are able to construct a scaling of the fluid Rossby number R{sub of}= ω-tilde /2Ω{sub ∗}∼ v-tilde /2Ω{sub ∗}R{sub ∗}, which we calibrate based on our 3D ASH simulations. We can use this scaling to map the behavior of differential rotation versus the global parameters of stellar mass and rotation rate. Finally, we isolate a region on this map (R {sub of} ≳ 1.5–2) where we posit that stars with an anti-solar differential rotation may exist in order to encourage observers to hunt for such targets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Predicting the Slowing of Stellar Differential Rotation by Instability-driven Turbulence

Abstract Differentially rotating stars and planets transport angular momentum (AM) internally due to turbulence at rates that have long been a challenge to predict reliably. We develop a self-consistent saturation theory, using a statistical closure approximation, for hydrodynamic turbulence driven by the axisymmetric Goldreich–Schubert–Fricke instability at the stellar equator with radial differential rotation. This instability arises when fast thermal diffusion eliminates the stabilizing effects of buoyancy forces in a system where a stabilizing entropy gradient dominates over the destabilizing AM gradient. Our turbulence closure invokes a dominant three-wave coupling between pairs of linearly unstable eigenmodes and a near-zero frequency, viscously damped eigenmode that features latitudinal jets. We derive turbulent transport rates of momentum and heat and provide them in analytic forms. Such formulae, free of tunable model parameters, are tested against direct numerical simulations; the comparison shows good agreement. They improve upon prior quasi-linear or “parasitic saturation” models containing a free parameter. Given model correspondences, we also extend this theory to heat and compositional transport for axisymmetric thermohaline-instability-driven turbulence in certain regimes.

Astronomy & Astrophysics↗
DOE OSTI2025

Joint Modeling of Quasar Variability and Accretion Disk Reprocessing Using Latent Stochastic Differential Equations

Quasars are bright active galactic nuclei powered by the accretion of matter around supermassive black holes at the center of galaxies. Their stochastic brightness variability depends on the physical properties of the accretion disk and black hole. The upcoming Rubin Observatory Legacy Survey of Space and Time (LSST) is expected to observe tens of millions of quasars, so there is a need for efficient techniques like machine learning that can handle the large volume of data. Quasar variability is believed to be driven by an X-ray corona, which is reprocessed by the accretion disk and emitted as UV/optical variability. We are the first to introduce an auto-differentiable simulation of the accretion disk and reprocessing. We use the simulation as a direct component of our neural network to jointly model the driving variability and reprocessing, trained with supervised learning on simulated LSST-like 10 yr quasar light curves. We encode the light curves using a transformer encoder, and the driving variability is reconstructed using latent stochastic differential equations, a physically motivated generative deep learning method that can model continuous-time stochastic dynamics. By embedding the physical processes of the driving signal and reprocessing into our network, we achieve a model that is more robust and interpretable. We demonstrate that our model outperforms a Gaussian process regression baseline and can infer accretion disk parameters and time delays between wave bands, even for out-of-distribution driving signals. Our approach provides a powerful framework that can be adapted to solve other inverse problems in multivariate time series.

Fagin, Joshua [City Univ. of New York (CUNY), NY (↗
DOE OSTI2021

Data from: Niche differentiation of bacterial versus archaeal soil nitrifiers induced by ammonium inhibition along a management gradient

Soil nitrification, mediated mainly by ammonia oxidizing archaea (AOA) and bacteria (AOB), converts ammonium (NH4+) to nitrite (NO2−) and thence nitrate (NO3−). To better understand ecological differences between AOA and AOB, we investigated the nitrification kinetics of AOA and AOB under eight replicated cropped and unmanaged ecosystems (including two fertilized natural systems) along a long-term management intensity gradient in the upper U.S. Midwest. For five of eight ecosystems, AOB but not AOA exhibited Haldane kinetics (inhibited by high NH4+ additions), especially in perennial and successional systems. In contrast, AOA predominantly exhibited Michaelis-Menten kinetics, suggesting greater resistance to high nitrogen inputs than AOB. These responses suggest the potential for NH4+-induced niche differentiation between AOA and AOB. Additionally, long-term fertilization significantly enhanced maximum nitrification rates (Vmax) in the early successional systems for both AOA and AOB, but not in the deciduous forest systems. This was likely due to pH suppression of nitrification in the acidic forest soils, corroborated by a positive correlation of Vmax with soil pH but not with amoA gene abundance. Results also demonstrated that soil nitrification potentials were relatively stable, as there were no seasonal differences. Overall, results suggest that (1) NH4+ inhibition of AOB but not AOA could be another factor contributing to niche differentiation between AOA and AOB in soil, and (2) nitrification by both AOA and AOB can be significantly promoted by long-term nitrogen inputs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Efficacy of a portable, moderate-resolution, fast-scanning differential mobility analyzer for ambient aerosol size distribution measurements

Abstract. Ambient aerosol size distributions obtained with a compact scanning mobility analyzer, the “Spider” differential mobility analyzer (DMA), are compared to those obtained with a conventional mobility analyzer, with specific attention to the effect of mobility resolution on the measured size distribution parameters. The Spider is a 12 cm diameter radial differential mobility analyzer that spans the 10–500 nm size range with 30 s mobility scans. It achieves its compact size by operating at a nominal mobility resolution R=3 (sheath flow = 0.9 L min−1; aerosol flow = 0.3 L min−1) in place of the higher ratio of sheath flow to aerosol flow commonly used. The question addressed here is whether the lower resolution is sufficient to capture key characteristics of ambient aerosol size distributions. The Spider, operated at R=3 with 30 s up- and downscans, was co-located with a TSI 3081 long-column mobility analyzer, operated at R=10 with a 360 s sampling duty cycle. Ambient aerosol data were collected over 26 consecutive days of continuous operation, in Pasadena, CA. Over the 17–500 nm size range, the two instruments exhibit excellent correlation in the total particle number concentrations and geometric mean diameters, with regression slopes of 1.13 and 1.00, respectively. Our results suggest that particle sizing at a lower resolution than typically employed may be sufficient to obtain key properties of ambient size distributions, at least for these two moments of the size distribution. Moreover, it enables better counting statistics, as the wider transfer function for a given aerosol flow rate results in a higher counting rate.

54 ENVIRONMENTAL SCIENCES↗