Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Covariance Data”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,153 records · Page 64

Modeling land surface processes over a mountainous rainforest in Costa Rica using CLM4.5 and CLM5

This study compares the performance of the Community Land Models (CLM4.5 and CLM5) against tower and ground measurements from a tropical montane rainforest in Costa Rica. The study site receives over 4000 mm of mean annual precipitation and has high daily levels of relative humidity. The measurement tower is equipped with eddy-covariance and vertical profile systems able to measure various micrometeorological variables, particularly in wet and complex terrain. In this work, results from point-scale simulations for both CLM4.5 and its updated version (CLM5) are compared to observed canopy flux and micrometeorological data. Both models failed to capture the effects of frequent rainfall events and mountainous topography on the variables of interest (temperatures, leaf wetness, and fluxes). Overall, CLM5 alleviates some errors in CLM4.5, but CLM5 still cannot precisely simulate a number of canopy processes for this forest. Soil, air, and canopy temperatures, as well as leaf wetness, remain too sensitive to incoming solar radiation rates despite updates to the model. As a result, daytime vapor flux and carbon flux are overestimated, and modeled temperature differences between day and night are higher than those observed. Slope effects appear in the measured average diurnal variations of surface albedo and carbon flux, but CLM5 cannot simulate these features. This study suggests that both CLMs still require further improvements concerning energy partitioning processes, such as leaf wetness process, photosynthesis model, and aerodynamic resistance model for wet and mountainous regions.

58 GEOSCIENCES↗

Determination of the auroral O/1S/ production sources from coordinated rocket and satellite measurements

This paper examines the auroral excitation mechanisms using data collected in a satellite/rocket/aurora coordination. Past theories considered that variation with altitude of the emission from the O(1S) state is the same as that from the N2(+) first negative band system. Actually, the source function for O(1S) does not covary with N2(+) but has a smaller scale height. The covariation is due to quenching of the source function by O(3P). The steadiness of the aurora in a flight of the Michigan Airglow Payload allows the profile to be differentiated for direct comparison with theory. It is concluded that when most of the relevant variables for production of the O(1S) state are simultaneously measured in aurora, the major source functions are dissociative recombination of O2(+) above 210 km to apogee at 250 km, and below 180 km the energy transfer from N2(A3sigma) in the oxygen atom quenching reaction.

Sharp, W. E.↗

Implementation of a Parallel Kalman Filter for Stratospheric Chemical Tracer Assimilation

A Kalman filter for the assimilation of long-lived atmospheric chemical constituents has been developed for two-dimensional transport models on isentropic surfaces over the globe. An important attribute of the Kalman filter is that it calculates error covariances of the constituent fields using the tracer dynamics. Consequently, the current Kalman-filter assimilation is a five-dimensional problem (coordinates of two points and time), and it can only be handled on computers with large memory and high floating point speed. In this paper, an implementation of the Kalman filter for distributed-memory, message-passing parallel computers is discussed. Two approaches were studied: an operator decomposition and a covariance decomposition. The latter was found to be more scalable than the former, and it possesses the property that the dynamical model does not need to be parallelized, which is of considerable practical advantage. This code is currently used to assimilate constituent data retrieved by limb sounders on the Upper Atmosphere Research Satellite. Tests of the code examined the variance transport and observability properties. Aspects of the parallel implementation, some timing results, and a brief discussion of the physical results will be presented.

Chang, Lang-Ping↗

Remote Sensing of Rain

The first problem addressed concerns passive-microwave rain retrievals. Most current approaches start by building off-line a cloud-model-derived database. Given data, the retrieval algorithms search the database for the microwave temperatures "closest" to the observed data, then after some fine-tuning (performed in different ways by different implementations) the rain is estimated to be that which corresponds to the selected (and fine-tuned) set of database temperatures. These approaches have three drawbacks: they cannot properly take into account the ambiguities which arise from the fact that several rain scenarios can produce the same observed temperatures; they are quite inefficient since they require manipulating a large database along with often complex "fine-tuning" procedures; and they cannot refine their estimates if additional data is available. This past year we have derived closed formulae relating observed microwave brightness temperatures, T(sub b), and the underlying rain rates, R: average T(sub b) =f (rain) and average rain = g (T(sub b)), along with the corresponding covariance matrices. These results are sufficient to describe the conditional probabilities p(R/T(sub b)) and p(T(sub b)/R) to second order. Progress has also been made towards deriving a robust description of the rain drop size distribution (DSD). The widespread approach consisting in parameterizing the DSD as a gamma-distribution in terms of the drop diameter D suffers from the facts that, in reality, the DSD is not a smooth function of D and that the largely arbitrary Gamma model imposes unintended behavior, which has implications on any quantities derived from the DSD model. We have therefore developed a non-parametric yet practical description of the DSD, which is particularly well-suited for use in remote-sensing applications. The diagram on the left shows a comparison between an actual DSD sample and the truncated non-parametric representation. One figure shows the relation between radar reflectivity and rain rate derived using this representation. Validation of the Tropical Rainfall Measuring Mission (TRMM) radar-radiometer combined R and DSD algorithm is underway. This algorithm was designed to make optimal use of the instantaneous reflectivity profiles measured by the TRMM radar and the microwave brightness temperatures measured by the TRMM passive radiometer. So far, it appears to be the most reliable TRMM rain algorithm.

