Change In Contact Angle Research Articles

Replaceable Dielectric Film for Low-Voltage and High-Performance Electrowetting-Based Digital Microfluidics.

Electrowetting-on-dielectric (EWOD) technology has been considered as a promising candidate for digital microfluidic (DMF) applications due to its outstanding flexibility and integrability. The dielectric layer with a hydrophobic surface is the key element of an EWOD device, determining its driving voltage, reliability, and lifetime. Hereby, inspired by the ionic-liquid-filled structuring polymer with high capacitance independent on thickness, namely ion gel (IG), we develop a polymer (P)-ion gel-amorphous fluoropolymer, namely, PIGAF, composite film as a replaceable hydrophobic dielectric layer for fabrication of a high-efficiency and stable EWOD-DMF device at relatively low voltage. The results show that the proposed EWOD devices using the PIGAF-based dielectric layer can achieve a large contact angle (θ) change of ∼50° and excellent reversibility with a contact angle hysteresis of ≤5° at a relatively low voltage of 30 Vrms. More importantly, the EWOD actuation voltage did not change obviously with the PIGAF film thickness in the range of several to tens of microns, enabling the thickness of the film to be adjusted according to the demand within a certain range while keeping the actuation voltage low. An EWOD-DMF device can be prepared by simply stacking a PIGAF film onto a PCB board, demonstrating stable droplet actuation (motion) at 30 Vrms and 1 kHz as well as a maximum moving velocity of 69 mm/s at 140 Vrms and 1 kHz. The PIGAF film was highly stable and reliable, maintaining excellent EWOD performance after multiple droplet manipulations (≥50 cycles) or long-term storage of 1 year. The proposed EWOD-DMF device has been demonstrated for digital chemical reactions and biomedical sensing applications.

Nanobubble-Induced Aggregation of Ultrafine Particles: A Molecular Dynamics Study.

Nanobubble-induced aggregation (NBIA) of fine and ultrafine particles in liquid is a promising method for enhancing floatation rates in mineral processing, cleaning contaminants from water, and reviving marine ecosystems. Although the current experimental techniques can measure the nanobubble capillary force between two surfaces with controlled approach speed, they are not capable of imaging NBIA dynamics of fine/ultrafine particles by real-time observation with nanoscale spatial resolution. In this work, we use molecular dynamics (MD) simulations to study dynamics of NBIA of Ag particles in a Lennard-Jones fluid system. The molecular-level modeling allows us to study microscopic details of NBIA dynamics that are inaccessible by current experimental means. Using MD simulations, we investigated the effects of NB size, surface wettability, surface roughness, and contact line pinning on NBIA dynamics. Our modeling results show that both concave NB bridges between two hydrophobic surfaces and convex NB bridges between two hydrophilic surfaces can result in an attractive nanobubble capillary force (NBCF) that causes the aggregation of Ag particles in liquids. The equilibrium separation between two fully aggregated particles can be well predicted by the improved capillary force model. We also observe that the change of contact angle after the contact line pinning occurs at the sharp edge of a particle, which slows the aggregation process. Our thermodynamics analysis shows that there is a critical contact angle below which the merged surface NBs will detach from the surface instead of causing aggregation. The prediction of the critical contact angle is corroborated by our MD simulation results.

Numerical study of thermocapillary migration of a droplet on an oleophilic track

The thermocapillary migration dynamics of a droplet with a non-circular footprint on an oleophilic track is studied numerically. Three-dimensional simulations showed that reducing the width of the oleophilic track improves the droplet migration velocity only upto a specific track width when compared to the change in velocity of a droplet with circular footprint. Further decrease in track width is seen to depreciate the enhancement in migration velocity. The spread length of the droplet on the oleophilic track is also observed to have a similar variation with track width. The two trends are related because the temperature difference in the droplet, which scales with spread length, determines the driving force governing the migration velocity characteristics. Further, a parametric study is performed to assess the migration behaviour on oleophilic tracks with change in substrate temperature gradient, droplet volume and contact angles within and outside the track. The influence of change in these parameters on the variation in enhancement in spread length and corresponding track width is also investigated. The study identified trends in the minimum track width that may be used to improve migration characteristics, as well as the order of enhancement that can be noticed with changes in each parameter.

