skip to main content


Title: Influence of Wetting on Viscous Fingering Via 2D Lattice Boltzmann Simulations
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
Award ID(s):
1918126
NSF-PAR ID:
10234482
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Transport in Porous Media
Volume:
138
Issue:
3
ISSN:
0169-3913
Format(s):
Medium: X Size: p. 511-538
Size(s):
p. 511-538
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The elliptic flow$$(v_2)$$(v2)of$${\textrm{D}}^{0}$$D0mesons from beauty-hadron decays (non-prompt$${\textrm{D}}^{0})$$D0)was measured in midcentral (30–50%) Pb–Pb collisions at a centre-of-mass energy per nucleon pair$$\sqrt{s_{\textrm{NN}}} = 5.02$$sNN=5.02 TeV with the ALICE detector at the LHC. The$${\textrm{D}}^{0}$$D0mesons were reconstructed at midrapidity$$(|y|<0.8)$$(|y|<0.8)from their hadronic decay$$\mathrm {D^0 \rightarrow K^-\uppi ^+}$$D0K-π+, in the transverse momentum interval$$2< p_{\textrm{T}} < 12$$2<pT<12 GeV/c. The result indicates a positive$$v_2$$v2for non-prompt$${{\textrm{D}}^{0}}$$D0mesons with a significance of 2.7$$\sigma $$σ. The non-prompt$${{\textrm{D}}^{0}}$$D0-meson$$v_2$$v2is lower than that of prompt non-strange D mesons with 3.2$$\sigma $$σsignificance in$$2< p_\textrm{T} < 8~\textrm{GeV}/c$$2<pT<8GeV/c, and compatible with the$$v_2$$v2of beauty-decay electrons. Theoretical calculations of beauty-quark transport in a hydrodynamically expanding medium describe the measurement within uncertainties.

     
    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

    We present a proof of concept for a spectrally selective thermal mid-IR source based on nanopatterned graphene (NPG) with a typical mobility of CVD-grown graphene (up to 3000$$\hbox {cm}^2\,\hbox {V}^{-1}\,\hbox {s}^{-1}$$cm2V-1s-1), ensuring scalability to large areas. For that, we solve the electrostatic problem of a conducting hyperboloid with an elliptical wormhole in the presence of anin-planeelectric field. The localized surface plasmons (LSPs) on the NPG sheet, partially hybridized with graphene phonons and surface phonons of the neighboring materials, allow for the control and tuning of the thermal emission spectrum in the wavelength regime from$$\lambda =3$$λ=3to 12$$\upmu$$μm by adjusting the size of and distance between the circular holes in a hexagonal or square lattice structure. Most importantly, the LSPs along with an optical cavity increase the emittance of graphene from about 2.3% for pristine graphene to 80% for NPG, thereby outperforming state-of-the-art pristine graphene light sources operating in the near-infrared by at least a factor of 100. According to our COMSOL calculations, a maximum emission power per area of$$11\times 10^3$$11×103W/$$\hbox {m}^2$$m2at$$T=2000$$T=2000K for a bias voltage of$$V=23$$V=23V is achieved by controlling the temperature of the hot electrons through the Joule heating. By generalizing Planck’s theory to any grey body and deriving the completely general nonlocal fluctuation-dissipation theorem with nonlocal response of surface plasmons in the random phase approximation, we show that the coherence length of the graphene plasmons and the thermally emitted photons can be as large as 13$$\upmu$$μm and 150$$\upmu$$μm, respectively, providing the opportunity to create phased arrays made of nanoantennas represented by the holes in NPG. The spatial phase variation of the coherence allows for beamsteering of the thermal emission in the range between$$12^\circ$$12and$$80^\circ$$80by tuning the Fermi energy between$$E_F=1.0$$EF=1.0eV and$$E_F=0.25$$EF=0.25eV through the gate voltage. Our analysis of the nonlocal hydrodynamic response leads to the conjecture that the diffusion length and viscosity in graphene are frequency-dependent. Using finite-difference time domain calculations, coupled mode theory, and RPA, we develop the model of a mid-IR light source based on NPG, which will pave the way to graphene-based optical mid-IR communication, mid-IR color displays, mid-IR spectroscopy, and virus detection.

     
    more » « less
  4. Abstract

    We propose a new observable for the measurement of the forward–backward asymmetry$$(A_{FB})$$(AFB)in Drell–Yan lepton production. At hadron colliders, the$$A_{FB}$$AFBdistribution is sensitive to both the electroweak (EW) fundamental parameter$$\sin ^{2} \theta _{W}$$sin2θW, the weak mixing angle, and the parton distribution functions (PDFs). Hence, the determination of$$\sin ^{2} \theta _{W}$$sin2θWand the updating of PDFs by directly using the same$$A_{FB}$$AFBspectrum are strongly correlated. This correlation would introduce large bias or uncertainty into both precise measurements of EW and PDF sectors. In this article, we show that the sensitivity of$$A_{FB}$$AFBon$$\sin ^{2} \theta _{W}$$sin2θWis dominated by its average value around theZpole region, while the shape (or gradient) of the$$A_{FB}$$AFBspectrum is insensitive to$$\sin ^{2} \theta _{W}$$sin2θWand contains important information on the PDF modeling. Accordingly, a new observable related to the gradient of the spectrum is introduced, and demonstrated to be able to significantly reduce the potential bias on the determination of$$\sin ^{2} \theta _{W}$$sin2θWwhen updating the PDFs using the same$$A_{FB}$$AFBdata.

     
    more » « less
  5. Abstract

    A well-known open problem of Meir and Moser asks if the squares of sidelength 1/nfor$$n\ge 2$$n2can be packed perfectly into a rectangle of area$$\sum _{n=2}^\infty n^{-2}=\pi ^2/6-1$$n=2n-2=π2/6-1. In this paper we show that for any$$1/21/2<t<1, and any$$n_0$$n0that is sufficiently large depending on t, the squares of sidelength$$n^{-t}$$n-tfor$$n\ge n_0$$nn0can be packed perfectly into a square of area$$\sum _{n=n_0}^\infty n^{-2t}$$n=n0n-2t. This was previously known (if one packs a rectangle instead of a square) for$$1/21/2<t2/3(in which case one can take$$n_0=1$$n0=1).

     
    more » « less