Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical modeling”

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 739 records · Page 41

A semiparametric latent factor model for large scale temporal data with heteroscedasticity

Large scale temporal data have flourished in a vast array of applications, and their sophisticated structures, especially the heteroscedasticity among subjects with inter- and intra-temporal dependence, have fueled a great demand for new statistical models. In this paper, with covariate information, we consider a flexible model for large scale temporal data with subject-specific heteroscedasticity. Formally, the model employs latent semiparametric factors to simultaneously account for the subject-specific heteroscedasticity and the contemporaneous and/or serial correlations. The subject-specific heteroscedasticity is modeled as the product of the unobserved factor process and subject’s covariate effect, which is further characterized via additive models. For estimation, we propose a two-step procedure. First, the latent factor process and nonparametric loading are recovered through projection-based methods, and following, we estimate the regression components by approaches motivated from the generalized least squares. By scrupulously examining the non-asymptotic rates for recovering the factor process and its loading, we show the consistency and efficiency of estimated regression coefficients in the absence of prior knowledge of latent factor process and subject’s covariate effect. Here, the statistical guarantees remain valid even for finite time points that makes our method particularly appealing when the subjects significantly outnumber the observation time points. Using comprehensive simulations, we demonstrate the finite sample performance of our method, which corroborates the theoretical findings. Finally, we apply our method to a data set of air quality and energy consumption collected at 129 monitoring sites in the United States in 2015.

97 MATHEMATICS AND COMPUTING↗

Cleanroom certification model

The Cleanroom software development methodology is designed to take the gamble out of product releases for both suppliers and receivers of the software. The ingredients of this procedure are a life cycle of executable product increments, representative statistical testing, and a standard estimate of the MTTF (Mean Time To Failure) of the product at the time of its release. A statistical approach to software product testing using randomly selected samples of test cases is considered. A statistical model is defined for the certification process which uses the timing data recorded during test. A reasonableness argument for this model is provided that uses previously published data on software product execution. Also included is a derivation of the certification model estimators and a comparison of the proposed least squares technique with the more commonly used maximum likelihood estimators.

Currit, P. A.↗

Stereo vision for planetary rovers - Stochastic modeling to near real-time implementation

JPL has achieved the first autonomous cross-country robotic traverses to use stereo vision, with all computing onboard the vehicle. This paper describes the stereo vision system, including the underlying statistical model and the details of the implementation. It is argued that the overall approach provides a unifying paradigm for practical domain-independent stereo ranging.

Matthies, Larry↗

Virtual Fit Assessment: Validation using Historical Spacesuit Fit Data

Virtual fit tests using 3D body scans have provided a cost-effective means to predictively assess spacesuit fit for the current and future astronaut population. However, fit is a complex issue influenced by physical interferences, subjective preference, and many other factors. Fit can substantially differ between the type of hardware, environmental conditions, and tasks being performed. Namely, fit changes across different contexts, such as 3D printed mockup evaluations, pressurized one-g suited test events, neutral buoyancy and other training events, and flight extravehicular activities (EVA). These challenges make it difficult to validate virtual fit frameworks to physical suit sizing. This study proposes a new method of validation using the wealth of historically archived suit fit data. Physically Assessed Fit (PAF) data was assessed from astronauts and test volunteers who wore the legacy Extravehicular Mobility Unit (EMU). The data was collected during the past decades at NASA and considered to be most reliable and dependable. Each participant was 3D scanned for body shape and spacesuit-critical dimensions were measured per NASA guidelines. The participants’ size preference for hard upper torso (HUT) assembly was retrieved from their pressurized suit fit check records. Fit assessments were also updated after one-g EVA training, neutral buoyancy training, and EVA flights. Virtually Assessed Fit (VAF) was done by overlaying the CAD model of a HUT with the 3D body scans of the PAF participants. Each body scan was iteratively adjusted for the position inside the HUT to minimize the suit-to-body contact interference while satisfying a set of prescribed requirements. Then the residual contact area, depth, and volume were quantified as VAF metrics. This process was repeated for the different HUT sizes, including medium, large, and extra-large. Statistical modeling is currently in progress and will be presented at the conference. A statistical classifier will be developed to predict a fit probability for different HUT sizes as a function of the corresponding person’s VAF metrics. The HUT size with the highest fit probability will constitute the most likely size selection for the person, and the prediction will be compared against the corresponding PAF data. The participant data will be randomly pre-grouped into either a model development or validation subset. While the model development subset will be used to build the probability model, the validation subset will be used to assess the model accuracy. Also, the locations and magnitudes of suit-body-contacts will be estimated from VAF. This information can identify the critical suit geometry and body shape features that influence suit fit. The variations found in PAF size selections by subjects with similar VAF metrics will allow for investigations into subjective preferences, for example, tight versus loose fit. Overall, this study is expected to provide a structured validation of a virtual suit fit framework, which has not been possible in the past. The outcome can also provide useful insights and potential limitations for interpreting virtual fit tests for future spacesuit designs and population accommodation.

