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 1,063 records · Page 59

Trojan horse method as an indirect approach to study resonant reactions in nuclear astrophysics

The primary goal of the Trojan horse method (THM) is to analyze resonant rearrangement reactions when the density of the resonance levels is low and statistical models cannot be applied. The main difficulty of the analysis is related with the facts that in the final state the THM reaction involves three particles and that the intermediate particle, which is transferred from the Trojan horse particle to the target nucleus to form a resonance state, is virtual. Another difficulty is associated with the Coulomb interaction between the particles, especially, taking into account that the goal of the THM is to study resonant rearrangement reactions at very low energies important for nuclear astrophysics. The exact theory of such reactions with three charged particles is very complicated and is not available. This is why different approximations are used to analyze THM reactions. In this review paper we describe a new approach based on a few-body formalism that provides a solid basis for deriving the THM reaction amplitude taking into account rescattering of the particles in the initial, intermediate and final states of the THM reaction. Since the THM uses a two-step reaction in which the first step is the transfer reaction populating a resonance state, we address the theory of the transfer reactions. The theory is based on the surface-integral approach and R-matrix formalism. We also discuss application of the THM to resonant reactions populating both resonances located on the second energy sheet and subthreshold resonances, which are subthreshold bound states located at negative energies close to thresholds. We consider the application of the THM to determine the astrophysical factors of resonant radiative-capture reactions at energies so low that direct measurements can hardly be performed due to the negligibly small penetrability factor in the entry channel of the reaction. Here, we elucidated the main ideas of the THM and outline necessary conditions to perform the THM experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraining the nuclear spin distribution using improved 197 Au neutron resonance parameters

New neutron transmission data at resonance energies using a 197 Au sample were measured using an early version of the Device for Indirect Capture Experiments on Radionuclides (DICER), which is under development at the Los Alamos Neutron Science Center (LANSCE). These data were combined with previous neutron transmission and capture data in a simultaneous R-matrix analysis to extract improved neutron resonance parameters for this nuclide. As a result, total radiation widths, Γ γ , were obtained for 33 J=1 and 44 J=2 197 Au+n resonances. Γ γ distributions for these two spins states were compared to distributions calculated according to the nuclear statistical model using published nuclear level density (NLD) and photon strength functions (PSF) measured using the Oslo technique. The calculated distributions were found to be narrower and the average values for the two spins states closer together than the data. The calculation can be brought into agreement with the data by substantial modifications to the spin distribution in 198 Au as a function of excitation energy. As far as we know, the spin distribution currently is otherwise poorly constrained. The modified spin distribution changes the shapes of the NLD and PSF extracted using the Oslo technique and so could have broad implications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurements of ${\Xi \left( 1530\right) ^{0}} $ and ${\overline{\Xi }\left( 1530\right) ^{0}} $ production in proton–proton interactions at $\sqrt{s_{NN}}$ = 17.3 $\text{ GeV }$ in the NA61/SHINE experiment

Double-differential yields of ${\Xi \left( 1530\right) ^{0}} $ and ${\overline{\Xi }\left( 1530\right) ^{0}} $ resonances produced in p+p interactions were measured at a laboratory beam momentum of 158 $\text{ GeV }\!/\!c$. This measurement is the first of its kind in p+p interactions below LHC energies. It was performed at the CERN SPS by the NA61/SHINE collaboration. Double-differential distributions in rapidity and transverse momentum were obtained from a sample of $26\times 10^6$ inelastic events. The spectra are extrapolated to full phase space resulting in mean multiplicity of ${\Xi \left( 1530\right) ^{0}} $ ($6.73 \pm 0.25\pm 0.67)\times 10^{-4}$ and ${\overline{\Xi }\left( 1530\right) ^{0}} $ ($2.71 \pm 0.18\pm 0.18)\times 10^{-4}$. The rapidity and transverse momentum spectra and mean multiplicities were compared to predictions of string-hadronic and statistical model calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Expanding the access of wearable silicone wristbands in community-engaged research through best practices in data analysis and integration

Wearable silicone wristbands are a rapidly growing exposure assessment technology that offer researchers the ability to study previously inaccessible cohorts and have the potential to provide a more comprehensive picture of chemical exposure within diverse communities. However, there are no established best practices for analyzing the data within a study or across multiple studies, thereby limiting impact and access of these data for larger meta-analyses. We utilize data from three studies, from over 600 wristbands worn by participants in New York City and Eugene, Oregon, to present a first-of-its-kind manuscript detailing wristband data properties. We further discuss and provide concrete examples of key areas and considerations in common statistical modeling methods where best practices must be established to enable meta-analyses and integration of data from multiple studies. Finally, we detail important and challenging aspects of machine learning, meta-analysis, and data integration that researchers will face in order to extend beyond the limited scope of individual studies focused on specific populations.

