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 145 records · Page 8

Validation of the DESI 2024 Lyα forest BAO analysis using synthetic datasets

The first year of data from the Dark Energy Spectroscopic Instrument (DESI) contains the largest set of Lyman-α (Lyα) forest spectra ever observed. This data, collected in the DESI Data Release 1 (DR1) sample, has been used to measure the Baryon Acoustic Oscillation (BAO) feature at redshift z = 2.33. In this work, we use a set of 150 synthetic realizations of DESI DR1 to validate the DESI 2024 Lyα forest BAO measurement presented in [1]. The synthetic data sets are based on Gaussian random fields using the log-normal approximation. We produce realistic synthetic DESI spectra that include all major contaminants affecting the Lyα forest. The synthetic data sets span a redshift range 1.8 < z < 3.8, and are analyzed using the same framework and pipeline used for the DESI 2024 Lyα forest BAO measurement. To measure BAO, we use both the Lyα auto-correlation and its cross-correlation with quasar positions. We use the mean of correlation functions from the set of DESI DR1 realizations to show that our model is able to recover unbiased measurements of the BAO position. We also fit each mock individually and study the population of BAO fits in order to validate BAO uncertainties and test our method for estimating the covariance matrix of the Lyα forest correlation functions. Finally, we discuss the implications of our results and identify the needs for the next generation of Lyα forest synthetic data sets, with the top priority being to simulate the effect of BAO broadening due to non-linear evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Higher-derivative relations between scalars and gluons

We extend the covariant color-kinematics duality introduced by Cheung and Mangan to effective field theories. We focus in particular on relations between the effective field theories of gluons only and of gluons coupled to bi-adjoint scalars. Maps are established between their respective equations of motion and between their tree-level scattering amplitudes. An additional rule for the replacement of flavor structures by kinematic factors realizes the map between higher-derivative amplitudes. As an example of new relations, the pure-gluon amplitudes of mass dimension up to eight, featuring insertions of the F 3 and F 4 operators which satisfy the traditional color-kinematics duality, can be generated at all multiplicities from just renormalizable amplitudes of gluons and bi-adjoint scalars. We also obtain closed-form expressions for the kinematic numerators of the dimension-six gluon effective field theory, which are valid in D space-time dimensions. Finally, we find strong evidence that this extended covariant color-kinematics duality relates the (DF) 2 +YM(+Φ 3 ) theories which, at low energies, generate infinite towers of operators satisfying the traditional color-kinematics duality, beyond aforementioned F 3 and F 4 ones.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NGEE Arctic Tram: Radiation and Environmental Measurements over Low- and High-Centered Polygon Vegetation Communities, Barrow, Alaska, 2014-2017

These data were generated from an observational platform (Tram), an automated cart carrying a suite of radiation and remote sensing instruments running on 68 meters of elevated track 1 to 1.5 meters above the surface of the vegetation communities of low- and high-centered polygon ground. The Tram is in the footprint of the NGEE Arctic/AmeriFlux eddy covariance flux tower (US-NGB) at the Barrow Environmental Observatory, Barrow, Alaska, and operated before, during, and after the growing seasons of 2014-2017. Dataset includes one user guide (*.pdf) with final and raw *.csv data files. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Generalized Augmented-State Covariance Analysis for Spaceflight

The use of linear covariance analysis techniques, also known as LinCov, has been used extensively for more than a half century for spaceflight applications. Originally, its primary purpose was to facilitate navigation analysis. For many past and current applications, the specific implementations only support navigation studies still. When the concept of an augmented-state linear covariance analysis approach was initially introduced that allowed for both navigation and trajectory dispersion analysis, the enhancement was motivated and primarily utilized to support navigation filter tuning and error budget analysis. Relatively few utilize this alternate augmented-state formulation of LinCov due to its additional complexity. The untapped potential of the augmented-state linear covariance analysis technique slowly unfolded in the past two-decades as its capability to rapidly and reliably capture the integrated closed-loop guidance, navigation, and control (GN&C) system performance became more apparent. Even with this dual purpose of generating insights to both navigation errors along with trajectory and delta-v dispersions, the core theoretical development had a heavy emphasis on the impacts of the navigation system and largely neglected the details of the actual guidance, targeting, and control systems. This paper extends the navigation-centric theoretical development by formulating a generalized augmented-state covariance analysis (GAUSCOV) technique that allows for the intricacies of a variety of targeting and control strategies along with ground planning and mission operations to be more formally included in assessing the impacts to spaceflight GN&C system performance.

Linear Covariance Analysis↗

Flexible nonstationary spatiotemporal modeling of high-frequency monitoring data

