Search NASA⌕ Search

SEARCH · Search NASA

Results for “wettability 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Application of unsupervised deep learning to image segmentation and in-situ contact angle measurements in a CO 2 -water-rock system

Rock surface wettability is a critical property that regulates multiphase flows in porous media, which can be quantified using the surface contact angle (CA). X-ray micro-computed tomography (μCT) provides an effective approach to in-situ measurements of surface CAs. However, the CA measurement accuracy depends significantly on the quality of CT image segmentation, which is the clustering of CT pixels into separate phases. Inspired by this, we developed a deep learning (DL)-based CA measurement workflow. Motivated by the recent tremendous progress in unsupervised learning techniques and aiming to avoid expensive manual data annotations, an unsupervised DL pipeline for CT image segmentation was proposed and implemented, which includes unsupervised model training and post-processing. The unsupervised model training was driven by a novel loss function constrained with feature similarity and spatial continuity and implemented by iterative forward and backward paths; the former clustered the pixel-wise feature vectors extracted by convolution neural networks, whereas the latter updated the parameters using gradient descent. An over-segmentation strategy was adopted for model training. The post-processing steps based on agglomerative hierarchical clustering (AHC) were implemented to further merge the over-segmented model output to the desired cluster number, which is intended to improve the efficiency of image segmentation. The developed unsupervised DL pipeline was compared with other commonly-used image segmentation methods using pixel-wise and physics-based evaluation metrics on a synthetic raw-image dataset, which had a known ground truth. The unsupervised DL pipeline showed the best performance. Next, the segmented images were input to an automatic CA measurement tool, and the results were validated by comparisons with manual measurements. The CA values from the manual and automatic measurements showed similar distributions and statistical properties. The automatic measurement demonstrated a wider spectrum because of the much larger number of measurement data points. The primary novelty of the unsupervised DL pipeline developed in this study lies in the novel loss function and the over-segmentation strategy associated with AHC post-processing. Finally, the workflow has been proven an efficient tool for pore-scale wettability characterization, which has a wide range of applications in fundamental studies of multiphase flows in natural porous media, which have critical implications to geological carbon sequestration, hydrocarbon energy recovery, and contaminant transport in groundwater.

42 ENGINEERING↗
DOE OSTI2025

Evaluating the Structural Response of Amphiphilic Monolayers to Environmental Stimuli

Amphiphilic monolayers composed of end groups with distinct polar and nonpolar functional groups offer rapid and reversible interfacial adaptation in response to environmental stimuli such as a change in interfacial medium polarity. We have synthesized and characterized a suite of monolayers with functional groups of competing polarity designed to reconfigure their interfacial chemical composition in response to solvent polarity. In these films, the end group is designed to be able to reorient and expose the functional groups that minimize the interfacial free energy between the film and the environment. Using a combination of spectroscopic, computational, and wettability studies, we have investigated the responsive interfacial behavior of different end groups upon exposure to environments with varying polarities. Contact angle measurements across a series of polar and dispersive probe liquids reveal trends that reflect the underlying molecular flexibility and composition. Vibrational sum frequency generation (SFG) spectroscopy and atomistic molecular dynamics (MD) simulations confirm solvent-driven reorientation of the end groups, with restructuring observed at the interface. To quantify these effects, we have developed a surface energy calculation model that incorporates solvent-induced surface rearrangements into the estimations. Our findings reveal a strong dependence of surface energy and switching behavior on the length and flexibility of the functionalities in the end group, which affects the exposure of certain interfacial compositions under different solvents. These results offer new insights into the design of adaptive monolayers and provide a framework for evaluating solvent-responsive surfaces.

functional groups↗
NASA NTRS1995

Evaporation of binary mixtures in microgravity

