Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface wetting”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Origin of Life on Mars: Suitability and Opportunities

Although the habitability of early Mars is now well established, its suitability for conditions favorable to an independent origin of life (OoL) has been less certain. With continued exploration, evidence has mounted for a widespread diversity of physical and chemical conditions on Mars that mimic those variously hypothesized as settings in which life first arose on Earth. Mars has also provided water, energy sources, CHNOPS elements, critical catalytic transition metal elements, as well as B, Mg, Ca, Na and K, all of which are elements associated with life as we know it. With its highly favorable sulfur abundance and land/ocean ratio, early wet Mars remains a prime candidate for its own OoL, in many respects superior to Earth. The relatively well-preserved ancient surface of planet Mars helps inform the range of possible analogous conditions during the now-obliterated history of early Earth. Continued exploration of Mars also contributes to the understanding of the opportunities for settings enabling an OoL on exoplanets. Favoring geochemical sediment samples for eventual return to Earth will enhance assessments of the likelihood of a Martian OoL.

transition elements↗
DOE OSTI2021

Surface Sampling Techniques for the Canister Deposition Field Demonstration

This report describes plans for dust sampling and analysis for the multi-year Canister Deposition Field Demonstration. The demonstration will use three commercial 32PTH2 NUHOMS welded stainless steel storage canisters, which will be stored at an ISFSI site in Advanced Horizontal Storage Modules. One canister will be unheated; the other two will have heaters to achieve canister surface temperatures that match, to the degree possible, spent nuclear fuel (SNF) loaded canisters with heat loads of 10 kW and 40 kW. Surface sampling campaigns will take place on a yearly or bi-yearly basis. The goal of the planned dust sampling and analysis is to determine important environmental parameters that impact the potential occurrence of stress corrosion cracking on SNF dry storage canisters. Specifically, the size, morphology, and composition of the deposited dust and salt particles will be quantified, as well as the soluble salt load per unit area and the rate of deposition, as a function of canister surface temperature, location, time, and orientation. Sampling locations on the canister surface will nominally include 25 locations, corresponding to 5 circumferential locations at each of the 5 longitudinal locations. At each sampling location, a 2x2 sampling grid (containing 4 sample cells) will be painted onto the metal surface. During each sampling campaign, two samples at each sampling location will be collected, in a specific routine to measure both periodic (yearly or bi-yearly) and cumulative deposition rates. For each sample, a wet and a dry sample will be collected. Wet samples will be analyzed to determine the composition of the soluble salt fraction and to estimate salt loading per unit area. Dry samples will be analyzed to assess particle size, morphology, mineralogy, and identity (e.g. for floral/faunal fragments). The data generated by this proposed sampling plan will provide detailed information on dust and salt aerosol deposits on spent nuclear fuel canister surfaces. The anticipated results include information regarding particle compositions, size distributions, and morphologies, in addition to particle deposition rates as a function of canister surface location, orientation, time, and temperature. The information gathered during the Canister Deposition Field Demonstration is critical for ongoing efforts to develop a detailed understanding of the potential for stress corrosion cracking on SNF dry storage canisters

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Catalytic Dehydrogenation of Ethane over Mn Oxide Supported on Zeolite Chabazite

A new class of Mn-containing zeolites prepared by incipient wetness impregnation (IWI) have been found to catalyze the ethane dehydrogenation reaction with high selectivity (98 %+). Preparation by IWI leads to the formation of Mn 2 O 3 nanoparticles on the external surface of the zeolite crystals and herein is shown that the primary active sites for the reaction are located on the surface of these particles. Propane dehydrogenation is also successfully catalyzed by this catalyst. Furthermore, other Mn-zeolites (MFI and BEA) also have high reactivity and selectivity towards light alkane dehydrogenation.

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

Revisiting Selenium Interactions with Pyrite: From Adsorption to Coprecipitation

