Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multivariate data”

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 631 records · Page 35

Measurement of inclusive and differential cross sections for single top quark production in association with a W boson in proton-proton collisions at $\sqrt{s}$ = 13 TeV

Measurements of the inclusive and normalised differential cross sections are presented for the production of single top quarks in association with a W boson in proton-proton collisions at a centre-of-mass energy of 13 TeV. The data used were recorded with the CMS detector at the LHC during 2016–2018, and correspond to an integrated luminosity of 138 fb -1 . Events containing one electron and one muon in the final state are analysed. For the inclusive measurement, a multivariate discriminant, exploiting the kinematic properties of the events is used to separate the signal from the dominant t$\overline{t}$ background. A cross section of 79.2 ± 0.9(stat)$^{+7.7}_{-8.0}$ (syst) ± 1.2 (lumi) pb is obtained, consistent with the predictions of the standard model. For the differential measurements, a fiducial region is defined according to the detector acceptance, and the requirement of exactly one jet coming from the fragmentation of a bottom quark. The resulting distributions are unfolded to particle level and agree with the predictions at next-to-leading order in perturbative quantum chromodynamics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

CsI(Tl) pulse shape discrimination with the Belle II electromagnetic calorimeter as a novel method to improve particle identification at electron–positron colliders

Here we describe the implementation and performance of CsI(Tl) pulse shape discrimination for the Belle II electromagnetic calorimeter, representing the first application of CsI(Tl) pulse shape discrimination for particle identification at an electron–positron collider. The pulse shape characterization algorithms applied by the Belle II calorimeter are described. Control samples of $γ, μ^+, π^±, K^±$ and $ρ/ \bar{ρ}$ are used to demonstrate the significant insight into the secondary particle composition of calorimeter clusters that is provided by CsI(Tl) pulse shape discrimination. Comparisons with simulation are presented and provide further validation for newly developed CsI(Tl) scintillation response simulation techniques, which when incorporated with GEANT4 simulations allow the particle dependent scintillation response of CsI(Tl) to be modelled. Comparisons between data and simulation also demonstrate that pulse shape discrimination can be a new tool to identify sources of improvement in the simulation of hadronic interactions in materials. The $K^0_L$ efficiency and photon-as-hadron fake-rate of a multivariate classifier that is trained to use pulse shape discrimination is presented and comparisons are made to a shower-shape based approach. CsI(Tl) pulse shape discrimination is shown to reduce the photon-as-hadron fake-rate by over a factor of 3 at photon energies of 0.2 GeV and over a factor 10 at photon energies of 1 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Learning the Temporal Effect in Infrared Thermal Videos With Long Short-Term Memory for Quality Prediction in Resistance Spot Welding

With the advances of sensing technology, in-situ infrared thermal videos can be collected from Resistance Spot Welding (RSW) processes. Each video records the formulation process of a weld nugget. The nugget evolution creates a “temporal effect” across the frames, which can be leveraged for real-time, nondestructive evaluation (NDE) of the weld quality. Currently, quality prediction with imaging data mainly focuses on optical feature extraction with Convolutional Neural Network (CNN) but does not make the most of such temporal effect. In this study, pixels corresponding to critical locations on the weld nugget surface are extracted from a video to form multivariate time series (MTS). Multivariate Adaptive Regression Splines (MARS) is used in MTS processing to remove noisy signals related to uninformative frames. A Stacked Long Short-Term Memory (LSTM) model is developed to learn from the processed MTS and then predicts weld nugget size and thickness in real-time NDE. Results from a case study on RSW of Boron steel demonstrates the improvement in prediction accuracy and computational time with the proposed method, as compared to CNN-based weld quality prediction.

Guo, Shenghan↗

Implementation of ILLIAC 4 algorithms for multispectral image interpretation

Research has focused on the design and partial implementation of a comprehensive ILLIAC software system for computer-assisted interpretation of multispectral earth resources data such as that now collected by the Earth Resources Technology Satellite. Research suggests generally that the ILLIAC 4 should be as much as two orders of magnitude more cost effective than serial processing computers for digital interpretation of ERTS imagery via multivariate statistical classification techniques. The potential of the ARPA Network as a mechanism for interfacing geographically-dispersed users to an ILLIAC 4 image processing facility is discussed.