Host-rock and caprock wettability during hydrogen drainage: Implications of hydrogen subsurface storage

Underground hydrogen storage (UHS) is an integral part of H2 economy value chain, essential for industrial-scale actualization of global decarbonization objectives. UHS in depleted hydrocarbon reservoirs is considered a safer and promising storage technique due to presence of large porous formation and impermeable seal, but its effectiveness depends on precise estimation of rock wettability, a crucial parameter in reservoir characterization, which controls pore-scale gas distribution, H2 containment safety and withdrawal efficiency. Several recent studies have provided contact angles data for host rock-water-H2 (quartz and sandstone) and caprock-water-H2 (mica and shale) measured through sessile drop method. However, the contact angle datasets for carbonate rock-water-H2 measured via captive bubble method, which can reflect the wettability of the rock during imbibition and drainage are largely unknown. The present work characterized the wettability of carbonate-water-H2 systems for various rock types, prepared from five different lithologies with different mineral assemblages. The contact angle θ was measured using the captive bubble method at two different temperatures of 303 K and 348 K, and three different pressures (3.44, 10.34, and 17.23 MPa). Experimental results showed that all rocks remained intrinsically strongly water-wet (θ ranged between 17° – 30°) at all experimental conditions. Furthermore, no significant change in contact angles occurred with changing temperature and pressure. For instance, at 17.23 MPa, contact angle of H2/brine on anhydrite (S-4) rock were measured as 19° and 20° at 303 K and 348 K respectively, suggesting that H2 remains the non-wetting phase with increasing storage depth and at warmer reservoirs. The study suggests that the pore-scale flow and fluid connectivity of H2 may not be influenced by changing wettability in pure storage rocks/caprocks. Thus, the wettability of carbonates storage rocks and caprocks may not be over-predicted by assuming strongly water-wet conditions for the non-contaminated rock surfaces during UHS.

Study on mechanical behavior of modified wood chips reinforced cementitious composites and reinforcing mechanism of wood chips from the perspective of energy

ABSTRACT In this study, the mechanical behavior of modified wood chips reinforced cementitious composites was tested and compared with the composites prepared by ordinary(untreated) wood chips and the control. The effect of wood chips content was also evaluated. Then the reinforcing mechanism of modified wood chips was discussed through force analysis, contact angle test and surface energy calculation. And the reasons for the change of contact angle and surface energy were explained by FT-IR and SEM. The results showed that the compressive strength of modified wood chips reinforced cementitious composites decreased, whereas the flexural strength, flexural toughness, toughness index, fracture energy, and ductility index increased significantly compared to the control. With the increase in wood chips content, the compressive strength gradually decreased, the flexural strength increased at first and then decreased, and flexural toughness, toughness index, fracture energy, and ductility index gradually increased. Further, all the mechanical properties of modified wood chips reinforced cementitious composites were superior to those of ordinary wood chips reinforced cementitious composites. The reinforcing mechanism of modified wood chips is mainly attributed to the energy storage and dissipation behavior, and the modification improves the interfacial compatibility between the wood chips and cement matrix. Also, the combined modification makes the wood chips hydrophobic.

Fundamental study on the growth process of interfacial microbubbles in an air-supersaturated solution

Interfacial microbubbles widely exist and play important roles in nature and industries, while their growth process remained unknown in many aspects. In this study, the interfacial microbubbles are created in an air-supersaturated solution utilizing a sudden drop of ambient pressure, and the effects of different influential parameters on the growth process of interfacial microbubbles were studied by a self-made observation system. The investigations confirmed that the growth process of interfacial microbubbles under all conditions had experienced three stages, including the floating stage, transition stage and expansion stage. The experimental results also show that the magnitude of depressurization amplitude, surface wettability, surface roughness, dissolved air content, and solution surface tension have significant effects on the nucleation and growth of interfacial microbubbles. In addition, the modified diffusion theoretical model of dissolved air was successfully used to predict the changes in contact angle and bubble radius for interfacial bubbles during the growth process.

