Search NASA⌕ Search

SEARCH · Search NASA

Results for “automatic differentiation”

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 289 records · Page 16

Emulation and detection of physical faults and cyber-attacks on building energy systems through real-time hardware-in-the-loop experiments

The increasing use of remote or mobile access, integrated wearable technologies, data exchange, and cloud-based data analytics in modern smart buildings is steering the building industry towards open communication technologies. The increased connectivity and accessibility could lead to more cyber-attacks in smart buildings. On the other hand, physical faults (e.g., HVAC -heating, ventilation, and air-conditioning faults) may have similar adverse impacts as those from the cyber-attacks on building energy systems, such as occupant discomfort, energy wastage, and equipment downtime. However, current physical behavior-based anomaly detection methods fail to differentiate between cyber-attacks and physical faults in building energy systems. Moreover, the challenge in collecting real-world threat data with ground truth has led researchers to rely on numerical models with user-defined assumptions, which may not accurately reflect real-world conditions due to the lack of in-situ experimental datasets. To address these challenges and gaps, this paper presents a flexible hardware-in-the-loop (HIL) testbed for generating cyber-attack and physical fault datasets and demonstrating threat detection algorithms in a real building automation system (BAS) environment. This testbed combines hardware (i.e., real BAS with local HVAC controllers and a physical network) with software (i.e., high-fidelity models to represent behaviors of building envelope and HVAC energy systems), enabling emulations of realistic threats. Five HIL experiments, including one baseline without any threats, two with physical faults, and two with cyber-attacks, were conducted to generate datasets containing detailed network traffic and system states. A joint classification framework, incorporating a network analyzer and a physical HVAC fault detector, was proposed to automatically detect cyber-physical abnormalities on BAS at both the network and the physical HVAC levels. The network analyzer comprises a conditional random fields (CRF) based command validator and a statistics-based detection strategy. The fault detector employs a weather and schedule-based pattern matching and feature-based principal component analysis (WPM-FPCA) method. Evaluation of the classification using four metrics from the multi-class confusion matrix revealed an average accuracy of 90.2%, recall of 89.7%, precision of 88.5% and F1-score of 89.2%. Finally, these results demonstrate that the proposed joint classification framework can effectively differentiate between specific types of cyber-attacks (e.g., device reinitialization attack, network Denial-of-Service attack) and physical faults (e.g., air handling unit operational fault, cooling coil valve stuck) in real time for improved building energy management.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Inferring the spatial and energy distribution of gamma-ray burst sources. 1: Methodology

We describe Bayesian methods for analyzing the distribution of gamma-ray burst peak photon fluxes and directions. These methods fit the differential distribution, and have the following advantages over rival methods: (1) they do not destroy information by binning or averaging the data (as do, say, chi squared, the averaged value of V/V(sub max), and angular moment analyses); (2) they straightforwardly handle uncertainties in the measured quantities; (3) they analyze the strength and direction information jointly; (4) they use information available about nondetections; and (5) they automatically identify and account for biases and selection effects given a precise description of the experiment. In these methods, the most important information needed about the instrument threshold is not its value at the times of burst triggers, as is used in the average value of V/V(sub max) analyses, but rather the value of the threshold at times when no trigger occurred. We show that this information can be summarized as an average detection efficiency that is similar to the product of the exposure and efficiency reported in the First Burst and Transient Source Experiment (BATSE) Burst (1B) Catalog, but significantly different from it at low fluxes. We also quantify an important bias that results from estimating the peak flux by scanning the burst to find the peak number of counts in a window of specified duration, as was done for the 1B Catalog. When the duration of the peak of the light curve is longer than the window duration, a simple flux estimate based on the peak counts significantly overestimates the peak flux in a nonlinear fashion that distorts the shape of the log(N)-log(P) distribution. This distortion also corrupts analyses of the V/V(sub max) distribution that use ratios of counts above background to estimate V/V(sub max). The Bayesian calculation specifies how to account for this bias. Implementation of the Bayesian approach requires some changes in the way burst data are reported that we describe in detail. Subsequent papers will report analyses of the 1B Catalog data using the methods described here.

Loredo, Thomas J.↗

Time and frequency domain analysis of sampled data controllers via mixed operation equations

