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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Effects of time-dependent deformation of shale on the integrity of a geological nuclear waste repository

Safety assessment of geological nuclear waste repositories is essential for the permanent disposal of spent nuclear fuels and high-level radioactive waste. The long-term integrity of the host rock as well as the engineered barrier system (e.g., bentonite buffer) surrounding nuclear waste canisters is of particular importance as decay heat from nuclear waste canisters, which will last over thousands of years, significantly disturbs the thermo-hydromechanical (THM) state of the repository. In this study, THM coupled simulations were carried out to investigate the effect of time-dependent deformation (i.e., creep) of shale on the long-term integrity of a generic subsurface nuclear waste repository. The Norton-Bailey creep model, which is also known as the Lemaitre-Menzel-Schreiner model, was employed to simulate the power-law type creep that is observed in shales. The TOUGH-FLAC simulator was employed for the THM coupled modeling of the repository. The objective of this study is to assess the effect of creep in different shales (i.e., high creep shale vs. low creep shale) on long-term stress and permeability changes in the repository. Results show potential advantages of constructing repositories in the high creep shale, as deviatoric stress levels in the formation decreased to zero in 100 years since the emplacement of nuclear waste canisters and the permeability also decreased to the undamaged, intrinsic levels in 10,000 years. Also, mean effective stress levels in the bentonite buffer increased by 100% in the high creep shale case relative to the low creep shale case at 10,000 years due to creep-induced contraction of the nuclear waste disposal tunnel. However, in earlier periods (e.g., 1000 years), the stress levels in the bentonite buffer were twice smaller in the high creep shale case than in the low creep shale case, which shows a tradeoff between the intermediate- (~1000 years) and long-term (>10,000 years) compaction levels in the bentonite buffer depending on the creep characteristics of the host shale.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2016

Geologic Modeling and Simulation Report for the Aquistore Project

In this update, a regional-scale model extends beyond the 34-square-kilometer PTRC 3-D seismic survey area. A local grid refinement (LGR) system near both the injection and observation wells was introduced for the history-matching and uncertainty analysis. Spinner log survey and pressure test data provided by partners were used and evaluated to adjust the nearwellbore local permeability in order to history-match the field pressure data. As of January 9, 2016, approximately 24,000 tonnes of CO 2 has been injected during a series of relatively short operating periods. However, the injection rate during these short operating periods increased near the end of year 2015 to a level of 300 to 500 tonnes per day, based on the quantity of CO 2 made available from the pipeline.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Optical analogues to the equatorial Kerr–Newman black hole

Optical analogues to black holes allow the investigation of general relativity in a laboratory setting. Previous works have considered analogues to Schwarzschild black holes in an isotropic coordinate system; the major drawback is that required material properties diverge at the horizon. We present the dielectric permittivity and permeability tensors that exactly reproduce the equatorial Kerr–Newman metric, as well as the gradient-index material that reproduces equatorial Kerr–Newman null geodesics. Importantly, the radial profile of the scalar refractive index is finite along all trajectories except at the point of rotation reversal for counter-rotating geodesics. Construction of these analogues is feasible with available ordinary materials. A finite-difference frequency-domain solver of Maxwell’s equations is used to simulate light trajectories around a variety of Kerr–Newman black holes. For reasonably sized experimental systems, ray tracing confirms that null geodesics can be well-approximated in the lab, even when allowing for imperfect construction and experimental error.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

DART-PFLOTRAN: An ensemble-based data assimilation system for estimating subsurface flow and transport model parameters

Ensemble-based Data Assimilation (EDA), based on the Monte Carlo approach, has been effectively applied to estimate model parameters through inverse modeling in subsurface flow and transport problems. However, implementation of EDA approach involves a complicated workflow that include setting up and executing ensemble forward model simulations, processing observations and model simulation results for parameter updates, and repeat for sequential or iterative EDA. To facilitate the management of such workflow and lower the barriers for adopting EDA-based parameter estimation in subsurface science, we develop a generic software frame-work linking the Data Assimilation Research Testbed (DART) with a massively parallel subsurface FLOw and TRANsport code PFLOTRAN. The new DART-PFLOTRAN leverages both the core data assimilation engines in DART and the computational power afforded by PFLOTRAN. In addition to the standard smoother and filtering options, DART-PFLOTRAN enables an iterative EDA workflow based on the Ensemble Smoother for Multiple Data Assimilation method (ES-MDA) to improve estimation accuracy for nonlinear forward problems. Here, we verify the implementation of ES-MDA in DART-PFLOTRAN using two synthetic cases designed to estimate static permeability and dynamic exchange fluxes across the riverbed, respectively, from continuous temperature measurements made across a depth profile. One-dimensional hydro-thermal simulations are performed in both cases to relate temperature responses with the parameters of interest. In the case of estimating dynamic parameters, we demonstrate the flexibility of DART-PFLOTRAN in automating sequential ES-MDA workflow, which will significantly reduce the time researchers spend on managing complex workflows in similar applications. Both studies yield accurate estimations of the parameters compared to their synthetic truth, while ES-MDA leads to more accurate estimation when a high level of nonlinearity exist between observed responses and unknown parameters. With a code base in Python and Fortran, DART-PFLOTRAN paves the way for applications in large-scale subsurface inverse modeling by automating the complex workflow of sequential ES-MDA that can be executed on various computing platforms.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

