Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Covariance Generation”

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 127 records · Page 7

True covariance simulation of the EUVE update filter

This paper presents a covariance analysis of the performance and sensitivity of the attitude determination Extended Kalman Filter (EKF) used by the On Board Computer (OBC) of the Extreme Ultra Violet Explorer (EUVE) spacecraft. The linearized dynamics and measurement equations of the error states are used in formulating the 'truth model' describing the real behavior of the systems involved. The 'design model' used by the OBC EKF is then obtained by reducing the order of the truth model. The covariance matrix of the EKF which uses the reduced order model is not the correct covariance of the EKF estimation error. A 'true covariance analysis' has to be carried out in order to evaluate the correct accuracy of the OBC generated estimates. The results of such analysis are presented which indicate both the performance and the sensitivity of the OBC EKF.

Bar-Itzhack, I. Y.↗

True covariance simulation of the EUVE (extreme ultra violet explorer) update filter

This paper presents a covariance analysis of the performance and sensitivity of the attitude determination Extended Kalman Filter (EKF) used by the onboard computer (OBC) of the Extreme Ultra Violet Explorer (EUVE) spacecraft. The linearized dynamics and measurement equations of the error states are used in formulating the 'truth model' describing the order of the systems involved. The 'design model' used by the OBC EKF is then obtained by reducing the order of the truth model. The covariance matrix of the EKF which uses the reduced order model is not the correct covariance of the EKF estimation error. A 'true covariance analysis' has to be carried out in order to evaluate the correct accuracy of the OBC generated estimates. The results of such analysis are presented which indicate both the performance and the sensitivity of the OBC EKF.

Bar-Itzhack, I. Y.↗

Macroscopic Violation of Three Cauchy-Schwarz Inequalities Using Correlated Light Beams From an Infra-Red Emitting Semiconductor Diode Array

We briefly review quantum mechanical and semi-classical descriptions of experiments which demonstrate the macroscopic violation of the three Cauchy-Schwarz inequalities: g(sup 2)(sub 11)(0) greater than or equal to 1; g(sup 2)(sub 11)(0) greater than or equal to g(sup 2)(sub 11)(t), (t approaches infinity); (the absolute value of g(sup 2)(sub 11)(0))(exp 2) less than or equal to g(sup 2)(sub 11)(0) g(sup 2)(sub 11)(0). Our measurements demonstrate the violation, at macroscopic intensities, of each of these inequalities. We show that their violation, although weak, can be demonstrated through photodetector current covariance measurements on correlated sub-Poissonian Poissonian, and super Poissonian light beams. Such beams are readily generated by a tandem array of infrared-emitting semiconductor junction diodes. Our measurements utilize an electrically coupled array of one or more infrared-emitting diodes, optically coupled to a detector array. The emitting array is operated in such a way as to generate highly correlated beams of variable photon Fano Factor. Because the measurements are made on time scales long compared with the first order coherence time and with detector areas large compared with the corresponding coherence areas, first order interference effects are negligible. The first and second inequalities are violated, as expected, when a sub-Poissonian light beam is split and the intensity fluctuations of the two split beams are measured by two photodetectors and subsequently cross-correlated. The third inequality is violated by bunched (as well as anti-bunched) beams of equal intensity provided the measured cross correlation coefficient exceeds (F - 1)/F, where F is the measured Fano Factor of each beam. We also investigate the violation for the case of unequal beams.

Edwards, P. J.↗

Impacts of Snow and Cloud Covers on Satellite-Derived PM 2.5 Levels

