Search NASA⌕ Search

SEARCH · Search NASA

Results for “Relative permeability”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Geomechanical Framework for Secure CO 2 Storage in Fractured Reservoirs and Caprocks for Sedimentary Basins in theMidwest United States

This report presents final technical results for the project Geomechanical Framework for Secure CO 2 Storage in Fractured Reservoirs and Caprocks for Sedimentary Basins in the Midwest United States (DE-FE0023330). The project was a three-year effort consisting of seven technical tasks focused on defining geomechanical factors for CO 2 storage applications in deep saline rock formations in Ohio and the Midwest United States, because geomechancial issues have been identified as a significant risk factor for large-scale CO 2 storage applications. A basin-scale stress-strain analysis was completed to describe the geomechanical setting for rock formations of Ordovician-Cambrian age in Ohio and adjacent areas of the Midwest United States in relation to geologic CO 2 storage applications. The tectonic setting, stress orientation-magnitude, and geomechanical and petrophysical parameters for CO 2 storage zones and caprocks in the region were cataloged. Ten geophysical image logs were analyzed for natural fractures, borehole breakouts, and drilling-induced fractures. The logs indicated mostly less than 10 fractures per 100 vertical feet in the borehole, with mostly N65E principal stress orientation through the section. Geophysical image logs and other logs were obtained for three wells located near the sites where specific models were developed for geomechanical simulations: Arches site in Boone County, Kentucky; Northern Appalachian Basin site in Chautauqua County, New York; and E-Central Appalachian Basin site in Tuscarawas County, Ohio. For these three wells, 9,700 feet of image logs were processed and interpreted to provide a systematic review of the distribution within each well of natural fractures, wellbore breakouts, faults, and drilling induced fractures. There were many borehole breakouts and drilling-induced tensile fractures but few natural fractures. Concentrated fractures were present at the Rome-basal sandstone and basal sandstone-Precambrian contacts at the Arches and East-Central Appalachian Basin sites. Geophysical logs were utilized to develop local-scale geologic models by determining geomechanical and petrophysical parameters within the geologic formations. These data were ported to coupled fluid-flow and reservoir geomechanics multi-phase CO 2 injection simulations. The models were developed to emphasize the geomechanical layers within the CO 2 storage zones and caprocks. A series of simulations were completed for each site to evaluate whether commercial-scale CO 2 could be safely injected into each site, given site-specific geologic and geomechanical controls. This involved analyzing the simulation results for the integrity of the caprock, intermediate, and reservoir zones, as well quantifying the areal uplift at the surface. Simulation results were also examined to ensure that the stress-stress perturbations were isolated within the subsurface, and that there was only limited upward migration of the CO 2 . Simulations showed capacity to inject more than 10 million metric tons of CO 2 in a single well at the Arches and East Central Appalachian Basin sites without excessive geomechanical risks. Low-permeability rock layers at the Northern Appalachian Basin study area well resulted in very low CO 2 injection capacity. Fracture models developed for the sites suggests that the sites have sparse fracture network in the deeper Cambrian rocks. However, there were indicators in image logs of a moderate fracture matrix in the Rose Run Sandstone at the Northern Appalachian Basin site. Dual permeability fracture matrix simulations suggest the much higher injection rates may be feasible in the fractured interval. Guidance was developed for geomechanical site characterization in the areas of geophysical logging, rock core testing, well testing, and site monitoring. The guidance demonstrates that there is a suitable array of options for addressing geomechanical issues at CO 2 storage sites. Finally, a review of Marcellus and Utica-Point Pleasant shale gas wells and CO 2 storage intervals indicates that these items are vertically separated, except for the Oriskany sandstone and Marcellus wells in southwest Pennsylvania and northern West Virginia. Together, project results present a more realistic portrayal of geomechanical risk factors related to CO 2 storage for existing and future coal-fired power plants in Ohio.

58 GEOSCIENCES↗
DOE OSTI2016

Play fairway analysis of geothermal resources across the state of Hawaii: 2. Resource probability mapping