Effect of sodium hypochlorite as a depressant for copper species in Cu-Mo flotation separation

One alternative to replace conventional depressants in Cu-Mo separation is using oxidizing agents to promote the formation of products on chalcopyrite that render its surface hydrophilic and thus favor its depression. This work provides a study to understand how sodium hypochlorite can act as a depressor for chalcopyrite in Cu-Mo flotation separation. The effect of NaClO in flotation tests was evaluated by means of electrochemical techniques, scanning electron microscopy (SEM), Fourier transform infrared analysis and contact angle measurements. The flotation results showed a selective depressant effect at a dosage of 1,020 g/t NaClO, conditioning time of 5 min, and pH 7. The mixed potential showed that the oxidation of chalcopyrite and molybdenite is possible through the addition of NaClO. The SEM observations showed changes in oxidized chalcopyrite morphology, which is supported by the decreasing contact angle change; no morphological changes were observed on molybdenite. FTIR analyses showed the characteristic peaks of the OH group after chalcopyrite was treated with NaClO; when molybdenite was treated with NaClO, the FTIR spectrum did not show changes. A scheme is proposed for the formation of oxidation species.

Binary Mixture Droplet Evaporation on Microstructured Decorated Surfaces and the Mixed Stick-Slip Modes.

The interactions between liquid droplets and solid surfaces during wetting and phase change are important to many applications and are related to the physicochemical properties of the substrate and the fluid. In this work, we investigate experimentally the evaporation of pure water, pure ethanol, and their binary mixture droplets, accessing a wide range of surface tensions, on hydrophobic micro-pillared surfaces varying the spacing between the pillars. Results show that on structured surfaces, droplets evaporate following three classical evaporative behaviors: constant contact radius/pinning, stick-slip, or mixed mode. In addition, we report two further droplet evaporation modes, which are a mixed stick-slip mode where the contact angle increases while the contact radius decreases in a stick-slip fashion and a mixed stick-slip mode where both the contact angle and the contact radius decrease in a stick-slip fashion. We name these evaporation modes not yet reported in the literature as the increasing and decreasing contact angle mixed stick-slip modes, respectively. The former ensues because the fluid surface tension increases as the most volatile fluid evaporates coupled to the presence of structures, whereas the latter is due to the presence of structures for either fluid. The duration of each evaporation mode is dissimilar and depends on the surface tension and on the spacing between structures. Pure water yields longer initial pinning times on all surfaces before stick-slip ensues, whereas for binary mixtures and pure ethanol, initial pinning ensues mainly on short spacing structures due to the different wetting regimes displayed. Meanwhile, mixed stick-slip modes ensue mainly for high ethanol concentrations and/or pure ethanol independent of the solid fraction and for low ethanol concentrations on large spacing. Contact line jumps, changes in contact angle and pinning forces are also presented and discussed. This investigation provides guidelines for tailoring the evaporation of a wide range of surface tension fluids on structured surfaces for inkjet printing, DNA patterning, or microfluidics applications.

Prediction of electrostatic properties of reservoir rock in low salinity water injection into carbonate reservoirs

