Search NASA⌕ Search

SEARCH · Search NASA

Results for “statistical models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,189 records · Page 66

Probabilistic Solar Power Forecasting Using Bayesian Model Averaging

There is rising interest in probabilistic forecasting to mitigate risks from solar power uncertainty, but the numerical weather prediction (NWP) ensembles readily available to system operators are often biased and underdispersed. We propose a Bayesian model averaging (BMA) post-processing method suitable for forecasting power from utility-scale photovoltaic (PV) plants at multiple time horizons up to at least the day-ahead timescale. BMA is a kernel dressing technique for NWP ensembles in which the forecast is a weighted sum of member-specific probability density functions. We tailor BMA for utility-scale PV forecasting by modeling power clipping at the AC inverter rating and advance the theory of BMA with a new beta kernel parameterization that accommodates theoretical constraints not previously addressed. BMA is demonstrated for a case study of 11 utility-scale PV plants in Texas, forecasting at hourly resolution for the complete year 2018. BMA's mixture-model approach mitigates underdispersion of the raw ensemble to significantly improve forecast calibration, while consistently outperforming an ensemble model output statistics (EMOS) parametric approach from the literature. At 4-hour lead time, the BMA post-processing achieves continuous ranked probability skill scores of 2--36% over the raw ensemble, with consistent performance at multiple lead times suitable for power system operations.

14 SOLAR ENERGY↗

Statistical mechanics in biology: how ubiquitous are long-range correlations?

The purpose of this opening talk is to describe examples of recent progress in applying statistical mechanics to biological systems. We first briefly review several biological systems, and then focus on the fractal features characterized by the long-range correlations found recently in DNA sequences containing non-coding material. We discuss the evidence supporting the finding that for sequences containing only coding regions, there are no long-range correlations. We also discuss the recent finding that the exponent alpha characterizing the long-range correlations increases with evolution, and we discuss two related models, the insertion model and the insertion-deletion model, that may account for the presence of long-range correlations. Finally, we summarize the analysis of long-term data on human heartbeats (up to 10(4) heart beats) that supports the possibility that the successive increments in the cardiac beat-to-beat intervals of healthy subjects display scale-invariant, long-range "anti-correlations" (a tendency to beat faster is balanced by a tendency to beat slower later on). In contrast, for a group of subjects with severe heart disease, long-range correlations vanish. This finding suggests that the classical theory of homeostasis, according to which stable physiological processes seek to maintain "constancy," should be extended to account for this type of dynamical, far from equilibrium, behavior.

Non-NASA Center↗

Optimizing Prediction Error for Time-dependent Solar Radiation Modeling

Numerical weather forecasting models and statistical methods have found wide use to help power companies estimate renewable output, but better methods are needed, particularly for extended forecasts. Machine learning approaches have been used here as well, but so far a major limitation is the ability to also predict the corresponding uncertainty in a forecast. Here we show that both can be done and demonstrate this using a long-term-short memory neural network where the difference between predicted and ground truth data are used to train a model for the corresponding forecast uncertainties.

97 MATHEMATICS AND COMPUTING↗

FACETS Dynamical Downscaling Simulations over North America by the CAM-MPAS Variable-Resolution Model

Comprehensive assessment of climate datasets created by statistical or dynamical models is important for effectively communicating model projection and associated uncertainty to stakeholders and decision-makers. The Department of Energy FACETS project aims to foster such communication through development of metrics and their demonstration on a hierarchy of downscaled climate datasets to quantify aspects of climate change projections that are credible, particularly for supporting decisions related to the energy-water-land nexus. As a part of this effort, we have produced a regional climate dataset using the Model for Prediction Across Scales coupled to the Community Atmosphere Model (CAM-MPAS). This global modeling framework is configured with variable-resolution meshes featuring higher resolutions over North America, as well as quasi-uniform resolution meshes across the globe. The variable-resolution configurations allow fine-scale features to be better resolved inside the refinement and interact with global-scale circulations. The dataset includes multiple uniform- (240km and 120km) and variable-resolution (200-50km, 100-25km, and 46-12km) simulations that are designed to be compatible with other regional climate simulations that contribute to the hierarchy of downscaled climate datasets of the project. Furthermore, the dataset consists of simulations for both the present-day (1989-2010) and future (2079-2100) climate and post-processing of the model output has been coordinated across the project for consistency to facilitate common analysis across the hierarchy of datasets. Altogether, this CAM-MPAS model dataset provides a unique opportunity to assess the influence of resolutions and modeling framework on model credibility and climate change projection.

