Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analytical sensitivity 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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Analysis of secondary organic aerosol simulation bias in the Community Earth System Model (CESM2.1)

Organic aerosol (OA) has been considered as one of the most important uncertainties in climate modeling due to the complexity in presenting its chemical production and depletion mechanisms. To better understand the capability of climate models and probe into the associated uncertainties in simulating OA, we evaluate the Community Earth System Model version 2.1 (CESM2.1) configured with the Community Atmosphere Model version 6 (CAM6) with comprehensive tropospheric and stratospheric chemistry representation (CAM6-Chem) through a long-term simulation (1988–2019) with observations collected from multiple datasets in the United States. We find that CESM generally reproduces the interannual variation and seasonal cycle of OA mass concentration at surface layer with a correlation of 0.40 compared to ground observations and systematically overestimates (69%) in summer and underestimates (-19%) in winter. Through a series of sensitivity simulations, we reveal that modeling bias is primarily related to the dominant fraction of monoterpene-formed secondary organic aerosol (SOA), and a strong positive correlation of 0.67 is found between monoterpene emission and modeling bias in the eastern US during summer. In terms of vertical profile, the model prominently underestimates OA and monoterpene concentrations by 37%–99% and 82%–99%, respectively, in the upper air (>500 m) as validated against aircraft observations. Our study suggests that the current volatility basis set (VBS) scheme applied in CESM might be parameterized with monoterpene SOA yields that are too high, which subsequently results in strong SOA production near the emission source area. We also find that the model has difficulty in reproducing the decreasing trend of surface OA in the southeastern US probably because of employing pure gas VBS to represent isoprene SOA which is in reality mainly formed through multiphase chemistry; thus, the influence of aerosol acidity and sulfate particle change on isoprene SOA formation has not been fully considered in the model. This study reveals the urgent need to improve the SOA modeling in climate models.

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

Version 4 CALIPSO Imaging Infrared Radiometer ice and liquid water cloud microphysical properties – Part I: The retrieval algorithms

Following the release of the version 4 Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) data products from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) mission, a new version (version 4; V4) of the CALIPSO Imaging Infrared Radiometer (IIR) Level 2 data products has been developed. The IIR Level 2 data products include cloud effective emissivities and cloud microphysical properties such as effective diameter and ice or liquid water path estimates. Dedicated retrievals for water clouds were added in V4, taking advantage of the high sensitivity of the IIR retrieval technique to small particle sizes. This paper (Part I) describes the improvements in the V4 algorithms compared to those used in the version 3 (V3) release, while results will be presented in a companion (Part II) paper. The IIR Level 2 algorithm has been modified in the V4 data release to improve the accuracy of the retrievals in clouds of very small (close to 0) and very large (close to 1) effective emissivities. To reduce biases at very small emissivities that were made evident in V3, the radiative transfer model used to compute clear-sky brightness temperatures over oceans has been updated and tuned for the simulations using Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) data to match IIR observations in clear-sky conditions. Furthermore, the clear-sky mask has been refined compared to V3 by taking advantage of additional information now available in the V4 CALIOP 5 km layer products used as an input to the IIR algorithm. After sea surface emissivity adjustments, observed and computed brightness temperatures differ by less than ±0.2 K at night for the three IIR channels centered at 08.65, 10.6, and 12.05 µm, and inter-channel biases are reduced from several tens of Kelvin in V3 to less than 0.1 K in V4. We have also improved retrievals in ice clouds having large emissivity by refining the determination of the radiative temperature needed for emissivity computation. The initial V3 estimate, namely the cloud centroid temperature derived from CALIOP, is corrected using a parameterized function of temperature difference between cloud base and top altitudes, cloud absorption optical depth, and CALIOP multiple scattering correction factor. As shown in Part II, this improvement reduces the low biases at large optical depths that were seen in V3 and increases the number of retrievals. As in V3, the IIR microphysical retrievals use the concept of microphysical indices applied to the pairs of IIR channels at 12.05 and 10.6 µm and at 12.05 and 08.65 µm. The V4 algorithm uses ice look-up tables (LUTs) built using two ice habit models from the recent “TAMUice2016” database, namely the single-hexagonal-column model and the eight-element column aggregate model, from which bulk properties are synthesized using a gamma size distribution. Four sets of effective diameters derived from a second approach are also reported in V4. Here, the LUTs are analytical functions relating microphysical index applied to IIR channels 12.05 and 10.6 µm and effective diameter as derived from in situ measurements at tropical and midlatitudes during the Tropical Composition, Cloud, and Climate Coupling (TC4) and Small Particles in Cirrus Science and Operations Plan (SPARTICUS) field experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Kinetic Modeling Analysis of Ar Addition to Atmospheric Pressure N 2 –H 2 Plasma for Plasma-Assisted Catalytic Synthesis of NH 3

