Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analysis of Variance”

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,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Remote sensing of Earth Terrain

The Mueller matrix and polarization covariance matrix for polarimetric radar systems were studied. The clutter is modelled by a layer of random permittivity, described by a three dimensional correlation function, with variance, and horizontal and vertical correlation lengths. This model is applied, using the wave theory with Born approximations carried to the second order, to find the backscattering elements of the polarimetric matrices. It is found that 8 out of 16 elements of the Mueller matrix are identically zero, corresponding to a covariance matrix with four zero elements. Theoretical predictions are matched with experimental data for vegetation fields. A general mixing formula is derived for discrete scatterers immersed in a host medium. The inclusion particles are assumed to be ellipsoidal. The electric field inside the scatterers is determined by quasistatic analysis, assuming the diameter of the inclusion particles to be much smaller than the wavelength. The results are applicable to general multiphase mixtures.

Kong, J. A.↗
NASA NTRS2018

Variability in Galactic Cosmic Radiation- Induced DNA Damage Response in Inbred Mice Is Modulated by Genetics

In radiation biology, the ability to predict cancer risk associated with exposure to low doses of high-LET (Linear Energy Transfer) ionizing radiation remains a challenge. Epidemiological methods lack the sensitivity and power to provide detailed risk estimates for cancer and ignore individual sensitivity. We have hypothesized that DNA repair capacity is the primary factor differentiating peoples radiation sensitivity. We previously showed in immortalized human cell lines that characterizing the dose and time dependence of p53-binding protein 1 (53BP1) foci formation in the nucleus following X-rays exposure is sufficient to predict DNA repair response to any other LET in the same cell line. We now tested this hypothesis across a population of mice with different genetic background. Fibroblast cells were extracted and cultivated from 76 individual mice from 15 different strains and exposed to HZE (high (H) atomic number (Z) and energy (E) galactic cosmic ray particles) particles and X-rays. Individual radiation sensitivities were investigated by high throughput measurement of DNA repair kinetics that evaluated 53bp1 foci numbers as a surrogate for DNA double-strand breaks at various times post-irradiation. Instead of just counting foci which can be hard to distinguish for high-LET or high doses, we also took into account the track structure of high-LET particles to compute the remaining number of unrepaired tracks as a function of time post-irradiation. As expected, the percentage of unrepaired track over a 48 hours follow-up period increased with LET. In addition, repair rate was modulated by genetics, with animals from the same strain showing small variance while large rate differences were observed between strains. Radiation strain sensitivity ranking was estimated based on repair rates from exposure to each LET evaluated in this work, and ranking for high-LET correlated better with ranking from high dose of X-ray, not low dose. At the in-vivo level, drops in T-cells and B-cells number measured 24 hours after 0.1 Gy (Gray) X-ray exposure, correlated with slower DNA repair kinetic in fibroblast cells of the same strains of mice. At the genomic level, mouse genome wide association (GWA) analysis identified seven significant genetic loci on chromosomes 2, 3, 7, 10, 11, 13 and 19 with different significance depending on the LET. Interestingly, for the two highest LET, a common locus on Chromosome 10 was identified with high enrichment for DNA repair associated genes.Overall, this work suggests that repair kinetics of primary skin fibroblasts is a good surrogate marker for in-vivo radiation sensitivities in other tissues and that this response is modulated by genetics. Our study also confirms that DNA repair kinetics following high doses of X-ray can be used to predict radiation sensitivity to high-LET.

High-LET Radiation↗
DOE OSTI2015

A stochastic approach to quantifying the blur with uncertainty estimation for high-energy X-ray imaging systems

