Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase velocity”

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 919 records · Page 51

Improving AHU Performance by Minimizing Approach Temperature, Reducing Air Maldistribution, and Efficiently Handling Sensible and Latent Loads (Phase 1 Interim Final Technical Report)

Residential air handling units (AHUs) have stayed the same in form and efficiency for the past 30+ years, with incremental improvements made to address safety, functionality, and energy-efficiency. The purpose of this research in Phase I, Topic 9a: Next Generation Residential Air Handlers, was to improve AHU performance by minimizing heat exchanger (HX) approach temperature, reducing air maldistribution, and developing alternative system configurations which more efficiently handle sensible and latent loads. In this research Optimized Thermal Systems (OTS) developed, modeled, and evaluated multiple alternative system concepts. A dual vapor compression system separate sensible and latent cooling (SSLC) concept was studied to inform work on alternative concepts and to show best-case performance benefit. System concepts included ejector enhanced vapor compression cycles, desiccant assisted dehumidification, dual evaporator SSLC, and alternative AHU HX configurations. A dual vapor compression system showed COP improvement of 20%, however, required additional components, increased unit size, and increased cost. Two types of ejector enhanced vapor compression cycles with dual evaporators improved system COP by 9 to 11%, and reduced AHU losses by as much as 18%, with design changes limited to the AHU, no unit physical size increase, and a moderate increase to system first cost. Desiccant assisted air-conditioning required increased air flow rate resulting in higher fan power and the desiccant wheel increased sensible heat load leading to increased compressor power and reduced system COP. Dual evaporator cycles were found to degrade performance due to increased expansion losses. Optimized single slab HX designs used in place of the traditional A-coil HX led to 44–49% reduction in aluminum, 47–60% less refrigerant charge, and improved HX velocity distribution.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Freestream turbulence effects on unsteady wind turbine loads and wakes: An IDDES study

We investigate numerically the effects of freestream turbulence on the unsteady aerodynamics and wakes of the National Renewable Energy Laboratory Phase VI wind turbine rotor for increasing wind speed. Turbulence is modeled using the Improved Delayed Detached-Eddy Simulation (IDDES) method. As a first step, a detailed mesh resolution study is conducted with the decaying freestream turbulence model at turbulence intensity of 0.5%. Our blade-resolved IDDES simulations show that grid-independent average torque and thrust results can be achieved with relatively coarse meshes, whereas dramatically higher mesh resolution is required for grid-independent results for power spectral densities of thrust force, especially in the deep-stall regime. Comparing the loads with the Shear-Stress Transport model demonstrates the superiority of IDDES in predicting massively separated flows. The aerodynamic performance and wake predictions with the decaying freestream turbulence model are compared with the synthetic freestream turbulence model. Both models predict nearly the same loads, spectral energy content, and wake characteristics. The properties of both the near- and far-wake regions are then examined. Furthermore, we show that separated boundary layers accelerate turbulent mixing and entrainment of the external flow, which results in faster wake recovery. The effect of increasing turbulence intensity to 6% is investigated using the synthetic freestream turbulence model. In contrast with the fully attached boundary layer, higher freestream turbulence in deep stall does not significantly affect the loads and vortex-shedding characteristics. However, the turbulent mixing in the wake is enhanced, which further hastens the recovery of the self-similar velocity profile. In general, increasing the wind speed at high turbulence intensity shifts the recovery farther upstream and increases the wake width.

17 WIND ENERGY↗

Development of simplified model for injection rate prediction of diesel injectors during transient and steady operation

Determining the injection rate of the diesel injector during the transient short injection operations is a critical matter to enhance the combustion process by improving air/fuel mixture formation. In this regard, numerous injection models have been suggested in the literature to predict the injection rate accurately. However, such algorithms have not precisely predicted the rate of injection (ROI) in transient short injection events. In addition, most of those models require many input parameters to characterize hydraulic and mechanical subsystems in the injector, which can make them inconvenient to tune for different injector specifications, in other words not suitable for universal usage. Therefore, our study aims to develop a simple injection model that can predict injection rate precisely in the short and long injection durations considering relatively fewer input parameters. The model traces the pressure variation inside the sac based on the bulk modulus theory to predict nozzle exit velocity, thus the injection rate of the injector. In the algorithm, the bulk modulus of fuel (E) and discharge coefficients of the nozzle hole and needle seat (C d,H and C d,N ) have been considered to vary rather than set constant during the injection process since such flow parameters would have a critical impact on the flow development at the transient stage, thus affecting the injection rate of diesel injectors. Here with this modelling strategy, the accuracy of the new model on ROI prediction was confirmed with acceptable model error for both short and long injection durations by validating against experimental results. To understand the effect of the new model scheme on ROI prediction, ROI results with the new model have been discussed by comparing it with the ROI results predicted with constant flow parameters.