Zero-dimensional kinetic modeling of atmospheric pressure Ar–N 2 –H 2 nonthermal plasma was carried out to gain mechanistic insights into plasma-assisted catalytic synthesis of ammonia. Ar dilution is a common technique for tailoring plasma discharge properties and has been shown to enhance NH 3 formation when added to N 2 –H 2 plasma. The kinetic model was developed for a coaxial dielectric barrier discharge quartz wool-packed bed reactor operating at near room temperature using a kHz-frequency plasma source. With 30% Ar mixed in a 1:1 N 2 –H 2 plasma at 760 Torr, we find that NH 3 production is dominated by Eley–Rideal (E-R) surface reactions, which heavily involve surface NH x species derived from N and H radicals in the gas phase, while the influence of excited N 2 molecules is negligible. This is contrary to the commonly proposed mechanism that excited N 2 molecules created by Penning excitation of N 2 by Ar(4s) and Ar(4p) play a significant role in assisting NH 3 formation. Our model shows that the enhanced NH 3 formation upon Ar dilution is unlikely due to the interactions between Ar and H species, as excited Ar atoms have a weak effect on H radical formation through H 2 dissociation compared to electrons. We find that excited Ar atoms contribute to 28% of the N radical production in the gas phase via N 2 dissociation, while the rest are dominated by electron-impact dissociation. Furthermore, Ar species play a negligible role in the product NH 3 dissociation. N 2 conversion sensitivity analyses were carried out for electron number density (ne) and reduced electric field (E/N), and contributions from Ar to gas-phase N radical production were quantified. Finally, the model can provide guidance on potential reasons for observing enhanced NH 3 formation upon Ar dilution in N 2 –H 2 plasma beyond changes in the discharge characteristics.

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

The Interplay of Pauli Repulsion, Electrostatics, and Field Inhomogeneity for Blueshifting and Redshifting Vibrational Probe Molecules

Many molecules' vibrational frequencies are sensitive to intermolecular electric fields, enabling them to probe the field in complex molecular environments. However, it is often unclear whether the probe is responding to the local electric field or other types of intermolecular interactions, inhibiting interpretation of the frequency and effectiveness as probes. This is especially true for molecules whose vibrational frequencies blueshift instead of the more typical redshift in hydrogen bonding configurations. In this work, we computationally investigate the causes of redshifting versus blueshifting over a range of vibrational reporters. First, we apply adiabatic energy decomposition analysis to a paradigmatic set of probes, finding that redshifting only occurs when electrostatic interactions are strong enough to overcome the dominant and large blueshifting contribution of Pauli repulsion. Furthermore, we demonstrate that field inhomogeneity can further shift the frequency of many probes substantially to either reinforce or counteract the shift expected from a homogeneous field. We find that redshifting is reinforced by electric field inhomogeneity, otherwise field inhomogeneity further weakens the electrostatic contribution relative to Pauli repulsion, leading to blueshifting. Further calculations indicate that the probe's response to field inhomogeneity can be understood by considering the mass of the atoms involved in the stretching mode and sign of the electric field. In explaining the interplay of different intermolecular interactions and field inhomogeneity for many probes, our results should enable the use and interpretation of spectroscopic probes and their connection to electric fields in more complex systems.

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

DOE Award No.: DE-FE0023919 Phase 4 Scientific/Technical Report