Bramer, Lisa M.↗

Density-Aware Clustering Based on Aggregated Heat Kernel and Its Transformation

Current spectral clustering algorithms suffer from the sensitivity to existing noise, and parameter scaling, and may not be aware of different density distributions across clusters. If these problems are left untreated, the consequent clustering results cannot accurately represent true data patterns, in particular, for complex real world datasets with heterogeneous densities. This paper aims to solve these problems by proposing a diffusion-based Aggregated Heat Kernel (AHK) to improve the clustering stability, and a Local Density Affinity Transformation (LDAT) to correct the bias originating from different cluster densities. AHK statistically\ models the heat diffusion traces along the entire time scale, so it ensures robustness during clustering process, while LDAT probabilistically reveals local density of each instance and suppresses the local density bias in the affinity matrix. Our proposed framework integrates these two techniques systematically. As a result, not only does it provide an advanced noise-resisting and density-aware spectral mapping to the original dataset, but also demonstrates the stability during the processing of tuning the scaling parameter (which usually controls the range of neighborhood). Furthermore, our framework works well with the majority of similarity kernels, which ensures its applicability to many types of data and problem domains. The systematic experiments on different applications show that our proposed algorithms outperform state-of-the-art clustering algorithms for the data with heterogeneous density distributions, and achieve robust clustering performance with respect to tuning the scaling parameter and handling various levels and types of noise.

97 MATHEMATICS AND COMPUTING↗

Bimolecular Chemistry in the Ultracold Regime

Advances in atomic, molecular, and optical physics techniques allowed the cooling of simple molecules down to the ultracold regime ([Formula: see text]1 mK) and opened opportunities to study chemical reactions with unprecedented levels of control. This review covers recent developments in studying bimolecular chemistry at ultralow temperatures. We begin with a brief overview of methods for producing, manipulating, and detecting ultracold molecules. We then survey experimental works that exploit the controllability of ultracold molecules to probe and modify their long-range interactions. Further combining the use of physical chemistry techniques such as mass spectrometry and ion imaging significantly improved the detection of ultracold reactions and enabled explorations of their dynamics in the short range. We discuss a series of studies on the reaction KRb + KRb → K 2 + Rb 2 initiated below 1 μK, including the direct observation of a long-lived complex, the demonstration of product rotational state control via conserved nuclear spins, and a test of the statistical model using the complete quantum state distribution of the products.

Chemistry↗

Code for the manuscript "Lagrangian Attention Tensor Networks for Velocity Gradient Statistical Mode

We disclose a python/pytorch implementation of the physics-informed machine learning algorithm described in "Lagrangian Attention Tensor Networks for Velocity Gradient Statistical Modeling", LA-UR-24-30678. Direct numerical simulation (DNS) of ubiquitous turbulence phenomena is computationally infeasible for realistic flows. As a result, reduced modeling for turbulent flows aim to reduce the number of resolved scales while retaining accurate representations of the small-scale physics. The dynamics of the velocity gradient tensor (VGT) is a key ingredient in reduced or subgrid turbulence models. The evolution equation for the VGT involves nonlocal terms, requiring closure modeling. This implementation of the novel methodology of Lagrangian Attention Tensor Networks (LATN), utilizes a structured representation of the history of the VGT to inform a physics-informed machine learning algorithm. This addition of structured memory terms is shown to outperform previous models when trained and evaluated on DNS data.

Livescu, Daniel [LANL]↗

Coding and noncoding variants in EBF3 are involved in HADDS and simplex autism

