Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational algorithms”

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,153 records · Page 64

Clutter Distributions for Tomographic Image Standardization in Ground-Penetrating Radar

Multistatic ground-penetrating radar (GPR) signals can be imaged tomographically to produce 3-D distributions of image intensities. In the absence of objects of interest, these intensities can be considered to be estimates of clutter. These clutter intensities spatially vary over several orders of magnitude and vary across different arrays, which makes a direct comparison of these raw intensities difficult. However, by gathering statistics on these intensities and their spatial variation, a variety of metrics can be determined. In this study, the clutter distribution is found to fit better to a two-parameter Weibull distribution than Gaussian or log-normal distributions. Based on the spatial variation of the two Weibull parameters, scale and shape, more information may be gleaned from these data. How well the GPR array is illuminating various parts of the ground, in depth and cross track, may be determined from the spatial variation of the Weibull scale parameter, which may in turn be used to estimate an effective attenuation coefficient in the soil. The transition in depth from clutter- to noise-limited conditions (which is one possible definition of GPR penetration depth) can be estimated from the spatial variation of the Weibull shape parameter. Lastly, the underlying clutter distributions also provide an opportunity to standardize image intensities to determine when a statistically significant deviation from background (clutter) has occurred, which is convenient for buried threat detection algorithm development that needs to be robust across multiple different arrays.

42 ENGINEERING↗

sourcePy

Pollutant source identification techniques (of which there are many variations) are either locked behind researchers writing their own code for each use case or GUI platforms that are easy to use but inflexible and opaque. The Python package sourcePy brings together many of the pollutant source identification algorithms, giving the user full control out of the box. It aims to create a platform for source identification experiments where the full analysis from beginning to end can be done in Python, with a level of specificity in design that isn't available in the GUI options. sourcePy provides users with a few key features: -A Python interface with HYSPLIT, which can be used to generate trajectories and concentration plumes -Several Python classes which standardize the preparation and processing of data related to source identification experiments -Example scripts and notebooks that allow even new python users to get started with their own experiments quickly -Visualization methods

