Search NASA⌕ Search

SEARCH · Search NASA

Results for “vertical 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 451 records · Page 25

Transition from wall modes to multimodality in liquid gallium magnetoconvection

Coupled laboratory-numerical experiments of Rayleigh-Bénard convection in liquid gallium subject to a vertical magnetic field are presented. The experiments are carried out in two cylindrical containers with diameter-to-height aspect ratio Γ = 1.0 and 2.0 at varying thermal forcing (Rayleigh numbers 10 5 ≲ Ra ≲ 10 8 ) and magnetic field strength (Chandrasekhar numbers 0 ≲ Ch ≲ 3 × 10 5 ). Laboratory measurements and numerical simulations confirm that magnetoconvection in our finite cylindrical tanks onsets via nondrifting wall-attached modes, in good agreement with asymptotic predictions for a semi-infinite domain. With increasing supercriticality, the experimental and numerical thermal measurements and the numerical velocity data reveal transitions between wall mode states with different azimuthal mode numbers and between wall-dominated convection to wall and interior multimodality. These transitions are also reflected in the heat transfer data, which combined with previous studies connect onset to supercritical turbulent behaviors in liquid metal magnetoconvection over a large parameter space. The gross heat transfer behaviors between magnetoconvection and rotating convection in liquid metals are compared and discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Integrating Experiments and Well Logs to Predict Caney Shale Static Mechanical Properties During Production with Supervised Machine Learning

Caney shale is one of the emerging oil reservoirs in Oklahoma. Understanding the impact of effective stress on its mechanical properties is critical for predicting hydraulic fracture geometry and overall hydrocarbon production. The objective of our study is to evaluate the impact of effective stress on the dynamic Young’s modulus using ultrasonic velocity measurements for Caney shale samples. A triaxial cell was utilized to measure ultrasonic (P-wave and S-wave) velocities for ten downhole Caney shale samples under various effective stresses to indirectly assess the impact of pore pressure change. The dynamic Young’s moduli estimated from these measurements were integrated with available conventional well logs (excluding sonic logs) and triaxial test results from Benge et al. (2021) to predict the static Young’s modulus using Random Forest (RF) and Extreme Gradient Boosting (XGBoost) models. The results showed that the estimated dynamic Young’s moduli from ultrasonic measurements were higher than the corresponding static Young’s modulus of cores from the same vertical well at similar depths. With increasing effective stress, the dynamic Young’s modulus increased for all samples. The estimated dynamic-to-static correction factor tended to be higher in zones of high neutron porosity (PHIN) and low density compared to other zones. Finally, SHapley Additive exPlanations (SHAP) for RF and XGBoost models identified depth, gamma ray (GR), and PHIN as key features for predicting the static Young’s modulus. This study enhances our understanding of the dynamic and static Young’s moduli for the Caney shale interval, as a function of effective stress and conventional well logs. The findings from this study can improve predictions of production throughout the well's lifespan by offering insights into the mechanical property degradation resulting from pore pressure depletion.

Kholy, Sherif M.↗

Tidal Energy Resource Characterization, Bottom Lander Measurements, Cook Inlet, AK, 2021

