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,369 records · Page 76

Polar motion from laser range measurements of GEOS-3

Using two-day arcs of GEOS-3 laser data, simultaneous solutions for pole position components, x sub p and y sub p, and orbit elements have been obtained for the period spanning 3 February to 6 March 1976 using three NASA Goddard Space Flight Center laser stations located near Washington, D.C. (STALAS) and on the islands of Bermuda and Grand Turk. The results are in general agreement with the BIH results. However, because of the locations of the laser sites, the x sub p solution is weaker than the y sub p solution. The x sub p and y sub p estimates were smoothed with a straight line by weighted least squares using the variance associated with the pole estimates as weights in order to reflect the effect of widely different data distributions. The smoothed y sub p differs by one meter with respect to the BIH smoothed values and the smoothed x sub p differs by about two meters. Spectral analysis of the results has identified frequencies associated with the orbital motion indicating the need for further improvements in the model of the physical system.

Schutz, B. E.↗

GRAVSAT error analysis by means of global spectral analyses of the marine geoid from SEASAT altimeter data

The geographical variability of short wavelength geoid power spectra (geoid roughness), was mapped for the world's oceans between latitudes 72 deg N and 72 deg S. A spectral analysis of SEASAT altimeter data, reduced to sea surface heights, was performed at 2 minute intervals for 15 consecutive days of the 3-day repeat orbit. The geoid roughness represented by these spectra for wavelengths shorter than about 220 km is separated from the total sea height variance and is displayed in the form of a global contour map. The global average geoid roughness is 32 cm RMS, varying from a high in excess of 2 m RMS near deep ocean trenches to a low of 2 cm RMS in the southeast Pacific near the East Pacific Rise. This average value agrees well with previous estimates based on gravimetry and GEOS-3 altimetry. In general, the smoothest areas in the marine geoid overlie relatively young seafloor adjacent mid-ocean spreading centers, where even short-wavelength topographic variations tend to be isostatically compensated.

Source record↗

Mechanisms for SAR imaging of ocean surface phenomena: Theory and experiment

Understanding the SAR response to surface wave is a central issue in the analysis of SAR ocean images. The imaging mechanism for gravity waves and the practical question of just which characteristics of the ocean wave field can be measured remotely using SAR were examined. Assessments of wave imaging theory are based primarily on comparisons of the directional wave height variance spectrum psi (K) measured by in situ buoys with estimates from SAR images. Other criteria are also recommended, e.g., the effects of focus adjustments. It is assumed that fluctuations in SAR image intensity are proportional to fluctuations in ocean surface height. If this were true, the Fourier power spectrum of a SAR image and corresponding surface measurements of psi would coincide. Differences between SAR estimates based on this hypothesis and buoy measurements of psi are then used to begin the assessment of rival wave imaging theories.

Vesecky, J. F.↗

An automatic frequency control loop using overlapping DFTs (Discrete Fourier Transforms)

An automatic frequency control (AFC) loop is introduced and analyzed in detail. The new scheme is a generalization of the well known Cross Product AFC loop that uses running overlapping discrete Fourier transforms (DFTs) to create a discriminator curve. Linear analysis is included and supported with computer simulations. The algorithm is tested in a low carrier to noise ratio (CNR) dynamic environment, and the probability of loss of lock is estimated via computer simulations. The algorithm discussed is a suboptimum tracking scheme with a larger frequency error variance compared to an optimum strategy, but offers simplicity of implementation and a very low operating threshold CNR. This technique can be applied during the carrier acquisition and re-acquisition process in the Advanced Receiver.

Aguirre, S.↗

Numerical Investigation and Optimization of a Flushwall Injector for Scramjet Applications at Hypervelocity Flow Conditions