This is the Phase 4 Report for the ‘Deepwater Methane Hydrate Characterization and Scientific Assessment or Genesis of Methane Hydrates in the Gulf of Mexico (GOM2)’ research project (DOE Award No. DE-FE0023919). The report summarizes activities from October 1, 2019 to September 30, 2020. The project is led by the University of Texas at Austin (UT). The project objective is to gain insight into the nature, formation, occurrence and physical properties of methane hydrate-bearing sediments for the purpose of methane hydrate resource appraisal through the planning and execution of drilling, coring, logging, testing and analytical activities that assess the geologic occurrence, regional context, and characteristics of marine methane hydrate deposits in the Gulf of Mexico outer continental shelf (OCS). We published a dedicated American Association of Petroleum Geologists Bulletin volume describing initial results from the UT-GOM2-1 expedition in Sept. 2020. This is part 1 of a multi-volume commitment by AAPG to this project. We further confirmed that the natural gas in hydrate at GC-955 was formed by primary microbial processes (>76.1 %). The in-situ effective permeability hydrate-bearing sandy silts at the GC-955 reservoir ranges from 0.1 md (1.0×10-16 m2) to 2.4 md (2.4×10-15m2) in cores with 83% to 93% hydrate saturation. The intrinsic permeability (the single phase permeability) is estimated from reconstituted samples to be ~12 md (1.2×10-14 m2) to ~41 md (4.1×10-14 m2). We used observation and models to interpret that the core degradation that is found in pressure cores is due to dissociation of the methane hydrate in the outer circumference of the core and dissolution of that methane into the fresh pore water that the core is stored with. We are designing approaches to minimize this core loss in the future. We spent an enormous amount of effort to further improve the ability of the pressure coring tool (the PCTB) to pressure seal correctly. We completed upgrading the upper section of the PCTB to address poor pressure. We successful tested the modifications at Geotek’s test facility in Salt Lake City (Bench Test II). We completed a Land Test of the PCTB at the Schlumberger Cameron Test and Training Facility (CTTF). The tool did not seal in 6 out of 7 tests and we clearly demonstrated that cuttings were wedging in the ball valve assembly, keeping the ball valve from sealing. We reproduced the failure mechanism observed during the land test at Salt Lake City and confirmed the sensitivity of the ball valve assembly to grit. Geotek designed and tested 9 modifications to address this issue and the PCTB is now 100% successfully sealing in the presence of grit. Our science expedition is scheduled for spring 2022 and we are fully focused on preparing for this. UT and Ohio State completed a Shallow Hazard Assessment report for each proposed UT-GOM2-2 drilling location, pursuant to 30 CFR 250.214(f) and 250.244 (f). The Shallow Hazard Reports will accompany the UT-GOM2-2 Exploration Plan that is submitted to BOEM, and completes the geological and geophysical analysis for UT-GOM2-2 permitting efforts. We updated the UT-GOM2-2 Operations Plan (Version 1). We completed the UT-GOM2-2 Science and Sample Distribution Plan (Version 1). We evaluated the scope, budget, and schedule that would result from using a commercial vessel. We developed detailed drilling schedule, mud volume, and resource estimates. We developed a vessel specification document and a well plan, and sent these documents to prospective vessel contractors.

03 NATURAL GAS↗
DOE OSTI2026

Titanium alloy response sensitivity to variations in spectral reconstructions of National Ignition Facility xenon line-emission x-ray sources

Thermomechanical shock experiments on the National Ignition Facility (NIF) aim to study high strain rate dynamic material response. In such experiments, the NIF laser is used to generate high fluence x-ray emission sources, which irradiate material samples of interest. Under sufficiently high x-ray energy deposition, thermomechanical impulses are generated in the materials. While it is known that the characteristics of x-ray generated impulses vary as a function of incident x-ray spectra, it remains unclear how spectral assumptions and uncertainties in NIF spectral reconstructions affect our interpretation of impulsive loading. Here, in this paper, we simulate the response of a standard titanium alloy baseline sample to synthetic analytically derived and measured NIF xenon line-emission x-ray sources with a radiation hydrodynamics code. We vary the source spectral characteristics based on different source reconstruction techniques to understand the resulting variation in baseline sample response and compare the simulated response with experimental results. We find that the response is highly sensitive to assumptions made about the spectral contents and that knowledge of spectral uncertainties bounds our understanding of the resulting material response. The results of this effort help to extend our ability to use baseline material samples to extract quantitative properties from x-ray experiments on the NIF.

Alloys↗
DOE OSTI2025

Synergistic effects of Pd single atoms and nanoclusters boosting SnO 2 gas sensing performance

Tin(IV) oxide-supported Pd is a promising heterogenous catalyst for CO oxidation relevant for environmental cleanup reactions. In this study, an atomically dispersed Pd catalyst on SnO 2 (ADC Pd/SnO 2 ) hybrid material is successfully synthesized via a straightforward wet chemistry method and is found to exhibit superior performance toward CO sensing. Ex situ EXAFS analysis confirms the formation of single Pd atoms and small Pd nanoclusters stabilized on the SnO 2 (110) surface. Further, the material exhibits high efficiency in generating adsorbed O 2 - as well as high activity in catalyzing CO oxidation at low temperatures, resulting in exceptional sensitivity and selectivity toward CO in comparison to pure SnO 2 and Pd nanoparticles loaded on SnO 2 respectively. In situ FTIR measurements unravel CO adsorption kinetics on ADC Pd/SnO 2 under reaction conditions, and a possible sensing mechanism is put forth in which CO is transformed into CO 2 by reaction with active oxygen species; and concurrently, carbon-related species (bicarbonates and carbonates) are formed and decomposed into CO 2 .

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

Systematically derived thermodynamic properties for alkane oxidation

