SEARCH · Engineering Papers

Results for “Parameter Identifiability”

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.

SourceResearch connectionFrom yearThrough year

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

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Regression-based reduced-order models to predict transient thermal output for enhanced geothermal systems

Reduced-order modeling is a promising approach, as many phenomena can be described by a few parameters/mechanisms. An advantage and attractive aspect of a reduced-order model is that it is computational inexpensive to evaluate when compared to running a high-fidelity numerical simulation. A reduced-order model takes couple of seconds to run on a laptop while a high-fidelity simulation may take couple of hours to run on a high-performance computing cluster. The goal of this paper is to assess the utility of regression-based reduced-order models (ROMs) developed from high-fidelity numerical simulations for predicting transient thermal power output for an enhanced geothermal reservoir while explicitly accounting for uncertainties in the subsurface system and site-specific details. Numerical simulations are performed based on equally spaced values in the specified range of model parameters. Key sensitive parameters are then identified from these simulations, which are fracture zone permeability, well/skin factor, bottom hole pressure, and injection flow rate. We found the fracture zone permeability to be the most sensitive parameter. The fracture zone permeability along with time, are used to build regression-based ROMs for the thermal power output. The ROMs are trained and validated using detailed physics-based numerical simulations. Finally, predictions from the ROMs are then compared with field data. We propose three different ROMs with different levels of model parsimony, each describing key and essential features of the power production curves. The coefficients in the proposed regression-based ROMs are developed by minimizing a non-linear least-squares misfit function using the Levenberg–Marquardt algorithm. The misfit function is based on the difference between numerical simulation data and reduced-order model. ROM-1 is constructed based on polynomials up to fourth order. ROM-1 is able to accurately reproduce the power output of numerical simulations for low values of permeabilities and certain features of the field-scale data. ROM-2 is a model with more analytical functions consisting of polynomials up to order eight, exponential functions and smooth approximations of Heaviside functions, and accurately describes the field-data. At higher permeabilities, ROM-2 reproduces numerical results better than ROM-1, however, there is a considerable deviation from numerical results at low fracture zone permeabilities. ROM-3 consists of polynomials up to order ten, and is developed by taking the best aspects of ROM-1 and ROM-2. ROM-1 is relatively parsimonious than ROM-2 and ROM-3, while ROM-2 overfits the data. ROM-3 on the other hand, provides a middle ground for model parsimony. Based on R 2 -values for training, validation, and prediction data sets we found that ROM-3 is better model than ROM-2 and ROM-1. For predicting thermal drawdown in EGS applications, where high fracture zone permeabilities (typically greater than 10 –15 m 2 ) are desired, ROM-2 and ROM-3 outperform ROM-1. As per computational time, all the ROMs are 10 4 times faster when compared to running a high-fidelity numerical simulation. In conclusion, this makes the proposed regression-based ROMs attractive for real-time EGS applications because they are fast and provide reasonably good predictions for thermal power output.

58 GEOSCIENCES↗

NASA NTRS2000

Effect of Discontinuities and Uncertainties on the Response and Failure of Composite Structures

The overall goal of this research was to assess the effect of discontinuities and uncertainties on the nonlinear response and failure of composite structures subjected to combined mechanical and thermal loads. The four key elements of the study were: (1) development of simple and efficient procedures for the accurate determination of transverse shear and transverse normal stresses in structural sandwiches as well as in unstiffened and stiffened composite panels and shells; (2) study the effects of transverse stresses on the response, damage initiation and propagation in composite and sandwich structures; (3) use of hierarchical sensitivity coefficients to identify the major parameters that affect the response and damage in each of the different levels in the hierarchy (micro-mechanical, layer, panel, subcomponent and component levels); and (4) application of fuzzy set techniques to identify the range and variation of possible responses. The computational models developed were used in conjunction with experiments, to understand the physical phenomena associated with the nonlinear response and failure of composite and sandwich structures. A toolkit was developed for use in conjunction with deterministic analysis programs to help the designer in assessing the effect of uncertainties in the different computational model parameters on the variability of the response quantities.

