Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random Effects”

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

Using iterative random forest to find geospatial environmental and Sociodemographic predictors of suicide attempts

Despite a recent global decrease in suicide rates, death by suicide has increased in the United States. It is therefore imperative to identify the risk factors associated with suicide attempts to combat this growing epidemic. In this study, we aim to identify potential risk factors of suicide attempt using geospatial features in an Artificial intelligence framework. We use iterative Random Forest, an explainable artificial intelligence method, to predict suicide attempts using data from the Million Veteran Program. This cohort incorporated 405,540 patients with 391,409 controls and 14,131 attempts. Our predictive model incorporates multiple climatic features at ZIP-code-level geospatial resolution. We additionally consider demographic features from the American Community Survey as well as the number of firearms and alcohol vendors per 10,000 people to assess the contributions of proximal environment, access to means, and restraint decrease to suicide attempts. In total 1,784 features were included in the predictive model. Our results show that geographic areas with higher concentrations of married males living with spouses are predictive of lower rates of suicide attempts, whereas geographic areas where males are more likely to live alone and to rent housing are predictive of higher rates of suicide attempts. We also identified climatic features that were associated with suicide attempt risk by age group. Additionally, we observed that firearms and alcohol vendors were associated with increased risk for suicide attempts irrespective of the age group examined, but that their effects were small in comparison to the top features. Taken together, our findings highlight the importance of social determinants and environmental factors in understanding suicide risk among veterans.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Photovoltage behaviour of p-Sb 2 S 3 photocathodes for hydrogen evolution: effect of n-In 2 S 3 passivation layers

The 1.76 eV band gap of antimony(iii) sulphide (Sb 2 S 3 ) makes this semiconductor material a promising light absorber for photoelectrochemical water splitting, but scalable fabrication approaches to efficient devices are still lacking. Here we show that compact Sb 2 S 3 films on FTO can be obtained by electrochemical growth from aqueous colloidal sulphur and antimony trichloride solutions, followed by mild annealing. These films can be converted into hydrogen evolution photocathodes after coating with In 2 S 3 passivation layers and the addition of Pt proton reduction co-catalysts. For the first time, vibrating Kelvin probe surface photovoltage (VKP-SPV) spectroscopy is used to observe the carrier dynamics in such photoelectrodes. While the bare Sb 2 S 3 films suffer from high surface recombination rates and poor electron extraction, the In 2 S 3 overlayer is found to raise the photovoltage and cathodic photocurrent density, due to passivation of surface defects and formation of a p–n heterojunction. In thick In 2 S 3 films, these benefits are offset by shading and slow electron transfer. Also, we find that O 2 strongly affects the band bending in the Sb 2 S 3 –air and In 2 S 3 –air junctions and their photovoltage. The optimised devices evolve H 2 at 77.5% Faradaic efficiency and with 0.084% applied bias photon-to-current efficiency (ABPE) at 0.12 V vs. RHE. The low ABPE value is attributed to Sb 2 S 3 sub-bandgap defects visible in SPV spectra, the random orientation of Sb 2 S 3 crystallites in the films, which inhibits charge transport, the absence of crystal facets of Sb 2 S 3 , and a detrimental Schottky junction at the FTO|Sb 2 S 3 interface.

