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,081 records · Page 60

An overview of UK C130 observations from the FIRE marine stratocumulus IFO

An overview of the characteristics of the ten flights made by the C130 during the Intensive Field Observations (IFO) are presented. An indication of data quality is given. Comparisons of the carbon dioxide radiation thermometer and the Rosemount suggest there may be potentially serious wetting problems with the latter, especially close to cloud base. The main features of the cloud layers sampled by the C130 are summarized, including microphysical parameters. Many of the results shown are derived from a first order analysis of turbulence quantities. The gross vertical variation of heat, vapor, and liquid water fluxes from case to case, together with the observed vertical velocity variance, a good indicator of convective activity, are discussed. Preliminary results suggest that the different cases can be usefully categorized using these data. The relationship between these and the measured boundary conditions is mentioned. Similarities and differences between the various cases are highlighted, and interpreted in terms of the researchers' current understanding of stratocumulus. Attention is given to particular features of the various flights which may deserve more intensive study.

Nicholls, S.↗

A geostatistical analysis of small-scale spatial variability in bacterial abundance and community structure in salt marsh creek bank sediments

Small-scale variations in bacterial abundance and community structure were examined in salt marsh sediments from Virginia's eastern shore. Samples were collected at 5 cm intervals (horizontally) along a 50 cm elevation gradient, over a 215 cm horizontal transect. For each sample, bacterial abundance was determined using acridine orange direct counts and community structure was analyzed using randomly amplified polymorphic DNA fingerprinting of whole-community DNA extracts. A geostatistical analysis was used to determine the degree of spatial autocorrelation among the samples, for each variable and each direction (horizontal and vertical). The proportion of variance in bacterial abundance that could be accounted for by the spatial model was quite high (vertical: 60%, horizontal: 73%); significant autocorrelation was found among samples separated by 25 cm in the vertical direction and up to 115 cm horizontally. In contrast, most of the variability in community structure was not accounted for by simply considering the spatial separation of samples (vertical: 11%, horizontal: 22%), and must reflect variability from other parameters (e.g., variation at other spatial scales, experimental error, or environmental heterogeneity). Microbial community patch size based upon overall similarity in community structure varied between 17 cm (vertical) and 35 cm (horizontal). Overall, variability due to horizontal position (distance from the creek bank) was much smaller than that due to vertical position (elevation) for both community properties assayed. This suggests that processes more correlated with elevation (e.g., drainage and redox potential) vary at a smaller scale (therefore producing smaller patch sizes) than processes controlled by distance from the creek bank. c2002 Federation of European Microbiological Societies. Published by Elsevier Science B.V. All rights reserved.

Non-NASA Center↗

Modified GMDH-NN algorithm and its application for global sensitivity analysis

Global sensitivity analysis (GSA) is a very useful tool to evaluate the influence of input variables in the whole distribution range. Sobol' method is the most commonly used among variance-based methods, which are efficient and popular GSA techniques. High dimensional model representation (HDMR) is a popular way to compute Sobol' indices, however, its drawbacks cannot be ignored. We show that modified GMDH-NN algorithm can calculate coefficients of metamodel efficiently, so this paper aims at combining it with HDMR and proposes GMDH-HDMR method. The new method shows higher precision and faster convergent rate. Several numerical and engineering examples are used to confirm its advantages. - Highlights: • The GMDH-NN is improved to construct the explicit polynomial model of optimal complexity by self-organization. • The paper aims at combining improved GMDH-NN with HDMR expansions and using it to compute Sobol' indices directly. • The method can be applied in uniform, normal and exponential distribution by using suitable orthogonal polynomials. • Engineering examples, e.g., electronic circuit models can be solved by the presented method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Performance Analysis of the Expected Track Length Estimator

This work investigates whether performance of the MCNP code is improved by the implementation of the expected track length estimator. The expected track length estimator is known to decrease variance when compared to the track length estimator in many classes of Monte Carlo particle transport simulations. However, the computational cost of computing an exponent for every particle collision within a tallied cell has previously been considered prohibitively expensive. Therefore, the goal of this work is to determine whether the variance improvements seen by the expected track length estimator in most cases makes up for the increase in transport time when implemented in the MCNP code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

