Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “maximum cut”

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 19 records

Comparing three generations of D-Wave quantum annealers for minor embedded combinatorial optimization problems

Abstract Quantum annealing (QA) is a novel type of analog computation that aims to use quantum mechanical fluctuations to search for optimal solutions of Ising problems. QA in the transverse Ising model, implemented on D-Wave quantum processing units, are available as cloud computing resources. In this study we report concise benchmarks across three generations of D-Wave quantum annealers, consisting of four different devices, for the NP-hard discrete combinatorial optimization problems unweighted maximum clique and unweighted maximum cut on random graphs. The Ising, or equivalently quadratic unconstrained binary optimization, formulation of these problems do not require auxiliary variables for order reduction, and their overall structure and weights are not highly variable, which makes these problems simple test cases to understand the sampling capability of current D-Wave quantum annealers. All-to-all minor embeddings of size 52, with relatively uniform chain lengths, are used for a direct comparison across the Chimera, Pegasus, and Zephyr device topologies. A grid-search over annealing times and the minor embedding chain strengths is performed in order to determine the level of reasonable performance for each device and problem type. Experiment metrics that are reported are approximation ratios for non-broken chain samples, chain break proportions, and time-to-solution for the maximum clique problem instances. How fairly the quantum annealers sample optimal maximum cliques, for instances which contain multiple maximum cliques, is quantified using entropy of the measured ground state distributions. The newest generation of quantum annealing hardware, which has a Zephyr hardware connectivity, performed the best overall with respect to approximation ratios and chain break frequencies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantum approximate multi-objective optimization

The goal of multi-objective optimization is to understand optimal trade-offs between competing objective functions by finding the Pareto front, that is, the set of all Pareto-optimal solutions, where no objective can be improved without degrading another one. Multi-objective optimization can be challenging classically, even if the corresponding single-objective optimization problems are efficiently solvable. Thus, multi-objective optimization represents a compelling problem class to analyze with quantum computers. Here we use a low-depth quantum approximate optimization algorithm to approximate the optimal Pareto front of certain multi-objective weighted maximum-cut problems. We demonstrate its performance on an IBM Quantum computer, as well as with matrix product state numerical simulation, and show its potential to outperform classical approaches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analyzing the Quantum Approximate Optimization Algorithm: Ansätze, Symmetries, and Lie Algebras

The quantum approximate optimization algorithm (QAOA) has been proposed as a method to obtain approximate solutions for combinatorial optimization tasks. In this work, we study the underlying algebraic properties of three QAOA ansätze for the maximum-cut problem on connected graphs, while focusing on the generated Lie algebras as well as their invariant subspaces. Specifically, we analyze the standard QAOA ansatz as well as the orbit and multiangle ansätze. We are able to fully characterize the Lie algebras of the multiangle ansatz across arbitrary connected graphs, finding that they only fall into one of just six families. Aside from the cycle and path graphs, the Lie dimensions for every graph are exponentially large in the system size, meaning that multiangle ansätze are extremely prone to exhibiting barren plateaus. Then, a similar quasi-graph-independent Lie-algebraic characterization beyond the multiangle ansatz is impeded as the circuit exhibits additional “hidden” symmetries besides those naturally arising from a certain parity-superselection operator and all automorphisms of the considered graph. Disregarding the “hidden” symmetries, we can upper bound the dimensions of the orbit and the standard Lie algebras, and the dimensions of the associated invariant subspaces are determined via explicit character formulas. To finish, we conjecture that (for most graphs) the standard Lie algebras have only components that are either exponential or that grow, at most, polynomially with the system size. This would imply that the QAOA is either prone to barren plateaus or classically simulable. More generally, our work provides a symmetry framework and tools to analyze any desired variational quantum algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The CMS Phase II Outer Tracker Upgrade at Fermilab

● Extracted threshold in electrons from S-curve measurement ● Estimated cut for noisy channels based on maximum module occupancy 10−6 ● ~40 channels offline masked ● Measured efficiencies and noise occupancies as a function of the threshold ● Source of inefficiency in the pixels due to punch through structure ● Inefficiencies in the strips after irradiation being investigated, the threshold used was not correct

