SEARCH · Engineering Papers

Results for “covariance constraints”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2026

NREL ASSIST Rhode Island / Thermodynamic retrievals TROPoe

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous observations from the NREL ASSIST-II (SN 11) infrared spectrometer. Observations are noise-filtered but not averaged in time to minimize errors due to non-uniform clouds. Additional input data in TROPoe are cloud base height from a collocated NREL Upgraded Galion lidar #1212-60 and temperature, RH, and pressure from collocated PNNL met station. The lidar beta is used to estimate the cloud base height through a modified algorithm by Newsom et al., 2019 used at ARM SGP. The full pipeline for running the retrieval is available at https://github.com/StefanoWind/TROPoe_processor. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) which provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see e.g. Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Upton, NY.

17 WIND ENERGY↗

DOE OSTI2021

Gaussian process regression constrained by boundary value problems

We develop a framework for Gaussian processes regression constrained by boundary value problems. The framework may be applied to infer the solution of a well-posed boundary value problem with a known second-order differential operator and boundary conditions, but for which only scattered observations of the source term are available. Scattered observations of the solution may also be used in the regression. The framework combines co-kriging with the linear transformation of a Gaussian process together with the use of kernels given by spectral expansions in eigenfunctions of the boundary value problem. Furthermore, it benefits from a reduced-rank property of covariance matrices. We demonstrate that the resulting framework yields more accurate and stable solution inference as compared to physics-informed Gaussian process regression without boundary condition constraints.

42 ENGINEERING↗

NASA NTRS2009

The Principle of Energetic Consistency

A basic result in estimation theory is that the minimum variance estimate of the dynamical state, given the observations, is the conditional mean estimate. This result holds independently of the specifics of any dynamical or observation nonlinearity or stochasticity, requiring only that the probability density function of the state, conditioned on the observations, has two moments. For nonlinear dynamics that conserve a total energy, this general result implies the principle of energetic consistency: if the dynamical variables are taken to be the natural energy variables, then the sum of the total energy of the conditional mean and the trace of the conditional covariance matrix (the total variance) is constant between observations. Ensemble Kalman filtering methods are designed to approximate the evolution of the conditional mean and covariance matrix. For them the principle of energetic consistency holds independently of ensemble size, even with covariance localization. However, full Kalman filter experiments with advection dynamics have shown that a small amount of numerical dissipation can cause a large, state-dependent loss of total variance, to the detriment of filter performance. The principle of energetic consistency offers a simple way to test whether this spurious loss of variance limits ensemble filter performance in full-blown applications. The classical second-moment closure (third-moment discard) equations also satisfy the principle of energetic consistency, independently of the rank of the conditional covariance matrix. Low-rank approximation of these equations offers an energetically consistent, computationally viable alternative to ensemble filtering. Current formulations of long-window, weak-constraint, four-dimensional variational methods are designed to approximate the conditional mode rather than the conditional mean. Thus they neglect the nonlinear bias term in the second-moment closure equation for the conditional mean. The principle of energetic consistency implies that, to precisely the extent that growing modes are important in data assimilation, this term is also important.

Cohn, Stephen E.↗

NASA NTRS2010

GLGM-3: A Degree-ISO Lunar Gravity Model from the Historical Tracking Data of NASA Moon Orbiters

In preparation for the radio science experiment of the Lunar Reconnaissance Orbiter (LRO) mission, we analyzed the available radio tracking data of previous NASA lunar orbiters. Our goal was to use these historical observations in combination with the new low-altitude data to be obtained by LRO. We performed Precision Orbit Determination on trajectory arcs from Lunar Orbiter 1 in 1966 to Lunar Prospector in 1998, using the GEODYN II program developed at NASA Goddard Space Flight Center. We then created a set of normal equations and solved for the coefficients of a spherical harmonics expansion of the lunar gravity potential up to degree and order 150. The GLGM-3 solution obtained with a global Kaula constraint (2.5 x 10(exp -4)/sq l) shows good agreement with model LP150Q from the Jet Propulsion Laboratory, especially over the nearside. The levels of data fit with both gravity models are very similar (Doppler RMS of approx.0.2 and approx. 1-2 mm/s in the nominal and extended phases, respectiVely). Orbit overlaps and uncertainties estimated from the covariance matrix also agree well. GLGM-3 shows better correlation with lunar topography and admittance over the nearside at high degrees of expansion (l > 100), particularly near the poles. We also present three companion solutions, obtained with the same data set but using alternate inversion strategies that modify the power law constraint and expectation of the individual spherical harmonics coefficients. We give a detailed discussion of the performance of this family of gravity field solutions in terms of observation fit, orbit quality, and geophysical consistency.

Mazarico, E.↗

DOE OSTI2021