Sakaguchi, Koichi↗

$\Sigma (1385)^{\pm }$ resonance production in Pb–Pb collisions at $\sqrt{s_{\textrm{NN}}} = 5.02$ TeV

Hadronic resonances are used to probe the hadron gas produced in the late stage of heavy-ion collisions since they decay on the same timescale, of the order of 1–10 fm/c, as the decoupling time of the system. In the hadron gas, (pseudo)elastic scatterings among the products of resonances that decayed before the kinetic freeze-out and regeneration processes counteract each other, the net effect depending on the resonance lifetime, the duration of the hadronic phase, and the hadronic cross sections at play. In this context, the Σ(1385) ± particle is of particular interest as models predict that regeneration dominates over rescattering despite its relatively short lifetime of about 5.5 fm/c. The first measurement of the Σ(1385) ± resonance production at midrapidity in Pb–Pb collisions at $\sqrt{s_{NN}}$=5.02 TeV with the ALICE detector is presented in this Letter. The resonances are reconstructed via their hadronic decay channel, Λπ, as a function of the transverse momentum (p T ) and the collision centrality. The results are discussed in comparison with the measured yield of pions and with expectations from the statistical hadronization model as well as commonly employed event generators, including PYTHIA8/Angantyr and EPOS3 coupled to the UrQMD hadronic cascade afterburner. None of the models can describe the data. For Σ(1385) ± , a similar behaviour as K*(892) 0 is observed in data unlike the predictions of EPOS3 with afterburner.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Theory of Alike Selectivity in Biological Channels

We introduce a statistical mechanical model of the selectivity filter that accounts for the interaction between ions within the channel and derive Eisenman equation of the filter selectivity directly from the condition of barrier-less conduction.

alike selectivity↗

Evaluation of the Maximum Soluble Water Recovery from Pretreated Urine Brines in Low Earth Orbit

The Urine Processor Assembly (UPA) on the International Space Station recovers water from pretreated urine by thin-film, vacuum compression distillation (VCD). A number of iterative improvements to the distiller assembly (DA) have enabled increases in on-orbit water recovery and improved distillate water quality at on-orbit recoveries up to 87% by volume relative to the volume of pretreated urine. The phosphoric acid-based pretreatment formulation has enabled these higher recoveries on-orbit since May of 2016. Application of simple empirical chemistry models have allowed water recovery of sulfuric-acid pretreated urine from the Russian Orbiting segment to be distilled with phosphoric acid pretreated urine from the U.S. Orbiting Segment. Returned liquid brine from brine filters in the UPA brine loop have enabled empirical measurement of gypsum solubility as a function of recovery, conductivity, or mass density in the real brines of ISS. Calculations of maximum soluble recovery, MSR, are based on preventing the precipitation of gypsum, which is the first mineral to precipitate for both types of pretreated urine. A statistical urine model of urine composition provides an estimate of worst-case urine from a crew of three during one year of operations on ISS. The MSR is defined based on the worst-case concentrations of urinary calcium and urinary sulfate that are pretreated, pooled, and distilled. The potential to increase water recoveries above 87% without gypsum precipitation are evaluated. The potential benefits and costs of the addition of a conductivity sensor or a Coriolis mass flow meter to the brine loop are evaluated.

Dean L Muirhead↗

Constraining the hadronic spectrum and repulsive interactions in a hadron resonance gas via fluctuations of conserved charges

