Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multivariate analysis”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Real-time elemental analysis of liquids for process monitoring using laser-induced breakdown spectroscopy with a liquid wheel sampling approach

This article presents an engineered sampling system that used a rotating wheel to form a thin liquid layer, permitting the use of laser-induced breakdown spectroscopy (LIBS) for in situ, real-time elemental impurity quantification during liquid processing. The sampling approach was demonstrated on eight elements from across the periodic table (Na, Al, K, Ca, Ti, Sr, Mo, and Yb). Univariate and multivariate calibrations were presented for each element. The average value for percent root mean square errors of cross-validation for the multivariate models was 3.64%, highlighting the method's strong prediction accuracy. Additionally, the limits of detection for each analyte were estimated from their univariate models: Na = 0.0532, Al = 18.5, K = 0.105, Ca = 0.273, Ti = 67.7, Sr = 0.640, Mo = 22.4, and Yb = 22.9 μg mL –1 . Finally, a test in which multivariate models were used to monitor a liquid system for 80 min was performed to investigate the real-time monitoring capabilities of this liquid LIBS sampling approach. Rigorous measurements were performed to effectively predict the absence and concentrations of multiple analytes as they were spiked and diluted. This demonstration showed the feasibility of using LIBS for real-time liquid quantification models with estimated precision ≤ 8.1%. Finally, the limitations of this approach and potential future improvements are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Search for Baryon Number Violation in Top Quark Production and Decay Using Proton-Proton Collisions at $\sqrt{𝑠}$ = 13 TeV

A search is presented for baryon number violating interactions in top quark production and decay. The analysis uses data from proton-proton collisions at a center-of-mass energy of 13 TeV, collected with the CMS detector at the LHC with an integrated luminosity of 138 fb −1 . Candidate events are selected by requiring two oppositely charged leptons (electrons or muons) and exactly one jet identified as originating from a bottom quark. Multivariate discriminants are used to separate the signal from the background. No significant deviation from the standard model prediction is observed. Upper limits are placed on the strength of baryon number violating couplings. For the first time the production of single top quarks via baryon number violating interactions is studied. This allows the search to set the most stringent constraints to date on the branching fraction of the top quark decay to a lepton, an up-type quark (𝑢 or 𝑐), and a down-type quark (𝑑, 𝑠, or 𝑏). The results improve the previous bounds by 3 to 6 orders of magnitude based on the fermion flavor combination of the baryon number violating interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Organic matter distribution in the icy environments of Taylor Valley, Antarctica

Glaciers can accumulate and release organic matter affecting the structure and function of associated terrestrial and aquatic ecosystems. Here we analyzed 18 ice cores collected from six locations in Taylor Valley (McMurdo Dry Valleys), Antarctica to determine the spatial abundance and quality of organic matter, and the spatial distribution of bacterial density and community structure from the terminus of the Taylor Glacier to the coast (McMurdo Sound). Our results showed that dissolved and particulate organic carbon (DOC and POC) concentrations in the ice core samples increased from the Taylor Glacier to McMurdo Sound, a pattern also shown by bacterial cell density. Fluorescence Excitation Emission Matrices Spectroscopy (EEMs) and multivariate parallel factor (PARAFAC) modeling identified one humic-like (C1) and one protein-like (C2) component in ice cores whose fluorescent intensities all increased from the Polar Plateau to the coast. The fluorescence index showed that the bioavailability of dissolved organic matter (DOM) also decreased from the Polar Plateau to the coast. Partial least squares path modeling analysis revealed that bacterial abundance was the main positive biotic factor influencing both the quantity and quality of organic matter. Marine aerosol influenced the spatial distribution of DOC more than katabatic winds in the ice cores. Certain bacterial taxa showed significant correlations with DOC and POC concentrations. Collectively, our results show the tight connectivity among organic matter spatial distribution, bacterial abundance and meteorology in the McMurdo Dry Valley ecosystem.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Anomaly Detection for Online Monitoring of Thermocouple Sensors in the Advanced Test Reactor

