Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Vectorization”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Generalizing mkFit and its Application to HL-LHC

mkFit is an implementation of the Kalman filter-based track reconstruction algorithm that exploits both thread- and data-level parallelism. In the past few years the project transitioned from the R&D phase to deployment in the Run-3 offline workflow of the CMS experiment. The CMS tracking performs a series of iterations, targeting reconstruction of tracks of increasing difficulty after removing hits associated to tracks found in previous iterations. mkFit has been adopted for several of the tracking iterations, which contribute to the majority of reconstructed tracks. When tested in the standard conditions for production jobs, speedups in track pattern recognition are on average of the order of 3.5x for the iterations where it is used (3-7x depending on the iteration). Multiple factors contribute to the observed speedups, including vectorization and a lightweight geometry description, as well as improved memory management and single precision. Efficient vectorization is achieved with both the icc and the gcc (default in CMSSW) compilers and relies on a dedicated library for small matrix operations, Matriplex, which has recently been released in a public repository. While the mkFit geometry description already featured levels of abstraction from the actual Phase-1 CMS tracker, several components of the implementations were still tied to that specific geometry. We have further generalized the geometry description and the configuration of the run-time parameters, in order to enable support for the Phase-2 upgraded tracker geometry for the HL-LHC and potentially other detector configurations. The implementation strategy and high-level code changes required for the HL-LHC geometry are presented. Speedups in track building from mkFit imply that track fitting becomes a comparably time consuming step of the tracking chain.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Extending Conformal Prediction to Hidden Markov Models with Exact Validity via de Finetti’s Theorem for Markov Chains

Conformal prediction is a widely used method to quantify uncertainty in settings where the data is independent and identically distributed (IID), or more generally, exchangeable. Conformal prediction takes in a pre-trained classifier and a calibration dataset as inputs, and returns a function which maps feature vectors to subsets of classes. The output of the returned function for a new feature vector is guaranteed to contain the true class with a pre-specified confidence. Despite its success and usefulness in IID settings, extending conformal prediction to non-exchangeable (e.g., Markovian) data in a manner that provably preserves all desirable theoretical properties has largely remained an open problem. As a solution, we extend conformal prediction to the setting of a Hidden Markov Model (HMM) with unknown parameters. The key idea behind the proposed method is to partition the non-exchangeable Markovian data from the HMM into exchangeable blocks by exploiting the de Finetti’s Theorem for Markov Chains discovered by Diaconis and Freedman (1980). The permutations of the exchangeable blocks are then viewed as randomizations of the observed Markovian data from the HMM. The proposed method provably retains all desirable theoretical guarantees offered by the classical conformal prediction framework and is general enough to be useful in many sequential prediction problems.

Nettasinghe, Don Buddhika Wijayantha↗
DOE OSTI2024

Constraining Cosmology with Simulation-based inference and Optical Galaxy Cluster Abundance

We test the robustness of simulation-based inference (SBI) in the context of cosmological parameter estimation from galaxy cluster counts and masses in simulated optical datasets. We construct ``simulations'' using analytical models for the galaxy cluster halo mass function (HMF) and for the observed richness (number of observed member galaxies) to train and test the SBI method. We compare the SBI parameter posterior samples to those from an MCMC analysis that uses the same analytical models to construct predictions of the observed data vector. The two methods exhibit comparable performance, with reliable constraints derived for the primary cosmological parameters, ($\Omega_m$ and $\sigma_8$), and richness-mass relation parameters. We also perform out-of-domain tests with observables constructed from galaxy cluster-sized halos in the Quijote simulations. Again, the SBI and MCMC results have comparable posteriors, with similar uncertainties and biases. Unsurprisingly, upon evaluating the SBI method on thousands of simulated data vectors that span the parameter space, SBI exhibits worsened posterior calibration metrics in the out-of-domain application. We note that such calibration tests with MCMC is less computationally feasible and highlight the potential use of SBI to stress-test limitations of analytical models, such as in the use for constructing models for inference with MCMC.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Systems and methods for architecture-independent binary code analysis

Binaries configured for execution within respective computing environments may be disassembled into architecture-specific intermediate (AIL) representations. The AIL representations may be converted into canonical intermediate language (CIL) representations. The CIL representations may comprise normalized, architecture-independent code configured to characterize functionality of respective components of a binary (e.g., respective functions or the like). Feature vectors may be extracted from the CIL representations. The feature vectors may be used to identify components of respective binaries, assign security classifications to the binaries, and/or the like.