MJO in Different Orbital Regimes: Role of the Mean State in the MJO’s Amplitude during Boreal Winter

Abstract The Madden–Julian oscillation (MJO) exhibits pronounced seasonality, with one of the key unanswered questions being the following: what controls the maximum in MJO precipitation variance in the Southern Hemisphere during boreal winter? In this study, we examine a set of global climate model simulations in which the eccentricity and precession of Earth’s orbit are altered to change the boreal winter mean state in an attempt to reveal the processes that are responsible for the MJO’s amplitude in the boreal winter. In response to the forced insolation changes, the north–south asymmetry in sea surface temperature is amplified in boreal fall, which intensifies the Hadley circulation in boreal winter. The stronger Hadley circulation yields higher mean precipitation and stronger mean lower-tropospheric westerlies in the southern part of the Indo-Pacific warm pool. The MJO precipitation variability increases significantly where the mean precipitation and lower-tropospheric westerlies strengthen. In the column-integrated moisture budget of the simulated MJO, only surface latent heat flux feedback shows a trend that is consistent with the MJO’s amplitude, suggesting an important role for the surface latent heat flux feedback in the MJO’s amplitude during the boreal winter. An analysis of the moisture–precipitation relationship in the simulations shows that the increase in the mean precipitation lowers the convective moisture adjustment time scale, leading to the increase in precipitation variance. Our results suggest that the mean-state precipitation plays a critical role in the maintenance mechanism of the MJO.

Meteorology & Atmospheric Sciences↗

Uncertainty quantification and sensitivity analysis applications to fuel performance modeling

Best-estimate fuel performance codes such as BISON currently under development at the Idaho National Laboratory, utilize empirical and mechanistic lower-length-scale informed correlations to predict fuel behavior under normal operating and accident reactor conditions. Traditionally, best-estimate results are presented using the correlations with no quantification of the uncertainty in the output metrics of interest. However, there are associated uncertainties in the input parameters and correlations used to determine the behavior of the fuel and cladding under irradiation. Therefore, it is important to perform uncertainty quantification and include confidence bounds on the output metrics that take into account the uncertainties in the inputs. In addition, sensitivity analyses can be performed to determine which input parameters have the greatest influence on the outputs. In this paper we couple the BISON fuel performance code to the DAKOTA uncertainty analysis software to analyze a representative fuel performance problem. The case studied in this paper is based upon rod 1 from the IFA-432 integral experiment performed at the Halden Reactor in Norway. The rodlet is representative of a BWR fuel rod. The input parameters uncertainties are broken into three separate categories including boundary condition uncertainties (e.g., power, coolant flow rate), manufacturing uncertainties (e.g., pellet diameter, cladding thickness), and model uncertainties (e.g., fuel thermal conductivity, fuel swelling). Utilizing DAKOTA, a variety of statistical analysis techniques are applied to quantify the uncertainty and sensitivity of the output metrics of interest. Specifically, we demonstrate the use of sampling methods, polynomial chaos expansions, surrogate models, and variance-based decomposition. The output metrics investigated in this study are the fuel centerline temperature, cladding surface temperature, fission gas released, and fuel rod diameter. The results highlighted the importance of the thermal conductivity of the fuel, the thermal conductivity of the cladding, the gap thickness, the fuel thermal expansion, and the system power factor as influential parameters.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

On-orbit frequency stability analysis of the GPS NAVSTAR-1 quartz clock and the NAVSTARs-6 and -8 rubidium clocks

An on-orbit frequency stability performance analysis of the GPS NAVSTAR-1 quartz clock and the NAVSTARs-6 and -8 rubidium clocks is presented. The clock offsets were obtained from measurements taken at the GPS monitor stations which use high performance cesium standards as a reference. Clock performance is characterized through the use of the Allan variance, which is evaluated for sample times of 15 minutes to two hours, and from one day to 10 days. The quartz and rubidium clocks' offsets were corrected for aging rate before computing the frequency stability. The effect of small errors in aging rate is presented for the NAVSTAR-8 rubidium clock's stability analysis. The analysis includes presentation of time and frequency residuals with respect to linear and quadratic models, which aid in obtaining aging rate values and identifying systematic and random effects. The frequency stability values were further processed with a time domain noise process analysis, which is used to classify random noise process and modulation type.

