Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “horizontal”

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 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

“Gray Zone” Simulations Using a Three-Dimensional Planetary Boundary Layer Parameterization in the Weather Research and Forecasting Model

Abstract Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δ x ) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δ x ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation. Significance Statement The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.

Juliano, Timothy W.↗
DOE OSTI2013

Slip and Dilation Tendency Analysis of the Patua Geothermal Area

Critically stressed fault segments have a relatively high likelihood of acting as fluid flow conduits (Sibson, 1994). As such, the tendency of a fault segment to slip (slip tendency; Ts; Morris et al., 1996) or to dilate (dilation tendency; Td; Ferrill et al., 1999) provides an indication of which faults or fault segments within a geothermal system are critically stressed and therefore likely to transmit geothermal fluids. The slip tendency of a surface is defined by the ratio of shear stress to normal stress on that surface: Ts = T / on (Morris et al., 1996). Dilation tendency is defined by the stress acting normal to a given surface: Td = (o1-on) / (o1-o3) (Ferrill et al., 1999). Slip and dilation were calculated using 3DStress (Southwest Research Institute). Slip and dilation tendency are both unitless ratios of the resolved stresses applied to the fault plane by ambient stress conditions. Values range from a maximum of 1, a fault plane ideally oriented to slip or dilate under ambient stress conditions to zero, a fault plane with no potential to slip or dilate. Slip and dilation tendency values were calculated for each fault in the focus study areas at, McGinness Hills, Neal Hot Springs, Patua, Salt Wells, San Emidio, and Tuscarora on fault traces. As dip is not well constrained or unknown for many faults mapped in within these we made these calculations using the dip for each fault that would yield the maximum slip tendency or dilation tendency. As such, these results should be viewed as maximum tendency of each fault to slip or dilate. The resulting along-fault and fault-to-fault variation in slip or dilation potential is a proxy for along fault and fault-to-fault variation in fluid flow conduit potential. Stress Magnitudes and directions Stress field variation within each focus area was approximated based on regional published data and the world stress database (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2010; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012; Moeck et al., 2010; Moos and Ronne, 2010 and Reinecker et al., 2005) as well as local stress information if applicable. For faults within these focus systems we applied either a normal faulting stress regime where the vertical stress (sv) is larger than the maximum horizontal stress (shmax) which is larger than the minimum horizontal stress (sv>shmax>shmin) or strike-slip faulting stress regime where the maximum horizontal stress (shmax) is larger than the vertical stress (sv) which is larger than the minimum horizontal stress (shmax >sv>shmin) depending on the general tectonic province of the system. Based on visual inspection of the limited stress magnitude data in the Great Basin we used magnitudes such that shmin/shmax = .527 and shmin/sv= .46, which are consistent with complete and partial stress field determinations from Desert Peak, Coso, the Fallon area and Dixie valley (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2011; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012). Slip and dilation tendency analysis for the Patua geothermal system was calculated based on faults mapped in the Hazen Quadrangle (Faulds et al., 2011). Patua lies near the margin between the Basin and Range province, which is characterized by west-northwest directed extension and the Walker Lane province, characterized by west-northwest directed dextral shear. As such, the Patua area likely has been affected by tectonic stress associated with either or both of stress regimes over geologic time. In order to characterize this stress variation we calculated slip tendency at Patua for both normal faulting and strike slip faulting stress regimes. Based on examination of regional and local stress data (as explained above) we applied at shmin direction of 105 to Patua. Whether the vertical stress (sv) magnitude is larger than th...

15 GEOTHERMAL ENERGY↗
DOE OSTI2013

Slip and Dilation Tendency Analysis of Neal Hot Springs Geothermal Area

