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 1,513 records · Page 84

Continuously differentiable PIC shape functions for triangular meshes

In this study, a new class of continuously-differentiable shape functions is developed and applied to two-dimensional electrostatic PIC simulation on an unstructured simplex (triangle) mesh. It is shown that troublesome aliasing instabilities are avoided for cold plasma simulation in which the Debye length is as small as 0.01 cell sizes. These new shape functions satisfy all requirements for PIC particle shape. They are non-negative, have compact support, and partition unity. They are given explicitly by cubic expressions in the usual triangle logical (areal) coordinates. The shape functions are not finite elements because their structure depends on the topology of the mesh, in particular, the number of triangles neighboring each mesh vertex. Nevertheless, they may be useful as approximations to solution of other problems in which continuity of derivatives is required or desired.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An improved iterative HDG approach for partial differential equations

Here we propose and analyze an iterative high-order hybridized discontinuous Galerkin (iHDG) discretization for linear partial differential equations. We improve our previous work [45] in several directions: 1) the improved iHDG approach converges in a finite number of iterations for the scalar transport equation; 2) it is unconditionally convergent for both the linearized shallow water system and the convection–diffusion equation; 3) it has improved stability and convergence rates; 4) we uncover a relationship between the number of iterations and time stepsize, solution order, meshsize and the equation parameters. This allows us to choose the time stepsize such that the number of iterations is approximately independent of the solution order and the meshsize; and 5) we provide both strong and weak scalings of the improved iHDG approach up to 16,384 cores. A connection between iHDG and time integration methods such as parareal and implicit/explicit methods are discussed. Extensive numerical results for linear (and nonlinear) PDEs are presented to verify the theoretical findings.

97 MATHEMATICS AND COMPUTING↗

Quantum mechanical closure of partial differential equations with symmetries

We develop a statistical framework for the dynamical closure of spatiotemporal dynamics governed by partial differential equations. Employing the mathematical framework of quantum mechanics to embed the original classical dynamics into a quantum mechanical representation, we use the space of quantum density operators to model the unresolved degrees of freedom of the original dynamics in a statistical sense, and the framework of quantum measurement to predict their contributions to the resolved dynamics. The embedded dynamics is discretized by a positivity preserving process, leading to a compressed representation that is invariant under the dynamical symmetries of the resolved dynamics. We present a data based formulation of the closure scheme and apply it to a closure problem for the shallow water equations. The numerical results demonstrate that our closure model can accurately predict the main features of the true dynamics, including for out of sample initial conditions.

Delay embedding↗

Evaluation of Respiratory Triggered Diffusion-Weighted MRI with Three b-Values Compared to ADC Map and Fast Spin Echo Heavily T2W in Differential Diagnosis of Hemangioma from Malignant Liver Lesions

The goal of this study was to evaluate the ability of diffusion-weighted imaging (DWI) with three different b-values compared to apparent diffusion coefficient (ADC) map and fast spin echo heavily T{sub 2} weighted (FSE-HT{sub 2}W) in differential diagnosis of hemangioma from malignant liver lesions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Load-Displacement Relation and Gap Distribution Between Rough Surfaces: Partial Differential Equations Approach

We develop a theoretical model to predict the load-displacement relation and probability density function for gaps between contacting rough surfaces. We derive partial differential equations from the previous model (Joe et al., 2018), and extend them to the non-adhesive contact problem. The predictions of the present theory are compared with numerical results using the Green's function molecular dynamics algorithm, and a good agreement is obtained.

adhesion↗

Characterization of the microstructures and the shape memory properties of the Fe-Mn-Si-Cr-Ni-C shape memory alloy after severe plastic deformation by differential speed rolling and subsequent annealing

Highlights: • The Fe-SMA was subjected to severe plastic deformation and then annealed. • Recovery strain increased with increasing grain size. • An optimum grain size existed for the largest recovery stress. • Grain reduction is useful for having a high recovery stress. - Abstract: The Fe-13.51Mn-4.82Si-8.32Cr-3.49Ni-0.15C shape memory alloy was subjected to severe plastic deformation by high-ratio differential speed rolling and subsequent annealing (between 873 and 1373 K) to find the factors that affect recovery stress and recovery strain. As the annealing temperature increased (i. e., as grain size increased), the recovery strain increased. This was attributed to decrease in the density of grain boundaries and the density of twins with increasing grain size. The recovery stress, however, showed a different behavior. It increased with annealing temperature up to 973 K, at which an ultrafine grain size near 1.6 μm was obtained, and then decreased with further increase in the temperature. The difference between the recovery strain and the plastic strain induced by recovery stress under constrained displacement during cooling was proposed to determine the amount of recovery stress. There was an optimum grain size for the largest difference between the recovery strain and the plastic strain. The current study shows that grain-size reduction is not good for having a high recovery strain but it is beneficial for having a high recovery stress in the Fe-Mn-Si shape memory alloys.