Mccaskill, T. B.↗

DETECTION OF PLASMA FLUCTUATIONS IN WHITE-LIGHT IMAGES OF THE OUTER SOLAR CORONA: INVESTIGATION OF THE SPATIAL AND TEMPORAL EVOLUTION

This work focuses on the first results from the identification and characterization of periodic plasma density fluctuations in the outer corona, observed in STEREO-A COR1 white-light image time series. A two-dimensional reconstruction of the spatial distribution and temporal evolution of the coronal fluctuation power has been performed over the whole plane of the sky, from 1.4 to 4.0 R{sub Sun }. The adopted diagnostic tool is based on wavelet transforms. This technique, with respect to the standard Fourier analysis, has the advantage of localizing non-persistent fluctuating features and exploring variations of the relating wavelet power in both space and time. The map of the variance of the coronal brightness clearly outlines intermittent spatially coherent fluctuating features, localized along, or adjacent to, the strongest magnetic field lines. In most cases, they do not correspond to the visible coronal structures in the brightness maps. The results obtained provide a scenario in which the solar corona shows quasi-periodic, non-stationary density variations characterized by a wide range of temporal and spatial scales and strongly confined by the magnetic field topology. In addition, structures fluctuating with larger power are larger in size and evolve more slowly. The characteristic periodicities of the fluctuations are comparable to their lifetimes. This suggests that plasma fluctuations lasting only one or two wave periods and initially characterized by a single dominant periodicity either rapidly decay into a turbulent mixed flow via nonlinear interactions with other plasma modes, or they are damped by thermal conduction. The periodic non-stationary coronal fluctuations outlined by the closed field lines at low and mid latitudes might be associated with the existence of slow standing magneto-acoustic waves excited by the convective supergranular motion. The fluctuating ray-like structures observed along open field lines appear to be linked either to the intermittent nature of the processes underlying the generation of magnetic reconnection in the polar regions or to the oscillatory transverse displacements of the coronal ray itself.

79 ASTRONOMY AND ASTROPHYSICS↗

Prostate-Specific Antigen at 4 to 5 Years After Low-Dose-Rate Prostate Brachytherapy Is a Strong Predictor of Disease-Free Survival

Purpose: To determine (1) the prognostic utility of prostate-specific antigen (PSA) concentration at 45 to 60 months (48mPSA) after low-dose-rate prostate brachytherapy (LDR-PB); (2) the predictors of 48mPSA; and (3) the prognostic utility of directional trends between PSA levels at 24, 36, and 48 months after LDR-PB. Methods and Materials: Between 1998 and 2008, 2223 patients with low- and intermediate-risk prostate cancer received LDR-PB monotherapy. A cohort of 1434 of these patients was identified with a documented 48mPSA and no evidence of disease relapse prior to the 48mPSA. In addition, a subset of this cohort (n=585) was identified with ≥72 months of follow-up and documented PSA values at both 24 and 36 months after implantation. Results: Median follow-up time was 76 months. Eight-year Kaplan-Meier disease-free survival (DFS) rates were 100% vs 73.4% for patients with 48mPSA ≤0.2 vs those with >0.2 ng/mL; 99.1% versus 53.8% for a 48mPSA threshold of ≤0.4 versus >0.4 ng/mL, respectively; and 97.3% versus 0% for a threshold of ≤1.0 versus >1.0 ng/mL, respectively. On multivariate analysis, the only factor predictive of DFS was 48mPSA (P<.0001). On subset analysis (n=585), 29 patients had a PSA rise (defined as >0.2 ng/mL) between 24 and 36 months, 24 patients had a rise between 36 and 48 months, and 11 patients had rises over both intervals. Failure rates in these patients were 52%, 79%, and 100%, respectively. On multivariate analysis, initial PSA, androgen deprivation therapy, and dose to 90% of the prostate significantly correlated with 48mPSA but together accounted for only ∼5% of its total variance. Conclusions: The 48mPSA after LDR-PB is highly predictive of long-term DFS. Patients with 48mPSA ≤0.4 ng/mL had a <1% risk of disease relapse at 8 years, whereas all patients with 48mPSA >1.0 ng/mL relapsed. Consecutive PSA rises of >0.2 ng/mL from 24 to 36 months and from 36 to 48 months were also highly predictive of subsequent failure.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Impact of Atmospheric Rivers on Surface Hydrological Processes in Western U.S. Watersheds

