Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential mode”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

CHEMNODE: CHEMICAL KINETICS SOLVER APPROACH BASED ON NEURAL ORDINARY DIFFERENTIAL EQUATIONS

This software contains an algorithm to model and predict the time evolution of chemical kinetics in reacting flow simulations using neural ordinary differential equations (NODEs). The approach works by using artificial neural networks to predict the chemical source terms. It achieves this by using forward-mode automatic differentiation and the Levenberg-Marquardt algorithm to adjust the neural network parameters, such that the discrepancies between the actual and predicted species profiles are minimized.

OWOYELE, OPEOLUWA↗
NASA NTRS1999

Demonstration of Automatically-Generated Adjoint Code for Use in Aerodynamic Shape Optimization

Gradient-based optimization requires accurate derivatives of the objective function and constraints. These gradients may have previously been obtained by manual differentiation of analysis codes, symbolic manipulators, finite-difference approximations, or existing automatic differentiation (AD) tools such as ADIFOR (Automatic Differentiation in FORTRAN). Each of these methods has certain deficiencies, particularly when applied to complex, coupled analyses with many design variables. Recently, a new AD tool called ADJIFOR (Automatic Adjoint Generation in FORTRAN), based upon ADIFOR, was developed and demonstrated. Whereas ADIFOR implements forward-mode (direct) differentiation throughout an analysis program to obtain exact derivatives via the chain rule of calculus, ADJIFOR implements the reverse-mode counterpart of the chain rule to obtain exact adjoint form derivatives from FORTRAN code. Automatically-generated adjoint versions of the widely-used CFL3D computational fluid dynamics (CFD) code and an algebraic wing grid generation code were obtained with just a few hours processing time using the ADJIFOR tool. The codes were verified for accuracy and were shown to compute the exact gradient of the wing lift-to-drag ratio, with respect to any number of shape parameters, in about the time required for 7 to 20 function evaluations. The codes have now been executed on various computers with typical memory and disk space for problems with up to 129 x 65 x 33 grid points, and for hundreds to thousands of independent variables. These adjoint codes are now used in a gradient-based aerodynamic shape optimization problem for a swept, tapered wing. For each design iteration, the optimization package constructs an approximate, linear optimization problem, based upon the current objective function, constraints, and gradient values. The optimizer subroutines are called within a design loop employing the approximate linear problem until an optimum shape is found, the design loop limit is reached, or no further design improvement is possible due to active design variable bounds and/or constraints. The resulting shape parameters are then used by the grid generation code to define a new wing surface and computational grid. The lift-to-drag ratio and its gradient are computed for the new design by the automatically-generated adjoint codes. Several optimization iterations may be required to find an optimum wing shape. Results from two sample cases will be discussed. The reader should note that this work primarily represents a demonstration of use of automatically- generated adjoint code within an aerodynamic shape optimization. As such, little significance is placed upon the actual optimization results, relative to the method for obtaining the results.

Green, Lawrence↗
DOE OSTI2024

Hutchinson Trace Estimation for high-dimensional and high-order Physics-Informed Neural Networks

Physics-Informed Neural Networks (PINNs) have proven effective in solving partial differential equations (PDEs), especially when some data are available by seamlessly blending data and physics. However, extending PINNs to high-dimensional and even high-order PDEs encounters significant challenges due to the computational cost associated with automatic differentiation in the residual loss function calculation. Herein, we address the limitations of PINNs in handling high-dimensional and high-order PDEs by introducing the Hutchinson Trace Estimation (HTE) method. Starting with the second-order high-dimensional PDEs, which are ubiquitous in scientific computing, HTE is applied to transform the calculation of the entire Hessian matrix into a Hessian vector product (HVP). This approach not only alleviates the computational bottleneck via Taylor-mode automatic differentiation but also significantly reduces memory consumption from the Hessian matrix to an HVP’s scalar output. We further showcase HTE’s convergence to the original PINN loss and its unbiased behavior under specific conditions. Comparisons with the Stochastic Dimension Gradient Descent (SDGD) highlight the distinct advantages of HTE, particularly in scenarios with significant variability and variance among dimensions. We further extend the application of HTE to higher-order and higher-dimensional PDEs, specifically addressing the biharmonic equation. By employing tensor-vector products (TVP), HTE efficiently computes the colossal tensor associated with the fourth-order high-dimensional biharmonic equation, saving memory and enabling rapid computation. The effectiveness of HTE is illustrated through experimental setups, demonstrating comparable convergence rates with SDGD under memory and speed constraints. Additionally, HTE proves valuable in accelerating the Gradient-Enhanced PINN (gPINN) version as well as the Biharmonic equation. Overall, HTE opens up a new capability in scientific machine learning for tackling high-order and high-dimensional PDEs.

