Search NASA⌕ Search

SEARCH · Search NASA

Results for “random variables”

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 667 records · Page 37

Automated classification of periodic variable stars detected by the wide-field infrared survey explorer

We describe a methodology to classify periodic variable stars identified using photometric time-series measurements constructed from the Wide-field Infrared Survey Explorer (WISE) full-mission single-exposure Source Databases. This will assist in the future construction of a WISE Variable Source Database that assigns variables to specific science classes as constrained by the WISE observing cadence with statistically meaningful classification probabilities. We have analyzed the WISE light curves of 8273 variable stars identified in previous optical variability surveys (MACHO, GCVS, and ASAS) and show that Fourier decomposition techniques can be extended into the mid-IR to assist with their classification. Combined with other periodic light-curve features, this sample is then used to train a machine-learned classifier based on the random forest (RF) method. Consistent with previous classification studies of variable stars in general, the RF machine-learned classifier is superior to other methods in terms of accuracy, robustness against outliers, and relative immunity to features that carry little or redundant class information. For the three most common classes identified by WISE: Algols, RR Lyrae, and W Ursae Majoris type variables, we obtain classification efficiencies of 80.7%, 82.7%, and 84.5% respectively using cross-validation analyses, with 95% confidence intervals of approximately ±2%. These accuracies are achieved at purity (or reliability) levels of 88.5%, 96.2%, and 87.8% respectively, similar to that achieved in previous automated classification studies of periodic variable stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Forest Potential Productivity Mapping by Linking Remote-Sensing-Derived Metrics to Site Variables

A fine-resolution region-wide map of forest site productivity is an essential need for effective large-scale forestry planning and management. In this study, we incorporated Sentinel-2 satellite data into an increment-based measure of forest productivity (biomass growth index (BGI)) derived from climate, lithology, soils, and topographic metrics to map improved BGI (iBGI) in parts of North American Acadian regions. Initially, several Sentinel-2 variables including nine single spectral bands and 12 spectral vegetation indices (SVIs) were used in combination with forest management variables to predict tree volume/ha and height using Random Forest. The results showed a 10–12 % increase in out of bag (OOB) r2when Sentinel-2 variables were included in the prediction of both volume and height together with BGI. Later, selected Sentinel-2 variables were used for biomass growth prediction in Maine, USA and New Brunswick, Canada using data from 7738 provincial permanent sample plots. The Sentinel-2 red-edge position (S2REP) index was identified as the most important variable over others to have known influence on site productivity. While a slight improvement in the iBGI accuracy occurred compared to the base BGI model (~2%), substantial changes to coefficients of other variables were evident and some site variables became less important when S2REP was included.

site productivity↗

Evaluating the impact of anatomical and physiological variability on human equivalent doses using PBPK models

Abstract Addressing human anatomical and physiological variability is a crucial component of human health risk assessment of chemicals. Experts have recommended probabilistic chemical risk assessment paradigms in which distributional adjustment factors are used to account for various sources of uncertainty and variability, including variability in the pharmacokinetic behavior of a given substance in different humans. In practice, convenient assumptions about the distribution forms of adjustment factors and human equivalent doses (HEDs) are often used. Parameters such as tissue volumes and blood flows are likewise often assumed to be lognormally or normally distributed without evaluating empirical data for consistency with these forms. In this work, we performed dosimetric extrapolations using physiologically based pharmacokinetic (PBPK) models for dichloromethane (DCM) and chloroform that incorporate uncertainty and variability to determine if the HEDs associated with such extrapolations are approximately lognormal and how they depend on the underlying distribution shapes chosen to represent model parameters. We accounted for uncertainty and variability in PBPK model parameters by randomly drawing their values from a variety of distribution types. We then performed reverse dosimetry to calculate HEDs based on animal points of departure for each set of sampled parameters. Corresponding samples of HEDs were tested to determine the impact of input parameter distributions on their central tendencies, extreme percentiles, and degree of conformance to lognormality. This work demonstrates that the measurable attributes of human variability should be considered more carefully and that generalized assumptions about parameter distribution shapes may lead to inaccurate estimates of extreme percentiles of HEDs.

Toxicology↗

Estimation of 3D Woven Design Sensitivities Using a Rapid Multiscale Analysis Technique

