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 883 records · Page 49

Upgrading Biogas through in situ Conversion of Carbon Dioxide to Biomethane in Anaerobic Digesters

Organic waste streams generated by wastewater treatment plants, agricultural operations, and food processing industries represent an important yet underutilized opportunity for renewable energy production in the United States. Through anaerobic digestion, these waste streams can produce biogas, a mixture primarily composed of methane (CH4) and carbon dioxide (CO2), that can be upgraded to pipeline-quality natural gas. However, most existing upgrading technologies remove CO2 from biogas rather than utilizing it, leaving a significant portion of the potential energy unused. This project investigates a novel biological upgrading approach that converts CO2 into additional CH4 by supplying hydrogen (H2) to specialized microorganisms capable of performing hydrogenotrophic methanation. The main challenges associated with biological biogas upgrading are related to hydrogen supply, gas-liquid mass transfer, and process stability. First, due to the high cost of hydrogen gas, it is preferable that H2 be produced on-site using renewable energy sources such as wind or solar power. Second, hydrogen has low solubility in liquids, which limits its availability to microorganisms and requires strategies to improve gas dissolution and transfer within the reactor. Third, process inhibition may occur as a result of increased pH caused by CO2 consumption or elevated H2 partial pressure, both of which can negatively affect methanogenic activity. Although research in these areas has advanced during the course of this project, these challenges have not yet been fully resolved. To date, the biological systems that have achieved the highest methane concentrations are typically ex-situ reactors, where operational conditions can be more easily controlled. For this reason, the findings of the present project remain highly relevant. The project goal was to develop an innovative system that can accomplish biogas upgrading via biological conversion of CO2 to CH4, in a novel hybrid approach that combines the advantages of both in-situ and ex-situ systems. The proposed system employs a three-phase upflow anaerobic bioreactor with H2 delivery through a gas-permeable membrane, enabling efficient hydrogen transfer and microbial conversion. Under optimized operating conditions, the system achieved 99% H2 consumption and 90% CO2 conversion. A subsequent gas cleaning stage was implemented to further improve gas quality and meet target purity standards. The upgraded gas composition reached 97.7% CH4, 2.2% CO2, and 0.97% O2, while H2S concentrations remained below detection limits. In addition, a flue gas-driven inorganic thermoelectric generator (TEG) system was designed and experimentally validated as a potential source of electricity for H2 production. The system consisted of six TEG modules connected in series and achieved an open-circuit voltage of 4.5 V and a maximum power output of 224 mW at a temperature difference of approximately 53.5 °C, demonstrating effective conversion of waste heat into electrical power under simulated flue gas conditions. Finally, a comprehensive techno-economic analysis was completed to evaluate the capital and operating costs associated with the proposed system. The results provide important insights to guide future scale-up, optimization, and potential deployment of integrated biological biogas upgrading technologies.

09 BIOMASS FUELS↗

Repurposing Drilling Control Diagnostics for Subsurface Edge Detection and Boundary Advisement During Planetary Drilling

Informed decision-making during lunar drilling and sampling missions will require data monitoring tools and specialized ground data systems. Accurate and updated situational awareness, with ongoing data monitoring, is critical for timely responses by to incoming science data. Traverse plans and scheduled activities may need to be flexibly changed in order to react to unexpected data or situations. Unlike (for example) Mars missions, the relative lightspeed closeness of the Moon allows for near-real-time ground processing of incoming mission and instrument data. An Apollo-class lunar regolith drill will in a sense “travel” a meter or two vertically at a given subsurface characterization site. As the drill penetrates into lunar regolith, it is likely to encounter a range of material densities, orientations, fracture toughness, and (perhaps) ice percentages. Lunar drill telemetry can provide science teams with a valuable first look into the subsurface structure, the regolith bulk properties, and constituents at each drilled site. Real-time AI-based recognition and reaction to downhole situations has been developed for automated deeper drilling on Mars and beyond. We can leverage the same knowledge bases and pattern-matching as areal-time interpreter of the subsurface, a situational awareness tool during drilling operations. We recently (Sept. 2019) demonstrated this AI drilling monitoring and analysis capability, in control of in-situ drilling and sampling operations, mounted on a KREX-2 rover in Chile’s Atacama Desert. Terrestrial automated drilling log analyses in oil exploration have used similar machine learning techniques in classifying and identifying features in drilling logs –but these typically are designed assuming a drilling fluid influencing downhole measurements and data (permeability, resistivity). Drilling models and existing AI software designed to detect and respond to drilling faults and hard materials can be repurposed, for near-real-time (ground-based) interpretation of drilling telemetry –a potentially valuable advisory tool for strata boundaries and changes in drilling parameters. On the Moon, this approach could be used to study the structure and to some extent the composition of lunar regolith vs. borehole depth, based on recognizable variations in fracture hardness, drilling energy and penetration rates while actively drilling. Since the early 2000s, a series of increasingly-capable real-time drilling telemetry interpretation and characterization software tools have been developed. These subsurface models and software tools have monitored the real-time drilling data received, and automatically identified changes in drill behavior (e.g., encountering a harder target layer, bit inclusions, drill choking due to infall downhole, and others) correlating these with subsurface structures and features. We discuss the mappings between drill borehole parameters, faults or events detected, and modeled changes in rock layer boundaries, in examples drawn from field testing at analog sites in an Arctic impact crater, Rio Tinto, and Chile’s Atacama Desert. These demonstrate how subsurface structural boundaries led to fault detections and responses by the software.

