Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface modeling”

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

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

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Effect of Ni Dopant on Furan Activation over Mo 2 C Surface: Insights from First-Principles-Based Microkinetic Modeling

Here, the role of Ni dopant in a molybdenum carbide (Mo 2 C) catalyst for butadiene and propyne production from furan is investigated by using a Mo 2 C(001) surface model and first-principles-based microkinetic modeling (MKM). A reaction mechanism including 18 reaction steps is generated to compute the turnover frequency (TOF) of butadiene and propyne from furan gas. The butadiene production is dominant with the TOF being 3-6 orders of magnitude higher than the TOF of propyne on a bare Mo 2 C surface due to an active C-O bond scission reaction compared to C-C bond cleavage reaction. The exothermic binding of ring-opened furan (CHCHCHCHO*) inhibits the activity of a bare Mo 2 C surface for butadiene and propyne production. Calculations predict that Mo 2 C surface modification with a Ni dopant energetically destabilizes the adsorption of CHCHCHCHO*, which facilitates the formation of butadiene and propyne. In the presence of Ni dopants on a Mo 2 C(001) surface, the selectivity toward the butadiene formation is decreased compared to propyne due to relatively favorable C-C bond scission as opposed to the C-O bond scission. This first-principles-based study provides mechanistic and kinetic insights into the role of dopants in cost-efficient Mo 2 C catalysts during the vapor phase conversion of furan to gaseous products.

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

A Bezier Curve Informed Melt Pool Geometry to Model Additive Manufacturing Microstructures Using SPPARKS

Additive manufacturing is a transformative technology with the potential to manufacture designs which traditional subtractive machining methods cannot. Additive manufacturing offers fast builds at near final desired geometry; however, material properties and variability from part to part remain a challenge for certification and qualification of metallic components. AM induced metallic microstructures are spatially heterogeneous and highly process dependent. Engineering properties such as strength and toughness are significantly affected by microstructure morphologies resulting from the manufacturing process Linking process parameters to microstructures and ultimately to the dynamic response of AM materials is critical to certifying and qualifying AM built parts and components and improving the performance of AM materials. The AM fabrication process is characterized by building parts layer by layer using a selective laser melt process guided by a computer. A laser selectively scans and melts metal according to a designated geometry. As the laser scans, metal melts, fuses, and solidifies forming the final geometry in a layerwise fashion. As the laser heat source moves away, the metal cools and solidifies forming metallic microstructures. This work describes a microstructure modeling application implemented in the SPPARKS kinetic Monte Carlo computational framework for simulating the resulting microstructures. The application uses Bzier curves and surfaces to model the melt pool surface and spatial temperature profile induced by moving the laser heat source; it simulates the melting and fusing of metal at the laser hot spot and microstructure formation and evolution when the laser moves away. The geometry of the melt pool is quite flexible and we explore effects of variances in model parameters on simulated microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electronic Structure of InAs and InSb Surfaces: Density Functional Theory and Angle-Resolved Photoemission Spectroscopy

The electronic structure of surfaces plays a key role in the properties of quantum devices. However, surfaces are also the most challenging to simulate and engineer. Here, in this work, the electronic structure of InAs(001), InAs(111), and InSb(110) surfaces is studied using a combination of density functional theory (DFT) and angle-resolved photoemission spectroscopy (ARPES). Large-scale first principles simulations are enabled by using DFT calculations with a machine-learned Hubbard U correction [npj Comput. Mater. 6, 180 (2020)]. To facilitate direct comparison with ARPES results, a “bulk unfolding” scheme is implemented by projecting the calculated band structure of a supercell surface slab model onto the bulk primitive cell. For all three surfaces, a good agreement is found between DFT calculations and ARPES. For InAs(001), the simulations clarify the effect of the surface reconstruction. Different reconstructions are found to produce distinctive surface states, which may be detected by ARPES with low photon energies. For InAs(111) and InSb(110), the simulations help elucidate the effect of oxidation. Owing to larger charge transfer from As to O than from Sb to O, oxidation of InAs(111) leads to significant band bending and produces an electron pocket, whereas oxidation of InSb(110) does not. The combined theoretical and experimental results may inform the design of quantum devices based on InAs and InSb semiconductors, for example, topological qubits utilizing the Majorana zero modes.