This study explores data-driven anomaly detection methods to analyze sensor fail- ures in the Advanced Gas Reactor (AGR) nuclear fuel irradiation experiments. Specifically, we examine failures of thermocouples (TCs), which are critical for mon- itoring and controlling in-reactor temperatures during operation. Failures were pri- marily observed during abrupt power transitions and manifested as sensor drop-outs, drifts, or unexplained behavior. We applied three time-series analysis techniques— rolling mean smoothing, matrix profile, and vector auto-regression (VAR)—to de- tect anomalies in TC data prior to failure events. The rolling mean method effec- tively highlighted deviations aligned with reported failures, while the matrix profile provided partial early warning but sometimes flagged normal fluctuations during power-down periods. VAR shows potential in capturing multivariate dependencies but requires further calibration. A rare case of TC drift was also documented, which did not result in failure, underscoring the challenge of building predictive models with sparse positive examples. Our findings demonstrate that traditional statistical tools can aid anomaly detection but have limited predictive power without richer training data. We propose future directions including synthetic data generation, real- time surrogate modeling, and multi-modal feature integration. This work provides a foundation for applying robust anomaly detection frameworks to mission-critical sensor systems in experimental settings.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Search for Higgsinos in final states with low-momentum lepton-track pairs at 13 TeV

We present a search for the pair production of Higgsinos in final states with large missing transverse momentum and either two reconstructed muons or a reconstructed lepton (muon or electron) and an isolated track. The analyzed data correspond to proton-proton collisions with an integrated luminosity of 137 fb −1 , collected by the CMS experiment at $\sqrt{𝑠}$ =13 TeV in 2016, 2017, and 2018. The signal scenario assumes four nearly mass degenerate Higgsino mass eigenstates: two neutralino states $\tilde{𝜒}^0_2$ and $\tilde{𝜒}^0_1$ with a small mass difference in the range 1–10 GeV and two chargino states $\tilde{𝜒}^±_1$ with an intermediate mass. The analysis focuses on the decay of the heavier neutralino into the lighter one and a virtual 𝑍 boson, which decays into two same-flavor leptons. The leptons have small transverse momentum and/or a small opening angle between the identified muons. An isolated track is used to recover events in which only one of the two leptons is identified. Multivariate discriminants are used to enhance the sensitivity by efficiently rejecting backgrounds from SM processes or misreconstructed tracks and/or leptons. The search explores a unique phase space and probes a previously unexplored region of the signal model parameter space. Mass differences between the two neutralinos are probed down to 1.5 GeV, assuming a Higgsino mass of 100 GeV. The maximum excluded Higgsino mass is 115 GeV.

Hayrapetyan, A. [Yerevan Physics Institute]↗
DOE OSTI2024

Long-Term Statistical Process Monitoring of an Ultrafiltration Water Treatment Process

As water treatment technology has improved, the amount of available process data has substantially increased, making real-time, data-driven fault detection a reality. One shortcoming of the fault detection literature is that methods are usually evaluated by comparing their performance on hand-picked, short-term case studies, which yields no insight into long-term performance. In this work, we first evaluate multiple statistical and machine learning approaches for detrending process data. Then, we evaluate the performance of a PCA-based fault detection approach, applied to the detrended data, to monitor influent water quality, filtrate quality, and membrane fouling of an ultrafiltration membrane system for indirect potable reuse. Based on two short case studies, the adaptive lasso detrending method is selected, and the performance of the multivariate approach is evaluated over more than a year. The method is tested for different sets of three critical tuning parameters, and we find that for long-term, autonomous monitoring to be successful, these parameters should be carefully evaluated. However, in comparison with industry standards of simpler, univariate monitoring or daily pressure decay tests, multivariate monitoring produces substantial benefits in long-term testing.

ammonia↗
DOE OSTI2023

Concordance

Dimension reduction techniques have long been an important topic in statistics, and active subspaces (AS) have received much attention this past decade in the computer experiments literature. The most common approach towards estimating the AS is to use Monte Carlo with numerical gradient evaluation. This software takes a new approach towards estimating the active subspace, in which Multivariate Adaptive Regression Splines (MARS) is used as a surrogate model, and relevant calculation are obtained analytically. Using a MARS surrogate has many advantages including improved scaling, better estimation of active subspaces in high dimensions and the ability to handle a large number of prior distributions in closed form. This software is also capable of conducting a so-called "concordance analysis", using a generalization of ASs which we refer to as "Co-Active" subspaces.

Rumsey, Kelin↗
DOE OSTI2024

Regularized Differentiation for Bioburden Density Estimation in Planetary Protection