Haile, Jedediah T.↗
DOE OSTI2024

Light Detecting and Ranging (LIDAR) for in-situ heliostat optical error assessment

A system and method for optical assessment of a heliostat includes obtaining a point cloud data representing an image of the heliostat; isolating the data; filtering and fitting the filtered heliostat data to a bounding box; translating the heliostat data to a plane to aid in segmentation; segmenting a plurality of facets of the heliostat fitting each of the segmented facets to a respective plane; generating normal vectors characterizing each of the plurality of facets; and calculating a canting angle associated with each respective facet of the plurality of facets. A heliostat with mirrored facets and a scanner are provided. The scanner captures point cloud data representing the heliostat, which is segmented for each facet. Normal vectors characterize the facets and a canting angle is calculated for the respective facet.

Small, Daniel E.↗
DOE OSTI2024

Methods and compositions for obtaining useful plant traits

The present invention provides methods for obtaining plants that exhibit useful traits by perturbation of plastid function in plant rootstocks and grafting the rootstocks to scions. Methods for identifying genetic loci that provide for useful traits in plants and plants produced with those loci are also provided. In addition, plants that exhibit the useful traits, parts of the plants including seeds, and products of the plants are provided as well as methods of using the plants. Recombinant DNA vectors and transgenic plants comprising those vectors that provide for plastid perturbation are also provided.

Mackenzie, Sally↗
DOE OSTI2024

Multiplex characterization of microbial traits using dual barcoded nucleic acid fragment expression library

Disclosed herein are barcoded expression libraries comprising a plurality of expression vectors, wherein each expression vector comprises a nucleic acid fragment flanked by a first barcode and a second barcode. Further disclosed herein are methods of making the barcoded expression libraries and methods of conducting functional analysis using the barcoded expression libraries.

Mutalik, Vivek K.↗
DOE OSTI2024

Uncertainty-refined image segmentation under domain shift

Digital image segmentation is provided. The method comprises training a neural network for image segmentation with a labeled training dataset from a first domain, wherein a subset of nodes in the neural net are dropped out during training. The neural network receives image data from a second, different domain. A vector of N values that sum to 1 is calculated for each image element, wherein each value represents an image segmentation class. A label is assigned to each image element according to the class with the highest value in the vector. Multiple inferences are performed with active dropout layers for each image element, and an uncertainty value is generated for each image element. Uncertainty is resolved according to expected characteristics. The label of any image element with an uncertainty above a threshold is replaced with a new label corresponding to a segmentation class based on domain knowledge.

Martinez, Carianne↗
DOE OSTI2022

Fully Distributed Multi-Material Magnetic Sensing Structures for Multiparameter DAS Applications

This dissertation demonstrates the first of its kind distributed magnetic field sensor based on a fiber optic distributed acoustic sensing (DAS) scheme. Ferromagnetic nickel and Metglas® were dispersed internally within a fiber optic preform and then drawn on an in-house fiber optic draw tower to lengths in the kilometers. Due to the close proximity of the ferromagnetic metals and fiber optic core, the magnetostrictive strain response of the ferromagnetic materials when exposed to a magnetic field would perturbate within the fiber cladding and transfer that strain, internally, to the fiber optic core. Strain resulting from the magnetostrictive effect allows the DAS based sensor to accurately translate strain into readable magnetic field data. Due to the high sensitivity seen in this sensor design, multiparameter sources, acoustic and magnetic fields, were tested and validated and a three dimensional magnetic-field vector sensor was proposed. Numerical analysis of the novel sensor design was first implemented using COMSOL Multiphysics, where inputs such as magnetostrictive element shape, size, distance, and number were first investigated. Upon optimizing system constraints, the sensor design was further modified such that single mode operation was consistent across multiple fiber draws while retaining high strain transfer from the ferromagnetic elements to the fiber optic core. Ferromagnetic material selection was evaluated as a function of the saturation magnetostriction constants and a total of 4 modules were used to fully characterize the complex physics involved in this sensor design. All fabrication and testing were performed in-house using a full scale 3-story fiber draw tower and custom environmental testing stations to imitate naturally occurring events such as magnetic or acoustic point sources. A unique stacking method was used to embed ferromagnetic nickel and Metglas® into a fiber optic preform which when combined with a custom fiber draw process resulted in consistent multi-material fibers drawn to lengths of 1-km. In-house testing facilities included different types of electromagnetic generators, in addition to a soil test bed, and an outdoor test bed which allowed 100 meters of fiber to be tested simultaneously. All tested sensors demonstrated high strain transfer capabilities on the order of 0.01-10 μϵ depending on the materials used, ferromagnetic rod number, and core to metal spacing. Due to the sensitivity of the system the difference between AC and DC was distinct, and directional magnetostriction was studied. Transverse and longitudinal magnetic wave propagation was controlled through a solenoid and rectangular Helmholtz coil, both built in-house. A three-dimensional magnetic field vector sensor was proposed due to the success of the magnetic field sensor, and a design was proposed and initially tested to validate direction as a function of field strength and distance. To summarize, this dissertation explores the first fully distributed magnetic field sensor using DAS based techniques and one of the first multi-material fiber draw processes which can produce consistent single mode fiber up to 1-km. Due to extensive FEA modeling, multiple iterations of the magnetic sensor were fully characterized and an equation describing the relationship between sensor design and strain transfer has been created and validated experimentally. Multi-parameter tests including acoustic and magnetic fields were implemented and an algorithm was developed to separate the mixed signals. Finally, a test was performed to demonstrate the feasibility of sensing magnetic fields directionally. Cumulative results demonstrate a high-quality sensor alternative to current designs which may surpass other magnetic sensors due to innate multi-parameter capabilities, in addition to the inexpensive production cost and extremely long operating lengths.