The motivation of this research is to obtain a better understanding of phase-change heat transfer within single and binary liquid meniscii, both in 1-g and 0-g environments. During phase 1 and part of phase 2, in a glass test cell with an inclined heated plate, 1-6 experiments on pentane with additions of decane up to 3% were conducted to determine the optimum concentration that will exhibit the maximum heat transfer and stability. During phase 2 emphasis was given to explore fundamental research issues and to ultimately develop a reliable capillary pumped loop (CPL) device for low gravity. In related experimental work, it was found that thermocapillary stresses near the contract line could result in a degraded wettability which ultimately could explain the observed failure of CPL devices in zero-gravity environment. Therefore, the current experimental effort investigates the effect of adding binary constituents in improving the thermocapillary characteristics near the contact line within the loop configuration. Achievements during second phase include: (1) Further enhancement of Central State University's Microgravity Laboratory by adding or improving upon capabilities of photography, video imaging, fluid visualization, and general experimental testing capabilities; (2) Experimental results for the inclined plate cell; (3) Modeling effort with a detailed scaling analysis; (4) Additional testing with a tube loop configuration to extend experimental work by Dickens, et al.; (5) Fabrication of a capillary loop to be tested using binary fluid (pentane/decane). The device that has been recently completed will be set up horizontally so that the effect of gravity on the performance is negligible. Testing will cover a wide range of parameters such as decane/pentane concentration, heat input value, heat input location (below or above meniscus), and loop temperature.

Girgis, Morris↗
DOE OSTI2026

Toward durable stacks: glass-ceramic sealants for intermediate-temperature protonic ceramic electrochemical systems

Protonic ceramic electrochemical cells (PCECs) are emerging as promising technologies for efficient energy conversion and hydrogen production because they operate at intermediate temperatures with improved efficiency and durability compared with conventional solid oxide electrochemical cells. However, the long-term reliability and commercialization of PCEC stacks remain strongly limited by the performance of sealants, which are required to maintain gas tightness, electrical insulation, and mechanical integrity under harsh thermal and chemical environments. Among various sealing approaches, glass-ceramic sealants are considered the most practical and scalable due to their excellent wettability, chemical tunability, and strong interfacial adhesion. This review provides a comprehensive overview of recent advances in glass-ceramic sealants for intermediate-temperature protonic ceramic electrochemical systems. The fundamental design principles of sealant compositions are first discussed, followed by recent developments in deposition methods, sintering strategies, surface treatments, and degradation monitoring techniques. Particular attention is given to the unique challenges associated with PCEC operating conditions, including hydrothermal degradation, interfacial reactions with barium-containing electrolytes, and thermal mismatch. Finally, future opportunities involving sustainable materials, multiscale modeling, additive manufacturing, and artificial intelligence-assisted sealant optimization are highlighted.

glass–ceramic sealants↗
DOE OSTI2024

Pore‐Scale Modeling of Carbon Dioxide and Hydrogen Transport During Geologic Gas Storage

Abstract Geologic storage of CO 2 and H 2 are climate‐positive techniques for meeting the energy transition. While similar formations could be considered for both gases, the flow dynamics could differ due to differences in their thermophysical properties. We conduct a rigorous pore‐scale study of water/CO 2 and water/H 2 systems at relevant reservoir conditions in a Bentheimer rock sample using the lattice Boltzmann method to quantify the effects of capillary, viscous, inertial, and wetting forces during gas invasion. At similar conditions, H 2 invasion is weaker compared to CO 2 due to unfavorable viscosity ratios. Increasing flow rate, however, increases the breakthrough saturation for both gas systems in the range of capillary numbers studied. At isolated conditions of flow rate, viscosity ratio, and wettability, local inertial effects are found to be critical and show consistent increase in the invaded gas saturation. The effect of inertial forces persits for both gases across all field conditions tested.

58 GEOSCIENCES↗
DOE OSTI2023

Large-scale physically accurate modelling of real proton exchange membrane fuel cell with deep learning

Proton exchange membrane fuel cells, consuming hydrogen and oxygen to generate clean electricity and water, suffer acute liquid water challenges. Accurate liquid water modelling is inherently challenging due to the multi-phase, multi-component, reactive dynamics within multi-scale, multi-layered porous media. In addition, currently inadequate imaging and modelling capabilities are limiting simulations to small areas (<1 mm 2 ) or simplified architectures. Herein, an advancement in water modelling is achieved using X-ray micro-computed tomography, deep learned super-resolution, multi-label segmentation, and direct multi-phase simulation. The resulting image is the most resolved domain (16 mm 2 with 700 nm voxel resolution) and the largest direct multi-phase flow simulation of a fuel cell. This generalisable approach unveils multi-scale water clustering and transport mechanisms over large dry and flooded areas in the gas diffusion layer and flow fields, paving the way for next generation proton exchange membrane fuel cells with optimised structures and wettabilities.