Reactive Transport Models with Geomechanics to Mitigate Risks of CO2 Utilization and Storage

Reactivity of carbon dioxide (CO 2 ), rocks and brine is important in a number of practical situations in carbon dioxide sequestration. Injectivity of CO 2 will be affected by near wellbore dissolution or precipitation. Natural fractures or faults containing specific minerals may reactivate leading to induced seismicity. In this project, we first examined if the reactions between CO 2 , brine and rocks affect the nature of the porous medium and properties including petrophysical properties in the timeframe of the injection operations. This was done by carrying out experiments at sequestration conditions (2000 psi for corefloods and 2400 psi for batch experiments, and 600°C) with three different types of rocks – sandstone, limestone and dolomite. Experiments were performed in batch mode and corefloods were conducted over a two-week period. Batch experiments were performed with samples of differing surface area to understand the impact of surface area on overall reaction rates. Toughreact, a reactive transport model was used to interpret and understand the experimental results. The role of iron in dissolution and precipitation reactions was observed to be significant. Iron containing minerals – siderite and ankerite dissolved resulting in changes in porosity and permeability. Corefloods and batch experiments revealed similar patterns. With the right cationic balance, there is a possibility of precipitation of iron bearing carbonates. The results indicate that during injection operations mineralogical changes may lead to injectivity enhancements near the wellbore and petrophysical changes elsewhere in the system. Limestone and dolomite cores showed consistent dissolution at the entrance of the core. The dissolution led to formation of wormholes and interconnected dissolution zones. Results indicate that near wellbore dissolution in these rock-types may lead to rock failure. Micro-CT images of the cores before and after the experiments revealed that an initial high-permeability pathway facilitated the formation of wormholes. The peak cation concentrations and general trends were matched using Toughreact. Batch reactor modeling showed that the geometric factors obtained using powder data that related effective surface area to the BET surface area had to be reduced for fractured samples and cores. This indicates that the available surface area in consolidated samples is lower than that deduced from powder experiments. Field-scale modeling of reactive transport and geomechanics was developed in parallel at Idaho National Laboratory. The model is able to take into account complex chemistry, and consider interactions of natural fractures and faults. Poroelastic geomechanical considerations are also included in the model.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2016

Wellbore Cement Porosity Evolution in Response to Mineral Alteration during CO 2 Flooding

Mineral reactions during CO 2 sequestration will change the pore-size distribution and pore surface characteristics, complicating permeability and storage security predictions. In this study, we report a small/wide angle scattering study of wellbore cement that has been exposed to carbon dioxide for three decades. We have constructed detailed contour maps that describe local porosity distributions and the mineralogy of the sample and relate these quantities to the carbon dioxide reaction front on the cement. We find that the initial bimodal distribution of pores in the cement, 1–2 and 10–20 nm, is affected differently during the course of carbonation reactions. Initial dissolution of cement phases occurs in the 10–20 nm pores and leads to the development of new pore spaces that are eventually sealed by CaCO 3 precipitation, leading to a loss of gel and capillary nanopores, smoother pore surfaces, and reduced porosity. This suggests that during extensive carbonation of wellbore cement, the cement becomes less permeable because of carbonate mineral precipitation within the pore space. Additionally, the loss of gel and capillary nanoporosities will reduce the reactivity of cement with CO 2 due to reactive surface area loss. Finally, this work demonstrates the importance of understanding not only changes in total porosity but also how the distribution of porosity evolves with reaction that affects permeability.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Universal Relationship Between Viscous and Inertial Permeability of Geologic Porous Media

