See discussions, stats, and author profiles for this publication at: [Link]

net/publication/305721466

Deducción de la fuerza de Stokes a partir de las ecuaciones de Navier-Stokes

Article · July 2016

CITATIONS READS

0 1,891

1 author:

Jonathan Estévez-Fernández
Complutense University of Madrid
117 PUBLICATIONS 0 CITATIONS

SEE PROFILE

All content following this page was uploaded by Jonathan Estévez-Fernández on 29 July 2016.

The user has requested enhancement of the downloaded file.

 PHYSICAL REVIEW D 73, 174412 (2016)

Deducción de la fuerza de Stokes a partir de las ecuaciones de Navier-Stokes

¨
Jonathan Estévez-Fernández
Facultad de Ciencias Fı́sicas, Universidad Complutense, 28040 Madrid, España
(Elaborado el 29 de julio de 2016)

En este artı́culo me propongo deducir analı́ticamente a partir de las ecuaciones de

Navier-Stokes la expresión para la fuerza de Stokes de una gota de agua de lluvia de
pequeño tamaño que cae lentamente (en aproximación de fluido de Stokes cuyo número
de Reynolds es mucho menor que la unidad) por el método de Hadamard-Rybczynski.
¨
jonestev@[Link]

La famosa fórmula de la fuerza de Stokes fue obte- Introducimos un sistema de coordenadas esféricas .r; ✓; '/
nida por primera vez en 1851 por Sir George Gabriel centrado en el centro de la gota de radio r0 . Las com-
Stokes (1819-1903) referida al movimiento de roza- ponentes de la velocidad del flujo de aire viscoso e
miento de una esfera (no necesariamente sólida) por incompresible alrededor de la gota v D .vr ; v✓ ; v' /
un fluido (movimiento de arrastre o drag). La fórmula tienen únicamente dependencia radial y polar, la com-
para el caso de una gota de agua de lluvia que cae por ponente azimutal es nula v' D 0 por la asimetrı́a del
la atmósfera es problema.
Las ecuaciones de Navier-Stokes para v W R4 ! R3
F D 6†r0 v0 e3 (i.1)
son
donde r0 > 0 es el radio de la esfera, > 0 el coefi- Ä @v
ciente de viscosidad (medido en poises o pascales por C .v r/v vD rp C f; (1.1a)
@t
segundo), v0 > 0 es la velocidad de caı́da (supuesta
lı́mite) y ei es el i-ésimo vector unitario coordenado rv D 0; (1.1b)
en un sistema de coordenadas cartesianas (coeficien-
tes métricos constantes en R3 ) ortogonales fe1 ; e2 ; e3 g. con la condición inicial v.r; 0/ D v0 , donde es el
Nótese que al ser F3 > 0 y F1 ; F2 D 0, la fuerza de coeficiente de viscosidad (que aparece en (i.1)), p es la
Stokes tiene la dirección del eje vertical e3 y sentido presión efectiva del fluido sobre la gota (donde se ha
positivo (hacia arriba) en todo momento t > 0. considerado la presión atmosférica p0 ) y f es la fuerza
externa. En nuestro caso, las ecuaciones de Navier-
Stokes (1.1) en aproximación de Stokes se comportan
para v W R3 ! R3 como
1. Ecuaciones de Navier-Stokes. (
v D rp; (1.2a)
Bajo la aproximación de Stokes (creeping flow re- rv D 0; (1.2b)
gime) el número de Reynolds ha de ser
con la condición v.r0 / D .vr ; v✓ ; 0/: Las propiedades
⇢r0 v0 de este fluido (la atmósfera) que verifica (1.2) son
NRe ⌘ ⌧ 1;
✏ Estacionario:
que se cumple para esferas de pequeño tamaño sumer- @v
D 0:
gidas en un gas o un lı́quido. La fórmula es usada para @t
hallar la velocidad lı́mite de caı́da de las gotas de agua ✏ Incompresible:
de lluvia o para determinar la carga del electrón en la
caı́da de gotas ionizadas de aceite en un condensador. rv D 0 , ⇢ D cte. 8r 2 U ⇢ R3

1
EST ÉVEZ PHYSICAL REVIEW D 73, 174412 (2016)

que es, de hecho, la ecuación (1.2b). Por tanto, la ecua- 1.1. Ecuaciones de Cauchy del momento
ción (1.2b) que en coordenadas esféricas se escribe
Sustituyendo en (1.4) obtenemos las ecuaciones
1 @ 1 @ 2
rv D .v✓ sin ✓ / C 2 .r vr / de Cauchy del momento (las presiones sobre la gota;
r sin ✓ @✓ ✓ ◆ r @r ✓ ◆ p' D 0)
1 @ 1@ 1 @ 1 @
D
r sin ✓ @✓ r @r
C 2
r @r sin ✓ @✓
D0
@p
‚@.Q /
D 2 ; (1.5a)
@r r sin ✓ @✓
donde se han introducido la función de corriente de
Stokes (stream function) .r; ✓ / de clase C 2 (para 1 @p @.Q /
D ; (1.5b)
verificar el lema de Schwarz) definido por r @✓ r sin ✓ @r

1 @ siendo p del clase C 2 que verifica el lema de Schwarz.

1 @
vr D 2 ; v✓ D Por este hecho, tomando derivadas parciales cruzadas
:
r sin ✓ @✓ r sin ✓ @r
respectivamente a ambos miembros y restando las ecua-
Analicemos ahora la ecuación de Cauchy del momento ciones obtenemos una EDP de cuarto orden
(1.2a) sabiendo que (identidad vectorial)
.@2r✓ @2✓ r /p D Q2 D 0: (1.6)
r ⇥ .r ⇥ v/ D [Link]/ v:
Para resolver la ec. (1.6) utilicemos el método de se-
Por (1.2b) o, equivalentemente (1.1b), y definiendo la paración de variables para .r; ✓ / D R.r/⇥.✓ /. Para
vorticidad como ! ⌘ r ⇥ v obtenemos analizar el comportamiento de ⇥.✓ / evaluemos una
condición de contorno
v D r ⇥ !: (1.3)
Z ✓ Z ✓
@
Sustituyendo (1.3) en (1.2a) obtenemos jr!1 D ds D vr r 2 sin s ds
0 @s 0
Z ✓
rp D r⇥ : ! (1.4) 1
D r2 v1 cos s sin s ds D r 2 v1 sin2 ✓ C g.r/;
0 2
Calculemos ahora la vorticidad: como el campo de
velocidades es este es el llamado método de Hadamard-Rybczynski.
Probemos con
v D vr .r; ✓ /er C v✓ .r; ✓ /e✓
D f .r/ sin2 ✓
su rotacional debe ser perpendicular al plano al que
pertenece v ya que v' D 0: que sustituyéndolo en (1.6)
ˇ ˇ h ⇤
ˇ er re✓ r sin ✓ e' ˇ 2 2
sin2 ✓ f 00 .r/
2
ˇ ˇ 0 D Q .f .r/ sin ✓/ D Q f .r/
! D r ⇥ v D ˇˇ@r @✓ @' ˇˇ D !' e' ; r2
ˇvr rv✓ r sin ✓ v' ˇ ⌘ Q.g.r/ sin ✓/ D h.r/ sin2 ✓
2

con donde
✓ ◆
1 @.rv✓ / @vr 1 2
!' D e' D Q h.r/ D g 00 .r/ D @2rr f 00 .r/ f .r/
r @r @✓ r sin ✓ r2
4 00 8 8
siendo Q un operador diferencial definido por D .f .4/ 2
f C 3f 0 f /:
r r r4
✓ ◆
2 sin ✓ 1 Nótese que la fórmula generalizada de Leibniz es
Q ⌘ @r C 2 @✓ @✓ :
r sin ✓
n ✓ ◆
X n
Hallemos el rotacional de la vorticidad .uv/.n/ D u.i / v .n i/
:
i
iD0
sin ✓!' r!'
✓ π
◆ ✓ »
◆ Es fácil ver que (1.6) es simplemente
1 1 1 1
r ⇥! D @✓ Q er @r Q e✓
r sin ✓ r r sin ✓ r 4 h.r/ D 0
1 1
D @✓ .Q /er C @r .Q /e✓ :
r sin ✓
2 r sin ✓

2
 TRAYECTORIA CICLOIDAL... PHYSICAL REVIEW D 73, 174412 (2016)

que es una ecuación diferencial ordinaria (EDO) de Las condiciones de contorno determinan las constantes:
Cauchy-Euler de cuarto orden en r la condición de vr y v✓ en la superficie determina c3 D
v0 =2 y c4 D 0, y la condición en el infinito determina
r 4 f .4/ .r/ 4r 2 f 00 .r/ C 8rf 0 .r/ 8f .r/ D 0: c1 D v0 r03 =4 y c2 D 3r0 v0 =4:
Realizando el cambio de variable independiente t D vr ; v✓ ! 0;
log jrj para r ¤ 0 y con u.t / D f .r/ se obtiene rDr0

y
u.4/ 6u000 C 7u00 C 6u0 8u D 0: (1.7)
vr ! v0 cos ✓; vr ! v0 sin ✓:
r!1 r!1
Nótese que la regla de la cadena para el cambio en las
Finalmente, la función de Stokes de corriente es
derivadas es para la primera
!
du df dr dr r03 3r0 r r2
u0 D D D f 0 .r/ D rf 0 .r/; .r; ✓ / D v0 C sin2 ✓;
dt dr dt dt 4r 4 2
para la segunda que sutituyéndolo en las ecuaciones de la velocidad
✓ ◆ resulta
du0 du0 =dr dr du0
u00 D D D !
dt dt=dr dr dr r03 3r0
D r f .r/ C rf .r/ D rf 0 .r/ C r 2 f 00 .r/;
0 00 vr D v0 C 1 cos ✓;
2r 3 2r
para la tercera !
✓ ◆ r03 3r0
du00 du00 =dr dr du00 v✓ D v0 3
C 1 sin ✓;
000
u D D D 4r 4r
dt dt=dr dr dr
y ✓ ◆
D r f 0 .r/ C rf 00 .r/ C 2rf 00 .r/ C r 2 f .3/ .r/ 3r0 v0
pD cos ✓:
D rf 0 .r/ C 3r 2 f 00 .r/ C r 3 f .3/ .r/; 2r 2
La fuerza debida al campo escalar de presiones es
y para la cuarta Z
!
du.3/ du.3/ =dr dr du.3/ Fpresión D p´ dSe3
u.4/ D D D S0
dt dt=dr dr dr
donde el diferencial de superficie de la gota esférica en
D r f 0 .r/ C rf 00 .r/ C 6rf 00 .r/ C 6r 2 f .3/ .r/ r D r0 es
C r 3 f .4/ .r/ dS D r02 sin ✓ d✓d'
D rf 0 .r/ C 7r 2 f 00 .r/ C 6r 3 f .3/ .r/ C r 4 f .4/ .r/: que integrando en la proyección vertical de la presión
normal (donde se sitúa la fuerza)
El proceso de cambio de variable para obtener (1.7)
Z 2† Z †
es proceder primero por las derivadas de orden mayor 2
Fpresión D r0 d' .p cos ✓/ sin ✓ d✓ e3
sabiendo que 0
Z 0†
r 4 f .4/ .r/ D u.4/ F .r/: D 3†r0 v0 sin ✓ cos2 ✓ d✓ e3 D 2†r0 v0 e3 :

La ec. (1.7) es una EDO lineal con coeficientes cons-

ù 0

2=3
tantes que se resuelve fácilmente gracias al polinomio
caracterı́stico Por otro lugar, ahora hay que integrar la presión tangen-
te a la superficie que es la tensión superficial definida
pA .r/ D 4 6 3 C7 2 C6 8 D . 2 1/. 2/. 4/ D 0 por un tensor tal que la componente que contribuye a
la fuerza de Stokes es
cuyas raı́ces son ¢ D f 1; 1; 2; 4g. La solución son
exponenciales en t que mediante el cambio de variable  ⇣v ⌘ 1
✓ 3r03 v0
anterior se transforma en ⌧r✓ D r@r C .@ ✓ vr / D sin ✓:
r r 2r 4
c1
f .r/ D C c2 r C c3 r 2 C c4 r 4 :
r

3
 TRAYECTORIA CICLOIDAL... PHYSICAL REVIEW D 73, 174412 (2016)

Luego de manera análoga Además, si igualamos la fuerza gravitatoria con el em-

Z puje de Arquı́medes (buoyancy force) y la fuerza de
Ftensión superficial D .⌧r✓ sin ✓ / dS e3 Stokes en el equilibrio
S0
Z †
D 3†r0 v0 sin3 ✓ d✓ e3 D 4†r0 v0 e3 : mgota g D 6†r0 v0 C ⇢aire gVgota
ö 0

4=3 resulta que la velocidad lı́mite de la gota de agua es

Sumando ambas fuerzas resulta la fórmula de la fuerza
de Stokes 2gr02
v1 ⌘ v0 D .⇢gota ⇢aire /:
Fpresión C Ftensión superficial D 2†r0 v0 e3 C 4†r0 v0 e3 9
D 6†r0 v0 e3 :

[1] G.K. Batchelor, (1967) An Introduction to fluid dynamics, Cambridge University Press. (p.229)
[2] [Link], Coltrin, Glarbor, (2003) Chemically reacting flow: theory and practice, John Wiley & Sons Co.
(apéndice B)
[3] U.H. Kurzweg, “Derivation of the Stokes drag formula”, University of Florida, Department of Mechanical
and Aerospace Engineering, noviembre (2003).

View publication stats