Satellite aerosol optical depth (AOD) has been widely employed to evaluate ground fine particle (PM 2.5 ) levels, whereas snow/cloud covers often lead to a large proportion of non-random missing AOD. As a result, the fully covered and unbiased PM 2.5 estimates will be hard to generate. Among the current approaches to deal with the data gap issue, few have considered the cloud-AOD relationship and none of them have considered the snow-AOD relationship. This study examined the impacts of snow and cloud covers on AOD and PM 2.5 and made full-coverage PM 2.5 predictions with the consideration of these impacts. To estimate the missing AOD, daily gap-filling models with snow/cloud fractions and meteorological covariates were developed using the random forest algorithm. By using these models in New York State, a daily AOD data set with a 1-km resolution was generated with a complete coverage. The“out-of-bag” R 2 of the gap-filling models averaged 0.93 with an interquartile range from 0.90 to 0.95. Subsequently, a random forest-based PM 2.5 prediction model with the gap-filled AOD and covariates was built to predict fully covered PM 2.5 estimates. A ten-fold cross-validation for the prediction model showed a good performance with an R 2 of 0.82. In the gap-filling models, the snow fraction was of higher significance in the snow season compared with the rest of the year. The prediction models fitted with/without the snow fraction also suggested the discernible changes in PM 2.5 patterns, further confirming the significance of this parameter. Compared with the methods without considering snow and cloud covers, our PM 2.5 prediction surfaces showed more spatial details and reflected small-scale terrain-driven PM 2.5 patterns. The proposed methods can be generalized to the areas with extensive snow/cloud covers and large proportions of missing satellite AOD for predicting PM 2.5 levels with high resolutions and complete coverage.

AOD↗

Evaluating Delta-V Dispersions Using Linear Covariance Techniques with Applications to Rendezvous and Docking

One of the attractive advantages of using linear covariance analysis (LinCov) is that it can accurately produce both navigation errors and trajectory dispersions in a fraction of the time comparable to the standard Monte Carlo analysis approach that is traditionally adopted. It also has the capability of generating delta-v dispersions which is perhaps the most influential performance metric since many aspects of a mission design revolve around the anticipated propellant usage. Accurately capturing this key performance criteria using linear covariance techniques has several subtle caveats raising reasonable doubts and a cautious sense of skepticism. This paper outlines multiple approaches of modeling delta-v usage in LinCov, highlights their limitations and advantages, and ultimately compares their corresponding results to the actual Monte Carlo performance. The theory is also extended to reliably account for an arbitrary thruster configuration and the resultant propellant consumptions due to both translational burns and attitude maneuvers. This provides near instant insight into the impacts of thruster layout design and different vehicle mass properties. These derived techniques are applied and verified using an NRHO rendezvous and docking scenario consistent with upcoming NASA Artemis missions for Orion.

Linear Covariance Analysis↗

Automated Analysis, Classification, and Display of Waveforms

A computer program partly automates the analysis, classification, and display of waveforms represented by digital samples. In the original application for which the program was developed, the raw waveform data to be analyzed by the program are acquired from space-shuttle auxiliary power units (APUs) at a sampling rate of 100 Hz. The program could also be modified for application to other waveforms -- for example, electrocardiograms. The program begins by performing principal-component analysis (PCA) of 50 normal-mode APU waveforms. Each waveform is segmented. A covariance matrix is formed by use of the segmented waveforms. Three eigenvectors corresponding to three principal components are calculated. To generate features, each waveform is then projected onto the eigenvectors. These features are displayed on a three-dimensional diagram, facilitating the visualization of the trend of APU operations.

Kwan, Chiman↗

Integrating in situ environmental covariates in an American lobster catch model to improve impact assessment

The installation and operation of floating offshore wind power is an integral component of societal transition to renewable energy generation where fixed bottom offshore wind is not possible. However, it will cause unique ecosystem changes. To disentangle the effects of offshore wind installations from the concurrent effects of climate change and the fishing practices on commercially significant resources, we must develop detailed characterizations of the resources before development occurs. In the Gulf of Maine, American lobster is the most commercially and culturally important fishery. At the time of writing, this is the largest fishery by value in North America. Our understanding of baseline localized parameters (such as catch per trap at the spatial scale of individual turbines) should be informed by relationships to environmental, biological, and survey-specific functional drivers of catch. A more mechanistic understanding of catch will allow for strategic adjustments to Post- Deployment fishery responses and ultimately, the development of research- and commercial-scale floating offshore wind development. Here, we used survey data from the New England Aqua Ventus Pre-Construction Commercial Trapping Survey to develop Generalized Additive Models describing seasonal catch per trap for legal and sublegal lobsters. We found fall catch to be nearly twice that of spring. Bottom temperature dynamics could be used to predict catch, and the Fall survey was associated with a warmer temperature regime. By using analytical tools that incorporate environmental heterogeneity, we developed monitoring methods from preconstruction baseline data that will be applicable over the post-construction operating period of an offshore wind farm.

BACI↗

Summary of light-element evaluation work: nuclear data uncertainty quantification