Specification of the mathematical equations required to define the dynamic response of a linear continuous plant, subject to sampled data control, is complicated by the fact that the digital components of the control system cannot be modeled via linear ordinary differential equations. This complication can be overcome by introducing two new mathematical operations; namely, the operation of zero order hold and digial delay. It is shown that by direct utilization of these operations, a set of linear mixed operation equations can be written and used to define the dynamic response characteristics of the controlled system. It also is shown how these linear mixed operation equations lead, in an automatable manner, directly to a set of finite difference equations which are in a format compatible with follow on time and frequency domain analysis methods.

Frisch, H. P.↗

Techniques for Automated Single-Star Photometry

Tennessee State University operates several automatic photoelectric telescopes (APTs) located at Fairborn Observatory in the Patagonia Mountains of southern Arizona. The APTs are dedicated to photometric monitoring programs that would be expensive and difficult to accomplish without the advantages of automation. I review the operation of two of the telescopes (a 0.75- and 0.80-meter (m) APT) and the quality-control techniques that result in the routine acquisition of single-star differential photometry with a precision of 0.001 mag for single observations and 0.0001 - 0.0002 mag for seasonal means. I also illustrate the capabilities of the APTS with sample results from a program to measure luminosity cycles in sun-like stars and a related program to search for the signatures of extrasolar planets around these stars.

Henry, Gregory W.↗

Numerical optimization in Hilbert space using inexact function and gradient evaluations

Trust region algorithms provide a robust iterative technique for solving non-convex unstrained optimization problems, but in many instances it is prohibitively expensive to compute high accuracy function and gradient values for the method. Of particular interest are inverse and parameter estimation problems, since function and gradient evaluations involve numerically solving large systems of differential equations. A global convergence theory is presented for trust region algorithms in which neither function nor gradient values are known exactly. The theory is formulated in a Hilbert space setting so that it can be applied to variational problems as well as the finite dimensional problems normally seen in trust region literature. The conditions concerning allowable error are remarkably relaxed: relative errors in the gradient error condition is automatically satisfied if the error is orthogonal to the gradient approximation. A technique for estimating gradient error and improving the approximation is also presented.

Carter, Richard G.↗

Automated Scoring of Morphological Changes in Images of Pentaerythritol Tetranitrate

Recent advances in characterization techniques that generate large datasets of material microstructure images require robust, automated image-processing. We applied an unsupervised anomaly detection method called feature anomaly detection system (FADS) to automatically detect and quantify microstructure changes in images of the explosive pentaerythritol tetranitrate (PETN) aged at various temperatures. We demonstrated the FADS approach on two-dimensional images extracted from computed tomography scans, but the same technique can be readily applied to other imaging modalities. FADS calculates anomaly scores on the basis of differences in filter activations of nominal and test data in pretrained convolutional neural networks. The FADS scores successfully differentiated between pristine PETN and PETN aged at a temperature where material coarsening occurred. Morphological metric analysis of segmented images verified observed trends in FADS scores as a function of aging temperature and aging time, specifically by calculating volume fractions, specific boundary lengths, two-point correlation functions, and local thicknesses. Here, the FADS technique has two important advantages compared to traditional morphological analysis: First, it uses grayscale images as input, rather than images that are segmented to separate the appropriate phases; and second, FADS scores capture any type of changes among image sets, rather than requiring prior knowledge or selection of a relevant set of metrics.

Accelerated aging↗

Parallel sorting algorithm classification: is manual instrumentation necessary?

Understanding parallel algorithms is crucial for accelerating scientific simulations on complex, distributed memory, high-performance computers. Modern algorithm classification approaches learn semantics directly from source code to differentiate between algorithms, however, accessing source code is not always possible. We can learn about parallel algorithms from observing their performance, as programs running the same algorithms and using the same hardware should exhibit similar performance characteristics. We present an approach to learn algorithm classes from parallel performance data directly in order to classify algorithms without access to the source code. We extend previous work to enable classifying parallel sorting algorithms using automatic instrumentation instead of requiring manual region annotations in the source code. In this work, we design and demonstrate a study for classification of parallel sorting algorithms using parallel performance data collected from automatic instrumentation, and evaluate the performance of our new methodology on classification. We leverage Caliper to collect the performance data, Thicket for our exploratory data analysis (EDA), and PyTorch and Scikit-learn to evaluate the effectiveness of random forests, support vector machines (SVMs), decision trees, neural networks, and logistic regressions on parallel performance data. Additionally, we study noise in parallel performance data, whether the removal of noise and pre-processing of the data is necessary to accurately classify parallel sorting algorithms, and determine the effectiveness of features created from performance data. In conclusion, we demonstrate classification accuracy for these five different models of up to 97.7% across four different parallel algorithm classes.