Han Kim↗

AutoBayes Program Synthesis System Users Manual

Program synthesis is the systematic, automatic construction of efficient executable code from high-level declarative specifications. AutoBayes is a fully automatic program synthesis system for the statistical data analysis domain; in particular, it solves parameter estimation problems. It has seen many successful applications at NASA and is currently being used, for example, to analyze simulation results for Orion. The input to AutoBayes is a concise description of a data analysis problem composed of a parameterized statistical model and a goal that is a probability term involving parameters and input data. The output is optimized and fully documented C/C++ code computing the values for those parameters that maximize the probability term. AutoBayes can solve many subproblems symbolically rather than having to rely on numeric approximation algorithms, thus yielding effective, efficient, and compact code. Statistical analysis is faster and more reliable, because effort can be focused on model development and validation rather than manual development of solution algorithms and code.

Schumann, Johann↗

MODELING THE ANOMALY OF SURFACE NUMBER DENSITIES OF GALAXIES ON THE GALACTIC EXTINCTION MAP DUE TO THEIR FIR EMISSION CONTAMINATION

The most widely used Galactic extinction map is constructed assuming that the observed far-infrared (FIR) fluxes come entirely from Galactic dust. According to the earlier suggestion by Yahata et al., we consider how FIR emission of galaxies affects the SFD map. We first compute the surface number density of Sloan Digital Sky Survey (SDSS) DR7 galaxies as a function of the r-band extinction, A {sub r,} {sub SFD}. We confirm that the surface densities of those galaxies positively correlate with A {sub r,} {sub SFD} for A {sub r,} {sub SFD} < 0.1, as first discovered by Yahata et al. for SDSS DR4 galaxies. Next we construct an analytical model to compute the surface density of galaxies, taking into account the contamination of their FIR emission. We adopt a log-normal probability distribution for the ratio of 100 μm and r-band luminosities of each galaxy, y ≡ (νL){sub 100} {sub μm}/(νL) {sub r}. Then we search for the mean and rms values of y that fit the observed anomaly, using the analytical model. The required values to reproduce the anomaly are roughly consistent with those measured from the stacking analysis of SDSS galaxies. Due to the limitation of our statistical modeling, we are not yet able to remove the FIR contamination of galaxies from the extinction map. Nevertheless, the agreement with the model prediction suggests that the FIR emission of galaxies is mainly responsible for the observed anomaly. Whereas the corresponding systematic error in the Galactic extinction map is 0.1-1 mmag, it is directly correlated with galaxy clustering and thus needs to be carefully examined in precision cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

Ground Motion and Intraevent Variability from 3D Deterministic Broadband (0-7.5 Hz) Simulations along a Nonplanar Strike-Slip Fault

In the article, we model deterministic broadband (-7.5 Hz) ground motion from an M w 7.1 bilateral strike-slip earthquake scenario with dynamic rupture propagation along a rough-fault topography embedded in a medium including small-scale velocity and density perturbations. Spectral accelerations (SAs) at periods 0.2-3 s and Arias intensity durations show a similar distance decay (at the level of 1-2 interevent standard deviations above the median) when compared to Next Generation Attenuation-West2 (NGA)-West2 ground-motion prediction equations (GMPEs) using a Q(f) power-law exponent of 0.6-0.8 above 1 Hz in models with a minimum V S of 750 m/s. With a trade-off from Q(f), the median ground motion is slightly increased by scattering from statistical models of small-scale heterogeneity with standard deviation (σ) of the perturbations at the lower end of the observed range (5%) but reduced by scattering attenuation at the upper end (10%) when using a realistic 3D background velocity model. The ground-motion variability is strongly affected by the addition of small-scale media heterogeneity, reducing otherwise large values of intraevent standard deviation closer to those of empirical observations. These simulations generally have intraevent standard deviations for SAs lower than the GMPEs for the modeled bandwidth, with an increasing trend with distance (most pronounced in low-to-moderate scattering media) near the level of observations at distances greater than 35 km from the fault. Durations for the models follow the same increasing trend with distance, in which σ~5% produces the best match to GMPE values. We find that a 3D background-velocity model reduces the pulse period into the expected range by breaking up coherent waves from directivity, generating a lognormal distribution of ground-motion residuals. These results indicate that a strongly heterogeneous medium is needed to produce realistic deterministic broadband ground motions. Finally, the addition of a thin surficial layer with low, frequency-independent Q in the model (with a minimum V S of 750 m/s) controls the high-frequency decay in energy, as measured by the parameter κ, that may be necessary to include as simulations continue to extend to higher frequencies.

