Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix factorization”

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

Estimating Eigenenergies from Quantum Dynamics: A Unified Noise-Resilient Measurement-Driven Approach

Ground state energy estimation in physical, chemical, and materials sciences is one of the most promising applications of quantum computing. In this work, we introduce a new hybrid approach that finds the eigenenergies by collecting real-time measurements and post-processing them using the machinery of dynamic mode decomposition (DMD). From the perspective of quantum dynamics, we establish that our approach can be formally understood as a stable variational method on the function space of observables available from a quantum many-body system. We also provide strong theoretical and numerical evidence that our method converges rapidly even in the presence of a large degree of perturbative noise, and show that the method bears an isomorphism to robust matrix factorization methods developed independently across various scientific communities. Our numerical benchmarks on spin and molecular systems demonstrate an accelerated convergence and a favorable resource reduction over state-of-the-art algorithms. The DMD-centric strategy can systematically mitigate noise and stands out as a leading hybrid quantum-classical eigensolver.

Shen, Yizhi↗
DOE OSTI2024

Accuracy optimized neural networks do not effectively model optic flow tuning in brain area MSTd

Accuracy-optimized convolutional neural networks (CNNs) have emerged as highly effective models at predicting neural responses in brain areas along the primate ventral stream, but it is largely unknown whether they effectively model neurons in the complementary primate dorsal stream. We explored how well CNNs model the optic flow tuning properties of neurons in dorsal area MSTd and we compared our results with the Non-Negative Matrix Factorization (NNMF) model, which successfully models many tuning properties of MSTd neurons. To better understand the role of computational properties in the NNMF model that give rise to optic flow tuning that resembles that of MSTd neurons, we created additional CNN model variants that implement key NNMF constraints – non-negative weights and sparse coding of optic flow. While the CNNs and NNMF models both accurately estimate the observer's self-motion from purely translational or rotational optic flow, NNMF and the CNNs with nonnegative weights yield substantially less accurate estimates than the other CNNs when tested on more complex optic flow that combines observer translation and rotation. Despite its poor accuracy, NNMF gives rise to tuning properties that align more closely with those observed in primate MSTd than any of the accuracy-optimized CNNs. This work offers a step toward a deeper understanding of the computational properties and constraints that describe the optic flow tuning of primate area MSTd.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Machine Learning for Geothermal Resource Exploration in the Tularosa Basin, New Mexico

Geothermal energy is considered an essential renewable resource to generate flexible electricity. Geothermal resource assessments conducted by the U.S. Geological Survey showed that the southwestern basins in the U.S. have a significant geothermal potential for meeting domestic electricity demand. Within these southwestern basins, play fairway analysis (PFA), funded by the U.S. Department of Energy’s (DOE) Geothermal Technologies Office, identified that the Tularosa Basin in New Mexico has significant geothermal potential. This short communication paper presents a machine learning (ML) methodology for curating and analyzing the PFA data from the DOE’s geothermal data repository. The proposed approach to identify potential geothermal sites in the Tularosa Basin is based on an unsupervised ML method called non-negative matrix factorization with custom k-means clustering. This methodology is available in our open-source ML framework, GeoThermalCloud (GTC). Using this GTC framework, we discover prospective geothermal locations and find key parameters defining these prospects. Our ML analysis found that these prospects are consistent with the existing Tularosa Basin’s PFA studies. This instills confidence in our GTC framework to accelerate geothermal exploration and resource development, which is generally time-consuming.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Neural Network Approaches for Mobile Spectroscopic Gamma-Ray Source Detection

Artificial neural networks (ANNs) for performing spectroscopic gamma-ray source identification have been previously introduced, primarily for applications in controlled laboratory settings. To understand the utility of these methods in scenarios and environments more relevant to nuclear safety and security, this work examines the use of ANNs for mobile detection, which involves highly variable gamma-ray background, low signal-to-noise ratio measurements, and low false alarm rates. Simulated data from a 2” × 4” × 16” NaI(Tl) detector are used in this work for demonstrating these concepts, and the minimum detectable activity (MDA) is used as a performance metric in assessing model performance.In addition to examining simultaneous detection and identification, binary spectral anomaly detection using autoencoders is introduced in this work, and benchmarked using detection methods based on Non-negative Matrix Factorization (NMF) and Principal Component Analysis (PCA). On average, the autoencoder provides a 12% and 23% improvement over NMF- and PCA-based detection methods, respectively. Additionally, source identification using ANNs is extended to leverage temporal dynamics by means of recurrent neural networks, and these time-dependent models outperform their time-independent counterparts by 17% for the analysis examined here. The paper concludes with a discussion on tradeoffs between the ANN-based approaches and the benchmark methods examined here.