We develop a new geostatistical method to combine evidence provided by diverse geological data sets and produce maps of geothermal resource probability. The application is to the State of Hawaii, and the data sets include the locations and ages of mapped volcanic centers, gravity and magnetotelluric measurements, groundwater temperature and geochemistry, ground surface deformation, seismicity, water table elevation, and groundwater recharge. Using the basic principles of Bayesian statistics, these data and expert knowledge about the effects and importance of the data are used to compute the probabilities of the primary resource qualities of elevated subsurface heat, reservoir permeability, and reservoir fluid content. The product of these marginal probabilities estimates the joint probability of all three qualities and hence the probability of a successful geothermal prospect at each map point. An analogous set of algorithms is used to quantify the confidence in the probability at each point. Not surprisingly, we find that successful geothermal prospects are most probable on the active volcanoes of Hawaii Island, including the area of Hawaii’s single geothermal energy plant. Probability decreases primarily with shield volcano age, being relatively moderate in select locations on Maui and Lanai, relatively low on Oahu, and minimal on Kauai. Future exploration efforts should consider these results as well as the practical, societal, and economic conditions that influence development viability. Here, the difficulties of interisland power transmission mean that even areas with moderate to low probabilities are worth investigating on islands with population centers.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

Propagation, proppant transport and the evolution of transport properties of hydraulic fractures

Hydraulic fracturing is a widely used method for well stimulation to enhance hydrocarbon recovery. Permeability, or fluid conductivity, of the hydraulic fracture is a key parameter to determine the fluid production rate, and is principally conditioned by fracture geometry and the distribution of the encased proppant. A numerical model is developed to describe proppant transport within a propagating blade-shaped fracture towards defining the fracture conductivity and reservoir production after fracture closure. Fracture propagation is formulated based on the PKN-formalism coupled with advective transport of an equivalent slurry representing a proppant-laden fluid. Empirical constitutive relations are incorporated to define rheology of the slurry, proppant transport with bulk slurry flow, proppant gravitational settling, and finally the transition from Poiseuille (fracture) flow to Darcy (proppant pack) flow. At the maximum extent of the fluid-driven fracture, as driving pressure is released, a fracture closure model is employed to follow the evolution of fracture conductivity with the decreasing fluid pressure. This model is capable of accommodating the mechanical response of the proppant pack, fracture closure of potentially contacting rough surfaces, proppant embedment into fracture walls, and most importantly flexural displacement of the unsupported spans of the fracture. Results show that reduced fluid viscosity increases the length of the resulting fracture, while rapid leak-off decreases it, with both characteristics minimizing fracture width over converse conditions. Proppant density and size do not significantly influence fracture propagation. Proppant settling ensues throughout fracture advance, and is accelerated by a lower viscosity fluid or greater proppant density or size, resulting in accumulation of a proppant bed at the fracture base. ‘Screen-out’ of proppant at the fracture tip can occur where the fracture aperture is only several times the diameter of the individual proppant particles. After fracture closure, proppant packs comprising larger particles exhibit higher conductivity. More importantly, high-conductivity flow channels are necessarily formed around proppant banks due to the flexural displacement of the fracture walls, which offer preferential flow pathways and significantly influence the distribution of fluid transport. Higher compacting stresses are observed around the edge of proppant banks, resulting in greater depths of proppant embedment into the fracture walls and/or an increased potential for proppant crushing.

Mechanics↗
DOE OSTI2016

Deep Borehole Field Test Research Activities at LBNL

The goal of the U.S. Department of Energy Used Fuel Disposition’s (UFD) Deep Borehole Field Test is to drill two 5 km large-diameter boreholes: a characterization borehole with a bottom-hole diameter of 8.5 inches and a field test borehole with a bottom-hole diameter of 17 inches. These boreholes will be used to demonstrate the ability to drill such holes in crystalline rocks, effectively characterize the bedrock repository system using geophysical, geochemical, and hydrological techniques, and emplace and retrieve test waste packages. These studies will be used to test the deep borehole disposal concept, which requires a hydrologically isolated environment characterized by low permeability, stable fluid density, reducing fluid chemistry conditions, and an effective borehole seal. During FY16, Lawrence Berkeley National Laboratory scientists conducted a number of research studies to support the UFD Deep Borehole Field Test effort. This work included providing supporting data for the Los Alamos National Laboratory geologic framework model for the proposed deep borehole site, conducting an analog study using an extensive suite of geoscience data and samples from a deep (2.5 km) research borehole in Sweden, conducting laboratory experiments and coupled process modeling related to borehole seals, and developing a suite of potential techniques that could be applied to the characterization and monitoring of the deep borehole environment. The results of these studies are presented in this report.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Proton Conducting Silicon Oxide Membranes as a Fluorine Free Alternative to Nafion for Low Temperature Water Electrolysis