Geochemistry & Geophysics↗

Seasonal variability of multiple leaf traits captured by leaf spectroscopy at two temperate deciduous forests

Understanding the temporal patterns of leaf traits is critical in determining the seasonality and magnitude of terrestrial carbon, water, and energy fluxes. However, we lack robust and efficient ways to monitor the temporal dynamics of leaf traits. Here we assessed the potential of leaf spectroscopy to predict and monitor leaf traits across their entire life cycle at different forest sites and light environments (sunlit vs. shaded) using a weekly sampled dataset across the entire growing season at two temperate deciduous forests. In addition, the dataset includes field measured leaf-level directional-hemispherical reflectance/transmittance together with seven important leaf traits [total chlorophyll (chlorophyll a and b), carotenoids, mass-based nitrogen concentration (N mass ), mass-based carbon concentration (C mass ), and leaf mass per area (LMA)]. All leaf traits varied significantly throughout the growing season, and displayed trait-specific temporal patterns. We used a Partial Least Square Regression (PLSR) modeling approach to estimate leaf traits from spectra, and found that PLSR was able to capture the variability across time, sites, and light environments of all leaf traits investigated (R 2 = 0.6–0.8 for temporal variability; R 2 = 0.3–0.7 for cross-site variability; R 2 = 0.4–0.8 for variability from light environments). We also tested alternative field sampling designs and found that for most leaf traits, biweekly leaf sampling throughout the growing season enabled accurate characterization of the seasonal patterns. Compared with the estimation of foliar pigments, the performance of N mass , C mass and LMA PLSR models improved more significantly with sampling frequency. Our results demonstrate that leaf spectra-trait relationships vary with time, and thus tracking the seasonality of leaf traits requires statistical models calibrated with data sampled throughout the growing season. In conclusion, our results have broad implications for future research that use vegetation spectra to infer leaf traits at different growing stages.

54 ENVIRONMENTAL SCIENCES↗

2008 GEM Modeling Challenge: Metrics Study of the Dst Index in Physics-Based Magnetosphere and Ring Current Models and in Statistical and Analytic Specifications

In this paper the metrics-based results of the Dst part of the 2008-2009 GEM Metrics Challenge are reported. The Metrics Challenge asked modelers to submit results for 4 geomagnetic storm events and 5 different types of observations that can be modeled by statistical or climatological or physics-based (e.g. MHD) models of the magnetosphere-ionosphere system. We present the results of over 25 model settings that were run at the Community Coordinated Modeling Center (CCMC) and at the institutions of various modelers for these events. To measure the performance of each of the models against the observations we use comparisons of one-hour averaged model data with the Dst index issued by the World Data Center for Geomagnetism, Kyoto, Japan, and direct comparison of one-minute model data with the one-minute Dst index calculated by the United States Geologic Survey (USGS).

Rastaetter, L.↗

Machine learning models for estimating contamination across different curbside collection strategies

Contaminated recyclables, which are frequently discarded as waste, pose a significant challenge to the implementation of a circular economy. These contaminated recyclables impede the circulation of resources, resulting in higher processing costs at material recovery facilities (MRFs). Over the past few decades, machine learning (ML) models such as linear regression (LR), support vector machine (SVM), and random forest (RF) have evolved to provide new methods for predicting inbound contamination rates in addition to traditional statistical models. In this study, we applied ML models to predict inbound contamination rates using demographic features from 15 counties in the U.S. with different curbside collection strategies. In general, we found that ML models outperformed linear mixed models. Specifically, SVM models had the highest performance (R 2 = 0.75; mean absolute error (MAE) = 0.06), which may be due to their ability to model nonlinear relationships between features and inbound contamination rates. Further, the key predictor was population, with poverty rate being positively correlated and median age negatively correlated with inbound contamination rates. To improve the management of contamination and enhance the implementation of a circular economy, better models are needed to understand and estimate inbound contamination rates as well as identify critical factors in the present and future.