02 PETROLEUM↗
DOE OSTI2026

$B \to π$, $B_{(s)} \to D_{(s)}$ from 2+1+1 Flavor Lattice QCD

We present a lattice-QCD calculation of the hadronic form factors for $B$-meson semileptonic decays computed using the highly improved staggered quark action for both valence and sea quarks on the MILC collaborations 2+1+1-flavor ensembles with lattice spacing ranging from 0.09 fm to 0.03 fm, many with physical pion masses On our finest ensembles, we compute the form factors directly at the physical $b$-quark mass. We discuss the computational setup and analysis strategies for two- and three-point correlation functions. For $B_{(s)} \to D_{(s)}$ we present preliminary results of chiral-continuum fits for the scalar and vector form factors. The goal of this project is a percent-level determination of the scalar and vector form factors to enable high-precision determinations of $|V_{ub}|$ and $|V_{cb}|$. This work fits into a broader program of lattice-QCD studies of weak $B$-meson decays by the Fermilab Lattice and MILC Collaborations.

Cassar, Nicholas [Colorado State U.]↗
DOE OSTI2026

Feature-Based PMU Event Classification under Variable PMU Participation and Overlapping Events

This paper is the basis for a presentation help at the 2026 Georgia Tech Fault & Disturbance Analysis Conference, which can be found at OSTI # 3168287 Paper Abstract—Phasor Measurement Units (PMUs) stream time synchronized, high-resolution measurements from the grid, enabling data-driven techniques for event detection and classification. Accurate event classification improves grid reliability and stability. Events can be detected by varying numbers of PMUs and exhibit different durations depending on the event type. This variability challenges standard classifiers that require uniform input sizes. Moreover, multiple events may coincide, which increases classification complexity. Standard classifiers assign each instance to the class with the highest predicted probability, whereas overlapping events may exhibit comparable probabilities across multiple classes. In this study, to handle data size variability, we extract a wide range of time–frequency domain features from all available PMUs for each event into a fixed-length vector, facilitating the application of standard machine learning classifiers, including Random Forest, XGBoost, LightGBM, Support Vector Machine, and Multilayer Perceptron. To account for overlapping events, a probabilistic post-processing step is applied. For a given data instance, if multiple predicted class probabilities exceed 30% and the differences between them are less than 10%, the event is assigned to multiple classes. Experiments using real-world PMU data demonstrate that the Random Forest and XGBoost models achieve the highest accuracy, while the proposed post-processing method yields perfect classification performance on external unseen test sets.

Nematirad, Reza [Danovo Energy Solutions]↗
DOE OSTI2026

Two-Tower Quantum Matrix Chain Multiplication: Trading Qubits for Depth