Ray, R. M.↗

High-resolution land/ice imaging using Seasat scatterometer measurements

A new method for obtaining high-resolution images (to 4 km) of the land backscatter from low-resolution Seasat-A scatterometer (SASS) measurements is introduced. The method utilizes the measurement cell overlap in multiple spacecraft passes over the region of interest and signal processing techniques to generate high-resolution images of the radar backscatter. The overlap is exploited to estimate the underlying high-resolution surface radar backscatter characteristics using a robust multivariate image reconstruction algorithm. The algorithm has been designed to operate in the high-noise environment typical of scatterometer measurements. The ultimate resolution obtainable is a function of the number of measurements and the measurement overlap. Sample results based on SASS data are provided.

Long, D. G.↗

Design and Implementation of a Parallel Multivariate Ensemble Kalman Filter for the Poseidon Ocean General Circulation Model

A multivariate ensemble Kalman filter (MvEnKF) implemented on a massively parallel computer architecture has been implemented for the Poseidon ocean circulation model and tested with a Pacific Basin model configuration. There are about two million prognostic state-vector variables. Parallelism for the data assimilation step is achieved by regionalization of the background-error covariances that are calculated from the phase-space distribution of the ensemble. Each processing element (PE) collects elements of a matrix measurement functional from nearby PEs. To avoid the introduction of spurious long-range covariances associated with finite ensemble sizes, the background-error covariances are given compact support by means of a Hadamard (element by element) product with a three-dimensional canonical correlation function. The methodology and the MvEnKF configuration are discussed. It is shown that the regionalization of the background covariances; has a negligible impact on the quality of the analyses. The parallel algorithm is very efficient for large numbers of observations but does not scale well beyond 100 PEs at the current model resolution. On a platform with distributed memory, memory rather than speed is the limiting factor.

Keppenne, Christian L.↗

Land Model Testbed: Accelerating Development, Benchmarking and Analysis of Land Surface Models

A Land Model Testbed (LMT), designed to provide a computational framework for systematically assessing model fidelity and supporting rapid development of complex multiscale models, offers a general-purpose workflow for conducting large ensemble simulations of multiple land surface models, post-processing large volumes of model output, and evaluating model results. It leverages existing tools for launching model simulations and the International Land Model Benchmarking (ILAMB) package for assessing model fidelity through comparison with best-available observational datasets. Increased complexity and proliferation of uncertain parameters in process representations in land surface models has driven the need for frequent and intensive testing and evaluating of models to quantify uncertainties and optimize parameters such that results are consistent with observations. The LMT described here meets these needs by providing tools to run thousands of ensemble simulations simultaneously and post-process their output files, by automating execution of an enhanced version of ILAMB with site-specific benchmarks and multivariate functional relationships, and by offering ensemble diagnostics and a customizable dashboard for displaying model performance metrics and associated graphics. We envision the LMT capabilities will serve as a foundational computational resource for a proposed user facility focused on terrestrial multiscale model--data integration.

Sreepathi, Sarat↗

Aggregating available soil water holding capacity data for crop yield models

The total amount of water available to plants that is held against gravity in a soil is usually estimated as the amount present at -0.03 MPa average water potential minus the amount present at -1.5 MPa water potential. This value, designated available water-holding capacity (AWHC), is a very important soil characteristic that is strongly and positively correlated to the inherent productivity of soils. In various applications, including assessing soil moisture status over large areas, it is necessary to group soil types or series as to their productivity. Current methods to classify AWHC of soils consider only total capacity of soil profiles and thus may group together soils which differ greatly in AWHC as a function of depth in the profile. A general approach for evaluating quantitatively the multidimensional nature of AWHC in soils is described. Data for 902 soil profiles, representing 184 soil series, in Indiana were obtained from the Soil Characterization Laboratory at Purdue University. The AWHC for each of ten 150-mm layers in each soil was established, based on soil texture and parent material. A multivariate clustering procedure was used to classify each soil profile into one of 4, 8, or 12 classes based upon ten-dimensional AWHC values. The optimum number of classes depends on the range of AWHC in the population of oil profiles analyzed and on the sensitivity of a crop to differences in distribution of water within the soil profile.

