Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geologic Modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Review: Pre-Darcy flows in low-permeability porous media

The widely used Darcy’s law specifies a linear relation between the Darcy velocity of fluid flow and the pressure gradient that drives the flow. However, studies have shown that Darcy velocity can exhibit a nonlinear dependence on the pressure gradient in low-permeability porous media such as clay and shale when the pressure gradient is adequately low. This phenomenon is referred to as low-velocity non-Darcian flow or pre-Darcy flow. This paper provides a comprehensive review of the theories, experimental data, and modeling methods for pre-Darcy flow in low-permeability porous media. The review begins by outlining the fundamental mechanisms underlying pre-Darcy flow that regulate the unique characteristics such as nonlinear dependence of the Darcy velocity on the pressure gradient and its relevance to fluid–rock interactions. The review then proceeds to present a thorough compilation of experimental investigations performed in various low-permeability geomaterials including tight sandstones, shales, and clays. Next, empirical and theoretical models and simulation methods that have been developed to fit and interpret experimental data are reviewed. Finally, the review underscores the challenges encountered in conducting and interpreting pre-Darcy flow experiments and suggests future research directions. By analyzing previous experimental investigations, this review aims to offer a valuable resource for researchers and practitioners seeking to enhance their understanding of fluid dynamics in low-permeability geomaterials. This provides insights into the application of pre-Darcy flow in numerous natural and engineered processes such as shale oil and gas recovery, contaminant transport in low-permeability aquifers, and geological disposal of nuclear waste.

58 GEOSCIENCES↗
DOE OSTI2020

Shear strength and permeability of the cement-casing interface

Here, the shear strength and hydraulic permeability of the interface between well cement and casing was investigated using a triaxial direct shear apparatus. For the first time, these experiments provide measurements under controlled stress conditions with fluid flow measurements along the interface. The low cohesion (1.1 ± 1.1 MPa) and the high friction angle (43.4 ± 2.0°) indicates that the shear strength of the interface is provided by friction. This implies that the state of stress of the cement is critical to well integrity. The hydraulic aperture of the undamaged cement-steel samples was 6.8 ± 1.0 microns. Shear damage to the interface caused a decrease (-20 %) in hydraulic aperture for samples aged up to 1 month, and an increase (+300 %) for samples cured for two years. We performed numerical simulations to estimate the leakage potential from a carbon storage operation. This model predicts negligible leakage amounts (47 tonnes) in a shear-damaged well for the modeled injection of ~1.26 million tonnes of CO 2 . Thus, our measurements indicate that the cement-casing interface is not a significant leakage pathway in its intact or damaged state, and that shear-driven failure scenarios for this interface are not a significant risk to CO 2 storage security.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Sand Migration Simulation during Gas Production from Gas Hydrate Reservoir at Kuparuk 7–11–12 site in the Prodhoe Bay Unit, Alaska

Uncontrolled sand production impedes continuous gas production from a hydrate reservoir as observed in the past field-scale gas production tests. Sand mobilization is strongly linked with sediment deformation and high pressure gradient. Sand production occurs when the mobilized sands reach the well. Throughout gas production from a hydrate reservoir, because of hydrate dissociation, deformation and pressure gradient change both in time and space. Therefore, it is necessary to consider the entire process to identify where sand is likely mobilized and how much mobilized sand could reach the well. Here, this study utilizes a coupled thermal, hydrological, chemical, and mechanical (thermo-hydro-chemo-mechanical) sand migration model to simulate a year long gas production from the hydrate reservoir at the Kuparuk 7–11–12 site in the Prodhoe Bay Unit, Alaska. It is found that sand would mainly come from the lower portion of the production zone where faster hydrate dissociation occurs. The relatively faster hydrate dissociation coupled with the fact that hydrate-bearing sediments and hydrate-free sediments have similar stiffness but different strengths causes complex stress transfer, which results in excessive shear deformation when the upper portion of the production zone starts to dissociate. This is evident not only near the well but also away from the well, leading to a large amount of sand mobilization. Another focus of the modeling study is evaluation of the effect of sand migration on gas production. The comparison of two extreme cases–with sandscreen and without sandscreen–suggests that migrated and settled solids around a sand control device can prevent pressure drop across the near-wellbore zone and lead to reduction of the gas production rate.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Earthquake-triggered submarine canyon flushing transfers young terrestrial and marine organic carbon into the deep sea

