arxiv: v1 [physics.flu-dyn] 21 Dec 2018

Transcription

1 1 axiv: v1 [physics.flu-dyn] 21 Dec 218 The cicula capillay jump Rajesh K. Bhagat 1, and P. F. Linden 2, 1 Depatment of Chemical Engineeing and Biotechnology, Univesity of Cambidge, Philippa Fawcett Dive, Cambidge, CB3 A 2 Depatment of Applied Mathematics and Theoetical Physics, Wilbefoce Road, Cambidge CB3 WA, UK Received xx; evised xx; accepted xx) In this pape we examine the flow poduced by the nomal impact of a lamina liquid jet on to an infinite plane. It is obseved expeimentally that afte impact the liquid speads adially ove the plane away fom the point of impact in a thin film. It is also obseved expeimentally that, at a finite adius, thee is an abupt incease in thickness of the film which has been identified as a hyaulic jump, and that this adius is independent of the oientation of the suface showing that gavity is unimpotant and that the flow is dominated by suface tension. We show that the application of consevation of momentum in the film flow subject to viscosity and suface tension alone pedicts a singulaity in the cuvatue of the liquid film and subsequently the depth of the film at a finite adius. This location is almost identical to the adius of the jump pedicted by consevation of enegy and agees with expeimental obsevations. 1. Intoduction In a ecent pape Bhagat et al. 218b) conducted expeiments which showed that in a thin liquid film, on scales typical in a kitchen sink, the jump poduced by the nomal impact of a lamina jet onto an infinite plane is independent of the oientation of the suface. Consequently, to fist ode, the flow and the location of the jump ae independent of gavity and detemined by suface tension alone. The adius of the jump, hee called a capillay jump to emphasise the dominance of suface tension ove gavity on these scales, was pedicted by applying the pinciple of consevation of enegy and balancing the etading foces of viscosity and suface tension. Hee we pedict the jump location using consevation of momentum applied to this adial flow, ignoing gavity completely. The pape is oganised as follows. The ole of suface tension at the liquid-gas inteface in a flowing, as distinct fom stationay liquid, is deived and applied to a adially speading thin film in 2. The application of consevation of momentum is both the adial and film-nomal diections is given in 3. Finally, some conclusions ae given in Theoy 2.1. uface tension foce Fluid motion is govened by the Navie-tokes equations which expess consevation of momentum in a fluid continuum. Fo flows with an inteface between two fluids, the Navie-tokes equations do not expess the suface tension foce acting on the inteface. This foce is intoduced as a nomal stess bounday condition at the inteface. Conside a suface with unit nomal n, bounded by a closed contou C with ac adess fo coespondence: kb29@cam.ac.uk

2 2 R.K. Bhagat, P.F. Linden Figue 1: chematic velocity pofiles in a flow with a) fee suface and b) a suface with non-zeo suface tension γ. The suface tension foce etads the flow nea the suface giving a non-zeo shea stess at the inteface 2.1). length l in the inteface between two immiscible fluids, taken hee to be the common case of a liquid and a gas denoted by the subscipts L and G, espectively, with constant suface tension γ. ince the suface tension foce acts in a diection pependicula n and the contou C, continuity of the nomal stess is expessed as T G T L ) nd + γ C dl n =, 2.1) whee T = pi + µ[ u + u) T ] is the total stess, with pessue p and velocity u, and µ is the viscosity of the fluid. Using a vecto identity and noting that is an abitay suface this gives the dynamic bounday condition T G T L ) n = γ s n)n, 2.2) whee s = [I nn]. is the suface gadient, elating the jump in the nomal stess to the cuvatue of the suface. Assuming the dynamic viscosity of the ai is negligible compaed to that of the liquid and denoting pessue in the ai as p G and in the liquid as p L, 2.1) can be witten as p G + p L )nd + µ n.[ u + u) T ]d + γ C dl n =, 2.3) whee u is the velocity in the liquid. In the case of a stationay liquid and gas, the middle tem of 2.3) is zeo and this equation gives the usual Laplace pessue in the liquid associated with the cuvatue of the suface. In the case of a flowing liquid the suface tension foce can be balanced by the viscous stesses at the suface, and fo a sufficiently fast flow this can be much lage than the Laplace pessue. This foce balance manifests itself though a non-zeo velocity shea at the suface as shown in figue Foce on an axisymmetic thin film Following Bush & Aistoff 23) fo an axisymmetic thin film on a plana suface we wite the equation of the suface in implicit fom J, z) = z h) =, 2.4)