drilling advisor↗

DEEPEN 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...

15 GEOTHERMAL ENERGY↗

Effect of iron and silica nanoparticles’ size on in vitro human skin binding and penetration

This in vitro skin study determined absorption, diffusion, and binding rates of four [14C]-labeled nanoparticles (NPs): 12 nm Fe 3 O 4 , 32 nm Fe 3 O 4 @SiO 2 , 33 nm SiO 2 , and 78 nm SiO 2 in each layer of human cadaver skin. In vitro microdialysis device and flow-through skin diffusion system were used to measure the binding affinity to the stratum corneum (SC) and permeability into/through skin layer of the four NPs with different physical–chemical properties, respectively, in short (30 min) and/or long (24 hours) exposures. Results show that NP size is an important factor affecting NP percutaneous absorption. The 12 nm Fe 3 O 4 NPs reached the SC and viable epidermis; 32 nm Fe 3 O 4 @SiO 2 core/shell NPs only reached SC. However, 33 nm and 78 nm silica NPs did not permeate SC. Similar patterns were observed for NP binding affinity to SC and dermatopharmacokinetic analysis using the tape stripping method. The binding affinity determination may be a useful method to efficiently screen skin penetration of NPs.

59 BASIC BIOLOGICAL SCIENCES↗

Multiscale Characterization of Structural Compositional and Textural Heterogeneity of Nano-porous Geomaterials

The purpose of the project was to perform multiscale characterization of low permeability rocks to determine the effect of physical and chemical heterogeneity on the poromechanical and flow responses of shales and carbonate rocks with a broad range of physical and chemical heterogeneity . An integrated multiscale imaging of shale and carbonate rocks from nanometer to centimeter scales include s dual focused ion beam - scanning electron microscopy (FIB - SEM) , micro computed tomography (micro - CT) , optical and confocal microscopy, and 2D and 3D energy dispersive spectroscopy (EDS). In addition, mineralogical mapping and backscattered imaging with nanoindentation testing advanced the quantitative evaluat ion of the relationship between material heterogeneity and mechanical behavior. T he spatial distribution of compositional heterogeneity, anisotropic bedding patterns, and mechanical anisotropy were employed as inputs for brittle fracture simulations using a phase field model . Comparison of experimental and numerical simulations reveal ed that proper incorporation of additional material information, such as bedding layer thickness and other geometrical attributes of the microstructures, can yield improvements on the numerical prediction of the mesoscale fracture patterns and hence the macroscopic effective toughness. Overall, a comprehensive framework to evaluate the relationship between mechanical response and micro-lithofacial features can allow us to make more accurate prediction of reservoir performance by developing a multi - scale understanding of poromechanical response to coupled chemical and mechanical interactions for subsurface energy related activities.

36 MATERIALS SCIENCE↗

Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies

