Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gradient information”

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 631 records · Page 35

Experiments in Neural-Network Control of a Free-Flying Space Robot

Four important generic issues are identified and addressed in some depth in this thesis as part of the development of an adaptive neural network based control system for an experimental free flying space robot prototype. The first issue concerns the importance of true system level design of the control system. A new hybrid strategy is developed here, in depth, for the beneficial integration of neural networks into the total control system. A second important issue in neural network control concerns incorporating a priori knowledge into the neural network. In many applications, it is possible to get a reasonably accurate controller using conventional means. If this prior information is used purposefully to provide a starting point for the optimizing capabilities of the neural network, it can provide much faster initial learning. In a step towards addressing this issue, a new generic Fully Connected Architecture (FCA) is developed for use with backpropagation. A third issue is that neural networks are commonly trained using a gradient based optimization method such as backpropagation; but many real world systems have Discrete Valued Functions (DVFs) that do not permit gradient based optimization. One example is the on-off thrusters that are common on spacecraft. A new technique is developed here that now extends backpropagation learning for use with DVFs. The fourth issue is that the speed of adaptation is often a limiting factor in the implementation of a neural network control system. This issue has been strongly resolved in the research by drawing on the above new contributions.

NEURAL NETS↗

Computational modeling of charge hopping dynamics along a disordered one-dimensional wire with energy gradients in quantum environments

This computational study investigates the effects of energy gradients on charge hopping dynamics along a one-dimensional chain of discrete sites coupled to quantum bath, which is modeled at the level of Pauli master equation (PME). This study also assesses the performance of different approximations for the hopping rates. Three different methods for solving the PME, a fourth order Runge–Kutta method, numerical diagonalization of the rate matrix followed by analytic propagation, and kinetic Monte Carlo simulation method, are tested and confirmed to produce virtually identical values of time dependent mean square displacement, diffusion constant, and mobility. Five different rate expressions, exact numerical evaluation of Fermi’s Golden Rule (FGR) rate, stationary phase interpolation (SPI) approximation, semiclassical approximation, classical Marcus rate, and Miller–Abrahams rate, are tested to help understand the effects of approximations in representing quantum environments in the presence of energy gradients. The results based on direct numerical evaluation of FGR rate exhibit transition from diffusive to non-diffusive behavior with the increase in the gradient and show that the charge transport in the quantum bath is more sensitive to the magnitude of the gradient and the disorder than in the classical bath. Among all the four approximations for the hopping rates, the SPI approximation is confirmed to work best overall. A comparison of two different methods to calculate the mobility identifies drift motion of the population distribution as the major source of non-diffusive behavior and provides more reliable information on the contribution of quantum bath.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Noise-Resilient Quantum Machine Learning for Stability Assessment of Power Systems

Transient stability assessment (TSA) is a cornerstone for resilient operations of todays interconnected power grids. This paper is a confluence of quantum computing, data science and machine learning to potentially address the power system TSA issue. Here, we devise a quantum TSA (QTSA) method to enable scalable and efficient data-driven transient stability prediction for bulk power systems, which is the first attempt to tackle the TSA issue with quantum computing. Our contributions are three-fold: 1) A high expressibility, low-depth (HELD) quantum circuit is designed for accurate and noise-resilient TSA; 2) A quantum natural gradient descent algorithm is developed for efficient HELD circuit training; 3) A systematical analysis on QTSAs performance under various quantum factors is per-formed. QTSA underpins a foundation of quantum-enabled and data-driven power grid stability analytics. It renders the intractable TSA straightforward and effortless in the Hilbert space, and therefore provides stability information for power system operations. Extensive experiments on quantum simulators and real quantum computers verify the accuracy, noise-resilience, scalability and universality of QTSA.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

On the temperature and velocity through the photosphere of a sunspot penumbra