3 The cicula capillay jump 3 Figue 2: chematic of the diffeential volume showing the slope of the thin liquid film which yields the nomal vecto n = J J = ẑ h ˆ, 2.5) 1 + h 2 ) 1/2 whee ˆ and ẑ ae unit vectos in the adial and wall-nomal diections, espectively, and h = dh/. We define the angle α as the tangent to the suface defined by h = tan α. 1 h Then cos α = and sin α =, and 2.5) can also be witten as 1+h 2 ) 1/2 1+h 2 ) 1/2 n = ẑ cos α ˆ sin α. 2.6) Conside the contol volume shown in figue 2. Hee dl = ˆθdθ and γ + dl n = γ + ˆθ ẑ cos α ˆ sin α)dθ, = γ 2π cos αˆ + 2π sin αẑ) +, = γ 2π cos αˆ + 2π sin αẑ), 2.7) + whee we have integated azimuthally fom to 2π. Noting that, since the adial length along the suface is cos α figue 2) so that d = 2π cos α, then fom 2.3) and 2.7) the adial and vetical components of nomal stess at the fee suface can be witten as + ) F γ, µ n [ u + u) T ] ˆ 2π ) + = P sin α2π cos α h cos α +2πγ cos α +, 2.8) and + ) F γ,z µ n [ u + u) T ] ẑ 2π ) + = P 2π+2πγ sin α +, cos α h 2.9) whee P = p L p G.

4 4 R.K. Bhagat, P.F. Linden As we will see below application of momentum consevation equies expessions fo the adial gadients of these foces and these ae given by df γ, and = d ) ) 2π µ n [ u + u) T ] ˆ cos α h = 2πP tan α + 2πγ d cos α) 2.1) df γ,z = d + ) ) ) 2π µ n [ u + u) T d sin α) ] ẑ = 2πP + 2πγ. cos α h 2.11) 3. Momentum consevation We now apply consevation of momentum to the axisymmetic flow speading adially fom the point of impact of the jet on the plane. We conside cylinical co-odinates and z in the the adial and wall-nomal jet-axial) coodinates, espectively, with the oigin at the point of impact of the jet on the plane, and u and w the associated velocity components, and we assume cicula symmety about the jet axis. In ode to detemine the velocity field we assume, following Watson 1964), that the flow is self simila and that the adial velocity can be expessed as the poduct of the adially vaying suface velocity u s ) and a function of the similaity vaiable η = z/h) in the fom u, z) = u s )fη), η z, η 1, 3.1) h) whee u s is the suface velocity and f) =, f1) = 1. Consevation of mass implies h udz = u s h 1 and C 1 1 fη)dη. Incompessibility gives fη)dη C 1 u s h = const., 3.2) Then, using 3.2) η wη) = = 1 1 u dz, d [u sh) η fη)dη]. 3.3) η w = u s h d fη)dη, = uh η = u s h ηfη). 3.4) which automatically satisfies the kinematic bounday conditions at the wall and inte-

5 The cicula capillay jump 5 face. We now use these expessions fo the velocity in equations expessing momentum consevation in the adial and wall-nomal diections Momentum balance in the adial diection ince the pessue in the gas at the suface is constant atmospheic) and the film is thin, the adial pessue gadient in the liquid is the Laplace pessue caused by the the cuvatue of the suface. In a thin film with h/r 1, whee R is the adius of the jump location, the cuvatue is small this tem can be ignoed. Then, in the absence of gavity, the momentum equation in the adial diection is ρ u u + w u z ) = 1 τ ) + τ z z, 3.5) whee τ = T + pi is the deviatoic stess tenso in cylinical coodinates). We use incompessibility u = ), integate 3.5) acoss the film fom the wall to the inteface and azimuthally fom to 2π, and substitute fo the velocity fom 3.1) and 3.4) and apply the suface tension bounday condition 2.1) and obtain d h ρu 2 d cos α) dz = P tan α + γ µ u z + w Noting that consevation of mass 3.2) implies u s h = const. we obtain C 2 [ρu s h du s d cos α) u ] = P tan α + γ µ z + w ) ). 3.6) 3.7) whee C 2 = 1 f 2 η)dη is a second integation constant aising fom the velocity pofile. Equation 3.7) can be witten in tems of the inteface slope C 2 [ρu s h du s ] = P 1 h + γ 1 + h 2 ) γ h h 1 + h 2 ) τ w, 3.8) 3/2 ) u whee τ w = µ z + w is the wall shea stess Momentum balance in wall-nomal diection We now apply consevation of momentum in the wall-nomal z diection in the diffeential contol volume shown in figue 2. ince the film is thin the pessue is independent of z and fo an axisymmetic flow ρu w w + ρw z = τ z) As befoe we integate acoss the film to obtain h ρu w dz + h ρw w z dz = h + τ zz z. 3.9) τ z ) dz z h. 3.1) ubstituting fo the velocity fom 3.1) and 3.4), and applying the suface bounday condition 2.1) gives ρu s h du sh ) +ρh 2 u 2 s ηf 2 η)dη ρh 2 u 2 s ηf 2 η)dη ρ u2 sh 2 η 2 f 2 η) 1 2 ηf 2 η)dη + ρw2 = df γ,z τz ) dz 2 z h. h 3.11)