Submarine canyons transfer substantial amounts of sediment and organic carbon (OC) into the deep ocean, nourishing deep-sea ecosystems and contributing to the global carbon cycle through OC burial and sequestration. Tracking lateral OC transport through submarine canyon systems is challenged by the deep-ocean setting, difficulties with constraining episodic depositional events, and the need to assess the composition and age of marine and terrestrial organic matter. We apply innovative parallel ramped pyrolysis oxidation accelerator mass spectrometry and pyrolysis-gas chromatography-mass spectrometry with isotope analyses to track OC age and sources in the 2016 Kaikoura earthquake-triggered, canyon-flushing event that deposited along >1300 km of a submarine canyon-channel system, offshore Aotearoa New Zealand. Specifically, these techniques allow us to determine the ages, sources, and partitioning of OC within the Kaikoura turbidite deposit and test hypotheses of how submarine canyon systems contribute to lateral OC flux and burial. Our results show that, despite considerable canyon floor erosion, substantial amounts of young OC were flushed into the deep sea, with relatively little (~2%) pre-Holocene OC contributions. Even without a direct connection between rivers and submarine canyons, most (~55%) of the OC in the Kaikoura event bed is from terrestrial sources. However, the deposit also contains substantial amounts (~22%) of marine-derived OC and ~23% of the material is of unassignable origin. Particle sorting imparts variability on the age and composition of OC within turbidite deposits and along the turbidity current flow path. Terrestrial-derived OC is preferentially older than marine-derived OC and concentrated in coarser particle sizes found more commonly at the deposit base and in proximal settings. Young, marine-derived OC is concentrated at the surface of the deposits and tends to be enriched in finer particle sizes. Such OC partitioning in turbidites supports the relevance of depositional models for predicting and quantifying distribution of OC in deep-sea deposits. Earthquake-triggered, canyon flushing events and resulting turbidites enhance OC burial efficiency and can sequester OC effectively, contributing an important carbon sink to the sedimentary carbon cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

On the Scaling and Growth Limit of Fluvial Dunes

Abstract Bedforms emerge in a variety of shapes and sizes when a granular erodible surface is subject to a strong enough shear flow, as observed in topographic data from submarine canyons, rivers, deserts, and planetary bodies. The two salient features of bedforms are the ability to collectively transport particles, and to generate form drag thus increasing flow resistance. These two mechanisms are in competition and contribute to force bedforms of increasing size migrating more slowly. In a dedicated large‐scale open channel flow facility, hierarchies of fluvial bedforms were generated and measured in equilibrium conditions. The corresponding scale‐dependent migration velocity and mass flux contributions were quantified in the wave number and frequency domains. Experimental results are paired with a validated set of theoretical models to demonstrate that ripples or dunes reach an equilibrium state when drag partitioning ensures high enough frictional drag to sustain the bedload transport of sediment, and low enough form drag to enable the migration of the largest bedform size. This mechanism is inferred to constrain the growth of bedforms when sediment supply is not the limiting factor.

58 GEOSCIENCES↗
DOE OSTI2021

Local summer insolation and greenhouse gas forcing drove warming and glacier retreat in New Zealand during the Holocene

Geological climate archives from the Holocene Epoch provide baseline information concerning natural climate variability. Temperate mountain glacier extent is sensitive to summer air temperature, thus geological records of past glacier length changes are a useful proxy for this climatic variable. Here we present a new cosmogenic 10 Be chronology of glacier length changes at Dart Glacier in the Southern Alps, New Zealand. Prominent moraines deposited 321 ± 44 yr ago (n = 11) and 7.8 ± 0.3 ka (n = 5) show glaciers during the Little Ice Age were less extensive than during the early Holocene. This pattern of net Holocene glacier retreat is consistent with emerging data from other catchments in New Zealand and across the southern mid-latitudes. Using the physical framework of a transient global climate model simulation, we suggest that cool summers in the early Holocene were promoted by the local summer insolation minimum, together with low atmospheric greenhouse gas concentrations, causing an early Holocene austral glacial maximum. An insolation-driven reduction in seasonality at southern mid-latitudes may reconcile differences between early Holocene temperature reconstructions where climate proxies have different seasonal sensitivities. Furthermore, we suggest that rising greenhouse gas concentrations after 7 ka caused regional-scale glacier retreat and appear to be the dominant driver of multi-millennial summer temperature trends in the southern mid-latitudes during the present interglacial.

