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 721 records · Page 40

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↗

Space Propulsion Design and Analysis

This software provides an improved methodology for predicting launcher base pressure and heat loads for RSRM (Reusable Solid Rocket Motor) launchers by accounting for complex anisotropic stress/strains and variable turbulent Prandtl and Schmidt numbers. A "building block" approach to turbulence model development, and validation has been applied for improved missile/launcher base region analysis. Modifications to existing kappa - epsilon turbulence models and application of scalar variance models are incorporated into a RANS-based method for aeropropulsive flow modeling, directly related to base flow methodology. (RANS stands for Reynolds-averaged Navier-Stokes.) The models are applied in a RANS solver framework and can improve analysis of other complex flow fields. The enhanced models provide a more accurate predictive capability for improving the design and analysis of RSRM launcher configuration. The kappa - epsilon model enhancements have been shown to improve the capability for predicting turbulence effects in base blow environments. The scalar variance models have been assessed over a wide range of flow configurations to improve prediction of turbulent scalar mixing.

Sinha, Neeraj↗

Machine learning and process-based modeling of spatiotemporal changes in active layer thickness across Alaska

Permafrost degradation poses a growing threat to infrastructure stability and ecosystem resilience in the rapidly warming Arctic. We investigated the spatiotemporal dynamics of active layer thickness (ALT) across Alaska by integrating field observations, environmental datasets, a physically based Stefan model, and machine learning (ML) techniques. Using weather projections from the Coupled Model Intercomparison Project Phase 6 under two Shared Socioeconomic Pathways (SSP 2-4.5 and SSP 5-8.5), we assessed ALT sensitivity to projected future weather conditions. The random forest (RF) model outperformed the Stefan approach in predicting ALT on the training dataset (R² = 0.84 vs. 0.53) but demonstrated lower generalizability on the test dataset (R² = 0.24 vs. 0.54). The root mean square error (RMSE) for the RF model for training and testing ranged from 14 to 22 cm, compared to 17 and 18 cm for the Stefan model. Variable importance analysis revealed that mean annual temperature and slope angle were the strongest predictors of ALT, accounting for 19% and 18% of the variance, respectively, followed by sediment transport index (14%) and stream power index (11%). Comparative analysis of baseline ALT predictions showed the Stefan model tended to project a thicker active layer (mean ± SD: 65 ± 16 cm), compared to the RF model (mean ± SD: 59 ± 8.8) cm). Both models indicated a latitudinal gradient in ALT, with shallower depths at higher latitudes. Projected ALT increases by 2100 were estimated at 3.3 ± 2.2 cm under SSP 2-4.5 and 5.9 ± 4.0 cm under SSP 5-8.5 for the ML model, whereas the Stefan model projected substantially larger increases of 13 ± 2.6 cm (SSP 2-4.5) and 28 ± 4.4 cm (SSP5-8.5). Spatial analysis showed the greatest ALT increases in northern Alaska, with relatively smaller changes in southern regions. These findings highlight the complex, multifactorial nature of ALT dynamics and the value of hybrid modeling approaches. As rising temperatures accelerate permafrost thaw, changes in ALT can disrupt ecosystems, damage infrastructures, and enhance the release of stored soil carbon, highlighting the urgent need for improved predictive capabilities to inform adaptation strategies in the Arctic.

Climate sciences↗

A statistical note on the redundancy of nine standard baroreflex parameters

An accepted method for measuring the responsiveness of the carotid-cardiac baroreflex to arterial pressure changes is to artificially stimulate the baroreceptors in the neck with a pressurized neck chamber. Nine physiological responses to this type of stimulation are quantified and used as indicators of the baroreflex response function. Thirty male humans between the ages of 27 and 46 underwent the carotid-cardiac baroreflex test. The data for the nine response parameters were analyzed by principle component factor analysis. The results indicated that 92.5 percent of the total variance across all nine parameters could be explained in four dimensions. The first two dimensions reflected location points for R-R interval and carotid distending pressure, respectively. The third factor was composed of measures reflecting the gain (responsiveness) of the reflex. The fourth dimension was the ratio of baseline R-R interval to the maximal R-R interval response during simulated hypertension. The data suggest that the analysis of all nine baroreflex parameters is likely to be redundant and researchers should account for these redundancies either in their analyses or conclusions.

Ludwig, David A.↗