An investigation utilizing Reynolds-averaged simulations (RAS) was performed in order to find optimal designs for an interdigitated flushwall injector suitable for scramjet applications at hypervelocity conditions. The flight Mach number, duct height, spanwise width, and injection angle were the design variables selected to maximize two objective functions: the thrust potential and combustion efficiency. A Latin hypercube sampling design-of-experiments method was used to select design points for RAS. A methodology was developed that automated building geometries and generating grids for each design. The ensuing RAS analysis generated the performance database from which the two objective functions of interest were computed using a one-dimensional performance utility. The data were fitted using four surrogate models: an artificial neural network (ANN) model, a cubic polynomial, a quadratic polynomial, and a Kriging model. Variance-based decomposition showed that both objective functions were primarily driven by changes in the duct height. Multiobjective design optimization was performed for all four surrogate models via a genetic algorithm method. Optimal solutions were obtained at the upper and lower bounds of the flight Mach number range. The Kriging model obtained an optimal solution set that predicted high values for both objective functions. Additionally, three challenge points were selected to assess the designs on the Pareto fronts. Further sampling among the designs of the Pareto fronts are required in order to lower the errors and perform more accurate surrogate-based optimization. sed optimization.

Shenoy, Rajiv R.↗

Data-Informed Evaluation Framework for Integrated Energy Systems: Insights from Power, Process Heat, and Hydrogen Production Applications

The multi-criteria decision analysis (MCDA) framework provides a systematic evaluation of the diverse preferences and performance metrics associated with alternative solutions. This approach is advantageous over a single-criterion methodology, which are only valid under conditions that assume ceteris paribus or an "apples-to-apples" comparison. However, selecting suitable technologies for integrated energy systems (IES) can be likened to an "apples-to-oranges" comparison, given the heterogeneous factors at stake. These factors include economics and performance parameters, geological compatibility, and environmental impacts. Consequently, past research has often employed a mixture of qualitative and quantitative criteria tailored to the specific interests of each study. While the method proves effective in handling the intricate interplay of criteria, the resulting rankings and scores can vary from study to study. This inconsistency is introduced from the use of subjectively defined thresholds and weights. As a result, decision-makers frequently find it challenging to establish clear connections between specific criteria and the resulting scores, as the transformation of criteria into ordinal scores results in a substantial loss of information. To address this challenge, we introduce a data-informed IES evaluation framework that offers comprehensive, interpretable, and traceable evaluations backed by quantifiable rationale. First, we identified key IES evaluation criteria from a decade of literature, focusing on relevant IES applications in power, process heat, and hydrogen production. We leveraged state-of-the-art cost estimates from the Idaho National Laboratory (INL) and technical data from 78 reactor designs from the International Atomic Energy Agency (IAEA) and the Organization for Economic Co-operation and Development - Nuclear Energy Agency (OECD-NEA). Lastly, we established thresholds by analyzing the mean, variance, root mean square, and slope of values across alternatives, categorizing the preferences of decision-makers into distinct utility functions, such as linear, saturating, exponential, and stepwise. Our approach yielded two main outcomes: (1) it provided consistent assessments across different stakeholder groups and (2) it visualized uncertainties in the decision-making context via comprehensive sensitivity analysis. To demonstrate the impact of our framework, we conducted case studies on 6 reactor designs (AP1000, NuScale, BWRX-300, Xe-100, eVinci, iMSR) for the three applications. Our data-driven framework proved to be highly effective in addressing heterogenous uncertainties faced by varied decision-makers? preferences and IES applications, as well as cost and technical estimates of advanced reactors.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Data-Informed Evaluation Framework for Integrated Energy Systems: Insights from Power, Process Heat, and Hydrogen Production Applications