Hale & Paris (T-2) have evaluated the 15 N system to provide nuclear cross section and covariance information for nuclear data uncertainty quantification. This memo gives a brief discussion of the R-matrix method used to generate this nuclear scattering and reaction data in the ENDF-6 format.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Influence of atmospherically induced random wave fronts on diffraction imagery - A computer simulation model for testing image reconstruction algorithms

This paper is devoted to the development of a two-dimensional computer-simulation model that is based on the rigid constraints of optical diffraction theory with careful attention paid to the generation of sample realizations of Gaussian-distributed, spatially random, isotropic wave fronts that have zero-mean and prescribed-covariance functions. Given a sample realization of the wave front, the corresponding centered point-spread function and optical-transfer function are evaluated. A detailed study is made of the statistics of random wave-front tilt, point-spread function, modulus squared of transfer function, and phase of transfer function.

Barakat, Richard↗

The unpolarized component in polarimetric radar observations of forested areas

A polarimetric radar measures the complete scattering matrix of a target. For the forest and other natural targets, the average return is, in general, partially polarized. Conditions on the second-order statistics of the scattering matrix under which the return is fully polarized are derived, and it is shown that these conditions are connected with the eigenvalues of the covariance matrix. The result of this analysis is a set of three quantities which indicate the ability of a scatterer to generate an unpolarized component. These quantities for forested areas are shown.

Durden, Stephen L.↗

Uncertainty Quantification of Inducer Natural Frequency using Conditional Assessment of Modeling and Modal Testing of Simpler Systems

The low pressure fuel pump inducer of the new Space Launch System RS25 core stage engine operates in a highly complex environment that substantially affects its modal characteristics. Some of the more important effects are fluid-added mass resulting from operation within a light liquid (Hydrogen), and the magnification of this effect due to tight tip clearance. Since higher order cavitation has been identified as a significant harmonic driver, knowledge of the natural frequency of potentially excitable modes is critical for safe operation, but this frequency cannot be measured during the severe operational environment. A comprehensive testing and analysis program has therefore been performed over the last four years to identify the nominal value and uncertainty of the frequency by modeling and testing four simpler structures which share some of the characteristics of the operational inducer. This testing was used to assess and adjust modeling techniques and excellent correlation was achieved. Identification of the uncertainty in the inducer frequency itself was still problematic, however. This difficulty led to an investigation of Bayesian uncertainty quantification techniques, and to the application of the relatively simple technique of Multi-Variate Normal conditional distributions to calculate the inducer natural frequency uncertainty. Assumptions on prior distributions of uncertainty of the fluid-added mass and tip clearance effect are initially applied to models of each of the simple structures and the inducer itself, and these uncertainties are propagated to generate natural frequencies using design of experiments. Simple response surfaces are then created from this data in order to calculate a Covariance Matrix relating all of these natural frequencies. Finally, the results from modal test of the simple structures are considered to be observations and used to calculate the conditional variance of the desired inducer frequencies. As this method is less rigorous than more complicated Bayesian methods reported in the literature, a conservative factor is applied to the result, but the resulting uncertainty is still significantly less than originally estimated and will greatly assist certification of the inducer for use in the engine.

Uncertainty Quantification↗

Uncertainty Quantification of Inducer Eigenvalues using Conditional Assessment of Models and Modal Test of Simpler Systems

The low pressure fuel pump inducer of the new Space Launch System RS25 core stage engine operates in a highly complex environment that substantially affects its modal characteristics. Some of the more important effects are fluid-added mass resulting from operation within a light liquid (Hydrogen), and the magnification of this effect due to tight tip clearance. Since higher order cavitation has been identified as a significant harmonic driver, knowledge of the natural frequency of potentially excitable modes is critical for safe operation, but this frequency cannot be measured during the severe operational environment. A comprehensive testing and analysis program has therefore been performed over the last four years to identify the nominal value and uncertainty of the frequency by modeling and testing four simpler structures which share some of the characteristics of the operational inducer. This testing was used to assess and adjust modeling techniques and excellent correlation was achieved. Identification of the uncertainty in the inducer frequency itself was still problematic, however. This difficulty led to an investigation of Bayesian uncertainty quantification techniques, and to the application of the relatively simple technique of Multi-Variate Normal conditional distributions to calculate the inducer natural frequency uncertainty. Assumptions on prior distributions of uncertainty of the fluid-added mass and tip clearance effect are initially applied to models of each of the simple structures and the inducer itself, and these uncertainties are propagated to generate natural frequencies using design of experiments. Simple response surfaces are then created from this data in order to calculate a Covariance Matrix relating all of these natural frequencies. Finally, the results from modal test of the simple structures are considered to be observations and used to calculate the conditional variance of the desired inducer frequencies. As this method is less rigorous than more complicated Bayesian methods reported in the literature, a conservative factor is applied to the result, but the resulting uncertainty is still significantly less than originally estimated and will greatly assist certification of the inducer for use in the engine.