These datasets are from tidal resource characterization measurements collected on the Terrasond High Energy Oceanographic Mooring (THEOM) from 1 July 2021 to 30 August 2021 (60 days) in Cook Inlet, Alaska. The lander was deployed at 60.7207031 N, 151.4294998 W in ~50 m of water. The dataset contains raw and processed data from the following two instruments: 1. A Nortek Signature 500 kHz acoustic Doppler current profiler (ADCP). Data were recorded in 4 Hz in the beam coordinate system from all 5 beams. Processed data has been averaged into 5 minutes bins and converted to the East-North-Up (ENU) coordinate system. 2. A Nortek Vector acoustic Doppler velocimeter (ADV). Data were recorded at 8 Hz in the beam coordinate system. Processed data has been averaged into 5 minutes bins and converted to the Streamwise - Cross-stream - Vertical (Principal) coordinate system. Turbulence statistics were calculated from 5-minute bins, with an FFT length equal to the bin length, and saved in the processed dataset. Data was read and analyzed using the DOLfYN (version 1.0.2) python package and saved in MATLAB (.mat) and netCDF (.nc) file formats. Files containing analyzed data (".b1") were standardized using the TSDAT (version 0.4.2) python package. NetCDF files can be opened using DOLfYN (e.g., `dat = dolfyn.load(''*.nc")`) or the xarray python package (e.g. `dat = xarray.open_dataset("*.nc"). All distances are in meters (e.g., depth, range, etc), and all velocities in m/s. See the DOLfYN documentation linked in the submission, and/or the Nortek documentation for additional details.

16 TIDAL AND WAVE POWER↗

Characterization of beam ion loss in high poloidal beta regime on EAST

A critical issue for achieving the integrated operation of steady-state long-pulse high-confinement (H-mode) plasmas on experimental advanced superconducting tokamak (EAST) is to improve beam ion population confinement during neutral beam injection (NBI). To study the characterization of beam ion loss and improve beam ion confinement, the steady-state long pulse scenario discharges were conducted on EAST (β p $\geqslant$ 2.0, β N $\geqslant$ 1.7, q 95 $\geqslant$ 6.7 and H 98y2 $\geqslant$ 1.1) with NBI heating. Based on neutron yield, the beam voltage and line-averaged electron density were adjusted from 50 kV to 60 kV and 4.4 × 10 19 m –3 to 5.0 × 10 19 m –3 , respectively. The results show that the dominant mechanisms of beam ion loss are shine-through loss, prompt loss, and stochastic ripple loss. The shine-through loss fraction is determined by initial velocity, flight time and entire beam path. The change in prompt loss fraction is caused by the change in the deposition of beam ions. The change in stochastic ripple loss fraction is caused by the change in the initial fraction of trapped-confined ions. Detailed physics shows that the prompt loss fraction during counter-Ip injections (~45%) is far larger than during co-Ip injections (~5%) due to the finite orbit width. The lost ions are mainly deposited on the lower divertor or below the midplane since the direction of magnetic drift is vertical down. The orbit types of prompt loss during counter-Ip injections are mainly trapped-lost and ctr-passing lost. To minimize the prompt loss fraction during counter-Ip injections, a reversed Ip configuration (rev-Ip) discharge #94758 was conducted. The result suggests that the beam ion wall load fraction during counter-Ip tangential injection (~3%) is far lower than that in normal Ip configuration (nor-Ip) discharge #94820. It is also found that the confinement of beam ion population in the counter-Ip injection #94758 was greatly improved when compared to #94820. Finally, this study can provide unique support for the improvement of beam ion population confinement and for the performance evaluation of the NBI system on EAST and future tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Designing a Particle Flow Control Apparatus

Flow control within a particle-based Concentrated Solar Power (CSP) system is essential in determining the heat transfer coefficient, and therefore, the power generation capability of these systems. There are three areas where particle flow control is significant: the receivers, storage tanks, and particle-sCO2 heat exchangers. The focus of this work is on designing a new mechanism to control the flow in the particle-sCO2 heat exchangers due to the simplicity and potential cost savings when compared to the other areas of interest. The goal is for this new design to have quicker response times in terms of particle flowrate than a slide gate or flow control valve, which are designs currently used. The design resembles that of a chuck mechanism within a drill where a rotation of the sleeve elicits movement of the jaws both vertically and horizontally to close the outlet area of the nozzle. Additionally, this design will utilize the current actuator that is already used within these heat exchangers to reduce the complexity of implementation. The jaws are designed to be closed at an angle of 76˚ which is just slightly steeper than the hopper leading to the mechanism. Furthermore, this design can be tuned to limit particle bridging and other particle flow phenomena that result in blockages. The prototypes were 3D printed out of polylactic acid (PLA) and scaled up to 100%, 200%, and 400% to be able to observe the velocity profiles of the mechanism more clearly. Experiments are performed with this prototype to compare the inlet and outlet mass flow rates at different configurations of the jaws. The particles used in these experiments are 0.3mm HSP 40/70 that are commonly used in particle-based CSP systems.

14 SOLAR ENERGY↗

DESIGNING A PARTICLE FLOW CONTROL APPARATUS

Flow control within a particle-based Concentrated Solar Power (CSP) system is essential in determining the heat transfer coefficient, and therefore, the power generation capability of these systems. There are three areas where particle flow control is significant: the receivers, storage tanks, and particle-sCO2 heat exchangers. The focus of this paper is on designing a new mechanism to control the flow in the particle-sCO2 heat exchangers due to the simplicity and potential cost savings when compared to the other areas of interest. The goal is for this new design to have quicker response times in terms of particle flowrate than a slide gate or flow control valve, which are designs currently used. The design resembles that of a chuck mechanism within a drill where a rotation of the sleeve elicits movement of the jaws both vertically and horizontally to close the outlet area of the nozzle. Additionally, this design will utilize the current actuator that is already used within these heat exchangers to reduce the complexity of implementation. The jaws are designed to be closed at an angle of 76 ̊ which is just slightly steeper than the hopper leading to the mechanism. Furthermore, this design can be tuned to limit particle bridging and other particle flow phenomena that result in blockages. The prototypes were 3D printed out of polylactic acid (PLA) and scaled up to 100%, 200%, and 400% to be able to observe the velocity profiles of the mechanism more clearly. Experiments are performed with this prototype to compare the inlet and outlet mass flow rates at different configurations of the jaws. The particles used in these experiments are 0.3mm HSP 40/70 that are commonly used in particle-based CSP systems.

14 SOLAR ENERGY↗

W-Band ARM Scanning Cloud Radar (WSACR) 2nd Generation CF-Radial Spectral Data, Vertically-Pointing Scan, Cross Mode (a1)

ARM's scanning cloud radars are fully coherent dual-frequency, dual-polarization Doppler radars mounted on a common scanning pedestal. Each pedestal includes a Ka-band radar (2kW peak power) and the deployment location determines whether the second radar is a W-band (WSACR; 1.7 kW peak power) or X-band (XSACR; 20 kW peak power). Beamwidths at Ka and W bands are roughly matched at 0.3 degrees. Due to the narrow antenna beamwidth, ARM’s scanning cloud radars use scanning strategies that are unlike typical weather radars. Rather than focusing on plan position indicator, or PPI, scans, the Ka-SACR uses range height indicator, or RHI, scans at numerous azimuths to obtain cloud volume data. Measurements collected with the W-SACR are copolar and cross-polar radar reflectivity, Doppler velocity, spectra width and spectra when not scanning, and linear depolarization ration.

54 ENVIRONMENTAL SCIENCES↗

Computational fluid dynamics analysis of char conversion in Sandia’s pressurized entrained flow reactor

Design and analysis of practical reactors utilizing solid feedstocks rely on reaction rate parameters that are typically generated in lab-scale reactors. Evaluation of the reaction rate information often relies on assumptions of uniform temperature, velocity, and species distributions in the reactor, in lieu of detailed measurements that provide local information. This assumption might be a source of substantial error, since reactor designs can impose significant inhomogeneities, leading to data misinterpretation. Spatially resolved reactor simulations help understand the key processes within the reactor and support the identification of severe variations of temperature, velocity, and species distributions. In this work, Sandia’s pressurized entrained flow reactor is modeled to identify inhomogeneities in the reaction zone. Tracer particles are tracked through the reactor to estimate the residence times and burnout ratio of introduced coal char particles in gasifying environments. The results reveal a complex mixing environment for the cool gas and particles entering the reactor along the centerline and the main high-speed hot gas reactor flow. Furthermore, the computational fluid dynamics (CFD) results show that flow asymmetries are introduced through the use of a horizontal gas pre-heating section that connects to the vertical reactor tube. Computed particle temperatures and residence times in the reactor differ substantially from the idealized plug flow conditions typically evoked in interpreting experimental measurements. Furthermore, experimental measurements and CFD analysis of heat flow through porous refractory insulation suggest that for the investigated conditions (1350 °C, <20 atm), the thermal conductivity of the insulation does not increase substantially with increasing pressure.

47 OTHER INSTRUMENTATION↗

W-Band ARM Scanning Cloud Radar (WSACR) 2nd Generation CF-Radial Spectral Data, Vertically-Pointing Scan, Co-Polarization Mode (a1)

ARM's scanning cloud radars are fully coherent dual-frequency, dual-polarization Doppler radars mounted on a common scanning pedestal. Each pedestal includes a Ka-band radar (2kW peak power) and the deployment location determines whether the second radar is a W-band (WSACR; 1.7 kW peak power) or X-band (XSACR; 20 kW peak power). Beamwidths at Ka and W bands are roughly matched at 0.3 degrees. Due to the narrow antenna beamwidth, ARM’s scanning cloud radars use scanning strategies that are unlike typical weather radars. Rather than focusing on plan position indicator, or PPI, scans, the Ka-SACR uses range height indicator, or RHI, scans at numerous azimuths to obtain cloud volume data. Measurements collected with the W-SACR are copolar and cross-polar radar reflectivity, Doppler velocity, spectra width and spectra when not scanning, and linear depolarization ration.

54 ENVIRONMENTAL SCIENCES↗

W-Band ARM Scanning Cloud Radar (WSACR) 2nd Generation CF-Radial Spectral Data, Vertically-Pointing Scan, Cross-Polarization Model (a1)

ARM's scanning cloud radars are fully coherent dual-frequency, dual-polarization Doppler radars mounted on a common scanning pedestal. Each pedestal includes a Ka-band radar (2kW peak power) and the deployment location determines whether the second radar is a W-band (WSACR; 1.7 kW peak power) or X-band (XSACR; 20 kW peak power). Beamwidths at Ka and W bands are roughly matched at 0.3 degrees. Due to the narrow antenna beamwidth, ARM’s scanning cloud radars use scanning strategies that are unlike typical weather radars. Rather than focusing on plan position indicator, or PPI, scans, the Ka-SACR uses range height indicator, or RHI, scans at numerous azimuths to obtain cloud volume data. Measurements collected with the W-SACR are copolar and cross-polar radar reflectivity, Doppler velocity, spectra width and spectra when not scanning, and linear depolarization ration.

54 ENVIRONMENTAL SCIENCES↗

Stepwise Dynamic Calibration of a Hydromechanical Simulation Using Time-Lapse Vertical Seismic Profile

This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse Vertical Seismic Profiling (VSP) integration. The selected study site is the 13-10A injector well within the ongoing CO2-EOR operation of the Farnsworth Field Unit. The Time-lapse VSP dataset carries the combined effects of fluid substitution and mean effective stress changes, thereby providing a dataset amenable for the calibration of production and injection-induced stress changes. The concept is similar to calibrating a reservoir simulation model in that the process honor real field data to set up an inverse problem. The solution optimizes the independent and impactful geomechanical parameters that replicate the observed time-lapse seismic velocity changes. This stress calibration is enabled by 4D geomechanical modeling and the VSP Integration workflow. This calibration benefits from extensive geological, geophysical and geomechanical characterization through 3D seismic data, geophysical well logs, and core assessed as part of the 1D MEM conducted on the 13-10A subject well. These data are used to develop a site-specific rock physics model. The Biot Gassmann workflow combines rock physics and reservoir simulation outputs to determine the fluid substitution contribution to seismic velocity change. Additionally, modeled seismic velocity attributed to mean effective stress are determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from the ultrasonic seismic velocity measurements on the extracted Morrow B core. A penalty function is then formed between the modeled seismic velocities and the observed time-lapse VSP dataset. Four independent and impactful geomechanical parameters have been determined. These are the bulk modulus and shear modulus for zero porosity and the shear and compressional seismic velocity to mean effective stress derivatives. The dataset of numerous coupled hydromechanical- geomechanical simulation realizations is built by combining variations of the four stated geomechanical parameters. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer are used to converge on the optimal geomechanical parameters. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling.

02 PETROLEUM↗

Spectral Element Modeling of Acoustic to Seismic Coupling Over Topography

Acoustic waves in the atmosphere are commonly recorded on seismometers as they couple into the ground. These signals, here called ground coupled airwaves, are not commonly considered in numerical modeling of infrasound propagation, which often assumes a rigid unmeshed boundary. Starting from an analytically-tractable spherical wave model, we analyze how the coupling of an acoustic wave into a planar elastic halfspace changes over a wide range of scenarios. We use energy admittance to quantify acoustic to seismic coupling over both a planar elastic halfspace and meshed topography. Our spectral element and analytic calculations have different maxima as a function of incidence angle, with very high admittance values for near-vertical incidence (maximum ≈78%). Energy admittance calculations at shallow incidence angles are much smaller (less than 1%). In simulations over the complex topography of Sakurajima Volcano, we attribute the variable spatial pattern of energy admittance to changes in earth parameters between each model. The observed pressure difference over the simulated 15 km region appears to be < 2%. While this value is relatively small, the cumulative addition over 100s of km and multiple acoustic bounce points may be significant. Acoustic to seismic coupling along the propagation path may bias long distance yield estimates, particularly when infrasound propagates over regions with steep topography or particularly slow seismic velocities, such as alluvial planes.

58 GEOSCIENCES↗

Dark matter raining on DUNE and other large volume detectors

Direct detection is a powerful means of searching for particle physics evidence of dark matter (DM) heavier than about a GeV with o(kiloton) volume, low-threshold detectors. In many scenarios, some fraction of the DM may be boosted to large velocities enhancing and generally modifying possible detection signatures. We investigate the scenario where 100% of the DM is boosted at the Earth due to new attractive long-range forces. This leads to two main improvements in detection capabilities: (1) the large boost allows for detectable signatures of DM well below a GeV at large-volume neutrino detectors, such as DUNE, Super-K, Hyper-K, and JUNO, as possible DM detectors, and (2) the flux at the Earth’s surface is enhanced by a focusing effect. In addition, the model leads to a significant anisotropy in the signal with the DM flowing dominantly vertically at the Earth’s surface instead of the typical approximately isotropic DM signal. We develop the theory behind this model and also calculate realistic constraints using a detailed GENIE simulation of the signal inside detectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Tidal Energy Resource Characterization, Velocity and Turbulence Measurements, Processed Data, Cook Inlet, AK, 2021

This submission contains processed datasets from a long-term deployment of 3 moorings and a transect survey of the proposed tidal energy site off the East Forelands in Cook Inlet, AK. The long-term mooring datasets were created from 8 instruments mounted on a Terrasond High Energy Oceanographic Mooring (THEOM) bottom lander and two Mid-Water Mooring (MWM) Stablemoor buoys from 1 July 2021 to 31 August 2021 (60 days). The west-most mooring (MWM1) was deployed at 60.720225 N, 151.436196 W in ~50 m of water. The middle mooring (THEOM) was deployed at 60.720703 N, 151.429500 W in ~52 m of water. The east-most buoy (MWM2) was deployed at 60.720081 N, 151.420896 W in ~50 m of water. Each Stablemoor carried three instruments: 1. A Nortek Vector acoustic Doppler velocimeter (ADV) mounted at the Stablemoor's nose. Data were recorded at 8 Hz on a 5 minute duty cycle every 20 minutes. Data was motion-corrected using the internal IMU and external ADCP bottom-track data and then bin-averaged into 4 minute bins and converted to the Principal (streamwise, cross-stream, vertical) coordinate system. (Note: 30 seconds were trimmed from the beginning and end of each 5 minute duty cycle to account for the filter end-effects from turning on and turning off the IMU.) 2. A down-looking Nortek Signature 1000 kHz acoustic Doppler current profiler (ADCP) mounted in the first Stablemoor instrument well. Data were recorded in 2 Hz with 5-beam burst and bottom-track enabled. Processed data has been averaged into 10 minute bins and converted into the Principal coordinate system. 3. An up-looking Nortek Signature 1000 kHz acoustic Doppler current profiler (ADCP) mounted in the second Stablemoor instrument well. Data were recorded at 4 Hz with 5 beam burst enabled. Processed data has been averaged into 10 minute bins and converted into the Principal coordinate system. Note: the down-facing ADCP on MWM1 failed on July 10th, 2021, only recording 9 days of data. Because ADV motion-correction required bottom track, the ADV from MWM1 also only has 9 days processed. Additionally, only 25 days of data were processed from the MWM2 ADV because it appeared to have been impacted by debris on 7/25. Two instruments were mounted on the THEOM (see MHKDR link further below for THEOM raw data): 4. A Nortek Vector acoustic Doppler velocimeter (ADV). Data were recorded at 8 Hz on a 5 minute duty cycle every 20 minutes. Data was bin-averaged into 5 minute bins, and converted to the Principal coordinate system. 5. A Nortek Signature 500 kHz acoustic Doppler current profiler (ADCP). Data were recorded in 4 Hz in the beam coordinate system from all 5 beams. Processed data has been averaged into 10 minutes bins and converted to the Principal coordinate system.

16 TIDAL AND WAVE POWER↗

A Moving-Wave Implementation in WRF to Study the Impact of Surface Water Waves on the Atmospheric Boundary Layer

Abstract While numerous modeling studies have focused on the interaction of ocean surface waves with the atmospheric boundary layer, most employ idealized waves that are either monochromatic or synthetically generated from a theoretical wave spectrum, and the atmospheric solvers are typically incompressible. To study wind–wave coupling in real-world scenarios, a model that can simulate both realistic meteorological and wave conditions is necessary. In this paper we describe the implementation of a moving bottom boundary condition into the Weather Research and Forecasting Model for large-eddy simulation applications. We first describe the moving bottom boundary conditions within WRF’s pressure-based vertical coordinate system. We then validate our code with idealized test cases that have analytical solutions, including flow over a monochromatic wave with and without viscosity. Finally, we present results from turbulent flows over a moving monochromatic wave with different wave ages, and demonstrate satisfactory agreement of the wave growth rate with results from the literature. We also compare atmospheric stress and wind parameters from two physically equivalent cases. The first specifies a wind moving in the same direction as a propagating wave, while the second involves a stationary wave with the wind adjusted such that the wind relative to the wave is the same as in the first case. Results indicate that the velocity and Reynolds stress profiles for the two cases match, further validating the moving bottom implementation.

17 WIND ENERGY↗

An intercomparison of wall fluxes in a turbulent thermal convection chamber: Direct numerical simulations and wall-modeled large-eddy simulations enhanced by machine learning

Thermal convection in a closed chamber is driven by a warm bottom, a cold top, and side walls at various temperatures. Although wall fluxes are the source of convection energy, accurately modeling these fluxes (i.e., the wall model) is challenging. In large-eddy simulations (LESs), many wall models are traditionally derived from the canonical boundary layer, which may be unsuitable for thermal convection bounded by both horizontal and vertical walls. This study conducts a model intercomparison of dry convection in a cubic-meter chamber using three direct numerical simulations (DNSs) and four LESs with different wall models. The LESs employ traditional wall models, a new wall model employing physics-aware neural networks, and a refined grid near the walls. The experiment involves four cases with varying sidewall temperatures. Our results show that LESs capture the main flow features and the trends of mean fluxes. The physics-aware neural networks and refined wall grids can improve the temporally averaged local fluxes when the large-scale circulation has a preferred direction. Even without the local improvement of wall fluxes, the LES flow quantities (temperature and velocities) can still largely match those in DNSs, provided the mean flux largely matches the DNSs. Additionally, DNSs reveal that a variation in corner treatments has minimal impacts on the flow quantities away from corners. Finally, LESs underestimate the mean fluxes of the entire wall due to their inability to resolve corner regions, but their mean flux away from the corner can better match DNS.

54 ENVIRONMENTAL SCIENCES↗

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

As part of the Midwest Regional Carbon Sequestration Partnership (MRCSP) Phase III project, a monitoring study was 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.

47 OTHER INSTRUMENTATION↗