Atmospheric rivers (ARs) can significantly modulate surface hydrological processes through the extreme precipitation they produce. Previously, studies have investigated the impact of ARs on individual hydrological processes, but a comprehensive evaluation of ARs’ impact on surface hydrology is lacking. This study used a high-resolution regional climate simulation to quantify the impact of ARs on surface hydrological processes across the western U.S. watersheds. The model performance is evaluated through extensive comparison against observations. Our analysis indicates that ARs produce heavy precipitation but suppress evapotranspiration. Snowpack ablates more during ARs, which is mainly driven by warmer air temperature and increased downward longwave radiation. At the 0oC to 10oC temperature range, ARs increase the probability of snow ablation from 0.33 to 0.57. Soil moisture responds similarly to AR and non-AR precipitation events. The runoff-to-precipitation ratio almost doubles due to modulation by heavy precipitation in the southwestern watersheds or intensification of snow ablation in the northwestern watersheds during AR events compared to non-AR events. From analysis of the relationship between the hydrological responses and different meteorological factors, precipitation, temperature and radiation were identified as the key drivers that distinguish the hydrologic responses between AR and non-AR events. Lastly, analysis of ARs and total runoff at annual scale and April 1st snowpack and winter precipitation shows that ARs explain 30% to 60% of the interannual variance and sharpen the seasonality of water resources availability in the west coast mountain watersheds.

54 ENVIRONMENTAL SCIENCES↗

Networked Gamma Radiation System for Tactical Deployment

A networked gamma radiation detection system with directional sensitivity and energy spectral data acquisition capability is being developed by the National Security Technologies, LLC, Remote Sensing Laboratory to support the close and intense tactical engagement of law enforcement who carry out counterterrorism missions. In the proposed design, three clusters of 2" × 4" × 16" sodium iodide crystals (4 each) with digiBASE-E (for list mode data collection) would be placed on the passenger side of a minivan. To enhance localization and facilitate rapid identification of isotopes, advanced smart real-time localization and radioisotope identification algorithms like WAVRAD (wavelet-assisted variance reduction for anomaly detection) and NSCRAD (nuisance-rejection spectral comparison ratio anomaly detection) will be incorporated. We will test a collection of algorithms and analysis that centers on the problem of radiation detection with a distributed sensor network. We will study the basic characteristics of a radiation sensor network and focus on the trade-offs between false positive alarm rates, true positive alarm rates, and time to detect multiple radiation sources in a large area. Empirical and simulation analyses of critical system parameters, such as number of sensors, sensor placement, and sensor response functions, will be examined. This networked system will provide an integrated radiation detection architecture and framework with (i) a large nationally recognized search database equivalent that would help generate a common operational picture in a major radiological crisis; (ii) a robust reach back connectivity for search data to be evaluated by home teams; and, finally, (iii) a possibility of integrating search data from multi-agency responders.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Accelerated Sensitivity Analysis in High-Dimensional Stochastic Reaction Networks