Key combustion properties, such as ignition delay time, show strong sensitivity to the thermochemistry of the main species in the standard radical oxidation pathway at low temperatures (600-1000K). Significant uncertainties persist in current estimates of thermodynamic properties, particularly for the larger species of relevance to common practical fuels. In this article, we use advanced computational schemes to evaluate thermodynamic properties for the fuel, fuel radical, peroxy, hydroperoxy-alkyl radical, and hydroperoxide species for a set of 17 fuels containing up to 9 carbon atoms, and with various degrees of branching in the alkane backbone. The procedure, termed STAR-1D, combines conformer sampling to find the minimum geometry, B2PLYP-D3/cc-pVTZ harmonic frequency evaluations, and omega B97X-D/cc-pVTZ one-dimensional torsional mapping. It includes two physically-based scaling routines: a frequency-dependent scaling to B2PLYP-D3/cc-pVTZ anharmonic frequencies and a scaling of one-dimensional omega B97X-D/ccpVTZ torsional profiles to reproduce the product of the B2PLYP-D3/cc-pVTZ frequencies in the harmonic limit. Substantive comparisons with existing experimental databases, together with careful examinations of key theoretical assumptions, are used to explore the accuracy of the predictions. These computationally intensive explorations of 195 species were facilitated by automated thermochemistry software. Further, high accuracy 0 K heats of formation from a separate study are used with the STAR-1D computations to generate NASA polynomial representations. In addition to their intrinsic value, the present results also provide a reliable database for the optimization of group additivity or machine learning schemes for scaling to larger combustion systems. Towards this end a complementary extensive conformational analysis is carried out for the medium sized species and the thermodynamic properties of the lowest energy hydrogen-bonded and non-hydrogen bonded conformers are contrasted for larger hydroperoxy-alkyl radical species.

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

Burning velocities of R-32/O 2 /N 2 mixtures: Experimental measurements and development of a validated detailed chemical kinetic model

This work entails characterizing the flammability of the refrigerant R-32 (CH 2 F 2 ) by both experimental measurements and modeling. Burning velocities S u were measured using a constant-volume spherical-flame method for R-32/O2/N 2 mixtures with O 2 /N 2 ratios ranging from 21% (synthetic air) to 40%, pressures of (1 to 3) bar, and equivalence ratios $\phi$ of (0.8 to 1.3). Based on a critical assessment of available data, and extended by our own calculations, a detailed chemical kinetic model was developed and key reactions determined using reaction path and sensitivity analyses. Initiation and combustion were identified as distinct kinetic regimes and burning velocities were found to be controlled by two primary reactions: unimolecular decomposition of CH 2 F 2 → CHF + HF and the subsequent reaction, CHF + O 2 → CHFO + O, the latter reaction initiating the radical chain propagating and branching by producing O atoms. Sensitive rate constants in the kinetic model were critically adjusted within their uncertainties and current knowledge bounds to best fit the experimental burning velocities. We found that rate constants in the model could be adjusted to match a given experimental S u for specific conditions (O 2 loading, P, T, $\phi$). This, however, then fixes predicted burning velocities for other all conditions within (3 to 4)% if physically realistic rate parameters are maintained. Thus, the entire set of experimental data is fit, not just to particular conditions. Relative random uncertainties in the experimental Su measurements were (4 to 6)%, but assumptions made for thermal radiation lost by the burned gas in the spherical-flame experiments add an additional systematic uncertainty. Systematic differences between the limiting cases of adiabatic (no thermal radiation lost) and optically-thin (all thermal radiation lost) varied significantly with conditions and ranged from (4 to 30)% at high to low velocities, respectively, translating into uncertainties of (2 to 15)% considering the average of two limiting cases. Comparison of experimental and kinetically modeled Su values suggests that the burned gas tends towards the optically-thin limit at the lowest pressures and fuel loadings and toward the adiabatic limit at the highest pressures and loadings. We tested and found support for this conclusion with a detailed analysis as a function of all the conditions (T, P, % O 2 , $\phi$). This behavior appears to transition from optically-thin to adiabatic as the density of the initial fuel increases, which results in increased CO 2 in the burned gas and thus increased absorption of the thermal radiation (consistent with the Beer-Lambert Law). The validated detailed model based on evaluated kinetics is shown to accurately predict burning velocities for R-32 O 2 /N 2 mixtures over a wide range of conditions and provides a reliable basis for extrapolation to other conditions.

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

Uniformly Ordered Binary Decision Algorithm for Benchmark Experiment Correlations in Whisper Validation

