Search NASA⌕ Search

SEARCH · Search NASA

Results for “permeability models”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An extended trajectory-mechanics approach for calculating two-phase flow paths

A technique originating in quantum dynamics is used to derive a trajectory-based, semi-analytical solution for two-phase flow. The partial differential equation governing the evolution of the aqueous phase is equivalent to a family of ordinary differential equations defined along a path through the porous medium. The trajectories may be found by solving the differential equations directly or by post-processing the output of a numerical solution to the full set of governing equations. The trajectories, which differ from conventional streamlines, are found to bend downward in response to gravitational forces. The curvature is more pronounced as the dip of the porous layer containing the flow increases. Subtle changes in the relative permeability curve can lead to significant variations in the trajectories. The ordinary differential equation for the trajectory provides an expression for the travel time along the path. The expression produces a semi-analytical approximation to the model parameter sensitivities, the partial derivatives of the travel times with respect to changes in the permeability model. The semi-analytical trajectory-based sensitivities generally agree with those computed using a numerical reservoir simulator and a perturbation approach. The sensitivities are useful in tomographic imaging algorithms designed to estimate the spatial variation in permeability within a porous medium using multiphase observations.

58 GEOSCIENCES↗
DOE OSTI2024

Anion Selectivities in Zwitterion Grafted Nanopores: Effect of Zwitterion Architecture

The separation of ions of similar charge is a crucial challenge in many applications, from water treatment to precious metal recovery. Membranes with cross-linked zwitterionic amphiphilic copolymer (ZAC-X) selective layers, which feature self-assembled, zwitterion-lined nanodomains for permeation, offer unique permselectivity between monovalent anions (e.g. $\mathrm{Cl^-}$/$\mathrm{F^-}$). This has motivated studies on the mechanisms of transport and selectivity in this family of materials. In this study, we conducted molecular dynamics simulations of aqueous salt solutions within zwitterion-functionalized nanopores to elucidate the influence of dipole orientation of the ZI ligands on anion diffusivities, partitioning, and permeabilities. Here, our model compares systems with contrasting ZI organization: surface-cation-anion ($\mathrm{S\text{-}ZI^{+}\text{-}ZI^{-}}$, Motif A) and surface-anion-cation ($\mathrm{S\text{-}ZI^{-}\text{-}ZI^{+}}$, Motif B). Our results reveal that Motif A exhibits less pronounced ion pairing due to a spatial separation in the radial profiles of cations and anions. Motif B demonstrates prominent ion pairing for smaller anions owing to their overlap with cation distributions. Further, our potential of mean force profiles reveal that anion partitioning increases with anion size in both ligand motifs, whereas Motif B exhibits significantly higher partitioning selectivity towards larger anions compared to Motif A. Our results for ion diffusivities show that the self-diffusivities of both anions and cations are lower for Motif B compared to Motif A. Such trends in anion partitioning and diffusivities can be explained by differences in the interactions and steric hindrance experienced by the anionic species in Motifs A and B. Finally, our results for anion permselectivity, obtained by combining partitioning and diffusivity, indicate that partitioning trends dominate over diffusivity trends. Consequently, anion permeability increases with anion size, and ligand Motif B yields much higher permselectivity towards larger anions compared to ligand Motif A.

anions↗
DOE OSTI2025

Bacterial microcompartment architectures as biomaterials for conversion of gaseous substrates

Bacterial microcompartments (BMCs) are protein shells encapsulating multiple enzymes of a metabolic pathway. Interpretations of early experiments on carboxysomes led to the narrative that transport of small gases (CO 2 , O 2 ) across the shell membrane is restricted. Since then, this notion has been largely contradicted by studies of engineered shells, although these shell constructs lack important proteins present in the native BMCs, altering the synthetic shells’ topology, surface and mechanical properties. Here, we discuss here an updated model of gas permeability that informs the design of engineered shells for catalysis on gas substrates and outline how nonshell suprastructures of BMC shell proteins could be used in formulating sustainable biomaterials for hydrogen generation via methane pyrolysis and for other greenhouse gas mitigations.