Bilton, Kyle J. (ORCID:0000000184553689)↗
DOE OSTI2024

New Measurements of the Lyα Forest Continuum and Effective Optical Depth with LyCAN and DESI Y1 Data

Abstract We present the Ly α Continuum Analysis Network (LyCAN), a convolutional neural network that predicts the unabsorbed quasar continuum within the rest-frame wavelength range of 1040–1600 Å based on the red side of the Ly α emission line (1216–1600 Å). We developed synthetic spectra based on a Gaussian mixture model representation of nonnegative matrix factorization (NMF) coefficients. These coefficients were derived from high-resolution, low-redshift ( z < 0.2) Hubble Space Telescope/Cosmic Origins Spectrograph (COS) quasar spectra. We supplemented this COS-based synthetic sample with an equal number of DESI Year 5 mock spectra. LyCAN performs extremely well on testing sets, achieving a median error in the forest region of 1.5% on the DESI mock sample, 2.0% on the COS-based synthetic sample, and 4.1% on the original COS spectra. LyCAN outperforms principal component analysis (PCA) and NMF-based prediction methods using the same training set by 40% or more. We predict the intrinsic continua of 83,635 DESI Year 1 spectra in the redshift range of 2.1 ≤ z ≤ 4.2 and perform an absolute measurement of the evolution of the effective optical depth. This is the largest sample employed to measure the optical depth evolution to date. We fit a power law of the form τ ( z ) = τ 0 ( 1 + z ) γ to our measurements and find τ 0 = (2.46 ± 0.14) × 10 −3 and γ = 3.62 ± 0.04. Our results show particular agreement with high-resolution, ground-based observations around z = 2, indicating that LyCAN is able to predict the quasar continuum in the forest region with only spectral information outside the forest.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Predictions of the glass transition temperature and viscosity of organic aerosols from volatility distributions

Abstract. Volatility and viscosity are important properties of organic aerosols (OA), affecting aerosol processes such as formation, evolution, and partitioning of OA. Volatility distributions of ambient OA particles have often been measured, while viscosity measurements are scarce. We have previously developed a method to estimate the glass transition temperature (Tg) of an organic compound containing carbon, hydrogen, and oxygen. Based on analysis of over 2400 organic compounds including oxygenated organic compounds, as well as nitrogen- and sulfur-containing organic compounds, we extend this method to include nitrogen- and sulfur-containing compounds based on elemental composition. In addition, parameterizations are developed to predict Tg as a function of volatility and the atomic oxygen-to-carbon ratio based on a negative correlation between Tg and volatility. This prediction method of Tg is applied to ambient observations of volatility distributions at 11 field sites. The predicted Tg values of OA under dry conditions vary mainly from 290 to 339 K and the predicted viscosities are consistent with the results of ambient particle-phase-state measurements in the southeastern US and the Amazonian rain forest. Reducing the uncertainties in measured volatility distributions would improve predictions of viscosity, especially at low relative humidity. We also predict the Tg of OA components identified via positive matrix factorization of aerosol mass spectrometer (AMS) data. The predicted viscosity of oxidized OA is consistent with previously reported viscosity of secondary organic aerosols (SOA) derived from α-pinene, toluene, isoprene epoxydiol (IEPOX), and diesel fuel. Comparison of the predicted viscosity based on the observed volatility distributions with the viscosity simulated by a chemical transport model implies that missing low volatility compounds in a global model can lead to underestimation of OA viscosity at some sites. The relation between volatility and viscosity can be applied in the molecular corridor or volatility basis set approaches to improve OA simulations in chemical transport models by consideration of effects of particle viscosity in OA formation and evolution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Source apportionment of soot particles and aqueous-phase processing of black carbon coatings in an urban environment