The multi-criteria decision analysis (MCDA) framework provides a systematic evaluation of the diverse preferences and performance metrics associated with alternative solutions. This approach is advantageous over a single-criterion methodology, which are only valid under conditions that assume ceteris paribus or an "apples-to-apples" comparison. However, selecting suitable technologies for integrated energy systems (IES) can be likened to an "apples-to-oranges" comparison, given the heterogeneous factors at stake. These factors include economics and performance parameters, geological compatibility, and environmental impacts. Consequently, past research has often employed a mixture of qualitative and quantitative criteria tailored to the specific interests of each study. While the method proves effective in handling the intricate interplay of criteria, the resulting rankings and scores can vary from study to study. This inconsistency is introduced from the use of subjectively defined thresholds and weights. As a result, decision-makers frequently find it challenging to establish clear connections between specific criteria and the resulting scores, as the transformation of criteria into ordinal scores results in a substantial loss of information. To address this challenge, we introduce a data-informed IES evaluation framework that offers comprehensive, interpretable, and traceable evaluations backed by quantifiable rationale. First, we identified key IES evaluation criteria from a decade of literature, focusing on relevant IES applications in power, process heat, and hydrogen production. We leveraged state-of-the-art cost estimates from the Idaho National Laboratory (INL) and technical data from 78 reactor designs from the International Atomic Energy Agency (IAEA) and the Organization for Economic Co-operation and Development - Nuclear Energy Agency (OECD-NEA). Lastly, we established thresholds by analyzing the mean, variance, root mean square, and slope of values across alternatives, categorizing the preferences of decision-makers into distinct utility functions, such as linear, saturating, exponential, and stepwise. Our approach yielded two main outcomes: (1) it provided consistent assessments across different stakeholder groups and (2) it visualized uncertainties in the decision-making context via comprehensive sensitivity analysis. To demonstrate the impact of our framework, we conducted case studies on 6 reactor designs (AP1000, NuScale, BWRX-300, Xe-100, eVinci, iMSR) for the three applications. Our data-driven framework proved to be highly effective in addressing heterogenous uncertainties faced by varied decision-makers? preferences and IES applications, as well as cost and technical estimates of advanced reactors.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Evaluation of Weakly Informed Priors for FLEX Data

On behalf of the U. S. Nuclear Regulatory Commission (NRC), Idaho National Laboratory (INL) has reviewed the development and use of “weakly informed prior” (WIP) probability distributions to obtain industry-wide failure probability and rate distributions for portable “FLEX” equipment, as documented in the Pressurized Water Reactor Owner’s Group (PWROG) report, FLEX Equipment Data Collection and Analysis (PWROG 18043-P) by MM Degonish. The review includes a description of the methods used by the PWROG analysts. It also includes an implementation of those methods for comparison purposes and an implementation of two constrained-noninformative (CN) distributions for further WIP comparisons. The report documents several issues that the review identified in the development of the WIP distributions. They relate to the quality and representativeness of the current portable “FLEX” component data, the choice of informed distributions for related installed components, and the lack of verification for the factors chosen to weaken the information from the installed components. Two errors occurred in the computations: the posterior distribution of the WIP was used to modify the WIP itself, resulting in a lack of independence between the WIP and the data; and an error was made in the conversion of WIP mean and variance values to beta distributions. The INL reviewers think the WIP distributions presented in the report are untenable for use in risk assessments at the present time because of these issues.

99 GENERAL AND MISCELLANEOUS↗

Unbiased Ratio Estimator for k-Eigenvalue in Monte Carlo Criticality Calculations