Existing sensitivity analysis approaches are not able to handle efficiently stochastic reaction networks with a large number of parameters and species, which are typical in the modeling and simulation of complex biochemical phenomena. In this paper, a two-step strategy for parametric sensitivity analysis for such systems is proposed, exploiting advantages and synergies between two recently proposed sensitivity analysis methodologies for stochastic dynamics. The first method performs sensitivity analysis of the stochastic dynamics by means of the Fisher Information Matrix on the underlying distribution of the trajectories; the second method is a reduced-variance, finite-difference, gradient-type sensitivity approach relying on stochastic coupling techniques for variance reduction. Here in this paper we demonstrate that these two methods can be combined and deployed together by means of a new sensitivity bound which incorporates the variance of the quantity of interest as well as the Fisher Information Matrix estimated from the first method. The first step of the proposed strategy labels sensitivities using the bound and screens out the insensitive parameters in a controlled manner. In the second step of the proposed strategy, a finite-difference method is applied only for the sensitivity estimation of the (potentially) sensitive parameters that have not been screened out in the first step. Results on an epidermal growth factor network with fifty parameters and on a protein homeostasis with eighty parameters demonstrate that the proposed strategy is able to quickly discover and discard the insensitive parameters and in the remaining potentially sensitive parameters it accurately estimates the sensitivities. The new sensitivity strategy can be several times faster than current state-of-the-art approaches that test all parameters, especially in “sloppy” systems. In particular, the computational acceleration is quantified by the ratio between the total number of parameters over the number of the sensitive parameters.

97 MATHEMATICS AND COMPUTING↗

SCALE 6.2.4 Validation: Radiation Shielding

For safe and reliable use of computer codes by the community, accuracy must be clearly evaluated. In particular, the nuclear reactor engineering and licensing field needs accurate tools for radiation shielding modeling. Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC) is one such tool, with built-in variance reduction methods distributed within the SCALE code, and its validity is demonstrated in this report for the released version 6.2.4. Representative benchmarks corresponding to shielding analysis are selected for the validation study. Typical experimental results analyzed from those benchmarks include neutron fluxes, detector count rates, detector energy response functions, neutron and gamma doses, foil neutron activation rates and activities, neutron leakage fluxes, and skyshine dose rates. Thousands of points of comparison between experiment and calculation are presented in this work. Other than rare outliers typically explained by either a lack of information or large uncertainties in the experiment conditions, material, or dimensions, MAVRIC agrees well with the experiment results. MAVRIC is also compared to Monte Carlo N-Particle (MCNP) calculations when available, and both codes generally produce good agreements within estimated uncertainties. The selected benchmarks are obtained from reliable sources such as the International Criticality Safety Benchmark Evaluation Project Handbook (ICSBEP Handbook), the Shielding Integral Benchmark Archive & Database (SINBAD), and other shielding validation work found in the literature. Additional datapoints and benchmarks will be added to future versions of this report to incrementally expand the shielding validation suite incrementally.

61 RADIATION PROTECTION AND DOSIMETRY↗

On the performance of digital phase locked loops in the threshold region

Extended Kalman filter algorithms are used to obtain a digital phase lock loop structure for demodulation of angle modulated signals. It is shown that the error variance equations obtained directly from this structure enable one to predict threshold if one retains higher frequency terms. This is in sharp contrast to the similar analysis of the analog phase lock loop, where the higher frequency terms are filtered out because of the low pass filter in the loop. Results are compared to actual simulation results and threshold region results obtained previously.

Hurst, G. T.↗

Trophic classification of Tennessee Valley area reservoirs derived from LANDSAT multispectral scanner data

LANDSAT MSS data from four different dates were extracted from computer tapes using a semiautomated digital data handling and analysis system. Reservoirs were extracted from the surrounding land matrix by using a Band 7 density level slice of 3; and descriptive statistics to include mean, variance, and ratio between bands for each of the four bands were calculated. Significant correlations ( 0.80) were identified between the MSS statistics and many trophic indicators from ground truth water quality data collected at 35 reservoirs in the greater Tennessee Valley region. Regression models were developed which gave significant estimates of each reservoir's trophic state as defined by its trophic state index and explained in all four LANDSAT frames at least 85 percent of the variability in the data. To illustrate the spatial variations within reservoirs as well as the relative variations between reservoirs, a table look up elliptical classification was used in conjunction with each reservoir's trophic state index to classify each reservoir on a pixel by pixel basis and produce color coded thematic representations.

Meinert, D. L.↗

Multitaper Spectral Analysis and Wavelet Denoising Applied to Helioseismic Data

