Ephemeral Antibubbles: Spatiotemporal Evolution from Direct Numerical Simulations Nairita Pal,1 Rashmi Ramadugu,2 Prasad Perlekar,2 and Rahul Pandit3 1Los Alamos National Laboratory, NM 87545, USA 2TIFR Center for Interdisciplinary Sciences, Tata Institute of Fundamental Research, Gopanpally, Hyderabad, 500046, India 3Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India (Dated: April 1, 2021) Antibubbles, which consist of a shell of a low-density fluid inside a high-density fluid, have several promising applications. We show, via extensive direct numerical simulations (DNSs), in both two and three dimensions (2D and 3D), that the spatiotemporal evolution of antibubbles can be described naturally by the coupled Cahn-Hilliard-Navier-Stokes (CHNS) equations for a binary fluid. Our DNSs capture elegantly the gravity-induced thinning and breakup of an antibubble via the time evolution of the Cahn-Hilliard scalar order parameter field φ, which varies continuously across interfaces, so we do not have to enforce complicated boundary conditions at the moving antibubble interfaces. To ensure that our results are robust, we supplement our CHNS simulations with sharp- interface Volume-of-Fluid (VoF) DNSs. We track the thickness of the antibubble and calculate the dependence of the lifetime of an antibubble on several parameters; we show that our DNS results agree with various experimental results; in particular, the velocity with which the arms of the antibubble retract after breakup scales as σ1=2, where σ is the surface tension. I. INTRODUCTION which have been used to study a variety of problems in binary-fluid flows [26{35], provide a minimal theo- Antibubbles, which comprise a shell of a low-density retical framework for studying the spatiotemporal evo- fluid inside a high-density fluid, have been known lution of antibubbles; in addition to a velocity field, for close to 90 years, since the work of Hughes and the CHNS system employs a phase field φ that distin- Hughes [1]. In contrast to an ordinary bubble, an an- guishes between the two fluid phases. We use the CHNS tibubble has two surfaces, which trap a certain volume equations [26, 27, 29{32] to study antibubbles in two- of fluid between them. Therefore, the contact area of an dimensions (2D) and in three dimensions (3D) by using antibubble is much larger than that of a bubble with the extensive DNSs. To complement our CHNS results, we same fluid volume; this property can be exploited for a also employ an alternative Volume-of-Fluid (VoF) numer- variety of chemical reactions. Furthermore, antibubbles ical scheme that is a sharp-interface method. are commonly used in clinical diagnostic imaging, sono- Our studies yield several interesting results that (a) poration (see, e.g., Ref. [2]), as agents for ultrasound- provide explanations for many experimental observations guided drug delivery [3{5], and for active leakage detec- and (b) suggest new experimental studies: Our results for tion [6]. Clearly, an antibubble is unstable in the pres- the spatiotemporal evolution of antibubbles, in 2D and ence of gravity, and, if the inner core of the antibubble 3D, show clearly how antibubble breakup occurs either is denser than the outer core, the antibubble rises under because of gravity-induced thinning or the puncturing of gravity; because of hydrostatic pressure in the outer core, its bottom boundary; the the collapsing antibubble then the fluid rises from the bottom to the top, resulting in a forms a rim that retracts with a velocity vrim. We un- thinning, and subsequent collapse, of the outer shell. Al- cover signatures of this collapse in Fourier-space spectra though there have been a number of experimental inves- of the Cahn-Hilliard field φ, the velocity, and the vortic- ity. We show that v σ1=2, where σ is the surface tigations of the spatiotemporal evolution of an antibub- rim ∼ ble [7{15] and drops [16{20] to name a few, over the past tension; and this power-law exponent is independent of few decades, theoretical studies of antibubble evolution the kinematic viscosity ν of the background fluid. We have been initiated only recently [15, 21{25] and they investigate the dependence of the scaled antibubble life- time τ1/τg on σ, ν, and the scaled outer radius of the have not, attempted, hitherto, to address the spatiotem- p poral evolution of antibubbles in detail. The number of antibubble R0=h0, in 2D; here, τg R0=Ag, with A arXiv:2103.16780v1 [physics.flu-dyn] 31 Mar 2021 ≡ experimental studies of the complete spatiotemporal evo- the Atwood number and g the acceleration due to grav- lution of antibubbles is also limited, partly because great ity, and h0 is the initial thickness of the antibubble shell. care has to be exercised to stabilise antibubbles; often, We compare our results with experiments and earlier the- surfactant molecules have to be introduced into the high- oretical studies. density liquid phase for such stabilization. The remainder of this paper is organized as follows: In We develop a natural, multiphase model for antibub- Sec. II we present the models and the numerical methods bles and demonstrate how we can use direct numerical that we use. In Sec. III we present the results our DNSs simulations (DNSs) to follow the spatiotemporal develop- in the following subsections: ment of ephemeral, but beautiful, antibubbles. We show • subsection III A on the spatiotemporal evolution of that the Cahn-Hilliard-Navier-Stokes (CHNS) equations, antibubbles (2D DNSs); 2 2 • subsection III B on the temporal evolution of the fluid-energy dissipation rate "(t) = ν !(x; t) x, h j j i Fourier spectra of φ and the velocity and vortic- where x denotes the average over space, and ! = u ity fields that are associated with antibubbles in is the fluidhi vorticity. We calculate the lifetime of anr× an- our 2D DNSs; tibubble falling under gravity in two ways: (a) from the energy time series; and (b) from the minimum thickness • subsection III C on the velocity of the retracting at the South pole of the antibubble. rim of the ruptured antibubble in our 2D DNSs; At time t = 0 we begin with the order-parameter pro- • subsection III D on the signatures of this rupture file [36] in the time dependence of the energy; 1 φ(x) = tanh ( x xc R0) • subsection III E 1 on the dependence of the − p2ξ j − j − antibubble-rupture time on the size of the antibub- 1 ble and the surface tension; + tanh ( x xc R1) 1:0; (4) p2ξ j − j − − • subsection III E 2 on the dependence of the where R0 and R1 are, respectively, the initial magnitudes antibubble-rupture time on the kinematic viscos- of the outer and inner radii of the antibubble, whose cen- ity; tre is initially at xc; the interface width ξ is measured by • subsection III F on the spatiotemporal evolution of the dimensionless Cahn number Ch = ξ=L. We use the antibubbles in our 3D DNSs. initial thickness of the antibubble shell h0 R0 R1 as a typical length scale. ≡ − We present conclusions and a discussion of our results in In our CHNS studies we use the Boussinesq approx- ρ + ρ Sec. IV. imation, wherein ρ(φ)D u 1 2 D u and [ρ(φ) t ≈ 2 t − ρ + ρ 1 2 ]g Aφg in Eq.(1); the surface-tension force II. MODEL AND NUMERICAL METHODS 2 ≈ − Fσ φ µ and the chemical potential µ(x; t) follow from≡ the − Cahn-Hilliardr free-energy functional (see the An antibubble consists of a shell of light fluid of density Supplemental Material [41]): F ρ1 inside a background heavy fluid of density ρ2; these fluids are immiscible and incompressible. The following µ = δ [φ]/δφ(x; t); CHNS equations provide a natural theoretical description FZ for such a binary-fluid mixture [36{39]: [φ] = Λ [(φ2 1)2=(4ξ2) + φ 2=2]dx; (5) F − jr j here, Λ is the energy density with which the two phases 2 ρ(φ)Dtu = p + &(φ) u + Fσ mix in the interfacial regime [36], ξ sets the scale of −∇ r p ρ1 + ρ2 2 2Λ + [ρ(φ) ]g αu; (1) the interface width , σ = 3ξ is the surface tension, − 2 − γ is the mobility [38] of the binary-fluid mixture; the D φ = γ 2µ and u = 0; (2) t r r · Cahn number Ch = ξ=L, where L is the linear size of our simulation domain, is a dimensionless measure of u In Eqs.(1)-(2), Dt (@t + ) is the convective deriva- the interfacial width; the Schmidt number Sc = ν=D. tive, u (u ; u ;≡ u ) is the· rfluid velocity, α is the co- ≡ x y z We assume that γ and ρ are independent of φ; we keep efficient of friction (we use this only in 2D), g is the γΛ the diffusivity D = 2 constant. We study the retrac- constant acceleration due to gravity, which points down- ξ tion dynamics of an antibubble by varying the Atwood wards, ρ(φ) ρ (1 + φ)=2 + ρ (1 φ)=2 is the density, 1 2 number A = (ρ ρ )=(ρ + ρ ) and the Bond num- &(φ) ν (1+≡φ)=2+ν (1 φ)=2 is the− dynamic viscosity, 2 1 2 1 1 2 ber Bo = Agρh2/σ−. The Bond number is a convenient φ(x; t≡) is the order-parameter− field at the point x and 0 dimensionless ratio of body forces (gravity), on the an- time t [with φ(x; t) > 0 in the background (majority) tibubble, and the surface or interfacial tension σ. We phase and φ(x; t) < 0 in the antibubble-shell (minority) p express times in multiples of τ R =Ag = 1:16.