Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variance estimation”

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

Recipes for when physics fails: recovering robust learning of physics informed neural networks

Abstract Physics-informed neural networks (PINNs) have been shown to be effective in solving partial differential equations by capturing the physics induced constraints as a part of the training loss function. This paper shows that a PINN can be sensitive to errors in training data and overfit itself in dynamically propagating these errors over the domain of the solution of the PDE. It also shows how physical regularizations based on continuity criteria and conservation laws fail to address this issue and rather introduce problems of their own causing the deep network to converge to a physics-obeying local minimum instead of the global minimum. We introduce Gaussian process (GP) based smoothing that recovers the performance of a PINN and promises a robust architecture against noise/errors in measurements. Additionally, we illustrate an inexpensive method of quantifying the evolution of uncertainty based on the variance estimation of GPs on boundary data. Robust PINN performance is also shown to be achievable by choice of sparse sets of inducing points based on sparsely induced GPs. We demonstrate the performance of our proposed methods and compare the results from existing benchmark models in literature for time-dependent Schrödinger and Burgers’ equations.

97 MATHEMATICS AND COMPUTING↗

New Probe of Cosmic Birefringence Using Galaxy Polarization and Shapes

We propose a novel statistical method to measure cosmic birefringence and demonstrate its power in probing parity violation due to axions. Exploiting an empirical correlation between the integrated radio polarization direction of a spiral galaxy and its apparent shape, we devise an unbiased minimum-variance estimator for the rotation angle, which should achieve an uncertainty of 5°–15° per galaxy. In conclusion, large galaxy samples from the forthcoming SKA continuum surveys, together with optical shape catalogs, promise a comparable or even lower noise power spectrum for the rotation angle than in the CMB Stage-IV (CMB-S4) experiment, with different systematics.

Axion-like particles↗

Key environmental and production factors for understanding variation in switchgrass chemical attributes

Abstract Switchgrass ( Panicum virgatum L.) is a promising feedstock for bioenergy and bioproducts; however, its inherent variability in chemical attributes creates challenges for uniform conversion efficiencies and product quality. It is necessary to understand the range of variation and factors (i.e., field management, environmental) influencing chemical attributes for process improvement and risk assessment. The objectives of this study were to (1) examine the impact of nitrogen fertilizer application rate, year, and location on switchgrass chemical attributes, (2) examine the relationships among chemical attributes, weather and soil data, and (3) develop models to predict chemical attributes using environmental factors. Switchgrass samples from a field study spanning four locations including upland cultivars, one location including a lowland cultivar, and between three and six harvest years were assessed for glucan, xylan, lignin, volatiles, carbon, nitrogen, and ash concentrations. Using variance estimation, location/cultivar, nitrogen application rate, and year explained 65%–96% of the variation for switchgrass chemical attributes. Location/cultivar × year interaction was a significant factor for all chemical attributes indicating environmental‐based influences. Nitrogen rate was less influential. Production variables and environmental conditions occurring during the switchgrass field trials were used to successfully predict chemical attributes using linear regression models. Upland switchgrass results highlight the complexity in plant responses to growing conditions because all production and environmental variables had strong relationships with one or more chemical attributes. Lowland switchgrass was limited to observations of year‐to‐year environmental variability and nitrogen application rate. All explanatory variable categories were important for lowland switchgrass models but stand age and precipitation relationships were particularly strong. The relationships found in this study can be used to understand spatial and temporal variation in switchgrass chemical attributes. The ability to predict chemical attributes critical for conversion processes in a geospatial/temporal manner would provide state‐of‐the‐art knowledge for risk assessment in the bioenergy and bioproducts industry.

09 BIOMASS FUELS↗

Neural Network-Based Electric Vehicle Range Prediction for Smart Charging Optimization

Range prediction is a standard feature in most modern road vehicles, allowing drivers to make informed decisions about when to refuel. Most vehicles make range predictions through data- or model-driven means, monitoring the average fuel consumption rate or using a tuned vehicle model to predict fuel consumption. The uncertainty of future driving conditions makes the range prediction problem challenging, particularly for less pervasive battery electric vehicles (BEV). Most contemporary machine learning-based methods attempt to forecast the battery SOC discharge profile to predict vehicle range. In this work, we propose a novel approach using two recurrent neural networks (RNNs) to predict the remaining range of BEVs and the minimum charge required to safely complete a trip. Each RNN has two outputs that can be used for statistical analysis to account for uncertainties; the first loss function leads to mean and variance estimation (MVE), while the second results in bounded interval estimation (BIE). These outputs of the proposed RNNs are then used to predict the probability of a vehicle completing a given trip without charging, or if charging is needed, the remaining range and minimum charging required to finish the trip with high probability. Training data was generated using a low-order physics model to estimate vehicle energy consumption from historical drive cycle data collected from medium-duty last-mile delivery vehicles. Here, the proposed method demonstrated high accuracy in the presence of day-to-day route variability, with the root-mean-square error (RMSE) below 6% for both RNN models.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Human Lacrimal Production Rate and Wetted Length of Modified Schirmer's Tear Test Strips