Many physical datasets are generated by collections of instruments that make measurements at regular time intervals. For such regular monitoring data, we extend the framework of half-spectral covariance functions to the case of nonstationarity in space and time and demonstrate that this method provides a natural and tractable way to incorporate complex behaviors into a covariance model. Further, we use this method with fully time-domain computations to obtain bona fide maximum likelihood estimators—as opposed to using Whittle-type likelihood approximations, for example—that can still be computed conveniently. Additionally, we apply this method to very high-frequency Doppler LIDAR vertical wind velocity measurements, demonstrating that the model can expressively capture the extreme nonstationarity of dynamics above and below the atmospheric boundary layer and, more importantly, the interaction of the process dynamics across it.

54 ENVIRONMENTAL SCIENCES↗

On the distribution of source coherency in subsonic jets

Recent experimental and theoretical studies on the sources of far field noise generated by subsonic circular jets have given convincing indications that a coherent wave motion does exist in the potential core of subsonic jets. Cross-covariance measurements using single and cross hot-wire probes in the potential core as well as the turbulent mixing region and microphone measurements in the far field have resulted in cross-correlation coefficients of sufficiently high magnitudes to indicate that significant, if not dominant, noise sources may be located in the potential core. The paper presents the analytical and experimental results of this research.

Seiner, J. M.↗

Refinement of Earth's gravity field with Topex GPS measurements

The NASA Ocean Topography Experiment satellite TOPEX will carry a microwave altimeter accurate to a few centimeters for the measurement of ocean height. The capability can be fully exploited only if TOPEX altitude can be independently determined to 15 cm or better. This in turn requires an accurate gravity model. The gravity will be tuned with selected nine 10-day arcs of laser ranging, which will be the baseline tracking data type, collected in the first six months of TOPEX flight. TOPEX will also carry onboard an experimental Global Positioning System (GPS) flight receiver capable of simultaneously observing six GPS satellites above its horizon to demonstrate the capability of GPS carrier phase and P-code pseudorange for precise determination of the TOPEX orbit. It was found that subdecimeter orbit accuracy can be achieved with a mere two-hour arc of GPS tracking data, provided that simultaneous measurements are also made at six of more ground tracking sites. The precision GPS data from TOPEX are also valuable for refining the gravity model. An efficient technique is presented for gravity tuning using GPS measurements. Unlike conventional global gravity tuning, this technique solves for far fewer gravity parameters in each filter run. These gravity parameters yield local gravity anomalies which can later be combined with the solutions over other parts of the earth to generate a global gravity map. No supercomputing power will be needed for such combining. The approaches used in this study are described and preliminary results of a covariance analysis presented.

Wu, Sien-Chong↗

Multidimensional Distributional Neural Network Output Demonstrated in Super‐Resolution of Surface Wind Speed

Accurate quantification of uncertainty in neural network predictions remains a central challenge for scientific applications involving high-dimensional, correlated data. While existing methods capture either aleatoric or epistemic uncertainty, few offer closed-form, multidimensional distributions that preserve spatial correlation while remaining computationally tractable. In this work, we present a framework for training neural networks with a multidimensional Gaussian loss, generating a closed-form predictive distribution over outputs informed by non-identically distributed training data. Our approach captures aleatoric uncertainty by iteratively estimating the means and covariance matrices, and is demonstrated on a super-resolution example out-of-training-sample. We leverage a Fourier representation of the covariance matrix to stabilize network training and preserve spatial correlation. We introduce a novel regularization strategy—referred to as information sharing—that interpolates between image-specific and global covariance estimates, enabling convergence of the super-resolution downscaling network trained on image-specific distributional loss functions. This framework allows for efficient sampling, explicit correlation modeling, and extensions to more complex distribution families all without disrupting prediction performance. We demonstrate the method on a surface wind speed downscaling task and discuss its broader applicability to uncertainty-aware prediction in scientific models.

17 WIND ENERGY↗

Quantitative radiography for determining density fluctuations in HED experiments

We have developed a method to extract density fluctuation measurements from x-ray radiographs of high-energy density (HED) instability growth and turbulence experiments. We use this information to calculate density fluctuation statistics for constraining the performance of turbulent mix models in HED systems. The density calculation combines image filtering, removal of systemic effects such as backlighter variation, calculation of transmission across multiple materials, and use of tracer materials to generate an approximate single-material density field. From the density map, we calculate both average density and a variance-like moment b (density-specific-volume covariance), which we compare to our models. We infer both quantities from a single image, which is significantly more information than the historic single scalar mix width measurements. We also develop a method of analyzing simulation outputs that incorporate both the density fluctuation metric from a turbulence model and the bulk material maps from the hydrodynamic code. This analysis helps address the question of how to initialize the simulations for best comparison to data from systems with large separations of scale in the mixing perturbation initial condition. We find that our data analysis method yields 1D average density and b curves with similar morphology and amplitudes as those from preliminary simulation comparisons.