Abstract. The impacts of soot particles on climate and human health depend on the concentration of black carbon (BC) as well as the thickness and composition of the coating material, i.e., organic and inorganic compounds internally mixed with BC. In this study, the size-resolved chemical composition of BC-containing aerosol was measured using a high-resolution soot particle aerosol mass spectrometer (SP-AMS) during wintertime in Fresno, California, a location influenced by abundant combustion emissions and frequent fog events. Concurrently, particle optical properties were measured to investigate the BC absorption enhancement. Positive matrix factorization (PMF) analysis was performed on the SP-AMS mass spectral measurements to explore the sources of soot particles and the atmospheric processes affecting the properties of BC coatings. The analysis revealed that residential woodburning and traffic are the dominant sources of soot particles. Alongside primary soot particles originating from biomass burning (BBOABC) and vehicles (HOABC) two distinct types of processed BC-containing aerosol were identified: fog-related oxidized organic aerosol (FOOABC) and winter-background OOABC (WOOABC). Both types of OOABC showed evidence of having undergone aqueous processing, albeit with differences. The concentration of FOOABC was substantially elevated during fog events, indicating the formation of aqueous secondary organic aerosol (aqSOA) within fog droplets. On the other hand, WOOABC was present at a relatively consistent concentration throughout the winter and was likely related to the formation of secondary organic aerosol (SOA) in both the gas phase and aerosol liquid water. By comparing the chemical properties and temporal variations in FOOABC and WOOABC, we gain insights into the key aging processes of BC aerosol. It was found that aqueous-phase reactions facilitated by fog droplets had a significant impact on the thickness and chemical composition of BC coatings, thereby affecting the light absorption and hygroscopic properties of soot particles. These findings underscore the important role of chemical reactions occurring within clouds and fogs and influencing the climate forcing of BC aerosol in the atmosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Photoaging of phenolic secondary organic aerosol in the aqueous phase: evolution of chemical and optical properties and effects of oxidants

Abstract. While gas-phase reactions are well established to have significant impacts on the mass concentration, chemical composition, and optical properties of secondary organic aerosol (SOA), the aqueous-phase aging of SOA remains poorly understood. In this study, we performed a series of long-duration photochemical aging experiments to investigate the evolution of the composition and light absorption of the aqueous SOA (aqSOA) from guaiacyl acetone (GA), a semivolatile phenolic carbonyl that is common in biomass burning smoke. The aqSOA was produced from reactions of GA with hydroxyl radical (•OH-aqSOA) or a triplet excited state of organic carbon (3C∗-aqSOA) and was then photoaged in water under conditions that simulate sunlight exposure in northern California for up to 48 h. The effects of increasing aqueous-phase •OH or 3C∗ concentration on the photoaging of the aqSOA were also studied. High-resolution aerosol mass spectrometry (HR-AMS) and UV–Vis spectroscopy were utilized to characterize the composition and the light absorptivity of the aqSOA and to track their changes during aging. Compared to •OH-aqSOA, the 3C∗-aqSOA is produced more rapidly and shows less oxidation, a greater abundance of oligomers, and higher light absorption. Prolonged photoaging promotes fragmentation and the formation of more volatile and less light-absorbing products. More than half of the initial aqSOA mass is lost, and substantial photobleaching occurs after 10.5 h of prolonged aging under simulated sunlight illumination for 3C∗-aqSOA and 48 h for •OH-aqSOA. By performing positive matrix factorization (PMF) analysis of the combined HR-AMS and UV–Vis spectral data, we resolved three generations of aqSOA with distinctly different chemical and optical properties. The first-generation aqSOA shows significant oligomer formation and enhanced light absorption at 340–400 nm. The second-generation aqSOA is enriched in functionalized GA species and has the highest mass absorption coefficients in 300–500 nm, while the third-generation aqSOA contains more fragmented products and is the least light absorbing. These results suggest that intermediately aged phenolic aqSOA is more light absorbing than other generations, and that the light absorptivity of phenolic aqSOA results from a competition between brown carbon (BrC) formation and photobleaching, which is dependent on aging time. Although photoaging generally increases the oxidation of aqSOA, a slightly decreased O/C of the •OH-aqSOA is observed after 48 h of prolonged photoaging with additional •OH exposure. This is likely due to greater fragmentation and evaporation of highly oxidized compounds. Increased oxidant concentration accelerates the transformation of aqSOA and promotes the decay of BrC chromophores, leading to faster mass reduction and photobleaching. In addition, compared with •OH, photoaging by 3C∗ produces more low-volatility functionalized products, which counterbalances part of the aqSOA mass loss due to fragmentation and evaporation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Chemical properties and single-particle mixing state of soot aerosol in Houston during the TRACER campaign