42 ENGINEERING↗

Powdermet CRADA for VELOCITI Voucher (CRADA Final Report)

NLR shall develop a model of an electric thermal energy storage (ETES) system in the System Advisor Model (SAM) framework that uses a salt and particle slurry developed by PowderMet as the thermal energy storage medium. The ETES shall consist of a recuperated sCO 2 Brayton cycle using 2-tank indirect thermal energy storage and an electric heater for heat input. The salt slurry uses chloride salt and operates at a maximum temperature of 500-720°C. The ETES uses dry air cooling. The salt slurry will be modeled as a sensible heat storage medium and as a latent heat storage medium. NLR will use the model to do parametric analysis of the factors affecting the performance (ex., round trip efficiency) of the ETES system. NLR will work with PowderMet to estimate the costs of system components and total system costs, with the levelized cost of electricity (LCOE) as the primary cost metric.

14 SOLAR ENERGY↗

Radiation and heat transport in divergent shock–bubble interactions

Shock–bubble interactions (SBIs) are important across a wide range of physical systems. In inertial confinement fusion, interactions between laser-driven shocks and micro-voids in both ablators and foam targets generate instabilities that are a major obstacle in achieving ignition. Experiments imaging the collapse of such voids at high energy densities (HED) are constrained by spatial and temporal resolution, making simulations a vital tool in understanding these systems. In this study, we benchmark several radiation and thermal transport models in the xRAGE hydrodynamic code against experimental images of a collapsing mesoscale void during the passage of a 300 GPa shock. We also quantitatively examine the role of transport physics in the evolution of the SBI. This allows us to understand the dynamics of the interaction at timescales shorter than experimental imaging framerates. We find that all radiation models examined reproduce empirical shock velocities within experimental error. Radiation transport is found to reduce shock pressures by providing an additional energy pathway in the ablation region, but this effect is small (∼1% of total shock pressure). Employing a flux-limited Spitzer model for heat conduction, we find that flux limiters between 0.03 and 0.10 produce agreement with experimental velocities, suggesting that the system is well-within the Spitzer regime. Higher heat conduction is found to lower temperatures in the ablated plasma and to prevent secondary shocks at the ablation front, resulting in weaker primary shocks. Finally, we confirm that the SBI-driven instabilities observed in the HED regime are baroclinically driven, as in the low energy case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High-Temperature Ceramic-Carbonate Dual-Phase Membrane Reactor for Pre-combustion Carbon Dioxide Capture (Final Scientific/Technical Report)