Curse of dimensionality↗
DOE OSTI2025

Quantitative phase retrieval and characterization of magnetic nanostructures via Lorentz (scanning) transmission electron microscopy

Magnetic materials phase reconstruction using Lorentz transmission electron microscopy (LTEM) measurements have traditionally been achieved using longstanding methods such as off-axis holography (OAH) fast-Fourier transform technique and the transport-of-intensity equation (TIE). The increase in access to processing power alongside the development of advanced algorithms have allowed for phase retrieval of nanoscale magnetic materials with greater efficacy and resolution. Specifically, reverse-mode automatic differentiation (RMAD) and the extended electron ptychography iterative engine (ePIE) are two recent developments of phase retrieval that can be applied to analyzing micro-to-nano- scale magnetic materials. This work evaluates phase retrieval using TIE, RMAD, and ePIE in simulations of Permalloy (Ni 80 Fe 20 ) nanoscale islands, or nanomagnets. Extending beyond simulations, we demonstrate total phase retrieval and image reconstructions of a NiFe nanowire using OAH and RMAD in LTEM and ePIE in Lorentz-mode-4D scanning transmission electron microscopy experiments and determine the saturation magnetization through corroborations with micromagnetic modeling. Finally, we demonstrate the efficacy of these methods in retrieving the total phase and highlight its use in characterizing and analyzing the proximity effect of the magnetic nanostructures.

Lorentz transmission electron microscopy↗
NASA NTRS1994

Mass analysis addition to the Differential Ion Flux Probe (DIFP) study

The objective of this study is to develop a technique to measure the characteristics of space plasmas under highly disturbed conditions; e.g., non-Maxwellian plasmas with strong drifting populations and plasmas contaminated by spacecraft outgassing. The approach, conducted in conjunction with current MSFC activities, is to extend the capabilities of the Differential Ion Flux Probe (DIFP) to include a high throughput mass measurement that does not require either high voltage or contamination sensitive devices such as channeltron electron multipliers or microchannel plates. This will significantly reduce the complexity and expense of instrument fabrication, testing, and integration of flight hardware compared to classical mass analyzers. The feasibility of the enhanced DIFP has been verified by using breadboard test models in a controlled plasma environment. The ability to manipulate particles through the instrument regardless of incident angle, energy, or ionic component has been amply demonstrated. The energy analysis mode is differential and leads directly to a time-of-flight mass measurement. With the new design, the DIFP will separate multiple ion streams and analyze each stream independently for ion flux intensity, velocity (including direction of motion), mass, and temperature (or energy distribution). In particular, such an instrument will be invaluable on follow-on electrodynamic TSS missions and, possibly, for environmental monitoring on the space station.

Wright, K. H., Jr.↗
DOE OSTI2022

Measurements of the Higgs boson inclusive and differential fiducial cross-sections in the diphoton decay channel with pp collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A measurement of inclusive and differential fiducial cross-sections for the production of the Higgs boson decaying into two photons is performed using 139 fb –1 of proton-proton collision data recorded at $\sqrt{s}$ = 13 TeV by the ATLAS experiment at the Large Hadron Collider. The inclusive cross-section times branching ratio, in a fiducial region closely matching the experimental selection, is measured to be 67 ± 6 fb, which is in agreement with the state-of-the-art Standard Model prediction of 64 ± 4 fb. Extrapolating this result to the full phase space and correcting for the branching ratio, the total cross-section for Higgs boson production is estimated to be 58 ± 6 pb. In addition, the cross-sections in four fiducial regions sensitive to various Higgs boson production modes and differential cross-sections as a function of either one or two of several observables are measured. All the measurements are found to be in agreement with the Standard Model predictions. The measured transverse momentum distribution of the Higgs boson is used as an indirect probe of the Yukawa coupling of the Higgs boson to the bottom and charm quarks. In addition, five differential cross-section measurements are used to constrain anomalous Higgs boson couplings to vector bosons in the Standard Model effective field theory framework.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Separable physics-informed DeepONet: Breaking the curse of dimensionality in physics-informed machine learning