25 ENERGY STORAGE↗
NASA NTRS1992

Melt infiltration of silicon carbide compacts. I - Study of infiltration dynamics

Countergravity, pressure-assisted infiltration with a 2014 Al alloy of suitably tamped porous compacts of platelet shaped single crystals of alpha (hexagonal) silicon carbide was used to measure particulate wettability and infiltration kinetics under dynamic conditions relevant to pressure casting of composites. A threshold pressure P(th) for ingression of the infiltrant was identified based on the experimental penetration length versus pressure profiles for a range of experimental variables which included infiltration pressure, infiltration time, SiC size and SiC surface chemistry. The results showed that P(th) decreased whereas the penetration length increased with increasing SiC size and infiltration time. Cu-coated SiC led to lower P(th) and larger penetration lengths compared to uncoated SiC under identical conditions. These observations have been discussed in the light of theoretical models of infiltration and the kinetics of wetting.

Asthana, Rajiv↗
DOE OSTI2024

Interdependence of Support Wettability ‐ Electrodeposition Rate‐ Sodium Metal Anode and SEI Microstructure

Abstract This study examines how current collector support chemistry (sodiophilic intermetallic Na 2 Te vs. sodiophobic baseline Cu) and electrodeposition rate affect microstructure of sodium metal and its solid electrolyte interphase (SEI). Capacity and current (6 mAh cm −2 , 0.5–3 mA cm −2 ) representative of commercially relevant mass loading in anode‐free sodium metal battery (AF‐SMBs) are analyzed. Synchrotron X‐ray nanotomography and grazing‐incidence wide‐angle X‐ray scattering (GIWAXS) are combined with cryogenic ion beam (cryo‐FIB) microscopy. Highlighted are major differences in film morphology, internal porosity, and crystallographic preferred orientation e.g. (110) vs. (100) and (211) with support and deposition rate. Within the SEI, sodium fluoride (NaF) is more prevalent with Te−Cu versus sodium hydride (NaH) and sodium hydroxide (NaOH) with baseline Cu. Due to competitive grain growth the preferred orientation of sodium crystallites depends on film thickness. Mesoscale modeling delineates the role of SEI (ionic conductivity, morphology) on electrodeposit growth and onset of electrochemical instability.

Lo, Chang‐An↗
DOE OSTI2022

Sensitivity Analysis of Fluid–Fluid Interfacial Area, Phase Saturation and Phase Connectivity on Relative Permeability Estimation Using Machine Learning Algorithms

Recent studies have shown that relative permeability can be modeled as a state function which is independent of flow direction and dependent upon phase saturation (S), phase connectivity (X), and fluid–fluid interfacial area (A). This study evaluates the impact of each of the three state parameters (S, X, and A) in the estimation of relative permeability. The relative importance of the three state parameters in four separate quadrants of S-X-A space was evaluated using a machine learning algorithm (out-of-bag predictor importance method). The results show that relative permeability is sensitive to all the three parameters, S, X, and A, with varying magnitudes in each of the four quadrants at a constant value of wettability. We observe that the wetting-phase relative permeability is most sensitive to saturation, while the non-wetting phase is most sensitive to phase connectivity. Although the least important, fluid–fluid interfacial area is still important to make the relative permeability a more exact state function.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Interdependence of Support Wettability - Electrodeposition Rate- Sodium Metal Anode and SEI Microstructure

Here, this study examines how current collector support chemistry (sodiophilic intermetallic Na 2 Te vs. sodiophobic baseline Cu) and electrodeposition rate affect microstructure of sodium metal and its solid electrolyte interphase (SEI). Capacity and current (6 mAh cm -2 , 0.5–3 mA cm -2 ) representative of commercially relevant mass loading in anode-free sodium metal battery (AF-SMBs) are analyzed. Synchrotron X-ray nanotomography and grazing-incidence wide-angle X-ray scattering (GIWAXS) are combined with cryogenic ion beam (cryo-FIB) microscopy. Highlighted are major differences in film morphology, internal porosity, and crystallographic preferred orientation e.g. (110) vs. (100) and (211) with support and deposition rate. Within the SEI, sodium fluoride (NaF) is more prevalent with Te-Cu versus sodium hydride (NaH) and sodium hydroxide (NaOH) with baseline Cu. Due to competitive grain growth the preferred orientation of sodium crystallites depends on film thickness. Mesoscale modeling delineates the role of SEI (ionic conductivity, morphology) on electrodeposit growth and onset of electrochemical instability.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Generalized Analysis of Dynamic Flow Fouling on Heat Transfer Surfaces

