Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Principal Component Analysis (PCA)”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

3-D Geologic Controls of Hydrothermal Fluid Flow at Brady geothermal field, Nevada, USA

In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a throughgoing network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density.

58 GEOSCIENCES↗
DOE OSTI2026

Effects of 9.5 years warming on SOC concentration and composition in bulk soil and density fractions

Original data of whole-soil warming experiment after 9.5 years at Blodgett Forest Research Station. The Blodgett Forest is a mixed coniferous temperate forest with Mediterranean climate. The annual air temperature is 12.5℃ and the annual precipitation is 1774 mm yr-1- The soil is mesic ultic Alfisol of granitic origin, equivalent to Dystric Cambisol according to The World Reference Base for Soil Resources (WRB) system. The soil is warmed down to 1 m at + 4℃ by vertically installed heating cables. At the time of soil sampling on 1 May 2023, the whole-soil warming experiment had been running for approximately 9.5 years, from January 2014 to May 2023. The dataset includes: - Bulk_EA: C, N content, δ13C, and CN ratio of bulk soil; - Density_fractionation: organic carbon concentration, δ13C, and C/N ratio of free light fraction (fLF), occluded ligh fraction (oLF), and heavy fraction (HF); - PCA_DRIFT_AUC: original data of area under the curve (AUC) values of eight carbon bond types integrated on diffuse reflectance infrared fourier transform spectroscopy for each soil sample and soil fraction, which are consequently used for principal component analysis (PCA); - DRIFTS_stability_index: the calculation of aliphatic C–H (3000–2800 cm-1) to aromatic C=C (1670–1600 cm-1) ratios for each bulk soil sample and soil fraction. All data are provided in CSV format and can be viewed using Microsoft Excel.

Climate change↗
DOE OSTI2023

Data-driven analysis and prediction of wastewater treatment plant performance: Insights and forecasting for sustainable operations

Here this study presents a comprehensive performance and forecasting analysis of the As-Samra wastewater treatment plant (WWTP) in Jordan, with two main objectives. Firstly, a thorough evaluation of the plant's performance is conducted. The analysis involves independently assessing historical operational conditions, plant production, and their statistical correlations using various statistical techniques. The second objective focuses on developing a data-driven forecasting approach to predict the plant's production one month in advance, using multiple machine learning models. The results highlight the effectiveness of principal component analysis (PCA) in simplifying operational data, revealing distinct operational clusters, and identifying seasonal production patterns while showing correlations between operational conditions and overall power production. The support vector machine (SVM) forecasting model emerged as the top performer, showcasing the potential of a hybrid forecasting approach. The findings offer valuable perspectives for enhancing operational efficiency, refining production planning, and ultimately improving the environmental impact of the plant.

42 ENGINEERING↗
DOE OSTI2022

Co-optimized machine-learned manifold models for large eddy simulation of turbulent combustion

Many modeling approaches in large eddy simulation (LES) of turbulent combustion employ a projection of the thermochemical state onto a low-dimensional manifold within state space to reduce the number of transported variables and hence computational cost. Flamelet-generated manifolds (FGM) is an example of a well-established, physics-based approach, but increasingly, principal component analysis (PCA) is being used as a data-driven method for generating manifold models. For both approaches, the nonlinear relationship between the location on the predefined manifold and the outputs of interest, such as reaction rates, can be tabulated or encoded in a neural network. This work proposes a new approach for manifold modeling that extends these existing approaches. A modified neural network structure simultaneously encodes the definition of the manifold variables, the nonlinear mapping, and the subfilter closure for LES. This allows all three of these aspects of the model to be co-optimized, generating a model from any source of combustion thermochemical state data. The manifold parameterizing variables are constrained to be linear combinations of species, as in FGM and PCA-based models, to aid in interpretability and implementation. For LES, subfilter variances of the manifold variables are also included as inputs. Two types of a priori analysis are performed to evaluate the new approach. In the first, the model is trained on data from one-dimensional premixed flames. In this case, the approach recovers the behavior of flamelet-based manifold approaches, and in fact slightly improves performance by identifying an optimized progress variable. The approach is also applied to data from direct numerical simulations of spherical ignition kernels in isotropic turbulence. For any specified manifold dimensionality, the new approach provides substantially lower prediction errors than a PCA-based model developed from the same data set. Additionally, the LES formulation of the new approach can provide accurate predictions for filtered reaction rates across a variety of filter widths.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Beyond PCA: Additional Dimension Reduction Techniques to Consider in the Development of Climate Fingerprints