58 GEOSCIENCES↗
DOE OSTI2020

Calculation of the Vertical Strata Load of Utility Tunnel Crossing Ground Fissure Zone

A ground fissure is a geological disaster in which the vertical dislocation of strata causes surface rupture. Ground fissures can cause extreme harm to the surface and underground buildings. Ground fissure activity can result in different settlement on the two sides of the strata, which will generate additional stress (pressure) that differs from the stress of the general stratum on underground structures across the ground fissure zone. It is essential to assess the effective stress of strata in the design of underground engineering structures across a ground fissure zone. The Xi’an ground fissure through a utility tunnel was focus of the research, and a physical model and data for oblique crossing of the 45° ground fissure were analyzed. A model of the utility tunnel structure was established, including the surrounding soil load as an active ground fissure environment. This model was used to calculate the vertical formation pressure of the overlying soil on the utility tunnel. A method to calculate the overlying load on the utility tunnel caused by ground fissure activity was proposed and compared with the calculation based on the A. Marston principle. The results showed that the ground fissure load calculation method based on the strata-holding effect can effectively calculate the earth pressure of the surrounding soil layer of the utility tunnel in the cross-ground fissure section. The results of this work provide guidance and reference value for the design of a utility tunnel in an area with the potential for a ground fissure.

Zhang, Dan↗
DOE OSTI2022

Exploring for Superhot Geothermal Targets in Magmatic Settings: Developing a Methodology

This paper presents preliminary results from a subset of work carried out as part of a multinational research project entitled DErisking Exploration for multiple geothermal Plays in magmatic ENvironments (DEEPEN). One objective of DEEPEN is to develop a customized approach to exploration for superhot geothermal plays in magmatic systems. This paper summarizes key geologic components, risk factors, and exploration methods for geothermal plays in magmatic settings based on a review and comparative analysis of international training sites. As part of a Play Fairway Analysis (PFA) approach to exploring for multiple play types in a single magmatic system, training data were compiled and weights assigned to various evidence layers. Two different approaches for weighting exploration datasets are described in this paper - one based on expert opinions and the other using statistical learning. Weights produced by both approaches will be input into a 3D PFA workflow that combines multiple exploration datasets to generate 3D geothermal favorability models, which will be applied to two international demonstration sites.

GEOTHERMAL ENERGY↗
DOE OSTI2025

A review of Geological Thermal Energy Storage for seasonal, grid-scale dispatching

Energy storage is essential for the decarbonization of the U.S. energy grid, especially with the increasing deployment of variable renewable energy sources like solar and wind. Geological thermal energy storage (GeoTES) has emerged as a promising long duration, grid scale solution, providing stability and security through flexible operations and valuable grid services. GeoTES utilizes subsurface reservoirs to store thermal energy for power generation and direct-use heating and cooling. This approach significantly enhances the use of low-temperature reservoirs, which would otherwise be unsuitable for geothermal power plants. It also aligns well with depleted oil and gas reservoirs, concentrating solar power, non-flexible renewables (photovoltaic and wind), and geothermal-related power cycles. Given the favorable marginal costs of GeoTES as storage duration increases, it becomes particularly competitive for seasonal, grid-scale dispatch, where few technologies are viable. This paper provides a comprehensive review of GeoTES systems and the research underpinning itsr development. This analysis begins by defining and categorizing the unique characteristics of thermal energy storage techniques, setting GeoTES apart from other technologies. The various components, configurations, subsurface characteristics, and modeling efforts that guide GeoTES development are then explored. Finally, challenges in GeoTES research, development, and deployment are discussed, along with mitigation strategies and lessons from related technologies. Beyond their economic benefits, GeoTES systems support grid resilience and decarbonize industrial processes. Their scalability, broad distribution, seasonal storage potential, and flexible dispatch capacity make GeoTES a valuable tool for expanding renewable energy deployment and addressing climate change.