Highly-refined finite element models of three-dimension (3D) woven composite systems currently require excessive computational demands that limit their use in sensitivity analysis, uncertainty quantification, and optimization. An alternative analysis methodology was developed using the NASA Multiscale Analysis Tool (NASMAT) where multiscale models of a 3D woven composite (including inter-tow matrix voids and constituent failure) can be completed on a single central processing unit(CPU)on the order of ~30 s. To develop inputs and validation data for the NASMAT model, coupon and acid-digesting testing and x-ray computed tomography were performed. The NASMAT inputs were parameterized using a set of 25 input variables and distributions. These inputs were randomly sampled to generate a total of 100,000 NASMAT analyses that could be used to understand the influence of different material and geometric properties on the warp and weft-direction stiffness and strength. These analyses (including pre/post-processing) were performed in less than eight hours on a 120 CPU cluster. The computational efficiency of the NASMAT model enabled a sensitivity analysis to be performed, and dominant input variables were able to be identified. Key results were consistent with theoretical and experimental observations for the specific 3D woven system studied in this work.

NASMAT↗

Estimation of 3D Woven Design Sensitivities Using a Rapid Multiscale Analysis Technique

Highly-refined finite element models of three-dimension (3D) woven composite systems currently require excessive computational demands that limit their use in sensitivity analysis, uncertainty quantification, and optimization. An alternative analysis methodology was developed using the NASA Multiscale Analysis Tool (NASMAT) where multiscale models of a 3D woven composite (including inter-tow matrix voids and constituent failure) can be completed on a single central processing unit (CPU) on the order of ~30s. To develop inputs and validation data for the NASMAT model, coupon and acid-digesting testing and x-ray computed tomography were performed. The NASMAT inputs were parameterized using a set of 25 input variables and distributions. These inputs were randomly sampled to generate a total of 100,000 NASMAT analyses that could be used to understand the influence of different material and geometric properties on the warp and weft-direction stiffness and strength. These analyses (including pre/post-processing) were performed in less than eight hours on a 120 CPU cluster. The computational efficiency of the NASMAT model enabled a sensitivity analysis to be performed, and dominant input variables were able to be identified. Key results were consistent with theoretical and experimental observations for the specific 3D woven system studied in this work.

NASMAT↗

Wave theory for microwave remote sensing of a half-space random medium with three-dimensional variations

The two-variable expansion technique is used to solve for the mean Green's functions from the Dyson equation under the nonlinear approximation for a half-space random medium with three-dimensional correlation functions. The Bethe-Salpeter equations are then solved under the ladder approximation. The radiative transfer equations, which have been applied extensively in the study of microwave remote sensing problems, are derived under these approximations. The limiting cases of large and small horizontal correlation lengths are discussed. It is found that there is only one propagation constant except for the case of large horizontal correlation lengths, in which there are two propagation constants. It is shown that boundary layer appears in first-order solutions and does not appear in zeroth-order solutions.

Tsang, L.↗

Random permutations fix a worst case for cyclic coordinate descent

Abstract Variants of the coordinate descent approach for minimizing a nonlinear function are distinguished in part by the order in which coordinates are considered for relaxation. Three common orderings are cyclic (CCD), in which we cycle through the components of $x$ in order; randomized (RCD), in which the component to update is selected randomly and independently at each iteration; and random-permutations cyclic (RPCD), which differs from CCD only in that a random permutation is applied to the variables at the start of each cycle. Known convergence guarantees are weaker for CCD and RPCD than for RCD, though in most practical cases, computational performance is similar among all these variants. There is a certain type of quadratic function for which CCD is significantly slower than for RCD; a recent paper by Sun & Ye (2016, Worst-case complexity of cyclic coordinate descent: $O(n^2)$ gap with randomized version. Technical Report. Stanford, CA: Department of Management Science and Engineering, Stanford University. arXiv:1604.07130) has explored the poor behavior of CCD on functions of this type. The RPCD approach performs well on these functions, even better than RCD in a certain regime. This paper explains the good behavior of RPCD with a tight analysis.

Lee, Ching-pei↗

Propagation in a shearing plasma. II - Turbulence and the frequency range of pulsar microstructure