Geochemical modeling along with chemical reactions is one of the challenges in modeling of low salinity water injection. The most important issue in the geochemical model is to determine the correct electrical charge distribution model and its tuning parameters. The composition of the rock as well as the candidate water used is effective in determining the type of model and its parameters, so that the tuning parameters are determined based on the history of zeta potential experiments. In this study, in order to determine the correct model of electrical charge distribution and its tuning parameters in carbonate rock samples, first, equilibrium samples of Candidate water with crushed rock are subjected to static zeta potential tests. Then, the diffuse electrical double layer model is used to determine the electrical charge of the rock/water and water/oil surfaces and to predict the zeta potential. In the following, by adjusting the tuning parameters of the model to match the prediction results of the model with the history of the laboratory data, the density of the carbonate rock surface, the equilibrium constants and the kinetics of the governing reactions are determined. The obtained results show that the range of error in zeta potential prediction by the model compared to the laboratory data is from 2 to 20%, which is within the acceptable range of the performance of electrical charge distribution models. Moreover, it could be observed that the error of prediction using DLM model is significantly less than the conventional models (CD-MUSIC and BSM) for different candidate water. Finally, the effect of calculated zeta potential changes is used to calculate the contact angle changes of low salinity water injection based on the coupling of DLVO theory and geochemical model. The results of the study prove that the prediction error is less than 5% compared to the results of the static wettability tests. Based on this, according to the good match between the model and the laboratory results, it is possible to determine the properties of surface sites in surface complexation models of carbonate samples using the proposed approach and the subsequent tuning data of the geochemical model.

Durable ultra-high hydrophobic oil absorbing materials based on biomass prepared via a simple low-cost strategy

Wood is renewable and biodegradable natural resource which has been widely used in construction industry and home furnishing; however, it absorbs moisture easily limiting its application. Herein, an ultra-hydrophobic wood was developed with acetylated β-cyclodextrin (Ac-β-CD). Due to the unique hydrophobic cavity structure of Ac-β-CD, a special micro/nanostructure was constructed on the wood cell wall achieving a maximum water contact angle of 150.5°. Also, it was stable at 137°, with a contact angle change rate of 1.16% during 10–100 s. The tensile strength of Ac-β-CD-wood increased by 58.76% due to the hydrogen bond network across the wood, and it also showed good performance in oil absorption capacity. Moreover, Ac-β-CD-wood can withstand more than 100 wipes while retaining its hydrophobic property. Therefore, this study provides a new eco-friendly strategy for producing biodegradable hydrophobic oil-absorbing wood with outstanding tensile strength and dimensional stability.

Improvement in Switching Characteristics and Bias Stability of Solution-Processed Zinc-Tin Oxide Thin Film Transistors via Simple Low-Pressure Thermal Annealing Treatment.

In this study, we used a low-pressure thermal annealing (LPTA) treatment to improve the switching characteristics and bias stability of zinc-tin oxide (ZTO) thin film transistors (TFTs). For this, we first fabricated the TFT and then applied the LPTA treatment at temperatures of 80 °C and 140 °C. The LPTA treatment reduced the number of defects in the bulk and interface of the ZTO TFTs. In addition, the changes in the water contact angle on the ZTO TFT surface indicated that the LPTA treatment reduced the surface defects. Hydrophobicity suppressed the off-current and instability under negative bias stress because of the limited absorption of moisture on the oxide surface. Moreover, the ratio of metal-oxygen bonds increased, while the ratio of oxygen-hydrogen bonds decreased. The reduced action of hydrogen as a shallow donor induced improvements in the on/off ratio (from 5.5 × 103 to 1.1 × 107) and subthreshold swing (8.63 to V·dec-1 and 0.73 V·dec-1), producing ZTO TFTs with excellent switching characteristics. In addition, device-to-device uniformity was significantly improved because of the reduced defects in the LPTA-treated ZTO TFTs.

Mechanism for the effects of surface chemical composition and crystal face on the wettability of α-quartz surface

Quartz is a ubiquitous matter in natural rock and industrial production, and the wetting property of quartz surface plays a critical role in determining numerous industrial processes. However, the key mechanism for the effects of surface chemical composition (based on surface pretreatment) and crystal face on the wettability of quartz surface remains vague. Herein, by applying an advanced atomic force microscopy (AFM), we measured the contact angles of liquids on three common crystal faces of α-quartz after different pretreatments. The time-of-flight secondary ion mass spectrometry (ToF-SIMS) analysis and molecular dynamics (MD) simulations were further conducted to uncover the microscopic mechanism involved. It is found that the surface pretreatment does significantly affect the wettability of α-quartz faces, mainly because of the variation in surface hydroxyl concentration after different surface pretreatments. The magnitude of change in contact angle after distinct pretreatments is also affected by the liquid property. After certain pretreatments, although there is no difference in the wettability of different faces, the inherent atomic arrangement of crystal faces still causes discrepancies in other interface properties, e.g., liquid structure, diffusion coefficient of liquid molecule. The reasonable pretreatments of α-quartz for a range of applications are recommended based on the current findings.