Algorithm Classification↗

Application and development of the Schwinger multichannel scattering theory and the partial differential equation theory of electron-molecule scattering

One version of the multichannel theory for electron-target scattering based on the Schwinger variational principle, the SMC method, requires the introduction of a projection parameter. The role of the projection parameter a is investigated and it is shown that the principal-value operator in the SMC equation is Hermitian regardless of the value of a as long as it is real and nonzero. In a basis that is properly orthonormalizable, the matrix representation of this operator is also Hermitian. The use of such basis is consistent with the Schwinger variational principle because the Lippmann-Schwinger equation automatically builds in the correct boundary conditions. Otherwise, an auxiliary condition needs to be introduced, and Takatsuka and McKoy's original value of a is one of the three possible ways to achieve Hermiticity. In all cases but one, a can be uncoupled from the Hermiticity condition and becomes a free parameter. An equation for a based on the variational stability of the scattering amplitude is derived; its solution has an interesting property that the scattering amplitude from a converged SMC calculation is independent of the choice of a even though the SMC operator itself is a-dependent. This property provides a sensitive test of the convergence of the calculation. For a static-exchange calculation, the convergence requirement only depends on the completeness of the one-electron basis, but for a general multichannel case, the a-invariance in the scattering amplitude requires both the one-electron basis and the N plus 1-electron basis to be complete. The role of a in the SMC equation and the convergence property are illustrated using two examples: e-CO elastic scattering in the static-exchange approximation, and a two-state treatment of the e-H2 Chi(sup 1)Sigma(sub g)(+) yields b(sup 3)Sigma(sub u)(+) excitation.

Weatherford, Charles A.↗

Numerical solution of potential flow about arbitrary 2-dimensional multiple bodies

A procedure for the finite-difference numerical solution of the lifting potential flow about any number of arbitrarily shaped bodies is given. The solution is based on a technique of automatic numerical generation of a curvilinear coordinate system having coordinate lines coincident with the contours of all bodies in the field, regardless of their shapes and number. The effects of all numerical parameters involved are analyzed and appropriate values are recommended. Comparisons with analytic solutions for single Karman-Trefftz airfoils and a circular cylinder pair show excellent agreement. The technique of application of the boundary-fitted coordinate systems to the numerical solution of partial differential equations is illustrated.

Thompson, J. F.↗

Machine Learning Approaches to Increasing Value of Spaceflight Omics Databases

The number of spaceflight bioscience mission opportunities is too small to allow all relevant biological and environmental parameters to be experimentally identified. Simulated spaceflight experiments in ground-based facilities (GBFs), such as clinostats, are each suitable only for particular investigations -- a rotating-wall vessel may be 'simulated microgravity' for cell differentiation (hours), but not DNA repair (seconds) -- and introduce confounding stimuli, such as motor vibration and fluid shear effects. This uncertainty over which biological mechanisms respond to a given form of simulated space radiation or gravity, as well as its side effects, limits our ability to baseline spaceflight data and validate mission science. Machine learning techniques autonomously identify relevant and interdependent factors in a data set given the set of desired metrics to be evaluated: to automatically identify related studies, compare data from related studies, or determine linkages between types of data in the same study. System-of-systems (SoS) machine learning models have the ability to deal with both sparse and heterogeneous data, such as that provided by the small and diverse number of space biosciences flight missions; however, they require appropriate user-defined metrics for any given data set. Although machine learning in bioinformatics is rapidly expanding, the need to combine spaceflight/GBF mission parameters with omics data is unique. This work characterizes the basic requirements for implementing the SoS approach through the System Map (SM) technique, a composite of a dynamic Bayesian network and Gaussian mixture model, in real-world repositories such as the GeneLab Data System and Life Sciences Data Archive. The three primary steps are metadata management for experimental description using open-source ontologies, defining similarity and consistency metrics, and generating testing and validation data sets. Such approaches to spaceflight and GBF omics data may soon enable unique insight into which measured phenomena correlate to biological mechanisms that are truly affected by spaceflight conditions; which are most likely to be confounded by other variables; and which are insufficiently characterized, significantly increasing existing and future science return from ISS and spaceflight missions.