Arseneau, Isaac [Oak Ridge National Laboratory (OR↗

Machine learning identifies novel signatures of antifungal drug resistance in Saccharomycotina yeasts

Antifungal drug resistance is a major challenge in fungal infection management. Numerous genomic changes are known to contribute to acquired drug resistance in clinical isolates of specific pathogens, but whether they broadly explain natural resistance across entire lineages is unknown. We leveraged genomic, ecological, and phenotypic trait data from naturally sampled strains from nearly all known species in subphylum Saccharomycotina to examine the evolution of resistance to eight antifungal drugs. The phylogenetic distribution of drug resistance varied by drug; fluconazole resistance was widespread, while 5-fluorocytosine resistance was rare, except in Lipomycetales. A random forest algorithm trained on genomic data predicted drug-resistant yeasts with 54–75% accuracy. Fluconazole resistance was consistently predicted with the highest accuracy (75.2%). Furthermore, fluconazole resistance prediction accuracy was similar between models trained on genome-wide variation in the presence and number of InterPro protein annotations across Saccharomycotina (75.2%) and those trained on amino acid sequence alignment data of Erg11, a protein known to be involved in fluconazole resistance (74.3-74.9%). Interestingly, the top Erg11 residues for predicting fluconazole resistance across Saccharomycotina do not overlap with, are not spatially close to, and are less conserved than those previously linked to resistance in clinical isolates of Candida albicans. In silico deep mutational scanning of the C. albicans Erg11 protein reveals that amino acid variants implicated in clinical cases of resistance are almost universally destabilizing while variants in our most informative residues are energetically more neutral, explaining why the latter are much more common than the former in natural populations. Importantly, previous experimental analyses of C. albicans Erg11 have shown that amino acid variation in our most informative residues, despite having never been directly implicated in clinical cases, can directly contribute to resistance. Our results suggest that studies of natural resistance in yeast species never encountered in the clinic will yield a fuller understanding of antifungal drug resistance.

Harrison, Marie-Claire [Vanderbilt Univ., Nashvill↗

Prediction and measurement of turbulent aerodynamic trailing edge flows

A viscous-inviscid interaction algorithm is developed for prediction of two-dimensional mean and fluctuating velocity distributions in the wake immediately downstream of an airfoil trailing edge. A composite pressure field is defined, and a Poisson equation solved for transverse pressure variations. A parabolized form of the time-averaged steady Navier-Stokes equations are solved in conjunction with a viscous-augmented two-dimensional inviscid potential flow analysis. A tensor constitutive equation is employed to predict Reynolds stress distributions from solutions of a turbulence kinetic energy two equation closure model. Numerical predictions compared favorably with detailed experimental data for mean and fluctuating velocities, and Reynolds shear stress distributions, in the trailing edge region of a NACA 63-012 airfoil.

Baker, A. J.↗

Monte Carlo simulations of precise timekeeping in the Milstar communication satellite system

The Milstar communications satellite system will provide secure antijam communication capabilities for DOD operations into the next century. In order to accomplish this task, the Milstar system will employ precise timekeeping on its satellites and at its ground control stations. The constellation will consist of four satellites in geosynchronous orbit, each carrying a set of four rubidium (Rb) atomic clocks. Several times a day, during normal operation, the Mission Control Element (MCE) will collect timing information from the constellation, and after several days use this information to update the time and frequency of the satellite clocks. The MCE will maintain precise time with a cesium (Cs) atomic clock, synchronized to UTC(USNO) via a GPS receiver. We have developed a Monte Carlo simulation of Milstar's space segment timekeeping. The simulation includes the effects of: uplink/downlink time transfer noise; satellite crosslink time transfer noise; satellite diurnal temperature variations; satellite and ground station atomic clock noise; and also quantization limits regarding satellite time and frequency corrections. The Monte Carlo simulation capability has proven to be an invaluable tool in assessing the performance characteristics of various timekeeping algorithms proposed for Milstar, and also in highlighting the timekeeping capabilities of the system. Here, we provide a brief overview of the basic Milstar timekeeping architecture as it is presently envisioned. We then describe the Monte Carlo simulation of space segment timekeeping, and provide examples of the simulation's efficacy in resolving timekeeping issues.

Camparo, James C.↗

Efficient Methods to Assimilate Satellite Retrievals Based on Information Content: Suboptimal Retrieval Assimilation - Part 2

One of the outstanding problems in data assimilation has been and continues to be how best to utilize satellite data while balancing the tradeoff between accuracy and computational cost. A number of weather prediction centers have recently achieved remarkable success in improving their forecast skill by changing the method by which satellite data are assimilated into the forecast model from the traditional approach of assimilating retrievals to the direct assimilation of radiances in a variational framework. The operational implementation of such a substantial change in methodology involves a great number of technical details, e.g., pertaining to quality control procedures, systematic error correction techniques, and tuning of the statistical parameters in the analysis algorithm. Although there are clear theoretical advantages to the direct radiance assimilation approach, it is not obvious at all to what extent the improvements that have been obtained so far can be attributed to the change in methodology, or to various technical aspects of the implementation. The issue is of interest because retrieval assimilation retains many practical and logistical advantages which may become even more significant in the near future when increasingly high-volume data sources become available. The central question we address here is: how much improvement can we expect from assimilating radiances rather than retrievals, all other things being equal? We compare the two approaches in a simplified one-dimensional theoretical framework, in which problems related to quality control and systematic error correction are conveniently absent. By assuming a perfect radiative transfer model and perfect knowledge of radiance and background error covariances, we are able to formulate a nonlinear local error analysis for each assimilation method. Direct radiance assimilation is optimal in this idealized context, while the traditional method of assimilating retrievals is suboptimal because it ignores the cross-covariances between background errors and retrieval errors. We show that interactive retrieval assimilation (where the same background used for assimilation is also used in the retrieval step) is equivalent to direct assimilation of radiances with suboptimal analysis weights. We illustrate and extend these theoretical arguments with several one-dimensional assimilation experiments, where we estimate vertical atmospheric profiles using simulated data from both the High-resolution InfraRed Sounder 2 (HIRS2) and the future Atmospheric InfraRed Sounder (AIRS).

Joiner, J.↗

Aeroelastic Stability of a Soft-Inplane Gimballed Tiltrotor Model in Hover

Soft-inplane rotor systems can significantly reduce the inplane rotor loads generated during the maneuvers of large tiltrotors, thereby reducing the strength requirements and the associated structural weight of the hub. Soft-inplane rotor systems, however, are subject to instabilities associated with ground resonance, and for tiltrotors this instability has increased complexity as compared to a conventional helicopter. Researchers at Langley Research Center and Bell Helicopter-Textron, Inc. have completed an initial study of a soft-inplane gimballed tiltrotor model subject to ground resonance conditions in hover. Parametric variations of the rotor collective pitch and blade root damping., and their associated effects on the model stability were examined. Also considered in the study was the effectiveness of an active swashplate and a generalized predictive control (GPC) algorithm for stability augmentation of the ground resonance conditions. Results of this study show that the ground resonance behavior of a gimballed soft-inplane tiltrotor can be significantly different from that of a classical soft-inplane helicopter rotor. The GPC-based active swashplate was successfully implemented, and served to significantly augment damping of the critical modes to an acceptable value.

Nixon, Mark W.↗

The New IERS Special Bureau for Loading (SBL)

Currently, the establishment of the International Earth Rotation Service (IERS) Special Bureau for Loading (SBL) is in progress as part of the IERS Global Geophysical Fluids Center (GGFC). The main purpose of the SBL is to provide reliable, consistent model predictions of loading signals that have been thoroughly tested and validated. The products will describe at least the surface deformation, gravity signal and geo-center variations due to the various surface loading processes in reference frames relevant for direct comparison with existing geodetic observing techniques. To achieve these goals, major scientific advances are required with respect to the Earth model, the theory and algorithms used to model deformations of the Earth as well as improvements in the observational data related to surface loading.

vanDam, Tonie↗

Effects of the Antenna Aperture on Remote Sensing of Sea Surface Salinity at L-Band

Remote sensing of sea surface salinity with sufficient accuracy to meet the needs of global oceanography is a challenging task. The global variability of the salinity signal in the open ocean is only a few Kelvin even at L-band and an accuracy on the order of 0.1K is desired to study the influence of salinity on ocean circulation and energy exchange with the atmosphere. On the other hand, resolution is not an issue for understanding the dynamics of the open ocean where scales of hundreds of km are not uncommon. This permits remote sensing with large antenna footprints and spatial averaging to reduce noise. However, antennas with large footprints introduce other problems. For example, the angle of incidence and hence the brightness temperature varies over the footprint. Similarly, the polarization of brightness temperature relative to the antenna ports changes. Studies have been conducted using antenna patterns representative of the antenna that will be flown on the Aquarius mission to examine these effects. Aquarius is a pushbroom style radiometer with three beams looking across track away from the sun. The beams are at incidences angles (at the spacecraft) of about 26.5, 34 and 40 degrees each with a half-power beam width of about 5.8 degrees. It is shown that the measured brightness temperature is biased relative to the value at boresight because of changes across the field of view. The bias can be as much as 4K and positive or negative depending on polarization. Polarization mixing because of the variations of the local plane of incidence across the footprint also occur and can result in biased polarimetric measurements. A bias in the third Stokes parameter of as much as 0.4K is possible. Such effects may affect algorithms that use the third Stokes parameter to correct for Faraday rotation. Another issue associated with the antenna is sun glint. This is an issue determined by surface roughness and antenna sidelobes. Examples will be given for the random component (glint) for the case of the Aquarius antenna beams. Fortunately, the Aquarius beams mostly look to the dark side of the day-night termination, but during some portions of the year they will see sun-lighted ocean. In this case, glint could be an issue for the inner-most beam.

Dinnat, Emmanuel P.↗

Hyperspectral Remote Sensing of Atmosphere and Surface Properties

Atmospheric Infrared Sounder (AIRS), Infrared Atmospheric Sounding Interferometer (IASI), and Cross-track Infrared Sounder (CrIS) are all hyper-spectral satellite sensors with thousands of spectral channels. Top of atmospheric radiance spectra measured by these sensors contain high information content on atmospheric, cloud, and surface properties. Exploring high information content contained in these high spectral resolution spectra is a challenging task due to computation e ort involved in modeling thousands of spectral channels. Usually, only very small fractions (4{10 percent) of the available channels are included in physical retrieval systems or numerical weather forecast (NWP) satellite data assimilations. We will describe a method of simultaneously retrieving atmospheric temperature, moisture, cloud, and surface properties using all available spectral channels without sacrificing computational speed. The essence of the method is to convert channel radiance spectra into super-channels by an Empirical Orthogonal Function (EOF) transformation. Because the EOFs are orthogonal to each other, about 100 super-channels are adequate to capture the information content of the radiance spectra. A Principal Component-based Radiative Transfer Model (PCRTM) developed at NASA Langley Research Center is used to calculate both the super-channel magnitudes and derivatives with respect to atmospheric profiles and other properties. There is no need to perform EOF transformations to convert super channels back to spectral space at each iteration step for a one-dimensional variational retrieval or a NWP data assimilation system. The PCRTM forward model is also capable of calculating radiative contributions due to multiple-layer clouds. The multiple scattering effects of the clouds are efficiently parameterized. A physical retrieval algorithm then performs an inversion of atmospheric, cloud, and surface properties in super channel domain directly therefore both reducing the computational need and preserving the information content of the IASI measurements. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance matrices and a priori information as constraints. One advantage of this approach is that it uses all information content from the hyper-spectral data so that the retrieval is less sensitive to instrument noise and eliminates the need for selecting a sub-set of the channels.

Liu, Xu↗

Autonomous Guidance Navigation and Control for Agile Quadrotors Using Polynomial Trajectory Planning and L1 Adaptive Control

We address the challenge to allow efficient autonomous flight in real world environments, both indoor and outdoor. We use a straight-line SE-SCP (Spherical Expansion and Sequential Convex Programming) [algorithm] to find an initial route through the environment and minimum snap trajectory generation using piecewise polynomials. Then, we implement an adaptive robust control able to address some robustness issues for quadrotors in outdoor flight, such as mass variation and wind disturbances. Coupling these techniques we allow high-speed and aggressive autonomous flight through obstacle-dense indoor environments, as well as address outdoor disturbances.

Landolfi, Mattia↗

Template Matching Used for Small Body Optical Navigation with Poorly Detailed Objects

Object and template matching becomes difficult when an image lacks detail. This is particularly worrisome when typical matching techniques, cross-correlation, log-polar mapping, and key point matching fail. Work herein describes a formulation that identifies objects of interest, estimates the affine transformation between a template object and scene using Principal Component Analysis (PCA), and provides a fit value for the objects and template incorporating Hu's Moments. The algorithm presented is tested on synthetic images and images obtained from the OSIRIS-REx mission while the spacecraft was approaching its target, Bennu. Results for the current formulation show that, with the presence of large-scale variations and rotation, the fitting scheme performs well when compared with other techniques.

Lyzhoft, Joshua R.↗

Eighteen Years of Aqua MODIS On-orbit Operation, Calibration, and Performance

Aqua MODIS has successfully operated for more than 18 years since its launch in May 2002and has generated numerous science products in support of studies of the Earth’s system and its changes via a set of geophysical and environmental parameters. On-orbit calibration and characterization activities have played a vital role in maintaining the quality of MODIS data products. In addition to data collected from sensor on-board calibrators (OBC), near-monthly lunar observations and select earth view targets have been used to monitor and characterize on-orbit changes in sensor responses and to derive and update the calibration look-up tables. In this paper, we provide an overview of Aqua MODIS on-orbit operations, calibration activities and approaches, and algorithm improvements and also illustrate the sensor on-orbit performance using examples derived from various calibration sources and targets. We will focus on issues identified from instrument operations and calibrations, such as solar diffuser degradation, electronic crosstalk, variations in the cold focal plane temperatures, and changes in response versus scan-angle. Also discussed in this paper are remaining challenges and future improvements.

Aqua↗

Successive Convexification for Powered Descent Guidance with Time-Varying Mass Properties

This paper extends the 6 degree-of-freedom successive convexification (SCvx) algorithm discussed in literature to account for time-varying mass properties in the representation of the attitude dynamics. Simulation results of a notional large-size lander performing the last minutes of Lunar powered descent with significant variations in the center-of-mass and moments of inertia over the duration of the trajectory are shown. Differences in the final optimal trajectories generated via SCvx are presented between the constant vs. variable mass property formulations for an in-plane and out-of-plane example.

Alex Hayes↗

Successive Convexification for Powered Descent Guidance with Time-Varying Mass Properties

This paper extends the 6 degree-of-freedom successive convexification (SCvx) algorithm discussed in literature to account for time-varying mass properties in the representation of the attitude dynamics. Simulation results of a notional large-size lander performing the last minutes of Lunar powered descent with significant variations in the center-of-mass and moments of inertia over the duration of the trajectory are shown. Differences in the final optimal trajectories generated via SCvx are presented between the constant vs. variable mass property formulations for an in-plane and out-of-plane example.

Alex Hayes↗

The Localized Active Space Method with Unitary Selective Coupled Cluster

Here, we introduce a hybrid quantum-classical algorithm, the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. Derived from the localized active space unitary coupled cluster (LAS-UCCSD) method, LAS-USCCSD first performs a classical LASSCF calculation, then selectively identifies the most important parameters (cluster amplitudes used to build the multireference UCC ansatz) for restoring interfragment interaction energy using this reduced set of parameters with the variational quantum eigensolver method. We benchmark LAS-USCCSD against LAS-UCCSD by calculating the total energies of (H 2 ) 2 , (H 2 ) 4 , and trans-butadiene, and the magnetic coupling constant for a bimetallic compound [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ . For these systems, we find that LAS-USCCSD reduces the number of required parameters and thus the circuit depth by at least 1 order of magnitude, an aspect which is important for the practical implementation of multireference hybrid quantum-classical algorithms like LAS-UCCSD on near-term quantum computers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interannual Variability of Tropical Rainfall as Seen From TRMM

Considerable uncertainty surrounds the issue of whether precipitation over the tropical oceans (30deg N/S) systematically changes with interannual sea-surface temperature (SST) anomalies that accompany El Nino (warm) and La Nina (cold) events. Although it is well documented that El Nino-Southern Oscillation (ENSO) events with marked SST changes over the tropical oceans produce significant regional changes in precipitation, water vapor, and radiative fluxes in the tropics, we still cannot yet adequately quantify the associated net integrated changes to water and heat balance over the entire tropical oceanic or land sectors. Resolving this uncertainty is important since precipitation and latent heat release variations over land and ocean sectors are key components of the tropical heat balance in its most aggregated form. Rainfall estimates from the Version 5 Tropical Rainfall Measuring Mission (TRMM) Precipitation Radar (PR) averaged over the tropical oceans have not solved this issue and, in fact, show marked differences with estimates from two TRMM Microwave Imager (TMI) passive microwave algorithms. In this paper we will focus on findings that uncertainties in microphysical assumptions necessitated by the single-frequency PR measurement pose difficulties for detecting climate-related precipitation signals. Recent work has shown that path-integrated attenuation derived from the effects of precipitation on the radar return from the ocean surface exhibits interannual variability that agrees closely with the TMI time series, yet the PR rainfall interannual variability (and attenuation derived predominantly from reflectivity) differs even in sign. We will explore these apparent inconsistencies and examine changes in new TRMM Version 6 retrievals. To place these results in a tropical water balance perspective we also examine interannual variations in evaporation over the tropical oceans made from TRMM and SSM/I (Special Sensor Microwave Imager) measurements of surface winds and humidity. Evaporation estimates from reanalysis and several global model experiments will also be compared to the TRMM findings and evaluated for consistency. The ability to detect regional shifts in freshwater flux over the oceans (equivalently, integrated moisture convergence) and moisture transport will be discussed.

Robertson, Franklin R.↗

Scalable semidefinite programming approach to variational embedding for quantum many-body problems

In quantum embedding theories, a quantum many-body system is divided into localized clusters of sites which are treated with an accurate ‘high-level’ theory and glued together self-consistently by a less accurate ‘low-level’ theory at the global scale. The recently introduced variational embedding approach for quantum many-body problems combines the insights of semidefinite relaxation and quantum embedding theory to provide a lower bound on the ground-state energy that improves as the cluster size is increased. The variational embedding method is formulated as a semidefinite program (SDP), which can suffer from poor computational scaling when treated with black-box solvers. Here, we exploit the interpretation of this SDP as an embedding method to develop an algorithm which alternates parallelizable local updates of the high-level quantities with updates that enforce the low-level global constraints. Moreover, we show how translation invariance in lattice systems can be exploited to reduce the complexity of projecting a key matrix to the positive semidefinite cone.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