Uncertainty Quantification↗

A method for examining ensemble averaging forms during the transition to turbulence in HED systems for application to RANS models

This paper discusses a strategy to initialize a two-dimensional (2D) Reynolds-averaged Navier–Stokes model [LANL's Besnard–Harlow–Rauenzahn (BHR) model] in order to describe an unsteady transitional Richtmyer–Meshkov (RM)-induced flow observed in on-going high-energy-density ensemble experiments performed on the OMEGA-EP facility. The experiments consist of a nominal single-mode perturbation (initial amplitude a 0 ≈ 10 and wavelength $λ$ = 100μm) with target-to-target variations in the surface roughness subjected to the RM instability with delayed Rayleigh–Taylor in a heavy-to-light configuration. Our strategy leverages high-resolution three-dimensional (3D) implicit large eddy simulations (ILES) simulations to initialize BHR-relevant parameters and subsequently validate the 2D BHR results against the 3D ILES simulations. A suite of five 3D ILES simulations corresponding to five experimental target profiles is undertaken to generate an ensemble dataset. Using ensemble averages from the 3D simulations to initialize the turbulent kinetic energy in the BHR model ( K 0 ) demonstrates the ability of the model to predict the time evolution of the interface as well as the density-specific-volume covariance, b . To quantify the sensitivity of the BHR results to the choice of K 0 and the initial turbulent length scale, S 0 , we execute a parameter sweep spanning four orders of magnitude for both S 0 and K 0 , generating a parameter space consisting of 26 simulations. The Pearson's correlation coefficient is used as a measure of discrepancy between the 2D BHR and 3D ILES simulations and reveals that the ranges 8≲S 0 ≲20 μm and 10 9 ≲K 0 ≲10 10 cm 2 /s 2 produce predictions that agree best with the 3D ILES results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sampling Functions from Gaussian Processes and Structured Covariance Gaussian Networks

When learning aerodynamic models from data, it is critical to incorporate estimates of model uncertainty. This motivates the design of probabilistic aerodynamic databases which can be sampled to generate physically and statistically plausible aerodynamic models. In this talk we discuss how to sample deterministic functions from two different kinds of probabilistic models and demonstrate their use. First, Gaussian Process Regressors (GPRs) are a widely used probabilistic kernel-based model which can be thought of as Gaussian distributions over functions. GPRs are generally trained by maximizing the marginal likelihood of seeing the training data over the kernel parameter space. Sample functions are easily generated by drawing points from the Gaussian distribution at desired input points. However, when the points are not known ahead of time, the classical sampling approach is not possible since successive function samples will generate different function realizations. We present an approach for sampling consistent function evaluations from a GPR over multiple samples. Second, we describe a neural network architecture which learns a conditional Gaussian distribution by maximizing the marginal likelihood at each point in the input space. We then discuss and compare several options for generating sample functions which match this distribution. Finally, we demonstrate the use of these probabilistic aerodynamic models in an atmospheric reentry simulation.

Gaussian process regression↗

Space Trajectories Error Analysis (STEAP) Programs. Volume 1: Analytic manual, update

Manual revisions are presented for the modified and expanded STEAP series. The STEAP 2 is composed of three independent but related programs: NOMAL for the generation of n-body nominal trajectories performing a number of deterministic guidance events; ERRAN for the linear error analysis and generalized covariance analysis along specific targeted trajectories; and SIMUL for testing the mathematical models used in the navigation and guidance process. The analytic manual provides general problem description, formulation, and solution and the detailed analysis of subroutines. The programmers' manual gives descriptions of the overall structure of the programs as well as the computational flow and analysis of the individual subroutines. The user's manual provides information on the input and output quantities of the programs. These are updates to N69-36472 and N69-36473.

Source record↗

The Goddard VLBI SINEX Files

The Goddard VLBI group generates SINEX files for individual 24-hr sessions and multiyear TRF solutions. Each 24-hr file includes the reduced normal matrix for the station positions, EOP, and EOP rates. A global Terrestrial Reference Frame (TRF) SINEX file has the reduced normal matrix for station positions and velocities but no EOP parameters. Both types of normal matrix are free from any datum constraints and are designed for combination solutions that will later apply the necessary datum constraints. The SINEX files also contain blocks that provide the constraints and covariance of well-defined 24-hr or global TRF solutions. These blocks are not yet part of the SINEX 2.0 specification. For the TRF solutions the Goddard VLBI group also generates the corresponding EOP time series. SINEX files and EOP time series are distributed through the IVS. The algorithms for generating the various blocks and the proper use of such SINEX files will be discussed.

Ma, Chopo↗

Improving E3SM Land Model Photosynthesis Parameterization via Satellite SIF, Machine Learning, and Surrogate Modeling

The parameterization of key photosynthesis parameters is one of the key uncertain sources in modeling ecosystem gross primary productivity (GPP). Solar-induced chlorophyll fluorescence (SIF) offers a good proxy for GPP since it marks the actual process of photosynthesis; while machine learning (ML) provides a robust approach to model the GPP-SIF relationship. Here, we trained the boosted regressing tree (BRT) and the Random Forest ML models with Greenhouse Gases Observing Satellite SIF data and in situ GPP observations from 49 eddy covariance towers. These trained ML GPP-SIF models were fed into the Energy Exascale Earth System Model (E3SM) Land Model (ELM) to generate ELM-simulated global SIF estimates, which were then benchmarked against satellite SIF observations with a surrogate modeling approach. Our results indicated good modeling performance of the ML-based GPP-SIF relationship. The ELM model when fed with the ML GPP-SIF models also can well predict the spatial-temporal variations in SIF. We also found high model accuracy for the surrogate modeling. Model parameter sensitivity analysis suggested that the fraction of leaf nitrogen in RuBisCO (flnr) is the most sensitive parameter to the SIF; other sensitive parameters include the Ball-Berry stomatal conductance slope (mbbopt) and the vcmax entropy (vcmaxse). The posterior uncertainty in simulated GPP was greatly reduced after benchmarking, and the model produced improved spatial patterns of mean GPP relative to FLUXCOM GPP. Our integrated approach provides a new avenue for improving land models and using remote-sensing SIF, which can be further improved in the future with more ground- and satellite-based observations.

54 ENVIRONMENTAL SCIENCES↗

Misaligned Wind‐Waves Behind Atmospheric Cold Fronts

Abstract Atmospheric fronts embedded in extratropical cyclones are high‐impact weather phenomena, contributing significantly to mid‐latitude winter precipitation. The three vital characteristics of the atmospheric fronts, high wind speeds, abrupt change in wind direction, and rapid translation, force the induced surface waves to be misaligned with winds exclusively behind the cold fronts. The effects of the misaligned waves under atmospheric cold fronts on air‐sea fluxes remain undocumented. Using the multi‐year in situ near‐surface observations and direct covariance flux measurements from the Pioneer Array off the coast of New England, we find that the majority of the passing cold fronts generate misaligned waves behind the cold front. Once generated, the waves remain misaligned, on average, for about 8 hr. The parameterized effect of misaligned waves in a fully coupled model significantly increases the roughness length (185%), drag coefficient (19%), and air‐sea momentum flux (11%). The increased surface drag reduces the wind speeds in the surface layer. The upward turbulent heat flux is weakly decreased by the misaligned waves because of the decrease in temperature and humidity scaling parameters being greater than the increase in friction velocity. The misaligned wave effect is not accurately represented in a commonly used wave‐based bulk flux algorithm. Yet, considering this effect in the current formulation improves the overall accuracy of parameterized momentum flux estimates. The results imply that better representing a directional wind‐wave coupling in the bulk formula of the numerical models may help improve the air‐sea interaction simulations under the passing atmospheric fronts in the mid‐latitudes.

Oceanography↗