Interactions of selenium (Se), a trace element bio-essential at low concentrations but highly toxic at high concentrations, with the most abundant sulfide mineral in the Earth's crust, namely pyrite, was investigated over a wide range of time scales. At the nanosecond scale, selenate Se(VI)O 4 2– adsorption onto the net pyrite surface is shown by ab-initio computations to proceed via the formation of a chemical bond between an oxyanion oxygen atom and a surface Fe atom, weakening the other Se-O bonds and reducing Se atom oxidation state. At the hour-to-day scale, adsorption and coprecipitation of selenate Se(VI)O 4 2– and selenite, Se(IV)O 3 2– , were investigated through wet chemical batch experiments at various pH values at different sulfide concentrations. Selenium removal from solution is slower and weaker for selenate than for selenite. After 24h, only 10% of selenate, against 60% of selenite (up to 100% in the presence of sulfide), is removed by the pyrite surface. Independently of its original oxidation state, adsorbed Se is completely reduced to elemental trigonal selenium via adsorption, precipitation or coprecipitation, as shown by XANES spectroscopy. Here, our EXAFS results, compared to published data on Se-rich pyrite, show a Se to S substitution within the pyrite structure. The reductive coprecipitation mechanism of selenium with pyrite represents valuable new insights for improving our understanding of modern and ancient biogeochemical cycles involving Se. In addition, several industries can benefit from direct applications of our findings, such as water treatment, green technologies and sustainable mining.

Adsorption↗
DOE OSTI2022

Evaporation of Water Nanodroplets on Heated Surfaces: Does Nano Matter?

While experiments and continuum models have provided a relatively good understanding of the evaporation of macroscopic water droplets, elucidating how sessile nanodroplets evaporate is an open question critical for advancing nanotechnological applications where nanodroplets can play an essential role. In this study, using molecular dynamics simulations, we find that evaporating nanodroplets, in contrast to their macroscopic counterparts, are not always in thermal equilibrium with the substrate and that the vapor concentration on the nanodroplet surface does not reach a steady state. As a result, the evaporative behavior of nanodroplets is significantly different. Regardless of hydrophobicity, nanodroplets do not follow conventional evaporation modes but instead exhibit dynamic wetting behavior characterized by huge, non-equilibrium, isovolumetric fluctuations in the contact angle and contact radius. For hydrophilic nanodroplets, the evaporation rate, controlled by the vapor concentration, decays exponentially over time. Hydrophobic nanodroplets follow stretched exponential kinetics arising from the slower thermalization with the substrate. The evaporative half-lifetime of the nanodroplets is directly related to the thermalization time scale and therefore increases monotonically with the hydrophobicity of the substrate. Finally, the evaporative flux profile along the nanodroplet surface is highly nonuniform but does not diverge at the contact line as the macroscopic continuum models predict.

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

Methane emission hotspots in a boreal forest-fen mosaic potentially linked to deep taliks