Numerical calculations are presented which explore the propagation of radio waves and pulses through a turbulent shearing plasma. A stochastic model is used to describe the turbulent velocity field, in which random fluctuations are superposed on a mean profile. Both Gaussian wave packets and shot-noise pulses polarized in the plane of shearing acquire quasi-periodic intensity modulations above the rest-frame plasma frequency. The frequency range over which these modulations appear is shown to depend on the spatial correlation length of the turbulent fluctuations and on their standard deviation from the mean velocity. The period of the modulations is variable from one realization of the random process to the next, and the average period is frequency dependent. The results of these calculations lend further support to a model for pulsar microstructure in which periodic micropulses are pure temporal modulations of the emitted radiation due to propagation effects in magnetospheric shearing regions. It is shown that turbulence on a scale of approximately 10 cm can produce these quasi-periodic modulations over a frequency range comparable to that of observed microstructure.

Harding, A. K.↗

A global reference atmospheric model for surface to orbital altitudes

An empirical atmospheric model has been developed which generates values for pressure, density, temperature and winds from surface levels to orbital altitudes. The output parameters consist of components for: (1) latitude, longitude, and altitude dependent monthly means; (2) quasi-biennial oscillations; and (3) random perturbations to partially simulate the variability due to synoptic, diurnal, planetary wave and gravity wave variations. The monthly mean models consist of: NASA's four dimensional worldwide model for height, latitude, and longitude dependent monthly means from the surface to 25 km; and a newly developed latitude-longitude dependent model which is an extension of the Groves latitude dependent model for the region between 25 and 90 km. The Jacchia 1970 model is used above 90 km and is faired with the modified Groves values between 90 and 115 km. Quasi-biennial and random variation perturbations are computed from parameters determined from various empirical studies, and are added to the monthly mean values.

Justus, C. G.↗

Signatures of aerosol-cloud interactions in GiOcean: a coupled global reanalysis with two-moment cloud microphysics

Aerosols influence the Earth's radiative balance through direct interactions with radiation and by affecting cloud properties. Anthropogenic aerosols have led to cooling during the industrial era through aerosol–cloud interactions (ACI), including aerosol effects on cloud microphysical properties and the subsequent adjustments. However, large uncertainties remain in Earth system models (ESMs) regarding the magnitude of this cooling. In part, ESMs substantially disagree on cloud properties, thermodynamics, the hydrological cycle, and general circulation. Reanalysis provides a useful avenue for exploring the impact of ACI on clouds and radiation because its atmosphere is forced to match realistic conditions through the assimilation of observations. Here, we explore the impact of ACI on clouds in the GiOcean reanalysis – the first to incorporate aerosol-cloud interactions. We contrast variables important for ACI between GiOcean and satellite observations and develop 2-dimensional lookup tables of ACI for both using a source-sink budget perspective to attribute the changes in cloud droplet number (Nd) and liquid water path (LWP) to aerosol and meteorology. A compositing analysis using lookup tables shows that GiOcean captures key aspects of aerosol–cloud–precipitation interactions, including (1) activation of aerosol into cloud droplets, (2) effective precipitation scavenging of Nd, (3) suppression of precipitation by high Nd in regions with heavy aerosol emissions. In contrast, satellite observations do not exhibit clear patterns for processes (2) and (3). Random Forest analysis shows that interannual variability in Nd and LWP over the Northern Hemisphere ocean in GiOcean is primarily driven by precipitation, consistent with satellite observations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Random-Cluster Correlation Inequalities for Gibbs Fields

In this note we prove a correlation inequality for local variables of a Gibbs field based on the connectivity in a random cluster representation of the non overlap configuration distribution of two independent copies of the field. As a consequence, we show that absence of a particular type of percolation (of Machta–Newman–Stein blue bonds) implies uniqueness of Gibbs distribution in EA Spin Glasses. In dimension two this could constitute a step towards a proof that the critical temperature is zero.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computational simulation of probabilistic lifetime strength for aerospace materials subjected to high temperature, mechanical fatigue, creep, and thermal fatigue