In this paper, we propose and investigate the performance of two novel shrinkage estimators for bioburden density estimation in planetary protection. The estimators are based on the regularized differentiation of a cumulative count of colony forming units collected throughout the data collecting session or the life cycle of the entire mission. The regularized differentiation recasts the problem of bioburden density estimation as a linear least squares problem. The least squares problem is then solved through regularization techniques, such as truncated singular value decomposition and penalized least squares. The regularization is necessary to avoid noise amplification during the differentiation of noisy data. The two regularization estimators are compared with four other commonly used estimators to simultaneously evaluate the means of multivariable independent Poisson distributions: the maximum likelihood, noninformative Bayes estimator with Jeffreys prior, Empirical Bayes using conjugate gamma-Poisson model with gamma parameters selected by method of moments, and the Clevenson-Zidek estimator. It is shown through computer-simulated data that the regularized differentiation based on ridge regression has the smallest mean-squared error among all estimators. The analysis of shrinkage mechanism implemented by regularized differentiation is performed, and it is shown that the regularized differentiation amounts to performing a weighted averaging of all the samples. The weights are determined by the regularization parameter automatically selected by the L-curve technique. Since the method of least squares makes no distributional assumptions about the data, it presents an attractive technique for bioburden density estimation when there are concerns about the misspecification of the distributional model. The paper concludes with the analysis of the bioburden data collected during InSight mission and directions for future work.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Modeling interval trendlines: Symbolic singular spectrum analysis for interval time series

Abstract In this article we propose an extension of singular spectrum analysis for interval‐valued time series. The proposed methods can be used to decompose and forecast the dynamics governing a set‐valued stochastic process. The resulting components on which the interval time series is decomposed can be understood as interval trendlines, cycles, or noise. Forecasting can be conducted through a linear recurrent method, and we devised generalizations of the decomposition method for the multivariate setting. The performance of the proposed methods is showcased in a simulation study. We apply the proposed methods so to track the dynamics governing the Argentina Stock Market (MERVAL) in real time, in a case study over a period of turbulence that led to discussions of the government of Argentina with the International Monetary Fund.

de Carvalho, Miguel↗
DOE OSTI2026

Measurement of the top quark pair production cross section in PbPb collisions at $\sqrt{s_{\mathrm{NN}}}=5.36$ TeV

The inclusive cross section for top quark pair ($\mathrm{t}\overline{\mathrm{t}}$) production in lead-lead (PbPb) collisions is reported for the first time at a center-of-mass energy per nucleon pair of 5.36 TeV. The analysis uses data corresponding to an integrated luminosity of 1.58 nb −1 collected by the CMS experiment at the CERN LHC in 2023. The $\mathrm{t}\overline{\mathrm{t}}$ production cross section, ${\sigma}_{\mathrm{t}\overline{\mathrm{t}}}={3.42}_{-0.51}^{+0.54}{\left(\mathrm{stat}\right)}_{-0.43}^{+0.50}\left(\mathrm{syst}\right)$ μb, is measured in dilepton final states using a fit to a multivariate discriminator that combines the decay electron and muon kinematic properties with the multiplicity of bottom quark jets. The result is consistent with perturbative quantum chromodynamics calculations at next-to-next-to-leading order (NNLO) accuracy employing several nuclear parton distribution functions. In addition, the Drell–Yan production cross section (σ DY ) for dilepton masses above 10 GeV and the ratio of $\mathrm{t}\overline{\mathrm{t}}$ to DY cross sections $\left({R}_{\mathrm{t}\overline{\mathrm{t}}/\mathrm{DY}}\right)$ are found to be compatible with the NNLO predictions. The observables ${\sigma}_{\mathrm{t}\overline{\mathrm{t}}}$, σ DY , and ${R}_{\mathrm{t}\overline{\mathrm{t}}/\mathrm{DY}}$ are measured separately for central and semicentral PbPb collisions to investigate for the first time the dependence of top quark production on the collision impact parameter.

Heavy Ion Experiments↗
DOE OSTI2025

Search for the production of a Higgs boson in association with a single top quark in pp collisions at $\sqrt{s}=13$ TeV with the ATLAS detector

A search for the production of a Higgs boson in association with a single top quark, tH, is presented. The analysis uses proton-proton collision data corresponding to an integrated luminosity of 140 fb −1 at a centre-of-mass energy of 13 TeV, collected by the ATLAS detector at the LHC. The search targets Higgs-boson decays into $b\bar{b}$, WW * , ZZ * , and ττ, accompanied by an isolated lepton (electron or muon) from the top-quark decay. Multivariate techniques are employed to enhance the separation between signal and background processes. The observed signal strength, μ tH , defined as the ratio between the measured cross-section and the predicted Standard Model value, is μ tH = 8.1 ± 2.6 (stat.) ± 2.0 (syst.). The significance of the observed (expected) signal above the background-only expectation is 2.8 (0.4) standard deviations. The corresponding observed (expected) upper limit at the 95% confidence level on the tH cross-section is found to be 13.9 (6.1) times the value predicted by the Standard Model. An interpretation with an inverted sign of the top-quark Yukawa coupling is performed, and the signal strength and corresponding limit are reported.