Permafrost thaw is transforming boreal forests into mosaics of wetlands and drier uplands. Topographic controls on hydrological and ecological conditions impact methane (CH 4 ) fluxes, contributing to uncertainty in local and regional CH 4 budgets and underlying drivers. The objective of this study was to explore CH 4 fluxes and their drivers in a transitioning boreal forest-fen ecosystem (Goldstream Valley, Alaska, USA). This landscape is characterized by thawing discontinuous permafrost and heterogeneous mosaics of fens, collapse-scar channels, and small mounds of permafrost soils. From a survey in July 2021, observed chamber CH4 fluxes included fen areas with intermediate to very high emissions (29.8–635.3 mg CH 4 m −2 d −1 ), clustered locations with CH 4 uptake (−2.11 to −0.7 mg CH 4 m −2 d −1 ), and three anomalous emission hotspots (342.4–772.4 mg CH 4 m −2 d −1 ) that were located near samples with lower emissions. Some surface and near-surface variables partially explained the spatial variation in CH 4 flux. Log-transformed CH 4 flux had a positive linear relationship with soil moisture at 20 cm depth ( R 2 = 0.31, p -value < 1e-5) and negative linear relationships with microtopography ( R 2 = 0.13, p -value < 0.006) and slope ( R 2 = 0.28, p -value < 2e-5). Methane emissions generally occurred in flat, wet, graminoid-dominated fens, whereas CH 4 uptake occurred on permafrost mounds dominated by feather mosses and woody vegetation. However, the CH 4 hotspots occurred on drier, slightly sloped locations with low or undetectable near-surface methanogen abundance, suggesting that CH 4 was produced in deeper soils. When the hotspot samples were omitted, log-transformed CH 4 flux had a positive linear relationship with near-surface methanogen abundance ( R 2 = 0.29, p -value = 0.0023), and stronger linear relationships with soil moisture, slope, and soil macronutrient concentrations. Our findings suggest that some CH 4 emission hotspots could arise from CH 4 in deep taliks. The inference that methanogenesis occurs in deep taliks was strengthened by the identification of intrapermafrost taliks across the study area using low-frequency geophysical induction. This study assesses surface spatial heterogeneity in the context of subsurface permafrost conditions and highlights the complexity of CH 4 flux patterns in transitioning forest-wetland ecosystems. To better inform regional CH 4 budgets, further research is needed to understand the spatial distribution of terrestrial CH 4 hotspots and to resolve their surface, near-surface, and subsurface drivers.

boreal↗
DOE OSTI2021

Optimal surface-tension isotropy in the Rothman-Keller color-gradient lattice Boltzmann method for multiphase flow

The Rothman-Keller color-gradient (CG) lattice Boltzmann method is a popular method to simulate two-phase flow because of its ability to deal with fluids with large viscosity contrasts and a wide range of interfacial tensions. Here, two fluids are labeled red and blue, and the gradient in the color difference is used to compute the effect of interfacial tension. It is well known that finite-difference errors in the color-gradient calculation lead to anisotropy of interfacial tension and errors such as spurious currents. Here, we investigate the accuracy of the CG calculation for interfaces between fluids with several radii of curvature and find that the standard CG calculations lead to significant inaccuracy. Specifically, we observe significant anisotropy of the color gradient of order 7% for high curvature of an interface such as when a pinchout occurs. We derive a second order accurate color gradient and find that the diagonal nearest neighbors can be weighted differently than in the usual color-gradient calculation such that anisotropy is minimized to a fraction of a percent. The optimal weights that minimize anisotropy for the smallest radius of curvature interface are found to be w = (0.298, 0.284, 0.275) for diagonal nearest neighbors for the cases of the interface smoothing parameter β = (0.5, 0.7, 0.99), somewhat higher than the w = 0.25 value derived by Leclaire et al. [Leclaire, Reggio, and Trepanier, Computers and Fluids 48, 98 (2011)] based on obtaining isotropic errors to second order. We find that use of these optimal w values yields over a factor of 10 decrease in anisotropy and over a factor of 30 decrease in mean anisotropy relative to using the standard w = 1 value. And we find a factor of about 2 decrease in the anisotropic error and up to factor 15 decrease in mean anisotropic error relative to the choice of w = 0.25 for small radius of curvature interfaces. The improved CG calculations will allow the method to be more reliably applied to studies of phenomenology and pore scale processes such as viscous and capillary fingering, and droplet formation where surface-tension isotropy of narrow fingers and small droplets plays a crucial role in correctly capturing phenomenology. We present an example illustrating how different phenomena can be captured using the improved color-gradient method. Namely, we present simulations of a wetting fluid invading a fluid filled pipe where the viscosity ratio of fluids is unity in which droplets form at the transition to fingering using the improved CG calculations that are not captured using the standard CG calculations. We present an explanation of why this is so which relates to anisotropy of the surface tension, which inhibits the pinchouts needed to form droplets.

58 GEOSCIENCES↗
DOE OSTI2020

Tin-based ionic chaperone phases to improve low temperature molten sodium–NaSICON interfaces