Using landscape typologies to model socioecological systems: Application to agriculture of the United States Gulf Coast

Agricultural enterprises are vulnerable to the effects of climate variability and change. Improved understanding of the determinants of vulnerability and adaptive capacity in agricultural systems is important for projecting and managing future climate risk. At present, three analytical tools dominate methodological approaches to understanding agroecological vulnerability to climate: process-based crop models, empirical crop models, and integrated assessment models. A common weakness of these approaches is their limited treatment of socio-economic conditions and human agency in modeling agroecological processes and outcomes. This study proposes a framework that uses spatial cluster analysis to generate regional socioecological typologies that capture geographic variance in regional agricultural production and enable attribution of that variance to climatic, topographic, edaphic, and socioeconomic components. This framework was applied to historical corn production (1986-2010) in the U.S. Gulf of Mexico region as a testbed. The results demonstrate that regional socioeconomic heterogeneity is an important driving force in human dominated ecosystems, which we hypothesize, is a function of the link between socioeconomic conditions and the adaptive capacity of agricultural systems. Meaningful representation of future agricultural responses to climate variability and change is contingent upon understanding interactions among biophysical conditions, socioeconomic conditions, and human agency their incorporation in predictive models.

58 GEOSCIENCES↗

Hyperspectral imaging for real-time waste materials characterization and recovery using endmember extraction and abundance detection

Hyperspectral imaging, combined with advanced spectral unmixing techniques and artificial intelligence, offers a powerful solution for improving material identification and classification. Here, this study evaluates the effectiveness of the pixel purity index and the sequential maximum angle convex cone algorithms in extracting and validating spectral signatures from pure samples of paper components (cellulose and lignin) and plastic (polypropylene). Principal-component analysis showed that both algorithms captured nearly all relevant variance for the tested materials. Spectral signatures were compared using the spectral angle mapper, revealing high similarity in the short-wave infrared region and greater variability in the visible near-infrared range. The methodology was then applied to a disposable coffee cup to detect and quantify mixed materials, accurately estimating material abundance and object area with less than 1% error. This approach enhances material classification, supporting product verification, quality control, and automated sorting for sustainable waste management and resource recovery.

36 MATERIALS SCIENCE↗

Characterizing warm atmospheric boundary layer over land by combining Raman and Doppler lidar measurements

PBL plays a critical role in the atmosphere by transferring heat, moisture, and momentum. The warm PBL has a distinct diurnal cycle including daytime convective mixing layer (ML) and nighttime residual layer developments. Thus, for PBL characterization and process study, simultaneous determinations of PBL height (PBLH) and ML height (MLH) are necessary. Here, new approaches are developed to provide reliable PBLH and MLH to characterize warm PBL evolution. The approaches use Raman lidar (RL) water vapor mixing ratio (WVMR) and Doppler lidar (DL) vertical velocity measurements at the Southern Great Plains (SGP) atmospheric observatory, which was established by the Atmospheric Radiation Measurement (ARM) user facility. Compared with widely used lidar aerosol measurements for PBLH, WVMR is a better trace for PBL vertical mixing. For PBLH, the approach classifies PBL water vapor structures into a few general patterns, then uses a slope method and dynamic threshold method to determine PBLH. For MLH, wavelet analysis is used to re-construct 2-D variance from DL vertical wind velocity measurements according to the turbulence eddy size to minimize the impacts of gravity wave and eddy size on variance calculations; then, a dynamic threshold method is used to determine MLH. Remotely-sensed PBLHs and MLHs are compared with radiosonde measurements based on the Richardson number method. Good agreements between them confirm that the proposed new algorithms are reliable for PBLH and MLH characterization. The algorithms are applied to warm seasons’ RL and ML measurements at the SGP site for five years to study warm season PBL structure and processes. The weekly composited diurnal evolutions of PBLHs and MLHs in warm climate were provided to illustrate diurnal and seasonal PBL evolutions. This reliable PBLH and MLH dataset will be valuable for PBL process study, model evolution, and PBL parameterization improvement.

54 ENVIRONMENTAL SCIENCES↗

Earth-based observations of Martian albedo changes

The large-scale albedo features of Mars, as observed with Earth-based telescopes, undergo secular changes, cyclic seasonal changes, seemingly random day-to-day changes, and cyclic diurnal changes. The suspension or redistribution of dust may be involved in all four. Some of the recent findings, based on extensive new Earth-based data, are at variance with classical ideas. Current methods of observation and analysis need to be pursued until all Martian seasons have been equally sampled.