Salimbeni, Nicolò↗

DESI peculiar velocity survey – Fundamental Plane

The Dark Energy Spectroscopic Instrument (DESI) peculiar velocity survey aims to measure the peculiar velocities of early- and late-type galaxies within the DESI footprint using both the Fundamental Plane and optical Tully–Fisher relations. Direct measurements of peculiar velocities can significantly improve constraints on the growth rate of structure, reducing uncertainty by a factor of approximately 2.5 at redshift 0.1 compared to the DESI Bright Galaxy Survey’s redshift space distortion measurements alone. We assess the quality of stellar velocity dispersion measurements from DESI spectroscopic data. These measurements, along with photometric data from the Legacy Survey, establish the Fundamental Plane relation and determine distances and peculiar velocities of early-type galaxies. During survey validation, we obtain spectra for 6698 unique early-type galaxies, up to a photometric redshift of 0.15. 64 per cent of observed galaxies (4267) have relative velocity dispersion errors below 10 per cent. This percentage increases to 75 per cent if we restrict our sample to galaxies with spectroscopic redshifts below 0.1. We use the measured central velocity dispersion, along with photometry from the DESI Legacy Imaging Surveys, to fit the Fundamental Plane parameters using a 3D Gaussian maximum likelihood algorithm that accounts for measurement uncertainties and selection cuts. In addition, we conduct zero-point calibration using the absolute distance measurements to the Coma cluster, leading to a value of the Hubble constant, H 0 = 76.05 ± 0.35 (statistical) ±0.49 (systematic Fundamental Plane) ±4.86 (statistical due to calibration) km s –1 Mpc –1 ⁠. This H 0 value is within 2σ of Planck cosmic microwave background results and within 1σ of other low-redshift distance indicator-based measurements.

cosmological parameters↗

In Search of Optimum Fresh-Cut Raw Material: Using Computer Vision Systems as a Sensory Screening Tool for Browning-Resistant Romaine Lettuce Accessions

The popularity of ready-to-eat (RTE) salads has prompted novel technology to prolong the shelf life of their ingredients. Fresh-cut romaine lettuce is widely used in RTE salads; however, its tendency to quickly discolor continues to be a challenge for the industry. Selecting the ideal lettuce accessions for use in RTE salads is essential to ensure maximum shelf life, and it is critical to have a practical way to assess and compare the quality of multiple lettuce accessions that are being considered for use in fresh-cut applications. Thus, in this work we aimed to determine whether a computer vision system (CVS) composed of image acquisition, processing, and analysis could be effective to detect visual quality differences among 16 accessions of fresh-cut romaine lettuce during postharvest storage. The CVS involved a post-capturing color correction, effective image segmentation, and calculation of a browning index, which was tested as a predictor of quality and shelf life of fresh-cut romaine lettuce. The results demonstrated that machine vision software can be implemented to replace or supplement the scoring of a trained panel and instrumental quality measurements. Overall visual quality, a key sensory parameter that determines food preferences and consumer behavior, was highly correlated with the browning index, with a Pearson correlation coefficient of −0.85. Other important sensory decision parameters were also strongly or moderately correlated with the browning index, with Pearson correlation coefficients of −0.84 for freshness, 0.79 for off odor, and 0.57 for browning. The ranking of the accessions according to quality acceptability from the sensory evaluation produced a similar pattern to those obtained with the CVS. This study revealed that multiple lettuce accessions can be effectively benchmarked for their performance as fresh-cut sources via a CVS-based method. Future opportunities and challenges in using machine vision image processing to predict consumer preferences for RTE salad greens is also discussed.

Agriculture↗

Disc dichotomy signature in the vertical distribution of [Mg/Fe] and the delayed gas infall scenario