Hadron-Hadron Scattering↗
DOE OSTI2025

Sensitive detection of structural dynamics using a statistical framework for comparative crystallography

Chemical and conformational changes are crucial to protein function and its pharmacological control. X-ray crystallography can reveal these changes in atomic detail, but standard analysis methods, which refine separate datasets, often overlook differences that are subtle or arise in only a subset of molecules. Direct comparison of crystallographic datasets is, in principle, more powerful, but systematic errors (“scales”) often mask changes in the crystallographic observables (“structure factors”). Machine learning algorithms that jointly estimate scales and structure factors can address this limitation. Here, we augment this approach with multivariate, structured priors derived from crystallographic theory, implemented in the variational deep learning framework Careless. Doing so strongly improves the detection of protein dynamics, element-specific anomalous signals, and the binding of drug candidates, offering a robust approach to comparative crystallography and, potentially, to detection of protein dynamics by other structure determination methods.

Hekstra, Doeke R. [Harvard Univ., Cambridge, MA (U↗
DOE OSTI2023

Contraceptive Sabotage and Contraceptive Use at the Time of Pregnancy: An Analysis of People with a Recent Live Birth in the United States

Contraceptive sabotage and other forms of intimate partner violence (IPV) can interfere with contraceptive use. We used 2012 to 2015 Pregnancy Risk Assessment Monitoring System data from 8,981 people residing in five states who reported that when they became pregnant, they were not trying to get pregnant. We assessed the relationships between ever experiencing contraceptive sabotage and physical IPV 12 months before pregnancy (both by the current partner) and contraceptive use at the time of pregnancy using multivariable logistic regression. We also assessed the joint associations between physical IPV 12 months before pregnancy and ever experienced contraceptive sabotage with contraceptive use at the time of pregnancy. Few people ever experienced contraceptive sabotage (1.8%; 95% confidence interval [CI]: 1.4, 2.3) or physical IPV 12 months before pregnancy (2.8%; 95% CI: 2.3, 3.3). In models adjusted for age, race/ethnicity, marital status, education, and state of residence, ever experiencing contraceptive sabotage was associated with contraceptive use at the time of pregnancy (adjusted odds ratio [aOR]: 1.73; 95% CI: 1.06, 2.82), but not with physical IPV 12 months before pregnancy (aOR: 0.69; 95% CI: 0.46, 1.02). When examining the joint association, compared to not ever experiencing contraceptive sabotage or physical IPV 12 months before pregnancy, ever experiencing contraceptive sabotage was significantly related to contraceptive use at the time of pregnancy (aOR: 1.72; 95% CI: 1.00, 2.95). However, it was not associated with experiencing physical IPV 12 months before pregnancy (aOR: 0.68; 95% CI: 0.45, 1.04) or with experiencing both contraceptive sabotage and physical IPV 12 months before pregnancy (aOR: 1.21; 95% CI: 0.42, 3.50), compared to not ever experiencing contraceptive sabotage or physical IPV 12 months before pregnancy. Our study highlights that current partner contraceptive sabotage may motivate those not trying to get pregnant to use contraception; however, all people in our sample still experienced a pregnancy.

Huber-Krum, Sarah↗
DOE OSTI2022

Predictability and empirical dynamics of fisheries time series in the North Pacific

Previous studies have documented a strong relationship between marine ecosystems and large-scale modes of sea surface height (SSH) and sea surface temperature (SST) variability in the North Pacific such as the Pacific Decadal Oscillation and the North Pacific Gyre Oscillation. In the central and western North Pacific along the Kuroshio-Oyashio Extension (KOE), the expression of these modes in SSH and SST is linked to the propagation of long oceanic Rossby waves, which extend the predictability of the climate system to ~3 years. Using a multivariate physical-biological linear inverse model (LIM) we explore the extent to which this physical predictability leads to multi-year prediction of dominant fishery indicators inferred from three datasets (i.e., estimated biomasses, landings, and catches). We find that despite the strong autocorrelation in the fish indicators, the LIM adds dynamical forecast skill beyond persistence up to 5-6 years. By performing a sensitivity analysis of the LIM forecast model, we find that two main factors are essential for extending the dynamical predictability of the fishery indicators beyond persistence. The first is the interaction of the fishery indicators with the SST/SSH of the North and tropical Pacific. The second is the empirical relationship among the fisheries time series. This latter component reflects stock-stock interactions as well as common technological and human socioeconomic factors that may influence multiple fisheries and are captured in the training of the LIM. These results suggest that empirical dynamical models and machine learning algorithms, such as the LIM, provide an alternative and promising approach for forecasting key ecological indicators beyond the skill of persistence.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Real-time confinement regime detection in fusion plasmas with convolutional neural networks and high-bandwidth edge fluctuation measurements

Abstract A real-time detection of the plasma confinement regime can enable new advanced plasma control capabilities for both the access to and sustainment of enhanced confinement regimes in fusion devices. For example, a real-time indication of the confinement regime can facilitate transition to the high-performing wide-pedestal (WP) quiescent H-mode, or avoid unwanted transitions to lower confinement regimes that may induce plasma termination. To demonstrate real-time confinement regime detection, we use the 2D beam emission spectroscopy (BES) diagnostic system to capture localized density fluctuations of long wavelength turbulent modes in the edge region at a 1 MHz sampling rate. BES data from 330 discharges in either L-mode, H-mode, quiescent H (QH)-mode, or WP QH-mode were collected from the DIII-D tokamak and curated to develop a high-quality database to train a deep-learning classification model for real-time confinement detection. We utilize the 6×8 spatial configuration with a time window of 1024 µ s and recast the input to obtain spectral-like features via fast Fourier transform preprocessing. We employ a shallow 3D convolutional neural network for the multivariate time-series classification task and utilize a softmax in the final dense layer to retrieve a probability distribution over the different confinement regimes. Our model classifies the global confinement state on 44 unseen test discharges with an average F 1 score of 0.94, using only ∼1 ms snippets of BES data at a time. This activity demonstrates the feasibility for real-time data analysis of fluctuation diagnostics in future devices such as ITER, where the need for reliable and advanced plasma control is urgent.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An Exploratory Approach Using Regression and Machine Learning in the Analysis of Mass Absorption Cross Section of Black Carbon Aerosols: Model Development and Evaluation

Mass absorption cross-section of black carbon (MAC BC ) describes the absorptive cross-section per unit mass of black carbon, and is, thus, an essential parameter to estimate the radiative forcing of black carbon. Many studies have sought to estimate MAC BC from a theoretical perspective, but these studies require the knowledge of a set of aerosol properties, which are difficult and/or labor-intensive to measure. We therefore investigate the ability of seven data analytical approaches (including different multivariate regressions, support vector machine, and neural networks) in predicting MAC BC for both ambient and biomass burning measurements. Our model utilizes multi-wavelength light absorption and scattering as well as the aerosol size distributions as input variables to predict MAC BC across different wavelengths. We assessed the applicability of the proposed approaches in estimating MAC BC using different statistical metrics (such as coefficient of determination (R 2 ), mean square error (MSE), fractional error, and fractional bias). Overall, the approaches used in this study can estimate MAC BC appropriately, but the prediction performance varies across approaches and atmospheric environments. Based on an uncertainty evaluation of our models and the empirical and theoretical approaches to predict MAC BC , we preliminarily put forth support vector machine (SVM) as a recommended data analytical technique for use. We provide an operational tool built with the approaches presented in this paper to facilitate this procedure for future users.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Downscaled CMIP5 projections of physical fire risk understate historical trends

Reliable projections of wildfire risk are important for multi-sector impacts analysis. Statistically downscaled and bias-corrected Earth system model ensemble products are routinely used to analyze regional physical wildfire risk, but evaluations of historical observed trends and variability are lacking. Here, we evaluate physical fire risk over the western United States using the Canadian Forest Fire Weather Index (FWI) by comparing model outputs from the Coupled Model Intercomparison Project Phase 5 (CMIP5), statistically downscaled via the Multivariate Adaptive Constructed Analogs (MACA) approach, against the observational target dataset gridMET, a gridded high-resolution surface meteorological product. We analyze multidecadal trends and interannual variability in seasonal average FWI for the historical period and future projections under two emissions scenarios, and we compare MACA-CMIP5 ensemble results with a simple time series model that generates historical and future projections of seasonal FWI based on bootstrapping observed historical trends and variability. Our findings indicate that MACA-CMIP5 accurately captures the magnitude and spatial patterns of seasonally averaged FWI but tends to underestimate historical decadal trends. We show that future increases in fire risk may be underestimated relative to the simple time series model that projects historical variability into the future. We also highlight that model biases in relative humidity contribute significantly to model-data differences. Our results underscore the importance of historical hindcasting exercises for informing broader multi-sector applications.

FWI↗