Non-invasive rust detection of steel plates determined through interfacial modulus

Initial methods to detect rust in pipelines have been conducted through invasive probes and sectioning off parts of the facility as the plant is running. These methods greatly increase the costs overall. The need for a feasible solution to this issue lies in the detection of rust formation through a non-invasive method. This study’s objective is to measure rust formation through droplet motion on the outer layer of pipelines. Multiple experiments are conducted using carbon steel sheets whose bottom layer has been exposed to acid for different durations of time. As rust formation in the metal is a voltaic phenomenon, it would mean that the acid corrosion of the bottom layer would adversely affect the top layer of the substrate. Consequentially, droplet motion and the droplet’s contour would change in different corrosive scenarios which we could then detect with novel parameters in our lab. One such parameter is the Interfacial Modulus (GS), which describes the initial resistance of the solid’s outer layer towards the liquid. We can understand this parameter with the aid of the novel device, known as the Centrifugal Adhesion Balance (CAB). As we cause the drop to slide across the substrate at constant normal force condition, we observe the difference in the contour of the drop as it slides across the substrate. The real-time change in contact angles at each edge of the drop, along with its change in external lateral force, causes a change in the GS values, which varies in different corrosive scenarios.

Effects of the Surface Roughness of Six Wood Species for Furniture Production on the Wettability and Bonding Quality of Coating

Wood surface roughness, surface free energy (SFE), wettability, and bonding quality for water-based acrylic coatings were investigated. The samples tested in this study included Pinus radiata, Pinus sylvestris, Larch, Hemp oak, Catalpa tree, and Camphor. Sandpaper with grits of 180, 240, 320, 400, and 500 was utilized to sand wood surfaces. The van OSS-Chaudhury-Good equation (vOCG) was used to calculate the SFE values. The modified model (M-D) was used to calculate the wettability based on the contact angle change rate (K value). The higher the K value, the faster the contact angle approaches equilibrium. A cross-cut test was used to evaluate the coating’s bonding quality. The anatomical structure of wood has an impact on the roughness of hardwood. The equilibrium contact angle is influenced by the wood species and sandpaper grit size. Sanding can make the surface of wood more wettable. Radiata pine that had been sanded to 180 grit had the highest SFE value. After finishing with waterborne acrylic, hardwood had a slightly better coating adhesion than softwood. Hemp oak wood had the lowest coating adhesion (0.6) and the highest K value (0.82). The best bonding quality (0.4) was supplied by the camphor wood with the lowest K value (0.13). Wettability in terms of K values was a good indication of determining the bonding quality of the water-based acrylic coatings.

Simulation and neural network analysis of spreading the oblique fuel jet impinging onto concave walls

A series of numerical simulations of oblique-jet impinging onto concave walls are carried out to gain insight into the formation process of the liquid fuel film in a liquid film cooling thrust chamber. Moreover, the influence of jet and wall surface parameters on spreading the fuel film is investigated. The volume-of-fluid method is employed to capture the spreading process of liquid film formation after an oblique jet impacts the curved surface. Furthermore, the neural network model is developed to analyze the effects of jet diameter (0.2 mm ≤ d ≤ 0.6 mm), jet velocity (15 m/s ≤ v ≤ 20 m/s), impingement angle (15°≤β ≤ 30°), equilibrium contact angle (45°≤θe≤75°), curvature radius (22.5 mm ≤ R ≤ 67.5 mm), and surface roughness (3.2 μm ≤ Ra ≤ 16μm) on the fuel film's spreading length and maximum width based on Radial Basis Function. The results show that the liquid fuel film gradually stabilizes from upstream to downstream during spreading and the circumferential spreading of the liquid film is at its greatest when all the forces in the circumferential direction are in equilibrium. The jet diameter has the most significant effect on film spreading among jet parameters. Doubling the jet diameter at least doubles the film width and length. The effect of jet angle and velocity on film spreading is also magnified at larger jet diameters. Compared to other surface parameters, a change in the wall contact angle has the most marked effect on the spreading of the film. Lastly, the wall radius of curvature between 40 mm and 50 mm is more favorable for obtaining larger film lengths.