Driven by environmental and health concerns related to per- and polyfluoroalkyl substances (PFAS), there has been growing interest in developing fluorine-free proton (H + ) exchange membrane (PEM) materials for fuel cells and water electrolyzers. In this study, we present a side-by-side comparison of the key transport properties of submicron thick, PFAS-free amorphous silicon dioxide (SiO 2 ) membranes to Nafion, a fluorinated polymer electrolyte membrane that represents the industry standard for PEM fuel cells and electrolyzers. Here, measurements of proton (H + ) conductivity (σ H + ), hydrogen (H 2 ) permeability (P H 2 ), and electrical resistivity (ρ e – ) were conducted using model thin films comprised of SiO 2 membranes deposited by atomic layer deposition (ALD). Although the H + conductivity of the SiO 2 membranes is 2–3 orders of magnitude lower than Nafion, the addition of phosphorus dopants (PO x ) improves H + conductivity such that the area specific membrane resistance of thin (<50 nm) PO x -doped SiO 2 membranes is more than an order of magnitude lower than Nafion-117. Importantly, the safe operation of such nanoscale membranes within a PEM electrolyzer is feasible thanks to the low H 2 permeability of dense SiO 2 -based membranes, which are predicted to limit H 2 crossover rates to acceptable levels for pressures up to ≈ 100 bar. As a proof-of-principle demonstration, a chip-scale water electrolyzer based on 100 nm thick PO x -SiO 2 membrane is shown to achieve a current density of 2 A cm –2 at a potential of 2.5 V. If this technology can be successfully scaled up, H + conducting oxide membranes offer an attractive PFAS-free alternative to Nafion for efficient and durable water electrolysis and fuel cell technologies.

Nafion↗
DOE OSTI2016

Field Testing Unvented Roofs with Asphalt Shingles in Cold and Hot-humid Climates

BSC TO5 Task 7.1 Field Testing Unvented Roofs with Asphalt Shingles in Cold Climates - Bolingbrook, IL 60440 BSC TO5 Task 7.1 Field Testing Unvented Roofs with Asphalt Shingles in Hot-Humid Climates - Friendswood, TX 77546 To understand the long-term moisture performance of unvented roof assemblies with fibrous insulation, the U.S. Department of Energy's Building America team Building Science Corporation built and monitored test roofs in two climates: a cold-climate test house in Bolingbrook, Illinois (Chicago area; zone 5A), and a hot-humid test house in Friendswood, Texas (Houston area, zone 2A). The Chicago-area test bed had seven parallel experimental rafter bays. The assemblies included a control vented compact (cathedral) roof, a dense-pack cellulose unvented roof, and an unvented roof with a "diffusion vent" (a strip of vapor-permeable gypsum sheathing at the ridge to allow drying). The interior finish was gypsum wallboard (GWB) with latex paint. The other four bays were top vent roof assemblies, which have a polypropylene breather mesh between the roof sheathing and the asphalt shingles to allow ventilation drying of the assembly from underneath the vapor-impermeable shingles. Two of the four "top vent" roofs were fiberglass batt and two were dense-pack cellulose roofs with either interior GWB or open to the interior. Wintertime moisture-related failure was accelerated by providing interior humidification (22.2 degree C [72 degree F]/50% relative humidity [RH]), which is known to be an extreme interior loading). Data were captured over 8 months, including a winter and the following spring and early summer.

Array↗
DOE OSTI2018

Characterization and property investigation of novel eco‐friendly agar/carrageenan/TiO 2 nanocomposite films

ABSTRACT Agar/kappa carrageenan (AG/CA) nanocomposite films are prepared using thermal casting method in the presence of various concentrations of nano‐TiO 2 (0.5, 1, 1.5, and 2 g/100 g). As nano‐TiO 2 content rises from 0 to 1 g/100 g, tensile strength (TS) increases from 30.8 to 60.6 MPa. However, TS decreases gradually as nano‐TiO 2 content further raises. TS of films containing 1 g/100 g nano‐TiO 2 increases by 96.8% when compared to those prepared without nano‐TiO 2 . The reason may be related to the formation of intermolecular hydrogen bonding and coordinate bonds between nano‐TiO 2 and matrix molecules. Elongation at break decreases gradually as nano‐TiO 2 concentration increases. The addition of nano‐TiO 2 improves barrier properties of films against ultraviolet light, while it decreases their transparency. Swelling and water solubility of the films containing 1 g/100 g nano‐TiO 2 decrease, further demonstrating that interaction forces are formed between nano‐TiO 2 and matrix molecules. The addition of nano‐TiO 2 significantly declines water vapor permeability. Overall, films containing 1 g/100 g nano‐TiO 2 show the best mechanical, optical, and water resistance properties. © 2018 Wiley Periodicals, Inc. J. Appl. Polym. Sci. 2019 , 136 , 47113.