In criticality calculations, the effective multiplication factor (k{sub eff}) is an important parameter for reactor analysis. It can be interpreted as the ratio of the neutron production rate to the loss rate in a reactor core. In Monte Carlo (MC) criticality calculations, k{sub eff} is obtained for each active cycle i as: k-i = P{sub i}/L{sub i} for I = 1 to I where P{sub i} and L{sub i} are the production rate and the loss rate at cycle i, respectively, and they are normalized by the nominal source size for cycle. Since all neutron particles in a cycle will be absorbed or leaks out at the end, the loss rate, L{sub i}, can be set as the source rate, which is the unity. Thus, in MCNP5, k{sub i} is simply calculated with P{sub i} by using three conventional estimators; collision, absorption, and track-length estimators. In domain decomposed MC criticality calculations based on the fission and surface source (FSS) iteration method, k{sub i} is calculated by a ratio of P{sub i} and L{sub i}, where both P{sub i} and L{sub i} are random variables, sampled by using the same particle histories. When the numerator and the denominator are correlated, the variance of the ratio can be smaller than that of the numerator. Then, the final estimate of k{sub eff} is provided as the sample mean of k{sub i}, which will be referred as 'simple' ratio estimator in this paper. In Ref., it is shown that this ratio estimator converges faster and fluctuates less than the conventional estimator. A similar ratio estimator is also used in Serpent 2, as an implicit estimate of k{sub eff}. However, this type of 'simple' ratio estimator is biased, which means that the expectation of the ratio is different from the ratio of the expectation of the numerator and that of the denominator. In this paper, an unbiased ratio estimator which was originally derived in the cytometry analysis is introduced to estimate k{sub eff} and is tested on a continuous-energy UO{sub 2} pin-cell problem. In this test problem, it is shown that the variance of the unbiased ratio estimator is smaller than those of the conventional estimators provided by MCNP5, while the simple ratio estimator can be significantly biased when the number of histories per cycle is too small. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Deep Analysis Net with Causal Embedding for Coal-fired Power Plant Fault Detection and Diagnosis (DANCE4CFDD)

Fault detection and diagnosis is critical to power plant operation to ensure attaining high reliability while reducing operation cost. As more renewable power is introduced to the power grid, traditional fossil power plants take on the extra burden of excessive load cycling to compensate the generation variability from renewable power. Such load cycling will pose more reliability challenges to power plant operation. There are a number of challenges faced by today’s asset health management system in coal- fired (or gas) power plants: 1) high-dimensional nonlinear interaction among multiple time series measurements; 2) high measurement variance induced by operational conditions/modes; 3) variation among asset types and plant configurations; and 4) a small number of faulty events to learn from. To cope with these challenges, today’s fielded asset health management systems rely heavily on manual efforts from domain experts and hand-crafted features or rules based on domain knowledge. Despite its role in plant reliability, such a practice is costly and hinders its scalability and sustainability, particularly when a plant undergoes modifications. The objective of this project is to develop a novel end-to-end AI learning system that is trainable (i.e., the AI representation of a complex system behavior can be directly learned from properly labeled data) for accurate fault detection and root cause analysis. The ability to create a fault detection model directly from time series could alleviate the efforts associated with today’s asset management solution development. In the course of this project, we have achieved the following: Created an AI model development environment incorporating state-of-the-art neural network architectures for rapid model development and evaluation; Developed novel learning strategies for training of fault detection model; Developed special-purpose neural network architecture embedded with variable association graph aiming for better interpretability; Developed a learning strategy to leverage a small number of faulty events for enhanced fault detection capability; Conducted detailed experimental study based on public benchmark datasets and demonstrated the effectiveness of the proposed solution; and Validated the developed system with data from both a coal-fired plant boiler dynamic simulation model and real-world coal-fired power plant covering multiple asset and fault types. Overall, the project attained a technology readiness level of TRL 5 from TRL 2 at the beginning of the project.

20 FOSSIL-FUELED POWER PLANTS↗

Uncertainty amplification due to density/refractive index gradients in background-oriented schlieren experiments

