Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical model”

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,369 records · Page 76

Predictive understanding of the surface tension and velocity of sound in ionic liquids using machine learning

Knowledge of the physical properties of ionic liquids (ILs), such as the surface tension and speed of sound, is important for both industrial and research applications. Unfortunately, technical challenges and costs limit exhaustive experimental screening efforts of ILs for these critical properties. Previous work has demonstrated that the use of quantum-mechanics-based thermochemical property prediction tools, such as the conductor-like screening model for real solvents, when combined with machine learning (ML) approaches, may provide an alternative pathway to guide the rapid screening and design of ILs for desired physiochemical properties. However, the question of which machine-learning approaches are most appropriate remains. In the present study, we examine how different ML architectures, ranging from tree-based approaches to feed-forward artificial neural networks, perform in generating nonlinear multivariate quantitative structure–property relationship models for the prediction of the temperature- and pressure-dependent surface tension of and speed of sound in ILs over a wide range of surface tensions (16.9–76.2 mN/m) and speeds of sound (1009.7–1992 m/s). The ML models are further interrogated using the powerful interpretation method, shapley additive explanations. We find that several different ML models provide high accuracy, according to traditional statistical metrics. The decision tree-based approaches appear to be the most accurate and precise, with extreme gradient-boosting trees and gradient-boosting trees being the best performers. However, our results also indicate that the promise of using machine-learning to gain deep insights into the underlying physics driving structure–property relationships in ILs may still be somewhat premature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Bayesian analysis of nuclear deformation properties with Skyrme energy functionals

In spite of numerous scientific and practical applications, there is still no comprehensive theoretical description of the nuclear fission process based solely on protons, neutrons and their interactions. Furthermore, the most advanced simulations of fission are currently carried out within nuclear density functional theory (DFT). In spite of being fully quantum-mechanical and rooted in the theory of nuclear forces, DFT still depends on a dozen or so parameters characterizing the energy functional. Calibrating these parameters on experimental data results in uncertainties that must be quantified for applications. This task is very challenging because of the high computational cost of DFT calculations for fission. In this paper, we use Gaussian processes to build emulators of DFT models in order to quantify and propagate statistical uncertainties of theoretical predictions for a range of nuclear deformations relevant to describing the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Micro on a macroscale: relating microbial-scale soil processes to global ecosystem function

ABSTRACT Soil microorganisms play a key role in driving major biogeochemical cycles and in global responses to climate change. However, understanding and predicting the behavior and function of these microorganisms remains a grand challenge for soil ecology due in part to the microscale complexity of soils. It is becoming increasingly clear that understanding the microbial perspective is vital to accurately predicting global processes. Here, we discuss the microbial perspective including the microbial habitat as it relates to measurement and modeling of ecosystem processes. We argue that clearly defining and quantifying the size, distribution and sphere of influence of microhabitats is crucial to managing microbial activity at the ecosystem scale. This can be achieved using controlled and hierarchical sampling designs. Model microbial systems can provide key data needed to integrate microhabitats into ecosystem models, while adapting soil sampling schemes and statistical methods can allow us to collect microbially-focused data. Quantifying soil processes, like biogeochemical cycles, from a microbial perspective will allow us to more accurately predict soil functions and address long-standing unknowns in soil ecology.

59 BASIC BIOLOGICAL SCIENCES↗

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: BAO and RSD measurements from the anisotropic power spectrum of the quasar sample between redshift 0.8 and 2.2

