SEARCH · Engineering Papers

Results for “capillary”

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

SourceResearch connectionFrom yearThrough year

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

At least 775 records · Page 43

All sources NASA NTRS DOE OSTI

DOE OSTI2021

From Langmuir isotherm to Brunauer–Emmett–Teller isotherm

Abstract Developed for multilayer adsorption, the Brunauer–Emmett–Teller ( BET) isotherm considers the adsorption of the first layer as an equilibrium chemical reaction between adsorbate molecules and adsorption sites and the adsorption of the second and subsequent layers as a condensation–evaporation process. Following the recent development of an activity‐based formulation for the Langmuir isotherm for monolayer adsorption, we present an activity‐based formulation for the BET isotherm in which species concentrations are replaced with species activities. Capturing the adsorbent surface heterogeneity for the adsorption of the first layer, the resulting thermodynamic BET isotherm is shown to accurately represent pure component adsorption isotherms over the relative pressure range of zero to unity or prior to the onset of capillary condensation. The thermodynamic BET isotherm should facilitate accurate estimation of monolayer adsorption capacity and the corresponding adsorbent surface area.

Vyawahare, Pradeep↗

DOE OSTI2024

Characterization of Porous Transport Layers Towards the Development of Efficient Proton Exchange Membrane Water Electrolysis

The current goals for implementing the hydrogen economy have highlighted a need to further optimize water-splitting technologies for clean hydrogen production. Proton exchange membrane water electrolysis (PEMWE) is a leading technology, but further optimizations of anode materials including the porous transport layer (PTL) and the adjacent catalyst layer (CL) are required to increase overall cell performance and reduce cost. This literature review describes advances in PTL development and characterization, highlighting early PTL characterization work and most common methods including capillary flow porometry and mercury intrusion porometry, optical imaging, neutron and x-ray radiography, and x-ray computed tomography. The article also discusses PTL protective coatings and their characterizations, focusing on platinum group metal (PGM)-based coatings, alternative non-PGM-based coatings, post-treated PTLs, and investigations into thin PGM-based coatings. Furthermore, it highlights the integration of the PTL and the adjacent CL along with associated characterization challenges. Lastly, this review discusses future developments in the characterization needed to improve PEMWE's performance and long-term durability are discussed.

08 HYDROGEN↗

DOE OSTI2022

Triptycene-Based Molecular Rods for Langmuir-Blodgett Monolayers

In this report we introduce fully modular synthesis leading to three representative examples of rigid molecular rods that are intended to form sturdy monolayers on various surfaces. These molecules contain two triptycene units that are designed to interlock into a compact "double-decker" structure. Two of the three final products provided suitable crystals for X-ray diffraction (analyzed on synchrotron), allowing deeper insight into packing in the 3-D crystal lattice. The acidity of all three compounds were determined by capillary electrophoresis, and the pK a values ranged between 2.06-2.53. All three rigid rods easily formed Langmuir-Blodgett monolayers (LBMs) on the water-air interfaces, with the area per molecule equal to 55-59 Å 2 /molecule, suggesting tight intermolecular packing. The thickness of all three films reached ~19 Å after transfer to a gold (111) surface, meaning that individual molecules are tilted maximally 38° from the axis perpendicular to the surface. The structure of one of these films on a gold (111) surface was visualized by AFM. These geometrically unique molecules represent promising platforms with a wide scope of applicability in the supramolecular architecture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Measuring Solvation Interactions of Deep Eutectic Solvents Formed by Metal Chlorides and an Imidazolium Salt by Inverse Gas Chromatography

Deep eutectic solvents (DESs) represent a class of solvents that offer a number of advantages including minimal toxicity, affordability, low vapor pressure, and simple, environmentally friendly preparation methods. They have found utility in areas, such as gas absorption, metal plating, and extractions. However, the relationship between their solvation properties and chemical composition remains poorly understood. In this study, a broad range of Type I DESs composed of metal chlorides and an imidazolium salt were prepared, employed as gas chromatographic stationary phases, and characterized using the Abraham solvation parameter model by inverse gas chromatography. The Abraham solvation parameter model allows for the study of DES solvation properties and the effects of varying structural components on system constants using a linear-free energy relationship. The DESs were investigated by systematically varying their composition, including the type of metal chloride and the molar ratio between the metal chloride and the imidazolium salt. The results show that hydrogen bond acidity, hydrogen bond basicity, and dipolarity/polarizability interactions are strongly influenced by the type of metal chloride within the DES and the ratio of metal chloride to imidazolium salt in the eutectic mixture. Furthermore, a column pretreatment procedure is presented that enables the effective coating of highly polar DESs onto open tubular capillary gas chromatography columns.