Slip and Dilation Tendency in focus areas Critically stressed fault segments have a relatively high likelihood of acting as fluid flow conduits (Sibson, 1994). As such, the tendency of a fault segment to slip (slip tendency; Ts; Morris et al., 1996) or to dilate (dilation tendency; Td; Ferrill et al., 1999) provides an indication of which faults or fault segments within a geothermal system are critically stressed and therefore likely to transmit geothermal fluids. The slip tendency of a surface is defined by the ratio of shear stress to normal stress on that surface: Ts = T / on (Morris et al., 1996). Dilation tendency is defined by the stress acting normal to a given surface: Td = (o1-on) / (o1-o3) (Ferrill et al., 1999). Slip and dilation were calculated using 3DStress (Southwest Research Institute). Slip and dilation tendency are both unitless ratios of the resolved stresses applied to the fault plane by ambient stress conditions. Values range from a maximum of 1, a fault plane ideally oriented to slip or dilate under ambient stress conditions to zero, a fault plane with no potential to slip or dilate. Slip and dilation tendency values were calculated for each fault in the focus study areas at, McGinness Hills, Neal Hot Springs, Patua, Salt Wells, San Emidio, and Tuscarora on fault traces. As dip is not well constrained or unknown for many faults mapped in within these we made these calculations using the dip for each fault that would yield the maximum slip tendency or dilation tendency. As such, these results should be viewed as maximum tendency of each fault to slip or dilate. The resulting along-fault and fault-to-fault variation in slip or dilation potential is a proxy for along fault and fault-to-fault variation in fluid flow conduit potential. Stress Magnitudes and directions Stress field variation within each focus area was approximated based on regional published data and the world stress database (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2010; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012; Moeck et al., 2010; Moos and Ronne, 2010 and Reinecker et al., 2005) as well as local stress information if applicable. For faults within these focus systems we applied either a normal faulting stress regime where the vertical stress (sv) is larger than the maximum horizontal stress (shmax) which is larger than the minimum horizontal stress (sv>shmax>shmin) or strike-slip faulting stress regime where the maximum horizontal stress (shmax) is larger than the vertical stress (sv) which is larger than the minimum horizontal stress (shmax >sv>shmin) depending on the general tectonic province of the system. Based on visual inspection of the limited stress magnitude data in the Great Basin we used magnitudes such that shmin/shmax = .527 and shmin/sv= .46, which are consistent with complete and partial stress field determinations from Desert Peak, Coso, the Fallon area and Dixie valley (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2011; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012). Based on inversion of fault kinematic data, Edwards (2013) interpreted that two discrete stress orientations are preserved at Neal Hot Springs. An older episode of east-west directed extension and a younger episode of southwest-northeast directed sinistral, oblique -normal extension. This interpretation is consistent with the evolution of Cenozoic tectonics in the region (Edwards, 2013). As such we applied a southwest-northeast (060) directed normal faulting stress regime, consistent with the younger extensional episode, to the Neal Hot Springs faults. Under these stress conditions northeast striking steeply dipping fault segments have the highest tendency to dilate and northeast striking 60 degrees dipping fault segments have the highest tendency to slip. Under these stress condi...

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Heavy Oil Polymer EOR in the Challenging Alaskan Arctic - It Works!

Under co-sponsorship of the US Department of Energy and Hilcorp Alaska LLC the first ever polymer field pilot commenced on 8/28/2018 in the Schrader Bluff heavy oil reservoir at the Milne Point Field on Alaska North Slope (ANS). The primary objective of the pilot is to prove the efficacy of polymer Enhanced Oil Recovery (EOR) to unlock the vast heavy oil resources on ANS. More than two and half years after startup, the polymer injection, supporting laboratory experiments and simulation studies steadily continue. The pilot started injecting hydrolyzed polyacrylamide (HPAM), at a concentration of 1,750 ppm to achieve a target viscosity of 45 cP, into the two horizontal injectors in the flood pattern. Production is monitored in the two horizontal producers. Based on laboratory measurements of polymer viscosity at reservoir conditions, the team decided to reduce polymer concentration to 1,200 ppm since July 2020 in an effort to control injection pressure and optimize polymer utilization. Quality control (QC) on the field ensures uniform polymer solution properties. Representative rock and fluid systems and test conditions are utilized in the corefloods on polymer retention, effect of injection water salinity, polymer loading, and their combinations on oil recovery. A history matched reservoir simulation model for forecasting oil recovery was developed on the basis of all the available field data. Field concerns related to the post-polymer breakthrough impact on flow assurance is addressed via specialized laboratory tests. Notwithstanding early operational disruptions and hydration issues, continuous polymer injection in both injectors has been achieved. To date, 950,000 lbs of polymer or 2 million barrels of polymer solution, equating to 13% of total pore volume (PV), has been placed in the flood pattern, serving as an effective indicator of adequate polymer injectivity. So far, the success of polymer EOR is evident from drastically reduced water cut in the producers, an estimated incremental 1,000 bopd over waterflood, and a favorable polymer utilization of 1.7 lbs/barrel of incremental oil. Polymer breakthrough was observed 26 months after the start of polymer injection. Main observations from corefloods are unusually high polymer retention values in some cases and a positive response to low salinity water. Although the heterogeneity in the flood pattern and exceptionally low water cut pose some challenges, persistent novel and justifiable simulation approaches have resulted in a robust history matched model. Experimental results on produced fluids treatment provide operational guidance for improved oil-water separation and mitigation of heater tube fouling. The ongoing success of the pilot is a key indicator of bringing the team ever closer to meeting the project’s primary objective. The encouraging results of the pilot is one of the drivers that has provided the impetus to apply polymer EOR throughout the Milne Point Field, which would increase oil recovery and extend the economic life of the Trans Alaska Pipeline System. The scientific knowledge, including the many lessons learned from this pilot also has referential value for other potential heavy oil EOR projects throughout the world. The collected data, operational lessons learned, and the overall success of the pilot are summarized in the paper.