Abstract. A high-resolution soot particle aerosol mass spectrometer (SP-AMS) was used to selectively measure refractory black carbon (rBC) and its associated coating material using both the ensemble size-resolved mass spectral mode and the event trigger single particle (ETSP) mode in Houston, Texas, in summer 2022. This study was conducted as part of the Department of Energy Atmospheric Radiation Measurement (ARM) program's TRacking Aerosol Convection interactions ExpeRiment (TRACER) field campaign. The study revealed an average (±1σ) rBC concentration of 103 ± 176 ng m−3. Additionally, the coatings on the BC particles were primarily composed of organics (59 %; 219 ± 260 ng m−3) and sulfate (26 %; 94 ± 55 ng m−3). Positive matrix factorization (PMF) analysis of the ensemble mass spectra of BC-containing particles resolved four distinct types of soot aerosol, including an oxidized organic aerosol (OOABC,PMF) factor associated with processed primary organic aerosol, an inorganic sulfate factor (SO4,BC,PMF), an oxidized rBC factor (O-BCPMF), and a mixed mineral dust–biomass burning aerosol factor with significant contribution from potassium (K-BBBC,PMF). Additionally, K-means clustering analysis of the single-particle mass spectra identified eight different clusters, including soot particles enriched in hydrocarbon-like organic aerosol (HOABC,ETSP), sulfate (SO4,BC,ETSP), two types of rBC, OOA (OOABC,ETSP), chloride (ClBC,ETSP), and nitrate (NO3,BC,ETSP). The single-particle measurements demonstrate substantial variation in BC coating thickness with coating-to-rBC mass ratios ranging from 0.1 to 100. The mixing state index (χ), which denotes the degree of homogeneity of the soot aerosol, varied from 4 % to 94 % with a median of 40 %, indicating that the aerosol population lies in between internal and external mixing but has large temporal and source type variability. In addition, a significant fraction of BC-containing particles, a majority enriched with oxidized organics and sulfate, exhibit sufficiently high κ values and diameters conducive to activation as cloud nuclei under atmospherically relevant supersaturation conditions. This finding bears significance in comprehending the aging processes of rBC-containing particles and their activation into cloud droplets. Our analysis highlights the complex nature of soot aerosol and underscores the need to comprehend its variability across different environments for accurate assessment of climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Observation-constrained kinetic modeling of isoprene SOA formation in the atmosphere

Isoprene has the largest global non-methane hydrocarbon emission, and the oxidation of isoprene plays a crucial role in the formation of secondary organic aerosol (SOA). Two primary processes are known to contribute to SOA formation from isoprene oxidation: (1) the reactive uptake of isoprene-derived epoxides on acidic or aqueous particle surfaces and (2) the absorptive gas–particle partitioning of low-volatility oxidation products. In this study, we developed a new multiphase condensed isoprene oxidation mechanism that includes these processes with key molecular intermediates and products. The new mechanism was applied to simulate isoprene gas-phase oxidation products and SOA formation from previously published chamber experiments under a variety of conditions and atmospheric observations during the Southern Oxidant and Aerosol Studies (SOAS) field campaign. Our results show that SOA formation from most of the chamber experiments is reasonably reproduced using our mechanism, except when the concentration ratios of initial nitric oxide to isoprene exceed ~ 2, the formed SOA is significantly underpredicted. The SOAS simulations also reasonably agree with the measurements regarding the diurnal pattern and concentrations of different product categories, while the total isoprene SOA remains underestimated. The molecular compositions of the modeled SOA indicate that multifunctional low-volatility products contribute to isoprene SOA more significantly than previously thought, with a median mass contribution of ~ 57 % to the total modeled isoprene SOA. However, this contribution is intricately intertwined with IEPOX-derived SOA (IEPOX: isoprene-derived epoxydiols), posing challenges for their differentiation using bulk aerosol composition analysis (e.g., the aerosol mass spectrometer with positive matrix factorization). Furthermore, the SOA from these pathways may vary greatly, mainly dependent on the volatility estimation and treatment of particle-phase processes (i.e., photolysis and hydrolysis). Our findings emphasize that the various pathways to produce these low-volatility species should be considered in models to more accurately predict isoprene SOA formation. The new condensed isoprene chemical mechanism can be further incorporated into regional-scale air quality models, such as the Community Multiscale Air Quality Modelling System (CMAQ), to assess isoprene SOA formation on a larger scale.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The sources and diurnal variations of submicron aerosols in a coastal–rural environment near Houston, US