Abstract Background Previous research in autism and other neurodevelopmental disorders (NDDs) has indicated an important contribution of protein-coding (coding) de novo variants (DNVs) within specific genes. The role of de novo noncoding variation has been observable as a general increase in genetic burden but has yet to be resolved to individual functional elements. In this study, we assessed whole-genome sequencing data in 2671 families with autism (discovery cohort of 516 families, replication cohort of 2155 families). We focused on DNVs in enhancers with characterized in vivo activity in the brain and identified an excess of DNVs in an enhancer named hs737. Results We adapted the fitDNM statistical model to work in noncoding regions and tested enhancers for excess of DNVs in families with autism. We found only one enhancer (hs737) with nominal significance in the discovery (p = 0.0172), replication (p = 2.5 × 10 −3 ), and combined dataset (p = 1.1 × 10 −4 ). Each individual with a DNV in hs737 had shared phenotypes including being male, intact cognitive function, and hypotonia or motor delay. Our in vitro assessment of the DNVs showed they all reduce enhancer activity in a neuronal cell line. By epigenomic analyses, we found that hs737 is brain-specific and targets the transcription factor gene EBF3 in human fetal brain. EBF3 is genome-wide significant for coding DNVs in NDDs (missense p = 8.12 × 10 −35 , loss-of-function p = 2.26 × 10 −13 ) and is widely expressed in the body. Through characterization of promoters bound by EBF3 in neuronal cells, we saw enrichment for binding to NDD genes (p = 7.43 × 10 −6 , OR = 1.87) involved in gene regulation. Individuals with coding DNVs have greater phenotypic severity (hypotonia, ataxia, and delayed development syndrome [HADDS]) in comparison to individuals with noncoding DNVs that have autism and hypotonia. Conclusions In this study, we identify DNVs in the hs737 enhancer in individuals with autism. Through multiple approaches, we find hs737 targets the gene EBF3 that is genome-wide significant in NDDs. By assessment of noncoding variation and the genes they affect, we are beginning to understand their impact on gene regulatory networks in NDDs.

59 BASIC BIOLOGICAL SCIENCES↗

General Bayesian Inference over the Stiefel Manifold via the Givens Transform

We introduce here an approach based on the Givens representation for posterior inference in statistical models with orthogonal matrix parameters, such as factor models and probabilistic principal component analysis (PPCA). We show how the Givens representation can be used to develop practical methods for transforming densities over the Stiefel manifold into densities over subsets of Euclidean space. We demonstrate how to deal with issues arising from the topology of the Stiefel manifold and how to inexpensively compute the change-of-measure terms. We introduce an auxiliary parameter approach that limits the impact of topological issues. We provide both analysis of our methods and numerical examples demonstrating the effectiveness of the approach. We also discuss how our Givens representation can be used to define general classes of distributions over the space of orthogonal matrices. We then give demonstrations on several examples showing how the Givens approach performs in practice in comparison with other methods.

97 MATHEMATICS AND COMPUTING↗

Godiva II—An Unmoderated Pulse-Irradiation Reactor

Design features of Godiva II, the improved pulsed reactor successor to Lady Godiva, are discussed together with characteristics of power excursions, and performance is compared with that of the original Godiva. Measurements of the wait time between stepwise reactivity insertion and the occurrence of a burst are presented and compared with theory based on a statistical model of fission chains. Analytical and numerical solutions of the reactor equations are developed to reproduce experimental data and extrapolate to higher energy release. Here, consideration is also given to perturbations arising from room-returned neutrons. Two different modes of operation are discussed and some design problems of Godiva-type pulsed reactors are briefly mentioned. Typical bursts are illustrated with peak powers up to 13,000 Mw and widths at half-maximum down to 35 μsec.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Limb-Enhancer Genie: An accessible resource of accurate enhancer predictions in the developing limb

Epigenomic mapping of enhancer-associated chromatin modifications facilitates the genome-wide discovery of tissue-specific enhancers in vivo. However, reliance on single chromatin marks leads to high rates of false-positive predictions. More sophisticated, integrative methods have been described, but commonly suffer from limited accessibility to the resulting predictions and reduced biological interpretability. Here we present the Limb-Enhancer Genie (LEG), a collection of highly accurate, genome-wide predictions of enhancers in the developing limb, available through a user-friendly online interface. We predict limb enhancers using a combination of > 50 published limb-specific datasets and clusters of evolutionarily conserved transcription factor binding sites, taking advantage of the patterns observed at previously in vivo validated elements. By combining different statistical models, our approach outperforms current state-of-the-art methods and provides interpretable measures of feature importance. Our results indicate that including a previously unappreciated score that quantifies tissue-specific nuclease accessibility significantly improves prediction performance. We demonstrate the utility of our approach through in vivo validation of newly predicted elements. Moreover, we describe general features that can guide the type of datasets to include when predicting tissue-specific enhancers genome-wide, while providing an accessible resource to the general biological community and facilitating the functional interpretation of genetic studies of limb malformations.

59 BASIC BIOLOGICAL SCIENCES↗

Assessing Inhalation Exposures Associated with Contamination Events in Water Distribution Systems