Several energy applications involve flow of untreated heat transfer fluids that contain trace amounts of mineral salts, such as calcium carbonate and calcium sulfate, which show reverse solubility with temperature. The gradual accretion of these salts onto the heat transfer surfaces causes fouling that degrades thermal performance over time. While fouling of bare surfaces has been extensively studied, fouling of nonwetting superhydrophobic or liquid-infused surfaces is less understood, especially under dynamic flow conditions. This study presents a systematic experimental analysis to elucidate dynamic flow fouling of calcium carbonate and calcium sulfate on smooth, superhydrophobic and liquid-infused surfaces in a generalized manner. Here, copper tubes with modified inner surface wettability are fabricated and subjected to forced convection at different flow rates and fouling salt concentrations. Fouling characteristics of the two salts on the different surfaces are examined through microstructural analysis, and quantified in terms of asymptotic fouling resistance that measures the total extent of fouling, and a threshold time, which correlates to the onset of significant fouling. Using Taguchi design of experiments, the study reports, for the first time, closed-form analytical relationships for the asymptotic fouling thermal resistance and threshold time on the Reynolds number, dimensionless concentration of fouling agent and dimensionless infused liquid viscosity that represents the different surface types in a unified manner. It is shown that the analytical model accurately predicts the fouling characteristics in up to about 97% of the data, based on which contour maps of optimum surface designs for minimizing fouling resistance are presented. The results provide information on designing surfaces for reduced fouling or for estimating the fouling characteristics for a given surface in applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ion irradiation effects on Cr-coated zircaloy-4 surface wettability and pool boiling critical heat flux.

The concept of coating the currently used nuclear fuel cladding (zirconium-based alloy, typically Zircaloy-4 or Zirc-4) with an oxidation preventive layer is a progressing Accident tolerant Fuel (ATF) candidate alloys. The coated Zirc-4-based alloys could be a solution to suppress undesirable fast reaction kinetics with high-temperature steam. Zirc-4 has been the most preferred cladding material in pressurized water reactors (PWRs). Chromium (Cr) based alloys as a coating material provides excellent corrosion protection and good strength and wear resistance. In this paper, we present the surface wettability measurements and pool boiling Critical Heat Flux (CHF) for Cr-coated Zirc-4 claddings pre- and post-exposure to an ion irradiation environment. The wettability measurements, including static contact angle (contact angle, θ) and average surface roughness (surface roughness, Ra), are introduced for samples of different coating thicknesses (5–30 μm thick). The coatings fabricated by the cold spray of Cr-Al particles to 10 mm × 10 mm × 1.95 mm Zirc-4 substrates. Post fabrication, a Pilgering (cold rolling) process, was applied to finalize the coating thickness and resulted in a significant reduction in surface roughness of initially fabricated rough surfaces. The process produced three distinguished samples 5-μm unpolished (as machined), 5-μm, and 30-μm polished (cold rolled). The measurements are presented for the three surfaces and bare Zirc-4 as a baseline surface. The contact angle analyses were implemented in theoretical models from the literature to predict pool boiling CHF. Pool boiling experiments were conducted to measure the pool boiling CHF values and compare them to the predicted values. Scanning Electron Microscope (SEM) images and Energy Dispersive X-ray Spectroscopy (EDS) analysis was performed to characterize the surfaces for better understanding and interpreting the results. The SEM images showed localized surface damage due to ion irradiation. No recognized change in the measured surface roughness due to ion irradiation. The contact angles of irradiated Cr-coated surfaces are consistently higher (10°) than pre-irradiated surfaces. Decreasing the Cr-coating layer thickness resulted in lower contact angle pre- and post- ion irradiation. The predicted pool boiling CHF using the Kandlikar model is in good agreement with the experimentally measured CHF values within ±12% for all samples.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Molecular-scale origins of wettability at petroleum–brine–carbonate interfaces