The deep operator network (DeepONet) has shown remarkable potential in solving partial differential equations (PDEs) by mapping between infinite-dimensional function spaces using labeled datasets. However, in scenarios lacking labeled data, the physics-informed DeepONet (PI-DeepONet) approach, which utilizes the residual loss of the governing PDE to optimize the network parameters, faces significant computational challenges, particularly due to the curse of dimensionality. This limitation has hindered its application to high-dimensional problems, making even standard 3D spatial with 1D temporal problems computationally prohibitive. Additionally, the computational requirement increases exponentially with the discretization density of the domain. Here, to address these challenges and enhance scalability for high-dimensional PDEs, we introduce the Separable physics-informed DeepONet (Sep-PI-DeepONet). This framework employs a factorization technique, utilizing sub-networks for individual one-dimensional coordinates, thereby reducing the number of forward passes and the size of the Jacobian matrix required for gradient computations. By incorporating forward-mode automatic differentiation (AD), we further optimize computational efficiency, achieving linear scaling of computational cost with discretization density and dimensionality, making our approach highly suitable for high-dimensional PDEs. We demonstrate the effectiveness of Sep-PI-DeepONet through three benchmark PDE models: the viscous Burgers’ equation, Biot’s consolidation theory, and a parameterized heat equation. Our framework maintains accuracy comparable to the conventional PI-DeepONet while reducing training time by two orders of magnitude. Notably, for the heat equation solved as a 4D problem, the conventional PI-DeepONet was computationally infeasible (estimated 289.35 h), while the Sep-PI-DeepONet completed training in just 2.5 h. These results underscore the potential of Sep-PI-DeepONet in efficiently solving complex, high-dimensional PDEs, marking a significant advancement in physics-informed machine learning.

Neural operator↗
DOE OSTI2022

Understanding the Impacts of Support–Polymer Interactions on the Dynamics of Poly(ethyleneimine) Confined in Mesoporous SBA-15

Supported amines are a promising class of CO 2 sorbents offering large uptake capacities and fast uptake rates. Additionally, among supported amines, poly(ethyleneimine) (PEI) physically impregnated in the mesopores of SBA-15 silica is widely used. Within these composite materials, the chain dynamics and morphologies of PEI strongly influence the CO 2 capture performance, yet little is known about chain and macromolecule mobility in confined pores. Here, we probe the impact of the support–PEI interactions on the dynamics and structures of PEI at the support interface and the corresponding impact on CO 2 uptake performance, which yields critical structure–property relationships. The pore walls of the support are grafted with organosilanes with different chemical end groups to differentiate interaction modes (spanning from strong attraction to repulsion) between the pore surface and PEI. Combinations of techniques, such as quasi-elastic neutron scattering (QENS), 1 H T 1 –T 2 relaxation correlation solid-state NMR, and molecular dynamics (MD) simulations, are used to comprehensively assess the physical properties of confined PEI. We hypothesized that PEI would have faster dynamics when subjected to less attractive or repulsive interactions. However, we discover that complex interfacial interactions resulted in complex structure–property relationships. Indeed, both the chain conformation of the surface-grafted chains and of the PEI around the surface influenced the chain mobility and CO 2 uptake performance. By coupling knowledge of the dynamics and distributions of PEI with CO 2 sorption performance and other characteristics, we determine that the macroscopic structures of the hybrid materials dictate the first rapid CO 2 uptake, and the rate of CO 2 sorption during the subsequent gradual uptake stage is determined by PEI chain motions that promote diffusive jumps of CO 2 through PEI-packed domains.

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

Augmenting subspace optimization methods with linear bandits

In this work, we consider the framework of methods for unconstrained minimization that are, in each iteration, restricted to a model that is only a valid approximation to the objective function on some affine subspace containing an incumbent point. These methods are of practical interest in computational settings where derivative information is either expensive or impossible to obtain. Recent attention has been paid in the literature to employing randomized matrix sketching for generating the affine subspaces within this framework. We consider a relatively straightforward, deterministic augmentation of such a generic subspace optimization method. In particular, we consider a sequential optimization framework where actions consist of one-dimensional linear subspaces and rewards consist of (approximations to) the magnitudes of directional derivatives computed in the direction of the action subspace. Reward maximization in this context is consistent with maximizing lower bounds on descent guaranteed by first-order Taylor models. This sequential optimization problem can be analysed through the lens of dynamic regret. We modify an existing linear upper confidence bound (UCB) bandit method and prove sublinear dynamic regret in the subspace optimization setting. We demonstrate the efficacy of employing this linear UCB method in a setting where forward-mode algorithmic differentiation can provide directional derivatives in arbitrary directions and in a derivative-free setting. For the derivative-free setting, we propose SS-POUNDers, an extension of the derivative-free optimization method POUNDers that employs the linear UCB mechanism to identify promising subspaces. Our numerical experiments suggest a preference, in either computational setting, for employing a linear UCB mechanism within a subspace optimization method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Efficient ptychographic phase retrieval via a matrix-free Levenberg-Marquardt algorithm