When a water distribution system (WDS) is contaminated, short-term inhalation exposures to airborne contaminants could occur as the result of domestic water use. The most important domestic sources of such exposures are likely to be showering and the use of aerosol-producing humidifiers, i.e., ultrasonic and impeller (cool-mist) units. A framework is presented for assessing the potential effects of short-term, system-wide inhalation exposures that could result from such activities during a contamination event. This framework utilizes available statistical models for showering frequency and duration, available exposure models for showering and humidifier use, and experimental results on both aerosol generation and the volatilization of chemicals during showering. New models for the times when showering occurs are developed using time-use data for the United States. Given a lack of similar models for how humidifiers are used, or the information needed to develop them, an analysis of the sensitivity of results to assumptions concerning humidifier use is presented. The framework is applied using network models for three actual WDSs. Simple models are developed for estimating upper bounds on the potential effects of system-wide inhalation exposures associated with showering and humidifier use. From a system-wide, population perspective, showering could result in significant inhalation doses of volatile chemical contaminants, and humidifier use could result in significant inhalation doses of microbial contaminants during a contamination event. From a system-wide perspective, showering is unlikely to be associated with significant doses of microbial contaminants. In conclusion, given the potential importance of humidifiers as a source of airborne contaminants during a contamination event, an improved understanding of the nature of humidifier use is warranted.

54 ENVIRONMENTAL SCIENCES↗

Using spatio-temporal graph neural networks to estimate fleet-wide photovoltaic performance degradation patterns

Accurate estimation of photovoltaic (PV) system performance is crucial for determining its feasibility as a power generation technology and financial asset. PV-based energy solutions offer a viable alternative to traditional energy resources due to their superior Levelized Cost of Energy (LCOE). A significant challenge in assessing the LCOE of PV systems lies in understanding the Performance Loss Rate (PLR) for large fleets of PV systems. Estimating the PLR of PV systems becomes increasingly important in the rapidly growing PV industry. Precise PLR estimation benefits PV users by providing real-time monitoring of PV module performance, while explainable PLR estimation assists PV manufacturers in studying and enhancing the performance of their products. However, traditional PLR estimation methods based on statistical models have notable drawbacks. Firstly, they require user knowledge and decision-making. Secondly, they fail to leverage spatial coherence for fleet-level analysis. Additionally, these methods inherently assume the linearity of degradation, which is not representative of real world degradation. To overcome these challenges, we propose a novel graph deep learning-based decomposition method called the Spatio-Temporal Graph Neural Network for fleet-level PLR estimation (PV-stGNN-PLR). PV-stGNN-PLR decomposes the power timeseries data into aging and fluctuation components, utilizing the aging component to estimate PLR. PV-stGNN-PLR exploits spatial and temporal coherence to derive PLR estimation for all systems in a fleet and imposes flatness and smoothness regularization in loss function to ensure the successful disentanglement between aging and fluctuation. We have evaluated PV-stGNN-PLR on three simulated PV datasets consisting of 100 inverters from 5 sites. Experimental results show that PV-stGNN-PLR obtains a reduction of 33.9% and 35.1% on average in Mean Absolute Percent Error (MAPE) and Euclidean Distance (ED) in PLR degradation pattern estimation compared to the state-of-the-art PLR estimation methods.

14 SOLAR ENERGY↗

QA/QC-ed BCI meteorological drivers

Time series raw datasets of Temperature, Solar Radiation, Relative Humidity, Rainfall, Barometric Pressure, Wind Speed, and Wind Direction for the Barro Colorado (PA-BCI) site for the period from 2003-01-01 to 2016-12-31 were downloaded from the STRI website. The QA/QC protocol of raw datasets included the detection and removal of NAs, outliers and bad data points, imputation of missing data, and creation of the time series with equidistant time stamps. The time series for different meteorological parameters were aligned into a single database to be used for modeling. Statistical QA/QC analysis was performed using a series of R libraries in Rstudio. The changes in the dataset were flagged. The file BCI_met_drivers_2003-2016_QAQC_summary_report.docx provides the details of the QA/QC data analysis.

54 ENVIRONMENTAL SCIENCES↗

Preface to the Focus Section on Machine Learning in Seismology

Machine learning (ML) is a collection of algorithms and statistical models that enable computers to extract relevant patterns and information from large data sets. Unlike physical modeling approaches, in which scientists develop theories based on physical laws to compare with real–world observations, ML approaches learn directly from data without explicitly reasoning about the underlying physical mechanisms. ML algorithms are often categorized into supervised and unsupervised learning (see fig. 2 in Kong et al., 2018). Supervised learning algorithms build a model from existing labeled input data with the goal of predicting the labels of new data.

58 GEOSCIENCES↗

Stan : A Probabilistic Programming Language

