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 703 records · Page 39

Landsat-TM identification of Amblyomma variegatum (Acari: Ixodidae) habitats in Guadeloupe

The feasibility of identifying specific habitats of the African bont tick, Amblyomma variegatum, from Landsat-TM images was investigated by comparing remotely sensed images of visible farms in Grande Terre (Guadeloupe) with field observations made in the same period of time (1986-1987). The different tick habitates could be separated using principal component analysis. The analysis clustered the sites by large and small variance of band values, and by vegetation and moisture indexes. It was found that herds in heterogeneous sites with large variances had more ticks than those in homogeneous or low variance sites. Within the heterogeneous sites, those with high vegetation and moisture indexes had more ticks than those with low values.

Hugh-Jones, M.↗

Fiber Bundle Morphisms as a Framework for Modeling Many-to-Many Maps

While it is not generally reflected in the `nice' datasets used for benchmarking machine learning algorithms, the real-world is full of processes that would be best described as many-to-many. That is, a single input can potentially yield many different outputs (whether due to noise, imperfect measurement, or intrinsic stochasticity in the process) and many different inputs can yield the same output (that is, the map is not injective). For example, imagine a sentiment analysis task where, due to linguistic ambiguity, a single statement can have a range of different sentiment interpretations while at the same time many distinct statements can represent the same sentiment. When modeling such a multivalued function $f: X \rightarrow Y$, it is frequently useful to be able to model the distribution on $f(x)$ for specific input $x$ as well as the distribution on fiber $f^{-1}(y)$ for specific output $y$. Such an analysis helps the user (i) better understand the variance intrinsic to the process they are studying and (ii) understand the range of specific input $x$ that can be used to achieve output $y$. Following existing work which used a fiber bundle framework to better model many-to-one processes we describe how morphisms of fiber bundles provide a template for building models which naturally capture the structure of many-to-many processes.

Coda, Elizabeth D.↗

The South Pole Telescope AGN Monitoring Campaign: First Release of SPTpol Bright AGN Light Curves

The South Pole Telescope (SPT) collaboration has recently embarked upon a campaign to monitor the brightness of a sample of active galactic nuclei (AGN), both in real time and in archival SPT data. The original design of the SPT was optimized for observations of the cosmic microwave background (CMB) at arc-minute and larger angular scales, and it has been used for this purpose for nearly twenty years, using three generations of CMB cameras. Recently it has been recognized that data from CMB experiments have the potential to be used for AGN monitoring. In this paper, we present the first public release of data from a full sample of SPT-monitored AGN, comprising 158 AGN light curves and associated data from the SPTpol camera, which was operational from 2012-2016. These light curves were created using observations from the SPTpol 500 deg$^{2}$ survey, in which the instrument was used to scan a 500 deg$^2$ patch of the sky several times per day with detectors sensitive to radiation in bands centered at 90 and 150 GHz. We provide a comprehensive description of the observations, the data processing methods, and the resulting light curve catalog. As an example of analyses that these data enable, we searched for a correlation between variability and spectral index, and we looked for ``bluer-when-brighter'' trends in the sample. Our analysis finds $> 10 σ$ correlation between fractional intrinsic variance and mean spectral index in the sample, but no significant evidence for bluer-when-brighter trends. The datasets from this study can be accessed through the SPT Treasury Record of AGN With Historical Activity and Time-Series or STRAWHAT catalog. This initial data release includes SPTpol light curves at 90 and 150 GHz, focusing on total intensity. In later updates, SPTpol polarization data and new observations from the SPT-3G instrument at 90, 150, and 220 GHz will be included.

Hood, J.C., II [Chicago U., KICP; Chicago U., Astr↗

A variance-based decomposition and global sensitivity index method for uncertainty quantification: Application to retrieved ice cloud properties

This study develops a novel uncertainty quantification (UQ) method for cloud microphysical property retrievals using variance-based decomposition and global sensitivity index. In this UQ framework, empirical orthogonal function (EOF) analysis is applied to the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) ground-based observations, which are the inputs for the cloud retrieval studied here. The principal components (PCs) in the EOF expansion are parameterized as random input variables, and hence, the input dimension is greatly reduced (up to a factor of 50), allowing large ensemble of random samplings. The EOF expansion improves the accuracy of the uncertainty estimation by taking into account the cross correlations in the input data profiles. Here, this method enables a probabilistic representation of a retrieval process by adding normally distributed perturbations into PCs of sample means of input data profiles within a time window. Therefore, it effectively facilitates objective validation of climate models against cloud retrievals under a probabilistic framework for rigorous statistical inferences. Moreover, the variance-based global sensitivity index analysis, part of this method, attributes the output uncertainties to each individual source, thus providing directions for improving retrieval algorithms and observation strategies. For demonstration, we apply this method to quantify the uncertainties of the ARM program's baseline cloud retrieval algorithm for an ice cloud case observed at the Southern Great Plains site on 9 March 2000.

54 ENVIRONMENTAL SCIENCES↗

Expendable Launch Vehicle Studies

Analytical support studies of expendable launch vehicles concentrates on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane linearized dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed. The design of a robust LQR controller based on the reduced order system is accomplished using the parameter perturbation technique. The ELV modeling and analysis team has spent the past three years working on the theoretical development and application of sensitivity analysis to solve eigenvalue problems associated with structural dynamics. Specific application areas include stochastic vibrations, viscously damped vibrations and nonlinear dynamics. In stochastic linear vibrations, sensitivity analysis methods have been developed which determined the expectation and variance of the response for given probability density functions of various stochastic parameters. Examples worked include beam vibrations with up to six stochastic parameters; and the eigensensitivity results differ little from those obtained, with much effect, from Monte Carlo techniques. For viscously damped vibrations, eigensensitivity analysis has given excellent results for homogeneous beams, modeled by either quadratic or quartic eigenvalue equations. Specific applications include Kevin and Maxwell-type viscoelastic beams.

Bainum, P. M.↗

Impact of three-dimensional hot-spot flow asymmetry on ion-temperature measurements in inertial confinement fusion experiments

Three-dimensional (3-D) implosion asymmetries lead to signficant variations in ion-temperature measurements in inertial confinement fusion experiments. We present an analytical method to generalize the physical properties of velocity variance in the Brysk ion-temperature model. This analysis provides a consistent explanation for the 3-D effects of inferred ion-temperature variations for various single modes and multimodes modeled by the deceleration-phase hydrocode DEC3D and the neutron transport code IRIS3D . The effect of the hot-spot ow asymmetry on variations in ion-temperature measurement is shown uniquely determined by a complete set of six hot-spot flow parameters. An approximated solution to the minimum inferred ion temperature is derived and shown to be close to the thermal ion temperature for low mode l = 1, which exhibits the largest anisotropic velocity variance in the single-mode spectrum. In conclusion, the isotropic velocity variance for low mode l = 2 is shown to result in the minimum inferred ion temperatures being well above the thermal ion temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Objective determination of image end-members in spectral mixture analysis

Spectral mixture analysis was shown to be a powerful, multifaceted tool for analysis of multi- and hyper-spectral data. The essence of the first phase of the approach is to determine a set of image end-members that best account for the spectral variance in an image cube within a constrained, linear least squares mixing model. The selection of the image end-member is usually achieved using a priori knowledge and successive trial and error solutions to refine the total number and physical location of the end-members. However, in many situations a more objective method of determining these essential components is desired. The problem of image end-member determination was approached objectively by using the inherent variance of the data. Unlike purely statistical methods such as factor analysis, this approach derives solutions that conform to a physically realistic model.

Tompkins, Stefanie↗

A Multiscale Approach to High-Resolution Ocean Profile Observations within a 4DVAR Analysis System

Most ocean data assimilation systems are tuned to process and assimilate observations to constrain features on the order of the mesoscale and larger. Typically this involves removal of observations or computing averaged observations. This procedure, while necessary, eliminates many observations from the analysis step and can reduce the overall effectiveness of a particular observing platform. Simply including these observations is not an option as doing so can produce an overdetermined, ill-conditioned problem that is more difficult to solve. An approach, presented here, aims to avoid such issues while at the same time increasing the number of observations within the assimilation. A two-step assimilation procedure with the four-dimensional variational data assimilation (4DVAR) system is adopted. The first step attempts to constrain the large-scale features by assimilating a set of super observations with appropriate background error correlation scales and error variances. The second step then attempts to correct smaller-scale features by assimilating the full observation set with shorter background error correlation scales and appropriate error variances; here the background state is taken as the analysis from the first step. Results using a real high-density observation set from underwater gliders in the region southeast of Iceland, collected during the 2017 Nordic Recognized Environmental Picture (NREP) experiment, will be shown using the Navy Coastal Ocean Model 4DVAR (NCOM-4DVAR).

54 ENVIRONMENTAL SCIENCES↗

Probabilistic Parameter Uncertainty Analysis of Single Input Single Output Control Systems

The current standards for handling uncertainty in control systems use interval bounds for definition of the uncertain parameters. This approach gives no information about the likelihood of system performance, but simply gives the response bounds. When used in design, current methods of m-analysis and can lead to overly conservative controller design. With these methods, worst case conditions are weighted equally with the most likely conditions. This research explores a unique approach for probabilistic analysis of control systems. Current reliability methods are examined showing the strong areas of each in handling probability. A hybrid method is developed using these reliability tools for efficiently propagating probabilistic uncertainty through classical control analysis problems. The method developed is applied to classical response analysis as well as analysis methods that explore the effects of the uncertain parameters on stability and performance metrics. The benefits of using this hybrid approach for calculating the mean and variance of responses cumulative distribution functions are shown. Results of the probabilistic analysis of a missile pitch control system, and a non-collocated mass spring system, show the added information provided by this hybrid analysis.

Smith, Brett A.↗

An evaluation with the Fourier Amplitude Sensitivity Test (FAST) of which land-surface parameters are of greatest importance in atmospheric modeling

Land-surface parameterizations based on a statistical-dynamical have been suggested recently to improve the representation of the surface forcing from heterogeneous land in atmospheric models. With this approach, land-surface characteristics are prescribed by probability density functions (PDFs) rather than single 'representative' values as in 'big-leaf' parameterizations. Yet the use of many PDFs results in an increased computational burden and requires the complex problem of representing covariances between PDFs to be addressed. In this study, a sensitivity analysis of a land-surface parameterization for atmospheric modeling was performed to evaluate the surface parameters most important to the variability of surface heat fluxes. The Fourier amplitude sensitivity test (FAST) used for this analysis determines the relative contribution of individual input parameters to the variance of energy fluxes resulting from a heterogeneous surface. By simultaneously varying all parameters according to their individual probability density functions, the number of computations needed is very much reduced by this technique. This analysis demonstrates that most of the variability of surface heat fluxes may be described by the distributions of relative stomatal conductance and surface roughness. Thus, the statistical-dynamical approach may be simplified by the use of only these two probability density functions.

Collins, Dan C.↗

Sensitivity Analysis of OECD Benchmark Tests in BISON

This report summarizes a NEAMS (Nuclear Energy Advanced Modeling and Simulation) project focused on sensitivity analysis of a fuels performance benchmark problem. The benchmark problem was defined by the Uncertainty Analysis in Modeling working group of the Nuclear Science Committee, part of the Nuclear Energy Agency of the Organization for Economic Cooperation and Development (OECD ). The benchmark problem involv ed steady - state behavior of a fuel pin in a Pressurized Water Reactor (PWR). The problem was created in the BISON Fuels Performance code. Dakota was used to generate and analyze 300 samples of 17 input parameters defining core boundary conditions, manuf acturing tolerances , and fuel properties. There were 24 responses of interest, including fuel centerline temperatures at a variety of locations and burnup levels, fission gas released, axial elongation of the fuel pin, etc. Pearson and Spearman correlatio n coefficients and Sobol' variance - based indices were used to perform the sensitivity analysis. This report summarizes the process and presents results from this study.

97 MATHEMATICS AND COMPUTING↗

Global sensitivity analysis in stochastic simulators of uncertain reaction networks

Stochastic models of chemical systems are often subjected to uncertainties in kinetic parameters in addition to the inherent random nature of their dynamics. Uncertainty quantification in such systems is generally achieved by means of sensitivity analyses in which one characterizes the variability with the uncertain kinetic parameters of the first statistical moments of model predictions. In this work, we propose an original global sensitivity analysis method where the parametric and inherent variability sources are both treated through Sobol’s decomposition of the variance into contributions from arbitrary subset of uncertain parameters and stochastic reaction channels. The conceptual development only assumes that the inherent and parametric sources are independent, and considers the Poisson processes in the random-time-change representation of the state dynamics as the fundamental objects governing the inherent stochasticity. Here, a sampling algorithm is proposed to perform the global sensitivity analysis, and to estimate the partial variances and sensitivity indices characterizing the importance of the various sources of variability and their interactions. The birth-death and Schlögl models are used to illustrate both the implementation of the algorithm and the richness of the proposed analysis method. The output of the proposed sensitivity analysis is also contrasted with a local derivative-based sensitivity analysis method classically used for this type of systems.

97 MATHEMATICS AND COMPUTING↗

Major to trace element imaging and analysis of iron age glasses using stage scanning in the analytical dual beam microscope (tandem)

Dark and clear silicate glasses formed during an iron age vitrification event ≈ 1500 years ago at the Broborg hillfort near Uppsala, Sweden have been analyzed using a scanning electron microscope equipped with a micro-X-ray fluorescence (μXRF) spectrometer. Correlated µXRF and electron beam-induced energy dispersive spectrometry (EDS) X-ray maps were collected via stage-scanning at constant velocity. This coupled procedure represents a new approach for the cultural heritage community to conduct analytical studies of archaeometric specimens composed of metal, ceramic, or mixed inorganic/organic materials, where major and trace element compositions are registered in space for areas up to the centimeter-length scale at micrometer-scale resolution. Overview images were used to select areas for EDS beam scan maps correlated with multispectral cathodoluminescence (CL) imaging and co-located quantitative EDS and μXRF point analysis. Fe, Ca, Mg, Ti, P, Mn, Zr, Zn, and Y are enriched in the dark glass, while Si, Al, K, Na, Ba, Sr, Rb, and Ga are enriched in the clear glass. Unmelted material is comprised predominately of quartz (SiO 2 ) along with trace apatite (Ca 5 (PO 4 ) 3 [Cl,OH]) and zircon (ZrSiO 4 ). Multivariate statistical analysis was used to measure the area fractions of high variance components while lower variance components represented phase mixtures. Differences between calculated melt viscosities for the glass compositions are consistent with field and laboratory observations. Coupled large area EDS and μXRF imaging shows significant promise for informed selection of higher spatial resolution and higher sensitivity follow-up studies, e.g., those performed using synchrotron analysis.

36 MATERIALS SCIENCE↗

Sensitivity analysis of thermal contact conductance modeling to inform MiniFuel irradiation capsule designs

The MiniFuel irradiation platform has been developed by Oak Ridge National Laboratory as a flexible, high-throughput separate effects testing capability within the High Flux Isotope Reactor (HFIR). Finite element thermal models are relied upon to design MiniFuel experiments to achieve a specific time-averaged irradiation temperature for experimental objectives. A previous study identified that uncertainty in the component heat generation rates and thermal contact conductance (TCC) model are the most significant contributors to predicted fuel temperature variance. To address both sources of uncertainty, this work performs sensitivity analysis on the TCC model to identify high-impact, high-uncertainty parameters that contribute to fuel temperature variance. The TCC model is analyzed in increasing detail, first using a standalone Python code, then again after coupling Python to the BISON fuel performance code. Furthermore, the parameters with the largest contributions to fuel temperature variance which can be reduced through design changes are identified as the initial subcapsule gas pressure, contact pressure between the fuel and dish, and the effective surface roughness of the interface. A set of design recommendations for future capsule designs has been established and applied to reduce the previously quantified average fuel temperature uncertainty ranges of ± 40 °C in the HFIR vertical experiment facilities (VXF) and ± 80 °C in the removable beryllium (RB) reflector to approximately ± 32 °C and ± 53 °C, respectively. This equates to a 21 % and 33 % reduction in the uncertainty range of the average fuel temperature for VXF and RB, respectively.

BISON↗

Vertically Resolved Analysis of the Madden-Julian Oscillation Highlights the Role of Convective Transport of Moist Static Energy

We simulate the Madden-Julian oscillation (MJO) over an aquaplanet with uniform surface temperature using the multiscale modeling framework (MMF) configuration of the Energy Exascale Earth System Model (E3SM-MMF). The model produces MJO-like features that have a similar spatial structure and propagation behavior to the observed MJO. To explore the processes involved in the propagation and maintenance of these MJO-like features, we perform a vertically resolved moist static energy (MSE) analysis for the MJO (Yao et al., 2022, https://doi.org/10.1175/jas-d-20-0254.1). Unlike the column-integrated MSE analysis, our method emphasizes the local production of MSE variance and quantifies how individual physical processes amplify and propagate the MJO's characteristic vertical structure. We find that radiation, convection, and boundary layer (BL) processes all contribute to maintaining the MJO, balanced by the large-scale MSE transport. Furthermore, large-scale dynamics, convection, and BL processes all contribute to the propagation of the MJO, while radiation slows the propagation. Additionally, we perform mechanism-denial experiments to examine the role of radiation and associated feedbacks in simulating the MJO. We find that the MJO can still self-emerge and maintain its characteristic structures without radiative feedbacks. This study highlights the role of convective MSE transport in the MJO dynamics, which was overlooked in the column-integrated MSE analysis.

54 ENVIRONMENTAL SCIENCES↗

Polarization bispectrum for measuring primordial magnetic fields

We examine the potential of polarization bispectra of the cosmic microwave background (CMB) to constrain primordial magnetic fields (PMFs). We compute all possible bispectra between temperature and polarization anisotropies sourced by PMFs and show that they are weakly correlated with well-known local-type and secondary ISW-lensing bispectra. From a Fisher analysis it is found that, owing to E-mode bispectra, in a cosmic-variance-limited experiment the expected uncertainty in the amplitude of magnetized bispectra is 80% improved in comparison with an analysis in terms of temperature auto-bispectrum alone. In the Planck or the proposed PRISM experiment cases, we will be able to measure PMFs with strength 2.6 or 2.2 nG. PMFs also generate bispectra involving B-mode polarization, due to tensor-mode dependence. We also find that the B-mode bispectrum can reduce the uncertainty more drastically and hence PMFs comparable to or less than 1 nG may be measured in a PRISM-like experiment.

79 ASTRONOMY AND ASTROPHYSICS↗

Reflectance Modeling

The main thrust of the work conducted is in the modeling of spectral responses of discontinuous conifer forest canopies, considered as collections of conical forms that are illuminated at an angle and cast shadows on contracting background. The three first-year objectives are: (1) modification of the variance-driven model to incorporate overlapping of shadows and crowns; (2) collection and analysis of photographic data to calibrate tree-spacing functions; and (3) collection of field data to determine height, spacing, and shape of conifers in open and closed stands in the Goosenest (California) test area. Modification of the invertible model uses a linear approximation to correct for the reduction in shadow and canopy area that occurs when shadows fall on crowns or tree crowns intersect. Comparison of the results of the approximation with a series of Monte Carlo simulations shows that the mean area covered by crowns or shadows is approximated quite accurately, but the variance calculated departs significantly from simulated values. The analysis of tree spacing patterns using air photos showed that spacing tended to be somewhat regular at a scale approximately equal to the average distance between trees, but at larger distance scales counts of trees in cells could be approximated acceptably by a Poisson (random) function.

Strahler, A. H.↗

Optimized mirror shape tuning using beam weightings based on distance, angle of incidence, reflectivity, and power

For glancing-incidence optical systems, such as short-wavelength optics used for nano-focusing, incorporating physical factors in the calculations used for shape optimization can improve performance. Wavefront metrology, including the measurement of a mirror's shape or slope, is routinely used as input for mirror figure optimization on mirrors that can be bent, actuated, positioned, or aligned. Modeling shows that when the incident power distribution, distance from focus, angle of incidence, and the spatially varying reflectivity are included in the optimization, higher Strehl ratios can be achieved. Following the works of Maréchal and Mahajan, optimization of the Strehl ratio (for peak intensity with a coherently illuminated system) occurs when the expectation value of the phase error's variance is minimized. We describe an optimization procedure based on regression analysis that incorporates these physical parameters. This approach is suitable for coherently illuminated systems of nearly diffraction-limited quality. Mathematically, this work is an enhancement of the methods commonly applied for ex situ alignment based on uniform weighting of all points on the surface (or a sub-region of the surface). It follows a similar approach to the optimization of apodized and non-uniformly illuminated optical systems. Significantly, it reaches a different conclusion than a more recent approach based on minimization of focal plane ray errors.

47 OTHER INSTRUMENTATION↗