Cosmic shear cosmology beyond two-point statistics: a combined peak count and correlation function analysis of DES-Y1

ABSTRACT We constrain cosmological parameters from a joint cosmic shear analysis of peak-counts and the two-point shear correlation functions, as measured from the Dark Energy Survey (DES-Y1). We find the structure growth parameter $S_8\equiv \sigma _8\sqrt{\Omega _{\rm m}/0.3} = 0.766^{+0.033}_{-0.038}$ which, at 4.8 per cent precision, provides one of the tightest constraints on S8 from the DES-Y1 weak lensing data. In our simulation-based method we determine the expected DES-Y1 peak-count signal for a range of cosmologies sampled in four w cold dark matter parameters (Ωm, σ8, h, w0). We also determine the joint covariance matrix with over 1000 realizations at our fiducial cosmology. With mock DES-Y1 data we calibrate the impact of photometric redshift and shear calibration uncertainty on the peak-count, marginalizing over these uncertainties in our cosmological analysis. Using dedicated training samples we show that our measurements are unaffected by mass resolution limits in the simulation, and that our constraints are robust against uncertainty in the effect of baryon feedback. Accurate modelling for the impact of intrinsic alignments on the tomographic peak-count remains a challenge, currently limiting our exploitation of cross-correlated peak counts between high and low redshift bins. We demonstrate that once calibrated, a fully tomographic joint peak-count and correlation functions analysis has the potential to reach a 3 per cent precision on S8 for DES-Y1. Our methodology can be adopted to model any statistic that is sensitive to the non-Gaussian information encoded in the shear field. In order to accelerate the development of these beyond-two-point cosmic shear studies, our simulations are made available to the community upon request.

Harnois-Déraps, Joachim↗

NASA NTRS1972

Computation of optimal output-feedback compensators for linear time-invariant systems

The control of linear time-invariant systems with respect to a quadratic performance criterion was considered, subject to the constraint that the control vector be a constant linear transformation of the output vector. The optimal feedback matrix, f*, was selected to optimize the expected performance, given the covariance of the initial state. It is first shown that the expected performance criterion can be expressed as the ratio of two multinomials in the element of f. This expression provides the basis for a feasible method of determining f* in the case of single-input single-output systems. A number of iterative algorithms are then proposed for the calculation of f* for multiple input-output systems. For two of these, monotone convergence is proved, but they involve the solution of nonlinear matrix equations at each iteration. Another is proposed involving the solution of Lyapunov equations at each iteration, and the gradual increase of the magnitude of a penalty function. Experience with this algorithm will be needed to determine whether or not it does, indeed, possess desirable convergence properties, and whether it can be used to determine the globally optimal f*.

Platzman, L. K.↗

NASA NTRS1974

Sequential parameter estimation using pseudoinverse

This correspondence presents a method for solution of the parameter estimation problem in an entirely sequential fashion. In this procedure each measurement is processed as it is received, without regard to whether it and the prior measurements have established a minimal observable data set, or without the assumption a priori of a covariance matrix in the manner of the Kalman filter. During the period prior to the establishment of the minimal data set the estimate maintains the character of the pseudoinverse constraint, that of minimum norm. In addition at each measurement, a parameter is available for deciding whether this measurement has increased the observability of the system.

Nayak, R. P.↗

DOE OSTI2024

Two-Dimensional Projected-Momentum Covariance Mapping for Coulomb Explosion Imaging

We introduce projected-momentum covariance mapping, an extension of recoil-frame covariance mapping for 2D ion imaging studies. By considering the two-dimensional projection of the ion momenta as recorded by the detector, one opens the door to a complex suite of analysis tools adapted from three-dimensional momentum imaging studies. This includes the use of different frames of reference to unravel the dynamics of fragmentation and the application of fragment momentum constraints to isolate specific fragmentation channels. The technique is demonstrated on data from a two-dimensional ion imaging study of the Coulomb explosion of the cis and trans isomers of 1,2-dichloroethene, following strong-field ionization by an intense near-infrared femtosecond laser pulse. Classical simulations are used to guide the interpretation of projected-momentum covariance maps. The results offer a detailed insight into the distinct Coulomb explosion dynamics for this pair of isomers and lay the groundwork for future time-resolved studies of photoisomerization dynamics in this molecular system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Measure This, Not That: Pareto Optimal Trade-Offs between Model-Based Information Content and Measurements Cost

The slides present a novel convex optimization formulation to compute the best set of measurements for multi-response dynamical systems with asynchronous time steps that maximize the Fisher information content subject to budget constraints. The trace (A-optimality) or determinant (D-optimality) of the Fisher Information Matrix (FIM) quantifies the information content. The framework supports arbitrary (positive semi-definite) variance and covariances between every pair of responses and their time steps.

Wang, Jialu↗

DOE OSTI2024