Fluid flow through geologic porous media is reflected by Darcy's law and its inertial and nonlinear extension, the Forchheimer equation. These relationships equate the product of the driving potential gradient and phenomenological coefficients representing momentum resistance and dissipation to flux. From decades of research, the coefficient of viscous permeability (k v ) in Darcy's law is largely predictable, but this is not the case for the coefficient of inertial permeability (k i ) in the Forchheimer equation. Synthesizing results from thousands of laboratory and field flow tests and pore-scale flow model results, we demonstrate that k i can be predicted from k v via the equation k i = 10 10 k v 3/2 across 12 and 20 orders of magnitude in k v and k i , respectively. Since it is related with k i , k v is thus sufficient for predicting flow across viscous to inertial regimes for most geologic porous media.

58 GEOSCIENCES↗

DOE OSTI2019

Measured moisture performance of sealed and insulated attics with permeable insulation in california homes

This study performed detailed thermal and moisture measurements in two homes in Fresno and Clovis, California with sealed attics insulated with glass fiber insulation to determine moisture performance in inland California climates. Data were recorded once every minute for more than a year in each home to allow detailed observation of time-varying moisture and thermal conditions. Moisture measurements were taken at multiple locations in each attic and included wood moisture, surface condensation and surface relative humidity. The results showed that there are strong solar-driven diurnal cycles in attic wood and air moisture content, as well as longer-term seasonal and annual variations. Mold growth was obdserved on the North roof deck in one of the two homes even in this sunny dry climate with attics built in compliance with California building codes that specifically address sealed and insulated attics with vapor-permeable insulation. We investigated moisture stratification in the attic and found that, although there is stratification during sunny daytime hours, this was not when wood moisture content, surface humidity or surface condensation are high. Therefore, other moisture transport mechanisms are acting to increase moisture levels near the roof peak. Finally, our measurements also showed that the attic with the highest measured moisture levels did not have any mold growth, while mold was visually observed in the attic whose measurements indicted little or no concern, indicating that the state-of-the-art in moisture measurement and mold growth prediction in building assemblies may be insufficient.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

Measured Moisture Performance of Sealed and Insulated Attics with Permeable Insulation In California Homes

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2018

Effect of environmental exposure on autogenous self-healing of cracked cement-based materials

Despite abundant literature on self-healing of cement-based materials, there is dearth of information on how environmental exposure affects self-healing. In this study, self-healing of cracks in cement mortar under different environmental exposure was investigated. Pre-cracked mortar specimens were submerged in water, while identical specimens were exposed to cyclic temperature and relative humidity. Change in crack width was examined using optical microscopy. SEM coupled with energy dispersive X-ray analysis was used to identify healing compounds. Mercury Intrusion Porosimetry, water absorption and permeability were employed to assess porosity. X-ray computed tomography was deployed to explore healing of internal cracks. No significant self-healing occurred in specimens exposed to cyclic T and RH. Although crack self-healing was identified in submerged specimens, X-ray μCT demonstrated that it was limited to exposed surface of specimens. This warrants further research to bridge the gap between laboratory findings on self-healing and actual field performance of ageing civil infrastructure.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Enhancing capillary trapping effectiveness through proper time scheduling of injection of supercritical CO 2 in heterogeneous formations

The impact of geologic heterogeneity on capillary trapping of supercritical CO 2 (scCO 2 ) has been recognized and appraised through laboratory experimentation and modeling. However, how different injection strategies can be optimized to improve capillary trapping has not received adequate attention. Here, we present a study based on stochastic analysis to show the impact of injection scheduling on capillary trapping of scCO 2 in heterogeneous geological formations. Improvement of trapping efficiency consists of maximizing the trapped volume within a predefined secure zone of the storage formation with the goal of avoiding plume intersection with potential leakage pathways. Using knowledge acquired from physical experiments in intermediate-scale tanks with controlled heterogeneity, we conduct numerical simulations following a stochastic approach to extend the observed outcomes to a series of equally probable permeability scenarios. Our simulations involve the same combination of surrogate fluids that are used in the experiments to mimic viscosity and density contrasts between scCO 2 and brine. To account for uncertainty of formation heterogeneity, several permeability fields with three different variances and two horizontal correlation lengths are generated. The same volume of scCO 2 is emplaced using four different modes of injection scheduling. Results suggest that injection strategies aimed at enhancing capillary trapping of scCO 2 and increasing the probability of constraining the plume within a secure zone are strictly related to the specific heterogeneity of the reservoir. Finally, this preliminary theoretical finding suggests the potential for maximizing secure capillary trapping through the proper selection of injection schedule to fit the formation heterogeneity.

58 GEOSCIENCES↗

DOE OSTI2019

Deep multiscale model learning

