Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Charge conservation and higher moments of charge fluctuations

Higher moments of distributions of net charge and baryon number in heavy-ion collisions have been proposed as signals of fundamental QCD phase transitions. In order to better understand background processes for these observables, models are presented which enable one to gauge the effects of local charge conservation, decays of resonances and clusters, Bose symmetrization, and volume fluctuations. Monte Carlo methods for generating samplings of particles consistent with local charge conservation are presented and are followed by a review of simple analytic models involving a single type of charge with a constant experimental efficiency. The main model consists of thermal emission superimposed onto a simple parametrization of collective flow, known as a blast wave, with emission being consistent with individual canonical ensembles. The spatial extent of local charge conservation is parameterized by the size and extent over which charge is conserved. Here, the sensitivity of third- and fourth-order moments, skewness and kurtosis, to these parameters, and to beam energy and baryon density is explored. Comparisons with STAR data show that a significant part of the observed non-Poissonian fluctuations in net-proton fluctuations are explained by charge and baryon-number conservation, but that measurements of the STAR collaboration for fluctuations of net electric charge significantly differ from expectations of the models presented here.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Kinetic modeling of neutral transport for a continuum gyrokinetic code

In this work we present the first-of-its-kind coupling of a continuum full- f gyrokinetic turbulence model with a 6D continuum model for kinetic neutrals, carried out using the Gkeyll code. Our objective is to improve the first-principle understanding of the role of neutrals in plasma fueling, detachment, and their interaction with edge plasma profiles and turbulence statistics. Our model includes only atomic hydrogen and incorporates electron-impact ionization, charge exchange, and wall recycling. These features have been successfully verified with analytical predictions and benchmarked with the DEGAS2 Monte Carlo neutral code. We carry out simulations for a scrape-off layer (SOL) with simplified geometry and National Spherical Torus Experiment parameters. We compare these results to a baseline simulation without neutrals and find that neutral interactions reduce the normalized density fluctuation levels and associated skewness and kurtosis, while increasing auto-correlation times. A flatter density profile is also observed, similar to the SOL density shoulder formation in experimental scenarios with high fueling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Characteristic momentum of Hydro + and a bound on the enhancement of the speed of sound near the QCD critical point

Near the critical point in the QCD phase diagram, hydrodynamics breaks down at a momentum where the frequency of the fastest hydrodynamic mode becomes comparable with the decay rate of the slowest nonhydrodynamic mode. Hydro+ was developed as a framework which extends the range of validity of hydrodynamics beyond that momentum value. This was achieved through coupling the hydrodynamic modes to the slowest nonhydrodynamic mode. In this work, analyzing the spectrum of linear perturbations in Hydro+, we find that a slow mode falls out of equilibrium if its momentum is greater than a characteristic momentum value. That characteristic momentum turns out to be set by the branch points of the dispersion relations. These branch points occur at the critical momenta of so-called spectral curves and are related to the radius of convergence of the derivative expansion. The existence of such a characteristic momentum scale suggests that a particular class of slow modes has no remarkable effect on the flow of the plasma. Based on these results, we show that there is an enhancement of the speed of sound due to the Hydro+ corrections compared to the known thermodynamic corrections near the critical point. This enhancement can not exceed a temperature-dependent upper bound which we derive near the critical point in the QCD phase diagram.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

SU-E-T-399: Determination of the Radiobiological Parameters That Describe the Dose-Response Relations of Xerostomia and Disgeusia From Head and Neck Radiotherapy