ABSTRACT We measure the clustering of quasars of the final data release (DR16) of eBOSS. The sample contains $343\, 708$ quasars between redshifts 0.8 ≤ z ≤ 2.2 over $4699\, \mathrm{deg}^2$. We calculate the Legendre multipoles (0,2,4) of the anisotropic power spectrum and perform a BAO and a Full-Shape (FS) analysis at the effective redshift zeff = 1.480. The errors include systematic errors that amount to 1/3 of the statistical error. The systematic errors comprise a modelling part studied using a blind N-body mock challenge and observational effects studied with approximate mocks to account for various types of redshift smearing and fibre collisions. For the BAO analysis, we measure the transverse comoving distance DM(zeff)/rdrag = 30.60 ± 0.90 and the Hubble distance DH(zeff)/rdrag = 13.34 ± 0.60. This agrees with the configuration space analysis, and the consensus yields: DM(zeff)/rdrag = 30.69 ± 0.80 and DH(zeff)/rdrag = 13.26 ± 0.55. In the FS analysis, we fit the power spectrum using a model based on Regularised Perturbation Theory, which includes redshift space distortions and the Alcock–Paczynski effect. The results are DM(zeff)/rdrag = 30.68 ± 0.90 and DH(zeff)/rdrag = 13.52 ± 0.51 and we constrain the linear growth rate of structure f(zeff)σ8(zeff) = 0.476 ± 0.047. Our results agree with the configuration space analysis. The consensus analysis of the eBOSS quasar sample yields: DM(zeff)/rdrag = 30.21 ± 0.79, DH(zeff)/rdrag = 3.23 ± 0.47, and f(zeff)σ8(zeff) = 0.462 ± 0.045 and is consistent with a flat ΛCDM cosmological model using Planck results.

Neveux, Richard↗

Atacama Cosmology Telescope: DR6 gravitational lensing and SDSS BOSS cross-correlation measurement and constraints on gravity with the 𝐸 𝐺 statistic

We derive new constraints on the 𝐸 𝐺 statistic as a test of gravity, combining the cosmic microwave background (CMB) lensing map estimated from Data Release 6 (DR6) of the Atacama Cosmology Telescope with Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (SDSS BOSS) CMASS and LOWZ galaxy data. We develop an analysis pipeline to measure the cross-correlation between CMB lensing maps and galaxy data, following a blinding policy and testing the approach through null and consistency checks. By testing the equivalence of the spatial and temporal gravitational potentials, the 𝐸 𝐺 statistic can distinguish Λ⁢ CDM from alternative models of gravity. We find 𝐸 𝐺 ⁡(𝑧 eff = 0.555) = 0.3⁢1$^{+0.06}_{−0.05}$ for Atacama Cosmology Telescope (ACT) and CMASS data at 68.28% confidence level, and 𝐸 𝐺 ⁡(𝑧 eff = 0.316) = 0.4⁢9$^{+0.14}_{−0.11}$ for the ACT and LOWZ. Systematic errors are estimated to be 3% and 4%, respectively. Including CMB lensing information from Planck PR4 results in 𝐸 𝐺 ⁡(𝑧 eff = 0.555) = 0.3⁢4$^{+0.05}_{−0.05}$ with CMASS and 𝐸 𝐺 ⁡(𝑧 eff = 0.316) = 0.4⁢3$^{+0.11}_{−0.09}$ with LOWZ. These are consistent with predictions for the Λ⁢ CDM model that best fits the Planck CMB anisotropy and SDSS BOSS baryon acoustic oscillations (BAO), where 𝐸$^{GR}_{𝐺⁡}$(𝑧 eff =0.555) =0.401 ± 0.005 for CMB lensing combined with CMASS and 𝐸$^{GR}_{𝐺}$⁡(𝑧 eff = 0.316) = 0.452 ± 0.005 combined with LOWZ. We also find 𝐸 𝐺 to be scale independent, with probability to exceed >5%, as predicted by general relativity. The methods developed in this work are also applicable to improved future analyses with upcoming spectroscopic galaxy samples and CMB lensing measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

Data assimilation empowered neural network parametrizations for subgrid processes in geophysical flows

