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 451 records · Page 25

A Scale-Invariant Model of Statistical Mechanics and Modified Forms of the First and the Second Laws of Thermodynamics

A scale-invariant statistical theory of fields is presented that leads to invariant definition of density, velocity, temperature, and pressure, The definition of Boltzmann constant is introduced as k(sub k) = m(sub k)v(sub k)c = 1.381 x 10(exp -23) J x K(exp -1), suggesting that the Kelvin absolute temperature scale is equivalent to a length scale. Two new state variables called the reversible heat Q(sub rev) = TS and the reversible work W(sub rev) = PV are introduced. The modified forms of the first and second law of thermodynamics are presented. The microscopic definition of heat (work) is presented as the kinetic energy due to the random (peculiar) translational, rotational, and pulsational motions. The Gibbs free energy of an element at scale Beta is identified as the total system energy at scale (Beta-1), thus leading to an invariant form of the first law of thermodynamics U(sub Beta) = Q(sub Beta) - W(sub Beta) +N(e3)U(sub Beta-1).

Sohrab, Siavash H.↗

Influence of the Convection Electric Field Models on Predicted Plasmapause Positions During Magnetic Storms

In the present work, we determine how three well documented models of the magnetospheric electric field, and two different mechanisms proposed for the formation of the plasmapause influence the radial distance, the shape and the evolution of the plasmapause during the geomagnetic storms of 28 October 2001 and of 17 April 2002. The convection electric field models considered are: Mcllwain's E51) electric field model, Volland-Stern's model and Weimer's statistical model compiled from low-Earth orbit satellite data. The mechanisms for the formation of the plasmapause to be tested are: (i) the MHD theory where the plasmapause should correspond to the last-closed- equipotential (LCE) or last-closed-streamline (LCS), if the E-field distribution is stationary or time-dependent respectively; (ii) the interchange mechanism where the plasmapause corresponds to streamlines tangent to a Zero-Parallel-Force surface where the field-aligned plasma distribution becomes convectively unstable during enhancements of the E-field intensity in the nightside local time sector. The results of the different time dependent simulations are compared with concomitant EUV observations when available. The plasmatails or plumes observed after both selected geomagnetic storms are predicted in all simulations and for all E-field models. However, their shapes are quite different depending on the E-field models and the mechanisms that are used. Despite the partial success of the simulations to reproduce plumes during magnetic storms and substorms, there remains a long way to go before the detailed structures observed in the EUV observations during periods of geomagnetic activity can be accounted for very precisely by the existing E-field models. Furthermore, it cannot be excluded that the mechanisms currently identified to explain the formation of "Carpenter's knee" during substorm events, will', have to be revised or complemented in the cases of geomagnetic storms.

Pierrard, V.↗

Assessing the Transferability of Statistical Predictive Models for Leaf Area Index Between Two Airborne Discrete Return LiDAR Sensor Designs Within Multiple Intensely Managed Loblolly Pine Forest Locations in the South-Eastern USA

Leaf area is an important forest structural variable which serves as the primary means of mass and energy exchange within vegetated ecosystems. The objective of the current study was to determine if leaf area index (LAI) could be estimated accurately and consistently in five intensively managed pine plantation forests using two multiple-return airborne LiDAR datasets. Field measurements of LAI were made using the LiCOR LAI2000 and LAI2200 instruments within 116 plots were established of varying size and within a variety of stand conditions (i.e. stand age, nutrient regime and stem density) in North Carolina and Virginia in 2008 and 2013. A number of common LiDAR return height and intensity distribution metrics were calculated (e.g. average return height), in addition to ten indices, with two additional variants, utilized in the surrounding literature which have been used to estimate LAI and fractional cover, were calculated from return heights and intensity, for each plot extent. Each of the indices was assessed for correlation with each other, and was used as independent variables in linear regression analysis with field LAI as the dependent variable. All LiDAR derived metrics were also entered into a forward stepwise linear regression. The results from each of the indices varied from an R2 of 0.33 (S.E. 0.87) to 0.89 (S.E. 0.36). Those indices calculated using ratios of all returns produced the strongest correlations, such as the Above and Below Ratio Index (ABRI) and Laser Penetration Index 1 (LPI1). The regression model produced from a combination of three metrics did not improve correlations greatly (R2 0.90; S.E. 0.35). The results indicate that LAI can be predicted over a range of intensively managed pine plantation forest environments accurately when using different LiDAR sensor designs. Those indices which incorporated counts of specific return numbers (e.g. first returns) or return intensity correlated poorly with field measurements. There were disparities between the number of different types of returns and intensity values when comparing the results from two LiDAR sensors, indicating that predictive models developed using such metrics are not transferable between datasets with different acquisition parameters. Each of the indices were significantly correlated with one another, with one exception (LAI proxy), in particular those indices calculated from all returns, which indicates similarities in information content for those indices. It can then be argued that LiDAR indices have reached a similar stage in development to those calculated from optical-spectral sensors, but which offer a number of advantages, such as the reduction or removal of saturation issues in areas of high biomass.