The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress in modeling and experimental approaches towards the characterization of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in deep clay/shale/argillaceous rock. International collaboration activities such as heater tests and postmortem analysis of samples recovered from these have elucidated key information regarding changes in the engineered barrier system (EBS) material exposed to years of thermal loads. Chemical and structural analyses of sampled bentonite material from such tests has as well as experiments conducted on these are key to the characterization of thermal effects affecting bentonite clay barrier performance and the extent of sacrificial zones in the EBS during the thermal period. Thermal, hydrologic, and chemical data collected from heater tests and laboratory experiments has been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches (e.g., continuum vs. discrete) to tackle issues related to flow and transport at various scales of the host-rock and EBS design concept. Consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the EBS design. This is particularly important for backfilled repository concepts where temperature plays a key role in the EBS behavior and long-term performance. This report describes multiple R&D efforts on disposal in argillaceous geologic media through development and application of coupled THMC process models, experimental studies on clay/metal/cement barrier and host-rock (argillite) material interactions, molecular dynamic (MD) simulations of water transport during (swelling) clay dehydration, first-principles studies of metaschoepite (UO 2 corrosion product) stability, and advances in thermodynamic plus surface complexation database development. Drift-scale URL experiments provides key data for testing hydrological-chemical (HC) model involving strong couplings of fluid mixing and barrier material chemical interactions. The THM modeling focuses on heater test experiments in argillite rock and gas migration in bentonite as part of international collaboration activities at underground research laboratories (URLs). In addition, field testing at an URL involves in situ analysis of fault slip behavior and fault permeability. Pore-scale modeling of gas bubble migration is also being investigated within the gas migration modeling effort. Interaction experiments on bentonite samples from heater test under ambient and elevated temperatures permit the evaluation of ion exchange, phase stability, and mineral transformation changes that could impact clay swelling. Advances in the development, testing, and implementation of a spent nuclear fuel (SNF) degradation model coupled with canister corrosion focus on the effects of hydrogen gas generation and its integration with Geologic Disposal Safety Assessment (GDSA). GDSA integration activities includes evaluation of groundwater chemistries in shale formations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DEEPEN: Final 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

Part of the DEEPEN (DE-risking Exploration of geothermal Plays in magmatic ENvironments) project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology, detailed in the journal article below, is based on the method proposed by Poux & O'brien (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., (2022). For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty.

15 GEOTHERMAL ENERGY↗

Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions

Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.

02 PETROLEUM↗

Tensile Properties of Hydrogels and of Snake Skin

Stimulus-responsive or 'smart' gels are of potential interest as sensors and actuators, in industrial separations, and as permeable delivery systems. In most applications, a certain degree of mechanical strength and toughness will be required, yet the large-strain behavior of gels has not been widely reported. Some exceptions include work on gelatin and other food gels, some characterization of soft gels applicable for in-vitro cell growth studies, and toughness determinations on commercial contact lens materials. In general, it can be anticipated that the gel stiffness will increase with increasing degree of crosslinking, but the tensile strength may go through a maximum. Gel properties can be tailored by varying not only the degree of crosslinking, but also the polymer concentration and the nature of the polymer backbone (e.g. its stiffness or solubility). Polypeptides provide an especially interesting case, where secondary structure affects trends in moduli and conformational transitions may accompany phase changes. A few papers on the tensile properties of responsive gels have begun to appear. The responsive hydrogel chosen for the present study, crosslinked hydroxypropylcellulose, shrinks over a rather narrow temperature range near 44 C. Some vertebrate skin is also subject to substantial strain. Among reptiles, the morphologies of the skin and scales show wide variations. Bauer et al. described the mechanical properties and histology of gecko skin; longitudinal tensile properties of snake skin were examined by Jayne with reference to locomotion. The present measurements focus on adaptations related to feeding, including the response of the skin to circumferential tension. Tensile properties will be related to interspecific and regional variation in skin structure and folding.

Hinkley, Jeffrey A.↗

Insights into scale translation of methane transport in nanopores

Accurate prediction of flow behavior in shale matrix is critical for efficient development of shale gas reservoirs. In these systems, the majority of pores are in the nano-size range. As a result, continuum-based approaches may not be appropriate to simulate flow in such systems. Molecular dynamics (MD) simulations are capable of capturing the relevant microscale physics. Their relatively high computational expense, however, restricts MD simulations to rather small systems and domains. This limitation creates a gap between computational need of macroscale systems and capabilities of MD simulations. The lattice Boltzmann method (LBM) is a suitable candidate to bridge this gap. In this work, the multiple-relaxation-time (MRT)-LBM is used to study methane transport in nano-size pores. Adsorption effects near solid boundaries, as well as non-ideal behavior of fluids, are accounted for via incorporating appropriate force terms in LBM. In this work, parameters associated with the force terms in the equation of state are studied in detail, and a workflow is proposed to determine optimal values of these parameters for gas flow in slit pores. Specifically, we establish these parameters such that the range of density values that the model is able to simulate is maximized. We demonstrate this workflow by simulating gas flow where velocity and density profiles from MD simulations are used as reference data. Results from LBM simulations are in good agreement with MD reference data for pores that are 4 nm in width or larger. Moreover, we propose a preconditioning scheme to improve the stability of LBM in dealing with complex geometries. The robustness of this scheme is demonstrated by simulating several roughness geometries. This work motivates the use of LBM in scale translation of the physics of mass transport in more complex permeable media.