In the past couple of years, there has been a proliferation in the use of machine learning approaches to represent subgrid-scale processes in geophysical flows with an aim to improve the forecasting capability and to accelerate numerical simulations of these flows. Despite its success for different types of flow, the online deployment of a data-driven closure model can cause instabilities and biases in modeling the overall effect of subgrid-scale processes, which in turn leads to inaccurate prediction. To tackle this issue, we exploit the data assimilation technique to correct the physics-based model coupled with the neural network as a surrogate for unresolved flow dynamics in multiscale systems. In particular, we use a set of neural network architectures to learn the correlation between resolved flow variables and the parametrizations of unresolved flow dynamics and formulate a data assimilation approach to correct the hybrid model during their online deployment. We illustrate our framework in a set of applications of the multiscale Lorenz 96 system for which the parametrization model for unresolved scales is exactly known, and the two-dimensional Kraichnan turbulence system for which the parametrization model for unresolved scales is not known a priori. Our analysis, therefore, comprises a predictive dynamical core empowered by (i) a data-driven closure model for subgrid-scale processes, (ii) a data assimilation approach for forecast error correction, and (iii) both data-driven closure and data assimilation procedures. We show significant improvement in the long-term prediction of the underlying chaotic dynamics with our framework compared to using only neural network parametrizations for future prediction. Moreover, we demonstrate that these data-driven parametrization models can handle the non-Gaussian statistics of subgrid-scale processes, and effectively improve the accuracy of outer data assimilation workflow loops in a modular nonintrusive way.

42 ENGINEERING↗

Continental Scale Hydrostratigraphy: Comparing Geologically Informed Data Products to Analytical Solutions

Abstract This study synthesizes two different methods for estimating hydraulic conductivity (K) at large scales. We derive analytical approaches that estimate K and apply them to the contiguous United States. We then compare these analytical approaches to three‐dimensional, national gridded K data products and three transmissivity (T) data products developed from publicly available sources. We evaluate these data products using multiple approaches: comparing their statistics qualitatively and quantitatively and with hydrologic model simulations. Some of these datasets were used as inputs for an integrated hydrologic model of the Upper Colorado River Basin and the comparison of the results with observations was used to further evaluate the K data products. Simulated average daily streamflow was compared to daily flow data from 10 USGS stream gages in the domain, and annually averaged simulated groundwater depths are compared to observations from nearly 2000 monitoring wells. We find streamflow predictions from analytically informed simulations to be similar in relative bias and Spearman's rho to the geologically informed simulations. R ‐squared values for groundwater depth predictions are close between the best performing analytically and geologically informed simulations at 0.68 and 0.70 respectively, with RMSE values under 10 m. We also show that the analytical approach derived by this study produces estimates of K that are similar in spatial distribution, standard deviation, mean value, and modeling performance to geologically‐informed estimates. The results of this work are used to inform a follow‐on study that tests additional data‐driven approaches in multiple basins within the contiguous United States.

54 ENVIRONMENTAL SCIENCES↗

Machine Learning for Improved Availability of the SNS Klystron High Voltage Converter Modulators

Beam availability has increased at the SNS, however, the targeted availability is greater than 95 %, while the SNS has failed to meet lower targets in the past. The HVCM used to power the linac klystrons have been one source of lost beam time and was chosen to explore using AI/ML techniques to improve reliability. Among the possibilities being explored are automating the tuning of HVCMs and predicting component failures such as capacitor aging, rectifier assemblies containing hundreds of diodes, and insulating oil degradation. The methodology pursued includes data cleaning, de-noising, post-analysis data labeling, and machine learning model development. We explore using Long Short-Term Memory and autoencoders for anomaly detection and prognostication used to schedule maintenance. We evaluate the use of model regularizers and constraints to improve the performance of the model and investigate methods to estimate the uncertainty of the models to provide a robust prediction with statistical interoperability. This paper describes the operational experience and known failures of the HVCMs and the proposed ML methodology and the preliminary results of training the AI/ML algorithms.

Pappas, G. C.↗

EMP-Resilient Electric Grid Transformer Analysis