Matrix chain multiplication -- computing $\mathcal{W} = M^{(0)}\cdots M^{(K-1)}$ where $M^{(k)} \in \mathbb{R}^{P_k \times P_{k+1}}$-- arises in scientific computing, machine learning, and graph analysis. Despite the importance of this problem, for chains of distinct matrices, the classical number of operations grows linearly with the chain length $K$ and polynomially in the matrix dimensions. We present \emph{Two-Tower Matrix Multiplication}, a quantum subroutine that encodes the product $\mathcal{W}$ of the $K$ matrices into a quantum state in circuit depth $\mathcal{O}(\max_{k} \mathrm{polylog} (P_k P_{k+1}))$, which is independent of~$K$ within the QRAM-based state-preparation model, whereas the qubit count is $\mathcal{O}\bigl(\sum_{k} \log P_k \bigr)$; the total gate count remains linear in $K$, so the gain is in the circuit depth. The construction interleaves state-preparation operators across two layers; within each layer, all operators act on disjoint registers and execute in parallel. This subroutine can be specialized for the chain-vector case, which computes the product of $K-1$ matrices applied to a vector. We prove the correctness of the subroutine for all $K$ and provide two implementations using the Qiskit and QCLAB frameworks. The subroutine is applicable to any downstream quantum algorithm that operates on a matrix encoded in the statevector, including norm estimation, graph-matrix powers, linear system solving, and quantum machine learning kernels.