Wettability control of carbonates is a central concept for enhanced petroleum recovery, but a mechanistic understanding of the associated molecular-scale chemical processes remains unclear. We directly probe the interface of calcium carbonate (calcite) with natural petroleum oil, synthetic petroleum analogues, and aqueous brines to understand the molecular scale behavior at this interface. The calcite–petroleum interface structure is similar whether or not calcite was previously exposed to an aqueous brine, and is characterized by an adsorbed interfacial layer, significant structural changes within the calcite surface, and increased surface roughness. No evidence for an often-assumed thin-brine wetting layer at the calcite–petroleum interface is observed. These features differ from those observed at the calcite–brine interface, and for parallel measurements using model synthetic petroleum mixtures (consisting of representative components, including dodecane, toluene, and asphaltene). Changes to the interface after petroleum displacement by aqueous brines are also discussed.

02 PETROLEUM↗
DOE OSTI2021

Technoeconomic Benefits of Film-Forming Amine Products Applied to Steam Surface Condensers

In a conventional Rankine cycle, the majority of power plants employ surface condensers that use pumped cooling water to reject heat from the cycle. In such cases, heat rejection occurs in a shell and tube heat exchanger by film-wise condensation of low-pressure steam on stainless steel, titanium, brass, or copper-nickel tubing. To improve the thermal performance of steam surface condensers, a replenishable film-forming substance (FFS) can be ap-plied to the condenser tubing to promote efficient dropwise condensation. Conventionally, film-forming amine prod-uct (FFAP) coatings protect boiler surfaces from oxidative corrosion, which substantially reduces the operation and maintenance costs. To quantify the technical and economic benefits of FFAP coatings applied to condenser tubing due to the promotion of dropwise condensation, a thermal resistance network model was established. Using a rep-resentative steam surface condenser, the improvements in thermal performance (overall heat transfer coefficient) and process parameters (net plant efficiency, cooling water flowrate, and turbine backpressure) were determined due to the enhancement in the condensation heat transfer coefficient. Experimentally measured condensation heat transfer coefficients for common condenser materials were compared with the modeling results and were found to be within attainable bounds. Lastly, the trend in total heat exchanger cost reduction is generalized to understand the trade-off between reduced surface area for heat rejection and increase in coating application costs for a re-plenishable coating system.

20 FOSSIL-FUELED POWER PLANTS↗
NASA NTRS1999

Enhanced Boiling on Micro-Configured Composite Surfaces Under Microgravity Conditions

In order to accommodate the growing thermal management needs of future space platforms, several two-phase active thermal control systems (ATCSs) have evolved and were included in the designs of space stations. Compared to the pumped single-phase liquid loops used in the conventional Space Transportation System and Spacelab, ATCSs offer significant benefits that may be realized by adopting a two-phase fluid-loop system. Alternately, dynamic power systems (DPSs), based on the Rankine cycle, seem inevitably to be required to supply the electrical power requirements of expanding space activities. Boiling heat transfer is one of the key technologies for both ATCSs and DPSs. Nucleate boiling near critical heat flux (CHF) can transport very large thermal loads with much smaller device size and much lower pumping power. However, boiling performance deteriorates in a reduced gravity environment and operation in the CHF regime is precarious because any slight overload will cause the heat transfer to suddenly move to the film boiling regime, which in turn, will result in burnout of the heat transfer surfaces. New materials, such as micro-configured metal-graphite composites, can provide a solution for boiling enhancement. It has been shown experimentally that this type of material manifests outstanding boiling heat transfer performance and their CHF is also extended to higher values. Due to the high thermal conductivity of graphite fiber (up to 1,200 W/m-K in the fiber direction), the composite surfaces are non-isothermal during the boiling process. The composite surfaces are believed to have a much wider safe operating region (a more uniform boiling curve in the CHF regime) because non-isothermal surfaces have been found to be less sensitive to variations of wall superheat in the CHF regime. The thermocapillary forces formed by the temperature difference between the fiber tips and the metal matrix play a more important role than the buoyancy in the bubble detachment, for the bubble detachment manifests itself by a necking process which should not be weakened by reduced gravity. In addition, the composite surfaces introduce no extra pressure drop, no fouling and do not impose significant primary or maintenance costs. All of these suggest that this type of composite is an ideal material for the challenge of accounting for both reliability and economy of the relevant components applied in the ATCSs, the DPSs and other devices in future space missions. The aim of the proposed work is to experimentally investigate high nucleate pool boiling performance on a micro-configured metal-graphite composite surface and to determine the mechanisms of the nucleate boiling heat transfer both experimentally and theoretically. Freon-113 and water will be used as the test liquids to investigate wettability effects on boiling characteristics. The Cu-Gr and Al-Gr composites with various volume fractions of graphite fibers will be tested to obtain the heat transfer characteristic data in the nucleate boiling region and in the CHF regime. In the experiments, the bubble emission and coalescence processes will be recorded by a video camera with a magnifying borescope probe immersed in the working fluid. The temperature profile in the thermal boundary layer on the composite surfaces will be measured by a group of micro thermocouples consisting of four ultra fine micro thermocouples. This instrument was developed and successfully used to measure the temperature profile of evaporating liquid thin layers by the proposers in a study performed at the NASA/Lewis Research Center. A two tier model to explain the nucleate boiling process and the performance enhancement on the composite surfaces has been suggested by the authors. According to the model, the thicknesses of the microlayer and the macrolayer underneath the bubbles and mushrooms, can be estimated by the geometry of the composite surface. The experimental results will be compared to the predictions from the model, and in turn, to revise and improve it.