The response of a high - voltage (HV) transformer to fast rise time voltages, such as that from a electromagnetic pulse (EMP) can result in interruption of power distribution and possibly system failure. To help identify these potential occurrences, it is necessary to develop a transformer model that not only captures the input/output response of the transformer but also the internal behavior. The model constructed should cover the frequency band of interest while capturing the internal physical and electrical characteristics. This broad-band, high-fidelity model would enable the prediction of unwanted effects through simulation. A proposed modeling scheme for a HV transformer is described in Part 1 of this report. Part 2 of this report details assessments of internal voltage and electrical field holdoff testing of transformer insulation dielectric breakdown, including comparison of low frequency (DC/60 Hz) holdoff to rise times characteristic of lightning (1 s) and EMP E1 transients (10 - 30 ns). This initial project is a path toward establishing electrical grid transformer failure criteria for EMP voltage transients. We developed modeling methods and measured breakdown electrical field statistical distributions for direct current, 60 Hz, lightning and EMP characteristic voltage rise times. Methods of nanosecond-scale capacitive discharge unit high voltage source development, suggestions for derating of 60 Hz insulation maximum electrical fields for EMP nanosecond pulse voltage withstand rating, and potential methods for increasing transformer resilience to such fast rise time pulses are described.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Multi-Differential Semi-Inclusive Charged-Current Muon Neutrino Cross Sections on Helium-4 in MINERvA and the Helium-to-CH Target Cross Section Ratio

A comprehensive international effort has been underway to elucidate the properties and behaviors of neutrinos. A major source of systematic uncertainties in studying neutrino-induced interactions comes from neutrino-nucleus cross-section models, highlighting the need for more precise statistical measurements. MINERvA, an on-axis neutrino-nucleus scattering experiment located at the Fermi National Accelerator Laboratory, was established to produce neutrino cross-section measurements with many different nuclei. The helium target is the lightest nucleus to be measured by MINERvA. Using the NuMI medium energy muon neutrino data set, we present preliminary results and summarize the extraction of the charged current (CC) muon neutrino - helium-4 semi-inclusive ($\nu_{\mu}$+$^{4}\textrm{He}$$\rightarrow$$\mu^{-}$+N$\pi$+M$\textit{p}$) multi-differential cross-section extraction as a function of transverse ($P_{T}$) and longitudinal ($P_{L}$) muon momentum with respect to the neutrino beamline. We define the final state topology as CC-N$\pi$M$\textit{p}$, with at least two reconstructed tracks: a muon and a combination of N protons and M pions, where $N + M > 0$. To probe the dependence on the size of the nucleus in neutrino-induced interaction phenomena we present a differential cross-section ratio of helium-4 to MINERvA's hydrocarbon target (\textrm{CH}) as a function of transverse muon momentum.

Nguyen, Christian↗

On the uncertainty of long-period return values of extreme daily precipitation

Methods for calculating return values of extreme precipitation and their uncertainty are compared using daily precipitation rates over the Western U.S. and Southwestern Canada from a large ensemble of climate model simulations. The roles of return-value estimation procedures and sample size in uncertainty are evaluated for various return periods. We compare two different generalized extreme value (GEV) parameter estimation techniques, namely L-moments and maximum likelihood (MLE), as well as empirical techniques. Even for very large datasets, confidence intervals calculated using GEV techniques are narrower than those calculated using empirical methods. Furthermore, the more efficient L-moments parameter estimation techniques result in narrower confidence intervals than MLE parameter estimation techniques at small sample sizes, but similar best estimates. It should be noted that we do not claim that either parameter fitting technique is better calibrated than the other to estimate long period return values. While a non-stationary MLE methodology is readily available to estimate GEV parameters, it is not for the L-moments method. Comparison of uncertainty quantification methods are found to yield significantly different estimates for small sample sizes but converge to similar results as sample size increases. Finally, practical recommendations about the length and size of climate model ensemble simulations and the choice of statistical methods to robustly estimate long period return values of extreme daily precipitation statistics and quantify their uncertainty.

54 ENVIRONMENTAL SCIENCES↗

SDSS J1058+5443: A Blue Quasar without Optical/NUV Broad Emission Lines