The objective of this article is to design novel multi-layer neural networks for multiscale simulations of flows taking into account the observed fine data and physical modeling concepts. Our approaches use deep learning techniques combined with local multiscale model reduction methodologies to predict flow dynamics. Using reduced-order model concepts is important for constructing robust deep learning architectures since the reduced-order models provide fewer degrees of freedom. We consider flow dynamics in porous media as multi-layer networks in this work. More precisely, the solution (e.g., pressures and saturation) at the time instant n+1 depends on the solution at the time instant n and input parameters, such as permeability fields, forcing terms, and initial conditions. One can regard the solution as a multi-layer network, where each layer, in general, is a nonlinear forward map and the number of layers relates to the internal time steps. We will rely on rigorous model reduction concepts to define unknowns and connections between layers. It is critical to use reduced-order models for this purpose, which will identify the regions of influence and the appropriate number of variables. Furthermore, due to the lack of available observed fine data, the reduced-order model can provide us sufficient inexpensive data as needed. The designed deep neural network will be trained using both coarse simulation data which is obtained from the reduced-order model and observed fine data. We will present the main ingredients of our approach and numerical examples. Numerical results show that using deep learning with data generated from multiscale models as well as available observed fine data, we can obtain an improved forward map which can better approximate the fine scale model.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Effects of natural zeolites on field-scale geologic noble gas transport

Improving predictive models for noble gas transport through natural materials at the field-scale is an essential component of improving US nuclear monitoring capabilities. Several field-scale experiments with a gas transport component have been conducted at the Nevada National Security Site (Non-Proliferation Experiment, Underground Nuclear Explosion Signatures Experiment). However, the models associated with these experiments have not treated zeolite minerals as gas adsorbing phases. This is significant as zeolites are a common alteration mineral with a high abundance at these field sites and are shown here to significantly fractionate noble gases during field-scale transport. This fractionation and associated retardation can complicate gas transport predictions by reducing the signal-to-noise ratio to the detector (e.g. mass spectrometers or radiation detectors) enough to mask the signal or make the data difficult to interpret. Omitting adsorption-related retardation data of noble gases in predictive gas transport models therefore results in systematic errors in model predictions where zeolites are present.Herein is presented noble gas adsorption data collected on zeolitized and non-zeolitized tuff. Experimental results were obtained using a unique piezometric adsorption system designed and built for this study. Data collected were then related to pure-phase mineral analyses conducted on clinoptilolite, mordenite, and quartz. These results quantify the adsorption capacity of materials present in field-scale systems, enabling the modeling of low-permeability rocks as significant sorption reservoirs vital to bulk transport predictions.

58 GEOSCIENCES↗

DOE OSTI2015

Permeation of "Hydromer" Film: An Elastomeric Hydrogen-Capturing Biopolymer.

This report analyzes the permeation resistance of a novel and proprietary polymer coating for hydrogen isotope resistance that was developed by New Mexico State University. Thermal gravimetric analysis and thermal desoprtion spectroscopy show the polymer is stable thermally to approximately 250 deg C. Deuterium gas-driven permeation experiments were conducted at Sandia to explore early evidence (obtained using Brunauer - Emmett - Teller) of the polymer's strong resistance to hydrogen. With a relatively small amount of the polymer in solution (0.15%), a decrease in diffusion by a factor of 2 is observed at 100 and 150 deg C. While there was very little reduction in permeability, the preliminary findings reported here are meant to demonstrate the sensitivity of Sandia's permeation measurements and are intended to motivate the future exploration of thicker barriers with greater polymer coverage.

36 MATERIALS SCIENCE↗

NASA NTRS2015

EM Modeling of Far-Field Radiation Patterns for Antennas on the GMA-TT UAV

To optimize communication with the Generic Modular Aircraft T-Tail (GMA-TT) unmanned aerial vehicle (UAV), electromagnetic (EM) simulations have been performed to predict the performance of two antenna types on the aircraft. Simulated far-field radiation patterns tell the amount of power radiated by the antennas and the aircraft together, taking into account blockage by the aircraft as well as radiation by conducting and dielectric portions of the aircraft. With a knowledge of the polarization and distance of the two communicating antennas, e.g. one on the UAV and one on the ground, and the transmitted signal strength, a calculation may be performed to find the strength of the signal travelling from one antenna to the other and to check that the transmitted signal meets the receiver system requirements for the designated range. In order to do this, the antenna frequency and polarization must be known for each antenna, in addition to its design and location. The permittivity, permeability, and geometry of the UAV components must also be known. The full-wave method of moments solution produces the appropriate dBi radiation pattern in which the received signal strength is calculated relative to that of an isotropic radiator.

Mackenzie, Anne I.↗

DOE OSTI2020

Exploring an ‘ideal hill': how lithology and transport mechanisms affect the possibility of a steady state during weathering and erosion