42 ENGINEERING↗
DOE OSTI2023

Community Based Data of Potentiometric Titration of Iron Oxides: Ferrihydrite (HFO), Goethite, Hematite, Magnetite

This data release includes experimental data of potentiometric titration for iron oxides. The data in the provided .csv files is not our own experimental data but have been compiled from the multiple literature sources. The master database is L-SCIE (LLNL Surface Complexation/Ion Exchange) database, and the provided .csv files are extracted data from L-SCIE. The .csv files were obtained by using the Lawrence Livermore National Laboratory Surface Complexation Database Converter (SCDC) code written in the R programming language (free licensing available at https://ipo.llnl.gov/technologies/software/llnl-surface-complexation-database-converter-scdc).The released data was used for developing a comprehensive community data-driven surface complexation modeling (SCM) framework for simulating potentiometric titration of mineral surfaces. Compiled community data for ferrihydrite, goethite, hematite, and magnetite are fit to produce representative protolysis constants that can reproduce potentiometric titration data collected from multiple literature sources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Re-evaluating cloud chamber constraints on depositional ice growth in cirrus clouds – Part 1: Model description and sensitivity tests

Abstract. Ice growth from vapor deposition is an important process for the evolution of cirrus clouds, but the physics of depositional ice growth at the low temperatures (<235 K) characteristic of the upper troposphere and lower stratosphere is not well understood. Surface attachment kinetics, generally parameterized as a deposition coefficient αD, control ice crystal habit and also may limit growth rates in certain cases, but significant discrepancies between experimental measurements have not been satisfactorily explained. Experiments on single ice crystals have previously indicated the deposition coefficient is a function of temperature and supersaturation, consistent with growth mechanisms controlled by the crystal's surface characteristics. Here we use observations from cloud chamber experiments in the Aerosol Interactions and Dynamics in the Atmosphere (AIDA) aerosol and cloud chamber to evaluate surface kinetic models in realistic cirrus conditions. These experiments have rapidly changing temperature, pressure, and ice supersaturation such that depositional ice growth may evolve from diffusion limited to surface kinetics limited over the course of a single experiment. In Part 1, we describe the adaptation of a Lagrangian parcel model with the Diffusion Surface Kinetics Ice Crystal Evolution (DiSKICE) model (Zhang and Harrington, 2014) to the AIDA chamber experiments. We compare the observed ice water content and saturation ratios to that derived under varying assumptions for ice surface growth mechanisms for experiments simulating ice clouds between 180 and 235 K and pressures between 150 and 300 hPa. We found that both heterogeneous and homogeneous nucleation experiments at higher temperatures (>205 K) could generally be modeled consistently with either a constant deposition coefficient or the DiSKICE model assuming growth on isometric crystals via abundant surface dislocations. Lower-temperature experiments showed more significant deviations from any depositional growth model, with different ice growth rates for heterogeneous and homogeneous nucleation experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Adsorption of REEs to Kaolinite via Ion Exchange and Surface Complexation as a Function of Water Chemistry

Rare earth elements (REEs) are critical components of modern technology behind renewable energy, transportation, and electronics but have a limited current supply. A substantial portion of global REE production relies on ion adsorption deposits. A high abundance of kaolinite in REE enrichment zones within these deposits suggests that kaolinite controls the subsurface migration of REEs. This study aimed to improve the current understanding of REE binding to kaolinite under varying water chemistry conditions. We conducted batch experiments with kaolinite (KGa-2) and three REEs (Nd, Dy, and Yb) at varying pH, electrolyte concentration, dissolved inorganic carbon (DIC), low molecular weight organic acids (citric and oxalic acids), and total REE concentration conditions. Increasing electrolyte concentration inhibits REE adsorption at pH < 7, suggesting that ion exchange contributes to adsorption at these pH values. DIC affects adsorption above pH 7–8 by forming strong aqueous complexes with heavy REEs. Citric acid decreases REE adsorption via aqueous complexation of REEs at pH > 5 but does not affect adsorption at pH < 5. The surface complexation model captures the main adsorption trends with two mechanisms: ion exchange on basal planes at pH < ∼6 and inner-sphere surface complexation to edge sites at pH > ∼6. Equilibrium constants for surface complexation increase in the order of Yb > Dy > Nd, indicating a higher strength of adsorption for heavy REEs. This study demonstrates how water chemistry conditions control the adsorption mechanisms that may determine the mobility of REEs in subsurface environments rich in kaolinite.

58 GEOSCIENCES↗
DOE OSTI2025

Adsorption of Rare Earth Elements to Kaolinite [dataset]

Rare earth elements (REEs) are critical components of modern technology behind renewable energy, transportation, and electronics but have a limited current supply. A substantial portion of global REE production relies on ion adsorption deposits. A high abundance of kaolinite in REE enrichment zones within these deposits suggests that kaolinite controls the subsurface migration of REEs. This study aimed to improve the current understanding of REE binding to kaolinite under varying water chemistry conditions. We conducted batch experiments with kaolinite (KGa-2) and three REEs (Nd, Dy, and Yb) at varying pH, electrolyte concentration, dissolved inorganic carbon (DIC), low molecular weight organic acids (citric and oxalic acids), and total REE concentration conditions. Increasing electrolyte concentration inhibits REE adsorption at pH 7, suggesting that ion exchange contributes to adsorption at these pH values. DIC affects adsorption above pH 7−8 by forming strong aqueous complexes with heavy REEs. Citric acid decreases REE adsorption via aqueous complexation of REEs at pH 5 but does not affect adsorption at pH 5. The surface complexation model captures the main adsorption trends with two mechanisms: ion exchange on basal planes at pH ∼6 and inner-sphere surface complexation to edge sites at pH ∼6. Equilibrium constants for surface complexation increase in the order of Yb Dy Nd, indicating a higher strength of adsorption for heavy REEs. This study demonstrates how water chemistry conditions control the adsorption mechanisms that may determine the mobility of REEs in subsurface environments rich in kaolinite.

58 GEOSCIENCES↗
DOE OSTI2021

Predictability and feedbacks of the ocean-soil-plant-atmosphere water cycle: deep learning water conductance in Earth System Model

This white paper responds to Focal Area 2. We seek to build predictive models of leaf and surface conductance of water by implementing deep learning (DL) data assimilation techniques. These new models would then be implemented in existing Land Surface Models (LSMs) and Earth System models (ESMs), generating novel water cycle feedbacks. In doing so, we would improve predictability of expected changes in land precipitation, soil moisture, and vegetation dynamics in the long-term, and the role of land cover on the impacts and feedbacks of extreme weather events in the short-term

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Uranium Retardation Capacity of Lithologies from the Negev Desert, Israel—Rock Characterization and Sorption Experiments

A series of batch experiments were performed to assess the uranium sorption capacity of four mineralogically distinct lithologies from the Negev Desert, Israel, to evaluate the suitability of a potential site for subsurface radioactive waste disposal. The rock specimens consisted of an organic-rich phosphorite, a bituminous marl, a chalk, and a sandstone. The sorption data for each lithology were fitted using a general composite surface complexation model (GC SCM) implemented in PHREEQC. Sorption data were also fitted by a non-mechanistic Langmuir sorption isotherm, which can be used as an alternative to the GC SCM to provide a more computationally efficient method for uranium sorption. This is because all the rocks tested have high pH/alkalinity/calcium buffering capacities that restrict groundwater chemistry variations, so that the use of a GC SCM is not advantageous. The mineralogy of the rocks points to several dominant sorption phases for uranyl (UO 2 2+ ), including apatite, organic carbon, clays, and iron-bearing phases. The surface complexation parameters based on literature values for the minerals identified overestimate the uranium sorption capacities, so that for our application, an empirical approach that makes direct use of the experimental data to estimate mineral-specific sorption parameters appears to be more practical for predicting uranium sorption.

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

Exploring topography-based methods for downscaling subgrid precipitation for use in Earth System Models

Topography exerts major control on land surface processes. To improve representation of topographic impacts on land surface processes, a new topography-based subgrid structure has been introduced to the Energy Exascale Earth System Model to represent the subgrid heterogeneity of surface elevation. To take advantage of the new subgrid structure for improving land surface modeling, this study explores four variations of topography-based methods for downscaling grid precipitation to the corresponding subgrids. In the first three methods, the deviation of the subgrid precipitation from the grid’s is equal to the grid precipitation multiplied by the ratio of the elevation difference between the subgrid and grid mean to a specified elevation equals to the grid elevation, the difference between the maximum and minimum subgrid elevation, and the maximum subgrid elevation, respectively. The second method limits the ratio to 0.5 to avoid extreme values on mountains and the third method accounts for the slope effect. The fourth method is similar to the third method except that the Froude Number is used to limit the numerator in a blocking regime of the ambient flow. The downscaled precipitation is evaluated using the PRISM precipitation data over the U.S. using statistical metrics. Results show that by accounting for topographic slope besides elevation, the third and fourth methods show clear advantages over the first and second methods. Furthermore, introducing the Froude Number in the fourth method improves downscaling skill and shows consistent advantages over the third method in areas with larger subgrid heterogeneity across different grids sizes.

Tesfa, Teklu K.↗
DOE OSTI2023

Effects of Size and Shape on the Tolerances for Misalignment and Probabilities for Successful Oriented Attachment of Nanoparticles

Oriented attachment (OA) of nanoparticles is an important pathway of crystal growth, but tools for quantitatively modeling OA are lacking. Here we present several simple models that relate the probability of achieving OA to basic geometric parameters such as particle size, shape, and lattice periodicity. A Moiré-domain model is applied to understand twist-misorientations between parallel surfaces, and it predicts that the range of twist angles yielding perfect OA is inversely related to the width of the contact area. This is confirmed and further developed using a surface functional model, which predicts how crystallographic registration forces drive the emergence of complex orientational energy landscapes. The energy landscapes are predicted to possess local minima that can trap particles in imperfect alignments, and these local minima become deeper and more numerous as the contact area increases, making OA more challenging for large particles. Further, a second set of models is presented to understand the sequence of events by which two crystallographic faces become co-planer after collision. We use a ‘central force approximation’ to quantitatively predict the odds of attaining coalignment between various faces when particles collide with random misalignments, and we show that in the absence of biasing forces, the probability of attaining alignment on a given face is roughly proportional to its solid angle as viewed from the center of the particle. The model predicts that OA is most favorable between well-faceted particles and becomes exceedingly unlikely for large spherical particles that express many microfacets.

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

In Situ Characterization of Surface Recombination in p-Si/SiO x Based Photoelectrochemical Cells

Time-resolved infrared (TRIR) and electrochemical impedance spectroscopy (EIS) were utilized to quantify surface states present at silicon oxide (SiO x ) protected crystalline p-Si electrolyte interfaces. The primary goal was to identify p-Si/SiO x photoelectrodes with both low surface recombination rates and efficient multi-electron transfer to an acceptor present in the external electrolyte. Three SiO x layers were investigated: native oxide (nOx), chemical oxide (cOx), and rapid thermal annealed (RTA) thermal oxide (tOx). Comparative study with [Ru(bpy) 3 ](PF 6 ) 2 as the electron acceptor indicated that tOx was most optimal with a small effective recombination rate, multi-electron transfer capability, and photovoltage of 500 ± 50 mV. A secondary goal was to analyze the surface recombination rates with the Shockley–Read–Hall (SRH) kinetic model. Two surface states were identified from this analysis, one closer to the CB edge (V t,1 ) and the other near the midgap (V t,2 ). EIS and SRH analyses revealed that a forming gas (5% H 2 /N 2 ) anneal (FGA) decreased surface recombination for tOx and nOx through a lower density of surface states. In the case of tOx, the infrared data indicated that V t,2 was completely removed. Here, the energetic positions of the band edges were correlated with the surface state density; low densities corresponded to more favorable potentials for inversion layer formation, which is expected to be most optimal for photocatalysis. Collectively this study indicates that the free carrier dynamics provided by TRIR represent a powerful in situ probe of the band edge and the surface state energetics in silicon based photoelectrochemical cells.

Electrochemical Impedance Spectroscopy↗
DOE OSTI2023

Resolving experimental biases in the interpretation of diffusion experiments with a user-friendly numerical reactive transport approach

The reactive transport code CrunchClay was used to derive effective diffusion coefficients (D e ), clay porosities (ε), and adsorption distribution coefficients (K D ) from through-diffusion data while considering accurately the influence of unavoidable experimental biases on the estimation of these diffusion parameters. These effects include the presence of filters holding the solid sample in place, the variations in concentration gradients across the diffusion cell due to sampling events, the impact of tubing/dead volumes on the estimation of diffusive fluxes and sample porosity, and the effects of O-ring-filter setups on the delivery of solutions to the clay packing. Doing so, the direct modeling of the measurements of (radio)tracer concentrations in reservoirs is more accurate than that of data converted directly into diffusive fluxes. While the above-mentioned effects have already been described individually in the literature, a consistent modeling approach addressing all these issues at the same time has never been described nor made easily available to the community. A graphical user interface, CrunchEase, was created, which supports the user by automating the creation of input files, the running of simulations, and the extraction and comparison of data and simulation results. While a classical model considering an effective diffusion coefficient, a porosity and a solid/solution distribution coefficient (D e –ε–K D ) may be implemented in any reactive transport code, the development of CrunchEase makes it easy to apply by experimentalists without a background in reactive transport modeling. CrunchEase makes it also possible to transition more easily from a D e –ε–K D modeling approach to a state-of-the-art process-based understanding modeling approach using the full capabilities of CrunchClay, which include surface complexation modeling and a multi-porosity description of the clay packing with charged diffuse layers.

58 GEOSCIENCES↗
DOE OSTI2021

Modelling the interlayer bond strength of 3D printed concrete with surface moisture

Highlights: • An analytical model for interlayer adhesion of 3DPC is proposed based on the amount of surface moisture on the interlayer. • For this mix, the model accurately (RMSE=2.5%) predicts the reduction in interlayer adhesion from 30% to 50%. • Internal curing with superabsorbent polymers increases the interlayer adhesion by 10% and flexural strength by 19%. • SAPs also improve the initial thixotropy by 49% and buildability by 25% compared to the reference. Providing additional water to the hydrating cementitious particles is essential to achieve the optimal mechanical performance of the low w/b concrete mixes preferred for 3D printing. This study incorporates superabsorbent polymers (SAP) and additional water in 3D printed concrete (3DPC) to promote the hydration process through delayed internal water release. The study shows that a retentive SAP modifies the rheological development by absorbing the pore fluid for a short period after printing. The absorption-induced stiffening increases the thixotropy and buildability by 49% and 25%, respectively. A retentive SAP increases the flexural strength and interlayer adhesion by 19% and 10%, respectively. This is due to the internal water release that promotes hydration. Evaporation of the interlayer moisture during the pass time has the opposite effect— evaporation reduces the interlayer adhesion. Based on this assumption, an analytical model is proposed. The model accounts for the pass time, bleeding, and the environmental evaporation rate to estimate the surface moisture and predict the lack of interlayer adhesion. In this study, the model accurately (RMSE = 2.5%) predicted an interlayer adhesion reduction from 30% to 50%. The interlayer adhesion results of other studies could also be predicted.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High reactivity of mesoporous CeO 2 to dissociate chemical warfare agent sarin

We perform spectroscopic surface characterization on mesoporous ceria (CeO 2 ) under exposure of sarin gas, also known as GB. We find mesoporous CeO 2 displays extraordinary reactivity towards GB dissociation under high vacuum. The spectroscopic results combined with density functional theory (DFT) calculations suggest that the cleavage of the P–F and P–OR bonds are both possible pathways for GB to dissociate on the CeO 2 surfaces. The modeling reveals differences in dissociation pathways of GB on the crystallographic (111)- and (110)-surfaces of ceria. Here we also report some dissociation of GB on mesoporous CeO 2 under ambient conditions. The finding that mesoporous CeO 2 can degrade GB under mild conditions is encouraging and implies the strong potential of ceria-based materials to be applied in combating chemical warfare agents related to the GB series.

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

A model for oblique impacts on material surfaces

Many practical situations of material damage, wear, and erosion involve collisions between small particles and surfaces at inclined angles. While there are many well-validated models of normal incidence impact situations, elastic-plastic models for oblique incidence impact events are lacking. Here the finite element method is used to predict the normal and tangential coefficient of restitution in oblique impacts for hard, elastic spheres impacting an elastic-perfectly plastic material surface. The proposed model covers various impact angles ranging from 0° to 45°, within a limiting impact velocity below which the effects of heating are negligible. The normal coefficient of restitution follows power-laws with respect to normalized values of the impact velocity. Interestingly, the tangential coefficient of restitution follows a linear relationship with impact velocity. Together, these results provide a semi-empirical set of equations predicting oblique impact rebounds (both velocity and trajectory) for a wide range of conditions and material properties, with which experimental results can be rapidly interpreted. As a result, Laser-Induced Particle Impact Test (LIPIT) data are also presented for aluminum particles impacting aluminum substrates, at impact angles of 25° and 40°; the results compare favorably with the model and validate the general use of such models for the analysis of experimental data.

Coefficient of restitution↗
DOE OSTI2026

Validated Reactive Force Field Quantifies MXene Interfacial Properties, Mechanics, and Thermal Transport

MXenes combine rich surface chemistry, mechanical strength, and high conductivity for a multitude of emerging applications. Predictive modeling supports accelerated materials designs and has been limited by the absence of validated and transferable force fields. Here, we introduce an interpretable, reactive INTERFACE force field (IFF and IFF-R) for Ti 3 C 2 T x MXenes that is trained based on chemical knowledge and achieves quantitative agreement with experiments across lattice parameters (<0.5%), density (<0.2%), liquid contact angles, Raman spectra, and the in-plane elastic modulus (∼320 GPa). The models cover surface terminations from hydroxyl (−OH) to fluorine (−F) groups and are extensible to other chemistries. We introduce pH-resolved surface chemistry and identify dopamine adsorption mechanisms at MXene–aqueous interfaces supported by QCM-D and UV–Vis experiments. The data reveal coplanar and perpendicular binding modes and concentration-dependent multilayer assembly. We predict previously inaccessible properties, including termination-dependent cleavage energies, interlayer shear moduli and dynamic shear failure, nanoindentation and brittle fracture, anisotropic in-plane and out-of-plane thermal conductivities, including the role of defects. Agreement with available experimental data is consistently close and exceeds DFT accuracy across the benchmark properties examined. The IFF/IFF-R model is compatible with CHARMM, AMBER, OPLS, and CVFF force fields for simulations of MXenes with diverse surface terminations, electrolyte interfaces, biointerfaces, and polymer composites without additional parameters. Parameter sets, 3D models, and analysis scripts are provided for community use. The validated, reactive, and transferable IFF framework facilitates predictive design of MXene-based films, membranes, sensing interfaces, and composites.

MXene↗
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↗