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 775 records · Page 43

Landscape‐scale lidar analysis of aboveground biomass distribution in secondary Brazilian Atlantic Forest

Abstract Secondary forests account for more than half of tropical forests and represent a growing carbon sink, but rates of biomass accumulation vary by a factor of two or more even among plots in the same landscape. To better understand the drivers of this variability, we used airborne lidar to measure forest canopy height and estimate biomass over 4529 ha at Serra do Conduru Park in Southern Bahia, Brazil. We measured trees in 30 georeferenced field plots (0.25‐ha each) to estimate biomass using allometry. Then we estimated aboveground biomass density ( ABD ) across the lidar study area using a statistical model developed from our field plots. This model related the 95th percentile of the distribution of lidar return heights to ABD . We overlaid this map of ABD on a Landsat‐derived forest age map to determine rates of biomass accumulation. We found rapid initial biomass regeneration (~6 Mg/ha yr), which slowed as forests aged. We also observed high variability in both height and biomass across the landscape within forests of similar age. Nevertheless, a regression model that accounted for spatial autocorrelation and included forest age, slope, and distance to roads or open areas explained 62 and 77 percent of the landscape variation in ABD and canopy height, respectively. Thus, while there is high spatial heterogeneity in forest recovery, and the drivers of this heterogeneity warrant further investigation, we suggest that a relatively simple set of predictor variables is sufficient to explain the majority of variance in both height and ABD in this landscape.

Becknell, Justin M.↗

Call Center Arrivals: When to Jointly Forecast Multiple Streams?

We consider call centers that have multiple (potentially inter‐dependent) demand arrival streams. Workforce management of such labor intensive service systems starts with forecasting future arrival demand. We investigate the question of whether and when to jointly forecast future arrivals of the multiple streams. We first develop a general statistical model to simultaneously forecast multi‐stream arrival rates. The model takes into account three types of inter‐stream dependence. We then show with analytical and simulation studies how the forecasting benefits of the multi‐stream forecasting model vary by the type, direction, and strength of inter‐stream dependence. In particular, we find that it is beneficial to simultaneously forecast multi‐stream arrivals (instead of separately forecasting each stream), when there exists inter‐stream lag dependence among daily arrival rates. Empirical studies, using two real call center datasets further demonstrate our findings, and provide operational insights into how one chooses forecasting models for multi‐stream arrivals.

Ye, Han↗

Stochastic analysis of surface metrology

The design and evaluation of the expected performance of optical systems require sophisticated and reliable information about the surface topography for planned optical elements before they are fabricated. Modern x-ray source facilities are reliant upon the availability of optics with unprecedented quality (surface slope accuracy <0.1 μrad). The problem is especially complex in the case of x-ray optics, particularly for the X-ray Surveyor under development and other missions. The high angular resolution and throughput of future x-ray space observatories requires hundreds of square meters of high-quality optics. The uniqueness of the optics and limited number of proficient vendors makes the fabrication extremely time consuming and expensive, mostly due to the limitations in accuracy and measurement rate of metrology used in fabrication. Here, we discuss improvements in metrology efficacy via comprehensive statistical analysis of a compact volume of metrology data. The data are considered stochastic, and a statistical model called invertible time-invariant linear filter (InTILF) is developed now for two-dimensional (2-D) surface profiles to provide compact description of the 2-D data in addition to one-dimensional data treated so far. The InTILF model captures stochastic patterns in the data and can be used as a quality metric and feedback to polishing processes, avoiding high-resolution metrology measurements over the entire optical surface. The modeling, implemented in our BeatMark™ software, allows simulating metrology data for optics made by the same vendor and technology. The data are vital for reliable specification for optical fabrication, to be exactly adequate for the required system performance.

42 ENGINEERING↗

Raman spectroscopy of oxygen carrier particles in harsh environments