Abstract Dimension reduction techniques are an essential part of the climate analyst’s toolkit. Due to the enormous scale of climate data, dimension reduction methods are used to identify major patterns of variability within climate dynamics, to create compelling and informative visualizations, and to quantify major named modes such as El Niño–Southern Oscillation. Principal components analysis (PCA), also known as the method of empirical orthogonal functions (EOFs), is the most commonly used form of dimension reduction, characterized by a remarkable confluence of attractive mathematical, statistical, and computational properties. Despite its ubiquity, PCA suffers from several difficulties relevant to climate science: high computational burden with large datasets, decreased statistical accuracy in high dimensions, and difficulties comparing across multiple datasets. In this paper, we introduce several variants of PCA that are likely to be of use in climate sciences and address these problems. Specifically, we introduce non-negative , sparse , and tensor PCA and demonstrate how each approach provides superior pattern recognition in climate data. We also discuss approaches to comparing PCA-family results within and across datasets in a domain-relevant manner. We demonstrate these approaches through an analysis of several runs of the E3SM climate model from 1991 to 1995, focusing on the simulated response to the Mt. Pinatubo eruption; our findings are consistent with a recently identified stratospheric warming fingerprint associated with this type of stratospheric aerosol injection.

Weylandt, Michael↗
DOE OSTI2021

Accurate Thermochemistry of Complex Lignin Structures via Density Functional Theory, Group Additivity, and Machine Learning

A molecular-level understanding of lignin structures and bond dissociation energies could facilitate depolymerization technologies. Still, this information is currently limited due to the lack of databases and the simplification of surrogate models. Here, substitution effects on seven common linkages in lignin polymers are systematically investigated. An automated reaction network generator is employed to create a database of structures. A new group additivity (GA) model based on principal component analysis (PCA) descriptors is introduced and trained on gas-phase density functional theory data of 4100 species at the M06-2X/6-311++G(d,p) level. Hydrogen bonds, local steric, and nonaromatic ring contributions are also incorporated. Lastly, we improve the accuracy of the group additivity model to reach the G4 theory by computing a data set of 770 species at this level and using a data fusion approach.

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

Molecular identification of wines using in situ liquid SIMS and PCA analysis

Composition analysis in wine is gaining increasing attention because it can provide information about the wine quality, source, and nutrition. In this work, in situ liquid secondary ion mass spectrometry (SIMS) was applied to 14 representative wines, including six wines manufactured by a manufacturer in Washington State, United States, four Cabernet Sauvignon wines, and four Chardonnay wines from other different manufacturers and locations. In situ liquid SIMS has the unique advantage of simultaneously examining both organic and inorganic compositions from liquid samples. Principal component analysis (PCA) of SIMS spectra showed that red and white wines can be clearly differentiated according to their aromatic and oxygen-contained organic species. Furthermore, the identities of different wines, especially the same variety of wines, can be enforced with a combination of both organic and inorganic species. Meanwhile, in situ liquid SIMS is sample-friendly, so liquid samples can be directly analyzed without any prior sample dilution or separation. Taken together, we demonstrate the great potential of in situ liquid SIMS in applications related to the molecular investigation of various liquid samples in food science.

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

Correlated conformational dynamics of the human GluN1-GluN2A type N-methyl-D-aspartate (NMDA) receptor

N-Methyl-D-aspartate receptors (NMDARs) are glutamate-gated ion channels found in the nerve cell membranes. As a result of overexcitation of NMDARs, neuronal death occurs and may lead to diseases such as epilepsy, stroke, Alzheimer’s disease, and Parkinson’s disease. In this study, human GluN1- GluN2A type NMDAR structure is modeled based on the X-ray structure of the Xenopus laevis template and missing loops are added by ab-initio loop modeling. The final structure is chosen according to two different model assessment scores. To be able to observe the structural changes upon ligand binding, glycine and glutamate molecules are docked into the corresponding binding sites of the receptor. Subsequently, molecular dynamics simulations of 1.3 μs are performed for both apo and ligand-bound structures. Structural parameters, which have been considered to show functionally important changes in previous NMDAR studies, are monitored as conformational rulers to understand the dynamics of the conformational changes. Moreover, principal component analysis (PCA) is performed for the equilibrated part of the simulations. Additionally, from these analyses, the differences in between apo and ligand-bound simulations can be summarized as the following: The girdle right at the beginning of the pore loop, which connects M2 and M3 helices of the ion channel, partially opens. Ligands act like an adhesive for the ligand-binding domain (LBD) by keeping the bi-lobed structure together and consequently this is reflected to the overall dynamics of the protein as an increased correlation of the LBD with especially the amino-terminal domain (ATD) of the protein.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Uncertainty Quantification for Smooth Functional Data with Application to Material Properties