Noor, Ahmed K.↗

DOE OSTI2021

Sensitivity analysis of in-pile critical heat flux experiments in TREAT for characterization of RIA power-transient effects

A reactivity-initiated accident (RIA) is one type of postulated design basis accident (DBA) that can cause a departure from nucleate boiling (DNB) event in pressurized-water reactors (PWRs). A DNB occurrence and its consequences depend on the thermophysical properties of the fuel components and coolant, characteristics of the transient energy insertion into the fuel rod, and the onset of the critical heat flux (CHF) phenomenon. To leverage the restart of the Transient Reactor Test (TREAT) Facility, an effort is currently underway to better understand the cladding-to-coolant heat transfer mechanisms and the CHF phenomenon under fast-transient irradiation conditions. This paper characterizes the impact of power transients on the thermal-hydraulic behavior of a TREAT Facility reactor heater rodlet CHF experiment to provide the priority of parameters that need to be investigated for an improved CHF model. Sobol sensitivity analysis methods and the Reactor Excursion and Leak Analysis Program (RELAP5-3D) code were used to identify key input parameters on the uncertainty in the prediction of peak outer- and inner-surface temperatures of the heater tube, as well as the time of the DNB event. A series of sensitivity analyses revealed the total energy deposition on the tube and the transient effects of power pulse had large impacts on the maximum temperatures. The CHF multiplier had the largest impact on the time occurrence of CHF. The overall results show the energy deposition rate in the tube is the most influencing factor to the manifestation of CHF and the resulting thermal-hydraulic behaviors of the tube. The multiplier for the CHF, which is interpreted as the predicted CHF value, has the largest Sobol indices for the time of the CHF in all cases, since it directly determines the occurrence of CHF. It is inferred that the uncertainties in the thermal-hydraulic behaviors of fuels increase with respect to the key parameters as the power pulse becomes broader, and an accurate estimation of the energy deposition rate is required to reduce the uncertainty in the evaluation of the integrity of fuel if the CHF is expected to occur near the peak power. Therefore, the outputs are expected to provide rigorous interpretation of ongoing in-pile CHF experiments in the TREAT Facility reactor regarding thermal-hydraulic behavior of the fuel system aiming for a new transient in-pile CHF model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1987

Photometric properties of lunar terrains derived from Hapke's equation

Hapke's bidirectional reflectance equations provide the most rigorous available description of photometric behavior in terms of physically meaningful parameters. The primary objective of this study was to derive Hapke parameters for the lunar surface from both disk-integrated and disk-resolved photometric data. Hapke's equation was fit to the disk-integrated phase curves and disk-resolved data for dark, average, and bright terrain classes using an iterative, nonlinear least squares algorithm described by Helfenstein. Parameters were initially determined from the disk-integrated data, and the result was applied as a first guess to the iterative solution of parameters for individual terrain classes. Plots are presented of the disk-resolved data normalized to corresponding brightnesses predicted from the disk-integrated solution under the same illumination and viewing geometries. Systematic trends in disk-resolved parameters can be identified. Values for single scattering albedo (w) of the dominantly anorthositic average and bright terrains are significantly larger than the value for the basaltic dark terrains (mare). Values for surge brightness parameters, particle phase function, and average topographic slope angle are also discussed.

Helfenstein, Paul↗

NASA NTRS1990

Attitude determination and calibration using a recursive maximum likelihood-based adaptive Kalman filter

An adaptive Kalman filter design that utilizes recursive maximum likelihood parameter identification is discussed. At the center of this design is the Kalman filter itself, which has the responsibility for attitude determination. At the same time, the identification algorithm is continually identifying the system parameters. The approach is applicable to nonlinear, as well as linear systems. This adaptive Kalman filter design has much potential for real time implementation, especially considering the fast clock speeds, cache memory and internal RAM available today. The recursive maximum likelihood algorithm is discussed in detail, with special attention directed towards its unique matrix formulation. The procedure for using the algorithm is described along with comments on how this algorithm interacts with the Kalman filter.