deep eutectic solvent↗

DOE OSTI2026

Continuous Counter‐Current Microfluidic Liquid–Liquid Extraction Achieved Using a Pair of Wettable Screen Meshes

Continuous counter‐current microfluidic liquid–liquid extraction performs separations by flowing immiscible liquids in opposing directions within a single flow channel. In principle, this flow arrangement enables a large number of theoretical separation units in a small footprint, without using interstage valving, pumping, and phase separation. Despite its potential for excellent separation performance, this microfluidic scheme rarely appears in literature due to the requirement for capillary forces to be greater than hydrodynamic forces for stable flow. We present a novel microfluidic device and flow approaches that overcome this force‐balance challenge, enabling stable, long‐duration continuous counter‐current flow. Additionally, we cover a suite of methodologies for quantifying the performance of the microfluidic device, revealing the number of theoretical equilibrium stages achieved. The enabling technologies include a woven mesh screen‐based microfluidic device architecture that is easily fabricated outside of a clean room, surface functionalization strategies to promote conjugate (organic/aqueous) wettability, flow approaches to eliminate bubbles and carryover, and computer‐aided flow automation with optical measurement of extraction performance. The reported experiments lasted for over 36 h, terminated only at experiment conclusion, where the device still exhibited good performance. Automated Raman spectroscopy was used for solute quantitation of the ternary system tert‐butanol in a toluene/water matrix, a ternary system that was specifically chosen to analyze the device's performance with a small solute partition ratio and to enable in‐line Raman measurements of solute concentrations in both phases. The microfluidic device possessed a 55 mm contact length and a 38.5 µL internal volume. During counter‐current flow, we observed approximately 37 equilibrium stages (37 ± 13) based on a best‐fit of the solute fraction remaining in the aqueous phase using a Kremser Group Method analysis.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Novel pore size-controlled, susceptibility matched, 3D-printed MRI phantoms

We report the design concept and fabrication of MRI phantoms, containing blocks of aligned microcapillaires that can be stacked into larger arrays to construct diameter distribution phantoms or fractured, to create a “powder-averaged” emulsion of randomly oriented blocks for vetting or calibrating advanced MRI methods, that is, diffusion tensor imaging, AxCaliber MRI, MAP-MRI, and multiple pulsed field gradient or double diffusion-encoded microstructure imaging methods. The goal was to create a susceptibility-matched microscopically anisotropic but macroscopically isotropic phantom with a ground truth diameter that could be used to vet advanced diffusion methods for diameter determination in fibrous tissues. Two-photon polymerization, a novel three-dimensional printing method is used to fabricate blocks of capillaries. Double diffusion encoding methods were employed and analyzed to estimate the expected MRI diameter. Susceptibility-matched microcapillary blocks or modules that can be assembled into large-scale MRI phantoms have been fabricated and measured using advanced diffusion methods, resulting in microscopic anisotropy and random orientation. This phantom can vet and calibrate various advanced MRI methods and multiple pulsed field gradient or diffusion-encoded microstructure imaging methods. We demonstrated that two double diffusion encoding methods underestimated the ground truth diameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Changes in water vapor adsorption and water film thickness in clayey materials as a function of relative humidity

Abstract Soil water, with adsorbed water being an important component, plays an important role in fluid flow and chemical movement in the unsaturated zone. The adsorbed water consists of water film and interlayer water, and its content is related to the relative humidity (RH). This study theoretically and experimentally focuses on the variation in adsorbed water content in clayey materials as the RH changes. Based on a slit‐pore model, three types of water (water film, interlayer water, and capillary water) could be present in the water vapor adsorption process. The variation of water film is obtained from the analysis of interfacial forces, and then the total water content can be quantitatively subdivided into these three types of water for different RH values. The slit‐pore adsorption model was used with experimental measurements (water vapor adsorption and N 2 physisorption) to quantify the amounts and interrelationships of these three different types of water to six clayey materials, namely three clay minerals (kaolinite, montmorillonite, and illite–smectite mixed layer) and three clay‐rich sediments from the Jianghan Plain in China. The main results are that (a) the volumes of water vapor adsorption are much greater than those of N 2 adsorption except for kaolinite; (b) interlayer water largely dominates the growth of total water content in montmorillonite and results in a concentration of 0.090 ml g −1 at 95% RH; and (c) the minimum thickness of water film is calculated to be 0.23 nm, and the maximum value is one‐third of the pore width by considering the interfacial forces.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Compressional wave velocity and effective stress in unsaturated soil: Potential application for monitoring moisture conditions in vadose zone sediments