Baum, W. A.↗

Analyses of mean and turbulent motion in the tropics with the use of unequally spaced data

Wind velocities from 25 km to 60 km over Ascension Island, Fort Sherman and Kwajalein for the period January 1970 to December 1971 are analyzed in order to achieve a better understanding of the mean flow, the eddy kinetic energy and the Eulerian time spectra of the eddy kinetic energy. Since the data are unequally spaced in time, techniques of one-dimensional covariance theory were utilized and an unequally spaced time series analysis was accomplished. The theoretical equations for two-dimensional analysis or wavenumber frequency analysis of unequally spaced data were developed. Analysis of the turbulent winds and the average seasonal variance and eddy kinetic energy of the turbulent winds indicated that maximum total variance and energy is associated with the east-west velocity component. This is particularly true for long period seasonal waves which dominate the total energy spectrum. Additionally, there is an energy shift for the east-west component into the longer period waves with altitude increasing from 30 km to 50 km.

Kao, S. K.↗

Long-Term Global Morphology of Gravity Wave Activity Using UARS Data

This is the first quarter's report on research to extract global gravity-wave data from satellite data and to model those observations synoptically. Preliminary analysis of global maps of extracted middle atmospheric temperature variance from the CRISTA instrument is presented, which appear to contain gravity-wave information. Corresponding simulations of global gravity-wave and mountain-wave activity during this mission period are described using global ray-tracing and mountain-wave models, and interesting similarities among simulated data and CRISTA data are noted. Climatological simulations of mesospheric gravity-wave activity using the HWM-03 wind-temperature climatology are also reported, for comparison with UARS MLS data. Preparatory work on modeling of gravity wave observations from space-based platforms and subsequent interpretation of the MLS gravity-wave product are also described. Preliminary interpretation and relation to the research objectives are provided, and further action for the next quarter's research is recommended.

Eckermann, Stephen D.↗

Advanced Electrocardiography Can Identify Occult Cardiomyopathy in Doberman Pinschers

Recently, multiple advanced resting electrocardiographic (A-ECG) techniques have improved the diagnostic value of short-duration ECG in detection of dilated cardiomyopathy (DCM) in humans. This study investigated whether 12-lead A-ECG recordings could accurately identify the occult phase of DCM in dogs. Short-duration (3-5 min) high-fidelity 12-lead ECG recordings were obtained from 31 privately-owned, clinically healthy Doberman Pinschers (5.4 +/- 1.7 years, 11/20 males/females). Dogs were divided into 2 groups: 1) 19 healthy dogs with normal echocardiographic M-mode measurements: left ventricular internal diameter in diastole (LVIDd . 47mm) and in systole (LVIDs . 38mm) and normal 24-hour ECG recordings (<50 ventricular premature complexes, VPCs); and 2) 12 dogs with occult DCM: 11/12 dogs had increased M-mode measurements (LVIDd . 49mm and/or LVIDs . 40mm) and 5/11 dogs had also >100 VPCs/24h; 1/12 dogs had only abnormal 24-hour ECG recordings (>100 VPCs/24h). ECG recordings were evaluated via custom software programs to calculate multiple parameters of high-frequency (HF) QRS ECG, heart rate variability, QT variability, waveform complexity and 3-D ECG. Student's t-tests determined 19 ECG parameters that were significantly different (P < 0.05) between groups. Principal component factor analysis identified a 5-factor model with 81.4% explained variance. QRS dipolar and non-dipolar voltages, Cornell voltage criteria and QRS waveform residuum were increased significantly (P < 0.05), whereas mean HF QRS amplitude was decreased significantly (P < 0.05) in dogs with occult DCM. For the 5 selected parameters the prediction of occult DCM was performed using a binary logistic regression model with Chi-square tested significance (P < 0.01). ROC analyses showed that the five selected ECG parameters could identify occult ECG with sensitivity 89% and specificity 83%. Results suggest that 12-lead A-ECG might improve diagnostic value of short-duration ECG in earlier detection of canine DCM as five selected ECG parameters can with reasonable accuracy identify occult DCM in Doberman Pinschers. Future extensive clinical studies need to clarify if 12-lead A-ECG could be useful as an additional screening test for canine DCM.