In this paper, the blue quasar SDSS J105816.19+544310.2 (=SDSS J1058+5443) at redshift 0.479 has been reported as the best true type 2 quasar candidate with the disappearance of central broad-line regions. There are no definite conclusions on the very existence of true type 2 active galactic nuclei (AGN), mainly due to detected optical broad emission lines in high-quality spectra of some previously classified true type 2 AGN candidates. Here, unlike previously reported true type 2 AGN candidates among narrow emission-line galaxies with weak AGN activities but strong stellar lights, the definitely blue quasar SDSS J1058+5443 can be well confirmed as a true type 2 quasar due to apparent quasar-shape blue continuum emissions but an apparent loss of both the optical broad Balmer emission lines and the near-UV (NUV) broad Mg II emission line. Based on different model functions and the F-test statistical technique, after considering blueshifted optical and UV Fe II emissions, there are no apparent broad optical Balmer emission lines and/or broad NUV Mg II lines, and the confidence level is smaller than 1σ in support of broad optical and NUV emission lines. Moreover, assuming the virialization assumption to broad-line emission clouds, the reconstructed broad emission lines strongly indicate that the probable intrinsic broad emission lines, if they exist, cannot be hidden or overwhelmed in the noise of the Sloan Digital Sky Survey spectrum of SDSS J1058+5443. Therefore, SDSS J1058+5443 is so far the best and most robust true type 2 quasar candidate, leading to the clear conclusion of the very existence of true type 2 AGN.

79 ASTRONOMY AND ASTROPHYSICS↗

Searching for Spectroscopic Signatures of Ongoing Quenching in SDSS Galaxies

In this paper, we estimate the “star formation change parameter,” SFR 79 , which characterizes the current SFR relative to the average during the last 800 Myr, for ~300,000 galaxies selected from the Sloan Digital Sky Survey (SDSS). The goals are to examine, in a much larger and independent sample, the trends previously reported in a sample of star-forming (SF) MaNGA galaxies and also to search for spectroscopic signatures of ongoing quenching in the so-called “Green Valley,” which is generally believed to contain galaxies that are migrating from the SF population to the quenched population of galaxies. Measuring SFR 79 for our large sample of SDSS galaxies, we first confirm the basic results of SF galaxies published by Wang & Lilly. We then discuss in detail the calibration and meaning of SFR 79 for galaxies that are well below the SFMS and establish the expected statistical signature of systematic ongoing quenching from modeling the z ~ 0 quenching rate of the SF population. We conclude that it is not possible at present to establish unambiguous observational evidence for systematic ongoing quenching processes, due to limitations both in the noise of the observational data, in particular in the measurements of H δ absorption, and in the calibration of SFR 79 , as well as biases introduced by the necessity of selecting objects with significant H α emission. We do however see plausible indications of ongoing quenching, which are quantitatively consistent with expectations from “growth+quenching” models of galaxy evolution and a typical e-folding timescale for quenching of ~500 Myr.

79 ASTRONOMY AND ASTROPHYSICS↗

Planetary Spin and Obliquity from Mergers

In planetary systems with sufficiently small inter-planet spacing, close encounters can lead to planetary collisions/mergers or ejections. Here, we study the spin property of the merger products of two giant planets in a statistical manner using numerical simulations and analytical modeling. Planetary collisions lead to rapidly rotating objects and a broad range of obliquities. We find that, under typical conditions for two-planet scatterings, the distributions of spin magnitude S and obliquity ${\theta }_{\mathrm{SL}}$ of the merger products have simple analytical forms: f S ∝ S and ${f}_{\cos {\theta }_{\mathrm{SL}}}\propto {(1-{\cos }^{2}{\theta }_{\mathrm{SL}})}^{-1/2}$. Through parameter studies, we determine the regime of validity for the analytical distributions of spin and obliquity. Since planetary mergers are a major outcome of planet–planet scatterings, observational search for the spin/obliquity signatures of exoplanets would provide important constraints on the dynamical history of planetary systems.

79 ASTRONOMY AND ASTROPHYSICS↗

First Muon Neutrino on Argon Cross Section Measurements With an Off-Axis Beam