Purpose: To estimate the radiobiological parameters that describe the doseresponse relations of xerostomia and disgeusia from head and neck cancer radiotherapy. To identify the organs that are best correlated with the manifestation of those clinical endpoints. Finally, to evaluate the goodnessof- fit by comparing the model predictions against the actual clinical results. Methods: In this study, 349 head and neck cancer patients were included. For each patient the dose volume histograms (DVH) of parotids (separate and combined), mandible, submandibular glands (separate and combined) and salivary glands were calculated. The follow-up of those patients was recorded at different times after the completion of the treatment (7 weeks, 3, 7, 12, 18 and 24 months). Acute and late xerostomia and acute disgeusia were the clinical endpoints examined. A maximum likelihood fitting was performed to calculate the best estimates of the parameters used by the relative seriality model. The statistical methods of the error distribution, the receiver operating characteristic (ROC) curve, the Pearson's test and the Akaike's information criterion were utilized to assess the goodness-of-fit and the agreement between the pattern of the radiobiological predictions with that of the clinical records. Results: The estimated values of the radiobiological parameters of salivary glands are D50 = 25.2 Gy, γ = 0.52, s = 0.001. The statistical analysis confirmed the clinical validity of those parameters (area under the ROC curve = 0.65 and AIC = 38.3). Conclusion: The analysis proved that the treatment outcome pattern of the patient material can be reproduced by the relative seriality model and the estimated radiobiological parameters. Salivary glands were found to have strong volume dependence (low relative seriality). Diminishing the biologically effective uniform dose to salivary glands below 30 Gy may significantly reduce the risk of complications to the patients irradiated for prostate cancer.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2004

Characterization of Transport Errors in Chemical Forecasts from a Global Tropospheric Chemical Transport Model

We propose a new methodology to characterize errors in the representation of transport processes in chemical transport models. We constrain the evaluation of a global three-dimensional chemical transport model (GEOS-CHEM) with an extended dataset of carbon monoxide (CO) concentrations obtained during the Transport and Chemical Evolution over the Pacific (TRACE-P) aircraft campaign. The TRACEP mission took place over the western Pacific, a region frequently impacted by continental outflow associated with different synoptic-scale weather systems (such as cold fronts) and deep convection, and thus provides a valuable dataset. for our analysis. Model simulations using both forecast and assimilated meteorology are examined. Background CO concentrations are computed as a function of latitude and altitude and subsequently subtracted from both the observed and the model datasets to focus on the ability of the model to simulate variability on a synoptic scale. Different sampling strategies (i.e., spatial displacement and smoothing) are applied along the flight tracks to search for systematic model biases. Statistical quantities such as correlation coefficient and centered root-mean-square difference are computed between the simulated and the observed fields and are further inter-compared using Taylor diagrams. We find no systematic bias in the model for the TRACE-P region when we consider the entire dataset (i.e., from the surface to 12 km ). This result indicates that the transport error in our model is globally unbiased, which has important implications for using the model to conduct inverse modeling studies. Using the First-Look assimilated meteorology only provides little improvement of the correlation, in comparison with the forecast meteorology. These general statements can be refined when the entire dataset is divided into different vertical domains, i.e., the lower troposphere (less than 2 km), the middle troposphere (2-6 km), and the upper troposphere (greater than 6 km). The best agreement between the observations and the model is found in the lower and middle troposphere. Downward displacements in the lower troposphere provide a better fit with the observed value, which could indicate a problem in the representation of boundary layer height in the model. Significant improvement is also found for downward and southward displacements in the upper troposphere. There are several potential sources of errors in our simulation of the continental outflow in the upper troposphere which could lead to such biases, including the location and/or the strength of deep convective cells as well as that of wildfires in Southeast Asia.

Bey, I.↗
NASA NTRS2019

X-ray Binary Luminosity Function Scaling Relations for Local Galaxies Based on Subgalactic Modeling