High temperature operation of molten sodium batteries impacts cost, reliability, and lifetime, and has limited the widespread adoption of these grid-scale energy storage technologies. Poor charge transfer and high interfacial resistance between molten sodium and solid-state electrolytes, however, prevents the operation of molten sodium batteries at low temperatures. In this study, in situ formation of tin-based chaperone phases on solid state NaSICON ion conductor surfaces is shown in this work to greatly improve charge transfer and lower interfacial resistance in sodium symmetric cells operated at 110 °C at current densities up to an aggressive 50 mA cm –2 . It is shown that static wetting testing, as measured by the contact angle of molten sodium on NaSICON, does not accurately predict battery performance due to the dynamic formation of a chaperone NaSn phase during cycling. This work demonstrates the promise of sodium intermetallic-forming coatings for the advancement of low temperature molten sodium batteries by improved mating of sodium–NaSICON surfaces and reduced interfacial resistance.

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

Diminishing seasonality of subtropical water availability in a warmer world dominated by soil moisture–atmosphere feedbacks

Global warming is expected to cause wet seasons to get wetter and dry seasons to get drier, which would have broad social and ecological implications. However, the extent to which this seasonal paradigm holds over land remains unclear. Here we examine seasonal changes in surface water availability (precipitation minus evaporation, P–E) from CMIP5 and CMIP6 projections. While the P–E seasonal cycle does broadly intensify over much of the land surface, ~20% of land area experiences a diminished seasonal cycle, mostly over subtropical regions and the Amazon. Using land–atmosphere coupling experiments, we demonstrate that 63% of the seasonality reduction is driven by seasonally varying soil moisture (SM) feedbacks on P–E. Declining SM reduces evapotranspiration and modulates circulation to enhance moisture convergence and increase P–E in the dry season but not in the wet season. Our results underscore the importance of SM–atmosphere feedbacks for seasonal water availability changes in a warmer climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Land management effects on wet aggregate stability and carbon content

Land management affects soil structure and many other soil properties and processes. Our objectives were to evaluate soil organic C (SOC), aggregate size distribution, aggregate-associated C, and soil structure as affected by long-term land management and slope. A chronosequence of 38 on-farm sites with low to high (5–18%) slopes was selected to evaluate 5–40 yr of management. The sites were classified as business as usual (BAU) cropland (BAU-Crop), BAU pasture (BAU-Past), newly established conservation reserve program (CRP) areas (CRP-New), and established CRP (CRP-Old). Soil samples were collected from the 0-to-5- and 5-to-15-cm depth increments and processed for soil property measurements including fractionation by wet sieving into five aggregate size classes (>2,000, 1,000–2,000, 500–1,000, 250–500, and 53–250 μm). Within the surface 5 cm, mean weight diameter (MWD) and geometric mean diameter (GMD) were used to characterize soil structural stability. The BAU-Past and CRP-Old sites had 79% more macroaggregates (>2,000, 1,000–2,000, and 500–1,000 μm), 123% higher MWD, 38% higher GMD, and 47% higher SOC than BAU-Crop or CRP-New sites. The 5-to-15-cm depth increment showed a similar but lower magnitude response. Aggregate-associated C was quantified using a constant soil mass that reflected aggregate size distribution to prevent overestimating C content. Lower-slope locations had more SOC, more macroaggregates, more C associated with macroaggregates, and higher GMD and MWD compared with high-slope locations across all management classifications and soil depths. The results support our hypothesis that the high-slop soils may benefits from specific management decisions than the lower-sloping soils as a function of landscape property. We recommend reestablishing grassland on sloping land that is susceptible to excessive soil erosion, although those practices will likely take a long time to restore soil structural stability and SOC content to precultivation levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The effect of interfacial phenomena on gas solubility measurements in molten salts