While over ninety years have passed since neutrinos were first predicted in 1930 by Wolfgang Pauli, much about them is still unknown [20]. How many types of neutrinos there are, how precisely neutrinos change flavor, and how neutrinos interact with various types of matter are largely unanswered questions that have implications for fields beyond particle physics, such as cosmology and astrophysics. The final question can be answered partially by differential cross section measurements, two of which constitute the primary result of this thesis. The NeUtrinos at the Main Injector (NuMI) Beam at Fermilab in Batavia, IL provides a source of muon neutrinos and muon antineutrinos to the MicroBooNE experiment, an 85 ton liquid argon time projection chamber (LArTPC). The technology allows for high precision particle tracking and energy reconstruction that can be used to make cross section measurements. Necessary work for making these measurements for which the author of this thesis made an essential contribution is presented in detail in References [21] [22] [11] [23]. The measurement presented in this thesis is made with 2.187e20 Protons on Target (POT), the standard metric for an amount of beam exposure, of data when the beam was in forward horn current mode. Single-differential cross section measurements in terms of muon candidate kinetic energy and direction are made, and their statistical fits with simulation are 5.8 (12 d.o.f.) and 11.4 (10 d.o.f.), respectively. This indicates that the muon kinematics for muon neutrino and muon antineutrino charged-current (CC) events primarily in the energy range [0-3] GeV can be trusted for oscillation experiments, including DUNE, because of the agreement to within the associated systematic and statistical uncertainties between the measurement and the model in every bin of both single-differential cross sections. This result adds to the corresponding result for the Booster Neutrino Beam (BNB) because it covers a completely different energy region than that measurement, a double-differential cross section in terms of muon momentum and direction [14]. This result contributes to reducing the uncertainty in near-far event rate comparisons in long baseline neutrino experiments and also adds to existing knowledge of the neutrino-nucleus interaction by describing the outgoing muon kinematics in interactions of muon neutrinos and muon antineutrinos. This result is the first-ever measurement of differential cross sections of muon neutrinos and muon antineutrinos in liquid argon made with an off-axis beam.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Systematic identification of cis -regulatory variants that cause gene expression differences in a yeast cross

Sequence variation in regulatory DNA alters gene expression and shapes genetically complex traits. However, the identification of individual, causal regulatory variants is challenging. Here, we used a massively parallel reporter assay to measure the cis-regulatory consequences of 5832 natural DNA variants in the promoters of 2503 genes in the yeast Saccharomyces cerevisiae. We identified 451 causal variants, which underlie genetic loci known to affect gene expression. Several promoters harbored multiple causal variants. In five promoters, pairs of variants showed non-additive, epistatic interactions. Causal variants were enriched at conserved nucleotides, tended to have low derived allele frequency, and were depleted from promoters of essential genes, which is consistent with the action of negative selection. Causal variants were also enriched for alterations in transcription factor binding sites. Models integrating these features provided modest, but statistically significant, ability to predict causal variants. This work revealed a complex molecular basis for cis-acting regulatory variation.

60 APPLIED LIFE SCIENCES↗

Machine Learning for Improved Availability of the SNS Klystron High Voltage Converter Modulators

Beam availability has increased at the SNS, however, the targeted availability is greater than 95 %, while the SNS has failed to meet lower targets in the past. The HVCM used to power the linac klystrons have been one source of lost beam time and was chosen to explore using AI/ML techniques to improve reliability. Among the possibilities being explored are automating the tuning of HVCMs and predicting component failures such as capacitor aging, rectifier assemblies containing hundreds of diodes, and insulating oil degradation. The methodology pursued includes data cleaning, de-noising, post-analysis data labeling, and machine learning model development. We explore using Long Short-Term Memory and autoencoders for anomaly detection and prognostication used to schedule maintenance. We evaluate the use of model regularizers and constraints to improve the performance of the model and investigate methods to estimate the uncertainty of the models to provide a robust prediction with statistical interoperability. This paper describes the operational experience and known failures of the HVCMs and the proposed ML methodology and the preliminary results of training the AI/ML algorithms.

Pappas, Chris↗