We present new Chandra constraints on the X-ray luminosity functions(XLFs)of X-ray binary(XRB)populations, as well as their scaling relations, for a sample of 38 nearby galaxies(D = 3.4–29 Mpc). Our galaxy sample is drawn primarily from the Spitzer Infrared Nearby Galaxies Survey(SINGS)and contains a wealth of Chandra(5.8 Ms total)and multiwavelength data, allowing for star formation rates(SFRs)and stellar masses(M⁎)to be measured on subgalactic scales. We divided the 2478 X-ray-detected sources into 21 subsamples in bins of specific SFR(sSFR ≡ SFR/M⁎)and constructed XLFs. To model the XLF dependence on sSFR, we fit a global XLF model, containing contributions from high-mass XRBs(HMXBs), low-mass XRBs(LMXBs), and background sources from the cosmic X-ray background that respectively scale with SFR, M⁎, and sky area. We find an HMXB XLF that is more complex in shape than previously reported and an LMXB XLF that likely varies with sSFR, potentially due to an age dependence. When applying our global model to XLF data for each individual galaxy, we discover a few galaxy XLFs that significantly deviate from our model beyond statistical scatter. Most notably, relatively low-metallicity galaxies have an excess of HMXBs above≈10^(38)erg/s, and elliptical galaxies that have relatively rich populations of globular clusters(GCs)show excesses of LMXBs compared to the global model. Additional modeling of how the XRB XLF depends on stellar age, metallicity, and GC specific frequency is required to sufficiently characterize the XLFs of galaxies.

Bret D. Lehmer↗
NASA NTRS2002

Scale Dependence of Statistics of Spatially Averaged Rain Rate Seen in TOGA COARE Comparison with Predictions from a Stochastic Model

A characteristic feature of rainfall statistics is that they in general depend on the space and time scales over which rain data are averaged. As a part of an earlier effort to determine the sampling error of satellite rain averages, a space-time model of rainfall statistics was developed to describe the statistics of gridded rain observed in GATE. The model allows one to compute the second moment statistics of space- and time-averaged rain rate which can be fitted to satellite or rain gauge data to determine the four model parameters appearing in the precipitation spectrum - an overall strength parameter, a characteristic length separating the long and short wavelength regimes and a characteristic relaxation time for decay of the autocorrelation of the instantaneous local rain rate and a certain 'fractal' power law exponent. For area-averaged instantaneous rain rate, this exponent governs the power law dependence of these statistics on the averaging length scale $L$ predicted by the model in the limit of small $L$. In particular, the variance of rain rate averaged over an $L \times L$ area exhibits a power law singularity as $L \rightarrow 0$. In the present work the model is used to investigate how the statistics of area-averaged rain rate over the tropical Western Pacific measured with ship borne radar during TOGA COARE (Tropical Ocean Global Atmosphere Coupled Ocean Atmospheric Response Experiment) and gridded on a 2 km grid depends on the size of the spatial averaging scale. Good agreement is found between the data and predictions from the model over a wide range of averaging length scales.

Kundu, Prasun K.↗
DOE OSTI2016

Performance Modeling of In Situ Rendering

With the push to exascale, in situ visualization and analysis will continue to play an important role in high performance computing. Tightly coupling in situ visualization with simulations constrains resources for both, and these constraints force a complex balance of trade-offs. A performance model that provides an a priori answer for the cost of using an in situ approach for a given task would assist in managing the trade-offs between simulation and visualization resources. In this work, we present new statistical performance models, based on algorithmic complexity, that accurately predict the run-time cost of a set of representative rendering algorithms, an essential in situ visualization task. To train and validate the models, we conduct a performance study of an MPI+X rendering infrastructure used in situ with three HPC simulation applications. We then explore feasibility issues using the model for selected in situ rendering questions.

Larsen, Matthew↗
DOE OSTI2013

Project of a Super Charm-Tau factory at the Budker Institute of Nuclear Physics in Novosibirsk

A project of a Super Charm-Tau factory is being developed at the Budker Institute of Nuclear Physics (Siberian Branch, Russian Academy of Sciences) in Novosibirsk. The electron-positron collider to be employed will operate at c.m. energies in the range between 2 and 5 GeV at an unprecedentedly high luminosity of 10{sup 35} cm{sup -2} s{sup -1} with a longitudinal electron polarization at the beam-interaction point. The main objective of experiments at the Super Charm-Tau factory is to study processes involving the production and properties of charmed quarks and tau leptons. A high luminosity of this setup will make it possible to obtain a statistical data sample that will be three to four orders of magnitude vaster than that from any other experiment performed thus far. Experiments at this setup are assumed to be sensitive to effects of new physics beyond the Standard Model. Investigations to be carried out at the Super-Charm-Tau factory will supplement future experiments at Super-B factories under construction in Italy and in Japan.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Light element ($$Z=1,2$$) production from spontaneous ternary fission of $$^{252}$$Cf