Aerosol properties were characterized at a rural site southwest of Houston from May to September 2022 during the intensive operation periods (IOPs) of the Tracking Aerosol Convection Interactions ExpeRiment (TRACER). Backward trajectory analysis reveals three major air mass types: marine air mass from the Gulf, urban air mass influenced by urban emissions, and regional air mass. Marine aerosols typically show a bimodal size distribution and have the lowest particle number and mass concentrations of PM 1 (particulate matter with an aerodynamic diameter of less than 1 µm), while aerosols from air masses strongly influenced by urban emissions exhibit the highest concentrations. Organic aerosol (OA) accounts for more than 50 % of PM 1 for urban and regional air masses, whereas sulfate is comparable to OA in marine air masses. Positive matrix factorization (PMF) analysis of aerosol mass spectra identifies 6 OA factors: hydrocarbon-like OA (HOA), OA from the oxidation of monoterpenes (MT-SOA), OA from the reactive uptake of isoprene epoxydiols by acidic sulfate particles (isoprene-SOA), oxygenated OA arising from shipping emissions (shipping-OOA), and two oxygenated OA factors with high O : C ratios (OOA1 and OOA2). OOA2 has the highest O : C ratio and exhibits elevated mass concentration in the afternoon. Similar diurnal variation of highly oxidized OA factors was commonly observed in the Houston area during previous studies and attributed to the SOA formation by photochemistry and mixing from aloft. Here, using air mass backward trajectories and a 1-D box model, we show the diurnal trend of OOA2 mass concentration is instead driven by changes in air mass arriving at the rural site. The air mass changes are likely caused by the shift between land breezes and sea/bay breezes. Within the same air mass type (e.g., either urban or marine air mass), OOA2 mass concentration is largely independent of wind direction and shows essentially no diurnal variation, suggesting OOA2 is related to aged OA with minimal influence by local emissions. This study helps identify the major sources of OA in the Houston region and highlights the impacts of both atmospheric chemistry and meteorology on aerosol properties in the coastal–rural environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Efficient Parallelization of Irregular Applications on GPU Architectures

With the enlarging computation capacity of general Graphics Processing Units (GPUs), leveraging GPUs to accelerate parallel applications has become a critical topic in academia and industry. However, a wide range of irregular applications with the computation-/memory-intensive nature cannot easily achieve high GPU utilization. The challenges mainly involve the following aspects: first, data dependence leads to coarse-grained kernel and inefficient parallelism; second, heavy GPU memory usage may cause frequent memory evictions and extra overhead of I/O; third, specific computation patterns produce memory redundancies; last, workload balance and data reusability conjunctly benefit the overall performance, but there may exist a dynamic trade-off between them. Targeting these challenges, this dissertation proposes multiple optimizations to accelerate two real-world applications: many-body correlation functions to simulate nuclear physics in a large-scale scientific system; the other is the eALS-based matrix factorization recommendation system. To accelerate the calculations of many-body correlation functions, this dissertation presents three frameworks in GPU memory management and multi-GPU scheduling. Firstly, an optimized systematic GPU memory management framework, MemHC, utilizes a series of new memory reduction designs in GPU memory allocation, CPU/GPU communications, and GPU memory oversubscription. Secondly, an enhanced multi-GPU scheduling framework, MICCO, particularly by taking both data dimension (e.g., data reuse and data eviction) and computation dimension into account. MICCO designs a heuristic scheduling algorithm and a machine learning-based regression model to generate the optimal settings of a proposed new concept to manage the trade-off. Thirdly, a locality-aware multi-GPU scheduling framework. This scheduler leverages pipeline batch generation with a looking-ahead strategy by building local dependency graphs for memory transfer reduction and better data reuse, achieving up to 79.92% memory cost reduction and 1.67x speedup. To parallelize the eALS-based recommendation system, this dissertation proposes an efficient CPU/GPU heterogeneous recommendation system, HEALS. HEALS employs newly designed architecture-adaptive data formats to achieve load balance and good data locality on CPU and GPU. To mitigate the data dependence, HEALS presents a CPU/GPU collaboration model for both task parallelism and data parallelism with multiple kernel computation optimizations. In summary, this dissertation efficiently accelerates two typical irregular applications on GPUs by building four frameworks, including CPU/GPU collaboration, GPU memory management, and multi-GPU scheduling.