Purpose: To assess and compare the wetting kinetics of sheathed and unsheathed Schirmer's tear test (STT) strips, and to determine the repeatability of 5-minute wetted length (WL) and basal tear production rate (BTPR). Methods:Seventeen subjects underwent two sheathed and unsheathed STTs each for both eyes on four visits on separate days. After administration of topical anesthetic, WLs were measured every 30 seconds for 5 minutes, and BTPRs were calculated for sheathed strips. Limits of agreement (LoA), difference-versus-mean plots (DVM), and the coefficient of repeatability (CR) assessed WL and BTPR repeatabilities. Variance estimates were used to calculate sample sizes for future study. Results: For the unsheathed STT, the mean (SD) difference in WLs between visits was 0.74 (5.05) mm, LoA were [-9.17, 10.64], and CR was 9.17 mm; for the sheathed STT, the mean (SD) intervisit difference was 0.16 (5.94) mm, LoA were [-11.49, 11.8], and CR was 10.53 mm. Eight of 48 sheathed STTs and 20 of 44 unsheathed STTs showed constant WL for the final 90 seconds of the test. The mean (SD) difference between repeated visits for BTPR was approximately 0.0 μL/min, LoA were [-1.82, 1.82], and CR was 1.91 μL/min. Conclusions: Repeatability of sheathed and unsheathed 5-minute WL and BTPR is inadequate for measuring within-subject changes, but is sufficient for group studies with moderate sample sizes. Constant WL for the fiant WL with the unsheathed STT can be explained by balancing evaporation and BTPR. Translational Relevance: Repeatability of the modified STT is evaluated clinically to establish quantitative BTPRs rather than inference from a strip WL.

59 BASIC BIOLOGICAL SCIENCES↗

Scalable Gaussian Processes, GPyTorch Application Benchmarking, and Targeted Adaptive Design (TAD) on ThetaGPU

We aim at showcasing the scalability of Gaussian Process (GP). The naive GP implementation scales cubically with data size, which can be prohibitive, so GP has not heretofore been considered suitable for very large-scale problem settings. We take advantage of GPyTorch, a library for scalable GPs built on top of PyTorch that incorporates GPU acceleration. With GPyTorch, one can achieve nearly linear scaling with structured kernel interpolation (SKI) and constant-time predictive covariances computation with LanczOs Variance Estimates (LOVE) while preserving accuracy. We also take advantage of the computational power of ThetaGPU, a supercomputer of Argonne Leadership Computing Facility (ALCF). In addition, we implement a scalable, GPU-ready version of Targeted Adaptive Design (TAD), a GP-based data-driven algorithm that efficiently searches the control space of an advanced manufacturing experiment for settings capable of producing a required design within a specified tolerance, despite the poorly known mapping from control settings to design. We finally show our benchmarking for GPyTorch and TAD performance on CPU vs. ThetaGPU and discuss the results and implications.

97 MATHEMATICS AND COMPUTING↗

The Sixth Data Release of the Radial Velocity Experiment (RAVE). I. Survey Description, Spectra, and Radial Velocities