The behavior of fission gases in molten fuel salt reactors governs activity transport from the reactor and can also affect the performance of the reactor itself. The gas solubility can be described thermodynamically by Henry’s law. However, the coupling of the condensed and gas phases depends on the interfacial area, which is difficult to measure or even to estimate. Surfaces of materials in the reactor will include disperse phases in the salt and porosity within the structural materials, covering a range of compositions and sizes. These attributes can affect measurements of fundamental properties such as gas solubility. Methods to obtain gas solubility, surface tension, interfacial energies, and bubble gas transport are reviewed. Recent data from manometric experiments are interpreted based on xenon sorption onto salt-wetted quartz.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

On the similarity of hillslope hydrologic function: a clustering approach based on groundwater changes

Abstract. Hillslope similarity is an active topic in hydrology because of its importance in improving our understanding of hydrologic processes and enabling comparisons and paired studies. In this study, we propose a holistic bottom-up hillslope clustering based on a region's integrative hydrodynamic response quantified by the seasonal changes in groundwater levels ΔP. The main advantage of the ΔP clustering is its ability to capture recharge and discharge processes. We test the performance of the ΔP clustering by comparing it to seven other common hillslope clustering approaches. These include clustering approaches based on the aridity index, topographic wetness index, elevation, land cover, and machine-learning that jointly integrate multiple data. We assess the ability of these clustering approaches to identify and categorize hillslopes with similar static characteristics, hydroclimate, land surface processes, and subsurface dynamics in a mountainous watershed – the East River – located in the headwaters of the Upper Colorado River Basin. The ΔP clustering performs very well in identifying hillslopes with six out of the nine characteristics studied. The variability among clusters as quantified by the coefficient of variation (0.2) is less in the ΔP and the machine learning approaches than in the others (> 0.3 for TWI, elevation, and land cover). We further demonstrate the robustness of the ΔP clustering by testing its ability to predict hillslope responses to wet and dry hydrologic conditions, of which it performs well when based on average conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Anomalous near-equilibrium capillary imbibition induced by nanoscale surface topography

Here, this work reports that the nanoscale surface topography induces a crossover from the classical Lucas–Washburn imbibition dynamics to a long-lived anomalously slow regime observed near thermodynamic equilibrium conditions. A compact analytical model considering the interplay between interfacial wetting, thermal motion, and the nanoscale physical topography accounts for experimental observations for the case of capillary rise of different nonvolatile liquids in simple glass tubes. Analytical predictions supported by topographic analysis via atomic force microscopy indicate that topographic features with a specific range of nanoscale dimensions determine the crossover condition and anomalous imbibition rate. Our findings have important implications for the scientific understanding and technical application of capillary imbibition and suggest strategies to control the adsorption of specific liquids in porous materials with complex surface topography.

Capillary driven imbibition↗
DOE OSTI2024

Molecular beam scattering of ammonia from a dodecane flat liquid jet

The evaporation and scattering of ND 3 from a dodecane flat liquid jet are investigated and the results are compared with previous studies on molecular beam scattering from liquid surfaces. Evaporation is well-described by a Maxwell–Boltzmann flux distribution with a cos θ angular distribution at the liquid temperature. Scattering experiments at E i = 28.8 kJ mol -1 over a range of deflection angles show evidence for impulsive scattering and thermal desorption. Further, at a deflection angle of 90°, the thermal desorption fraction is 0.49, which is higher than that of other molecules previously scattered from dodecane and consistent with work performed on NH 3 scattering from a squalane-wetted wheel. ND 3 scattering from dodecane results in super-specular scattering, as seen in previous experiments on dodecane. The impulsive scattering channel is fitted to a “soft-sphere” model, yielding an effective surface mass of 55 amu and an internal excitation of 5.08 kJ mol -1 . Overall, impulsively scattered ND 3 behaves similarly to other small molecules scattered from dodecane.

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

Catalyst Deactivation Modes of PdO/γ-Al 2 O 3 Catalysts for Lean Methane Oxidation