Multi-treatments based on polydimethylsiloxane and metal-organic framework wrapped with graphene oxide for achieving long-term corrosion and fouling protection: experimental and density functional theory aspects

In this study, ZIF-8 metal–organic framework decorated on surface of graphene oxide nanosheets and functionalized with silane moieties (F-GO@MOF) is employed to enhance the anti-corrosion, anti-bacterial, and anti-fouling properties of epoxy coatings modified with polydimethylsiloxane (epoxy-PDMS). The epoxy-PDMS/F-GO@MOF sample has revealed super-hydrophobic behavior with least water contact angle changes (Δθ ≈ 5° while Δθ ≈ 17° is recorded for epoxy-PDMS sample) during 20 days of immersion, twice pull-off adhesion strength (6.8 MPa compared to 3.7 MPa for epoxy-PDMS coating), and significantly higher impedance modulus (8.82 × 108 Ω·cm2 which is 3 times more than the epoxy-PDMS sample) during 60 days of immersion in NaCl electrolyte. The anti-bacterial test against Pseudomonas sp. and anti-fouling monitoring in Qingdao port have shown that fewer bacteria and marine organisms are attached on the epoxy-PDMS/F-GO@MOF coating, respectively, indicating advanced fouling and bacterial resistance of the nanocomposite sample. Based on the theoretical studies by density functional theory (DFT) method, the binding of the components in the GO@MOF complex is notably large which could inhibit the penetration of corrosive agents toward the steel substrate. In addition, the stabilization of the GO@MOF nanohybrid is mainly controlled through van der Waals (vdW) forces.

Experimental Investigation of Tube-in-Tube Nanocomposite Coated Heat Exchanger

In this paper, heat transfer performance is investigated for plain and nanoparticle coated tube-in-tube heat exchangers. Four types of tubes, i.e. bare copper tube, bare aluminium tube, Cu-Al2O3 nanoparticle coated tube and anodized aluminium tubes are used for performing the experimental investigation. The coating thickness of Cu-Al2O3 nanocomposite surface and anodized aluminium tube varies from 10, 25 and 30 micrometres to 15, 20 and 30 micrometres respectively. The surface is found to be hydrophilic in nature as the contact angle changes from 79.82deg to 55.47deg. The prepared surfaces are characterized by FESEM, EDS and FTIR. By adjusting the hot and cold fluids relative mass flow rates, one may calculate not only the overall heat transfer coefficient but also the efficiency of the tube-in-tube heat exchanger. The Cu-Al2O3 nanocomposite coated surface has the highest overall heat transfer coefficient and efficiency, followed by an anodized aluminium surface, bare copper surface and bare aluminium surface. The aluminium anodic oxide (AAO) surface also exhibits an increased heat transfer coefficient, but to a lesser extent than the nanocomposite coating.

Arylazopyrazole-Modified Thiolactone Acrylate Copolymer Brushes for Tuneable and Photoresponsive Wettability of Glass Surfaces.