54 ENVIRONMENTAL SCIENCES↗

Structure and Reactivity of Binuclear Cu Active Sites in Cu-CHA Zeolites for Stoichiometric Partial Methane Oxidation to Methanol

Aluminosilicate zeolites exchanged with copper ions facilitate partial methane oxidation (PMO) to methanol in stoichiometric oxidation and reduction cycles, yet the identities of active Cu sites and details of the reaction mechanism remain debated. Here, we use the high-symmetry chabazite (CHA) zeolite framework as a model support to probe the relationship between bulk composition, Cu speciation, and response to various oxidizing and reducing treatments. Density functional theory and first-principles thermodynamics combined with statistical models reveal that Cu speciation and composition depend strongly on Al configuration and external gas conditions. Cu-CHA samples were synthesized to survey broad regions of Si/Al and Cu/Al composition space and framework Al proximity. Characterization by in situ X-ray absorption and UV–visible spectroscopy during exposure to different oxidation conditions reveal that the extent of Cu oxidation is sensitive to activation conditions and thus that both kinetic and thermodynamic factors influence Cu oxidizability in a given material. Similar characterizations during CO reduction reveal that CO titrates Cu 2+ in amounts suggesting the presence of both O- and O 2 -bridged species. In contrast, CH 4 and autoreduction (He) treatments reduce similar but smaller numbers of Cu sites than CO, implicating O 2 -bridged Cu dimers as a potential common intermediate in the former reduction pathways. Furthermore, a systematic increase in methanol yields (per Cu) in stoichiometric PMO cycles increase with the fraction of binuclear O x -bridged Cu sites suggests these species as active sites, as depicted in an updated PMO reaction mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advanced data analysis in inertial confinement fusion and high energy density physics

Bayesian analysis enables flexible and rigorous definition of statistical model assumptions with well-characterized propagation of uncertainties and resulting inferences for single-shot, repeated, or even cross-platform data. This approach has a strong history of application to a variety of problems in physical sciences ranging from inference of particle mass from multi-source high-energy particle data to analysis of black-hole characteristics from gravitational wave observations. The recent adoption of Bayesian statistics for analysis and design of high-energy density physics (HEDP) and inertial confinement fusion (ICF) experiments has provided invaluable gains in expert understanding and experiment performance. In this Review, we discuss the basic theory and practical application of the Bayesian statistics framework. We highlight a variety of studies from the HEDP and ICF literature, demonstrating the power of this technique. Due to the computational complexity of multi-physics models needed to analyze HEDP and ICF experiments, Bayesian inference is often not computationally tractable. Two sections are devoted to a review of statistical approximations, efficient inference algorithms, and data-driven methods, such as deep-learning and dimensionality reduction, which play a significant role in enabling use of the Bayesian framework. We provide additional discussion of various applications of Bayesian and machine learning methods that appear to be sparse in the HEDP and ICF literature constituting possible next steps for the community. We conclude by highlighting community needs, the resolution of which will improve trust in data-driven methods that have proven critical for accelerating the design and discovery cycle in many application areas.

47 OTHER INSTRUMENTATION↗

Reconstruction of implanted marker trajectories from cone-beam CT projection images using interdimensional correlation modeling