Gentry, Diana↗

Field-Programmable Gate Array Computer in Structural Analysis: An Initial Exploration

This paper reports on an initial assessment of using a Field-Programmable Gate Array (FPGA) computational device as a new tool for solving structural mechanics problems. A FPGA is an assemblage of binary gates arranged in logical blocks that are interconnected via software in a manner dependent on the algorithm being implemented and can be reprogrammed thousands of times per second. In effect, this creates a computer specialized for the problem that automatically exploits all the potential for parallel computing intrinsic in an algorithm. This inherent parallelism is the most important feature of the FPGA computational environment. It is therefore important that if a problem offers a choice of different solution algorithms, an algorithm of a higher degree of inherent parallelism should be selected. It is found that in structural analysis, an 'analog computer' style of programming, which solves problems by direct simulation of the terms in the governing differential equations, yields a more favorable solution algorithm than current solution methods. This style of programming is facilitated by a 'drag-and-drop' graphic programming language that is supplied with the particular type of FPGA computer reported in this paper. Simple examples in structural dynamics and statics illustrate the solution approach used. The FPGA system also allows linear scalability in computing capability. As the problem grows, the number of FPGA chips can be increased with no loss of computing efficiency due to data flow or algorithmic latency that occurs when a single problem is distributed among many conventional processors that operate in parallel. This initial assessment finds the FPGA hardware and software to be in their infancy in regard to the user conveniences; however, they have enormous potential for shrinking the elapsed time of structural analysis solutions if programmed with algorithms that exhibit inherent parallelism and linear scalability. This potential warrants further development of FPGA-tailored algorithms for structural analysis.

Singleterry, Robert C., Jr.↗

Design and flight test of a differential GPS/inertial navigation system for approach/landing guidance

NASA-Langley has conducted a flight test program evaluating a differential GPS/inertial navigation system's (DGPS/INS) utility as an approach/landing aid. The DGPS/INS airborne and ground components are based on off-the-shelf transport aircraft avionics, namely a global positioning/inertial reference unit (GPIRU) and two GPS sensor units (GPSSUs). Systematic GPS errors are measured by the ground GPSSU and transmitted to the aircraft GPIRU, allowing the errors to be eliminated or greatly reduced in the airborne equipment. Over 120 landings were flown; 36 of these were fully automatic DGPS/INS landings.

Vallot, Lawrence↗

Adaptive computational methods for SSME internal flow analysis

Adaptive finite element methods for the analysis of classes of problems in compressible and incompressible flow of interest in SSME (space shuttle main engine) analysis and design are described. The general objective of the adaptive methods is to improve and to quantify the quality of numerical solutions to the governing partial differential equations of fluid dynamics in two-dimensional cases. There are several different families of adaptive schemes that can be used to improve the quality of solutions in complex flow simulations. Among these are: (1) r-methods (node-redistribution or moving mesh methods) in which a fixed number of nodal points is allowed to migrate to points in the mesh where high error is detected; (2) h-methods, in which the mesh size h is automatically refined to reduce local error; and (3) p-methods, in which the local degree p of the finite element approximation is increased to reduce local error. Two of the three basic techniques have been studied in this project: an r-method for steady Euler equations in two dimensions and a p-method for transient, laminar, viscous incompressible flow. Numerical results are presented. A brief introduction to residual methods of a-posterior error estimation is also given and some pertinent conclusions of the study are listed.

Oden, J. T.↗

Flight demonstration of integrated airport surface automation concepts

A flight demonstration was conducted to address airport surface movement area capacity issues by providing pilots with enhanced situational awareness information. The demonstration showed an integration of several technologies to government and industry representatives. These technologies consisted of an electronic moving map display in the cockpit, a Differential Global Positioning System (DGPS) receiver, a high speed VHF data link, an ASDE-3 radar, and the Airport Movement Area Safety System (AMASS). Aircraft identification was presented to an air traffic controller on AMASS. The onboard electronic map included the display of taxi routes, hold instructions, and clearances, which were sent to the aircraft via data link by the controller. The map also displayed the positions of other traffic and warning information, which were sent to the aircraft automatically from the ASDE-3/AMASS system. This paper describes the flight demonstration in detail, along with preliminary results.