This document outlines a method for processing functional output (i.e., curves) for the ultimate purpose of sampling curves under specified input conditions for use in modeling and simulation uncertainty quantification (UQ) studies. A set of benchmark curves sufficiently representative of the relevant scenario(s) being simulated are provided to the process and formatted as described in Section 1. Principal Component Analysis (PCA) is utilized to discover the components of uncertainty in the benchmark curves and is outlined in Section 2. Section 3 describes the application of uncertainty quantification to the PCA results for the purpose of sampling curves to be used in UQ analysis. Section 4 applies these techniques to an example benchmark dataset. Concluding remarks are provided in the final section.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Elemental and isotopic signatures of Asteroid Ryugu support three early Solar System reservoirs

Understanding the number and locations of different reservoirs present in the early Solar System is crucial to understanding the Solar System’s origin and evolution. Previous work has suggested that three unique isotopic reservoirs existed in the early Solar System but subsequent works have challenged that idea. Here we present elemental abundances along with Ca, Ti, Cr, Fe, Ni, and Zn isotopic data from primitive material returned by the Japan Aerospace Exploration Agency’s (JAXA) Hayabusa2 mission to asteroid (162173) Ryugu to make inferences on the Solar System’s early architecture. Data from Ryugu particle A0208 are consistent with a close genetic heritage between Ryugu and CI chondrites. Here, we employ principal component analysis (PCA) on these Ryugu and published meteorite data to demonstrate that Ryugu and CI chondrites are distinct from other known astromaterials, strongly supporting the existence of a third major isotopic reservoir in the early Solar System.

Isotopes↗
DOE OSTI2023

Machine-learning-assisted analysis of transition metal dichalcogenide thin-film growth

In situ reflective high-energy electron diffraction (RHEED) is widely used to monitor the surface crystalline state during thin-film growth by molecular beam epitaxy (MBE) and pulsed laser deposition. With the recent development of machine learning (ML), ML-assisted analysis of RHEED videos aids in interpreting the complete RHEED data of oxide thin films. The quantitative analysis of RHEED data allows us to characterize and categorize the growth modes step by step, and extract hidden knowledge of the epitaxial film growth process. In this study, we employed the ML-assisted RHEED analysis method to investigate the growth of 2D thin films of transition metal dichalcogenides (ReSe2) on graphene substrates by MBE. Principal component analysis (PCA) and K-means clustering were used to separate statistically important patterns and visualize the trend of pattern evolution without any notable loss of information. Using the modified PCA, we could monitor the diffraction intensity of solely the ReSe2 layers by filtering out the substrate contribution. These findings demonstrate that ML analysis can be successfully employed to examine and understand the film-growth dynamics of 2D materials. Further, the ML-based method can pave the way for the development of advanced real-time monitoring and autonomous material synthesis techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling MTS pyrolysis and SiC deposition kinetics using principal component analysis and neural networks

Accurate chemical kinetics modeling is crucial for improving the efficiency of chemical processing and synthesis of ceramic matrix composites. Detailed kinetic models are computationally expensive due to the large number of transported chemical species, while the simplified physics-based models, such as single-step global mechanisms, are efficient but often overlook key chemical intermediates and pathways. Recent deep learning approaches promise accurate and cost-effective models. Yet, they require additional closures for the transported nonlinear latent variables, complicating integration with existing solvers. In this work, we develop a hybrid linear—nonlinear reduced model for silicon carbide deposition from methyltrichlorosilane precursor by combining principal component analysis (PCA) and autoencoder (AE) neural network (NN) approaches. PCA is used to identify a smaller set of linear transport variables, enabling direct reuse of conventional transport solvers. NNs then reconstruct the full chemical state from these reduced variables. We demonstrate the method on a chemical vapor deposition reactor—comprising a gas-phase pyrolysis plug flow reactor and a heterogeneous surface reactor—over a wide range of temperatures, pressures, and residence times. Our PCA–AE model achieves high accuracy with only five transported scalars, achieving an eightfold cost reduction compared to detailed mechanisms, in both a priori (using data from the test set only) and a posteriori (coupled with a differential equation solver). In conclusion, notable errors arise primarily near training domain boundaries and for long residence times, indicating the need for domain shift indicators and better long-horizon predictions in future reduced chemistry model development.