15 - GEOTHERMAL ENERGY↗
DOE OSTI2023

Derivation of heterogeneous material distributions and their sensitivity to HM-coupled two-phase flow models exemplified with the LASGIT experiment

Abstract Advective gas transport in bentonite, a possible buffer material in repositories for radioactive materials, is difficult to simulate in numerical continuum models, partly due to the complicated microstructure of bentonite. To generate reliable models of repositories nevertheless, spatially distributed heterogeneous material properties can be used to allow localization of gas flow. In this study, a pore-size-dependent stochastic approach of the gas entry pressure is derived from Mercury Intrusion Porosimetry, which is used to replicate measurements from the LASGIT experiment. In addition, three benchmark tests are simulated to investigate the dependence of heterogeneous distributions of material properties on the mesh discretization, the temporal dependence, and the coupling between the processes influenced by the heterogeneous parameters. The numerical modeling results of the LASGIT experiment show that the onset of gas flow into the system and the subsequent increase in pressure and stress can be well reproduced using heterogeneous distributions. Compared to a model with homogeneous material properties, heterogeneous distributions may allow the generation of dilatancy-controlled microfractures—an important feature with regard to the advective gas flow in bentonites. However, it can be observed that the heterogeneous distributions in LASGIT are less significant, as technical gaps or differences in material types could have a greater impact.

Environmental Sciences & Ecology↗
DOE OSTI2022

Paw-Net: Stacking ensemble deep learning for segmenting scanning electron microscopy images of fine-grained shale samples

Segmentation of scanning electron microscopy (SEM) images is critical yet time-consuming for geological analyses, as it needs to differentiate the boundaries for different mineral objects to facilitate subsequent analyses, such as porosity calculation. Recently, various machine learning methods, especially convolutional neural networks (CNNs), have been explored to segment SEM images of fine-grained shale samples. However, we found that general CNNs do not yield optimal performance due to insufficient training data and imbalanced objects in SEM images. This work has revised the U-Net architecture, a popular approach for biomedical image analyses, by incorporating a loss function that addresses the imbalance issue. Furthermore, we used the ensemble learning method to train multiple models and combined the results to improve the overall performance of segmentation. We prepared 2162 sub-images from raw SEM images in our experiments and divided them into training, validation, and testing datasets. The overall results show that our method improves the average Intersection over Union (IOU) of mineral objects from 0.49 to 0.58, compared to the original U-Net model. Our method can clearly distinguish each object from others with boundaries, even in highly imbalanced images. Training our models takes less than three minutes using a single GPU, while manual labeling can take up to three hours for each image. Furthermore, the method helps geoscientists gain insights quickly and effectively by building neural network models from a small dataset of SEM images.

58 GEOSCIENCES↗
DOE OSTI2020

Integrated hydrogeophysical modelling and data assimilation for geoelectrical leak detection

Time-lapse electrical resistivity tomography (ERT) measurements provide indirect observations of hydrological processes in the Earth's shallow subsurface at high spatial and temporal resolution. ERT has been used in the past decades to detect leaks and monitor the evolution of associated contaminant plumes. Specifically, inverted resistivity images allow visualization of the dynamic changes in the structure of the plume. However, existing methods do not allow the direct estimation of leak parameters (e.g. leak rate, location, etc.) and their uncertainties. We propose an ensemble-based data assimilation framework that evaluates proposed hydrological models against observed time-lapse ERT measurements without directly inverting for the resistivities. Each proposed hydrological model is run through the parallel coupled hydro-geophysical simulation code PFLOTRAN-E4D to obtain simulated ERT measurements. The ensemble of model proposals is then updated using an iterative ensemble smoother. In this paper, we demonstrate the proposed framework on synthetic and field ERT data from controlled tracer injection experiments. Our results show that the approach allows joint identification of contaminant source location, initial release time, and solute loading from the cross-borehole time-lapse ERT data, alongside with an assessment of uncertainties in these estimates. We demonstrate a reduction in site-wide uncertainty by comparing the prior and posterior plume mass discharges at a selected image plane. This framework is particularly attractive to sites that have previously undergone extensive geological investigation (e.g., nuclear sites). It is well suited to complement ERT imaging and we discuss practical issues in its application to field problems.