Jones, Denise R.↗

Flight Demonstration of Integrated Airport Surface Automation Concepts

A flight demonstration was conducted to address airport surface movement area capacity issues by providing pilots with enhanced situational awareness information. The demonstration showed an integration of several technologies to government and industry representatives. These technologies consisted of an electronic moving map display in the cockpit, a Differential Global Positioning System (DGPS) receiver, a high speed VHF data link, an ASDE-3 radar, and the Airport Movement Area Safety System (AMASS). Aircraft identification was presented to an air traffic controller on AMASS. The onboard electronic map included the display of taxi routes, hold instructions, and clearances, which were sent to the aircraft via data link by the controller. The map also displayed the positions of other traffic and warning information, which were sent to the aircraft automatically from the ASDE-3/AMASS system. This paper describes the flight demonstration in detail, along with preliminary results.

Jones, Denise R.↗

Use of piloted simulation for studies of fighter departure/spin susceptibility

The NASA-Langley Research Center has incorporated into its stall/spin research program on military airplanes the use of piloted, fixed-base simulation to complement the existing matrix of unique research testing techniques. The piloted simulations of fighter stall/departure flight dynamics are conducted on the Langley Differential Maneuvering Simulator (DMS). The paper reviews the objectives of the simulation research, presents the rationale underlying the simulation methods and procedures used in the evaluation of airplane characteristics, and discuss in detail the evaluation steps used to assess fighter stall/departure characteristics. Simulation results are presented to illustrate the flight-dynamics phenomena dealt with. The use of the piloted simulation methods and procedures described has been found very effective in identifying stability and control problem areas and in developing automatic control concepts to alleviate many of these problems. A good level of correlation between simulated flight dynamics and flight-test results has been obtained over the many fighter configurations studied in the simulator.

Gilbert, W. P.↗

Point to Point Control of the Hydrogen Mixer

We continue to build on the theoretical modelling of a liquid hydrogen (LH2) and gas hydrogen (GH2) mixer subsystem. The mixer described in this work is responsible for combining high pressure LH2 and GH2 to produce a hydrogen flow that meets certain thermodynamic properties. The desired properties are maintained by precise control of the LH2 and GH2 flows. The mixer is modelled as a general multi-flow lumped volume for single constituent fluids using density and internal energy as states. The set of nonlinear differential equations is modelled in the SIMULINK environment including a table look-up feature of the fluid thermodynamic properties. A small signal (linear) model is developed based on the nonlinear model and simulated as well. Pulse disturbances are introduced to the valve positions and the quality of the linear model is ascertained by comparing its behavior against the nonlinear model simulations. Valve control strategies that simulate an operator-in-the-loop scenario are then explored demonstrating the need for automatic feedback control. Finally, optimal single-output and multi-output Proportional/Integral controllers are designed based on the linear model and applied to the nonlinear model with excellent results to track simultaneous, constant setpoint changes in desired exit flow, exit temperature, and mixer pressure, as well as to reject unmeasureable but bounded additive step perturbations in the valve positions.

Barbieri, Enrique↗

Small Signal Point-to-Point Tracking of a Propellant Mixer

This paper addresses some theoretical modelling and control issues for a mixing chamber used in rocket engine testing at NASA Stennis Space Center. The mixer is responsible for combining high pressure LH2 and GH2 to produce a hydrogen flow that meets certain thermodynamic properties before it is fed into a test article. The desired properties are maintained by precise control of the LH2 and GH2 flows. The mixer is modelled as a general multi-flow lumped volume for single constituent fluids using density and internal energy as states. The set of nonlinear differential equations is modelled in the SIMULINK environment including a table look-up feature of the fluid thermodynamic properties. a small-signal (linear) model is developed based on the nonlinear model and simulated as well. Pulse disturbances are introduced to the valve positions and the quality of the linear model is ascertained by comparing its behavior against the nonlinear model simulations. Valve control strategies that simulate an operator-in-the-loop scenario are then explored demonstrating the need for automatic feedback control. Finally, classical optimal single-output and multi-output Proportional/Integral controllers are designed based on the linear model and applied to the nonlinear model with excellent results to track simultaneous, constant setpoint changes in desired exit flow, exit temperature, and mixer pressure, as well as to reject unmeasurable but bounded additive step perturbations in the valve positions.

Barbieri, Enrique↗