Zhang, Nengli↗
DOE OSTI2020

High Pressure Infiltration/Expulsion of Aqueous NaCl in Planar Hydrophobic Nanopores

Pressure-driven permeation of water in a poorly wettable material results in a conversion of mechanical work into surface free energy representing a new form of energy storage, or absorption. When water is replaced by a concentrated electrolyte solution, the storage capacity of a nanoporous medium becomes comparable to high-end supercapacitors. The addition of salt can also reduce the hysteresis of the infiltration/expulsion cycle. Our molecular simulations provide a theoretical perspective into the mechanisms involved in the process, and underlying structures and interactions in compressed nanoconfined solutions. Specifically, we consider aqueous NaCl in planar confinements of widths of 1.0 nm and 1.64 nm and pressures of up to 3 kbar. Open ensemble Monte Carlo simulations utilizing fractional exchanges of molecules for efficient additions/removal of ions have been utilized in conjunction with pressure-dependent chemical potentials to model bulk phases under pressure. Confinements open to these pressurized bulk, aqueous electrolyte phases show reversibility at narrow pore sizes and irreversibility in wider ones, consistent with experiment, as well as strong hysteresis at both pore size. The addition of salt results in significant increases in the solid/liquid interfacial tension in narrower pores and associated infiltration and expulsion pressures. Furthermore, these changes are consistent with strong desalination effects at the lower pore size observed irrespective of external pressure and initial concentration.

25 ENERGY STORAGE↗
DOE OSTI2020

Two-Phase Fluid Flow Properties of Rough Fractures With Heterogeneous Wettability: Analysis With Lattice Boltzmann Simulations

Fractures are conduits for fluid flow in low-permeability geological formations. Multiphase flow properties of fractures are important in natural processes and in engineering applications such as the evaluation of the sealing capacity of caprocks and productivity of hydrocarbon-bearing tight rocks. Investigations of flow and transport through fractures typically focus on the effects of fracture geometric and mechanical factors such as aperture, roughness, and compressibility. The wettability of the fracture surfaces and its influence on microscale interfacial phenomena and macroscale effective transport properties are seldom studied. In this study, we investigated the effect of heterogeneous wetting properties on the displacement of water by supercritical CO2 through a series of lattice Boltzmann method simulations. The results show the evolution of the CO2 plume within a fracture is controlled by both the roughness of the aperture field and the wetting distribution. We combined these factors into a capillary pressure map that can be related to the macroscopic flow behavior of the fracture. We observed that heterogeneous wetting distributions promote the residual trapping of water where lower capillary pressures allowed for isolated water pockets in higher capillary pressure zones. Analysis of fracture unsteady relative permeability shows the effect of wetting on permeability evolution and provides support for the viscous-coupling relative permeability model. Finally, analysis of the steady-state relative permeability and saturation demonstrates a strong correlation between permeability and the standard deviation of the capillary pressure field. Thus, characterizing the distribution of wetting properties of fractures is crucial to understanding multiphase fracture flow and transport properties.

58 GEOSCIENCES↗