Sumnall, Matthew↗

LAVA Contribution to HiLiftPW-3

Launch Ascent and Vehicle Aerodynamics (LAVA) (Computational Fluid Dynamics (CFD) solver) contribution to HiLiftPW-3 (3rd AIAA CFD High Lift Prediction Workshop). Including HL-CRM (NASA High Lift Common Research Model) and JSM (Joint Statistical Modeling) case studies.

Contribution↗

Ellicott City Disasters III: Building a Real-Time Statistical Flood Model for Improving Early Warning Systems in Ellicott City, Maryland

As flood events in the United States grow in frequency and intensity, the uses of applied remote sensing analyses are increasingly necessary for effective flood monitoring and warning systems. The NASA DEVELOP Ellicott City Disasters III project investigated the use of machine learning for applications in flood risk detection to support the improvement of early warning systems in Ellicott City, Maryland. To strengthen the efforts of the Howard County Office of Emergency Management (OEM) in building a more robust flood monitoring system, this term built on the predictive capability of the long-short term memory (LSTM) model created by the second term of this DEVELOP project to create a Sequentially Trained Real-time EstimAted Model (STREAM). Enhancements to the model included the integration of both real-time and predicted weather products from the National Weather Service to increase predictive capacity. These weather products were supplemented by stream gauge data from the OEM as well as real-time radar products. The resultant flood risk model was trained to evaluate input variables and predict stage height in Ellicott City in real time. The model, upgraded to predict stage height up to 8 hours in advance, was incorporated into an online dashboard in a user-friendly interface. The project demonstrated the potential for integration of open data and NASA Earth observations into a flood risk forecasting tool capable of informing real-time decision-making.

Erika Munshi↗

Bhutan Agriculture: Developing a Crop Mask for Rice and Creating a Data Collection Protocol Utilizing Remotely Sensed Data in Bhutan

Rice cultivation in Bhutan has been increasingly threatened by deteriorating soil health and outbreaks of diseases and pests associated with the global change in climate patterns. Field surveys, which the national government of Bhutan has relied on to monitor remote agricultural lands, are becoming increasingly overwhelmed by growing threats to agricultural health. To address these concerns, NASA DEVELOP partnered with the Department of Agriculture of Bhutan, the Bhutan Foundation, and the Ugyen Wangchuck Institute of Conservation and Environmental Research (UWICER) and worked to increase the government of Bhutan’s agricultural monitoring capacity. Utilizing Earth observations including Landsat 8 Operational Land Imager (OLI), Sentinel-1 C-band Synthetic Aperture Radar (C-SAR), Shuttle Radar Topography Mission (SRTM), and Planet imagery, the DEVELOP team worked with NASA SERVIR and created a sampling protocol to identify rice plantations and supplement field surveys for more efficient agriculture monitoring. The analysis focused on districts Paro, Punakha, Samtse, Sarpang, Trongsa, Zhemgang, Wangdue Phodrang, and Samdrup Jongkhar in the year 2020 during the period of transplantation (June) to harvesting of rice (November). The team provided the partners with a sampling protocol for integrating NASA Earth observations into their crop monitoring methods, as well as a crop mask for rice identification and to aid crop management. The crop mask for rice was developed using the Random Forest (RF) classifier for the eight districts of Bhutan. Visually, the random forest model has proved to be more accurate and precise than the classification and Regression Tree model. Statistically, the Random Forest model was 91.8% accurate in identifying rice in Bhutan.

Yeshey Seldon↗

Efficient First-Order Algorithms for Large-Scale, Non-Smooth Maximum Entropy Models with Application to Wildfire Science