We investigate the structure in depth of a sunspot penumbra by means of the inversion code of the radiative transfer equation proposed by Ruiz Cobo & del Toro Iniesta (1992), applied to a set of filtergrams of a sunspot, scanning the Fe I line at 5576.1 A, with a sampling interval of 30 mA, from -120 to 120 mA from line center (data previously analyzed by Title et al. 1993). The temperature structure of this penumbra is obtained for each of the 801 pixels selected (0.32 sec x 0.32 sec). On the average, the temperatures seem to decrease as we move inward, but the differences are of the order of the rms values (approximately equal 100-200 K) at a given distance to sunspot center. The outer parts of the penumbra have also a bigger curvature in the T versus log tau(sub 5) relation than the inner parts. We realize, however, that these differences might be influenced by possible stray light effects. Compared to the quiet Sun, penumbral temperatures are cooler at deep layers and hotter at high layers. A mean penumbral model atmosphere is presented. The asymmetries observed in the intensity profile (the line is magnetically insensitive) are deduced to be produced by strong gradients of the line-of-sight velocity that sharply vary spatially along slices of almost constant distance to sunspot center. These variations suggest that such gradients are not only needed to explain the broadband circular polarization observed in sunspots (see Sanchez Almeida & Lites 1992) but are a main characteristic of the fine-scale penumbra. The results are compatible with an Evershed flow present everywhere, but its gradient with depth turns out to vary so that the flow seems to be mainly concentrated in some penumbral fibrils when studied through Dopplergrams. Finally, as by-products of this study, we put constraints to the practical usefulness of the Eddington-Barbier relation, and we explain the values of the Fourier Dopplergrams to be carrying information of layers around the centroid of the generalized response function of Dopplergrams to velocity fluctuations.

Del Toro Iniesta, J. C.↗

The Derivation of the Gradient of the Acoustic Pressure on a Moving Surface for Application to the Fast Scattering Code (FSC)

In this paper we develop an analytic expression for calculation of the the acoustic pressure from a rotating blade on a moving surface for application to the Fast Scattering Code (FSC). The analytic result is intended to be used in the helicopter noise prediction code PSU-WOPWOP. One of the goals of the derivation is obtaining a result that will not use any more information than are needed for the prediction of the thickness and loading noise. The result derived here achieves this goal and its incorporation in PSU-WOPWOP is straight forward and attainable.

Farassat, F.↗

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↗

Fluid flow and heat convection studies for actively cooled airframes

This report details progress made on the jet impingement - liquid crystal - digital imaging experiment. With the design phase complete, the experiment is currently in the construction phase. In order to reach this phase two design related issues were resolved. The first issue was to determine NASP leading edge active cooling design parameters. Meetings were arranged with personnel at SAIC International, Torrance, CA in order to obtain recent publications that characterized expected leading edge heat fluxes as well as other details of NASP operating conditions. The information in these publications was used to estimate minimum and maximum jet Reynolds numbers needed to accomplish the required leading edge cooling, and to determine the parameters of the experiment. The details of this analysis are shown in Appendix A. One of the concerns for the NASP design is that of thermal stress due to large surface temperature gradients. Using a series of circular jets to cool the leading edge will cause a non-uniform temperature distribution and potentially large thermal stresses. Therefore it was decided to explore the feasibility of using a slot jet to cool the leading edge. The literature contains many investigations into circular jet heat transfer but few investigations of slot jet heat transfer. The first experiments will be done on circular jets impinging on a fiat plate and results compared to previously published data to establish the accuracy of the method. Subsequent experiments will be slot jets impinging on full scale models of the NASP leading edge. Table 1 shows the range of parameters to be explored. Next a preliminary design of the experiment was done. Previous papers which used a similar experimental technique were studied and elements of those experiments adapted to the jet impingement study. Trade-off studies were conducted to determine which design was the least expensive, easy to construct, and easy to use. Once the final design was settled, vendors were contacted to verify that equipment could be obtained to meet our specifications. Much of the equipment required to complete the construction of the experiment has been ordered or received. The material status list is shown in Appendix B.

Mills, A. F.↗

Large gradients via correlation in random parameterized quantum circuits