Context. Analysis of the Apache Point Observatory Galactic Evolution Experiment project (APOGEE) data suggests the existence of a clear distinction between two sequences of disc stars in the [α/Fe] versus [Fe/H] abundance ratio space, known as the high- and low-α sequence, respectively. This dichotomy also emerges from an analysis of the vertical distribution of the [α/Fe] abundance ratio. Aims. We aim to test whether the revised two-infall chemical evolution models designed to reproduce the low- and high-α sequences in the [α/Fe] versus [Fe/H] ratios in the solar neighbourhood are also capable of predicting the disc bimodality observed in the vertical distribution of [Mg/Fe] in APOGEE DR16 data. Methods. Along with the chemical composition of the simple stellar populations born at different Galactic times predicted by our reference chemical evolution models in the solar vicinity, we provide their maximum vertical height above the Galactic plane |zmax| computed assuming the relation between the vertical action and stellar age in APOGEE thin-disc stars. Result. The vertical distribution of the [Mg/Fe] abundance ratio predicted by the reference chemical evolution models is in agreement with that observed when combining the APOGEE DR16 data (chemical abundances) with the astroNN catalogue (stellar ages, orbital parameters) for stars younger than 8 Gyr (only low-α sequence stars). Including the high-α disc component, the dichotomy in the vertical [Mg/Fe] abundance distribution is reproduced considering the observational cut in the Galactic height of |z|< 2 kpc. However, our model predicts an overly flat (almost constant) growth of the maximum vertical height |zmax| quantity as a function of [Mg/Fe] for high-α objects in contrast with the median values from APOGEE data. Possible explanations for such a tension are that: (i) the APOGEE sample with |z|< 2 kpc is more likely than ours to be contaminated by halo stars, causing the median values to be kinematically hotter, and (ii) external perturbations – such as minor mergers – that the Milky Way experienced in the past could have heated up the disc, and the heating of the orbits cannot be modeled by only scattering processes. Assuming a disc dissection based on chemistry for APOGEE-DR16 stars (|z|< 2 kpc), the observed |zmax| distributions for high-α and low-α sequences are in good agreement with our model predictions if we consider the errors in the vertical action estimates in the calculation. Moreover, a better agreement between predicted and observed stellar distributions at different Galactic vertical heights is achieved if asteroseismic ages are included as a constraint in the best-fit model calculations. Conclusions. The signature of a delayed gas infall episode, which gives rise to a hiatus in the star formation history of the Galaxy, are imprinted both in the [Mg/Fe] versus [Fe/H] relation and in vertical distribution of [Mg/Fe] abundances in the solar vicinity.

79 ASTRONOMY AND ASTROPHYSICS↗

The Effects of Plasma Pre-Cracking on Drilling of Hard Rocks: A Single Insert Cutting Experiment

Abstract This paper studies the effects of plasma-induced cracks on rock cutting to support the concept of a plasma-integrated drag bit for accelerated geothermal drilling through hard rocks. For this, a single polycrystalline diamond compact (PDC) drilling technique is used in cutting granite to compare thrust and cutting forces between plasma-treated and untreated rock samples. The cracks are produced using underwater plasma at 80 J per pulse. This energy level does not produce visible damage to the rock. The cutting tests are conducted at a cutting speed of 12.7 m/min and four feed rates of 0.127, 0.201, 0.267, and 0.414 mm/s to represent actual drilling scenarios. The results show a general trend of reduction in both thrust and cutting forces for these feed rates, but the magnitude of reduction highly depends on the feed rate. The maximum force reduction of around 50% is found at the 0.267 mm/s feed rate with statistical significance. Cases with a higher force reduction are also found to have rougher surface topography, which indicates more excessive fracturing and, thus, a cracks-accelerated material removal process. This study provides initial evidence of using underwater plasma to increase the downhole drilling rate of hard rocks.

Engineering↗

Five Year Comparison of Mixing Height Determinations at the Savannah River Site