Maximum entropy (MaxEnt) models are a class of statistical models that use the maximum entropy principle to estimate probability distributions from data. Due to the size of modern data sets, MaxEnt models need efficient optimization algorithms to scale well for big data applications. State-of-the-art algorithms for MaxEnt models, however, were not originally designed to handle big data sets; these algorithms either rely on technical devices that may yield unreliable numerical results, scale poorly, or require smoothness assumptions that many practical MaxEnt models lack. In this paper, we present novel optimization algorithms that overcome the shortcomings of state-of-the-art algorithms for training large-scale, non-smooth MaxEnt models. Our proposed first-order algorithms leverage the Kullback–Leibler divergence to train large-scale and non-smooth MaxEnt models efficiently. For MaxEnt models with discrete probability distribution of n elements built from samples, each containing m features, the stepsize parameter estimation and iterations in our algorithms scale on the order of O(mn) operations and can be trivially parallelized. Moreover, the strong ℓ1 convexity of the Kullback–Leibler divergence allows for larger stepsize parameters, thereby speeding up the convergence rate of our algorithms. To illustrate the efficiency of our novel algorithms, we consider the problem of estimating probabilities of fire occurrences as a function of ecological features in the Western US MTBS-Interagency wildfire data set. Our numerical results show that our algorithms outperform the state of the art by one order of magnitude and yield results that agree with physical models of wildfire occurrence and previous statistical analyses of wildfire drivers.

Physics↗

Atmospheric turbulence simulation techniques with application to flight analysis

Statistical modeling of atmospheric turbulence is discussed. The statistical properties of atmospheric turbulence, in particular the probability distribution, the spectra, and the coherence are reviewed. Different atmospheric turbulence simulation models are investigated, and appropriate statistical analyses are carried out to verify their validity. The models for simulation are incorporated into a computer model of aircraft flight dynamics. Statistical results of computer simulated landings for an aircraft having characteristics of a DC-8 are reported for the different turbulence simulation techniques. The significance of various degrees of sophistication in the turbulence simulation techniques on the landing performance of the aircraft is discussed.

Wang, S. T.↗

Mott-Anderson freeze-out and the strange matter “horn”

We discuss the √s-dependence of the K {sup +}/π{sup +} ratio in heavy-ion collisions (the “horn” effect) within a Mott-Anderson localization model for chemical freeze-out. The different response of pion and kaon radii to the hot and dense hadronic medium results in different freeze-out conditions. We demonstrate within a simple model that this circumstance enhances the “horn” effect relative to statistical models with universal chemical freeze-out.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

State densities of heavy nuclei in the static-path plus random-phase approximation

We report that nuclear state densities are important inputs to statistical models of compound-nucleus reactions. State densities are often calculated with self-consistent mean-field approximations that do not include important correlations and must be augmented with empirical collective enhancement factors. Here, we benchmark the static-path plus random-phase approximation (SPA + RPA) to the state density in a chain of samarium isotopes 148–155 Sm against exact results (up to statistical errors) obtained with the shell-model Monte Carlo (SMMC) method. The SPA + RPA method incorporates all static fluctuations beyond the mean field together with small-amplitude quantal fluctuations around each static fluctuation. Using a pairing plus quadrupole interaction, we show that the SPA + RPA state densities agree well with the exact SMMC densities for both the even- and odd-mass isotopes. For the even-mass isotopes, we also compare our results with mean-field state densities calculated with the finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation. We find that the SPA + RPA repairs the deficiencies of the mean-field approximation associated with broken rotational symmetry in deformed nuclei and with the violation of particle-number conservation in the pairing condensate. In particular, in deformed nuclei the SPA + RPA reproduces the rotational enhancement of the state density relative to the mean-field state density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Inverse problems-based maximum likelihood estimation of ground reflectivity for selected regions of interest from stripmap SAR data [Regularized maximum likelihood estimation of ground reflectivity from stripmap SAR data]