Scaling of variational quantum algorithms to large problem sizes requires efficient optimization of random parameterized quantum circuits. For such circuits with uncorrelated parameters, the presence of exponentially vanishing gradients in cost function landscapes is an obstacle to optimization by gradient descent methods. In this work, we prove that reducing the dimensionality of the parameter space by utilizing circuit modules containing spatially or temporally correlated gate layers can allow one to circumvent the vanishing gradient phenomenon. Here, examples are drawn from random separable circuits and asymptotically optimal variational versions of Grover's algorithm based on the quantum alternating operator ansatz. In the latter scenario, our bounds on cost function variation imply a transition between vanishing gradients and efficient trainability as the number of layers is increased toward $\mathcal{O}\left({2}^{n/2}\right)$, the optimal oracle complexity of quantum unstructured search.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Gas sampling method for determining pollutant concentrations in the flame zone of two swirl-can combustor modules

A gas sampling probe and traversing mechanism were developed to obtain detailed measurements of gaseous pollutant concentrations in the primary and mixing regions of combustors in order to better understand how pollutants are formed. The gas sampling probe was actuated by a three-degree-of-freedom traversing mechanism and the samples obtained were analyzed by an on-line gas analysis system. The pollutants in the flame zone of two different swirl-can combustor modules were measured at an inlet-air temperature of 590 K, pressure of 6 atmospheres, and reference velocities of 23 and 30 meters per second at a fuel-air ratio of 0.02. Typical results show large spatial gradients in the gaseous pollutant concentration close to the swirl-can module. Average concentrations of unburned hydrocarbons and carbon monoxide decrease rapidly in the downstream wake regions of each module. By careful and detailed probing, the effect of various module design features on pollutant formation can be assessed. The techniques presently developed seem adequate to obtain the desired information.

Duerr, R. A.↗

Meteorite spectroscopy and characterization of asteroid surface materials

The analysis of visible and near-infrared reflectance spectra is the primary means to determine surface mineralogy and petrology of individual asteroids. These individual studies provide the data to investigate the broader relationships between the asteroids and meteorites and between asteroids at different heliocentric distances. The main purpose is to improve the understanding of the origin, evolution, and inter-relationships of the asteroids; of their relationships to the meteorites; and of the processes active and the conditions present in the early inner solar system. Empirical information from the study of asteroids and the meteorites is essential to the adequate development and testing of the theoretical models for the accretion of the terrestrial planets, and for their early post-accretionary evolution. The recent results are outined in the following sections: (1) asteroid igneous processes, and (2) spinel-bearing asteroids and the nebular compositional gradient.

Gaffey, Michael J.↗

Characterization of radiation-induced damage in high performance polymers by electron paramagnetic resonance imaging spectroscopy

The potential for long-term human activity beyond the Earth's protective magnetosphere is limited in part by the lack of detailed information on the effectiveness and performance of existing structural materials to shield the crew and spacecraft from highly penetrating space radiations. The two radiations of greatest concern are high energy protons emitted during solar flares and galactic cosmic rays which are energetic ions ranging from protons to highly oxidized iron. Although the interactions of such high-energy radiations with matter are not completely understood at this time, the effects of the incident radiation are clearly expected to include the formation of paramagnetic spin centers via ionization and bond-scission reactions in the molecular matrices of structural materials. Since this type of radiation damage is readily characterized by Electron Paramagnetic Resonance (EPR) spectroscopy, the NASA Langley Research Center EPR system was repaired and brought on-line during the 1991 ASEE term. A major goal of the 1992 ASEE term was to adapt the existing core of the LaRC EPR system to meet the requirements for EPR Imaging--a powerful new technique which provides detailed information on the internal structure of materials by mapping the spatial distribution of unpaired spin density in bulk media. Major impetus for this adaptation arises from the fact that information derived from EPRI complements other methods such as scanning electron microscopy which primarily characterize surface phenomena. The modification of the EPR system has been initiated by the construction of specially designed, counterwound Helmholtz coils which will be mounted on the main EPR electromagnet. The specifications of the coils have been set to achieve a static linear magnetic field gradient of 10 gauss/mm/amp along the principal (Z) axis of the Zeeman field. Construction is also in progress of a paramagnetic standard in which the spin distribution is known in all three dimensions. This sample will be used to assess the linearity of the magnetic field gradient and to ensure authentic image reconstruction. A second major task was to secure the computer capability to enable image reconstruction from projection data generated by the magnetic field gradients. To this end, commercially available and public domain software packages which perform inverse Fourier Transform and convoluted (filtered) back projection functions are being integrated into the existing EPR data processing system.