autoencoder neural networks↗
DOE OSTI2025

Advancements in Constitutive Model Calibration: Leveraging the Power of Full‐Field DIC Measurements and In Situ Load Path Selection for Reliable Parameter Inference

Accurate material characterization and model calibration are essential for computationally supported high-consequence engineering decisions. Historically, characterization and calibration methods (1) use simplified test specimen geometries and global data, (2) cannot guarantee that sufficient characterization data are collected for a specific model of interest, (3) use deterministic methods that provide best-fit parameter values with no uncertainty quantification, and (4) are sequential, inflexible, and time-consuming. This work brings together several recent advancements into an improved workflow called interlaced characterization and calibration (ICC) that advances the state-of-the-art in constitutive model calibration. The ICC paradigm (1) employs tools to efficiently use full-field data to calibrate high-fidelity material models, (2) aligns the data needed with the data collected by adopting an optimal experimental design protocol, (3) quantifies parameter uncertainty through Bayesian inference and (4) incorporates these advancements into a quasi real-time feedback loop. The ICC framework is demonstrated here on the calibration of a material model using simulated full-field data for an aluminium cruciform specimen being deformed biaxially. The cruciform is actively driven through the myopically preferred load path using Bayesian optimal experimental design, which selects load steps that yield the maximum expected information gain (EIG). Principal component analysis (PCA) is performed on the model predictions of full-field displacements, and fast surrogate models are built to approximate the input-output relationships of the expensive finite element model. Furthermore, the tools developed and demonstrated here show that high-fidelity constitutive models can be efficiently and reliably calibrated with quantified uncertainty, thus supporting credible decision-making and potentially increasing the agility of solid mechanics modelling by enabling utilization of computational simulations at earlier stages of the design cycle.

Bayesian optimal experimental design↗
DOE OSTI2023

Rheological Properties of Small-Molecular Liquids at High Shear Strain Rates

Molecular-scale understanding of rheological properties of small-molecular liquids and polymers is critical to optimizing their performance in practical applications such as lubrication and hydraulic fracking. We combine nonequilibrium molecular dynamics simulations with two unsupervised machine learning methods: principal component analysis (PCA) and t-distributed stochastic neighbor embedding (t-SNE), to extract the correlation between the rheological properties and molecular structure of squalane sheared at high strain rates (10 6 –10 10 s -1 ) for which substantial shear thinning is observed under pressures P ϵ 0.1–955 MPa at 293 K. Intramolecular atom pair orientation tensors of 435 × 6 dimensions and the intermolecular atom pair orientation tensors of 61 × 6 dimensions are reduced and visualized using PCA and t-SNE to assess the changes in the orientation order during the shear thinning of squalane. Dimension reduction of intramolecular orientation tensors at low pressures P = 0.1,100 MPa reveals a strong correlation between changes in strain rate and the orientation of the side-backbone atom pairs, end-backbone atom pairs, short backbone-backbone atom pairs, and long backbone-backbone atom pairs associated with a squalane molecule. At high pressures P ≥ 400 MPa, the orientation tensors are better classified by these different pair types rather than strain rate, signaling an overall limited evolution of intramolecular orientation with changes in strain rate. Dimension reduction also finds no clear evidence of the link between shear thinning at high pressures and changes in the intermolecular orientation. The alignment of squalane molecules is found to be saturated over the entire range of rates during which squalane exhibits substantial shear thinning at high pressures.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Transient uncertainty quantification and Global Sensitivity Analysis of the open-source Molten Chloride Reactor Experiment (MCRE) using GP-PCA surrogate models

Uncertainties in the thermophysical properties of molten salts impact both the steady-state and transient behavior of Molten Salt Reactors (MSRs). In this work, we aim to quantify the influence of such uncertainties on the transient operation of the Molten Chloride Reactor Experiment (MCRE), utilizing the open-source specifications provided for this reactor. Seven representative transient scenarios are considered. For each scenario, we evaluate the impact of thermophysical property uncertainties on four key multiphysics model output variables of interest (VoIs): maximum power density, maximum fuel temperature, maximum reflector temperature, and average fuel velocity magnitude. In addition, we perform a Global Sensitivity Analysis (GSA) by computing Sobol’ indices for the uncertain input parameters to determine their contribution to the variability of each VoI. Conducting GSA is computationally intensive due to the large number of required evaluations of the high-fidelity multiphysics model. To mitigate this cost, we develop a surrogate modeling framework that combines Gaussian Process (GP) regression with Principal Component Analysis (PCA), enabling efficient sample generation for the GSA. Our results show that for energy-related VoIs, thermal conductivity is the dominant contributor to uncertainty. In contrast, for flow-related VoIs, density and dynamic viscosity are the primary sources of uncertainty. The specific heat of the fuel salt was found to play a secondary role in the transient analyses.