Wang, Qihan↗
DOE OSTI2021

Butterfly Factorization Via Randomized Matrix-Vector Multiplications

This paper presents an adaptive randomized algorithm for computing the butterfly factorization of an m × n matrix with m ≈ n provided that both the matrix and its transpose can be rapidly applied to arbitrary vectors. The resulting factorization is composed of O(log n) sparse factors, each containing O(n) nonzero entries. The factorization can be attained using O(n 3/2 log n) computation and O(n log n) memory resources. Furthermore, the proposed algorithm can be implemented in parallel and can apply to matrices with strong or weak admissibility conditions arising from surface integral equation solvers as well as multi-frontal-based finite-difference, finite-element, or finite-volume solvers. A distributed-memory parallel implementation of the algorithm demonstrates excellent scaling behavior.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Spectroscopic factors, overlaps, and isospin symmetry from an R -matrix point of view

Background: Spectroscopic factors, overlaps, and isospin symmetry are often used in conjunction with single-particle wave functions for the phenomenological analysis of nuclear structure and reactions. Many differing prescriptions for connecting these quantities to physically relevant asymptotic normalization constants or widths are available in the literature, but their relationship and degree of validity are not always clear. Purpose: This paper derives relationships among the above quantities of interest using well-defined methodology and starting assumptions. Method: R -matrix theory is used as the primary tool to interoperate among the quantities of interest to this work. Particular attention is paid to effects arising from beyond the nuclear surface, where isospin symmetry is strongly violated. Results: Relationships between the quantities of interest are derived. Example applications of these methods to mirror levels in nucleon + 12 C, nucleon + 16 O, and nucleon + 26 Al are presented. A new approach to multilevel mirror symmetry is derived and applied to the first three 2 + states of 18 O and 18 Ne. Conclusions: In this work, the relationship between the quantities of interest is clarified and certain procedures are recommended. It is found that the asymptotic normalization constant of the second 2 + state in 18 Ne deduced from the mirror state in 18 O is significantly larger than found in previous work. This finding has the effect of increasing the 17 F($\textit{p,γ}$) 18 Ne reaction rate in novae.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Materials structure–property factorization for identification of synergistic phase interactions in complex solar fuels photoanodes

Abstract Properties can be tailored by tuning composition in high-order composition spaces. For spaces with complex phase behavior, modeling the properties as a function of composition and phase distribution remains a formidable challenge. We present materials structure–property factorization (MSPF) as an approach to automate modeling of such data and identify synergistic phase interactions. MSPF is an interpretable machine learning algorithm that couples phase mapping via Deep Reasoning Networks (DRNets) to matrix factorization-based modeling of the representative properties of each phase in a dataset. MSPF is demonstrated for Bi–Cu–V oxide photoanodes for solar fuel generation, which contains 25 different phase combinations and correspondingly exhibits complex composition-structure-photoactivity relationships. Comparing the measured photoactivity to a learned model for non-interacting phases, synergistic phase interactions are identified to guide further photoactivity optimization and understanding. MSPF identifies synergistic interactions of a BiVO 4 -like phase with both Cu 2 V 2 O 7 -like and CuV 2 O 6 -like phases, creating avenues for understanding complex photoelectrocatalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Application of Atmospheric Gases and Particulate Matter to the Assessment of Urban Heat Island

Background: Urban heat island (UHI), where built areas are warmer compared to non-urban regions, increases human related diseases and mortality. A key challenge in UHI analysis is the designation of sites as urban or suburban/rural; however, the growing complexity of green spaces in urban areas and the predominance of the transportation sector in nonurban areas creates a dilemma for distinct delineation. Objectives: This study aims to utilize the variability of atmospheric components such as particulate matter (PM), inorganic gases, and volatile organic compounds (VOCs) as direct tracers of the degree of urbanization for ground-based measurements to fully comprehend UHI in convoluted regions with indistinct delineation of urban and nonurban environments. Methods: Atmospheric gases and aerosols were used as direct tracers of urbanization for UHI analysis. Inorganic gases and particulate matter were monitored in two sites in a southeastern US city with varying degrees of urbanization. VOCs were analyzed using a proton transfer reaction time-of-flight mass spectrometer. Results: The more-urbanized site exhibited warmer night conditions and elevated total oxidant levels, leading to the formation of nanometer-sized particles. Machine learning analysis revealed similar atmospheric pollutant profiles for both sites, suggesting comparable sources and variability. Biogenic VOCs were enhanced at the less-urbanized site; however, levels of anthropogenic aromatic VOCs were comparable for both sites. A comprehensive mass spectra analysis revealed distinct molecular backbones per site that further affirmed the applicability of VOCs as indicators of urbanization. Conclusion: This study concludes that VOCs provide more direct and accurate information than typical inorganic gases and PM parameters for characterizing the degree of urbanization. Further exploration of VOCs can enhance our understanding of UHI dynamics and its interaction with vegetation in urban green spaces.