When performing a validation exercise for determining the upper subcritical limit of a nuclear criticality safety application, an analyst should select and perform a statistical analysis on a population of benchmark experiments that are neutronically similar to the application. The size of this population should be sufficiently large such that the statistical analysis has a high degree of confidence that the bias plus bias uncertainty (calculational margin) has been accurately quantified. A complication arises because many benchmark experiments share common components, leading to correlations in their measured effective multiplication factors. Correlations between benchmark experiments within the population reduces its predictive power. This motivates the need for methods that consider benchmark experiment correlations and ensure adequate statistical significance of results. The Whisper code is a statistical analysis pack- age that incorporates nuclear data sensitivity coefficients from MCNP to assess benchmark experiment similarity and then performs an extreme-value analysis to estimate the bias plus bias uncertainty. The original methodology in Whisper does not consider the effect of benchmark experiment correlations when making this estimation, and this summary proposes the uniformly ordered binary decision algorithm to address this shortcoming. The original methodology in Whisper computes similarity coefficients ck for an application compared to all benchmark experiments in its library and develops weighting factors for a selected population proportional to the ck values. The methodology can be interpreted as statistically emulating a validation exercise for a particular application where the weighting factors may be viewed as the likelihood that an analyst would include a particular benchmark experiment within the population. The effective sample size of the population is the expected or mean number of benchmark experiments in the population. The uniformly ordered binary decision algorithm identifies clusters of correlated benchmark experiments within the population and then computes adjusted weighting factors based on the magnitude of the correlation coefficients within the cluster to compute a reduced effective sample size accounting for the lower information content because of correlations. Benchmark experiments within the cluster are ordered randomly with equal probability and probabilistic decisions are made as to whether a benchmark. experiment within the cluster should treated as redundant with a previous one; if two redundant benchmark experiments are included, then the conservative worst case bias plus bias uncertainty is used and the pair is counted as a single benchmark experiment in the population. Results are provided for HEU solutions in a research version of the Whisper software using benchmark experiment correlations provided by DICE, the Database for the International Criticality Safety Benchmark Evaluation Project (ICSBEP). These show that there can be a significant increase in the bias plus bias uncertainty because the effective sample size is reduced, and therefore the algorithm, needing to meet sample size requirements, expands the benchmark experiment population by accepting less similar benchmark experiments that would have otherwise not been included.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Combined Nitrogen and Phosphorous Recovery via Electrochemical Technology Integration into Municipal Wastewater Treatment Plants

This project developed electroN-P, an electrochemical technology for recovering nitrogen and phosphorus as a fertilizer product from anaerobic digester centrate, targeting an advancement from TRL 4 to 6. The reactor was scaled 220× from a 250 mL batch cell to a continuous 4-channel system treating 55 L of wastewater, recovering >80% of phosphorus. Energy consumption under constant-voltage operation was lower than the embedded energy of conventional fertilizers, while constant-current operation produced a cost-competitive product with a smaller reactor footprint. TEA and LCA outcomes were highly sensitive to the magnesium source; alternative magnesium salt configurations projected reduced costs and embedded energy compared to sacrificial rods, but at the expense of significantly longer operation times. The recovered fertilizer performed comparably to diammonium phosphate and triple superphosphate in soil and plant trials. Integrating the technology into a whole-plant model reduced aeration energy by nearly 50%, lowering the levelized cost of water treatment from $0.668/m 3 to $0.648/m 3 . Results support targeting commercialization at smaller (≤1 MGD) facilities with high-strength digester streams. Further investigation into the transport and corrosion kinetics governing sacrificial anode wear is recommended to support continuous, longer-term operation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗
DOE OSTI2025

Real-Time Automated pH Control within Batch Processes Relying on Raman pH Measurement

Nuclear fission is an energy source that can provide consistent power with very low associated carbon emissions. However, management of the used nuclear fuel is an important aspect of the application of nuclear power. Recycling of useful components from used fuel is an attractive option, but this involves chemical processing of the fuel. Possible chemical separation technologies that might be used in this regard are sensitive to solution pH. Raman spectroscopy is a promising technique for monitoring the pH of solutions in real time. Classical pH probes are too fragile to be used in the harsh environments encountered in nuclear fuel processing. Raman probes are robust and can withstand these harsh environments to track pH. Coupled with chemometric analysis, the demonstration of the use of Raman spectroscopy to track and predict the pH in carboxylate-buffered systems is made possible. Utilizing this spectroscopy in conjunction with Programmable Logic Controllers mimics industrial control systems used in many modern industrial settings. This showcases a pragmatic approach toward leveraging Raman spectroscopy and chemometric model outputs as inputs for a real-time control system. The model to predict pH created by chemometrics proved to be successful in tracking pH. The optimal pH for TALSPEAK extraction of lanthanides and actinides from aqueous solution is known to proceed in a narrow pH range of around pH = 2.8 ± 0.1. This study uses Raman optical monitoring and automated control to return and maintain solution pH within this range after acid or base perturbations move the solution pH well outside this region. Root-mean-square errors show that pH changes measured using Raman spectroscopy on the batch process solution are reliably measured and used to automatically correct and maintain solution pH. Measurement of solution pH tracks favorably with electrochemical pH probe comparison measurements. As a result, the ability to showcase Raman spectroscopy paired with chemometrics analysis acts as a durable, better alternative data source compared to traditional pH probes to optimize the separation efficiency in the used nuclear fuel processing.

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