Bacterial microcompartment↗
DOE OSTI2021

Porosity-Permeability Evolution During Simultaneous Mineral Dissolution and Precipitation

Permeability in reactive porous media evolves in complex ways that are hard to predict. Macroscopic empirical equations are often used to estimate permeability evolution but fail to reflect the impact of pore scale reactions on permeability. This work aims to understand the evolution of permeability in systems where mineral dissolution and precipitation occur simultaneously. Pore network models are used for permeability simulation of a Paluxy sandstone using pore and pore-throat size distributions from the analysis of X-ray CT images. Pore and pore-throat radii are increased and decreased to reflect the effects of dissolution and precipitation reactions. Varying spatial distributions of reactions are simulated and the resulting porosity and permeability values compared with commonly used macroscopic equations. It is observed that when dissolution occurs at the inlet and precipitation at the outlet, porosity increases while permeability decreases. Similar results are also observed in the respective opposite scenario where precipitation occurs at the inlet and dissolution at the outlet. When dissolution and precipitation reactions are randomly distributed throughout the network, porosity increases with little change in permeability that is not captured by the empirical equations. When reactions are controlled by pore and pore-throat sizes, such as dissolution occurring in pores and pore-throats of small size and precipitation in pores and pore-throats of larger size, porosity and permeability decrease. When precipitation occurs in pores and pore-throats of smaller size and dissolution in pores and pore-throats of larger size, porosity and permeability increase and is described relatively well by the Verma-Pruess porosity-permeability equation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Effect of solid-electrolyte pellet density on failure of solid-state batteries

Abstract Despite the potentially higher energy density and improved safety of solid-state batteries (SSBs) relative to Li-ion batteries, failure due to Li-filament penetration of the solid electrolyte and subsequent short circuit remains a critical issue. Herein, we show that Li-filament growth is suppressed in solid-electrolyte pellets with a relative density beyond ~95%. Below this threshold value, however, the battery shorts more easily as the density increases due to faster Li-filament growth within the percolating pores in the pellet. The microstructural properties (e.g., pore size, connectivity, porosity, and tortuosity) of $$75\%{{{{{\rm{L}}}}}}{{{{{{\rm{i}}}}}}}_{2}{{{{{\rm{S}}}}}}-25\%{{{{{{\rm{P}}}}}}}_{2}{{{{{{\rm{S}}}}}}}_{5}$$ 75 % L i 2 S − 25 % P 2 S 5 with various relative densities are quantified using focused ion beam–scanning electron microscopy tomography and permeability tests. Furthermore, modeling results provide details on the Li-filament growth inside pores ranging from 0.2 to 2 μm in size. Our findings improve the understanding of the failure modes of SSBs and provide guidelines for the design of dendrite-free SSBs.

25 ENERGY STORAGE↗
DOE OSTI2021

Modelling the development of capillary pressure in freshly 3D-printed concrete elements

3D concrete printing is a promising technology recently developed to automate construction. Since no formwork is used in this technology to support and protect fresh concrete, there are two aspects which considerably accelerate the development of capillary pressure in 3D-printed concrete in comparison to conventionally placed concrete: i) high stiffness of 3D-printed needed to provide sufficient buildability, and ii) very early and fast evaporation of pore water. Accelerated development of capillary pressure may lead to severe plastic shrinkage cracking in 3D-printed elements and, hence, need to be mitigated. This investigation aims at providing a poromechanical model for capillary pressure development in 3D-printed elements. To simulate the development of capillary pressure and plastic shrinkage, environmental factors, material properties, and element geometry need to be considered as a whole. The model inputs – coefficient of permeability, static bulk modulus, air entry pressure and chemical shrinkage rate – were determined experimentally. The model was validated for two fine-grained concretes. Both 3D-printed materials yielded faster capillary pressure increase in comparison to cast concrete, while partial substitution of cement with silica fume further accelerated the capillary pressure development. Furthermore, due to the lower permeability of the mixture containing silica fume, the gradient of capillary pressure between 3D-printed layers increased, as did the gradient of plastic shrinkage.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Relative permeability as a stationary process: Energy fluctuations in immiscible displacement