36 MATERIALS SCIENCE↗

Shape memory and superelasticity of nanograined Ti-51.2 at.% Ni alloy processed by severe plastic deformation via high-ratio differential speed rolling

A novel method of producing nanograined Ni-rich superelastic NiTi alloys in sheet form was proposed using a combination of severe plastic deformation via high-ratio differential speed rolling (HRDSR) and post-deformation annealing. The HRDSR-processed microstructure was composed of heavily deformed austenite and martensite grains, and amorphous phases. After annealing at 673 K, the severely deformed microstructure with no functional properties evolved to the nanograined structure (20–70 nm) composed of austenite and martensite nanograins and sub-nanograins through static recovery or continuous static recrystallization process. The nanograined microstructure had a high resistance to martensitic transformation upon cooling and slip deformation during straining. As a result, the HRDSR-processed alloy annealed at 673K exhibited superior superelasticity compared to the alloys with coarse grains. At the higher annealing temperature of 873 K, the micron-sized recrystallized grains with low dislocations developed through discontinuous static recrystallization process. In this case, deformation during straining was governed by the detwinning of twinned martensite, and as a result, shape memory effect was more significantly pronounced than superelasticity.

36 MATERIALS SCIENCE↗

Study of the dislocation activity in a Mg–Y alloy by differential aperture X-ray microscopy

Differential aperture X-ray microscopy (DAXM) was employed in this study to reveal the development of geometrically necessary dislocations (GNDs) when a Mg-5 wt% Y alloy is subjected to bending. The characters of the GNDs were determined from the streak direction of the diffraction peaks. Five grains were analyzed voxel by voxel for the GND content. Among those voxels with peak streaking, 77% contain basal or non-basal < a > GNDs, while the other 23% contain < c + a > GNDs. Furthermore, those < c + a > dislocations were nucleated to accommodate local stress concentration. Even though < c + a > dislocations are less mobile, they make important contribution to the high ductility of Mg-Y alloys.

36 MATERIALS SCIENCE↗

Differentiating supported platinum single atoms, clusters and nanoparticles by styrene hydrogenation

Supported metal catalysts often consist of metal sites ranging from nanoparticles to subnanometer clusters and single atoms. It remains a necessity to differentiate these sites to guide design of optimal catalysts. Here we report a simple method to assess the distribution of metal active sites in catalyst samples. The method takes the advantage of the structure sensitivity of styrene hydrogenation over titania supported platinum (Pt) catalysts with Pt aggregates varied from single atom to ~1.40 nm nanoparticles. The physicochemical properties were characterized by STEM, XPS, XANES, H2-TPR and CO-chemisorption measurements. The reactivity of Pt sites was quantified by styrene hydrogenation at ambient conditions. The nanometer-sized Pt clusters have significantly higher activity than Pt nanoparticles, sub-nanometer clusters or isolated single atoms. Additionally, the relationship between activity and structural/electronic properties of Pt sites influenced by particle sizes was discussed. Similar relationship was found in the carbon supported Pt catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring the interplay of resilience and energy consumption for a task-based partial differential equations preconditioner

Here, we discuss algorithm-based resilience to silent data corruptions (SDCs) in a task-based domain-decomposition preconditioner for partial differential equations (PDEs). The algorithm exploits a reformulation of the PDE as a sampling problem, followed by a solution update through data manipulation that is resilient to SDCs. The implementation is based on a server-client model where all state information is held by the servers, while clients are designed solely as computational units. Scalability tests run up to ~51K cores show a parallel efficiency greater than 90%. We use a 2D elliptic PDE and a fault model based on random single and double bit-flip to demonstrate the resilience of the application to synthetically injected SDC. We discuss two fault scenarios: one based on the corruption of all data of a target task, and the other involving the corruption of a single data point. We show that for our application, given the test problem considered, a four-fold increase in the number of faults only yields a 2% change in the overhead to overcome their presence, from 7% to 9%. We then discuss potential savings in energy consumption via dynamic voltage/frequency scaling, and its interplay with fault-rates, and application overhead.