Air quality sensor↗
DOE OSTI2023

High-dimensional data analytics in civil engineering: A review on matrix and tensor decomposition

Recent developments in sensing and monitoring techniques have led to the generation of high-dimensional data in the field of civil engineering. High-dimensional data analytics methods have thus been developed to interpret such complex data. Among the different high-dimensional data analytics techniques, matrix and tensor decomposition methods have acquired a notable interest in the civil engineering community over the past decade. Due to their unique ability to deal with highly redundant and correlated data, these methods are establishing themselves as promising and efficient tools to analyze high-dimensional data in the civil engineering arena. In this paper, high-dimensional data is referred to as a data set in which the number of features is comparable or larger than the number of observations. This review paper aims to summarize the applications of matrix and tensor decomposition methods in civil engineering over the last decade. The survey begins with a general overview of matrix and tensor decomposition followed by highlighting their significance in the field. Afterward, various applications of these high-dimensional data analytics methods in civil engineering are presented, while the advantages offered by these methods are discussed. Lastly, challenges and potential research avenues for employing matrix and tensor decomposition and future emerging trends for their novel use are highlighted.

42 ENGINEERING↗
DOE OSTI2023

Analysis of gas storage and transport in Eagle Ford shale using pressure pulse decay measurements with He, Kr and CO 2

Detailed characterization of gas transport and storage in low-permeability and organic-rich shales is associated with an array of challenges due to the complex morphology of the pore space, representing a broad range of pore sizes, combined with the heterogeneous fabric of the shale matrix. These factors, and their interplay during gas transport and sorption, complicate a) the analysis of shale samples at the laboratory scale (~ up to a ft. In length), and b) the prediction of natural gas production and carbon sequestration potential at larger scales. Here, in this work, tandem experiments with inert (helium – He) and adsorbing (krypton - Kr and carbon dioxide - CO 2 ) gases were performed and analyzed to develop an efficient workflow for characterizing and modeling transport and sorption at the laboratory scale. In particular, pressure pulse-decay (PPD) measurements were conducted on an Eagle Ford shale core sample at room temperature using He, Kr, and CO 2 . PPD measurements with He (a non-sorbing gas) were employed to probe the overall porosity, including natural fractures, microcracks, mesopores, and micropores. A triple-porosity model (TPM) was adopted to interpret the gas transport in the shale sample: The pore space is represented by three interacting continua, including macropores (larger fractures), mesopores (including microcracks), and micropores. To facilitate the application of the TPM, a modified analytical approach is introduced to extract effective transport parameters (in terms of the characteristic time for transport at relevant porosity levels) directly from the PPD measurements with He. To validate the modeling approach, model parameters extracted from two He PPD measurements are demonstrated to provide excellent agreement with a 3rd He PPD measurement performed at a higher pressure. The effective transport parameters, extracted from the He experiments, were subsequently converted for application to the Kr and CO 2 PPD experiments, by accounting for relevant transport modes and differences in fluid properties. Excess adsorption isotherms were extracted from the equilibrium pressures of the Kr and CO 2 PPD experiments using the He pore volume as a baseline. These adsorption isotherms were then integrated into the TPM to predict combined gas transport and sorption for Kr and CO 2 , with transport coefficients translated from the He measurements. The predictions for Kr and CO 2 are demonstrated to be in excellent agreement with the experimental observations. This, in turn, demonstrates that the proposed analytical approach provides for an effective characterization of mass transfer rates in shales, that can be applied directly in a TPM representation of mass transfer and sorption.

58 GEOSCIENCES↗