Relative permeability is commonly used to model immiscible fluid flow through porous materials. In this work, we derive the relative permeability relationship from conservation of energy, assuming that the system to be non-ergodic at large length scales and relying on averaging in both space and time to homogenize the behavior. Explicit criteria are obtained to define stationary conditions: (1) there can be no net change for extensive measures of the system state over the time averaging interval; (2) the net energy inputs into the system are zero, meaning that the net rate of work done on the system must balance with the heat removed; and (3) there is no net work performed due to the contribution of internal energy fluctuations. Results are then evaluated based on direct numerical simulation. Dynamic connectivity is observed during steady-state flow, which is quantitatively assessed based the Euler characteristic. We show that even during steady-state flow at low capillary number ([Formula: see text]), typical flow processes will explore multiple connectivity states. The residence time for each connectivity state is captured based on the time-and-space average. The distribution for energy fluctuations is shown to be multi-modal and non-Gaussian when terms are considered independently. However, we demonstrate that their sum is zero. Given an appropriate choice of the thermodynamic driving force, we show that the conventional relative permeability relationship is sufficient to model the energy dissipation in systems with complex pore-scale dynamics that routinely alter the structure of fluid connected pathways.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Microstructure and water absorption of ancient concrete from Pompeii: An integrated synchrotron microtomography and neutron radiography characterization

There is renewed interest in using advanced techniques to characterize ancient Roman concrete due to its exceptional durability and low-carbon footprint. In the present work, samples were drilled from the “Hospitium” in Pompeii and were analyzed by synchrotron microtomography (μCT) and neutron radiography to study how the microstructure, including the presence of induced cracks, affects their water adsorption. The water distribution and absorptivity were quantified by neutron radiography. The 3D crack propagation, pore size distribution and orientation, tortuosity, and connectivity were analyzed from μCT results using advanced imaging methods. Porosity was also measured by mercury intrusion porosimetry (MIP) as a reference. Ductile fracture patterns were observed once cracks were introduced. Compared to Portland cement mortar/concrete, the Pompeii samples had relatively high porosity, low connectivity, and a similar coefficient of capillary penetration. In addition, permeability was predicted from models based on percolation theory and pore structure data to evaluate the fluid transport properties. Understanding the microstructure of ancient Pompeii concrete is important because it could inspire the development of modern concrete with high durability.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Discovering Blind Geothermal Systems in the Great Basin Region: An Integrated Geologic and Geophysical Approach for Establishing Geothermal Play Fairways: All Phases

