Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Natural Probabilities”

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 271 records · Page 15

Hydrogen Plus Other Alternative Fuels Risk Assessment Models (HyRAM+) (Technical Reference Manual V.4)

The HyRAM+ software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen, methane, and propane infrastructure and transportation systems. HyRAM+ is designed to facilitate the use of state-of-the-art science and engineering models to conduct robust, repeatable assessments of safety, hazards, and risk. HyRAM+ includes generic probabilities for equipment failures, probabilistic models for the impact of heat flux on humans and structures, and experimentally validated first-order models of release and flame physics. HyRAM+ integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, and characterizing hazards (thermal effects from jet fires, overpressure effects from delayed ignition), and assessing impact on people and structures. HyRAM+ is developed at Sandia National Laboratories to support the development and revision of national and international codes and standards. HyRAM+ is a research software in active development and thus the models and data may change. This report will be updated at appropriate developmental intervals. This document provides a description of the methodology and models contained in HyRAM+ version 4.0. The most significant change for HyRAM+ version 4.0 from HyRAM version 3.1 is the incorporation of other alternative fuels, namely methane (as a proxy for natural gas) and propane into the toolkit. This change necessitated significant changes to the installable graphical user interface as well as changes to the back-end Python models. A second major change is the inclusion of physics models for the overpressure associated with the delayed ignition of an unconfined jet/plume of flammable gas.

02 PETROLEUM↗

Enhancing quantum annealing accuracy through replication-based error mitigation *

Abstract Quantum annealers like those manufactured by D-Wave Systems are designed to find high quality solutions to optimization problems that are typically hard for classical computers. They utilize quantum effects like tunneling to evolve toward low-energy states representing solutions to optimization problems. However, their analog nature and limited control functionalities present challenges to correcting or mitigating hardware errors. As quantum computing advances towards applications, effective error suppression is an important research goal. We propose a new approach called replication based mitigation (RBM) based on parallel quantum annealing (QA). In RBM, physical qubits representing the same logical qubit are dispersed across different copies of the problem embedded in the hardware. This mitigates hardware biases, is compatible with limited qubit connectivity in current annealers, and is well-suited for currently available noisy intermediate-scale quantum annealers. Our experimental analysis shows that RBM provides solution quality on par with previous methods while being more flexible and compatible with a wider range of hardware connectivity patterns. In comparisons against standard QA without error mitigation on larger problem instances that could not be handled by previous methods, RBM consistently gets better energies and ground state probabilities across parameterized problem sets.

Djidjev, Hristo N. (ORCID:0000000192868824)↗

Optimizing Materials to Boost the Valorization of CO 2 : Tuning Cobalt–Cobalt Interactions on In 2 O 3 -Based Photothermal Catalysts

The valorization of CO 2 is an important challenge within the current panorama, since this molecule is probably the main contributor to climate change. In this study, the synthesis of materials based on a nanostructured batonnet-type indium oxide is carried out. In them, different amounts of Co are introduced, varying between 2 and 8% mol. It is verified that the most active sample in the transformation of carbon dioxide to carbon monoxide contains 6 mol %. of Co. This sample’s activity under dual excitation exceeds the thermal counterpart by more than 30%. After carrying out a complete physical and chemical characterization with the help of X-ray absorption spectroscopy and other techniques, it is shown that catalysts with amounts of cobalt equal to or below 4 mol % contain isolated single-atom species, while those with higher amounts of metal display a Co–Co interaction which triggers the evolution of the samples under reaction conditions. The optimum control of this Co–Co interaction and the nature of the final cobalt-containing species determine dual photothermal catalytic properties. In conclusion, this work establishes a structure–activity relationship to interpret the catalytic behavior of highly dispersed subnanometric cobalt species, and thus an avenue to optimize the photothermal valorization of carbon dioxide.

36 MATERIALS SCIENCE↗

Probing for the host galaxies of the fast X-ray transients XRT 000519 and XRT 110103