Air quality dispersion modeling is performed for the Savannah River Site (SRS) to demonstrate compliance with applicable regulations. The AMS/EPA Regulatory Model (AERMOD) modeling system is an EPA recommended model for air quality applications with a data preprocessor (AERMET) to incorporate meteorological data collected on site. AERMET parameterizes or calculates meteorological variables that are not directly measured onsite. One of the parameters estimated by AERMET is the atmospheric mixing height. While the mixing height is not currently a measurement input into AERMET, SRS has the capability to measure the local mixing height. The Savannah River National Laboratory (SRNL) operates a Vaisala CL31 Lidar Ceilometer which estimates mixing height from aerosol backscatter. This study compares the parameterized mixing height from AERMET to the ceilometer estimated mixing height for the current regulatory period at SRS incorporating data from 2015-2019. Results from this study showed the average daily minimum values (morning) from AERMET were an order of magnitude lower than the commonly used Holzworth (1972) method and the ceilometer estimated mixing heights. Additionally, on average, the ceilometer exhibited a daily maximum mixing height value that occurred 1-3 hours later than the AERMET estimated maximum. This difference is likely due to the nighttime atmospheric mixing height assumptions and calculations used by AERMET. The AERMET algorithm cuts off mixing height growth at sunset while the ceilometer data show ongoing evening convection typical of the southeastern United States. These results suggest that the AERMET parametrization scheme assumptions may not be representative of a forested landscape and evening convection which could account for more mixing overnight. The results obtained in this study are significant for air dispersion modeling applications for regulatory purposes and worker safety. Mixing height can impact model estimated pollutant concentrations. A greater mixing height will provide more volume for pollutant dispersion. This report documents efforts to quantify the dependence of mixing height inputs toward a conservative estimated pollutant concentration.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Analysis of Vector Particle-In-Cell (VPIC) memory usage optimizations on cutting-edge computer architectures

Vector Particle-In-Cell (VPIC) is one of the fastest plasma simulation codes in the world, with particle numbers ranging from one trillion on the first petascale system, Roadrunner, to ten trillion particles on the more recent Blue Waters supercomputer. As supercomputers continue to grow rapidly in size, so too does the gap between computing capability and memory capability. Current memory systems limit VPIC simulations greatly as the maximum number of particles that can be simulated directly depends on the available memory. In this study, we present a suite of VPIC memory optimizations (i.e., particle weight, half-precision, and fixed-point optimizations) that enable a significant increase in the number of particles in VPIC simulations. Here, we assess the optimizations’ impact on memory and runtime performance for a suite of cutting-edge computer architectures such has the NVIDIA V100 GPU, the IBM Power9, and the Fujitsu A64FX architectures. Our optimizations enable a 31.25% reduction in memory usage and up to 40% increase in the number of particles. This paper extends our work on developing particle storage format optimizations Tan et al.

97 MATHEMATICS AND COMPUTING↗

Development of a Superconducting Adaptive Gap Undulator Prototype for NSLS-II

The photon flux and brightness of synchrotron radiation, crucial parameters for any light source, vary significantly depending on the type of source employed. Among the 23 Insertion Device (ID) sources at the National Synchrotron Light Source II (NSLS-II) at Brookhaven National Lab (BNL), the 12-year-old 3m-long In-Vacuum Undulator (IVU20) stands out for its superior performance, although it no longer represents the cutting edge of technology. Recently, there has been a shift in focus towards developing next-generation sources, particularly Superconducting Undulators (SCUs), characterized by smaller gaps, shorter periods, and maximum lengths. However, despite ongoing research and development efforts, SCUs have yet to surpass their predecessors, the Cryogenic Permanent Magnet Undulators (CPMUs), in terms of performance. This is largely attributed to the limitations posed by traditional superconducting wire, as well as challenges in the design of the magnetic structure and vacuum chamber. In this paper, we aim to overcome such limitations through the development of a unique prototype Superconducting Adaptive Gap Undulator (SC-AGU) magnet core and vacuum chamber design. This paper will outline a novel technical approach aimed at constructing a compact prototype magnet array utilizing state-of-the-art superconducting wire technology. This approach provides a more efficient magnetic structure, allowing for enhanced magnetic field strength and stability.

36 MATERIALS SCIENCE↗

Resource Characterization to Estimate Potential for Electricity Co-Production at Blackburn Oil Field, Nevada