6 6 R.K. Bhagat, P.F. Linden ubstituting 2.11) into 3.11) yields ρu s h du sh ) 1 ηf 2 η)dη = P d sin α) τz ) + γ dz cos α z h, = γp sin α + γ cos α dα τz ) dz z h, = P 1 + h 2 ) 1/2 γh γh h 2 ) 1/2 1 + h 3/2 ) τz ) dz z h. In the thin liquid film upsteam of the hyaulic jump, the inteface slope emains small and we will ignoe the highe ode tems in dh. Applying this appoximation and e-aanging 3.12) gives an expession fo the cuvatue of the film h ρu 2 sh 1 Finally, substituting ρu s h dus [C 3 ρu 2 sh γ]h + C 3 C 2 h 1 ηf 2 η)dη γ) + h ρu s h du s ηf 2 η)dη γ) = τz ) P dz z h. fom the adial momentum balance 3.8) gives 3.12) ) P h 1 + γ 1 + h 2 ) γ h h 1 + h 2 ) τ w C 2 γ 3/2 C ) τz ) = P dz z h, whee C 3 = 1 ηf 2 η)dη. Ignoing highe ode tems in h and e-aanging gives an expession fo the cuvatue of the film h [ C3 h C 2 τ w + 1 C3 C 2 )γ] P τ z) dz = C 3 ρu 2 sh1 We 1 ) whee the Webe numbe We is defined by z h 3.14) We C 3ρu 2 sh. 3.15) γ Consequently we pedict a singulaity in the cuvatue of the film at a citical adius whee the film thickness is such that We = Relation to enegy consevation Revisiting 2.8) and ecognising that upsteam of the hyaulic jump, the liquid film is almost flat which implies that α 2.8) can be witten as µ n.[ u + u) T ]d.ˆ γ dl n.ˆ = 2πγ 2πγ. 4.1) + h C Consequently, in a contol volume appoach the foce due to suface tension, which appeas as a nomal stess bounday condition, can be incopoated as a suface foce

7 The cicula capillay jump 7 on the cicumfeence of the contol volume consistent with the analysis in Bhagat et al. 218b) which applied consevation of enegy to this contol volume. The nomal stess bounday condition at the fee suface also implies that the inteio of the fluid expeiences the effect of suface tension though a change in the velocity pofile. In ou pevious analysis Bhagat et al. 218b)) we did not conside this and simply used Watson s similaity velocity pofile Watson 1964) which assumed zeo stess u z = ) at the fee suface. The suface velocity obtained using Watson s similaity velocity pofile is, theefoe, not the coect suface velocity see figue 1). Futhemoe, tevens & Webb 1993) measued the velocity pofile in the liquid film and obseved the expected lowe suface velocity. Consequently, it is appopiate to associate the flux of suface enegy with the aveage velocity athe than the suface velocity obtained fom a bounday laye velocity pofile see Bhagat et al. 218a)). 5. Conclusions Applying consevation of adial and wall-nomal momentum to the flow in an expanding axisymmmetic thin film shows that the cuvatue of the film is singula at a finite adius detemined by a citical value of the Webe numbe. This singulaity aises fom the wall-nomal momentum consevation which implies that d2 h wheeas dh 2 )2 at a finite adius. Using consevation of enegy Bhagat et al. 218b) showed that the adial velocity gadient is also singula, in this case dh, at a citical Webe numbe which is numeically slightly diffeent. This adius was identified in expeiments as a jump in the flow depth to a thicke and slowe flow downsteam, and excellent quantitative ageement was found in the pedicted and obseved values of the jump adius. Tanslating these esults into pedictions of the jump adius, consevation of enegy pedicts that the jump adius R is given by ) 1/4 R 1 C 2 = R f )2π) 3 C1 3 =.275, 5.1) whee R = Q3/4 ρ 1/4 ν 1/4 γ 1/4. On the othe hand consevation of momentum pedicts R R = 1 f )2π) 3 C 3 C 3 1 ) 1/4 = ) The numeical values ae obtained fom Watson s similaity pofile which, as we discuss above, do not stictly apply, and the pedictions ae both smalle than the expeimentally measued values of.289 ±.15. ince C 1 is the aea unde the cuve fη) then this will be smalle fo the eal pofile figue 1), leading to a lage pediction. Also, since η 1, C 2 > C 3, the jump adius estimate fom momentum consevation is always smalle than that obtained fom enegy consevation. Acknowledgement RKB would like to thank t John s College and the Commonwealth cholaship Commission. RKB would also like to thank his PhD supeviso Pof Ian Wilson fo his advice and suppot.

8 8 R.K. Bhagat, P.F. Linden REFERENCE Bhagat, R. K., Jha, N. K., Linden, P. F. & Wilson, D I. 218a On the oigin of the cicula hyaulic jump: a diffeential analysis. papescontibution 819 final.pdf IBN Bhagat, R. K., Jha, N. K., Linden, P. F. & Wilson, D I. 218b On the oigin of the cicula hyaulic jump in a thin liquid film. Jounal of Fluid Mechanics 851, R5. Bush, J. W. M. & Aistoff, J. M. 23 The influence of suface tension on the cicula hyaulic jump. Jounal of Fluid Mechanics 489, tevens, J. & Webb, B. W Measuements of flow stuctue in the adial laye of impinging fee-suface liquid jets. Intenational Jounal of Heat and Mass Tansfe 36 15), Watson, E. J The adial spead of a liquid jet ove a hoizontal plane. Jounal of Fluid Mechanics 2 3),