Over the past few years, ~30 extragalactic fast X-ray transients (FXRTs) have been discovered, mainly in Chandra and XMM-Newton data. Their nature remains unclear, with proposed origins, including a double neutron star merger, a tidal disruption event involving an intermediate-mass black hole and a white dwarf, or a supernova shock breakout. A decisive differentiation between these three promising mechanisms for their origin requires an understanding of the FXRT energetics, environments, and/or host properties. We present optical observations obtained with the Very Large Telescope for the FXRTs XRT 000519 and XRT 110103 and Gran Telescopio Canarias observations for XRT 000519 designed to search for host galaxies of these FXRTs. In the g s , r s , and R-band images, we detect an extended source on the north-west side of the $\sim \, 1^{\prime \prime }$ (68 per cent confidence) error circle of the X-ray position of XRT 000519 with a Kron magnitude of g s = 26.29 ± 0.09 (AB magnitude). We discuss the XRT 000519 association with the probable host candidate for various possible distances, and we conclude that if XRT 000519 is associated with the host candidate a supernova shock breakout scenario is likely excluded. No host galaxy is found near XRT 110103 down to a limiting magnitude of R > 25.8.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of Parameters for the Particle Size Distribution of TATB

Laser light scattering (LLS), manual counting of scanning electron microscopy (SEM) images, and reduction of SEM images using ImageJ (open source program) to determine Feret (caliper) diameters were applied to determine particle size distribution (PSD) of four preparations of 1,3,5‐triamino‐2,4,6‐trinitrobenzene (TATB) and yttria‐stabilized zirconia (YSZ), an SEM certified standard. Mie theory was used to reduce the LLS data. The spherical nature of the YSZ made it a good candidate for LLS. Variations in n , the refractive index, and iκ, the imaginary component, produced very little change in the PSD. However, changing the carrier liquid from H 2 O to a 40% aqueous sucrose solution, thereby changing the carrier refractive index, n 0 , substantially affected the PSD. The Mie complex refractive indices for the YSZ were n = 2.200, iκ = 0.100, with a 40% aqueous sucrose solution, n 0 = 1.400. The triclinic crystal structure of TATB made refractive index determinations more difficult, so a study was conducted varying Mie parameters and comparing them to the same data reduced using the Fraunhofer theory. Changing the n and κ parameters produced PSD with a small concentration of particles less than 1 µm in size or none in this range. SEM images, Feret data, manual counting, and Fraunhofer data reduction indicate particles less than 1 µm are probably < 5% in concentration. The final selection of Mie parameters for TATB was n = 2.283, iκ = 0.1, and suspension medium, n 0 = 1.330. Finally, computations, using density functional theory produced similar parameters.

Feret diameter↗

Solving independent set problems with photonic quantum circuits

An independent set (IS) is a set of vertices in a graph such that no edge connects any two vertices. In adiabatic quantum computation [E. Farhi, et al ., Science 292, 472–475 (2001); A. Das, B. K. Chakrabarti, Rev. Mod. Phys. 80, 1061–1081 (2008)], a given graph G ( V , E ) can be naturally mapped onto a many-body Hamiltonian H IS G ( V , E ) , with edges E being the two-body interactions between adjacent vertices V . Thus, solving the IS problem is equivalent to finding all the computational basis ground states of H IS G ( V , E ) . Very recently, non-Abelian adiabatic mixing (NAAM) has been proposed to address this task, exploiting an emergent non-Abelian gauge symmetry of H IS G ( V , E ) [B. Wu, H. Yu, F. Wilczek, Phys. Rev. A 101, 012318 (2020)]. Here, we solve a representative IS problem G ( 8 , 7 ) by simulating the NAAM digitally using a linear optical quantum network, consisting of three C-Phase gates, four deterministic two-qubit gate arrays (DGA), and ten single rotation gates. The maximum IS has been successfully identified with sufficient Trotterization steps and a carefully chosen evolution path. Remarkably, we find IS with a total probability of 0.875(16), among which the nontrivial ones have a considerable weight of about 31.4%. Our experiment demonstrates the potential advantage of NAAM for solving IS-equivalent problems.

Science & Technology - Other Topics↗

BeyondPlanck I. Global Bayesian analysis of the Planck Low Frequency Instrument data