Photoswitches have long been employed in coatings for surfaces and substrates to harness light as a versatile stimulus to induce responsive behavior. We previously demonstrated the viability of arylazopyrazole (AAP) as a photoswitch in self-assembled monolayers (SAMs) on silicon and glass surfaces for photoresponsive wetting applications. We now aim to transfer the excellent photophysical properties of AAPs to polymer brush coatings. Compared to SAMs, polymer brushes offer enhanced stability and an increase of the thickness and density of the functional organic layer. In this work, we present thiolactone acrylate copolymer brushes which can be post-modified with AAP amines as well as hydrophobic acrylates, making use of the unique chemistry of the thiolactones. This strategy enables photoresponsive wetting with a tuneable range of contact angle change on glass substrates. We show the successful synthesis of thiolactone hydroxyethyl acrylate copolymer brushes by means of surface-initiated atom-transfer radical polymerization with the option to either prepare homogeneous brushes or to prepare micrometer-sized brush patterns by microcontact printing. The polymer brushes were analyzed by atomic force microscopy, time-of-flight secondary ion spectrometry, and X-ray photoelectron spectroscopy. Photoresponsive behavior imparted to the brushes by means of post-modification with AAP is monitored by UV/vis spectroscopy, and wetting behavior of homogeneous brushes is measured by static and dynamic contact angle measurements. The brushes show an average change in static contact angle of around 13° between E and Z isomer of the AAP photoswitch for at least five cycles, while the range of contact angle change can be fine-tuned between 53.5°/66.5° (E/Z) and 81.5°/94.8° (E/Z) by post-modification with hydrophobic acrylates.

Automatic measurement of three-phase contact angles in pore throats based on digital images

With the help of digital image processing technology, an automatic measurement method for the three-phase contact angles in the pore throats of the microfluidic model was established using the microfluidic water flooding experiment videos as the data source. The results of the new method were verified through comparing with the manual measurement data. On this basis, the dynamic changes of the three-phase contact angles under flow conditions were clarified by the contact angles probability density curve and mean value change curve. The results show that, for water-wetting rocks, the mean value of the contact angles is acute angle during the early stage of the water flooding process, and it increases with the displacement time and becomes obtuse angle in the middle-late stage of displacement as the dominant force of oil phase gradually changes from viscous force to capillary force. The droplet flow in the remaining oil occurs in the central part of the pore throats, without three-phase contact angle. The contact angles for the porous flow and the columnar flow change slightly during the displacement and present as obtuse angles in view of mean values, which makes the remaining oil poorly movable and thus hard to be recovered. The mean value of the contact angle for the cluster flow tends to increase in the flooding process, which makes the remaining oil more difficult to be recovered. The contact angles for the membrane flow are mainly obtuse angles and reach the highest mean value in the late stage of displacement, which makes the remaining oil most difficult to be recovered. After displacement, the remaining oils under different flow regimes are just subjected to capillary force, with obtuse contact angles, and the wettability of the pore throat walls in the microfluidic model tends to be oil-wet under the action of crude oil.

Numerical Experiments for Surfactant Infiltration in the Vadose Zone to Demonstrate Concentration-Dependent Capillarity, Viscosity, and Sorption Characteristics

Surfactants (i.e., solutes that reduce the surface tension of water) exist in the subsurface either naturally or are introduced to the subsurface due to anthropogenic activities (e.g., agricultural purposes, environmental remediation strategies). Surfactant-induced changes in surface tension, contact angle, density, and viscosity alter the water retention and conduction properties of the vadose zone. This research numerically investigates the effects of surfactants in the vadose zone by comparing the flow and transport of three different surfactant solutions, namely butanol, ethanol, and Triton X-100. For each surfactant case, surfactant-specific concentration-dependent surface tension, contact angle, density, and viscosity relationships were incorporated by modifying a finite element unsaturated flow and transport code. The modified code was used to simulate surfactant infiltration in the vadose zone at residual state under intermittent boundary conditions. The modelling results show that all three surfactant solutions led to unique and noteworthy differences in comparison to the infiltration of pure water containing a conservative tracer. Results indicate that surfactant infiltrations led to complex patterns with reduced vertical movement and enhanced horizontal spreading, which are a function of concentration-dependent surface tension, density, contact angle, viscosity and sorption characteristics. The findings of this research will help understanding the effects of surfactant presence in the subsurface on unsaturated flow and its possible links to future environmental problems.