Antonioli, Giacomo [Pisa U.] (ORCID:00090000668703↗
DOE OSTI2022

Cyber-Attack Identification of Synchrophasor Data Via VMD and Multifusion SVM

A large amount of synchrophasor data in the wide area measurement system (WAMS) needs to be collected and transmitted to the phasor data concentrator, thereby increasing the possibility of being attacked by hackers. The attacked data are therefore hidden into the normal synchrophasor data so that the synchrophasor data based application will be affected. To remedy this problem, an identification framework is proposed to detect the data cyber-attack in WAMS utilizing variational mode decomposition (VMD) and multifusion support vector machine (MSVM). First, VMD is used to transform the attacked data into multiple modal components. Thereafter, a novel MSVM is employed to classify the deterministic features using the proposed linear combined multikernel (LCM). Further, this LCM can fuse multiple types of features, including the time, frequency, and statistical domains of the synchrophasor data. Utilizing the actual data from FNET/GridEye, different experiments are conducted under multiple attack strengths and types. The results demonstrate that the identification framework has higher precision and robustness compared with other conventional classifiers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Direct Imaging of Antiferromagnet-Ferromagnet Phase Transition in van der Waals Antiferromagnet CrSBr

The advent of van der Waals (vdW) ferromagnetic (FM) and antiferromagnetic (AFM) materials offers unprecedented opportunities for spintronics and magneto-optic devices. Combining magnetic Kerr microscopy and density functional theory calculations, the AFM-FM transition is investigated and a surprising abnormal magneto-optic anisotropy in vdW CrSBr associated with different magnetic phases (FM, AFM, or paramagnetic state) is discovered. This unique magneto-optic property leads to different anisotropic optical reflectivity from different magnetic states, permitting direct imaging of the AFM Néel vector orientation and the dynamic process of the AFM-FM transition within a magnetic field. Using Kerr microscopy, not only the domain nucleation and propagation process is imaged but also the intermediate spin-flop state in the AFM-FM transition is identified. In conclusion, the unique magneto-optic property and clear identification of the dynamics process of the AFM-FM phase transition in CrSBr demonstrate the promise of vdW magnetic materials for future spintronic technology.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Anatomy of an agricultural antagonist: Feeding complex structure and function of three xylem sap‐feeding insects illuminated with synchrotron‐based 3D imaging

Abstract Many insects feed on xylem or phloem sap of vascular plants. Although physical damage to the plant is minimal, the process of insect feeding can transmit lethal viruses and bacterial pathogens. Disparities between insect‐mediated pathogen transmission efficiency have been identified among xylem sap‐feeding insects; however, the mechanistic drivers of these trends are unclear. Identifying and understanding the structural factors and associated integrated functional components that may ultimately determine these disparities are critical for managing plant diseases. Here, we applied synchrotron‐based X‐ray microcomputed tomography to digitally reconstruct the morphology of three xylem sap‐feeding insect vectors of plant pathogens: Graphocephala atropunctata (blue‐green sharpshooter; Hemiptera, Cicadellidae) and Homalodisca vitripennis (glassy‐winged sharpshooter; Hemiptera, Cicadellidae), and the spittlebug Philaenus spumarius (meadow spittlebug; Hemiptera, Aphrophoridae). The application of this technique revealed previously undescribed anatomical features of these organisms, such as key components of the salivary complex. The visualization of the 3D structure of the precibarial valve led to new insights into the mechanism of how this structure functions. Morphological disparities with functional implications between taxa were highlighted as well, including the morphology and volume of the cibarial dilator musculature responsible for extracting xylem sap, which has implications for force application capabilities. These morphological insights will be used to target analyses illuminating functional differences in feeding behavior.

3D imaging↗
DOE OSTI2023

Measurement of e + e − → ωπ+π− cross section at $\sqrt{s}$ = 2.000 to 3.080 GeV

A partial wave analysis on the process e⁺e⁻ → ωπ⁺π⁻ is performed using 647 pb -1 of data sample collected by using the BESIII detector operating at the BEPCII storage ring at center-of-mass (c.m.) energies from 2.000 GeV to 3.080 GeV. The Born cross section of the e⁺e⁻ → ωπ⁺π⁻ process is measured, with precision improved by a factor of 3 compared to that of previous studies. A structure near 2.25 GeV is observed in the energy-dependent cross sections of e⁺e⁻ → ωπ⁺π⁻ and ωπ 0 π 0 with a statistical significance of 7.6σ, and its determined mass and width are 2232 ± 19 ± 27 MeV/c 2 and 93 ± 53 ± 20 MeV, respectively, where the first and second uncertainties are statistical and systematic, respectively. By analyzing the cross sections of subprocesses e⁺e⁻ → ωf 0 (500), ωf 0 (980), ωf 0 (1370), ωf 2 (1270), and b1(1235)π, a structure, with mass M = 2200 ± 11 ± 17 MeV/c 2 and width Γ = 74 ± 20 ± 24 MeV, is observed with a combined statistical significance of 7.9σ. The measured resonance parameters will help to reveal the nature of vector states around 2.25 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

The deconstruction of flavor in the privately democratic Higgs sector

The Standard Model (SM) of particle physics fails to explain the observed hierarchy in fermion masses or the origin of fermion-flavor structure. We construct a model to explain these observations in the quark sector. We introduce a spectrum of new particles consisting of six of each — massive singlet vector-like quarks (VLQs), singlet scalars, and SU(2)-doublet scalars. SM quark masses are generated when the neutral components of the SU(2)-doublet scalars acquire non-zero vacuum expectation values (VEVs). We impose global symmetries to ensure that Yukawa couplings stay roughly flavor diagonal and democratic (of the same order), as well as to suppress tree-level flavor-changing neutral currents. Quark-mass hierarchy then follows from a hierarchy in scalar VEVs. The singlet scalars also acquire weak-scale VEVs. Together with the VLQs, they act as messengers between different generations of quarks in the SM. These messenger particles are responsible for generating the elements of the Cabibbo-Kobayashi-Masakawa (CKM) matrix which depend on the ratios of the singlet VEVs and VLQ masses. Constructed this way, the CKM matrix is found to be independent of the SM fermion masses. Using the measured values of the CKM matrix elements and assuming order-one couplings, we derive constraints on the masses of the VLQs and discuss prospects for probing our model in the near future.

CKM Parameters↗
DOE OSTI2024

Measurement and interpretation of same-sign W boson pair production in association with two jets in $pp$ collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

This paper presents the measurement of fiducial and differential cross sections for both the inclusive and electroweak production of a same-sign W-boson pair in association with two jets (W ± W ± jj) using 139 fb –1 of proton-proton collision data recorded at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV by the ATLAS detector at the Large Hadron Collider. The analysis is performed by selecting two same-charge leptons, electron or muon, and at least two jets with large invariant mass and a large rapidity difference. The measured fiducial cross sections for electroweak and inclusive W ± W ± jj production are 2.92 ± 0.22 (stat.) ± 0.19 (syst.) fb and 3.38 ± 0.22 (stat.) ± 0.19 (syst.) fb, respectively, in agreement with Standard Model predictions. The measurements are used to constrain anomalous quartic gauge couplings by extracting 95% confidence level intervals on dimension-8 operators. A search for doubly charged Higgs bosons H ±± that are produced in vector-boson fusion processes and decay into a same-sign W boson pair is performed. The largest deviation from the Standard Model occurs for an H ±± mass near 450 GeV, with a global significance of 2.5 standard deviations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