Compressional seismic (P-wave) velocities were measured experimentally at different water potentials and confining pressure states using a variety of sediments collected from the US Department of Energy Hanford Site in southeastern Washington State. P-wave velocity was measured in variably saturated sediments using ultrasonic piezoelectric transducers and water potential was measured with a heat dissipation sensor. We propose a model to relate changes in P-wave velocity to soil stress conditions using an effective stress formulation that includes capillary stress as well as adsorptive stress. We show that compressional wave velocities can be related to water potential and confining stress using a model with five fitted parameters. A single set of parameters provides an excellent fit to data acquired from all of the sediment samples. Under very dry conditions, corresponding to low water potentials, water flow occurs in adsorptive thin films associated with very low unsaturated hydraulic conductivities. The results presented here suggest the utility of seismic methods for establishing low water potential and flow conditions within variably saturated Hanford formation sediments and providing feedback on the performance of mitigation efforts such as soil desiccation and surface infiltration barriers. Seismic methods can provide minimally invasive, three-dimensional information on relevant flow conditions in vadose zone sediments as well as changes that occur over time.

Linneman, Dorothy C.↗

DOE OSTI2022

Parameter estimation from spontaneous imbibition into volcanic tuff

Two-phase fluid flow properties underlie quantitative prediction of water and gas movement, but constraining these properties typically requires multiple time-consuming laboratory methods. The estimation of two-phase flow properties (van Genuchten parameters, porosity, and intrinsic permeability) is illustrated in cores of vitric nonwelded volcanic tuff using Bayesian parameter estimation that fits numerical models to observations from spontaneous imbibition experiments. The uniqueness and correlation of the estimated parameters is explored using different modeling assumptions and subsets of the observed data. The resulting estimation process is sensitive to both moisture retention and relative permeability functions, thereby offering a comprehensive method for constraining both functions. The data collected during this relatively simple laboratory experiment, used in conjunction with a numerical model and a global optimizer, result in a viable approach for augmenting more traditional capillary pressure data obtained from hanging water column, membrane plate extractor, or mercury intrusion methods. This method may be useful when imbibition rather than drainage parameters are sought, when larger samples (e.g., including heterogeneity or fractures) need to be tested that cannot be accommodated in more traditional methods, or when in educational laboratory settings.

58 GEOSCIENCES↗

DOE OSTI2020

Diverse Thaumarchaeota Dominate Subsurface Ammonia-oxidizing Communities in Semi-arid Floodplains in the Western United States

Subsurface microbial communities mediate biogeochemical transformations that drive both local and ecosystem-level cycling of essential elements, including nitrogen. However, their study has been largely limited to the deep ocean, terrestrial mines, caves, and topsoils (< 30 cm). Here, we present regional insights into the microbial ecology of aerobic ammonia oxidation within the terrestrial subsurface of five semi-arid riparian sites spanning a 900-km N-S transect. We sampled sediments, profiled communities to depths of ≤10 m, and compared them to reveal trends regionally within and surrounding the Upper Colorado River Basin (CRB). The diversity and abundance of ammonia-oxidizing microbial communities were evaluated in the context of subsurface geochemistry by applying a combination of amoA (encoding ammonia monooxygenase subunit A) gene sequencing, quantitative PCR, and geochemical techniques. Analysis of 898 amoA sequences from ammonia-oxidizing archaea (AOA) and bacteria (AOB) revealed extensive ecosystem-scale diversity, including archaeal amoA sequences from four of the five major AOA lineages currently found worldwide as well as distinct AOA ecotypes associated with naturally reduced zones (NRZs) and hydrogeochemical zones (unsaturated, capillary fringe, and saturated). Overall, AOA outnumber AOB by 2- to 5000-fold over this regional scale, suggesting that AOA may play a prominent biogeochemical role in nitrification within terrestrial subsurface sediments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Spatially Quasi-Periodic Water Waves of Infinite Depth