Estimates of solar normal mode frequencies from helioseismic observations can be improved by using Multitaper Spectral Analysis (MTSA) to estimate spectra from the time series, then using wavelet denoising of the log spectra. MTSA leads to a power spectrum estimate with reduced variance and better leakage properties than the conventional periodogram. Under the assumption of stationarity and mild regularity conditions, the log multitaper spectrum has a statistical distribution that is approximately Gaussian, so wavelet denoising is asymptotically an optimal method to reduce the noise in the estimated spectra. We find that a single m-upsilon spectrum benefits greatly from MTSA followed by wavelet denoising, and that wavelet denoising by itself can be used to improve m-averaged spectra. We compare estimates using two different 5-taper estimates (Stepian and sine tapers) and the periodogram estimate, for GONG time series at selected angular degrees l. We compare those three spectra with and without wavelet-denoising, both visually, and in terms of the mode parameters estimated from the pre-processed spectra using the GONG peak-fitting algorithm. The two multitaper estimates give equivalent results. The number of modes fitted well by the GONG algorithm is 20% to 60% larger (depending on l and the temporal frequency) when applied to the multitaper estimates than when applied to the periodogram. The estimated mode parameters (frequency, amplitude and width) are comparable for the three power spectrum estimates, except for modes with very small mode widths (a few frequency bins), where the multitaper spectra broadened the modest compared with the periodogram. We tested the influence of the number of tapers used and found that narrow modes at low n values are broadened to the extent that they can no longer be fit if the number of tapers is too large. For helioseismic time series of this length and temporal resolution, the optimal number of tapers is less than 10.

Komm, R. W.↗

Atmospheric Boundary Layer Wind Data During the Period January 1, 1998 Through January 31, 1999 at the Dallas-Fort Worth Airport: Quality Assessment - Volume 1

The quality of the Aircraft Vortex Spacing System (AVOSS) is critically dependent on representative wind profiles in the atmospheric boundary layer. These winds observed from a number of sensor systems around the Dallas-Fort Worth airport were combined into single vertical wind profiles by an algorithm developed and implemented by MIT Lincoln Laboratory. This process, called the AVOSS Winds Analysis System (AWAS), is used by AVOSS for wake corridor predictions. During times when AWAS solutions were available, the quality of the resultant wind profiles and variance was judged from a series of plots combining all sensor observations and AWAS profiles during the period 1200 to 0400 UTC daily. First, input data was evaluated for continuity and consistency from criteria established. Next, the degree of agreement among all wind sensor systems was noted and cases of disagreement identified. Finally, the resultant AWAS solution was compared to the quality-assessed input data. When profiles differed by a specified amount from valid sensor consensus winds, times and altitudes were flagged. Volume one documents the process and quality of input sensor data. Volume two documents the data processing/sorting process and provides the resultant flagged files.

Zak, J. Allen↗

Networked Gamma Radiation Detection System for Rapid Tactical Response

A networked gamma radiation detection system with directional sensitivity and energy spectral data acquisition capability is being developed by the National Security Technologies, LLC, to support the close and intense tactical engagement of law enforcement who carry out counterterrorism missions. In the proposed design, three clusters of 2' × 4' × 16' sodium iodide crystals (4 each) with digiBASE-E (for list mode data collection) would be placed on the passenger side of a minivan. To enhance localization and facilitate rapid identification of isotopes, advanced smart realtime localization and radioisotope identification algorithms like WAVRAD (wavelet-assisted variance reduction for anomaly detection) and NSCRAD (nuisance-rejection spectral comparison ratio anomaly detection) will be incorporated. We are testing a collection of algorithms and analysis on the problem of radiation detection with a distributed sensor network, examining such a network’s associated characteristics, and also the trade-offs between false positive alarm rates, true positive alarm rates, and time to detect multiple radiation sources in a large area. Empirical and simulation analyses of critical system parameters, such as number of sensors, sensor placement, and sensor response functions, are being examined. This networked system will provide an integrated radiation detection architecture and framework with (i) a large nationally recognized search database equivalent that would help generate a common operational picture in a major radiological crisis; (ii) a robust reach back connectivity for search data to be evaluated by home teams; and, finally, (iii) a possibility of integrating search data from multi-agency responders.

false alarm rate↗