Haddad, Ziad S.↗

On the Choice of Variable for Atmospheric Moisture Analysis

The implications of using different control variables for the analysis of moisture observations in a global atmospheric data assimilation system are investigated. A moisture analysis based on either mixing ratio or specific humidity is prone to large extrapolation errors, due to the high variability in space and time of these parameters and to the difficulties in modeling their error covariances. Using the logarithm of specific humidity does not alleviate these problems, and has the further disadvantage that very dry background estimates cannot be effectively corrected by observations. Relative humidity is a better choice from a statistical point of view, because this field is spatially and temporally more coherent and error statistics are therefore easier to obtain. If, however, the analysis is designed to preserve relative humidity in the absence of moisture observations, then the analyzed specific humidity field depends entirely on analyzed temperature changes. If the model has a cool bias in the stratosphere this will lead to an unstable accumulation of excess moisture there. A pseudo-relative humidity can be defined by scaling the mixing ratio by the background saturation mixing ratio. A univariate pseudo-relative humidity analysis will preserve the specific humidity field in the absence of moisture observations. A pseudorelative humidity analysis is shown to be equivalent to a mixing ratio analysis with flow-dependent covariances. In the presence of multivariate (temperature-moisture) observations it produces analyzed relative humidity values that are nearly identical to those produced by a relative humidity analysis. Based on a time series analysis of radiosonde observed-minus-background differences it appears to be more justifiable to neglect specific humidity-temperature correlations (in a univariate pseudo-relative humidity analysis) than to neglect relative humidity-temperature correlations (in a univariate relative humidity analysis). A pseudo-relative humidity analysis is easily implemented in an existing moisture analysis system, by simply scaling observed-minus background moisture residuals prior to solving the analysis equation, and rescaling the analyzed increments afterward.

Dee, Dick P.↗

Observations of the Interaction and Transport of Fine Mode Aerosols With Cloud and/or Fog in Northeast Asia From Aerosol Robotic Network and Satellite Remote Sensing

Analysis of Sun photometer measured and satellite retrieved aerosol optical depth (AOD) datahas shown that major aerosol pollution events with very highfine mode AOD (>1.0 in midvisible) in theChina/Korea/Japan region are often observed to be associated with significant cloud cover. This makesremote sensing of these events difficult even for high temporal resolution Sun photometer measurements.Possible physical mechanisms for these events that have high AOD include a combination of aerosolhumidification, cloud processing, and meteorological covariation with atmospheric stability andconvergence. The new development of Aerosol Robotic Network Version 3 Level 2 AOD with improved cloudscreening algorithms now allow for unprecedented ability to monitor these extremefine mode pollutionevents. Further, the spectral deconvolution algorithm (SDA) applied to Level 1 data (L1; no cloud screening)provides an even more comprehensive assessment offine mode AOD than L2 in current and previous dataversions. Studying the 2012 winter-summer period, comparisons of Aerosol Robotic Network L1 SDA dailyaveragefine mode AOD data showed that Moderate Resolution Imaging Spectroradiometer satellite remotesensing of AOD often did not retrieve and/or identify some of the highestfine mode AOD events in thisregion. Also, compared to models that include data assimilation of satellite retrieved AOD, the L1 SDAfinemode AOD was significantly higher in magnitude, particularly for the highest AOD events that were oftenassociated with significant cloudiness.

Eck, T. F.↗

High-Speed Schlieren Analysis of Retropropulsion Jet in Mach 10 Flow

High-speed schlieren imaging of a 5-in.-diam, 70° sphere-cone model with a single centerline 15° half-angle nozzle was performed in a Mach 10 flow. Image sequences captured at 100 kHz were obtained with the nozzle plugged (fully blunted model), no nozzle flow, and with nozzle flow over a range of supply pressures. The freestream unit Reynolds number was kept constant for all runs at3.38×106m−1. Time-average information pertaining to the location, width, and stand-off distance of the bow shock, interface, terminal shock, and triple point as a function of jet pressure ratio are presented based on analysis of covariance images from each run. Measurement of jet plume boundary shape as a function of jet pressure ratio is also presented. Analyses of frequency content associated with the jet structure are shown, as are results of the proper orthogonal decomposition of these data. For all runs with nozzle flow, a dominant fundamental frequency of 2 kHz was observed. Finally, a discussion of bow shock unsteadiness resulting from interaction with freestream disturbances for the blunted and no nozzle flow is provided, as is a discussion on the formation of weak waves emanating from the bow shock resulting from this interaction.