Arizona State University, in collaboration with University of South Carolina, worked on a project aimed at development of a new high temperature, high pressure CO 2 perm-selective membrane reactor for water-gas-shift reaction (WGS) with simulated gasifier syngas to produce a high concentration H 2 stream with CO 2 capture. The membrane reactor is made of a CO 2 semi-permeable ceramic-carbonate dual-phase (CCDP) membrane with high CO 2 perm-selectivity/permeance and thermal/mechanical stability for application in WGS reaction. The objectives of this project were to (1) synthesize the chemically/thermally stable tubular CCDP membranes with CO 2 permeance and selectivity (with respect to H 2 , CO or H 2 O) larger than 6.5×10-7 mol/m2·s·Pa and 500, respectively; (2) establish CCDP membrane reactor setup and study high pressure CO 2 permeation and WGS reaction with CO 2 capture using the setup; and (3) identify conditions for WGS in the CCDP membrane reactor that produce CO 2 and H 2 streams with purity of >99% and >90% respectively at CO conversion >95% and overall carbon capture >90%. The work in this project included both membrane development and membrane reactor process study. The membrane development efforts were focused on investigating a H 2 S resistant and highly oxygen-ionic conducting metal oxide material and membrane for CO 2 separation, fabrication of tubular samaria-doped-ceria/molten-carbonate CCDP membrane with high mechanical strength, and experimental and modeling study of high-pressure CO 2 permeation of the CCDP membranes. Mathematical models were developed to describe WGS in the CCDP membrane reactor without a catalyst or packed with a commercial high temperature WGS catalyst. Experiments on WGS in the CCDP membrane reactor with the commercial WGS catalyst, guided by the model analysis, were performed to identify optimum conditions for achieving the CO conversion, carbon capture, and the purity of the H 2 and CO 2 streams mentioned above. At 30 atm feed pressure, 750°C operation temperature, space velocity of 250 h-1, and with steam sweep, a single-stage CCDP membrane reactor with average CO 2 permeation flux of 0.5 cm3(STP)/min.cm2 can achieve CO 2 conversion of 95% and overall carbon capture of 94%, and produce CO 2 and H 2 streams with dry-based purity of >99% and 92% respectively. The project also included process design and techno-economic analysis (TEA) for a CCDP membrane reactor process for WGS reaction with CO 2 capture for a 550 MW coal-fired IGCC power plant, and its comparison with the conventional fixed-bed reactor system for WGS with follow-up CO 2 capture by an amine absorption process. The target performance for the reactor for WGS with CO 2 capture includes CO conversion >95%, hydrogen stream purity >90%, CO 2 stream purity >95%, and total carbon capture >90%. The CCDP membrane developed in this project can achieve the performance target, without subsequent CO 2 capture process at the optimum conditions identified in this project. The outcome of the process design and TEA analysis shows that the membrane reactor for WGS with in-situ CO 2 capture has an operating cost about 40% lower than that for the conventional fixed-bed reactor with a separate amine absorption process for CO 2 capture. However, the capital cost of the membrane reactor process is about twice that of the conventional process because of the higher cost of the CCDP membrane. Modeling analysis shows that a membrane reactor using a CCDP membrane with higher CO 2 permeance (about three times the current value) can deliver the targeted performance for WGS reaction with CO 2 capture at a much higher space velocity and lower membrane surface area to catalyst volume ratio, leading to a smaller catalyst amount and/or membrane area and hence significantly reduced membrane reactor capital costs.

20 FOSSIL-FUELED POWER PLANTS↗

Date Release Report for Large Surface Explosion Coupling Experiment (LSECE) Nevada National Security Site