The U.S. Department of Energy estimates that an annual average of 25 billion barrels of hot water are produced from oil and gas wells within the United States. The thermal energy available in the co-produced water stream is usually discarded, as the produced waters are considered an inconvenience by the operators and are disposed of using injection wells. However, utilizing organic Rankine cycle (ORC) generators, a vast amount of thermal energy can be captured and converted into electricity (albeit at relatively low efficiency due to the low temperatures). The National Renewable Energy Laboratory (NREL), in collaboration with Transitional Energy and Grant Canyon Oil & Gas, evaluated the feasibility of geothermal co-production of electricity by utilizing existing oil wells in Blackburn oil field in Nevada. The once prolific Blackburn oil field is located in Pine Valley, approximately 45 miles east-southeast of Elko, Nevada. Currently, the wells targeting the highly fractured Devonian Nevada dolomite reservoir are operating at a water cut ratio of more than 99%, with individual fluid (oil and water) production rates reaching 7.4 L/s (4,021 BBL/day). Analysis of publicly available data showed that the combination of the suitable wells' maximum historical production rates reached 22.90 L/s. The production from these wells occurs naturally and the wells are choked (and even shut down) by the operator to mitigate excessive water production, indicating a strong reservoir recharge and future opportunity to increase the water production for geothermal electricity generation. The main goal of this study was to evaluate the productivity of the existing wells, the performance of the reservoir, the surface network, and the operational constraints in order to achieve 1 MWe of electricity production from the field's water production. Utilizing the GEOPHIRES tool, we have determined that a twofold to threefold increase in the total fluid production, compared to the historical production under artificial restraint (choke), is required to reach a 1-MWe net target output for a low-temperature ORC system with air-cooled condensers. Lower flow rates would be required when utilizing water-based condensers instead of air-cooled condensers. However, that would require a constant supply of cold water, which may be challenging given the arid environment of the project site.

Blackburn↗

Develop Accurate Techniques for Passive SiC Temperature Monitoring of Miniature Samples for Cross-Cutting Applications

Passive thermometry is critically important because most fuels and materials irradiation experiments are not instrumented, and it is necessary to understand the irradiation temperature to properly interpret any post-irradiation examination data, including evolving properties and/or microstructures. The standard passive thermometry approach uses continuous dilatometry to evaluate changes in the instantaneous coefficient of thermal expansion during post-irradiation thermal annealing. This approach has limitations in terms of sample size (minimum length requirements) and the maximum irradiation temperature that can be accurately determined, which is limited by the reduced swelling (and therefore recovery) following higher temperature irradiation and limitations on the furnaces used with push-rod dilatometers. This work evaluates two new proposed techniques for post-irradiation evaluation of passive SiC temperature monitors: differential scanning calorimetry (DSC) and Raman spectroscopy. DSC is an extremely sensitive technique that can be used for any specimen geometry and is capable of higher temperature operation. Raman spectroscopy is similar in that it is a surface technique capable of examining extremely small samples (submillimeter), can be used with a heated stage up to 1,500°C (planned for future work), and is capable of mapping local irradiation temperatures throughout a sample. Existing SiC samples that were previously irradiated over a wide range of temperatures were cut into multiple pieces to allow for annealing studies using multiple different techniques: dilatometry, DSC, and Raman spectroscopy. This approach mitigates the concern that samples analyzed using one technique may have a slightly different irradiation history than those analyzed using a different technique. Recovery was clearly observed during annealing using both DSC and dilatometry. In some cases, a direct comparison could not be made due to some of the DSC runs accidentally including material from multiple specimens and issues with using an alternative DSC sample holder for the highest temperature annealing studies. Nevertheless, one trend was clear: the DSC runs resulted in higher irradiation temperatures compared to those of the dilatometry runs. Part of this could be attributed to the higher temperature ramp rates used during the DSC runs, which are often preferred to reduce noise in the measurements. By comparison, dilatometry has previously been shown to produce better data at lower ramp rates. Future work should further investigate the ideal ramp rate for both techniques to produce consistent results. Additional work should evaluate the best holder material to use for DSC runs exceeding 1,000°C to provide reliable data while preventing interactions between SiC and the holder.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Develop Accurate Techniques for Passive SiC Temperature Monitoring of Miniature Samples for Cross-Cutting Applications

