UNIVERSIDAD EIA

Fenómenos de Transporte
Perfiles de velocidad y esfuerzos viscosos

Objetivo
Aplicar el estudio realizado en la clase del cálculo de perfiles de velocidad y esfuerzos viscosos de fluidos simples.

Conocimientos previos
Para resolver situaciones que involucran cálculo de perfiles de velocidad y esfuerzos viscosos de fluidos, el
estudiante debe tener solvencia en la manipulación de:

• la ecuación de continuidad

∂ρ
= − (∇ · ρv) (1)
∂t
• la ecuación de la cantidad de movimiento

∂ ∂ ∂ ∂
ρvi δ i = − ρvj vi δ i − τji δ i − pδ i + ρgi δ i (2)
∂t ∂xj ∂xj ∂xi

• Condiciones lı́mite (Condiciones de frontera)

1. En la interfaz sólido-lı́quido, la velocidad del fluido es igual a la velocidad con la que se mueve el
sólido. Se aplica a la componente tangencial y normal de las velocidad. La igualdad de las compo-
nentes tangenciales se denomina ”condición sin deslizamiento”.

2. En el plano interfacial lı́quido-lı́quido de x constante, las componentes tangenciales de velocidad vy

y vz son continuas a través de la interfase (Condición sin deslizamiento) ası́ como las componentes
de esfuerzo molecular.

3. En un plano interfacial lı́quido-gas de x constante, las componentes del tensor de esfuerzo τxy y τxz
se toman como cero, siempre que el gradiente de velocidad del lado del gas no sea demasiado grande.
Esto es razonable, ya que las viscosidades de los lı́quidos son mucho mayores que las de los gases.

Repaso algunos conocimientos previos

Para tener éxito en el desarrollo de la actividad que se propone es fundamental que el estudiante haga un
repaso profundo de los perfiles de velocidad y esfuerzos viscosos que se calcularon en la clase.

Planteamiento del ejercicio

Considere un fluido que se mueve entre dos placas, como se muestra en la figura 1. El fluido tiene una viscosidad
µ = 0.8 Pa s y conductividad térmica kf = 0.145W/mK . La placa inferior se mueve hacia la izquierda con
una velocidad v1 = 5m/s y la placa superior se mueve con una velocidad uniforme v2 = 3v1 . La distancia
entre las placas es h = 5mm . Suponga la conductividad térmica y la viscosidad constantes en todo momento
y además que el fluido es incompresible. Las placas tienen un ancho (perpendicular a la hoja) w = 10cm y un
largo l = 50cm. Nota: Todos los desarrollos se deben hacer con el sistema de referencia que se
muestra en la figura 1.

1. Para cuando se alcance el estado estacionario, determinar los perfiles de la velocidad y de los esfuerzos
viscosos del fluido.

2. Determinar la fuerza necesaria para mantener la rapidez de la tabla superior.

Fluido entre las placas

Figure 1: Fluido entre un par de láminas horizontales.

Solución:

1. Como el fluido es incompresible la ecuación de continuidad queda:

0 = ∇·v
0 0
∂vx ∂vy
7 ∂vz7
⇒0 = +  + 
∂x ∂y ∂z
| {z  }
se mueve en x
∂vx
∴0 = (3)
∂x

Ahora, en la ecuación de cantidad de movimiento,

0
∂  ∂ ∂ ∂p
ρvi =− ρvj vi − τji − + ρgi (4)
|∂t{z } ∂xj ∂xj ∂xi
Es.estacionario

Haciendo i = x y j = x, y, z la ecuación (4) toma la forma:

 
∂  *0 ∂ ∂ 0
 0
0=−  x x + ρ
ρv v v
y vx +
> ρ
v
z vx 
> 
|∂x {z } ∂y ∂z
 
Ec.continuidad
 
∂  0 ∂ ∂ * 0−
∂p 0
− τxx
*+ τyx +  τzx
 + ρ
g
x
>
∂x  ∂y ∂z ∂x
 
∂ ∂vx ∂p
⇒0=− −µ −
∂y ∂y ∂x

∂2 ∂p
∴0=µ vx − (5)
∂y 2 ∂x
 Haciendo i = y y j = x, y, z la ecuación (4) toma la forma:

 
0
 ∂ *
 ∂ 0 ∂ 0
0=− x
ρv v y + ρ
v
y vy +
> ρ
v
z vy 
> 
  ∂x
| {z } ∂y ∂z
Ec.continuidad
 
∂ ∂ 0 ∂ ∂p
− τxy +  τyy
*+ * 0−
τzy + ρgy
∂x ∂y ∂z  ∂y
    0
∂ ∂vx ∂p ∂ ∂v
x
>
 ∂p
⇒0=− −µ − + ρgy = µ − + ρgy
∂x ∂y ∂y ∂y  ∂x ∂y
| {z }
Ec.continuidad

∂p
∴ = ρgy (6)
∂y
Al resolver la ecuación diferencial (6) se obtiene:

p = ρgy y + C1 (7)
lo cual indica que la presión solo depende de y, por ende la ecuación (5) queda:

∂2
0=µ vx
∂y 2

⇒ v x = C2 y + C3 (8)

Para determinar las constantes de integración C2 y C3 se utilizan las condiciones de frontera de Dirichlet,
luego:

Para y = 0; vx = −v1 ⇒ C3 = −v1

4v1
Para y = h; vx = v2 ⇒ C2 = h
 
4v1
∴ vx = y − v1 (9)
h

Por otra parte, se sabe que:

   
∂vx 4v1
τxy = τyx = −µ ∴ τyx = −µ
∂y h

2. Para determinar la fuerza necesaria para mantener la rapidez de la tabla superior se usa el perfil de
cizalladura, esto es:
 
Fx 4v1
τyx = ⇒ Fx = Aτyx ∴ |Fx | = µA
A h
Reemplazando los valores ,

4 × 0.8Pa s × 5m/s × 0.05m2

 
|Fx | = = 160N
0.005m