03 NATURAL GAS↗

The impact of mineral reactive surface area variation on simulated mineral reactions and reaction rates

Reactive transport modeling is an essential tool to simulate complex geochemical reactions in porous media that can impact formation properties including porosity and permeability. However, simulating these reactions is challenging due to uncertainties in model parameters, particularly mineral surface areas. Imaging has emerged as a powerful means of estimating model parameters including porosity and mineral abundance, accessibility and accessible surface area. However, these parameters, particularly mineral accessible surface area, vary with image resolution. This work aims to enhance understanding of the impact of image resolution and other means of estimating mineral reactive surface area on simulated mineral reactions and reaction rates. Mineral surface areas calculated from images with resolutions of 0.34 μm and 5.71 μm were used to simulate mineral reactions in the context of geologic CO 2 sequestration in the Paluxy formation at the continuum scale. Additional simulations were carried out using BET surface areas collected from the literature and geometric surface areas. Simulations were run for 7300 days and mineral volume fractions and effluent ion concentrations tracked and compared. Variations in mineral surface areas measured from images are within 1 order of magnitude and yield similar simulation results, indicating the impact of image resolution on simulated reactions and reaction rates is minimum for the resolutions and sample considered. In comparison, surface areas obtained from BET and geometric approaches are 1–4 orders of magnitude higher than image-obtained surface areas and result in greater simulated reaction rates and extents. Minerals with high reaction rates (calcite and siderite) are most impacted by surface area values at short times where simulated mineral volume fractions at longer times agree relatively well, even for simulations with several orders of magnitude variation in surface area. Phases with lower reaction rates, such as K-feldspar and muscovite, are predominantly impacted over longer times where variations in surface areas impact reaction extents and porosity evolution. Overall, variations in surface areas due to image resolution are small and result in little variation in simulated reactions and reaction rates while there are significant variations in simulation results when other surface area estimates are used. These variations, however, depend on the reactivity of the mineral phase where surface areas of fast-reacting phases largely impact simulations at short (hours to days) timescales and surface areas of slower-reactive phases impact longer term simulations (years).

58 GEOSCIENCES↗

Prototyping and Testing a New Volumetric Curvature Tool for Modeling Reservoir Compartments and Leakage Pathways in the Arbuckle Saline Aquifer: Reducing Uncertainty in CO 2 Storage and Permanence