Cu-, Co-, and Zn-Based Metal–Organic Framework-Derived Nanoporous Ion Emitters for Picogram Level Analysis of Actinides

Thermal ionization mass spectrometry (TIMS) is often regarded as the preferred technique for trace-level isotopic analysis of actinides owing to its high sensitivity and absence of carry-over effects. However, actinide sample utilization efficiency (SUE) is typically low (<0.05%) without the use of activators or specialized loading approaches which can yield SUEs of >5%. To this effect, we investigate a series of metal–organic framework (MOF)-based nanoporous ion emitters (nano-PIEs) that show enhanced ionization of actinides when used in TIMS loading. We study the impact of physical and chemical properties of MOFs on TIMS SUEs using two families of MOFs that can be synthesized under similar reaction conditions. The structural and chemical properties of these MOFs can be systematically modified one at a time while keeping other features the same. This allows us to strategically investigate their impact on SUEs. The first family of MOFs considered in this study is Zeolitic Imidazole Frameworks (ZIFs) which are made using 2 methyl-imidazole as an organic linker with zinc (ZIF-8) and cobalt (ZIF-67) as metal centers. Additionally, the effect of morphology was also studied using Zn-based ZIF-L with a 2-dimensional structure. The second family of MOFs was synthesized using benzene-1,3,5-tricarboxylate (BTC) with copper (Cu-BTC) and zinc (Zn-BTC) as the metal center. Among the MOFs tested, Cu-BTC showed the highest SUE with an average SUE of 0.27 ± 0.15%, followed closely by Zn-BTC (0.24 ± 0.10%), ZIF-8 (0.17 ± 0.10%), and trailed by the other MOFs. When the MOFs were pyrolyzed in N 2 before loading, an apparent increase in the SUE was observed with ZIF-8 and Cu-BTC showing average SUEs of 0.25 ± 0.08% and 0.34 ± 0.13%, respectively, with the highest measured SUE of 0.53% for pyrolyzed Cu-BTC. In conclusion, this observed increase in SUE by up to an order of magnitude compared to bare filaments demonstrates the efficacy and potential of MOF-derived nano-PIEs for TIMS application.

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

Standoff Detection of Chemical Plumes Using Swept-ECQCL Remote Detection Technology: Results from EMRTC Complex Terrain Dual Tracer Experiment

As part of the Complex Terrain Dual Tracer Experiments performed in support of the LYNM program, a swept-wavelength external cavity quantum cascade laser (swept-ECQCL) remote detection system was deployed to measure chemical plume propagation with high temporal resolution. The swept-ECQCL remote detection technology was initially developed at PNNL, and is now under further development at Opticslah, LLC under SBIR funding from DOE (DOE SBIR Phase I Project DE-SC0019855, Remote Detection Technologies, Program Manager Chris Ramos). The swept-ECQCL technology provides a single-frequency infrared laser source which is rapidly scanned over a large wavelength range to measure the infrared absorption of gases, liquids, or solids. Analysis of the measured infrared absorption spectrum is used to detect, identify, and quantify the chemical species which are present. While most previous laser-based techniques used small wavelength tuning ranges to measure only a few isolated spectral lines of a few species, the swept-ECQCL provides the large tuning range needed to measure multiple species simultaneously, and to detect broad absorption features from large molecules or condensed phase materials. In addition, the swept-ECQCL systems have high enough spectral resolution and sensitivity to measure narrow lines of small molecules and determine isotope ratios. The remote detection system for gas plumes is based on the swept-ECQCL technology, and consists of a swept-ECQCL source, infrared detector, control electronics, and beam transmitting and alignment optics. The remote detection system directs the infrared output of the swept-ECQCL to a target retro reflector typically located at ~10-10,000 m distance. The swept-ECQCL wavelength is varied continuously over a large wavelength range (>1000 nm) at rates up to 1000 Hz. The swept ECQCL light reflected from the remote target is collected and focused onto an infrared photodetector. By measuring the detector signal continuously as the swept-ECQCL wavelength is varied, a time-series of absorption spectra are measured for the gases along the beam line of-sight. Based on analysis and fitting of the measured absorption spectra, the time-dependent chemical concentrations of species are determined. Multiple species and mixtures can be detected simultaneously, for all species with absorption features above the sensor noise floor. The sensor output is thus a continuous, high-speed record of identified chemicals along with their concentrations and isotope ratios of interest. The high-speed operation enables low-noise operation in conditions of atmospheric turbulence and allows tracking of fluctuating concentrations in transient chemical plumes as needed for reliable remote detection applications. All components operate with minimal cooling requirements and no cryogens are needed. In addition, the output infrared light is invisible with an intensity below the maximum permissible exposure (MPE) threshold of 100 mW/cm2 for these infrared wavelengths and thus eye safe. The measurement concept for plume detection using the swept-ECQCL remote detection system is as follows. The ECQCL beam is directed toward a remote retro-reflector, which reflects the beam back to the ECQCL system where it is focused onto an infrared photodetector. The line-of-sight from the ECQCL system to the retro-reflector defines the measurement path. The detected infrared signal is recorded continuously and detects the absorption of molecular species along the measurement path. Before a plume release, the system measures absorption from atmospheric constituents, which in this LWIR band consist primarily of lines from H 2 O and CO 2 , with additional weaker absorption from N 2 O, NH 3 , and O 3 possibly contributing. The system then measures any changes in absorption due to the tracer compounds when they enter the measurement path. The shape of absorption features is used to distinguish and identify the different tracer compounds, and the strength of the absorption is used to determine the path-integrated concentration – also called the column density.

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