Hou, X.↗
DOE OSTI2023

Utah FORGE Core Material Friction Experiments

The following folders contain 3 friction experiments on Utah FORGE core material from depths of 6735 - 8533ft. The experiment conditions range from 100% relative humidity at room temperature to fluid-saturated at 110C. Experiments were done on both intact samples in a single-direct shear configuration, double-direct shear, and gouge.

15 GEOTHERMAL ENERGY↗
NASA NTRS1974

Progress of the LASL dry hot rock geothermal energy project

The possibilities and problems of extracting energy from geothermal reservoirs which do not spontaneously yield useful amounts of steam or hot water are discussed. The system for accomplishing this which is being developed first is a pressurized-water circulation loop intended for use in relatively impermeable hot rock. It will consist of two holes connected through the hot rock by a very large hydraulic fracture and connected at the surface through the primary heat exchanger of an energy utilization system. Preliminary experiments in a hole 2576 ft (0.7852 km) deep, extending about 470 ft (143 m) into the Precambrian basement rock underlying the Jemez Plateau of north-central New Mexico, revealed no unexpected difficulties in drilling or hydraulically fracturing such rock at a temperature of approximately 100 C, and demonstrated a permeability low enough so that it appeared probable that pressurized water could be contained by the basement rock. Similar experiments are in progress in a second hole, now 6701 ft (2.043 km) deep, about 1.5 miles (2.4 km) south of the first one.

Smith, M. C.↗
DOE OSTI2021

Probing the Damage Zone at Parkfield

Rocks are heterogeneous materials that exhibit nonlinear elastic (anelastic) behavior at scales ranging from the laboratory to Earth. In the laboratory, typical, complex relationships exist between stress and strain that include hysteresis, finite relaxation times, strain rate, and history dependence. These behaviors are linked to important characteristics such as stress, porosity, permeability, material integrity, and material failure. In this work, we adopted a “pump-probe” type experiment common in laboratory studies, using solid earth tides as the low-frequency pump and empirical Green's function as the high-frequency probe. By probing the velocity at different points in the pump cycle, we constrained important information about the strain-modulus relationship. Near the San Andreas Fault, we observed strongly nonlinear elastic behavior that characterizes the damage zone. We also constrained important aspects of hysteretic behavior that are related to damage properties and possibly pore pressure. Away from the fault, the nonlinear behavior is diminished.

58 GEOSCIENCES↗
DOE OSTI2017

Conformal Lithium Fluoride Protection Layer on Three-Dimensional Lithium by Nonhazardous Gaseous Reagent Freon

Research on lithium (Li) metal chemistry has been rapidly gaining momentum nowadays not only because of the appealing high theoretical capacity, but also its indispensable role in the next-generation Li–S and Li–air batteries. However, two root problems of Li metal, namely high reactivity and infinite relative volume change during cycling, bring about numerous other challenges that impede its practical applications. In the past, extensive studies have targeted these two root causes by either improving interfacial stability or constructing a stable host. However, efficient surface passivation on three-dimensional (3D) Li is still absent. Here in this paper, we develop a conformal LiF coating technique on Li surface with commercial Freon R134a as the reagent. In contrast to solid/liquid reagents, gaseous Freon exhibits not only nontoxicity and well-controlled reactivity, but also much better permeability that enables a uniform LiF coating even on 3D Li. By applying a LiF coating onto 3D layered Li-reduced graphene oxide (Li-rGO) electrodes, highly reduced side reactions and enhanced cycling stability without overpotential augment for over 200 cycles were proven in symmetric cells. Furthermore, Li–S cells with LiF protected Li-rGO exhibit significantly improved cyclability and Coulombic efficiency, while excellent rate capability (~800 mAh g -1 at 2 C) can still be retained.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

New insights into the drainage of inundated ice-wedge polygons using fundamental hydrologic principles

