skip to main content


Title: Optimal Wetting Angles in Lattice Boltzmann Simulations of Viscous Fingering
Abstract We conduct pore-scale simulations of two-phase flow using the 2D Rothman–Keller colour gradient lattice Boltzmann method to study the effect of wettability on saturation at breakthrough (sweep) when the injected fluid first passes through the right boundary of the model. We performed a suite of 189 simulations in which a “red” fluid is injected at the left side of a 2D porous model that is initially saturated with a “blue” fluid spanning viscosity ratios $$M = \nu _\mathrm{r}/\nu _\mathrm{b} \in [0.001,100]$$ M = ν r / ν b ∈ [ 0.001 , 100 ] and wetting angles $$\theta _\mathrm{w} \in [ 0^\circ ,180^\circ ]$$ θ w ∈ [ 0 ∘ , 180 ∘ ] . As expected, at low-viscosity ratios $$M=\nu _\mathrm{r}/\nu _\mathrm{b} \ll 1$$ M = ν r / ν b ≪ 1 we observe viscous fingering in which narrow tendrils of the red fluid span the model, and for high-viscosity ratios $$M \gg 1$$ M ≫ 1 , we observe stable displacement. The viscous finger morphology is affected by the wetting angle with a tendency for more rounded fingers when the injected fluid is wetting. However, rather than the expected result of increased saturation with increasing wettability, we observe a complex saturation landscape at breakthrough as a function of viscosity ratio and wetting angle that contains hills and valleys with specific wetting angles at given viscosity ratios that maximize sweep. This unexpected result that sweep does not necessarily increase with wettability has major implications to enhanced oil recovery and suggests that the dynamics of multiphase flow in porous media has a complex relationship with the geometry of the medium and the hydrodynamical parameters.  more » « less
Award ID(s):
1918126
NSF-PAR ID:
10253872
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Transport in Porous Media
Volume:
136
Issue:
3
ISSN:
0169-3913
Page Range / eLocation ID:
831 to 842
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present simulations of two-phase flow using the Rothman and Keller colour gradient Lattice Boltzmann method to study viscous fingering when a “red fluid” invades a porous model initially filled with a “blue” fluid with different viscosity. We conducted eleven suites of 81 numerical experiments totalling 891 simulations, where each suite had a different random realization of the porous model and spanned viscosity ratios in the range$$M\in [0.01,100]$$M[0.01,100]and wetting angles in the range$$\theta _w\in [180^\circ ,0^\circ ]$$θw[180,0]to allow us to study the effect of these parameters on the fluid-displacement morphology and saturation at breakthrough (sweep). Although sweep often increased with wettability, this was not always so and the sweep phase space landscape, defined as the difference in saturation at a given wetting angle relative to saturation for the non-wetting case, had hills, ridges and valleys. At low viscosity ratios, flow at breakthrough is localized through narrow fingers that span the model. After breakthrough, the flow field continues to evolve and the saturation continues to increase albeit at a reduced rate, and eventually exceeds 90% for both non-wetting and wetting cases. The existence of a complicated sweep phase space at breakthrough, and continued post-breakthrough evolution suggests the hydrodynamics and sweep is a complicated function of wetting angle, viscosity ratio and time, which has major potential implications to Enhanced Oil Recovery by water flooding, and hence, on estimates of global oil reserves. Validation of these results via experiments is required to ensure they translate to field studies.

     
    more » « less
  2. Abstract

    Hemiwicking is the phenomena where a liquid wets a textured surface beyond its intrinsic wetting length due to capillary action and imbibition. In this work, we derive a simple analytical model for hemiwicking in micropillar arrays. The model is based on the combined effects of capillary action dictated by interfacial and intermolecular pressures gradients within the curved liquid meniscus and fluid drag from the pillars at ultra-low Reynolds numbers$${\boldsymbol{(}}{{\bf{10}}}^{{\boldsymbol{-}}{\bf{7}}}{\boldsymbol{\lesssim }}{\bf{Re}}{\boldsymbol{\lesssim }}{{\bf{10}}}^{{\boldsymbol{-}}{\bf{3}}}{\boldsymbol{)}}$$(107Re103). Fluid drag is conceptualized via a critical Reynolds number:$${\bf{Re}}{\boldsymbol{=}}\frac{{{\bf{v}}}_{{\bf{0}}}{{\bf{x}}}_{{\bf{0}}}}{{\boldsymbol{\nu }}}$$Re=v0x0ν, wherev0corresponds to the maximum wetting speed on a flat, dry surface andx0is the extension length of the liquid meniscus that drives the bulk fluid toward the adsorbed thin-film region. The model is validated with wicking experiments on different hemiwicking surfaces in conjunction withv0andx0measurements using Water$${\boldsymbol{(}}{{\bf{v}}}_{{\bf{0}}}{\boldsymbol{\approx }}{\bf{2}}\,{\bf{m}}{\boldsymbol{/}}{\bf{s}}{\boldsymbol{,}}\,{\bf{25}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{\lesssim }}{{\bf{x}}}_{{\bf{0}}}{\boldsymbol{\lesssim }}{\bf{28}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{)}}$$(v02m/s,25µmx028µm), viscous FC-70$${\boldsymbol{(}}{{\boldsymbol{v}}}_{{\bf{0}}}{\boldsymbol{\approx }}{\bf{0.3}}\,{\bf{m}}{\boldsymbol{/}}{\bf{s}}{\boldsymbol{,}}\,{\bf{18.6}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{\lesssim }}{{\boldsymbol{x}}}_{{\bf{0}}}{\boldsymbol{\lesssim }}{\bf{38.6}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{)}}$$(v00.3m/s,18.6µmx038.6µm)and lower viscosity Ethanol$${\boldsymbol{(}}{{\boldsymbol{v}}}_{{\bf{0}}}{\boldsymbol{\approx }}{\bf{1.2}}\,{\bf{m}}{\boldsymbol{/}}{\bf{s}}{\boldsymbol{,}}\,{\bf{11.8}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{\lesssim }}{{\bf{x}}}_{{\bf{0}}}{\boldsymbol{\lesssim }}{\bf{33.3}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{)}}$$(v01.2m/s,11.8µmx033.3µm).

     
    more » « less
  3. Abstract

    Two-dimensional electron systems subjected to high transverse magnetic fields can exhibit Fractional Quantum Hall Effects (FQHE). In the GaAs/AlGaAs 2D electron system, a double degeneracy of Landau levels due to electron-spin, is removed by a small Zeeman spin splitting,$$g \mu _B B$$gμBB, comparable to the correlation energy. Then, a change of the Zeeman splitting relative to the correlation energy can lead to a re-ordering between spin polarized, partially polarized, and unpolarized many body ground states at a constant filling factor. We show here that tuning the spin energy can produce fractionally quantized Hall effect transitions that include both a change in$$\nu$$νfor the$$R_{xx}$$Rxxminimum, e.g., from$$\nu = 11/7$$ν=11/7to$$\nu = 8/5$$ν=8/5, and a corresponding change in the$$R_{xy}$$Rxy, e.g., from$$R_{xy}/R_{K} = (11/7)^{-1}$$Rxy/RK=(11/7)-1to$$R_{xy}/R_{K} = (8/5)^{-1}$$Rxy/RK=(8/5)-1, with increasing tilt angle. Further, we exhibit a striking size dependence in the tilt angle interval for the vanishing of the$$\nu = 4/3$$ν=4/3and$$\nu = 7/5$$ν=7/5resistance minima, including “avoided crossing” type lineshape characteristics, and observable shifts of$$R_{xy}$$Rxyat the$$R_{xx}$$Rxxminima- the latter occurring for$$\nu = 4/3, 7/5$$ν=4/3,7/5and the 10/7. The results demonstrate both size dependence and the possibility, not just of competition between different spin polarized states at the same$$\nu$$νand$$R_{xy}$$Rxy, but also the tilt- or Zeeman-energy-dependent- crossover between distinct FQHE associated with different Hall resistances.

     
    more » « less
  4. McCartney, J.S. ; Tomac, I. (Ed.)
    Multiphase flow patterns in porous media largely depend on the properties of the fluids and interfaces such as viscosity, surface tension, and contact angle. Microorganisms in soils change the fluid and interfacial properties, and thus can alter multiphase fluid flow in porous media. This study investigates the impact of motile bacterium Escherichia coli ( E. coli ) on fluid displacement patterns in a microfluidic chip. The fluid displacement is observed during the saturation and the desaturation processes of the microfluidic chip with and without E.coli suspension. Time-lapse photography results show that the presence of E.coli alters the displacement patterns during the wetting and drying process and changes the residual saturation of the chip. Although studies of the impacts of motility on interfacial properties remain elusive, these results bring the expectation to the manipulation of multiphase transport in porous media and the adaptive control of industrial and environmental flow processes using active particles. 
    more » « less
  5. Abstract

    This paper reports the first integration of laser‐etched polycrystalline diamond microchannels with template‐fabricated microporous copper for extreme convective boiling in a composite heat sink for power electronics and energy conversion. Diamond offers the highest thermal conductivity near room temperature, and enables aggressive heat spreading along triangular channel walls with 1:1 aspect ratio. Conformally coated porous copper with thickness 25 µm and 5 µm pore size optimizes fluid and heat transport for convective boiling within the diamond channels. Data reported here include 1280 W cm−2of heat removal from 0.7 cm2surface area with temperature rise beyond fluid saturation less than 21 K, corresponding to 6.3 × 105W m−2K−1. This heat sink has the potential to dissipate much larger localized heat loads with small temperature nonuniformity (5 kW cm−2over 200 µm × 200 µm with <3 K temperature difference). A microfluidic manifold assures uniform distribution of liquid over the heat sink surface with negligible pumping power requirements (e.g., <1.4 × 10−4of the thermal power dissipated). This breakthrough integration of functional materials and the resulting experimental data set a very high bar for microfluidic heat removal.

     
    more » « less