Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sample complexity analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,441 records · Page 80

Managing Uncertainty Due to a Fundamental Error Source Arising from Scatterer Distribution Complexity in Radar Remote Sensing of Precipitation

The assumption that cloud and rain drops are spatially distributed according to a Poisson distribution within a scattering volume probed by a radar being used to estimate precipitation has represented bedrock theory in establishing 'rules of the game' for pulse averaging--the process needed to beat down noise to an acceptable level in the measurement of radar reflectivity factor. Based on relatively recent observations of 'realistic' spatial distributions of hydrometeor scatterers in a cloudy atmosphere motivates a renewed examination of the consequences of using a too simplified assumption underlying volume scattering--particularly in regards to the standard pulse averaging rule. Our investigation addresses two extremes, simple to complex, insofar as allowed for complexities in an underlying scatterer distribution. It is demonstrated that as the spatial distribution ranges from Poisson (a narrow distribution) to multi-fractal (much broader distribution), uncertainty in a measurement increases if the rule for pulse averaging goes unchanged from its Poisson distribution reference county. [A bounded cascade is used for the multi-fractal distribution, a regularly observed distribution vis-a-vis cloud liquid water content.] The resultant measurement uncertainty leads to a fundamental source of error in the estimation of rain rate from radar measurements, one that has been disregarded since the early 1950s when radar sets first began to be used for rainfall measuring. It is shown how this source of error can be 'managed'--under the assumption that number of data analysis experiments would be carried out, experiments involving pulse-by-pulse measurements obtained from a radar set modified to output individual pulses of reflectivity factor. For practical applications, a new parameter called normalized k-sample intensity invariance is developed to enable defining the required pulse average count according to a preferred degree of uncertainty.

Smith, Eric A.↗

Hydrogen Analysis by Gas Chromatography–Mass Spectrometry

The detection of hydrogen in complex gas mixtures is essential for many applications. Conventional approaches such as gas chromatography (GC) with thermal conductivity detection (TCD) and residual gas analyzers (RGAs) face significant limitations: TCD exhibits poor response when helium is used as a carrier gas, whereas RGAs lack chromatographic separation, preventing reliable quantification of hydrogen because of interference from other species. Traditional GC methods rely on dual-column configurations with packed or molecular-sieve porous layer open tubular (PLOT) columns to separate hydrogen from O2, N2, CO, CO2, CH4, and other hydrocarbons, increasing system complexity and limiting compatibility with mass spectrometry (MS) detection. In this work, we developed and validated a robust GC–MS method capable of directly detecting and quantifying hydrogen using electron ionization (EI) without dopants, reagent gases, or ion–molecule reaction schemes. By integrating a modified EI source and a cryogenically cooled single capillary column configuration, we achieved baseline separation of hydrogen from all major permanent gases and hydrocarbons in a refinery gas mixture using helium as the carrier gas. The method demonstrated excellent linearity, high sensitivity, and exceptional reproducibility. Adjustable sample-loop volumes and split ratios enabled optimization of peak shape and signal-to-noise performance for trace-level and percent-level hydrogen concentrations. Beyond hydrogen quantitation, the method provides simultaneous compositional profiling of other gases in a mixture in a single run, making it valuable for a wide range of tasks.

Lobodin, Vlad [ORNL]↗

Probabilistic Structural Analysis Program

NASA/NESSUS 6.2c is a general-purpose, probabilistic analysis program that computes probability of failure and probabilistic sensitivity measures of engineered systems. Because NASA/NESSUS uses highly computationally efficient and accurate analysis techniques, probabilistic solutions can be obtained even for extremely large and complex models. Once the probabilistic response is quantified, the results can be used to support risk-informed decisions regarding reliability for safety-critical and one-of-a-kind systems, as well as for maintaining a level of quality while reducing manufacturing costs for larger-quantity products. NASA/NESSUS has been successfully applied to a diverse range of problems in aerospace, gas turbine engines, biomechanics, pipelines, defense, weaponry, and infrastructure. This program combines state-of-the-art probabilistic algorithms with general-purpose structural analysis and lifting methods to compute the probabilistic response and reliability of engineered structures. Uncertainties in load, material properties, geometry, boundary conditions, and initial conditions can be simulated. The structural analysis methods include non-linear finite-element methods, heat-transfer analysis, polymer/ceramic matrix composite analysis, monolithic (conventional metallic) materials life-prediction methodologies, boundary element methods, and user-written subroutines. Several probabilistic algorithms are available such as the advanced mean value method and the adaptive importance sampling method. NASA/NESSUS 6.2c is structured in a modular format with 15 elements.