We describe the BeyondPlanck project in terms of motivation, methodology and main products, and provide a guide to a set of companion papers that describe each result in fuller detail. Building directly on experience from ESA's Planck mission, we implement a complete end-to-end Bayesian analysis framework for the Planck Low Frequency Instrument (LFI) observations. The primary product is a joint posterior distribution P(omega|d), where omega represents the set of all free instrumental (gain, correlated noise, bandpass etc.), astrophysical (synchrotron, free-free, thermal dust emission etc.), and cosmological (CMB map, power spectrum etc.) parameters. Some notable advantages of this approach are seamless end-to-end propagation of uncertainties; accurate modeling of both astrophysical and instrumental effects in the most natural basis for each uncertain quantity; optimized computational costs with little or no need for intermediate human interaction between various analysis steps; and a complete overview of the entire analysis process within one single framework. As a practical demonstration of this framework, we focus in particular on low-l CMB polarization reconstruction, paying special attention to the LFI 44 GHz channel. We find evidence of significant residual systematic effects that are still not accounted for in the current processing, but must be addressed in future work. These include a break-down of the 1/f correlated noise model at 30 and 44 GHz, and scan-aligned stripes in the Southern Galactic hemisphere at 44 GHz. On the Northern hemisphere, however, we find that all results are consistent with the LCDM model, and we constrain the reionization optical depth to tau = 0.067 +/- 0.016, with a low-resolution chi-squared probability-to-exceed of 16%. The marginal CMB dipole amplitude is 3359.5 +/- 1.9 uK. (Abridged.)

Andersen, KJ↗

Emergent constraints on future methane emissions from global wetlands

Future methane (CH 4 ) emissions from natural wetlands are predicted to increase due to global warming, leading to positive feedback on climate change. However, the magnitude of this increase remains highly uncertain. Here we present novel ensemble simulations of seven state-of-the-art terrestrial biosphere models to estimate wetland CH 4 emissions (eCH 4 ) during the twenty-first century. Our estimates suggest that for every 1 °C increase in global land surface temperature, there is a 24 ± 10 Tg CH 4 yr −1 increase in eCH 4 . We also identify an emergent relationship between contemporary temperature dependence and projected eCH 4 . When constrained by 163 site-year eddy-covariance measurements of eCH 4 , we show that wetland emissions can increase by 50–60% by the 2090s relative to the 2010s under a high-warming scenario. The projected decadal increase in eCH 4 from the 2010–2019 baseline to the 2030s would very likely (90% probability) offset an amount equivalent in scale to 8–10% of anthropogenic eCH 4 at the 2020 level, comparable to the reductions committed under the Global Methane Pledge. However, the constraint is dominated by mid- and high-latitude observations, with limited tropical coverage, and uncertainties in projected wetland inundation contribute substantially to uncertainty in eCH 4 . Our findings reduce the uncertainty in projected wetland methane–climate feedback and highlight its potential impacts on methane mitigation efforts to slow global warming.