Most geothermal resources in the Great Basin region of the western USA are blind, and thus the discovery of new commercial-grade systems requires synthesis of favorable characteristics for geothermal activity. The geothermal play fairway concept involves integration of multiple parameters indicative of geothermal activity to identify promising areas for new development. This project integrated multiple datasets to apply the play fairway concept and assess geothermal potential in a large region of the Great Basin in Nevada. It is therefore referred to as the Nevada play fairway project. This project was a strong collaborative effort between several organizations, led by the Nevada Bureau of Mines and Geology at the University of Nevada, Reno, but with key support from the U.S. Geological Survey, ATLAS Geosciences, Inc,, Hi-Q Geophysical, Inc., Lawrence Berkeley National Laboratory, Utah Geological Survey, and Innovative Geothermal Ltd. In Budget Period 1 of this project, available data for nine geologic, geochemical, and geophysical parameters were initially synthesized to produce a new detailed geothermal potential map of 96,000 km2 from west-central to eastern Nevada. These parameters were grouped into subsets and individually weighted to delineate rankings for local permeability, intermediate permeability, regional permeability, and thermal potential, which collectively defined geothermal play fairways (i.e., most likely locations for significant geothermal fluid flow). This initial work was aimed at reducing the risks in regional exploration and therefore facilitating discovery of new commercial-grade systems in blind settings, as well as in areas with surface expressions of geothermal activity. Budget Period 2 of the project involved detailed analysis of some of the most promising areas identified in Phase 1. Twenty-four highly prospective areas, including both known undeveloped systems and previously undiscovered potential blind systems, were identified for further analysis. After reconnaissance of these areas, five of the most promising sites were selected for detailed studies. Multiple techniques were employed in the detailed studies, including geologic mapping, shallow temperature surveys, gravity surveys, Lidar, geochemical studies, seismic reflection analysis, and 3D modeling. The goal of the detailed studies was to identify specific areas with the highest likelihood for high permeability and thermal fluids, such that drill sites could be targeted. Three main sets of predictive maps were generated for each detailed study area: 1) play fairway maps, 2) play fairway error maps, and 3) direct evidence maps. Local- and intermediate-scale permeability models were revised to reflect results of the detailed geologic, geophysical, and geochemical analyses. Budget Period 3 of the project involved more detailed geophysical analyses and temperature-gradient (TG) drilling in southeastern Gabbs Valley and northern Granite Springs Valley, deemed the two most promising sites, with the goal of providing preliminary validation of the play fairway methodology. In southeastern Gabbs Valley, the collocation of a favorable structural setting (displacement transfer zone and fault intersections), Quaternary faults, intersecting and terminating gravity gradients, magnetic low, shallow (2 m) temperature anomaly, low resistivity anomaly, and promising geothermometry from nearby water wells provided evidence for a blind system. Drilling of six TG holes defines an apparent geothermal system at this locality with temperatures as high as 124°C at 152 m. This system is blind, with no surface hot springs, fumaroles, or paleo-geothermal deposits. For northern Granite Springs Valley, a favorable structural setting (termination of a major Quaternary normal fault), terminating gravity gradient, magnetic gradient, newly discovered sinter deposits, nearby warm water wells, previously drilled TG holes in the vicinity, and promising geothermometry suggest a hidden system. Drilling of six new TG holes yields temperatures of ~96°C at ~250 m, suggesting the presence of a geothermal system. Major lessons learned in the course of this project include: 1) initially identified sites commonly include multiple favorable structural settings at a finer scale; 2) promising sites in Cenozoic basins cannot be recognized without detailed geophysical surveys; and 3) play fairway analysis should be refined as the exploration program vectors into the most promising sites and finer-scale data are acquired. In addition to producing copious amounts of data, this project resulted in 16 published papers, 10 abstracts, more than 40 presentations across the U.S. and abroad (including several keynote addresses), 2 Masters theses, and 7 media reports.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Induced microearthquakes predict permeability creation in the brittle crust

Predicting the evolution of permeability accurately during stimulation at the reservoir scale and at the resolution of individual fractures is essential to characterize the fluid transport and the reactive/heat-transfer characteristics of reservoirs where stress exerts significant control. Here, we develop a hybrid machine learning (ML) model to visualize in situ permeability evolution for an intermediate-scale (~10 m) hydraulic stimulation experiment. This model includes an ML model that was trained using the well history of flow rate and wellhead pressure and MEQ data from the first three stimulation episodes to predict average permeability from the statistical features of the MEQs alone for later episodes. Moreover, a physics-inspired model is integrated to estimate in situ fracture permeability spatially. This method relates fracture permeability to fracture dilation and scales dilation to the equivalent MEQ magnitude, according to laboratory observations. The seismic data are then applied to define incremental changes in permeability in both space and time. Our results confirm the excellent agreement between the ground truth and model-predicted permeability evolution. The resulting permeability map defines and quantifies flow paths in the reservoir with the averaged permeability comparing favorably with the ground truth of permeability.

induced seismicity↗
DOE OSTI2023

Morrow1 Model - No Capillary Pressure

Farnsworth Unit (FWU) CO2EOR Eclipse compositional model: Morrow1 model uses the Morrow1 relative permeability curve with no capillary pressure assigned.

capillary pressure↗
DOE OSTI2023

Morrow1 Model - with Capillary Pressure