Seubert, C. E.↗

Assessing the capability of EOS sensors in measuring ocean-atmosphere moisture exchange

As part of the Science Synergism Studies to identify interdisciplinary Scientific studies, which could be addressed by the Environmental Observing System (EOS), the techniques being developed to measure ocean-atmosphere moisture exchanges using satellite sensors were evaluated. Studies required to use sensors proposed for EOS were examined. A method has been developed to compute the moisture flux using the wind speed, sea surface temperature, and preciptable water measured by satellite sensors. It relies on a statistical model which predicts surface-level humidity from precipitable water. The Scanning Multichannel Microwave Radiometer (SMMR) measures all 3 parameters and was found to be sensitive to the annual cycle and large interannual variations such as the 1982 to 1983 El Nino. There are systematic differences between geophysical parameters measured by Nimbus/SMMR and in situ measurements. After quadratic trends and crosstalks were removed from the parameters through multivariate regressions, the latent heat fluxes computed from SMMR agree with those computed from ship reports to within 30 W/sq m. The poor quality of ship reports may be the cause of a portion of this scatter. Similar results are found using SEASAT/SMMR data. When the scatterometer winds were used instead of the SMMR winds, the difference between the satellite fluxes and the ship fluxes was reduced.

Liu, W. T.↗

Model-agnostic search for dijet resonances with anomalous jet substructure in proton–proton collisions at $\sqrt{s}$ = 13 TeV

This paper presents a model-agnostic search for narrow resonances in the dijet final state in the mass range 1.8-6 TeV. The signal is assumed to produce jets with substructure atypical of jets initiated by light quarks or gluons, with minimal additional assumptions. Search regions are obtained by utilizing multivariate machine-learning methods to select jets with anomalous substructure. A collection of complementary anomaly detection methods - based on unsupervised, weakly supervised, and semisupervised algorithms - are used in order to maximize the sensitivity to unknown new physics signatures. These algorithms are applied to data corresponding to an integrated luminosity of 138 fb -1 , recorded by the CMS experiment at the LHC, at a center-of-mass energy of 13 TeV. No significant excesses above background expectations are seen. Exclusion limits are derived on the production cross section of benchmark signal models varying in resonance mass, jet mass, and jet substructure. Many of these signatures have not been previously sought, making several of the limits reported on the corresponding benchmark models the first ever. When compared to benchmark inclusive and substructure-based search strategies, the anomaly detection methods are found to significantly enhance the sensitivity to a variety of models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for pair production of third-generation leptoquarks decaying into a bottom quark and a $\tau $-lepton with the ATLAS detector

Abstract A search for pair-produced scalar or vector leptoquarks decaying into a b -quark and a $$\tau $$ τ -lepton is presented using the full LHC Run 2 (2015–2018) data sample of 139 fb $$^{-1}$$ - 1 collected with the ATLAS detector in proton–proton collisions at a centre-of-mass energy of $$\sqrt{s} =13$$ s = 13 TeV. Events in which at least one $$\tau $$ τ -lepton decays hadronically are considered, and multivariate discriminants are used to extract the signals. No significant deviations from the Standard Model expectation are observed and 95% confidence-level upper limits on the production cross-section are derived as a function of leptoquark mass and branching ratio $$\mathcal {B}$$ B into a $$\tau $$ τ -lepton and b -quark. For scalar leptoquarks, masses below 1460 GeV are excluded assuming $$\mathcal {B}=100$$ B = 100 %, while for vector leptoquarks the corresponding limit is 1650 GeV (1910 GeV) in the minimal-coupling (Yang–Mills) scenario.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Bayesian approach to microwave precipitation profile retrieval