Suleman, Naushadalli K.↗

HfTa_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-tantalum (Hf-Ta). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ta. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice sites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,040 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗

Development of a molecularly informed biogeochemical framework for reactive transport modeling of subsurface carbon inventories, transformations and fluxes (Final Report)

The overall objective of the project was to combine new molecular-level characterization strategies with soil carbon flux measurements to develop and evaluate model representations of subsurface carbon cycling. We expanded on existing studies in the East River watershed, Colorado, in collaboration with the Berkeley Lab Watershed Function Scientific Focus Area (SFA), SLAC Groundwater Quality SFA, and Rocky Mountain Biological Laboratory (RMBL) to develop an elevation and vegetation gradient that is now the subject of long-term monitoring by the USGS. To achieve the overall objective, we combined field studies of soil respiration with laboratory analyses, ranging from spectroscopy to incubation studies. This combination of techniques enabled us to develop a new understanding of the drivers of high-elevation soil respiration. In the process, we developed three new modeling approaches that improve our ability to conceptualize and ultimately to represent soil respiration in numerical models. The first approach is a plot-scale transient inverse model that can be used to determine in situ CO 2 production rates from measured concentration profiles and surface fluxes. Application of this method revealed the importance of plant phenology and deep CO 2 production in moderating CO 2 fluxes to the atmosphere. The second modeling approach is a molecular-scale tool that enables spectroscopic and elemental data for carbon speciation to be transformed into functional group abundances, or the ‘SOC-fga model’. This method uniquely enables carbon speciation to be tracked within a reactive transport framework to partition carbon among different pathways and storage zones within the soil. The third modeling tool builds strongly on the previous approaches and captures the microbial processes driving heterotrophic respiration. This ‘dormancy model’ allows the native soil microbial population to respond transiently to the presence or absence of water in order to catalyze carbon respiration. This approach was also compared to the simpler and more widely used first-order model using two experimental datasets with different temporal and spatial resolutions. Our results illustrate that the simpler first order model provides a robust and efficient representation of deep (>1 meter) soil respiration, but that shallow soils, where most respiration occurs, require explicit representation of moisture-dependent activation and dormancy rates. By assessing soil organic carbon turnover at multiple scales, we see a complex array of controls emerge. Atte scale of a hillslope, spatial heterogeneity in soil respiration rates dominates and is uncorrelated with instantaneous soil moisture and plant community. At the profile scale, the balance between plant inputs and water availability is the dominant control. At the microbial to molecular scale, physiological processes associated with carbon use and carbon speciation are important controls. These scale-dependent controls emphasize the need for new modeling approaches that examine their interactions and hierarchies.

54 ENVIRONMENTAL SCIENCES↗

Reinforcement Learning-based Output Structured Feedback for Distributed Multi-Area Power System Frequency Control

Load frequency control (LFC) is a key factor to maintain the stable frequency in multi-area power systems. As the modern power systems evolve from centralized to decentralized paradigm, LFC needs to consider the decentralized scheme that considers limited information from the information-exchange graph for the generator control of each interconnected area. This paper aims to solve a data-driven constrained LQR problem with mean-variance risk constraints and output structured feedback, and applies this framework to solve the LFC problem in multi-area power systems. By reformulating the constrained optimization problem into a minimax problem, the stochastic gradient descent max-oracle (SGDmax) algorithm with zero-order policy gradient (ZOPG) is adopted to find the optimal feedback gain from the learning, while guaranteeing the convergence. In addition, to improve the adaptation of the proposed learning method to new or varying models, we construct an emulator grid that approximates the dynamics of a physical grid and performs training based on this model. Once the feedback gain is obtained from the emulator grid, it is applied to the physical grid with a robustness test to check whether the controller from the approximated emulator applies to the actual system. Numerical tests show that the obtained feedback controller can successfully control the frequency of each area, while mitigating the uncertainty from the loads, with reliable robustness that ensures the adaptability of the obtained feedback gain to the actual physical grid.