Spiljak, M.↗

FIGS: spectral fitting constraints on the star formation history of massive galaxies since the cosmic noon

We constrain the stellar population properties of a sample of 52 massive galaxies – with stellar mass log (M(s)/M(⊙)) ≳ 10.5 – over the redshift range 0.5 < z < 2 by use of observer-frame optical and near-infrared slitless spectra from Hubble Space Telescope’s ACS and WFC3 grisms. The deep exposures (∼100 ks) allow us to target individual spectra of massive galaxies to F160W = 22.5 AB. Our spectral fitting approach uses a set of six base models adapted to the redshift and spectral resolution of each observation, and fits the weights of the base models, including potential dust attenuation, via a Markov Chain Monte Carlo method. Our sample comprises a mixed distribution of quiescent (19) and star-forming galaxies (33). We quantify the width of the age distribution (Δt) that is found to dominate the variance of the retrieved parameters according to principal component analysis. The population parameters follow the expected trend towards older ages with increasing mass, and Δt appears to weakly anticorrelate with stellar mass, suggesting a more efficient star formation at the massive end. As expected, the redshift dependence of the relative stellar age (measured in units of the age of the Universe at the source) in the quiescent sample rejects the hypothesis of a single burst (aka monolithic collapse). Radial colour gradients within each galaxy are also explored, finding a wider scatter in the star-forming subsample, but no conclusive trend with respect to the population parameters.

Ignacio Ferreras↗

Objective determination of image end-members in spectral mixture analysis of AVIRIS data

Spectral mixture analysis has been shown to be a powerful, multifaceted tool for analysis of multi- and hyper-spectral data. Applications of AVIRIS data have ranged from mapping soils and bedrock to ecosystem studies. During the first phase of the approach, a set of end-members are selected from an image cube (image end-members) that best account for its spectral variance within a constrained, linear least squares mixing model. These image end-members are usually selected 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. We approach the problem of image end-member determination 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↗

The meandering Gulf Stream as seen by the Geosat altimeter - Surface transport, position, and velocity variance from 73 deg to 46 deg W

Results are presented of an analysis of the surface geostrophic velocity field for the Gulf Stream region for the position, structure, and surface transport of the Gulf Stream for 2.5 yr of the Geosat altimeter Exact Repeat Mission. Synthetic data using a Gaussian velocity profile were generated and fit to the sea surface residual heights to create a synthetic mean sea surface height field and profiles of absolute geostrophic currents. An analysis of the model parameters and the actual geostrophic velocity profiles revealed two different flow regimes for the Gulf Stream connected by a narrow transition region coincident with the New England Seamount Chain. The upstream region was found to exhibit relatively straight Gulf Stream paths, long Eulerian time scales, and eastward propagating meanders. The downstream region had more large meanders, no consistent propagation direction, and shorter Eulerian time scales. A 25-percent reduction in surface transport occurred in the transition region, with a corresponding reduction in current speed and no change in Gulf Stream width.

Kelly, Kathryn A.↗

Southern Idaho Ecological Conservation: Investigating the Impact of Targeted Grazing to Improve Wetland Habitat in the Sterling Wildlife Management Area

Wetland ecosystems are vital for biodiversity conservation and ecosystem services. The Sterling Wildlife Management Area in Bingham County, Idaho, has management concerns about decadent and accumulated vegetation growth encroaching on wetland habitat, which presents challenges for wildlife, decreases biodiversity, and limits public access. Targeted grazing has been proposed as a sustainable alternative to chemical herbicides or burning. Land managers introduced targeted cattle grazing in January 2021 to reduce biomass. NASA DEVELOP partnered with the Idaho Department of Fish and Game to determine the impact of grazing using NASA Earth observations from Landsat 8 Operational Land Imager (OLI) in Google Earth Engine (GEE). Images were processed with TerrSet’s Land Change Modeler and ArcGIS Pro’s Change Detection Wizard to understand land changes following grazing. A Normalized Difference Vegetation Index (NDVI) analysis was performed to assess impacts on vegetation productivity and compare variance in biomass before and after grazing. A Normalized Difference Water Index (NDWI) was used to compare changes in the wetland and its vegetation content to evaluate the suitability of the area for migratory birds post-grazing. Results showed a decrease in the vegetation index and an increase in the water index postgrazing. The DEVELOP team’s analysis suggests that grazing helps break down thick, senesced vegetation and increase soil moisture. Providing a workflow model will aid partners in continuing to monitor this management area and other management areas across the state.