The phase retrieval problem, where one aims to recover a complex-valued image from far-field intensity measurements, is a classic problem encountered in a range of imaging applications. Modern phase retrieval approaches usually rely on gradient descent methods in a nonlinear minimization framework. Calculating closed-form gradients for use in these methods is tedious work, and formulating second order derivatives is even more laborious. Additionally, second order techniques often require the storage and inversion of large matrices of partial derivatives, with memory requirements that can be prohibitive for data-rich imaging modalities. We use a reverse-mode automatic differentiation (AD) framework to implement an efficient matrix-free version of the Levenberg-Marquardt (LM) algorithm, a longstanding method that finds popular use in nonlinear least-square minimization problems but which has seen little use in phase retrieval. Furthermore, we extend the basic LM algorithm so that it can be applied for more general constrained optimization problems (including phase retrieval problems) beyond just the least-square applications. Since we use AD, we only need to specify the physics-based forward model for a specific imaging application; the first and second-order derivative terms are calculated automatically through matrix-vector products, without explicitly forming the large Jacobian or Gauss-Newton matrices typically required for the LM method. We demonstrate that this algorithm can be used to solve both the unconstrained ptychographic object retrieval problem and the constrained “blind” ptychographic object and probe retrieval problems, under the popular Gaussian noise model as well as the Poisson noise model. We compare this algorithm to state-of-the-art first order ptychographic reconstruction methods to demonstrate empirically that this method outperforms best-in-class first-order methods: it provides excellent convergence guarantees with (in many cases) a superlinear rate of convergence, all with a computational cost comparable to, or lower than, the tested first-order algorithms.

36 MATERIALS SCIENCE↗
NASA NTRS1978

The airborne Laser Absorption Spectrometer - A new instrument of remote measurement of atmospheric trace gases

The Laser Absorption Spectrometer is a portable instrument developed by JPL for remote measurement of trace gases from an aircraft platform. It contains two carbon dioxide lasers, two optical heterodyne receivers, appropriate optics to aim the lasers at the ground and detect the backscattered energy, and signal processing and recording electronics. Operating in the differential-absorption mode, it is possible to monitor one atmospheric gas at a time and record the data in real time. The system can presently measure ozone, ethylene, water vapor, and chlorofluoromethanes with high sensitivity. Airborne measurements were made in early 1977 from the NASA/JPL twin-engine Beechcraft and in May 1977 from the NASA Convair 990 during the ASSESS-II Shuttle Simulation Study. These flights resulted in measurements of ozone concentrations in the lower troposphere which were compared with ground-based values provided by the Air Pollution Control District. This paper describes the details of the instrument and results of the airborne measurements.

Shumate, M. S.↗
NASA NTRS1986

Analytical determination of propeller performance degradation due to ice accretion

A computer code has been developed which is capable of computing propeller performance for clean, glaze, or rime iced propeller configurations, thereby providing a mechanism for determining the degree of performance degradation which results from a given icing encounter. The inviscid, incompressible flow field at each specified propeller radial location is first computed using the Theodorsen transformation method of conformal mapping. A droplet trajectory computation then calculates droplet impingement points and airfoil collection efficiency for each radial location, at which point several user-selectable empirical correlations are available for determining the aerodynamic penalities which arise due to the ice accretion. Propeller performance is finally computed using strip analysis for either the clean or iced propeller. In the iced mode, the differential thrust and torque coefficient equations are modified by the drag and lift coefficient increments due to ice to obtain the appropriate iced values. Comparison with available experimental propeller icing data shows good agreement in several cases. The code's capability to properly predict iced thrust coefficient, power coefficient, and propeller efficiency is shown to be dependent on the choice of empirical correlation employed as well as proper specification of radial icing extent.

Miller, T. L.↗
NASA NTRS1992

Emerging solid-state laser technology by lidar/DIAL remote sensing

Significant progress has been made in recent years in the development of new, solid-state laser sources. This talk will present an overview of some of the new developments in solid-state lasers, and their application toward lidar/DIAL measurements of the atmosphere. Newly emerging lasers such as Ho:YAG, Tm:YAG, OPO, and Ti:Sapphire will be covered, along with the spectroscopic parameters required for differential operational modes of atmospheric remote sensing including Doppler-Windshear lidar, Tunable laser detection of water/CO2, and broad linewidth OPO's for open path detection of pollutant hydrocarbon gases. Additional considerations of emerging laser technology for lidar/DIAL will also be covered.

Killinger, Dennis↗
NASA NTRS2006

The Potential of Water Vapor & Precipitation Estimation with a Differential-frequency Radar

