∫° Nonlinear Capillary Oscillation of a Charged Drop
Total Page:16
File Type:pdf, Size:1020Kb
Technical Physics, Vol. 45, No. 8, 2000, pp. 1001–1008. Translated from Zhurnal Tekhnicheskoœ Fiziki, Vol. 70, No. 8, 2000, pp. 45–52. Original Russian Text Copyright © 2000 by Belonozhko, Grigor’ev. GASES AND LIQUIDS Nonlinear Capillary Oscillation of a Charged Drop D. F. Belonozhko and A. I. Grigor’ev Yaroslavl State University, Sovetskaya ul. 14, Yaroslavl, 150000 Russia E-mail: [email protected] Received May 18, 1999; in final form October 14, 1999 Abstract—In the quadratic approximation with respect to the amplitudes of capillary oscillation and velocity field of the liquid moving inside a charged drop of a perfectly conducting fluid, it is shown that the liquid drop oscillates about a weakly prolate form. This refines the result obtained in the linear theory developed by Lord Rayleigh, who predicted oscillation about a spherical form. The extent of elongation is proportional to the initial amplitude of the principal mode and increases with the intrinsic charge carried by the drop. An estimate is obtained for the characteristic time of instability development for a critically charged drop. © 2000 MAIK “Nauka/Interperiodica”. Studies of capillary oscillation and stability of a 2. Consider a drop of an inviscid, perfectly conduct- charged drop are motivated by numerous applications ing liquid of density ρ characterized by the surface ten- [1]. The problem was reviewed in [1–7] and references sion γ and carrying a charge Q. The linear theory of therein. However, most theoretical analyses are based capillary oscillation postulates that there exists a time on the linearized system of fluid-dynamics equations. moment that can be treated as t = 0, when the drop Only recent studies have captured the nonlinear nature geometry is described by the second mode of linear of the phenomenon and provided essentially new infor- capillary oscillation of a small finite amplitude ε and mation about the mechanism of instability of a highly the velocity field is zero. The initial drop volume is charged drop and its capillary oscillation [2–7]. equal to that of a spherical drop of radius R. The prob- lem is to determine the axially symmetric oscillations 1. A charged spherical drop with an intrinsic charge of the drop in the case of an axially symmetric potential Q greater than a certain critical value becomes unsta- velocity field in the drop. ble, because the electrical repulsive forces exceed the We use the dimensionless variables in which ρ ≡ 1, surface tension forces. At the close of the nineteenth γ ≡ ≡ century, Lord Rayleigh developed a linear model of this 1, and R 1. The mathematical model of the prob- instability [8]. Since then, studies of the instability of a lem is then written in spherical coordinates with the ori- highly charged drop and its generalizations have grown gin at the center of the drop as follows: into a broad area of research having important technical ∆ψ ==0; U — ⋅ ψ, (1) applications. ∆Φ ==0; E–—⋅Φ, (2) From the perspective of classical physics, the ther- mal motion of molecules gives rise to drop oscillations r ∞: |—Φ| 0, (3) of amplitude ~kT/γ about the equilibrium spherical r = 0: |—ψ| < ∞, (4) form, where k is Boltzmann’s constant, T is the absolute temperature, and γ is the surface tension of the liquid. r=1+ξΘ(),t Each mode is characterized by a specific critical value 1∂Φ –------ ∫------- dS =Q; Sr(),,Θϕ ≡0≤≤Θπ (5) of surface charge, and a supercritical charge leads to the 4π°∂n onset of instability. The stability of the drop with 0≤ϕ<2π, respect to its intrinsic charge Q is generally character- ≡ 2 πγ 3 ized by the so-called Rayleigh parameter W Q /4 R , r==1+ξΘ(),t: Φ const, (6) where R is the drop radius. The second mode is the most unstable, and the corresponding critical value of ∂ψ ∂ξ 1 ∂ξ ∂ψ ------- = ------ + ---- ------- ------- , (7) the Rayleigh parameter is W = 4. 2∂Θ∂Θ c ∂r∂tr In this paper, we solve Rayleigh’s problem in the 2 quadratic approximation with respect to the amplitudes ∂ψ 1 2 E ∆F– ------- –--- ()— ⋅ ψ +F==Fγ;F------ , (8) of the velocity field and surface perturbation and deter- ∂t 2 E E 8π mine the geometric form about which a subcritically ξ ε ()Θ charged liquid drop oscillates. t==0: r 1 ++* P2 cos , (9) 1063-7842/00/4508-1001 $20.00 © 2000 MAIK “Nauka/Interperiodica” 1002 BELONOZHKO, GRIGOR’EV In the zeroth approximation, the nontrivial relations 1 2 4 P (µ) = --(3µ – 1), dv = --π; are (2), (3), (5), and (8). They contain the electric 2 2 ∫ 3 potential, the capillary pressure, and the constant com- V (10) ponent of the difference of pressures inside and outside ≤ ≤ ε ε ( (Θ)) the drop. Problem (1)–(11) reduces to one for the static 0 r 1 + + P2 cos , Φ ψ ≡ ξ ≡ V(r, Θ, t) ≡ * field 0(r), with the known functions 0 0 and 0 0: 0 ≤ Θ ≤ π, 0 ≤ ϕ < 2π, 2Φ Φ d 0 d 0 ψ = 0. (11) ----------- = 0, E1 = –---------er, dr2 dr This is a boundary value problem with the free sur- face r = 1 + ξ(Θ, t) and three unknown functions: the 2π π Φ potential Φ = Φ(r, Θ, t) of the electric field outside the 1 d 0 (Θ) Θ ϕ ψ ψ Θ –----π-- ∫ ∫--------- sin d d = Q, (16) drop, the potential = (r, , t) of the velocity field 4 dr r = 1 inside the drop, and the surface deviation ξ = ξ(Θ, t) 0 0 from a regular sphere with r = 1. Φ ∞ d 0 The potentials Φ and ψ satisfy the Laplace equa- r : --------- 0, tions (1) and (2) under the boundary and initial condi- dr tions (3)–(8) and (9)–(11), respectively. where er is the unit radial vector. Both kinematic and dynamic boundary conditions Φ The solution of the problem is 0 = Q/r. Using this for the velocity field U, (7) and (8), contain nonlinear solution and formulas (1B), (2C), and (5D)–(10D) from terms. Condition (4) implies that U is bounded at the ∆F the Appendices, the first- and second order problems center of the drop. In (8), is the difference between are easily formulated. The first order problem is written the constant pressure components inside and outside as the drop, F is the electric pressure on the drop surface E ∆ψ ⋅ ψ [9, 10], and Fγ is the capillary pressure under the 1 = 0, U1 = — 1, (17) curved drop surface [11]. ∆Φ = 0, E = –— ⋅ Φ , (18) It follows from (3) that the electric field E must be 1 1 1 ∞ | Φ | bounded at infinity. The continuity of the tangential r : — 1 0, (19) components of E(r, t) at the charged surface of a con- | ψ | ∞ ductor means that the surface is equipotential (see (6)) r = 0: — 1 < , (20) [10]. The condition for the normal components of E(r, t) 2π π at the interface is written in (5) in an integral form, Φ 1 d 1 (Θ) Θ ϕ where n is the outward unit normal to the drop surface. r = 1: –----π-- ∫ ∫--------- sin d d = 0, (21) ξ 4 dr The quantity ∗ in (9) is determined by (10), which 0 0 means that the volumes of the initial drop and a spher- Φ – ξ Q = 0, (22) ical drop of unit radius are equal. 1 1 3. In dimensionless variables, the amplitude ε of an ∂ψ ∂ξ ---------1 = --------1, (23) initial deviation from spherical form is measured in ∂r ∂t units of the spherical drop radius. We treat ε as a small parameter and seek a solution to the problem as a series ∂ψ ∂Φ ε 1 Q 1 ξ ξ ∆ ξ expansion in integer powers of , omitting the terms of ---∂------ + ----π-- ---∂------ + 2 1Q = 2 1 + Ω 1, (24) order higher than two: t 4 r 3 ξ ε ( (Θ)) Φ Φ Φ Φ (ε ) t = 0: 1 = P2 cos ; = 0 + 1 + 2 + O , (12) 1 2 (25) ψ ψ ψ (ε3) P (cos(Θ)) = --(3cos (Θ) – 1), = 1 + 2 + O , (13) 2 2 ξ ξ ξ ξ (ε3) ψ r = 1 + , = 1 + 2 + O , (14) 1 = 0. (26) Φ ∼ ψ ∼ ξ ∼ O(εm); m = 1, 2; Φ = O(1). (15) The problem for the second-order quantities is writ- m m m 0 ten as follows: The functions Φ, ψ, and ξ are assumed to be of the ∆ψ ⋅ ψ same order as their partial derivatives. Under these 2 = 0, U2 = — 2, (27) assumptions, since the Laplace operator is linear, we ∆Φ ⋅ Φ can use the expansions obtained in the Appendices to 2 = 0, E2 = –— 2, (28) split the problem (1)–(11) into problems of zeroth, first, ∞ | Φ | Φ Φ r : — · 2 0, (29) and second order [11] for the seven functions 0, m, ψ ξ | ψ | ∞ m, and m, where m = 1 and 2. r = 0: — · 2 < , (30) TECHNICAL PHYSICS Vol. 45 No. 8 2000 NONLINEAR CAPILLARY OSCILLATION OF A CHARGED DROP 1003 π π − Θ 2 2 m(m + 1)Pm(cos )); and formulate the Cauchy prob- 1 dΦ ∂Φ ∂Φ r = 1: –------ ---------2 + 2ξ ---------1 + ξ ---------1 lem for a homogeneous system of ordinary differential π ∫ ∫ 1 ∂ 1 2 4 dr r ∂r equations for the unknown functions Cm(t), Fm(t), and 0 0 (31) Z (t). By finding these unknowns and substituting them ∂ξ ∂Φ m 1 1 Θ Θ ϕ into (37)–(39), we can write the stable solutions to the – -∂---Θ------∂---Θ---- sin d d = 0, first-order problem as ∂Φ ξ = εcos(ω t)P (cosΘ), (40) Φ ξ ξ 1 ξ2 1 2 2 2 – 2Q = – 1--------- + 1 Q, (32) ∂r εω ψ = –--------2 sin(ω t)P (cosΘ), (41) ∂ψ ∂2ψ ∂ξ ∂ξ ∂ψ 1 2 2 2 ---------2 + ξ -----------1 = --------2 + --------1 ---------1, (33) ∂r 1 2 ∂t ∂Θ ∂Θ Φ ε (ω ) ( Θ) ∂r 1 = Qcos 2t P2 cos , (42) ∂ψ ∂Φ 2 2 Q 2 ξ ξ ∆ ξ Q ---∂------ + ----π-- ---∂------ + 2 2Q = 2 2 + Ω 2 ω = 8 – 2W, W = ------ < 4.