SEARCH · Search NASA
Results for “PRESSURE DISTRIBUTION - AIRFOILS”
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.
Irrigation system
A gravity irrigation system includes a distribution piping having apertures to distribute water to a field, and a valve located upstream of the distribution piping. The valve limits a pressure of the water being delivered to the distribution piping. The system also includes a sump to receive the water at a lowest elevation of the field, a depth sensor disposed within the sump, and a return pump disposed at least partially within the sump to move the water to an elevated portion of the field. The system also includes a motor to drive the return pump, and a power source coupled to a variable frequency drive that powers the motor and controls a motor speed proportionately to an indication of the depth sensor. The system also includes a transfer piping to bring the water from the return pump to a check valve and to the distribution piping.
Gaseous hydrogen storage system with cryogenic supply
A gaseous hydrogen storage and distribution system with a cryogenic supply and a method for the cryogenic conversion of liquid hydrogen into high-pressure gaseous hydrogen are provided. The gaseous hydrogen storage and distribution system includes pressuring liquid hydrogen from a cryogenic tank using a low pressure liquid pump before vaporization within a relatively small vaporizer. The resulting high pressure gaseous hydrogen is transferred to a plurality of storage tanks at ambient temperature according to a desired fill sequence. The high pressure hydrogen gas is subsequently distributed from the storage tanks through a hydrogen fueling dispenser according to a desired dispensing sequence. The present system and method provide improvements in operational safety, eliminates the use of high pressure gas compressor, and minimizes boiling off and ventilation losses at a reduced cost when compared to existing thermal compression storage systems.
Impact of air distribution on dew point evaporative cooler thermal performance
A new geometry of the dew-point indirect evaporative cooler structure was proposed which enables uniform distribution of water in the working channels and application in traditional air handling units. This arrangement requires a complicated structure of the product channel in which the air is cooled. Further, the airflow direction changes and complicated structure it is associated with local losses affecting the total pressure loss. High-pressure drop limits the use of this dew point evaporative cooler as an alternative source of cooling. Therefore, a verified method based on numerical fluid dynamics (CFD) was used to determine the pressure drop and air distribution in the exchanger. The combination of the CFD method and modified epsilon-NTU allows the thermal performance of the device including the uniformity of air distribution, cooling capacity, and coefficient of performance to be determined. Finally, the appropriate exchanger dimensions, i.e., the channel height, or final distribution regulation may be found depending on the exchanger application.
Characterization of Surrogate Molten Salt Reactor Aerosol Streams
Measuring the aerosol evolution from MSRs is important for monitoring the off-gas system of the reactor and is particularly important for detecting off-normal conditions. In a molten salt reactor (MSR) accident scenario, an aerosol release would be a major factor in the source term. This aerosol stream would likely be generated from a breach in the cover gas system, which causes particles produced from fission itself to escape, or from a salt spill that produces aerosols through splashing and secondary reactions. The particle size of the produced aerosols is anticipated to vary greatly and range from 0.01 to 10 µm. The transport of these aerosols would be dependent on the particle size. A better understanding of aerosol generation, size, and monitoring methods are needed to inform estimation and mitigation of potential aerosolized source terms from MSRs. While salt spill experiments are being performed at Argonne National Laboratory, the development of aerosol characterization and monitoring methods are being developed at Oak Ridge National Laboratory. To generate prototypic aerosols for use in testing monitoring instruments and mitigation methods, a surrogate aerosol stream was produced with a Collison nebulizer, and the particle size distributions were measured with a cascade impactor. The results demonstrated that by changing the nebulizer pressure, the aerosol particle size distribution can be adjusted to best match the region of interest for experiments with higher pressures, driving the particle size down. However, nearly all aerosols formed exceeded 1 µm in diameter, providing a lower bound for the surrogate aerosol stream. In addition to verifying the applicability of this surrogate aerosol stream, this work has shown that a laser induced breakdown spectroscopy monitoring system that is under development is resilient to changes in particle sizes, increasing its robustness for off-gas monitoring.
Experimental and Computational Studies on Superhard Material Rhenium Diboride under Ultrahigh Pressures
An emerging class of superhard materials for extreme environment applications are compounds formed by heavy transition metals with light elements. In this work, ultrahigh pressure experiments on transition metal rhenium diboride (ReB 2 ) were carried out in a diamond anvil cell under isothermal and non-hydrostatic compression. Two independent high-pressure experiments were carried out on ReB 2 for the first time up to a pressure of 241 GPa (volume compression V/V 0 = 0.731 ± 0.004), with platinum as an internal pressure standard in X-ray diffraction studies. The hexagonal phase of ReB 2 was stable under highest pressure, and the anisotropy between the a-axis and c-axis compression increases with pressure to 241 GPa. The measured equation of state (EOS) above the yield stress of ReB 2 is well represented by the bulk modulus K 0 = 364 GPa and its first pressure derivative K 0 ´ = 3.53. Corresponding density-functional-theory (DFT) simulations of the EOS and elastic constants agreed well with the experimental data. DFT results indicated that ReB 2 becomes more ductile with enhanced tendency towards metallic bonding under compression. The DFT results also showed strong crystal anisotropy up to the maximum pressure under study. The pressure-enhanced electron density distribution along the Re and B bond direction renders the material highly incompressible along the c-axis. Our study helps to establish the fundamental basis for anisotropic compression of ReB 2 under ultrahigh pressures
The oxygen partial pressure in solid oxide electrolysis cells with multilayer electrolytes
Not provided.
TGA Study of the Reaction Kinetics of a CuFeMnAlO 4+δ Oxygen Carrier for Methane Chemical Looping Combustion under High Pressure and Atmospheric Pressure Conditions
Reduction and oxidation kinetics of a CuFeMnAlO 4+δ oxygen carrier for methane were studied using an atmospheric thermogravimetric analyzer (TGA) and a pressurized TGA. Important parameters such as temperature (500–900 °C), particle size distribution (150–650 μm), and partial pressures of reactants (methane and oxygen) were evaluated. Results showed that the reduction of the oxygen carrier was not affected by particle size distribution studied in this work, and that above 750 °C, the effect of the temperature on oxygen carrier reduction started to diminish. The reaction time approaches an asymptotic limit as methane partial pressure increases. For oxygen carrier oxidation, it was found that the rate of oxidation increased with increasing oxygen concentration. The effect of the oxygen concentration was more noticeable when the oxygen concentration was below 7.5%. Because of the inherent feature of the oxygen carrier, its chemical looping with oxygen uncoupling property was also assessed. A phase boundary-controlled shrinking sphere model was used to model the test data from the TGAs. Modeling results showed that the oxygen carrier reduction had a reaction order of 1 with an activation energy of ∼46 kJ/mol. Oxygen carrier oxidation demonstrated a reaction order of 1.1 with an activation energy of ∼55 kJ/mol. Additionally, the differences of the experimental results between the atmospheric pressure TGA and the pressurized TGA were discussed in detail and the limitations of the pressurized TGA were highlighted.
Theoretical and Phenomenological Studies of the Nucleon Structure in High Energy Reactions
The understanding of the internal structure of the proton and other strongly interacting particles is at the forefront of modern nuclear physics research. Generalized Parton Distribution Functions (GPDs) are a powerful tool to advance the understanding of hadron structure. In addition to the information about the one- dimensional collinear momentum distributions of partons (quarks, anti-quarks, and gluons) known from studies of high energy deep-inelastic reactions, GPDs also carry information on the distribution of partons in the transverse plane, and allow us in this way to access the three-dimensional structure of the nucleon. GPDs can be studied in hard-exclusive reactions and contain also information on the energy-momentum tensor form factors which will allow us to gain insights on quantities like pressure or angular momentum distribution inside the nucleon. The goal of this thesis is to deepen our understanding of the three-dimensional structure of the nucleon. We investigate energy-momentum tensor form factors and densities, and all leading-twist GPDs in the bag model. This quark model provides a consistent theoretical framework to investigate many general concepts that have recently attracted interest, and allows one to study insightful limits like the large-Nc limit, heavy-quark limit or the non-relativistic limit. Another important aspect of this thesis is the discussion of the monopole and quadrupole contributions to the angular momentum density. Finally, the description of pseudoscalar meson production in exclusive processes in terms of GPDs is discussed in the Goloskokov-Kroll model and implemented in the PARTONS framework, a software development project which will provide direct support for experiments at the Jefferson National Lab and the future Electron-Ion Collider.
Heat Transfer Experiments of a 1st Stage Blade Cascade for Supercritical CO2 Oxy-Combustion Turbine Application
The results of internally cooled 1st stage blade (S1B) cascade testing in a supercritical CO2 environment is presented. The turbine blade design has been previously established for the end application of an oxy-combustion turbine operating in the Allam-Fetvedt cycle with turbine inlet conditions of 305 bar and 1150°C. The internally cooled blade features leading edge (LE) region impingement cooling, mid-section ribbed serpentine passages, and a pin-finned trailing edge (TE) region before cooling ejection holes. The geometry for the tested blade cascade has a cooled central blade with un-cooled blades on either side to match flowpath areas of the actual turbine. The flowpath reuses internal components previously employed for mid-section region ribbed serpentine passage experiments that established Nusselt number enhancement ratios over a range of Reynolds numbers from 100,000-400,000. New components include flow conditioning plates upstream and downstream of the blade cascade to adequately represent the flow field and blade external heat transfer coefficient profiles for the actual turbine. The cooled central blade utilizes uniform crystal temperature sensors (UCTS) with six sensors each on the blade pressure and suction surfaces distributed radially and from LE to TE. The post-processed UCTS quantified the maximum wall temperature seen at each installed sensor location. The test procedure consisted of establishing supercritical CO2 cooling flow temperature and flow rate and maintaining it throughout the test. The flow rate aims to match that for the actual in-service turbine blade design and is maintained through an orifice restriction to keep the pressure differential between internal cooling flow and external hot flow nearly constant. For the sCO2 flow path external to the blade, temperatures were ramped throughout the test via control of the test loop’s natural gas burner heater. The maximum temperature seen was 468°C and held constant for a duration of 10 minutes at which the blade metal temperature was predicted to be at its maximum before ramping down. For the turbine blade design for service inlet conditions, external flow path computational fluid dynamics (CFD) results and an internal cooling 1-D thermal and hydraulic flow network model using experimentally validated correlations served as thermal finite element (FE) boundary conditions to predict blade metal temperatures. These predicted temperatures were subsequently utilized in a structural FE model to predict blade life ratings dictated by Haynes 282 creep strength data, having a strong dependence on temperature. The boundary conditions experienced during testing are used in the same workflow and compared to the experimental results, with the goal of validating the analysis methodology and providing insight on the uncertainty in local metal temperature predictions.
Modeling Oxygen Partial Pressure in Solid Oxide Electrolysis Cells: The Microstructure Effect
Oxygen partial pressure is an important thermodynamic state variable that affects both the performance and degradation of solid oxide electrolysis cells. In this work, a 3D model developed from Virkar’s 1D model has been applied to reconstructed and synthetic microstructures of Ni-YSZ-GDC-LSCF cell. The effect of the microstructures, including the thickness of the YSZ and GDC layer, and the compositions of the hydrogen and oxygen electrode, on the distribution of oxygen partial pressure has been investigated. The results show that the maximum oxygen partial pressure may occur on the interface between oxygen electrode and GDC layer or between GDC layer and YSZ layer depending on the rate of oxygen ion exchange between GDC and YSZ. A thicker GDC layer lowers the maximum oxygen partial pressure in the cell, while a thicker YSZ layer lowers the maximum oxygen partial pressure in the hydrogen electrode. In addition, the Ni:YSZ ratio and porosity also affects the maximum partial pressure. These findings provide insights on mitigating degradation in solid oxide electrolysis cell by tuning the microstructures.
Deformation of pores in response to uniaxial and hydrostatic stress cycling in Marcellus Shale: Implications for gas recovery
One of the main challenges during gas production from shale reservoirs is low recovery rate. One contributing factor to this outcome is an insufficient understanding of pore systems, especially pore behaviour following changes in reservoir conditions or resulting from gas production practices. Because the pressure in the producing well can be controlled, understanding the effects of pressure variation on the pore size distribution and methane trapping is necessary to help design optimal conditions to improve the gas recovery rate. This work is the first systematic study of sub-millimeter pore deformation in shale caused by uniaxial and hydrostatic stress up to 100 MPa. Overmature samples from the Middle Devonian Marcellus Shale were analyzed using neutron scattering (SANS and USANS) techniques to interpret the response of nanopores to stress cycling of magnitude and duration compatible with the hydraulic fracturing treatments. Experiments reported here are performed at a series of uniaxial pressures up to 100 MPa and at hydrostatic pressures of deuterated methane 0 and 50 MPa. Since at the original depth of the shale samples’ burial of 2184 m the hydrostatic pressure is approximately 27 MPa and the lithostatic pressure is about 55 MPa, the experimental conditions reasonably well simulate the reservoir pressure regime. Our SANS and USANS results show that different pore sizes are affected by uniaxial stress in different ways. Specifically, in the pore size range from 1 nm to 800 nm, a decrease of pore density with pressure is observed, with the most depleted being mesopores of about 100 nm in diameter. The observed decrease is likely related to deformation of kerogen, followed by a loss of pore nano-volume, as well as methane trapped in the micropores. For pores larger than 5 μm, USANS data suggest that the negative trend is reversed at above 74 MPa, and the number density of large macropores may increase with increased stress even above the original value. The increased number of macropores at high pressure may create new interconnected conduits for gas migration, resulting in a better recovery rate. Another important finding of this study is an irreversible rearrangement of pore size distribution taking place after pressure cycling. Furthermore, this irreversible reorganization of pore size distribution should be taken into account during management of well production to maximize recovery rate.
Ripple-Type Control for Enhancing Resilience of Networked Physical Systems
Distributed control agents have been advocated as an effective means for improving the resiliency of our physical infrastructures under unexpected events. Purely local control has been shown to be insufficient, centralized optimal resource allocation approaches can be slow. In this context, we put forth a hybrid low-communication saturation-driven protocol for the coordination of control agents that are distributed over a physical system and are allowed to communicate with peers over a “hotline” communication network. According to this protocol, agents act on local readings unless their control resources have been depleted, in which case they send a beacon for assistance to peer agents. Our ripple-type scheme triggers communication locally only for the agents with saturated resources and it is proved to converge. Moreover, under a monotonicity assumption on the underlying physical law coupling control outputs to inputs, the devised control is proved to converge to a configuration satisfying safe operational constraints. The assumption is shown to hold for voltage control in electric power systems and pressure control in water distribution networks. Numerical tests corroborate the efficacy of the novel scheme.
Ripple-Type Control for Enhancing Resilience of Networked Physical Systems: Preprint
Distributed control agents have been advocated as an effective means for improving the resiliency of our physical in-frastructures under unexpected events. While purely local control has been shown to be insufficient, centralized optimal resource allocation approaches can be slow. In this context, we put forth a hybrid low-communication saturation-driven protocol for the coordination of control agents that are distributed over a physicalsystem and are allowed to communicate with peers over a 'hot-line' communication network. According to this protocol, agents act upon on local readings unless their control resources have been depleted, in which case they send a beacon for assistance to peer agents. Our ripple-type scheme triggers communication locally only for the agents with saturated resources, and is proved to converge. Moreover, under a monotonicity assumption on the underlying physical law coupling control outputs to inputs, the devised control is proved to converge to a configuration satisfying safe operational constraints. The assumption is shown to hold for voltage control in electric power systems and pressure control in water distribution networks. Numerical tests on both networks corroborate the efficacy of the novel scheme.
Shedding light on shadow generalized parton distributions
The feasibility of extracting generalized parton distributions (GPDs) from deeply virtual Compton scattering (DVCS) data has recently been questioned because of the existence of an infinite set of so-called “shadow GPDs” (SGPDs). These SGPDs depend on the process and manifest as multiple solutions (at a fixed scale 𝑄 2 ) to the inverse problem that needs to be solved to infer GPDs from DVCS data. SGPDs therefore pose a significant challenge for extracting GPDs from DVCS data. With this motivation we study the extent to which QCD evolution can provide constraints on SGPDs. This is possible because the known classes of SGPDs begin to contribute to observables after evolution, and can then be constrained (at the input scale 𝑄$^{2}_{0}$) by data that has a finite 𝑄 2 range. The impact that SGPDs could have on determining the total angular momentum, pressure and sheer force distributions, and tomography is also discussed. We find that 𝑄 2 evolution, given empirical data over a wide range of skewness (𝜉 ) and 𝑄 2 scales, constrains the SGPDs studied at low 𝑥 and low 𝜉. If this finding generalizes to all SGPDs—and we have not been able to find any SGPDs with qualitatively different behavior—then this potentially makes extraction of GPDs from DVCS possible over a limited range of 𝑥 and 𝜉.
A Novel Technique for Distributed Measuring of Strain in the Vadose Zone
Characterizing water content and pressure changes in the vadose zone is important to understanding a variety of geologic processes, ranging from permeability, to evapotranspiration, and aquifer recharge. Changes in water content, or pressure, cause strain in the solid porous medium of the vadose zone, and it is possible that measuring those strains could be used a characterization tool. The Coherence-length-gated Microwave Photonics Interferometry (CMPI) technique measures strain at high resolution along many intervals defined by pairs of reflectors distributed along an optical fiber. This technique has recently been developed at Clemson University, and it has the spatial, and temporal resolution to characterize the strains that are expected to occur with hydrologic changes in the vadose zone. However, the technique has never been used to measure strain in porous media, so its capabilities remain uncertain. The objective of this thesis is to evaluate the ability of using CMPI to measure strain changes in the vadose zone. The research approach consists of conducting laboratory tests using a column filled with sand that was subjected to changes in water content and pressure. An optical fiber with CMPI reflectors was deployed in a high surface area ribbon and used to measure strain along the axis of the column. The column was made from 8-inch, Schedule 40 PVC pipe (75 cm tall, 20 cm inner diameter) and filled with medium-grained sand (K = 2.5x10-6 m/s, porosity = 0.28, Coefficient of Uniformity = 2, van Genuchten (n = 4.34 and = 0.00039 m-1), Young’s Modulus = 18 – 47 MPa). Five pressure and four temperature sensors attached to probes inserted into the wall. The optical fiber includes five reflectors spaced 10 cm apart on the inside and outside of the column along the vertical axis. Strain measurements are made between pairs of reflectors. The reflectors are created in 250-micron-diameter acrylate-coated single mode Corning SMF28e+ optical fiber using a femtosecond laser. The optical fiber was laminated between two pieces of polyester film creating a large surface area to transfer strain from the porous media to the optical fiber. The experiments were conducted by filling the column from the bottom or infiltrating water from the top. The water was allowed to equilibrate to room temperature prior to each test in order to limit thermoelastic strain. Five injection tests with a rate of 250 ml/min and three infiltration tests at varying rates were conducted, and the results show patterns of strain and pressure are generally similar. Hydrologic conditions define three zones based on the pressure magnitude and distribution. 1.) Ambient Zone where the pressure heads are quasi-static and the pressure gradient is roughly unity (head gradient of zero). This is the uppermost zone and is characterized by negative pressures. 2.) Transition Zone where the pressures increase from ambient to zero, are changing relatively rapidly and the pressure gradient is relatively steep (pressure head gradients of 2). The Transition zone is 10 to 15 cm thick. 3.) Positive Pressure Zone where the pressure is positive, the rate of change is slower than in the transition zone and the gradient is flatter (pressure head gradient 1.1 to 1.2). This is the lowest zone in the column. Injection of water causes the pressure to increase and the three zones to move upward. The Transition zone moves at a rate of approximately 0.0002 m/s +/- .00005, according to analyses of pressure profiles. This is the same as the average velocity of the water calculated as volumetric flux/effective porosity =5.4x10-5 m/s/0.28. Strain signals in the range of 10s were observed with a noise level of generally less than 0.1 (signal to noise ratio of greater than 100) during injection and drainage. The spatial and temporal distributions of strain caused by injection depend on the location of the moving pressure zones. The locations of the different pressure zones were determined from pressure profiles at different times and these data were transferred to strain time series. This showed that strain in the Ambient pressure zone is either unchanged or slightly compressive, whereas strain in Transition zone is tensile and roughly proportional to the pressure change, and strain ranges from tensile to compressive in the Positive Pressure zone. In general, the strain is variable at the top of the Positive Pressure zone, but it appears to be consistently tensile lower in the zone. The spatial distribution of strain causes three distinct stages in the strain times series measured between pairs of CMPI reflectors. The strain is unchanged or slightly compressive during Stage 1, it increases (positive strain is tensile) approximately linearly at a rate of 0.02 to 0.03 μ/s during the Stage 2, and then flattens out in Stage 3. The strain increases again during Stage 4 of the time series. The pressure also changes in stages. It is unchanged during Stage 1 and then increases relatively rapidly at a rate of between 4 and 10 Pa/s during Stage 2 and continues to increase during Stages 3 and 4, but at a rate slightly slower than during Stage 2 (from 2.5 to 3.5 Pa/s). The rates of pressure and strain change are consistent with basic analyses. The rate of pressurization is similar to the calculated velocity of flow in the saturated zone, and this is similar to the measured velocity of the strain increase. The ratio of the pressurization rate and strain rate during Stage 2 (2.5 (Pa/s) /0.03 (μ/s)) is 80 MPa, which is approximately equal to the upper range expected considering the uniaxial compression and the Young’s Modulus of the sand. An interesting effect occurs when the upper surface of the sand becomes saturated. Significant compression (several 10s of ) occurs as the pressure and saturation increase at the upper surface. Compression occurs throughout the column, but the effect is greatest at the top of the column. The rate of compression is fastest slightly before ponding occurs, but it slows markedly and nearly stops when water starts to accumulate at the surface (ponding). This effect reverses when the surface of the sand is drained, resulting in tension throughout the column. The compression caused by this effect can be as large or larger than the tensile strain that accumulated during filling. This effect was unexpected because increasing pressure is normally associated with tensile strain. Nevertheless, this effect was observed consistently in all tests where the pressure changes at the upper surface, including tests where water was injected from below or infiltrated from above. This effect behaves as if the pore pressure at the upper surface of the sand exerts a normal force on the boundary (increasing pore pressure exerts a downward compression on the boundary). The results of these laboratory experiments indicate that CMPI can measure strain caused by fluid pressure changes in the vadose zone with a signal to noise ratio of 100 or more. Injection and drainage cause a strain signal that is complex, but repeatable. The magnitudes of the strain signal are consistent with magnitudes that are expected based on poroelastic calculations using independently measured properties of the sand. These results indicate that the CMPI technique with an optical fiber laminated in a polyester ribbon generates data that represent the strain distribution during hydrologic processes of imbibition and drainage in variably saturated sand. This suggests that distributed strain measurements using CMPI could be a viable approach for evaluating processes in the vadose zone, laying the groundwork for future field implementation.
Formulation and Experimental Validation of an Agricultural Implement-Only MPR System for Maximum Compatibility with Existing Agricultural Tractors
Tightening emissions regulations and rising fuel costs have driven a desire across many industries for more efficient actuation systems. This is particularly true of the agricultural sector. An extremely common arrangement in this sector is the tractor and implement pairing, in which actuators on an implement are powered by a hydraulic supply system on the towing tractor. This arrangement complicates the development of energy efficient hydraulic systems, as many new system designs require modification of both machines to reap efficiency benefits. Past work by the authors’ team has demonstrated great potentials for Multi-Pressure-Rail (MPR) technology involving both the tractor and implement subsystems. However, applicability of this MPR technology in a more realistic scenario where only one vehicle is equipped with such technology was not addressed. This work proposes an implementation of the MPR technology to an agricultural planter that allows significant savings, while only modifying the implement machine. This is done by manipulating the load sense network of a stock tractor to set system pressures to those required by the MPR system. This greatly reduces the barrier to implementation of MPR technology in agriculture. The work begins by outlining the reference machine for the system, then reviews the MPR system working principle. After this, the proposed expansion to the MPR concept is laid out and applied to the reference system. Finally, experimental validation is carried out, demonstrating up to a 35% reduction in system power consumption when paired with a state of the art, double-LS System tractor, and 15% with a single-LS System tractor.
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.