Kelly, D. A.↗

NASA NTRS2017

Monte Carlo Computational Modeling of Atomic Oxygen Interactions

Computational modeling of the erosion of polymers caused by atomic oxygen in low Earth orbit (LEO) is useful for determining areas of concern for spacecraft environment durability. Successful modeling requires that the characteristics of the environment such as atomic oxygen energy distribution, flux, and angular distribution be properly represented in the model. Thus whether the atomic oxygen is arriving normal to or inclined to a surface and whether it arrives in a consistent direction or is sweeping across the surface such as in the case of polymeric solar array blankets is important to determine durability. When atomic oxygen impacts a polymer surface it can react removing a certain volume per incident atom (called the erosion yield), recombine, or be ejected as an active oxygen atom to potentially either react with other polymer atoms or exit into space. Scattered atoms can also have a lower energy as a result of partial or total thermal accommodation. Many solutions to polymer durability in LEO involve protective thin films of metal oxides such as SiO2 to prevent atomic oxygen erosion. Such protective films also have their own interaction characteristics. A Monte Carlo computational model has been developed which takes into account the various types of atomic oxygen arrival and how it reacts with a representative polymer (polyimide Kapton H) and how it reacts at defect sites in an oxide protective coating, such as SiO2 on that polymer. Although this model was initially intended to determine atomic oxygen erosion behavior at defect sites for the International Space Station solar arrays, it has been used to predict atomic oxygen erosion or oxidation behavior on many other spacecraft components including erosion of polymeric joints, durability of solar array blanket box covers, and scattering of atomic oxygen into telescopes and microwave cavities where oxidation of critical component surfaces can take place. The computational model is a two dimensional model which has the capability to tune the interactions of how the atomic oxygen reacts, scatters, or recombines on polymer or nonreactive surfaces. In addition to the specification of atomic oxygen arrival details, a total of 15 atomic oxygen interaction parameters have been identified as necessary to properly simulate observed interactions and resulting polymer erosion that have been observed in LEO. The tuning of the Monte Carlo model has been accomplished by adjusting interaction parameters so the erosion patterns produced by the model match those from several actual LEO space experiments. Surface texturing in LEO can also be predicted by the model. Such comparison of space tests with ground laboratory experiments have enabled confidence in ground laboratory lifetime prediction of protected polymers. Results of Monte Carlo tuning, examples of surface texturing and undercutting erosion prediction, and several examples of how the model can be used to predict other LEO and Mars orbital space results are presented.

Monte Carlo Modeling↗

NASA NTRS2017

Monte Carlo Computational Modeling of Atomic Oxygen Interactions

Atomic Oxygen↗

DOE OSTI2024

An open-source framework for balancing computational speed and fidelity in production cost models

Studies of bulk power system operations need to incorporate uncertainty and sensitivity analyses, especially around exposure to weather and climate variability and extremes, but this remains a computational modeling challenge. Commercial production cost models (PCMs) have shorter runtimes, but also important limitations (opacity, license restrictions) that do not fully support stochastic simulation. Open-source PCMs represent a potential solution. They allow for multiple, simultaneous runs in high-performance computing environments and offer flexibility in model parameterization. Yet, developers must balance computational speed (i.e. runtime) with model fidelity (i.e. accuracy). In this paper, we present Grid Operations (GO), a framework for instantiating open-source, scale-adaptive PCMs. GO allows users to search across parameter spaces to identify model versions that appropriately balance computational speed and fidelity based on experimental needs and resource limits. Results provide generalizable insights on how to navigate the fidelity and computational speed tradeoff through parameter selection. We show that models with coarser network topologies can accurately mimic market operations, sometimes better than higher-resolution models. It is thus possible to conduct large simulation experiments that characterize operational risks related to climate and weather extremes while maintaining sufficient model accuracy.

42 ENGINEERING↗

DOE OSTI2019

Sensitivity analysis of the effect of wind characteristics and turbine properties on wind turbine loads