One of the primary causes of blur in a high-energy X-ray imaging system is the shape and extent of the radiation source, or ‘spot’. It is important to be able to quantify the size of the spot as it provides a lower bound on the recoverable resolution for a radiograph, and penumbral imaging methods – which involve the analysis of blur caused by a structured aperture – can be used to obtain the spot’s spatial profile. We present a Bayesian approach for estimating the spot shape that, unlike variational methods, is robust to the initial choice of parameters. The posterior is obtained from a normal likelihood, which was constructed from a weighted least squares approximation to a Poisson noise model, and prior assumptions that enforce both smoothness and non-negativity constraints. A Markov chain Monte Carlo algorithm is used to obtain samples from the target posterior, and the reconstruction and uncertainty estimates are the computed mean and variance of the samples, respectively. Lastly, synthetic data-sets are used to demonstrate accurate reconstruction, while real data taken with high-energy X-ray imaging systems are used to demonstrate applicability and feasibility.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Evapotranspiration simulations in ISIMIP2a—Evaluation of spatio-temporal characteristics with a comprehensive ensemble of independent datasets

Actual land evapotranspiration (ET) is a key component of the global hydrological cycle and an essential variable determining the evolution of hydrological extreme events under different climate change scenarios. However, recently available ET products show persistent uncertainties that are impeding a precise attribution of human-induced climate change. Here, we aim at comparing a range of independent global monthly land ET estimates with historical model simulations from the global water, agriculture, and biomes sectors participating in the second phase of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP2a). Among the independent estimates, we use the EartH2Observe Tier-1 dataset (E2O), two commonly used reanalyses, a pre-compiled ensemble product (LandFlux-EVAL), and an updated collection of recently published datasets that algorithmically derive ET from observations or observations-based estimates (diagnostic datasets). A cluster analysis is applied in order to identify spatiotemporal differences among all datasets and to thus identify factors that dominate overall uncertainties. The clustering is controlled by several factors including the model choice, the meteorological forcing used to drive the assessed models, the data category (models participating in the different sectors of ISIMIP2a, E2O models, diagnostic estimates, reanalysis-based estimates or composite products), the ET scheme, and the number of soil layers in the models. By using these factors to explain spatial and spatio-temporal variabilities in ET, we and that the model choice mostly dominates (24%{40% of variance explained), except for spatio-temporal patterns of total ET, where the forcing explains the largest fraction of the variance (29%). The most dominant clusters of datasets are further compared with individual diagnostic and reanalysis-based estimates to assess their representation of selected heat waves and droughts in the Great Plains, Central Europe and Western Russia. Although most of the ET estimates capture these extreme events, the generally large spread among the entire ensemble indicates substantial uncertainties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Impacts of the Pacific–Japan and Circumglobal Teleconnection Patterns on the Interdecadal Variability of the East Asian Summer Monsoon

Based on the Twentieth Century Reanalysis (20CR) dataset, the dominant modes of interdecadal variability of the East Asian summer monsoon (EASM) are investigated through a multivariate empirical orthogonal function analysis (MV-EOF). The first mode (EA1) is characterized by an anomalous cyclone centered over Taiwan and an anomalous anticyclone centered over the Bohai Sea. These phenomena are part of the meridional wave–like teleconnection pattern propagating poleward from the southern tropical western North Pacific (WNP), referred to as the interdecadal Pacific–Japan (PJ) pattern. The interdecadal PJ pattern is driven by negative anomalous convective heating over the southern tropical WNP, which is associated with the interdecadal Pacific oscillation (IPO) and the interdecadal Indian Ocean basin mode (IOBM). The amplitude of the EA1 and its contribution to the total variance of the EASM decrease remarkably after the 1960s. The second MV-EOF mode (EA2) is characterized by cyclone anomalies extending from northeastern China to Japan, which are part of a circumglobal wave train. Given the spatial scale of the wave train in the zonal direction (wavenumber 5), as well as the fact that it possesses barotropic structures and propagates along the Northern Hemispheric jet stream, it is referred to herein as the interdecadal circumglobal teleconnection (CGT) pattern. The interdecadal CGT pattern is associated with the forcing from the Atlantic multidecadal oscillation (AMO). Though the interdecadal PJ and CGT patterns are derived from the 20CR dataset, they are carefully verified through comparisons with various observational and reanalysis datasets from different perspectives.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1983

LANDSAT-D thematic mapper image dimensionality reduction and geometric correction accuracy

When principal component analysis of a subscene of a section of the Sacramento River showed lower correlation among the TM spectral components that were observed for the MSS spectral components, principal component analysis was applied to a LANDSAT 2 MSS subscene of the same area for comparison purposes. Correlation coefficient matrices indicate the pairwise similarity and correlation of the data for the spectral components. The principal components transformation matrix, indicates the weights applied to the original components to generate the transformed components. The first two TM components can be described as visible and near infrared. For the MSS data, the first transformed component is roughly the average of the four original components. The second transformed component is roughly the difference between the visible and infrared components. Tables show that 97.0% of the variance in an MSS image is contained in only two transformed components.

Ford, G. E.↗
NASA NTRS2011

Assessing the Urban Heat Island Effect Across Biomes in the Continental USA Using Landsat and MODIS

Impervious surface area (ISA) from the Landsat TM and land surface temperature (LST) from MODIS averaged over three annual cycles (2003-2005) are used in a spatial analysis to assess the urban heat island (UHI) skin temperature amplitude and its relationship to development intensity, size, and ecological setting for 38 of the most populous cities in the continental United States. Development intensity zones based on %ISA are defined across urban gradients and used to stratify sampling of LST and NDVI. We find that ecological context significantly influences the amplitude of summer daytime UHI (urban - rural temperature difference) with the largest 8 C (average) for cities built in mixed forest biomes. For all cities ISA is the primary driver for increase in temperature explaining 70% of the total variance. Annually, urban areas are warmer than the non-urban fringe by 2.9 C, except in biomes with arid and semiarid climates. The average amplitude of the UHI is asymmetric with a 4.3 C difference in summer and 1.3 C in winter. In desert environments, UHI's point to a possible heat sink effect. Results show that the urban heat island amplitude increases with city size and is seasonally asymmetric for a large number of cities across most biomes. The implications are that for urban areas developed within forested ecosystems the summertime UHI can be quite high relative to the wintertime UHI suggesting that the residential energy consumption required for summer cooling is likely to increase with urban growth within those biomes.

Imhoff, Marc L.↗
DOE OSTI2015

Inconsistent Investment and Consumption Problems

In a traditional Black–Scholes market we develop a verification theorem for a general class of investment and consumption problems where the standard dynamic programming principle does not hold. The theorem is an extension of the standard Hamilton–Jacobi–Bellman equation in the form of a system of non-linear differential equations. We derive the optimal investment and consumption strategy for a mean-variance investor without pre-commitment endowed with labor income. In the case of constant risk aversion it turns out that the optimal amount of money to invest in stocks is independent of wealth. The optimal consumption strategy is given as a deterministic bang-bang strategy. In order to have a more realistic model we allow the risk aversion to be time and state dependent. Of special interest is the case were the risk aversion is inversely proportional to present wealth plus the financial value of future labor income net of consumption. Using the verification theorem we give a detailed analysis of this problem. It turns out that the optimal amount of money to invest in stocks is given by a linear function of wealth plus the financial value of future labor income net of consumption. The optimal consumption strategy is again given as a deterministic bang-bang strategy. We also calculate, for a general time and state dependent risk aversion function, the optimal investment and consumption strategy for a mean-standard deviation investor without pre-commitment. In that case, it turns out that it is optimal to take no risk at all.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Statistical Uncertainty in Paleoclimate Proxy Reconstructions

A quantitative analysis of any environment older than the instrumental record relies on proxies. Uncertainties associated with proxy reconstructions are often underestimated, which can lead to artificial conflict between different proxies, and between data and models. In this paper, using ordinary least squares linear regression as a common example, we describe a simple, robust and generalizable method for quantifying uncertainty in proxy reconstructions. We highlight the primary controls on the magnitude of uncertainty, and compare this simple estimate to equivalent estimates from Bayesian, nonparametric and fiducial statistical frameworks. We discuss when it may be possible to reduce uncertainties, and conclude that the unexplained variance in the calibration must always feature in the uncertainty in the reconstruction. This directs future research toward explaining as much of the variance in the calibration data as possible. We also advocate for a “data-forward” approach, that clearly decouples the presentation of proxy data from plausible environmental inferences.

58 GEOSCIENCES↗
DOE OSTI2013

Differentiation of Microbial Species and Strains in Coculture Biofilms by Multivariate Analysis of Laser Desorption Postionization Mass Spectra

7.87 to 10.5 eV vacuum ultraviolet (VUV) photon energies were used in laser desorption postionization mass spectrometry (LDPI-MS) to analyze biofilms comprised of binary cultures of interacting microorganisms. The effect of photon energy was examined using both tunable synchrotron and laser sources of VUV radiation. Principal components analysis (PCA) was applied to the MS data to differentiate species in Escherichia coli-Saccharomyces cerevisiae coculture biofilms. PCA of LDPI-MS also differentiated individual E. coli strains in a biofilm comprised of two interacting gene deletion strains, even though these strains differed from the wild type K-12 strain by no more than four gene deletions each out of approximately 2000 genes. PCA treatment of 7.87 eV LDPI-MS data separated the E. coli strains into three distinct groups two ?pure? groups and a mixed region. Furthermore, the ?pure? regions of the E. coli cocultures showed greater variance by PCA when analyzed by 7.87 eV photon energies than by 10.5 eV radiation. Comparison of the 7.87 and 10.5 eV data is consistent with the expectation that the lower photon energy selects a subset of low ionization energy analytes while 10.5 eV is more inclusive, detecting a wider range of analytes. These two VUV photon energies therefore give different spreads via PCA and their respective use in LDPI-MS constitute an additional experimental parameter to differentiate strains and species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

SPT-3G D1: Compton-$y$ maps using data from the SPT-3G and Planck surveys

We present thermal Sunyaev-Zel'dovich (tSZ) Compton-$y$ parameter maps constructed from two years (2019-2020) of observations with the South Pole Telescope (SPT) third-generation camera, SPT-3G, combined with data from the Planck satellite. Using a linear combination (LC) pipeline, we obtain a suite of reconstructions that explore different trade-offs between statistical sensitivity and suppression of astrophysical contaminants, including minimum-variance, CMB-deprojected, and CIB-deprojected $y$-maps. We validate these maps through different statistical techniques such as auto- and cross-power spectra with large-scale structure tracers as well as stacking on cluster locations. These tests are used to understand the balance between noise and astrophysical foreground residuals (such as the CIB) in combination with the recovery of the tSZ signal for different maps. For example, results from stacking at the location of clusters confirm the robustness of the recovered tSZ signal over the $\sim 1500\: {\rm deg}^2$ SPT-3G survey field used in this analysis. The high-resolution and low-noise maps produced here provide an important cosmological tool for future studies, including measurements of the Compton-$y$ map power spectrum, cross-correlations with other tracers of the large-scale structure, detailed modeling of cluster pressure profiles, and study of the thermodynamic state of the baryons in the Universe.

Maniyar, A. S. [Harvard-Smithsonian Ctr. Astrophys↗
NASA NTRS1985

Probabilistic Structural Analysis of SSME Turbopump Blades: Probabilistic Geometry Effects

A probabilistic study was initiated to evaluate the precisions of the geometric and material properties tolerances on the structural response of turbopump blades. To complete this study, a number of important probabilistic variables were identified which are conceived to affect the structural response of the blade. In addition, a methodology was developed to statistically quantify the influence of these probabilistic variables in an optimized way. The identified variables include random geometric and material properties perturbations, different loadings and a probabilistic combination of these loadings. Influences of these probabilistic variables are planned to be quantified by evaluating the blade structural response. Studies of the geometric perturbations were conducted for a flat plate geometry as well as for a space shuttle main engine blade geometry using a special purpose code which uses the finite element approach. Analyses indicate that the variances of the perturbations about given mean values have significant influence on the response.

Nagpal, V. K.↗
NASA NTRS2010

A Framework for Propagation of Uncertainties in the Kepler Data Analysis Pipeline

The Kepler space telescope is designed to detect Earth-like planets around Sun-like stars using transit photometry by simultaneously observing 100,000 stellar targets nearly continuously over a three and a half year period. The 96-megapixel focal plane consists of 42 charge-coupled devices (CCD) each containing two 1024 x 1100 pixel arrays. Cross-correlations between calibrated pixels are introduced by common calibrations performed on each CCD requiring downstream data products access to the calibrated pixel covariance matrix in order to properly estimate uncertainties. The prohibitively large covariance matrices corresponding to the ~75,000 calibrated pixels per CCD preclude calculating and storing the covariance in standard lock-step fashion. We present a novel framework used to implement standard propagation of uncertainties (POU) in the Kepler Science Operations Center (SOC) data processing pipeline. The POU framework captures the variance of the raw pixel data and the kernel of each subsequent calibration transformation allowing the full covariance matrix of any subset of calibrated pixels to be recalled on-the-fly at any step in the calibration process. Singular value decomposition (SVD) is used to compress and low-pass filter the raw uncertainty data as well as any data dependent kernels. The combination of POU framework and SVD compression provide downstream consumers of the calibrated pixel data access to the full covariance matrix of any subset of the calibrated pixels traceable to pixel level measurement uncertainties without having to store, retrieve and operate on prohibitively large covariance matrices. We describe the POU Framework and SVD compression scheme and its implementation in the Kepler SOC pipeline.

Clarke, Bruce D.↗
DOE OSTI2019

Methane Production Pathway Regulated Proximally by Substrate Availability and Distally by Temperature in a High-Latitude Mire Complex

Projected 21st century changes in high-latitude climate are expected to have significant impacts on permafrost thaw, which could cause substantial increases in emissions to the atmosphere of carbon dioxide (CO 2 ) and methane (CH 4 , which has a global warming potential 28 times larger than CO 2 over a 100-year horizon). However, predicted CH 4 emission rates are very uncertain due to difficulties in modeling complex interactions among hydrological, thermal, biogeochemical, and plant processes. Methanogenic production pathways (i.e., acetoclastic [AM] and hydrogenotrophic [HM]) and the magnitude of CH 4 emissions may both change as permafrost thaws, but a mechanistic analysis of controls on such shifts in CH 4 dynamics is lacking. In this study, we reproduced observed shifts in CH 4 emissions and production pathways with a comprehensive biogeochemical model (ecosys) at the Stordalen Mire in subarctic Sweden. Our results demonstrate that soil temperature changes differently affect AM and HM substrate availability, which regulates magnitudes of AM, HM, and thereby net CH 4 emissions. We predict very large landscape-scale, vertical, and temporal variations in the modeled HM fraction, highlighting that measurement strategies for metrics that compare CH 4 production pathways could benefit from model informed scale of temporal and spatial variance. Finally, our findings suggest that the warming and wetting trends projected in northern peatlands could enhance peatland AM fraction and CH 4 emissions even without further permafrost degradation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

SDSS-IV MaNGA: Exploring the Local Scaling Relations for N/O

We present, for the first time, the relationship between local stellar mass surface density, Σ * , and N/O derived from SDSS-IV MaNGA data, using a sample of 792,765 high signal-to-noise ratio star-forming spaxels. Using a combination of phenomenological modeling and partial correlation analysis, we find that Σ * alone is insufficient to predict the N/O in MaNGA spaxels and that there is an additional dependence on the local star formation rate surface density, Σ SFR . This effect is a factor of 3 stronger than the dependence of 12+log(O/H) on Σ SFR . Surprisingly, we find that the local N/O scaling relations also depend on the total galaxy stellar mass at fixed Σ * and the galaxy size at fixed stellar mass. We find that more compact galaxies are more nitrogen rich, even when Σ * and Σ SFR are controlled for. We show that ~50% of the variance of N/O is explained by the total stellar mass and size. Thus, the evolution of nitrogen in galaxies is set by more than just local effects and does not simply track the buildup of oxygen in galaxies. The precise form of the N/O–O/H relation is therefore sensitive to the sample of galaxies from which it is derived. This result casts doubt on the universal applicability of nitrogen-based strong-line metallicity indicators derived in the local universe.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

In-flight radiometric calibration of the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS)

A reflectance-based method was used to provide an analysis of the in-flight radiometric performance of AVIRIS. Field spectral reflectance measurements of the surface and extinction measurements of the atmosphere using solar radiation were used as input to atmospheric radiative transfer calculations. Five separate codes were used in the analysis. Four include multiple scattering, and the computed radiances from these for flight conditions were in good agreement. Code-generated radiances were compared with AVIRIS-predicted radiances based on two laboratory calibrations (pre- and post-season of flight) for a uniform highly reflecting natural dry lake target. For one spectrometer (C), the pre- and post-season calibration factors were found to give identical results, and to be in agreement with the atmospheric models that include multiple scattering. This positive result validates the field and laboratory calibration technique. Results for the other spectrometers (A, B and D) were widely at variance with the models no matter which calibration factors were used. Potential causes of these discrepancies are discussed.

Conel, James E.↗
DOE OSTI2023

BEYONDPLANCK IV. On end-to-end simulations in CMB analysis — Bayesian versus frequentist statistics

End-to-end simulations play a key role in the analysis of any high-sensitivity cosmic microwave background (CMB) experiment, providing high-fidelity systematic error propagation capabilities that are unmatched by any other means. In this paper, we address an important issue regarding such simulations, namely, how to define the inputs in terms of sky model and instrument parameters. These may either be taken as a constrained realization derived from the data or as a random realization independent from the data. We refer to these as posterior and prior simulations, respectively. We show that the two options lead to significantly different correlation structures, as prior simulations (contrary to posterior simulations) effectively include cosmic variance, but they exclude realization-specific correlations from non-linear degeneracies. Consequently, they quantify fundamentally different types of uncertainties. We argue that as a result, they also have different and complementary scientific uses, even if this dichotomy is not absolute. In particular, posterior simulations are in general more convenient for parameter estimation studies, while prior simulations are generally more convenient for model testing. Before BEYONDPLANCK, most pipelines used a mix of constrained and random inputs and applied the same hybrid simulations for all applications, even though the statistical justification for this is not always evident. BEYONDPLANCK represents the first end-to-end CMB simulation framework that is able to generate both types of simulations and these new capabilities have brought this topic to the forefront. The BEYONDPLANCK posterior simulations and their uses are described extensively in a suite of companion papers. In this work, we consider one important applications of the corresponding prior simulations, namely, code validation. Specifically, we generated a set of one-year LFI 30 GHz prior simulations with known inputs and we used these to validate the core low-level BEYONDPLANCK algorithms dealing with gain estimation, correlated noise estimation, and mapmaking.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1985

Mesospheric Winds and Tides at Arecibo/puerto Rico (18 Deg N, 67 Deg W) November 20-23, 1981

The SOUSY-VHF-Radar was operated in the period November 20 to 23, 1981 during a continuous run at the Arecibo Observatory. The data used for this analysis were recorded from 06 AST on November 20 until 06 AST on November 23. Meridional winds were measured at odd hours and zonal winds at even hours. The occurrence of mesospheric VHF radar echoes was limited as usual to the daylight hours between 06 and 18 AST. The height limits, were variable from day to day, they were roughly 62 to 90 km. The intermittency of echoes in time and height, as well as radio interference, made it fairly difficult to deduce tidal information for this time period. About 10 to 45 single velocity estimates could be used to calculate an hourly mean. The variance of the displayed resulting in data, obtained by weighted averaging over the 3 day period, is about similar to the mean. The winds were extremely low during the observation period. Velocities of the prevailing wind, observed during an earlier campaign were 5 to 10 times larger than in the period of November 20 to 23, 1981. The amplitudes of the diurnal tide are comparable.

Rottger, J.↗