None, None↗
DOE OSTI2021

Analysis of Completion Design Impact on Cluster Efficiency and Pressure-Based Well Communication in HFTS-2 Delaware Basin

The Hydraulic Fracturing Test Site 2 (HFTS-2) is a joint industry project in the Delaware basin to advance hydraulic fracturing understanding and improve productivity in shale reservoirs. The project integrates multi-disciplinary approaches to evaluate different completion designs, well spacing, inter-well communication, and stimulated rock volume, among other factors. This paper focused on two major areas related to hydraulic fracture performance. First, an analysis of different completion designs on cluster efficiency based on near well Distributed Acoustic Sensing (DAS). Second, an evaluation of well and completion designs on inter-well fracture driven interactions (FDIs) based on downhole pressure monitoring. Fluid/sand distribution and cluster efficiency analyses were based on near wellbore DAS data collected from two adjacent horizontal wells completed in two different landing zones in the Wolfcamp formation. These wells had different completion designs aiming to evaluate the effect of normal vs. extended stage lengths, perforation hole tapering and limited entry. Standard deviation from ideal fluid/sand distribution and waterfall plots were used to evaluate cluster efficiency for each design and stage. Inter-well FDIs analysis was conducted among the horizontal wells and a vertical monitor well. One horizontal well served as the monitor well while the other horizontal well was being treated. The vertical well was instrumented with downhole pressure and temperature gauges to aid monitoring fracture height growth. The pressure response during and after fracturing was characterized based on maximum pressure increase value and slope. Pressure response vs. FDIs trigger factors such distance, cluster efficiency and stage fluid volume were also analyzed. Based on the different completion and perforation designs tested, DAS analysis suggests that limited entry design worked best. Extreme limited entry showed the potential of high perforation erosion and reduced cluster efficiency. Here, the limited entry and tapered perforation design demonstrated potential to improve the cluster efficiency for extended stage lengths. Pressure monitoring across formation units proved to be critical to understand fracture interactions and fracture vertical growth. Pressure communication across different formation units during hydraulic fracturing operation indicate fractures grew upwards during Wolfcamp wells fracturing. However, this pressure communication dissipated over time. High intensity FDIs were recorded when the frac stages Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/1-21URTC/D011S018R001/2477607/urtec-2021-5289-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5289 were closer to the pressure gauge location in the monitor wells. Some of these stages that produced high intensity FDIs also had high fluid volume per cluster and low cluster efficiency. The multi-disciplinary and high-quality data collected from HFTS-2 helped to further understand why completion approaches such as limited entry and tapered perforation design are successful in improving cluster efficiency. The DAS data combined with downhole high-resolution pressure measurements also helped to quantify the effect of lower cluster efficiency data on the incidence and intensity of FDIs.

58 GEOSCIENCES↗
DOE OSTI2019

First Ever Polymer Flood Field Pilot - A Game Changer to Enhance the Recovery of Heavy Oils on Alaska’s North Slope

