24.12.2013 Views

The force on a bubble, drop, or particle in arbitrary time-dependent ...

The force on a bubble, drop, or particle in arbitrary time-dependent ...

The force on a bubble, drop, or particle in arbitrary time-dependent ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

T*)dS<br />

J'+ VfJV,')-lV fVd Vt* 9*I<br />

-PVduWJ (t = -RFU*- QU~), (42)<br />

where Us(t) [=U(t)-U'(t)] is the slip velocity of the<br />

<strong>drop</strong>. Here we have used the fact that the first two <strong>in</strong>tegrals<br />

<strong>on</strong> the rhs of (26) may be simplified by not<strong>in</strong>g<br />

fsd(nT)dS = -RFU and Sfd(n * = 0. <str<strong>on</strong>g>The</str<strong>on</strong>g> goal<br />

now is to estimate the c<strong>on</strong>tributi<strong>on</strong>s from the two <strong>in</strong>tegrals<br />

<strong>in</strong> (42) with the c<strong>on</strong>diti<strong>on</strong> that the Reynolds number<br />

(Re=aU/v) f<strong>or</strong> the fluid <strong>in</strong>side and outside the <strong>drop</strong>,<br />

based <strong>on</strong> the <strong>drop</strong>'s slip velocity, is small but f<strong>in</strong>ite. (F<strong>or</strong> a<br />

n<strong>on</strong>spherical body, a denotes the characteristic <strong>drop</strong> dimensi<strong>on</strong>;<br />

otherwise it is the <strong>drop</strong> radius.) In so do<strong>in</strong>g, we<br />

will obta<strong>in</strong> an expressi<strong>on</strong> f<strong>or</strong> the hydrodynamic <str<strong>on</strong>g>f<strong>or</strong>ce</str<strong>on</strong>g> act<strong>in</strong>g<br />

<strong>on</strong> the <strong>drop</strong> to O(Re) f<strong>or</strong> <strong>arbitrary</strong> <strong>time</strong>-<strong>dependent</strong><br />

moti<strong>on</strong>.<br />

First note the follow<strong>in</strong>g equalities f<strong>or</strong> the fluid exteri<strong>or</strong><br />

to the <strong>drop</strong>:<br />

-A-U5-Vu'+u' -VU') =-Vp'-V<br />

And if we def<strong>in</strong>e<br />

Pa<br />

2 u'=V a'.<br />

(43)<br />

we can note the follow<strong>in</strong>g equalities f<strong>or</strong> the fluid <strong>in</strong>side the<br />

<strong>drop</strong>:<br />

* **(a1u*' Vu*±* .*)<br />

=-vp*t +,*v2u*' =V O,*-p*U-, H-<br />

(45)<br />

where we have applied the c<strong>on</strong>diti<strong>on</strong> that<br />

-(p*/p)VpO =p*Uw (t). Us<strong>in</strong>g f' and f*' to signify the<br />

<strong>in</strong>ertial terms from the first equalities of (43) and (45), the<br />

hydrodynamic <str<strong>on</strong>g>f<strong>or</strong>ce</str<strong>on</strong>g> can now be expressed as<br />

FRj-pVdU''(t)=-RFU (t)- ff *MdV<br />

Vf<br />

fVd f*' (M* -)dV. (46)<br />

Here we have used arguments similar to (30) to show<br />

Now s<strong>in</strong>ce the boundary c<strong>on</strong>diti<strong>on</strong>s f<strong>or</strong> the "primed"<br />

fields are<br />

and<br />

n. [[L(Vu'+Vu' Ț ) _*(Vu*t+Vu*, T )] * (I-nn) =0,<br />

u =u*', n-'nu*= J <strong>on</strong> Sd<br />

(48a)<br />

u'-0, p'-*O as r-,co (48b)<br />

it can be seen that the two volume <strong>in</strong>tegrals <strong>in</strong> (46) represent<br />

the <strong>in</strong>ertial c<strong>or</strong>recti<strong>on</strong>s to the steady Stokes drag f<strong>or</strong><br />

a <strong>drop</strong> translat<strong>in</strong>g with velocity Us <strong>in</strong> a quiescent fluid. In<br />

additi<strong>on</strong>, other than the presence of the <strong>in</strong>tegral over the<br />

volume of the <strong>drop</strong>, (46) is identical to the expressi<strong>on</strong> f<strong>or</strong><br />

a solid <strong>particle</strong>. Thus, with appropriate modificati<strong>on</strong>s, we<br />

can make use of the results f<strong>or</strong> solid <strong>particle</strong>s from LB. We<br />

will summarize the general ideas from that paper to show<br />

the similarities and differences with the current derivati<strong>on</strong>.<br />

<str<strong>on</strong>g>The</str<strong>on</strong>g> <strong>in</strong>terested reader is referred to the <strong>or</strong>ig<strong>in</strong>al w<strong>or</strong>k f<strong>or</strong><br />

further details.<br />

F<strong>or</strong> small Reynolds number and sh<strong>or</strong>t <strong>time</strong> scale moti<strong>on</strong><br />

(i-s < v/U,2 where r, is the <strong>time</strong> scale f<strong>or</strong> the change <strong>in</strong><br />

the <strong>drop</strong>s slip velocity), the flow is governed, to lead<strong>in</strong>g<br />

<strong>or</strong>der <strong>in</strong> Re, by the unsteady Stokes equati<strong>on</strong>s throughout<br />

the fluid doma<strong>in</strong>:<br />

au'<br />

This approximati<strong>on</strong> is appropriate f<strong>or</strong> the flow <strong>in</strong>side as<br />

well as outside the <strong>drop</strong>. <str<strong>on</strong>g>The</str<strong>on</strong>g> c<strong>on</strong>vective terms of the<br />

Navier-Stokes equati<strong>on</strong>s, u* Vu and Us- Vu, are everywhere<br />

smaller than the viscous <strong>or</strong> the unsteady <strong>in</strong>ertial<br />

terms, because the v<strong>or</strong>ticity produced at the surface of the<br />

<strong>drop</strong> has not diffused out to the Oseen distance, vIUc,<br />

where c<strong>on</strong>vecti<strong>on</strong> becomes imp<strong>or</strong>tant as a transp<strong>or</strong>t mechanism.<br />

Under these c<strong>on</strong>diti<strong>on</strong>s, the c<strong>on</strong>tributi<strong>on</strong>s from the<br />

c<strong>on</strong>vective terms are obta<strong>in</strong>ed solely from a regular perturbati<strong>on</strong><br />

analysis.<br />

On the other hand, f<strong>or</strong> l<strong>on</strong>g <strong>time</strong> scale moti<strong>on</strong><br />

(r->v/U2) the flow <strong>in</strong> the near-field regi<strong>on</strong> (f<strong>or</strong> length<br />

scales sh<strong>or</strong>ter than the Oseen distance v/vc) is governed by<br />

the steady Stokes equati<strong>on</strong>s, while that <strong>in</strong> the far field [def<strong>in</strong>ed<br />

by distances from the <strong>drop</strong> of O(v/Uc) <strong>or</strong> greater] is<br />

determ<strong>in</strong>ed by the unsteady Oseen equati<strong>on</strong>s to lead<strong>in</strong>g<br />

<strong>or</strong>der:<br />

-Vp' +tV~uI~(-j 7 -U. CVu') +F 't8(X),<br />

(49)<br />

.fVd(V-a-*)<br />

(Ml*-~I)dV- f~df*t. (M1*1)dV<br />

= P*toQ(t) .J~ (MA*-I)dV<br />

= P*to (t). TfV V [(M1*-I)r~dV<br />

~P*iOO (t). .fSd n-[(1M*-I)r~dS=0.<br />

2109 Phys. Fluids A, Vol. 5, No. 9, September 1993<br />

V -u'= O. (50)<br />

Here, the boundary c<strong>on</strong>diti<strong>on</strong>s at the <strong>drop</strong> surface are replaced<br />

by the presence of the <str<strong>on</strong>g>f<strong>or</strong>ce</str<strong>on</strong>g> m<strong>on</strong>opole <strong>in</strong> the govern<strong>in</strong>g<br />

equati<strong>on</strong>; that is, to lead<strong>in</strong>g <strong>or</strong>der <strong>in</strong> the far-field<br />

regi<strong>on</strong>, the <strong>particle</strong> appears as a po<strong>in</strong>t-<str<strong>on</strong>g>f<strong>or</strong>ce</str<strong>on</strong>g> disturbance of<br />

magnitude the pseudosteady Stokes drag FIS<br />

[= -R k, UJ(t)]. Also, <strong>in</strong> this case, diffusi<strong>on</strong> and c<strong>on</strong>vecti<strong>on</strong><br />

are of equal imp<strong>or</strong>tance <strong>in</strong> the transp<strong>or</strong>t of v<strong>or</strong>ticity.<br />

F<strong>or</strong> moti<strong>on</strong> of <strong>arbitrary</strong> <strong>time</strong> scale, the unsteady Stokes<br />

(47) equati<strong>on</strong>s describe the flow to lead<strong>in</strong>g <strong>or</strong>der everywhere,<br />

except <strong>in</strong> the far field where the unsteady Oseen equati<strong>on</strong>s<br />

P. M. Lovalenti and J. F. Brady 2109<br />

Downloaded 13 Jan 2006 to 131.215.225.172. Redistributi<strong>on</strong> subject to AIP license <strong>or</strong> copyright, see http://pof.aip.<strong>or</strong>g/pof/copyright.jsp

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!