Abstract The yields of light elements ( $$Z=1,2$$ Z = 1 , 2 ) obtained from spontaneous ternary fission of $$^{252}$$ 252 Cf are treated within a nonequilibrium approach, and the contribution of unstable nuclei and excited bound states is taken into account. These light cluster yields may be used to probe dense matter, and to infer in-medium corrections such as Pauli blocking which is determined by the nucleon density. Continuum correlations are calculated from scattering phase shifts using the Beth-Uhlenbeck formula, and the effect of medium modification is estimated. The relevant distribution is reconstructed from the measured yields of isotopes. This describes the state of the nucleon system at scission and cluster formation, using only three Lagrange parameters which are the nonequilibrium counterparts of the temperature and chemical potentials, as defined in thermodynamic equilibrium. We concluded that a simple nuclear statistical equilibrium model neglecting continuum correlations and medium effects is not able to describe the measured distribution of H and He isotopes. Moreover, the freeze-out concept may serve as an important ingredient to the nonequilibrium approach using the relevant statistical operator concept.

Röpke, G.↗
NASA NTRS2023

Spacesuit Fit and Mobility Assessments by Digital Human Modeling

Spacesuit Fit and Mobility Assessments by Digital Human Modeling K. Han Kim (Leidos, Inc.) Elizabeth A. Benson (KBR, Inc.) Sudhakar L. Rajulu (NASA Johnson Space Center) Spacesuits are required to accommodate safe operations for astronauts across gender and a wide variety of body shapes and sizes. This goal has been of particular importance given the increasing diversity of NASA crewmembers for upcoming Missions. While testing with design prototypes is a critical step for spacesuit development, iterative mockup design, fabrication, and human subject tests can be extremely costly and time consuming. Moreover, testing with a limited subject pool has often raised questions for validity, as test subjects need to represent the entire astronaut population, not only of the past or current, but also the future. This study is aimed at demonstrating how digital human modeling (DHM) tools have been built and directly supported NASA spacesuit developments. With DHM, the computer aided design (CAD) model of a spacesuit was integrated with human body models. Two use scenarios are presented in this paper, namely fit and mobility. For fit assessments, the suit-to-body contact and compression patterns were estimated using 3D human body scan models virtually wearing a spacesuit model. A statistical fit classifier was made from the contact patterns and applied against a large database of body scans (N=2,500). With this technique, the NASA reference design spacesuit Exploration EVA Mobility Unit (xEMU)was verified to accommodate 90% of the astronaut-like population, with the critical dimensions covering 1stto 99thpercentiles of the target body measurements. A similar technique assessed XEMU mobility. The maximum reach envelopes were considered, within which the work objects and critical hardware interfaces should be located for safe and ergonomic operations. While the reach envelope geometry varies significantly with the suit wearer’s body size and strength, the existing test data did not include the subjects critically required to define suit mobility requirements, such as very small females or large males. Using the xEMU virtual model kinematically simulated and permuted for a hypothetical wearer, however, the existing data were statistically transformed and scaled. This method enabled for a parametric estimation of the reach envelopes from the 1stpercentile female or 99thpercentile male. The outcome was successfully incorporated for requirement developments. With the new DHM tools, human integration of the spacesuit was structurally simulated and predictively assessed, which would have been otherwise impossible. Also, the needs for iterative mockup and subject tests were significantly reduced, which resulted in time and cost saving. Additional work is in progress to integrate additional vehicles, tools, and hardware with DHM framework.

K Han Kim↗
DOE OSTI2020

Propagation of statistical uncertainties of Skyrme mass models to simulations of r -process nucleosynthesis