Purpose: Cone-beam CT (CBCT) is a widely used imaging modality for image-guided radiotherapy. Most vendors provide CBCT systems that are mounted on a linac gantry. Thus, CBCT can be used to estimate the actual 3-dimensional (3D) position of moving respiratory targets in the thoracic/abdominal region using 2D projection images. The authors have developed a method for estimating the 3D trajectory of respiratory-induced target motion from CBCT projection images using interdimensional correlation modeling. Methods: Because the superior–inferior (SI) motion of a target can be easily analyzed on projection images of a gantry-mounted CBCT system, the authors investigated the interdimensional correlation of the SI motion with left–right and anterior–posterior (AP) movements while the gantry is rotating. A simple linear model and a state-augmented model were implemented and applied to the interdimensional correlation analysis, and their performance was compared. The parameters of the interdimensional correlation models were determined by least-square estimation of the 2D error between the actual and estimated projected target position. The method was validated using 160 3D tumor trajectories from 46 thoracic/abdominal cancer patients obtained during CyberKnife treatment. The authors’ simulations assumed two application scenarios: (1) retrospective estimation for the purpose of moving tumor setup used just after volumetric matching with CBCT; and (2) on-the-fly estimation for the purpose of real-time target position estimation during gating or tracking delivery, either for full-rotation volumetric-modulated arc therapy (VMAT) in 60 s or a stationary six-field intensity-modulated radiation therapy (IMRT) with a beam delivery time of 20 s. Results: For the retrospective CBCT simulations, the mean 3D root-mean-square error (RMSE) for all 4893 trajectory segments was 0.41 mm (simple linear model) and 0.35 mm (state-augmented model). In the on-the-fly simulations, prior projections over more than 60° appear to be necessary for reliable estimations. The mean 3D RMSE during beam delivery after the simple linear model had established with a prior 90° projection data was 0.42 mm for VMAT and 0.45 mm for IMRT. Conclusions: The proposed method does not require any internal/external correlation or statistical modeling to estimate the target trajectory and can be used for both retrospective image-guided radiotherapy with CBCT projection images and real-time target position monitoring for respiratory gating or tracking.

60 APPLIED LIFE SCIENCES↗

HYBRD (High Resolution HYBrid Regional Downscaling) Model: Input data and Code

The HYBRD (HYBrid Regional Downscaling) model is a high-resolution urban land downscaling model that can be used to downscale intermediate urban land use and land cover (LULC) products into a high-resolution (30-meters). HYBRD uses a sequential hybrid process, combining statistical models with cellular-automata-based spatial algorithms. This repository contains all the necessary model code and inputs needed to successfully run HYBRD for Los Angeles, California. The repo also contains example outputs of each model step, except the final simulated raster outputs. Examples of simulated raster outputs for multiple scenarios for Los Angeles are available at DOI: 10.57931/2575233. Please refer to Related Works below.

Land↗

System and Software Reliability (C103)

Within the last decade better reliability models (hardware. software, system) than those currently used have been theorized and developed but not implemented in practice. Previous research on software reliability has shown that while some existing software reliability models are practical, they are no accurate enough. New paradigms of development (e.g. OO) have appeared and associated reliability models have been proposed posed but not investigated. Hardware models have been extensively investigated but not integrated into a system framework. System reliability modeling is the weakest of the three. NASA engineers need better methods and tools to demonstrate that the products meet NASA requirements for reliability measurement. For the new models for the software component of the last decade, there is a great need to bring them into a form that they can be used on software intensive systems. The Statistical Modeling and Estimation of Reliability Functions for Systems (SMERFS'3) tool is an existing vehicle that may be used to incorporate these new modeling advances. Adapting some existing software reliability modeling changes to accommodate major changes in software development technology may also show substantial improvement in prediction accuracy. With some additional research, the next step is to identify and investigate system reliability. System reliability models could then be incorporated in a tool such as SMERFS'3. This tool with better models would greatly add value in assess in GSFC projects.

Wallace, Dolores↗

Infusing Statistical Thinking into the NASA Quesst Community Test Campaign

Statistical thinking permeates many important decisions as NASA plans its Quesst mission, which will culminate in a series of community overflights using the X-59 aircraft to demonstrate low-noise supersonic flight. Month-long longitudinal surveys will be deployed to assess human perception and annoyance to this new acoustic phenomenon. NASA works with a large contractor team to develop systems and methodologies to estimate noise doses, to test and field socio-acoustic surveys, and to study the relationship between the two quantities, dose and response, through appropriate choices of statistical models. This latter dose-response relationship will serve as an important tool as national and international noise regulators debate whether overland supersonic flights could be permitted once again within permissible noise limits. In this presentation we highlight several areas where statistical thinking has come into play, including issues of sampling, classification and data fusion, and analysis of longitudinal survey data that are subject to rare events and the consequences of measurement error. We note several operational constraints that shape the appeal or feasibility of some decisions on statistical approaches, and we identify several important remaining questions to be addressed.

Bayesian model↗