This DOE-funded project evaluates the utility of seismic volumetric curvature (VC) for predicting stratal and structural architecture diagnostic of paleokarst reservoirs. Of special interest are applications geared toward carbon capture, utilization, and storage (CCUS). VC has been championed for identifying faults (offset <¼ λ) that cannot be imaged by conventional 3-D seismic attributes such as coherence. The objective of this research was to evaluate VC-techniques for reducing uncertainties in reservoir compartmentalization studies and seal risk assessments especially for saline aquifers. A 2000-ft horizontal lateral was purposefully drilled across VC-imaged lineaments—interpreted to record a fractured and a fault-bounded doline—to physically confirm their presence. The 15-mi² study area is located in southeastern Bemis-Shutts Field, which is situated along the crest of the Central Kansas Uplift (CKU) in Ellis County, Kansas. The uppermost Arbuckle (200+ ft) has extensive paleokarst including collapsed paleocaverns and dolines related to exceedingly prolonged pre-Simpson (Sauk–Tippecanoe) and/or pre-Pennsylvanian subaerial exposure. A lateral borehole was successfully drilled across the full extent (~1100 ft) of a VC-inferred paleokarst doline. Triple combo (GR-neutron/density-resistivity), full-wave sonic, and borehole micro-imager logs were successfully run to TD on drill-pipe. Results from the formation evaluation reveal breccias (e.g., crackle, mosaic, chaotic), fractures, faults, vugs (1-6"), and unaffected host strata consistent with the pre-spud interpretation. Well-rounded pebbles were also observed on the image log. VC-inferred lineaments coincide with 20–80-ft wide intervals of high GR values (100+ API), matrix-rich breccias, and faults. To further demonstrate their utility, VC attributes are integrated into a geocellular modeling workflow: 1) to constrain the structural model; 2) to generate facies probability grids, and; 3) to collocate petrophysical models to separate-vug rock fabrics along solution-enlarged fault and fracture systems. Simulation-based studies demonstrate a potential alternative field development model for developing CO 2 storage sites that target carbonate reservoirs overprinted by paleokarst. Simulation results for this complex reservoir indicate that individual fault blocks could function as discrete containers for CO 2 storage thereby reducing the risk of plume migration outside the legally defined extent of the permitted storage site. Vertically extensive, anastomosing, solution-enlarged fault/fracture systems — infilled by clay-rich sediments — would operate as non-to-low permeability vertical "curtains" that restrict CO 2 movement beyond the confines of the CO 2 storage site. Such a location could be developed in a checker-board fashion with CO 2 injection operations occurring in one block and surveillance operations occurring in the adjacent block. Such naturally partitioned reservoirs may be ideal candidates for reducing risks associated with CO 2 plume breakthrough.

02 PETROLEUM↗

High-Reliability Waveguide Vacuum/Pressure Window

The NASA Deep Space Network (DSN) uses commercial waveguide windows on the output waveguide of Ka-band (32 GHz) low-noise amplifiers. Mechanical failure of these windows resulted in an unacceptable loss in tracking time. To address this issue, a new Ka-band WR-28 waveguide window has been designed, fabricated, and tested. The window uses a slab of low-loss, low-dielectric constant foam that is bonded into a 1/2-wave-thick waveguide/flange. The foam is a commercially available, rigid, closed-cell polymethacrylimide. It has excellent electrical properties with a dielectric constant of 1.04, and a loss tangent of 0.01. It is relatively strong with a tensile strength of 1 MPa. The material is virtually impermeable to helium. The finished window exhibits a leak rate of less than 3x10(exp -3)cu cm/s with helium. The material is also chemically resistant and can be cleaned with acetone. The window is constructed by fabricating a window body by brazing a short length of WR-28 copper waveguide into a standard rectangular flange, and machining the resulting part to a thickness of 4.6 mm. The foam is machined to a rectangular shape with a dimension of 7.06x3.53 mm. The foam is bonded into the body with a two-part epoxy. After curing, the excess glue and foam are knife-trimmed by hand. The finished window has a loss of less than 0.08 dB (2%) and a return loss of greater than 25 dB at 32 GHz. This meets the requirements for the DSN application. The window is usable for most applications over the entire 26-to-40-GHz waveguide band. The window return loss can be tuned to a required frequency by var y in g the thickness of the window slightly. Most standard waveguide windows use a thin membrane of material bonded into a recess in a waveguide flange, or sandwiched between two flanges with a polymer seal. Designs using the recessed window are prone to mechanical failure over time due to constraints on the dimensions of the recess that allow the bond to fail. Designs using the sandwich method are often permeable to helium, which prohibits the use of helium leak detection. At the time of this reporting, 40 windows have been produced. Twelve are in operation with a combined operating time of over 30,000 hours without a failure.

Britcliffe, Michael J.↗

The non-Newtonian heat and mass transport of He 2 in porous media used for vapor-liquid phase separation

This investigation of vapor-liquid phase separation (VLPS) of He 2 is related to long-term storage of cryogenic liquid. The VLPS system utilizes porous plugs in order to generate thermomechanical (thermo-osmotic) force which in turn prevents liquid from flowing out of the cryo-vessel (e.g., Infrared Astronomical Satellite). An apparatus was built and VLPS data were collected for a 2 and a 10 micrometer sintered stainless steel plug and a 5 to 15 micrometer sintered bronze plug. The VLPS data obtained at high temperature were in the nonlinear turbulent regime. At low temperature, the Stokes regime was approached. A turbulent flow model was developed, which provides a phenomenological description of the VLPS data. According to the model, most of the phase separation data are in the turbulent regime. The model is based on concepts of the Gorter-Mellink transport involving the mutual friction known from the zero net mass flow (ZNMF) studies. The latter had to be modified to obtain agreement with the present experimental VLPS evidence. In contrast to the well-known ZNMF mode, the VLPS results require a geometry dependent constant (Gorter-Mellink constant). A theoretical interpretation of the phenomenological equation for the VLPS data obtained, is based on modelling of the dynamics of quantized vortices proposed by Vinen. In extending Vinen's model to the VLPS transport of He 2 in porous media, a correlation between the K*(GM) and K(p) was obtained which permits an interpretation of the present findings. As K(p) is crucial, various methods were introduced to measure the permeability of the porous media at low temperatures. Good agreement was found between the room temperature and the low temperature K(p)-value of the plugs.

Yuan, S. W. K.↗

Effects of ingested nanocellulose on intestinal microbiota and homeostasis in Wistar Han rats

Micron scale crystalline and fibrillar cellulose materials are “generally regarded as safe” (GRAS) as binders and thickeners in food products. Yet, partly due to a lack of relevant toxicological data, nanocellulose (NC) materials, which have unique properties that can be exploited to improve food quality and safety, have yet to receive FDA approval as food ingredients. Results of in vitro and in vivo toxicological studies of ingested NC, detailed in a recent companion report, revealed minimal acute in vitro cytotoxicity, and no toxicity in a subacute rat gavage model, suggesting that NC materials are non-hazardous. However, ingested materials may also modulate gut microbial populations, or alter aspects of intestinal function not evaluated in standard totoxicity testing, which could have important health implications. Here, we report the results of studies of the effects of ingested cellulose nanofibrils (CNF) on the fecal microbiome and metabolome, intestinal (ileal) epithelial cell mRNA expression of cell junction genes, and ileal cytokine production. Feces, plasma, and ileum samples were collected from Wistar Han rats before and after five weeks of biweekly gavages of water or cream, with or without 1% w/w CNF. Analysis of fecal microbial populations revealed that CNF caused changes in both genus and species diversity, with specific effects on species that produce short chain fatty acids, and that have been associated with increased IgA production and elevation of insulin levels. Fecal metabolomic analysis revealed relatively few effects of CNF, with significant changes in the levels of only ten metabolites out of 366 measured. Exposure to CNF also caused differential regulation of mRNA expression of several genes involved in epithelial cell junctions, and increased production of cytokines that both promote and inhibit proliferation of CD8 T cells. These perturbations in intestinal homeostasis contrast somewhat with the absence of in vitro and in vivo toxicity observed previously but would appear to represent minor effects. Additionally, the gut microbiome is highly sensitive to ingested substances, and such disturbances of the intestinal microbial ecosystem do not necessarily represent or predict meaningful pathology. Further studies, including chronic feeding studies, are needed to assess the real health implications, if any, of these changes.

59 BASIC BIOLOGICAL SCIENCES↗

Establishing a silica gel zone in well annulus and evaluating its performance in blocking vertical water flow

Wells are often constructed for monitoring purposes with relatively long screen lengths (e.g., >10 m). Vertical water flows can occur within the artificial or natural filterpack annulus that surrounds the screened interval, bypassing packer assemblies installed inside the wellbore. Attempts to isolate discrete vertical zones during groundwater sampling are unsuccessful when annular flow occurs and lead to remedy decisions based on biased or incorrect interpretations. Blocking vertical annular water flow and contaminant transport will help obtain more accurate concentrations of contaminants from sampling in targeted depth intervals. The application of silica gels formed from the injected colloidal silica CS suspensions is a novel approach to minimize or prevent movement of vertical movement of groundwater in the surrounding filterpack annulus. In this work, we tested the feasibility of injecting CS suspensions to target locations and developed a modified CS formulation that is injectable and prevents gravity sinking. We studied the distribution and penetration of silica gel at laboratory scale in mock well annulus with surrounding formations. We evaluated the performance of the silica gel in blocking vertical water flow in the annulus and in minimizing chemical transport through the gel zone. CS suspension formulations have been defined that are ready for injection, stay in target locations, and form gel within desired time frames. Injection of CS suspensions achieved uniform distribution in a well annulus filter pack, fully occupied the annulus pore space, and penetrated the formation surrounding the filter packer with a sufficient distance to create a hydraulic annular seal when the injection was applied at a sufficient rate. The depth of penetration into the formation was dependent on the permeability contrast between the filter pack and the surrounding formation. Silica gel that formed in the annulus blocked vertical water flow and stopped the chemical transport through the gel zone. In conclusion, this research reveals that using CS suspension injection and sequential gelation (CS-GEL) is a promising technology for blocking vertical water flow and chemical transport through the filter pack in targeted zones within the annulus of long-screened well systems.

Colloidal silica suspension↗

Ma'adim Vallis Estuarine Delta in Elysium Basin and Its Relevance as a Landing Site for Exobiology Exploration on Mars

The debouche of Ma'adim Vallis in the Elysium Basin generated a transitional transported sediment structure, which planimetric shape is controlled by the enclosing topography of a deep reentrant gulf of the Basin into the highland. We defined it as an estuarine delta. The location and the importance of this estuarine delta is supported by the theoretical model of graded profile constructed for Ma'adim Vallis, and by two approaches: (i) the reconstruction of Ma'adim Vallis downstream course from Gusev to Elysium Basin, and (ii) the survey of the sediment deposit in the alleged estuary. The longitudinal graded profile of Ma'adim Vallis finds its base-level in the Elysium Basin, at a about 1000 m elevation, which is in agreement with the observed Basin shoreline. This model is supported by observational evidence of flow between the northern rim of Gusev crater, and the Elysium Basin shoreline. This downstream course of Ma'adim Vallis can be divided into three hydrogeologic regions. into three hydrogeologic regions. (a) The first region is a flooded plain (Zephiria Mensae), consisting in chaotic terrain formed by highland rocks, and disintegrated lava of the western flank of Apollinaris. Morphologic indicators of the flood process are: (1) the sediment deposit over the Gusev crater northern rim that reflects the overspilling of the crater-lake water through a 40-km wide gap provided by an ancient impact crater, (2) the tear-drop shaped feature on the northeastern flank of Apollinaris Patera, and (3) the chaotic terrain that suggest the emergence of ground water generated by the seepage of the crater lake through high-permeable broken rampart material. This underground water circulation sustained by the hydrostatic pressure of the crater-lake has likely generated a hydrothermal system in the volcanic environment of Apollinaris Patera. The stratigraphy of the flooded area is identified as Hesperian age, with occurrences of Noachian hilly individual features, and as Amazonian flooded plain and chaotic material, (b) The second region is located on the western flank of Apollinaris Patera. It is surrounded by relics of deep valleys that suggest a former downstream course of Ma'adim Vallis. The geologic setting of this region (Lucus Planum) is interpreted to be an Amazonian formation composed by the middle and lower members of the Medusae Formation., c) The third region corresponds to the convergence of the west and east branches of Ma'adirn Vallis into a deep re-entrant wide gulf that penetrates about 100 km into the highland. This topographic depression is delineated by the 1000 in elevation contour. This gulf has formed an estuarine configuration centered at 3S/190W within the Elysium Basin. This configuration has favored the formation of a estuarine sedimentary delta, because of topographically controlled lateral migration. This estuarine structure is strongly dominated by the incoming supply of Ma'adim Vallis fluvial sediment extracted from Zephiria Mensae and Lucus Planum. The obtuse-angle geometry of the estuary increases the sedimentation rate, which is higher than in the course of the channel. The sediment deposition process is governed by the estuarine water circulation. The inflowing loaded fluvial water enters the estuary as a bottom current, and mixes with the relatively less-loaded water of the receiving basin. When they mixed. the inflowing fluvial material, and the landward basin circulating water generate an accumulation of highly-diversified estuarine deposit stratification. This accumulation of material is mostly centered in the transitional zone of the delta. The sediment trapping efficiency of the estuary is function of the energy balance between the inflowing fluvial water, and the ingoing basin current. The submergence of the delta by the rising of the water-level increases the estuary water-depth, and consequently the sediment entrapment is favored. The locus of sediment accumulation moves landward in the zone of inflowing fluvial water. This results in the rising of the channel base-level, thus in the increase of the length of the longitudinal graded-profile. The sediment deposit facies of the zone A shows a generally smooth surface. The longitudinal deposit is bordered by alluvial terraces that reflect the variations of the channel level. The waning of the Elysium Basin caused the erosion of the Basin estuarine zone by small channels, this episode being characterized by dissected tear-drop shaped mesa-like morphologies in the delta. Our estuarine delta model predicts a lithostratigraphic depositional sequence associated with the water submergence and the transgression of Elysium Basin. The thickness of the estuarine sediment corresponds to the Elysium Basin levels changes relatively to the bed floor of the estuary, The depositional sequence of Ma'adim Vallis are described: (1) a pro-current filled region (A), where fluvial are longitudinally accumulated by the inflowing water, (2) inverse current from Elysium Basin (B), where fluvial and lacustrine sediments are accumulated, and (3) zone of current equilibrium (C), where the sediments are distributed as a shoreline at the boundary of the estuarine delta. The estuary sedimentology dynamics collects and keeps the record of the geologic unit material crossed by Ma'adim Vallis, and those of the lakebed deposit of Elysium Basin. The predicted mixed stratigraphic sequence from fluvial and lacustrine sediment makes this site an exceptional environment to concentrate potential multi-origin biologic records. We envision four possible strategies to explore this sedimentologic record: (1) longitudinal surface and subsurface traverses in region A to investigate outcrop levees, (2) exploration of the mesa walls in region B, (3) deep drilling hole lodging of the sequential deposits in the zones A and B, and (4) surface and subsurface exploration of the shoreline delta. The expected results for each of these strategies are: (1) in the deepest layers of region A are predicted frequent and abundant coarse material, sandy lenses lamination grading downward from sand to cobbles. Volcanic debris from the Noachian crustal Plateau unit material, hydrothermal altered rocks, carbonates, Hesperian and possibly Amazonian volcanic material, from Apollinaris Patera, altered rocks and carbonates from Zephiria Mensae are expected. As a favorable environment for inception of life, possible biological records are expected in transported rock, (2) At the surface, and subsurface (<=100 m), large deposits sandy to silted material from Elysium paleolake basin mixed with fine-grained sediments from Ma'adim Vallis are expected mostly in the upstream part of region B, (3) on the shoreline of the estuarine delta, abundant fine material from Elysium paleolake basin (evaporites, carbonates), mostly Amazonian in age are expected. The Ma'adini estuary is a favorable landing site for all the above mentioned science aspects, and .for its location. The site lies near the equator, which is favorable for the rover solar power supply, and at 1000m elevation, which is a favorable configuration for the descent system braking. Another advantage is the extent of the area of high scientific interest (33,000 sq km), which provides a good ellipse, and potential long study traverses.