Here, we theoretically analyze the effect of density/refractive index gradients on the measurement precision of background-oriented schlieren (BOS) experiments by deriving the Cramer–Rao lower bound (CRLB) for the 2D centroid estimation process. A model is derived for the diffraction limited image of a dot viewed through a medium containing density gradients that includes the effect of the experimental parameters such as the magnification and f-number. It is shown using the model that nonlinearities in the density gradient field lead to blurring of the dot image. This blurring amplifies the effect of image noise on the centroid estimation process, leading to an increase in the CRLB and a decrease in the measurement precision. The ratio of position uncertainties of a dot in the reference and gradient images is shown to be a function of the ratio of the dot diameters and dot intensities. We termed this parameter the amplification ratio (A F ), and a methodology for reporting position uncertainties in tracking-based BOS measurements is proposed. The theoretical predictions of the dot position estimation variance from the CRLB are compared to ray tracing simulations, and agreement is obtained. The uncertainty amplification is also demonstrated on experimental BOS images of flow induced by a spark discharge, where it is seen that regions of high amplification ratio correspond to regions of density gradients. This analysis elucidates the dependence of the position uncertainty on density and refractive index gradient-induced distortion parameters, provides a methodology for accounting its effect on uncertainty quantification and provides a framework for optimizing experiment design.

42 ENGINEERING↗

Indra: a public computationally accessible suite of cosmological N -body simulations

ABSTRACT Indra is a suite of large-volume cosmological N-body simulations with the goal of providing excellent statistics of the large-scale features of the distribution of dark matter. Each of the 384 simulations is computed with the same cosmological parameters and different initial phases, with 10243 dark matter particles in a box of length 1 h−1 Gpc, 64 snapshots of particle data and halo catalogues, and 505 time-steps of the Fourier modes of the density field, amounting to almost a petabyte of data. All of the Indra data are immediately available for analysis via the SciServer science platform, which provides interactive and batch computing modes, personal data storage, and other hosted data sets such as the Millennium simulations and many astronomical surveys. We present the Indra simulations, describe the data products and how to access them, and measure ensemble averages, variances, and covariances of the matter power spectrum, the matter correlation function, and the halo mass function to demonstrate the types of computations that Indra enables. We hope that Indra will be both a resource for large-scale structure research and a demonstration of how to make very large data sets public and computationally accessible.

Falck, Bridget↗

Wireless Sensor Networks for Improved Snow Water Equivalent and Runoff Estimates

We leverage the frontiers of the Internet of Things technology in a recently developed end-to-end wireless sensor network (WSN) system that samples, collects, stores, and displays mountain hydrology measurements in near real-time. At the core of the system lies an ultra-low power, radio channel-hoping, and self-organizing mesh that allows for remote autonomous sampling of snow. Such properties, combined with a rugged weather-sealed design of the devices and multi-level data replication, provides reliable real-time data at spatial and temporal scales previously impractical to achieve in mountain environments. The system was deployed at three 1 km 2 sites across the North Fork of the Feather River basin with a cluster of 12 sensor nodes for each location. Measurements show that existing operational autonomous systems are non-representative spatially, with biases that can reach up to 50%. A comparison between a wet and dry year showed that snow depths exhibit strong multi-scale inter-year spatial stationarity with major rank conservation. Temporally dense analysis using elastic net regression shows that dominant features at the sub-km 2 scale are site-dependent and differ from the watershed scale. Newly introduced explanatory variables, based on the nearest neighbor with a Landsat assimilated historical product, consistently explained up to 90% of the variance in the watershed-scale SWE for both years. At two WSN sites, lagged cross-correlation of snowmelt with stream flow measurements showed a significant improvement of up to 100% compared with existing systems, suggesting that WSNs can be instrumental in improving runoff forecasting and water management.

47 OTHER INSTRUMENTATION↗

EV Profile Capture

NextGen Profiles' EV profile capture efforts aimed to explore the variance in performance and evaluate how different operational conditions influence production EV charging behavior. Data were collected at a frequency of 10 Hz from both the EV and EVSE during each charge session. These charge session parameters were then entered into a time-series database for further analysis. The data were gathered under different operational conditions to examine the effects of various factors such as battery state of charge, battery temperature, vehicle condition, smart charge management, and EVSE limitations. The EV profile capture dataset includes extensive high-power charging data from 16 different EVs—comprising light-, medium-, and heavy-duty vehicles—along with EVSE from various suppliers. To protect confidentiality, the EV and EVSE metadata are anonymized, and the publicly released datasets are aggregated to 0.1-Hz frequency.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Roughness of the marine geoid from Seasat altimetry