change detection↗

Risk Importance Ranking of Fire Data Parameters to Enhance Fire PRA Model Realism

Fire is historically and analytically a significant contributor to nuclear power plant risk. The level of fire risk and the methods, tools and data for modeling this risk is highly debated by experts. One area of debate is the input data used in fire modeling and how to deal with this data’s high uncertainty. This report outlines initial work performed for determining the key parameters causing this uncertainty and how it propagates into nuclear power plant models. This research paves the way for the development of methods to reduce fire data uncertainty used in modeling. The Nuclear Regulatory Commission has mandated that nuclear power plants perform fire risk modeling. However, there are several issues with the current risk modeling implementation that affect the results. Approved modeling methods can be overly conservative and often do not match plant experience. Also, the data used in the modeling can have high uncertainties and is influenced by expert judgement. To evaluate input data uncertainty, researchers performed an initial review of several fire experiments done at Sandia National Laboratories. Uncertainties for fire data can come from many sources, such as experiment design constraints, environmental conditions, or other plant-specific aspects. There are many different significant and insignificant parameters driving the uncertainty. Additionally, the uncertainty of the different input data used in the fire modeling could have a significant or insignificant effect on the entire plant risk. A four-step methodology was developed to perform Integrated Probabilistic Risk Assessment Importance Ranking. A demonstration case using these steps was set up and three of the four steps were completed in fiscal year (FY) 2019 and the fourth step done FY 2020. These steps are: 1. The qualitative analysis of potential sources was conducted with the following items identified for the demonstration. • Maximum heat release rate • Time to maximum heat release rate • Duration of max heat release rate • Time to decay • Thermal conductivity of concrete • Specific heat of concrete • Density of concrete • Cable jacket thickness 2. A quantitative characterization of dominant sources of uncertainty was performed. A list of distributions and determined values of the dominant sources is shown in Appendix A. 3. A quantitative screening of the potential sources of uncertainty using Morris Elementary Effects Analysis was completed. An experimental model using the physics-based fire modeling tool Fire Dynamics Simulator was developed and coupled with the Risk Analysis Virtual Environment. The Morris analysis identified at least two parameters that can be eliminated as significant contributors (specific heat of concrete and cable jacket thickness). 4. Global importance measure (Global IM) analysis to generate a comprehensive ranking based on their influence on the plant risk. In this research, a moment-independent Global IM is used since it can address (a) uncertainty in the input parameters of the fire model, (b) uncertainty in the risk outputs, and (c) non-linearity and interactions among input parameters in the fire model, more accurately than the correlation-based and variance-based global methods. The observations from the research showed that, depending on the initial and boundary conditions of the fire scenarios, fire-induced damage could have a very small probability and could be dominated by the tail of the uncertainty distribution; hence, the accuracy of the correlation-based and variance-based methods is questionable. The moment-independent Global IM analysis in this research provides a better understanding of how experimental uncertainty data affects industry’s plant models and where improvements in that data will have the largest benefit for improving fire modeling accuracy in causing core damage. Among the five unscreened parameters obtained from the Morris EE analysis, the Global IM analysis results for the case study indicated that max heat release rate and fire location are the most important parameters. The report also outlines benefits of using a unified computational platform that integrates the underlying simulations (e.g., a fire progression model), quantitative screening (using the Morris EE method), and the Global IM analysis. A unified platform can (i) facilitate the ranking of input parameters considering multiple key fire scenarios simultaneously, rather than considering one scenario at a time, (ii) contribute to more explicit and accurate treatment of dependencies at multiple levels of Fire PRA, (iii) facilitate the sampling-based uncertainty quantification for Fire PRA, and (iv) help generating both “industry-wide” and "plant-specific" ranking of uncertainty sources in Fire PRA. Future research should be done to include additional parameters such as detection/suppression or cable fire spread. Adding a PRA software such as SAPHIRE to the RAVEN platform would help with plant model integration and improve treatment of fire-induced dependency. The I-PRA risk importance ranking methodology offered in this report can provide valuable information for efficiently (a) enhancing the realism of Fire PRA for existing plants and (b) supporting the development of Dynamic Fire PRA for advanced reactors and new plants.

97 MATHEMATICS AND COMPUTING↗