Brett F Bathel↗

Event-by-event neutron–photon multiplicity correlations in 252Cf(sf)

Excited nuclear fragments are emitted during nuclear fission. The de-excitation of these fission fragments takes place as sequential emission of neutrons followed by photons. A correlation between neutron and photon multiplicities accompanying fission is thus expected. Fission event generators based on established statistical nuclear physics models predict a negative event-by-event correlation in neutron-photon multiplicity. A survey of published experimental results of an event-by-event covariance between the neutron - photon multiplicities emitted following the spontaneous fission of 252 Cf is presented. Analytic unfolding expressions are developed in this work to determine the bias introduced by background sources, particle misclassification, pulse pileup, and inelastic photon production. The published experimental data are re-analyzed using these unfolding techniques and are found to be in qualitative agreement with the predictions of model-based calculations. Finally, in particular, we have concluded that there exists a significant event-by-event neutron-photon emission competition following the spontaneous fission of 252 Cf.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Physiological Controls on Carbon Fluxes and Biomass Production in Miscanthus: Insights From a Process‐Based Agroecosystem Model

Biomass crops serve as essential feedstocks for renewable energy and bioproducts and play a critical role in achieving lower emissions in the transportation sector. However, dedicated perennial biomass crops such as Miscanthus × giganteus (Miscanthus) remain underrepresented in process-based agroecosystem models, limiting robust evaluation of their economic and environmental performance. In this study, we developed a data-constrained representation of the sterile triploid Miscanthus (IL clone) within the process-based model ecosys, integrating global sensitivity analysis, ensemble simulation, and parameter calibration. Planting, harvesting, and fertilization practices consistent with field management were incorporated, and phenology was constrained using PhenoCam-derived Green Chromatic Coordinate (GCC) data. Using the Morris global sensitivity analysis method, we identified 11 key physiological parameters governing plant carbon, water, and nutrient relations, particularly processes associated with CO 2 assimilation. We then conducted ensemble simulations by perturbing these parameters and calibrated the model against eddy covariance fluxes and field-measured biomass. Building on the calibrated operating state, parameter-response analyses show that different photosynthetic processes influence productivity in different ways. Protein allocation determines whether productivity increases toward a higher level, whereas electron transport capacity controls additional gains once protein allocation approaches saturation. These findings demonstrate that parameter importance depends on physiological context and on which photosynthetic processes remain limiting. Calibration and validation against observations show that ecosys can reliably reproduce carbon and water fluxes, as well as both above- and belowground biomass, with post-calibration GPP R 2 improving from 0.67 to 0.95 during the calibration period and remaining high during validation (R 2 = 0.95). These results provide a mechanistic foundation for regional simulations and sustainable bioenergy assessments.

ecosys↗

AmeriFlux FLUXNET-1F US-Wkg Walnut Gulch Kendall Grasslands

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-Wkg Walnut Gulch Kendall Grasslands. This is the FLUXNET version of the carbon flux data for the site US-Wkg Walnut Gulch Kendall Grasslands produced by applying the standard ONEFlux (1F) software. Site Description - This site is located in a small, intensively-studied, experimental watershed within USDA-ARS's Walnut Gulch Experimental Watershed. Eddy covariance measurements of energy, water and CO2 fluxes began in the spring of 2004, though meteorological (including Bowen ratio) and hydrological measurements are available much further back.

Scott, Russell↗

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↗

Estimation of Kalman filter model parameters from an ensemble of tests

A methodology for estimating initial mean and covariance parameters in a Kalman filter model from an ensemble of nonidentical tests is presented. In addition, the problem of estimating time constants and process noise levels is addressed. Practical problems such as developing and validating inertial instrument error models from laboratory test data or developing error models of individual phases of a test are generally considered.

Gibbs, B. P.↗

Simultaneous stochastic inversion for geomagnetic main field and secular variation. I - A large-scale inverse problem

The method of stochastic inversion is extended to the simultaneous inversion of both main field and secular variation. In the present method, the time dependency is represented by an expansion in Legendre polynomials, resulting in a simple diagonal form for the a priori covariance matrix. The efficient preconditioned Broyden-Fletcher-Goldfarb-Shanno algorithm is used to solve the large system of equations resulting from expansion of the field spatially to spherical harmonic degree 14 and temporally to degree 8. Application of the method to observatory data spanning the 1900-1980 period results in a data fit of better than 30 nT, while providing temporally and spatially smoothly varying models of the magnetic field at the core-mantle boundary.