The development pace of Alaska’s vast, 20-25 billion barrels, heavy oil resources has been very slow due to high development costs and low oil recovery using conventional waterflood, and the impracticality of deploying thermal methods due to the presence of continuous permafrost. Although, polymer flooding has attracted attention and has become a promising EOR technique in heavy oil reservoirs due to the extensive application of horizontal wells and advancement of polymer flooding technology, no field tests have been performed to date in Alaska’s underdeveloped heavy oil reservoirs. The overall objective of this research is to perform a field experiment to validate the use of polymer flooding for extracting heavy oil in Alaska’s challenging environment. Two pre-existing pairs of horizontal injection and production wells in an isolated fault block of the Schrader Bluff heavy oil reservoir at the Milne Point Field are currently being used for the field experiment. Hydrolyzed polyacrylamide (HPAM) polymer injection started on August 28, 2018 at 600 ppm (4 cP viscosity) concentration ramping up to 1,800 ppm (45 cP viscosity) over a three week time period, and has been maintained at an average concentration of ~1,800 ppm. Current injection rates in the two horizontal injectors are ~2,200 and 600 bwpd. Laboratory experiments to determine the polymer retention, optimum water salinity, synergistic effects of water salinity and polymer, and handling of produced fluids, in support of the field experiment, are currently ongoing. Similarly, reservoir simulation of coreflood behavior and history match of previous waterfloods to predict polymer flood performance in the project area are also conducted in parallel. The field data and scientific knowledge that have been collected since the start of the injection indicates that the field pilot is performing as predicted. To date, no unexpected injectivity issues or polymer breakthrough have been encountered, and the two horizontal producers are showing positive response to the polymer injection, resulting in incremental increase in oil production rate. Since the research is still in its early stages, selected field, laboratory and simulation results are presented and discussed to highlight the integrative approach adopted in this first ever polymer flood field pilot in Alaska.

02 PETROLEUM↗
DOE OSTI2020

First Ever Polymer Flood Field Pilot to Enhance the Recovery of Heavy Oils on Alaska North Slope – Producer Responses and Operational Lessons Learned

The US Department of Energy and Hilcorp Alaska, LLC are jointly cosponsoring a field pilot project entitled “First Ever Field Pilot on Alaska’s North Slope to Validate the Use of Polymer Floods for Heavy Oil Enhanced Oil Recovery (EOR).” The objective of this pilot is to perform a field experiment to validate the use of an advanced polymer flooding technology to unlock the vast heavy oil resources on Alaska North Slope (ANS). Two pairs of horizontal injection and production wells have been deployed in an isolated fault block of the Schrader Bluff heavy oil reservoir at the Milne Point Field to conduct a polymer flood pilot which was designed to acquire field performance data to optimize polymer flood design in the Schrader Bluff heavy oil reservoirs on ANS. Hydrolyzed polyacrylamide (HPAM) polymer injection started on August 28, 2018 using a custom made polymer mixing and pumping unit. Polymer solution properties such as concentration, viscosity, and filter ratio, have been closely monitored and adjusted as necessary. Operational procedures have been improved to ensure adequate injectivity and polymer propagation through the formation. Oil and water production is tested frequently to monitor water cut performance, and produced water is sampled and analyzed for polymer content. Tracers were pumped pre and post polymer injection startup to compare polymer flood breakthrough timing and sweep efficiency. Polymer injectivity of the horizontal wells drilled in the Schrader Bluff NB sand is found to be sufficient to replace the reservoir production voidage although there were some declines as high viscosity polymer swept the near wellbore region. Production data show significant reduction in water cut and the oil production rate has stabilized. No polymer production has been confirmed from the two horizontal producers after 22 months of polymer injection into the two supporting horizontal injectors. Tracer data show very low sweep efficiencies with waterflood alone, whereas significant improvement has been achieved with polymer flood. Many operational lessons have been learned and improvements have been made regarding polymer mixing, pumping, and onsite quality control, particularly in a challenging Arctic environment. Adequate polymer injection rate can be achieved with diligent monitoring and quality assurance of the properties of the injected polymer solution.

Ning, Samson↗
DOE OSTI2016

Data Report for Monitoring at Six West Virginia Marcellus Shale Development Sites using NETL’s Mobile Air Monitoring Laboratory (July–November 2012)

The West Virginia Department of Environmental Protection’s Office of Oil and Gas was directed according to the Natural Gas Horizontal Well Control Act of December 14, 2011 (West Virginia Code §22-6A) to conduct studies of horizontal well drilling activities related to air quality. The planned study, “Noise, Light, Dust, Volatile Organic Compounds Related to Well Location Restrictions,” required determination of the effectiveness of a 625 ft minimum set-back from the center of the pad of a horizontal well drilling site to the nearest occupied dwelling. An investigation was conducted at seven drilling sites by West Virginia University (WVU) and the National Energy Technology Laboratory (NETL) to collect data on dust, hydrocarbon compounds and on noise, radiation, and light levels. NETL’s role in this study was to collect measurements of ambient pollutant concentrations at six of the seven selected sites using NETL’s Mobile Air Monitoring Laboratory. The trailer-based laboratory was situated a distance of 492–1,312 ft from each well pad, on which activities included well pad construction, vertical drilling, horizontal drilling, hydraulic fracturing, and flaring, with the objective of evaluating the air quality impact of each activity for 1–4 weeks per site. Measured pollutants included volatile organic compounds (VOCs), coarse and fine particulate matter (PM 10 and PM 2.5 , respectively), ozone, methane (CH 4 ), carbon dioxide (CO 2 ), carbon isotopes of CH 4 and CO 2 , organic carbon (OC), elemental carbon (EC), oxides of nitrogen (NOx), and sulfur dioxide (SO 2 ).