Proper wind turbine design relies on the ability to accurately predict ultimate and fatigue loads of turbines. The load analysis process requires precise knowledge of the expected wind-inflow conditions as well as turbine structural and aerodynamic properties. However, uncertainty in most parameters is inevitable. It is therefore important to understand the impact such uncertainties have on the resulting loads. The goal of this work is to assess which input parameters have the greatest influence on turbine power, fatigue loads, and ultimate loads during normal turbine operation. An elementary effects sensitivity analysis is performed to identify the most sensitive parameters. Separate case studies are performed on (1) wind-inflow conditions and (2) turbine structural and aerodynamic properties, both cases using the National Renewable Energy Laboratory 5 MW baseline wind turbine. The Veers model was used to generate synthetic International Electrotechnical Commission (IEC) Kaimal turbulence spectrum inflow. The focus is on individual parameter sensitivity, though interactions between parameters are considered. The results of this work show that for wind-inflow conditions, turbulence in the primary wind direction and shear are the most sensitive parameters for turbine loads, which is expected. Secondary parameters of importance are identified as veer, u -direction integral length, and u components of the IEC coherence model, as well as the exponent. For the turbine properties, the most sensitive parameters are yaw misalignment and outboard lift coefficient distribution; secondary parameters of importance are inboard lift distribution, blade-twist distribution, and blade mass imbalance. This information can be used to help establish uncertainty bars around the predictions of engineering models during validation efforts, and provide insight to probabilistic design methods and site-suitability analyses.

17 WIND ENERGY↗

DOE OSTI2021

Validation of Jezebel Reactivity Coefficients and Sensitivity Analysis

Nuclear data validation is often performed today using criticality measurements. The gold standard for criticality measurements is the International Criticality Safety Benchmark Experiment Project (ICSBEP). The validation specifically focuses on the effective multiplication factor ($k_{eff}$). $K_{eff}$ is a relatively easy parameter to infer and has reduced uncertainty due to being at or above critical. However, while $k_{eff}$ is the most documented parameter and its uncertainties and sensitivities have been evaluated in great detail, it cannot be used as a standalone metric to determine inaccuracies in nuclear data (e.g., cross section data, PFNS, nu), which is based on theory, physics, and differential measurements. The Experiments Underpinned by Computational Learning for Improvements in nuclear Data (EUCLID) project aims to identify compensating errors in specific isotope nuclear data by optimally designing experiments that are, or are not sensitive to a suite of measurement parameters beyond $k_{eff}$. By identifying parameters that are sensitive to each other, oppositely sensitive, or have substantial magnitude differences in sensitivity, experiments can be designed to constrain questionable nuclear data. One sensitivity that is of particular interest to this project includes the sensitivity of reactivity coefficients. Reactivity coefficients compare reactivity, which is related to $k_{eff}$ at two different states therefore being sensitive to small changes in the system. The most common type of reactivity coefficient measurements is comparison to void for a small sample within the assembly. It is key that the sample sizes are small enough to not affect the flux of the full system. Reactivity coefficients were evaluated for many early experiments to better understand transport corrected cross sections. In fact, ICSBEP includes reactivity coefficient results as “Supplemental Measurements” in appendices for a handful of older benchmarks. One of those benchmarks is Jezebel, the bare Pu critical assembly. This paper compares new simulations of reactivity coefficients for Jezebel, and explores the sensitivity of reactivity coefficients to small changes in nuclear data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Model form and sensitivity analysis of CALPHAD-based nucleation models in b-stabilized Ti alloys