42 - ENGINEERING↗
DOE OSTI2021

Structural Polymorphism of Chitin and Chitosan in Fungal Cell Walls From Solid-State NMR and Principal Component Analysis

Chitin is a major carbohydrate component of the fungal cell wall and a promising target for novel antifungal agents. However, it is technically challenging to characterize the structure of this polymer in native cell walls. Here, we recorded and compared 13 C chemical shifts of chitin using isotopically enriched cells of six Aspergillus, Rhizopus, and Candida strains, with data interpretation assisted by principal component analysis (PCA) and linear discriminant analysis (LDA) methods. The structure of chitin is found to be intrinsically heterogeneous, with peak multiplicity detected in each sample and distinct fingerprints observed across fungal species. Fungal chitin exhibits partial similarity to the model structures of α- and γ-allomorphs; therefore, chitin structure is not significantly affected by interactions with other cell wall components. Addition of antifungal drugs and salts did not significantly perturb the chemical shifts, revealing the structural resistance of chitin to external stress. In addition, the structure of the deacetylated form, chitosan, was found to resemble a relaxed two-fold helix conformation. This study provides high-resolution information on the structure of chitin and chitosan in their cellular contexts. The method is applicable to the analysis of other complex carbohydrates and polymer composites.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Combining ToF‐SIMS and Multivariate Analysis to Resolve Active Sites on Ni‐Based HER Catalysts

Unambiguous identification of active sites in heterogeneous catalysis remains a major challenge, particularly for materials with ultrathin, chemically mixed surface layers. Here, we demonstrate a generalizable approach that combines time-of-flight secondary ion mass spectrometry (ToF-SIMS) with multivariate statistical analysis (principal component analysis [PCA] and multivariate curve resolution [MCR]) to resolve catalytically relevant motifs at the nanoscale. Using Ni electrodes as a model system, PCA distinguished hydroxide-enriched domains from oxide- and metal-rich regions, while MCR decomposed depth profiles and 3D images into hydroxide, oxide, and metallic layers with nanometer resolution. A unique secondary-ion fragment, NiO 3 H 3 − (m/z 108.94), emerged as a marker of hydroxide-rich environments and correlated with hydrogen evolution reaction (HER) activity across a series of Ni electrodes. Complementary density functional theory (DFT) calculations revealed that Ni(OH) 2 clusters adjacent to metallic Ni offer the most favorable water dissociation energetics, establishing the structural origin of the marker. While illustrated here for Ni-based HER, this workflow provides a broadly applicable framework to isolate and rank near-surface patterns that govern catalytic activity, thereby extending ToF-SIMS from a qualitative probe to a predictive tool for active site identification.

HER active sites↗
DOE OSTI2022

Decomposition of Irganox 1010 in plastic bonded explosives

Abstract Degradation pathways of Irganox 1010 in aged plastic bonded explosive (PBX) 9501 were investigated using ultrahigh performance liquid chromatography coupled to quadrupole time of flight mass spectrometry (UHPLC‐QTOF). Using a targeted approach, a total of 44 Irganox 1010 decomposition products were discovered. These decomposition products were formed through hydrolysis, scission, and/or oxidation of Irganox 1010. The hydrolytic decomposition of Irganox is a straightforward process resulting in the cleavage of the ester group(s) while oxidation and scission are more complicated and can happen at multiple locations on the Irganox 1010 molecule. Moreover, due to the symmetric nature of Irganox 1010, multiple decomposition reactions can occur. Indeed some decomposition products exhibited hydrolysis, oxidation, and scission. In order to probe any trends in the aged PBX 9501 samples, principal component analysis (PCA) was implemented. The greatest chemical differences between the aged PBX samples was hydrolysis of the ester functional groups on Irganox 1010. Despite the negative connotations of hydrolysis, the Irganox 1010 decomposition products are still able to function as a radical scavenger in PBX 9501 as intended.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