Quantum reference frames from top-down crossed products

All physical observations are made relative to a reference frame, which is a system in its own right. If the system of interest admits a group symmetry, the reference frame observing it must transform commensurately under the group to ensure the covariance of the combined system. We point out that the crossed product is a way to realize quantum reference frames from the bottom-up; adjoining a quantum reference frame and imposing constraints generates a crossed product algebra. We provide a top-down specification of crossed product algebras and show that one cannot obtain inequivalent quantum reference frames using this approach. As a remedy, we define an abstract algebra associated to the system and symmetry group built out of relational crossed product algebras associated with different choices of quantum reference frames. We term this object the G -framed algebra, and show how potentially inequivalent frames are realized within this object. We comment on this algebra’s analog of the classical Gribov problem in gauge theory, its importance in gravity where we show that it is relevant for semiclassical de Sitter and potentially beyond the semiclassical limit, and its utility for understanding the frame dependence of physical notions like observables, density states, and entropies. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

NASA NTRS1990

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↗

NASA NTRS1986

Higher-order structure of rRNA

A comparative search for phylogenetically covarying basepair replacements within potential helices has been the only reliable method to determine the correct secondary structure of the 3 rRNAs, 5S, 16S, and 23S. The analysis of 16S from a wide phylogenetic spectrum, that includes various branches of the eubacteria, archaebacteria, eucaryotes, in addition to the mitochondria and chloroplast, is beginning to reveal the constraints on the secondary structures of these rRNAs. Based on the success of this analysis, and the assumption that higher order structure will also be phylogenetically conserved, a comparative search was initiated for positions that show co-variation not involved in secondary structure helices. From a list of potential higher order interactions within 16S rRNA, two higher-order interactions are presented. The first of these interactions involves positions 570 and 866. Based on the extent of phylogenetic covariation between these positions while maintaining Watson-Crick pairing, this higher-order interaction is considered proven. The other interaction involves a minimum of six positions between the 1400 and 1500 regions of the 16S rRNA. Although these patterns of covariation are not as striking as the 570/866 interaction, the fact that they all exist in an anti-parallel fashion and that experimental methods previously implicated these two regions of the molecule in tRNA function suggests that these interactions be given serious consideration.

Gutell, R. R.↗

DOE OSTI2021

Implications of PREX-2 on the electric dipole polarizability of neutron-rich nuclei

Background: The recent announcement by the PREX collaboration of an unanticipated thick neutron skin in 208 Pb ($R^{208}_{skin}$) has challenged the understanding of neutron-rich matter in the vicinity of nuclear saturation density. Whereas earlier constraints indicate that the symmetry energy is relatively soft, the PREX-2 result seems to suggest the opposite. Purpose: To confront constraints on the symmetry energy obtained from measurements of the electric dipole polarizability against those informed by the PREX-2 measurement of $R^{208}_{skin}$ and by the correlations that it entails. Methods: Here, covariant energy density functionals informed by the properties of finite nuclei are used to compute the electric dipole response of 48 Ca, 68 Ni, 132 Sn, and 208 Pb. The set of functionals used in this work are consistent with experimental data, yet are flexible enough in that they span a wide range of values of $R^{208}_{skin}$ Results: It is found that theoretical predictions of the electric dipole polarizability that are consistent with the PREX-2 measurement systematically overestimate the corresponding values extracted from the direct measurements of the distribution of electric dipole strength. Conclusions: The neutron skin thickness of 208 Pb extracted from parity violating electron scattering and the electric dipole polarizability measured in photoabsorption experiments are two of the cleanest experimental tools used to constrain the symmetry energy around nuclear saturation density. However, the recent value of $R^{208}_{skin}$ that suggests a fairly stiff symmetry energy stands in stark contrast to the conclusions derived from the electric dipole polarizability. At present, I offer no solution to this dilemma.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2002

Site distribution and aliasing effects in the inversion for load coefficients and geocenter motion from GPS data

Precise GPS measurements of elastic relative site displacements due to surface mass loading offer important constraints on global surface mass transport. We investigate effects of site distribution and aliasing by higher- degree (n ≥ 2) loading terms on inversion of GPS data for n = 1 load coefficients and geocenter motion. Covariance and simulation analyses are conducted to assess the sensitivity of the inversion to aliasing and mismodeling errors and possible uncertainties in the n = 1 load coefficient determination. We found that the use of center-of-figure approximation in the inverse formulation could cause 10– 15% errors in the inverted load coefficients. n = 1 load estimates may be contaminated significantly by unknown higher-degree terms, depending on the load scenario and the GPS site distribution. The uncertainty in n = 1 zonal load estimate is at the level of 80–95% for two load scenarios.

Webb, Frank H.↗

NASA NTRSDate not supplied