de Araújo, Moisés A. [University of California, Da↗
DOE OSTI2023

Community response of soil microorganisms to combined contamination of polycyclic aromatic hydrocarbons and potentially toxic elements in a typical coking plant

Both polycyclic aromatic hydrocarbons (PAHs) and potentially toxic elements (PTEs) of coking industries impose negative effects on the stability of soil ecosystem. Soil microbes are regarded as an essential moderator of biochemical processes and soil remediation, while their responses to PAHs-PTEs combined contamination are largely unknown. In the present study, soil microbial diversity and community composition in the typical coking plant under the chronic co-exposure of PAHs and PTEs were investigated and microbial interaction networks were built to reveal microbial co-occurrence patterns. The results indicated that the concentrations of PAHs in the soil inside the coking plant were significantly higher than those outside the plant. The mean concentration of ∑16PAHs was 2894.4 ng·g −1 , which is 5.58 times higher than that outside the plant. The average Hg concentration inside the coking plant was 22 times higher than the background value of Hebei province. The soil fungal community inside the coking plant showed lower richness compared with that of outside community, and there are significant difference in the bacterial and fungal community composition between inside and outside of coking plant ( p < 0.01). Predicted contribution of different environmental factors to each dominant species based on random forest identified 20 and 25 biomarkers in bacteria and fungi, respectively, that were highly sensitive to coking plant soil in operation, such as Betaproteobacteria，Sordariomycetes and Dothideomycetes . Bacterial and fungal communities were shaped by the soil chemical properties (pH), PTEs (Hg), and PAHs together in the coking plant soils. Furthermore, the bacterial and fungal interaction patterns were investigated separately or jointly by intradomain and interdomain networks. Competition is the main strategy based on the co-exclusion pattern in fungal community, and the competitive relationship inside the coking plant is more complex than that outside the plant. In contrast, cooperation is the dominant strategy in bacterial networks based on the co-occurrence pattern. The present study provided insights into microbial response strategies and the interactions between bacteria and fungi under long-term combined contamination.

Shen, Qihui↗
DOE OSTI2024

Machine Learning-Assisted Stability Boundary Determination of Multiport Autonomous Reconfigurable Solar Power Plants

The multiport autonomous reconfigurable solar (MARS) power plant is a promising solution to integrate renewable resources and energy storage systems into the alternating current (ac) power grid and an high-voltage direct current (HVdc) link. In the MARS system, various input power sources are connected to the individual submodules (SMs) through direct current (dc)–dc converters. However, the presence of external power sources can result in unbalanced capacitor voltages of SMs, thereby violating stability constraints under multiple/diverse operating conditions. This article aims to address the gap by accurately determining the stability boundary of the MARS system. As such, a novel machine learning (ML)-assisted energy balancing control (EBC) criterion is proposed. Further, in conjunction with a refined EBC, this approach ensures balanced capacitor voltages across various types of SMs, significantly enhancing the overall system efficiency. The proposed EBC criterion effectively controls EBC activation and deactivation, achieving remarkable accuracy. Both power systems computer aided design (PSCAD)/electromagnetic transients including direct current (EMTDC) simulations and control hardware-in-the-loop (cHIL) tests are conducted to validate the feasibility and efficiency of the proposed method. By combining the EBC and ML-assisted EBC criterion, efficient energy management is achieved for systems featuring multiple input power sources, such as MARS. This approach enables the system to fully exploit its potential across an expanded operational range while upholding high-efficiency standards.

14 SOLAR ENERGY↗
DOE OSTI2022

Enhancing scalability of a matrix-free eigensolver for studying many-body localization

We propose several techniques to enhance the parallel scalability of a matrix-free eigensolver designed for studying many-body localization (MBL) of quantum spin chain models with nearest-neighbor interactions and on-site disorder. This type of problem is computationally challenging because the dimension of the associated Hamiltonian matrix grows exponentially with respect to the number of spins L, and we need to average over different realizations of the random disorder to obtain relevant statistical behavior. For each disorder realization, we need to compute eigenvalues from different regions of the spectrum and their corresponding eigenvectors. In previous work, the interior eigenstates for a single eigenvalue problem are computed via the shift-and-invert Lanczos algorithm. Due to the extremely high memory footprint of the LU factorizations, this technique is not well suited for large L’s. For example, we need thousands of compute nodes on modern high performance computing infrastructures to go beyond L = 24. The matrix-free approach does not suffer from this memory bottleneck, however, its scalability is limited by a computation and communication load imbalance. To reduce this imbalance and to significantly enhance the scalability of the matrix-free eigensolver, we reorder the matrix and leverage the consistent space runtime, CSPACER. We also show its efficiency in managing irregular communication patterns at scale compared to optimized MPI non-blocking two-sided and one-sided RMA implementation variants. This effort enables us to study MBL for spin chains with a larger number of spins. The efficiency and effectiveness of the proposed algorithm is demonstrated by computing eigenstates on a massively parallel many-core high performance computer.

METIS↗
DOE OSTI2023

Gamma Spectrometry Code Rodeo for Uranium Enrichment—FY 2022 Report

In the first two quarters of FY22, data acquisition continued at ORNL and LLNL using uranium sources of known enrichments. This was an FY21 task which could not be completed in FY21 because of problems encountered with the ORNL M400 CZT in Q4 of FY21, and the subsequent repairs. The detector was received back from H3D in the first of September 2021 , and the measurements resumed . Measurements using the repaired detector were completed in Q1 of FY22. The spectra were distributed by ORNL to the analyzing labs. Analysis results were received in Q2 of FY2022. The results from the various codes were intercompared and an ANOVA analysis was performed. Random and systematic uncertainties were established for each code. The ANOVA results and discussions were included in a revised version of FY21Annual Report issued in March 2022. A paper was presented at the INMM 2022Annual Conference, with the analysis results from the various isotopic codes, and the ANOVA table with random and systematic uncertainties. The Project Work Plan (PWP) for FY22 included a task to perform field testing of the M400 CZT and the analysis codes using UF 6 cylinder measurements at the Framatome Fuel Fabrication Facility in Richland, WA. PNNL was the lead for the field testing task. PNNL drafted a Field Test Plan, and refined it based on comments received from the team. PNNL coordinated with Framatome facility, the logistics of carrying out the field testing . A collimator and shield made out of TFlex (tungsten impregnated polymer) was designed and professionally manufactured. The collimators were used in the field test measurements. The measurements at Framatome were completed on April 21, 2022. A total of 34 Type 30B cylinders were measured using three M400 detectors (PNNL, LLNL, and ORNL detectors). Measurements using M400 were taken at two locations on the side of each cylinder and from the end-on bottom location. Additionally, HPGe measurements were taken at the end-on location to establish ground truth. Cylinder wall thickness measurements were also made at all three locations. To gain a better understanding of the effect of background from surrounding cylinders, the same five cylinders measured individually in low background locations were measured again in the cylinder storage yards. Due to inclement weather, manufacturer delays, shipping delays, and equipment failure, the measurement campaign spanned twice as long compared to the original timeline. Gamma-ray spectra from M400 and HPGe detectors, along with the cylinder data and photographs were organized and shared with the collaborators for further analysis. Spectra were analyzed by the participating laboratories. FY22 PWP also consists of tasks related to plutonium source measurements, adapting the codes to analyze plutonium spectra, and inter-comparison of the results from various codes. Plutonium spectra are being acquired at LANL, ORNL, and LLNL. LANL, SNL and LLNL are in the process of modifying FRAM, GADRAS, and CZTU, respectively. The plutonium related tasks will be completed in Q2 of FY23.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

The Impacts of Single-Scattering and Microphysical Properties of Ice Particles Smaller Than 100 µm on the Bulk Radiative Properties of Tropical Cirrus

There are large uncertainties in the single-scattering (i.e., morphologies) and microphysical (i.e., concentrations) properties of ice particles whose size are less than ~100 µm. Insufficient resolutions of the most advanced cloud probes (e.g., cloud particle imager) cannot resolve the micrometer-scale morphologies of small ice particles. Further, the shattering of large ice particles on probes’ inlets or tips causes uncertainties in the measurement of the concentrations of small ice particles. These uncertainties have large impacts on the single-scattering and microphysical properties of small ice particles that are utilized to quantify the bulk radiative properties of cirrus. In this study, the impacts of uncertainties in the morphologies and concentrations of small ice particles on the bulk radiative properties of tropical cirrus were calculated using measurements acquired during the Tropical Warm Pool-International Cloud Experiment. Five different models (i.e., budding Buckyball, Chebyshev particle, droxtal, Gaussian random sphere, and sphere) that represent the shapes of small ice particles were used to calculate the single-scattering properties. The bulk radiative properties, average phase-function ($\overline{P_{11}}$), and average asymmetry parameter ($\overline{g}$) were computed by combining the measured size/habit distributions and the calculated single-scattering properties of ice particles. The impacts of the selection of varying morphologies of small particles on the bulk radiative properties were quantified. For these calculations, the possible range of the concentrations of small ice particles which depend on the degree of shattered large particles were also used. The impacts of varying the single-scattering properties of small ice particles on the bulk radiative properties were the largest in the upper parts of cirrus (T < –60 °C), while they were the smallest in the lower parts of cirrus (–45 < T < –30 °C). The impacts of uncertainties in the concentrations of small ice particles on the bulk radiative properties were largest in the lower parts of cirrus (–45 < T < –30 °C), whereas they were smallest in the upper parts of cirrus (T < –60 °C). The effect of shattering was maximum in the lower parts of cirrus, whilst it was minimum in the upper parts of cirrus. The combined impacts of uncertainties in the single-scattering (i.e., morphologies) and microphysical (i.e., concentrations) properties of small ice particles revealed variations of up to 11.2% (127.1%; 67.3%) of the integrated intensity in the forward (sideward; backward) angles in $\overline{P_{11}}$ and a corresponding change in $\overline{g}$ by up to 12.61%.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Data trustworthiness signatures for nuclear reactor dynamics simulation

With the increased reliance on digitization in industrial control systems, the need for effective monitoring techniques has risen dramatically. Specifically, there is now a growing concern about the so-called false data injection (FDI) attacks. These attacks aim to alter the raw sensors’ data to cause malicious outcomes. Any serious FDI algorithm is based on an intimate knowledge of the system and its associated physics models, which renders conventional outlier/anomaly detection techniques almost obsolete in the face of such attacks. Thus, a critical need has emerged to develop a new class of defense methods that are capable of detecting FDI attacks under the assumption that the attacker has a strong familiarity with the system and its physics modeling. This class of defense methods are denoted by model-based defenses which are premised on the assumption that the attacker, while having a good understanding of the system, does not have full privileged access to all proprietary data and historical records of operation. However, (s)he is assumed to be capable of learning system behavior using self-learning techniques during an initial lie-in-wait period. To defend against this scenario, we propose a new model-based randomized window algorithm that searches time-series data for signatures that can serve as classifiers between normal and FDI scenarios. The classifiers are based on the correlations between the dominant degrees of freedom (DOFs) and the less-dominant DOFs (expected to be very sensitive to the system details that are unknown to the attacker). For demonstration, RELAP5 models are employed to calculate representative nuclear reactor behavior during a number of transient scenarios. Finally, falsified data are injected into the RELAP5-simulated behavior, and the proposed signature-identification algorithm is employed to detect the injected data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Improved quality metrics for association and reproducibility in chromatin accessibility data using mutual information

Correlation metrics are widely utilized in genomics analysis and often implemented with little regard to assumptions of normality, homoscedasticity, and independence of values. This is especially true when comparing values between replicated sequencing experiments that probe chromatin accessibility, such as assays for transposase-accessible chromatin via sequencing (ATAC-seq). Such data can possess several regions across the human genome with little to no sequencing depth and are thus non-normal with a large portion of zero values. Despite distributed use in the epigenomics field, few studies have evaluated and benchmarked how correlation and association statistics behave across ATAC-seq experiments with known differences or the effects of removing specific outliers from the data. Here, we developed a computational simulation of ATAC-seq data to elucidate the behavior of correlation statistics and to compare their accuracy under set conditions of reproducibility. Using these simulations, we monitored the behavior of several correlation statistics, including the Pearson’s R and Spearman’s ρ coefficients as well as Kendall’s τ and Top–Down correlation. We also test the behavior of association measures, including the coefficient of determination R 2 , Kendall’s W, and normalized mutual information. Our experiments reveal an insensitivity of most statistics, including Spear man’s ρ, Kendall’s τ, and Kendall’s W, to increasing differences between simulated ATAC-seq replicates. The removal of co-zeros (regions lacking mapped sequenced reads) between simulated experiments greatly improves the estimates of correlation and association. After removing co-zeros, the R 2 coefficient and normalized mutual information display the best performance, having a closer one-to-one relationship with the known portion of shared, enhanced loci between simulated replicates. When comparing values between experimental ATAC-seq data using a random forest model, mutual information best predicts ATAC-seq replicate relationships. Collectively, this study demonstrates how measures of correlation and association can behave in epigenomics experiments. We provide improved strategies for quantifying relationships in these increasingly prevalent and important chromatin accessibility assays.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Dual-unitary shadow tomography

We introduce a classical shadow tomography scheme based on dual-unitary brick-wall circuits termed "dual-unitary shadow tomography" (DUST). For this we study operator spreading and Pauli weight dynamics in one-dimensional qubit systems, evolved by random two-local dual-unitary gates arranged in a brick-wall structure, ending with a final measurement layer. We do this by deriving general constraints on the Pauli weight transfer matrix and specializing to the case of dual-unitarity. We first show that dual-unitaries must have a minimal amount of entropy production. Remarkably, we find that operator spreading in these circuits have a rich structure resembling that of relativistic quantum field theories, with massless chiral excitations that can decay or fuse into each other, which we call left- or right-movers. We develop a mean-field description of the Pauli weight in terms of $\rho(x,t)$, which represents the probability of having nontrivial support at site $x$ and depth $t$ starting from a fixed weight distribution. We develop an equation of state for $\rho(x,t)$, and simulate it numerically using Monte Carlo simulations. Lastly, we demonstrate that the fast-thermalizing properties of dual-unitary circuits make them better at predicting large operators than shallow brick-wall Clifford circuits. Our results are robust to finite-size effects due to the chirality of dual-unitary brick-wall circuits.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The Sixth Data Release of the Radial Velocity Experiment (Rave). II. Stellar Atmospheric Parameters, Chemical Abundances, and Distances

We present part 2 of the sixth and final Data Release (DR6) of the Radial Velocity Experiment (Rave), a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in Earth’s southern hemisphere. The Rave medium-resolution spectra (R ∼ 7500) cover the Ca triplet region (8410–8795 Å) and span the complete time frame from the start of Rave observations on 2003 April 12 to their completion on 2013 April 4. In the second of two publications, we present the data products derived from 518,387 observations of 451,783 unique stars using a suite of advanced reduction pipelines focusing on stellar atmospheric parameters, in particular purely spectroscopically derived stellar atmospheric parameters (T{sub eff}, logg, and the overall metallicity), enhanced stellar atmospheric parameters inferred via a Bayesian pipeline using Gaia DR2 astrometric priors, and asteroseismically calibrated stellar atmospheric parameters for giant stars based on asteroseismic observations for 699 K2 stars. In addition, we provide abundances of the elements Fe, Al, and Ni, as well as an overall [α/Fe] ratio obtained using a new pipeline based on the GAUGUIN optimization method that is able to deal with variable signal-to-noise ratios. The Rave DR6 catalogs are cross-matched with relevant astrometric and photometric catalogs, and are complemented by orbital parameters and effective temperatures based on the infrared flux method. The data can be accessed via the Rave website (http://rave-survey.org) or the Vizier database.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Assessing the design of integrated methane sensing networks

Abstract While methane is the second largest contributor to global warming after carbon dioxide, it has a larger warming effect over a much shorter lifetime. Despite accelerated technological efforts to radically reduce global carbon dioxide emissions, rapid reductions in methane emissions are needed to limit near-term warming. Being primarily emitted as a byproduct from agricultural activities and energy extraction, methane is currently monitored via bottom–up (i.e. activity level) or top–down (via airborne or satellite retrievals) approaches. However, significant methane leaks remain undetected and emission rates are challenging to characterize with current monitoring frameworks. In this paper, we study the design of a layered monitoring approach that combines bottom–up and top–down approaches as an integrated sensing network. By recognizing that varying meteorological conditions and emission rates impact the efficacy of bottom–up monitoring, we develop a probabilistic approach to optimal sensor placement in its bottom–up network. Subsequently, we derive an inverse Bayesian framework to quantify the improvement that a design-optimized integrated framework has on emission-rate quantifications and their uncertainties. We find that under realistic meteorological conditions, the overall error in estimating the true emission rates is approximately 1.3 times higher, with their uncertainties being approximately 2.4 times higher, when using a randomized network over an optimized network, highlighting the importance of optimizing the design of integrated methane sensing networks. Further, we find that optimized networks can improve scenario coverage fractions by more than a factor of 2 over experimentally-studied networks, and identify a budget threshold beyond which the rate of optimized-network coverage improvement exhibits diminishing returns, suggesting that strategic sensor placement is also crucial for maximizing network efficiency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Simulations of axion minihalos

The axion, motivated as a solution to the strong CP problem, is also a viable dark matter candidate. The axion field takes random values in causally disconnected regions if the symmetry breaking that establishes the particle occurs after inflation, leading to white-noise density fluctuations at low wave numbers and forming dense minihalos with sub-planetary masses subsequently. There have been two recent proposals that appear capable of testing this scenario, namely using pulsar timing arrays and studying cosmological microlensing caustics. Motivated by these proposals, we use N-body simulations to study the formation of substructures from white-noise density fluctuations. The density profiles of our relaxed axion minihalos can be described by the Navarro–Frenk–White profile, and the minihalos’ concentration number agrees well with a simple, physically-motivated model. We develop a semi-analytic formula to fit the mass function from our simulation, which agrees broadly at different redshifts and only differs at factor of two level from classic halo mass functions. This analytic mass function allows us to consider uncertainties in the post-inflation axion scenario, as well as extrapolate our high-redshift simulations results to the present. Our work estimates the present-day abundance of axion substructures, as is necessary for predicting their effect on cosmological microlensing caustics and pulsar timing. Our calculations suggest that if pulsar timing and microlensing probes can reach recent sensitivity forecasts, they may be sensitive to the post-inflation axion dark matter scenario, even when accounting for uncertainties pertaining to axion strings. For pulsar timing, the most significant caveat is whether axion minihalos are disrupted by stars, which our estimates show is mildly important at the most relevant masses. In conclusion, as our gravitational simulations are scale invariant, the results can be extended to models where the dark matter is comprised of other axion-like particles and even clusters of primordial black holes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Gas generation and intramolecular isotope study in laboratory pyrolysis of the Springfield coal from the Illinois Basin

Position-specific (PS) isotopes of propane have been proposed as a potential geochemical tool to decipher various geological processes (e.g., thermal cracking, biodegradation, H exchange) in natural reservoirs. The limited studies have been conducted on the PS isotopes of propane from the pyrolysis gases from marine shales, and natural gases sourced from lacustrine and marine kerogens, but little is known on gases produced from the humic kerogen. Here, this study investigated the PS δ 13 C of propane in the closed-system pyrolysis of the Springfield coal, Illinois Basin, Indiana, at 310 to 470 °C (Easy %Ro: 0.76 to 3.07). The C kinetic isotope effect (KIE) of CH 4 produced in both this study and previous low-temperature pyrolysis of the same coal indicates the cleavage of C—O bonds is the main generation pathway at the early kerogen cracking stage, followed by the breakdown of alkyl groups. At the wet-gas cracking stage, C 3 H 8 production from thermally stable compounds has a significant influence on the bulk and position-specific C KIE in the pyrolysis of marine Woodford kerogen and Springfield coal. According to the PS δ 13 C of propane, the central site is likely more enriched in 13 C and the δ 13 C of the terminal site is relatively heterogeneous within the propyl group attached to different functional groups of the gas-prone kerogen. Our findings based on the pyrolysis experiments and natural gas samples indicate thermal cracking and biodegradation appear to alter the δ 13 C cen values more significantly than the δ 13 C ter values of propane. The larger magnitude of ΔC c-t in the kMC simulations (Peterson et al., 2018) compared with those from the marine shale and coal possibly implies the non-random distributions of 13 C of propane precursors in the kerogens. As a new dimension of intramolecular isotopic information of propane, the PS δ 13 C values can contribute to fingerprinting the gas origins and identifying the various geological processes (e.g., kerogen cracking, wet-gas cracking, microbial activities) in sedimentary basins.

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

Computational Imaging for Intelligence in Highly Scattering Aerosols (Final Report)

Natural and man-made degraded visual environments pose major threats to national security. The random scattering and absorption of light by tiny particles suspended in the air reduces situational awareness and causes unacceptable down-time for critical systems and operations. To improve the situation, we have developed several approaches to interpret the information contained within scattered light to enhance sensing and imaging in scattering media. These approaches were tested at the Sandia National Laboratory Fog Chamber facility and with tabletop fog chambers. Computationally efficient light transport models were developed and leveraged for computational sensing. The models are based on a weak angular dependence approximation to the Boltzmann or radiative transfer equation that appears to be applicable in both the moderate and highly scattering regimes. After the new model was experimentally validated, statistical approaches for detection, localization, and imaging of objects hidden in fog were developed and demonstrated. A binary hypothesis test and the Neyman-Pearson lemma provided the highest theoretically possible probability of detection for a specified false alarm rate and signal-to-noise ratio. Maximum likelihood estimation allowed estimation of the fog optical properties as well as the position, size, and reflection coefficient of an object in fog. A computational dehazing approach was implemented to reduce the effects of scatter on images, making object features more readily discernible. We have developed, characterized, and deployed a new Tabletop Fog Chamber capable of repeatably generating multiple unique fog-analogues for optical testing in degraded visual environments. We characterized this chamber using both optical and microphysical techniques. In doing so we have explored the ability of droplet nucleation theory to describe the aerosols generated within the chamber, as well as Mie scattering theory to describe the attenuation of light by said aerosols, and correlated the aerosol microphysics to optical properties such as transmission and meteorological optical range (MOR). This chamber has proved highly valuable and has supported multiple efforts inclusive to and exclusive of this LDRD project to test optics in degraded visual environments. Circularly polarized light has been found to maintain its polarization state better than linearly polarized light when propagating through fog. This was demonstrated experimentally in both the visible and short-wave infrared (SWIR) by imaging targets made of different commercially available retroreflective films. It was found that active circularly polarized imaging can increase contrast and range compared to linearly polarized imaging. We have completed an initial investigation of the capability for machine learning methods to reduce the effects of light scattering when imaging through fog. Previously acquired experimental long-wave images were used to train an autoencoder denoising architecture. Overfitting was found to be a problem because of lack of variability in the object type in this data set. The lessons learned were used to collect a well labeled dataset with much more variability using the Tabletop Fog Chamber that will be available for future studies. We have developed several new sensing methods using speckle intensity correlations. First, the ability to image moving objects in fog was shown, establishing that our unique speckle imaging method can be implemented in dynamic scattering media. Second, the speckle decorrelation over time was found to be sensitive to fog composition, implying extensions to fog characterization. Third, the ability to distinguish macroscopically identical objects on a far-subwavelength scale was demonstrated, suggesting numerous applications ranging from nanoscale defect detection to security. Fourth, we have shown the capability to simultaneously image and localize hidden objects, allowing the speckle imaging method to be effective without prior object positional information. Finally, an interferometric effect was presented that illustrates a new approach for analyzing speckle intensity correlations that may lead to more effective ways to localize and image moving objects. All of these results represent significant developments that challenge the limits of the application of speckle imaging and open important application spaces. A theory was developed and simulations were performed to assess the potential transverse resolution benefit of relative motion in structured illumination for radar systems. Results for a simplified radar system model indicate that significant resolution benefits are possible using data from scanning a structured beam over the target, with the use of appropriate signal processing.

58 GEOSCIENCES↗
DOE OSTI2023

Machine learning and deep learning for mineralogy interpretation and CO 2 saturation estimation in geological carbon Storage: A case study in the Illinois Basin

Carbon capture and storage (CCS) is a promising approach to simultaneously maintaining energy security and reducing carbon dioxide (CO 2 ) emissions under the current energy portfolio that is dominated by fossil fuel energy. Pre-injection formation characterization and post-injection CO 2 monitoring are two critical tasks to guarantee storage efficiency in CCS. The CCS projects in the Illinois Basin, the first large-scale CO 2 injection into saline aquifers in the United States, employed conventional and the latest pulsed neutron logging (PNL) tools for mineralogy interpretation and CO 2 saturation estimation, which provide valuable references for future CCS projects. Because of the inherent fuzziness of petrophysical measurements and complex subsurface heterogeneity, interpreting well-logging data is time-consuming, and its accuracy can be user-biased. In recent years, data-driven methods have been widely used to capture the non-linear patterns between input features and interpretation results. This work applied and evaluated four commonly used machine learning (ML) models, including ridge regression (RR), random forest (RF), gradient boosting regression (GBR), support vector regression (SVR), and one deep learning (DL) model, the artificial neural network (ANN). We optimized the hyperparameters of the four ML models and the DL model using the simulated annealing algorithm and the grid search strategy, respectively. The input features of the mineralogy interpretation models were eleven conventional well-logging parameters, and the label data (i.e., ground truth) were the porosity and volumetric fractions of six minerals, including quartz, feldspar, dolomite, calcite, clay, and iron minerals. The results demonstrated that the GBR and RF models were superior in predicting volumetric fractions of minerals and porosity; label data with low coefficient of variation (CV) values tended to yield better performance. For CO 2 saturation estimation, the RF was the best-performing model, followed by SVR, ANN, GBR, and RR. Furthermore, we conducted feature importance ranking using the permutation importance algorithm and found that the formation sigma and well pressure were the most important features in this study. In conclusion, the study of CCS projects in the Illinois Basin bridges the gap between the limited knowledge and understanding of geological carbon storage and the increasing demand for reliable, cost-effective, and sustainable energy solutions.

58 GEOSCIENCES↗
DOE OSTI2025

Sign Problem in Tensor-Network Contraction

We investigate how the computational difficulty of contracting tensor networks depends on the sign structure of the tensor entries. Using results from computational complexity, we observe that the approximate contraction of tensor networks with only positive entries has lower computational complexity as compared to tensor networks with general real or complex entries. This raises the question of how this transition in computational complexity manifests itself in the hardness of different tensor-network-contraction schemes. We pursue this question by studying random tensor networks with varying bias toward positive entries. First, we consider contraction via Monte Carlo sampling and find that the transition from hard to easy occurs when the tensor entries become predominantly positive; this can be understood as a tensor-network manifestation of the well-known negative-sign problem in quantum Monte Carlo. Second, we analyze the commonly used contraction based on boundary tensor networks. The performance of this scheme is governed by the number of correlations in contiguous parts of the tensor network (which by analogy can be thought of as entanglement). Remarkably, we find that the transition from hard to easy—i.e., from a volume-law to a boundary-law scaling of entanglement—already occurs for a slight bias of the tensor entries toward a positive mean, scaling inversely with the bond dimension D , and thus the problem becomes easy the earlier the larger D occurs. This is in contrast both to expectations and to the behavior found in Monte Carlo contraction, where the hardness at fixed bias increases with the bond dimension. To provide insight into this early breakdown of computational hardness and the accompanying entanglement transition, we construct an effective classical statistical-mechanical model that predicts a transition at a bias of the tensor entries of 1 / D , confirming our observations. We conclude by investigating the computational difficulty of computing expectation values of tensor-network wave functions (projected entangled-pair states, PEPSs) and find that in this setting, the complexity of entanglement-based contraction always remains low. We explain this by providing a local transformation that maps PEPS expectation values to a positive-valued tensor network. This not only provides insight into the origin of the observed boundary-law entanglement scaling but also suggests new approaches toward PEPS contraction based on positive decompositions. Published by the American Physical Society 2025

Chen, Jielun (ORCID:0000000178411545)↗
DOE OSTI2026

Heat exposure and maternal stress: evidence from the GRAPHS pregnancy cohort in Ghana

Heat exposure has been linked to psychosocial stress, an established antecedent of perinatal depression; however, evidence on heat-related stress during pregnancy in sub-Saharan Africa remains limited. We analyzed psychosocial stress scores and covariate data from the Ghana Randomized Air Pollution and Health Study, linking daily maximum and minimum shaded wet bulb globe temperature (WBGT) metrics to participants’ stress scores derived from the Crisis in Family Systems-Revised Life Events Questionnaire. We evaluated associations using ordinal logistic regression of pregnancy-average and trimester-average exposures and distributed lag non-linear models (DLNMs) to assess time-varying associations across gestation. Higher average maximum WBGT exposure across pregnancy was associated with increased odds of higher psychosocial stress; each 1 °C increase in maximum WBGT was associated with 64% higher odds of belonging to a higher stress category (OR = 1.64; 95% CI = 1.17–2.31; p = 0.0040). In trimester-average models, higher first-trimester maximum WBGT was also associated with higher stress (OR = 1.44; 95% CI = 1.15–1.81; p = 0.0014). DLNMs suggested that relatively cooler daily maximum WBGT values (25th percentile) were associated with decreased odds of stress in early pregnancy, whereas extreme daily maximum WBGT values (99th percentile) showed a pattern consistent with increased odds of stress in mid-to-late gestation. These findings highlight gestational windows in which heat exposure may influence stress, emphasizing the need for further research into underlying mechanisms and effective interventions to protect maternal mental health in heat-vulnerable settings.

White, Lewis [Columbia University] (ORCID:00090005↗