Abstract We formulate the two-dimensional gravity-capillary water wave equations in a spatially quasi-periodic setting and present a numerical study of solutions of the initial value problem. We propose a Fourier pseudo-spectral discretization of the equations of motion in which one-dimensional quasi-periodic functions are represented by two-dimensional periodic functions on a torus. We adopt a conformal mapping formulation and employ a quasi-periodic version of the Hilbert transform to determine the normal velocity of the free surface. Two methods of time-stepping the initial value problem are proposed, an explicit Runge–Kutta (ERK) method and an exponential time-differencing (ETD) scheme. The ETD approach makes use of the small-scale decomposition to eliminate stiffness due to surface tension. We perform a convergence study to compare the accuracy and efficiency of the methods on a traveling wave test problem. We also present an example of a periodic wave profile containing vertical tangent lines that is set in motion with a quasi-periodic velocity potential. As time evolves, each wave peak evolves differently, and only some of them overturn. Beyond water waves, we argue that spatial quasi-periodicity is a natural setting to study the dynamics of linear and nonlinear waves, offering a third option to the usual modeling assumption that solutions either evolve on a periodic domain or decay at infinity.

Wilkening, Jon (ORCID:0000000327827596)↗

DOE OSTI2022

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗

DOE OSTI2025

Thermal Response of a Lithium Vapor Divertor to Cyclical Operation

The lithium vapor divertor concept is being developed as a method to achieve detached divertor conditions in a tokamak while minimizing impurity radiation losses from the core plasma. SOLPS-ITER modeling has previously been used to identify some of the geometric constraints and required lithium evaporation rate of a lithium vapor divertor in a medium-sized tokamak during steady-state operation. Here an updated conceptual design based on these operating requirements is introduced and the thermal response of the system is modeled during cyclical operation, consistent with operation in a short-pulse tokamak. Controllability of the temperature of the lithium capillary porous system (CPS) is achieved by adopting a design where there is no line-of-sight for radiation from the plasma to reach the heated CPS surface. Operational strategies to minimize the amount of lithium evaporated between plasma discharges while achieving steady evaporation rates during plasma discharges are discussed and modeled here. The optimal feedforward control strategy demonstrated in this work is to ramp up the temperature of the evaporator as quickly as possible immediately before a plasma discharge and then reduce the heating to match the desired steady-state net evaporation rate just before the plasma discharge begins, allowing the thermal inertia of the system to stabilize the evaporation rate during the first second of the plasma discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Lithium Divertor Targets and Walls for the ASTER Liquid Stellarator Reactor, Distributed Divertor

Stellarators may have advantages for certain liquid metal options as Plasma Facing Components (PFC) for divertor targets and first walls due to the wide range of possible magnetic configurations, which additionally are free of disruptions and fast field variations. In a previous work (V. Queral et al., IEEE Trans. Plasma Sci. 52, 2024), a concept of stellarator reactor (ASTER-CP) based on swirling Li-molten salts and liquid lithium floating on the molten salt as PFC was presented. The divertor matters were not studied then and, thus, they are being studied and experimentally tested now. The ASTER-CP reactor concept, the initial liquid metal experiments and potential concepts for the ASTER-CP divertor and first wall are reported. Concerning the experiments, several small scale experiments of galinstan in a small rotating cylinder under magnetic field have been produced, including one experiment with high viscosity galinstan-mixture for increased thickness of layer. An experiment of floating lithium on the molten salt LiCl-PbCl2 gave fast volatilization/decomposition of the molten salt. Particularly for divertors, the traditional free-flow, Capillary Porous Systems and ‘divertorlets’ have been studied for application to ASTER-CP. Surface waves (hot spots), lack of enough surface fluid turbulence and excessive fluid speed are the main issues found in fast free-flow. The perhaps original concept of Distributed Divertor and Equi-power Surface is tentatively proposed and studied, taking advantage of stellarator fields and low recycling regime.

distributed divertor↗

DOE OSTI2021

Lattice Boltzmann Modeling of Drying of Porous Media Considering Contact Angle Hysteresis