Abstract. The pathways and timing of drainage from the inundated centers of ice-wedge polygons in a warming climate have important implications for carbon flushing, advective heat transport, and transitions from methane to carbon dioxide dominated emissions. Here, we expand on previous research using a recently developed analytical model of drainage from a low-centered polygon. Specifically, we perform (1) a calibration to field data identifying necessary model refinements and (2) a rigorous model sensitivity analysis that expands on previously published indications of polygon drainage characteristics. This research provides intuition on inundated polygon drainage by presenting the first in-depth analysis of drainage within a polygon based on hydrogeological first principles. We verify a recently developed analytical solution of polygon drainage through a calibration to a season of field measurements. Due to the parsimony of the model, providing the potential that it could fail, we identify the minimum necessary refinements that allow the model to match water levels measured in a low-centered polygon. We find that (1) the measured precipitation must be increased by a factor of around 2.2, and (2) the vertical soil hydraulic conductivity must decrease with increasing thaw depth. Model refinement (1) accounts for runoff from rims into the ice-wedge polygon pond during precipitation events and possible rain gauge undercatch, while refinement (2) accounts for the decreasing permeability of deeper soil layers. The calibration to field measurements supports the validity of the model, indicating that it is able to represent ice-wedge polygon drainage dynamics. We then use the analytical solution in non-dimensional form to provide a baseline for the effects of polygon aspect ratios (radius to thaw depth) and coefficient of hydraulic conductivity anisotropy (horizontal to vertical hydraulic conductivity) on drainage pathways and temporal depletion of ponded water from inundated ice-wedge polygon centers. By varying the polygon aspect ratio, we evaluate the relative effect of polygon size (width), inter-annual increases in active-layer thickness, and seasonal increases in thaw depth on drainage. The results of our sensitivity analysis rigorously confirm a previous analysis indicating that most drainage through the active layer occurs along an annular region of the polygon center near the rims. This has important implications for transport of nutrients (such as dissolved organic carbon) and advection of heat towards ice-wedge tops. We also provide a comprehensive investigation of the effect of polygon aspect ratio and anisotropy on drainage timing and patterns, expanding on previously published research. Our results indicate that polygons with large aspect ratios and high anisotropy will have the most distributed drainage, while polygons with large aspect ratios and low anisotropy will have their drainage most focused near their periphery and will drain most slowly. Polygons with small aspect ratios and high anisotropy will drain most quickly. These results, based on parametric investigation of idealized scenarios, provide a baseline for further research considering the geometric and hydraulic complexities of ice-wedge polygons.

Harp, Dylan R. (ORCID:0000000197778000)↗
DOE OSTI2015

Particle tracking approach for transport in three-dimensional discrete fracture networks: Particle tracking in 3-D DFNs

The discrete fracture network (DFN) model is a method to mimic discrete pathways for fluid flow through a fractured low-permeable rock mass, and may be combined with particle tracking simulations to address solute transport. However, experience has shown that it is challenging to obtain accurate transport results in three-dimensional DFNs because of the high computational burden and difficulty in constructing a high-quality unstructured computational mesh on simulated fractures. We present a new particle tracking capability, which is adapted to control volume (Voronoi polygons) flow solutions on unstructured grids (Delaunay triangulations) on three-dimensional DFNs. The locally mass-conserving finite-volume approach eliminates mass balance-related problems during particle tracking. The scalar fluxes calculated for each control volume face by the flow solver are used to reconstruct a Darcy velocity at each control volume centroid. The groundwater velocities can then be continuously interpolated to any point in the domain of interest. The control volumes at fracture intersections are split into four pieces, and the velocity is reconstructed independently on each piece, which results in multiple groundwater velocities at the intersection, one for each fracture on each side of the intersection line. This technique enables detailed particle transport representation through a complex DFN structure. Verified for small DFNs, the new simulation capability enables numerical experiments on advective transport in large DFNs to be performed. As a result, we demonstrate this particle transport approach on a DFN model using parameters similar to those of crystalline rock at a proposed geologic repository for spent nuclear fuel in Forsmark, Sweden.

58 GEOSCIENCES↗
DOE OSTI2019

Achieving high strength and ductility in traditionally brittle soft magnetic intermetallics via additive manufacturing

Intermetallic alloys possess exceptional soft magnetic properties, including high permeability, low coercivity, and high saturation induction, but exhibit poor mechanical properties that make them impractical to bulk process and use at ideal compositions. We used laser-based Additive Manufacturing to process traditionally brittle Fe–Co and Fe–Si alloys in bulk form without macroscopic defects and at near-ideal compositions for electromagnetic applications. The binary Fe–50Co, as a model material, demonstrated simultaneous high strength (600–700 MPa) and high ductility (35%) in tension, corresponding to a ~300% increase in strength and an order-of-magnitude improvement in ductility relative to conventionally processed material. Atomic-scale toughening and strengthening mechanisms, based on engineered multiscale microstructures, are proposed to explain the unusual combination of mechanical properties. Here, this work presents an instance in which metal Additive Manufacturing processes are enabling, rather than limiting, the development of higher-performance alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Iron-modified biochar-based bilayer permeable reactive barrier for Cr(VI) removal