A multichannel passive microwave precipitation retrieval algorithm is developed. Bayes theorem is used to combine statistical information from numerical cloud models with forward radiative transfer modeling. A multivariate lognormal prior probability distribution contains the covariance information about hydrometeor distribution that resolves the nonuniqueness inherent in the inversion process. Hydrometeor profiles are retrieved by maximizing the posterior probability density for each vector of observations. The hydrometeor profile retrieval method is tested with data from the Advanced Microwave Precipitation Radiometer (10, 19, 37, and 85 GHz) of convection over ocean and land in Florida. The CP-2 multiparameter radar data are used to verify the retrieved profiles. The results show that the method can retrieve approximate hydrometeor profiles, with larger errors over land than water. There is considerably greater accuracy in the retrieval of integrated hydrometeor contents than of profiles. Many of the retrieval errors are traced to problems with the cloud model microphysical information, and future improvements to the algorithm are suggested.

Evans, K. Franklin↗

Transfer Learning-Based Independent Component Analysis

Understanding the underlying component structure is crucial for multivariate signal analysis. Among all the techniques that try to learn the latent structure, independent component analysis (ICA) is one of the most important and popular methods, which aims to extract independent components from multivariate signals and enables further analysis. For example, in electroencephalogram (EEG) analysis, artifacts filtering and disease detection are conducted based on the independent components of the signals. One critical challenge in existing ICA approaches is that the component extraction accuracy may degrade when the available data of a unit are limited. To address this issue, this paper proposes a transfer learning-based ICA method by innovatively transferring component distribution from a source domain, so that accurate component extraction results can be achieved even when only limited data are available in the target domain. To the best of our knowledge, this is the first work that leverages transfer learning to improve ICA accuracy with limited available data. In particular, we first extract all the independent components from the source domain by maximizing the log-likelihood function with a Newton-like method on a smooth manifold. Then for the target domain, the component with the largest negentropy is extracted in each round. To effectively leverage the knowledge from the source domain and to prevent the negative transfer, we try to find a component in the source domain that matches the component we are extracting. The probability density function of the matched component will then be used to improve the component extraction accuracy if such matched component can be found; otherwise, no knowledge will be transferred. Finally, numerical simulations and a case study with electrocardiogram (ECG) data are conducted, showing the effectiveness of the proposed method in transferring knowledge and reducing negative transfer.

42 ENGINEERING↗

Applying Two-Dimensional Correlation Spectroscopy and Principal Component Analysis to Understand How Temperature Affects the Neptunium(V) Absorption Spectrum

The visible-near infrared (Vis-NIR) electronic absorption spectrum of neptunium(V) (NpO 2 + ) comprises numerous f-f electronic transitions with mostly undocumented temperature dependencies. The effect of temperature on the absorption spectrum of the pentavalent neptunyl dioxocation (NpO 2 + ) is an important factor to consider with spectrophotometric applications but has often been overlooked. Optical Vis-NIR absorption spectra (400–1700 nm) of Np(V) (0.017–0.89 M) in 1 M nitric acid were evaluated with varying temperatures (T = 10–80 °C). The intensity, position, and overall shape of the bands were sensitive to interactions with the solvent and coordination environment. Numerous temperature-induced isosbestic points were identified resulting from dynamic, overlapping peak shifts. Spectral variations were characterized using principal component analysis (PCA) and 2D correlation spectroscopy (COS). 2D-COS revealed that the absorption band near 1095 nm likely consists of two bands centered near 1087 and 1096 nm, which cannot be explained by current computational methods. 2D-COS analysis also provided an unambiguous assignment of unresolved peaks in the visible region for comparison with computational predictions. PCA was used to identify nonlinearity in the spectral response at elevated Np(V) concentrations ≥ 0.5 M. This unique experimental data and interpretation will foster a deeper understanding of the absorption spectra for complex actinyl ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A survey on the structured singular value

The structured singular value, U, is an important linear algebra tool to study a class of matrix perturbation problems. It is useful for analyzing the robustness of stability and performance of uncertain, (nominally) linear systems. Computation of (M) is difficult, and usually, upper and lower bounds are all that can be reliably computed. Upper bounds give conservative estimates of the sizes of allowable perturbations. The maximum singular value of a matrix M is an upper bound for (M). As an upper bound, it can be improved by finding a transformations to the data (i.e. M) which do not change the structured singular value, but do reduce the maximum singular value. Typically, upper bound algorithms involve searches over sets of transformations to yield the tightest bound. Lower bound algorithms are intelligent searches for minimum-norm solutions to multivariable polynomial equations, and are based on various optimality conditions that hold at the global (and, unfortunately, some local) minima. The current methods to compute both of these types of bounds are reviewed. Theoretical justification and extensive numerical experience with the various algorithms are covered.