47 OTHER INSTRUMENTATION↗

Approaches to Evaluating Probability of Collision Uncertainty

While the two-dimensional probability of collision (Pc) calculation has served as the main input to conjunction analysis risk assessment for over a decade, it has done this mostly as a point estimate, with relatively little effort made to produce confidence intervals on the Pc value based on the uncertainties in the inputs. The present effort seeks to try to carry these uncertainties through the calculation in order to generate a probability density of Pc results rather than a single average value. Methods for assessing uncertainty in the primary and secondary objects' physical sizes and state estimate covariances, as well as a resampling approach to reveal the natural variability in the calculation, are presented; and an initial proposal for operationally-useful display and interpretation of these data for a particular conjunction is given.

PROBABILITY↗

Structured Covariance Gaussian Networks for Orion Crew Module Aerodynamic Uncertainty Quantification

In this paper we propose a new approach for nonlinear regression and uncertainty quantification. The method is based on a pair of neural networks which parameterize mean and dense covariance functions of a multivariate Gaussian process, trained together to maximize the log-likelihood of observing the given data. The covariance matrix is made positive definite at every input by construction. We also propose a sampling approach that produces viable surrogate function realizations from the Gaussian process. We call the proposed model a Structured Covariance Gaussian Network (SCGN). We illustrate the use of SCGNs for learning an aerodynamic response surface with built-in uncertainty for the Orion crew module. We find that SCGN provides an efficient and systematic way to learn nonlinear functional relationships and dense covariances. We compare results to a baseline Gaussian process regressor and observe that the SCGN provides comparable uncertainty descriptions with improved scalability to dataset size. The sample functions generated by SCGN are fast to evaluate online and are therefore convenient for use in trajectory simulations. These results suggest that SCGN may be a viable computational method for aerodynamic uncertainty quantification.

machine learning↗

Structured Covariance Gaussian Networks for Orion Crew Module Aerodynamic Uncertainty Quantification

In this paper we propose a new approach for nonlinear regression and uncertainty quantification. The method is based on a pair of neural networks which parameterize mean and dense covariance functions of a multivariate Gaussian process, trained together to maximize the log-likelihood of observing the given data. The covariance matrix is made positive definite at every input by construction. We also propose a sampling approach that produces viable surrogate function realizations from the Gaussian process. We call the proposed model a Structured Covariance Gaussian Network (SCGN). We illustrate the use of SCGNs for learning an aerodynamic response surface with built-in uncertainty for the Orion crew module. We find that SCGN provides an efficient and systematic way to learn nonlinear functional relationships and dense covariances. We compare results to a baseline Gaussian process regressor and observe that the SCGN provides comparable uncertainty descriptions with improved scalability to dataset size. The sample functions generated by SCGN are fast to evaluate online and are therefore convenient for use in trajectory simulations. These results suggest that SCGN may be a viable computational method for aerodynamic uncertainty quantification.

machine learning↗

Non-asymptotic analysis of ensemble Kalman updates: effective dimension and localization

Many modern algorithms for inverse problems and data assimilation rely on ensemble Kalman updates to blend prior predictions with observed data. Ensemble Kalman methods often perform well with a small ensemble size, which is essential in applications where generating each particle is costly. This paper develops a non-asymptotic analysis of ensemble Kalman updates, which rigorously explains why a small ensemble size suffices if the prior covariance has moderate effective dimension due to fast spectrum decay or approximate sparsity. Here, we present our theory in a unified framework, comparing everal implementations of ensemble Kalman updates that use perturbed observations, square root filtering and localization. As part of our analysis, we develop new dimension-free covariance estimation bounds for approximately sparse matrices that may be of independent interest.

Mathematics↗

Data and Code for: Observation-constrained agroecosystem model inversion reveals continental-scale variation of winter wheat traits