PdO/γ-Al 2 O 3 catalysts are one of the most active catalytic components for the complete oxidation of methane. Under reaction conditions, especially in a wet feed, the catalysts suffer severe performance degradation. This study establishes a series of testing protocols to systematically investigate the causes of catalyst deactivation under methane oxidation reaction conditions. Four distinct catalyst deactivation modes are identified. Two of the deactivation modes are directly related to H 2 O, either from the feed gas or as a part of the reaction products, with one (Mode 2) being attributed to the formation of surface hydroxyl groups and the other (Mode 3) to the competitive adsorption of H 2 O on the catalysts. The impact of the two deactivation modes is acute and severe but reversible. In contrast, the other two deactivation modes are gradual and persistent but irreversible. Both modes are induced by CH 4 oxidation reaction, with the impact of a wet feed (Mode 4) being substantially more severe than that of a dry feed (Mode 1). The major cause of the irreversible catalyst deactivation is attributed to surface reconstruction of PdO nanoparticles, which behaves as a passivation layer lowering the number of coordinately unsaturated Pd sites for CH 4 activation. Although the passivation layer is relatively stable against thermal or hydrothermal treatment, it is not completely inert. Formation and partial regeneration of the passivation layer is a highly dynamic process and heavily depends on the reaction temperature: a lower reaction temperature (≤ 450 ℃) can lead to quicker catalyst deactivation; but a higher reaction temperature (between 500 – 550 ℃) can result in a greater extent of catalyst deactivation.

PdO/γ-Al2O3↗
DOE OSTI2025

Surface modification of cathode material enhances electrochemical performance in dry-processed Li-ion battery electrodes

The transition to electric vehicles (EVs) is pivotal for achieving energy security and integrating grid stability, with lithium-ion batteries (LIBs) playing a central role in this transformation. However, the conventional wet electrode manufacturing relying on N-methyl-2-pyrrolidone (NMP) solvent is energy intensive and costly. Dry processing (DP) has emerged as a promising alternative, eliminating solvents and using polytetrafluoroethylene (PTFE) binder for electrode fabrications. Despite its advantages, DP faces a critical challenge: poor interfacial adhesion between the hydrophobic PTFE binder and the hydrophilic cathode active material (CAM), particularly LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811), which undermines electrode performance. Here, to address this, we introduced a novel vapor-phase trimethoxymethylsilane (TMMS) coating to hydrophobize the CAM surface, enhancing compatibility with PTFE. This surface modification significantly enhances binder – CAM interactions, enabling uniform mixing and robust electrode integrity without damaging the CAM particles. Our findings advance the feasibility of environmentally sustainable and cost-effective dry processing, representing a significant step toward sustainable battery manufacturing.

Choi, Junbin [Oak Ridge National Laboratory (ORNL)↗
DOE OSTI2026

Testing of spark plasma sintered porous tungsten under neon glow discharge cleaning conditions in LTX-β

This paper demonstrates that tungsten (W) based powder reconstituted Plasma Facing Components (PFCs) can be treated in situ in a fusion reactor to remove W oxide and carbon (C) contamination. Doing so should ease the challenge of using these materials in a Capillary Porous System (CPS) with lithium (Li) by enabling better wetting and less contamination of the Li by the underlying CPS. Most powder reconstituted materials including 3D printed and sintered PFCs suffer from a high surface contamination from oxides and surface C, which complicates their use with liquid Li, a primary PFC candidate. Spark plasma sintered porous W samples were fabricated to be used as a CPS with liquid Li. The samples were characterized in terms of morphology and surface chemistry. Analysis confirms a high C and oxygen (O) contamination. We present the results of exposing this type of CPS to Glow Discharge Cleaning (GDC) cycles in the Lithium Tokamak Experiment-β (LTX-β). The sample was exposed to neon (Ne) GDC in the midplane of the low-field side of LTX-β and analyzed in vacuo with Temperature Programmed Desorption (TPD) and Secondary Ion Mass Spectrometry (SIMS) to investigate the effects the Ne GDC had on the chemical composition of the sample. The combination of Ne GDC with rapid heating as done in TPD was successful in reducing the W oxides and removing the C contamination.

Capillary porous system↗
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↗