Robust Cislunar Trajectory Optimization Via Midcourse Correction and Optical Navigation Scheduling

This paper presents a new approach to optimal trajectory design that considers uncertainties in the system, referred to herein as robust trajectory optimization. This approach assumes an existing reference trajectory and optimizes the locations of midcourse correction burns and utilization of onboard navigation sensors to minimize dispersions in ∆v or final position. Navigation errors, maneuver execution errors, orbit insertion errors, and environmental modeling errors are considered. The application in this paper is cislunar flight with the goal of injecting into a Near-Rectilinear Halo Orbit for rendezvous with a target vehicle. Two complementary optimization problems are proposed. One problem minimizes the total ∆v dispersion subject to a final position dispersion constraint. The other problem minimizes the final position dispersion subject to a total ∆v dispersion constraint. The results from each optimization problem are shown for a complete mission profile.

Linear Covariance Analysis↗

NASA NTRSDate not supplied

Generalized Linear Targeting For Cislunar Flight

An important element of Artemis and NASA’s campaign to explore the Moon is the autonomous onboard two-level targeter (TLT) used during all cislunar flight phases. The function of the TLT is to autonomously recompute the burn targets for the upcoming burn (or multiple burns) in response to navigation and vehicle dispersion providing a solution that meets all of the trajectory constraints. Although the TLT has been utilized previously as a ground-based planning tool, and flown onboard during the Artemis I mission, it’s complexity and iterative nature make is difficult to incorporate into and support rapid analyses such as robust optimal trajectory design applications where speed is essential. In this paper, a set of generalized linear targeting algorithms that mimics many of the properties of the TLT is derived. The generalized algorithms can handle single or multiple impulsive maneuvers, with multiple constraints at multiple fixed or variable times. A linear targeting algorithm for finite burn maneuvers is also derived. The generalized linear targeting algorithms are exceptionally fast and easy to implement in Monte Carlo analysis, linear covariance (LinCov) analysis, and robust optimal trajectory design. Several cislunar flight examples are provided.

Linear Covariance Analysis↗

DOE OSTI2024

Testing the ΛCDM Cosmological Model with Forthcoming Measurements of the Cosmic Microwave Background with SPT-3G

We forecast constraints on cosmological parameters enabled by three surveys conducted with SPT-3G, the third-generation camera on the South Pole Telescope. The surveys cover separate regions of 1500, 2650, and 6000 deg$^{2}$ to different depths, in total observing 25% of the sky. These regions will be measured to white noise levels of roughly 2.5, 9, and 12μK-armin, respectively, in cosmic microwave background (CMB) temperature units at 150 GHz by the end of 2024. The survey also includes measurements at 95 and 220 GHz, which have noise levels a factor of ∼1.2 and 3.5 times higher than 150 GHz, respectively, with each band having a polarization noise level ∼2times higher than the temperature noise. We use a novel approach to obtain the covariance matrices for jointly and optimally estimated gravitational lensing potential band powers and unlensed CMB temperature and polarization band powers. We demonstrate the ability to test the ΛCDM model via the consistency of cosmological parameters constrained independently from SPT-3G and Planck data, and consider the improvement in constraints on ΛCDM extension parameters from a joint analysis of SPT-3G and Planck data. The ΛCDM cosmological parameters are typically constrained with uncertainties up to ∼2 times smaller with SPT-3G data, compared to Planck, with the two data sets measuring significantly different angular scales and polarization levels, providing additional tests of the standard cosmological model.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Response of Ecosystem Productivity to High Vapor Pressure Deficit and Low Soil Moisture: Lessons Learned From the Global Eddy-Covariance Observations

Although there is mounting concern about how high vapor pressure deficit (VPD) and low soil moisture (SM) affect ecosystem productivity, their relative importance is still under debate. Here, we comprehensively quantified the relative impacts of these two factors on ecosystem gross primary production (GPP) using observations from a global network of eddy-covariance towers and two approaches (sensitivity analysis and linear regression model). Both approaches agree that a higher percentage of sites experience GPP reduction from high VPD than from low SM over the growing season. However, the constraint of high VPD and low SM on GPP reduction is tightly linked with climates and plant functional types. Humid and mesic ecosystems including forests and grasslands are dominated by VPD, while the semi-arid and arid ecosystems including shrublands and savannas are dominated by SM. The varying dominant role of these two factors on GPP is closely related to plant stomatal behavior, as predicted by a stomatal conductance model. Additionally, we highlight the non-linear impact of SM on GPP during droughts and the possible underestimation of the SM effects for deep-rooted plants when only using surface-layer SM. Our results shed light on a better understanding of the impacts of VPD and SM on vegetation productivity, with important implications for modeling the response and feedback of ecosystem dynamics to current and future climates.

54 ENVIRONMENTAL SCIENCES↗