Entangled Photon Spectroscopy

Conspectus The enhanced interest in quantum-related phenomena has provided new opportunities for chemists to push the limits of detection and analysis of chemical processes. As some have called this the second quantum revolution, a time has come to apply the rules learned from previous research in quantum phenomena toward new methods and technologies important to chemists. While there has been great interest recently in quantum information science (QIS), the quest to understand how nonclassical states of light interact with matter has been ongoing for more than two decades. Our entry into this field started around this time with the use of materials to produce nonclassical states of light. Here, the process of multiphoton absorption led to photon-number squeezed states of light, where the photon statistics are sub-Poissonian. In addition to the great interest in generating squeezed states of light, there was also interest in the formation of entangled states of light. While much of the effort is still in foundational physics, there are numerous new avenues as to how quantum entanglement can be applied to spectroscopy, imaging, and sensing. These opportunities could have a large impact on the chemical community for a broad spectrum of applications. Here, we discuss the use of entangled (or quantum) light for spectroscopy as well as applications in microscopy and interferometry. The potential benefits of the use of quantum light are discussed in detail. From the first experiments in porphyrin dendrimer systems by Dr. Dong-Ik Lee in our group to the measurements of the entangled two photon absorption cross sections of biological systems such as flavoproteins, the usefulness of entangled light for spectroscopy has been illustrated. These early measurements led the way to more advanced measurements of the unique characteristics of both entangled light and the entangled photon absorption cross-section, which provides new control knobs for manipulating excited states in molecules. The first reports of fluorescence-induced entangled processes were in organic chromophores where the entangled photon cross-section was measured. These results would later have widespread impact in applications such as entangled two-photon microscopy. From our design, construction and implementation of a quantum entangled photon excited microscope, important imaging capabilities were achieved at an unprecedented low excitation intensity of 10 7 photons/s, which is 6 orders of magnitude lower than the excitation level for the classical two-photon image. New reports have also illustrated an advantage of nonclassical light in Raman imaging as well. From a standpoint of more precise measurements, the use of entangled photons in quantum interferometry may offer new opportunities for chemistry research. Experiments that combine molecular spectroscopy and quantum interferometry, by utilizing the correlations of entangled photons in a Hong–Ou–Mandel (HOM) interferometer, have been carried out. The initial experiment showed that the HOM signal is sensitive to the presence of a resonant organic sample placed in one arm of the interferometer. In addition, parameters such as the dephasing time have been obtained with the opportunity for even more advanced phenomenology in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Uranium-Containing and Thorium-Containing Anions Studied by Photoelectron Spectroscopy