This repository contains the simulation outputs and processing scripts associated with the study of winter wheat traits across the United States, utilizing the Ecosys agroecosystem model. The dataset includes model results for both rainfed and irrigated winter wheat systems, supporting the findings presented in the manuscript titled "Observation-constrained agroecosystem model inversion reveals continental-scale variation of winter wheat traits." Data includes the original Ecosys simulation outputs (archived in .db format within the compressed .zip files) and extracted analysis data (stored in .pkl files for efficient processing). Python code for data processing and figure generation is provided in a Jupyter notebook. External Observational Datasets should refer to the following official repositories for the input and validation data used in this study. The eddy covariance data from the AmeriFlux network (https://ameriflux.lbl.gov/). Climate-forcing data of NLDAS-2 from NASA LDAS (https://ldas.gsfc.nasa.gov/nldas/nldas-2-forcing-data). Soil data from the Gridded Soil Survey Geographic Database (gSSURGO), available at (https://www.nrcs.usda.gov/resources/data-and-reports/gridded-soil-survey-geographic-gssurgo-database). Crop yields, planting and harvest dates from the USDA public databases (https://quickstats.nass.usda.gov/; https://webapp.rma.usda.gov/apps/actuarialinformationbrowser/CropCriteria.aspx). Satellite-derived SLOPE GPP data from ORNL DAAC (https://daac.ornl.gov/cgi-bin/dsviewer.pl?ds_id=1786). Land use and crop progress information from the USDA Crop Data Layer and Crop Progress and Condition Gridded Layers (https://www.nass.usda.gov/Research_and_Science/). The Ecosys model code is available online at https://github.com/jinyun1tang/ECOSYS.

Wheat↗

Theory of electromagnetic transmission structures. I - Relativistic foundation and network formalisms

A new theorem on a class of four-dimensional skew-symmetric tensors is demonstrated. Coupled with the relativistic covariant form of Maxwell's equations, this theorem consolidates the classifications of guided waves by combining the three types - TE, TM, TEM - under a uniform condition applied to the generating four-potential which is Lorentz invariant. Each type corresponds to a potential of which a pair of the four components vanishes in a particular frame. Through appropriate normalization conditions, the resulting time-domain equations for the field amplitudes are readily reduced to modified telegraphist equations, which in turn lead to distributed network representations for each of the three types. The ambiguity of distributed network formalisms in general is elucidated and the concept of network parameter densities such as traditionally employed in TEM transmission line theory is questioned.

Gabriel, G. J.↗

Insights Into the Aerodynamic Versus Radiometric Surface Temperature Debate in Thermal-Based Evaporation Modeling

Global evaporation monitoring from Earth observation thermal infrared satellite missions is historically challenged due to the unavailability of any direct measurements of aerodynamic temperature. State-of-the-art one-source evaporation models use remotely sensed radiometric surface temperature as a substitute for the aerodynamic temperature and apply empirical corrections to accommodate for their inequality. This introduces substantial uncertainty in operational drought mapping over complex landscapes. By employing a non-parametric model, we show that evaporation can be directly retrieved from thermal satellite data without the need of any empirical correction. Independent evaluation of evaporation in a broad spectrum of biome and aridity yielded statistically significant results when compared with eddy covariance observations. While our simplified model provides a new perspective to advance spatio-temporal evaporation mapping from any thermal remote sensing mission, the direct retrieval of aerodynamic temperature also generates the highly required insight on the critical role of biophysical interactions in global evaporation research.

54 ENVIRONMENTAL SCIENCES↗

Beyond the visible: Accounting for ultraviolet and far‐red radiation in vegetation productivity and surface energy budgets

Photosynthetically active radiation (PAR) is typically defined as light with a wavelength within 400–700 nm. However, ultra-violet (UV) radiation within 280–400 nm and far-red (FR) radiation within 700–750 nm can also excite photosystems, though not as efficiently as PAR. Vegetation and land surface models (LSMs) typically do not explicitly account for UV's contribution to energy budgets or photosynthesis, nor FR's contribution to photosynthesis. However, whether neglecting UV and FR has significant impacts remains unknown. Here, we explored how canopy radiative transfer (RT) and photosynthesis are impacted when explicitly implementing UV in the canopy RT model and accounting for UV and FR in the photosynthesis models within a next-generation LSM that can simulate hyperspectral canopy RT. We validated our improvements using photosynthesis measurements from plants under different light sources and intensities and surface reflection from an eddy-covariance tower. Our model simulations suggested that at the whole plant level, after accounting for UV and FR explicitly, chlorophyll content, leaf area index (LAI), clumping index, and solar radiation all impact the modeling of gross primary productivity (GPP). At the global scale, mean annual GPP within a grid would increase by up to 7.3% and the increase is proportional to LAI; globally integrated GPP increases by 4.6 PgC year −1 (3.8% of the GPP without accounting for UV + FR). Further, using PAR to proxy UV could overestimate surface albedo by more than 0.1, particularly in the boreal forests. Our results highlight the importance of improving UV and FR in canopy RT and photosynthesis modeling and the necessity to implement hyperspectral or multispectral canopy RT schemes in future vegetation and LSMs.

energy budget↗

COBRA:COMPUTED-TOMOGRAPHY BASED RANDOM-FIELD APPROXIMATION

SF-25-115 COBRA (COmputed-tomography Based Random-field Approximation) is a Python application for generating statistically equivalent random fields from CT-scan imagery. It leverages Karhunen–Loève expansions to model microstructural variability, enabling users to: Preprocess CT scans (filtering and Gaussian transformation); Fit covariance kernels fromempirical data; Solve eigenproblems to obtain KL modes; Sample random fields onsistent with fitted statistics; Postprocess samples back into the physical domain.

Hu, Tianchen↗