SEARCH · Engineering Papers
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.
An enrichment of sulfidic sediments from Zodletone spring was sequenced as a metagenome. Draft genomes representing Cloacimonadota, Deltabacterota, Firmicutes, and Patescibacteria were binned and annotated and will aid functional genomics and cultivation efforts.
The Caney Shale is emerging as a target for hydrocarbon production, creating opportunities to study rock mechanical origins of observed/anticipated challenges to effective and sustained stimulation. Here we examine five subunits within the Caney, two of which are dubbed “ductile” and three of which are dubbed “reservoir” rock types. The “ductile” versus “reservoir” identification is initially based on elastic properties ascertained from well logs. By comparing core-based mechanical tests, it is found nominally ductile and reservoir layers do not differ in terms of brittleness, but instead the nominally ductile layers tend to be weaker and more prone to creep deformation compared to the nominally more brittle reservoir layers. This difference in mechanical properties is shown to generate higher susceptibility to proppant embedment in the weaker and more creep-prone layers. Specifically, over a 5-year period, the creep is expected to have little impact on reservoir layers but is predicted to reduce propped fracture aperture by a factor of 2 in the ductile layers.
This report describes the methodology and results of a state cost-effectiveness analysis of ASHRAE Standard 90.1-2016. Moving to the Standard 90.1-2016 edition from Standard 90.1-2013 is found to be cost-effective for the state. Annual energy cost savings, added construction costs, and life-cycle costs are all described and presented here for various building types and climate zones.
This report describes the methodology and results of a state cost-effectiveness analysis of ASHRAE Standard 90.1-2019. Moving to the Standard 90.1-2019 edition from Standard 90.1-2016 is found to be cost-effective for the state. Annual energy cost savings, added construction costs, and life-cycle costs are all described and presented here for various building types and climate zones.
This report is part of a publication series that focuses on approximately 100 different local geographies, or "clusters." Each report provides commercial and multifamily building characteristic and energy data for a local geography, with the intention of helping policymakers at the city, county and state levels better understand building energy use. Specifically, this report breaks down the building stock in the counties shown in Figure 1 by building type, size, energy consumption, and emissions.
Seismology is the study of earthquakes. It involves looking at the origin and propagation of seismic waves through the Earth to study the interior.
This report is part of the second phase of a publication series focusing on approximately 100 different local geographies, or "clusters." Each report provides characteristic features and energy data for commercial buildings in a specific area to help policy makers at the city, county, and state level better understand building energy use and emissions. This report breaks down the energy consumption and emissions of the building stock in the counties shown in Figure 2 by building type, building size, end use, energy consumption, emissions, and segment.
This report is an addendum to a publication series that focuses on approximately 100 different local geographies, or "clusters". This addendum expands the report series to include large multifamily building characteristics as well as energy and emission data for each local geography. The intention of this addendum is to help policymakers at the city, county, and state levels better understand building energy use and emissions in large multifamily buildings.
Explore the source record for details and available documents.
ABSTRACT: The ultimate aim of hydraulic fracturing is to have a long and conductive flow path that extends from the wellbore into the formation. The effective fracture length is part of a hydraulically propped fracture which contributes to production. The difficulty in achieving economical production targets from shale reservoirs is at the forefront in many exploration companies. Fracture conductivity loss is related to; proppant embedment under depletion, proppant crushing, damage as a result of fracturing fluid, fines migration and proppant-pack permeability-damage are some of the factors that contribute to production decline after hydraulic fracturing in shale reservoirs. The Caney Shale is a calcareous organic-rich mudrock. Various studies have investigated the effect that clay on shale well productivity, however, there is currently no literature on the Caney shale in relation to horizontal wells; all available literature exists in vertical wells as well as on formations of the Caney that are shallow in comparison to an emerging play which is twice the depth. In this paper we investigate stress-dependent fracture conductivity of proppant-filled fractures and proppant embedment in Caney shale through laboratory and modeling studies. API fracture conductivity tests were conducted using 2% KCl on five locations within the Caney shale that consisted of selecting three brittle(reservoir) zones and two ductile zones. Confining pressures range from 1,000 psi to 12,000 psi at 210°F. Conductivity, permeability as well as embedment were measured during the test. Our experimental results have confirmed that improved fracture conductivity is attributed to; proppant size, the increase in porosity of the proppant pack, closure pressure changes and the reduction in fracture conductivity are a function of many factors such as fracture closure stress. The findings from this study could help the stimulation design by providing new insights into the critical factors that are to be determined to facilitate the choice of proppants as well as fracturing fluids for long term production and recovery from shale reservoirs. 1 INTRODUCTION The development of low permeability formations, like shales, has been aided by hydraulic fracturing of horizontal wells (Radonjic et al., 2020). Hydraulic fracturing fluid is injected at a high pressure to induce tensile fractures that can link to and stimulate natural fractures (Katende et al., 2021a,b). Preserving adequate conductivity in hydraulic fractures over the life of the wells is required for economic production; nevertheless, conserving such conductivity can be difficult in some circumstances, particularly in soft, clay-rich formations (Wang et al., 2021). Proppant particles help to keep the fractures open when the pumping stops and the fracturing fluid returns to the wellbore, producing one or more propped hydraulic fractures of varying length, breadth, and height (Katende et al., 2021a). The proppant pack within the hydraulic fracture boosts well output by providing a greater permeability flowpath for hydrocarbons (Duenckel et al., 2016). Proppant in the fracture is under complicated stress conditions, and the interplay between the rock formation and the proppant pack has a significant impact on proppant-pack permeability (Karazincir et al., 2019). Proppant may be embedded (Katende et al., 2021a) in the rock or crushed into small pieces if the proppant size and strength characteristics are not specified appropriately, resulting in a loss in proppant-pack permeability and fracture aperture, and consequently a fall in well output.
Explore the source record for details and available documents.
SPICE (the SPARTAN and IMPROVE Comparison Experiment) aims to obtain and quantify comparisons between aerosol PM 2.5 mass concentration measurements from the University of Oklahoma (OU) Surface Particulate Matter Network (SPARTAN) station and the U.S. Environmental Protection Agency (EPA) Interagency Monitoring of Protected Visual Environments (IMPROVE) station hosted by the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility at ARM’s Southern Great Plains (SGP) observatory in Oklahoma. Aerosols—both natural and anthropogenic—affect human populations in multiple ways. Much of ARM’s research focuses on how aerosols influence weather and climate through their optical and radiative effects, as well as their impacts on clouds and precipitation. However, aerosols also pose direct risks to humans and other organisms through inhalation, with the severity of health impacts depending on particle size, chemical composition, and duration of exposure. The IMPROVE network was established by the U.S. EPA to monitor air quality, including visible clarity as well as total and chemically speciated aerosol mass concentrations. The ARM SGP site hosts the IMPROVE SOGP station. Separately, the SPARTAN network operates a globally distributed set of stations similar to IMPROVE but with an emphasis on remote deployment and semi-autonomous operation for use beyond the borders of the United States (IMPROVE only operates within the U.S.). The University of Oklahoma operates a SPARTAN station. To establish confidence in the OU SPARTAN instrumentation and measurement protocol relative to the EPA-certified IMPROVE station, the OU SPARTAN station is currently deployed at SGP in close proximity to the IMPROVE SOGP station. The SPICE campaign was envisioned as a contiguous calendar-year effort for 2025 to capture seasonal variation in mass loading as well as composition. However, independent of the SPICE campaign, the SPARTAN network adopted a new filter construction part-way through the year, interrupting our contiguous data set. Thus, to obtain a contiguous data set with a uniform consistent configuration, SPICE desires an extension through 2026.
This study explores the updraft and downdraft properties of mature stage mesoscale convective systems (MCSs) in terms of draft core width, shape, intensity, and mass flux characteristics. The observations use extended radar wind profiler (RWP) and surveillance radar datasets from the U.S. Department of Energy Atmospheric Radiation Measurement program for midlatitude (Oklahoma, USA) and tropical (Amazon, Brazil) sites. MCS drafts behave qualitatively similar to previous aircraft and RWP cloud summaries. The Oklahoma MCSs indicate larger and more intense convective up- and downdraft cores, and greater mass flux than Amazon MCS counterparts. However, similar size-intensity relationships and draft vertical profile behaviors are observed for both regions. Additional similarities include weak positive correlations between core intensity and core width (correlation coefficient r ~ 0.5), and increases in draft intensity with altitude. A model-observational intercomparison for draft properties (core width, intensity, mass flux) is also performed to illustrate the potential usefulness of statistical observed draft characterizations. Idealized simulations with the Weather Research and Forecasting model aligned with midlatitude MCS conditions are performed at model grid spacings (?x) that range from 4 km to 250 m. It is shown that the simulations performed at ?x = 250 m at similar mature MCS lifecycle stages are those that exhibit draft intensity, width, mass flux, and shape parameter performances best matching with observed properties.
Abstract This paper evaluates the representation of turbulence and its effect on transport and dispersion within multiscale and microscale-only simulations in an urban environment. These simulations, run using the Weather Research and Forecasting Model with the addition of an immersed boundary method, predict transport and mixing during a controlled tracer release from the Joint Urban 2003 field campaign in Oklahoma City, Oklahoma. This work extends the results of a recent study through analysis of turbulence kinetic energy and turbulence spectra and their role in accurately simulating wind speed, direction, and tracer concentration. The significance and role of surface heat fluxes and use of the cell perturbation method in the numerical simulation setup are also examined. Our previous study detailed the model development necessary for our multiscale simulations, examined model skill at predicting wind speeds and tracer concentrations, and demonstrated that dynamic downscaling from mesoscale to microscale through a sequence of nested simulations can improve predictions of transport and dispersion relative to a microscale-only simulation forced by idealized meteorology. Here, predictions are compared with observations to assess qualitative agreement and statistical model skill at predicting wind speed, wind direction, tracer concentration, and turbulent kinetic energy at locations throughout the city. We also investigate the scale distribution of turbulence and the associated impact on model skill, particularly for predictions of transport and dispersion. Our results show that downscaled large-scale turbulence, which is unique to the multiscale simulations, significantly improves predictions of tracer concentrations in this complex urban environment. Significance Statement Simulations of atmospheric transport and mixing in urban environments have many applications, including pollution modeling for urban planning or informing emergency response following a hazardous release. These applications include phenomena with spatial scales spanning from millimeters to kilometers. Most simulations resolve flow only within the urban area of interest, omitting larger scales of turbulence and regional influences. This study examines a method that resolves both the small and large-scale flow features. We evaluate simulation accuracy by comparing predictions with observations from an experiment involving the release of a tracer gas in Oklahoma City, Oklahoma, with emphasis on correctly modeling turbulent fluctuations. Our results demonstrate the importance of resolving large-scale flow features when predicting transport and dispersion in urban environments.
Seismicity at The Geysers geothermal field and in Oklahoma is heavily influenced by industrial activities related to energy production, though the mechanism in which earthquakes are induced or triggered is different. At The Geysers, much of the seismicity is linked to thermoelastic stresses caused by injecting cold water into hot rocks, whereas in Oklahoma, the seismicity is linked to a reduction of confining stress on faults due to increasing pore pressure resulting from wastewater injections. Here, in this study, we show that these contrasting conditions are also evident in tidally triggered earthquakes. At The Geysers, earthquakes preferentially occur during maximum extensional stress, which does not occur at the same time as maximum shear stress on optimally oriented faults in the regional stress field. In Oklahoma, earthquakes preferentially occur during maximum shear stress on optimally oriented faults rather than maximum extensional stress. The magnitude of tidal extensional stress is naturally much greater than tidal shear stress. However, in a fluid-saturated environment, pore pressure responds to changes in volume, which can counteract or reduce the effect of the applied stress. The difference in behavior at these two sites is indicative of the level of coupling between applied stress and pore pressure, corresponding to unsaturated conditions at The Geysers and high pore pressure in Oklahoma.
The lower atmosphere contributes directly to the modulation of weather and climate. Understanding the importance of this part of the atmosphere, scientists have worked to improve the representation of the atmospheric boundary layer in numerical prediction tools. Such work depends upon information from observing systems, including remote sensors, weather balloons, and research aircraft. Recent years have seen significant advances in uncrewed aircraft systems (UAS) for atmospheric research. Those efforts have provided new perspectives on atmospheric and surface conditions, particularly from smaller UAS (sUAS) platforms. This field campaign included collection of data using several sUAS, and comparing the sUAS data to those from different observing facilities at the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility’s Southern Great Plains (SGP) observatory in Oklahoma. Over the course of one week, the Black Swift Technologies and University of Nebraska teams conducted flights for the sole purpose of system evaluation and intercomparison. In addition, these teams conducted joint flight operations at the Marshall, Oklahoma Mesonet site alongside teams from the University of Colorado, University of Oklahoma, and Oklahoma State University to collect additional data for side-by-side comparison. Platforms operated by Black Swift Technologies and the University of Nebraska under this project included one fixed-wing and two rotary-wing platforms, as well as an instrumented surface vehicle. These include the Black Swift Technologies S0 UAS, the University of Nebraska M600 UAS, a University of Nebraska-operated Meteomatics Meteodrone UAS, and the University of Nebraska CoMeT (Combined Mesonet and Tracker) vehicle. Additional details on all of these platforms can be found in de Boer et al. (in prep). In total, 95 sUAS flights were conducted for a combined total of 18.6 flight hours. The S0 was operated both to the south and north of the SGP 60 m tower, while the M600 and Meteodrone were operated directly to the east of the tower. Weather conditions were generally good, with moderate winds. Flights were conducted both to follow the extra radiosondes that were launched as part of this campaign, as well as to conduct extended statistical sampling at tower instrument heights. The later flights also provided data for statistical platform intercomparison.
This DOE-funded joint project had the over-arching goal of understanding the reasons for observed size distributions of ice particles in cold clouds. Its approach involves the use of cloud models and field observations by radar and aircraft of real storms. The project addresses three major research objectives: (1) Utilize a novel polarimetric radar technique to retrieve size distributions and amounts of ice from radar data collected during the previous DOE ARM field campaigns; (2) Evaluate the results of ice microphysical retrievals using available in situ aircraft measurements and other remote sensors (e.g., vertically pointing cloud radars and wind profilers); (3) Improve the treatment of microphysical processes leading to cloud glaciation and ice multiplication in numerical cloud models. The study was performed by three research organizations: University of Oklahoma, The Hebrew University of Jerusalem, Israel, and Lund University, Sweden. The Israel and Swedish partners were funded by the sub-awards from the University of Oklahoma. Cloud modeling studies have been performed by the research teams at The Hebrew University of Jerusalem (HUJ) and Lund University to identify the origins of high concentrations of cloud ice in areas of high ice water content (HIWC). The Hebrew University Cloud Model (HUCM) with full spectral bin microphysics and the Lund University aerosol-cloud (AC) model with a hybrid bin / bulk microphysics scheme complementing HUCM were utilized for simulations. Both research teams had particular focus on secondary ice production (SIP) as one of the possible sources of enhanced ice concentration. The HUJ group suggested a novel concept of ice multiplication during droplet freezing. It is assumed that splintering and droplet fragmentation during droplet freezing takes place because of dendritic growth within a supercooled drop. The resulting simulations of SIP generated small ice in concentrations exceeding hundreds per liter similar to what was observed in the HIWC regions of the tropical storms. The Lund team explored the SIP mechanisms such as breakup of ice particles due to ice-ice collisions and ice sublimation that are expected to dominate the continental storms. They also quantified the impact of homogeneous nucleation of cloud droplets on the total number concentration of ice at very low temperatures near the tops of the clouds. Additionally, the Lund AC model is able to simulate the effect of aerosols of various types (including biological) on the cloud life cycle and the corresponding ice production. The HUCM / AC model was used to simulate one of the “golden” cases of the DOE MC3E campaign on 20 May 2011. The model output was converted into the fields of polarimetric radar variables using the polarimetric radar forward operator developed at the University of Oklahoma and compared with radar observations and in situ microphysical measurements onboard research aircraft. It was demonstrated that specific differential phase KDP is the best radar parameter to identify the HIWC areas and quantify the corresponding ice parameters. For the first time, the shape of the vertical profile of KDP was realistically reproduced by the cloud model with the KDP maximum in the dendritic growth layer (DGL) centered at the -15°C isotherm. The University of Oklahoma team has developed a methodology for polarimetric radar retrievals of such microphysical parameters of ice as ice water content (IWC), mean volume diameter (Dm), and total number concentration (Nt) of ice particles. These retrievals have been validated using in situ aircraft measurements during 6 field campaigns and proved to be quite robust and reliable. This allowed to build the first climatology of the vertical profiles of polarimetric radar variables and retrieved microphysical parameters for the three types of weather systems: continental MCSs, maritime MCSs, and tropical cyclones / hurricanes (Hu and Ryzhkov 2022). The data were collected by a multitude of the WSR-88D radars in 13 continental and 10 maritime MCSs and 11 landfalling hurricanes. The HIWC areas were identified within the examined storms and the corresponding “HIWC statistics” was compared with the “background” one without HIWC. An overarching conclusion of the study is that maritime tropical storms (MCSs and hurricanes) are characterized by smaller size ice in higher concentration compared to the continental MCSs. High ice water content in the HIWC areas is primarily caused by a strong jump in a number concentration of ice particles rather than the increase of their size compared to the “background” environment. This may point to the homogeneous nucleation of excessive amounts of supercooled droplets and / or secondary ice production as the possible origins of HIWC. Such a climatology provides a good observational reference for the modelers to evaluate the performance of their models. As an example, the in-depth analysis of the 20 May 2011 MC3E case shows that the advanced cloud models developed in the course of this study still tend to underestimate the number concentration of ice in the HIWC areas although they succeed in reproducing realistically looking vertical profiles of IWC and Nt. The results of the project research are summarized in 13 journal papers.
Abstract Boundary layer convergence lines (CLs) are highly effective at deep-convection initiation (DCI), suggesting that their associated updraft properties differ from those of more widespread turbulent updrafts in the planetary boundary layer (PBL). This study exploits observations at the Atmospheric Radiation Measurement Southern Great Plains (ARM SGP) observatory in Oklahoma from 2011 to 2016 to quantify CL properties and their relation to turbulent PBL eddies preceding CL arrival. Two independent methods for estimating CL properties are developed at two locations in the SGP region, both relying on the assumption of a 2D circulation in the CL-normal plane but using different combinations of instruments. The first (the radar method) relies mainly on scanning radar data and is applied to 61 CLs passing near a high-resolution scanning radar based in Nardin, Oklahoma, while the second (the surface method) relies mainly on surface wind data and is applied to 68 CLs crossing the SGP facility in nearby Lamont, Oklahoma. Mean daytime (1000–1900 LST) CL width (∼2 km) and convergence magnitude (∼0.003 s −1 ) are similar for both methods, and mean daytime CL depth is ∼0.75 km. The two methods disagree at night (0000–1000 and 1900–2400 LST), where the surface method estimates wider and weaker CLs than the radar method. This difference may stem from the radar beam overshooting the shallow, highly stable nocturnal PBL. The largest CL updrafts are slightly wider (∼20%) and stronger (∼40%) than the largest PBL updrafts in the pre-CL period, generating 50%–100% larger updraft mass fluxes over most of the PBL depth. Significance Statement Deep convection is commonly initiated by boundary layer convergence lines (CLs), which are associated with intense surface-based wind convergence and strong updrafts that may lift air to saturation. Although CLs form regularly, they are far less common than ordinary, short-lived turbulent thermals in the daytime boundary layer. To better understand why CLs are so effective at deep-convection initiation, we observationally quantify their morphologies and strengths and compare these properties to those of surrounding turbulent updrafts. Perhaps surprisingly, the CLs are found to exhibit only slightly larger scales and strengths as the turbulent updrafts. Although these marginal increases help to explain the preference for storms to initiate along CLs, they likely are not the whole story.