Stan is a probabilistic programming language for specifying statistical models. A Stan program imperatively defines a log probability function over parameters conditioned on specified data and constants. As of version 2.14.0, Stan provides full Bayesian inference for continuous-variable models through Markov chain Monte Carlo methods such as the No-U-Turn sampler, an adaptive form of Hamiltonian Monte Carlo sampling. Penalized maximum likelihood estimates are calculated using optimization methods such as the limited memory Broyden-Fletcher-Goldfarb-Shanno algorithm. Stan is also a platform for computing log densities and their gradients and Hessians, which can be used in alternative algorithms such as variational Bayes, expectation propagation, and marginal inference using approximate integration. To this end, Stan is set up so that the densities, gradients, and Hessians, along with intermediate quantities of the algorithm such as acceptance probabilities, are easily accessible. Stan can also be called from the command line using the cmdstan package, through R using the rstan package, and through Python using the pystan package. All three interfaces support sampling and optimization-based inference with diagnostics and posterior analysis. rstan and pystan also provide access to log probabilities, gradients, Hessians, parameter transforms, and specialized plotting.

97 MATHEMATICS AND COMPUTING↗

The Installation of Direct Water-Cooling Systems to Reduce Cooling Energy Requirements for High-Performance Computing Centers

A large cluster of High-Performance Computing (HPC) equipment at the Lawrence Livermore National Laboratory in California was retrofitted with an Asetek cooling system. The Asetek system is a hybrid scheme with water-cooled cold plates on high-heat-producing components in the information technology (IT) equipment, and with the remainder of the heat being removed by conventional air-cooling systems. In order to determine energy savings of the Asetek system, data were gathered and analyzed two ways: using top-down statistical models, and bottom-up engineering models. The cluster, “Cabernet”, rejected its heat into a facilities cooling water loop which in turn rejected the heat into the same chilled water system serving the computer-room air handlers (CRAHs) that provided the air-based cooling for the room. Because the “before” and “after” cases both reject their heat into the chilled water system, the only savings is due to reduction in CRAH fan power. The top-down analysis showed a 4% overall energy savings for the data center (power usage effectiveness (PUE) —the ratio of total data center energy to IT energy— dropped from 1.60 to 1.53, lower is better); the bottom-up analysis showed a 3% overall energy savings (PUE from 1.70 to 1.66) and an 11% savings for the Cabernet system by itself (partial PUE of 1.51). Greater savings, on the order of 15-20%, would be possible if the chilled water system was not used for rejecting the heat from the Asetek system. About 37% of the heat from the Cab system was rejected to the cooling water, lower than at other installations.

Earni, Shankar↗

Fundamental Understanding of Ambient and High-Temperature Plasticity Phenomena in Structural Materials in Advanced Reactors

The goal of this research project is to develop the methods and tools necessary to link unit processes analyzed using atomistic simulations involving interaction of vacancies and interstitials with dislocations, as well as dislocation mediation at sessile junctions and interfaces as affected by radiation, with cooperative influence on higher-length scale behavior of polycrystals. These tools and methods are necessary to design and enhance radiation-induced damage-tolerant alloys. The project will achieve this goal by applying atomistic simulations to characterize unit processes of: 1. Dislocation nucleation, absorption, and desorption at interfaces 2. Vacancy production, radiation-induced segregation of substitutional Cr at defect clusters (point defect sinks) in BCC Fe-Cr ferritic/martensitic steels 3. Investigation of interaction of interstitials and vacancies with impurities (V, Nb, Ta, Mo, W, Al, Si, P, S) 4. Time evolution of swelling (cluster growth) phenomena of irradiated materials 5. Energetics and kinetics of dislocation bypass of defects formed by interstitial clustering and formation of prismatic loops, informing statistical models of continuum character with regard to processes of dislocation glide, vacancy agglomeration and swelling, climb and cross slip This project will consider the Fe, Fe-C, and Fe-Cr ferritic/martensitic material system, accounting for magnetism by choosing appropriate interatomic potentials and validating with first principles calculations. For these alloys, the rate of swelling and creep enhancement is considerably lower than that of face-centered cubic (FCC) alloys and of austenitic Fe-Cr-Mo alloys. The team will confirm mechanisms, validate simulations at various time and length scales, and improve the veracity of computational models. The proposed research's feasibility is supported by recent modeling of radiation effects in metals and alloys, interfacial dislocation transfer reactions in nano-twinned copper, and dislocation reactions at general boundaries, along with extensive modeling cooperative effects of dislocation interactions and migration in crystals and polycrystals using continuum models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