We simultaneously incorporate two common extensions of the hadron resonance gas model, namely the addition of extra, unconfirmed resonances to the particle list and the excluded volume repulsive interactions. We emphasize the complementary nature of these two extensions and identify combinations of conserved charge susceptibilities that allow to constrain them separately. In particular, ratios of second-order susceptibilities like $\chi_{11}^{BQ}/\chi_2^B$ and $\chi_{11}^{BS}/\chi_2^B$ are sensitive only to the baryon spectrum, while fourth-to-second order ratios like $\chi_4^B/\chi_2^B$, $\chi_{31}^{BS}/\chi_{11}^{BS}$, or $\chi_{31}^{BQ}/\chi_{11}^{BQ}$ are mainly determined by repulsive interactions. Analysis of the available lattice results suggests the presence of both the extra states in the baryon-strangeness sector and the repulsive baryonic interaction, with indications that hyperons have a smaller repulsive core than non-strange baryons. Furthermore, the modified hadron resonance gas model presented here significantly improves the description of lattice QCD susceptibilities at chemical freeze-out and can be used for the analysis of event-by-event fluctuations in heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Towards Inner Magnetosphere Particle and Field Models

We report on the initial stages of an effort to construct comprehensive empirical models of the plasma and fields in the inner magnetosphere. The models are based not only on archival data (e.g., SCATHA) but also on current spacecraft mission data (e.g., POLAR). In this component of the effort, we incorporate the recently archived data provided by the SCATRA (Spacecraft Charging AT High Altitude) satellite. The SCAT14A satellite was in a near-geostationary orbit and was in operation for more than a decade. In this paper, we focus on the SCATRA plasma and magnetic field data from approximately the first two years of operation. The time-series data are binned according to spatial location and geomagnetic activity. Examples of statistical, empirical models from this initial effort are presented; even the simply-constructed preliminary models reveal such subtle features as the previously reported warm (greater than 100 eV) plasma density enhancement near the magnetic equator.

Spence, H. E.↗

A statistical equilibrium analysis of interstellar CO toward zeta Ophiuchi as recorded by the Goddard High Resolution Spectrograph

High-resolution and high signal-to-noise ratio echelle observations made by the Goddard High Resolution Spectrograph (GHRS), aboard the Hubble Space Telescope, of molecular (12)C(16)O and (13)C(16)O in the absorbing gas toward zeta Oph are reported. A statistical equilibrium model, incorporating a different way of dealing with the radiation field in the interstellar clouds, is presented. Based upon previous observations, two velocity components of the gas at V(sub LSR) = - 0.79 and 0.39 km/s, corresponding to solar velocity = -14.6 and -13.4 km/s, together with three different sets of velocity dispersions are adopted in our models. The ratios of the computed emergent intensities (I(sub em)) resulting from the emission lines in the lowest transitions (J = 1 goes to 0, 2 goes to 1, 3 goes to 2) of ground state of CO are presented. The ratios for the two velocity components, (I(sub em)(1 goes to 0)(sub 1))/(I(sub em)(1 goes to 0)(sub 2)), and (I(sub em)(2 goes to 1)(sub 1))/I(sub em)(1 goes to 0)(sub 1)) are comparable to the observed corrected antenna temperature ratios. From the solutions to the statistical equilibrium equations including the six lowest rotational transitions of CO and from multiple Voigt profile fits, the observed (12)CO and (13)CO UV abroption features are successfully reproduced. Profile fitting and double cloud curve-of-growth analysis using the GHRS and Copernicus data indicate that the abundance of CO is (1.8 +/- 0.2) x 10(exp 15)/sq cm. The most recently measured wavelength of the (6 goes to 0) (13)CO band by Haridass & Huber (R(0) = 1370.616 A) shows that the tentative (13)CO features in the GHRS data are unquestionably identified. If there is no blending from any other unknown interstellar species, the isotope ratio of (12)CO/(13)CO is approximately 82 +/- 25. This result is close to the terrestrial value of 90 and implies that the CO is not substantially fractionated, and (12)C(16)O experiences little preferential self-shielding through the selective isotopic photodissociation process. This conclusion differs strongly with that of Sheffer et al. Other implications of the results are discussed.

Lyu, Cheng-Hsuan↗

Sensitivity of Simulated Global Ocean Carbon Flux Estimates to Forcing by Reanalysis Products