58 GEOSCIENCES↗
DOE OSTI2020

Modeling Holocene Peatland Carbon Accumulation in North America

Peatlands are a large carbon reservoir. Yet the quantification of their carbon stock still has a large uncertainty due to lacking observational data and well-tested peatland biogeochemistry models. Here, a process-based peatland model was calibrated using long-term peat carbon accumulation data at multiple sites in North America. The model was then applied to quantify the peat carbon accumulation rates and stocks within North America over the last 12,000 years. We estimated that 85–174 Pg carbon was accumulated in North American peatlands over the study period including 0.37–0.76 Pg carbon in subtropical peatlands. During the period from 10,000 to 8,000 years ago, the warmer and wetter conditions might have played an important role in stimulating peat carbon accumulation by enhancing plant photosynthesis. Enhanced peat decomposition due to warming slowed the carbon accumulation through the rest of the Holocene. While recent modeling studies indicate that the northern peatlands will continue to act as a carbon sink in this century, our studies suggest that future enhanced peat decomposition accompanied by peatland areal changes induced by permafrost degradation and other disturbances shall confound the sink and source analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Low 13 C- 13 C abundances in abiotic ethane

Distinguishing biotic compounds from abiotic ones is important in resource geology, biogeochemistry, and the search for life in the universe. Stable isotopes have traditionally been used to discriminate the origins of organic materials, with particular focus on hydrocarbons. However, despite extensive efforts, unequivocal distinction of abiotic hydrocarbons remains challenging. Recent development of clumped-isotope analysis provides more robust information because it is independent of the stable isotopic composition of the starting material. Here, we report data from a 13 C- 13 C clumped-isotope analysis of ethane and demonstrate that the abiotically-synthesized ethane shows distinctively low 13 C- 13 C abundances compared to thermogenic ethane. A collision frequency model predicts the observed low 13 C- 13 C abundances (anti-clumping) in ethane produced from methyl radical recombination. In contrast, thermogenic ethane presumably exhibits near stochastic 13 C- 13 C distribution inherited from the biological precursor, which undergoes C-C bond cleavage/recombination during metabolism. Further, we find an exceptionally high 13 C- 13 C signature in ethane remaining after microbial oxidation. In summary, the approach distinguishes between thermogenic, microbially altered, and abiotic hydrocarbons. The 13 C- 13 C signature can provide an important step forward for discrimination of the origin of organic molecules on Earth and in extra-terrestrial environments.

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

Grain boundary facilitated dissolution of nanocrystalline NpO 2 (s) from legacy waste processing

Dissolution of actinide dioxides, including neptunium dioxide (NpO 2 (s)), is paramount for the prediction of the environmental fate of nuclear materials. Quantifying dissolution rates, as well as understanding qualitative dissolution mechanisms, informs performance assessment for geologic disposal of spent nuclear fuel and management of legacy radioactive waste. The aim of this research was to measure the dissolution rate of nanocrystalline NpO 2 (s), produced through legacy nuclear waste processing, under oxidizing conditions, as well as to characterize surface alteration to the material. The solid phase was characterized using electron microscopy techniques (SEM/STEM) and X-ray photoelectron spectroscopy (XPS), indicating preferential dissolution of Np-hydroxide contained in the grain boundaries of NpO 2 (s) and fragmentation of grains from the matrix. The oxidative dissolution was monitored over 40 weeks, yielding a two-step kinetic dissolution model involving hydration of NpO 2 (s) and subsequent oxidation and dissolution of the hydroxide phase. Here, the proposed dissolution models for nanocrystalline NpO 2 (s) suggest that microstructural features such as grain boundaries are key factors affecting dissolution, including release of colloidal particles, and ultimately, environmental fate and transport of nuclear materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Glacial geology of the Hudson Mountains, Amundsen Sea sector, West Antarctica