Bloxham, Jeremy↗

Optical and X-ray observations of the low-mass X-ray binary EXO 0748-676

Optical and X-ray observations of EXO-676 in late March 1989 are presented. Our optical observations provide some support for the previously observed correlation between the mean optical brightness and light curve morphology. Unexpectedly, the mean X-ray and optical flux levels during this period do not reflect similar system states. The optical counterpart is found to be in an intermediate to low state, while X-ray data imply a bright (high) state. The changed relationship between optical and X-ray fluxes is evidence showing that EXO 0748-676 has possibly evolved. We fail to find correlated variability in simultaneous X-ray and optical observations. The lack of covariability is attributed to the limited simultaneous coverage of the source and/or significant geometric modulation in the optical light curve.

Thomas, Brian↗

Uncertainty Quantification for a Multi-phase Carbon Equation of State Model

Many physics models have tunable parameters that are calibrated by matching the model output to experimental or calculated data. However, given that calibration data often contain uncertainty and that different model parameter sets might result in a very similar simulated output for a finite calibration data set, it is advantageous to provide an ensemble of parameter sets that are consistent with the calibration data. Uncertainty quantification (UQ) provides a means to generate such an ensemble in a statistically rigorous fashion. In this work, we perform UQ for a multi-phase equation of state (EOS) model for carbon containing the diamond, graphite, and liquid phases. We use a Bayesian framework for the UQ and introduce a novel strategy for including phase diagram information in the calibration. The method is highly general and accurately reproduces the calibration data without any material-specific prior knowledge of the EOS model parameters.

36 MATERIALS SCIENCE↗

Model-Based Detection Scheme for Spoofed Sensor Data in Grid-Connected Inverters

In this article, a novel method to detect spoofed sensor data (SSD) for grid-connected inverters is proposed. First, a Kalman filter is used to estimate the inverter status in the proposed method. Then, the impact of SSD on the residual between estimation and control reference is analyzed. By leveraging the standard error of residual, an adaptive cumulative sum (CUSUM) chart is developed to achieve fast SSD detection and improve robustness. In addition, the proposed SSD detection method is verified by the hardware experiment. A comparative experiment with a conventional method also demonstrates the rapidity and robustness of the proposed method. Lastly, the performance of the proposed method under different noise covariance matrices is analyzed. Comprehensive analysis shows that the proposed method meets the detection requirement in IEEE 1547 Standard.

42 ENGINEERING↗

Simultaneous inference of equation of state parameters and unknown data errors with uncertainty quantification via hierarchical Bayesian posterior maximization

Equations of state (EOSs) are a key component in running hydrodynamic simulations as they relate the thermodynamic states for the material. The Davis reactants EOS is commonly used for modeling high explosives (HEs), and the EOS model parameters are calibrated using material specific data. The calibrations are often performed with uncertainty quantification via Bayesian inference to account for uncertainty in the data and generate ensembles of likely parameters. However, there are relatively few HE data sets to use for calibration and many are historical and lack error information. In this work, we simultaneously calibrate the Davis reactants EOS model parameters and unknown data error terms for the high explosive PBX 9501. To quantify the uncertainty in the models and the data, we use a Bayesian framework for the calibration and compute the hierarchical Bayesian posterior distribution with both a posteriori maximization approach and Markov Chain Monte Carlo. In general, we find that, given our assumptions, the two approaches result in similar calibrated parameters, posterior covariance matrices, and insights about the parameters but that the posterior maximization requires far less computational resources.

97 MATHEMATICS AND COMPUTING↗

Giant Graviton Expansion from Bubbling Geometry: Discreteness from Quantized Geometry

The superconformal index of half-BPS states in N = 4 supersymmetric Yang-Mills with gauge group U⁡(N) admits an expansion in terms of giant gravitons, J N (q) = J ∞ ⁡(q)⁢Σ$^{∞}_{m=0}$ q m⁢N ⁢ J^ m ⁡(q), where m is the number of giant gravitons and J ∞ ⁡(q) is the graviton index. The expansion can be viewed as the implementation of trace relations for finite N. We derive this expansion directly in supergravity from the class of half-BPS solutions due to Lin, Lunin, and Maldacena in type IIB supergravity. The moduli space of these configurations can be quantized using covariant quantization methods. We show how this quantization leads to the precise expression for the expansion in terms of giant gravitons. Our proposal provides a derivation of the giant graviton expansion directly in terms of quantized supergravity degrees of freedom, and it recovers discrete data via quantum geometries that are classically nonsmooth.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