Reanalysis products from MERRA, NCEP2, NCEP1, and ECMWF were used to force an established ocean biogeochemical model to estimate air-sea carbon fluxes (FCO2) and partial pressure of carbon dioxide (pCO2) in the global oceans. Global air-sea carbon fluxes and pCO2 were relatively insensitive to the choice of forcing reanalysis. All global FCO2 estimates from the model forced by the four different reanalyses were within 20% of in situ estimates (MERRA and NCEP1 were within 7%), and all models exhibited statistically significant positive correlations with in situ estimates across the 12 major oceanographic basins. Global pCO2 estimates were within 1% of in situ estimates with ECMWF being the outlier at 0.6%. Basin correlations were similar to FCO2. There were, however, substantial departures among basin estimates from the different reanalysis forcings. The high latitudes and tropics had the largest ranges in estimated fluxes among the reanalyses. Regional pCO2 differences among the reanalysis forcings were muted relative to the FCO2 results. No individual reanalysis was uniformly better or worse in the major oceanographic basins. The results provide information on the characterization of uncertainty in ocean carbon models due to choice of reanalysis forcing.

MERRA↗

Constraining the destruction rate of K 40 in stellar nucleosynthesis through the study of the Ar 40 ( p , n ) K 40 reaction

Background: We present that K 40 plays a significant role in the radiogenic heating of Earth-like exoplanets, which can affect the development of a habitable environment on their surfaces. The initial amount of K 40 in the interior of these planets depends on the composition of the interstellar clouds from which they formed. Within this context, nuclear reactions that regulate the production of K 40 during stellar evolution can play a critical role. Purpose: In this study, we constrain for the first time the astrophysical reaction rate of K 40 ( n , p ) Ar 40 , which is responsible for the destruction of K 40 during stellar nucleosynthesis. We provide to the nuclear physics community high-resolution data on the cross section and angular distribution of the Ar 40 ( p , n ) K 40 reaction. These are important to various applications involving Ar 40 . The associated reaction rate of the Ar 40 ( p , n ) K 40 process addresses a reaction rate gap in the Joint Institute for Nuclear Astrophysics REACLIB database in the region of intermediate-mass isotopes. Methods: We performed differential cross-section measurements on the Ar 40 ( p , n ) K 40 reaction, for six energies in the center-of-mass system between 3.2 and 4.0 MeV and various angles between 0 ° and 135 ° . The experiment took place at the Edwards Accelerator Laboratory at Ohio University using the beam swinger target location and a standard neutron time-of-flight technique. We extracted total and partial cross sections by integrating the double differential cross sections we measured. Results: The total and partial cross sections varied with energy due to the contribution from isobaric analog states and Ericson type fluctuations. The energy-averaged neutron angular distributions were symmetrical relative to 90 ° . Based on the experimental data, local transmission coefficients were extracted and were used to calculate the astrophysical reaction rates of Ar 40 ( p , n ) K 40 and K 40 ( n , p ) Ar 40 reactions. The new rates were found to vary significantly from the theoretical rates in the REACLIB library. We implemented the new rates in network calculations to study nucleosynthesis via the slow neutron capture process, and we found that the produced abundance of K 40 is reduced by up to 10% compared to calculations with the library rates. At the same time, the above result removes a significant portion of the previous theoretical uncertainty on the K 40 yields from stellar evolution calculations. Conclusions: Our results support a destruction rate of K 40 in massive stars via the K 40 ( n , p ) Ar 40 reaction that is larger compared to previous estimates. The rate of K 40 destruction via the K 40 ( n , p ) Ar 40 reaction now has a dramatically reduced uncertainty based on our measurement. Lastly, this result directly affects the predicted stellar yields of K 40 from nucleosynthesis, which is a critical input parameter for the galactic chemical evolution models that are currently employed for the study of significant properties of exoplanets.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Machine learning prediction of electron density and temperature from He I line ratios

We propose to utilize machine learning to predict the electron density, ne, and temperature, T e , from He I line intensity ratios. In this approach, training data consist of measured He I line ratios as input and ne and T e measured using other diagnostic(s) as desired output, which is a Langmuir probe in our study. Support vector machine regression analysis is, then, performed with the training data to develop a predictive model for n e and T e , separately. It is confirmed that n e and T e predicted using the developed models agree well with those from the Langmuir probe in the ranges of 0.28 × 10 18 ≤ n e (m -3 ) ≤ 3.8 × 10 18 and 3.2 ≤ T e (eV) ≤ 7.5. The developed models are, further, examined with an evaluation data, which are not included in the training data, and are found to well reproduce absolute values and radial profiles of probe-measured n e and T e .

47 OTHER INSTRUMENTATION↗