Drying of porous media is governed by a combination of evaporation and movement of the liquid phase within the porous structure. Contact angle hysteresis induced by surface roughness is shown to influence multi-phase flows, such as contact line motion of droplet, phase distribution during drainage and coffee ring formed after droplet drying in constant contact radius mode. However, the influence of contact angle hysteresis on liquid drying in porous media is still an unanswered question. Lattice Boltzmann model (LBM) is an advanced numerical approach increasingly used to study phase change problems including drying. In this paper, based on a geometric formulation scheme to prescribe contact angle, we implement a contact angle hysteresis model within the framework of a two-phase pseudopotential LBM. The capability and accuracy of prescribing and automatically measuring contact angles over a large range are tested and validated by simulating droplets sitting on flat and curved surfaces. Afterward, the proposed contact angle hysteresis model is validated by modeling droplet drying on flat and curved surfaces. Then, drying of two connected capillary tubes is studied, considering the influence of different contact angle hysteresis ranges on drying dynamics. Finally, the model is applied to study drying of a dual-porosity porous medium, where phase distribution and drying rate are compared with and without contact angle hysteresis. The proposed model is shown to be capable of dealing with different contact angle hysteresis ranges accurately and of capturing the physical mechanisms during drying in different porous media including flat and curved geometries.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Experimental investigation of the effect of pore size distribution on nano-particle capture efficiency within ceramic particulate filters

The effect of the pore size distribution on size-resolved filtration efficiency was investigated for two ceramic particulate filters using particulate matter (PM) generated by a spark-ignition direct-injection engine fueled with gasoline. The cordierite filter tested had a porosity of 43%, a median pore diameter (d_50) of 12 µm, and a wide pore size distribution with a lognormal standard deviation (s’) of 0.3. The aluminum titanate filter had very similar porosity, d_50, and thickness, but significantly narrower pore size distribution (s’ = 0.1). The testing of two filters under identical experimental conditions enabled the impact of the pore size distribution on filtration performance to be evaluated. Filtration experiments were performed focusing on just the filter wall, starting from a clean filter until the transition to cake filtration (filtration efficiency > 99%). Time-resolved particle size distribution measurements were used to evaluate the progression of filtration performance and estimate trapped mass within the filter. The aluminum titanate filter, with a narrow pore size distribution, exhibited significantly better diffusion capture efficiency. The negative impact of higher flow velocity on diffusion capture efficiency was more pronounced for a narrower pore size distribution. Flow distributions measured using capillary-flow porometry were used to develop a cylindrical pore flow model to understand the impact of the differences in pore size distribution on observed trends in diffusion capture efficiency within a clean filter. The model helps demonstrate how large pores on the high end of a wider pore size distribution can carry a large fraction of the flow, with less efficient capture of small particles by diffusion. The experimental results and data demonstrate that the bubble point diameter and width of the pore size distribution significantly influence diffusion capture efficiencies for filters with very similar median pore diameter, porosity, and thickness.

aerosol filtration, Aftertreatment, ceramic exhaus↗

DOE OSTI2025

Interface morphology and dislocation-mediated processes during rapid solidification of thin films

Rapid solidification experiments have, in recent years, revealed a wealth of new microstructural phenomena that suggest a strong connection between the kinetics of solidification and the crystalline structures that emerge as a result. In this work, we investigate the interplay between interface morphology and defect-mediated processes during rapid solidification conditions using a Phase Field Crystal (PFC) model, enabling us to simultaneously and efficiently explore the physics of solidification and elasto-plasticity in the formalism of a single-field theory. We predict that there are two mechanisms by which dislocations emitted directly from the solid–liquid interface induce orientation gradients as well as the formation of subgrain boundaries within a single solidifying cell. We relate these mechanisms to the morphology of the moving solid–liquid interface and identify a suitable control parameter in the PFC model with which we can go between said morphologies by effectively changing the relative strength of the capillary length and kinetic coefficients of the solid–liquid interface. Thus, we are able to provide mechanistic explanations for several microstructural features (with an emphasis on orientation gradients and subgrain boundaries) observed during the rapid solidification of pure materials. We also provide a simple explanation for the formation of “jagged” subgrain boundaries, which is consistent with our experimental observations in rapidly solidified samples of Aluminum, whose mechanisms have thus far been unknown.

Interface morphology↗

DOE OSTI2020

The role of the spatial heterogeneity and correlation length of surface wettability on two-phase flow in a CO 2 -water-rock system

This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO 2 (scCO 2 )-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO 2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO 2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 × 800 × 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO 2 ). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO 2 and showed that the variations in the relative permeability curves for both water and scCO 2 increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. Furthermore, the impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.

58 GEOSCIENCES↗