In this study, we derive a comprehensive forward model for the data collected by stripmap synthetic aperture radar (SAR) that is linear in the ground reflectivity parameters. It is also shown that if the noise model is additive, then the forward model fits into the linear statistical model framework, and the ground reflectivity parameters can be estimated by statistical methods. We derive the maximum likelihood (ML) estimates for the ground reflectivity parameters in the case of additive white Gaussian noise. Furthermore, we show that obtaining the ML estimates of the ground reflectivity requires two steps. The first step amounts to a cross-correlation of the data with a model of the data acquisition parameters, and it is shown that this step has essentially the same processing as the so-called convolution back-projection algorithm. The second step is a complete system inversion that is capable of mitigating the sidelobes of the spatially variant impulse responses remaining after the correlation processing. We also state the Cramer-Rao lower bound (CRLB) for the ML ground reflectivity estimates.We show that the CRLB is linked to the SAR system parameters, the flight path of the SAR sensor, and the image reconstruction grid.We demonstrate the ML image formation and the CRLB bound for synthetically generated data.

47 OTHER INSTRUMENTATION↗

Regional land use studies

Remote sensing technology and data from instrumented satellites and high altitude aircraft are proposed for mapping land use on a current national basis, for monitoring changes and trends, and for creating statistical models which can be manipulated to demonstrate the probable effects of proposed land use and of environmental changes over large areas. Both Apollo spacecraft and aircraft photography were used; the spacecraft pictures delineated the cropland and urban boundaries more clearly. A computer model is also proposed for statistical analysis and for printing out updated maps automatically; this model will include a data bank which can be updated rapidly with changes detected by the computer.

Place, J. L.↗

Strong one-neutron emission from two-neutron unbound states in β decays of the r -process nuclei Ga 86 , 87

β-delayed one-neutron and two-neutron branching ratios (P$_{1n}$ and P$_{2n}$) have been measured in the decay of A = 84 to 87 Ga isotopes at the Radioactive-Isotope Beam Factory (RIBF) at the RIKEN Nishina Center using a high-efficiency array of 3 He neutron counters (BRIKEN). Two-neutron emission was observed in the decay of 84, 85, 87 Ga for the first time and the branching ratios were measured to be P$_{2n}$ = 1.6 (2)%, 1.3 (2)%, and 10.2 ( 28 ) stat ( 5 ) sys %, respectively. One-neutron branching ratio of 87 Ga (P$_{1n}$= 81 (9) stat (8) sys %) and half-life of 29(4) ms were measured for the first time. The branching ratios of 86 Ga were also measured to be P$_{1n}$ = 74 (2) stat (8) sys % and 16.2 (9) stat ( 6 ) sys % with better precision than a previous study. The observation that P$_{1n}$ > P$_{2n}$ for both 86,87 Ga was unexpected and is interpreted as a signature of dominating one-neutron emission from the two-neutron unbound excited states in 86, 87 Ge. In order to interpret the experimental results, shell-model and Hauser-Feshbach statistical model calculations of delayed particle and γ-ray emission probabilities were performed. This model framework reproduces the experimental results. The shell model alone predicts P$_{2n}$ significantly larger than P$_{1n}$ for the 87 Ga decay, and it is necessary to invoke a statistical description to successfully explain the observation that P$_{1n}$ > P$_{2n}$. Our new results demonstrate the relevance and importance of a statistical description of neutron emission for the prediction of the decay properties of multineutron emitters and that it must be included in the r -process modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SIMPLIFIED PREDICTIVE MODELS FOR CO₂ SEQUESTRATION PERFORMANCE ASSESSMENT RESEARCH TOPICAL REPORT ON TASK #3 STATISTICAL LEARNING BASED MODELS

We compare two approaches for building a statistical proxy model (metamodel) for CO₂ geologic sequestration from the results of full-physics compositional simulations. The first approach involves a classical Box-Behnken or Augmented Pairs experimental design with a quadratic polynomial response surface. The second approach used a space-filling maxmin Latin Hypercube sampling or maximum entropy design with the choice of five different meta-modeling techniques: quadratic polynomial, kriging with constant and quadratic trend terms, multivariate adaptive regression spline (MARS) and additivity and variance stabilization (AVAS). Simulations results for CO₂ injection into a reservoir-caprock system with 9 design variables (and 97 samples) were used to generate the data for developing the proxy models. The fitted models were validated with using an independent data set and a cross-validation approach for three different performance metrics: total storage efficiency, CO₂ plume radius and average reservoir pressure. The Box-Behnken–quadratic polynomial metamodel performed the best, followed closely by the maximin LHS–kriging metamodel.

Mishra, Srikanta↗