97 MATHEMATICS AND COMPUTING↗

Observation of a multiplicity dependence in the $\mathcal{p}_T$-differential charm baryon-to-meson ratios in proton–proton collisions at $\sqrt{s}$ = 13 TeV

The production of prompt D 0 , D$^+_s$, and Λ$^+_c$ hadrons, and their ratios, D$^+_s$/D 0 and Λ$^+_c$/D 0 , are measured in proton–proton collisions at $\sqrt{s}$ = 13 TeV at midrapidity (|y| < 0.5) with the ALICE detector at the LHC. The measurements are performed as a function of the charm-hadron transverse momentum (p T ) in intervals of charged-particle multiplicity, measured with two multiplicity estimators covering different pseudorapidity regions. While the strange to non-strange D$^+_s$/D 0 ratio indicates no significant multiplicity dependence, the baryon-to-meson p T -differential Λ$^+_c$/D 0 ratio shows a multiplicity-dependent enhancement, with a significance of 5.3σ for 1 < p T < 12 GeV/c, comparing the highest multiplicity interval with respect to the lowest one. The measurements are compared with a theoretical model that explains the multiplicity dependence by a canonical treatment of quantum charges in the statistical hadronisation approach, and with predictions from event generators that implement colour reconnection mechanisms beyond the leading colour approximation to model the hadronisation process. The Λ$^+_c$/D 0 ratios as a function of p T present a similar shape and magnitude as the Λ/K$^0_s$ ratios in comparable multiplicity intervals, suggesting a potential common mechanism for light- and charm-hadron formation, with analogous multiplicity dependence. The p T -integrated ratios, extrapolated down to p T = 0, do not show a significant dependence on multiplicity within the uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurement of inclusive and differential cross sections for W + W − production in proton-proton collisions at $\sqrt{s} = 13$ TeV

Measurements at $\sqrt{s} = 13$ TeV of the opposite-sign W boson pair production cross section in proton-proton collisions are presented. The data used in this study were collected with the CMS detector at the CERN LHC in 2022, and correspond to an integrated luminosity of 34.8 fb -1 . Events are selected by requiring one electron and one muon of opposite charge. A maximum likelihood fit is performed on signal- and background-enriched data categories defined by the flavor and charge of the leptons, the number of jets, and number of jets originating from b quarks. The overall sensitivity is significantly better than that of previous results with a similar integrated luminosity. The improvement comes from a more refined control of experimental uncertainties and an improved fit strategy. An inclusive W + W - production cross section of 125.7 ± 5.6 pb is measured, in agreement with standard model predictions. Cross sections are also reported in a fiducial region close to that of the detector acceptance, both inclusively and differentially, as a function of the jet multiplicity in the event. For the first time in proton-proton collisions, WW events with zero, one, and at least two jets are studied simultaneously and compared with recent theoretical predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

First measurement of differential cross sections for muon neutrino charged current interactions on argon with a two-proton final state using the MicroBooNE detector

We present the first measurement of differential cross sections for charged-current muon neutrino interactions on argon with one muon, two protons, and no pions in the final state. These final states are dominated by two-nucleon knockout interactions, which are complicated to model and for which there is currently limited information about the characteristics of these interactions in existing neutrino-nucleus scattering data. Detailed investigations of two-nucleon knockout are vital to support upcoming experiments exploring the nature of the neutrino. Among the different kinematic quantities measured, the opening angle between the two protons, the angle between the total proton momentum and the muon, and the total transverse momentum of the final state system are most sensitive to the underlying physics processes as embodied in various theoretical models.

2-particle 2-hole states↗

Toward an event-level analysis of hadron structure using differential programming

Reconstructing the internal properties of hadrons in terms of fundamental quark and gluon de- grees of freedom is a central goal in nuclear and particle physics. This effort lies at the core of major experimental programs, such as the Jefferson Lab 12 GeV program and the upcoming Electron-Ion Collider. A primary challenge is the inherent inverse problem: converting large-scale observational data from collision events into the fundamental QCD-defined densities that characterize the micro- scopic structure of hadronic systems. Recent advances in AI and machine learning have opened new avenues for addressing this challenge using deep learning techniques. A particularly promising direction is the integration of complex theoretical calculations and experimental simulations into a unified framework capable of reconstructing these densities directly from event-level information. In this document, we introduce a key algorithm called LOITS, which enables differentiable program- ming within such a framework, facilitating the use of AI/ML techniques to solve the inverse problem of QCF reconstruction at the event level.