Farnsworth Unit (FWU) CO2EOR Eclipse compositional model: Morrow1 model uses the Morrow1 relative permeability curve with capillary pressure assigned homogeneously

Capillary Pressure↗
DOE OSTI2020

Relating land surface, water table, and weathering fronts with a conceptual valve model for headwater catchments

Knowing little about how porosity and permeability are distributed at depth, we commonly develop models of groundwater by treating the subsurface as a homogeneous black box even though porosity and permeability vary with depth. One reason for this depth variation is that infiltrating meteoric water reacts with minerals to affect porosity in localized zones called reaction fronts. We are beginning to learn to map and model these fronts beneath headwater catchments and show how they are distributed. The subsurface landscapes defined by these fronts lie subparallel to the earth-air interface but with lower relief. They can be situated above, below, or at the water table. These subsurface landscapes of reaction are important because porosity developed from weathering can control subsurface water storage. In addition, porosity often changes at the weathering fronts, and when this affects permeability significantly, the front can act like a valve that re-orients water flowing through the subsurface. We explore controls on the positions of reaction fronts under headwater landscapes by accounting for the timescales of erosion, chemical equilibration, and solute transport. One strong control on the landscape of subsurface reaction is the land surface geometry, which is in turn a function of the erosion rate. In addition, the reaction fronts, like the water table, are strongly affected by the lithology and water infiltration rate. We hypothesize that relationships among the land surface, reaction fronts, and the water table are controlled by feedbacks that can push landscapes toward an “ideal hill”. In this steady state, reaction-front valves partition water volumes into shallow and deep flowpaths. Furthermore, these flows dissolve low- and high-solubility minerals, respectively, allowing their reaction fronts to advance at the erosion rate. This conceptualization could inform better models of subsurface porosity and permeability, replacing the black box.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2021

Estimating Watershed Subsurface Permeability From Stream Discharge Data Using Deep Neural Networks

Subsurface permeability is a key parameter in watershed models that controls the contribution from the subsurface flow to stream flows. Since the permeability is difficult and expensive to measure directly at the spatial extent and resolution required by fully distributed watershed models, estimation through inverse modeling has had a long history in subsurface hydrology. The wide availability of stream surface flow data, compared to groundwater monitoring data, provides a new data source to infer soil and geologic properties using integrated surface and subsurface hydrologic models. As most of the existing methods have shown difficulty in dealing with highly nonlinear inverse problems, we explore the use of deep neural networks for inversion owing to their successes in mapping complex, highly nonlinear relationships. We train various deep neural network (DNN) models with different architectures to predict subsurface permeability from stream discharge hydrograph at the watershed outlet. The training data are obtained from ensemble simulations of hydrographs corresponding to an permeability ensemble using a fully-distributed, integrated surface-subsurface hydrologic model. The trained model is then applied to estimate the permeability of the real watershed using its observed hydrograph at the outlet. Our study demonstrates that the permeabilities of the soil and geologic facies that make significant contributions to the outlet discharge can be more accurately estimated from the discharge data. Their estimations are also more robust with observation errors. Compared to the traditional ensemble smoother method, DNNs show stronger performance in capturing the nonlinear relationship between permeability and stream hydrograph to accurately estimate permeability. Our study sheds new light on the value of the emerging deep learning methods in assisting integrated watershed modeling by improving parameter estimation, which will eventually reduce the uncertainty in predictive watershed models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Flow-Driven Branching in a Frangible Porous Medium

Channel formation and branching is widely seen in physical systems where movement of fluid through a porous structure causes the spatiotemporal evolution of the medium. We provide a simple theoretical framework that embodies this feedback mechanism in a multiphase model for flow through a frangible porous medium with a dynamic permeability. Numerical simulations of the model show the emergence of branched networks whose topology is determined by the geometry of external flow forcing. Furthermore, this allows us to delineate the conditions under which splitting and/or coalescing branched network formation is favored, with potential implications for both understanding and controlling branching in soft frangible media.

36 MATERIALS SCIENCE↗