We introduce a model of chemical reaction within hills to explore how evolving porosity (and by inference, permeability) affects flow pathways and weathering. The model consists of hydrologic and reactive-transport equations that describe alteration of ferrous minerals and feldspar. These reactions were chosen because previous work emphasized that oxygen- and acid-driven weathering affects porosity differently in felsic and mafic rocks. A parameter controlling the order of the fronts is presented. In the absence of erosion, the two reaction fronts move at different velocities and the relative depths depend on geochemical conditions and starting composition. In turn, the fronts and associated changes in porosity drastically affect both the vertical and lateral velocities of water flow. For these instances, estimates of weathering advance rates based on simple models that posit unidirectional constant-velocity advection do not apply. In the model hills, weathering advance rates diminish with time as the Darcy velocities decrease with depth. The system can thus attain a dynamical steady state at any erosion rate where the regolith thickness is constant in time and velocities of both fronts become equal to one another and to the erosion rate. The slower the advection velocities in a system, the faster it attains a steady state. For example, a massive rock with relatively fast-dissolving minerals such as diabase -- where solute transport across the reaction front mainly occurs by diffusion -- can reach a steady state more quickly than granitoid rocks in which advection contributes to solute transport. The attainment of a steady state is controlled by coupling between weathering and hydrologic processes that force water to flow horizontally above reaction fronts where permeability changes rapidly with depth and acts as a partial barrier to fluid flow.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2017

The cyclophilin D/Drp1 axis regulates mitochondrial fission contributing to oxidative stress-induced mitochondrial dysfunctions in SH-SY5Y cells

Oxidative stress plays a central role in the pathogenesis of various neurodegenerative diseases. Increasing evidences have demonstrated that structural abnormalities in mitochondria are involved in oxidative stress related nerve cell damage. And Drp1 plays a critical role in mitochondrial dynamic imbalance insulted by oxidative stress-derived mitochondria. However, the status of mitochondrial fusion and fission pathway and its relationship with mitochondrial properties such as mitochondrial membrane permeability transition pore (mPTP) have not been fully elucidated. Here, we demonstrated for the first time the role of Cyclophilin D (CypD), a crucial component for mPTP formation, in the regulation of mitochondrial dynamics in oxidative stress treated nerve cell. We observed that CypD-mediated phosphorylation of Drp1 and subsequently augmented Drp1 recruitment to mitochondria and shifts mitochondrial dynamics toward excessive fission, which contributes to the mitochondrial structural and functional dysfunctions in oxidative stress-treated nerve cells. CypD depletion or over expression accompanies mitochondrial dynamics/functions recovery or aggravation separately. We also demonstrated first time the link between the CypD to mitochondrial dynamics. Our data offer new insights into the mechanism of mitochondrial dynamics which contribute to the mitochondrial dysfunctions, specifically the role of CypD in Drp1-mediated mitochondrial fission. The protective effect of CsA, or other molecules affecting the function of CypD hold promise as a potential novel therapeutic strategy for governing oxidative stress pathology via mitochondrial pathways. - Highlights: • Demonstrated first time the link between the mPTP to mitochondrial dynamics. • The role of Cyclophilin D in the regulation of Drp1-mediated mitochondrial fission. • CsA as a potential target for governing oxidative stress related neuropathology.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Gas Production from Hot Water Circulation through Hydraulic Fractures in Methane Hydrate-Bearing Sediments: THC-Coupled Simulation of Production Mechanisms

Methane hydrates, widely found in permafrost and deep marine sediments, have great potential as a future energy 6 source. Conventional production schemes perform poorly for challenging hydrate reservoirs with low permeability. We propose an 7 efficient production scheme by combining hydraulic fracturing from horizontal wells and hot water circulation through fractures. A 8 fully coupled thermo-hydro-chemical (THC) model is developed to simulate the key physical processes during gas production from 9 a hydrate reservoir representative of typical geological settings in Shenhu, South China Sea. We found that the gas production 10 process has two distinct stages divided by thermal breakthrough: a relatively short prebreakthrough stage and a postbreakthrough 11 stage yielding stable gas production. Heat advection along and near the hydraulic fracture dominates the prebreakthrough stage, 12 whereas conduction-driven thermal recovery in the volume around fractures dominates the postbreakthrough stage. We identified 13 that the steady-state injection temperature has a strong effect on the performance of the proposed scheme while the fluid mass 14 circulation rate has a moderate impact beyond a threshold. The proposed scheme proves to be efficient and robust over a range of 15 reservoir conditions with respect to initial hydrate saturation and intrinsic permeability, including their spatial heterogeneities, 16 thereby offering a promising solution for challenging reservoir conditions.

gas hydrate, fracture flow, gas production, numeri↗