03 NATURAL GAS↗
DOE OSTI2016

Evaluation of Erosion of the Dummy “EE” Plate 19 in YA Type ATR Fuel Element During Reactor PALM Cycles

Advanced Test Reactor (ATR) Cycle 153B-1 was a 14-day, high-power, powered axial locator mechanism (PALM) operating cycle that completed on April 12, 2013. Cycle 153B-1 was a typical operating cycle for the ATR, and did not result in any unusual plant transients. ATR was started up and shut down as scheduled. The PALM drive physically moves the selected experiments into and out of the core to simulate reactor startup and heat up, and shutdown and cooldown transients, while the reactor remains in steady-state conditions. However, after the cycle was over, when the fuel elements were removed from the core and inspected, several thousand flow-assisted erosion pits and “horseshoeing” defects were readily observed on the surface of the several YA-type fuel elements (these are aluminum “dummy” plates that contain no fuel). In order to understand these erosion phenomena, a thermal-hydraulic model of coolant channel 20 on a YA-M fuel element was generated. The boundaries of the model were the aluminum EE plate of a YA-M fuel element and a beryllium reflector block with 13 horizontal saw cuts which represented regions of zero flow. The heat generated in fuel plates 1 through 18 was modeled to be passing through the aluminum EE plate. The coolant channel 20 width was set at 0.058 in. (58 mils). It was established that the horizontal saw cuts had a significant effect on the temperature of the coolant. The flow, which was expected to vary linearly with gradual heating of the coolant as it passed through the channel, was extremely turbulent. The temperature rise, which was expected to be a smooth “S” curve, was represented by a series temperature rise “humps,” which occurred at each horizontal saw cut in the beryllium reflector block. Each of the 13 saw cuts had a chamfered edge which resulted in the coolant flow being re-directed as a jet across the coolant channel into the surface of the EE plate, which explained the temperature rise and the observed scalloping and pitting degradation on the YA-M fuel elements. In the case of scalloping (horseshoeing) a surprising similarity of that defect to those appearing on aluminum plate rolled in over-lubrication conditions, were established. In turn, this made us think that the principal feature responsible for the appearance of these defects, was horizontal cuts in the beryllium reflector block created to arrest the propagation of large vertical crack(s) in Be in PALM cycles with higher overall fluence. This assumption was fully confirmed by the results of thermo-hydraulic simulations. The neutronics data for these modeling experiments were provided using advanced irradiation simulations (MCNP, HELIOS). In the case of pitting erosion the following corrective measures were proposed based upon the results of JMatPro v.8.2 modeling (TTT- and CCT-diagrams): change the fabrication process by adding blister anneal before program anneal, immediately after cold rolling of AA6061plate. This step will allow achieving complete recrystallization, eliminating of strengthening due to metastable precipitates, and reduce the possibility of forming sharp microstructural features upon the surface.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Retrieval of Boundary Layer 3D Cloud Properties Using Scanning Cloud Radar and 3D Radiative Transfer

Retrievals of cloud optical and microphysical properties for boundary layer clouds, including those widely used by ASR investigators, frequently assume that clouds are sufficiently horizontally homogeneous that scattering and absorption (at all wavelengths) can be treated using one dimensional (1D) radiative transfer, and that differences in the field-of-view of different sensors are unimportant. Unfortunately, most boundary layer clouds are far from horizontally homogeneous, and numerous theoretical and observational studies show that the assumption of horizontal homogeneity leads to significant errors. The introduction of scanning cloud and precipitation radars at the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program sites presents opportunities to move beyond the horizontally homogeneous assumption. The primary objective of this project was to develop a 3D retrieval for warm-phase (liquid only) boundary layer cloud microphysical properties, and to assess errors in current 1D (non-scanning) approaches. Specific research activities also involved examination of the diurnal cycle of hydrometeors as viewed by ARM cloud radar, and continued assessment of precipitation impacts on retrievals of cloud liquid water path using passive microwaves.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗
DOE OSTI2021