Uncertainties in nuclear models have a major impact on simulations that aim at understanding the origin of heavy elements in the universe through the rapid neutron capture process ( r process) of nucleosynthesis. Within the framework of the nuclear density functional theory, we use results of Bayesian statistical analysis to propagate uncertainties in the parameters of energy density functionals to the predicted r -process abundance pattern, by way not only of the nuclear masses but also through the influence of the masses on β -decay and neutron capture rates. We point out the importance of the nonequilibrium end stage of the r process in determining the width of the resulting abundance pattern uncertainty bands. We additionally make the first identifications of specific parameters of Skyrme-like energy density functionals which show tentative correlations with particular aspects of the r -process abundance pattern. While previous studies have explored the reduction in the abundance pattern uncertainties due to anticipated new measurements of neutron-rich nuclei, here we point out that an even larger reduction will occur when these new measurements are used to reduce the uncertainty of model predictions of masses, which are then propagated through to the abundance pattern. We make a quantitative prediction for how large this reduction will be.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Windowed multipole representation of R -matrix cross sections

Nuclear cross sections are basic inputs to any nuclear computation. Campaigns of experiments are fitted with the parametric R-matrix model of quantum nuclear interactions, and the resulting cross sections are documented—both pointwise and as resonance parameters (with uncertainties)—in standard evaluated nuclear data libraries (ENDF, JEFF, BROND, JENDL, CENDL, TENDL): these constitute our common knowledge of fundamental low-energy nuclear cross sections. In the past decade, a collaborative effort has been deployed to establish a new nuclear cross-section library format—the Windowed Multipole Library—with the goal of considerably reducing the computational cost of cross-section calculations in nuclear transport simulations. This work lays the theoretical foundations underpinning these efforts. From general R-matrix scattering theory, we derive the windowed multipole representation of nuclear cross sections. Though physically and mathematically equivalent to R-matrix cross sections, the windowed multipole representation is particularly well suited for subsequent temperature treatment of angle-integrated cross sections, in particular Doppler broadening, which is the averaging of cross sections over the thermal motion of the target atoms. Doppler broadening is of critical importance in neutron transport applications, as it ensures the stability of many nuclear reactors (negative thermal reactivity). Yet, Doppler broadening of nuclear cross sections has been a considerable bottleneck for nuclear transport computations, often requiring memory-costly pretabulations. We show that the windowed multipole representation can perform accurate Doppler broadening analytically (up to the first reaction threshold), from which we derive cross-section temperature derivatives to any order—all computable on the fly (without precalculations stored in memory). Furthermore, we here establish a way of converting the R-matrix resonance parameters uncertainty (covariance matrices) into windowed multipole parameters uncertainty. We show that generating stochastic nuclear cross sections by sampling from the resulting windowed multipole covariance matrix can reproduce the cross-section uncertainty in the original nuclear data file. The windowed multipole representation is therefore a novel nuclear physics formalism able to generate Doppler broadened stochastic nuclear cross sections on the fly, unlocking breakthrough computational gains for nuclear computations. Through this foundational paper, we hope to make the windowed multipole representation accessible, reproducible, and usable for the nuclear physics community, as well as provide the theoretical basis for future research on expanding its capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Simple, empirical approach to predict neutron capture cross sections from nuclear masses

Here, neutron capture cross sections are essential to understanding the astrophysical s and r processes, the modeling of nuclear reactor design and performance, and for a wide variety of nuclear forensics applications. Often, cross sections are needed for nuclei where experimental measurements are difficult. Enormous effort, over many decades, has gone into attempting to develop sophisticated statistical reaction models to predict these cross sections. Such work has met with some success but is often unable to reproduce measured cross sections to better than 40%, and has limited predictive power, with predictions from different models rapidly differing by an order of magnitude a few nucleons from the last measurement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Challenges in quantifying unparameterized spatial uncertainties in deep geologic repositories for nuclear waste