Passive thermometry is critically important because most fuels and materials irradiation experiments are not instrumented, and it is necessary to understand the irradiation temperature to properly interpret any post-irradiation examination data, including evolving properties and/or microstructures. The standard passive thermometry approach uses continuous dilatometry to evaluate changes in the instantaneous coefficient of thermal expansion during post-irradiation thermal annealing. This approach has limitations in terms of sample size (minimum length requirements) and the maximum irradiation temperature that can be accurately determined, which is limited by the reduced swelling (and therefore recovery) following higher temperature irradiation and limitations on the furnaces used with push-rod dilatometers. This work evaluates two new proposed techniques for post-irradiation evaluation of passive SiC temperature monitors: differential scanning calorimetry (DSC) and Raman spectroscopy. DSC is an extremely sensitive technique that can be used for any specimen geometry and is capable of higher temperature operation. Raman spectroscopy is similar in that it is a surface technique capable of examining extremely small samples (submillimeter), can be used with a heated stage up to 1,500°C (planned for future work), and is capable of mapping local irradiation temperatures throughout a sample. Existing SiC samples that were previously irradiated over a wide range of temperatures were cut into multiple pieces to allow for annealing studies using multiple different techniques: dilatometry, DSC, and Raman spectroscopy. This approach mitigates the concern that samples analyzed using one technique may have a slightly different irradiation history than those analyzed using a different technique. Recovery was clearly observed during annealing using both DSC and dilatometry. In some cases, a direct comparison could not be made due to some of the DSC runs accidentally including material from multiple specimens and issues with using an alternative DSC sample holder for the highest temperature annealing studies. Nevertheless, one trend was clear: the DSC runs resulted in higher irradiation temperatures compared to those of the dilatometry runs. Part of this could be attributed to the higher temperature ramp rates used during the DSC runs, which are often preferred to reduce noise in the measurements. By comparison, dilatometry has previously been shown to produce better data at lower ramp rates. Future work should further investigate the ideal ramp rate for both techniques to produce consistent results. Additional work should evaluate the best holder material to use for DSC runs exceeding 1,000°C to provide reliable data while preventing interactions between SiC and the holder.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Energy Dependence Analysis of Neutrino-Nucleus Interactions in Scintillator for MINERvA

This thesis is on the intrinsic energy dependence in neutrino-nucleus interaction models, especially quasielastic scattering and 2p2h process, and the implications for understanding data from MINERvA, NOvA, and DUNE. Neutrino-nucleus interactions are modeled using fve structure functions, W1, W2, W3, W4 and W5, and each structure function contributes diferently to the energy dependence once it enters the cross-section calculations. The energy independent term in the calculations contains only W2 and the next lower order term (of the order 1/E) contains W3 also. In this thesis, we show that W2 is good enough to describe MINERvA data and energy dependence coming from the cross-section directly is small. However, W3 will be signifcant for lower energy experiments such as NOvA and DUNE (especially at and below the oscillation maximum). We have also developed a method to extract structure functions directly from MINERvA data though we find that systematic uncertainties present in the data are too big to successfully determine 1/E Term. In addition to this, we have applied muon angle cuts on MINERvA data and MC, which have their own energy dependent effects. An apparent energy dependent discrepancy comes from known shortcomings especially in the pion production model. Finally, because the energy dependence analysis has similarities to the low-$\nu$ method, we have evaluated discrepancies between the data and MC in light of diferent fux constraint techniques. The results suggest including MINERvA s low-$\nu$ method constraint is an improvement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Improved light yield and growth of large-volume ultrafast single crystal scintillators Cs 2 ZnCl 4 and Cs 3 ZnCl 5