Modeling Environmental Effects on Ventilated Spent Fuel Storage Systems

This report describes newly developed external environment wind effects models of spent nuclear fuel (SNF) dry storage systems. The primary purpose of these wind effects models is to better understand particle deposition on SNF canisters in the context of chloride-induced stress corrosion cracking. The goal of this effort is to further the understanding and improve the wind effects model of the Nuclear Horizontal Modular Storage (NUHOMS ® ) Advanced Horizontal Storage Module–High Seismic (AHSM-HS) storage system in support of the planned Canister Deposition Field Demonstration (CDFD) project (Durbin et al. 2021, Suffield et al. 2021). The steps in this wind effects investigation begin with a validation case of comparing experimental data with a STAR-CCM+ computational fluid dynamics (CFD) model of the Belowground Vertical Dry Cask Simulator (BVDCS) with the external environment explicitly modeled. Next, a test case is performed with the AHSM-HS models comparing solar-loading strategies for a standalone and wind effects model in STAR-CCM+ and a standalone model in ANSYS Mechanical Parametric Design Language (APDL). The final portion of this study compares results of the validation exercise with applications of wind effects models for two horizontal storage systems, a site specific NUHOMS ® horizontal storage module and a NUHOMS ® AHSM-HS.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Experimental Measurements of Soot and Soot Precursors of Sustainable Aviation Fuels

Sustainable Aviation Fuels (SAFs) are promising for decarbonizing the challenging-to-electrify aerospace sector and reducing yearly greenhouse gas emissions by thousands of tons. Developing predictive detailed chemical kinetic models for SAF (and traditional) fuel mixtures necessitates a robust experimental database to elucidate the intricate pyrolysis and oxidation chemistry of SAF fuel blend components. Two critical components of SAF surrogate mixtures are iso-dodecane isomers and trans-decalin. This report characterizes the structure of three slightly sooting, laminar, non-premixed, nitrogen-diluted Planar Mixing Layer Flames (PMLFs) fueled by either pure ethylene, ethylene doped with 1500ppm of either 2,2,4,6,6-pentamethyl-heptane or trans-decalin. The flames have the same stoichiometric mixture fraction and total mole fraction of fuels in the fuel stream (X F,F = X C2H4,F + X Dopant,F = 0.26), resulting in approximately the same maximum temperature (T max ≈ 1800 K) to highlight the effect of doping on soot and its gaseous precursors. The PMLF is unconfined and established between adjacent planar jets of the fuel and oxidizer streams surrounded by annular shielding nitrogen. The flow is stabilized by exhausting the hot buoyant flame products through a slot in a downstream plate onto which the cold flow impinges. Importantly, any horizontal cross-section of the PMLF has a self-similar thermochemical structure of a reactive boundary layer that grows at increasing Height Above the Burner (HAB). Therefore, the horizontal cross-sections of a PMLF can be modeled as a One-Dimensional Counterflow Flame (1D-CF) with very low strain rates which are unachievable in traditional CFs and yield a structure that is several millimeters thick. In brief, the PMLF is ideal for performing spatially resolved measurements with minimal effects of the sampling-induced perturbations, to develop and validate detailed chemical kinetic models on time scales of a few tens of milliseconds that are relevant in aviation applications. The horizontal profiles of C0-C18 gaseous species at HAB = 50 mm and soot volume fraction at 25 and 50 mm are measured using capillary sampling followed by GC-MS analyses and Laser Induced Emission Spectroscopy (LIES), respectively. Additionally, the Elastic Laser Light Scattering (E-LLS) coefficient is measured at HAB = 50 mm and 25 mm to determine the strain rates of the 1D-CFs equivalent to the characterized PMLF horizontal cross sections and the profile of the E-LLS equivalent diameter of soot. The comparison of the results in the base PMLF fueled only by ethylene with those in either doped PMLF in which 1500 ppm of ethylene are replaced with 2,2,4,6,6-pentamethyl-heptane or trans-decalin, reveals that the doping increases the mole fraction of several Polycyclic Aromatic Hydrocarbons (PAHs) at both HABs by a factor of approximately 2.5, and the soot volume fraction at both HABs by factors of ~1.5 for 2,2,4,6,6-pentamethyl-heptane and ~1.25 for trans-decalin doping, respectively. Conversely, doping causes a near doubling of the soot equivalent diameter only at HAB = 50 mm, without significant effects at 25 mm. The experimental results partially validate the tested state-of-the-art detailed chemical kinetic model and also point toward the further improvement of its predictive capabilities.