Generation of Fragment Angular Momentum in Fission

A recent analysis of experimental data [J. Wilson et al., Nature (London) 590, 566 (2021)] found that the angular momenta of nuclear fission fragments are uncorrelated. Based on this finding, the authors concluded that the spins are therefore determined only after scission has occurred. We show here that the nucleon-exchange mechanism, as implemented in the well-established event-by-event fission model freya, while agitating collective rotational modes in which the two spins are highly correlated, nevertheless leads to fragment spins that are largely uncorrelated. We find that this counterexample invalidates the conclusion in [J. Wilson et al.] that uncorrelated spins must necessarily have been generated after scission (a potentious conclusion that would rule out all models that generate the fragment spins prior to scission). Furthermore, it was reported [J. Wilson et al.] that the mass dependence of the average fragment spin has a sawtooth structure. We demonstrate that such a behavior naturally emerges when shell and deformation effects are included in the moments of inertia of the fragments at scission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modified Brink-Axel hypothesis for astrophysical Gamow-Teller transitions

Weak interaction charged current transition strengths from highly excited nuclear states are fundamental ingredients for accurate modeling of compact object composition and dynamics, but are difficult to obtain either from experiment or theory. For lack of alternatives, calculations have often fallen back upon a generalized Brink-Axel hypothesis, that is, assuming the strength function (transition probability) is independent of the initial nuclear state but depends only upon the transition energy and the weak interaction properties of the parent nucleus ground state. Here we present numerical evidence for a modified “local” Brink-Axel hypothesis for Gamow-Teller transitions for pf-shell nuclei relevant to astrophysical applications. Specifically, while the original Brink-Axel hypothesis does not hold globally, strength functions from initial states nearby in energy are similar within statistical fluctuations. This agrees with previous work on strength function moments. Using this modified hypothesis, we can tackle strength functions at previously intractable initial energies, using semiconverged initial states at arbitrary excitation energy. Our work provides a well-founded method for computing accurate thermal weak transition rates for medium-mass nuclei at temperatures occurring in stellar cores near collapse. We finish by comparing results to previous calculations of astrophysical rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dynamic density functional theory of polymers with salt in electric fields

Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Remote sensing of earth terrain

Two monographs and 85 journal and conference papers on remote sensing of earth terrain have been published, sponsored by NASA Contract NAG5-270. A multivariate K-distribution is proposed to model the statistics of fully polarimetric data from earth terrain with polarizations HH, HV, VH, and VV. In this approach, correlated polarizations of radar signals, as characterized by a covariance matrix, are treated as the sum of N n-dimensional random vectors; N obeys the negative binomial distribution with a parameter alpha and mean bar N. Subsequently, and n-dimensional K-distribution, with either zero or non-zero mean, is developed in the limit of infinite bar N or illuminated area. The probability density function (PDF) of the K-distributed vector normalized by its Euclidean norm is independent of the parameter alpha and is the same as that derived from a zero-mean Gaussian-distributed random vector. The above model is well supported by experimental data provided by MIT Lincoln Laboratory and the Jet Propulsion Laboratory in the form of polarimetric measurements.

Kong, J. A.↗

Modeling Buried Object Brightness and Visibility for Ground Penetrating Radar

Comparing the observed brightness of various buried objects is a straightforward way to characterize the performance of a ground penetrating radar (GPR) system. However, a limitation arises. A simple comparison of buried object brightness values does not disentangle the effects of the GPR system itself from the system's operating environment and the objects being observed. Therefore, with brightness values exhibiting an unknown synthesis of systemic, environmental, and object factors, GPR system analysis becomes a convoluted affair. In this work, we use an experimentally collected dataset of over 25,000 object observations from five different multistatic radar arrays to develop models of buried object brightness and control for these various effects. Additionally, our modeling efforts provide a means for quantifying the relative brightness of GPR systems, the objects they detect, and the physical properties of those objects that influence observed brightness. To evaluate the models' performance on new object observations, we repeatedly simulate fitting them to half the dataset and predicting the observed brightness values of the unseen half. In addition, we introduce a method for estimating the probability that individual observations constitute a visible object, which aids in failure analysis, performance characterization, and dataset cleaning.

97 MATHEMATICS AND COMPUTING↗