Una masa semi infinita de un líquido de 𝜌 y 𝜇 constantes está limitada, por un lado, mediante una superficie plana (el
plano
𝑥𝑧 ) inicialmente, tanto el fluido como la superficie sólida están en reposo; pero en el instante 𝑡 = 0 la superficie solida se
pone en movimiento en la dirección x positiva con una velocidad 𝑉, tal como se indica en la Figura. Se desea conocer la
velocidad en función de 𝑦 y 𝑡 . En la dirección x no existe gradiente de presión ni fuerza de gravedad y se supone que el
flujo es laminar.
t<0
Fluido en reposo
(Placa inferior estatica)
y
x
t=0
Fluido en reposo
(Placa inferior inicia movimiento
y con velocidad constante V)
x V
t > 0 (Pequeño)
Vx (y,t) Fluido en movimiento
(Formación de flujo no estacionario, perfil
y de velocidad dependiente del tiempo)
x V
Vx (y) t > 0 (Grande)
Fluido en movimiento
(Formación de flujo estacionario, perfil de
y velocidad independiente del tiempo)
x V
Fig. 1 Flujo en las proximidades de una pared que se pone bruscamente en movimiento
Presunciones para el ejercicio:
1. Estado no estacionario.
2. Fluido incompresible con 𝜌 y 𝜇 constante.
3. Flujo laminar.
4. VY = VZ = 0
5. VX = VX (𝑦, 𝑡)
6. 𝑃𝑋 = 𝑝𝑌 = 𝑃𝑍 = 0
7. Sistema horizontal. Gravedad no afecta.
A partir de la ecuación de continuidad en coordenadas rectangulares:
𝜕𝜌 𝜕 𝜕 𝜕
+ (𝜌𝑣𝑥 ) + (𝜌𝑣𝑦 ) + (𝜌𝑣𝑧 ) = 0
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
2 5 4 4
0=0
Con el uso de las ecuaciones de Navier-Stokes en términos de la velocidad para la coordenada x que es la dirección en la
que se efectúa el movimiento tenemos.
𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑃 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑥
𝜌( + 𝑣𝑥 + 𝑣𝑦 + 𝑣𝑧 )=− +𝜇( 2 + + ) + 𝜌𝑔𝑥
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑥 𝜕𝑦 2 𝜕𝑧 2
5 4 4 6 5 5 7
𝜕𝑣𝑥 𝜕 2 𝑣𝑥
=𝜈
𝜕𝑡 𝜕𝑦 2
Donde 𝜈 = 𝜇⁄𝜌
Obtenemos una ecuación diferencial parcial que debe cumplir las siguientes dos condiciones de frontera (segundo orden
para el eje x), y condición inicial (primer orden para la variable t):
Codición Inicial: para 𝑡 = 0, 𝛿(𝑡) = 0 para todo valor de 𝑦
Codición de frontera #1: para 𝑦 = 0, 𝑣𝑥 = 𝑉 para todo valor de 𝑡 > 0
Codición de frontera #2: para 𝑦 = 𝛿, 𝑣𝑥 = 0 para todo valor de 𝑡 > 0
δ(t3)
Vx (y,t)
Vx (y,t) δ(t2)
t creciente t creciente
δ(t1)
y y
x V x V
Fig. 2a. Solución verdadera admitida Fig. 2b. Aproximación de la capa limite en función del
tiempo
Para dar solución a la anterior ecuación diferencial parcial, aplicaremos un procedimiento aproximado, en el cual se prevé
que la verdadera solución de problema será algo semejante a las curvas representadas en la Fig. 2a, pero supondremos que
la solución puede representarse satisfactoriamente mediante el gráfico de la Fig. 2b. Es decir, suponemos que en cualquier
instante 𝑡, hay un espesor 𝛿(𝑡) de capa límite, más allá del cual no existe flujo, de forma que los efectos de la viscosidad se
reducen a la región 0 ≤ 𝑦 ≤ 𝛿(𝑡). Se admite, finalmente, que los perfiles de velocidad adimensional son semejantes a lo
largo del tiempo; esta idea se expresa analíticamente de la siguiente forma
1 𝜕𝑣𝑥 1 𝜕 2 𝑣𝑥
=𝜈
𝑉 𝜕𝑡 𝑉 𝜕𝑦 2
Puesto que 𝑉 es constante, puede entrar dentro del operador derivada.
𝜕(𝑣𝑥 /𝑉) 𝜕 2 (𝑣𝑥 /𝑉)
=𝜈 (1)
𝜕𝑡 𝜕𝑦 2
Con 𝜙(𝜂) = 𝑣𝑥 ⁄𝑉 (Función 𝜙 con variable independiente adimensional 𝜂)
𝜕(𝑣𝑥 /𝑉) 𝜕𝜙 𝜕𝜂
= (2)
𝜕𝑡 𝜕𝜂 𝜕𝑡
En el anterior paso, se aplico regla de la cadena.
Ahora, si 𝜂 = 𝑦⁄𝛿(𝑡) tenemos que:
𝜕𝜂 𝜕 𝑦 𝜕 1 1 𝑑𝛿
= ( ) = 𝑦 ( ) = 𝑦 (− 2 )
𝜕𝑡 𝜕𝑡 𝛿 𝜕𝑡 𝛿 𝛿 𝑑𝑡
𝑦 no depende de 𝑡 Derivada de 1/𝛿 con respecto a 𝑡, y
(Sale del operador derivada) aplicación de regla de la cadena
Lo cual puede ser reescrito de la siguiente manera:
𝜕𝜂 𝑦 1 𝑑𝛿 𝜂 𝑑𝛿
=− =− (3)
𝜕𝑡 𝛿 𝛿 𝑑𝑡 𝛿 𝑑𝑡
𝑦
=𝜂
𝛿
Remplazando (3) en (2)
𝜕(𝑣𝑥 /𝑉) 𝜂 𝑑𝛿 𝜕𝜙
=− (4)
𝜕𝑡 𝛿 𝑑𝑡 𝜕𝜂
De la misma manera con 𝜙(𝜂) = 𝑣𝑥 ⁄𝑉 (Función 𝜙 con variable independiente adimensional 𝜂)
𝜕 2 (𝑣𝑥 /𝑉) 𝜕 𝜕(𝑣𝑥 /𝑉) 𝜕 𝜕𝜙 𝜕𝜂 𝜕 𝜕𝜙 𝜕𝜂 𝜕 𝜕𝜂 𝜕𝜙 𝜕 2 𝜙 𝜕𝜂 𝜕𝜂 𝜕 2 𝜂 𝜕𝜙
2
= ( )= ( ) = ( ( )) + ( ( )) = ( 2 ( )) + ( 2)
𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝜂 𝜕𝑦 𝜕𝑦 𝜕𝜂 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝜂 𝜕𝜂 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝜂
Regla de la Derivada de un Regla de la
cadena producto cadena
Con lo cual:
𝜕 2 (𝑣𝑥 /𝑉) 𝜕 2 𝜙 𝜕𝜂 2 𝜕 2 𝜂 𝜕𝜙
2
= 2 ( ) + ( 2) (5)
𝜕𝑦 𝜕𝜂 𝜕𝑦 𝜕𝑦 𝜕𝜂
Si 𝜂 = 𝑦⁄𝛿(𝑡) tenemos que:
𝜕𝜂 𝜕 𝑦 1 𝜕 1 1
= ( ) = ( ) (𝑦) = (1) = (6)
𝜕𝑦 𝜕𝑦 𝛿 𝛿 𝜕𝑦 𝛿 𝛿
𝛿 no depende de 𝑦 Derivada de 𝑦 con respecto a 𝑦
(Sale del operador derivada)
Y la segunda derivada teniendo en cuenta la ecuación (6)
𝜕2𝜂 𝜕 𝜕𝜂 𝜕 1
= ( )= ( )=0 (7)
𝜕𝑦 2 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝛿
Termino que no depende de 𝑦
(derivada de una constante)
Remplazando (6) y (7) en (5)
𝜕 2 (𝑣𝑥 /𝑉) 𝜕 2 𝜙 1 2 𝜕𝜙 𝜕 2 𝜙 1 2
2 = 2 ( ) + (0) = ( )
𝜕𝑦 𝜕𝜂 𝛿 𝜕𝜂 𝜕𝜂2 𝛿
Con lo cual se obtiene
𝜕 2 (𝑣𝑥 /𝑉) 1 𝜕2𝜙
2
= 2 2 (8)
𝜕𝑦 𝛿 𝜕𝜂
Remplazando (8) y (4) en (1)
𝜂 𝑑𝛿 𝜕𝜙 1 𝜕2𝜙
− =𝜈 2 2
𝛿 𝑑𝑡 𝜕𝜂 𝛿 𝜕𝜂
𝜕𝜙 𝑑𝛿 𝜕2𝜙
−𝛿𝜂 =𝜈 2
𝜕𝜂 𝑑𝑡 𝜕𝜂
En esta ecuación 𝜙 es solo función de 𝜂, por lo cual las derivadas parciales se convierten en derivadas totales.
𝑑𝜙 𝑑𝛿 𝑑2 𝜙
−𝜂 𝛿 =𝜈 2 (9)
𝑑𝜂 𝑑𝑡 𝑑𝜂
Con las siguientes condiciones de frontera:
Codición de frontera #1: para 𝑦 = 0 ⇒ η = 0, 𝑣𝑥 = 𝑉 ⇒ 𝜙 = 1, para todo valor de 𝑡 > 0
Codición de frontera #2: para 𝑦 = δ ⇒ η = 1, 𝑣𝑥 = 0 ⇒ 𝜙 = 0, para todo valor de 𝑡 > 0
Multiplicando toda la ecuación (9) por 𝑑𝜂
𝑑𝜙 𝑑𝛿 𝑑2 𝜙
−𝜂 𝑑𝜂𝛿 = 𝜈 2 𝑑𝜂
𝑑𝜂 𝑑𝑡 𝑑𝜂
Integrando respecto a 𝜂
1 1 2
𝑑𝜙 𝑑𝛿 𝑑 𝜙
(− ∫ 𝜂 𝑑𝜂 ) 𝛿 = 𝜈 (∫ 2 𝑑𝜂 )
(10)
0 𝑑𝜂 𝑑𝑡 𝑜 𝑑𝜂
Teniendo en cuenta que:
1 1 1 1
𝑑𝜙
−∫ 𝜂 𝑑𝜂 = − [(𝜂𝜙)|10 − ∫ 𝜙𝑑𝜂 ] = − [(1)(0) − (0)(1) − ∫ 𝜙𝑑𝜂 ] = + ∫ 𝜙𝑑𝜂 = 𝑀
0 𝑑𝜂 0 0 0
Integración por partes, con: Evaluación del producto
teniendo en cuenta las
𝑢 = 𝜂 ⇒ 𝑑𝑢 = 𝑑𝜂
condiciones limites
𝑑𝜙
𝑑𝑣 = 𝑑𝜂 ⇒ 𝑣 = 𝜙
𝑑𝜂
1 1
𝑑𝜙
𝑀 = −∫ 𝜂 𝑑𝜂 = + ∫ 𝜙𝑑𝜂 (11)
0 𝑑𝜂 0
Además
1
𝑑2 𝜙 𝑑𝜙 1
𝑁=∫ 𝑑𝜂 = | (12)
𝑜 𝑑𝜂2 𝑑𝜂 0
Remplazando (11) y (12) en (10)
𝑑𝛿
𝑀𝛿 = 𝜈𝑁 (13)
𝑑𝑡
Separando variables.
𝑀𝛿𝑑𝛿 = 𝜈𝑁𝑑𝑡
𝛿2
𝑀 = 𝜈𝑁𝑡 + 𝐶1
2
Codición Inicial: para 𝑡 = 0, 𝛿(𝑡) = 0 para todo valor de 𝑦
Por lo cual se obtiene que 𝐶1 = 0, Remplazando 𝐶1 y despejando 𝛿
𝛿(𝑡) = √2(𝑁/𝑀)𝜈𝑡 (14)
De esta manera podemos determinar una ecuación para el espesor de la capa límite para cualquiera función 𝜙(𝜂).
Elegimos ahora una función específica, cualquiera que sea que cumpla las condiciones de frontera.
Codición de frontera #1: para 𝑦 = 0 ⇒ η = 0, 𝑣𝑥 = 𝑉 ⇒ 𝜙 = 1, para todo valor de 𝑡 > 0
Codición de frontera #2: para 𝑦 = δ ⇒ η = 1, 𝑣𝑥 = 0 ⇒ 𝜙 = 0, para todo valor de 𝑡 > 0
Por ejemplo
3 1
𝜙(𝜂) = 1 − 𝜂 + 𝜂3 (15)
2 2
Entonces:
1 1 1
𝑑𝜙 3 1 3
𝑀 = −∫ 𝜂 𝑑𝜂 = + ∫ 𝜙𝑑𝜂 = ∫ (1 − 𝜂 + 𝜂3 ) 𝑑𝜂 =
0 𝑑𝜂 0 0 2 2 8
1
𝑑2 𝜙 𝑑𝜙 1 3 3 3
𝑁=∫ 𝑑𝜂 = | = (− + 𝜂2 ) |10 =
𝑜 𝑑𝜂2 𝑑𝜂 0 2 2 2
Remplazando valores de 𝑁 y 𝑀 en (14)
3/2
𝛿(𝑡) = √2 ( ) 𝜈𝑡 = 2√2𝜈𝑡 (16)
3/8
Por lo cual.
𝑦 𝑦
𝜂= =
𝛿(𝑡) 2√2𝜈𝑡
Remplazando 𝜂 en (15)
𝑣𝑥 3 𝑦 1 𝑦 3
𝜙(𝜂) = =1− ( )+ ( ) (17)
𝑉 2 2√2𝜈𝑡 2 2√2𝜈𝑡
Este resultado a partir del método aproximado, se compara con la solución exacta, cuyo desarrollo en serie, es:
∞
𝑣𝑥 2 𝑦 2 23 𝑦 3 2(−1)𝑛 22𝑛−1 𝑦 2𝑛−1
= 1 − 2√ ( )+ √ ( ) +∑ √ ( ) (18)
𝑉 𝜋 2√2𝜈𝑡 3 𝜋 2√2𝜈𝑡 2𝑛 − 1 𝜋 2√2𝜈𝑡
𝑛=3
Por lo tanto, la función 𝜙(𝑛) seleccionada conduce a los 3 primeros términos de la solución exacta, con modificaciones en
la elección de las constantes para una mejor aproximación.