Due to their reported fast decay times, Cs 2 ZnCl 4 and Cs 3 ZnCl 5 are promising candidates for detection of gamma rays and X-rays in high count rate and fast timing applications. In this work, we show that single crystals with better optical quality than previously demonstrated – and larger in size – can be grown via the vertical Bridgman method. Highly transparent Ø7 mm crystals of undoped Cs 2 ZnCl 4 and Cs 3 ZnCl 5 are grown and measured to have light yields surpassing those previously reported, achieving 1980 ± 100 ph/MeV and 1460 ± 70 ph/MeV at 662 keV – a 55% and 232% improvement, respectively. We observe single-component scintillation decay times for both Cs 2 ZnCl 4 (1.66 ns) and Cs 3 ZnCl 5 (0.82 ns) and radioluminescence emission with maximum intensity at ~290 nm. Scalability of these materials is also evaluated based on growth of Ø22 mm crystals. Minimal cracking is observed, and the fast decay times are maintained at this size. Coincidence time resolution of 3 × 3 × 5 mm 3 and 7 × 7 × 10 mm 3 pixels cut from Ø22 mm Cs 2 ZnCl 4 are measured to be 148 ± 1 ps FWHM and 175 ± 1 ps FWHM, respectively. Here, the improved performance and ability to be fabricated in large sizes now place Cs 2 ZnCl 4 and Cs 3 ZnCl 5 on the map as potential contenders for radiation detection applications where BaF 2 – the most commonly used ultrafast inorganic scintillator – is typically considered.

36 MATERIALS SCIENCE↗

Air quality impacts from the development of unconventional oil and gas well pads: Air toxics and other volatile organic compounds

Unconventional oil and natural gas development (UOGD) has expanded rapidly across the United States in recent decades and raised concerns about associated air quality impacts. While significant effort has been made to quantify methane emissions, relatively few observations have been made of Volatile Organic Compounds (VOCs), especially during drilling and completion of new wells. Extensive air monitoring during development of several large, multi-well pads in Broomfield, Colorado, in the Denver-Julesburg Basin, provides a novel opportunity to examine changes in local air toxics and other VOC concentrations during well drilling and completions and production. These operations offer an especially useful case to study as several management practices were implemented to reduce emissions (e.g., electrified, grid-powered drill rigs and closed loop fluid handling systems to reduce truck traffic and limit fluid handling on the pad). With simultaneous measurements of methane and 50 VOCs from October 2018 to December 2022 at as many as 19 sites near well pads, in adjacent neighborhoods, and at a more distant reference location, we identify impacts from each phase of well development and production. Use of weekly, time-integrated canisters, a Proton Transfer Reaction Mass Spectrometer (PTR-MS), continuous photoionization detectors (PID) to trigger canister collection upon detection of VOC-rich plumes, and an instrumented vehicle, provided a powerful suite of measurements to characterize both transient plumes and longer-term changes in air quality. Prior to the start of well development, VOC gradients were small across Broomfield. Once drilling commenced, concentrations of oil and gas (O&G) related VOCs, including alkanes and aromatics, increased around active well pads. Concentration increases were clearly apparent during certain operations, including drilling, coil tubing/millout operations, and production tubing installation. Emissions of C 8 –C 10 n-alkanes during drilling operations highlighted the importance of VOC emissions from synthetic drilling mud chosen to reduce odor impacts. More than 90 samples were collected of transient plumes. Using composition measurements, meteorological data, and information about well pad activities, these plumes were connected with specific UOGD operations including drilling, flowback, and production equipment maintenance. The chemical signatures of these plumes differed by operation type (e.g., C 8 –C 10 n-alkanes constituted a larger fraction of measured VOCs in drilling-related plumes). Concentrations of individual, oil and gas-related VOCs in these plumes were often several orders of magnitude higher than in background air, with maximum ethane and benzene concentrations of 79,600 and 819 ppbv, respectively. Because these plumes typically impact a monitoring site for just several minutes, they are easily missed by slower-responding instruments. Study measurements highlight future emission mitigation opportunities during UOGD operations, including better control of emissions from shakers that separate drill cuttings from drilling mud, production separator maintenance operations, and periodic emptying of sand cans during flowback operations.

54 ENVIRONMENTAL SCIENCES↗