An in-depth knowledge of actinide chemistry is fundamental to many aspects of nuclear science and technology, including the synthesis and processing of materials and the remediation of waste disposal sites. Among the actinides, the chemical bonding behaviors of actinium and thorium resemble those of the transition metals; the 5f-electrons of protactinium, uranium, neptunium, and plutonium often play important roles in their bonding; and among the still heavier elements, their bonding tends to mimic the lanthanide elements in terms of electron shielding and their f-electron contributions. Bonding that involves 5f-electrons, however, is especially important, in part because of the significance of uranium and plutonium, but also because these elements are among the few where f-electron participation in bonding is relatively common. This work focused on studying uranium-containing and thorium-containing anions in the gas phase using negative ion photoelectron spectroscopy. Since this technique directly probed valence electrons, it was uniquely positioned to address open questions regarding molecular bonding and electron configurations. A particularly important issue concerned how bonding in actinide-containing molecules was affected by modifications to their actinide atoms’ environment, i.e., due to their interaction with ligands. A closely related question was how actinide atoms’ suborbitals were qualitatively reordered and their energies quantitatively shifted as a result of their ligated environments. These were especially relevant issues in regard to uranium due to it having multiple possible oxidation states (OS) and the potential for 5f electron participation in bonding. The effects of ligands on oxidation states and 5f-orbital energies in uranium bonding was expected to be pronounced. Both ligands and excess electrons were seen as probes of actinide atoms within actinide-containing molecules. Our strategy for advancing knowledge of chemical bonding in the actinide-containing species utilized the synergy between experiments and theory, where in some cases experimental results validated theory and where in others computational results assisted in interpreting experiments. Calculations on actinide systems are terrifically challenging due to large spin-orbit interactions, relativistic effects, and just the sheer number of electrons involved. Even in the simplest species, e.g., U and U2, the most sophisticated, modern calculations carried out by the most experienced theorists often only approximate experimentally-measured values, such as electron affinities. For theory to provide confident predictions that can be used to solve real problems it needed an iterative and ultimately corrective mechanism by which its methods can develop further. Experiments can be used to identify when theory has failed; whereupon the subsequent process of using the experiment-theory interplay can be used to find the cause of the failure. Upon fixing it in one case, different test species can be proposed and studied by the experiment-theory combination to determine whether the problem has been corrected. Thus, experiments not only measure the values of molecular properties, they also provide navigational 3 beacons that keep computations off the reefs in an otherwise dark sea with few reference points. Experimental measurements in the actinide field are not only important, they are in actuality essential to computational progress. While it was not always possible to compare the theoreticallydetermined quantity of interest directly with the same experimentally-measured observable, it was usually possible to compare consequential properties that are both calculable and measurable. In the work completed here electron affinities and electronic state spacings were often sensitive consequential parameters. Reasonable agreement between measured and computational values signaled that a calculation that was very likely to be on-track. We had established collaborative relationships with five computational groups, all of which have expertise in computational actinide chemistry. Their PI’s are L. Cheng, D. Dixon, L. Gagliardi, K. Peterson, and B. Vlaisavljevich. Our close interaction with our theory partners led to us suggesting systems to them and them to us. This reciprocal interaction between our experimental and their computational results was among the most important strengths of this work and was a thread woven throughout. Even though anion photoelectron spectroscopic studies are conducted on anions, much of the information that they provide, pertains to the electronic structure of the neutral counterparts of those anions; among these are electron affinities and electronically excited state spacings. Our experimental tools included several specialized ion sources for forming the anionic species of interest, a mass spectrometer for identifying and mass-selecting them, and an anion photoelectron spectrometer for determining their electron affinities (EA) and characterizing the electronic states of the selected anions’ neutral counterparts. Anion photoelectron spectroscopy is conducted by crossing a mass-selected beam of anions with a fixed-frequency laser beam and energy-analyzing the resultant photodetached electrons. The photodetachment process is governed by the energyconserving relationship: hν = EBE + EKE, where hν is the photon’s energy, EBE is the electron binding (photodetachment transition) energy, and EKE is the electron’s kinetic energy. In our apparatus mass-selection is accomplished via time-of-flight mass spectrometry (TOF-MS), electron energy analysis is achieved with either a magnetic bottle or by velocity mapped imaging. Photodetachment of electrons from anions is implemented via either Nd:YAG or excimer lasers. The photodetachment transition energy, i.e., the EBE, between the ground vibrational and electronic state of an anion and the ground vibrational and electronic state of that anion’s neutral counterpart is the adiabatic electron affinity (EA) of that neutral molecule. Likewise, photodetachment transitions between the ground vibrational and electronic state of an anion and the various electronically-excited states of that anion’s corresponding neutral map the electronic spectrum of that neutral species, i.e., the spectral spacings in the photoelectron spectrum are a mirror image of the neutral’s electronic spectrum. It was, of course, crucial to be able to form the anionic species of interest. There, we had a particularly broad field of anion sources from which to choose. These included several variants of pulsed laser vaporization (LV), laser photoemission, infrared desorption plus photoemission, pulsed arc discharge (PACIS), electrospray ionization (ESI), and Rydberg electron transfer (RET). Each of these anion sources were readily combined with, i.e., connected to, the anion photoelectron spectroscopic portion of our apparatus as described above. Among the sources that utilize lasers, visible light for LV sources as well as IR for desorption sources are provided by Nd:YAG lasers. Ultraviolet photons are provided by both Nd:YAG and excimer lasers, whereas the excitation wavelengths for RET experiments come from two Nd:YAG-pumped dye lasers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