Grin, E. A.↗

Carbon Dioxide Shuttling Thermochemical Storage Using Strontium Carbonate

Phase I concludes with significant progress made towards the SunShot ELEMENTS goals of high energy density, high power density, and high temperature by virtue of a SrO/SrCO 3 based material. A detailed exploration of sintering inhibitors has been conducted and relatively stable materials supported by YSZ or SrZO 3 have been identified as the leading candidates. In 15 cycle runs using a 3 hour carbonation duration, several materials demonstrated energy densities of roughly 1500 MJ/m 3 or greater. The peak power density for the most productive materials consistently exceeded 40 MW/m 3 —an order of magnitude greater than the SOPO milestone. The team currently has a material demonstrating nearly 1000 MJ/m 3 after 100 abbreviated (1 hour carbonation) cycles. A subsequent 8 hour carbonation after the 100 cycle test exhibited over 1500 MJ/m 3 , which is evidence that the material still has capacity for high storage albeit with slower kinetics. Kinetic carbonation experiments have shown three distinct periods: induction, kinetically-controlled, and finally a diffusion-controlled period. In contrast to thermodynamic equilibrium prediction, higher carbonation temperatures lead to greater conversions over a 1 hour periods, as diffusion of CO 2 is more rapid at higher temperatures. A polynomial expression was fit to describe the temperature dependence of the linear kinetically-controlled regime, which does not obey a traditional Arrhenius relationship. Temperature and CO 2 partial pressure effects on the induction period were also investigated. The CO 2 partial pressure has a strong effect on the reaction progress at high temperatures but is insignificant at temperatures under 900°C. Tomography data for porous SrO/SrCO 3 structures at initial stage and after multiple carbonation/decomposition cycles have been obtained. Both 2D slices and 3D reconstructed representations have been obtained from the raw data. Porosity and surface area have been computed based on tomography. The two parameters, permeability and Forchheimer constants, in the Darcy’s law for linear and non-linear range have been obtained from the fluid dynamics simulation. The structure is nearly isotropic from the hydrodynamics point of view. Thermal conductivity at elevated temperature has been obtained from the 3-D pore-scale simulations. The computed thermal conductivities of the structure compare well with an existing correlation. Computational models have been developed for solving the energy and mass transport equations with chemical reaction on the surface and inside the particles. The chemical kinetics or rate law is needed on the particle scale; however it is generally measured using a lumped model based on the average of reacted as well as unreacted mass. Many simulations have been conducted on a 2-D model to determine the effects of those pore scale parameters on the heat release based on the measured data. The simulation results using 2-D model show that average volumetric heat flow rate of 4 MW/m 3 in the first 5 minutes can be obtained. The insight gained from the 2-D simulation will be useful for the next stage simulation based on realistic 3-D structures. A probabilistic analysis has been conducted on the baseline system of recovering heat from the hot CO 2 evolved from the endothermic reaction and compressing the cooled CO 2 into a storage tank at 20-21 bar. Using a Monte Carlo simulation using 2,000,000 samples and both uniform and triangular probability distributions of four variables (heat exchanger pinch points, compressor efficiency, and power block efficiency), exergetic efficiency of 90% appears to be very likely.

14 SOLAR ENERGY↗