Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “VECTOR 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The development of Gibbs's dyadic and implications for the gradient of a vector field

In this paper, we review the history of the dyadic as developed by Gibbs. This mathematical construct appeared in the second part of Gibbs's pamphlet on vector analysis (published in 1884), and it represented the first known development of a Cartesian theory of tensors. Gibbs made extensive use of the dyadic to express his theory of linear vector functions, that is, functions that acted on vectors and mapped them to new vectors. The dyadic proved to be a capable vehicle in Gibbs's hands, and his theory for dyadics (which we would now call second-order Cartesian tensors) was relatively advanced. The theory detailed notions such as the decomposition of vectors and conditions under which a tensor would have an inverse. While Gibbs's theory for linear operators expressed by dyadics was robust, it did not seem to garner the attention that the more conventional vector analysis (published in the first half of his pamphlet in 1881) did. Perhaps in part because of the general unfamiliarity with the dyadic, two distinct and conflicting definitions of the gradient of a vector field have arisen in the literature. The details of these differences in notation, possible reasons for the difference, and a potential resolution are proposed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

ZMPY3D: accelerating protein structure volume analysis through vectorized 3D Zernike moments and Python-based GPU integration

Abstract Motivation Volumetric 3D object analyses are being applied in research fields such as structural bioinformatics, biophysics, and structural biology, with potential integration of artificial intelligence/machine learning (AI/ML) techniques. One such method, 3D Zernike moments, has proven valuable in analyzing protein structures (e.g., protein fold classification, protein–protein interaction analysis, and molecular dynamics simulations). Their compactness and efficiency make them amenable to large-scale analyses. Established methods for deriving 3D Zernike moments, however, can be inefficient, particularly when higher order terms are required, hindering broader applications. As the volume of experimental and computationally-predicted protein structure information continues to increase, structural biology has become a “big data” science requiring more efficient analysis tools. Results This application note presents a Python-based software package, ZMPY3D, to accelerate computation of 3D Zernike moments by vectorizing the mathematical formulae and using graphical processing units (GPUs). The package offers popular GPU-supported libraries such as CuPy and TensorFlow together with NumPy implementations, aiming to improve computational efficiency, adaptability, and flexibility in future algorithm development. The ZMPY3D package can be installed via PyPI, and the source code is available from GitHub. Volumetric-based protein 3D structural similarity scores and transform matrix of superposition functionalities have both been implemented, creating a powerful computational tool that will allow the research community to amalgamate 3D Zernike moments with existing AI/ML tools, to advance research and education in protein structure bioinformatics. Availability and implementation ZMPY3D, implemented in Python, is available on GitHub (https://github.com/tawssie/ZMPY3D) and PyPI, released under the GPL License.

Lai, Jhih-Siang (ORCID:0000000156775890)↗
DOE OSTI2025

Temperature and dose effects on dislocation loops in self-ion irradiated high-purity iron

Body-centered cubic (BCC) Fe-based alloys are promising candidate materials for advanced nuclear reactors. However, a detailed understanding of irradiation induced dislocation loop microstructure development remains unresolved. It is a widespread belief that 〈001〉 loops become increasingly favorable over ½〈111〉 loops as irradiation temperature rises above ∼300 C. Unfortunately, the temperature effects on 〈001〉 loop have been primarily examined in in-situ irradiation on TEM thin foils but poorly explored on bulk Fe due to exceedingly limited experimental studies on bulk specimens, raising concerns about the potential influence of TEM thin foil artifacts on observed results. Here, in this study, we conducted experiments on ultra-high purity BCC Fe specimens irradiated with 6.7–8 MeV Fe ions over a wide temperature range on bulk samples. We investigated the effects of temperature (T irr = 250–500 °C), dose rate (10⁻⁵ to 10⁻³ dpa/s), and dose (0.35 to 3.5 dpa) on the formation and evolution of 〈001〉 and ½〈111〉 loops as well as cavity (void) formation. Post-irradiation Burgers vector analysis via g•b method on dislocation segments and loops revealed that 〈001〉 loop fraction does not show a monotonic positive correlation with irradiation temperature. Combined with previous and current theoretical as well as experimental findings, we explore the temperature effects on all existing models of 〈001〉 loop formation. We conclude that the prevailing reports regarding the dominance of 〈001〉 loops in Fe at elevated temperatures are mainly attributable to the loss of glissile ½〈111〉 clusters in TEM thin foil experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Temperature and dose rate effects on dislocation loop formation and decoration in self-ion irradiated high-purity Fe-Cr model alloys

Body-centered cubic (BCC) Fe-based alloys are promising candidate materials for advanced nuclear reactors. A widespread perception is that the addition of Cr favors the formation of ½ <111> dislocation loops. However, the Cr effects on dislocation loop formation/evolution remain unclear due to the limited experimental studies. In this study, we irradiate ultra-high purity bulk BCC Fe-xCr (x = 3, 5, and 8 wt%) with 6.7-8 MeV Fe ions over a wide range of temperatures (Tirr = 250-500 °C) and dose rates (10⁻⁵ to 10⁻³ dpa/s) to two midrange doses (0.35 and 3.5 dpa). The major irradiation-induced microstructures are loop decoration along non-screw dislocation lines at low temperatures and high dose rates as well as self-assembled loops at elevated temperatures. Post-irradiation Burgers vector analysis via g•b method on loops and electron tomography reveals that <001> preponderance is observable in all irradiated Fe-Cr samples at elevated temperature, resulting from stronger elastic energy originating from <001> edge.

Li, Yao [University of Tennessee (UT)]↗
DOE OSTI2024

Assessing decision boundaries under uncertainty

In order to make design decisions, engineers may seek to identify regions of the design domain that are acceptable in a computationally efficient manner. A design is typically considered acceptable if its reliability with respect to parametric uncertainty exceeds the designer’s desired level of confidence. Despite major advancements in reliability estimation and in design classification via decision boundary estimation, the current literature still lacks a design classification strategy that incorporates parametric uncertainty and desired design confidence. To address this gap, this paper offers a novel interpretation of the acceptance region by defining the decision boundary as the hypersurface which isolates the designs that exceed a user-defined level of confidence given parametric uncertainty. This work addresses the construction of this novel decision boundary using computationally efficient algorithms that were developed for reliability analysis and decision boundary estimation. The approach proposed in this paper is verified on two physical examples from structural and thermal analysis using Support Vector Machines and Efficient Global Optimization-based contour estimation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Swept-Lookback Deflectometry for High Performance Concentrating Solar Power Optical Metrology

This report describes an initial investigation into a proposed solution to the important problem of performing a detailed evaluation of heliostat optical performance, in situ in a heliostat field. Our approach is to place digital cameras in a position near the receiver where they look back toward the heliostat mirrors. The pixels of each camera sensor identify a set of small cells on the mirror surface, each corresponding to a “mixel.” By either passing reflected sunbeam over the camera or passing the camera through the reflected sunbeam, the cameras intercept sunlight reflected from each mixel. We then analyze the recorded video data to determine times when each mixel transitions from dark to light, and then back to dark. We then use these transitions to construct vectors from the camera to the mixel, and then from the mixel to the edge of the Sun at that moment. We then compute the surface normal at the mixel, which bisects the angle between these vectors. Performing this analysis for all mixels in the mirror yields a high-resolution map of slope across the mirror surface. We have implemented most of this process, successfully collecting data for an example heliostat facet and computing a preliminary estimated slope map. However, more work remains to complete this calculation, since certain factors and transformations are not yet included. Our observations so far support our hypothesis that such a system is possible, but we have not yet completed our quantitative evaluation of the concept.

14 SOLAR ENERGY↗
DOE OSTI2025

Complete definition of 𝑁 →Δ transition generalized parton distributions

We revisit the definition of the leading-twist chiral-even generalized parton distributions (GPDs) for 𝑁 →Δ baryon transitions. We identify and address deficiencies in previous definitions of the transition GPDs inspired by the transition form factors of the vector and axial vector currents. Through systematic analysis of all possible covariant structures, respecting discrete symmetries and the baryon spinor equations of motion, we derive complete sets of independent structures for the transition matrix elements of the vector and axial vector partonic operators. They contain additional structures proportional to the light-cone vector, corresponding to transition GPDs of vanishing first moment, which were not included in previous parametrizations. Their presence is confirmed independently by the light-front multipole expansion and the cross-channel SO(3) partial-wave analysis of the transition matrix elements. Our analysis provides a complete definition of the 𝑁 →Δ transition GPDs for use in theoretical and phenomenological studies.

baryons↗
DOE OSTI2026

Combination of searches for heavy vector boson resonances in proton-proton collisions at $\sqrt{s}=13$ TeV

A combined statistical analysis of searches for heavy vector boson resonances decaying into pairs of W, Z, or Higgs bosons, as well as into quark pairs $\left(\mathrm{q}\overline{\mathrm{q}},\mathrm{b}\overline{\mathrm{b}},\mathrm{t}\overline{\mathrm{t}},\mathrm{t}\overline{\mathrm{b}}\right)$ or lepton pairs ℓ + ℓ – , ℓ$\bar{v}$, with ℓ = e, μ, τ, is presented. The results are based on proton-proton collision data at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb −1 , collected by the CMS experiment from 2016 to 2018. No significant deviation from the expectations of the standard model is observed. The results are interpreted in the simplified heavy vector triplet (HVT) framework, setting 95% confidence level upper limits on the production cross sections and on the coupling strengths of the HVT bosons to standard model particles. The results exclude HVT resonances with masses below 5.5 TeV in a weakly coupled scenario, below 4.8 TeV in a strongly coupled scenario, and up to 2.0 TeV in the case of production via vector boson fusion. The combination provides the most stringent constraints to date on new phenomena predicted by the HVT model.

beyond Standard Model↗
DOE OSTI2025

SigTime: Learning and Visually Explaining Time Series Signatures

Understanding and distinguishing temporal patterns in time series data is essential for scientific discovery and decision-making. For example, in biomedical research, uncovering meaningful patterns in physiological signals can improve diagnosis, risk assessment, and patient outcomes. However, existing methods for time series pattern discovery face major challenges, including high computational complexity, limited interpretability, and difficulty in capturing meaningful temporal structures. Here, to address these gaps, we introduce a novel learning framework that jointly trains two Transformer models using complementary time series representations: shapelet-based representations to capture localized temporal structures and traditional feature engineering to encode statistical properties. The learned shapelets serve as interpretable signatures that differentiate time series across classification labels. Additionally, we develop a visual analytics system—SigTime—with coordinated views to facilitate exploration of time series signatures from multiple perspectives, aiding in useful insights generation. We quantitatively evaluate our learning framework on eight publicly available datasets and one proprietary clinical dataset. Additionally, we demonstrate the effectiveness of our system through two usage scenarios along with the domain experts: one involving public ECG data and the other focused on preterm labor analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

PLC Vulnerabilities and Mitigations

Programmable Logic Controllers (PLCs) are used extensively in many high-importance industrial and nonindustrial settings including controlling elevators, manufacturing machines, and utility facilities such as power and natural gas, however cybersecurity protection for them has been neglected. Within recent years, PLCs have been put under more security scrutiny and experts have advocated for changes from the addition of protocol encryption and network segmentation to intrusion detection systems on the PLCs themselves. While PLC security is critical, a large portion of the PLCs available today will never receive these changes due to being legacy or the difficulty of overhauling the security on existing systems. Due to the infeasibility of applying many recommended security measures towards currently available machines, we aim to provide realistic and affordable best practices for hardening PLCs. We will first conduct security analysis and consider attack vectors within our target PLC. Once we’ve analyzed the device’s security, we will evaluate a variety of mitigation methods and create guidelines to effectively reduce the threat posed by PLC attacks with minimal disruption to operations.

42 ENGINEERING↗
DOE OSTI2025

Validation of the DESI DR2 measurements of baryon acoustic oscillations from galaxies and quasars

The Dark Energy Spectroscopic Instrument (DESI) Data Release 2 (DR2) galaxy and quasar clustering data represents a significant expansion of data from Data Release 1 (DR1), providing improved statistical precision in baryon acoustic oscillation (BAO) constraints across multiple tracers, including bright galaxies, luminous red galaxies, emission line galaxies, and quasars. In this paper, we validate the BAO analysis of DR2. We present the results of robustness tests on the blinded DR2 data and, after unblinding, consistency checks on the unblinded DR2 data. All results are compared with those obtained from a suite of mock catalogs that replicate the selection and clustering properties of the DR2 sample. We confirm the consistency of DR2 BAO measurements with DR1 while achieving a reduction in statistical uncertainties due to the increased survey volume and completeness. The combined BAO precision, including both statistical and systematic errors, improves from ∼0.52% in DR1 to 0.30% in DR2—a factor of 1.7 gain. We assess the impact of analysis choices, including different data vectors (correlation function vs power spectrum), modeling approaches and systematics treatments, and an assumption of the Gaussian likelihood, finding that our BAO constraints are stable across these variations and assumptions with a few minor refinements to the baseline setup of the DR1 BAO analysis. We summarize a series of pre-unblinding tests that confirmed the readiness of our analysis pipeline, the final systematic errors, and the DR2 BAO analysis baseline. The successful completion of these tests led to the unblinding of the DR2 BAO measurements, ultimately leading to the DESI DR2 cosmological analysis, with their implications for the expansion history of the Universe and the nature of dark energy presented in the DESI key paper (companion paper).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Machine-learning-enabled on-the-fly analysis of RHEED patterns during thin film deposition by molecular beam epitaxy

Thin film deposition is a fundamental technology for the discovery, optimization, and manufacturing of functional materials. Deposition by molecular beam epitaxy (MBE) typically employs reflection high-energy electron diffraction (RHEED) as a real-time in situ probe of the growing film. However, the state-of-the-art for RHEED analysis during deposition requires human observation. Here, we present an approach using machine learning (ML) methods to monitor, analyze, and interpret RHEED images on-the-fly during thin film deposition. In the analysis workflow, RHEED pattern images are collected at one frame per second and featurized using a pretrained deep convolutional neural network. The feature vectors are then statistically analyzed to identify changepoints; these changepoints can be related to changes in the deposition mode from initial film nucleation to a transition regime, smooth film deposition, and in some cases, an additional transition to a rough, islanded deposition regime. The feature vectors are additionally analyzed via graph analysis and community classification. The graph is quantified as a stabilization plot, and we show that inflection points in the stabilization plot correspond to changes in the growth regime. The full RHEED analysis workflow is termed RHAAPsody and includes data transfer and output to a visual dashboard. We demonstrate the functionality of RHAAPsody by analyzing the precaptured RHEED images from epitaxial depositions of anatase TiO2 on SrTiO3(001) and show that the analysis workflow can be executed in less than 1 s. Our approach shows promise as one component of ML-enabled real-time feedback control of the MBE deposition process.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Machine Learning–Augmented Laser-Induced Breakdown Spectroscopy for Spectral Discrimination of Iron Oxalates

Enhanced characterization and phase identification of post-PUREX Pu Oxalates (PuOXA) are pivotal for nonproliferation and pre-detonation nuclear forensics. Despite significant advances in the characterization of PuO 2 samples, little is known about the impact of both the chemical structure and oxidation states of PuOXA (i.e., Pu(III) and Pu(IV)) have on optical emission signatures. Here, we demonstrate the analytical capabilities of laser-induced breakdown spectroscopy (LIBS) applied to Fe(II) and Fe(III) oxalate samples as surrogates for PuOXA, highlighting the discriminating features in the LIBS emission spectra arising from differences in the oxidation states within mixed FeOXA samples. We report the enhancement of spectral feature selection using Principal Component Analysis (PCA), which enables the analytical superiority of machine learning algorithms such as Linear Discriminant Analysis (LDA), Quadratic Discriminant Analysis (QDA), Partial Least Squares Regression (PLSR), Support Vector Regression (SVR), and Random Forest Regression (RFR) over conventional univariate techniques for phase discrimination and chemometric analysis. Cluster analysis revealed how both matrix effects and laser ablation influence cluster separability by introducing spectral artifacts that misdirect the maximization of variance. PCA-selected emission lines were used in the regression models, demonstrating that both univariate and multivariate linear regression models (i.e., PLSR and SVR) can achieve acceptable performance, with machine learning models outperforming conventional calibration regressions. Furthermore, the application of non-linearly activated PCA-selected emission lines illustrates how simplifying the data while retaining captured variance enables the use of less complex and more computationally efficient models. Furthermore, this is particularly evident in the underperformance of RFR, which suffers from increased computational costs and overfitting owing to its high complexity.

Oxalates↗
DOE OSTI2024

NSTXU Diagnostic Disruption Dynamic Loading Represented by Response Spectra

This article presents the results of transient dynamic simulations of loads due to disruption eddy currents on the NSTXU vacuum vessel. Dynamic loading at diagnostic mounting locations is expressed as response spectra derived from the time history results of the dynamic structural simulations of a variety of disruption scenarios. The disruption simulations draw on a history of the project assessments of worst case disruptions for specific components. Major efforts to assess disruption loading have included the vacuum vessel which is the major structural support for the machine, as well as the passive plates (PPs), high harmonic fast wave (HHFW) antenna, and centerstack casing. Each one of these efforts included transient electromagnetic simulations producing time-dependent eddy current Lorentz loads (and in some cases halo loads) which then were applied to time-dependent structural dynamic analyses intended to obtain the proper dynamic amplification factors. In some instances, the EM model and structural model were identical allowing direct transfer of EM forces to the structural model. In other cases, the EM and structural model were not identical and the vector potential (VP) transfer method was used. The results files from these analyses were available (or re-run) to post process in ANSYS Classic time history postprocessor. In conclusion, the ANSYS command is used to create response spectra from time history data at desired points on the vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Probing the Higgs boson CP properties in vector-boson fusion production in the H → τ + τ − channel with the ATLAS detector

The CP properties of the Higgs boson are studied in the vector-boson fusion production mode. The analysis exploits the decay mode of the Higgs boson into two τ-leptons using 140 fb −1 of proton-proton collision data at $\sqrt{s}$ = 13 TeV collected by the ATLAS experiment at the Large Hadron Collider. Results are obtained using the Optimal Observable method. CP-violating interactions between the Higgs boson and electroweak gauge bosons are considered in the effective field theory framework, with the interaction strength described in the HISZ basis by $\tilde{d}$, and in the Warsaw basis by c H$\tilde{W}$ , c H$\tilde{B}$ , and c H$\tilde{W}$B . No deviations relative to the Standard Model are observed, and limits are obtained on the strength parameters. The $\tilde{d}$ parameter is constrained to the interval [−0.012, 0.044] at the 95% confidence level while c H$\tilde{W}$ is constrained to [−0.24, 0.83], when considering both linear and quadratic effects of physics beyond the Standard Model.

Hadron-Hadron Scattering↗
DOE OSTI2024

Three-dimensional continuum point cloud method for large deformation and its verification

This study presents a strong form based meshfree collocation method, which is named Continuum Point Cloud Method, to solve nonlinear field equations derived from classical mechanics for deformed bodies in three-dimensional Euclidean space. The method and its implementation are benchmarked against a nonlinear vector field using manufactured solutions. The analysis of mechanical fields firstly focuses on the study of St. Venant Kirchhoff and compressible neo-Hookean materials. Results for various initial boundary value problems are presented, including benchmark cases involving unidirectional tension and simple shear. Subsequently, the study concludes with an analysis of a displacement-controlled simulation of a compressible neo-Hookean material, specifically a bar that is pulled to 50% of its original length and rotated 90°. The pure tension case yields a 1.5% error in displacement between computed and expected values and a combined tension and torsion loading case provides further insight into material behavior under complex loading conditions. The resulting normal axial and transverse stress-strain curves are also presented. Lastly, the consistency and robustness of the proposed nonlinear numerical schemes are successfully demonstrated through various numerical experiments.

Compressible neo-Hookean materials↗