Pai, Shantaram S.↗

Analysis of Exposure during Transportation of Waste to the Nevada National Security Site (NNSS) - 20445

Safety is a priority for the U.S. Department of Energy (DOE) Environmental Management (EM) Nevada Program. This commitment to safety is integrated into requirements and work procedures at DOE sites across the U.S. that generate low-level radioactive waste, mixed low-level radioactive waste, and classified waste that is transported to the Nevada National Security Site (NNSS) for permanent disposal. Generators of this waste are responsible for ensuring it is packaged and transported in a safe and compliant manner. These generators and their contracted commercial carriers must comply with U.S. Department of Transportation (DOT) regulations; the Nevada National Security Site Waste Acceptance Criteria (NNSSWAC); and other applicable federal, state, and local regulations and requirements. This includes packaging and shipping the waste to meet DOT regulatory radiation dose limits to ensure public safety during transportation. NNSS waste generators have a solid safety record for complying with these stringent regulations as demonstrated by the fact that no shipments exceeding DOT dose rate limits were received at the NNSS Area 5 Radioactive Waste Management Complex (RWMC) during FY 2014 through FY 2018. A 2005 study conducted by the Desert Research Institute validated shipments to the NNSS comply with DOT radiation dose limits. In response to recent Nevada stakeholder requests for an updated study, an analysis was conducted to verify and provide evidence that the public continues to be safe from exposure to radioactive waste shipment to the NNSS. Radiation exposure from radioactive waste shipments in transit to the NNSS was evaluated. Three studies were performed: 1. Confirm compliance of waste shipments to the DOT radiation limits using dose rate surveys taken at the generator site and at the disposal site. 2. Estimate the dose to the public from waste shipments using the dosimetry results for disposal site personnel by comparing time of exposure for various individuals/scenarios. 3. Estimate the truck driver dose using the receipt dose rate surveys and time of exposure, and from dosimetry results. Study 1 confirmed compliance with the DOT dose rate limits for a sample of waste shipments (123 of 1,038) received during fiscal year (FY) 2017. Study 2 showed that the average annual radiation equivalent dose for disposal site personnel during FY 2017 was approximately 6.7% [0.415 mSv (41.5 millirem)] of the annual average dose that a person living in the U.S. receives from all sources [6.2 mSv (620 mrem)]. Maximum dose to a member of the public, a person resting/sleeping in an adjacent truck at a rest area to a waste shipment, was 0.0166 mSv (1.66 mrem), which is less than the dose received from a person taking a round trip from the east coast to the west coast of the continental U.S. on an airplane [0.07 mSv (7 mrem) per round trip]. Study 3 showed that the truck driver who made the most shipments of radioactive waste to the NNSS, who drove the most hours (estimated time of exposure 2,520 hours), during FY 2017 had a total dose of 0.15 mSv (15 mrem), which was obtained from their personal dosimeter. The annual dose recorded for disposal site personnel, the dose estimates for the public, and the truck driver dose were all a small percentage of the average radiation dose a person receives from all sources (e.g., background and medical radiation). At these levels, risk should only be described qualitatively. Radiation exposure from radioactive waste shipments to the NNSS were shown to be compliant with the DOT dose rate limits, safe, and very low (i.e., within background levels). (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Large Deviations Analysis of Certain Qualitative Properties of Parallel Tempering and Infinite Swapping Algorithms

Parallel tempering, or replica exchange, is a popular method for simulating complex systems. The idea is to run parallel simulations at different temperatures, and at a given swap rate exchange configurations between the parallel simulations. From the perspective of large deviations it is optimal to let the swap rate tend to infinity and it is possible to construct a corresponding simulation scheme, known as infinite swapping. In this paper we propose a novel use of large deviations for empirical measures for a more detailed analysis of the infinite swapping limit in the setting of continuous time jump Markov processes. Using the large deviations rate function and associated stochastic control problems we consider a diagnostic based on temperature assignments, which can be easily computed during a simulation. Here we show that the convergence of this diagnostic to its a priori known limit is a necessary condition for the convergence of infinite swapping. Finally, the rate function is also used to investigate the impact of asymmetries in the underlying potential landscape, and where in the state space poor sampling is most likely to occur.

97 MATHEMATICS AND COMPUTING↗

Development of an Ion Mobility Spectrometry-Orbitrap Mass Spectrometer Platform

Complex samples benefit from multidimensional measurements where higher resolution enables more complete characterization of biological and environmental systems. To address this challenge, we developed a drift tube-based ion mobility spectrometry-Orbitrap mass spectrometer (IMS-Orbitrap MS) platform. To circumvent the time scale disparity between the fast IMS separation and the much slower Orbitrap MS acquisition, we utilized a dual gate and pseudorandom sequences to multiplexed injection of ions and allowing operation in signal averaging (SA), single multiplexing (SM) and double multiplexing (DM) IMS modes to optimize the signal-to-noise ratio of the measurements. For the SM measurements, a previously developed algorithm was used to reconstruct the IMS data. A new algorithm was developed for the DM analyses involving a two-step process that first recovers the SM data and then decodes the SM data. The algorithm also performs multiple refining procedures in order to minimize demultiplexing artifacts. The new IMS-Orbitrap MS platform was demonstrated by the analysis of proteomic and petroleum samples, where the integration of IMS and high mass resolution proved essential for accurate assignment of molecular formulae.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling Stochastic Variability in Multiband Time-series Data

In preparation for the era of time-domain astronomy with upcoming large-scale surveys, we propose a state-space representation of a multivariate damped random walk process as a tool to analyze irregularly-spaced multifilter light curves with heteroscedastic measurement errors. We adopt a computationally efficient and scalable Kalman filtering approach to evaluate the likelihood function, leading to maximum O(k 3 n) complexity, where k is the number of available bands and n is the number of unique observation times across the k bands. This is a significant computational advantage over a commonly used univariate Gaussian process that can stack up all multiband light curves in one vector with maximum O(k 3 n 3 ) complexity. Using such efficient likelihood computation, we provide both maximum likelihood estimates and Bayesian posterior samples of the model parameters. Three numerical illustrations are presented: (i) analyzing simulated five-band light curves for a comparison with independent single-band fits; (ii) analyzing five-band light curves of a quasar obtained from the Sloan Digital Sky Survey Stripe 82 to estimate short-term variability and timescale; (iii) analyzing gravitationally lensed g- and r-band light curves of Q0957+561 to infer the time delay. Two R packages, Rdrw and timedelay, are publicly available to fit the proposed models.

79 ASTRONOMY AND ASTROPHYSICS↗

Multi-Sensor Aerosol Products Sampling System

Global and local properties of atmospheric aerosols have been extensively observed and measured using both spaceborne and ground-based instruments, especially during the last decade. Unique properties retrieved by the different instruments contribute to an unprecedented availability of the most complete set of complimentary aerosol measurements ever acquired. However, some of these measurements remain underutilized, largely due to the complexities involved in analyzing them synergistically. To characterize the inconsistencies and bridge the gap that exists between the sensors, we have established a Multi-sensor Aerosol Products Sampling System (MAPSS), which consistently samples and generates the spatial statistics (mean, standard deviation, direction and rate of spatial variation, and spatial correlation coefficient) of aerosol products from multiple spacebome sensors, including MODIS (on Terra and Aqua), MISR, OMI, POLDER, CALIOP, and SeaWiFS. Samples of satellite aerosol products are extracted over Aerosol Robotic Network (AERONET) locations as well as over other locations of interest such as those with available ground-based aerosol observations. In this way, MAPSS enables a direct cross-characterization and data integration between Level-2 aerosol observations from multiple sensors. In addition, the available well-characterized co-located ground-based data provides the basis for the integrated validation of these products. This paper explains the sampling methodology and concepts used in MAPSS, and demonstrates specific examples of using MAPSS for an integrated analysis of multiple aerosol products.

Petrenko, M.↗

A statistical and simulation-informed model for estimating permeability from pore size distribution in saturated geomaterials

Accurate permeability estimation is essential across subsurface engineering applications but remains challenging due to the complex pore structures of natural geomaterials. Traditional empirical methods and simplified theoretical models often inadequately capture the role of pore size distribution and connectivity. Here, this study develops a statistical and simulation-informed permeability model that collapses pore-scale complexity into a compact scaling of the form k = αϕμ d 2 , where ϕ is porosity, μ d is mean pore size, and α is a weakly varying coefficient. By combining pore network simulations with statistical analysis of unimodal and bimodal pore size distributions, we identify three key findings: (i) permeability is much more sensitive to mean pore size than to porosity; (ii) across extensive datasets, the ratio σ d /μ d (standard deviation to mean) clusters around a characteristic value ∼0.4, allowing the effects of the full pore size distribution to be represented by μ d and a narrowly varying α ≈ 0.05; and (iii) for bimodal systems, there exists a critical fraction of small pores ∼0.78 above which flow becomes small-pore dominated, enabling the definition of an effective flow-controlling pore population and facilitating simplified permeability estimation for such systems. The resulting model, which requires only porosity and a representative mean pore size as inputs, is validated against comprehensive experimental datasets (>1700 samples) spanning diverse soils and rocks and achieves good predictive accuracy. Overall, this work provides a physically grounded yet practically simple permeability estimator suitable for subsurface engineering, environmental protection, and resource management applications.

Permeability↗

Effect of annealing on Curie temperature and phase transition in La{sub 0.55}Sr{sub 0.08}Mn{sub 0.37}O{sub 3} epitaxial films grown on SrTiO{sub 3} (100) substrates by reactive radio frequency magnetron sputtering

Mn-poor LaSrMnO{sub 3} (LSMO) epitaxial films were grown on SrTiO{sub 3} (100) substrates by radio frequency magnetron sputtering in an argon and oxygen gas mix, and then the samples were annealed in air at various temperatures (T{sub a}). 2 theta-chi X-ray diffraction mapping, nano-beam diffraction analysis through transmission electron microscopy, and electron back scatter diffraction through scanning electron microscopy revealed that the crystal symmetry of the LSMO films changed from monoclinic/orthorhombic to rhombohedral on annealing in air. Curie temperature (T{sub C}) of the LSMO films was found to increase with increasing T{sub a}, and become higher than the room temperature at T{sub a} ≥ 861 °C, indicating that the cause of these changes was the filling of oxygen and the transition of the crystal symmetry into rhombohedral. - Highlights: •Mn-poor LSMO changed from monoclinic/orthorhombic to rhombohedral by oxygen supply. •Mn-poor LSMO was increased T{sub C} by changed crystal symmetry, and it showed T{sub C} above RT. •Annealed in air effectively supplied O{sub 2} more than O{sub 2} gas during sputtering •EBSD is useful to evaluate crystal symmetry of complex oxide film from the substrate.

36 MATERIALS SCIENCE↗

Groundwater Monitoring Report, U.S. Department of Energy Y-12 National Security Complex, Oak Ridge, Tennessee

This report contains the groundwater and surface water monitoring data obtained during calendar year (CY) 2019 at the U.S. Department of Energy (DOE) Y-12 National Security Complex (Y-12) on the DOE Oak Ridge Reservation (ORR) in Oak Ridge, Tennessee. The monitoring data were obtained from wells, springs, and surface water sampling locations in three hydrogeologic regimes at Y-12. The Bear Creek Hydrogeologic Regime (Bear Creek Regime) encompasses a section of Bear Creek Valley (BCV) between the west end of Y-12 and the west end of the Bear Creek Watershed (directions are in reference to the Y-12 grid system, shown as Plant North. The Upper East Fork Poplar Creek Hydrogeologic Regime (East Fork Regime) encompasses the Y-12 industrial facilities and support structures in BCV. The Chestnut Ridge Hydrogeologic Regime (Chestnut Ridge Regime) encompasses a section of Chestnut Ridge directly south of Y-12. Background information in Section 2 of this report outlines the hydrogeologic framework for groundwater and surface water quality monitoring at Y-12 and includes an overview of the groundwater contamination in each hydrogeologic regime. Section 3 provides details regarding the groundwater and surface water sampling and analysis activities implemented under the Y-12 GWPP, including sampling locations and frequency, sample collection and handling, field measurements and laboratory analytes, quality assurance (QA)/quality control (QC) sampling, data management, and data quality assessment (DQA). However, the equivalent QA/QC or DQA information for the groundwater and surface water data associated with the monitoring programs implemented by UCOR are not included in this report and instead are deferred to referenced programmatic plans and reports issued by OREM and UCOR. Section 4 of this report presents a summary evaluation of the CY 2019 monitoring data with regard to the respective objectives of surveillance monitoring and exit pathway/perimeter monitoring. The evaluation is based primarily on the analytical results for the following principal groundwater contaminants at Y-12: nitrate, uranium, gross alpha activity, gross beta activity, and volatile organic compounds (VOCs). Section 5 summarizes the most significant findings with respect to the principal contaminants along with recommendations for any proposed changes to the ongoing groundwater and surface water quality monitoring performed under the Y-12 GWPP. Technical reports and plans cited in the narrative sections of the report are listed in Section 6. Narrative sections of this report reference several appendices. Figures (maps and diagrams) and data tables (excluding data summary tables incorporated in the narrative sections) are in Appendix A and Appendix B, respectively. Appendix C contains construction details for each well sampled during CY 2019 by either the Y-12 GWPP or UCOR, along with schematic diagrams for wells equipped with Westbay™ multiport sampling equipment or Barcad® pump systems. Appendix D supports the background summary discussion in Section 2 and provides more detailed information about the hydrogeologic framework for groundwater and surface water monitoring at Y-12, including the primary sources of groundwater contamination in each hydrogeologic regime. Results for all field measurements and laboratory analyses obtained by the Y-12 GWPP and UCOR are presented in Appendix E, which also includes the sample numbers for the QA/QC samples associated with groundwater and surface water monitoring performed by the Y-12 GWPP.

54 ENVIRONMENTAL SCIENCES↗

Evaluating the Utility of Permethylated Polysaccharide Solution NMR Data for Characterization of Insoluble Plant Cell Wall Polysaccharides

Plant cell wall polysaccharide analysis encompasses the utilization of a variety of analytical tools, including gas and liquid chromatography, mass spectrometry (MS), and nuclear magnetic resonance (NMR) spectroscopy. These methods provide complementary data, which enable confident structural proposals of the many complex polysaccharide structures that exist in the complex matrices of plant cell walls. However, cell walls contain fractions of varying solubilities, and a few techniques are available that can analyze all fractions simultaneously. We have discovered that permethylation affords the complete dissolution of both soluble and insoluble polysaccharide fractions of plant cell walls in organic solvents such as chloroform or acetonitrile, which can then be analyzed by a number of analytical techniques including MS and NMR. In this work, NMR structure analysis of 10 permethylated polysaccharide standards was undertaken to generate chemical shift data providing insights into spectral changes that result from permethylation of polysaccharide residues. This information is of especial relevance to the structure analysis of insoluble polysaccharide materials that otherwise are not easily investigated by solution-state NMR methodologies. The preassigned NMR chemical shift data is shown to be vital for NMR structure analysis of minor polysaccharide components of plant cell walls that are particularly difficult to assign by NMR correlation data alone. With the assigned chemical shift data, we analyzed the permethylated samples of destarched, alcohol-insoluble residues of switchgrass and poplar by two-dimensional NMR spectral profiling. Thus, we identified, in addition to the major polysaccharide components, two minor polysaccharides, namely, <5% 3-linked arabinoxylan (switchgrass) and <2% glucomannan (poplar). In particular, the position of the arabinose residue in the arabinoxylan of the switchgrass sample was confidently assigned based on chemical shift values, which are highly sensitive to local chemical environments. Furthermore, the high resolution afforded by the 1H NMR spectra of the permethylated switchgrass and poplar samples allowed facile relative quantitative analysis of their polysaccharide composition, utilizing only a few milligrams of the cell wall material. As a result, the concepts herein developed will thus facilitate NMR structure analysis of insoluble plant cell wall polysaccharides, more so of minor cell wall components that are especially challenging to analyze with current methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High Semi-Volatile Organic Aerosol Contributions Associated with Ammonium Nitrate in the Coastal Urban Environment

Organic components often contribute 50% or more of the submicron aerosol mass in coastal urban environments, but their partitioning between the gas and particle phases is controlled by a complex mixture of unidentified organic compounds that are poorly constrained by observations. This study compares daily filter organic functional groups (OFG) with online organic mass fragments from La Jolla, California, as part of the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE), quantifying the contributions of four types of non-volatile (NV) organic emission sources to the submicron composition. Daily filters retained only 0.79–0.98 µg/m3 NV submicron organic mass concentration, even though 1.8–1.9 µg/m3 non-refractory (NR) submicron organic mass concentration was measured online. The 62–64% of measured NR submicron organic mass concentration that exceeded what remained on the filters after 23-hr sampling is interpreted as semi-volatile, consistent with the moderate correlation of the NR NV difference to NR ammonium, NR nitrate, and biomass burning-related NV and NR organic factors. The association between semi-volatile organic components and ammonium nitrate likely results from both co-emission and co-evaporation. Size-resolved filter analysis showed that NV organic mass concentration accounts for 68% of NR organic mass concentration for 0.5–1 µm dry diameter but account for 9.0% for 0.18–0.5 µm dry diameter showing the differences in volatility between particle modes. Importantly, volatility of organic components was size-dependent Information Classification: General and associated with ammonium nitrate and biomass burning, providing guidance for constraining atmospheric aerosol properties in global models.

Pelayo, Christian↗

Irradiation of Sintered Mo Disks with Presence of Organic Impurities in He Gas Flow

Argonne National Laboratory (Argonne) is assisting NorthStar Medical Technologies in the development of a domestic supply of 99 Mo. Specifically, the present study focuses on the production of 99 Mo-feed solution used by the RadioGenix™ 99m Tc generator. During the target-irradiation phase of production, impurities can potentially be introduced into the feed, and can lead to disturbance of ligand- 99m Tc complexation chemistry and contamination of the final radiopharmaceutical that directly interacts with the patient. To address this issue, Argonne performed irradiations and chemical processing to identify whether the potential contamination of He flow with hydrocarbon oil during irradiation affects the radiochemical purity of the final K 2 MoO 4 (K 2 TcO 4 ) in 5M KOH solution. To mimic the conditions of real irradiation at NorthStar, Argonne used its electron linear accelerator and Van de Graaff facilities, heated the target and oil source to >800°C, and controlled oxygen in the presence of He during irradiation. Following irradiation, scanning electron microscopy (SEM) and carbon analysis (CA) were used to detect carbon contamination on the solid targets. The radiochemical purity of the dissolved targets was studied via thin-layer paper chromatography (TLC). As a result of these experiments, small regions of the surface of some irradiated disks were found to be high in carbon, but the total carbon content was still negligible in comparison to the reference sample, which had experienced no irradiation or contact with oil. The summarized results from the SEM, CA, and TLC tests lead to the conclusion that even an excess of oil and heating during irradiation do not affect the radiochemical purity of the final product.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Explaining and predicting human behavior and social dynamics in simulated virtual worlds: reproducibility, generalizability, and robustness of causal discovery methods

Ground Truth program was designed to evaluate social science modeling approaches using simulation test beds with ground truth intentionally and systematically embedded to understand and model complex Human Domain systems and their dynamics Lazer et al. (Science 369:1060–1062, 2020). Our multidisciplinary team of data scientists, statisticians, experts in Artificial Intelligence (AI) and visual analytics had a unique role on the program to investigate accuracy, reproducibility, generalizability, and robustness of the state-of-the-art (SOTA) causal structure learning approaches applied to fully observed and sampled simulated data across virtual worlds. In addition, we analyzed the feasibility of using machine learning models to predict future social behavior with and without causal knowledge explicitly embedded. In this paper, we first present our causal modeling approach to discover the causal structure of four virtual worlds produced by the simulation teams—Urban Life, Financial Governance, Disaster and Geopolitical Conflict. Our approach adapts the state-of-the-art causal discovery (including ensemble models), machine learning, data analytics, and visualization techniques to allow a human-machine team to reverse-engineer the true causal relations from sampled and fully observed data. We next present our reproducibility analysis of two research methods team’s performance using a range of causal discovery models applied to both sampled and fully observed data, and analyze their effectiveness and limitations. We further investigate the generalizability and robustness to sampling of the SOTA causal discovery approaches on additional simulated datasets with known ground truth. Our results reveal the limitations of existing causal modeling approaches when applied to large-scale, noisy, high-dimensional data with unobserved variables and unknown relationships between them. We show that the SOTA causal models explored in our experiments are not designed to take advantage from vasts amounts of data and have difficulty recovering ground truth when latent confounders are present; they do not generalize well across simulation scenarios and are not robust to sampling; they are vulnerable to data and modeling assumptions, and therefore, the results are hard to reproduce. Finally, when we outline lessons learned and provide recommendations to improve models for causal discovery and prediction of human social behavior from observational data, we highlight the importance of learning data to knowledge representations or transformations to improve causal discovery and describe the benefit of causal feature selection for predictive and prescriptive modeling.

97 MATHEMATICS AND COMPUTING↗

Electron Spin Resonance (ESR) studies of returned comet nucleus samples

The most important objective of the Comet Nucleus Sample Returm Mission is to return samples which could reflect formation conditions and evolutionary processes in the early solar nebula. It is expected that the returned samples will consist of fine-grained silicate materials mixed with ices composed of simple molecules such as H2O, NH3, CH4 as well as organics and/or more complex compounds. Because of the exposure to ionizing radiation from cosmic-ray, gamma-ray, and solar wind protons at low temperature, free radicals are expected to be formed and trapped in the solid ice matrices. The kind of trapped radical species together with their concentration and thermal stability can be used as a dosimeter as well as a geothermometer to determine thermal and radiation histories as well as outgassing and other possible alternation effects since the nucleus material was formed. Since free radicals that are known to contain unpaired electrons are all paramagnetic in nature, they can be readily detected and characterized in their native form by the Electron Spin Resonance (ESR) method. In fact, ESR has been shown to be a non-destructive, highly sensitive tool for the detection and characterization of paramagnetic, ferromagnetic, and radiation damage centers in terrestrial and extraterrestrial geological samples. The potential use of ESR as an effective method in the study of returned comet nucleus samples, in particular, in the analysis of fine-grained solid state icy samples is discussed.

Tsay, Fun-Dow↗

Development of Multiresolution Capabilities for the Holistic Energy Resource Optimization Network (HERON) tool A progress update

INL researchers work on technoeconomic analyses for integrated energy systems (IES) using the Framework for Optimization of ResourCes and Economics (FORCE). Within FORCE, researchers use the Holistic Energy Resource Optimization Network (HERON) tool to conduct optimization of grid portfolios under uncertain market conditions. These optimizations determine optimal capacities for all IES components and strategies for resource dispatch which maximize some economic metric (e.g., net present value). Resource dispatch occurs on finer timescales (typically hours) and thus are asked to respond to a given time series (e.g. hourly load demand profiles for a grid, or pre-determined electricity prices). Volatile and complex bidding dynamics as well as poorly forecasted weather events within deregulated markets add uncertainty to the time series; FORCE can address this uncertainty by training a reduced order model on historical time series and generate unique synthetic time series which represent individual scenarios or realizations of the market. The IES configuration can be simulated under these different sampled realizations and a stochastic optimization is conducted which optimizes the expected value of the desired economic metric. The training of a synthetic time series generator is limited by the chosen time resolution; dynamics can occur on different time scales. Seasonal demand trends can dominate faster dynamical events (such as power outages from certain sectors or severe weather events) which might not get captured correctly by the trained model. In this report, we investigate different ways of addressing the training and generation of time series on multiple time scales using three main algorithms: wavelet decomposition, dynamic mode decomposition, and generative adversarial networks for time series. We demonstrate a time series analysis that yields information on not just the frequency space but also temporal space: where a fast Fourier transform can provide what frequencies dominate, the new algorithms can provide when the frequencies dominate as well. These analyses can help improve IES optimization by allowing researchers to couple simulations at different timescales when it is most needed - seasonal, day-ahead, and real time optimization - with greater computational efficiency. Future work will include implementation of a subset of the proposed algorithms into the FORCE toolset and application of these analyses into multiple timescale optimization.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Optical vibrational spectroscopic signatures of ammonium diuranate process parameters

Ammonium diuranate (ADU) is commonly encountered in the nuclear fuel cycle; however, previous investigations have shown that ADU is a complex mixture of distinct compounds. Moreover, production parameters are known to heavily influence the composition of the resulting ADU. Here, we examine four samples of ADU prepared at Oak Ridge National Laboratory (ORNL), and one sample of ADU made at Pacific Northwest National Laboratory (PNNL), with the goal of further characterizing and elucidating the effect of processing parameters such as stir rate, strike direction, and temperature on material composition. Process parameters during ADU precipitation at ORNL and PNNL were well documented, and we relate process variables to optical vibrational spectroscopic signatures observed using Raman and infrared (IR) spectroscopy. In addition, powder X-ray diffraction (PXRD) reveals differences in the solid-phase composition of ADU precipitates, but we find that the primary phase is similar to the uranyl oxyhydroxyhydrate mineral metaschoepite. Despite the significant phase contributions of a metaschoepite-like phase, spectroscopic evidence of both nitrate and ammonium are observed for all samples. To gain a more holistic understanding of spectroscopic features of process parameters in ADU, principal component analysis (PCA) is employed and results in observable signatures that relate to the stir rate used during synthesis. These results provide further information about the process-dependence of ADU precipitate composition.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