Kwon, Kyung-bin↗

New results and techniques in space radio astronomy.

The methods and results of early space radioastronomy experiments are reviewed, with emphasis on the RAE 1 spacecraft which was designed specifically and exclusively for radio astronomical studies. The RAE 1 carries two gravity-gradient-stabilized 229-m traveling-wave V-antennas, a 37-m dipole antenna, and a number of radiometer systems to provide measurements over the 0.2 to 9.2 MHz frequency range with a time resolution of 0.5 sec and an absolute accuracy of plus or minus 25%. Observations of solar bursts at frequencies down to 0.2 MHz provide new information on the density, plasma velocity, and dynamics of coronal streamers out to distances greater than 50 solar radii. New information on the distribution of the ionized component of the interstellar medium is being obtained from galactic continuum background maps at frequencies around 4 MHz. Cosmic noise background spectra measured down to 0.5 MHz produce new estimates on the interstellar flux of cosmic rays, on magnetic fields in the galactic halo, and on distant extragalactic radio sources.

Alexander, J. K.↗

30 years of adaptive neural networks - Perceptron, Madaline, and backpropagation

Fundamental developments in feedforward artificial neural networks from the past thirty years are reviewed. The history, origination, operating characteristics, and basic theory of several supervised neural-network training algorithms (including the perceptron rule, the least-mean-square algorithm, three Madaline rules, and the backpropagation technique) are described. The concept underlying these iterative adaptation algorithms is the minimal disturbance principle, which suggests that during training it is advisable to inject new information into a network in a manner that disturbs stored information to the smallest extent possible. The two principal kinds of online rules that have developed for altering the weights of a network are examined for both single-threshold elements and multielement networks. They are error-correction rules, which alter the weights of a network to correct error in the output response to the present input pattern, and gradient rules, which alter the weights of a network during each pattern presentation by gradient descent with the objective of reducing mean-square error (averaged over all training patterns).

Widrow, Bernard↗

Rheology, tectonics, and the structure of the Venus lithosphere

Given the absence of ground truth information on seismic structure, heat flow, and rock strength, or short wavelength gravity or magnetic data for Venus, information on the thermal, mechanical and compositional nature of the shallow interior must be obtained by indirect methods. Using pre-Magellan data, theoretical models constrained by the depths of impact craters and the length scales of tectonic features yielded estimates on the thickness of Venus' brittle-elastic lithosphere and the allowable range of crustal thickness and surface thermal gradient. The purpose of this study is to revisit the question of the shallow structure of Venus based on Magellan observations of the surface and recent experiments that address Venus' crustal rheology.

Zuber, M. T.↗

Solar wind velocity and temperature in the outer heliosphere

At the end of 1992, the Pioneer 10, Pioneer 11, and Voyager 2 spacecraft were at heliocentric distances of 56.0, 37.3, and 39.0 AU and heliographic latitudes of 3.3 deg N, 17.4 deg N, and 8.6 deg S, respectively. Pioneer 11 and Voyager 2 are at similar celestial longitudes, while Pioneer 10 is on the opposite side of the Sun. All three spacecraft have working plasma analyzers, so intercomparison of data from these spacecraft provides important information about the global character of the solar wind in the outer heliosphere. The averaged solar wind speed continued to exhibit its well-known variation with solar cycle: Even at heliocentric distances greater than 50 AU, the average speed is highest during the declining phase of the solar cycle and lowest near solar minimum. There was a strong latitudinal gradient in solar wind speed between 3 deg and 17 deg N during the last solar minimum, but this gradient has since disappeared. The solar wind temperature declined with increasing heliocentric distance out to a heliocentric distance of at least 20 AU; this decline appeared to continue at larger heliocentric distances, but temperatures in the outer heliosphere were suprisingly high. While Pioneer 10 and Voyager 2 observed comparable solar wind temperatures, the temperature at Pioneer 11 was significantly higher, which suggests the existence of a large-scale variation of temperature with heliographic longitude. There was also some suggestion that solar wind temperatures were higher near solar minimum.

Gazis, P. R.↗