The results of a fourth year effort of a research program conducted for NASA-LeRC by The University of Texas at San Antonio (UTSA) are presented. The research included on-going development of methodology that provides probabilistic lifetime strength of aerospace materials via computational simulation. A probabilistic material strength degradation model, in the form of a randomized multifactor interaction equation, is postulated for strength degradation of structural components of aerospace propulsion systems subjected to a number of effects or primitive variables. These primitive variables may include high temperature, fatigue, or creep. In most cases, strength is reduced as a result of the action of a variable. This multifactor interaction strength degradation equation was randomized and is included in the computer program, PROMISC. Also included in the research is the development of methodology to calibrate the above-described constitutive equation using actual experimental materials data together with regression analysis of that data, thereby predicting values for the empirical material constants for each effect or primitive variable. This regression methodology is included in the computer program, PROMISC. Actual experimental materials data were obtained from industry and the open literature for materials typically for applications in aerospace propulsion system components. Material data for Inconel 718 was analyzed using the developed methodology.

Boyce, Lola↗

Quantifying the Variability of “Fixed-Width” Buffers on Harvested Lands in Western Oregon and Washington

In contemporary forest management, buffers of unharvested trees are left along streams to protect riparian and aquatic ecosystems. Buffer regulations often focus on specific minimum width requirements, which aid in straightforward regulation and application, but minimum widths also suggest buffered edges are uniform and contain little variability. Conceptual papers suggest that alternative buffer configurations may offer greater flexibility in landscape-level protection, increase forest complexity, and enhance aquatic and riparian biodiversity. However, before considering alternatives to fixed-width buffers, it is necessary to quantify the inherent variability in current buffer practices present on the landscape. In this study, we used aerial imagery to quantify variability of buffer widths in two hundred randomly selected recently harvested units on managed land in Oregon and Washington with both fish-bearing and non-fish-bearing sections of stream. Wider buffers on larger streams had a greater magnitude of variability, but when normalized by stream size, variation was greater in smaller streams, and overall, variation ranged from 25% to 50% of the mean width. Despite local variability, buffer widths rarely fell below 9.14 m (30 ft). The variation quantified here provides an initial measure of variability to inform future management, given emerging interest in variable retention buffers.

58 GEOSCIENCES↗

Polarization Angle Swings in Blazars Detected in the Millimeter-wave with the South Pole Telescope

We present the first systematic search for electric vector position angle (EVPA) swings in the millimeter-wave (mm-wave) emission of blazars, using five years of observations from the South Pole Telescope SPT-3G camera at 95 and 150 GHz, and investigate their connection to gamma-ray flares. Of the 168 bright sources in the ~1500 square degrees SPT-3G Main Field, eight have sufficient polarization signal-to-noise for reliable EVPA measurement, four of which have continuous Fermi LAT gamma-ray detections. We detect EVPA swings in all four gamma-ray-active blazars and in none of the remaining four, consistent with the established connection between EVPA swings and high-energy emission seen at optical wavelengths. The observed swing amplitudes and rotation rates are smaller than those found in optical studies, consistent with mm-wave EVPA variability being slower than at shorter wavelengths. Random walk simulations of the polarization angle using a multi-cell model fail to reproduce both the number, amplitude, and duration of the observed swings, suggesting that this mechanism is insufficient to explain the observed swing population. Analysis of EVPA variability on one-week timescales is consistent with the anti-correlation between polarization degree and EVPA rotation rate previously observed at optical wavelengths. Of the 14 detected swings, eight are within 30 days of a gamma-ray flare. While many individual swing-flare associations are found to have a very low probability of happening by chance, the full ensembles of 95 and 150 GHz swing time lags with respect to gamma-ray flares are found to be consistent with random coincidence.