10 SYNTHETIC FUELS↗
DOE OSTI2024

Hydrogen, Methane, Brine Flow Behavior, and Saturation in Sandstone Cores During H 2 and CH 4 Injection and Displacement

Large-scale underground hydrogen storage (UHS) is a critical component in the emerging hydrogen economy. Knowledge of multiphase flow behavior involving hydrogen in storage reservoir formations is crucial to characterizing hydrogen transport properties and essential for the deliverability and storage operations of UHS. There are still many gaps in fully understanding hydrogen–methane–brine multiphase phase flow that require further investigation. In this work, H 2 and CH 4 were injected through brine-saturated sandstone cores using a tri-axial core holder system fitted with flow rate meters and pressure transducers, while the effluent gas concentrations were analyzed using an online micro gas chromatograph. Brine displacement, permeability, and gas breakthrough curves were measured. We studied the flow behavior of hydrogen and methane in sandstone cores through testing brine displacement by gas injection and comparing the hydrogen displacement of methane with the methane displacement of hydrogen. We also tested the differences between horizontal and vertical flow in brine displacement. The results showed that brine displacement was more efficient in a core with higher permeability and porosity, resulting in a higher initial gas saturation. A higher gas injection rate brought about faster gas breakthrough measured by pore volume and sharper concentration curves. Hydrogen did not exhibit abnormal flow in the sandstone when the flow was horizontal and downward vertical. Gas overriding was observed in brine displacements when the flow was horizontal, with hydrogen showing this behavior more profoundly compared to methane. Downward vertical gas injection induced higher efficiency brine displacement compared to horizontal displacement and resulted in a higher initial gas saturation in the sandstone cores. These findings address critical knowledge gaps regarding gas flow patterns and displacement behaviors during hydrogen injection and recovery phases in UHS facilities using methane as the cushion gas. The insights from this research offer valuable guidance for optimizing UHS systems, ensuring operational efficiency, and advancing sustainable energy solutions in alignment with decarbonization goals.

08 HYDROGEN↗
DOE OSTI2014

Technical Note: On the use of nudging for aerosol–climate model intercomparison studies

Nudging as an assimilation technique has seen increased use in recent years in the development and evaluation of climate models. Constraining the simulated wind and temperature fields using global weather reanalysis facilitates more straightforward comparison between simulation and observation, and reduces uncertainties associated with natural variabilities of the large-scale circulation. On the other hand, the forcing introduced by nudging can be strong enough to change the basic characteristics of the model climate. In the paper we show that for the Community Atmosphere Model version 5 (CAM5), due to the systematic temperature bias in the standard model and the sensitivity of simulated ice formation to anthropogenic aerosol concentration, nudging towards reanalysis results in substantial reductions in the ice cloud amount and the impact of anthropogenic aerosols on long-wave cloud forcing. In order to reduce discrepancies between the nudged and unconstrained simulations, and meanwhile take the advantages of nudging, two alternative experimentation methods are evaluated. The first one constrains only the horizontal winds. The second method nudges both winds and temperature, but replaces the long-term climatology of the reanalysis by that of the model. Results show that both methods lead to substantially improved agreement with the free-running model in terms of the top-of-atmosphere radiation budget and cloud ice amount. The wind-only nudging is more convenient to apply, and provides higher correlations of the wind fields, geopotential height and specific humidity between simulation and reanalysis. Results from both CAM5 and a second aerosol–climate model ECHAM6-HAM2 also indicate that compared to the wind-and-temperature nudging, constraining only winds leads to better agreement with the free-running model in terms of the estimated shortwave cloud forcing and the simulated convective activities. This suggests nudging the horizontal winds but not temperature is a good strategy for the investigation of aerosol indirect effects since it provides well-constrained meteorology without strongly perturbing the model's mean climate.

Zhang, K.↗
DOE OSTI2022

Further improvement and evaluation of nudging in the E3SM Atmosphere Model version 1 (EAMv1): simulations of the mean climate, weather events, and anthropogenic aerosol effects