Packard, Andy↗

Temporal Analysis and Scene Change Detection in Multispectral Overhead Imagery

Scene change detection can be a tedious and time consuming process especially when concerning large geographical areas, and the process can be even more cumbersome when analyzing changes in an area over large spans of time. Developing a useful way to help analysts recognize at what points in time significant changes to a scene have occurred can allow them to better focus their efforts in characterizing events. Applications include: Facility monitoring, Construction chronology, Monitoring of vehicle/aircraft activity, Characterization of larger sequences of events. In large areas exceeding hundreds to thousands of square kilometers in size, it can be difficult localizing when scene changes have occurred. Analysts can spend hours going through imagery to try to identify new construction, monitor facility activities, monitor vehicle movement, etc. where the object of interest may only be a few square meters. Our goal is to help cut down this time by giving analysts change maps with hot spots of change, allowing them to focus on regions that have experienced actual change in time frames they're interested in. Additionally, by combining these change maps into layers within a data cube, analysts can examine the change maps from a temporal perspective, allowing events to be characterized over spans of time. By opening the data cube in an imaging software capable of separating the layers, we can analyze the change maps sequentially, allowing us to examine scene changes occurring over time. As an example, we examined overhead imagery from Planet Labs of what appears to be a parking lot on Fort Irwin over the course of a year using ENVI, a geospatial satellite imagery analysis software. Using ENVI, we generate a graph of changes over time, and notice a particular segment near the end of our analysis window where no changes are detected. Examination of the actual satellite imagery reveals that during this time span, the parking lot was empty. This could be due to facility shutdown for maintenance or upgrades, or possibly even total workforce/vehicle fleet movement. Information like this could help analysts better characterize events, as well as to help create clearer timelines in larger sequences of events. Workflow steps: - Collect multiple maps of the same AOI (Area of Interest) during a time span of interest; - Generate change maps from AOI maps; - Generate data cube from change maps. An analyst can use the data cube to help inspect an AOI for activities within a time span of interest. If an event of interest is discovered, the analyst can then refer to the maps corresponding to the appropriate dates and times in the data cube to see precisely what is transpiring. The biggest objective being worked on is improving the change detection methodology employed. We currently use PCA-EM (Principal Component Analysis with Expectation Maximization), but we are currently focusing on implementing IR-MAD (Iteratively Reweighted Multivariate Alteration Detection) to be used in conjunction with PCA-EM in an effort to decrease false positivity and noise in the change maps we generate.

42 ENGINEERING↗

Time-series forecasting using manifold learning, radial basis function interpolation, and geometric harmonics

We address a three-tier numerical framework based on nonlinear manifold learning for the forecasting of high-dimensional time series, relaxing the “curse of dimensionality” related to the training phase of surrogate/machine learning models. At the first step, we embed the high-dimensional time series into a reduced low-dimensional space using nonlinear manifold learning (local linear embedding and parsimonious diffusion maps). Then, we construct reduced-order surrogate models on the manifold (here, for our illustrations, we used multivariate autoregressive and Gaussian process regression models) to forecast the embedded dynamics. Finally, we solve the pre-image problem, thus lifting the embedded time series back to the original high-dimensional space using radial basis function interpolation and geometric harmonics. The proposed numerical data-driven scheme can also be applied as a reduced-order model procedure for the numerical solution/propagation of the (transient) dynamics of partial differential equations (PDEs). In conclusion, we assess the performance of the proposed scheme via three different families of problems: (a) the forecasting of synthetic time series generated by three simplistic linear and weakly nonlinear stochastic models resembling electroencephalography signals, (b) the prediction/propagation of the solution profiles of a linear parabolic PDE and the Brusselator model (a set of two nonlinear parabolic PDEs), and (c) the forecasting of a real-world data set containing daily time series of ten key foreign exchange rates spanning the time period 3 September 2001–29 October 2020.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