The DTRA sponsored 2020-2022 Large Surface Explosion Coupling Experiment (LSECE) consists of two large ground surface chemical explosions, data collection, analysis and modeling carried out in 2020-2022. The LSECE chemical explosions were carried out at the site of prior NNSA sponsored buried chemical explosions that were part of the Source Physics Experiment (SPE) Phase II in Dry Alluvium Geology (DAG). This report described the data collected under LSECE that is being made publicly available. The prior buried explosion DAG data are described in a separate report (https://www.osti.gov/biblio/1825534) and are also publicly available. The LSECE explosions, data analysis and numerical modeling work sponsored by DTRA were intended to address three objectives: 1) Generate seismo-acoustic data to test and improve numerical models of explosion energy coupling in dry alluvium geology; 2) To improve seismo-acoustic yield estimation techniques as a function of depth and medium properties; 3) To study and improve acoustic propagation modeling under two different atmospheric conditions. The LSECE data that were collected have expanded the prior DAG buried explosion dataset to include surface chemical explosions recorded on a common set of seismo-acoustic stations. They have allowed a more detailed study of above and below surface seismo-acoustic energy coupling with applications to explosion monitoring and assessment. The two approximately 1-ton TNT equivalent yield LSECE explosions were detonated at different times of the day to explore the effects of the different atmospheric conditions. The first chemical explosion “Artemis” was conducted before dawn when temperature inversions were present. The second chemical explosion “Apollo” was conducted on a sunny afternoon when the temperature gradient was more linear. The LSECE chemical explosion data were also collected across a variety of different sensor types to allow evaluation of their effectiveness in recovering useful information. The LSECE instrumentation included fiber optic or Distributed Acoustic Sensing (DAS), a dense array or Large-N array of seismometers, borehole accelerometers, DAS and a velocity meter, airborne acoustic instruments and a variety of visual and remote sensing data.

58 GEOSCIENCES↗

Parameterization of electron attachment rate constants for impurities in LArTPC detectors

The ability of free electrons to drift long distances at high velocities in pure liquid argon under an applied electric field has been exploited for the past forty years to implement detectors with increasingly larger volumes for high energy physics research. The attachment of free electrons to impurities in the LAr is an important limit on the free instrumented volume of these extremely large detectors, and impurity concentrations as small as 100 ppt can reduce their resolution and efficiency. In this paper, we summarize the electron attachment rate constants as a function of the applied electric field, for common impurities in LArTPCs, obtained from data in the literature. Here, we further provide analytical functions to parameterize the data, which are useful to compare with new measurements, to model and analyze the performance of existing detectors, and to predict the performance of new detectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Distributed Acoustic Sensing (DAS) Seismic Monitoring of CO 2 Injected for Enhanced Oil Recovery in Northern Michigan (Volume II.C)

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. This report describes the monitoring study conducted to assess the effectiveness of DAS (Distributed Acoustic Sensing) -based VSP (Vertical Seismic Profiling) technology for delineating CO 2 injected into the Silurian-age pinnacle reefs in northern Michigan, the host rocks for the MRCSP Phase III demonstration project. The DAS VSP study was conducted in the Chester 16 reef, one of several reefs in Otsego County Michigan that is operated by Core Energy, LLC of Traverse City, Michigan. Time-lapse DAS VSP was implemented at the Chester 16 reef to attempt to detect approximately 85,000 tonnes of CO 2 injected into the A-1 Carbonate and Brown Niagaran Formations. A baseline survey was conducted in February 2017 prior to injecting CO 2 and a repeat survey was conducted in August 2018. During the interim period between the baseline and repeat surveys, CO 2 was injected into the Chester 16 reef via the 6-16 injection well without production (withdrawal) of fluids from the reef. A grid of 181 source positions consisting of 44 vibrator positions, plus 137 dynamite shot locations, was used to give approximately continuous spatial coverage of the injection zone (A-1 Carbonate and upper Brown Niagaran) in the area between the two wells.

01 COAL, LIGNITE, AND PEAT↗

Real-time damage characterization for GFRCs using high-speed synchrotron X-ray phase contrast imaging

Here, we report the application of high-speed synchrotron X-ray phase contrast imaging (PCI) in real-time damage characterization for glass fiber reinforced composites (GFRCs) subjected to dynamic loading. Dynamic single-edge notched bending (DSENB) experiments on pre-notched S-2 GFRCs were performed on a modified Kolsky compression bar. During loading, the synchrotron X-ray beam penetrated through the composite specimen from the side to detect damage evolution inside the material. Entire dynamic events were recorded by a high-speed camera as image sequences.0° and 90° unidirectional and cross-ply composites were investigated. An optical imaging technique was also employed to capture similar dynamic events in comparison with the radiographic imaging. It is demonstrated that high-speed X-ray PCI had sufficient phase contrast to characterize a crack initiation at a 20-μm spatial resolution within 920 ns and track the crack geometry during propagation, thereby providing reliable data to quantify the dynamic damage resistance of GFRCs. Furthermore, being capable of recognizing microscopic damage-related features at a sub-10-μm resolution, high-speed X-ray PCI provided fundamental material failure mechanisms to reveal the essential of macroscale structural failure of composites. It can also track the damage evolution inside and between individual plies of laminated composites. However, current high-speed X-ray PCI technique only supports in-situ observation and the high timing and spatial resolutions are limited within a field of view of ~2.5 mm in square, preventing its application in the three-dimensional and larger-area damage detection for GFRC structures.

42 ENGINEERING↗

Astrometry in two-photon interferometry using an Earth rotation fringe scan

Optical interferometers may not require a phase-stable optical link between the stations if instead sources of quantum-mechanically entangled pairs could be provided to them, enabling long baselines. We developed a new variation of this idea, proposing that photons from two different astronomical sources could be interfered at two decoupled stations. Interference products can then be calculated in post-processing or requiring only a slow, classical connection between stations. In this work, we investigated practical feasibility of this approach. We developed a Bayesian analysis method for the earth rotation fringe scanning technique and showed that in the limit of high signal-to-noise ratio it reproduced the results from a simple Fisher matrix analysis. Here we identify candidate stair pairs in the northern hemisphere, where this technique could be applied. With two telescopes with an effective collecting area of ~ 2 m 2 , we could detect fringing and measure the astrometric separation of the sources at ~100 μas precision in a few hours of observations, in agreement with previous estimates.

47 OTHER INSTRUMENTATION↗

Big Sky Regional Carbon Sequestration Partnership (Phase III Final Scientific/Technical Report)

The Big Sky Carbon Sequestration Partnership (BSCSP) pursued a Phase III demonstration project at Kevin Dome in north central Montana. Kevin Dome covers approximately 700 square miles and is a naturally occurring CO 2 reservoir that is flanked by oil and gas fields. The carbon dioxide (CO 2 ) is in the upper Devonian Duperow (carbonate) formation and does not reach the spill point of the dome; therefore, the dome has potential as a CO 2 sequestration reservoir, a CO 2 supply, or as both if anthropogenic sources and enhanced oil recovery (EOR) operations are associated with the dome. Kevin Dome could potentially act as a buffer to continue accepting anthropogenic CO 2 when EOR flooding operations are interrupted or completed. The project objective was to produce one million tonnes of CO 2 from the gas cap of Kevin Dome, pipe it laterally, inject, and re-store it in the brine leg of the same formation to test the hub / buffer storage concept. This was to be accomplished by drilling up to five production wells, building a short pipeline and compression facilities, and drilling an injection well and several monitoring wells. BSCSP commenced outreach and site characterization activities including acquiring baseline data for near-surface insurance monitoring, acquiring 3-dimensional, 9-component surface seismic over the project area, drilling two test wells (one in the production area and one in the injection area), coring key intervals, and performing comprehensive logging. Well tests of those wells revealed two barriers to the project. The production (Danielson 33-17) well was perforated in multiple zones but failed to produce any significant CO 2 . This was despite being drilled in the near vicinity of a historic well that had produced 3700 MCF per day in a drill stem test. Modeling indicated that this was likely due to a phase change during production causing a temperature drop resulting in hydrate and/or water ice formation that clogged the formation. Tests of the injection zone (Wallewein 22-1) well indicated total dissolved solids (TDS) slightly below the EPA required 10,000 parts per million (ppm) for a Class VI underground injection control permit. While the project was initiated before Class VI rules were promulgated, and this was an experimental project (seemingly qualified for a Class V permit), the Environmental Protection Agency (EPA) indicated that injection would require a Class VI permit. The low salinity result was unexpected as contours plotted based on regional formation water quality data indicated an expected TDS above 20,000 ppm, and wells between the recharge zone and the Wallewein well tested above 10,000 ppm. Faced with the inability to obtain an injection permit, the demonstration project could not proceed. However, the project had generated valuable samples and data on a large natural analog including 32 sq. mi. of 3-D, 9-C seismic, 430 ft. of carbonate core covering seven different depositional environments taken from areas with, and without the presence of CO 2 , 30 ft. of core of two caprocks, a tight carbonate and an anhydrite, a full set of modern logs on both wells, and well tests. DOE decided to re-scope the project around completing studies utilizing this data. This report covers both the initial scope and the re-scope (Task / Section numbers preceded with an R). While the report covers a wide range of project activities, highlights of this work include: Development of a geostatic model using neural nets to match well logs to facies and using multi-waveform seismic to inform reservoir heterogeneity; Unique mechanical testing of permeability – stress relationship in two caprock materials; Development of full waveform inversion to generate a high resolution velocity model; Model development for dual permeability (fracture and matrix) systems to better account for matrix-matrix interactions; Joint seismic wave inversion (including the first quadr-joint inversion) exhibiting better imaging of a challenging reservoir zone in stiff rock; Core flow and core flood results on a reactive carbonate; and Innovative laboratory measurements of seismic response of fractured core as a function of fluid fill.

54 ENVIRONMENTAL SCIENCES↗

UNH TDP - Concurrent Measurements of Inflow, Power Performance, and Loads for a Grid-Synchronized Vertical Axis Cross-Flow Turbine Operating in a Tidal Estuary

This data was collected between October 12 and December 15 of 2021 at the University of New Hampshire (UNH) and Atlantic Marine Energy Center (AMEC) turbine deployment platform (TDP). This data set includes over 29 days of grid connected turbine operation during this 65 day time frame. The priority for this measurement campaign was to collect data while the turbine was electrically connected to the grid by means of a rectifier and inverter. The Fall_2021_UNH_Measurement_Timeline.png highlights when each instrument was functioning and the Fall_2021_UNH_Test_Log.jpg indicates the four main regions for analysis available from this measurement campaign. The TDP is a floating structure moored on the Portsmouth facing side of Memorial Bridge pier #2, which spans the Piscataqua River between Portsmouth, NH and Kittery, ME. The Piscataqua River connects the Great Bay Estuary to the Gulf of Maine and the river currents are dominated by tidal forcing with water velocities exceeding 2.5 m/s during spring ebb tides at this site which were previously characterized by Kaelin Chancey (Assessment Of The Localized Flow And Tidal Energy Conversion System At An Estuarine Bridge - UNH MS Thesis 2019). The turbine under test was a modified New Energy Corporation (Calgary, CA) model EVG-025 4-blade H-Darrius type vertical axis cross flow turbine that rotates in the clockwise direction with a rotor diameter of 3.2m and blade length of 1.7m. The hydro-foil profile was a NACA 0021 with a 10 inch chord length and a blade preset pitch angle of +4deg with a positive angle corresponding with the toe in direction. The standard EVG-025 has a rotor diameter of 3.4m and its rated power output is 25kW at 3 m/s. The rotor diameter was reduced to accommodate the size of the existing TDP moon-pool. This project was pursued to quantify device performance for cross flow turbines operating in a marine environment. Accurate physical models, to characterize cross flow turbine performance, require real operational data sets due to the complexity of blade fluid interactions. This data can help support model development which will help predict turbine performance when analyzing perspective project locations in the future. Instrumentation was deployed to measure; water speed/direction, electrical power output, turbine shaft speed, turbine thrust force, and platform motion. Concurrent measurements of these parameters allow for correlations (cause and affect) to be inferred, allowing for characterization of device performance over a range of operating conditions. Water currents were measured using Acoustic Doppler Current Profilers (ADCP's) and Acoustic Doppler Velocimeters (ADV's) directly upstream and downstream of the turbine for inflow, wake and turbulence measurements. Electrical power output was measured using the Voltsys rectifier and the Shark power meter. Shaft speed was calculated based on the Voltsys measurements of the permanent magnet three phase generator AC generation frequency, coupled directly to the cross flow turbine under test (i.e., no gear box). Platform motions were captured using a Yost IMU (inertial measurement unit). Turbine thrust loading was measured using a reaction arm about the turbine deployment platform spanning beam, where two bi-directional load cells were connected to the system via a pinned connection. This submission includes zipped folders for each instrument containing quality controlled (QC'd) data in daily .csv files for the relevant duration specific to each instrument, along with separate .csv file that contains the units for each variable. Some instrument daily files are quite large and can pose a challenge for a visual spreadsheet editor to open. A processing software like MATLAB or Python is recommended. Note the degree of QC varied between each instrument due to time constraints. Particular time and attention was given to perform quality control tests on the acoustic based instruments that are particularly suscep...

16 TIDAL AND WAVE POWER↗

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects by Incorporating Machine Learning Techniques

The overall objective of this project is to develop a framework for reliable characterization and prediction of the state of stress in the overburden and underburden (including the basement) in CO 2 storage reservoirs using machine learning and integrated geomechanics and geophysical methods. Specifically, we propose to develop workflow encompassing of technologies and/or methods to predict stress and pressure changes due to CO 2 injection in an active tertiary recovery site and their impacts on subtle fault activation, fractures and occurrence of microseismic events and compare responses to field observations. In this project, we anticipate using dataset from the Farnsworth field Unit (FWU) which is operated by Purdure Petroleum. A novel elastic-waveform VSP inversion technique will be used to estimate high-resolution spatial and temporal changes of elastic moduli in CO 2 storage reservoirs, which will be combined with velocity-stress relationship derived from laboratory tests to obtain subsurface pressure and stress. Clustered microseismic data will be jointly inverted for improved focal mechanisms. Least-squares reverse-time migration of microseismic waveform data will be performed to directly image fracture/fault zones. Additionally, a deep neural network machine learning technique with convolutional and recurrent layers will be used for learning the spectro-temporal structures in microseismic waveforms. The results of this geotechnical data analysis will be integrated to develop a high-resolution 3D mechanical earth model extending from the overburden sealing formations to the underburden including the basement. Mechanical properties will be derived through integration of mechanical logs, tests, available results from chemo-mechanical laboratory tests, and elastic inversion of seismic data using a combination of Bayesian and stochastic methods as well as machine learning technique. Failure features (faults/fractures) will be represented and/or modeled based on seismic and core data analysis. A transient hydrodynamic-geomechanical model will be developed through coupling with the calibrated FWU reservoir simulation model. The full physics coupled model will be used to train a reduced order proxy model using machine learning algorithm for estimating stress which will then be used with appropriate constitutive relationships and forward seismological models to simulate pressure changes and induced microseismicity. An advanced optimization framework will be developed to perform a history match to minimize error between field observations and simulated. The history matched proxy model will be verified against the full-physics equivalent. The field observations that will be used in the coupled model calibration process include pressure/stress inverted from VSP, moment magnitude from microseismic analysis, real time downhole pressure measurements, production and injection data. Parameter sensitivity and uncertainty analysis will be performed to characterize the impact of model parameter uncertainty on stress estimates. The proposed project will have significant impact on future field implementation of the proposed technology. Because the project field site is an ongoing CO 2 EOR development, the value of the new technology will be demonstrated in an operational context and evaluated as a viable risk mitigation strategy. Cost/benefit will be evaluated together with the various commercial incentives for CO 2 sequestration available to oil and gas operators. The extensive available dataset and ongoing data acquisition under the SWP Phase III work plan provides flexibility for investigation of multiple approaches and reduces technical risk.

58 GEOSCIENCES↗

A method to determine the electric field of liquid argon time projection chambers using a UV laser system and its application in MicroBooNE

Liquid argon time projection chambers (LArTPCs) are now a standard detector technology for making accelerator neutrino measurements, due to their high material density, precise tracking, and calorimetric capabilities. An electric field (E-field) is required in such detectors to 2020 JINST 15 P07010 drift ionization electrons to the anode where they are collected. The E-field of a TPC is often approximated to be uniform between the anode and the cathode planes. However, significant distortions can appear from effects such as mechanical deformations, electrode failures, or the accumulation of space charge generated by cosmic rays. The latter effect is particularly relevant for detectors placed near the Earth’s surface and with large drift distances and long drift time. To determine the E-field in situ, an ultraviolet (UV) laser system is installed in the MicroBooNE experiment at Fermi National Accelerator Laboratory. The purpose of this system is to provide precise measurements of the E-field, and to make it possible to correct for 3D spatial distortions due to E-field non-uniformities. Here we describe the methodology developed for deriving spatial distortions, the drift velocity and the E-field from UV-laser measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Quasi-Steady Simulation of Glaze Ice Accretion and Heat Transfer in the Supercooled Large Droplet Regime in Atmospheric Turbulent Air Flow

Abstract To avoid the use of a computationally intensive full unsteady simulation while providing an accurate solution, a quasi-steady simulation has been performed to study glaze and mixed ice accretion in the supercooled large droplet (SLD) regime in turbulent flow of atmospheric air. We have attempted to find the minimum time-step necessary to adequately simulate the icing process on an airfoil surface. Based on node displacement, a mesh morphing scheme has been adopted in the computations to account for the moving boundaries that are caused by the continuous ice buildup. We have modeled ice accretion on an airfoil surface for a time period of 232 s using several time steps. At each time-step we solved the steady-state conservation equations for the air and droplet phases and then used the results as initial conditions for the next time position. The magnitude of the time-step ranged from a least accurate two-shot simulation (where the time-step is 116 s) to a most accurate 2320-shot simulation (where the time-step is 0.1 s). In-between these two extreme time steps, we have performed a three-shot simulation (where the time-step is 77.33 s), a four-shot simulation (where the time-step is 58 s), a six-shot simulation (where the time-step is 38.67 s), a 46-shot simulation (where the time-step is 5 s), and a 232-shot simulation (where the time-step is 1 s). We have done so to find out the degree of accuracy (or inaccuracy) of the multishot simulation approach and to find out the appropriate time-step needed for a successful and valid quasi-steady simulation. A valid quasi-steady simulation needs to use a time-step that is small enough to reproduce the full time-dependent solution within a very small error band. We have found that both the 1 s and the 0.1 s time steps produce virtually identical results. This is the primary litmus test that proves the validity of the quasi-steady-state assumption. The results adopted in this paper are thus all based on the more conservative 0.1 s time-step. In the process of performing the simulations, remeshing was required in order to maintain the grid density in zones of high curvature to be able to capture the full physics in those zones. After successfully modeling glaze ice accretion over the airfoil surface using the 0.1 s quasi-steady simulation approach, the effects of supercooled large droplets (SLDs) impacting the surface have been examined and presented in terms of the variation of the local collection efficiency, the water film thickness, and the heat transfer rate. Examination of the variation of the angle that the ice horn makes with the airfoil chord line demonstrated a 20% improvement in angle prediction when the time-step is reduced from 116 s to 0.1 s. The analysis also reveals a 12% or 8.5% increase in the maximum collection efficiency, βmax, depending on whether the value of the liquid water content (LWC) has been doubled from 0.5 g/m3 to 1 g/m3, or the value of the freestream velocity has been doubled from 75 m/s to 150 m/s, respectively. Because of the need to monitor the local collection efficiency and convective heat fluxes at each shot (regardless of the number of shots employed), the approach adopted here was found to be effective in successively and successfully reproducing the curvature of the glaze ice horn.

Engineering↗

Symmetrization of Strong Hydrogen Bond under High Pressure in Bihydroxide-Ion-Containing NaCu 2 (SO 4 ) 2 ·H 3 O 2 Revealed by Experimental Charge Density, Single-Crystal Electron Diffraction, and Neutron Diffraction Studies

In minerals and inorganic compounds, strong hydrogen bonding can lead to the formation of complex ionic species such as the H 3 O 2 – bihydroxide anion and Zundel cation H 5 O 2 + . We studied [NaCu 2 (SO 4 ) 2 ·H 3 O 2 ] natrochalcite, which contains bihydroxide anions and undergoes hydrogen bond symmetrization at the lowest pressure reported so far among inorganic compounds. Hydrogen bond symmetrization leads to changes in the bulk modulus, seismic wave velocities, and proton mobility and plays a primary role in high-temperature superconductivity, but its characteristics are not well understood due to a lack of systematic studies and limitations of experimental methods sensitive to this subtle change. In this work, we applied experimental charge density analysis based on in situ single-crystal X-ray diffraction data, along with the single-crystal neutron and electron diffraction experiments, to probe the behavior of hydrogen atoms during the hydrogen bond symmetrization process under high-pressure conditions. On the way to the symmetrical H-bonding, natrochalcite undergoes a series of complex redistributions of electron density, which we trace with multipole refinement and detailed analysis of changes in the Laplacian of electron density values. Additionally, we deconvoluted the equation of state (volume of the unit cell vs pressure relation) into the atomic equation of states describing dependencies of atomic charges or volumes vs pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Photoinitiated Reactions of Molecules and Radicals in Molecular Beams

The UV photochemistry of organic molecules is a fundamental process that governs reactions in the atmosphere, synthetic chemistry, processing of organic aerosols, and biological damage in living tissues. In many organic molecules, the ensuing evolution involves pathways that are in competition, giving rise to different products, isomerization, coupling to other electronic states, and secondary reactive collisions. Because photodissociation is usually fast (picoseconds to microseconds), these processes are far from equilibrium and are controlled by kinetic competition and dynamical forces. The work accomplished was focused primarily on the photochemistry of alpha-keto carboxylic acids, a group of acids produced from natural sources, which are implicated in aerosol formation and biological processes. In the atmosphere, they are destroyed mainly by solar radiation. Studying their photochemistry has been surprisingly difficult because of the complexity of their excited electronic states, and the effect of collisions and secondary reactions. The second project investigated the lifetime of excited electronic states of pyrazine and picoline and had both experimental and theoretical aspects. The work focused on a model molecule, pyruvic acid (PA), because it was hypothesized that it should be a good source of the unstable carbene, methylhydroxycarbene (MHC). PA has an internal hydrogen bond that controls the evolution of its decomposition. The decomposition was studied following excitation to two of its lowest excited states reached by laser irradiation at 351 and 193 nm. The photodissociation dynamics in the absence and presence of collisions was monitored and compared. Understanding the UV photochemistry requires the use of complementary experimental approaches. Two methods were enlisted, which together generated a comprehensive and detailed set of results: (i) The time-sliced velocity map imaging (SVMI) instrument at USC was exploited to determine kinetic energy release of fragments, fast dissociation timescales, and internal state distributions of fragments for which Resonance Enhanced Multiphoton Ionization (REMPI) schemes exist; and (ii) The multiplexed photoionization mass spectrometer (MPIMS) setup developed at the Sandia Combustion Research Facility was used for product discovery, achieved by exploiting tunable narrowband VUV radiation at the Advanced Light Source (ALS), and to follow in real time their subsequent unimolecular and bimolecular reactions. The experimental conditions ranged from collisionless molecular beams to study nascent products to flow reactors of variable pressures to study subsequent bimolecular reactions of the products. The goals stated above were successfully accomplished. By exciting PA to its lowest excited state, MHC was identified as the only primary product and its isomerization to vinylalcohol and acetaldehyde was directly observed in real time. Moreover, we reported the first bimolecular reaction of MHC, which was with acetaldehyde, and identified its reaction product. This paves the way for the study of other reactions of this important carbene intermediate. When excitation was carried out at 193 nm, which imparted much higher energy to PA, many more products were directly observed in real time, some deriving from three-body dissociation. Quantitative branching ratios and secondary reactions of radical products were also determined and analyzed. The studies of pyrazine and picoline focused on the second excited state of these molecules and showed (via ionization studies) that their excited states were very short lived (<100 fs) because of efficient couplings to lower electronic states via vibronic coupling. Theoretical work on modeling the absorption spectrum and decay mechanisms of the excited states are in progress in collaboration with theoreticians.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