Braga, Kevin [College of William and Mary, William↗

Reproducing curvature effects due to differential diffusion in tabulated chemistry for premixed flames

Tabulated chemistry reduces the cost of modeling reacting flows by transporting fewer scalars. It also has the benefit of introducing fewer terms that must be closed for Reynolds-averaged Navier-Stokes and large eddy simulation modeling. In this study of premixed lean hydrogen–air flames, an extension to an existing tabulated chemistry model was developed by including differential diffusion and thermal diffusion of multiple species, using a mixture-averaged diffusion model. The transport equations for the tabulation variables (a progress variable and a mixture fraction-like variable) have been formulated, and curvature effects are intrinsic to the transport equations derived from first principles. The mathematical derivation is general, such that other fuel/air mixtures could be examined. The new tabulated chemistry model was then validated using a range of configurations, including flames with curvature, extinction, and hot spots. The model validation indicated excellent agreement with detailed chemistry simulations for lean premixed hydrogen–air mixtures. Conditional means of the progress variable source term were shown to agree well with detailed chemistry simulations of two-dimensional freely propagating flames and three-dimensional turbulent flames. Finally, progress variable source terms three times higher than one-dimensional flat flame results were identified and validated with detailed chemistry, and super-adiabatic hot spots were also predicted in the turbulent flame configuration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Distributed Automatic Differentiation for Ptychography

Synchrotron radiation light source facilities are leading the way to ultrahigh resolution X-ray imaging. High resolution imaging is essential to understanding the fundamental structure and interaction of materials at the smallest length scale possible. Diffraction based methods achieve nanoscale imaging by replacing traditional objective lenses by pixelated area detectors and computational image reconstruction. Among these methods, ptychography is quickly becoming the standard for sub-30 nanometer imaging of extended samples, but at the expense of increasingly high data rates and volumes. This paper presents a new distributed algorithm for solving the ptychographic image reconstruction problem based on automatic differentiation. Input datasets are subdivided between multiple graphics processing units (GPUs); each subset of the problem is then solved either entirely independent of other subsets (asynchronously) or through sharing gradient information with other GPUs (synchronously). The algorithm was evaluated on simulated and real data acquired at the Advanced Photon Source, scaling up to 192 GPUs. The synchronous variant of our method outperformed an existing multi-GPU implementation in terms of accuracy while running at a comparable execution time.

distributed algorithms↗

Oil price states and drivers: An analysis of the second-month spot-futures price differential

Oil remains a dominant component of global energy use, and its price, characterized by frequent changes and an ever-present potential for large swings, continues to be a focus of industry participants, policymakers and analysts attention. Here, this study examines the behavior of future spot oil prices using a continuous-time hidden Markov model (HMM) and daily price data spanning years 2007 to 2024. We identify six states in the second-month WTI spot-futures price differential and assess the roles of eleven futures price, quantity, financial market, and geopolitical/volatility variables in each state. The model yields several insights into the workings of the oil market and the relative roles of these drivers. We find support for several theoretical and empirical findings in the oil market literature, including the role of inventory, volatility/risk, and market responses to contango/backwardation in futures markets. A novel finding is that “normal contango” conditions represent a significant portion of next-day states in our in-sample data. Under the most volatile normal contango state, many of the oil market drivers differ markedly in coefficient signs and magnitudes from those in other states. The resulting model also performed well out-of-sample and would, in addition to understanding the impact of market drivers, be useful for short-term forecasting. Overall, the findings highlight the highly non-linear, regime-dependent interactions of the oil price and its drivers, emphasizing the importance of detailed information to market stakeholders.

Oladosu, Gbadebo A. [Oak Ridge National Laboratory↗

Thermodynamic analysis of a MEMS based differential scanning calorimeter model

Thin film based differential scanning calorimeter (DSC) was developed and utilized for the study of phase-transition kinetics in protein samples. In this paper, we proposed a MEMS based DSC model. The distance between the sample and reference cell was optimized based on the thermal analysis result. The influence of the response time to the shape of the DSC thermograms and the transition temperature was also studied systematically. Meanwhile, Lumry-Eyring model was used to investigate the relationship between the scanning rate and the transition temperature. The modeling result was verified by the lysozyme denaturation tests. This study provided a pathway for the optimization of MEMS based DSC array and the correction of the DSC thermograms.

47 OTHER INSTRUMENTATION↗