Spatially heterogeneous uncertainties are prevalent in geophysical modeling applications, such as probabilistic post-closure performance assessment (PA) of deep geologic repositories for nuclear waste. Such uncertainties are often highly influential to model outputs, so it is desirable to identify the most important mechanisms by which they influence model predictions. However, these uncertainties are often unparameterized in the sense that there is no set of parameters that can be specified to yield a particular realization of the uncertainty. Additionally, the uncertainty is not intrinsically endowed with a parameterization that captures a realization’s mechanistic influence on model outputs. Therefore, in this work we present a novel methodology to develop and assess a set of proxy variables that aim to represent this influence. We show how they can be computed, downselected, and incorporated into the construction of statistical surrogate models mapping model inputs to outputs. We present our methodology in the context of a motivating application problem in deep geologic repository PA and discuss the challenges in capturing the effects of these spatial heterogeneities in uncertainty analyses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

Dynamic structural disorder in supported nanoscale catalysts

We investigate the origin and physical effects of “dynamic structural disorder” (DSD) in supported nano-scale catalysts. DSD refers to the intrinsic fluctuating, inhomogeneous structure of such nano-scale systems. In contrast to bulk materials, nano-scale systems exhibit substantial fluctuations in structure, charge, temperature, and other quantities, as well as large surface effects. The DSD is driven largely by the stochastic librational motion of the center of mass and fluxional bonding at the nanoparticle surface due to thermal coupling with the substrate. Our approach for calculating and understanding DSD is based on a combination of real-time density functional theory/molecular dynamics simulations, transient coupled-oscillator models, and statistical mechanics. This approach treats thermal and dynamic effects over multiple time-scales, and includes bond-stretching and -bending vibrations, and transient tethering to the substrate at longer ps time-scales. Potential effects on the catalytic properties of these clusters are briefly explored. Model calculations of molecule-cluster interactions and molecular dissociation reaction paths are presented in which the reactant molecules are adsorbed on the surface of dynamically sampled clusters. This model suggests that DSD can affect both the prefactors and distribution of energy barriers in reaction rates, and thus can significantly affect catalytic activity at the nano-scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

Fluctuation statistics of outgoing longwave radiation in a general circulation model and in satellite data

Fluctuation statistics of Northern-Hemisphere-winter outgoing longwave radiation (OLR) obtained from data collected by Nimbus 7 and NOAA AVHRR are compared with those computed by a general circulation model (GCM). It was found that the GCM OLR has more variance than does the satellite, but that the GCM clear-sky OLR has less variance than the satellite OLR. The variance of the GCM outgoing full-sky OLR was found to be forced by clouds. While the standard deviation of the GCM full OLR is larger than that of the satellite OLR, the geographical patterns of both are similar. The autocorrelations of the modeled and the observed OLR also display similar geographical distributions, and they are in turn correlated with the dynamics.

Charlock, Thomas P.↗
DOE OSTI2022

Are Terrestrial Biosphere Models Fit for Simulating the Global Land Carbon Sink?

The Global Carbon Project estimates that the terrestrial biosphere has absorbed about one-third of anthropogenic CO 2 emissions during the 1959–2019 period. This sink-estimate is produced by an ensemble of terrestrial biosphere models and is consistent with the land uptake inferred from the residual of emissions and ocean uptake. The purpose of our study is to understand how well terrestrial biosphere models reproduce the processes that drive the terrestrial carbon sink. One challenge is to decide what level of agreement between model output and observation-based reference data is adequate considering that reference data are prone to uncertainties. To define such a level of agreement, we compute benchmark scores that quantify the similarity between independently derived reference data sets using multiple statistical metrics. Models are considered to perform well if their model scores reach benchmark scores. Our results show that reference data can differ considerably, causing benchmark scores to be low. Model scores are often of similar magnitude as benchmark scores, implying that model performance is reasonable given how different reference data are. While model performance is encouraging, ample potential for improvements remains, including a reduction in a positive leaf area index bias, improved representations of processes that govern soil organic carbon in high latitudes, and an assessment of causes that drive the inter-model spread of gross primary productivity in boreal regions and humid tropics. The success of future model development will increasingly depend on our capacity to reduce and account for observational uncertainties.

58 GEOSCIENCES↗