Fossil fuel combustion processes generate CO2, a greenhouse gas of considerable concern. To reduce this emission, the National Energy Technology Laboratory (NETL) is evaluating alternative combustion approaches, including chemical looping combustion (CLC). This technique generates relatively pure CO2, suitable for subsequent capture and storage. CLC uses oxygen-carrier particles (OCPs) such as iron oxides, copper oxides, calcium sulfates, etc. to provide oxygen for the combustion process. Optimization of the overall combustion process requires knowledge of the oxidation state (e.g., content of Fe2O3 vs. Fe3O4) of the OCPs during the different stages of the CLC process. Unfortunately, the ability to make on-line measurements of the oxidation state of OCPs in harsh environments is lacking and new sensors need to be developed. We are evaluating non-contact, stand-off Raman spectroscopy to determine the relative concentrations of the oxidized and reduced forms of OCPs at temperatures between 800 °C and 1000 °C, and pressures of about 10 atm. Using cw and pulsed Raman spectroscopy, in combination with a pressurized high-temperature sample chamber, we have optimized the operating parameters such as laser wavelength, laser intensity, collection optic design, focal spot size, etc. and measured Raman spectra of various OCP materials at high temperatures. To extract from the Raman spectra relevant information such as the concentration ratio of a material in different oxidation states, the measured data needs to be processed, and statistical modeling and multivariate calibration need to be performed.

20 FOSSIL-FUELED POWER PLANTS↗

Genome-Wide Analysis of RNA Decay in the Cyanobacterium Synechococcus sp. Strain PCC 7002

ABSTRACT RNA degradation is an important process that influences the ultimate concentration of individual proteins inside cells. While the main enzymes that facilitate this process have been identified, global maps of RNA turnover are available for only a few species. Even in these cases, there are few sequence elements that are known to enhance or destabilize a native transcript; even fewer confer the same effect when added to a heterologous transcript. To address this knowledge gap, we assayed genome-wide RNA degradation in the cyanobacterium Synechococcus sp. strain PCC 7002 by collecting total RNA samples after stopping nascent transcription with rifampin. We quantified the abundance of each position in the transcriptome as a function of time using RNA-sequencing data and later analyzed the global mRNA decay map using machine learning principles. Half-lives, calculated on a per-ORF (open reading frame) basis, were extremely short, with a median half-life of only 0.97 min. Despite extremely rapid turnover of most mRNA, transcripts encoding proteins involved in photosynthesis were both highly expressed and highly stable. Upon inspection of these stable transcripts, we identified an enriched motif in the 3′ untranslated region (UTR) that had similarity to Rho-independent terminators. We built statistical models for half-life prediction and used them to systematically identify sequence motifs in both 5′ and 3′ UTRs that correlate with stabilized transcripts. We found that transcripts linked to a terminator containing a poly(U) tract had a longer half-life than both those without a poly(U) tract and those without a terminator. IMPORTANCE RNA degradation is an important process that affects the final concentration of individual mRNAs, affecting protein expression and cellular physiology. Studies of how RNA is degraded increase our knowledge of this fundamental process as well as enable the creation of genetic tools to manipulate RNA stability. By studying global transcript turnover, we searched for sequence elements that correlated with transcript (in)stability and used these sequences to guide tool design. This study probes global RNA turnover in a cyanobacterium, Synechococcus sp. strain PCC 7002, that both has a unique array of RNases that facilitate RNA degradation and is an industrially relevant strain that could be used to convert CO 2 and sunlight into useful products.

59 BASIC BIOLOGICAL SCIENCES↗

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.↗

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↗

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↗

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↗

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↗

SAVY-4000 O-Ring and Filter Lifetime Extension Report

The 3-year accelerated aging study of the SAVY-4000 O-ring detailed in our previous Lifetime Extension Report showed very little evidence of significant degradation to samples subjected to aggressive elevated temperature and radiation conditions. Whole container thermal aging studies followed by helium leakage testing and compression set measurements were used to establish a failure criterion estimate for O-ring compression set of ≥ 65 %, which led to an O-ring lifetime estimate of greater than 65 years at 80°C. Herein, we describe a statistical model of the 3-year accelerated aging study and present additional results gathered in FY17 that further support the lifetime extension of the O-ring to 40 years at 80°C. Because we are choosing an extremely conservative lifetime value (40 years) compared to the aging study results (>65 years), the O-ring lifetime could be extended further, especially if supported by surveillance results collected over the years to come.

36 MATERIALS SCIENCE↗