Simpson, A. [Chicago U., Astron. Astrophys. Ctr.; ↗

Revisiting stochastic variability of AGNs with structure functions

Discrepancies between reported structure function (SF) slopes and their overall flatness as compared to the expectations from the damped random walk (DRW) model, which generally well describes the variability of active galactic nuclei (AGNs), have triggered us to study this problem in detail. We review common AGN variability observables and identify their most common problems. Equipped with this knowledge, we study ∼9000 r-band AGN light curves from Stripe 82 of the Sloan Digital Sky Survey, using SFs described by stochastic processes with the power exponential covariance matrix of the signal. We model the “subensemble” SFs in the redshift–absolute magnitude bins with the full SF equation (including the turnover and the noise part) and a single power law (SPL; in the “red noise regime” after subtracting the noise term). The distribution of full-equation SF (SPL) slopes peaks at γ=0.55±0.08 (0.52 ± 0.06) and is consistent with the DRW model. There is a hint of a weak correlation of γ with the luminosity and a lack of correlation with the black hole mass. The typical decorrelation timescale in the optical is τ=0.97±0.46 year. The SF amplitude at one year obtained from the SPL fitting is SF{sub 0}=0.22±0.06 mag and is overestimated because the SF is already at the turnover part, so the true value is SF{sub 0}=0.20±0.06 mag. The asymptotic variability is SF{sub ∞}=0.25±0.06 mag. It is strongly anticorrelated with both the luminosity and the Eddington ratio and is correlated with the black hole mass. The reliability of these results is fortified with Monte Carlo simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantum Adiabatic Optimization and Combinatorial Landscapes

In this paper we analyze the performance of the Quantum Adiabatic Evolution (QAE) algorithm on a variant of Satisfiability problem for an ensemble of random graphs parametrized by the ratio of clauses to variables, gamma = M / N. We introduce a set of macroscopic parameters (landscapes) and put forward an ansatz of universality for random bit flips. We then formulate the problem of finding the smallest eigenvalue and the excitation gap as a statistical mechanics problem. We use the so-called annealing approximation with a refinement that a finite set of macroscopic variables (verses only energy) is used, and are able to show the existence of a dynamic threshold gamma = gammad, beyond which QAE should take an exponentially long time to find a solution. We compare the results for extended and simplified sets of landscapes and provide numerical evidence in support of our universality ansatz.

Smelyanskiy, V. N.↗

Capacitated p -hub approach for park-and-ride facility location problem under nested logit demand function: polyhedral approaches

By generalizing the unconstrained p-hub approach for the park-and-ride (P&R) facility location problem under the multinomial logit demand function, the capacitated p-hub approach for the problem under the nested logit demand function captures a broader range of real-world cases. To solve this problem optimally, we introduce a mixed-integer linear program and accelerate its solution by enhancing the branch-and-cut procedure. To address the problem at a large scale, we introduce two other polyhedral approaches: variable neighborhood search (VNS) and adaptive randomized rounding (ARR). Downtown areas in Seoul have a high modal share of public transportation and congested road traffic, yet P&R has not been widely implemented. Therefore, we apply the ARR procedure to solve a real-world problem using traffic and geographic data from the Seoul metropolitan area. ARR performs better than VNS and addresses real-world cases. The solutions obtained by ARR present a phased expansion plan that encourages policymakers to start installing a small number of P&Rs immediately.

Capacitated p-hub approach↗

Temperature and Water Levels Collectively Regulate Methane Emissions From Subtropical Freshwater Wetlands

Abstract Wetlands are the largest and most climate‐sensitive natural sources of methane. Accurately estimating wetland methane emissions involves reconciling inversion (“top‐down”) and process‐based (“bottom‐up”) models within the global methane budget. However, estimates from these two model types are inherently interdependent and often reveal substantial discrepancies. To enhance the reliability of both approaches, we need a comprehensive understanding of wetland methane emissions and an independent high‐resolution long‐term flux data set. Here, we employed a data‐driven random forest approach to identify key variables influencing methane emissions from subtropical freshwater wetlands in the Southeastern United States. The model‐estimated monthly mean methane fluxes fit well with measured methane fluxes ( R 2 = 0.67) at four representative FLUXNET‐CH4 wetland sites across the region. Variable importance analysis highlighted the sensitivity of subtropical freshwater wetland methane emissions to variations in both temperature and water levels. High temperatures facilitate methanogenesis by enhancing microbial activities, while elevated water levels maintain anaerobic conditions necessary for methane production. Notably, the response of methane emissions to water level fluctuations is contingent on temperature conditions, and vice versa. Moreover, we constructed the first high‐spatial‐resolution (∼1 km × 1 km) and long‐term (1982–2010) gridded regional wetland methane flux product for the Southeastern United States, estimating annual methane emissions from subtropical freshwater wetlands in the region at 4.93 ± 0.11 Tg CH 4 yr −1 for 1982–2010. This new benchmark product holds promise for validating and parameterizing uncertain wetland methane emission processes in bottom‐up models and provides improved prior information for top‐down models.

He, Keqi [Earth and Climate Sciences Nicholas Scho↗