The Radial Velocity Experiment (Rave) is a magnitude-limited (9 < I < 12) spectroscopic survey of Galactic stars randomly selected in Earth’s southern hemisphere. The Rave medium-resolution spectra (R ~ 7500) cover the Ca-triplet region (8410–8795 Å). The sixth and final data release (DR6) is based on 518,387 observations of 451,783 unique stars. Rave observations were taken between 2003 April 12 and 2013 April 4. Here we present the genesis, setup, and data reduction of Rave as well as wavelength-calibrated and flux-normalized spectra and error spectra for all observations in Rave DR6. Furthermore, we present derived spectral classification and radial velocities for the Rave targets, complemented by cross-matches with Gaia DR2 and other relevant catalogs. A comparison between internal error estimates, variances derived from stars with more than one observing epoch, and a comparison with radial velocities of Gaia DR2 reveals consistently that 68% of the objects have a velocity accuracy better than 1.4 km s{sup –1}, while 95% of the objects have radial velocities better than 4.0 km s{sup –1}. Stellar atmospheric parameters, abundances and distances are presented in a subsequent publication. The data can be accessed via the Rave website (http://rave-survey.org) or the Vizier database.

79 ASTRONOMY AND ASTROPHYSICS↗

HYPHY: Deep Generative Conditional Posterior Mapping of Hydrodynamical Physics

Generating large-volume hydrodynamical simulations for cosmological observables is a computationally demanding task necessary for next-generation observations. In this work, we construct a novel fully convolutional variational autoencoder (VAE) to synthesize hydrodynamic fields conditioned on dark matter fields from N -body simulations. After training the model on a single hydrodynamical simulation, we are able to probabilistically map new dark-matter-only simulations to corresponding full hydrodynamical outputs. By sampling over the latent space of our VAE, we can generate posterior samples and study the variance of the mapping. We find that our reconstructed field provides an accurate representation of the target hydrodynamical fields as well as reasonable variance estimates. This approach has promise for the rapid generation of mocks as well as for implementation in a full inverse model of observed data.

79 ASTRONOMY AND ASTROPHYSICS↗

Guidance analysis.

Unmanned spacecraft guidance policy and orbit determination, using linear perturbation theory and minimum variance estimation

SPACECRAFT GUIDANCE↗

Map synchronization in optical communication systems

The time synchronization problem in an optical communication system is approached as a problem of estimating the arrival time (delay variable) of a known transmitted field. Maximum aposteriori (MAP) estimation procedures are used to generate optimal estimators, with emphasis placed on their interpretation as a practical system device, Estimation variances are used to aid in the design of the transmitter signals for best synchronization. Extension is made to systems that perform separate acquisition and tracking operations during synchronization. The closely allied problem of maintaining timing during pulse position modulation is also considered. The results have obvious application to optical radar and ranging systems, as well as the time synchronization problem.

Gagliardi, R. M.↗

Intermediate energy proton-deuteron elastic scattering

A fully symmetrized multiple scattering series is considered for the description of proton-deuteron elastic scattering. An off-shell continuation of the experimentally known twobody amplitudes that retains the exchange symmeteries required for the calculation is presented. The one boson exchange terms of the two body amplitudes are evaluated exactly in this off-shell prescription. The first two terms of the multiple scattering series are calculated explicitly whereas multiple scattering effects are obtained as minimum variance estimates from the 146-MeV data of Postma and Wilson. The multiple scattering corrections indeed consist of low order partial waves as suggested by Sloan based on model studies with separable interactions. The Hamada-Johnston wave function is shown consistent with the data for internucleon distances greater than about 0.84 fm.

Wilson, J. W.↗

ALTKAL: An optimum linear filter for GEOS-3 altimeter data

ALTKAL is a computer program designed to smooth sea surface height data obtained from the GEOS 3 altimeter, and to produce minimum variance estimates of sea surface height and sea surface slopes, along with their standard derivations. The program operates by processing the data through a Kalman filter in both the forward and backward directions, and optimally combining the results. The sea surface height signal is considered to have a geoid signal, modeled by a third order Gauss-Markov process, corrupted by additive white noise. The governing parameters for the signal and noise processes are the signal correlation length and the signal-to-noise ratio. Mathematical derivations of the filtering and smoothing algorithms are presented. The smoother characteristics are illustrated by giving the frequency response, the data weighting sequence and the transfer function of a realistic steady-state smoother example. Based on nominal estimates for geoidal undulation amplitude and correlation length, standard deviations for the estimated sea surface height and slope are 12 cm and 3 arc seconds, respectively.

Fang, B. T.↗

Parameter identification of linear discrete stochastic systems with time delays

An identification algorithm that uses the maximum likelihood technique to identify the unknown time delays, plant parameters, and noise covariances of linear discrete stochastic systems is presented. Cases of additive white noise and colored measurement noises are considered. The likelihood function is evaluated using either a minimum-variance (Kalman) filter or a minimal-order observer. The Kalman filter is used in the identification algorithm to provide minimum-variance estimates. The minimal-order observer is a lower-dimensional and computationally simpler filter, and is advantageous especially for systems with long delays. It provides a less optimal solution to the minimum-mean-square state estimation problem. The colored-noise observer algorithm has the disadvantage of having to compute an extra error covariance matrix of lower order.

Wong, E. C.↗

Evaluation of corn/soybeans separability using Thematic Mapper and Thematic Mapper Simulator data

Multitemporal Thematic Mapper, Thematic Mapper Simulator, and detailed ground truth data were collected for a 9- by 11-km sample segment in Webster County, IA, in the summer of 1982. Three dates were acquired each with Thematic Mapper Simulator (June 7, June 23, and July 31) and Thematic Mapper (August 2, September 3, and October 21). The Thematic Mapper Simulator data were converted to equivalent TM count values using TM and TMS calibration data and model based estimates of atmospheric effects. The July 31, TMS image was compared to the August 2, TM image to verify the conversion process. A quantitative measure of proportion estimation variance (Fisher information) was used to evaluate the corn/soybeans separability for each TM band as a function of time during the growing season. The additional bands in the middle infrared allowed corn and soybeans to be separated much earlier than was possible with the visible and near-infrared bands alone. Using the TM and TMS data, temporal profiles of the TM principal components were developed. The greenness and brightness exhibited behavior similar to MSS greenness and brightness for corn and soybeans.

Pitts, D. E.↗

Improved importance sampling technique for efficient simulation of digital communication systems

A new, improved importance sampling (IIS) approach to simulation is considered. Some basic concepts of IS are introduced, and detailed evolutions of simulation estimation variances for Monte Carlo (MC) and IS simulations are given. The general results obtained from these evolutions are applied to the specific previously known conventional importance sampling (CIS) technique and the new IIS technique. The derivation for a linear system with no signal random memory is considered in some detail. For the CIS technique, the optimum input scaling parameter is found, while for the IIS technique, the optimum translation parameter is found. The results are generalized to a linear system with memory and signals. Specific numerical and simulation results are given which show the advantages of CIS over MC and IIS over CIS for simulations of digital communications systems.

Lu, Dingqing↗