The Hudson Mountains are situated in the eastern Amundsen Sea sector of the West Antarctic Ice Sheet, adjacent to Pine Island Glacier. They form a volcanic field of 17 stratovolcanoes and parasitic vents, preserved as nunataks. Two former tributaries of Pine Island Glacier (Larter and Lucchitta glaciers) flow through the mountains. Here we present a detailed study of the glacial geology of the area. We describe field observations and measurements of geomorphological features from 15 of the nunataks, meltwater ponds found on the surface of three nunataks and supraglacial features (ice dolines) from two sites near the present grounding line. Together these provide constraints on the past ice sheet extent, flow pathways and thermal regime, and enhance our understanding of the present hydrological regime – all of which are important as context for the observed modern ice sheet behaviour. We find evidence suggesting that all nunataks in the Hudson Mountains were covered by ice during the Last Glacial Maximum (defined here as 26.5-19 ka) and have since deglaciated. Faceted and polished erratic cobbles and boulders of exotic lithologies (syenites, alkali granites, granites, granodiorites, tonalites and gabbros) are numerous and perched on nunatak surfaces. A marked difference between the dominant erratic lithologies on nunataks adjacent to Pine Island Glacier (granite) and Lucchitta Glacier (granodiorite-tonalite) indicates that the ice sheet was transporting clasts from at least two distinct upstream source regions. The similarity in degree of weathering suggests, however, that all the erratics were transported by one phase of (warm-based) glaciation; their presence on or close to the summits of all except one nunatak indicates that the ice sheet during that time was at least 700 m thicker than present. These results are consistent with ice sheet model simulations which suggest that all nunataks in the Hudson Mountains were completely submerged by the Last Glacial Maximum ice sheet.

58 GEOSCIENCES↗
DOE OSTI2021

Electrical Properties of Carbon Dioxide Hydrate: Implications for Monitoring CO 2 in the Gas Hydrate Stability Zone

CO 2 and CH 4 clathrate hydrates are of keen interest for energy and carbon cycle considerations. While both typically form on Earth as cubic structure I (sI), we find that pure CO 2 hydrate exhibits over an order of magnitude higher electrical conductivity ( σ ) than pure CH 4 hydrate at geologically relevant temperatures. In this study, the conductivity was obtained from frequency-dependent impedance ( Z ) measurements made on polycrystalline CO 2 hydrate (CO 2 ·6.0 ± 0.2H 2 O by methods here) with 25% gas-filled porosity, compared with CH 4 hydrate (CH 4 ·5.9H 2 O) formed and measured in the same apparatus and exhibiting closely matching grain characteristics. The conductivity of CO 2 hydrate is 6.5 × 10 –4 S/m at 273K with an activation energy ( E a ) of 46.5 kJ/mol at 260–281 K, compared with ~5 × 10 –5 S/m and 34.8 kJ/m for CH 4 hydrate. Equivalent circuit modeling indicates that different pathways govern conduction in CO 2 versus CH 4 hydrate. Results show promise for use of electromagnetic methods in monitoring CO 2 hydrate formation in certain natural settings or in CO 2 /CH 4 exchange efforts.

58 GEOSCIENCES↗
DOE OSTI2021

Diffusion-to-Imbibition Transition in Water Sorption in Nanoporous Media: Theoretical Studies

The ability to predict multiphase fluid transport in nanoporous rocks such as shales is critical for many geoscience applications, for example unconventional hydrocarbon production, geologic carbon sequestration, and nuclear waste disposal. When the pore sizes approach nanoscales, the impact of the molecular interaction forces between fluids and solids becomes increasingly important. These forces can alter macroscopic fluid phase behavior and control transport. Recent experimental studies have shown that capillary condensation and subsequent imbibition of liquid water can occur in hydrophilic nanoporous media even if the vapor phase is at a critical relative humidity (rhcrit) well below vapor saturation. This study presents a theoretical investigation of the processes controlling adsorption, capillary condensation and imbibition in nanoporous media, using the square-gradient classical density functional theory. The proposed theoretical model explicitly includes the relevant interaction forces among fluids and solids in macroscopic porous media. Application of the model to a relative-humidity-controlled water adsorption experiment is presented to demonstrate the impact of water-pore wall attractive forces on multiphase water behavior in a hydrophilic silicon nanoporous medium. Here, the model represents well the measured time-dependent evolution of the water imbibition front inside the nanoporous medium and also explains the diffusion-like water transport regimes observed at rh < rh crit and the imbibition-like flow regimes observed at rh > rh crit . The study furthermore gives an insight on hysteresis phenomenon in adsorption and desorption isotherms.

58 GEOSCIENCES↗