In this study, iron (Fe)-modified biochar (FeBC) has been developed to remove hexavalent chromium (Cr(VI)) from groundwater and is suitable for use in permeable reactive barriers (PRBs). However, Cr(VI) removal behavior and chemical processes in FeBC-based PRBs are not fully understood, and the potential for Fe release has not been addressed. In this study, three FeBC-based PRBs were assessed in column experiments for 563 days with respect to their ability to remove Cr(VI). Bilayer column filled with FeBC+limestone and BC+limestone in two separate layers (FeBC_Ca_BC) showed the best performance in terms of Cr(VI) removal with a low treatment cost. The corrosion of FeBC was mainly related to pH and Cr(VI) concentration rather than flow rate. Leached Fe was attenuated by BC and limestone and reutilized in FeBC_Ca_BC. Cr(VI) was reduced to Cr(III) and then adsorbed or precipitated on the biochars. Cr and Fe formed inner-sphere complexes and then transformed from double corner sharing to edge sharing. During the reaction, Cr penetrated from the surface to the interior of the biochars and became a more stable species. This study provides evidence of the effectiveness of a new combination of biochars for Cr(VI) removal and insights into the reaction mechanisms.

42 ENGINEERING↗
DOE OSTI2015

Design and installation of a ferromagnetic wall in tokamak geometry

Low-activation ferritic steels are leading material candidates for use in next-generation fusion development experiments such as a prospective component test facility and DEMO power reactor. Understanding the interaction of plasmas with a ferromagnetic wall will provide crucial physics for these facilities. In order to study ferromagnetic effects in toroidal geometry, a ferritic wall upgrade was designed and installed in the High Beta Tokamak–Extended Pulse (HBT-EP). Several material options were investigated based on conductivity, magnetic permeability, vacuum compatibility, and other criteria, and the material of choice (high-cobalt steel) is characterized. Installation was accomplished quickly, with minimal impact on existing diagnostics and overall machine performance, and initial results demonstrate the effects of the ferritic wall on plasma stability.

42 ENGINEERING↗
NASA NTRS1994

Sensible heat flux - Radiometric surface temperature relationship for eight semiarid areas

Measurements of sensible heat flux, radiometric surface temperature, air temperature, and wind speed made at eight semiarid rangeland sites were used to investigate the sensible heat flux-aerodynamic resistance relationship. The individual sites covered a wide range of vegetation (0.1-4 m tall) and cover (3%-95% bare soil) conditions. Mean values of k/B, a quantity related to the resistance of heat versus momentum transfer at the surface, for the individual sites varied between 3.5 and 12.5. A preliminary test of the utility of an excess resistance based on the mean value of k/B showed that the difference between the mean estimated and measured sensible heat fluxes varied +/- 60 W/sq m for the eight semiarid sites. For the eight sites the values of k/B were plotted against the roughness Reynolds number. The plot showed considerable scatter with values ranging between and beyond the theoretical curves for bluff rough and permeable rough surfaces.

Stewart, J. B.↗
DOE OSTI2017

Delineating Facies Spatial Distribution by Integrating Ensemble Data Assimilation and Indicator Geostatistics with Level Set Transformation.

A new approach is developed to delineate the spatial distribution of discrete facies (geological units that have unique distributions of hydraulic, physical, and/or chemical properties) conditioned not only on direct data (measurements directly related to facies properties, e.g., grain size distribution obtained from borehole samples) but also on indirect data (observations indirectly related to facies distribution, e.g., hydraulic head and tracer concentration). Our method integrates for the first time ensemble data assimilation with traditional transition probability-based geostatistics. The concept of level set is introduced to build shape parameterization that allows transformation between discrete facies indicators and continuous random variables. The spatial structure of different facies is simulated by indicator models using conditioning points selected adaptively during the iterative process of data assimilation. To evaluate the new method, a two-dimensional semi-synthetic example is designed to estimate the spatial distribution and permeability of two distinct facies from transient head data induced by pumping tests. The example demonstrates that our new method adequately captures the spatial pattern of facies distribution by imposing spatial continuity through conditioning points. The new method also reproduces the overall response in hydraulic head field with better accuracy compared to data assimilation with no constraints on spatial continuity on facies.

97 MATHEMATICS AND COMPUTING↗