In the presence of rain, the radar return powers from a three-frequency radar, with center frequency at 22.235 GHz and upper and lower frequencies chosen with equal water vapor absorption coefficients, can be used to estimate water vapor density and parameters of the precipitation. A linear combination of differential measurements between the center and lower frequencies on one hand and the upper and lower frequencies on the other provide an estimate of differential water vapor absorption. Conversely, the difference in radar reflectivity factors (in dB) between the upper and lower frequencies is independent of water vapor absorption and can be used to estimate the median mass diameter of the hydrometeors. For a down-looking radar, path-integrated estimates of water vapor absorption may be possible under rain-free as well as raining conditions by using the surface returns at the three frequencies. Cross-talk or interference between the precipitation and water vapor estimates depends on the frequency separation of the channels as well as on the phase state and the median mass diameter of the hydrometeors. Simulations of the retrieval of water vapor absorption show that the largest source of variability arises from the variance in the measured radar return powers while the largest biases occur in the mixed-phase region. Use of high pulse repetition frequencies and signal whitening methods may be needed to obtain the large number of independent samples required. Measurements over a fractional bandwidth, defined as the ratio of the difference between the upper and lower frequencies to the center frequency, up to about 0.2 should be passible in a differential frequency mode, where a single transceiver and antenna are used. Difficulties in frequency allocation may require alternative choices of frequency where the water vapor absorptions at the low and high frequencies are unequal. We consider the degradation in the retrieval accuracy when the frequencies are not optimum.

Meneghini, Robert↗
DOE OSTI2023

Apparatus and methods of detecting transient discharge modes and/or closed loop control of pulsed systems and method employing same

A system and method for differentiating between different modes of pulsed electrical discharges via of an amplitude to time (ATC) conversion circuit is described. A bipolar ATC circuit is used to add together the positive and negative portions of an attenuated and filtered signal derived either from the voltage or current of a pulse. Alternatively, a unipolar ATC circuit may be employed. The resulting processed signal is compared against a reference voltage to generate an output signal that is active for the amount of time that the processed signal exceeds the reference voltage. Discharge mode is determined based on three factors: did a pulse occur, if a pulse occurred when did the pulse start relative to the original pulse event, and what is the duty cycle of the pulse. Subsequent pulse generated may be controlled accordingly.

Fitzpatrick, Joseph F.↗
NASA NTRS1992

Convective instability in differentially rotating disks

A normal mode analysis for nonaxisymmetric perturbations in a thin, differentially rotating disk with a vertical structure that is isothermal and convectively unstable is performed. The vertical gravity is assumed to be external and constant. The perturbation scale is assumed to be much shorter than the radius of the disk but comparable to or less than the thickness. The initial value problem is formulated in shearing coordinates. Dispersion relations are obtained for the three limiting cases of zero shear, axisymmetric perturbations, and small radial wavelengths. The full effects of shear are studied by integrating numerically the initial value problem. Nonaxisymmetric local Fourier modes are found to have a radial wavenumber that increases linearly with time in proportion to the shear times the azimuthal wavenumber. While Coriolis forces exert stabilizing effects on the convective modes, reducing their growth rate and the range of unstable wavelengths, shear has destablizing effects inasmuch as it reduces the epicyclic frequency at a given angular velocity. In a Keplerian disk, perturbations with azimuthal wavelengths about 2 times smaller than vertical wavelengths grow exponentially.

Ryu, Dongsu↗
NASA NTRS1973

An experimental evaluation of metallic diaphragms for positive fuel expulsion in the Atmosphere Explorer, hydrazine propulsion subsystem

Four Arde conospheroid metallic diaphragms were tested at NASA to evaluate their capability for use in the orbit adjust propulsion subsystem (OAPS) of the Atmosphere Explorer spacecraft (AE's C, D, and E). The diaphragms will be used for positive propellant expulsion and spacecraft center of mass (c.m.) control. A leak-free cycle life capability of nine reversals was demonstrated. The diaphragms rolled smoothly from ring to ring in a predictable manner on the first reversal. Varying amounts of diaphragm cocking and ring skipping were observed on subsequent reversals. The diaphragm pressure differential did not exceed 10 psid during any reversal. Cycle life capability, reversal mode, and pressure differential were not affected by sudden reversals, environmental tests, or 18,000 partial reversals. An expulsion efficiency of approximately 97 percent was demonstrated. The results of these tests show that metallic diaphragms can be used as an effective means of positive fuel expulsion; however, to achieve spacecraft c.m. control, the diaphragm must not be reversed prior to flight.

Woodruff, W. L.↗