Accurate prediction of α-phase nucleation and growth in β-stabilized titanium alloys is crucial for designing heat treatments to optimize mechanical properties in additively manufactured lightweight components. Ideally, predictions of nucleation and growth would incorporate both top-down observations of past experimental heat treatments and bottom-up modeling of phase transformations; however, the appropriate method of combining these information sources is not self-evident. Combining top-down and bottom-up information requires a unified form of model that can connect between spatiotemporal scales, as well as sets of fitting parameters that can be identified by each data source. The selection of which parameters to fit to which data source can be made based on expert opinion, or by performing a sensitivity analysis. In solid-solid nucleation, direct observation of the nucleation and growth process is challenging. Most data on the heat treatment-controlled phase transformations are not in-situ. To predict the process and outcome of the nucleation, growth and coarsening of precipitates, theoretical models of the nucleation pathway are used to bridge the gap. Many sources of uncertainty affect the modeling of this nucleation process. It can be influenced by small variations in the thermomechanical processing history, chemical composition, and initial microstructure. If molecular dynamics (MD) simulations are used to determine thermodynamic quantities and inform CALPHAD modeling, additional uncertainty can be introduced and accounted for using Bayesian methods. Top-down uncertainties require additional steps to quantify. The influence of nucleation model form on the sensitivity of predictions to input parameters and physical conditions is the focus of this study. Classical nucleation theory (CNT) allows modeling to formulate the nucleation as homogeneous or, more commonly, heterogeneous. Non-classical nucleation models are also increasingly explored as a means of reconciling top-down and bottom-up data. In this study, the sensitivity of the intragranular nucleation of α in a β-annealed, slow-cooled aging (BASCA) heat treatment of β-stabilized Ti5553 alloy is explored using CNT and both heterogeneous and homogeneous assumptions. The Kampmann-Wagner Numerical model of precipitate nucleation and growth is employed. Using open-source tools (pyCalphad and thermodynamic modeling of TiMo as a surrogate system, a sensitivity analysis is performed to measure variations in key parameters, including chemical driving force, interfacial energy, and diffusivity, as they relate to predictions of precipitate number density. The inclusion of top-down and bottom-up data in selection of nucleation model form is discussed.

Rodriguez Negron, A. M.↗

DOE OSTI2015

SU-E-T-249: Determining the Sensitivity of Beam Profile Parameters for Detecting Energy Changes in Flattening Filter-Free Beams

Purpose: To identify the beam profile parameters that can be used to detect energy changes in a flattening filter-free photon beams. Methods: Flattening filter-free beam profiles (inline, crossline, and diagonals) were measured for multiple field sizes (25×25cm and 10×10cm) at 6MV on a clinical system (Truebeam, Varian Medical Systems Palo Alto CA). Profiles were acquired for baseline energy and detuned beams by changing the bending magnet current (BMC), above and below baseline. The following profile parameters were measured: flatness (off-axis ratio at 80% of field size), symmetry, uniformity, slope, and the off-axis ratio (OAR) at several off-axis distances. Tolerance values were determined from repeated measurements. Each parameter was evaluated for sensitivity to the induced beam changes, and the minimum detectable BMC change was calculated for each parameter by calculating the change in BMC that would Result in a change in the parameter above the measurement tolerance. Results: Tolerance values for the parameters were-Flatness≤0.1%; Symmetry≤0.4%; Uniformity≤0.01%; Slope≤ 0.001%/mm. The measurements made with a field size of 25cm and a depth of d=1.5cm showed the greatest sensitivity to bending magnet current variations. Uniformity had the highest sensitivity, able to detect a change in BMC of BMC=0.02A. The OARs and slope were sensitive to the magnitude and direction of BMC change. The sensitivity in the flatness parameter was BMC=0.04A; slope was sensitive to BMC=0.05A. The sensitivity decreased for OARs measured closer to central axis-BMC(8cm)=0.23A; BMC(5cm)=0.47A; BMC(2cm)=1.35A. Symmetry was not sensitive to changes in BMC. Conclusion: These tests allow for better QA of FFF beams by setting tolerance levels to beam parameter baseline values which reflect variations in machine calibration. Uniformity is most sensitive to BMC changes, while OARs provide information about magnitude and direction of miscalibration. Research funding provided by Varian Medical Systems. Dr. Sasa Mutic receives compensation for providing patient safety training services from Varian Medical Systems, the sponsor of this study.

60 APPLIED LIFE SCIENCES↗

NASA NTRS2000

Uncertainty Analysis of Stiffened Composite Panels

A study is made of the variability in the nonlinear response of three stiffened composite panels associated with variations in their geometric and material parameters. The three panels have a cylindrical skin with either four or five T-shaped stiffeners. Two of the panels have a notch and the third panel has a circular cutout. Hierarchical sensitivity analysis is used to identify the major parameters at the micro-mechanical, layer, laminate and sub-component levels. The major parameters are then taken to be fuzzy parameters, and a fuzzy set analysis is used to determine the range of variation of the response associated with pre-selected variations in the major parameters.

Noor, Ahmed K.↗

NASA NTRS2006

Estimation of Stability and Control Derivatives of an F-15

A technique for real-time estimation of stability and control derivatives (derivatives of moment coefficients with respect to control-surface deflection angles) was used to support a flight demonstration of a concept of an indirect-adaptive intelligent flight control system (IFCS). Traditionally, parameter identification, including estimation of stability and control derivatives, is done post-flight. However, for the indirect-adaptive IFCS concept, parameter identification is required during flight so that the system can modify control laws for a damaged aircraft. The flight demonstration was carried out on a highly modified F-15 airplane (see Figure 1). The main objective was to estimate the stability and control derivatives of the airplane in nearly real time. A secondary goal was to develop a system to automatically assess the quality of the results, so as to be able to tell a learning neural network which data to use. Parameter estimation was performed by use of Fourier-transform regression (FTR) a technique developed at NASA Langley Research Center. FTR is an equation- error technique that operates in the frequency domain. Data are put into the frequency domain by use of a recursive Fourier transform for a discrete frequency set. This calculation simplifies many subsequent calculations, removes biases, and automatically filters out data beyond the chosen frequency range. FTR as applied here was tailored to work with pilot inputs, which produce correlated surface positions that prevent accurate parameter estimates, by replacing half the derivatives with predicted values. FTR was also set up to work only on a recent window of data, to accommodate changes in flight condition. A system of confidence measures was developed to identify quality-parameter estimates that a learning neural network could use. This system judged the estimates primarily on the basis of their estimated variances and of the level of aircraft response. The resulting FTR system was implemented in the Simulink software system and auto-coded in the C programming language for use on the Airborne Research Test System (ARTS II) computer installed in the F-15 airplane. The Simulink model was also used in a control room that utilizes the Ring Buffered Network Bus hardware and software, making it possible to evaluate test points during flights. In-flight parameter estimation was done for piloted and automated maneuvers, primarily at three test conditions. Figure 2 shows results for pitching moment due to symmetric stabilator actuations for a series of three pitch doublet maneuvers (in a doublet maneuver, a command to change attitude in a given direction by a given amount is followed immediately by a command to change attitude in the opposite direction by the same amount). A time window of 5 seconds was used. The portions of the curves shown in red are those that passed the confidence tests. The technique showed good convergence for most derivatives for both kinds of maneuvers - typically within a few seconds. The confidence tests were marginally successful, and it would be necessary to refine them for use in an IFCS.

Smith, Mark↗

DOE OSTI2024

Parameter development and characterization of laser powder directed energy deposition of Nb – Alloy C103 for thin wall geometries

This work focuses on the parameter development and microstructural characterization of Nb-based alloy C103 for thin wall structures produced via laser powder – directed energy deposition (LP – DED). Laser power and scanning speeds were varied as part of a design of experiments to identify adequate print parameters. Combinations were evaluated for relative density, porosity, and geometrical accuracy. A combination of a laser power of 1420 W and scanning speed of 14 mm/s resulted in a relative density >99%, exhibited a consistent weld bead profile and was used for further microstructural evaluation. A mix of optical and scanning electron microscopy (SEM) revealed small, slightly elongated grains along the edges and large epitaxial grains in the central region along the Build – Transverse view. Scanning electron microscopy (SEM) and energy dispersive spectroscopy (EDS) analysis revealed Hf rich columnar cells in the center that transition into an evenly spaced cellular structure towards the edges of the built sample. Electron backscatter diffraction (EBSD) scans show a sharp [001] texture when looking at the cross section of the sample along the build direction. The Build – Scan view revealed a zig – zag pattern that follows the back – and – forth deposition strategy that was used. Finally, microhardness measurements were taken in the as – built (AB) and stress relieved (SR) conditions to baseline preliminary mechanical properties. The AB condition exhibited a large amount of scatter in the data and averages up to 11% larger than the SR condition. The reduction in scatter upon applying the SR cycle are indicative of a large concertation of dislocations present in the AB condition.

36 MATERIALS SCIENCE↗

NASA NTRS2017

Using Advanced Mass Spectrometry Techniques to Fully Characterize Atmospheric Organic Carbon: Current Capabilities and Remaining Gaps

Organic compounds in the atmosphere vary widely in their molecular composition and chemical properties, so no single instrument can reasonably measure the entire range of ambient compounds. Over the past decade, a new generation of in-situ, field-deployable mass spectrometers has dramatically improved our ability to detect, identify, and quantify these organic compounds, but no systematic approach has been developed to assess the extent to which currently available tools capture the entire space of chemical identity and properties that is expected in the atmosphere. Reduced-parameter frameworks that have been developed to describe atmospheric mixtures are exploited here to characterize the range of chemical properties accessed by a suite of instruments. Multiple chemical spaces (e.g. oxidation state of carbon vs. volatility, and oxygen number vs. carbon number) were populated with ions measured by several mass spectrometers, with gas- and particle-phase α-pinene oxidation products serving as the test mixture of organic compounds. Few gaps are observed in the coverage of the parameter spaces by the instruments employed in this work, though the full extent to which comprehensive measurement was achieved is difficult to assess due to uncertainty in the composition of the mixture. Overlaps between individual ions and regions in parameter space were identified, both between gas- and particle-phase measurements, and within each phase. These overlaps were conservatively found to account for little (<10%) of the measured mass. However, challenges in identifying overlaps and in accurately converting molecular formulas into chemical properties (such as volatility or reactivity) highlight a continued need to incorporate structural information into atmospheric measurements.

G. Isaacman-VanWertz↗

DOE OSTI2018

Analyzing the impact of reaction models on the production of hydrocarbons from thermally upgraded oil shales

Reaction parameters significantly affect oil production from shales by means of heating and in-situ upgrading. In this work, we perform numerical simulations of two chemical reaction models, which are mainly used in the research of kerogen pyrolysis and subsequent hydrocarbon decomposition in organic-rich porous media. They are the Braun and Burnham model and Wellington model. In these forward numerical simulations, we present the influence of the two reaction models on hydrocarbon production. The Braun and Burnham reaction model shows more vigorous kerogen and subsequent decomposition reactions and more hydrocarbon production than the Wellington model. A local sensitivity analysis identifies the reaction parameters with the highest influence on productivity, and the most sensitive outputs. A data-worth analysis identifies the most valuable observation data to be measured for the best prediction of total hydrocarbon production. We find that the most valuable observation data is the cumulative production of heavy oil in the Braun and Burnham model and of light oil in the Wellington model, respectively. Finally, once we determine the maximum allowable prediction uncertainty and the expected measurement uncertainty, the observation data to be measured for the minimization of prediction uncertainty can be obtained.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2021

Intelligent Agents for the Optimization of Atomic Layer Deposition

Atomic layer deposition (ALD) is a highly controllable thin film synthesis approach with applications in computing, energy, and separations. The flexibility of ALD means that it can access a massive chemical catalogue; however, this chemical and process diversity results in significant challenges in determining processing parameters that result in stable and uniform film growth with minimal precursor consumption. In situ measurements of the ALD growth per cycle (GPC) can accelerate process development but it still requires expert intuition and time-consuming trial and error to identify acceptable processing parameters. This procedure is made more difficult by the presence of experimental noise in the GPC values and the complexity of ALD surface chemistries. Here, a need exists for efficient optimization approaches capable of autonomously determining processing conditions resulting in optimal ALD film growth. In this work, we present the development of three optimization strategies and compare their performance in optimizing four simulated ALD processes. Furthermore, the effect of noise in the GPC measurements on optimization convergence is studied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