Abstract. A previous study on the use of nudging in E3SM Atmosphere Model version 1 (EAMv1) had an unresolved issue; i.e., a simulation nudged to EAMv1's own meteorology showed non-negligible deviations from the free-running baseline simulation over some of the subtropical marine stratocumulus and trade cumulus regions. Here, we demonstrate that the deviations can be substantially reduced by (1) changing where the nudging tendency is calculated in the time integration loop of a nudged EAM simulation so as to improve consistency with the free-running baseline and (2) increasing the frequency of the constraining data so as to better capture strong sub-diurnal variations. The fact that modification (2) improves the climate representativeness of the nudged simulations has motivated us to investigate whether the use of newer reanalysis products with higher data frequency can help improve nudged hindcast simulations by better capturing the observed weather events. To answer this question, we present simulations conducted at EAMv1's standard horizontal resolution (approximately 1∘) with nudging towards 6-hourly ERA-Interim reanalysis or 6-hourly, 3-hourly, or hourly ERA5 reanalysis. These simulations are evaluated against the climatology of free-running EAMv1 simulations as well as reanalyses, satellite retrievals, and in situ measurements from the Atmospheric Radiation Measurement user facility. For the 1∘ EAMv1 simulations, we recommend using the relocated nudging tendency calculation and the ERA5 reanalysis at 3-hourly or higher frequency. Simulations used for estimating the anthropogenic aerosol effects often use nudging to help discern signal from noise. The sensitivity of such estimates to the configuration of nudging is investigated in EAMv1, again using the standard 1∘ horizontal resolution. We find that, when estimating the global mean effects, the frequency of constraining data has relatively small impacts, while the choice of nudged variables can change the results substantially. The nudging of air temperature (in addition to horizontal winds) has two non-negligible effects. First, when the constraining data come from reanalysis, the nudging-induced mean bias correction can cause significant changes in the simulated clouds and hence substantially different estimates of the aerosol effects. The impact of the mean bias correction on ice cloud formation has been noted in previous studies and is also seen in EAMv1. For applications like ours, where the preferred configurations of nudging are those capable of providing results consistent with the multi-year free-running simulations, the consequence of the mean bias correction is undesirable. The second important impact of temperature nudging is a significant suppression of adjustments to aerosol forcing, which also causes changes in the estimated aerosol effects. This effect can be seen in simulations nudged to either reanalysis or EAM's own meteorology. These results suggest that nudging horizontal winds but not temperature is a better choice for estimating the anthropogenic aerosol effects.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1973

Transient boiling heat transfer in saturated liquid nitrogen and F113 at standard and zero gravity

Transient and steady state nucleate boiling in saturated LN2 and F113 at standard and near zero gravity conditions were investigated for the horizontal up, vertical and horizontal down orientations of the heating surface. Two distinct regimes of heat transfer mechanisms were observed during the interval from the step increase of power input to the onset of nucleate boiling: the conduction and convection dominated regimes. The time duration in each regime was considerably shorter with LN2 than with F113, and decreased as heat flux increased, as gravity was reduced, and as the orientation was changed from horizontal up to horizontal down. In transient boiling, boiling initiates at a single point following the step increase in power, and then spreads over the surface. The delay time for the inception of boiling at the first site, and the velocity of spread of boiling varies depending upon the heat flux, orientation, body force, surface roughness and liquid properties, and are a consequence of changes in boundary layer temperature levels associated with changes in natural convection. Following the step increase in power input, surface temperature overshoot and undershoot occur before the steady state boiling temperature level is established.

Oker, E.↗
NASA NTRS1975

Lightning properties in Florida thunderstorms from video tape records

From video-tape recordings of lightning flashes to ground in convective daytime storms at the NASA Kennedy Space Center, Florida, and near Ocala, Florida, statistical data are presented for the flash-time duration, the number of strokes per flash, the time between all strokes, and the time between strokes having spatially separate channels. A daytime video tape 'photograph' of a stepped leader is reproduced. Finally, statistical data are presented for 13 long horizontal lightning discharges occurring along a line of nocturnal thunderstorms near Gainesville, Florida. These discharges all propagated horizontally near the freezing level (4.5 km) at an apparent average velocity between 5600 and 11,000 m/s for a typical distance of 2 to 3 km. Eleven of the 13 horizontal discharges were preceded by or occurred within 16.7 ms of a glow in the cloud roughly 0.5 km above the horizontal channel.

Brantley, R. D.↗