Zhang, Zhen [Chinese Academy of Sciences (CAS), Be↗

Draft genome of multiple resistance donor plant Sinapis alba: An insight into SSRs, annotations and phylogenetics

Sinapis alba is a wild member of the Brassicaceae family reported to possess genetic resistance against major biotic and abiotic stresses of oilseed brassicas. However, the resistance nature of S. alba was not exploited generously due to the unavailability of usable genome sequences in public databases. Therefore, the present study was conducted to assemble the first draft genome from raw whole genome shotgun sequences with annotation and develop simple sequence repeat markers for molecular genetics and marker-assisted breeding. Results The raw genome sequences had 96x coverage on the Illumina platform with 170 Gbp data. The developed assembly by SOAPdenovo2 has ~459 Mbp genome size covered in 403,423 contigs with an average size of 1138.04 bp. The assembly was BLASTX with Arabidopsis thaliana which showed 32.9% positive hits between both plants. The top hit species distribution analysis showed the highest similarity with A. thaliana. A total of 809,597 GO level annotations were recorded after BLASTX results, and 34,012 sequences were annotated with different enzyme codes grouped under seven classes. The gene prediction tool AUGUSTUS identified 113,107 probable genes with an average size of 684 bp. The biochemical pathway annotation assigned 16,119 potential genes to 152 KEGG maps and 1751 enzyme codes. The development of potential SSRs from the de-novo assembly yielded 70731 unique primer pairs. Out of 159 randomly selected SSR markers for validation, 149 successfully amplified in S. alba. However, 10 SSR markers did not amplify during the validation experiment. Conclusion The annotated genome assembly with a large number of SSRs was developed in the present study. To the best of our knowledge, this is the first report of S. alba genome assembly development, annotation, and SSRs mining to date. The data presented here will be a very important resource for future crop improvement programs, especially for resistant breeding.

59 BASIC BIOLOGICAL SCIENCES↗

A Machine Learning-based Reliability Evaluation Model for Integrated Power-Gas Systems

This article proposes a hybrid machine learning method for the reliability evaluation of integrated power-gas systems (IPGS) under the uncertain component failure probability distributions. The Random Forest (RF) method is designed to select important features to solve the insufficient quantity of data and the curse of dimensionality problems. The Extreme Gradient Boosting (XGBoost) regression algorithm is developed to quantify the relationship between the uncertain parameters and reliability metrics. Moreover, a ten-fold cross-validation method is employed to further improve the accuracy of the regression model. Simulation results on three test systems show that the proposed method can achieve high accuracy for the reliability evaluation.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Multistep Coulomb excitation of Ni 64 : Shape coexistence and nature of low-spin excitations

Here, the structure of 64 Ni, the heaviest stable Ni isotope, has been investigated via high-statistics, multistep safe Coulomb excitation to search for shape coexistence, a phenomenon recently observed in neutron-rich 66 Ni and 70 Ni as well as in doubly magic, N = 40, 68 Ni. The study was motivated by recent, state-of-the-art Monte Carlo shell-model calculations (MCSM), where a Hamiltonian with effective interactions incorporating the monopole tensor force predicts the existence of shape coexistence, also in the lower-mass 62,64 Ni isotopes. A set of transition and static E2 matrix elements for both yrast and near-yrast structures was extracted from the differential Coulomb excitation cross sections. From comparisons between the new results and MCSM as well as other shell-model calculations, a clearer picture of the structure of 64 Ni emerges. Specifically, the low-spin states are shown to be dominated by proton and neutron excitations mainly within the fp shell, with minimal contribution from the g 9/2 shape-driving neutron orbital. The agreement between experimental data and MCSM results indicates a small oblate deformation for the $0$$^{+}_{2}$ level and a spherical shape for the $0$$^{+}_{3}$ state. In addition, the small upper limit determined for the B(E2) probability of a transition associated with the decay of the recently observed 3463-keV, $0$$^{+}_{4}$ state agrees with its proposed assignment to a prolate shape, herewith providing first evidence for triple shape coexistence in a stable Ni isotope.

59 ≤ A ≤ 89↗

Modern Monte Carlo methods for efficient uncertainty quantification and propagation: A survey

Uncertainty quantification (UQ) includes the characterization, integration, and propagation of uncertainties that result from stochastic variations and a lack of knowledge or data in the natural world. Monte Carlo (MC) method is a sampling-based approach that has widely used for quantification and propagation of uncertainties. However, the standard MC method is often time-consuming if the simulation-based model is computationally intensive. This article gives an overview of modern MC methods to address the existing challenges of the standard MC in the context of UQ. Specifically, multilevel Monte Carlo (MLMC) extending the concept of control variates achieves a significant reduction of the computational cost by performing most evaluations with low accuracy and corresponding low cost, and relatively few evaluations at high accuracy and corresponding high cost. Multifidelity Monte Carlo (MFMC) accelerates the convergence of standard Monte Carlo by generalizing the control variates with different models having varying fidelities and varying computational costs. Multimodel Monte Carlo method (MMMC), having a different setting of MLMC and MFMC, aims to address the issue of UQ and propagation when data for characterizing probability distributions are limited. Multimodel inference combined with importance sampling is proposed for quantifying and efficiently propagating the uncertainties resulting from small data sets. All of these three modern MC methods achieve a significant improvement of computational efficiency for probabilistic UQ, particularly uncertainty propagation. Finally, an algorithm summary and the corresponding code implementation are provided for each of the modern MC methods. The extension and application of these methods are discussed in detail.

97 MATHEMATICS AND COMPUTING↗

Structure of Pb(Fe 2/3 W 1/3 )O 3 single crystals with partial cation order

Despite intensive studies on the complex perovskite Pb(Fe 2/3 W 1/3 )O 3 (PFWO) relaxor, understanding the exact nature of its multifunctional properties has remained a challenge for decades. In this work we report a comprehensive structural study of the PFWO single crystals using a combination of synchrotron X-ray diffraction and high-resolution electron microscopy. The set of {h + ½, k + ½, l + ½} superlattice reflections was observed for the first time based on single-crystal synchrotron X-ray experiments (100–450 K) and transmission electron microscopy investigations, which indicates some kind of B-cation ordering in PFWO which had been thought to be totally disordered. It was found that (1) the crystal structure of PFWO should be described by a partly ordered cubic perovskite (i.e. Fm – 3m ), (2) the weak ferromagnetic properties and excess magnetic moment of PFWO can be understood based on non-random distribution of Fe cations between the 4a and 4b sites, and (3) the Pb displacement disorder is present in this material and the cations are probably displaced along the <100> directions. The X-ray diffraction results of this investigation show that partial cation ordering indeed exists in PFWO, which makes it necessary to revisit the generally accepted interpretations of the results obtained up to date. In agreement with X-ray diffraction study the main results of TEM study include: (1) a long range order that can be described with the Fm – 3m symmetry is reliably detected, (2) the coherence length of that long range order is in the order of 1–2 nm and (3) no remarkable chemical inhomogeneity is found in the tested PFWO crystal, excluding the possibility of a compositional ordering arising from substitutional defects in the perovskite structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interpreting machine learning prediction of fire emissions and comparison with FireMIP process-based models

Annual burned areas in the United States have increased 2-fold during the past decades. With more large fires resulting in more emissions of fine particulate matter, an accurate prediction of fire emissions is critical for quantifying the impacts of fires on air quality, human health, and climate. This study aims to construct a machine learning (ML) model with game-theory interpretation to predict monthly fire emissions over the contiguous US (CONUS) and to understand the controlling factors of fire emissions. The optimized ML model is used to diagnose the process-based models in the Fire Modeling Intercomparison Project (FireMIP) to inform future development. Results show promising performance for the ML model, Community Land Model (CLM), and Joint UK Land Environment Simulator-Interactive Fire And Emission Algorithm For Natural Environments (JULES-INFERNO) in reproducing the spatial distributions, seasonality, and interannual variability of fire emissions over the CONUS. Regional analysis shows that only the ML model and CLM simulate the realistic interannual variability of fire emissions for most of the subregions (r >0.95 for ML and r =0.14~0.70 for CLM), except for Mediterranean California, where all the models perform poorly (r =0.74 for ML and r <0.30 for the FireMIP models). Regarding seasonality, most models capture the peak emission in July over the western US. However, all models except for the ML model fail to reproduce the bimodal peaks in July and October over Mediterranean California, which may be explained by the smaller wind speeds of the atmospheric forcing data during Santa Ana wind events and limitations in model parameterizations for capturing the effects of Santa Ana winds on fire activity. Furthermore, most models struggle to capture the spring peak in emissions in the southeastern US, probably due to underrepresentation of human effects and the influences of winter dryness on fires in the models. As for extreme events, both the ML model and CLM successfully reproduce the frequency map of extreme emission occurrence but overestimate the number of months with extremely large fire emissions. Comparing the fire PM 2.5 emissions from the ML model with process-based fire models highlights their strengths and uncertainties for regional analysis and prediction and provides useful insights into future directions for model improvements.

54 ENVIRONMENTAL SCIENCES↗

Targeted engineering of camelina and pennycress seeds for ultrahigh accumulation of acetyl-TAG

Acetyl-TAG (3-acetyl-1,2-diacylglycerol), unique triacylglycerols (TAG) possessing an acetate group at the sn -3 position, exhibit valuable properties, such as reduced viscosity and freezing points. Previous attempts to engineer acetyl-TAG production in oilseed crops did not achieve the high levels found in naturally producing Euonymus seeds. Here, we demonstrate the successful generation of camelina and pennycress transgenic lines accumulating nearly pure acetyl-TAG at 93 mol% and 98 mol%, respectively. These ultrahigh acetyl-TAG synthesizing lines were created using gene-edited FATTY ACID ELONGASE1 ( FAE1 ) mutant lines as an improved genetic background to increase levels of acetyl-CoA available for acetyl-TAG synthesis mediated by the expression of EfDAcT, a high-activity diacylglycerol acetyltransferase isolated from Euonymus fortunei . Combining EfDAcT expression with suppression of the competing TAG-synthesizing enzyme DGAT1 further enhanced acetyl-TAG accumulation. These ultrahigh levels of acetyl-TAG exceed those in earlier engineered oilseeds and are equivalent or greater than those in Euonymus seeds. Imaging of lipid localization in transgenic seeds revealed that the low amounts of residual TAG were mostly confined to the embryonic axis. Similar spatial distributions of specific TAG and acetyl-TAG molecular species, as well as their probable diacylglycerol (DAG) precursors, provide additional evidence that acetyl-TAG and TAG are both synthesized from the same tissue-specific DAG pools. Remarkably, this ultrahigh production of acetyl-TAG in transgenic seeds exhibited minimal negative effects on seed properties, highlighting the potential for production of designer oils required for economical biofuel industries.

09 BIOMASS FUELS↗

Muon-induced fission of actinide nuclei

A negative muon captured by an actinide cascades down through the muonic atomic levels; deeply bound transitions can proceed via inverse internal conversion, depositing the muonic transition energy directly into the nucleus and, when the deposited energy exceeds the fission barrier, inducing prompt fission. Because the muon mean lifetime exceeds the saddle-to-scission timescale by orders of magnitude, the muon can survive the entire fission event as a 1⁢𝑠 spectator and ultimately attach to one or both of the emerging fragments. Its postscission attachment probability to the light fragment, 𝑃 𝐿 , can be used as a direct electromagnetic probe of fission dynamics on a timescale of 10 −21 s. In previous work, we introduced a three-dimensional lattice solution of the time-dependent Dirac equation coupled to the electromagnetic field generated by a fissioning nucleus and reported 𝑃 𝐿 for several actinides at a single dissipation strength. In this work, we extend that framework to a systematic survey of 232 Th , 238 U , and 240 Pu and implement a more realistic fission model which incorporates dynamic pairing correlations. We find that 𝑃 𝐿 falls steeply with the fragment charge asymmetry, a robust structural fingerprint of the fissioning system, while its dependence on nuclear dissipation is secondary and sensitive to the phenomenological friction prescription. These results establish 𝑃 𝐿 as a clean electromagnetic probe of fragment charge asymmetry and motivate a self-consistent, coordinate- and time-dependent treatment of nuclear dissipation as the natural next step.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Cumulative short-lived photofission product yields for nuclear forensic application

Experiments have long been carried out with linear accelerator (linac) produced bremsstrahlung X-rays to study high-energy photon-induced fission, known as “photofission”. One benefit of photofission, when compared with neutron-induced fission, is the ability to investigate fundamental physics of fission at excitation energies lower than what is possible with neutrons [1]. Absorption of an incident neutron results in a compound nucleus with an excitation energy at least equal to the neutron binding energy; thus, the lowest energy region of the fission barrier cannot be studied. Fission product yield distributions are dependent on the incident particle energy and the target nucleus mass. They are generally represented by a double-humped curve, with low and high mass “peaks” and an intermediate mass “valley”. As incident particle energy increases, the excitation energy of the nucleus increases, thereby increasing the overall yield of isotopes in the valley with respect to the two maxima or peaks of the fission product curve [3]. Photofission is most probable for incident photons around 14 MeV due to the giant dipole resonance. Induced photofission around this energy results in a fission product yield distribution that resembles 14 MeV DT neutron fission more than thermal or fast (500 keV) neutron fission. In post-detonation nuclear forensics, the peak-to-valley ratio of the fission product yield curve could be used for characterizing device type by providing an indication of the neutron energy spectrum that had induced fission (e.g., 500 keV vs. 14 MeV). The goal of this work is to improve nuclear data for isotopes of interest in nuclear forensics and special nuclear material (SNM) detection, specifically the fission product yield data via photofission experiments for the isotopes of 238U and 232Th. Photofission has been proposed as a less costly and more flexible production method for the valley isotopes that are representative of DT neutron fission product yield distributions, which would benefit nuclear forensics exercises in testing for analysis techniques. Natural uranium and thorium targets were irradiated with bremsstrahlung photons (nominal endpoint energies of 8, 14, and 20 MeV), generated via impingement of high-energy electrons on a tungsten electron-to-photon radiator. Gamma-ray spectroscopy was used to observe and identify the decay of the resulting fission products. The photofission yields of 102 fission products were measured for the set of the 238U and 232Th targets at the three bremsstrahlung X-ray endpoint energies 8, 14, and 20 MeV. These fragments included masses ranging from A=84 to A=144, with half-lives as short as 1.07 seconds.

07 ISOTOPE AND RADIATION SOURCES↗

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↗