The geographical variability of short wavelength geoid power spectra (geoid roughness) has been mapped for the world's oceans between latitudes 72 deg N and 72 deg S. A spectral analysis of Seasat altimeter data, reduced to sea surface heights, has been performed at 2-min intervals for 15 consecutive days of the 3-day repeat orbit. The geoid roughness represented by these spectra for wavelengths shorter than about 220 km is separated from the total sea height variance and is displayed in the form of a global contour map. The global average geoid roughness is 32 cm RMS, varying from a high in excess of 2 m RMS near deep ocean trenches to a low of 2 cm RMS in the southeast Pacific near the east Pacific rise. This average value agrees well with previous estimates based on gravimetry and GEOS 3 altimetry. In general, the smoothest areas in the marine geoid overlie relatively young sea floor adjacent mid-ocean spreading centers, where even short wavelength topographic variations tend to be isostatically compensated.

Brown, R. D.↗

Climatology of gravity waves over Poker Flat, Alaska for 1983

An analysis of short-period wind fluctuations over Poker Flat, Alaska, obtained using the Poker Flat mesosphere-stratosphere-troposphere radar is presented. Results are shown for the troposphere and lower stratosphere as well as for the upper mesosphere and lower thermosphere. Contours depict various levels of wind variance (m2s-2). These results pertain only to wind fluctuation periods lying between one and six hours. These particular fluctuations are generally considered to arise primarily from atmospheric gravity waves. Insofar as this is true, the figure thus describes a general climatology of gravity waves at high latitudes.

Balsley, B. B.↗

Analytic Steady-State Accuracy of a Spacecraft Attitude Estimator

This paper extends Farrenkopf's analysis of a single-axis spacecraft attitude estimator using gyro and angle sensor data to include the angle output white noise of a rate-integrating gyro. Analytic expressions are derived for the steady-state pre-update and post-update angle and drift bias variances and for the state update equations. It is shown that only part of the state update resulting from the angle sensor measurement is propagated to future times.

Markley, F. Landis↗

Status and Plans for the WCRP/GEWEX Global Precipitation Climatology Project (GPCP)

Status and plans for GPCP are presented along with scientific findings from the current data set. Global and large regional rainfall variations and possible long-term changes are examined using the 26-year (1979-2004) monthly dataset from the Global Precipitation Climatology Project (GPCP). One emphasis is to discriminate among the variations due to ENSO, volcanic events and possible long-term changes. Although the global change of precipitation in the data set is near zero, the data set does indicate an upward trend (0.13 mm/day/25yr) and a downward trend (-0.06 mm/day/25yr) over tropical oceans and lands (25S-25N), respectively. This corresponds to a 4% increase (ocean) and 2% decrease (land) during this time period. Simple techniques are derived to attempt to eliminate variations due to ENSO and major volcanic eruptions in the Tropics. Using only annual values two "volcano years" are determined by examining ocean-land coupled variations in precipitation related to ENSO and other phenomena. The outlier years coincide with Pinatubo and El Chicon eruptions. The ENSO signal is reduced by deriving mean ocean and land values for El Nino, La Nina and neutral conditions based on Nino 3.4 SST and normalizing the annual ocean and land precipitation to the neutral set of cases. The impact of the two major volcanic eruptions over the past 25 years is estimated to be about a 5% reduction in tropical rainfall. The modified data set (with ENSO and volcano effect at least partially removed) retains the same approximate linear change slopes over the data set period, but with reduced variance leading to significance tests with results in the 90-95% range. Intercomparisons between the GPCP, SSM/I (1988-2004), and TRMM (1998-2004) satellite rainfall products and alternate gauge analyses over land are made to attempt to increase or decrease confidence in the changes seen in the GPCP analysis.

Adkerm Robert F.↗